Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite strains”

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
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A comparative study of calibration techniques for finite strain elastoplasticity: Numerically-exact sensitivities for FEMU and VFM

Accurate identification of material parameters is crucial for predictive modeling in computational mechanics. Here, the two primary approaches in the experimental mechanics community for calibration from full-field digital image correlation data are known as finite element model updating (FEMU) and the virtual fields method (VFM). In VFM, the objective function is a squared mismatch between internal and external virtual work or power. In FEMU, the objective function quantifies the weighted mismatch between model predictions and corresponding experimentally measured quantities of interest. It is minimized by iteratively updating the parameters of an FE model. While FEMU is seen as more flexible, VFM is commonly used instead of FEMU due to its considerably greater computational expense. However, comparisons between the two methods usually involve approximations of gradients or sensitivities with finite difference schemes, thereby making direct assessments difficult. Hence, in this study, we compare VFM and FEMU in the context of numerically-exact sensitivities obtained through local sensitivity analyses and the application of automatic differentiation software. To this end, we conduct a series of test cases to assess both methods under practical challenges using a finite strain elastoplasticity model.

Automatic differentiation↗
DOE OSTI2024

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2020

Primal interface debonding formulation for finite strain isotropic plasticity

In this work, a framework is developed for modeling ductile damage of nonlinear materials whose plastic deformation is characterized using rate independent classical plasticity. This method relies on the assumption that the free energy can be decomposed into elastic, plastic and damage parts. A thermodynamically consistent method is derived which satisfies the second law of thermodynamics in the Clausius–Duhem inequality form. The dissipation associated with plasticity takes place in the domain only, while damage dissipation is localized to the interface. The method is developed using Variational Multiscale ideas to obtain definitions of the interface fluxes within a primal formulation analogous to the Discontinuous Galerkin method, which ensures weakly vanishing interface gap prior to reaching a damage initiation criterion. The local nonlinear problem to calculate both plastic deformation gradient and damage variable follows an incremental approach similar to classical plasticity return mapping algorithm. This elastoplastic damage formulation is developed for material undergoing finite strain, and it naturally accommodates a trapezoidal traction separation law (TSL) whose shape can be varied to model either ductile interface behavior or brittle interface behavior. The formulation's performance is assessed through modeling a patch test and a compact tension specimen.

42 ENGINEERING↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2025

Calibrating a finite-strain phase-field model of fracture for bonded granular materials with uncertainty quantification

To study the mechanical behavior of mock high explosives, an experimental and simulation program was developed to calibrate, with quantified uncertainty, a material model of the bonded granular material Idoxuridine and nitroplasticized Estane-5703. This paper reports on the efficacy of such a framework as a generalizable methodology for calibrating material models against experimental data with uncertainty quantification. Additionally, this paper studies the effect of two manufacturing temperatures and three initial granular configurations on the unconfined compressive behavior of the resulting bonded granular materials. In each of these cases, the same calibration framework was used; in that, hundreds of high-fidelity direct numerical simulations using a new, graphics processing unit-enabled, high-performance finite element method software, Ratel, were run to calibrate a finite-strain phase-field fracture model against experimental data. It was found that manufacturing temperature influenced the elastic response of the mock high explosives, with higher temperatures yielding a stiffer response. By contrast, it was found that the initial configuration of the grains had a negligible impact on the overall behavior of the mock high explosives though it remains possible that local damage accumulation within the specimens could be altered by the initial configurations. Overall, the calibration framework was successful at creating well-calibrated models, showing its usefulness as an engineering and scientific tool.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Critical Scaling of Solid Fragmentation at Quasistatic and Finite Strain Rates

Here, using two-dimensional simulations of sheared, brittle solids, we characterize the resulting fragmentation and explore its underlying critical nature. Under quasistatic loading, a power-law distribution of fragment masses emerges after fracture which grows with increasing strain. With increasing strain rate, the maximum size of a grain decreases and a shallower distribution is produced. We propose a scaling theory for distributions based on a fractal scaling of the largest mass with system size in the quasistatic limit or with a correlation length that diverges as a power of rate in the finite-rate limit. Critical exponents are measured using finite-size scaling techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A unified finite strain gradient-enhanced micropolar continuum approach for modeling quasi-brittle failure of cohesive-frictional materials

