Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic strain engineering”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Phase-field modeling of orientation-dependent crack growth in ductile single crystals with anisotropic elasticity

Crack growth in ductile single crystals (DuSCs) is orientation dependent due to the anisotropies of crystal plasticity and elastic tensor. This study develops a phase-field model incorporating both crystal plasticity and crack growth and proposes a general method to decompose the elastic energy into compressive and tensile parts to prevent crack growth under compression in the phase-field description. The phase-field model, in combination with three Euler angles, is employed to simulate orientation-dependent crack growth in DuSCs. The contributions from crystal plasticity and anisotropic elasticity are compared, and the former is found to dominate in the anisotropy of crack growth in copper single crystals. Furthermore, the simulation results demonstrate that crystal orientation strongly affects the heterogeneous distribution of plastic strain and the interaction between plastic strain and crack growth. High-throughput phase-field simulations are performed with exhaustive crystal orientations, and the results are explained based on the anisotropy of the Taylor factor.

Computational Solid Mechanics↗
DOE OSTI2024

Reversible two-way tuning of thermal conductivity in an end-linked star-shaped thermoset

Polymeric thermal switches that can reversibly tune and significantly enhance their thermal conductivities are desirable for diverse applications in electronics, aerospace, automotives, and medicine; however, they are rarely achieved. Here, we report a polymer-based thermal switch consisting of an end-linked star-shaped thermoset with two independent thermal conductivity tuning mechanisms—strain and temperature modulation—that rapidly, reversibly, and cyclically modulate thermal conductivity. The end-linked star-shaped thermoset exhibits a strain-modulated thermal conductivity enhancement up to 11.5 at a fixed temperature of 60 °C (increasing from 0.15 to 2.1 W m –1 K –1 ). Additionally, it demonstrates a temperature-modulated thermal conductivity tuning ratio up to 2.3 at a fixed stretch of 2.5 (increasing from 0.17 to 0.39 W m –1 K –1 ). When combined, these two effects collectively enable the end-linked star-shaped thermoset to achieve a thermal conductivity tuning ratio up to 14.2. Moreover, the end-linked star-shaped thermoset demonstrates reversible tuning for over 1000 cycles. The reversible two-way tuning of thermal conductivity is attributed to the synergy of aligned amorphous chains, oriented crystalline domains, and increased crystallinity by elastically deforming the end-linked star-shaped thermoset.

42 ENGINEERING↗
DOE OSTI2021

Anomalous elastic behavior of tantalum at high pressures: Experimental and theoretical studies

Elastic wave velocities of polycrystalline tantalum have been measured up to 13.6 GPa at room temperature using ultrasonic interferometry technique in a multi-anvil apparatus. The bulk and shear moduli, as well as their pressure derivatives are obtained using Eulerian 3rd order finite strain equations, yielding K S0 = 193.9(29) GPa, G 0 = 69.3(10) GPa, K S0' = 3.18(5) and G 0' = 0.88(1). Here the Steinberg-Guinan yield stress model is examined using the current shear modulus and its pressure derivative, where a slightly lower A 0 = 0.0127 GPa –1 compared to previous studies is obtained. First-principles calculations on the elasticity of tantalum have also been performed up to 200 GPa. Anomalous softening of the elastic shear constant C 44 and shear wave velocity V S is observed between 80– 180 GPa, whereas the anisotropy shows a softening-stiffening behavior with increasing pressure and reaches a minimum at around 150 GPa. Such anomalies might commonly originate from the nesting of Fermi surface for Group VB metals.

36 MATERIALS SCIENCE↗
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 OSTI2020

Localized Strain Analysis of Ce- and Mg-Treated Cast Iron under Uniaxial Compression

Cast iron exhibits a wide range of mechanical properties, depending on its microstructural features. The microstructure of cast iron consists of several microconstituents with different elastic-plastic behavior, making the strain non-uniform across the bulk material. To understand the individual effects of these microconstituents on the overall mechanical behavior, local strain analysis using digital image correlation analysis was carried out. Samples with two different compositions (varying cerium, magnesium and silicon) were processed at different solidification velocities in a Bridgman furnace. Sections of the directionally solidified samples were loaded under uniaxial compression to measure global and local strain behavior. Despite the variability of the microstructure, the stress–strain curves obtained by digital image correlation (DIC) were found to react in a well-controlled way to changes in solidification velocity. It was observed that high-strain failure (greater than 15%) was accompanied by local straining of the softer ferritic phase, but during low-strain failure, local straining was not prominent. Higher nodularities, due to higher solidification velocities, raised the compressive strength without affecting the toughness significantly. Higher percentages of carbides led to higher compressive strengths with corresponding losses in ductility. The continuity of the matrix was also found to play an important role in the behavior during compression.

