Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture propagation”

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 523 records · Page 29

Numerical and experimental analysis of mechanically induced failure in electric vehicle battery modules

Mitigating thermal runaway and cell-to-cell propagation is essential for improving the safety of electric and hybrid vehicles. Enhancing digital twin capabilities to predict battery mechanical abuse is particularly critical for automotive and aerospace applications, where crashworthiness is a key concern. Understanding failure conditions and propagation in battery modules during mechanical abuse is complex due to interactions between structural deformation, heat transfer, electrochemical processes, exothermic reactions and mechanical fracture. While prior studies have focused on modeling cell-level behavior, extending these models to module or pack level is necessary for a system level understating of electric vehicle safety. This study develops coupled large deformation finite element models that simultaneously solve for electrochemistry, material failure, internal short circuit and thermal runaway propagation. The models account for mechanical and thermal interactions between lithium-ion cells and other battery components while the contact interfaces are evolving with time. Model-predicted voltage, temperature and force responses are compared with experimental data for validation. The results demonstrate that the approach captures key failure mechanisms, including thermal propagation through heat transfer, electrical propagation from short circuits in parallel-connected cells, and mechanical propagation via penetration and crack formation. These findings show that computational models are valuable tools for understanding battery module failure and providing insight that can reduce the need for extensive experimental testing.

25 ENERGY STORAGE↗

Microfracture in high temperature metal matrix laminates

Computational simulation procedures are described to evaluate the composite microfracture behavior, establish the hierarchy/sequence of fracture modes, and the influence of compliant layers and partial debonding on composite properties and microfracture initiation. These procedures are based upon three-dimensional finite element analysis and composite micromechanics equations. Typical results for the effects of compliant layers and partial debonding, microfracture initiation, and propagation and the thermomechanical cyclic loading on a SiC/Ti15 composite system are presented and discussed. The results show that interfacial debonding follows fiber or matrix fracture, and the thermomechanical cyclic loading severely degrades the composite integrity.

Mital, Subodh K.↗

Effects of stitching on fracture toughness of uniweave textile graphite/epoxy laminates

The effects of through-the-thickness stitching on impact damage resistance, impact damage tolerance, and Mode 1 and Mode 2 fracture toughness of textile graphite/epoxy laminates were studied experimentally. Graphite/epoxy laminates were fabricated from AS4 graphite uniweave textiles and 3501-6 epoxy using Resin Transfer Molding. The cloths were stitched with Kevlar(tm) and glass yarns before resin infusion. Delamination was implanted during processing to simulate impact damage. Sublaminate buckling tests were performed in a novel fixture to measure Compression After Impact (CAI) strength of stitched laminates. The results show that CAI strength can be improved up to 400% by through-the-thickness stitching. Double Cantilever Beam tests were performed to study the effect of stitching on Mode 1 fracture toughness G(sub 1c). It was found that G(sub 1c) increased 30 times for a low stitching density of 16 stitches/sq in. Mode 2 fracture toughness was measured by testing the stitched beams in End Notch Flexure tests. Unlike in the unstitiched beams, crack propagation in the stitched beams was steady. The current formulas for ENF tests were not found suitable for determining G(sub 2C) for stitched beams. Hence two new methods were developed - one based on crack area measured from ultrasonic C-scanning and the other based on equivalent crack area measured from the residual stiffness of the specimen. The G(sub 2c) was found to be at least 5-15 times higher for the stitched laminates. The mechanisms by which stitching increases the CAI strength and fracture toughness are discussed.

Sankar, Bhavani V.↗

Cohesive phase-field chemo-mechanical simulations of inter- and trans- granular fractures in polycrystalline NMC cathodes via image-based 3D reconstruction

The optimal design and durable utilization of lithium-ion batteries necessitates an objective modeling approach to understand fracture and failure mechanisms. This paper presents a comprehensive chemo-mechanical modeling study focused on elucidating fracture-induced damage and degradation phenomena in the polycrystalline Li $\mathcal{x}$ Ni 0.5 Mn 0.3 Co 0.2 O 2 (NMC532) cathode. An innovative approach that utilizes image-based reconstructed 3D geometry as finite element (FE) mesh input is employed to enhance the precision in capturing the convoluted architecture and morphological features. For accurately representing the intricate crack configurations within the polycrystalline system, we adopted the cohesive phase-field fracture (CPF) model. Through the integration of advanced image-based geometry reconstruction technique and the promising CPF modeling approach, lithium (de)intercalation induced crack evolution (e.g., nucleation, propagation, branching and diverse modes including inter-/trans-(intra-) granular patterns) and the resulting chemical degradation can be precisely captured, which is also compared and validated with numerical predictions using a continuum damage model. In particular, this model predicts fracture induced degradation under varying fracture properties of grain boundaries and charging rates; the conclusion that NMC particles comprised of larger grains are predicted to have less degradation than those with smaller grains can also be drawn. This comprehensive analysis provides valuable insights into the fracture and degradation within polycrystalline NMC cathodes.

