Search NASA⌕ Search

SEARCH · Search NASA

Results for “strain gradient elasticity”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Higher-order effects of initial deformation on the vibrations of crystal plates

A system of plate equations for the thickness-shear and flexural vibrations superposed on large initial deflection due to bending is derived; in the stress-strain relations, the terms associated with the fourth-order elastic stiffness coefficients are retained. An explicit formula for the change in the fundamental cut-off thickness shear frequency is obtained, and the effects of the terms associated with the fourth-order constants appear to be significant for large gradients of the rotation angles.

Markenscoff, X.↗
DOE OSTI2021

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗
NASA NTRS2016

A Continuum Damage Mechanics Model to Predict Kink-Band Propagation Using Deformation Gradient Tensor Decomposition

A new model is proposed that represents the kinematics of kink-band formation and propagation within the framework of a mesoscale continuum damage mechanics (CDM) model. The model uses the recently proposed deformation gradient decomposition approach to represent a kink band as a displacement jump via a cohesive interface that is embedded in an elastic bulk material. The model is capable of representing the combination of matrix failure in the frame of a misaligned fiber and instability due to shear nonlinearity. In contrast to conventional linear or bilinear strain softening laws used in most mesoscale CDM models for longitudinal compression, the constitutive response of the proposed model includes features predicted by detailed micromechanical models. These features include: 1) the rotational kinematics of the kink band, 2) an instability when the peak load is reached, and 3) a nonzero plateau stress under large strains.

Bergan, Andrew C.↗
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↗
NASA NTRS1999

Determination of Creep Behavior of Thermal Barrier Coatings Under Laser Imposed High Thermal and Stress Gradient Conditions

A laser sintering/creep technique has been established to determine the creep behavior of thermal barrier coatings under steady-state high heat flux conditions. For a plasma sprayed zirconia-8 wt. % yttria coating, a significant primary creep strain and a low apparent creep activation energy were observed. Possible creep mechanisms involved include stress induced mechanical sliding and temperature and stress enhanced cation diffusion through the splat and grain boundaries. The elastic modulus evolution, stress response, and total accumulated creep strain variation across the ceramic coating are simulated using a finite difference approach. The modeled creep response is consistent with experimental observations.

Zhu, Dongming↗
NASA NTRS1987

Thermal stress analysis of space shuttle orbiter wing skin panel and thermal protection system

Preflight thermal stress analysis of the space shuttle orbiter wing skin panel and the thermal protection system (TPS) was performed. The heated skin panel analyzed was rectangular in shape and contained a small square cool region at its center. The wing skin immediately outside the cool region was found to be close to the state of elastic instability in the chordwise direction based on the conservative temperature distribution. The wing skin was found to be quite stable in the spanwise direction. The potential wing skin thermal instability was not severe enough to tear apart the strain isolation pad (SIP) layer. Also, the preflight thermal stress analysis was performed on the TPS tile under the most severe temperature gradient during the simulated reentry heating. The tensile thermal stress induced in the TPS tile was found to be much lower than the tensile strength of the TPS material. The thermal bending of the TPS tile was not severe enough to cause tearing of the SIP layer.

Ko, William L.↗
NASA NTRS1995

Numerical simulation of the nonlinear response of composite plates under combined thermal and acoustic loading

A time-domain study of the random response of a laminated plate subjected to combined acoustic and thermal loads is carried out. The features of this problem also include given uniform static inplane forces. The formulation takes into consideration a possible initial imperfection in the flatness of the plate. High decibel sound pressure levels along with high thermal gradients across thickness drive the plate response into nonlinear regimes. This calls for the analysis to use von Karman large deflection strain-displacement relationships. A finite element model that combines the von Karman strains with the first-order shear deformation plate theory is developed. The development of the analytical model can accommodate an anisotropic composite laminate built up of uniformly thick layers of orthotropic, linearly elastic laminae. The global system of finite element equations is then reduced to a modal system of equations. Numerical simulation using a single-step algorithm in the time-domain is then carried out to solve for the modal coordinates. Nonlinear algebraic equations within each time-step are solved by the Newton-Raphson method. The random gaussian filtered white noise load is generated using Monte Carlo simulation. The acoustic pressure distribution over the plate is capable of accounting for a grazing incidence wavefront. Numerical results are presented to study a variety of cases.

