Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear 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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Buckling Behavior of Compression-Loaded Quasi-Isotropic Curved Panels with a Circular Cutout

Results from a numerical and experimental study of the response of compression-loaded quasi-isotropic curved panels with a centrally located circular cutout are presented. The numerical results were obtained by using a geometrically nonlinear finite element analysis code. The effects of cutout size, panel curvature and initial geo- metric imperfections on the overall response of compression-loaded panels are described. In addition, results are presented from a numerical parametric study that indicate the effects of elastic circumferential edge restraints on the prebuckling and buckling response of a selected panel and these numerical results are compared to experimentally measured results. These restraints are used to identify the effects of circumferential edge restraints that are introduced by the test fixture that was used in the present study. It is shown that circumferential edge restraints can introduce substantial nonlinear prebuckling deformations into shallow compression-loaded curved panels that can results in a significant increase in buckling load.

Hilburger, Mark W.↗
DOE OSTI2026

Algebraic Multigrid with Filtering: An Efficient Preconditioner for Interior Point Methods in Large-Scale Contact Mechanics Optimization

Large-scale contact mechanics simulations are crucial in many engineering fields such as structural design and manufacturing. In the frictionless case, contact can be modeled by minimizing an energy functional; however, these problems are often nonlinear, nonconvex, and increasingly difficult to solve as mesh resolution increases. In this work, we employ a Newton-based interior-point (IP) filter line-search method, an effective approach for large-scale constrained optimization. While this method converges rapidly, each iteration requires solving a large saddle-point linear system that becomes ill-conditioned as the optimization process converges, largely due to IP treatment of the contact constraints. Such ill-conditioning can hinder solver scalability and increase iteration counts with mesh refinement. Here, to address this, we introduce a novel preconditioner, algebraic multigrid with filtering (AMGF), tailored to the Schur complement of the saddle-point system. Building on the classical AMG solver, commonly used for elasticity, we augment it with a specialized subspace correction that filters near null space components introduced by contact interface constraints. Through theoretical analysis and numerical experiments on a range of linear and nonlinear contact problems, we demonstrate that the proposed solver achieves mesh independent convergence and maintains robustness against the ill-conditioning that notoriously plagues IP methods. These results indicate that AMGF makes contact mechanics simulations more tractable and broadens the applicability of Newton-based IP methods in challenging engineering scenarios. More broadly, AMGF is well suited for problems, optimization or otherwise, where solver performance is limited by a low-dimensional subspace, such as those arising from localized constraints, interface conditions, or model heterogeneities. This makes the method widely applicable beyond contact mechanics and constrained optimization.

Mathematics and Computing↗
DOE OSTI2023

Deep material network via a quilting strategy: visualization for explainability and recursive training for improved accuracy

Recent developments integrating micromechanics and neural networks offer promising paths for rapid predictions of the response of heterogeneous materials with similar accuracy as direct numerical simulations. The deep material network is one such approaches, featuring a multi-layer network and micromechanics building blocks trained on anisotropic linear elastic properties. Once trained, the network acts as a reduced-order model, which can extrapolate the material’s behavior to more general constitutive laws, including nonlinear behaviors, without the need to be retrained. However, current training methods initialize network parameters randomly, incurring inevitable training and calibration errors. Here, we introduce a way to visualize the network parameters as an analogous unit cell and use this visualization to “quilt” patches of shallower networks to initialize deeper networks for a recursive training strategy. The result is an improvement in the accuracy and calibration performance of the network and an intuitive visual representation of the network for better explainability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Assessing Dynamic Time Warping Techniques for Discriminating Seismic Sources at Local and Regional Distances

Effective monitoring of seismic explosions and hazard assessment relies heavily on the accurate discrimination of underground seismic sources. This study investigates the application of novel nonlinear alignment techniques, specifically Dynamic Time Warping (DTW), for event-type discrimination at regional and local distances. Building on prior research that used DTW and Elastic Shape Analysis (ESA) in discrimination at regional distances, we evaluate the performance of recently developed variants of DTW, including a method that employs Pearson cross-correlation as a measure of warping distance and a time distortion coefficient that quantifies the type and degree of time distortion between signals. By analyzing observational datasets that include different source types, we assess the performance of these approaches for realistic monitoring scenarios. Specifically, we consider a dataset recorded at regional distances in the Korean Peninsula and a local-distance subset from the Unconstrained Utah Event Bulletin catalog to evaluate DTW-based discrimination across multiple distance scales. Additionally, we introduce the maximum cross-correlations of warped waveforms as a similarity metric for event classification. Through hierarchical cluster analysis and dendrogram interpretation, we present our findings, highlighting the strengths and limitations of these techniques in seismic event classification.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
NASA NTRS1987

Simultaneous analysis and design for eigenvalue maximization

