SEARCH · Engineering Papers

Results for “fracture modelling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Signature of transition to supershear rupture speed in the coseismic off-fault damage zone

Most earthquake ruptures propagate at speeds below the shear wave velocity within the crust, but in some rare cases, ruptures reach supershear speeds. The physics underlying the transition of natural subshear earthquakes to supershear ones is currently not fully understood. Most observational studies of supershear earthquakes have focused on determining which fault segments sustain fully grown supershear ruptures. Experimentally cross-validated numerical models have identified some of the key ingredients required to trigger a transition to supershear speed. However, the conditions for such a transition in nature are still unclear, including the precise location of this transition. In this work, we provide theoretical and numerical insights to identify the precise location of such a transition in nature. We use fracture mechanics arguments with multiple numerical models to identify the signature of supershear transition in coseismic off-fault damage. Furthermore, we then cross-validate this signature with high-resolution observations of fault zone width and early aftershock distributions. We confirm that the location of the transition from subshear to supershear speed is characterized by a decrease in the width of the coseismic off-fault damage zone. We thus help refine the precise location of such a transition for natural supershear earthquakes.

42 ENGINEERING↗

NASA NTRS1987

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, John P.↗

NASA NTRS1987

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

Gyekenyesi, J. P.↗

NASA NTRS2006

Dependence of Microelastic-plastic Nonlinearity of Martensitic Stainless Steel on Fatigue Damage Accumulation

Self-organized substructural arrangements of dislocations formed in wavy slip metals during cyclic stress-induced fatigue produce substantial changes in the material microelastic-plastic nonlinearity, a quantitative measure of which is the nonlinearity parameter Beta extracted from acoustic harmonic generation measurements. The contributions to Beta from the substructural evolution of dislocations and crack growth for fatigued martensitic 410Cb stainless steel are calculated from the Cantrell model as a function of percent full fatigue life to fracture. A wave interaction factor f(sub WI) is introduced into the model to account experimentally for the relative volume of material fatigue damage included in the volume of material swept out by an interrogating acoustic wave. For cyclic stress-controlled loading at 551 MPa and f(sub WI) = 0.013 the model predicts a monotonic increase in Beta from dislocation substructures of almost 100 percent from the virgin state to roughly 95 percent full life. Negligible contributions from cracks are predicted in this range of fatigue life. However, over the last five percent of fatigue life the model predicts a rapid monotonic increase of Beta by several thousand percent that is dominated by crack growth. The theoretical predictions are in good agreement with experimental measurements of 410Cb stainless steel samples fatigued in uniaxial, stress-controlled cyclic loading at 551 MPa from zero to full tensile load with a measured f(sub WI) of 0.013.

Cantrell, John H.↗

DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗

NASA NTRS1975

The effect of adhesive layer elasticity on the fracture mechanics of a blister test specimen

An analytical model of a blister type specimen for evaluating adhesive bond strength was developed. Plate theory with shear deformation was used to model the deformation of the plate, and elastic deformation of the adhesive layer is taken into account. It is shown that the inclusion of the elastic deformation of the adhesive layer can have a significant influence in the energy balance calculations of fracture mechanics.

Updike, D. P.↗

NASA NTRS1976

A reliability-based cost effective fail-safe design procedure

The authors have developed a methodology for cost-effective fatigue design of structures subject to random fatigue loading. A stochastic model for fatigue crack propagation under random loading has been discussed. Fracture mechanics is then used to estimate the parameters of the model and the residual strength of structures with cracks. The stochastic model and residual strength variations have been used to develop procedures for estimating the probability of failure and its changes with inspection frequency. This information on reliability is then used to construct an objective function in terms of either a total weight function or cost function. A procedure for selecting the design variables, subject to constraints, by optimizing the objective function has been illustrated by examples. In particular, optimum design of stiffened panel has been discussed.

Hanagud, S.↗

NASA NTRS1993

Plains tectonism on Venus: The deformation belts of Lavinia Planitia

High-resolution radar images from the Magellan spacecraft have revealed the first details of the morphology of the Lavinia Planitia region of Venus. A number of geologic units can be distinguished, including volcanic plains units with a range of ages. Transecting these plains over much of the Lavinia region are two types of generally orthogonal features that we interpret to be compressional wrinkle ridges and extensional grooves. The dominant tectonic features of Lavinia are broad elevated belts of intense deformation that transect the plains with complex geometry. They are many tens to a few hundred kilometers wide, as much as 1000 km long, and elevated hundreds of meters above the surrounding plains. Two classes of deformation belts are seen in the Lavinia region. 'Ridge belts' are composed of parallel ridges, each a few hundred meters in elevation, that we interpret to be folds. Typical fold spacings are 5-10 km. 'Fracture belts' are dominated instead by intense faulting, with faults in some instances paired to form narrow grabens. There is also some evidence for modest amounts of horizontal shear distributed across both ridge and fracture belts. Crosscutting relationships among the belts show there to be a range in belt ages. In western Lavinia, in particular, many ridge and fracture belts appear to bear a relationship to the much smaller wrinkle ridges and grooves on the surrounding plains: ridge morphology tends to dominate belts that lie more nearly parallel to local plains wrinkle ridges, and fracture morphology tends to dominate belts that lie more nearly parallel to local plains grooves. We use simple models to explore the formation of ridge and fracture belts. We show that convective motions in the mantle can couple to the crust to cause horizontal stresses of a magnitude sufficient to induce the formation of deformation belts like those observed in Lavinia. We also use the small-scale wavelengths of deformation observed within individual ridge belts to place an approximate lower limit on the venusian thermal gradient in the Lavinia region at the time of deformation.