Mei, Chuh↗
DOE OSTI2023

Operando neutron diffraction reveals mechanisms for controlled strain evolution in 3D printing

Abstract Residual stresses affect the performance and reliability of most manufactured goods and are prevalent in casting, welding, and additive manufacturing (AM, 3D printing). Residual stresses are associated with plastic strain gradients accrued due to transient thermal stress. Complex thermal conditions in AM produce similarly complex residual stress patterns. However, measuring real-time effects of processing on stress evolution is not possible with conventional techniques. Here we use operando neutron diffraction to characterize transient phase transformations and lattice strain evolution during AM of a low-temperature transformation steel. Combining diffraction, infrared and simulation data reveals that elastic and plastic strain distributions are controlled by motion of the face-centered cubic and body-centered cubic phase boundary. Our results provide a new pathway to design residual stress states and property distributions within additively manufactured components. These findings will enable control of residual stress distributions for advantages such as improved fatigue life or resistance to stress-corrosion cracking.

36 MATERIALS SCIENCE↗
NASA NTRS1988

MHOST: An efficient finite element program for inelastic analysis of solids and structures

An efficient finite element program for 3-D inelastic analysis of gas turbine hot section components was constructed and validated. A novel mixed iterative solution strategy is derived from the augmented Hu-Washizu variational principle in order to nodally interpolate coordinates, displacements, deformation, strains, stresses and material properties. A series of increasingly sophisticated material models incorporated in MHOST include elasticity, secant plasticity, infinitesimal and finite deformation plasticity, creep and unified viscoplastic constitutive model proposed by Walker. A library of high performance elements is built into this computer program utilizing the concepts of selective reduced integrations and independent strain interpolations. A family of efficient solution algorithms is implemented in MHOST for linear and nonlinear equation solution including the classical Newton-Raphson, modified, quasi and secant Newton methods with optional line search and the conjugate gradient method.

Nakazawa, S.↗
NASA NTRS1995

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗
NASA NTRS2009

Use of Raman Spectroscopy and Delta Volume Growth from Void Collapse to Assess Overwrap Stress Gradients Compromising the Reliability of Large Kevlar/Epoxy COPVs

Composite Overwrapped Pressure Vessels (COPVs) are frequently used for storing pressurized gases aboard spacecraft and aircraft when weight saving is desirable compared to all-metal versions. Failure mechanisms in fibrous COPVs and variability in lifetime can be very different from their metallic counterparts; in the former, catastrophic stress-rupture can occur with virtually no warning, whereas in latter, a leak before burst design philosophy can be implemented. Qualification and certification typically requires only one burst test on a production sample (possibly after several pressure cycles) and the vessel need only meet a design burst strength (the maximum operating pressure divided by a knockdown factor). Typically there is no requirement to assess variability in burst strength or lifetime, much less determine production and materials processing parameters important to control of such variability. Characterizing such variability and its source is crucial to models for calculating required reliability over a given lifetime (e.g. R = 0.9999 for 15 years). In this paper we present a case study of how lack of control of certain process parameters in COPV manufacturing can result in variations among vessels and between production runs that can greatly increase uncertainty and reduce reliability. The vessels considered are 40-inch ( NASA Glenn Research center, Cleveland, OH, 44135 29,500 in3 ) spherical COPVs with a 0.74 in. thick Kevlar49/epoxy overwrap and with a titanium liner of which 34 were originally produced. Two burst tests were eventually performed that unexpectedly differed by almost 5%, and were 10% lower than anticipated from burst tests on 26-inch sister vessels similar in every detail. A major observation from measurements made during proof testing (autofrettage) of the 40-inch vessels was that permanent volume growth from liner yielding varied by a factor of more than two (150 in3 to 360 in3 ), which suggests large differences in the residual stress gradient through their overwraps. This resulted in large uncertainty in true fiber stress ratio (fiber stress at operating pressure divided by fiber stress at burst) which governs lifetime. The vessels were originally designed with tight safety margins, so it became crucial to develop a non-destructive evaluation (NDE) technique to directly measure the overwrap residual stress state of each vessel, and to identify those vessels at highest risk of having poor reliability. This paper describes a Raman Spectroscopy technique for measuring certain patterns of fluctuation in fiber elastic strains over the outside vessel surface (where all but one wrap is exposed at certain locations) that are shown to directly correlate to increased fiber stress ratios and reduced reliability.

