Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear structural analysis”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Nonlinear dynamics, fractals, cardiac physiology and sudden death

The authors propose a diametrically opposite viewpoint to the generally accepted tendency of equating healthy function with order and disease with chaos. With regard to the question of sudden cardiac death and chaos, it is suggested that certain features of dynamical chaos related to fractal structure and fractal dynamics may be important organizing principles in normal physiology and that certain pathologies, including ventricular fibrillation, represent a class of 'pathological periodicities'. Some laboratory work bearing on the relation of nonlinear analysis to physiological and pathophysiological data is briefly reviewed, with tentative theories and models described in reference to the mechanism of ventricular fibrillation.

Goldberger, Ary L.↗
NASA NTRS2023

Seasonality of the Migrating Semidiurnal Tide in the Tropical Upper Mesosphere and Lower Thermosphere and its Thermodynamic and Momentum Budget

This work uses the Specified Dynamics-Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X) to determine and explain the seasonality of the migrating semidiurnal tide (SW2) components of tropical upper mesosphere and lower thermosphere (UMLT) temperature, zonal wind, and meridional wind. This work also quantifies aliasing due to SW2 in satellite-based tidal estimates. Results show that during equinox seasons, the vertical profiles of tropical UMLT temperature SW2 and zonal-wind SW2’s amplitudes have a double-peak structure while they, along with meridional-wind SW2, have a single-peak structure in their amplitudes in June solstice. Hough mode reconstruction reveals that a linear combination of five SW2 Hough modes cannot fully reproduce these features. Tendency analysis reveals that for temperature, the adiabatic term, nonlinear advection term, and linear advection term are important. For the winds, the classical terms, nonlinear advection term, linear advection term, and gravity wave drag are important. Results of our alias analysis then indicate that SW2 can induce an ∼60% alias in zonal-mean and DW1 components calculated from sampling like that of the Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics satellite and the Aura satellite. This work concludes that in situ generation by wave–wave interaction and/or by gravity waves plays significant roles in the seasonality of tropical UMLT temperature SW2, zonal-wind SW2, and meridional-wind SW2. The alias analysis further adds that one cannot simply assume that SW2 in the tropical UMLT is negligible.

Tides↗
DOE OSTI2025

Nuclear Materials Packaging, Transportation, and Systems Analysis Group Software Quality Assurance Plan: ANSYS Mechanical Finite Element Analysis Software Version 2023R1

ANSYS Inc. develops and markets engineering simulation software and services used in the aerospace, automotive, manufacturing, electronics, biomedical, energy, defense, and many other industries. ANSYS is dedicated to engineering simulation and is the world’s leading software provider. ANSYS was founded in 1970 and is headquartered in Canonsburg, Pennsylvania. ANSYS provides an engineering analysis tool combining structural, thermal, computational fluid dynamics, acoustic, and electromagnetic simulation capabilities. ANSYS has two main programs, which use the same solvers: (1) Mechanical APDL (ANSYS Design Parametric Language), a Fortran-based coding platform, and (2) ANSYS Workbench, which uses a graphical user interface to aid in finite element analysis implementation. This plan covers both APDL and Workbench. The ANSYS computer program is a large-scale, multipurpose finite element program that can be used to solve several classes of engineering analyses. The analysis capabilities of ANSYS include the ability to solve static and dynamic structural analyses, steady-state and transient heat transfer problems, mode-frequency and buckling eigenvalue problems, static or time-varying magnetic analyses, and various types of field and coupled-field applications. The program contains many special features that allow nonlinearities or secondary effects such as plasticity, large strain, hyperelasticity, creep, swelling, large deflections, contact, stress stiffening, temperature dependency, material anisotropy, and radiation to be included in the solution. As ANSYS has been developed, other special capabilities such as substructuring, submodeling, random vibration, kinetostatics, kinetodynamics, free convection fluid analysis, acoustics, magnetics, piezoelectrics, coupled-field analysis, and design optimization have been added to the program. These capabilities contribute further to making ANSYS a multipurpose analysis tool for varied engineering disciplines. The ANSYS program has been in commercial use for over 50 years and has been used extensively in the aerospace, automotive, construction, electronic, energy services, manufacturing, nuclear, plastics, oil, and steel industries. Additionally, many consulting firms and hundreds of universities have used ANSYS for analysis, research, and educational purposes. ANSYS is recognized worldwide as one of the most widely used and capable programs of its type. Ansys design analysis software is the first created within a quality system with ISO 9001 certification, the internationally accepted quality standard. Product development, testing, maintenance and support processes also meet the United States Nuclear Regulatory Commission's quality requirements, as they have for nearly four decades. The Quality Assurance Service Agreement is suitable for the customers working in the nuclear industry who need to meet specific federal regulations including 10CRF50 Appendix B and provisions of 10CFR21. ANSYS has retained its original International Organization for Standardization (ISO) 9001 accreditation certificate since1995-05-04, It’s current certificate is valid until 2027-05-29.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Exponential Decay in the Sensitivity Analysis of Nonlinear Dynamic Programming

