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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A predictive modeling tool for damage analysis and design of hydrogen storage composite pressure vessels

In this study, a predictive modeling tool is developed for damage analysis and design of hydrogen (H 2 ) storage composite pressure vessels. It integrates micromechanics of matrix cracking into a continuum damage mechanics (CDM) description for damage evolution, and three-dimensional (3D) finite element (FE) modeling of the vessel structural response. At the scale of the composite layer (mesoscale), the temperature-dependent stiffness reduction law in terms of the damage variable for transverse matrix cracking is computed using an Eshelby-Mori-Tanaka approach for the initial composite thermoelastic properties and a self-consistent model for the stiffness reduction as a function of the damage variable. While transverse matrix cracking obeying a damage evolution relation can progressively evolve from an initiation to a saturation state, fiber failure is predicted by a micromechanical fiber rupture criterion that accounts for the fiber strength and matrix stress. The implementation of this integrated multiscale modeling model into a 3D FE formulation enables damage analysis and design of H 2 storage composite pressure vessels. The developed tool is illustrated through 3D damage analyses of a cryogenically compressed H 2 storage vessel model subjected to thermomechanical loadings to investigate effects of the helical layer fiber orientation and loading scenario on damage development, vessel integrity and burst pressure.

08 HYDROGEN↗
DOE OSTI2021

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Relationship between the contact force strength and numerical inaccuracies in piecewise-smooth systems

This work studies the different types of behavior and inaccuracies that can occur when contact is not adequately accounted for in a dynamical system with freeplay, as the strength of the contact stiffness increases. The MATLAB® ode45 time integration solver, with the built-in Event Location capability, is first validated using past experimental data from a forced Duffing oscillator with freeplay. Next, numerical results utilizing event location are compared to results neglecting event location in order to highlight possible numerical errors and effects on multistable dynamical responses. Inaccuracies tend to occur in two different ways. First, neglecting event location can affect the boundaries between basins of attraction. Second, neglecting event location has little effect on the behaviors of the attractor solutions themselves besides merely resembling poorly converged solutions. Errors are less pronounced at the limits of soft or hard contact stiffness. Finally, this study shows the importance of accurately solving piecewise-smooth systems and the existing correlation between the strength of the contact force and possible numerical inaccuracies.

42 ENGINEERING↗
DOE OSTI2022

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Dislocation loop evolution and radiation hardening in nickel-based concentrated solid solution alloys

Effects of chemical composition, ion irradiation dose and temperature on unfaulting of irradiation induced Frank dislocation loops to perfect loops in two nickel based single-phase solid solution alloys, Ni–20Fe and NiFe–20Cr, have been studied. The fraction of Frank loops decreases with irradiation dose from 7.2 to 38.4 dpa at 500°C, but with more Frank loops remaining in the ternary alloy. However, perfect loops and dislocation networks become the dominant features of defects at 580°C in both alloys. The results indicate a thermally assisted loop unfaulting process that may be hindered by more sluggish defect motion in the alloy with more chemical components. Nano-indentation with both continuous stiffness method and single indentation method are used to measure radiation hardening. Loop unfaulting in both alloys irradiated at 580°C reduced radiation hardening while significant hardening is observed after irradiation at 500°C. The quasi-static single indentation method exhibits lower hardness results compared to continuous stiffness method, because dislocations induced from the cyclic loading in the latter method get relaxed and stabilized, resulting in higher resistance to the indenter.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tailoring ductile-phase toughened tungsten hierarchical microstructures for plasma-facing materials

Biological materials, such as bones, nacre, etc., have been found to show attractive and unique combinations of stiffness, strength and fracture toughness. Research has been conducted in various engineering areas to mimic the naturally hierarchical microstructures or nanostructures of these materials to produce materials with optimum stiffness, high strength and fracture toughness for their intended structural applications. In the same objective, this work applies a multiscale microstructural approach recently developed to investigate the deformation and fracture behavior of ductile phase toughened tungsten (W) materials such as tungsten-nickel/iron (W-Ni/Fe) composites that possess lamellar-like and hierarchical “brick-and-mortar” (BAM) microstructures. First, the approach is used to simulate tensile loading of W-Ni/Fe specimens cut out from hot-rolled W-Ni/Fe plates which possess lamellar-like microstructures. The finite element model of the gage section of the specimen consists of a W-Ni/Fe dual-phase microstructural domain in which the constitutive behaviors of W and Ni/Fe phases are described by an elastic-plastic damage model. The predicted material stress-strain response and crack pattern development as a function of loading are then compared to the corresponding experimental results to determine the constitutive model parameters for W and Ni/Fe. Subsequently, the model parameters are used to analyze W-Ni/Fe BAM microstructures that are artificially created to investigate the effects of microstructural features and morphologies on the composite stress-strain response, damage, and fracture patterns. The modeling shows that regular BAM microstructures that experience bridging mechanism combined with crack penetration across W-phase regions exhibit significantly higher strengths than the rather random lamellar-like microstructures. More significantly, however adjusting the brick’s length-to-height ratio, the BAM microstructure can be designed to allow a more distributed damage zone which leads to increased strength, ductility and fracture energy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Thermo-Elastic Properties and Damping of U-6wt%Nb

