Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear 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 919 records · Page 51

Nature of the wiggle instability of galactic spiral shocks

Gas in disk galaxies interacts nonlinearly with an underlying stellar spiral potential to form galactic spiral shocks. While numerical simulations typically show that spiral shocks are unstable to wiggle instability (WI) even in the absence of magnetic fields and self-gravity, its physical nature has remained uncertain. To clarify the mechanism behind the WI, we conduct a normal-mode linear stability analysis and nonlinear simulations assuming that the disk is isothermal and infinitesimally thin. We find that the WI is physical, originating from the generation of potential vorticity at a deformed shock front, rather than Kelvin-Helmholtz instabilities as previously thought. Since gas in galaxy rotation periodically passes through the shocks multiple times, the potential vorticity can accumulate successively, setting up a normal mode that grows exponentially with time. Eigenfunctions of the WI decay exponentially downstream from the shock front. Both shock compression of acoustic waves and a discontinuity of shear across the shock stabilize the WI. The wavelength and growth time of the WI depend on the arm strength quite sensitively. When the stellar-arm forcing is moderate at 5%, the wavelength of the most unstable mode is about 0.07 times the arm-to-arm spacing, with the growth rate comparable to the orbital angular frequency, which is found to be in good agreement with the results of numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

A unified perspective on robot control - The energy Lyapunov function approach

A unified framework for the stability analysis of robot tracking control is presented. By using an energy-motivated Lyapunov function candidate, the closed-loop stability is shown for a large family of control laws sharing a common structure of proportional and derivative feedback and a model-based feedforward. The feedforward can be zero, partial or complete linearized dynamics, partial or complete nonlinear dynamics, or linearized or nonlinear dynamics with parameter adaptation. As result, the dichotomous approaches to the robot control problem based on the open-loop linearization and nonlinear Lyapunov analysis are both included in this treatment. Furthermore, quantitative estimates of the trade-offs between different schemes in terms of the tracking performance, steady state error, domain of convergence, realtime computation load and required a prior model information are derived.

Wen, John T.↗

Minimum-mass design of filamentary composite panels under combined loads: Design procedure based on a rigorous buckling analysis

A procedure is presented for designing uniaxially stiffened panels made of composite material and subjected to combined inplane loads. The procedure uses a rigorous buckling analysis and nonlinear mathematical programing techniques. Design studies carried out with the procedure consider hat-stiffened and corrugated panels made of graphite-epoxy material. Combined longitudinal compression and shear and combined longitudinal and transverse compression are the loadings used in the studies. The capability to tailor the buckling response of a panel is also explored. Finally, the adequacy of another, simpler, analysis-design procedure is examined.

Stroud, W. J.↗

Structural analysis and sizing of stiffened, metal matrix composite panels for hypersonic vehicles

The present method for strength and stability analyses of stiffened, fiber-reinforced composite panels to be used in hypersonic vehicle structures is of great generality, and can be linked with planar finite-element analysis (FEA). Nonlinear temperature and load-dependent material data for each laminate are used to 'build-up' the stiffened panel's membrane, bending, and membrane-bending coupling stiffness terms, as well as thermal coefficients. The resulting, FEA-solved thermomechanical forces and moments are used to calculate strain at any location in the panel; this allows an effective ply-by-ply orthotropic strength analysis to be conducted, together with orthotropic instability checks for each laminated segment of the cross-section.

Collier, Craig S.↗

Analysis of gain scheduled control for nonlinear plants

Gain scheduling has proven to be a successful design methodology in many engineering applications. In the absence of a sound theoretical analysis, these designs come with no guarantees of the robustness, performance, or even nominal stability of the overall gain-scheduled design. An analysis is presented for two types of nonlinear gain-scheduled control systems: (1) scheduling on a reference trajectory and (2) scheduling on the plant output. Conditions which guarantee stability, robustness, and performance properties of the global gain schedule designs are given. These conditions confirm and formalize popular notions regarding gain-scheduled designs, such as that the scheduling variable should vary slowly and capture the plant's nonlinearities.

Shamma, Jeff S.↗

