SEARCH · Engineering Papers

Results for “Linear elasticity”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗

DOE OSTI2025

Modelling the brittle failure of graphite induced by the controlled impact of runaway electrons in DIII-D

The thermo-mechanical response of an ATJ graphite sample to controlled runaway electron (RE) dissipation, realized in DIII-D, is modelled with a novel work-flow that features the RE orbit code KORC, the Monte Carlo particle transport code Geant4 and the finite element multiphysics software COMSOL. KORC provides the RE striking positions and momenta, Geant4 calculates the volumetric energy deposition and COMSOL simulates the thermoelastic response. Brittle failure is predicted according to the maximum normal stress criterion, which is suitable for ATJ graphite owing to its linear elastic behavior up to fracture and its isotropic mechanical properties. Measurements of the conducted energy, damage topology, explosion timing and blown-off material volume, impose a number of empirical constraints that suffice to distinguish between different RE impact scenarios and to identify RE parameters which provide the best match to the observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Asymmetric equilibrium core structures of pyramidal-II < c + a > dislocations in ten hexagonal-close-packed metals

The structures of pyramidal-II < c + a > dislocations, one of the most important defects in structural hexagonal-close-packed (HCP) metals, have not been fully characterized for many of the HCP metals in use today. Here, we employ ab initio informed phase-field dislocation dynamics to determine the minimum energy structure of pyramidal {1¯1¯22} < 11¯23 > dislocations in ten HCP metals, including Be, Co, Mg, Re, Ti, Zn, Cd, Hf, Y, and Zr. As input for the simulations, we calculate, using first-principles density functional theory, the {1¯1¯22} generalized stacking fault energy (GSFE) curves for all ten metals. From these calculations, it is found that magnetism in Co is necessary for achieving a local minimum in the GSFE curve. We observe in simulations that edge and screw character dislocations split into two partials separated by a low-energy intrinsic stacking fault. The splitting distance is shown to scale inversely with the local minimum energy normalized by the product of its shear modulus and Burgers vector. Interestingly, some HCP metals exhibit an asymmetric structure, with either unequal partial Burgers vectors or widths, in contrast to the symmetric configuration expected from linear elastic dislocation theory. We explain these structures by properties of the local maxima in their GSFE curves. Metals with larger degrees of elastic anisotropy result in dislocations with larger splitting distances than would be expected under the commonly used assumption of elastic isotropy. Furthermore, these findings on the sizes and asymmetry in the structures of pyramidal-II < c + a > dislocations are fundamental to understanding how these dislocations glide and interact or react with other defects when these metals are mechanically strained.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Viscoelastic Material Characterization and Modeling of Photovoltaic Module Packaging Materials for Direct Finite-Element Method Input

Numerical tools, such as the finite-element method, are increasingly used to design and evaluate the photovoltaic (PV) modules, providing for the reduction of development time and improved performance and reliability. However, high-fidelity material models are necessary to accurately model the complex structural behavior of the involved packaging materials. A common simplification used in recent years is to model the polymer materials (i.e., encapsulant and backsheet) as linear elastic, which will lead to inaccurate results. Therefore, in this work, we present a thorough characterization of the time- and temperature-dependent mechanical response of predominant PV module encapsulant and backsheet materials. Based on this material characterization, we developed and experimentally validated generalized Maxwell models to describe each material's viscoelastic response. In addition, we included measurements of the coefficient of thermal expansion and presented all material models in such a fashion for direct input into commercial finite-element method modeling software.

FEM↗

DOE OSTI2020

Fracture Failure Analysis for MQXFA Magnet Aluminum Shells

The High-Luminosity LHC Accelerator Upgrade Project (HL-LHC AUP) is approaching the production phase of the US-contributed Q1 and Q3 Interaction Region Quadrupoles (MQXFA). The structures for the MQXFA prototypes were design and inspected by the US-LARP (LHC Accelerator Research Program), AUP developed criteria, which will be used for the pre-series structures. As the first two full-length prototypes with 4.2 m magnetic length, MQXFAP1 and MQXFAP2, were designed and assembled at Lawrence Berkeley National Laboratory (LBNL), and tested at Brookhaven National Laboratory (BNL). The end aluminum short shell of MQXFAP2 was fractured along the shell length during the test, and tests were stopped. Analytical and Finite Element analysis were performed in light of the graded procedure defined in the Structure Design Criteria to investigate the fracture failure for MQXFAP2. In this paper, we report the fracture analysis of the current shell design, including the elasto-plastic simulations with submodel technique, and calculations with Linear Elastic Fracture Mechanics (LEFM). Test material properties are also presented. The results of this analysis explain why the end shell of MQXFAP2 failed, and suggest fillets on the end shell notches to meet the margin specified in the Structural Design Criteria.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