Squyres, Steven W.↗

NASA NTRS2008

Mixed Mode Fracture of Plasma Sprayed Thermal Barrier Coatings: Effects of Anisotropy and Heterogeneity

The combined mode I-mode II fracture behavior of anisotropic ZrO2-8wt%Y2O3 thermal barrier coatings was determined in asymmetric flexure loading at both ambient and elevated temperatures. A fracture envelope of KI versus KII was determined for the coating material at ambient and elevated temperatures. Propagation angles of fracture as a function of KI/KII were also determined. The mixed-mode fracture behavior of the microsplat coating material was modeled using Finite Element approach to account for anisotropy and micro cracked structures, and predicted in terms of fracture envelope and propagation angle using mixed-mode fracture theories.

Zhu, Dongming↗

NASA NTRS1972

Fracture Behavior in Nylon 6 Fibers

Electron paramagnetic resonance (EPR) techniques are used to determine the number of free radicals produced during deformation leading to fracture of nylon 6 fibers. A reaction rate molecular model is proposed to explain some of the deformation and bond rupture behavior leading to fracture. High-strength polymer fibers are assumed to consist of a sandwich structure of disordered and ordered regions along the fiber axis. In the disordered or critical flaw regions, tie chains connecting the ordered or crystalline block regions are assumed to have a statistical distribution in length. These chains are, therefore, subjected to different stresses. The effective length distribution was determined by EPR. The probability of bond rupture was assumed to be controlled by reaction-rate theory with a stress-aided activation energy and behavior of various loadings determined by numerical techniques. The model is successfully correlated with experimental stress, strain, and bond rupture results for creep, constant rate loadings, cyclic stress, stress relaxation and step strain tests at room temperature.

Lloyd, B. A.↗

NASA NTRS2016

A Fragment-Cloud Model for Breakup of Asteroids with Varied Internal Structures

As an asteroid descends toward Earth, it deposits energy in the atmosphere through aerodynamic drag and ablation. Asteroid impact risk assessments rely on energy deposition estimates to predict blast overpressures and ground damage that may result from an airburst, such as the one that occurred over Chelyabinsk, Russia in 2013. The rates and altitudes at which energy is deposited along the entry trajectory depend upon how the bolide fragments, which in turn depends upon its internal structure and composition. In this work, we have developed an analytic asteroid fragmentation model to assess the atmospheric energy deposition of asteroids with a range of structures and compositions. The modeling approach combines successive fragmentation of larger independent pieces with aggregate debris clouds released with each fragmentation event. The model can vary the number and masses of fragments produced, the amount of mass released as debris clouds, the size-strength scaling used to increase the robustness of smaller fragments, and other parameters. The initial asteroid body can be seeded with a distribution of independent fragment sizes amid a remaining debris mass to represent loose rubble pile conglomerations, can be given an outer regolith later, or can be defined as a coherent or fractured monolith. This approach enables the model to represent a range of breakup behaviors and reproduce detailed energy deposition features such as multiple flares due to successive burst events, high-altitude regolith blow-off, or initial disruption of rubble piles followed by more energetic breakup of the constituent boulders. These capabilities provide a means to investigate sensitivities of ground damage to potential variations in asteroid structure.

Asteroid↗

NASA NTRS1999

Three-Dimensional Gear Crack Propagation Studied

