Search NASA⌕ Search

Engineering topics

Cheryl A Rose

Publications and source records attributed to Cheryl A Rose.

Progressive Damage Analysis of Post-Buckled Stiffened Panels under Static Compressive Loading

The validation of progressive damage analysis methods requires experimental data that captures the initiation of damage, its evolution with continued loading, and the morphology of damage throughout its evolution. In support of the NASA Advanced Composites Project, several validation specimens of increasing complexity have been tested and characterized in order to provide the data required to assess the capabilities of several progressive damage analysis methods for fiber-reinforced composite materials. The analytical capabilities of interest for this validation program include networks of interacting matrix cracks and delaminations in post-buckled stiffened structures. A series of single stringer compression specimens was tested and characterized at NASA Langley Research Center. The experimental data from these tests are compared with progressive damage finite element analysis results. Test/analysis comparisons are made in terms of load/displacement histories, damage morphology, and damage size versus load metrics.

Frank A Leone↗

TPSAS-NF1676L-32227-DND

The ability of a recently proposed cohesive fatigue model to predict the initiation and propagation of skin/stiffener separation was evaluated using a three-point bend (3PB) configuration. The 3PB specimen consists of a quasi-isotropic laminate representing the skin of a fuselage. A shorter doubler representing the flange of a stiffener is co-cured at the center of the specimen. A simplified modeling technique was applied, in which the complex damage mechanisms involved in the separation, such as delamination migration, matrix cracking, bridging, and crack delving are represented by a mixed-mode trilinear cohesive law at the interface between the skin and the doubler. The quasi-static properties of the trilinear law were characterized. The main advantage of the proposed fatigue model is that it relies on quasi-static properties as well as engineering assumptions to compose the shape of the S-N diagram. The model does not require Paris law data, and it is capable of predicting the initiation of separation, the rate of propagation, as well as the effects of crack resistance curves and the stress ratio.

Carlos G Davila↗

Initiation and Propagation of Skin/Stiffener Separation in Postbuckled Structures Subjected to Cyclic Loads

A methodology is proposed to evaluate experimentally and numerically the initiation and propagation of skin/stiffener separation in postbuckled panels subjected to cyclic loads. The response mechanisms that characterize skin/stiffener separation were investigated using a three point bend (3PB) test specimen, which is a simple specimen consisting of a laminate representing the skin and a doubler representing the stiffener flange. The 3PB tests were performed under quasi-static and cyclic loads. A fatigue cohesive model was used to analyze the initiation and propagation of skin/stiffener separation in the 3PB. The results indicate that the fatigue cohesive model with model parameters obtained from quasi-static fracture tests can provide accurate predictions of initiation from a pristine state, and the rate of fatigue crack propagation for a variety of mode mixities, load levels, and stress ratios.

Carlos G Davila↗

Experimental Observations of Damage States in Unnotched and Notched 3D Orthogonal Woven Coupons Loaded in Tension

This paper describes a series of unnotched and notched tensile tests on two 3D orthogonal woven architectures. Dog bone specimens were used for the unnotched tests, and open-hole specimens were used for the notched tests. For each architecture and test type, two specimens were loaded under monotonic displacement control to failure and two specimens were incrementally loaded to failure in five steps. X-ray computed tomography (CT) was used to characterize the woven architecture and defects on each specimen prior to testing. Digital image correlation (DIC) images were obtained from both sides of the specimen throughout the loading and unloading. Acoustic emission (AE) data were also continuously acquired throughout the loading. For the incremental tests, X-ray CT was used to document the damage progression after each load increment. The damage progression was analyzed based on the results from these three techniques. The results for both architectures and test specimens followed a pattern where matrix cracks formed at the surface perpendicular to the tensile load direction. New surface cracks formed continuously throughout the loading and cracks extended in depth and width. In general, fiber failure was not detected until specimen failure.

Wade C Jackson↗

Challenges in Scaling Analysis from Coupons to Elements in 3D Woven Composites for a Launch Vehicle End-Ring Structure

The limitations of structural design with 3D woven composites based on homogenized material models characterized using coupon test data, representative unit cell analysis, or a combination of the two are discussed in the context of a space launch system (SLS)-scale end ring. Selected building block test configurations at the coupon, element, and subcomponent scales are described in detail. Highly instrumented tests and X-ray computed tomography inspections are used to characterize the structural behavior and associated damage progression with increasing load. It is shown that a single homogenized material constitutive relationship, applied at coupon and element scale, is not sufficient to model the structural response observed in both tests. In fact, substantially different stress-strain laws are required to predict the response of the coupon and element scale specimens, despite similar loading conditions. A model with discrete representation of the binder tows and cohesive interfaces to model matrix cracking shows much better agreement with the test data at the two length scales. Thus, evidence is provided showing that predictive structural analysis of 3D woven composites requires discrete representation of the tows, matrix cracking morphology, and their interactions in the structural configuration. With accurate knowledge of the element scale behavior, engineering approximations readily enable analysis scale-up to the subcomponent level.

Andrew C Bergan↗