While thermal expansion data exists for quenched (as well as aged) U-6wt%Nb, there is wide variation in the reported room temperature elastic moduli. To better understand the room temperature elastic behavior and to address the complete absence of data on the temperature dependence of the elastic moduli, room temperature and in-situ elevated temperature resonant ultrasonic spectroscopy (RUS) was performed along with x-ray diffraction (XRD) and dilatometry using a thermomechanical analyzer (TMA). An in-situ small- and wide-angle x-ray scattering (SAXS/WAXS) experiment was performed on a companion sample to help interpret the results. The room temperature polycrystalline dynamic moduli were measured to be E = 95.4 ± 4.5 GPa and G = 35.2 ± 0.2 GPa, respectively. As the temperature is raised, the stiffness slowly decreases, consistent with the empirical rule proposed by Varshni, up to the point of precipitation of the equilibrium α-U phase, which is stiffer than the martensitic phases. Homogenization theory of composites can be used to rationalize the observed response at high temperatures and after high-temperature exposures. Further, aging of the material at low temperatures (≤200°C) does not affect the linear elastic stiffness, but does impact the damping behavior which can be measured through RUS. This change in damping provides another perspective on the microstructure changes induced by low-temperature aging which also result in significant strengthening.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Formulation and calibration of two-dimensional constitutive models for composite structures based on panel tests

Steel Plate Concrete (SC) composite members have been widely adopted because of its cost-efficiency and enhanced structural behavior. While researchers have attempted to study its in-plane shear behavior in the past twenty years, very limited number of large-scale pure shear tests were performed due to the challenge of experimental set-up and the availability of facilities. In this paper, a series of uniaxial loading tests and two full-scale pure shear panel tests of SC members were reported, on which the “mechanics-based Membrane Model of SC elements (MM-SC)” is developed. The MM-SC model is based on the fixed-angle crack formulation and the smeared-crack formulation, in which the experimental-based uniaxial constitutive models are implemented, considering the local buckling of faceplate, the tension stiffening of steel plate, the strength degradation of cracked concrete and the confinement effect of concrete. The proposed MM-SC model is subsequently incorporated into the object-oriented software OpenSEES. Finally, the simulation results of proposed model well predict the SC test observations in terms of critical branch points and structural behaviors, including initial stiffness, cracking strength, post-crack stiffness, yield strength, maximum strength, and failure modes.

42 ENGINEERING↗
DOE OSTI2021

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗
DOE OSTI2023

One-pot ternary sequential reactions for photopatterned gradient multimaterials

Seamless multimaterial construction is a common motif in animal physiology. Such continuous mechanical gradients remain challenging to reproduce in engineered systems, as current resin chemistries typically result in a single fixed set of properties. As an alternative to single-property materials, we introduce a thiol-ene-epoxy-based photothermal reaction scheme that produces multimaterials by altering the polymer microstructure within a single resin. In this system, the photodosage during the first stage of processing dictates the extent of conversion for each subsequent reaction. As a result, our photochemistry can exhibit a diverse range of soft (Young’s modulus, E ~ 400 kPa; elongation, dL/L0 ~ 300%) and stiff (E ~ 1.6 GPa; dL/L 0 ~ 3%) mechanical properties. Furthermore, we pattern photostable and mechanically robust modulus gradients (d[E r, stiff /E r , soft] /dx > 1,000 mm -1 ) that exceed those found in squid beaks and human knee entheses. We demonstrate the ability to build intricate multimaterial architectures including a soft, wearable braille display.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A constitutive model for elastomers tailored by ionic bonds and entanglements

