Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stiffness”

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

Simulation-driven design optimization of reaction injection molding (RIM) process for polydicyclopentadiene (pDCPD): Minimizing cycle time, defects, and warpage

Replacing metal components in trucks, trailers, and buses with lightweight polymer composites is challenging due to high temperatures and complex manufacturing. The Reaction Injection Molding (RIM) process using Dicyclopentadiene (DCPD) resin offers a solution by producing robust parts with excellent stiffness, impact strength, and resistance properties. Simulations are essential for optimizing this process, predicting defects, and improving quality. However, most commercial software is tailored for thermoplastics, requiring thermoset users to generate their own datasets. In this study, a material data card for DCPD was developed to perform RIM simulations. Design of Experiments (DOE) was used to identify key factors affecting filling, curing, and warpage, aiming to minimize cycle time and defects. The simulations explored varying injection gate parameters (size, location, number) and process conditions (mold/resin temperature, injection/curing pressure). Results showed that gate design significantly impacts filling behavior and defects. A single central gate provided balanced flow with fewer defects, while two corner gates led to more defects. Additionally, lower injection pressure increased filling time, while higher mold temperature accelerated curing but led to more warpage. In conclusion, this optimization framework aims to enhance DCPD part performance and promote sustainable manufacturing by reducing waste and energy consumption.

42 ENGINEERING↗
DOE OSTI2022

Magnetoelastic properties of multiferroic hexagonal ErMnO 3

Here, the strength and dynamics of magnetoelastic coupling through the paramagnetic (PM) – antiferromagnetic (AFM) – ferrimagnetic (FIM) transitions in multiferroic hexagonal ErMnO 3 have been investigated by Resonant Ultrasound Spectroscopy. Elastic stiffening by up to 2% below the PM – AFM transition at 80 K arises from biquadratic coupling between strain and the magnetic order parameter with relaxation times longer than ~10 -6 s for the response of spins to changes in strain. In contrast with YMnO 3 , the PM – AFM transition in ErMnO 3 is accompanied by a peak in acoustic loss immediately below the Néel point which is interpreted in terms of strain relaxation accompanying ordering of spins of Er 3+ at 4b sites. Changes in the magnetic ordering scheme at the AFM – FIM transition near 3 K are accompanied by elastic softening of ~0.03 %. During poling of the low temperature ferrimagnetic structure round magnetic hysteresis loops, small changes in elastic stiffness which arise due to the contribution of piezomagnetic and/or piezoelectric moduli are detected. Contributions of piezoelectric moduli to acoustic resonance frequencies also permit changes in the configuration of ferroelectric domains to be detected in response both to cycling through this transition and to application of a magnetic field. A peak in acoustic loss in the vicinity of 250 K is attributed to strain-mediated pinning/freezing of some aspect of the domain microstructure with an activation energy of ~0.25–0.3 eV. A return to the original elastic properties on heating to temperatures above ~250 K is interpreted in terms of backswitching of domains to the configuration they had at the start. These observations confirm the existence of subtle variations in magnetoelastic coupling behaviour relating to both the magnetic order parameters and magnetic domain structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spherical Indentation and Implementation of S 3 /P for yield stress determination of brittle materials

Here, a mathematically transparent and robust experimental method has been developed to estimate the yield stress of brittle materials through the analysis of depth-sensing spherical indentation. Employing Hertzian contact mechanics, an elastically invariant ratio based on the simple equation S 3 / P = 6RE r 2 , (where S and P are contact stiffness and indentation load, respectively) has been derived that enables more accurate and confident determination of the transition from elastic to inelastic deformation; a transition that the yield stress dictates and represents. Using two diamond spheres with radii of 3.2 and 8.6 μm, the indentation test method and analyses are applied to two vitreous silicates: Corning's HPFS 7980® fused silica and Vitro's Starphire® soda lime silicate. The estimated yield strengths are 8.15 GPa ± 2.5% for the fused silica and 6.1 GPa ± 3.3% for the soda-lime silicate, and both were independent of indenter radius. Verification of this new experimental method is demonstrated with an as-drawn titanium by showing equivalence of measured yield stress by its spherical indentation and that from uniaxial compression testing. This method will enable easier and more confident estimation of yield stress in brittle materials - a property that historically has been elusive to measure for these materials using common laboratory mechanical test methods.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗
DOE OSTI2026

