Search NASA⌕ Search

SEARCH · Search NASA

Results for “fracture simulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Numerical fracture simulation of compact and bend specimens

An elastic-plastic finite-element analysis with a critical crack-tip-opening displacement criterion was used to simulate fracture of various size compact and bend specimens made of HY-130 steel. From the calculated load-crack-extension and load-displacement curves, J-resistance (J-R) curves were determined by several methods. The simulated 3-R curves were insensitive to specimen size up to maximum load but were sensitive to specimen configuration for crack extensions greater than 10 percent of the initial uncracked ligament length.

Shivakumar, K. N.↗
NASA NTRS2011

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting, during flight, the residual strength of aircraft structures that sustain discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. Two ductile fracture simulations are presented to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data does, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗
NASA NTRS2011

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting in real time the residual strength of flight structures with discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. A residual strength test of a metallic, integrally-stiffened panel is simulated to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data would, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high-fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗
NASA NTRS1992

Finite-element analyses and fracture simulation in thin-sheet aluminum alloy

A two-dimensional, elastic-plastic finite-element analysis was used with a critical crack-tip-opening angle (CTOA) fracture criterion to model stable crack growth in thin-sheet 2024-T3 aluminum alloy under monotonic loading after precracking at different cyclic stress levels. Tests were conducted on three types of specimens: middle-crack, three-hole-crack and blunt-notch tensile specimens. An experiment technique was developed to measure CTOA during crack growth initiation and stable tearing using a high-resolution video camera and recorder. Crack front shapes were also measured during initiation and stable tearing using a fatigue marker-load technique. Three-dimensional elastic-plastic finite-element analyses of these crack shapes for stationary cracks were conducted to study the crack-front opening displacements. Predicted load against crack extension on middle-crack tension specimens agreed well with test results even for large-scale plastic deformations. The analyses were able to predict the effects of specimen size and precracking stress history on stable tearing. Predicted load against load-line displacements agreed well with test results up to maximum load bu the analyses tended to overpredict displacements as crack grew beyond the maximum load under displacement-controlled conditions. During the initiation phase, the measured CTOA values were high but decreased and remained nearly constant after a small amount of stable tearing. The constant value of CTOA agree well with the calculated value from the finite-element analysis. The larger CTOA values measured at the sheet surface during the initiation phase may be associated with the crack tunneling observed in the tests. Three-dimensional analyses for nonstraight crack fronts predicted much higher displacements near the free surface than in the interior.

Newman, J. C., Jr.↗
NASA NTRS1989

Numerical fracture simulation of bend specimens using a CTOD criterion

Stable crack growth in 50-mm-thick side-grooved three-point bend specimens of HY-130 steel is modeled numerically. A two-dimensional elastoplastic FEM formulation incorporating a CTOD fracture criterion is developed analytically, and the theoretical predictions are compared with experimental data in graphs. Good general agreement is obtained, except for the crack extension, which was typically overpredicted. The value of the J integral at crack initiation is found to be 83 kJ/sq m.

Shivakumar, K. N.↗
NASA NTRS1995

Three-dimensional CTOA and constraint effects during stable tearing in a thin-sheet material

A small strain theory, three-dimensional elastic-plastic finite element analysis was used to simulate fracture in thin sheet 2024-T3 aluminum alloy in the T-L orientation. Both straight and tunneled cracks were modeled. The tunneled crack front shapes as a function of applied stress were obtained from the fracture surface of tested specimens. The stable crack growth behavior was measured at the specimen surface as a function of applied stress. The fracture simulation modeled the crack tunneling and extension as a function of applied stress. The results indicated that the global constraint factor, alpha(sub g), initially dropped during stable crack growth. After peak applied stress was achieved, alpha(sub g) began to increase slightly. The effect of crack front shape on alpha(sub g) was small, but the crack front shape did greatly influence the local constraint and through-thickness crack-tip opening angle (CTOA) behavior. The surface values of CTOA for the tunneled crack front model agreed well with experimental measurements, showing the same initial decrease from high values during the initial 3mm of crack growth at the specimen's surface. At the same time, the interior CTOA values increased from low angles. After the initial stable tearing region, the CTOA was constant through the thickness. The three-dimensional analysis appears to confirm the potential of CTOA as a two-dimensional fracture criterion.

Dawicke, D. S.↗
NASA NTRS1997

Probabilistic Evaluation of Bolted Joints in Polymer Matrix Composites

Computational methods are described to probabilistically simulate fracture in bolted composite structures. Progressive fracture is simulated via an innovative approach independent of stress intensity factors and fracture toughness. The effect on structure damage of design variable uncertainties is quantified. The Fast Probability Integrator is used to assess the scatter in the composite structure response before and after damage. Sensitivity of the response to design variables is evaluated. The methods are demonstrated for bolted joint polymer matrix composite panels under end loads. The effects of fabrication process are included in the simulation of damage in the bolted panel. The results show that the most effective way to reduce the end displacement at fracture is to control the load and ply thickness.

