Search NASA⌕ Search

SEARCH · Search NASA

Results for “large 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 19 records
All sourcesNASA NTRSDOE OSTI
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 NTRS1972

Large-deformation modal coordinates for nonrigid vehicle dynamics

The derivation of minimum-dimension sets of discrete-coordinate and hybrid-coordinate equations of motion of a system consisting of an arbitrary number of hinge-connected rigid bodies assembled in tree topology is presented. These equations are useful for the simulation of dynamical systems that can be idealized as tree-like arrangements of substructures, with each substructure consisting of either a rigid body or a collection of elastically interconnected rigid bodies restricted to small relative rotations at each connection. Thus, some of the substructures represent elastic bodies subjected to small strains or local deformations, but possibly large gross deformations, in the hybrid formulation, distributed coordinates referred to herein as large-deformation modal coordinates, are used for the deformations of these substructures. The equations are in a form suitable for incorporation into one or more computer programs to be used as multipurpose tools in the simulation of spacecraft and other complex electromechanical systems.

Likins, P. W.↗
DOE OSTI2025

A Simple, Scalable Large Deformation Solid Mechanics Implementation in the MOOSE Framework

This article describes a large deformation solid mechanics solver implemented as part of the freely available and open source MOOSE finite element simulation framework. The article documents the choices made in developing the solid mechanics framework and describes novel formulations for the gradient operator and constitutive modeling framework made to simplify implementations of different coordinate systems, stabilized gradient operators, and different constitutive model inputs and outputs. In the process, the article describes a new formulation that casts objective integration of the Cauchy stress as a linear transformation of the small stress rate. Finally, the article presents key implementation details and examines the parallel efficiency of the solid mechanics solver implemented in MOOSE. The implementation retains a good weak scaling efficiency beyond 1,000 parallel processes. The article includes a discussion of the factors limiting the parallel efficiency of implicit, large deformation solid mechanics codes on current high-performance computers, with the main current limitation being the scalability of the algebraic multigrid methods used to solve the linearized equilibrium equations.

Applied computing → Computer-aided design↗
NASA NTRS1988

An internal variable constitutive model for the large deformation of metals at high temperatures

The advent of large deformation finite element methodologies is beginning to permit the numerical simulation of hot working processes whose design until recently has been based on prior industrial experience. Proper application of such finite element techniques requires realistic constitutive equations which more accurately model material behavior during hot working. A simple constitutive model for hot working is the single scalar internal variable model for isotropic thermal elastoplasticity proposed by Anand. The model is recalled and the specific scalar functions, for the equivalent plastic strain rate and the evolution equation for the internal variable, presented are slight modifications of those proposed by Anand. The modified functions are better able to represent high temperature material behavior. The monotonic constant true strain rate and strain rate jump compression experiments on a 2 percent silicon iron is briefly described. The model is implemented in the general purpose finite element program ABAQUS.

Brown, Stuart↗
NASA NTRS1986

Large Deformation Dynamic Bending of Composite Beams

Studies were conducted on the large deformation response of composite beams subjected to a dynamic axial load. The beams were loaded with a moderate eccentricity to promote bending. The study was primarily experimental but some finite element results were obtained. Both the deformation and the failure of the beams were of interest. The static response of the beams was also studied to determine potential differences between the static and dynamic failure. Twelve different laminate types were tested. The beams were loaded dynamically with a gravity driven impactor traveling at 19.6 ft/sec and quasi-static tests were conducted on identical beams in a displacement controlled manner. For laminates of practical interest, the failure modes under static and dynamic loadings were identical. Failure in most of the laminate types occurred in a single event involving 40% to 50% of the plies. However, failure in laminates with 30 deg or 15 deg off-axis plies occured in several events. All laminates exhibited bimodular elastic properties. Using empirically determined flexural properties, a finite element analysis was reasonably accurate in predicting the static and dynamic deformation response.

Derian, E. J.↗
NASA NTRS1986

Large Deformation Dynamic Bending of Composite Beams

Studies were conducted on the large deformation response of composite beams subjected to a dynamic axial load. The beams were loaded with a moderate eccentricity to promote bending. The study was primarily experimental but some finite element results were obtained. Both the deformation and the failure of the beams were of interest. The static response of the beams was also studied to determine potential differences between the static and dynamic failure. Twelve different laminate types were tested. The beams tested were 23 in. by 2 in. and generally 30 plies thick. The beams were loaded dynamically with a gravity-driven impactor traveling at 19.6 ft/sec and quasi-static tests were conducted on identical beams in a displacement controlled manner. For laminates of practical interest, the failure modes under static and dynamic loadings were identical. Failure in most of the laminate types occurred in a single event involving 40% to 50% of the plies. However, failure in laminates with 300 or 150 off-axis plies occurred in several events. All laminates exhibited bimodular elastic properties. The compressive flexural moduli in some laminates was measured to be 1/2 the tensile flexural modulus. No simple relationship could be found among the measured ultimate failure strains of the different laminate types. Using empirically determined flexural properties, a finite element analysis was reasonably accurate in predicting the static and dynamic deformation response.

Derian, E. J.↗
NASA NTRS1991

Large Deformation Behavior of Long Shallow Cylindrical Composite Panels