Nonlinear elasticity and damping govern ultrafast dynamics in click beetles

Many small animals use springs and latches to overcome the mechanical power output limitations of their muscles. Click beetles use springs and latches to bend their bodies at the thoracic hinge and then unbend extremely quickly, resulting in a clicking motion. When unconstrained, this quick clicking motion results in a jump. While the jumping motion has been studied in depth, the physical mechanisms enabling fast unbending have not. Here, we first identify and quantify the phases of the clicking motion: latching, loading, and energy release. We detail the motion kinematics and investigate the governing dynamics (forces) of the energy release. We use high-speed synchrotron X-ray imaging to observe and analyze the motion of the hinge’s internal structures of four Elater abruptus specimens. We show evidence that soft cuticle in the hinge contributes to the spring mechanism through rapid recoil. Using spectral analysis and nonlinear system identification, we determine the equation of motion and model the beetle as a nonlinear single-degree-of-freedom oscillator. Quadratic damping and snap-through buckling are identified to be the dominant damping and elastic forces, respectively, driving the angular position during the energy release phase. The methods used in this study provide experimental and analytical guidelines for the analysis of extreme motion, starting from motion observation to identifying the forces causing the movement. Here, the tools demonstrated here can be applied to other organisms to enhance our understanding of the energy storage and release strategies small animals use to achieve extreme accelerations repeatedly.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling the effects of exercise during 100% oxygen prebreathe on the risk of hypobaric decompression sickness

BACKGROUND: Several previous studies indicated that exercise during prebreathe with 100% O2 decreased the incidence of hypobaric decompression sickness (DCS). We report a meta-analysis of these investigations combined with a new study in our laboratory to develop a statistical model as a predictive tool for DCS. HYPOTHESIS: Exercise during prebreathe increases N2 elimination in a theoretical 360-min half-time compartment decreasing the incidence of DCS. METHODS: A dose-response probability tissue ratio (TR) model with 95% confidence limits was created for two groups, prebreathe with exercise (n = 113) and resting prebreathe (n = 113), using nonlinear regression analysis with maximum likelihood optimization. RESULTS: The model predicted that prebreathe exercise would reduce the residual N2 in a 360-min half-time compartment to a level analogous to that in a 180-min compartment. This finding supported the hypothesis. The incidence of DCS for the exercise prebreathe group was significantly decreased (Chi-Square = 17.1, p < 0.0001) from the resting prebreathe group. CONCLUSIONS: The results suggested that exercise during prebreathe increases tissue perfusion and N2 elimination approximately 2-fold and markedly lowers the risk of DCS. Based on the model, the prebreathe duration may be reduced from 240 min to a predicted 91 min for the protocol in our study, but this remains to be verified. The model provides a useful planning tool to develop and test appropriate prebreathe exercise protocols and to predict DCS risks for astronauts.

NASA Program Environmental Health↗

Failure of wooden sandwich beam reinforced with glass/epoxy faces

The mechanical properties and the failure of wooden beam strengthened with two faces from glass/epoxy composite and a wooden beam without strengthening was studied. Stresses and deflections on both beams, which are imposed in three point bending loading. On the idealized geometry of the specimens with detailed nonlinear orthotropic analysis was performed with a finite elements program. The failure study of the wooden beams was performed, applying the criterion of Tsai-Hill. The shear strength of the adhesive was taken into account. All the specimens were tested with three point bending loading and the experimental results were compared to those of the theoretical approach with the finite elements analysis. Comparing the results, the advantage of strengthened wooden beam against the simple wooden beam becomes obvious. Theoretical predictions were in good agreement with experimental results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A critical examination of stresses in an elastic single lap joint

The results of an approximate nonlinear finite-element analysis of a single lap joint are presented and compared with the results of a linear finite-element analysis, and the geometric nonlinear effects caused by the load-path eccentricity on the adhesive stress distributions are determined. The results from finite-element, Goland-Reissner, and photoelastic analyses show that for a single lap joint the effect of the geometric nonlinear behavior of the joint has a sizable effect on the stresses in the adhesive. The Goland-Reissner analysis is sufficiently accurate in the prediction of stresses along the midsurface of the adhesive bond to be used for qualitative evaluation of the influence of geometric or material parametric variations. Detailed stress distributions in both the adherend and adhesive obtained from the finite-element analysis are presented to provide a basis for comparison with other solution techniques.