In this report, we study the sensitivity of discrete-time dynamic programs with nonlinear dynamics and objective to perturbations in the initial conditions and reference parameters. Under uniform controllability and boundedness assumptions for the problem data, we prove that the directional derivative of the optimal state and control at time $k$, $x^*_k$, and $u^*_k$, with respect to the reference signal at time $i$, $d_i$, will have exponential decay in terms of $|k-i|$ with a decay rate $\rho$ independent of the temporal horizon length. The key technical step is to prove that a version of the convexification approach proposed by Verschueren et al. can be applied to the KKT conditions and results in a convex quadratic program with uniformly bounded data. In turn, Riccati techniques can be further employed to obtain the sensitivity result, borne from the observation that the directional derivatives are solutions of quadratic programs with structure similar to the KKT conditions themselves. We validate our findings with numerical experiments on a small nonlinear, nonconvex, dynamic program.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. In this work, a set of nonlinear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The nonlinear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP β increases, the EP β threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the nonlinear destabilization of the 3/4 and 2/3 energetic particle modes, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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.↗
NASA NTRS2023

An Enabling Platform for Achieving Multiscale Multiphysics Analysis of Multiphase Materials

This paper will address current NASA GRC research efforts utilizing the NASA Multiscale Analysis Tool (NASMAT) which address technical gaps in three of the Vision 2040 key discipline areas. NASMAT is a versatile platform for performing computationally efficient multiscale analyses of heterogeneous materials and is available free through the NASA software catalog. It offers the user flexibility to define an arbitrary number of length scales (levels) where a variety of micromechanics theories can be implemented at each level. Micromechanics theories can be selected to balance accuracy and computational efficiency and range from analytical (Mori-Tanaka) to several semi-analytical (method of cells) formulations. The resulting anisotropic, evolving nonlinear, thermomechanical constitutive model can also be coupled with external software and used to perform multiscale analyses of more complex structures. A recent application to model a complex, three-dimensional woven composite, with a particular emphasis placed on multiscale measurements utilized to enhance the quality of the multiscale analysis will be discussed. Since typical NASMAT analyses can be completed on the order of seconds to minutes, a second example will demonstrate the ability to generate large quantities of data useful for sensitivity analysis, uncertainty quantification, or machine learning applications. Current progress on implementing NASMAT within a multiscale digital thread/digital twin framework will also be addressed. These examples will demonstrate that NASMAT is an enabling platform for achieving the goals in Vision 2040.

Multiscale↗
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 NTRS1997

A Nonlinear Programming Perspective on Sensitivity Calculations for Systems Governed by State Equations

This paper discusses the calculation of sensitivities. or derivatives, for optimization problems involving systems governed by differential equations and other state relations. The subject is examined from the point of view of nonlinear programming, beginning with the analytical structure of the first and second derivatives associated with such problems and the relation of these derivatives to implicit differentiation and equality constrained optimization. We also outline an error analysis of the analytical formulae and compare the results with similar results for finite-difference estimates of derivatives. We then attend to an investigation of the nature of the adjoint method and the adjoint equations and their relation to directions of steepest descent. We illustrate the points discussed with an optimization problem in which the variables are the coefficients in a differential operator.

Lewis, Robert Michael↗
NASA NTRS2003

Balloon Borne Arc-Second Pointer Feasibility Study

