Search NASA⌕ Search

SEARCH · Search NASA

Results for “energy release rate”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Dynamic delamination buckling in composite laminates under impact loading: Computational simulation

A unique dynamic delamination buckling and delamination propagation analysis capability has been developed and incorporated into a finite element computer program. This capability consists of the following: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading, and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. This computer program is used to predict both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. It is shown that delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate.

Grady, Joseph E.↗

Dynamic delamination buckling in composite laminates under impact loading - Computational simulation

A unique dynamic delamination buckling and delamination propagation analysis capability has been developed and incorporated into a finite element computer program. This capability consists of the following: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading, and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. This computer program is used to predict both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. It is shown that delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate.

Grady, Joseph E.↗

A theory for the fracture of thin plates subjected to bending and twisting moments

Stress fields near the tip of a through crack in an elastic plate under bending and twisting moments are reviewed assuming both Kirchhoff and Reissner plate theories. The crack tip displacement and rotation fields based on the Reissner theory are calculated. These results are used to calculate the J-integral (energy release rate) for both Kirchhoff and Reissner plate theories. Invoking Simmonds and Duva's (1981) result that the value of the J-integral based on either theory is the same for thin plates, a universal relationship between the Kirchhoff theory stress intensity factors and the Reissner theory stress intensity factors is obtained for thin plates. Calculation of Kirchhoff theory stress intensity factors from finite elements based on energy release rate is illustrated. It is proposed that, for thin plates, fracture toughness and crack growth rates be correlated with the Kirchhoff theory stress intensity factors.

Hui, C. Y.↗

Analyses for Debonding of Stitched Composite Sandwich Structures Using Improved Constitutive Models

A fracture mechanics analysis based on strain energy release rates is used to study the effect of stitching in bonded sandwich beam configurations. Finite elements are used to model the configurations. The stitches were modeled as discrete nonlinear spring elements with a compliance determined by experiment. The constitutive models were developed using the results of flatwise tension tests from sandwich material rather than monolithic material. The analyses show that increasing stitch stiffness, stitch density and debond length decrease strain energy release rates for a fixed applied load.

Glaessgen, E. H.↗

Fracture-based shape optimization built upon the topological derivative

In Silva et al. (2011) and Alidoost et al. (2020), the authors developed an approximation of the energy release rate field associated with a small edge or surface crack at any boundary location and with any orientation using the topological derivative. The approximation is computationally attractive because it requires only a single analysis on the non-cracked domain in contrast with conventional boundary-element and finite-element-based methods, which require a separate and costlier analysis for each crack length-location-orientation combination. Here, a shape optimization scheme for fracture-resistant structures is developed using the energy release rate approximation. In the gradient-based optimization scheme, the domain and its boundary are defined implicitly using level-set functions. The level-set functions of arbitrary geometries are constructed using Boolean operations from the level-set functions of simple primitives. This geometrical representation has the dual advantage of (i) allowing shapes to intersect and/or separate during the optimization and (ii) simplifying the computation of the shape sensitivities.

42 ENGINEERING↗

Sizing Single Cantilever Beam Specimens for Characterizing Facesheet/Core Peel Debonding in Sandwich Structure

This paper details part of an effort focused on the development of a standardized facesheet/core peel debonding test procedure. The purpose of the test is to characterize facesheet/core peel in sandwich structure, accomplished through the measurement of the critical strain energy release rate associated with the debonding process. The specific test method selected for the standardized test procedure utilizes a single cantilever beam (SCB) specimen configuration. The objective of the current work is to develop a method for establishing SCB specimen dimensions. This is achieved by imposing specific limitations on specimen dimensions, with the objectives of promoting a linear elastic specimen response, and simplifying the data reduction method required for computing the critical strain energy release rate associated with debonding. The sizing method is also designed to be suitable for incorporation into a standardized test protocol. Preliminary application of the resulting sizing method yields practical specimen dimensions.