Here, in this work, a novel framework for modeling quasi-brittle crack propagation and shear band dominated failure of cohesive-frictional materials like concrete, mortar, rock, tough ceramics, energetic materials, but also granular materials like sands or powders in terms of a unified continuum approach is proposed. It is based on a combination of the gradient-enhanced continuum with gradients of internal variables for representing quasi-brittle cracking, and the micropolar continuum, accounting for the deformation of the microstructure. For developing the gradient-enhanced micropolar framework, the set of balance equations and the kinematic relations are derived, and the constitutive relations are established in a general manner. The framework is formulated in a geometrically exact setting, based on the thermodynamically sound theory of hyperelasto-plasticity, and the numerical implementation by means of the finite element method is discussed. For assessing the approach, realizations of this new approach in terms of constitutive models for particular materials are developed. They are applied to numerical benchmark examples, investigating various loading conditions, and the obtained results are validated by means of a comparison with experiments from the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Constraint and size effects in confined layer plasticity

Plane strain finite element analyses are used to model the experiments for a thin metal layer confined between elastic solids. Here, the thin metal layer undergoes elastic-plastic deformations, here modelled by strain gradient plasticity, while the rest of the solid deforms only elastically. Plane strain is assumed and finite strains are accounted for. The microscopic boundary condition at the band interface is initially taken to be zero plastic strain increment, but plastic straining at the interface is allowed when the plastic strain gradient at the interface has exceeded a critical value. When the narrow metal band is inclined 45° relative to the direction of compression of the solid simple shear develops in the band and predictions are rather similar to those of a one-dimensional analysis for an infinite band. For a band at 90° the elastic parts give significant constraint on plastic flow, the stress fields are non-uniform along the band and the average nominal compressive stress shows strong sensitivity to the layer thickness and to the width of the specimen. Remarkably good qualitative, and in some cases quantitative, agreement is found with the experimental observations. Considering also other band angles, both the average nominal shear and compressive stress show variations with the angle.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sound velocities and thermal equation of state of fcc -iron-nickel alloys at high pressure and high temperature: Implications for the cores of Moon and several planets

Fcc-Fe-Ni alloy is believed to be the most dominant solid constitute of moderate-sized terrestrial planetary cores. Investigating the physical properties, especially the density and sound velocity of Fe-Ni alloys and comparing them with seismic observations is an indispensable approach to constructing compositional models for planetary interiors. In this study, we conducted sound velocity measurements on Fe-Ni alloys with 10 wt.% and 20 wt.% Ni up to ∼13.5 GPa and 1073 K, using the ultrasonic interferometry technique in a multi-anvil apparatus in conjunction with synchrotron radiation. By fitting the experimental data to finite strain equations, the bulk and shear moduli and their pressure and temperature derivatives are derived, yielding K S0 =145.8(14) GPa, G 0 = 73.2(7) GPa, K S0 ’ = 5.89(24), G 0 ’ = 2.89(8), (∂K S /∂T) P = -0.0181(12) GPa/K and (∂G/∂T) P = -0.0393(10) GPa/K for fcc-Fe 80 Ni 20 . An examination of the density-velocity relationship shows that compressional wave velocity is insensitive to temperature within the current pressure and temperature range, while shear wave velocity exhibits a large reduction with increasing temperature. Here, extrapolation of the sound velocities following the finite strain theories suggests that much slower Vs should be expected at pressure and temperature conditions corresponding to those of the lunar core. Possible core density and velocity profiles for other moderate planets and satellites, such as Mars, Mercury, and Ganymede are also calculated.

Equation of state↗
DOE OSTI2021

Sound velocities of iron-nickel (Fe90Ni10) alloy up to 8 GPa and 773 K: The effect of nickel on the elastic properties of bcc-iron at high P-T

