Search NASA⌕ Search

SEARCH · Search NASA

Results for “Elastic-Plasticity”

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.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Elastic-plastic conductor damage evaluation at over 0.4% strain using a high-stress REBCO coil

Abstract Recent reports on screening current stress simulations of high-field REBCO magnets frequently present peak stresses over 1 GPa. However, this result is probably an unrealistic artifact of purely elastic calculations, considering the macroscopic yield and fracture stresses of approximately 900 MPa and less than 1.1 GPa for Hastelloy substrate-coated conductors. Here, we evaluate elastic-plastic conductor damage at over 0.4% strain using a high-stress REBCO coil exposed to a high field to explore this elastic-plastic regime. The coil was located off-center in a low-temperature superconductor magnet so as to induce a significant screening current in the enhanced radial field. Voltage taps, a Hall sensor, and two strain gauges were used for the instrumentation. We obtained strains exceeding 0.4% near the outward edge during the coil current charge from 350 A to 390 A, where the coil was exposed to external axial and radial magnetic fields of 13 T and 0.5 T. Post mortem results showed wavy plastic deformation, electrical damage, and REBCO defects. An elastic-plastic simulation reproduced the measured strains and predicted that ∼1 GPa stress is sufficient to induce ∼0.9% strain, thus validating our initial concerns with purely elastic models. This paper provides our experimental and simulation results.

Bang, Jeseok (ORCID:0000000247616334)↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2023

Microstructure characterization and elastic-plastic self-consistent simulation studies of anisotropic deformation of β -tin

Uniaxial compression tests of 99.9% pure polycrystalline $β$-tin (Sn) were conducted at various strain rates (10 -3 /s, 10 -1 /s and 1/s) and temperatures (294 K and 193 K) to understand the effects of these variables on the stress–strain response, microstructure evolution, and stress relaxation behavior. Multiple Sn specimens were subjected to complex compressive loading/unloading/reloading paths at different strain rates. Specimens initially compressed at higher strain rates showed strain rate-dependent texture evolution and more pronounced relaxation upon reloading, as compared to those pre-strained at lower rates at room temperature (RT). Compression tests conducted at low temperature (193 K) revealed increased strength, similar to the enhanced strength observed when the strain rate was increased at RT. Neutron diffraction was employed to characterize the initial and final bulk textures of the RT specimens. Electron backscatter diffraction was utilized to examine the crystallographic grain orientation and morphology, thereby identifying the signatures of dislocation-mediated deformation, grain refinement, recrystallization, and twinning behaviors. Elastic-plastic self-consistent simulations were performed to investigate the deformation modes responsible for the strain rate-dependent macroscopic stress–strain response and texture evolution. Effects of crystallographic orientation on stress relaxation behavior was also examined. In conclusion, the model predictions are in reasonable agreement with experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Non-linear elastic-plastic behavior of the invert glass lithium phosphorous oxynitride (LiPON)

Micro and nano-scale mechanical behavior of network/ionic glasses is dictated by their composition and the number of constraints per structural unit. From this perspective, the glass LiPON presents the opportunity for enhancement of the microscale ductility due to its rather unconstrained orthophosphate structure. Here, we use instrumented nanoindentation with different tip geometries to investigate the mechanical response of LiPON glass. The results reveal that the elastic modulus of LiPON is not constant and depends on pressure. With the method utilizing spherical nanoindentation and continuous stiffness measurement (CSM) we determine the yield point of LiPON. We propose the Drucker-Prager type of yield criterion for LiPON and estimate its yield stress in compression as 2.4 GPa. There exists another stress threshold, however, around 490 MPa, at which the elastic deformation becomes non-linear, and this can be mistaken for the yield stress.

Constitutive model↗
DOE OSTI2023

FY23 Simulation of Elastic-Plastic Failure Propagation

This memo summarizes the simulation of ductile failure propagation work conducted under the ASC project “V&V of Ductile Failure” conducted during FY 23. Physically, the failure propagation consists of crack propagation in the material. In the numerical setting—specifically in a finite element model—propagation can be accomplished through element death when critical conditions occur locally at an element that is then deleted from the simulation. The validation of the finite element models is evaluated by direct comparison between the experimental and simulation results regarding the rate of crack growth and its influence on the load-deflection response of the specimens tested. This work considers two geometries that display stable crack propagation under displacement-controlled conditions. The first geometry consists of hat specimens loaded in compression with nominally identical geometries but made with three different materials: Steel A286, Al 7075-T651 and 304L stainless steel. The three materials represent a range of ductility values that affect the response and crack propagation within the specimen. The crack induced propagates under an essentially mode-II type of deformation. The second geometry consists of a pre-cracked 304L stainless steel compact tension test specimen loaded so as to induce a mode-I deformation at the crack.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Brittle failure analysis and modeling of high-burnup PWR fuel cladding alloys