Ratcliffe, James G.↗

A Shell/3D Modeling Technique for the Analyses of Delaminated Composite Laminates

A shell/3D modeling technique was developed for which a local three-dimensional solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provided a kinematically compatible interface between the local three-dimensional model and the global structural model which has been meshed with plate or shell finite elements. Double Cantilever Beam (DCB), End Notched Flexure (ENF), and Single Leg Bending (SLB) specimens were modeled using the shell/3D technique to study the feasibility for pure mode I (DCB), mode II (ENF) and mixed mode I/II (SLB) cases. Mixed mode strain energy release rate distributions were computed across the width of the specimens using the virtual crack closure technique. Specimens with a unidirectional layup and with a multidirectional layup where the delamination is located between two non-zero degree plies were simulated. For a local three-dimensional model, extending to a minimum of about three specimen thicknesses on either side of the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential for reducing the model size, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗

A Shell/3D Modeling Technique for the Analysis of Delaminated Composite Laminates

A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a shell finite element model. Multi-point constraints provided a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with shell finite elements. Double Cantilever Beam, End Notched Flexure, and Single Leg Bending specimens were analyzed first using full 3D finite element models to obtain reference solutions. Mixed mode strain energy release rate distributions were computed using the virtual crack closure technique. The analyses were repeated using the shell/3D technique to study the feasibility for pure mode I, mode II and mixed mode I/II cases. Specimens with a unidirectional layup and with a multidirectional layup were simulated. For a local 3D model, extending to a minimum of about three specimen thicknesses on either side of the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures the shell/3D modeling technique offers a great potential for reducing the model size, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗

Modified Mode-I Cracked Sandwich Beam (CSB) Fracture Test

Five composite sandwich panels were fabricated using vacuum assisted resin transfer molding (VARTM). Four of these panels had E-glass/vinylester facesheets and one had carbon/epoxy facesheets. The sandwich panels had different density PVC foam cores. The four E-glass panels had core densities of 80, 100, 130, 200 kg/cu m. The sandwich with carbon/epoxy 3 facesheets had a core with density of 100 kg/cu m. Fracture tests were conducted using a modified Cracked Sandwich Beam (CSB) test configuration. Load displacement curves were obtained for loading and unloading of the specimens during crack growth. Various increments of crack growth were monitored. Critical Strain Energy Release Rates (SERR) were determined from the tests using the area method. The critical values of SERR can be considered the fracture toughness of the sandwich material. The fracture toughness ranged 367 J/sq m to 1350 J/sq m over the range of core densities. These results are compared to the Mode-I fracture toughness of the PVC foam core materials and values obtained for foam-cored sandwiches using the TSD specimen. Finite-element analyses (FEA) were performed for the test configuration and Strain Energy Release Rates were calculated using the Virtual Crack Closure Technique (VCCT). The SERR values determined from the FEA were scaled to the fracture loads, or critical loads, obtained from the modified CSB tests. These critical loads were in close agreement with the test values.

Smith, S. A.↗

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗

Finite element analysis of the end notched flexure specimen for measuring Mode II fracture toughness

The paper presents a finite element analysis of the end-notched flexure (ENF) test specimen for Mode II interlaminar fracture testing of composite materials. Virtual crack closure and compliance techniques employed to calculate strain energy release rates from linear elastic two-dimensional analysis indicate that the ENF specimen is a pure Mode II fracture test within the constraints of small deflection theory. Furthermore, the ENF fracture specimen is shown to be relatively insensitive to process-induced cracks, offset from the laminate midplane. Frictional effects are investigated by including the contact problem in the finite element model. A parametric study investigating the influence of delamination length, span, thickness, and material properties assessed the accuracy of beam theory expressions for compliance and strain energy release rate, GII. Finite element results indicate that data reduction schemes based upon beam theory underestimate GII by approximately 20-40 percent for typical unidirectional graphite fiber composite test specimen geometries. Consequently, an improved data reduction scheme is proposed.