Cooper, P. A.↗

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exact nonlinear excitations in double-degenerate plasmas

In this work, we use the conventional hydrodynamics formalism and incorporate the Chew-Goldberger-Low double-adiabatic theory to evaluate the nonlinear electrostatic ion excitations in double-degenerate (electron spin-orbit degenerate) magnetized quantum plasmas. Based on the Sagdeev pseudopotential method, an exact general pseudopotential is calculated which leads to the allowed Mach-number range criteria for such localized density structures in an anisotropic magnetized plasma. We employ the criteria on the Mach-number range for diverse magnetized quantum plasma with different equations of state. It is remarked that various plasma fractional parameters such as the system dimensionality, ion-temperature, relativistic-degeneracy, Zeeman-energy, and plasma composition are involved in the stability of an obliquely propagating nonlinear ion-acoustic wave in a double-degenerate quantum plasma. Current study is most appropriate for nonlinear wave analysis in dense astrophysical magnetized plasma environments such as white-dwarfs and neutron-star crusts where the strong magnetic fields can be present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High Performance-Multidisciplinary Simulation for Regional Scale Earthquake Hazard and Risk Assessments

Moderate-to-large magnitude earthquakes represent a significant risk around the world and a major issue across the Department of Energy (DOE) mission interests, ranging from DOE's one-of-a-kind mission critical facilities to national energy systems (nuclear power plants, electric/gas distribution systems, and renewable energy production facilities). As a result, for the DOE there is a strong motivation in promoting and developing a risk-informed and performance-based approach for both the assessment of existing facilities and the design of new structures. DOE National Labs have historically made multidisciplinary contributions to the solution of challenging seismic hazard and risk problems and have been at the forefront for the development of design standards for major facilities. The current approach for estimating the ground motion hazard for critical facilities is based upon a Probabilistic Seismic Hazard Analysis (PSHA), whereby measured ground motion data are collected and homogenized into a probabilistic framework to estimate the hazard for future events. However, recent studies have highlighted the site-specific nature of earthquakes and shown the limitations of a purely empirical approach. Moreover, the design of infrastructure systems has been moving toward a performance-based approach, in which specific performance goals are specified at different levels of excitation and for different limit states. Such an approach requires the availability of robust and reliable computational tools for the nonlinear dynamic analysis of complex structural systems. Ongoing research at Lawrence Berkeley National Laboratory (LBNL) and Lawrence Livermore National Laboratory (LLNL) aims at developing an advanced simulation environment at the exascale for a coupled assessment of earthquake hazard (ground motions) and earthquake risk (structural demands)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A comparative study of diffraction of shallow-water waves by high-level IGN and GN equations

This work is on the nonlinear diffraction analysis of shallow-water waves, impinging on submerged obstacles, by two related theories, namely the classical Green–Naghdi (GN) equations and the Irrotational Green–Naghdi (IGN) equations, both sets of equations being at high levels and derived for incompressible and inviscid flows. Recently, the high-level Green–Naghdi equations have been applied to some wave transformation problems. The high-level IGN equations have also been used in the last decade to study certain wave propagation problems. However, past works on these theories used different numerical methods to solve these nonlinear and unsteady sets of differential equations and at different levels. Moreover, different physical problems have been solved in the past. Therefore, it has not been possible to understand the differences produced by these two sets of theories and their range of applicability so far. We are thus motivated to make a direct comparison of the results produced by these theories by use of the same numerical method to solve physically the same wave diffraction problems. We focus on comparing these two theories by using similar codes; only the equations used are different but other parts of the codes, such as the wave-maker, damping zone, discretion method, matrix solver, etc., are exactly the same. This way, we eliminate many potential sources of differences that could be produced by the solution of different equations. The physical problems include the presence of various submerged obstacles that can be used for example as breakwaters or to represent the continental shelf. A numerical wave tank is created by placing a wavemaker on one end and a wave absorbing beach on the other. The nonlinear and unsteady sets of differential equations are solved by the finite-difference method. The results are compared with different equations as well as with the available experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The 3D inelastic analysis methods for hot section components

