Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

In-situ strain behavior and BISON simulations of Zircaloy cladding subjected to temperature cycling separate-effects tests in a steam environment

Understanding fuel system performance during anticipated transients without scram (ATWSs) in boiling water reactors (BWRs) is necessary for refining current and future safety limits. High-fidelity material models and simulations are fundamental to rigorous assessment of zirconium-based cladding performance. However, experimental thermomechanical data during simulated ATWSs to validate these modes are limited. To provide relevant in-situ data, Zircaloy-4 cladding was subjected to cyclic heating in a steam environment to simulate an out-of-pile BWR ATWS. Digital image correlation was used to capture the cladding strain behavior in-situ for comparison against simulations using the BISON finite element code. Conventional high-temperature models were compared using multiple schemes to gain a better understanding of the applicability of three BISON models to BWR ATWS: (1) the default combination of creep models in BISON, (2) the high-temperature Erbacher model alone, and (3) the low-temperature Limback-Andersson model alone. The cases run with the Limback-Andersson model alone produced the lowest root mean square error (RMSE). The lowest RMSE for the Limback-Andersson model alone was 0.659%, and the highest RMSE reported was 4.22%. A data gap within the model in the temperature regime of interest was also identified, and to account for this gap, the current model in BISON is linearly interpolated between two separate datasets. In conclusion, this evaluation highlights the need to either develop a new model or to improve the existing model to capture transient creep effects resulting from a cyclic temperature transient.

Anticipated transient without SCRAM↗
DOE OSTI2023

Grayscale Digital Light Processing Gradient Printing for Stress Concentration Reduction and Material Toughness Enhancement

Avoiding stress concentrations is essential to achieve robust parts since failure tends to originate at such concentrations. With recent advances in multimaterial additive manufacturing, it is possible to alter the stress (or strain) distribution by adjusting the material properties in selected locations. Here, we investigate the use of grayscale digital light processing (g-DLP) 3D printing to create modulus gradients around areas of high stress. These gradients prevent failure by redistributing high stresses (or strains) to the neighboring material. The improved material distributions are calculated using finite element analysis. The much-enhanced properties are demonstrated experimentally for thin plates with circular, triangular, and elliptical holes. In conclusion, this work suggests that multimaterial additive manufacturing techniques like g-DLP printing provide a unique opportunity to create tougher engineering materials and parts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Efficient mapping between void shapes and stress fields using Deep Convolutional Neural Networks with sparse data

Establishing fast and accurate structure-to-property relationships is an important component in the design and discovery of advanced materials. Physics-based simulation models like the finite element method (FEM) are often used to predict deformation, stress, and strain fields as a function of material microstructure in material and structural systems. Such models may be computationally expensive and time intensive if the underlying physics of the system is complex. This limits their application to solve inverse design problems and identify structures that maximize performance. In such scenarios, surrogate models are employed to make the forward mapping computationally efficient to evaluate. However, the high dimensionality of the input microstructure and the output field of interest often renders such surrogate models inefficient, especially when dealing with sparse data. Deep convolutional neural network (CNN) based surrogate models have shown great promise in handling such high-dimensional problems. In this paper, a single ellipsoidal void structure under a uniaxial tensile load represented by a linear elastic, high-dimensional and expensive-to-query, FEM model. We consider two deep CNN architectures, a modified convolutional autoencoder framework with a fully connected bottleneck and a UNet CNN, and compare their accuracy in predicting the von Mises stress field for any given input void shape in the FEM model. Additionally, a sensitivity analysis study is performed using the two approaches, where the variation in the prediction accuracy on unseen test data is studied through numerical experiments by varying the number of training samples from 20 to 100.

surrogate modeling; convolutional neural networks;↗
DOE OSTI2022

The Mechanical Response of High Burnup 17x17 PWR Fuel Rods Under Bending

