Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic deformation”

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

Tensile Modeling PVC Gels for Electrohydraulic Actuators

Polyvinyl chloride (PVC)-dibutyl adipate (DBA) gels are a fascinating dielectric elastomer actuator showing promise in soft robotics. When actuated with high voltages, the gel deforms towards the anode. A recent application of PVC gels in electrohydraulic actuators motivates elastic and hyperelastic constitutive relationships for tensile loading modes. PVC gels with plasticizer-to-polymer weight ratios of 2:1, 4:1, 6:1, and 8:1 w/w were evaluated. PVC gels exhibit a linear elastic region up to 25% strain. The elastic modulus decreased with increasing plasticizer content from 288.8 kPa, 56.1 kPa, 24.7 kPa, to 11 kPa. Poisson’s ratio also decreased with increasing plasticizer content from 0.42, 0.43, 0.39, to 0.35. We suggest that the decrease in polymer concentration facilitates a weakly interconnected polymer network susceptible to chain slippage that hinders the network response, thus lowering Poisson’s ratio. Our work suggests that PVC gels can be treated as isotropic and incompressible for large strains and hyperelastic modeling; however, highly plasticized gels tend to act less incompressible at small strains. The power scaling law between the elastic modulus and plasticizer weight ratio showed high agreement, making the elastic modulus deterministic for any plasticizer content. The Neo–Hookean, Mooney–Rivlin, Yeoh, Gent, Ogden, and extended tube hyperelastic constitutive models are investigated. The Yeoh model shows the highest feasibility when evaluated up to 3.5 stretch, showing a maximum normalized root-mean-square-error of 6.85%. Together, these findings establish a constitutive basis for PVC-DBA gels, incorporating small strain elasticity, large strain non-linear behavior, and network analysis while providing suggestive insight into the network structure required for accurately modeling the EPIC.

Polymer Science↗
DOE OSTI2020

Lattice strains and diffraction elastic constants of cubic polycrystals

In situ synchrotron X-ray and neutron diffraction experiments provide a powerful approach to measure lattice strains in bulk polycrystalline materials. They are being increasingly used for quantitative characterization of microscale deformation within and between grains and phases. Here we use a self-consistent micromechanics model to obtain a general analytic solution of the grain-level lattice strains and diffraction elastic constants for a broad class of elastically isotropic polycrystals with cubic crystal symmetry. This analytic solution reveals a direct linear relationship between the reciprocal of the elastic diffraction constant and the orientation index parameter along the direction of any diffraction vector, including tensile loading and transverse directions. The straightforward numerical implementation of this solution provides diffraction elastic constants for 26 representative cubic polycrystals. Finally, analytic solutions of this kind can serve to benchmark in situ diffraction measurements of lattice strains and also facilitate high-throughput studies of microscale stresses and diffraction elastic constants in polycrystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea-Ice Deformation Fields

As the sea-ice modeling community is shifting to advanced numerical frameworks, developing new sea-ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea-ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state-of-the-art sea-ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea-ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea-ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous-plastic (and the elastic-viscous-plastic variant), elastic-anisotropic-plastic, and Maxwell-elasto-brittle rheologies are compared in a single study. We find that both plastic and brittle sea-ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi-fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea-ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Energy based fracture initiation criterion for strain-crystallizing rubber-like materials with pre-existing cracks

Fracture prediction is indispensable for polymers, like rubbers, which have a broad range of applications mainly due to their high extensibility. The phenomenon known as strain-induced crystallization further contributes to the fracture toughness of certain rubbers. Here in this study, a criterion based on internal bond energy, incorporating the effects of crystallization, is proposed to predict fracture initiation in rubber-like materials with pre-existing cracks. First, a multi-scale mechanical model is developed for characterizing the behavior of rubber when subjected to both uniaxial and biaxial deformation states. At the microscale, both the amorphous and crystalline chain segments are modeled as elastic in order to consider the energy contribution by the molecular bond distortions. This internal energy is considered along with the entropic and crystalline free energy for each chain. In the chain model, the effects of loading condition and the relative orientation of a chain on its crystallinity are taken into account. At the macroscopic scale, an existing crystallinity distribution function is adapted and a mixed finite element formulation with an augmented Lagrangian multiplier is utilized to impose the incompressibility constraint. A non-affine maximal advance path constraint based homogenization model is utilized for bridging the two scales. Its potential to account for anisotropy in the stretched network compels the model to be preferable due to its physical significance, for the purpose of fracture modeling. The rigidity of the crystallites is accounted for by proposing a crystallite distortion energy in addition to the critical bond dissociation energy, for fracture initiation to occur. The model is validated by comparison with existing experimental results for both crystallizing and non-crystallizing rubbers. In addition to its potential to predict the material behavior when subjected to uniaxial and biaxial loading, the capability of the model to quantitatively estimate the effect of crystallization on fracture initiation is also verified.