25 ENERGY STORAGE↗

Prior-to-failure extension of flaws under monotonic and pulsating loadings: Inelastic fatigue

An equation governing the prior to failure crack propagation is proposed. For a rate-sensitive solid containing two-dimensional crack and subject to the tensile mode of fracture, the differential equations are integrated numerically for the loads increasing monotonically in time. The resulting integral curves sigma = sigma(l) and l = l(t), i.e. load vs. crack length and length vs. time, indicate that the growth of cracks in the subcritical range is strongly rate dependent. The fatigue growth, viewed as a sequence of slow growth periods, is simulated on an EAI 380 analogue computer. The fourth power law proposed by Paris is confirmed only within a certain range of high-cycle fatigue propagation and for a rate-insensitive solid. For a more pronounced rate dependency, induced by the viscosity of a solid and/or in the proximity of the final instability point, the growth is markedly enhanced.

Wnuk, M. P.↗

Rayleigh wave studies in lunar rocks.

An ultrasonic surface wave technique described by the author (1971) is used to verify a hypothesis that links the seismic wave propagation velocities in lunar crust, much too low as compared to those on earth, to the extensive fracturing of lunar rock in the absence of liquids and gases which changed drastically the elastic and inelastic properties of lunar rock. Measurements on lunar rock samples and synthetic analogs suggest that the presence of microfractures have influence on both the wave velocity and Q factor in lunar rocks.

Tittmann, B. R.↗

Composite strengthening

The mechanical behavior of unidirectionally reinforced metals is examined, with particular attention to fabrication techniques for artificial composites and eutectic alloys and to principles of fiber reinforcement. The properties of artificial composites are discussed in terms of strength of fiber composites, strength of ribbon-reinforced composites, crack initiation, crack propagation, and creep behavior. The properties of eutectic composites are examined relative to tensile strength, compressive strength, fracture, high-temperature strength, and fatigue. In the case of artificial composites, parallelism of fibers, good bonding between fibers and matrix, and freedom of fibers from damage are all necessary to ensure superior performance. For many eutectic systems there are stringent boundary conditions relative to melt purity and superheat, atmosphere control, temperature gradient, and growth rate in order to provide near-perfect alignment of the reinforcements with a minimum of growth defects.

Stoloff, N. S.↗

GENOA-PFA: Progressive Fracture in Composites Simulated Computationally

GENOA-PFA is a commercial version of the Composite Durability Structural Analysis (CODSTRAN) computer program that simulates the progression of damage ultimately leading to fracture in polymer-matrix-composite (PMC) material structures under various loading and environmental conditions. GENOA-PFA offers several capabilities not available in other programs developed for this purpose, making it preferable for use in analyzing the durability and damage tolerance of complex PMC structures in which the fiber reinforcements occur in two- and three-dimensional weaves and braids. GENOA-PFA implements a progressive-fracture methodology based on the idea that a structure fails when flaws that may initially be small (even microscopic) grow and/or coalesce to a critical dimension where the structure no longer has an adequate safety margin to avoid catastrophic global fracture. Damage is considered to progress through five stages: (1) initiation, (2) growth, (3) accumulation (coalescence of propagating flaws), (4) stable propagation (up to the critical dimension), and (5) unstable or very rapid propagation (beyond the critical dimension) to catastrophic failure. The computational simulation of progressive failure involves formal procedures for identifying the five different stages of damage and for relating the amount of damage at each stage to the overall behavior of the deteriorating structure. In GENOA-PFA, mathematical modeling of the composite physical behavior involves an integration of simulations at multiple, hierarchical scales ranging from the macroscopic (lamina, laminate, and structure) to the microscopic (fiber, matrix, and fiber/matrix interface), as shown in the figure. The code includes algorithms to simulate the progression of damage from various source defects, including (1) through-the-thickness cracks and (2) voids with edge, pocket, internal, or mixed-mode delaminations.

Murthy, Pappu L. N.↗

VCCT with Progressive Nodal Release for Simulating Mixed-Mode Delamination: Formulation, Algorithmic Improvements and Implications