For many years scientists have been utilizing stratospheric balloons as low-cost platforms on which to conduct space science experiments. A major hurdle in extending the range of experiments for which these vehicles are useful has been the imposition of the gondola dynamics on the accuracy with which an instrument can be kept pointed at a celestial target. A significant number of scientists have sought the ability to point their instruments with jitter in the arc-second range. This paper presents the design and analysis of a stratospheric balloon borne pointing system that is able to meet this requirement. The foundation for a high fidelity controller simulation is presented. The flexibility of the flight train is represented through generalized modal analysis. A multiple controller scheme is introduced for coarse and fine pointing. Coarse azimuth pointing is accomplished by an established pointing system, with extensive flight history, residing above the gondola structure. A pitch-yaw gimbal mount is used for fine pointing, providing orthogonal axes when nominally on target. Fine pointing actuation is from direct drive dc motors, eliminating backlash problems. An analysis of friction nonlinearities and a demonstration of the necessity in eliminating static fiction are provided. A unique bearing hub design is introduced that eliminates static fiction from the system dynamics. A control scheme involving linear accelerometers for enhanced disturbance rejection is also presented. Results from a linear analysis of the total system and the high fidelity simulation are given. This paper establishes that the proposed control strategy can be made robustly stable with significant design margins. Also demonstrated is the efficacy of the proposed system in rejecting disturbances larger than those considered realistic. Finally, we see that sub arc-second pointing stability can be achieved for a large instrument pointing at an inertial target.

Ward, Philip R.↗
NASA NTRS2002

Packaging Technology for SiC High Temperature Circuits Operable up to 500 Degrees Centigrade

New high temperature low power 8-pin packages have been fabricated using commercial fabrication service. These packages are made of aluminum nitride and 96 percent alumina with Au metallization. The new design of these packages provides the chips inside with EM shielding. Wirebond geometry control has been achieved for precise mechanical tests. Au wirebond samples with 45 degree heel-angle have been tested using wireloop test module. The geometry control improves the consistency of measurement of the wireloop breaking point.Also reported on is a parametric study of the thermomechanical reliability of a Au thick-film based SiC die-attach assembly using nonlinear finite element analysis (FEA) was conducted to optimize the die-attach thermo-mechanical performance for operation at temperatures from room temperature to 500 degrees Centigrade. This parametric study centered on material selection, structure design and process control.

Chen, Lian-Yu↗
NASA NTRS2006

Demonstration of a Balloon Borne Arc-Second Pointer Design

Many designs for utilizing stratospheric balloons as low-cost platforms on which to conduct space science experiments have been proposed throughout the years. A major hurdle in extending the range of experiments for which these vehicles are useful has been the imposition of the gondola dynamics on the accuracy with which an instrument can be kept pointed at a celestial target. A significant number of scientists have sought the ability to point their instruments with jitter in the arc-second range. This paper presents the design and analysis of a stratospheric balloon borne pointing system that is able to meet this requirement. The test results of a demonstration prototype of the design with similar ability are also presented. Discussion of a high fidelity controller simulation for design analysis is presented. The flexibility of the flight train is represented through generalized modal analysis. A multiple controller scheme is utilized for coarse and fine pointing. Coarse azimuth pointing is accomplished by an established pointing system, with extensive flight history, residing above the gondola structure. A pitch-yaw gimbal mount is used for fine pointing, providing orthogonal axes when nominally on target. Fine pointing actuation is from direct drive dc motors, eliminating backlash problems. An analysis of friction nonlinearities and a demonstration of the necessity in eliminating static friction are provided. A unique bearing hub design is introduced that eliminates static friction from the system dynamics. A control scheme involving linear accelerometers for enhanced disturbance rejection is also presented. Results from a linear analysis of the total system and the high fidelity simulation are given. Results from a generalized demonstration prototype are presented. Commercial off-the-shelf (COTS) hardware was used to demonstrate the efficacy and performance of the pointer design for a mock instrument. Sub-arcsecond pointing ability from a ground hang test setup is shown from the testing results. This paper establishes that the proposed control strategy can be made robustly stable with significant design margins. Also demonstrated is the efficacy of the proposed system in rejecting disturbances larger than those considered realistic. The system is implemented and demonstrates sub arc second pointing ability using COTS hardware. Finally, we see that sub arc-second pointing stability can be achieved for a large instrument pointing at an inertial target.

DeWeese, Keith D.↗
DOE OSTI2020

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Global nonlinear optimization of spacecraft protective structures design

The global optimization of protective structural designs for spacecraft subject to hypervelocity meteoroid and space debris impacts is presented. This nonlinear problem is first formulated for weight minimization of the space station core module configuration using the Nysmith impact predictor. Next, the equivalence and uniqueness of local and global optima is shown using properties of convexity. This analysis results in a new feasibility condition for this problem. The solution existence is then shown, followed by a comparison of optimization techniques. Finally, a sensitivity analysis is presented to determine the effects of variations in the systemic parameters on optimal design. The results show that global optimization of this problem is unique and may be achieved by a number of methods, provided the feasibility condition is satisfied. Furthermore, module structural design thicknesses and weight increase with increasing projectile velocity and diameter and decrease with increasing separation between bumper and wall for the Nysmith predictor.