A physically based mechanical model for Mullins effect in thermoplastic polyurethanes

Despite decades of research, connecting the chemical and physical structure of thermoplastic polyurethanes to their mechanical properties remains highly challenging. Of particular note are their large-deformation and rate-dependent behaviors, which vary greatly with molecular chemistry, including the type and relative content of soft and hard segments. In this work, we develop a physically motivated mechanical theory for predicting the behavior of thermoplastic polyurethanes. The theory incorporates a representation of microstructural evolution during mechanical deformation, which captures the signatures of stress softening over cyclic loading (commonly referred to as the Mullins effect). There are only eight physically motivated fitting parameters, including a direct dependence on the hard segment fraction. The model predicts that increasing the hard segment fraction leads to higher stiffness and greater energy dissipation, in quantitative agreement with published experimental data. Furthermore, we provide a comprehensive analysis of the model and validate its predictions across several independent datasets focused on mechanical characterization. Direct comparisons to experimental data demonstrate its predictive capability on the effect of loading rate, cyclic deformations, and applied tension or compression. Altogether, this work establishes a predictive framework that connects polymer chemistry and microstructure to emergent mechanical behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Method to account for natural fracture induced elastic anisotropy in geomechanical characterization of shale gas reservoirs

Shale has been usually recognized as a transverse isotropic (TI) medium in conventional geomechanical log interpretation due to its laminated nature. However, when natural fractures exist in the shale rock, additional elastic anisotropy is introduced, converting laminated Shale to an orthorhombic (OB) medium. Previous studies illustrate that neglecting the natural fracture induced anisotropy in shale geomechanical log interpretation could lead to inaccurate evaluations of elastic moduli and in-situ stresses. In this paper, a new method is developed to account for the natural fracture induced anisotropy in geomechanical log interpretation based upon the TI acoustic model developed by the author and a characterization technique of elastic wave anisotropy (Sayers, 1991). The new OB model incorporates the four acoustic log data inputs and five modeling constraints in a nonlinear optimization algorithm to solve for the nine independent stiffness coefficients of an OB rock, and further to solve for the geomechanical properties and in-situ stress profiles in an OB formation. The new method was validated with a Marcellus Gas Shale field case. Both the new OB model and the conventional TI model were applied to interpret the minimum horizontal stress profile for the same formation. By comparing the results, the OB model is more robust than the TI from two aspects. First, the average stress magnitude predicted by the OB model is closer to the one measured by the Diagnostic Fracture Injection Test (DFIT). Second, the OB model predicts a more obvious stress barrier between the lower Marcellus and upper Onondaga Limestone than the TI model does. Finally, the predicted stress barrier is consistent with the observation of the microseismic events of a horizontal well drilled and completed nearby, which reveals that no hydraulic fracture propagates downward through the bottom boundary of Marcellus Shale into the underlying Onondaga Limestone.

03 NATURAL GAS↗
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 OSTI2020

Defect evolution in Ni and solid-solution alloys of NiFe and NiFeCoCr under ion irradiation at 16 and 300 K