ductile failure↗
DOE OSTI2026

Compressive Response and Energy Absorption of Additively Manufactured Elastomers with Varied Simple Cubic Architectures

Additive manufacturing, and particularly the vat photopolymerization process, enables the fabrication of complex geometries at high resolution and small length scales, making it well-suited for fabricating cellular structures (e.g., foams and lattices). Among these, elastomeric cellular structures are of growing interest due to their tunable compliance and energy dissipation. However, comprehensive data on the compressive behavior of these structures remains limited, especially for investigating the structure-property effects from changing the density and distribution of material within the cellular structure. This study explores how the mechanical response of polyurethane-based simple cubic structures changes when varying volume fraction, unit cell length, and unit cell patterning, which have not been systematically investigated previously in additively manufactured elastomers. Increasing volume fraction from 10% to 50% yielded significant changes in compressive stress–strain performance (decreasing strain at 0.5 MPa by 41.6% and increasing energy absorption density by 3962.5%). Although changing the unit cell length between 2.5 and 7 mm in ~30 mm parts did not result in statistically different stress–strain responses, modifying the configuration of struts of different thicknesses across designs with 30% volume fraction altered the stress–strain behavior (differences of 12.5% in strain at 0.5 MPa and 109.4% for energy absorption density). Power law relationships were developed to understand the interactions between volume fraction, unit cell length, and elastic modulus, and experimental data showed strong fits (R 2 > 0.91). These findings enhance the understanding of how multiple structural design aspects influence the performance of elastomeric cellular materials, providing a foundation for informing strategic design of tailorable materials for diverse mechanical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Flexoelectricity and New Phenomena

It is easy to miss the scientific implications of our recent work on Triboelectricity. Everyone knows that rubbing and contact can produce static electricity; less appreciated is that the thermodynamic driver has been an open question since static electricity was first observed by Thales of Miletus around 585 BC. People think they understand it, for instance one common explanation that can still be found in the current literature is that differences in the work function drives charge transfer, often called the Volta-Helmholtz hypothesis. As summarized in 1967 by Harper, this fails to explain many experimental observations, for instance that charging can occur when two pieces of the same material are rubbed against each other. We are the first to place triboelectricity on a solid foundation rooted in quantum mechanics – the flexoelectric effect. We were able to explain a significant number of previously unexplained phenomena: 1) Bipolar tribocurrents associated with stick-slip, due to the change in sign of the strain gradients. 2) A one-third power scaling of tribocurrents with indentation force. 3) Tribocharging when two identical materials are used, these being due to local variations in the asperities so there are usually local potential differences. 4) Inhomogeneous charging of insulators, related to the statistical nature of asperities. 5) An experimentally observed reversal in the sign of charge transfer for negative and positive curvature, which is related to a change in the sign of the strain gradient. Exploiting our DOE prior funded work on flexoelectricity, we obtained semi-quantitative matching to existing experimental measurements of the surface charge in triboelectric experiments. The work has been well received in the literature. The work has been the focus of a number of popular science press articles, and also formed the basis for a Podcast for children 6-10 “The Rise and Fall of Static Man” posted in December 2019 by NPR as part of their “Wow in the World” series. I was also briefly interviewed by the Chicago PBS station in January 2020. This work are significant for a wide range of energy applications; to quote from an independent source: "Triboelectric power has plenty of potential, says Wenzhuo Wu, an assistant professor of engineering at Purdue. If the basics of static electricity are better understood, we could maximize the efficiency of wind or wave power generators, Wu says. The body's own movement could be used to power internal medical devices. Imagine being able to create a roof shaped to harness the power of a raindrop — the friction of the rain passing over the surface — to generate triboelectricity, powering the building below it." This is the start of new science, some of which we already partially understand such as the role of band bending in charge transfer. We need to understand charge transfer combining elasticity, quantum mechanics, band bending and defect states. These directly involve several of the DOE Grand Challenges: "How do we control material processes at the level of electrons? How do remarkable properties of matter emerge from complex correlations of the atomic or electronic constituents and how can we control these properties? How do we characterize and control matter away—especially very far away— from equilibrium?" I will argue that this work truly falls into the class of disruptive science; it is not just a simple extension, linear science. Not everyone will accept the approach. Since we explain far more about triboelectricity than anyone before, the preponderance of evidence supports the model. The feedback I have received is that many agree with the work, to quote: "The model makes sense, says Michael McAlpine, a professor of engineering at the University of Minnesota. "It's such a simple explanation, I was surprised I didn't put my finger on that," McAlpine says." The proposal received strong reviews. It was also publicized on the Department of Energy Web Site.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Bending-induced isostructural transitions in ultrathin layers of van der Waals ferrielectrics