A simultaneous analysis and design approach to the maximization of buckling or vibration eigenvalues is presented. Both unimodal and bimodal optimum solutions are considered. A discretization of the structure and response is used to obtain sets of nonlinear algebraic equations which are solved numerically. The formulation is applied to the solution of the optimum design of a beam-column supported by an elastic foundation for maximum buckling load. Two cases are considered: the optimum design of the column for a given foundation and the optimum design of the foundation for a given column. Results are compared to published solutions when possible.

Plaut, Raymond H.↗
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.↗
NASA NTRS1991

Experimental and analytical analysis of stress-strain behavior in a (90/0 deg)2s, SiC/Ti-15-3 laminate

The nonlinear stress strain behavior of 90 degree/0 degree sub 2s, SiC/Ti-15-3 composite laminate was numerically investigated with a finite element, unit cell approach. Tensile stress-strain curves from room temperature experiments depicted three distinct regions of deformation, and these regions were predicted by finite element analysis. The first region of behavior, which was linear elastic, occurred at low applied stresses. As applied stresses increased, fiber/matrix debonding in the 90 degree plies caused a break in the stress-strain curve and initiated a second linear region. In this second region, matrix plasticity in the 90 degree plies developed. The third region, which was typified by nonlinear, stress-strain behavior occr red at high stresses. In this region, the onset of matrix plasticity in the 0 degree plies stiffened the laminate in the direction transverse to the applied load. Metallographic sections confirmed the existence of matrix plasticity in specific areas of the structure. Finite element analysis also predicted these locations of matrix slip.

Lerch, Bradley A.↗
DOE OSTI2021

Unexpected Scaling of Peak Acceleration for a Yielding Mass-Spring System Subjected to a Triangular Base Acceleration Pulse

The use of bounding scenarios is a common practice that greatly simplifies the design and qualification of structures. However, this approach implicitly assumes that the quantities of interest increase monotonically with the input to the structure, which is not necessarily true for nonlinear structures. This paper surveys the literature for observations of nonmonotonic behavior of nonlinear systems and finds such observations in both the earthquake engineering and applied mechanics literature. Numerical simulations of a single degree-of-freedom mass-spring system with an elastic–plastic spring subjected to a triangular base acceleration pulse are then presented, and it is shown that the relative acceleration of this system scales nonmonotonically with the input magnitude in some cases. The equation of motion for this system is solved symbolically and an approximate expression for the relative acceleration is developed, which qualitatively agrees with the nonmonotonic behavior seen in the numerical results. The nonmonotonicity is investigated and found to be a result of dynamics excited by the discontinuous derivative of the base acceleration pulse, the magnitude of which scales nonmonotonically with the input magnitude due to the fact that the first yield of the spring occurs earlier as the input magnitude is increased. The relevance of this finding within the context of defining bounding scenarios is discussed, and it is recommended that modeling be used to perform a survey of the full range of possible inputs prior to defining bounding scenarios.

42 ENGINEERING↗
NASA NTRS1966

Waves in materials with memory. V.

One-dimensional acceleration waves and higher order waves propagating in general nonlinear Maxwellian materials with fading memory

VISCOELASTICITY↗
NASA NTRS1972

Finite element analysis of elasto-plastic soils. Report no. 4: Finite element analysis of elasto-plastic frictional materials for application to lunar earth sciences

The behavior of finite element models employing different constitutive relations to describe the stress-strain behavior of soils is investigated. Three models, which assume small strain theory is applicable, include a nondilatant, a dilatant and a strain hardening constitutive relation. Two models are formulated using large strain theory and include a hyperbolic and a Tresca elastic perfectly plastic constitutive relation. These finite element models are used to analyze retaining walls and footings. Methods of improving the finite element solutions are investigated. For nonlinear problems better solutions can be obtained by using smaller load increment sizes and more iterations per load increment than by increasing the number of elements. Suitable methods of treating tension stresses and stresses which exceed the yield criteria are discussed.

Marr, W. A., Jr.↗
NASA NTRS1973

Some conclusions regarding the aeroelastic stability of hingeless helicopter blades in hover and in forward flight

In this paper results and conclusions obtained from the study of the aeroelastic instability of hingeless helicopter blades are presented. First, the large amplitude coupled flap-lag equations of motion of a hingeless elastic helicopter blade are solved using an asymptotic expansion procedure in multiple time scales. Both hover and forward flight cases are considered. Stability boundaries and amplitudes of nonlinear response are obtained. From these, the importance of the nonlinear coupling and the effect of the periodic coefficients is determined. Next, using a system of linearized coupled flap-lag-pitch equations in hover, various divergence mechanisms for hingeless blades are shown. Finally, the flutter boundaries for coupled flap-lag-pitch are obtained. The effect of the torsional degree of freedom on the flap-lag type of instability is investigated. Similarly the effect of lag on the flap-pitch type of instability is considered. In addition, the effect of various blade parameters on the stability boundaries is shown.