Over the past decade or two, the concept has emerged of using multiple types of weak interactions simultaneously to enhance the mechanical properties of elastomers. These weak interactions include physical entanglements, hydrogen bonds, metal-coordination bonds, dynamic covalent bonds, and ionic bonds. The combination of entanglements and ionic bonding has been minimally explored and is particularly exciting because of the broad application space for polyelectrolytes. In this work, a constitutive model framework is developed to describe the response of elastomers with both ionic bonds and entanglements. Here we formulate a micromechanical model that couples together chain stretching, ionic bond slipping, and entanglement evolution. The ionic bonds provide toughness by enabling plastic deformation in comparison to covalently crosslinked material and add strength compared to a linear polymer. Evolution of the entanglement density is taken as a key mechanism that can govern stiffness, toughness, and self-recovery in elastomers. The model is used to match bulk polyelectrolytes with different fractions of ionic components under a variety of loading histories. The variations in material parameters are then used to help understand the relative importance of different governing mechanisms in the bulk polymers. We show that the theoretical framework can explain our experimental uniaxial tensile experimental results for polyelectrolytes. This model can help to design better material with high stiffness and toughness. We expect that our model can be extended to explain the mechanical behavior of other polyelectrolytes and other soft materials with a wide range of dynamic bonds.

42 ENGINEERING↗
DOE OSTI2025

Elucidating texture and grain morphology contributions to the micromechanical response of additively manufactured Inconel 625

Microstructural variation of additively manufactured (AM) metal components in comparison to wrought counterparts makes certification for critical applications a challenge. Microscale simulations leveraging modern computational tools may be used to supplement testing of AM microstructures, thus accelerating certification by reducing the number of experiments needed. However, as micromechanical response is closely tied to critical properties like fatigue-life and fracture, utilization of these simulations with macroscale experimental data alone is insufficient. One means to attain microscale experimental data is in situ diffraction data collected from synchrotron X-ray sources. In this work, such data were collected during in situ compression of AM Inconel 625 superalloy. Interpretation of experimental results was assisted by massive (8M element) complementary micromechanical simulations performed on sets of virtual microstructures generated using cellular automata. Together, micromechanical data from diffraction experiments and simulations were used to probe the effects of textured “track” microstructures generated during laser powder bed fusion and directional strength-to-stiffness on micromechanical response. Though fiber-averaged directional strength-to-stiffness ratios were expected to dominate given the high elastic anisotropy of the material, the combination of small variations in texture and specific grain configurations unique to AM microstructures lead to significant variability in micromechanical response after yield. The findings emphasize the importance of high-fidelity microstructural representation that captures key texture components and AM-specific morphology for property prediction of AM metals.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effect of fiber sizing and glass fiber laminate hybridization on vibration damping and mechanical properties of banana fiber reinforced polypropylene composites

Modern automotive applications demand lightweight, multifunctional materials to reach mileage goals and natural fiber reinforced composites (NFRCs) are one of the classes of materials proposed as a solution. NFRCs exhibit good vibration damping properties and have low density, but are often limited by processing challenges, poor-fiber matrix compatibility and variable performance. Herein, we investigate non-woven wet-lay of comingled banana fiber (BF), recycled glass fiber (rGF), and polypropylene (PP) fibers to in situ sizing and preparation of composite feedstocks for compression molding. BF and rGF hybrids were prepared by stacking rGF layers during compression molding to produce composites with various fiber ratios. The effect of fiber content, in-situ sizing and ratio of BF to rGF on tensile, flexural and vibration damping performance are investigated. Key results are the significant increase in tensile strength by in situ sizing (40 % sized at 60 wt% BF) and in flexural modulus (+58 % sized at 60 wt% BF) and flexural strength (+41 % sized 60 wt% BF) compared to the unsized equivalent. For BF-rGF hybrid composites with40 wt% total fiber content, flexural strength and modulus were improved by 51 % and 231 % respectively for a 1:1 ratio BF:rGF compared to BF reinforced system. Lastly, identifying the cross-over point where damping and stiffness are optimized for a hybrid composite. These findings demonstrate that these composites can be used as alternative to synthetic fiber or mineral filled composites in automotive applications, particularly where weight reduction, vibration damping and stiffness are desired.

Banana fiber↗
DOE OSTI2025

Wood template-supported phase change material composites for durable and form-stable thermal energy storage in buildings