As a part of the US Department of Energy Office of Nuclear Energy High Burnup Spent Fuel Data Project, Oak Ridge National Laboratory has performed 4pt-bend tests on a series of 17x17 PWR fuel rods that were discharged with rod average burnups between 50 and 60 GWd/MTU. The rods include four different kinds of fuel cladding: standard Zircaloy-4 (Zirc-4), low-tin (LT) Zirc-4, ZIRLO, and M5. Six rods of each type (combining Zirc-4 and LT Zirc-4 as one type) were tested to failure at room temperature and 200 oC. Three of each type experienced a full-length heat treatment cycle from RT to 400 oC to RT with an 8 hr soak at 400 oC to simulate the expected limiting vacuum drying condition, which potentially induces hydride reorientation and causes some annealing of irradiation damage in the cladding. Overall, the fuel rods show significant capacity to deflect in bending and exhibit non-linear deformation before failing. The calculated failure strains range from 1.5% to 3.8%. The simulated vacuum drying heat treatment had the largest effect on fuel rods with ZIRLO and M5 cladding. In both cases the heat treatment caused increases in the plastic strain at failure, suggesting some annealing of irradiation damage is occurring. The finite element modeling of the 4pt bending test assume the pellet is a collection of pellet wedge pieces. The results for the Zircaloy-4 rods suggest that the pellets are bonded to the cladding and friction controls the relative motion of the pellet wedge pieces.

Cantonwine, Paul↗
DOE OSTI2022

3D Bragg Coherent Diffraction Imaging of Extended Nanowires: Defect Formation in Highly Strained InGaAs Quantum Wells

InGaAs quantum wells embedded in GaAs nanowires can serve as compact near-infrared emitters for direct integration onto Si complementary metal oxide semiconductor technology. While the core–shell geometry in principle allows for a greater tuning of composition and emission, especially farther into the infrared, the practical limits of elastic strain accommodation in quantum wells on multifaceted nanowires have not been established. One barrier to progress is the difficulty of directly comparing the emission characteristics and the precise microstructure of a single nanowire. Here, in this study, we report an approach to correlating quantum well morphology, strain, defects, and emission to understand the limits of elastic strain accommodation in nanowire quantum wells specific to their geometry. We realize full 3D Bragg coherent diffraction imaging (BCDI) of intact quantum wells on vertically oriented epitaxial nanowires, which enables direct correlation with single-nanowire photoluminescence. By growing In 0.2 Ga 0.8 As quantum wells of distinct thicknesses on different facets of the same nanowire, we identified the critical thickness at which defects are nucleated. A correlation with a traditional transmission electron microscopy analysis confirms that BCDI can image the extended structure of defects. Finite element simulations of electron and hole states explain the emission characteristics arising from strained and partially relaxed regions. This approach, imaging the 3D strain and microstructure of intact nanowire core–shell structures with application-relevant dimensions, can aid the development of predictive models that enable the design of new compact infrared emitters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Meso-scale modelling of compressive fracture in concrete with irregularly shaped aggregates

This paper presents a meso-scale modelling framework to investigate the fracture process in concrete subjected to uniaxial and biaxial compression accounting for its mesostructural characteristics. 3D mesostructure of concrete consisting of coarse aggregates, mortar and interfacial transition zone between them was developed using an in-house code based on the Voronoi tessellation and splining method, which enables to generate the realistic-look aggregates with controllable structural features such as content, location, size and shape. Based on the generated 3D mesostructure, the concrete damage plasticity approach was employed to simulate the compressive fracture behaviour of concrete in terms of crack morphology and stress-strain response against the shape parameters of aggregate. Results indicate that the shape of aggregate has a negligible effect on compressive strength of concrete, which is highly associated with the random location and size distribution of aggregate. The aggregate irregularity has a significant influence on crack initiation and growth of concrete.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Prediction of high strain polar phases in antiferroelectric PbZrO 3 from multiscale approach

PbZrO 3 is regarded as the first antiferroelectric and currently under intense reexamination. We propose a multiscale approach that combines classical and first principles density functional theory based simulations to explore polar phases in this material. Application of this approach to PbZrO 3 predicts three polar phases, none of which is the common $R3c$. The lowest energy polar phase - $Cc$ - is metastable at zero Kelvin but can be stabilized by application of electric field. Its structural and polarization response to the electric field is in good agreement with experimental data from the literature. The other two phases are $Ima2$ and $I4cm$ of orthorhombic and tetragonal symmetry, respectively, and predicted to stabilize at finite temperatures and under application of larger electric fields. These two phases exhibit very large strain which is of technological significance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