Friedmann, P.↗
NASA NTRS1973

Liapunov stability analysis of hybrid dynamical systems in the neighborhood of nontrivial equilibrium

This paper is concerned with the stability of a hybrid dynamical system in the neighborhood of a nontrivial equilibrium, where the system consists of one rigid part and n elastic members. The body moves in a central-force field with its mass center describing a circular orbit. The nontrivial equilibrium is defined by steady rotation of the system at an angular velocity equal to the orbital velocity, with the elastic members being in deformed state. A Liapunov stability analysis is performed by assuming small perturbations about the nontrivial equilibrium, where the latter is generally defined by nonlinear differential equations. The theory is applied to a gravity-gradient stabilized satellite with flexible appendages.

Meirovitch, L.↗
NASA NTRS1989

Computing Stress, Stability, And Vibration Of Shells

BOSOR4 computer program developed as comprehensive program for analysis of stress, stability, and vibration of complex, branched shells of revolution made of elastic materials. Used to analyze prismatic shells and panels. Performs large-deflection axisymmetric stress analysis, small-deflection nonsymmetric stress analysis, modal vibration analysis with axisymmetric nonlinear prestress included, and buckling analysis with axisymmetric or nonsymmetric prestress. One main advantage, provisions for such realistic engineering details as eccentric load paths, internal supports, arbitrary branching conditions, and "library" of wall constructions. Based on finite-difference energy method and offers very rapid convergence with increasing numbers of mesh points. Written in FORTRAN 77.

Bushnell, David↗
NASA NTRS2002

Nonlinear Viscoelastic Mechanics of Cross-linked Rubbers

The paper develops a general theory for finite rubber viscoelasticity, and specifies it in the form, convenient for solving problems important for rubber, tire and space industries. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general nonlinear theory has been developed for arbitrary nonisothermal deformations of viscoelastic solids. In this theory, the constitutive equations are presented as the sum of known equilibrium (rubber elastic) and non-equilibrium (liquid polymer viscoelastic) terms. These equations are then simplified using several modeling arguments. Stability constraints for the proposed constitutive equations are also discussed. It is shown that only strong ellipticity criteria are applicable for assessing stability of the equations governing viscoelastic solids.

Freed, Alan D.↗
DOE OSTI2021

Elastic-instability–enabled locomotion

Locomotion of an organism interacting with an environment is the consequence of a symmetry-breaking action in space-time. Here we show a minimal instantiation of this principle using a thin circular sheet, actuated symmetrically by a pneumatic source, using pressure to change shape nonlinearly via a spontaneous buckling instability. This leads to a polarized, bilaterally symmetric cone that can walk on land and swim in water. In either mode of locomotion, the emergence of shape asymmetry in the sheet leads to an asymmetric interaction with the environment that generates movement––via anisotropic friction on land, and via directed inertial forces in water. Scaling laws for the speed of the sheet of the actuator as a function of its size, shape, and the frequency of actuation are consistent with our observations. The presence of easily controllable reversible modes of buckling deformation further allows for a change in the direction of locomotion in open arenas and the ability to squeeze through confined environments––both of which we demonstrate using simple experiments. Our simple approach of harnessing elastic instabilities in soft structures to drive locomotion enables the design of novel shape-changing robots and other bioinspired machines at multiple scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Jump phenomena

The paper considers jump phenomena composed of large amplitude responses of nonlinear systems caused by small amplitude disturbances. Physical problems where large jumps in the solution amplitude are important features of the response are described, including snap buckling of elastic shells, chemical reactions leading to combustion and explosion, and long-term climatic changes of the earth's atmosphere. A new method of rational functions was then developed which consists of representing the solutions of the jump problems as rational functions of the small disturbance parameter; this method can solve jump problems explicitly.

Reiss, E. L.↗
NASA NTRS1984

Calculated effects of varying Reynolds Number and dynamic pressure on flexible wings at transonic speeds

A computational method is described that includes the effects of static aeroelastic wing deflections in steady transonic aerodynamic calculations. This method, known as the Transonic Aero-elastic Program System (TAPS), interacts a 3D transonic computer code with boundary layer and a linear finite element structural analysis codes to calculate wing pressures and deflections. The nonlinear nature of the transonic flow makes it necessary to couple the aerodynamic and structures codes in an iterative manner. TAPS has been arranged in a modular fashion so that different aerodynamic or structures programs may be used with a minimum of coding changes required. Results obtained using two different aerodynamic codes in TAPS are given, and those results are correlated with experimental data.

Campbell, R. L.↗