Using Landau-Ginzburg-Devonshire (LGD) phenomenological approach we analyze the bending-induced re-distribution of electric polarization and field, elastic stresses and strains inside ultrathin layers of van der Waals ferrielectrics. We consider a CuInP2S6 (CIPS) thin layer with fixed edges and suspended central part, the bending of which is induced by external forces. The unique aspect of CIPS is the existence of two ferrielectric states, FI1 and FI2, corresponding to big and small polarization values, which arise due to the specific four-well potential of the eighth-order LGD functional. When the CIPS layer is flat, the single-domain FI1 state is stable in the central part of the layer, and the FI2 states are stable near the fixed edges. With an increase of the layer bending below the critical value, the sizes of the FI2 states near the fixed edges decreases, and the size of the FI1 region increases. When the bending exceeds the critical value, the edge FI2 states disappear being substituted by the FI1 state, but they appear abruptly near the inflection regions and expand as the bending increases. The bending-induced isostructural FI1-FI2 transition is specific for the bended van der Waals ferrielectrics described by the eighth (or higher) order LGD functional with consideration of linear and nonlinear electrostriction couplings. Furthermore, the isostructural transition, which is revealed in the vicinity of room temperature, can significantly reduce the coercive voltage of ferroelectric polarization reversal in CIPS nanoflakes, allowing for the curvature-engineering control of various flexible nanodevices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Revealing the unusual rate-dependent mechanical behaviors of nematic liquid crystal elastomers

Liquid crystal elastomers (LCEs) exhibit unique mechanical properties of soft elasticity and enhanced energy dissipation with rate dependency. They are potentially transformative materials for applications in mechanical impact mitigation and vibration isolation. However, previous studies have primarily focused on the mechanics of LCEs under equilibrium and quasistatic loading conditions. Critical knowledge gaps exist in understanding their rate-dependent behaviors, which are a complex mixture of traditional network viscoelasticity and the soft elastic behaviors with changes in the mesogen orientation and order parameter. Together, these inelastic mechanisms lead to unusual rate-dependent energy absorption responses of LCEs. In this work, we developed a viscoelastic constitutive theory for monodomain nematic LCEs to investigate how multiple underlying sources of inelasticity manifest in the rate-dependent and dissipative behaviors of monodomain LCEs. Here, the theoretical modeling framework combines the neo-classical network theory with evolution rules for the mesogen orientation and order parameter with conventional viscoelasticity. The model is calibrated with uniaxial tension and compression data spanning six decades of strain rates. The established 3D constitutive model enables general loading predictions taking the initial mesogen orientation and order parameter as inputs. Additionally, parametric studies were performed to further understand the rate dependence of monodomain LCEs in relation to their energy absorption characteristics. Based on the parametric studies, particularly loading scenarios are identified as conditions where LCEs outperform conventional elastomers regarding energy absorption.

42 ENGINEERING↗
DOE OSTI2021

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗
DOE OSTI2022

Non-destructive depth-resolved characterization of residual strain fields in high electron mobility transistors using differential aperture x-ray microscopy

