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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Embedded smoothed particle hydrodynamics
In this paper, a hybrid method for embedding a meshfree smoothed particle hydrodynamics into mesh-based ALE hydro schemes is presented. The method relies on the embedded mesh framework presented in Puso et al. (2015, 2012) and utilizes SPH (Owen, 2010) in the solid foreground; when an appropriate damage or failure model is used with SPH, it becomes a natural way to model material fracture and cracking in high-explosive driven experiments. Overall, we provide various examples which validate the coupling between SPH and the embedded mesh method.
Three-dimensional phase-field modeling of porosity dependent intergranular fracture in UO2
Not Available
Phase-field simulations of fission gas bubbles in high burnup UO2 during steady-state and LOCA transient conditions
To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized this is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory's Marmot application for phase-field simulation of nuclear materials. Simulations of bubble growth in the high burnup region during steady-state conditions showed that initially overpressurized bubbles decreased in pressure during growth, but still remained above equilibrium pressure. During a LOCA transient, simulations of bubbles in the high burnup region showed that bubble size did not change significantly, and the pressure increase due to the transient was calculated and passed to a phase-field model of fracture.
Numerical manifold method modeling of coupled processes in fractured geological media at multiple scales
Not Available
Review of High‐Speed Digital Image Correlation: Advancements and Good Practices
This paper reviews the current state of the art in high‐speed (HS) and ultrahigh‐speed (UHS) digital image correlation (DIC) techniques, emphasizing their critical role in experimental research across various scientific domains. HS and UHS DIC have evolved significantly, driven by advancements in camera systems, image processing algorithms and experimental methodologies. These developments have opened new avenues for capturing and analysing dynamic events with unprecedented temporal and spatial detail, but not without introducing challenges such as optical distortions, motion blur and lighting issues that can affect measurement quality. This review advocates for standardized reporting practices in HS/UHS DIC methodologies to improve reproducibility and reliability across studies, drawing on guidelines from the International Digital Image Correlation Society (iDICs). Through a comprehensive analysis of over 150 articles, this review identifies key advancements in imaging technology and their application in six research domains: material characterization, test development, fracture mechanics, model validation, ballistic and explosive phenomena assessment and measurement uncertainties. Distinctions between two‐dimensional (2D) and stereo‐DIC applications are explored, offering insights into their practical implementation and trade‐offs. Good practices for HS/UHS DIC applications are proposed along with suggestions for future directions for this evolving field, highlighting the indispensable role of technological innovation in expanding the capabilities of optical metrology.
The Convergence of Heat, Groundwater & Fracture Permeability: Innovative Play Fairway Modelling Applied to the Tularosa Basin (Final Report)
This report summarizes the activities and key findings of the project team occurring during Phase 1 (August 2014 – October 2015) of the Tularosa Basin Geothermal Play Fairway Analysis Project.
Local-nonlocal coupling in Emu/PDMS
A technique called the splice method for coupling local to peridynamic subregions of a body is described. The method relies on ghost nodes, whose values of displacement are interpolated from nearby physical nodes, to make each subregion visible to the other. In each time step, the nodes in each subregion treat the nodes in the other subregion as boundary conditions. Adaptively changing the subregions is possible through the creation and deletion of ghost nodes. Example problems in 2D and 3D illustrate how the method is used to perform multiscale modeling of fracture and impact events within a larger structure.
Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code
Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.
Modeling of diffusion in complex fractured porous media by using hierarchical material properties.
Abstract not provided.
Apollo 14 and 16 Active Seismic Experiments, and Apollo 17 Lunar Seismic Profiling
Seismic refraction experiments were conducted on the moon by Apollo astronauts during missions 14, 16, and 17. Seismic velocities of 104, 108, 92, 114 and 100 m/sec were inferred for the lunar regolith at the Apollo 12, 14, 15, 16, and 17 landing sites, respectively. These data indicate that fragmentation and comminution caused by meteoroid impacts has produced a layer of remarkably uniform seismic properties moonwide. Brecciation and high porosity are the probable causes of the very low velocities observed in the lunar regolith. Apollo 17 seismic data revealed that the seismic velocity increases very rapidly with depth to 4.7 km/sec at a depth of 1.4 km. Such a large velocity change is suggestive of compositional and textural changes and is compatible with a model of fractured basaltic flows overlying anorthositic breccias. 'Thermal' moonquakes were also detected at the Apollo 17 site, becoming increasingly frequent after sunrise and reaching a maximum at sunset. The source of these quakes could possibly be landsliding.
Status of the Ford program to evaluate ceramics for stator applications in automotive gas turbine engines
The paper reviews the progress of the major technical tasks of the DOE/NASA/Ford program Evaluation of Ceramics for Stator Applications in Automotive Gas Turbine Engines: reliability prediction, stator fabrication, material characterization, and stator evaluation. A fast fracture reliability model was prepared for a one-piece ceramic stator. Periodic inspection results are presented.
Possible role of crustal flexure in the initial detachment of extensional allochthons
The existence of low-angle normal faults indicates that the ratio of shear stress (TAU) to normal stress (sigma sub N) needed to cause slip on faults is substantially less than would be predicted based on experimental data. Because the tensional strength of rock at a large scale is exceedingly low, the upper plate of a low-angle normal fault cannot be pulled down the fault ramp, but must be driven down it by its own weight. The active or recently active Sevier Desert detachment fault in western Utah dips regionally at 12 deg. The radio of shear stress to normal stress due to the weight of the upper plate on a 12 deg dipping fault surface is 0.2. In contrast, laboratory experiments indicate that slip on fracture surfaces occurs with almost all rock types when (tau/sigma) reaches values of 0.6 to 0.85, corresponding to normal-fault dips of 30 deg to 40 deg. Seismological data indicate that low deviatoric stresses are associated with movement on faults of other geometrics and are not unique to low-angle normal faults. It thus appears that approximately planar fault zones with surface areas of hundreds to thousands of square kilometers have different mechanical properties than would be predicted based on laboratory studies of fractured rock. Modeling of stresses indicates that deviatoric stresses greater than 5kb exist and are sustainable in oceanic lithosphere, and that failure occurs when tau/sigma sub N approaches 0.6 to 0.85.
Sublaminate analysis of interlaminar fracture in composites. I - Analytical model
A simple analysis method based on a transverse shear-deformation theory and a sublaminate approach is utilized to analyze a mixed-mode edge delamination specimen. The analysis provides closed-form expressions for the interlaminar shear stresses ahead of the crack, the total strain-energy release rate, and the strain-energy release-rate components. The parameters controlling the behavior are identified. The effect of specimen stacking sequence and delamination interface on the strain energy release rate components is investigated. Results are compared with a finite-element simulation for reference. The simple nature of the method makes it suitable for preliminary design analyses which require a large number of configurations to be evaluated quickly and economically.
Matrix cracking in composite laminates with resin-rich interlaminar layers
Fracture mechanics analysis and test data for a toughened composite material that has a resin-rich interlaminar layer (RIL) were used to investigate in situ strength. Exposure to a range of environmental conditions was considered. A parametric analysis study was performed to judge the effects of laminate and material variables. A finite thickness effect, indicating an interaction between ply group thickness and effective flaw size, was found dominant. The magnitude of the effect was directly related to RIL stiffness. In situ strength was found to decrease with decreasing RIL stiffness. This work indicates the need to use a fracture mechanics model of actual lamina microstructure and heterogeneous properties to predict in situ strength in materials with RIL.
Damage and strength of composite materials: Trends, predictions, and challenges
Research on damage mechanisms and ultimate strength of composite materials relevant to scaling issues will be addressed in this viewgraph presentation. The use of fracture mechanics and Weibull statistics to predict scaling effects for the onset of isolated damage mechanisms will be highlighted. The ability of simple fracture mechanics models to predict trends that are useful in parametric or preliminary designs studies will be reviewed. The limitations of these simple models for complex loading conditions will also be noted. The difficulty in developing generic criteria for the growth of these mechanisms needed in progressive damage models to predict strength will be addressed. A specific example for a problem where failure is a direct consequence of progressive delamination will be explored. A damage threshold/fail-safety concept for addressing composite damage tolerance will be discussed.
Crack Turning Mechanics of Composite Wing Skin Panels
The safety of future composite wing skin integral stiffener panels requires a full understanding of failure mechanisms of these damage tolerance critical structures under both in-plane and bending loads. Of primary interest is to derive mathematical models using fracture mechanics in anisotropic cracked plate structures, to assess the crack turning mechanisms, and thereby to enhance the residual strength in the integral stiffener composite structures. The use of fracture mechanics to assess the failure behavior in a cracked structure requires the identification of critical fracture parameters which govern the severity of stress and deformation field ahead of the flaw, and which can be evaluated using information obtained from the flaw tip. In the three-year grant, the crack-tip fields under plane deformation, crack-tip fields for anisotropic plates and anisotropic shells have been obtained. In addition, methods for determining the stress intensity factors, energy release rate, and the T-stresses have been proposed and verified. The research accomplishments can be summarized as follows: (1) Under plane deformation in anisotropic solids, the asymptotic crack-tip fields have been obtained using Stroh formalism; (2) The T-stress and the coefficient of the second term for sigma(sub y), g(sub 32), have been obtained using path-independent integral, the J-integral and Betti's reciprocal theorem together with auxiliary fields; (3) With experimental data performed by NASA, analyses indicated that the mode-I critical stress intensity factor K(sub Q) provides a satisfactory characterization of fracture initiation for a given laminate thickness, provided the failure is fiber-dominated and crack extends in a self-similar manner; (4) The high constraint specimens, especially for CT specimens, due to large T-stress and large magnitude of negative g(sub 32) term may be expected to inhibit the crack extension in the same plane and promote crack turning; (5) Crack turning out of crack plane in generally anisotropic solids under plane deformation has been studied; (6) The role of T-stress and the higher-order term of sigma(sub y) on the crack turning and stability of the kinked crack has been quantified; (7) Asymptotic crack-tip fields including the effect of transverse shear deformation (Reissner plate theory) in an anisotropic plate under bending, twisting moments, and transverse shear loads has been presented; (8) The expression of the path-independent J-integral in terms of the generalized stress and strain has been derived; (9) Asymptotic crack-tip fields including the effect of transverse shear deformation (Reissner shallow shell theory) in a general anisotropic shell has been developed; (10) The Stroh formalism was used to characterize the crack tip fields in shells up to the second term and the energy release rate was expressed in a very compact form.
Failure Criteria for FRP Laminates in Plane Stress
A new set of six failure criteria for fiber reinforced polymer laminates is described. Derived from Dvorak's fracture mechanics analyses of cracked plies and from Puck's action plane concept, the physically-based criteria, denoted LaRC03, predict matrix and fiber failure accurately without requiring curve-fitting parameters. For matrix failure under transverse compression, the fracture plane is calculated by maximizing the Mohr-Coulomb effective stresses. A criterion for fiber kinking is obtained by calculating the fiber misalignment under load, and applying the matrix failure criterion in the coordinate frame of the misalignment. Fracture mechanics models of matrix cracks are used to develop a criterion for matrix in tension and to calculate the associated in-situ strengths. The LaRC03 criteria are applied to a few examples to predict failure load envelopes and to predict the failure mode for each region of the envelope. The analysis results are compared to the predictions using other available failure criteria and with experimental results. Predictions obtained with LaRC03 correlate well with the experimental results.