36 MATERIALS SCIENCE↗
DOE OSTI2026

High-temperature mechanical properties of a γ′-strengthened Co-based superalloy designed for additive manufacturing

GammaPrint®-700 is a recently developed high γ′ (∼70% volume fraction) CoNi-based superalloy designed to combine high-temperature mechanical performance with laser powder bed fusion processability. Room-temperature yield strength ranged from 610 to 658 MPa and increased to 661 MPa (longitudinal) and 730 MPa (transverse) at 760 °C. The creep behavior compared favorably to high-γ′ Ni-based superalloys manufactured via laser powder bed fusion such as Incoloy® 939 and Inconel® 738LC. In-situ neutron diffraction data measured during creep revealed that plastic deformation was largely localized to the γ phase at 760 °C, allowing the γ’ phase to elastically compensate and maintain creep strain resistance. Furthermore, at 900 °C, this load sharing behavior was weakened, contributing to accelerated creep rate and rupture.

Creep↗
DOE OSTI2021

Pseudoelastic deformation in Mo-based refractory multi-principal element alloys

Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal-element alloys. We present phase and property analysis of quinary (MoW) x Zr y (TaTi) 1-x-y refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular dynamics simulations. Both CALPHAD and density-functional theory analysis of phase stability indicates a Mo-W-rich region of this quinary has a stable single-phase body-centered-cubic structure. We report first quinary composition from Mo-W-Ta-Ti-Zr family of alloy with pseudo-elastic behavior, i.e., hysteresis in stress-strain. Our analysis shows that only Mo-W-rich compositions of Mo-W-Ta-Ti-Zr, i.e., Mo+W, show reproducible hysteresis in stress-strain responsible for pseudo-elastic behavior. The (MoW) 85 Zr 7.5 (TaTi) 7.5 was down-selected based on temperature-dependent phase diagram analysis and molecular dynamics simulations predicted elastic behavior that reveals twinning-assisted pseudoelastic behavior. While mostly unexplored in body-centered-cubic crystals, twinning is a fundamental deformation mechanism that competes against dislocation slip in crystalline solids. Furthermore, this alloy shows identical cyclic deformation characteristics during uniaxial < 100 > loading, i.e., the pseudoelasticity is isotropic in loading direction. Additionally, a temperature increase from 77 to 1,500 K enhances the elastic strain recovery in load-unload cycles, offering possibly control to tune the pseudoelastic behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sound Speeds of Solids from Ultrasonic Pulse Receiver Measurements

A common method of determining elastic material properties is utilizing an ultrasonic pulse receiver. This is a non-destructive test (NDT) causing no plastic deformation of the material. It requires a pulse generator, oscilloscope, and two transducers to measure the sound velocity of a material. Both longitudinal and shear sound velocities may be obtained by utilizing the appropriate transducer. Acoustic nondestructive testing methods are often used in manufacturing certification processes and for material defect detection in industry. These methods can also be utilized to directly measure the elastic moduli of a material without damaging the material. Overall, it is a low-cost experiment with minimal preparation, and is easy to use.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The revised elastic field of an edge dislocation

We modify the traditional linear-elastic field of an edge dislocation. The modifications stem from symmetry and energy requirements imposed in terms of embedded (deformed) coordinates. These requirements are satisfied if a line force is added to the dislocation field. The modification to the stress field is expressed by coefficients that are a function only of Poisson’s ratio. Qualitatively, the field of the dislocation is increased in the glide direction and decreased in the climb direction. There are small changes to the field of a screw dislocation that do not entail the addition of a line force. Effects of nonlinearity, anisotropic elasticity, and image forces are briefly considered.

Hirth, J. P.↗
DOE OSTI2022

Comparison of three state-of-the-art crystal plasticity based deformation twinning models for magnesium alloys

