Search NASA⌕ Search

Engineering topics

Carlos G Davila

Publications and source records attributed to Carlos G Davila.

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↗

TPSAS-NF1676L-15315-DND

Some advantages over mechanical joints: - Lower structural weight - Lower manufacturing cost - Improve damage tolerance However: - Lack of failure criteria and material/structural validated models - Lack of an accurate and reliable quality control

Donato Girolamo↗

TPSAS-NF1676L-25778-DND

Testing required to design and certify aircraft structures has increased dramatically with the growth in use of composite materials, and the development time and cost have become an issue as industry looks towards the next generation of aircraft. To help address this challenge, NASA has formed an industry-spanning consortium with the goal of cutting by 30% today's 10-year timeline for developing and certifying composite aircraft structures. This goal is to be achieved by infusing leading edge, physics-based tools to improve the prediction of the response and failure of composite aerospace structures. The plan for the public-private partnership of NASA, the FAA, airframe manufacturers, aircraft engine manufacturers, and several universities is to conduct joint research under the five-year Advanced Composites Program (ACP). The emphasis of the ACP project is the development of damage prediction and non-destructive inspection methods capabilities that are sufficiently reliable to enable a reduction in the element and subcomponent testing necessary for design and certification. ACP also is expected to create rapid design tools to speed preliminary design and improve safety and performance. A critical element of this research concerns the residual strength of postbuckled stiffened panels with impact damage. This presentation will discuss different aspects of ongoing research on the prediction of structural collapse due to skin/stiffener separation. In particular, the presentation will examine issues that must be addressed in the analysis, such as R-curve toughening, delamination migration, bridging, as well as some unexpected effects of ply orientation.

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↗

Scalability of Cohesive Fatigue Analyses Using Explicit Solvers

A cohesive fatigue law has been integrated into a constitutive material model compatible with an explicit finite element solver. The cohesive fatigue model response is based on engineering approximations of the endurance limit and the Goodman diagram. This approach can predict stress-life diagrams for crack initiation, the Paris law regime, and transient effects of crack initiation and stable tearing. Simplified cyclic loading is utilized so that the applied load(or displacement) corresponds to the peak load of a fatigue cycle. Loads are held constant during fatigue while damage develops with increasing solution increments. An automatically-calculated ratio of fatigue cycles per solution increment controls the rate of damage growth, ensuring that damage growth is modeled with a sufficient minimum number of increments and damage growth advances to a minimum desired extent within the explicit analysis step time. The compatibility with an explicit finite element solver enables the analysis of structures that are computationally intractable for implicit finite element solvers. Scalability studies are conducted for geometrically nonlinear problems involving fiber-reinforced composite structures that exhibit fatigue damage growth of interacting matrix cracks and delaminations

Frank A Leone↗

Modal Test and Analysis Correlation of Wind Tunnel Blades for the National Transonic Facility

NASA has a significant maintenance challenge with the aging wind tunnel infrastructure across the agency . Wind tunnel blades, such as the ones in the National Transonic Facility at NASA Langley Research Center, have been subjected to damage from impacts and fatigue. Accurate analytical models of these blades are necessary to either recertify replacement blades to be built or for redesign of these blades. In this paper an upgraded finite element model of these wind tunnel blades with the capability of representing damage is presented. An eigenvalue analysis of the structural model is performed, and the boundary conditions are adjusted to simulate the natural frequencies obtained by experimental testing of a set of twenty five blades in 1992. Additional numerical studies of the effect of delamination damage on the natural frequencies of the blade are conducted using the tuned finite element model.

wind tunnel↗

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↗

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↗

Softening the Gap between Wöhler and Paris – New Approaches for Fatigue Analysis –

Fatigue analysis tools can vary across industries. For example, automotive engineers often use the Wöhler (S-N) approach to design for safe-life, while aerospace engineers prioritize damage tolerance and inspection intervals, relying instead on crack growth models such as Paris’ law. Although both approaches may deal with the control of cracks in similar materials, their analysis tools and material characterizations are fundamentally distinct. This divide mirrors the classic split between stress-based strength analysis and linear elastic fracture mechanics. However, modern nonlinear models that incorporate material softening, such as cohesive laws, blur this boundary and capture fracture behaviors across scales. This presentation describes the CF23 fatigue model, which uses cohesive softening to link S-N crack initiation with crack propagation rates. CF23 spans the full fatigue spectrum, from initial propagation transients to steady-state growth and threshold conditions, offering a unified framework that bridges Wöhler and Paris-based methodologies. Example applications include fatigue crack propagation transients in adhesive interfaces and skin/stiffener separation.

cohesive elements↗