Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Finite-element models of continental extension

Numerical models of the initial deformation of extending continental lithosphere, computed to investigate the control of preexisting thermal and mechanical heterogeneities on the style of deformation, are presented. The finite element method is used to calculate deformation with a viscoelastic-plastic model for the lithosphere. Comparisons of the results of analytic models and finite-element models using this method show that good results may be obtained by the numerical technique, even with elements containing both brittle and viscoelastic sampling points. It is shown that the gross style of initial extensional deformation is controlled by the depth and width of the initial heterogeneity which localizes deformation.

Lynch, H. David↗
NASA NTRS1972

Finite element analysis of large transient elastic-plastic deformations of simple structures, with application to the engine rotor fragment containment/deflection problem

Assumed-displacement versions of the finite-element method are developed to predict large-deformation elastic-plastic transient deformations of structures. Both the conventional and a new improved finite-element variational formulation are derived. These formulations are then developed in detail for straight-beam and curved-beam elements undergoing (1) Bernoulli-Euler-Kirchhoff or (2) Timoshenko deformation behavior, in one plane. For each of these categories, several types of assumed-displacement finite elements are developed, and transient response predictions are compared with available exact solutions for small-deflection, linear-elastic transient responses. The present finite-element predictions for large-deflection elastic-plastic transient responses are evaluated via several beam and ring examples for which experimental measurements of transient strains and large transient deformations and independent finite-difference predictions are available.

Wu, R. W.↗
NASA NTRSDate not supplied

Alloy Evaluation and Flow Forming Process Modeling for Net Shape Aerospace Structures

Over the past decade, NASA Langley Research Center (LaRC) has led several manufacturing demonstration projects exploiting flow forming technology. The work has resulted in the commercial-scale manufacture of 10-ft. diameter single-piece, integrally stiffened cylinders. Near-net-shape flow forming offers simplified manufacturing schedules and cost savings through reduced part count. Reduction or elimination of machining, welding, and/or riveting can also lead to significant performance gains. NASA LaRC recently established a flow forming research facility to investigate new alloys and stiffener geometries for fuselage and launch vehicle cryotank applications. Candidate aluminum alloys and heat treatment combinations have been characterized through advanced mechanical testing and microscopy to maximize workability during flow forming trials. Elasto-plastic deformation simulations of the forming process have been performed using the finite element software DEFORM and correlated with forming trial results. The overall objective is to optimize structural performance through a combination of innovative materials, processes, and designs.

single-piece↗
NASA NTRS1983

Analyses of large quasistatic deformations of inelastic bodies by a new hybrid-stress finite element algorithm

A new hybrid-stress finite element algorithm, suitable for analyses of large, quasistatic, inelastic deformations, is presented. The algorithm is base upon a generalization of de Veubeke's complementary energy principle. The principal variables in the formulation are the nominal stress rate and spin, and thg resulting finite element equations are discrete versions of the equations of compatibility and angular momentum balance. The algorithm produces true rates, time derivatives, as opposed to 'increments'. There results a complete separation of the boundary value problem (for stress rate and velocity) and the initial value problem (for total stress and deformation); hence, their numerical treatments are essentially independent. After a fairly comprehensive discussion of the numerical treatment of the boundary value problem, we launch into a detailed examination of the numerical treatment of the initial value problem, covering the topics of efficiency, stability and objectivity. The paper is closed with a set of examples, finite homogeneous deformation problems, which serve to bring out important aspects of the algorithm.

Reed, K. W.↗
NASA NTRS1989

Nonlinear finite element formulation for the large displacement analysis in multibody system dynamics

A total Lagrangian finite element formulation for the deformable bodies in multibody mechanical systems that undergo finite relative rotations is developed. The deformable bodies are discretized using finite element methods. The shape functions that are used to describe the displacement field are required to include the rigid body modes that describe only large translational displacements. This does not impose any limitations on the technique because most commonly used shape functions satisfy this requirement. The configuration of an element is defined using four sets of coordinate systems: Body, Element, Intermediate element, Global. The body coordinate system serves as a unique standard for the assembly of the elements forming the deformable body. The element coordinate system is rigidly attached to the element and therefore it translates and rotates with the element. The intermediate element coordinate system, whose axes are initially parallel to the element axes, has an origin which is rigidly attached to the origin of the body coordinate system and is used to conveniently describe the configuration of the element in undeformed state with respect to the body coordinate system.

Rismantab-Sany, J.↗
NASA NTRS1995

Thermocapillary convection in a rectangular cavity with a deformable interface

A finite-volume method and a boundary element technique are used to compute two-dimensional (2D) thermocapillary convection in a rectangular cavity. The free surface is deformable and the deviations from a flat interface, h, are assumed small in the finite-volume calculations as appropriate for small Capillary (Ca). Two asymptotic approaches are employed; the first is an expansion valid as Ca tends towards 0 and the second assumes h tends towards 0, retaining Ca explicitly as a parameter. On the other hand, the boundary elements approach can be used with O(1) surface deformations. These three different formulations are used to calculate thermocapillary motions in fluids with small Prandtl numbers of about 0.01, Ca less than 0.05, aspect ratios (width/height) 1, 2 and 4, and various values of the Maragoni number (Ma). The same solutions are calculated with these different approaches and are found in good agreement for values of Ca up to 0.05. All solutions calculated are steady, which is both in agreement and disagreement with recently published results by different authors employing different numerical techniques.