To simulate the propagation of cracks using the Virtual Crack Closure Technique (VCCT), intermediate crack positions between existing node pairs can be modeled by progressively releasing the nodes using fracture mechanics principles. This approach has been implemented in commercial software and can be applied to simulate delamination growth. Recent research has provided insight into the synchronized mixed mode nodal release, extending and formalizing the methodology for delamination growth. This paper documents the development of the Jacobians associated with delamination propagation under quasi-static mixed-mode conditions, leading to improvements in convergence of the implicit nonlinear solvers used.

Mabson, G. E.↗

VCCT with Progressive Nodal Release for Simulating Mixed-Mode Delamination: Formulation, Algorithmic Improvements and Implications

To simulate the propagation of cracks using the Virtual Crack Closure Technique (VCCT), intermediate crack positions between existing node pairs can be modeled by progressively releasing the nodes using fracture mechanics principles. This approach has been implemented in commercial software and can be applied to simulate delamination growth. Recent research has provided insight into the synchronized mixed mode nodal release, extending and formalizing the methodology for delamination growth. This paper documents the development of the Jacobians associated with delamination propagation under quasi-static mixed-mode conditions, leading to improvements in convergence of the implicit nonlinear solvers used.

G E Mabson↗

Shock recompaction of spall damage

Spall fracture is a high strain-rate damage phenomenon associated with shock or impulsive loading events. When a material that has been subjected to shock compression is allowed to release, rarefaction waves propagate into the sample and reduce the internal stress to zero. If multiple rarefaction waves intersect, they generate tension which, if sufficient, can nucleate voids in the material. It has been observed in several works investigating spall fracture that although the shock-wave profiles suggested spall occurred, imaging of the recovered sample revealed no voids or cracks. In this study, we aim to determine whether a second shock event could recompact existing spall damage, and if so, what form does the microstructure at the recompaction interface have? Through a series of gas-gun flyer-plate impact experiments, we demonstrate that modest shock stresses of 2 GPa–3 GPa are enough to both fully compact a damaged copper target back to a state of zero porosity and, furthermore, drive recrystallization of the interface such that there is a new bond formed where the free surfaces were brought together.

42 ENGINEERING↗

Prior-to-failure extension of flaws under monotonic and pulsating loadings.

An equation governing the prior-to-failure crack propagation is proposed. For a rate-sensitive solid containing a two-dimensional crack and subject to the tensile mode of fracture the differential equations are integrated numerically for loads increasing monotonically in time. The resulting integral curves of load vs crack length and length vs time indicate that the growth of cracks in the subcritical range is strongly rate-dependent. The fatigue growth, viewed as a sequence of slow growth periods, is simulated on an EAI 380 analog computer. The fourth-power law proposed by Paris is confirmed only within certain range of high-cycle fatigue propagation and for a rate-insensitive solid. Otherwise, i.e., for a more pronounced rate dependence induced by viscosity of a solid and/or in the proximity of the final instability point, the growth is markedly enhanced.

Wnuk, M. P.↗

The delineation and interpretation of the earth's gravity field

In an attempt to understand the mechanical interaction of a growing lithosphere containing fracture zones with small and large scale mantle convection, which gives rise to geoid anomalies in oceanic regions, a series of fluid dynamical experiments is in progress to investigate: (1) the influence of lithosphere structure, fluid depth and viscosity field on the onset, scale, and evolution of sublithospheric convection; (2) the role of this convection in determining the rate of growth of lithosphere, especially in light of the flattening of the lithosphere bathymetry and heat flow at late times; and (3) combining the results of both numerical and laboratory experiments to decide the dominate factors in producing geoid anomalies in oceanic regions through the thermo-mechanical interaction of the lithosphere and subjacent mantle. The clear existence of small scale convection associated with a downward propagating solidification front (i.e., the lithosphere) and a larger scale flow associated with a discontinuous upward heat flux (i.e., a fracture zone) has been shown. The flows exist simultaneously and each may have a significant role in deciding the thermal evolution of the lithosphere and in understanding the relation of shallow mantle convection to deep mantle convection. This overall process is reflected in the geoid, gravity, and topographic anomalies in the north-central Pacific. These highly correlated fields of intermediate wavelength (approx. 200 to 2000 km) show isostatic compensation by a thin lithosphere for shorter (less than or equal to approx. 500 km), but not the longer, wavelengths. The ultimate, dynamic origin of this class of anomalies is being investigated.

Marsh, Bruce D.↗