Kezirian, Michael T.↗
NASA NTRS2005

Effects of Control Mode and R-Ratio on the Fatigue Behavior of a Metal Matrix Composite

Composite Because of their high specific stiffness and strength at elevated temperatures, continuously reinforced metal matrix composites (MMC's) are under consideration for a future generation of aeropropulsion systems. Since components in aeropropulsion systems experience substantial cyclic thermal and mechanical loads, the fatigue behavior of MMC's is of great interest. Almost without exception, previous investigations of the fatigue behavior of MMC's have been conducted in a tension-tension, load-controlled mode. This has been due to the fact that available material is typically less than 2.5-mm thick and, therefore, unable to withstand high compressive loads without buckling. Since one possible use of MMC's is in aircraft skins, this type of testing mode may be appropriate. However, unlike aircraft skins, most engine components are thick. In addition, the transient thermal gradients experienced in an aircraft engine will impose tension-compression loading on engine components, requiring designers to understand how the MMC will behave under fully reversed loading conditions. The increased thickness of the MMC may also affect the fatigue life. Traditionally, low-cycle fatigue (LCF) tests on MMC's have been performed in load control. For monolithic alloys, low-cycle fatigue tests are more typically performed in strain control. Two reasons justify this choice: (1) the critical volume from which cracks initiate and grow is generally small and elastically constrained by the larger surrounding volume of material, and (2) load-controlled, low-cycle fatigue tests of monolithics invariably lead to unconstrained ratcheting and localized necking--an undesired material response because the failure mechanism is far more severe than, and unrelated to, the fatigue mechanism being studied. It is unknown if this is the proper approach to composite testing. However, there is a lack of strain-controlled data on which to base any decisions. Consequently, this study addresses the isothermal, LCF behavior of a [0]_32 MMC tested under strain- and load-controlled conditions for both zero-tension and tension-compression loading conditions. These tests were run at 427 C on thick specimens of SiC-reinforced Ti-15-3. For the fully-reversed tests, no difference was observed in the lives between the load- and strain-controlled tests. However, for the zero-tension tests, the strain-controlled tests had longer lives by a factor of 3 in comparison to the load-controlled tests. This was due to the fact that under strain-control the specimens cyclically softened, reducing the cracking potential. In contrast, the load-controlled tests ratcheted toward larger tensile strains leading to an eventual overload of the fibers. Fatigue tests revealed that specimens tested under fully-reversed conditions had lives approximately an order of magnitude longer than for those specimens tested under zero tension. When examined on a strain-range basis, the fully reversed specimens had similar, but still shorter lives than those of the unreinforced matrix material. However, the composite had a strain limitation at short lives because of the limited strain capacity of the brittle ceramic fiber. The composite also suffered at very high lives because of the lack of an apparent fatigue limit in comparison to the unreinforced matrix. The value of adding fibers to the matrix is apparent when the fatigue lives are plotted as a function of stress range. Here, the composite is far superior to the unreinforced matrix because of the additional load-carrying capacity of the fibers.

Source record↗
DOE OSTI2025

Prediction of carbon nanostructure mechanical properties and the role of defects using machine learning

Graphene-based nanostructures hold immense potential as strong and lightweight materials, however, their mechanical properties such as modulus and strength are difficult to fully exploit due to challenges in atomic-scale engineering. This study presents a database of over 2,000 pristine and defective nanoscale CNT bundles and other graphitic assemblies, inspired by microscopy, with associated stress–strain curves from reactive molecular dynamics (MD) simulations using the reactive INTERFACE force field (IFF-R). These 3D structures, containing up to 80,000 atoms, enable detailed analyses of structure-stiffness-failure relationships. By leveraging the database and physics- and chemistry-informed machine learning (ML), accurate predictions of elastic moduli and tensile strength are demonstrated at speeds 1,000 to 10,000 times faster than efficient MD simulations. Hierarchical Graph Neural Networks with Spatial Information (HS-GNNs) are introduced, which integrate chemistry knowledge. HS-GNNs as well as extreme gradient boosted trees (XGBoost) achieve forecasts of mechanical properties of arbitrary carbon nanostructures with only 3 to 6% mean relative error. The reliability equals experimental accuracy and is up to 20 times higher than other ML methods. Predictions maintain 8 to 18% accuracy for large CNT bundles, CNT junctions, and carbon fiber cross-sections outside the training distribution. The physics- and chemistry-informed HS-GNN works remarkably well for data outside the training range while XGBoost works well with limited training data inside the training range. The carbon nanostructure database is designed for integration with multimodal experimental and simulation data, scalable beyond 100 nm size, and extendable to chemically similar compounds and broader property ranges. The ML approaches have potential for applications in structural materials, nanoelectronics, and carbon-based catalysts.

Winetrout, Jordan J.↗
DOE OSTI2022

Mechanical Properties of Functionally Graded Concrete Lining for Deep Underground Structures

With the mining depth of coal resources increasing, the thickness of traditional lining for deep mines becomes large. The bearing capacity of the outer lining cannot be fully utilized, so a new type of functionally graded lining structure (FGL) with radial Young’s modulus varying in gradient is proposed. In this study, through theoretical analysis and numerical simulation, the mechanical properties of the functionally graded lining were studied, including the characteristics of the elastic working state and the ultimate bearing state. The influences of the structural parameters of the functionally graded lining, including the inner radius, the thickness, Young’s modulus, the compressive strength of the concrete, Poisson’s ratio, and the number of layers on the mechanical properties were analyzed. The calculation formula of the ultimate bearing capacity of the lining and the calculation formula of the maximum tangential strain at the time of lining failure were put forward. The accuracy of the formula was verified by comparing with the numerical simulation results. The research results provide the basis for the construction of the design theory of a functionally graded lining structure and have great engineering significance for the construction of the kilometer-deep mines in the future.

Li, De-Chun↗
NASA NTRS2006

Porous Media Approach for Modeling Closed Cell Foam

In order to minimize boil off of the liquid oxygen and liquid hydrogen and to prevent the formation of ice on its exterior surface, the Space Shuttle External Tank (ET) is insulated using various low-density, closed-cell polymeric foams. Improved analysis methods for these foam materials are needed to predict the foam structural response and to help identify the foam fracture behavior in order to help minimize foam shedding occurrences. This presentation describes a continuum based approach to modeling the foam thermo-mechanical behavior that accounts for the cellular nature of the material and explicitly addresses the effect of the internal cell gas pressure. A porous media approach is implemented in a finite element frame work to model the mechanical behavior of the closed cell foam. The ABAQUS general purpose finite element program is used to simulate the continuum behavior of the foam. The soil mechanics element is implemented to account for the cell internal pressure and its effect on the stress and strain fields. The pressure variation inside the closed cells is calculated using the ideal gas laws. The soil mechanics element is compatible with an orthotropic materials model to capture the different behavior between the rise and in-plane directions of the foam. The porous media approach is applied to model the foam thermal strain and calculate the foam effective coefficient of thermal expansion. The calculated foam coefficients of thermal expansion were able to simulate the measured thermal strain during heat up from cryogenic temperature to room temperature in vacuum. The porous media approach was applied to an insulated substrate with one inch foam and compared to a simple elastic solution without pore pressure. The porous media approach is also applied to model the foam mechanical behavior during subscale laboratory experiments. In this test, a foam layer sprayed on a metal substrate is subjected to a temperature variation while the metal substrate is stretched to simulate the structural response of the tank during operation. The thermal expansion mismatch between the foam and the metal substrate and the thermal gradient in the foam layer causes high tensile stresses near the metal/foam interface that can lead to delamination.

Ghosn, Louis J.↗
DOE OSTI2025

Nanoscale Polarization-Dependent Young’s Modulus of Ferroelectric BaTiO 3 Single Crystals

Flexoelectric fields produced by strong strain gradients at the nanoscale couple to ferroelectric polarization, promoting changes in the mechanical properties of ferroelectric materials as a function of the direction of the ferroelectric polarization. In this work, we calculate the asymmetry in the Young’s modulus found in oppositely polarized out-of-plane domains of BaTiO 3 single crystals by means of contact resonance AFM, and we evaluate its impact on the electromechanical response as measured by piezoresponse force microscopy, both using band excitation modes. Furthermore, we analyze the electromechanical response of the different BaTiO 3 domains using k-means to build up their mean elastic and electromechanical features and quantify the flexoelectrically induced modulation of the Young’s modulus. Finally, we discuss the use of contact Kelvin probe force microscope measurements to decouple electrostatic artifacts from purely piezoelectric ones when flexoelectric coupling appears.

BaTiO 3 single crystals↗
DOE OSTI2020

Needle-Like Ferroelastic Domains in Individual Ferroelectric Nanoparticles

Superior structural, physical and electronic properties make ferroelectric nanocrystals essential in enabling a range of next-generation devices. Ferroelectric responses are determined by crystal structure and domain morphology. The ability to reversibly displace, create and annihilate elastic domains is critical to device applications. Whereas electric-field control has been demonstrated for ferroelectric 180° surface domain walls and vortices, similar control of ferroelastic domains and domain boundaries within individual nanocrystals remains challenging. Using controlled external compressive and tensile axial stress, deterministic and reversible control of highly mobile ferroelastic domains and axial polarization in three dimensions is demonstrated in this study. While many studies exist on ferroelastic domains in thin films and bulk, little is known about ferroelastic interactions at the single nanoparticle level, especially involving domain boundaries. Through combining Bragg coherent X-ray diffractive imaging and Landau theory, strain gradients in individual BaTiO3 nanocrystals are shown to stabilize needle-like ferroelastic twin domains. These domains are highly labile under applied axial stress, producing a locally enhanced electric polarization mediated by a ferroelectric phase transition. The efficacy of Bragg coherent X-ray diffractive imaging in studying in operando domains in three-dimensions is demonstrated, while synergy with theory provides a paradigm for domain boundary engineering and potential for nanoscale functional devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Emergent chirality in a polar meron to skyrmion phase transition

Polar skyrmions are predicted to emerge from the interplay of elastic, electrostatic and gradient energies, in contrast to the key role of the anti-symmetric Dzyalozhinskii-Moriya interaction in magnetic skyrmions. Here, we explore the reversible transition from a skyrmion state (topological charge of –1) to a two-dimensional, tetratic lattice of merons (with topological charge of –1/2) upon varying the temperature and elastic boundary conditions in [(PbTiO 3 ) 16 /(SrTiO 3 ) 16 ] 8 membranes. This topological phase transition is accompanied by a change in chirality, from zero-net chirality (in meronic phase) to net-handedness (in skyrmionic phase). We show how scanning electron diffraction provides a robust measure of the local polarization simultaneously with the strain state at sub-nm resolution, while also directly mapping the chirality of each skyrmion. Using this, we demonstrate strain as a crucial order parameter to drive isotropic-to-anisotropic structural transitions of chiral polar skyrmions to non-chiral merons, validated with X-ray reciprocal space mapping and phase-field simulations.

42 ENGINEERING↗