A Novel Multiscale Physics-Based Progressive Damage and Failure Modeling Tool for Advanced Composite Structures.

A novel, multiscale mechanics model for predicting the evolution of damage and failure in continuous fiber-reinforced laminates was developed. The thermodynamically-based work potential internal state variable (ISV) theory, Schapery theory (ST), is utilized to model matrix microdamage at the lamina level within a finite element method (FEM) setting. Failure due to transverse cracking and fiber breakage is modeled at the microscale within a repeating unit cell (RUC) using the semi-analytical generalized method of cells (GMC). A multiscale procedure is employed to link the microscale GMC calculations to the macroscale at every integration point in the FEM model. Micromechanics calculations are precluded if the macroscale damage is below some nominal value, increasing the overall computational efficiency of the multiscale scheme. Computational results and predicted failure modes are compared to experimental data of two center-notched, carbon fiber/epoxy panels containing different stacking sequences. A novel, single-scale extension of ST, the enhanced Schapery theory (EST), is also presented. Three additional ISVs are introduced to account for failure via matrix transverse cracking (mode I and mode II) and fiber breakage (mode I only). These ISVs incorporate a characteristic finite element length scale, and are directly related to the fracture toughnesses of the material. In doing so, the pathological mesh dependency, resulting from the failure degradation scheme that was used in the previous multiscale model is eliminated; however, the explicit influence of the fiber-matrix architecture is lost. The EST model is evaluated against the same center-notched panel data. Finally, a mesh objective, smeared crack band model is implemented into the high-fidelity generalized method of cells (HFGMC) micromechanics theory. This failure model utilizes local fields to resolve the orientation of the crack band locally within the subcells of the RUC. The capabilities of the model are demonstrated using an RUC containing multiple randomly oriented fibers subjected to transverse tension and compression. The results of the model are compared to experimental data, and it is concluded that the newly developed model is viable for mesh objective, multiscale simulations.Ph.D.Mechanical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/91483/1/ejpineda_1.pd

Evolving Material Porosity on an Additive Manufacturing Simulation with the Generalized Method of Cells

The effect of material porosity on final part distortion and residual stresses in a selective laser sintering manufacturing simulation is presented here. A time-dependent thermomechanical model is used with the open-source FEA software CalculiX. Effective homogenized material properties for Inconel 625 are precomputed using NASAs Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC). The evolving porosity of the material is estimated with each pass of the laser beam during simulation runtime. A comparison with a homogenous model and the evolving model shows that the evolving porous model predicts larger distortions with greater residual stresses

Hierarchical Coupling of Molecular Dynamics and Micromechanics to Predict the Elastic Properties of Three-Phase and Four-Phase Silicon Carbide Composites

The results obtained from previously conducted molecular dynamics analysis of silicon carbide (-SiC (6H, 4H, & 2H-SiC), -SiC (3C SiC)), silicon and boron nitride, were utilized as inputs in the MAC/GMC micromechanics software to model and evaluate the elastic properties of three-phase SiC/BN/SiC and four-phase SiC/BN/Si/SiC composites. This method of analysis eliminates the need for back-calculation of the apparent properties of the base constituents from the measured ceramic matrix composites properties. The multiscale models are validated against the available data in literature

Enhanced Schapery Theory Software Development for Modeling Failure of Fiber-Reinforced Laminates