ALGEBRA v.1.27

The ALGEBRA program allows the user to manipulate data from a finite element analysis before it is plotted. The finite element output data is in the form of variable values (e.g., stress, strain, and velocity components) in an EXODUS II database. The ALGEBRA program evaluates user-supplied functions of the data and writes the results to an output EXODUS II database that can be read by plot programs.

Sjaardema, Gregory↗
DOE OSTI2021

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Finite domain solution of a KGD hydraulic fracture in the viscosity-dominated regime

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future.

Domain-based method↗
DOE OSTI2021

Plasticity and fracture behavior of Inconel 625 manufactured by laser powder bed fusion: Comparison between as-built and stress relieved conditions

In this work, the influence of stress relief on the plasticity and fracture behavior of Inconel 625 fabricated through laser powder bed fusion additive manufacturing (AM) was investigated. The as-built versus stress relieved microstructures were compared, showing similar grain structures but the presence of ~10 vol % δ phase in the stress relieved condition, and no δ phase in the as-built condition. Mechanical tests under plane strain tension were performed on the stress relieved samples, and an anisotropic plasticity model was calibrated and validated using finite element simulations. Uniaxial and notched tension tests were performed on both as-built and stress relieved samples to probe the effect of stress relief on stress state- and direction-dependent fracture behavior. It was found that on average, the fracture strain of the stress relieved samples along the build direction was 30% higher than that along the perpendicular build direction in the stress state range studied, and the stress relief heat treatment resulted in a 45% decrease in fracture strain. The fracture strain in stress relieved samples was more strongly dependent on stress state than in as-built samples.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling of surface mechanical behaviors of soft elastic solids: theory and examples

Surfaces of soft solids can have significant surface stress, extensional modulus and bending stiffness. Previous theoretical studies have usually examined cases in which both the surface stress and bending stiffness are constant, assuming small deformation. In this work we consider a general formulation in which the surface can support large deformation and carry both surface stresses and surface bending moments. We show that the large deformation theory can be reduced to the classical linear theory (Shuttleworth equation). We obtain exact solutions for problems of an inflated cylindrical shell and bending of a plate with a finite thickness. Our analysis illustrates the different manners in which surface stiffening and surface bending stabilize these structures. We discuss how the complex surface constitutive behaviors affect the stress field of the bulk. Our calculation provides insights into effects of strain-dependent surface stress and surface bending in the large deformation regime, and can be used as a model to implement surface finite elements to study large deformation of complex structures.

36 MATERIALS SCIENCE↗
DOE OSTI2026

ACT University Collaboration Proposal: Role of manufacturing defects on material failure under dynamic loading for developing enhanced failure models and theories (Final Report)

Research on the tensile behavior of additively manufactured 316L stainless steel coupon specimens at increasing strain rates was conducted over the past year at Penn State. Dynamic loading rates into the 1000 1/s loading rates were performed on specimens with nearly 1.32” gage lengths. Stress-strain plots show ductile and plastic behavior well beyond 12% strain with necked specimens, having smaller effective gage lengths showing up to 56% ultimate strain. Additional tests performed on compact tension specimens helped with simulation work to understand the deformation behavior. Using finite elements, it was possible to determine the feasibility of a comprehensive experimental validation study towards an improved model for AM material failure – specifically, the Bai-Wierzbicki approach which accounts for lode angle and triaxiality. A non-significant number of tests are projected for eight specific failure nodes, with additional replicates to provide strain-rate capabilities to the existing formulation. One such correction factor is explored for quasi-static, notched specimens. Finally, stress intensity factor was explored using a set of compact-tension specimens which would also provide useful validation data for any simulation work.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A continuum framework for coupled solid deformation–fluid flow through anisotropic elastoplastic porous media