Chamis, C. C.↗
NASA NTRSDate not supplied

Simulation of Dynamic Crushing Behavior of Simulated Carbon/Epoxy Stanchions Using MAT58 and MAT213 in LS-DYNA

Carbon fiber-reinforced polymer (CFRP) composites have shown promise as a material for structures designed to absorb energy in crush-style impact loading. In complementing the lightweight and tailorable characteristics of CFRP structures, the benefits of impact energy absorption are particularly interesting in aerospace vehicles that require excellent crashworthiness performance. However, simulating the behavior of CFRP structures in dynamic, crushing-style impact loading is challenging because of the many damage and failure modes that are essential to capture in the model. For the present investigation, CFRP stanchions were tested using a crash sled experimental test rig. The stanchions were designed by the Composite Materials Handbook-17 (CMH-17) Crashworthiness Working Group for the purpose of comparing experimental crash sled tests to corresponding simulations of the tests. The stanchions are C-channel shaped and represent a geometry common in the interior of aerospace vehicle structures. Explicit simulations in LSDYNA were performed using a well-established composite material model (MAT58) and a next-generation material model (MAT213). Fully integrated shell elements were utilized instead of 3D solid elements to limit computation time. In all models, each of the sixteen plies were represented with individual layers of shell elements with tiebreak contact between each adjacent ply pair to simulate interlaminar fracture. Simulating the crushing of the stanchions occurred in two phases. First, the material models were calibrated using flat specimens that were manufactured with the same layup as the stanchions. While the original goal was calibration of material-related properties, meshdependent behavior was observed in simulations with either material model, and an unstructured mesh was selected to remediate undesirable mesh-dependent failure modes. Additionally, for both the MAT58 and MAT213 models, it was found that either the crush force or the failure mode could be modeled accurately, but no set of parameters could be identified to attain both results in the same model. Once satisfactory calibration was achieved, the same material parameters were applied to the stanchion crush simulations. The stanchion simulations showed that MAT213 more accurately predicted the experimentally determined crush force, and both material models predicted key aspects of the experimentally observed failure modes.

Composite Materials↗
NASA NTRS2006

Computer Simulation of Fracture in Aerogels

Aerogels are of interest to the aerospace community primarily for their thermal properties, notably their low thermal conductivities. While the gels are typically fragile, recent advances in the application of conformal polymer layers to these gels has made them potentially useful as lightweight structural materials as well. In this work, we investigate the strength and fracture behavior of silica aerogels using a molecular statics-based computer simulation technique. The gels' structure is simulated via a Diffusion Limited Cluster Aggregation (DLCA) algorithm, which produces fractal structures representing experimentally observed aggregates of so-called secondary particles, themselves composed of amorphous silica primary particles an order of magnitude smaller. We have performed multi-length-scale simulations of fracture in silica aerogels, in which the interaction b e e n two secondary particles is assumed to be described by a Morse pair potential parameterized such that the potential range is much smaller than the secondary particle size. These Morse parameters are obtained by atomistic simulation of models of the experimentally-observed amorphous silica "bridges," with the fracture behavior of these bridges modeled via molecular statics using a Morse/Coulomb potential for silica. We consider the energetics of the fracture, and compare qualitative features of low-and high-density gel fracture.

Good, Brian S.↗
NASA NTRS1984

Computer simulation of fracture using long range pair potentials

The behavior of a crack in iron is studied using computer simulation. A Morse interatomic potential is used which is of a longer ranged nature than the potentials used in previous studies. With this potential, blunting of the crack tip by spontaneous dislocation nucleation occurs under most conditions. The presence of hydrogen appears to inhibit this blunting and thus to encourage brittle fracture. In the absence of hydrogen, a temperature dependent brittle-ductile transition is observed for some loads. This behavior is quite different from that observed when shorter ranged potentials are used, and it appears to be in somewhat better agreement with experimental results. This agreement with experiment is not conclusive, however, indicating the need for further work to determine more accurate potentials.

Mullins, M.↗
NASA NTRS2009

An Overview of the State of the Art in Atomistic and Multiscale Simulation of Fracture

The emerging field of nanomechanics is providing a new focus in the study of the mechanics of materials, particularly in simulating fundamental atomic mechanisms involved in the initiation and evolution of damage. Simulating fundamental material processes using first principles in physics strongly motivates the formulation of computational multiscale methods to link macroscopic failure to the underlying atomic processes from which all material behavior originates. This report gives an overview of the state of the art in applying concurrent and sequential multiscale methods to analyze damage and failure mechanisms across length scales.