Experimental Observations of Damage States in Unnotched and Notched 3D Orthogonal Woven Coupons Loaded in Tension

This paper describes a series of unnotched and notched tensile tests on two 3D orthogonal woven architectures loaded in the warp direction. Dogbone specimens were used for the unnotched tests, and open-hole specimens were used for the notched tests. For each architecture and test type, two specimens were loaded under monotonic displacement control to failure, and two specimens were incrementally loaded to failure in five steps. X-ray computed tomography (CT) was used to characterize the woven architecture and defects on each specimen prior to testing. Digital image correlation (DIC) images were obtained from both sides of the specimen throughout the loading and unloading. For the incremental tests, X-ray CT was used to document the damage progression after each load increment. The results for both architectures and test specimens followed a similar pattern where matrix cracks formed perpendicular to the load direction and were primarily located between the stacks of warp tows. Two-piece failure path always occurred perpendicular to the loading direction, along the through-thickness segments of the Z-tows. The combination of stress concentrations due to the specimen geometry, weave geometry, and cut tows created stress concentrations that appear to initiate failure.

Woven materials↗

Interlaminar R-Curve Effects in Thermoplastic Composites: Experiments and Model Correlations

The interlaminar fracture behavior of thermoplastic composites under Mode-I and Mode-II static loading is investigated in this research. The material system under consideration is a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK)†. Characterization experiments were conducted to measure the load-displacement response, peak load, and propagation of delamination from a pre-implanted Kapton® insert. The Mode-I response was characterized using a double cantilever beam (DCB) test according to the ASTM D5528-22 standard. It was observed that the propagation exhibited significant resistance to delamination growth, i.e., an R-Curve effect. The test data indicated that the steady-state fracture toughness, G SS I , was approximately twice the mean initiation fracture toughness, G Ic . The Mode-II response of an end notch flexure (ENF) specimen tested according to the ASTM D7905-19 standard exhibited an unexpected stable behavior upon attaining peak load. Prior to attaining peak load, the Mode-II response exhibited some non-linearity, which may be attributed to nonlinear bulk material response, e.g., hyper-elasticity, viscoelasticity, plasticity, etc., and potentially friction at the load/support rollers and R-Curve effects. Finite element models based on the Cohesive Zone Method (CZM) and the Virtual Crack Closure Technique (VCCT) were used to predict the DCB and ENF specimen responses. DCB analysis results obtained using G SS I predicted propagation consistent with experimental data. However, results obtained with the as characterized G Ic were inaccurate in predicting delamination growth onset and peak load. ENF analysis results obtained with the as-characterized Mode-II fracture toughness, G IIc , were inaccurate in predicting the peak load and propagation. Hence, R-Curve effects were included in both DCB and ENF simulations, which enabled improved predictions. The experimentally determined Mode-I R-Curve was used for the DCB simulations, whereas an R-Curve was derived for Mode-II assuming that the change in compliance of the specimen was associated with crack growth only.

Prabhakar M Rao↗

Low-Velocity Impact Analyses of Composite Panels Using Cohesive Zone Modeling and Continuum Damage Mechanics

A parametric study was conducted to gage the sensitivity of simulations of low-velocity impacts on laminated composite panels to various finite element model and material model features. The material model used in this parametric study was previously verified and validated for conducting progressive damage analyses of post-buckled composite structures subjected to quasi-static loading conditions. Based on the results of this parametric study, updates have been made to a progressive damage analysis material model and the best practices for its application. Experimental results are compared to computational results obtained with and without the updates to the models. Both thermoset and thermoplastic fiber-reinforced materials are considered. The resulting enhancements to the material model are presented in detail.

Frank A Leone↗

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗

Analysis of Crack Propagation Transients in Accelerated Fatigue Characterization of Adhesives

The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.

fracture mechanics; delamination; J-integral↗

Calculating Mode Mixity in Interlaminar Fracture of Composite Laminates

To accurately predict delamination under mixed-mode loading, cohesive elements must partition the energies dissipated in Mode I and Mode II correctly. The ability of cohesive models to predict mode mixity was evaluated by performing analyses of the single leg bending (SLB) test, which exhibits a variable mode mixity with crack length. Accordingly, two numerical integration methods were developed to compute the energy partitioning: (i) a local calculation of the accumulated work of tractions at each point along the crack path, and (ii) a global J-integral “snapshot” calculation along a contour around the process zone. In the case of a composite laminate with properties characterized by short fracture process zone (FPZ), the results of these two methods correlate well with each other and with a linear elastic fracture mechanics solution. However, FPZscan be long due to fiber bridging or due to the intrinsic toughness of the material. Then, the results of mode-mixity calculations obtained with various methods vary widely from each other. It is argued that the global J-integral calculation is the most accurate interpretation of the energy release rate and mode mixity for either short or long process zones. This solution will provide a more accurate characterization of resistance curves and mode mixity in experimental results, and it will also guide the development of more accurate mixed-mode cohesive laws for tougher materials characterized by strong R-curve responses, such as thermo-plastic composites.

cohesive zone analysis↗