We present a continuum framework for coupled solid deformation–fluid flow in anisotropic elastoplastic porous media. A thermodynamic formulation of the coupled processes gives rise to an anisotropic Biot tensor that is a function of drained elastic tangent moduli tensor of the solid skeleton and the intrinsic bulk modulus of the solid constituent. Two effective stress measures emerge from the formulation, namely, σ', which is energy-conjugate to the elastic strain, and σ'', which is energy-conjugate to the plastic strain. For the special case of transverse isotropy that is commonly encountered in natural rocks, the Biot tensor can be expressed in terms of its normal and tangential components to the bedding plane, along with a microstructure tensor. Apart from its thermodynamic consistency, an advantage of this new formulation is that standard mixed finite element formulation can be employed to discretize the domain and solve initial boundary-value problems. We conduct plane strain simulations of coupled solid deformation–fluid flow in a transversely isotropic porous medium to demonstrate the impacts of material anisotropy, stress history, and the Biot tensor on the system response.

42 ENGINEERING↗
DOE OSTI2026

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗
DOE OSTI2025

Machine learning-enabled multiscale modeling of mechanical deformation of aluminum and Al-SiC nanocomposites

A machine learning-enabled multiscale framework is developed for modeling the mechanical response of both pure metal and nanoparticle-reinforced metal matrix nanocomposites (MMNCs). Using aluminum–silicon carbide (Al-SiC) as an example MMNC, atomistic simulations reveal three distinct deformation mechanisms (i.e., defect-free, dislocation-based, and interface separation) governed by the interfaces between the Al matrix and SiC nanoparticles. As compared with single crystal Al, the lattice undergoes a more abrupt failure once the dislocation network becomes extensive and void nucleation initiates, whereas in Al-SiC, nanoparticle interfaces enable a more gradual progression of damage. These mechanisms are captured through a combined classification-regression neural network surrogate model that bridges atomic-scale insights with continuum-scale finite element analysis. Machine learning-enabled multiscale modeling of pure Al accurately predicted strain localization and confirmed by in-situ scanning electron microscopic tensile testing on perforated Al specimens. This study underscores the promise of integrating physics-informed machine learning with hierarchical modeling to capture the interface dominated phenomena and guide the design of advanced MMNCs.

Al-SiC↗
DOE OSTI2023

A materials data framework and dataset for elastomeric foam impact mitigating materials

The availability of materials data for impact-mitigating materials has lagged behind applications-based data. For example, data describing on-field helmeted impacts are available, whereas material behaviors for the constituent impact-mitigating materials used in helmet designs lack open datasets. Here, we describe a new FAIR (findable, accessible, interoperable, reusable) data framework with structural and mechanical response data for one example elastic impact protection foam. The continuum-scale behavior of foams emerges from the interplay of polymer properties, internal gas, and geometric structure. This behavior is rate and temperature sensitive, therefore, describing structure-property characteristics requires data collected across several types of instruments. Data included are from structure imaging via micro-computed tomography, finite deformation mechanical measurements from universal test systems with full-field displacement and strain, and visco-thermo-elastic properties from dynamic mechanical analysis. These data facilitate modeling and design efforts in foam mechanics, e.g., homogenization, direct numerical simulation, or phenomenological fitting. The data framework is implemented using data services and software from the Materials Data Facility of the Center for Hierarchical Materials Design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Impact of Plate Thickness and Joint Geometry on Residual Stresses in 347H Stainless Steel Welds

Weldments of 347H stainless steel are potentially susceptible to stress relaxation cracking at elevated service temperatures. Mitigation of stress relaxation cracking susceptibility within a multipass weld requires a good understanding of welding practices and manufacturing techniques to control high tensile residual stresses. Here in this work, the dependence of residual stress distribution in 347H stainless steel on base plate thickness, joint geometry design, and preheating condition was systematically investigated by using three-dimensional finite element models. The finite element models were validated through good agreement between neutron diffraction measurements and calculated elastic strains. The single-V-groove welds with and without a preheating step all produced similar peak von Mises residual stresses, above 450 MPa, within both the fusion zone and heat-affected zone (HAZ). In plates thicker than 0.5 in. (12.7 mm), high tensile residual stress could be observed in a relatively large area, from the middle of the plate thickness to underneath the top surface. A double-V groove shifted the high tensile stress area to the middle thickness of the weld. A single-J-groove weld was able to confine the residual stress to a very small region near the middle thickness within the fusion zone and suppressed the von Mises residual stress within the HAZ to below 400 MPa.

14 SOLAR ENERGY↗