Mundrane, M.↗
NASA NTRS1992

Analysis of Large Quasistatic Deformations of Inelastic Solids by a New Stress Based Finite Element Method

A new hybrid stress finite element algorithm suitable for analyses of large quasistatic deformation of inelastic solids is presented. Principal variables in the formulation are the nominal stress rate and spin. The finite element equations which result are discrete versions of the equations of compatibility and angular momentum balance. Consistent reformulation of the constitutive equation and accurate and stable time integration of the stress are discussed at length. Examples which bring out the feasibility and performance of the algorithm conclude the work.

Reed, Kenneth W.↗
NASA NTRS1985

Thermodynamically consistent constitutive equations for nonisothermal large strain, elasto-plastic, creep behavior

The paper is concerned with the development of constitutive relations for large nonisothermal elastic-viscoplastic deformations for metals. The kinematics of elastic-plastic deformation, valid for finite strains and rotations, is presented. The resulting elastic-plastic uncoupled equations for the deformation rate combined with use of the incremental elasticity law permits a precise and purely deductive development of elastic-viscoplastic theory. It is shown that a phenomenological thermodynamic theory in which the elastic deformation and the temperature are state variables, including few internal variables, can be utilized to construct elastic-viscoplastic constitutive equations, which are appropriate for metals. The limiting case of inviscid plasticity is examined.

Riff, R.↗
NASA NTRS2011

A Deployable 4 Meter 180 to 680 GHz Antenna for the Scanning Microwave Limb Sounder

The Scanning Microwave Limb Sounder (SMLS) is a space-borne heterodyne radiometer which will measure pressure, temperature and atmospheric constituents from thermal emission between 180 and 680 GHz. SMLS, planned for the Global Atmospheric Composition Mission of the NRC Decadal Survey, uses a novel toric Cassegrain antenna to perform both elevation and azimuth scanning. These provide better horizontal and temporal resolution and coverage than were possible with elevation-only scanning at typical Low-Earth orbit spacing in the two previous MLS satellite instruments. Development of the SMLS antenna was the focus of a 2006 Small Business Innovative Research (SBIR) program whose phase II culminated in the fabrication and thermal stability testing of a composite demonstration model of the SMLS primary reflector. This reflector has the full 4m height and 1/3 the width planned for flight. An Instrument Incubator Program (IIP) titled "A deployable 4 Meter 180 to 680 GHz antenna for the Scanning Microwave Limb Sounder" continues development of the SMLS antenna with the study of 5 topics: 1) detailed mathematical modeling of the antenna patterns from which we simulate geophysical parameter retrievals in order to establish FOV performance requirements; 2) thorough correlation of finite element model predictions with measurements made on the SBIR reflector. We will again measure deformations of this reflector, under more flight-like thermal gradients, using higher precision metrology techniques available in a new large-aperture facility at JPL; 3) fabrication of a full-width primary reflector whose asbuilt surface figure will better meet the figure requirements of SMLS than did the SBIR reflector; 4) integration of the primary with other reflectors, and with residual front ends built in a 2007 IIP, in a breadboard antenna; and finally 5) RF testing of the breadboard on a Near Field Range at JPL. We report on significant progress in 3 areas of the current IIP: development of the mathematical model to predict SMLS antenna patterns and their application in a preliminary set of geophysical retrievals; the correlation between surface deformation predicted by finite element models and measurement in the 2009 isothermal stability tests of the SBIR, with implications for the thermal gradient testing to be performed at JPL; and aspects of the conceptual design of the full-width primary reflector to be fabricated and tested in the 2nd and 3rd years of the IIP.

heterodyne radiometer↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2022

Validation of material models for puncture of 7075-T651 aluminum plate

Plate puncture simulations are challenging computational tasks that require advanced material models including high strain rate and thermal-mechanical effects on both deformation and failure, plus finite element techniques capable of representing large deformations and material failure. The focus of this work is on the material issues, which require large sets of experiments, flexible material models and challenging calibration procedures. In this study, we consider the puncture of 12.7 mm thick, 7075-T651 aluminum alloy plates by a cylindrical punch with a hemispherical nose and diameter of 12.7 mm. The plasticity and ductile failure models were isotropic with calibration data obtained from uniaxial tension tests at different temperatures and strain rates plus quasi-static notched tension tests and shear-dominated tests described here. Sixteen puncture experiments were conducted to identify the threshold penetration energy, mode of puncture and punch acceleration during impact, The punch was mounted on a 139 kg mass and dropped on the plates with different impact speeds. Since the mass was the same in all tests, the quantity of interest was the impact speed. The axis and velocity of the punch were perpendicular to the plate surface. The mean threshold punch speed was 3.05 m/s, and the mode of failure was plugging by thermal-mechanical shear banding accompanied by scabbing fragments. Application of the material models in simulations of the tests yielded accurate estimates of the threshold puncture speed and of the mode of failure. Time histories of the punch acceleration compared well between simulation and test. Remarkably, the success of the simulations occurred in spite of even the smallest element used being larger than the width of the shear bands.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
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.↗
NASA NTRS1983