In magnesium alloys, deformation twinning and its interactions with dislocation slip are responsible for a sigmoidal shape stress–strain behavior and an asymmetrical tension–compression yield strength in magnesium alloys. The sensitivity of twinning to the underlying microstructure renders the crystal plasticity method the most commonly adopted modeling approach for magnesium-twinning. This paper compares three state-of-the-art crystal plasticity-based twinning models from the literature, namely the elastic-viscoplastic self-consistent twinning-detwinning (EVPSC-TDT) model, crystal plasticity finite element model based on enhanced predominate twin reorientation approach (CPFE-ePTR), and the crystal plasticity finite element model based on “discrete twinning” approach (CPFE-DT). A polycrystalline microstructure is simulated with all three methods to compare the resulting stress–strain curves and lattice strains to those from the experimentally measured data. All three methods showed the capability of capturing the experimental results with different levels of accuracy. Additionally, the EVPSC-TDT method avoids solving the finite element matrices and showed the highest computational efficiency. The CPFE-ePTR model shows a higher accuracy in capturing the lattice strain. The CPFE-DT relies on high-resolution finite element mesh and is much slower than the other two methods, but it captured the local deformation concentration and stress reversal phenomena near the twin band, which was not possible with the other two methods. Based on the comparison, guidance for the selection of the appropriate model based on the specific modeling target is provided in this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamic behavior of chemically tunable mechano-responsive hydrogels

Here, using theory and simulation, we model the mechanical behavior of gels that encompass loops and dangling chain ends. If the loops remain folded and dangling ends are chemically inert, then these topological features just serve as defects. If, however, the loops unfold to expose the hidden (“cryptic”) binding sites and the ends of the dangling chains are reactive, these moieties can form bonds that improve the gel's mechanical properties. For gels with a lower critical solubility temperature (LCST), we systematically switch on the possible unfolding and binding events. To quantify the resulting effects, we derive equations for the gel's equilibrium and dynamic elastic moduli. We also use a finite element approach to simulate the gel's response to deformation and validate the analytic calculations. Herein, we show that the equilibrium moduli are highly sensitive to the presence of unfolding and binding transitions. The dynamical moduli are sensitive not only to these structural changes, but also to the frequency of deformation. For example, when reactive ends bind to exposed cryptic sites at T = 29 °C and relatively high frequency, the storage shear modulus is 119% greater than the corresponding equilibrium value, while the storage Young's modulus is 109% greater than at equilibrium. These findings provide guidelines for tuning the chemical reactivity of loops and dangling ends and the frequency of deformation to tailor the mechano-responsive behavior of polymer networks.

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

Model of ramp compression of diamond from ab initio simulations

Ramp compression experiments characterize high-pressure states of matter at temperatures well below those present in shock compression. However, because temperature is typically not directly measured during ramp compression, it is uncertain how much heating occurs under these shock-free conditions. Here, we performed a series of ab initio simulations on carbon in order to match the density-stress measurements of Smith et al. [Smith et al., Nature (London) 511, 330 (2014)]. We considered isotropically as well as uniaxially compressed solid carbon in the diamond and BC8 phases, with and without defects, as well as liquid carbon. Our idealized model ascribes heating during ramp compression to an initially uniaxially compressed cell transforming isochorically into an isotropically (hydrostatic equivalent) compressed state having lower internal energy, hence higher temperature so as to conserve energy. Multiple such heating events can occur during a single ramp experiment, leading to higher temperatures than with isentropic compression. Comparison with experiments shows that heating alone does not explain the equation of state measurements on diamond, instead implying that a significant uniaxial stress component remains present at high compression. The temperature predictions of our ramp compression model remain to be verified with future laboratory measurements.

58 GEOSCIENCES↗
DOE OSTI2022

Cavitation in a soft porous material

Abstract We study the collapse and expansion of a cavitation bubble in a deformable porous medium. We develop a continuum-scale model that couples compressible fluid flow in the pore network with the elastic response of a solid skeleton. Under the assumption of spherical symmetry, our model can be reduced to an ordinary differential equation that extends the Rayleigh–Plesset equation to bubbles in soft porous media. The extended Rayleigh–Plesset equation reveals that finite-size effects lead to the breakdown of the universal scaling relation between bubble radius and time that holds in the infinite-size limit. Our data indicate that the deformability of the porous medium slows down the collapse and expansion processes, a result with important consequences for wide-ranging phenomena, from drug delivery to spore dispersion.

Leng, Yu↗
DOE OSTI2024

Ab initio translationally invariant nucleon-nucleus optical potentials