The aim of this research is the development of methods for predicting mechanical behavior and identification of limiting conditions to prevent brittle failure of high-burnup (HBU) pressure water reactor (PWR) fuel cladding alloys. A finite element (FE) model of the ring compression test (RCT) was created to analyze the failure behavior of zirconium-based alloys with radial hydrides during the RCT. An elastic-plastic material model describes the zirconium alloy. The stress-strain curve needed for the elastic-plastic material model was derived by inverse finite element analyses. Cohesive zone modeling is used to reproduce sudden load drops during RCT loading. Based on the failure mechanism in non-irradiated ZIRLO (R) claddings, a micro-mechanical model was developed that distinguishes between brittle failure along hydrides and ductile failure of the zirconium matrix. Two different cohesive laws representing these types of failure are present in the same cohesive interface. The key differences between these constitutive laws are the cohesive strength, the stress at which damage initiates, and the cohesive energy, which is the damage energy dissipated by the cohesive zone. Statistically generated matrix-hydride distributions were mapped onto the cohesive elements and simulations with focus on the first load drop were performed. Computational results are in good agreement with the RCT results conducted on high-burnup M5 (R) samples. It could be shown that crack initiation and propagation strongly depend on the specific configuration of hydrides and matrix material in the fracture area.

Simbruner, Kai↗
DOE OSTI2024

A model for oblique impacts on material surfaces

Many practical situations of material damage, wear, and erosion involve collisions between small particles and surfaces at inclined angles. While there are many well-validated models of normal incidence impact situations, elastic-plastic models for oblique incidence impact events are lacking. Here the finite element method is used to predict the normal and tangential coefficient of restitution in oblique impacts for hard, elastic spheres impacting an elastic-perfectly plastic material surface. The proposed model covers various impact angles ranging from 0° to 45°, within a limiting impact velocity below which the effects of heating are negligible. The normal coefficient of restitution follows power-laws with respect to normalized values of the impact velocity. Interestingly, the tangential coefficient of restitution follows a linear relationship with impact velocity. Together, these results provide a semi-empirical set of equations predicting oblique impact rebounds (both velocity and trajectory) for a wide range of conditions and material properties, with which experimental results can be rapidly interpreted. As a result, Laser-Induced Particle Impact Test (LIPIT) data are also presented for aluminum particles impacting aluminum substrates, at impact angles of 25° and 40°; the results compare favorably with the model and validate the general use of such models for the analysis of experimental data.

Coefficient of restitution↗
DOE OSTI2023

A parametric finite element study for determining burst strength of thin and thick-walled pressure vessels

To accurately predict the burst strength of both thin and thick-walled pressure vessels (PVs), a parametric study of PV burst strength was performed for a wide range of vessel geometries and materials using elastic-plastic finite element analysis (FEA). A valid FEA model was established through a detailed study of 2D versus 3D FEA models, the critical stress failure criterion versus the limit load criteria, and the thick-wall effect on the FEA simulations. Here, the results show that the stresses and strains at the mean diameter, rather than outside diameter, determines a more accurate burst strength for both thin and thick-walled PVs. On this basis, a parametrized FEA script using the ABAQUS Python application programming interface (API) was used to create a large database of PV burst strengths for a variety of vessel geometries and materials, demonstrating that Python scripting is a powerful technique for performing parametric studies or generating large databases. From the FEA results, using the regression method, a new burst pressure model was developed as a function of the vessel geometry (D/t ratio) and material properties (UTS and n). As validated by a large number of full-scale burst test data, the proposed burst model can very accurately predict the burst strength for both thin and thick-walled PVs.

42 ENGINEERING↗
DOE OSTI2025

Burst pressure models and validations for thick-walled pipelines containing corrosion defects

Corrosion is one major threat to pipeline integrity. Over the past decades, many corrosion models have been developed for determining the remaining strength of corroded pipelines, including ASME B31.G, Modified B31.G, LPC, PCORRC and their modified models. All these corrosion models are applicable only to large diameter, thin-walled pipelines with a diameter to wall thickness ratio D/t ≥ 20. In practice, many pipelines have a small diameter and thick wall with a D/t ratio < 20, and thus an adequate corrosion model is needed for assessing remaining strength for corroded thick-walled pipelines. This paper briefly reviews the theoretical burst pressure models for defect-free thin and thick-walled pipelines and four representative corrosion assessment models for thin-walled corroded pipelines. On this basis, two modified corrosion models are proposed to thick-walled pipelines in terms of the average shear stress yield theory. To verify the proposed corrosion models, comprehensive validations are performed. Numerical validations include the elastic-plastic finite element analysis to determine burst pressure for pipelines without and with corrosion defects and the model evaluation using a large dataset of available FEA results of burst pressure for machined defects. Experimental validations include a set of burst pressure tests for defect-free thick-walled pipes with different thicknesses and the model evaluation using one large burst dataset for machined defects with flat bottoms and another large dataset for real corrosion defects with curved river bottom profiles. Both numerical and experimental validations show that the proposed corrosion models can more accurately predict the remaining strength for corroded thin and thick-walled pipelines.

Pipeline↗
DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ExactPack