Gillespie, J. W., Jr.↗

Characterization of the Edge Crack Torsion (ECT) Test for Mode III Fracture Toughness Measurement of Laminated Composites

The edge crack torsion (ECT) test is designed to initiate mode III delamination growth in composite laminates. An ECT specimen is a rectangular laminate, containing an edge delamination at the laminate mid-plane. Torsion load is applied to the specimens, resulting in relative transverse shear sliding of the delaminated faces. The test data reduction schemes are intended to yield initiation values of critical mode III strain energy release rate, G(sub IIIc), that are constant with delamination length. The test has undergone several design changes during its development. The objective of this paper was to determine the suitability of the current ECT test design as a mode III fracture test. To this end, ECT tests were conducted on specimens manufactured from IM7/8552 and specimens made from S2/8552 tape laminates. Several specimens, each with different delamination lengths are tested. Detailed, three-dimensional finite element analyses of the specimens were performed. The analysis results were used to calculate the distribution of mode I, mode II, and mode III strain energy release rate along the delamination front. The results indicated that mode III-dominated delamination growth would be initiated from the specimen center. However, in specimens of both material types, the measured values of G(sub IIIc) exhibited significant dependence on delamination length. Furthermore, there was a large amount of scatter in the data. Load-displacement response of the specimens exhibited significant deviation from linearity before specimen failure. X-radiographs of a sample of specimens revealed that damage was initiated in the specimens prior to failure. Further inspection of the failure surfaces is required to identify the damage and determine that mode III delamination is initiated in the specimens.

Ratcliffe, James G.↗

Mode 1 and Mode 2 Analysis of Graphite/Epoxy Composites Using Double Cantilever Beam and End-Notched Flexure Tests

The critical strain energy release rates associated with debonding of the adhesive bondlines in graphite/epoxy IM6/3501-6 interlaminar fracture specimens were investigated. Two panels were manufactured for this investigation; however, panel two was layed-up incorrectly. As a result, data collected from Panel Two serves no real purpose in this investigation. Double Cantilever Beam (DCB) specimens were used to determine the opening Mode I interlaminar fracture toughness, G1(sub c), of uni-directional fiber re-inforced composites. The five specimens tested from Panel One had an average value of 946.42J/sq m for G1(sub c) with an acceptable coefficient of variation. The critical strain energy release rate, G2(sub c), for initiation of delamination under inplane shear loading was investigated using the End-Notched Flexure (ENF) Test. Four specimens were tested from Panel One and an average value of 584.98J/sq m for G2(sub c) was calculated. Calculations from the DCB and ENF test results for Panel One represent typical values of G1(sub c) and G2(sub c) for the adhesive debonding in the material studied in this investigation.

Hufnagel, Kathleen P.↗

Analysis of Composite Skin-Stiffener Debond Specimens Using a Shell/3D Modeling Technique and Submodeling

The application of a shell/3D modeling technique for the simulation of skin/stringer debond in a specimen subjected to tension and three-point bending was studied. The global structure was modeled with shell elements. A local three-dimensional model, extending to about three specimen thicknesses on either side of the delamination front was used to model the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from shell/3D simulations were in good agreement with results obtained from full solid models. The good correlation of the results demonstrated the effectiveness of the shell/3D modeling technique for the investigation of skin/stiffener separation due to delamination in the adherents. In addition, the application of the submodeling technique for the simulation of skin/stringer debond was also studied. Global models made of shell elements and solid elements were studied. Solid elements were used for local submodels, which extended between three and six specimen thicknesses on either side of the delamination front to model the details of the damaged section. Computed total strain energy release rates and mixed-mode ratios obtained from the simulations using the submodeling technique were not in agreement with results obtained from full solid models.

OBrien, T. Kevin↗

Fracture toughness computational simulation of general delaminations in fiber composites