A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults

Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off-fault damage over a wide range of spatial and temporal scales. Here, in this paper, we develop a phase-field model for quasi-dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off-fault damage processes with a single formulation. Built on a recently developed phase-field framework for shear fractures with frictional contact, the proposed formulation incorporates rate- and state-dependent friction, radiation damping, and their impacts on fault mechanics and off-fault damage. We show that the numerical results of the phase-field model are consistent with those obtained from well-verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase-field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

58 GEOSCIENCES↗

Molecular dynamics simulation of propagating cracks

Steady state crack propagation is investigated numerically using a model consisting of 236 free atoms in two (010) planes of bcc alpha iron. The continuum region is modeled using the finite element method with 175 nodes and 288 elements. The model shows clear (010) plane fracture to the edge of the discrete region at moderate loads. Analysis of the results obtained indicates that models of this type can provide realistic simulation of steady state crack propagation.

Mullins, M.↗

Predicting The Propagation Of Cracks

Advanced Crack Propagation Predictive Analysis Program, FLAGR04, developed as aid in predicting growth of preexisting flaws and cracks in structural components. Provides fracture-mechanics analyst with computerized method of evaluation of "safe crack-growth life" capabilities of structural components. Also used to evaluate tolerance-to-damage aspects of given structural design. Predicts growth of crack by use of two-dimensional model independently predicting growth of flaw in two directions based on calculation of stress-intensity factors. Input to FLAGR04 consists of initial definition of crack, type of rate solution, type and geometry of flaw, properties of material, data on load spectrum, load-stress functions, and design-limit stress levels. Output includes echo of input with any error or warning message and history of propagation of crack. Written in FORTRAN IV.

Tianlai, HU↗

Crack Growth Prediction Methodology for Multi-Site Damage: Layered Analysis and Growth During Plasticity

A finite element program has been developed to perform quasi-static, elastic-plastic crack growth simulations. The model provides a general framework for mixed-mode I/II elastic-plastic fracture analysis using small strain assumptions and plane stress, plane strain, and axisymmetric finite elements. Cracks are modeled explicitly in the mesh. As the cracks propagate, automatic remeshing algorithms delete the mesh local to the crack tip, extend the crack, and build a new mesh around the new tip. State variable mapping algorithms transfer stresses and displacements from the old mesh to the new mesh. The von Mises material model is implemented in the context of a non-linear Newton solution scheme. The fracture criterion is the critical crack tip opening displacement, and crack direction is predicted by the maximum tensile stress criterion at the crack tip. The implementation can accommodate multiple curving and interacting cracks. An additional fracture algorithm based on nodal release can be used to simulate fracture along a horizontal plane of symmetry. A core of plane strain elements can be used with the nodal release algorithm to simulate the triaxial state of stress near the crack tip. Verification and validation studies compare analysis results with experimental data and published three-dimensional analysis results. Fracture predictions using nodal release for compact tension, middle-crack tension, and multi-site damage test specimens produced accurate results for residual strength and link-up loads. Curving crack predictions using remeshing/mapping were compared with experimental data for an Arcan mixed-mode specimen. Loading angles from 0 degrees to 90 degrees were analyzed. The maximum tensile stress criterion was able to predict the crack direction and path for all loading angles in which the material failed in tension. Residual strength was also accurately predicted for these cases.

James, Mark Anthony↗

Multiscale Modeling of Intergranular Fracture in Aluminum: Constitutive Relation For Interface Debonding

Intergranular fracture is a dominant mode of failure in ultrafine grained materials. In the present study, the atomistic mechanisms of grain-boundary debonding during intergranular fracture in aluminum are modeled using a coupled molecular dynamics finite element simulation. Using a statistical mechanics approach, a cohesive-zone law in the form of a traction-displacement constitutive relationship, characterizing the load transfer across the plane of a growing edge crack, is extracted from atomistic simulations and then recast in a form suitable for inclusion within a continuum finite element model. The cohesive-zone law derived by the presented technique is free of finite size effects and is statistically representative for describing the interfacial debonding of a grain boundary (GB) interface examined at atomic length scales. By incorporating the cohesive-zone law in cohesive-zone finite elements, the debonding of a GB interface can be simulated in a coupled continuum-atomistic model, in which a crack starts in the continuum environment, smoothly penetrates the continuum-atomistic interface, and continues its propagation in the atomistic environment. This study is a step towards relating atomistically derived decohesion laws to macroscopic predictions of fracture and constructing multiscale models for nanocrystalline and ultrafine grained materials.

Yamakov, V.↗