A two-year program to develop advanced 3D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades, and vanes is described. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulation models were developed: an eight-noded midsurface shell element; a nine-noded midsurface shell element; and a twenty-noded isoparametric solid element. A separate computer program has been developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.

Mcknight, R. L.↗

On 3D inelastic analysis methods for hot section components

Accomplishments are described for the 2-year program, to develop advanced 3-D inelastic structural stress analysis methods and solution strategies for more accurate and cost effective analysis of combustors, turbine blades and vanes. The approach was to develop a matrix of formulation elements and constitutive models. Three constitutive models were developed in conjunction with optimized iterating techniques, accelerators, and convergence criteria within a framework of dynamic time incrementing. Three formulations models were developed; an eight-noded mid-surface shell element, a nine-noded mid-surface shell element and a twenty-noded isoparametric solid element. A separate computer program was developed for each combination of constitutive model-formulation model. Each program provides a functional stand alone capability for performing cyclic nonlinear structural analysis. In addition, the analysis capabilities incorporated into each program can be abstracted in subroutine form for incorporation into other codes or to form new combinations.

Mcknight, R. L.↗

Space Shuttle Main Engine structural analysis and data reduction/evaluation. Volume 2: High pressure oxidizer turbo-pump turbine end bearing analysis

The high-pressure oxidizer turbopump (HPOTP) consists of two centrifugal pumps, on a common shaft, that are directly driven by a hot-gas turbine. Pump shaft axial thrust is balanced in that the double-entry main inducer/impeller is inherently balanced and the thrusts of the preburner pump and turbine are nearly equal but opposite. Residual shaft thrust is controlled by a self-compensating, non-rubbing, balance piston. Shaft hang-up must be avoided if the balance piston is to perform properly. One potential cause of shaft hang-up is contact between the Phase 2 bearing support and axial spring cartridge of the HPOTP main pump housing. The status of the bearing support/axial spring cartridge interface is investigated under current loading conditions. An ANSYS version 4.3, three-dimensional, finite element model was generated on Lockheed's VAX 11/785 computer. A nonlinear thermal analysis was then executed on the Marshall Space Flight Center Engineering Analysis Data System (EADS). These thermal results were then applied along with the interference fit and bolt preloads to the model as load conditions for a static analysis to determine the gap status of the bearing support/axial spring cartridge interface. For possible further analysis of the local regions of HPOTP main pump housing assembly, detailed ANSYS submodels were generated using I-DEAS Geomod and Supertab (Appendix A).

Sisk, Gregory A.↗

Discrete Time McKean–Vlasov Control Problem: A Dynamic Programming Approach

We consider the stochastic optimal control problem of nonlinear mean-field systems in discrete time. We reformulate the problem into a deterministic control problem with marginal distribution as controlled state variable, and prove that dynamic programming principle holds in its general form. We apply our method for solving explicitly the mean-variance portfolio selection and the multivariate linear-quadratic McKean–Vlasov control problem.

97 MATHEMATICS AND COMPUTING↗

Analytical and finite element solutions of some problems using a viscoplastic model

The feasibility of using a viscoplastic model developed by Robinson to perform a nonlinear structural analysis was investigated. The paper presents analytical solutions for three classical problems: (1) a pressurized thick-walled cylinder, (2) a thin rotating disk, and (3) a pressurized thick-walled sphere. The analytical expressions derived for the stress and the strain rates for these components are general in nature as they consider both the mechanical and thermal loadings to be time-dependent. A computer program VISTAN (VIscoplastic STructural ANalyzer) was developed to obtain the stress and strain distributions. The finite element solutions for these problems are also presented. The numerical results pertaining to isothermal loading conditions are provided. The results obtained demonstrate the feasibility of using the viscoplastic model developed by Robinson to perform nonlinear structural analyses.

Arya, V. K.↗