Mog, R. A.↗
NASA NTRS1982

Effects of internal heat transfer on the structure of self-similar blast waves

An analysis of the problem of self-similar, nonadiabatic blast waves, where both conduction and radiation are allowed to take place, show the problem to be reducible to the integration of a system of six coupled nonlinear ordinary differential equations. Consideration of these equations shows that although radiation tends to produce uniform fields through temperature gradient attenuation, all the energy carried by radiation is deposited on the front and the bounding shock becomes increasingly overdriven. When conduction is taken into account, the distribution of gasdynamic parameters in blast waves in the case of Rosseland diffusion radiation is more uniform than in the case of the Planck emission radiation.

Ghoniem, A. F.↗
NASA NTRS1980

Finite-element nonlinear transient response computer programs PLATE 1 and CIVM-PLATE 1 for the analysis of panels subjected to impulse or impact loads

Two computer programs are described for predicting the transient large deflection elastic viscoplastic responses of thin single layer, initially flat unstiffened or integrally stiffened, Kirchhoff-Lov ductile metal panels. The PLATE 1 program pertains to structural responses produced by prescribed externally applied transient loading or prescribed initial velocity distributions. The collision imparted velocity method PLATE 1 program concerns structural responses produced by impact of an idealized nondeformable fragment. Finite elements are used to represent the structure in both programs. Strain hardening and strain rate effects of initially isotropic material are considered.

Spilker, R. L.↗
DOE OSTI2026

Multiresolution Quantum Chemistry: Nonlinear Response Properties at the Basis Set Limit

We benchmark the accuracy of Dunning correlation-consistent Gaussian basis sets for computing frequencydependent second-order hyperpolarizabilities relevant to second-harmonic generation (SHG), using multiresolution analysis (MRA) as a reference. Basis set errors are analyzed using a unit-sphere representation of the effective hyperpolarizability vector, enabling direct assessment of directional error structure. We introduce a relative RMS total error metric that integrates directional deviations over the unit sphere and complement it with signed projection errors that distinguish over- and underestimation. Unsupervised clustering based on these signed directional metrics reveals four distinct convergence behaviors across a set of 68 molecules. Unitsphere visualizations of representative systems show that basis set errors are often highly anisotropic and localized along specific bond directions, even when global error measures appear small. Doubly augmented basis sets consistently outperform singly augmented ones, and core-polarization functions are required for uniform convergence in second-row systems. Overall, this work demonstrates that directional analysis combined with clustering provides a robust framework for understanding basis set convergence in nonlinear optical response properties.

Basis sets↗
DOE OSTI2025

Heterogeneous Mixtures of Dictionary Functions to Approximate Subspace Invariance in Koopman Operators: Why Deep Koopman Operators Work

Abstract Koopman operators model nonlinear dynamics as a linear dynamic system acting on a nonlinear function as the state. This nonstandard state is often called a Koopman observable and is usually approximated numerically by a superposition of functions drawn from a dictionary . In a widely used algorithm, extended dynamic mode decomposition (EDMD), the dictionary functions are drawn from a fixed class of functions. Deep learning combined with EDMD has been used to learn novel dictionary functions in an algorithm called deep dynamic mode decomposition (deepDMD). The learned representation both (1) accurately models and (2) scales well with the dimension of the original nonlinear system. In this paper, we analyze the learned dictionaries from deepDMD and explore the theoretical basis for their strong performance. We explore State-Inclusive Logistic Lifting (SILL) dictionary functions to approximate Koopman observables. Error analysis of these dictionary functions show they satisfy a property of subspace approximation, which we define as uniform finite approximate closure. Typically, a Koopman dictionary’s nonlinear functions are homogeneous. In this paper, we discover that structured mixing of heterogeneous dictionary functions drawn from different classes of nonlinear functions achieve the same accuracy and dimensional scaling as the deep-learning-based deepDMD algorithm Yeung et al. ( In: 2019 American Control Conference (ACC), 2019). We specifically show this by building a heterogeneous dictionary comprised of SILL functions and conjunctive radial basis functions (RBFs). This mixed dictionary achieves similar accuracy and dimensional scaling to deepDMD with an order of magnitude reduction in parameters, while maintaining geometric interpretability. These results strengthen the viability of dictionary-based Koopman models to solving high-dimensional nonlinear learning problems.

Johnson, Charles A.↗