Sound velocities of iron and iron-based alloys at high pressure and high temperature are crucial for understanding the composition and structure of Earth’s and other telluric planetary cores. In this study, we performed ultrasonic interferometric measurements of both compressional (ν P ) and shear (ν S ) velocities on a polycrystalline body-centered-cubic (bcc)-Fe 90 Ni 10 up to 8 GPa and 773 K. The elastic moduli and their pressure and temperature derivatives are derived from least-square fits to third-order finite strain equations, yielding K S0 = 154.2(8) GPa, G 0 = 73.2(2) GPa, K' S0 = 4.6(2), G' 0 = 1.5(1), ∂K S /∂T = –0.028(1) GPa/K, and ∂G/∂T = –0.023(1) GPa/K. Furthermore, a comparison with literature data on bcc-Fe suggests that nickel not only decreases both P and S wave velocities but also weakens the temperature effects on the elastic moduli of Fe-Ni alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reassessment of Birch's Law on hcp‐Fe From Ultrasonic Sound Velocity Measurement and Implications on the Velocity Profiles of Earth's Inner Core

Here, we performed in situ X-ray diffraction and ultrasonic sound velocity measurements on hcp-Fe up to 15 GPa, 873 K in a multi-anvil apparatus. The elastic moduli and their pressure and temperature derivatives were determined by fitting the velocity and density data to the third-order finite strain equations, yielding K S0 = 169.0(57) GPa, K$^{′}_{S0}$ = 5.4(6), (∂K S /∂T) P = −0.031(3) GPa/K, G 0 = 104.5(27) GPa, G$^{′}_{0}$ = 1.7(2), (∂G/∂T) P = −0.060(2) GPa/K. Within the experimental P-T range, we find significant temperature effect on the density-velocity (ρ − V P or S ) relations and caution the use of the temperature-independent Birch's law for extrapolation to core conditions. Furthermore, temperature-induced velocity decrease is more significant in V S than in V P , offering a possible explanation for the high Poisson's ratio in the core. Extrapolations based on our results together with previous experimental data suggest that V P of hcp-Fe aligns with PREM at Earth's core conditions, while Vs and density are approximately 10% and 2.7% higher than PREM, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exact solutions of burst pressure for thick-walled cylinders in power-law strain hardening steels

Exact solutions of burst pressure for defect-free, thick-walled cylindrical pressure vessels with capped ends are developed using the flow theory of plasticity in terms of the Tresca, von Mises, and Zhu-Leis yield criteria. The material is assumed to obey a power-law strain hardening rule, and the finite strain theory is adopted to describe large plastic deformation. On this basis, internal pressure is obtained as a complex polynomial function of effective strains on the inside and outside surfaces of the pressure vessel with the use of the Bernoulli numbers. At burst failure, the effective strains and stresses on the inside and outside surfaces are first determined, and then three flow solutions of burst pressure are obtained as a power series function of the diameter ratio (D o /D i ), strain hardening exponent (n), and ultimate tensile strength (UTS). The power series solution is confirmed to agree well with the exact flow solution of burst pressure in a closed-form for each yield criterion, and the results show that the von Mises flow solution is an upper bound prediction, the Tresca flow solution is a lower bound prediction, and the Zhu-Leis flow solution is an intermediate prediction of burst pressure. In conclusion, two sets of full-scale burst test data are then utilized to evaluate and validate the proposed flow solutions of burst pressure for both thin-walled pipes and thick-walled cylindrical pressure vessels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Topology optimization of compliant mechanisms under transient thermal conditions

This work considers multi-material topology optimization of compliant mechanisms under transient thermal and quasi-static mechanical conditions wherein thermally actuated devices are optimized for different operating conditions. The materials are modeled using finite strain thermo-hyperelasticity and a two way coupling between the energy balance and equilibrium equations is investigated. The design updates are generated from the gradient-based method of moving asymptotes optimizer and the sensitivities are computed using the time dependent adjoint sensitivity analysis. Results show the impact of designing for short versus long actuation times.

42 ENGINEERING↗