Analyses of large quasistatic deformations of inelastic bodies by a new hybrid-stress finite element algorithm - Applications

A new hybrid-stress finite element algorithm suitable for analyzing large quasistatic deformations of inelastic solids is presented and its feasibility and performance are demonstrated with examples. The algorithm provides extremely accurate bifurcation analysis which is stable with respect to variation in the finite element mesh, so long as the same type of element is used in every mesh. When the mesh element is varied, the result changes in a predictable manner. The method does not necessarily lead to an upper or lower bound for the critical load. An explicit forward gradient scheme is used to improve stability and is shown to be useful also for elongation-dominated deformations. The application of the method to the onset of necking in plane extension and to deformation and stress in plane extension of an elasticoviscous fluid with an array of cylindrical voids is given in detail.

Reed, K. W.↗
NASA NTRS2015

Face Sheet/Core Disbond Growth in Honeycomb Sandwich Panels Subjected to Ground-Air-Ground Pressurization and In-Plane Loading

Typical damage modes in light honeycomb sandwich structures include face sheet/core disbonding and core fracture, both of which can pose a threat to the structural integrity of a component. These damage modes are of particular interest to aviation certification authorities since several in-service occurrences, such as rudder structural failure and other control surface malfunctions, have been attributed to face sheet/core disbonding. Extensive studies have shown that face sheet/core disbonding and core fracture can lead to damage propagation caused by internal pressure changes in the core. The increasing use of composite sandwich construction in aircraft applications makes it vitally important to understand the effect of ground-air-ground (GAG) cycles and conditions such as maneuver and gust loads on face sheet/core disbonding. The objective of the present study was to use a fracture mechanics based approach developed earlier to evaluate the loading at the disbond front caused by ground-air-ground pressurization and in-plane loading. A honeycomb sandwich panel containing a circular disbond at one face sheet/core interface was modeled with three-dimensional (3D) solid finite elements. The disbond was modeled as a discrete discontinuity and the strain energy release rate along the disbond front was computed using the Virtual Crack Closure Technique (VCCT). Special attention was paid to the pressure-deformation coupling which can decrease the pressure load within the disbonded sandwich section significantly when the structure is highly deformed. The commercial finite element analysis software, Abaqus/Standard, was used for the analyses. The recursive pressure-deformation coupling problem was solved by representing the entrapped air in the honeycomb cells as filled cavities in Abaqus/Standard. The results show that disbond size, face sheet thickness and core thickness are important parameters that determine crack tip loading at the disbond front. Further, the pressure-deformation coupling was found to have an important load decreasing effect [6]. In this paper, a detailed problem description is provided first. Second, the analysis methodology is presented. The fracture mechanics approach used is described and the specifics of the finite element model, including the fluid-filled cavities, are introduced. Third, the initial model verification and validation are discussed. Fourth, the findings from a closely related earlier study [6] are summarized. These findings provided the basis for the current investigation. Fifth, an aircraft ascent scenario from 0 to 12192 m (0 to 40000 ft) is considered and the resulting crack tip loading at the disbond front is determined. In-plane loading to simulate maneuvers and gust conditions are also considered. Sixth, the results are shown for a curved panel, which was used to simulate potential fuselage applications. Finally, a brief summary of observations is presented and recommendations for improvement are provided.

Chen, Zhi M.↗
NASA NTRS1999

Non-Linear Finite Element Modeling of THUNDER Piezoelectric Actuators

A NASTRAN non-linear finite element model has been developed for predicting the dome heights of THUNDER (THin Layer UNimorph Ferroelectric DrivER) piezoelectric actuators. To analytically validate the finite element model, a comparison was made with a non-linear plate solution using Von Karmen's approximation. A 500 volt input was used to examine the actuator deformation. The NASTRAN finite element model was also compared with experimental results. Four groups of specimens were fabricated and tested. Four different input voltages, which included 120, 160, 200, and 240 Vp-p with a 0 volts offset, were used for this comparison.

Taleghani, Barmac K.↗
DOE OSTI2021

Impact-induced twinning in a magnesium alloy under different stress conditions

We investigate deformation twinning in a highly textured magnesium alloy plate of Mg-3.1%Al-0.9%Zn-0.4%Mn under edge-on impact with in situ, synchrotron-based, ultrafast X-ray diffraction measurements, and corresponding stress states are simulated with the finite element method. Deformation twinning and its anisotropy under the triaxial stress condition are explained by a statistical analysis of resolved shear stress. The critical resolved shear stress criterion is applicable under complicated stress conditions induced by high strain rate impact loading. Three typical simple stress conditions are further explored as verification and application cases: uniaxial-stress, uniaxial-strain, and plane-stress. Extension twinning in the magnesium alloy is prone to occur for impact loading applied perpendicular to the crystallographic c-axis, regardless of the exact stress conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