A procedure is described to computationally simulate composite laminate fracture toughness in terms of strain energy release rate. It is also used to evaluate the degradation in laminate structural integrity in terms of displacements, loss in stiffness, loss in vibration frequencies and loss in buckling resistance. Specific laminates are selected for detail studies in order to demonstrate the generality of the procedure. These laminates had center delaminations, off-center delaminations, and pocket delaminations (center and off-center) at the free-edge and center delaminations at the interior. The laminates had two different thicknesses and were made from three different materials. The results obtained are presented in graphical form to illustrate the effects of delamination on the laminate structural integrity and on the laminate strain energy release rate (composite fracture toughness).

Wilt, T. E.↗

Fracture toughness computational simulation of general delaminations in fiber composites

A procedure is described to computationally simulate composite laminate fracture toughness in terms of strain energy release rate (SERR). It is also used to evaluate the degradation in laminate structural integrity in terms of displacements, loss in stiffness, loss in vibration frequencies and loss in buckling resistance. Specific laminates are selected for detail studies in order to demonstrate the generality of the procedure. These laminates had center delaminations, off-center delaminations, and pocket delaminations (center and off-center) at the free-edge and center delaminations at the interior. The laminates had two different thicknesses and were made from three different materials. The results obtained are presented in graphical form to illustrate the effects of delamination on the laminate structural integrity and on the laminate strain energy release rate (composite fracture toughness).

Murthy, P. L. N.↗

Fracture Mechanics of Thin, Cracked Plates Under Tension, Bending and Out-of-Plane Shear Loading

Cracks in the skin of aircraft fuselages or other shell structures can be subjected to very complex stress states, resulting in mixed-mode fracture conditions. For example, a crack running along a stringer in a pressurized fuselage will be subject to the usual in-plane tension stresses (Mode-I) along with out-of-plane tearing stresses (Mode-III like). Crack growth and initiation in this case is correlated not only with the tensile or Mode-I stress intensity factor, K(sub I), but depends on a combination of parameters and on the history of crack growth. The stresses at the tip of a crack in a plate or shell are typically described in terms of either the small deflection Kirchhoff plate theory. However, real applications involve large deflections. We show, using the von-Karman theory, that the crack tip stress field derived on the basis of the small deflection theory is still valid for large deflections. We then give examples demonstrating the exact calculation of energy release rates and stress intensity factors for cracked plates loaded to large deflections. The crack tip fields calculated using the plate theories are an approximation to the actual three dimensional fields. Using three dimensional finite element analyses we have explored the relationship between the three dimensional elasticity theory and two dimensional plate theory results. The results show that for out-of-plane shear loading the three dimensional and Kirchhoff theory results coincide at distance greater than h/2 from the crack tip, where h/2 is the plate thickness. Inside this region, the distribution of stresses through the thickness can be very different from the plate theory predictions. We have also explored how the energy release rate varies as a function of crack length to plate thickness using the different theories. This is important in the implementation of fracture prediction methods using finite element analysis. Our experiments show that under certain conditions, during fatigue crack growth, the presence of out-of-plane shear loads induces a great deal of contact and friction on the crack surfaces, dramatically reducing crack growth rate. A series of experiments and a proposed computational approach for accounting for the friction is discussed.

Zehnder, Alan T.↗

A Shell/3D Modeling Technique for Delaminations in Composite Laminates

A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provide a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with plate or shell finite elements. For simple double cantilever beam (DCB), end notched flexure (ENF), and single leg bending (SLB) specimens, mixed mode energy release rate distributions were computed across the width from nonlinear finite element analyses using the virtual crack closure technique. The analyses served to test the accuracy of the shell/3D technique for the pure mode I case (DCB), mode II case (ENF) and a mixed mode I/II case (SLB). Specimens with a unidirectional layup where the delamination is located between two 0 plies, as well as a multidirectional layup where the delamination is located between two non-zero degree plies, were simulated. For a local 3D model extending to a minimum of about three specimen thicknesses in front of and behind the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Krueger, Ronald↗