Gears used in current helicopters and turboprops are designed for light weight, high margins of safety, and high reliability. However, unexpected gear failures may occur even with adequate tooth design. To design an extremely safe system, the designer must ask and address the question, "What happens when a failure occurs?" With gear-tooth bending fatigue, tooth or rim fractures may occur. A crack that propagates through a rim will be catastrophic, leading to disengagement of the rotor or propeller, loss of an aircraft, and possible fatalities. This failure mode should be avoided. A crack that propagates through a tooth may or may not be catastrophic, depending on the design and operating conditions. Also, early warning of this failure mode may be possible because of advances in modern diagnostic systems. One concept proposed to address bending fatigue fracture from a safety aspect is a splittooth gear design. The prime objective of this design would be to control crack propagation in a desired direction such that at least half of the tooth would remain operational should a bending failure occur. A study at the NASA Lewis Research Center analytically validated the crack-propagation failsafe characteristics of a split-tooth gear. It used a specially developed three-dimensional crack analysis program that was based on boundary element modeling and principles of linear elastic fracture mechanics. Crack shapes as well as the crack-propagation life were predicted on the basis of the calculated stress intensity factors, mixed-mode crack-propagation trajectory theories, and fatigue crack-growth theories. The preceding figures show the effect of the location of initial cracks on crack propagation. Initial cracks in the fillet of the teeth produced stress intensity factors of greater magnitude (and thus, greater crack growth rates) than those in the root or groove areas of the teeth. Crack growth was simulated in a case study to evaluate crack-propagation paths. Tooth fracture was predicted from the crackgrowth simulation for an initial crack in the tooth fillet region. This was the desired failure mode for an ultrasafe design. Lastly, tooth loads on the uncracked mesh of the split-tooth design were up to five times greater than those on the cracked mesh if equal deflections of the cracked and uncracked teeth were considered. This effect needs to be considered in the design of a split-tooth configuration. This work was done in-house at Lewis in support of the National Rotorcraft Technology Center project, Ultra-Safe Gear Design, with the Boeing Defense and Space Group. The crack-propagation package was developed by the Cornell Fracture Group at Cornell University. The reported results, which are the initial findings of Lewis gear-crackpropagation research for the Rotorcraft Base Program, will be further investigated to develop generalized gear design guidelines.

Lewicki, David G.↗

NASA NTRS1989

R-Curve Instability Calculations Of Crack Growth

Report discusses use of instability method of calculation and R-curve mathematical models to analyze growth of cracks in fracture-mechanics specimens. In case of single material and structure, such analysis sometimes simple enough to be done on pocket calculator. Where microcomputer or larger computer available, comprehensive program includes libraries of driving-force equations for various configurations and R-curve mathematical models for different materials. Author concludes instability method simple and effective and model equations studied all viable in sense at lease one of them should fit almost any applicable set of crack-growth data. Method and models constitute powerful mathematical tools for analysis of fractures.

Orange, Thomas W.↗

DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1977

Some fundamental aspects of the fatigue and fracture response of composite materials

A simple analytical model, based on the concept that the constraining plies surrounding cracks through the thickness of a given ply or plies control the spacing of those cracks by controlling the surrounding stress state, appears to agree well with experimental data. Such a model can be used to evaluate critically various constraint situations including the prediction of the ability of a given constraint layer to arrest flaw development. The analytical model, which uses an equilibrium element approach in the present case, is sensitive to lamina and laminate elastic properties, flaw size, stacking sequence, and an interlaminar stress transfer distance which is estimated from experimental observations. The analysis is easily extended to study flaw growth and plastic region formation.

Reifsnider, K. L.↗

NASA NTRS1989

Analysis of interlaminar fracture in composites under combined loading

Delamination is a predominant failure mode in continuous fiber reinforced laminated composite structures. One type of delamination is the transverse crack tip delamination which originates at the tip of transverse matrix cracks. An analytical model based on the sublaminate approach and fracture mechanics is developed to study the growth of such delaminations. Plane strain conditions are assumed and estimates are provided for the total strain energy release rate as well as the mode 1 and mode 2 contribution. The energy release rate estimates are used in combination with a simple failure law to predict critical delamination growth strains and stresses. These predictions are compared with experimental data on T300/934 Graphite Epoxy (+ or - 25/90 nano seconds) laminates in the range n=.5 to 8. A good agreement is demonstrated for the range of n where the experimental observations indicate transverse crack tip delamination to be the predominant failure mode.

Armanios, Erian A.↗

NASA NTRS2004

Fracture Mechanics Analysis of an Annular Crack in a Three-concentric-cylinder Composite Model

A boundary-value problem governing a three-phase concentric-cylinder model was analytically modeled to analyze annular interfacial crack problems with Love s strain functions in order to find the stress intensity factors (SIFs) and strain energy release rates (SERRs) at the tips of an interface crack in a nonhomogeneous medium. The complex form of a singular integral equation (SIE) of the second kind was formulated using Bessel s functions in the Fourier domain, and the SIF and total SERR were calculated using Jacoby polynomials. For the validity of the SIF equations to be established, the SIE of the three-concentric-cylinder model was reduced to the SIE for a two-concentric-cylinder model, and the results were compared with the previous results of Erdogan. A preliminary set of parametric studies was carried out to show the effect of interphase properties on the SERR. The method presented here provides insight about the effect of interphase properties on the crack driving force.

Kuguoglu, Latife H.↗

DOE OSTI2022

EGS Collab: Modeling and Simulation Working Group Teleconference Series (129-130)

This submission contains the presentation slides and recordings from EGS Collab Modeling and Simulation Working Group (MSWG) teleconferences number 129 through 130. These teleconferences served three objectives for the project: 1) share simulation results, 2) communicate field activities and results to the simulation teams, and 3) hold open scientific discussions on EGS topics.

15 GEOTHERMAL ENERGY↗