Search NASA⌕ Search

SEARCH · Search NASA

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗
DOE OSTI2020

Shock compression of molybdenum single crystals to 110 GPa: Elastic–plastic deformation and crystal anisotropy

To explore the role of crystal anisotropy on the elastic-plastic deformation of BCC single crystals at high shock stresses, molybdenum (Mo) single crystals were shock compressed along [100], [111], and [110] orientations at elastic impact stresses between 20 and 110 GPa. Laser interferometry was used to measure shock wave velocities and particle velocity histories. Along the [100] and [111] orientations, elastic-plastic waves (two wave structure) were observed up to 110 GPa. Along the [110] orientation, the two wave structure was observed only up to 90 GPa. The measured elastic wave amplitudes were analyzed to determine crystal anisotropy effects, impact stress dependence, and the activated slip systems on the Hugoniot elastic limit. The results from our work have provided insight into the role of crystal anisotropy on the elastic-plastic deformation under shock compression at high stresses.

42 ENGINEERING↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
NASA NTRS1984

Finite elastic-plastic deformation of polycrystalline metals

Applying Hill's self-consistent method to finite elastic-plastic deformations, the overall moduli of polycrystalline solids are estimated. The model predicts a Bauschinger effect, hardening, and formation of vertex or corner on the yield surface for both microscopically non-hardening and hardening crystals. The changes in the instantaneous moduli with deformation are examined, and their asymptotic behavior, especially in relation to possible localization of deformations, is discussed. An interesting conclusion is that small second-order quantities, such as shape changes of grains and residual stresses (measured relative to the crystal elastic moduli), have a first-order effect on the overall response, as they lead to a loss of the overall stability by localized deformation. The predicted incipience of localization for a uniaxial deformation in two dimensions depends on the initial yield strain, but the orientation of localization is slightly less than 45 deg with respect to the tensile direction, although the numerical instability makes it very difficult to estimate this direction accurately.

Iwakuma, T.↗
NASA NTRS1971

Analysis of finite deformations of elastic solids by the finite element method.

Finite element applications, particularly to analyses of finite deformations in elastic solids, are reviewed, along with the difficulties encountered in the formulation of certain problems and in their numerical solution. Various approaches are discussed for overcoming these and other difficulties. A computer program designed for finite elasticity problems is described, and several numerical examples are presented.

Oden, J. T.↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonlinear deformation and elasticity of BCC refractory metals and alloys

Application of isotropic pressure or uniaxial strain alters the elastic properties of materials; sufficiently large strains can drive structural transformations. Linear elasticity describes stability against infinitesimal strains, while nonlinear elasticity describes the response to finite deformations. Here, it was previously shown that uniaxial strain along [100] drives refractory metals and alloys towards mechanical instabilities. These include an extensional instability, and a symmetry-breaking orthorhombic distortion caused by a Jahn-Teller-Peierls instability that splays the cubic lattice vectors. Here we analyze these transitions in depth. Eigenvalues and eigenvectors of the Wallace tensor identify and classify linear instabilities in the presence of strain. We show that both instabilities are discontinuous, leading to discrete jumps in the lattice parameters. We provide physical intuition for the instabilities by analyzing the changes in first-principles energy, stress, bond lengths, and angles upon application of strain. Electronic band structure calculations show differential occupation of bonding and antibonding orbitals, driven by the changing bond lengths and leading to the structural transformations. Strain thresholds for these instabilities depend on the valence electron count.

36 MATERIALS SCIENCE↗
NASA NTRS1975

Finite-element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is ideally suited to isotropically hardening Prandtl-Reuss materials. Further, the formulation is given in a manner which allows any conventional finite element program, for 'small strain' elastic-plastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. The paper closes with a unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures. Further, a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain, and the inadequacies of some of these are commented upon.

Mcmeeking, R. M.↗
NASA NTRS1974

Finite element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is suited to isotropically hardening Prandtl-Reuss materials. The formulation is given in a manner which allows any conventional finite element program, for "small strain" elasticplastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. A unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures, and a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain.

Mcmeeking, R. M.↗
NASA NTRS2003

Integrated Approach to the Dynamics and Control of Maneuvering Flexible Aircraft

This work uses a fundamental approach to the problem of simulating the flight of flexible aircraft. To this end, it integrates into a single formulation the pertinent disciplines, namely, analytical dynamics, structural dynamics, aerodynamics, and controls. It considers both the rigid body motions of the aircraft, three translations (forward motion, sideslip and plunge) and three rotations (roll, pitch and yaw), and the elastic deformations of every point of the aircraft, as well as the aerodynamic, propulsion, gravity and control forces. The equations of motion are expressed in a form ideally suited for computer processing. A perturbation approach yields a flight dynamics problem for the motions of a quasi-rigid aircraft and an 'extended aeroelasticity' problem for the elastic deformations and perturbations in the rigid body motions, with the solution of the first problem entering as an input into the second problem. The control forces for the flight dynamics problem are obtained by an 'inverse' process and the feedback controls for the extended aeroservoelasticity problem are determined by the LQG theory. A numerical example presents time simulations of rigid body perturbations and elastic deformations about 1) a steady level flight and 2) a level steady turn maneuver.