Saether, Erik↗
NASA NTRS2003

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗
NASA NTRS1996

Probabilistic Simulation of Progressive Fracture in Bolted-Joint Composite Laminates

This report describes computational methods to probabilistically simulate fracture in bolted composite structures. An innovative approach that is independent of stress intensity factors and fracture toughness was used to simulate progressive fracture. The effect of design variable uncertainties on structural damage was also quantified. A fast probability integrator assessed the scatter in the composite structure response before and after damage. Then the sensitivity of the response to design variables was computed. General-purpose methods, which are applicable to bolted joints in all types of structures and in all fracture processes-from damage initiation to unstable propagation and global structure collapse-were used. These methods were demonstrated for a bolted joint of a polymer matrix composite panel under edge loads. The effects of the fabrication process were included in the simulation of damage in the bolted panel. Results showed that the most effective way to reduce end displacement at fracture is to control both the load and the ply thickness. The cumulative probability for longitudinal stress in all plies was most sensitive to the load; in the 0 deg. plies it was very sensitive to ply thickness. The cumulative probability for transverse stress was most sensitive to the matrix coefficient of thermal expansion. In addition, fiber volume ratio and fiber transverse modulus both contributed significantly to the cumulative probability for the transverse stresses in all the plies.

Minnetyan, L.↗
NASA NTRS2004

Residual Strength Characterization of a Curved Integrally-Stiffened Panel

Over the years, Finite-element fracture simulation methodology has been very well established at NASA Langley to predict the residual strength of damaged aircraft structures. This methodology has been experimentally verified at NASA Langley for structures ranging from laboratory coupons up to full-scale built-up structural components with single and multiple-site damage cracking. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic-plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intension of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels. In this regard, a series of fracture tests were conducted on curved aluminum-alloy integrally-stiffened panels. These curved panels were subjected to uniaxial tension and pressure loading. During the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (psi(sub c) using three-dimensional code (ZIP3D) and the plane-strain core height (h(sub c)) using two-dimensional code (STAGS).

Seshadri, B. R.↗
NASA NTRS1990

Impact fracture experiments simulating interstellar grain-grain collisions

Oxide and silicate grains condensing during the early phases of the formation of the solar system or in the outflow of stars are exposed to high partial pressures of the low-z elements H, C, N and O and their simple gaseous compounds. Though refractory minerals are nominally anhydrous and non-carbonate, if they crystallize in the presence of H2O, N2 and CO or CO2 gases, they dissolve traces of the gaseous components. The question arises: How does the presence of dissolved gases or gas components manifest itself when grain-grain collisions occur. What are the gases emitted when grains are shattered during a collision event. Researchers report on fracture experiments in ultrahigh vacuum (UHV, approximately less than 10 to the -8th power mbar) designed to measure (by means of a quadrupole mass spectrometer, QMS, with microns to ms time resolution) the emission of gases and vapors during and after impact (up to 1.5 sec). Two terrestrial materials were chosen which represent structural and compositional extremes: olivine (San Carlos, AZ), a densely packed Mg-Fe(2+) silicate from the upper mantle, available as 6 to 12 mm single crystals, and obsidian (Oregon), a structurally open, alkaline-SiO2-rich volcanic glass. In the olivine crystals OH- groups have been identified spectroscopically, as well as H2 molecules. Obsidian is a water-rich glass containing OH- besides H2O molecules. Olivine from the mantle often contains CO2, either as CO2-rich fluid in fluid inclusions or structurally dissolved or both. By analogy to synthetic glasses CO2 in the obsidian may be present in form of CO2 molecules in voids of molecular dimensions, or as carbonate anions, CO3(2-). No organic molecules have been detected spectroscopically in either material. Results indicate that refractory oxide/silicates which contain dissolved traces of the H2O and CO/CO2 components but no spectroscopically detectable traces of organics may release complex H-C-O (possibly H-C-N-O) molecules upon fracture, plus metal vapor. This points: (1) at complex reaction mechanisms between dissolved H2O, CO/CO2 (and N2) components within the mineral structure or during fracture, and (2) at the possibility that similar emission processes occur following grain-grain collisions in interstellar dust clouds.

Freund, Friedemann↗
NASA NTRS1986

Computational simulation of progressive fracture in fiber composites

Computational methods for simulating and predicting progressive fracture in fiber composite structures are presented. These methods are integrated into a computer code of modular form. The modules include composite mechanics, finite element analysis, and fracture criteria. The code is used to computationally simulate progressive fracture in composite laminates with and without defects. The simulation tracks the fracture progression in terms of modes initiating fracture, damage growth, and imminent global (catastrophic) laminate fracture.

Chamis, C. C.↗
NASA NTRS1988

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