Single-phase concentrated solid-solution alloys (SP-CSAs) have shown unique chemical complexity at the levels of electrons and atoms, and their defect evolution is expected to be different from conventional dilute alloys. Single crystals of Ni, NiFe and NiFeCoCr are chosen as model systems to understand the chemical complexity on defect formation and damage accumulation in SP-CSAs under ion irradiation. The high-quality crystals were irradiated at 16 and 300 K to different ion fluences, to form irradiated region with little to heavy damages. Additionally, the ion-induced damage was determined using Rutherford backscattering spectrometry technique along a channeling direction (RBS/C) and the level of lattice damage in irradiated Ni and SP-CSAs was quantified from Monte Carlo (MC) simulations. The results are interpreted using the Multi Step Damage Accumulation model to reveal material damage accumulation kinetics. Key findings of the study are that in case of room temperature irradiations the damage level measured for complex alloys at the highest irradiation fluence of 2 × 10 15 cm -2 (~3 dpa) is significantly higher than that obtained for pure nickel samples and suggest two-step damage accumulation process with a defect transformation taking place at a fluence of about 1.5 × 10 15 cm -2 . Furthermore, structural and damage kinetic differences clearly imply that, with increasing degree of chemical complexity and high solid-solution strengthening effects from Ni to NiFe and to NiFeCoCr, the enhanced lattice stiffness resists to randomization of atomic configurations and inhibits the growth of extended defects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Micro-cantilever beam experiments and modeling in porous polycrystalline UO 2

