Search NASA⌕ Search

Engineering topics

Prabhakar Rao

Publications and source records attributed to Prabhakar Rao.

Progressive Damage and Failure Analysis of Thermoplastic Composites in Low Velocity Impact Using MAT299

As part of the NASA Hi-Rate Composite Aircraft Manufacturing (HiCAM) Project, state-of-the-art progressive damage and failure analysis (PDFA) tools developed for use with thermosets are being evaluated for use in modeling alternative material systems, like thermoplastics. Experimental low-velocity impact data of a thermoplastic material system, AS4D/PEKK-FC, is presented and includes characterization of the impact damage mechanisms as well as associated load and displacement data. Following the presentation of experimental results, two simulation approaches using the PDFA tool MAT299 in the commercial off-the-shelf finite element software LS-DYNA are employed to predict damage area and force-displacement responses of thermoplastic panels subjected to various impact energies. The first modeling method uses solid elements with a high-density mesh and a ply-by-ply modeling approach similar to previously published work for thermosets. This method has the capability of capturing individual crack development, progression, and delamination on a per ply basis. The second modeling method uses TSHELL elements that have in-plane dimensions that are an order of magnitude larger than the elements used in the solid element approach and reduces the number of elements through the thickness of the laminate. The second method produces a lower-fidelity model with reduced run times that is incapable of monitoring every delamination plane and damage within each ply. Results from both simulation methods are compared to experiment, and limitations of the methods are discussed.

Composite material↗

Progressive Damage and Failure Analysis of Thermoplastic Composites in Low Velocity Impact Using MAT299

As part of the NASA Hi-Rate Composite Aircraft Manufacturing (HiCAM) Project, state-of-the-art progressive damage and failure analysis (PDFA) tools developed for use with thermosets are being evaluated for use in modeling alternative material systems, like thermoplastics. Experimental low-velocity impact data of a thermoplastic material system, AS4D/PEKK-FC, is presented and includes characterization of the impact damage mechanisms as well as associated load and displacement data. Following the presentation of experimental results, two simulation approaches using the PDFA tool MAT299 in the commercial off-the-shelf finite element software LS-DYNA are employed to predict damage area and force-displacement responses of thermoplastic panels subjected to various impact energies. The first modeling method uses solid elements with a high-density mesh and a ply-by-ply modeling approach similar to previously published work for thermosets. This method has the capability of capturing individual crack development, progression, and delamination on a per ply basis. The second modeling method uses TSHELL elements that have in-plane dimensions that are an order of magnitude larger than the elements used in the solid element approach and reduces the number of elements through the thickness of the laminate. The second method produces a lower-fidelity model with reduced run times that is incapable of monitoring every delamination plane and damage within each ply. Results from both simulation methods are compared to experiment, and limitations of the methods are discussed.

Composite material↗

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 Rao↗

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↗