Localized residual stress and elastic strain concentrations in microelectronic devices often affect the electronic performance, resistance to thermomechanical damage, and, likely, radiation tolerance. A primary challenge for the characterization of these concentrations is that they exist over sub-μm length-scales, precluding their characterization by more traditional residual stress measurement techniques. Here, in this study, we demonstrate the use of synchrotron x-ray-based differential aperture x-ray microscopy (DAXM) as a viable, non-destructive means to characterize these stress and strain concentrations in a depth-resolved manner. DAXM is used to map two-dimensional strain fields between the source and the drain in a gallium nitride (GaN) layer within high electron mobility transistors (HEMTs) with sub-μm spatial resolution. Strain fields at various positions in both pristine and irradiated HEMT specimens are presented in addition to a preliminary stress analysis to estimate the distribution of various stress components within the GaN layer. γ-irradiation is found to significantly reduce the lattice plane spacing in the GaN along the sample normal direction, which is attributed to radiation damage in transistor components bonded to the GaN during irradiation.

42 ENGINEERING↗
DOE OSTI2023

Type 2 diabetes impairs annulus fibrosus fiber deformation and rotation under disc compression in the University of California Davis type 2 diabetes mellitus (UCD-T2DM) rat model

Abstract Understanding the biomechanical behavior of the intervertebral disc is crucial for studying disease mechanisms and developing tissue engineering strategies for managing disc degeneration. We used synchrotron small-angle X-ray scattering to investigate how changes to collagen behavior contribute to alterations in the disc’s ability to resist compression. Coccygeal motion segments from 6-month-old lean Sprague-Dawley rats ( n=7) and diabetic obese University of California Davis type 2 diabetes mellitus (UCD-T2DM) rats ( n=6, diabetic for 68±7 days) were compressed during simultaneous synchrotron scanning to measure collagen strain at the nanoscale (beamline 7.3.3 of the Advanced Light Source). After compression, the annulus fibrosus was assayed for nonenzymatic cross-links. In discs from lean rats, resistance to compression involved two main energy-dissipation mechanisms at the nanoscale: (1) rotation of the two groups of collagen fibrils forming the annulus fibrosus and (2) straightening (uncrimping) and stretching of the collagen fibrils. In discs from diabetic rats, both mechanisms were significantly impaired. Specifically, diabetes reduced fibril rotation by 31% and reduced collagen fibril strain by 30% (compared to lean discs). The stiffening of collagen fibrils in the discs from diabetic rats was consistent with a 31% higher concentration of nonenzymatic cross-links and with evidence of earlier onset plastic deformations such as fibril sliding and fibril–matrix delamination. These findings suggest that fibril reorientation, stretching, and straightening are key deformation mechanisms that facilitate whole-disc compression, and that type 2 diabetes impairs these efficient and low-energy elastic deformation mechanisms, thereby altering whole-disc behavior and inducing the earlier onset of plastic deformation.

Rosenberg, James L.↗
DOE OSTI2022

Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods

Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.

42 ENGINEERING↗
DOE OSTI2021

Anisotropic elastoplastic response of double-porosity media

Here we present a continuum framework to simulate fluid flow through anisotropic elastoplastic media with double porosity. Two effective stress measures σ' and δ" emerge from the thermodynamic formulation, which are energy-conjugate to the elastic and plastic components of strain, respectively. Both effective stress measures can be expressed as a combination of the total Cauchy stress σ and the average pore pressure $\bar{p}$ in the two pore scales. In the effective stress for elasticity, $\bar{p}$ is scaled with a rank-2 Biot tensor, whereas the effective stress for plasticity follows the Terzaghi form in which $\bar{p}$ is scaled by the Kronecker delta. The Biot tensor and storage coefficients are derived as functions of elasticity parameters and porosities. A mixed finite element formulation is introduced to discretize the domain and solve initial boundary-value problems. A stabilization scheme is employed on equal-order interpolation for both displacement and pressure fields throughout the entire range of drainage responses. Numerical simulations reproduce the hydromechanical response of Opalinus shale in one-dimensional consolidation tests throughout the range of primary and secondary consolidation under different external loads. Numerical simulations of the consolidation of a rectangular domain subjected to a strip load demonstrate the efficacy of the proposed stabilization scheme, as well as illustrate the impacts of stress history, mass transfer, and different pore systems on the hydromechanical response.

42 ENGINEERING↗
DOE OSTI2025

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗
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↗