Here, to reduce and shift peak energy loads in buildings, phase change materials (PCMs) with high transition enthalpies and transition temperatures near human thermal comfort are desirable for thermal energy storage (TES). Traditional solid/liquid PCMs suffer from leakage during thermal cycling, requiring encapsulation that lowers heat storage capacity. Wood templates (WTs), with porous and hierarchical structures, provide natural encapsulation scaffolds for PCM containment. Polyethylene glycol (PEG) is a compatible PCM with WTs, but when infiltrated alone, ~30 wt% of PEG binds to wood with no detectable phase change, limiting TES efficiency. To address this, we developed a method to (1) improve the form stability of balsa and pine-based composites (BWT + PCMs and PWT + PCMs) and (2) reduce inactive PCM within the composites to ~10 wt%. BWT + PCMs exhibit transition properties upwards of 114.2 J/g at 25.4 degrees C, with no degradation after 1000 thermal cycles, and similar stiffness compared to raw balsa. Meanwhile, PWT + PCMs exhibit 25 % higher storage efficiency with the addition of poly(ethylene glycol) diacrylate compared to solely PEG-infiltrated PWT. PWT + PCM retains 84 % of raw pine's mechanical stiffness, sufficient for light-duty construction. Our shape-stabilized WT + PCM composites enrich the functionality of wood materials as both ideal TES material candidates and light-duty building construction applications.

25 ENERGY STORAGE↗
DOE OSTI2025

Towards an intuitive application of WEC control co-design

A simple co-design example in a reduced parameter space is presented for an oscillating flap device. Initially, the WEC geometry and mass properties are considered along with drivetrain gear ratio, inertia, motor constant and stiffness under both PI and optimal control. This parameter space is reduced to those to which performance is most sensitive for a fixed geometry. The gear ratio, drivetrain stiffness, and flap mass are found to be the most impactful design criteria as they can create orders of magnitude variations in power performance. The performance of the optimized system is compared with several sub-optimal variants in terms of electrical and mechanical power capture, transmission coefficients, and transducer power gain. Notably, though substantial power capture improvements are demonstrated when an optimal controller is employed, this power capture remains sensitive to appropriate selections of drivetrain and flap design parameters, implying that control co-design procedures remain necessary for high-performing WECs. Furthermore, a number of practical caveats and extensions to the presented co-design methodology are suggested, including the characterization of system static friction, especially in the presence of high gear ratios.

Control co-design↗
DOE OSTI2023

Structural deformation and mechanical response of CrS 2 , CrSe 2 and Janus CrSSe

In the framework of density functional theory (DFT), we investigate the structural deformation, and mechanical behavior of the Janus CrSSe, which has out-of-plane structural asymmetry, with conventional transition metal dichalcogenides (TMDs) CrS 2 and CrSe 2 . The Janus CrSSe could be a potential candidate for machinable optoelectronic and piezoelectric applications. Here, we predict that these compounds are chemically, mechanically, and dynamically stable with the covalent bond between the TM(Cr) and chalcogen(X=S, Se) atoms. Due to the influence of tensile strain, the Cr–X bond length of each monolayers increases and the thickness decreases. Interestingly, the in-plane stiffness, shear and layer moduli, Poisson’s ratio, ultimate bi/uni-axial stress of Janus CrSSe are in between the values of CrS 2 and CrSe 2 monolayers. Similar to TMDs, the orientation-dependent in-plane stiffness and Poisson’s ratio demonstrates the isotropic behavior in Janus CrSSe. Furthermore, it can sustain a larger value of uni/bi-axial tensile strain with the critical strain equivalent to CrX 2 monolayers. By applying higher-order strain, we have also found average elastic–plastic behavior as expected. These findings demonstrate that the Janus CrSSe monolayer is a mechanically stable and ductile compound that maintains the hybrid behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A gravo-aeroelastically scaled wind turbine rotor at field-prototype scale with strict structural requirements

A new sub-scale field-prototype design solution is developed to realize the dynamics, structural response, and distributed loads (gravitational, aerodynamic, centrifugal) that are characteristic of a full-scale large, modern wind turbine rotor. Prior work in sub-scale wind turbine testing has focused on matching aerodynamic/aero-elastic characteristics of full-scale rotors at wind tunnel scale. However, large-scale rotor designs must expand beyond this limited set of scaling parameters for cost-effective prototyping and meet strict requirements for structural safety for field testing. Here, the challenge lies in producing a structural design meeting two competing objectives: novel scaling objectives that prescribe the sub-scale blade to have low mass and stiffness; and traditional structural safety objectives that drive the design to have higher stiffness and mass. A 20% gravo-aeroelastically scaled wind turbine blade is developed successfully that satisfies these competing objectives. First, it achieved close agreement for non-dimensional tip deflection and flap-wise blade frequency (both within 2.1%) with a blade mass distribution constrained to produce target gravitational and centrifugal loads. Second, the entire blade structure was optimized to ensure a safe, manufacturable solution meeting strict strength requirements for a testing site that can experience up to 45 m/s wind gusts. The prototype-scale blade was fabricated and successfully proof-load tested.

17 WIND ENERGY↗