On the bulk compaction of brittle granular materials, Part II: Representing interparticle comminution by K COMP

An applied and empirical compaction fracture resistance parameter, K COMP , was developed to represent the compression-induced interparticle-comminution response of granular brittle materials (e.g., ceramic or glass particles). The development of K COMP is an outcome from Part I of this three-paper series. The K COMP represents a continuum response of compaction fracture resistance and was adapted from established Griffith linear elastic fracture mechanics theory. From that, the K COMP relates macroscopically applied compressive stress to a corresponding inverse-square-root representative particle size, where the latter is estimated from an entire particle size distribution or specific surface area measurement. In conclusion, the K COMP value was not constant over a wide range of stresses for all examined brittle granular materials, and this is indicative of a change in the dominant mode of permanent deformation at higher compaction stresses.

brittle↗

DOE OSTI2022

Calibrating constitutive models with full‐field data via physics informed neural networks

Abstract The calibration of solid constitutive models with full‐field experimental data is a long‐standing challenge, especially in materials that undergo large deformations. In this paper, we propose a physics‐informed deep‐learning framework for the discovery of hyperelastic constitutive model parameterizations given full‐field surface displacement data and global force‐displacement data. Contrary to the majority of recent literature in this field, we work with the weak form of the governing equations rather than the strong form to impose physical constraints upon the neural network predictions. The approach presented in this paper is computationally efficient, suitable for irregular geometric domains, and readily ingests displacement data without the need for interpolation onto a computational grid. A selection of canonical hyperelastic material models suitable for different material classes is considered including the Neo–Hookean, Gent, and Blatz–Ko constitutive models as exemplars for general non‐linear elastic behaviour, elastomer behaviour with finite strain lock‐up, and compressible foam behaviour, respectively. We demonstrate that physics informed machine learning is an enabling technology and may shift the paradigm of how full‐field experimental data are utilized to calibrate constitutive models under finite deformations.

Hamel, Craig M.↗

DOE OSTI2020

A Generalized Stress Inversion Approach With Application to Residual Stress Estimation

In this work develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.

42 ENGINEERING↗

DOE OSTI2022

Designing Cladded Components for High Temperature Nuclear Service: Part I—Analysis Methods

Use of corrosion-resistant cladding can greatly extend the design life of structural components in many advanced reactor systems. However, there are currently no ASME design rules for cladded components to guard against creep-fatigue failure and ratcheting strain accumulation in elevated temperature nuclear service. This paper, presented in two parts, addresses this gap by proposing a design strategy for cladded components that do not require long-term testing of clad materials. The proposed approach relies on approximate design analysis methods for two types of clad materials—a soft clad that creeps much faster than and has lower yield stress than the class A base material and a hard clad that creeps much slower than and has higher yield stress than the class A base material. Part I discusses the approximate analysis strategies for the clad materials—treat a soft clad as perfectly compliant and a hard clad as linear elastic—and Part II develops a complete set of design rules for each of the two types of cladded components. In conclusion, finite element analyses of representative high temperature reactor components show that the proposed design analysis methods can bound the design quantities in soft cladded components and approximate the design quantities in hard cladded components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Designing Cladded Components for High Temperature Nuclear Service: Part II—Design Rules

The challenge of using existing ASME Section III, Division 5, class A metallic materials for the construction of structural components of advanced reactors with corrosive coolants could be mitigated by allowing designers to use cladding to protect the base material from corrosion. However, the existing Section III, Division 5 rules provide no guidance on the evaluation of strain accumulation and creep-fatigue damage in cladded components. The availability of design rules for cladded components that do not require long-term clad material testing could promote the application of the cladding approach to accelerate the deployment schedule of these advanced reactor systems. To avoid long-term properties for the clad materials Part I of this work proposes two approximate design analysis methods for two types of clad materials—soft clads that creep much faster and have lower yield stress than the class A base material, and hard clads that creep much slower and have higher yield stress than the class A base material. The proposed analysis methods approximate the response of a soft clad by treating it as perfectly compliant and of a hard clad by treating it as linear elastic. Based on these approximate design analysis strategies this Part II develops a complete set of design rules for class A components cladded with either soft or hard clad materials. In conclusion, Part II discusses the reasoning behind the proposed design rules and uses example finite element analyses of representative reactor components to illustrate the use of these design methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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 OSTI2024

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗

DOE OSTI2020

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

49 EE - Wind and Water Power Program - Wind (EE-4W↗

DOE OSTI2022

Analyses of internal structures and defects in materials using physics-informed neural networks

Characterizing internal structures and defects in materials is a challenging task, often requiring solutions to inverse problems with unknown topology, geometry, material properties, and nonlinear deformation. Here, we present a general framework based on physics-informed neural networks for identifying unknown geometric and material parameters. By using a mesh-free method, we parameterize the geometry of the material using a differentiable and trainable method that can identify multiple structural features. We validate this approach for materials with internal voids/inclusions using constitutive models that encompass the spectrum of linear elasticity, hyperelasticity, and plasticity. We predict the size, shape, and location of the internal void/inclusion as well as the elastic modulus of the inclusion. Our general framework can be applied to other inverse problems in different applications that involve unknown material properties and highly deformable geometries, targeting material characterization, quality assurance, and structural design.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Reducing Operator Complexity of Galerkin Coarse-grid Operators with Machine Learning

Here, we propose a data-driven and machine-learning-based approach to compute non-Galerkin coarse-grid operators in multigrid (MG) methods, addressing the well-known issue of increasing operator complexity. Guided by the MG theory on spectrally equivalent coarse-grid operators, we have developed novel machine learning algorithms that utilize neural networks combined with smooth test vectors from multigrid eigenvalue problems. The proposed method demonstrates promise in reducing the complexity of coarse-grid operators while maintaining overall MG convergence for solving parametric partial differential equation problems. Numerical experiments on anisotropic rotated Laplacian and linear elasticity problems are provided to showcase the performance and comparison with existing methods for computing non-Galerkin coarse-grid operators.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Dynamic Response of a Semiactive Suspension System with Hysteretic Nonlinear Energy Sink Based on Random Excitation by means of Computer Simulation

This paper aims to investigate the property and behavior of the hysteretic nonlinear energy sink (HNES) coupled to a half vehicle system which is a nine-degree-of-freedom, nonlinear, and semiactive suspension system in order to improve the ride comfort and increase the stability in shock mitigation by using the computer simulation method. The HNES model is a semiactive suspension device, which comprises the famous Bouc–Wen (B-W) model employed to describe the force produced by both the purely hysteretic spring and linear elastic spring of potentially negative stiffness connected in parallel, for the half vehicle system. Nine nonlinear motion equations of the half vehicle system are derived in terms of the seven displacements and the two dimensionless hysteretic variables, which are integrated numerically by employing the direct time integration method for studying both the variables of vertical displacements, velocities, accelerations, chassis pitch angle, and the ride comfort and driver safety, respectively, based on the bump and random road inputs of the pseudoexcitation method as excitation signal. Simulation results show that, compared with the HNES model and the magnetorheological (MR) model coupled to the half vehicle system, the ride comfort and stability have been evidently improved. A successful validation process has been performed, which indicated that both the ride comfort and driver safety properties of the HNES model coupled to half vehicle significantly improved.

Chen, Hui↗

DOE OSTI2024

Tardigrade, Version 1.1x

Tardigrade is a collection of functions, libraries, and stand-alone code which enables the simulation of materials using the higher-order micromorphic framework. The tools present are the interface to INL's MOOSE code which is used as the underlying FEA solver, several constitutive equations (linear elasticity, a pressure-sensitive elasto-plastic model, crystal plasticity, and additional models for geo-materials) useful for demonstration purposes and the simulation of some materials, the micromorphic filter which converts DNS information to micromorphic stress-deformation quantities, and the overlap coupling framework which enables multi-scale simulations of a coupled higher-order macroscale and a classically defined microscale.

Miller, Nathan↗

DOE OSTI2023

Canonical and Noncanonical Hamiltonian Operator Inference

SAND2024-08577O Canonical and Noncanonical Hamiltonian Operator Inference is research code that applies old and new model reduction algorithms to well-known and publicly available "toy" problems. It applies operator inference (OpInf) methods for model reduction to several benchmark hyperbolic partial differential equations: linear wave, Korteweg-de Vries, Benjamin-Bona-Mahoney, and linear elasticity. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