An exact solution is presented for the large deformation response of a simply supported orthotropic cylindrical panel subjected to a uniform line load along a cylinder generator. The cross section of the cylinder is circular and deformations up to the fully snapped through position are investigated. The orthotropic axes are parallel to the generator and circumferential directions. The governing equations are derived using laminated plate theory, nonlinear strain-displacement relations, and applying variational principles. The response is investigated for the case of a panel loaded exactly at midspan and for a panel with the load offset from midspan. The mathematical formulation is one dimensional in the circumferential coordinate. Solutions are obtained in closed-form. An experimental apparatus was designed to load the panels. Experimental results of displacement controlled tests performed on graphite-epoxy curved panels are compared with analytical predictions.

Carper, Douglas M.↗
NASA NTRS1983

Large deformation analysis of layered composite shells

A three-dimensional element based on the total Lagrangian description of large deformations of a layered anisotropic composite medium is developed, validated, and used to analyze layered composite shells. The element contains the following features: geometric nonlinearity, dynamic (transient) behavior, and arbitrary lamination scheme and lamina properties. Numerical results of nonlinear bending, natural vibration, and transient response are presented to illustrate the capabilities of the element.

Chao, W. C.↗
DOE OSTI2025

A cohesive zone treatment for the material point method involving problems of large deformation and damage

A new algorithm is described that permits the use of cohesive zones in the material point method for problems involving large deformation and fracture. In contrast to previous cohesive zone implementations, this method does not utilize massless surface-element particles. Instead, cohesive tractions are computed using the shape function mappings from a reference grid configuration in combination with explicitly defined particle surface normals and surface positions. These normals and relative surface positions are updated each time step according to particle deformation. The tractions are converted to cohesive forces using the nodal areas and mapped back to particles using the same reference shape function mappings. These forces are then remapped by conventional particle-to-grid interpolation as external forces using the current-configuration shape-function mappings. This allows highly compliant cohesive zones to function over jump displacements larger than a grid cell. Upon damage, these interfaces can revert to conventional multi-field contact surfaces. This approach is general and readily applies to two and three dimensions as well as being compatible with damage-field gradient partitioning offering exceptional computational flexibility. The framework for this method enables other capabilities, such as improved contact precision using explicitly defined surface normals and positions, and a method to mitigate spurious material damage at weak discontinuities between stiff brittle materials and soft or compliant materials.

Cohesive zone↗
NASA NTRS1993

Influence of transverse-shear and large-deformation effects on the low-speed impact response of laminated composite plates

An analytical procedure is presented for determining the transient response of simply supported, rectangular laminated composite plates subjected to impact loads from airgun-propelled or dropped-weight impactors. A first-order shear-deformation theory is included in the analysis to represent properly any local short-wave-length transient bending response. The impact force is modeled as a locally distributed load with a cosine-cosine distribution. A double Fourier series expansion and the Timoshenko small-increment method are used to determine the contact force, out-of-plane deflections, and in-plane strains and stresses at any plate location due to an impact force at any plate location. The results of experimental and analytical studies are compared for quasi-isotropic laminates. The results indicate that using the appropriate local force distribution for the locally loaded area and including transverse-shear-deformation effects in the laminated plate response analysis are important. The applicability of the present analytical procedure based on small deformation theory is investigated by comparing analytical and experimental results for combinations of quasi-isotropic laminate thicknesses and impact energy levels. The results of this study indicate that large-deformation effects influence the response of both 24- and 32-ply laminated plates, and that a geometrically nonlinear analysis is required for predicting the response accurately.

Ambur, Damodar R.↗
NASA NTRS1984

Hybrid stress finite elements for large deformations of inelastic solids

A new hybrid stress finite element algorithm, based on a generalization of Fraeijs de Veubeke's complementary energy principle is presented. Analyses of large quasistatic deformation of inelastic solids (hypoelastic, plastic, viscoplastic) are within its capability. Principle variables in the formulation are the nominal stress rate and spin. A brief account is given of the boundary value problem in these variables, and the 'equivalent' variational principle. The finite element equation, along with initial positions and stresses, comprise an initial value problem. Factors affecting the choice of time integration schemes are discussed. Results found by application of the new algorithm are compared to those obtained by a velocity based finite element algorithm.

Reed, K. W.↗
NASA NTRS1973

Hybrid-coordinate spacecraft dynamics using large-deformation modal coordinates

For a vehicle amenable to idealization as a set of n + 1 rigid bodies interconnected by n line hinges (implying tree topology), a set of minimum-dimension matrix equations of motion is recorded in terms of variables including the relative rotations of contiguous bodies. When a substructure or portion of the vehicle is subjected to large gross deformations but is restricted to 'small' strains or local deformations, the n-body model is (for sufficiently large n) assumed to experience 'small' relative rotations of contiguous bodies within the substructure, and the equations of motion are linearized in these variables. The application of coordinate transformations and truncations to the linearized subset of the system equations of motion then produces a hybrid-coordinate formulation of the problem involving both discrete and distributed coordinates.

Likins, P. W.↗
DOE OSTI2026

Large deformations of Tr(Φ 3 ) and the world at infinity

The amplitudes of the non-linear sigma model can be obtained from those of Tr(Φ 3 ) theory by sending the kinematic (Mandelstam) variables to infinity in a certain direction. In this paper we characterize the behavior of Tr(Φ 3 ) amplitudes under a general class of large kinematic shifts called g-vector shifts. The objects that live in this world at infinity retain certain key amplitude-like properties, most notably factorization, and admit descriptions in terms of polytopes, but they are not generally amplitudes of any cognizable theory. We identify particular g-vector shifts that lead at infinity to mixed amplitudes involving two pions and any number of scalars, allowing us to provide polytopal descriptions of these amplitudes.

effective field theories↗