We combine the ab initio symmetry-adapted no-core shell model (SA-NCSM) with the single-particle Green's function approach to construct optical potentials rooted in first principles. Specifically, we show that total cross sections and phase shifts for neutron elastic scattering from a 4 He target with projectile energies between 0.5 and 10 MeV closely reproduce the experiment. In addition, we discuss an important new development that resolves a long-standing issue with spurious center-of-mass motion in the Green's function formalism for many-body approaches. Furthermore, the new development opens a path for first-principle predictions of cross sections for elastic scattering of single-nucleon projectiles, nucleon capture, and deuteron breakup reactions, feasible for a broad range of open-shell spherical and deformed nuclei in the SA-NCSM approach.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electronic density of states as the descriptor of elastic bond strength, ductility, and local lattice distortion in BCC refractory alloys

Although electronic density of states (DOS) is fundamental to materials properties, its general relationship to mechanical properties of alloys is not well established. In this paper, using density functional theory (DFT) calculations, we show that the electronic occupancy at the Fermi level, N(E f ), obtained from DOS is a key descriptor of alloy strength and ductility. Our comprehensive analysis of numerous body centered cubic (BCC) refractory high entropy alloys (RHEAs) shows an overwhelming correlation that low N(E f ) indicates strong bonds that have high stiffness resulting in high elastic constants. High bond stiffness indicates presence of covalent nature of bonds that are directional in nature resulting in resistance to deformation leading to high bulk (B) and shear (G) moduli. Consequently, N(E f ) provides a direct correlation to the tendency of alloy ductility evidenced in the Pugh ratio (G/B). As stiffer bonds result in lower local lattice distortion (LLD), N(E f ) are LLD are also found to be corelated which opens up a correlation to solid solution strengthening and yield strength. Thus, this work unveils fundamental correlations between N(E f ) and (1) elastic bond strength, (2) ductility, and (3) LLD. These correlations open opportunities for the design of high strength high ductile RHEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Isotropic seismic radiation from rock damage and dilatancy

SUMMARY We present semi-analytical and numerical results on isotropic radiation expected to be generated by seismic sources with damaged dilatant materials that sustain during failure an overall shear deformation. The results are derived for inclusions with different stiffness components, including those responsible for shear-volumetric coupling in dilatant material, and different geometries, in a surrounding elastic solid. The obtained relations indicate that the isotropic seismic potency and seismic moment components are proportional to the shear potency components multiplied by the aspect ratio between the thickness and length of the source zone. For a typical rupture zone with a thin tabular geometry, the expected small aspect ratio explains the general seismological observations of small isotropic radiation. However, isotropic radiation induced by dynamic dilation for short and thick sources is more pronounced and may reach for spherical sources sustaining shear deformation 18 per cent of the shear components. Sources that include also tensile components and explosions can produce larger isotropic radiation. The results provide rigorous relations that can be used to estimate different source properties in numerical simulations, experiments and field studies.

Geochemistry & Geophysics↗
DOE OSTI2023

Highly Deformable Rigid Glassy Conjugated Polymeric Thin Films

Wearable devices benefit from the use of stretchable conjugated polymers (CPs). Traditionally, the design of stretchable CPs is based on the assumption that a low elastic modulus (E) is crucial for achieving high stretchability. However, this research, which analyzes the mechanical properties of 65 CP thin films, challenges this notion. It is discovered that softness alone does not determine stretchability; rather, it is the degree of entanglement that is critical. This means that rigid CPs can also exhibit high stretchability, contradicting conventional wisdom. To inverstigate further, the mechanical behavior, electrical properties, and deformation mechanism of two model CPs: a glassy poly(3-butylthiophene-2,5-diyl) (P3BT) with an E of 2.2 GPa and a viscoelastic poly(3-octylthiophene-2,5-diyl) (P3OT) with an E of 86 MPa, are studied. Ex situ transmission X-ray scattering and polarized UV–vis spectroscopy revealed that only the initial strain (i.e., <20%) exhibits different chain alignment mechanisms between two polymers, while both rigid and soft P3ATs showed similarly behavior at larger strains. In conclusion, by challenging the conventional design metric of low E for high stretchability and highlighting the importance of entanglement, it is hoped to broaden the range of CPs available for use in wearable devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Attaining regularization length insensitivity in phase-field models of ductile failure

In this work, a cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

42 ENGINEERING↗
DOE OSTI2025

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