Waszak, Martin R.↗
NASA NTRS1993

Stress stiffening and approximate equations in flexible multibody dynamics

A useful model for open chains of flexible bodies undergoing large rigid body motions, but small elastic deformations, is one in which the equations of motion are linearized in the small elastic deformations and deformation rates. For slow rigid body motions, the correctly linearized, or consistent, set of equations can be compared to prematurely linearized, or inconsistent, equations and to 'oversimplified,' or ruthless, equations through the use of open loop dynamic simulations. It has been shown that the inconsistent model should never be used, while the ruthless model should be used whenever possible. The consistent and inconsistent models differ by stress stiffening terms. These are due to zeroth-order stresses effecting virtual work via nonlinear strain-displacement terms. In this paper we examine in detail the nature of these stress stiffening terms and conclude that they are significant only when the associated zeroth-order stresses approach 'buckling' stresses. Finally it is emphasized that when the stress stiffening terms are negligible the ruthlessly linearized equations should be used.

Padilla, Carlos E.↗
DOE OSTI2020

Exploiting ultra-large linear elasticity over a wide temperature range in nanocrystalline NiTi alloy

Many shape memory alloys can support large recoverable strains of a few percent by reversible stress-induced martensite transformation, yet they behave non-linear within a narrow operating temperature range. Developing the bulk metallic materials with ultra-large linear elasticity over a wide temperature range has proven to be difficult. In this work, a material design concept was proposed, that is true elastic deformation and reversible twinning-detwinning deformation run in parallel to overcome this challenge. By engineering the residual internal stress to realize the concurrency of true elastic deformation and twinning-detwinning deformation, a bulk nanocrystalline NiTi that possesses an ultra-large linear elastic strain up to 5.1 % and a high yield stress of 2.16 GPa over a wide temperature range of 270 °C was developed. This study offers a new avenue for developing the metallic materials with ultra-large linear elasticity over a wide temperature range of 270 °C (from 70 °C to -197 °C).

36 MATERIALS SCIENCE↗
DOE OSTI2020

In Situ Synchrotron X-ray Diffraction Investigations of the Nonlinear Deformation Behavior of a Low Modulus β -Type Ti36Nb5Zr Alloy

The low modulus β-type Ti alloys usually have peculiar deformation behaviors due to their low phase stability. However, the study of the underlying mechanisms is challenging since some physical mechanisms are fully reversible after the release of the load. In this paper, the deformation behavior of a low modulus β-type Ti36Nb5Zr alloy was investigated with the aid of in situ synchrotron X-ray diffraction (SXRD) during tensile loading. The evolution of lattice strains and relative integrated diffraction peak intensities of both the β and α” phases were analyzed to determine the characteristics of the potential deformation mechanisms. Upon loading, the α” diffraction spots appeared at specific azimuth angles of the two-dimensional SXRD patterns due to the <110> fiber texture of original β grains and the selection of favorable martensitic variants. The nonlinear deformation behavior originated from a reversible stress-induced martensitic transformation (SIMT). However, the SIMT contributed a little to the large recoverable strain of over 2.0%, which was dominated by the elastic deformation of the β phase. Various deformation mechanisms were activated successively at different applied strains, including elastic deformation, SIMT and plastic deformation. Our investigations provide in-depth understandings of the deformation mechanisms in β-type Ti alloys with low elastic modulus.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Identification and control of structures in space

The derivation of the equations of motion for the Spacecraft Control Laboratory Experiment (SCOLE) is reported and the equations of motion of a similar structure orbiting the earth are also derived. The structure is assumed to undergo large rigid-body maneuvers and small elastic deformations. A perturbation approach is proposed whereby the quantities defining the rigid-body maneuver are assumed to be relatively large, with the elastic deformations and deviations from the rigid-body maneuver being relatively small. The perturbation equations have the form of linear equations with time-dependent coefficients. An active control technique can then be formulated to permit maneuvering of the spacecraft and simultaneously suppressing the elastic vibration.

Meirovitch, L.↗
NASA NTRS1985

Maneuvering of flexible spacecraft with application to SCOLE

This paper is concerned with the derivation of the equations of motion for the Spacecraft Control Laboratory Experiment (SCOLE). For future reference, the equations of motion of a similar structure orbiting the earth are also derived. The structure is assumed to undergo large rigid-body maneuvers and small elastic deformations. A perturbation approach is presented where the quantities defining the rigid-body maneuver are assumed to be relatively large, with the elastic deformations and deviations from the rigid-body maneuver being relatively small. The perturbation equations have the form of linear, non-self-adjoint equations with time-dependent coefficients. An active control technique can then be formulated to permit maneuvering of the spacecraft and simultaneously suppressing the elastic vibration.

Meirovitch, L.↗