Progressive damage and failure analysis (PDFA) tools are needed to predict the nonlinear response of advanced fiber-reinforced composite structures. Predictive tools should incorporate the underlying physics of the damage and failure mechanisms observed in the composite, and should utilize as few input parameters as possible. The purpose of the Enhanced Schapery Theory (EST) was to create a PDFA tool that operates in conjunction with a commercially available finite element (FE) code (Abaqus). The tool captures the physics of the damage and failure mechanisms that result in the nonlinear behavior of the material, and the failure methodology employed yields numerical results that are relatively insensitive to changes in the FE mesh. The EST code is written in Fortran and compiled into a static library that is linked to Abaqus. A Fortran Abaqus UMAT material subroutine is used to facilitate the communication between Abaqus and EST. A clear distinction between damage and failure is imposed. Damage mechanisms result in pre-peak nonlinearity in the stress strain curve. Four internal state variables (ISVs) are utilized to control the damage and failure degradation. All damage is said to result from matrix microdamage, and a single ISV marks the micro-damage evolution as it is used to degrade the transverse and shear moduli of the lamina using a set of experimentally obtainable matrix microdamage functions. Three separate failure ISVs are used to incorporate failure due to fiber breakage, mode I matrix cracking, and mode II matrix cracking. Failure initiation is determined using a failure criterion, and the evolution of these ISVs is controlled by a set of traction-separation laws. The traction separation laws are postulated such that the area under the curves is equal to the fracture toughness of the material associated with the corresponding failure mechanism. A characteristic finite element length is used to transform the traction-separation laws into stress-strain laws. The ISV evolution equations are derived in a thermodynamically consistent manner by invoking the stationary principle on the total work of the system with respect to each ISV. A novel feature is the inclusion of both pre-peak damage and appropriately scaled, post-peak strain softening failure. Also, the characteristic elements used in the failure degradation scheme are calculated using the element nodal coordinates, rather than simply the square root of the area of the element

Numerical Implementation of a Multiple-ISV Thermodynamically-Based Work Potential Theory for Modeling Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments

Numerical Predictions of Damage and Failure in Carbon Fiber Reinforced Laminates Using a Thermodynamically-Based Work Potential Theory

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, referred to as enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Consistent characteristic lengths are introduced into the formulation to govern the evolution of the failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs are derived. The theory is implemented into a commercial finite element code. The model is verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared against the experimental results

Buckling Testing and Analysis of Honeycomb Sandwich Panel Arc Segments of a Full-Scale Fairing Barrel Part 4: Six-ply Out-of-Autoclave Facesheets

Four honeycomb sandwich panel types, representing 1/16th arc segments of a 10-m diameter barrel section of the Heavy Lift Launch Vehicle (HLLV), were manufactured and tested under the NASA Composites for Exploration program and the NASA Constellation Ares V program. Two configurations were chosen for the panels: 6-ply facesheets with 1.125 in. honeycomb core and 8-ply facesheets with 1.000 in. honeycomb core. Additionally, two separate carbon fiber/epoxy material systems were chosen for the facesheets: in-autoclave IM7/977-3 and out-of-autoclave T40-800b/5320-1. Smaller 3 ft. by 5 ft. panels were cut from the 1/16th barrel sections. These panels were tested under compressive loading at the NASA Langley Research Center (LaRC). Furthermore, linear eigenvalue and geometrically nonlinear finite element analyses were performed to predict the compressive response of each 3 ft. by 5 ft. panel. This manuscript summarizes the experimental and analytical modeling efforts pertaining to the panels composed of 6-ply, T40-800b/5320-1 facesheets (referred to as Panels D). To improve the robustness of the geometrically nonlinear finite element model, measured surface imperfections were included in the geometry of the model. Both the linear and nonlinear models yield good qualitative and quantitative predictions. Additionally, it was correctly predicted that the panel would fail in buckling prior to failing in strength. Furthermore, three-dimensional (3D) effects on the compressive response of the panel were studied

An Efficient Semi-Analytical Framework for Micromechanical Modeling of Transverse Cracks in Fiber-reinforced Composites

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/83576/1/AIAA-2010-2693-128.pd

Computational Implementation of a Thermodynamically Based Work Potential Model for Progressive Microdamage and Transverse Cracking in Fiber-reinforced Laminates

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/83581/1/AIAA-2010-2815-444.pd

The EST Model for Predicting Progressive Damage and Failure of Open Hole Bending Specimens

Progressive damage and failure in open hole composite laminate coupons subjected to flexural loading is modeled using Enhanced Schapery Theory (EST). Previous studies have demonstrated that EST can accurately predict the strength of open hole coupons under remote tensile and compressive loading states. This homogenized modeling approach uses single composite shell elements to represent the entire laminate in the thickness direction and significantly reduces computational cost. Therefore, when delaminations are not of concern or are active in the post-peak regime, the version of EST presented here is a good engineering tool for predicting deformation response. Standard coupon level tests provides all the input data needed for the model and they are interpreted in conjunction with finite element (FE) based simulations. Open hole bending test results of three different IM7/8552 carbon fiber composite layups agree well with EST predictions. The model is able to accurately capture the curvature change and deformation localization in the specimen at and during the post catastrophic load drop event