Understanding the impact of microstructure on the thermo-mechanical behavior of oxide nuclear fuels is vital to predicting their performance through multiscale models. Evaluating the mechanical properties at the sub-grain length scale is key to developing these multiscale models. In this work, 3D finite element (FE) models were constructed to simulate the micrometer-scale bending of micro-cantilever beams fabricated using porous polycrystalline uranium dioxide (UO 2 ) and tested at room temperature. Here, the results showed that the porosity and elastic anisotropy of individual grains can play a significant role in determining the effective mechanical properties of the material deduced from the tests. Specifically, the porosity had a non-negligible effect, given that the pore size was of the same order of magnitude as the dimensions of the micro-beams. Correlations between load-deflection data, pore location, and elastic properties (effective Young's modulus) were investigated using UO 2 micro-beam FE models, where pore clusters were included and placed at different locations along the length of the beam. Results indicated that the presence of pore clusters near the substrate, i.e., the clamp of the micro-cantilever beam, has the strongest effect on the load-deflection behavior, with the porosity leading to a reduction of stiffness that is the largest for any location of the pore clusters. Furthermore, it was also found that pore clusters located towards the middle of the span and close to the end of the beam have a comparatively small effect on the load-deflection behavior. Therefore, it is concluded that accurate estimates of Young's modulus can be obtained from micro-cantilever experiments after accounting for porosity on the one third of the beam length close to the clamp. This, in turn, provides an avenue to improve microscale experiments and their analysis in porous, anisotropic elastic materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Anisotropic temperature-dependent elastic constants and thermal conductivities of TRISO particle coatings

Tristructural isotropic (TRISO) particles show great promise as a candidate fuel for use in several next-generation high-temperature nuclear reactor designs due to their structural integrity and fuel performance at high temperatures and burnups. Computational codes exist that can simulate TRISO fuel performance characteristics and failure probabilities under extreme conditions which require knowledge of the TRISO coatings’ thermophysical properties. The thermophysical descriptions of the TRISO particle's layers (i.e., buffer, pyrolytic carbon, and silicon carbide) currently used in fuel performance codes, however, assume that many of these properties are constant with respect to temperature or texture. In this paper, we obtain the full elastic stiffness tensors of the carbon and silicon carbide layers, which have transversely isotropic symmetry. Additionally, the calculated elastic properties of each of the coating layers are in remarkable agreement with the current models at room temperature and correct orientations. Additionally, the textured 3C-SiC layer was found to exhibit novel auxetic behavior above 1500 °C. The anisotropic temperature-dependent thermal conductivities of the carbon layers were calculated using acoustical Grüneisen-Debye theory which are in excellent agreement with current models at room temperature and correct orientations. These texture- and temperature-dependent relationships can be incorporated into the thermophysical description of TRISO particles in order to more accurately model fuel performance and failure probabilities under extreme conditions in forthcoming high-fidelity computational simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fracture mechanics approach to TRISO fuel particle failure analysis

Weibull stress-based methods for failure probability assessment have been developed and analyzed to assess the integrity of tristructural isotropic (TRISO) fuel particles during fuel life cycles and accident operating conditions. While simple, these methods entail a number of drawbacks when stress concentrates near crack tips, including finite element mesh size dependency when the Weibull stress is averaged over the finite element domain. Fracture mechanics approaches involving the use of interaction integrals eliminate this lack of mesh convergence and produce consistent fracture predictions. In this work, we use an interaction integral approach to computing stress intensity factors for a crack in the inner pyrolytic carbon layer perpendicular to the silicon carbide layer, which is simplified representation of a failure mode in TRISO particles. Further, the interface between these two TRISO layers has been shown to become porous, which we simulate by considering a transition of mechanical properties over such porous length, i.e. the layers are modeled as a functionally graded material. Aspects such as porosity and thermal and irradiation eigenstrains are considered in computing the stress intensity factor from a fracture mechanics approach and compared with the known Weibull stress failure approach. The methodology introduced in this paper enables a more general fracture probability assessment in TRISO particles and eliminates the need to identify best suited parameters when using local or averaged stress-based failure criterion. Finally, the numerical sensitivity studies show how parameters such as the porous transition zone length, the material stiffness, and creep affect the probability of TRISO fuel particle failure.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Probing the elastic modulus and hardness of superionic boron cluster solid electrolytes

Polyhedral borane salts are highly tunable compounds with the potential to be used as solid electrolytes due to their high Li and Na superionic conductivity and relatively wide electrochemical stability window. In considering their application to all solid-state batteries, their mechanical properties play a critical role in their operation due to stresses the solid electrolyte is subjected to during battery operation. Density functional theory (DFT) calculations have provided an initial assessment of bulk and elastic moduli for some selected boron cluster solid electrolytes, but direct measurements are still scarce. Here, in this paper, we report the elastic moduli and hardnesses of LiCB 11 H 12 , LiCB 9 H 10 , NaCB 11 H 12 , and NaCB 9 H 10 for the first time using nanoindentation continuous stiffness method (CSM) measurements under inert atmosphere. Experimental modulus values ranging 8.8–12.3 GPa and hardness values ranging 0.15–0.38 GPa are lower than that of other inorganic solid-state electrolytes such as oxide- and sulfide-based solid electrolytes. DFT calculations on the expected moduli of these compounds are also presented and discussed. Analysis of the indentation plasticity index reveals that these compounds have higher plasticity index compared to other common oxide and sulfide solid electrolytes. According to the calculated Pugh's ratio, all compounds in this study except LiCB 11 H 12 are considered ductile.

25 ENERGY STORAGE↗
DOE OSTI2023

Experimental and numerical study of the shielding effect of two tandem rough cylinders in flow-induce oscillation

The shielding effect of the downstream cylinder in flow induced oscillation (FIO) of two cylinders arranged in tandem is studied experimentally and numerically at Reynolds number 30,000 to 120,000. Both cylinders are in one degree-of-freedom, transverse-oscillations, and have turbulence stimulation in the form of selective surface roughness to expand FIO beyond vortex-induced vibration (VIV) into galloping. Shielding of the downstream cylinder has a negative effect on harnessing hydrokinetic energy. Further, to study its effect and mechanics, selective cases are studied both numerically and experimentally and discussed to demonstrate the shielding effect on the downstream cylinder and understand its cause. The main conclusions are: (1) The shielding effect for the downstream cylinder shows a strong relation to the damping ratio. As the damping ratio increases, the shielding effect is mitigated. Additionally, the oscillation of the rear cylinder becomes stable and shows stable frequency. (2) In the VIV region, as the stiffness and natural frequency increase, the shielding effect decreases substantially. (3) In the VIV region, the vorticity of the vortices shedding from both the upper and the lower sides of the downstream cylinder does not accumulate enough due to the attraction by the vortices shed from the upstream cylinder, thus resulting in partial suppression of the oscillation on the downstream cylinder. (4) In the galloping region, the shielding effect for the downstream cylinder depends on whether the vorticity near the downstream cylinder is strengthened by the vortices generated by the shear layers of the upstream cylinder or weakened.

42 ENGINEERING↗
DOE OSTI2022

Effect of laser scanning speed on the microstructure, phase transformation and mechanical property of NiTi alloys fabricated by $\mathrm{LPBF}$

This study investigated the effects of laser scanning speed on the microstructure, phase transformation and properties of NiTi alloys fabricated by laser powder bed fusion (LPBF). In this study, the contributions of metallurgical factors under different scanning speeds, such as Ni evaporation, Ni4Ti3 precipitation, dislocations and internal stress, to the transformation temperature and transformation latent heat were clarified through specially designed experiments. Ni evaporation is found to have the most profound effect, followed by precipitation. Increasing scanning speed is found to reduce Ni loss, thus cause less increase in the transformation temperature and transformation heat of the LPBF-NiTi alloys. Increasing scanning speed also increases the microstructure non-uniformity and thus widens the transformation temperature interval. The orientations of residual stress exhibit strong crystallographic stiffness dependence. The LPBF-NiTi alloys with different scanning speeds all exhibited high strains (>13.4%) and excellent shape memory effect. A LPBF-NiTi honeycomb structure exhibited 96% shape recovery rate after a 60% pre-compressive deformation. Besides, there is an optimum scanning speed for minimum porosity and smallest average pore size. However, the pore structure is found to have weak influence on the tensile behaviour of the LPBF-NiTi alloys, possibly due to the high defect tolerance of the martensitic transformation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Property optimized energy absorber for automotive bumpers utilizing multi-material and structural design strategies

This study proposes a novel design for automotive bumper using optimized lattice structures and multi-materials to balance low-speed collision and high-speed pedestrian impact performance. Different blends of 20 % carbon fiber-reinforced acrylonitrile butadiene styrene with thermoplastic polyurethane were used to tailor material properties. The energy absorber features lattice structures with customized mechanical responses, created by varying the incline angle θ from 0 to 180°. We conducted 576 finite element simulations on a half-scale model to optimize energy absorption and stiffness, leading to 66 optimized designs that met both low-speed and high-speed impact criteria. Two sub-scale optimized energy absorbers with different peak forces—both meeting low-speed impact requirements—were 3D printed and validated through drop-weight testing. The one with lower peak stress demonstrated a more compliant response, exhibiting approximately 90 % lower initial peak force and an increase in energy absorption of around 33 % (from 24 J to 32 J). Compared to the baseline triangular lattice, the optimized absorber increased energy absorption by 68 % from (19 J to 32 J) and reduced peak stress by 70 %. It also showed near-complete recovery with minimal fractures, making it suitable for repeated use. This design improves safety while offering a lightweight, durable, and cost-effective bumper system.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Rethinking production of machine tool bases: Polymer additive manufacturing and concrete

Cast iron and steel weldments are the most common machine tool base elements. However, both construction methods have associated disadvantages for domestic machine tool manufacturers. Here, this paper documents the investigation of an alternative method for machine tool base production using concrete to fill an additively manufactured polymer mold, where the motion components are attached to the concrete base after the initial concrete curing. Modal testing results for a three-axis, vertical spindle prototype indicate high damping and stiffness can be achieved using the concrete base construction. Advantages are reduced cost and lead time compared to traditional methods.

Additive manufacturing↗
DOE OSTI2022

A machining digital twin for hybrid manufacturing

Hybrid manufacturing consisting of metal additively manufactured preforms and computer numerical control (CNC) machining has been established to be an effective method for high material use rates. However, hybrid manufacturing introduces unique challenges. Near-net shape designs are typically selected, which result in a smaller margin for part placement within the stock and stringent requirements for work coordinate system identification. Additionally, less stock material reduces the preform stiffness, which limits the material removal rates during machining. This paper demonstrates a digital twin for CNC machining of a wire arc additively manufactured preform that implements: 1) structured light scanning for stock model identification and tool path generation; 2) a fused filament fabrication apparatus to attach temporary fiducials and scan targets to the preform that enable coordinate system definition for both the CAM and CNC machine; 3) preform and tool tip frequency response function measurements to enable stable milling parameter selection; and 4) post-manufacturing measurements of geometry, surface finish, and structural dynamics to confirm designer intent. These efforts define key components of the machining digital twin for hybrid manufacturing.

36 MATERIALS SCIENCE↗