A new version of the verification tool ExactPack, 1.7.0, includes three new solvers: A Radshocks solver, Riemann solver and an Elastic-Plastic Piston problem solver. The code has been refactored to work with Python 3. All Fortran code has been removed and solvers implemented in Fortran have been re-implemented in Python where necessary. All unit tests have been converted to use the Pytest format and all the GUI, CLI and verification analysis code has been removed.

Thrussell, Jasper↗
DOE OSTI2024

Fracture resistance of vintage cast iron in gaseous hydrogen

In an effort to decarbonize legacy energy systems, several projects around the world are exploring alternatives to natural gas. Gaseous hydrogen is proposed as a carbon-free fuel to displace natural gas in existing legacy natural gas distribution systems, some of which continue to operate after 100 years (or more) in service. These systems, particularly in older industrial centers, contain cast iron pipe. However, the fracture resistance of most metals is degraded in gaseous hydrogen environments. This study evaluated the fracture resistance of several legacy cast iron pipe materials while exposed to gaseous hydrogen. Measurements were performed in three environments: air, a gas mixture with hydrogen partial pressure of 1 bar and pure hydrogen with a partial pressure of 34 bar. Although cast iron is generally considered a low ductility metal, elastic-plastic fracture methods are needed to assess the fracture resistance of the relatively small specimens that can be extracted from legacy pipe. Hydrogen reduced the fracture resistance of these cast iron materials by 10-40%. In air, the fracture resistance was determined to be as high as 21 MPa m 1/2 , whereas in gaseous hydrogen at pressure of about 1 bar the fracture resistance was as low as 13 MPa m 1/2 . Additional modest degradation of the fracture resistance was assessed at higher partial pressure (as low as 12 MPa m 1/2 ).

08 HYDROGEN↗
DOE OSTI2024

Calibration of Al7075 with Plate-Puncture Predictions

The following details calibration of a material model for Al7075-T6511. This aluminum alloy is commonly used across a host of engineering applications. Owing to its widespread prevalence, there is great benefit in improving simulation predictions for this alloy. In the present effort, a calibration is performed of its elastic-plastic response accounting for both rate and temperature dependence. The calibration is informed by a series of tests that include specimens of different geometries tested at different rates and temperatures. All specimens are derived from the same barstock, 3.5 inches in diameter. The fitted model itself uses an anisotropic, Hill yield surface coupled with a Johnson-Cook hardening model. Failure predictions are had by means of a modified Wilkins failure criterion. Following calibration of the material model, a validation exercise is performed against platepuncture experiments. These experiments include multiple probe shapes, probe diameters, and plate thicknesses. The puncture experiments are replicated in simulation with mesh studies performed to assess uncertainty. Key quantities of interest, notably the absorbed energy up to failure, are compared between simulation and experiment providing a means to assess the suitability of the calibration in puncture simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Develop a Fast Analysis Solver for Welding Sequence Optimization

During the shipbuilding manufacturing process, materials are exposed to significant stresses, as induced both thermally and mechanically, that alter the intended design and significantly affect the production schedule, labor hours (fitting, welding, rework, etc.), and material structural performance. The type and magnitude of deformation of a given structure depends on many factors such as the material, thickness and quality of components, the process heat input, preheat and inter-pass temperatures, type and size of welds, welding sequence and direction, location, sequence, and degree of fixturing. Numerical simulations using finite element analysis (FEA) have long been used to analyze welding-induced structural distortion. For large assemblies, transient thermal elastic-plastic analysis (TEPA) can take days or weeks to run, and optimization of welding sequence is not feasible. Simplified analysis methods were developed to reduce computational time. However, it is challenging to use these techniques to fully optimize welding sequencing because of their applied simplifications in modeling weld details. A fast analysis solver that could be used by the shipbuilding industry is being developed for optimizing welding sequences by taking full advantage of modern GPU-based HPC hardware and incorporating patented acceleration schemes. The accelerated processing factors are up to 2200 times greater for large, multi-pass welded structures.

Yang, Yu-Ping↗
DOE OSTI2023

Robust deep learning framework for constitutive relations modeling

Modeling the full-range deformation behaviors of materials under complex loading and materials conditions is a significant challenge for constitutive relations (CRs) modeling. Here, we propose a general encoder-decoder deep learning framework that can model high-dimensional stress-strain data and complex loading histories with robustness and universal capability. The framework employs an encoder to project high-dimensional input information (e.g., loading history, loading conditions, and materials information) to a lower-dimensional hidden space and a decoder to map the hidden representation to the stress of interest. We evaluated various encoder architectures, including gated recurrent unit (GRU), GRU with attention, temporal convolutional network (TCN), and the Transformer encoder, on two complex stress-strain datasets that were designed to include a wide range of complex loading histories and loading conditions. All architectures achieved excellent test results with an root-mean-square error (RMSE) below 1 MPa. Additionally, we analyzed the capability of the different architectures to make predictions on out-of-domain applications, with an uncertainty estimation based on deep ensembles. The proposed approach provides a robust alternative to empirical/semi-empirical models for CRs modeling, offering the potential for more accurate and efficient materials design and optimization.

36 MATERIALS SCIENCE↗