Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 775 records · Page 43

Recent advances in reduction methods for nonlinear problems

Status and some recent developments in the application of reduction methods to nonlinear structural mechanics problems are summarized. The aspects of reduction methods discussed herein include: (1) selection of basis vectors in nonlinear static and dynamic problems, (2) application of reduction methods in nonlinear static analysis of structures subjected to prescribed edge displacements, and (3) use of reduction methods in conjunction with mixed finite element models. Numerical examples are presented to demonstrate the effectiveness of reduction methods in nonlinear problems. Also, a number of research areas which have high potential for application of reduction methods are identified.

Noor, A. K.↗

A finite element analysis of fatigue crack closure

Experiments have shown that fatigue cracks close at positive loads during constant-amplitude load cycling. The crack-closure phenomenon is caused by residual plastic deformations remaining in the wake of an advancing crack tip. The present paper is concerned with the application of a two-dimensional, nonlinear, finite-element analysis for predicting crack-closure and crack-opening stresses during cyclic loading. A two-dimensional finite-element computer program, which accounts for both elastic-plastic material behavior and changing boundary conditions associated with crack extension and intermittent contact of the crack surfaces under cyclic loading, has been developed. An efficient technique to account for changing boundary conditions was also incorporated into the nonlinear analysis program. This program was subsequently used to study crack extension and crack closure under constant-amplitude and two-level block loading. The calculated crack-closure and crack-opening stresses were qualitatively consistent with experimental observations.

Newman, J. C., Jr.↗

Test and Analysis of the 8-foot Diameter Cylindrical Sandwich Composite Test Article CTA8.2B: As part of the NASA Engineering and Safety Center Shell Buckling Knockdown Factor Project

The buckling response of thin-walled cylindrical structures have long been shown to be sensitive to various imperfections. That is, imperfections related to load introduction, boundary conditions, material property variation, and geometric deviations can all contribute to experimentally obtained buckling loads being lower than the buckling loads predicted for perfect representations of the same thin-walled cylindrical structure. To account for the influence of geometry-based imperfections of a thin-walled cylindrical structure on buckling response, a design factor (also known as a buckling knockdown factor) is used during the design phase. Guidelines for buckling knockdown factors are readily available in NASA SP-8007, Buckling of Thin Walled Circular Cylinders, which was updated in 2020. Prior to this recent update, NASA SP-8007 had not been updated since 1969, and since that time, computational analysis, experimental testing, and manufacturing methods have significantly improved. The NASA Engineering and Safety Center (NESC) Shell Buckling Knockdown Factor Project (SBKF, NESC Assessment 07-010-E) had the goal of developing buckling design recommendations for select classes of metallic and composite launch-vehicle structures [3]. Specifically, the test and analysis results from the fourth SBKF composite test article SBKF-P3-CYL-CTA8.2B, which is referred to as CTA8.2B, are described in this paper. This test was the fourth in a series of four tests on sandwich composite cylinders that will be used to experimentally validate analysis methods, which in turn can be used to develop new analysis-based shell buckling design guidelines for typical sandwich composite launch vehicle cylindrical structures. Publications describing previous test validation with analysis are available in Refs. CTA8.2B was designed to occupy one corner of the design space with approximately equal bending stiffnesses in the axial and transverse directions and have a larger ratio of radius to effective shell thickness. Each test article had a unique design and was designed to occupy a specific region of the design space. CTA8.2B was an 8-foot diameter honeycomb-core sandwich cylinder that was fabricated and tested at the NASA Marshall Space Flight Center (MSFC). The primary objectives of this test were to interrogate the structural capability of a composite test article, and to verify the test article design and analysis approach for cylinders subjected to axial compression and combined axial and bending loads. First, descriptions of the test article design and test are given in Section 2.0, and modeling and analysis methods used in support of the test article design and testing activities are described briefly in Section 3.0. Then, select test results are presented and compared to predicted results in Section 4.0, and concluding remarks are presented in Section 5.0. Finally, drawings of CTA8.2B and references used to support the test and analysis are provided in the Appendix.

Nonlinear analysis↗

Equivalent linearization for fatigue life estimates of a nonlinear structure

An analysis is presented of the suitability of the method of equivalent linearization for estimating the fatigue life of a nonlinear structure. Comparisons are made of the fatigue life of a nonlinear plate as predicted using conventional equivalent linearization and three other more accurate methods. The excitation of the plate is assumed to be Gaussian white noise and the plate response is modeled using a single resonant mode. The methods used for comparison consist of numerical simulation, a probabalistic formulation, and a modification of equivalent linearization which avoids the usual assumption that the response process is Gaussian. Remarkably close agreement is obtained between all four methods, even for cases where the response is significantly linear.

Miles, R. N.↗

Equivalent linearization for fatigue life estimates of a nonlinear structure

An analysis is presented of the suitability of the method of equivalent linearization for estimating the fatigue life of a nonlinear structure. Comparisons are made of the fatigue life of a nonlinear plate as predicted using conventional equivalent linearization and three other more accurate methods. The excitation of the plate is assumed to be Gaussian white noise and the plate response is modeled using a single resonant mode. The methods used for comparison consist of numerical simulation, a probabilistic formulation, and a modification of equivalent linearization which avoids the usual assumption that the response process is Gaussian. Remarkably close agreement is obtained between all four methods, even for cases where the response is significantly linear.

Miles, R. N.↗

Study of stability of large maneuvers of airplanes

A predictive method of nonlinear system analysis is used to investigate airplane stability and dynamic response during rolling maneuvers. The maneuver roll-rate is not assumed to be constant, and the airplane motion is represented by a set of coupled nonlinear differential equations. The general rolling maneuver is kinematically specified by its roll-rate variation p(t). A method for relating the airplane dynamic response to p(t) is developed. The method provides analytical expressions for the motion variables in terms of the maneuver descriptor p(t). A parameterized family of rolling maneuvers is considered, for which the method is used to predict specific dynamic response information, such as the dependence of the peak angle-of-attack excursion on the maneuver parameters. The stability and motion of the airplane in response to an arbitrary actuation of aileron input is considered. Analytical expressions relating motion variables to aileron input are obtained. Explicit analytical bounds on the motion variables are derived. A stability criterion which guarantees nondivergence of motion in response to aileron actuation is presented.

Haddad, E. K.↗

Transient analysis techniques in performing impact and crash dynamic studies

Because of the emphasis being placed on crashworthiness as a design requirement, increasing demands are being made by various organizations to analyze a wide range of complex structures that must perform safely when subjected to severe impact loads, such as those generated in a crash event. The ultimate goal of crashworthiness design and analysis is to produce vehicles with the ability to reduce the dynamic forces experienced by the occupants to specified levels, while maintaining a survivable envelope around them during a specified crash event. DYCAST is a nonlinear structural dynamic finite element computer code that started from the plans systems of a finite element program for static nonlinear structural analysis. The essential features of DYCAST are outlined.

Pifko, A. B.↗

Nonlinear effects in thermal stress analysis of a solid propellant rocket motor

Direct characterization procedures were used to determine the relaxation modulus as a function of time, temperature, and state of strain. Using the quasi-elastic method of linearviscoelasticity, these properties were employed in a finite element computer code to analyze a thick-walled, nonlinear viscoelastic cylinder in the state of plane strain bonded to a thin (but stiff) elastic casing and subjected to slow thermal cooling. The viscoelastic solution is then expressed as a sequence of elastic finite element solutions. The strain-dependent character of the relaxation modulus is included by replacing the single relaxation curve used in the linear viscoelastic theory by a family of relaxation functions obtained at various strain levels. These functions may be regarded as a collection of stress histories or responses to specific loads (in this case, step strains) with which the cooldown solution is made to agree by iterations on the modulus and strain level.

Francis, E. C.↗

TAP 1: A Finite Element Program for Steady-State Thermal Analysis of Convectively Cooled Structures

The program has a finite element library of six elements: two conduction/convection elements to model heat transfer in a solid, two convection elements to model heat transfer in a fluid, and two integrated conduction/convection elements to represent combined heat transfer in tubular and plate/fin fluid passages. Nonlinear thermal analysis due to temperature dependent thermal parameters is performed using the Newton-Raphson iteration method. Program output includes nodal temperatures and element heat fluxes. Pressure drops in fluid passages may be computed as an option. A companion plotting program for displaying the finite element model and predicted temperature distributions is presented. User instructions and sample problems are presented in appendixes.

Thornton, E. A.↗

Precision attitude determination for multimission spacecraft

Attitude determination algorithms for a multimission spacecraft are derived and their performance analyzed. The attitude determination system is composed of a strapdown Inertial Reference Unit (IRU), two fixed head star trackers, and an onboard computer. IRU data is processed to maintain real-time knowledge of spacecraft attitude relative to an inertial reference frame. Star tracker data is processed using Kalman filtering techniques to estimate and correct the attitude determination errors and the gyro drift compensation errors. The results of a star availability analysis for stellar, solar and earth pointing missions are presented. Linear covariance analysis techniques are used to evaluate nominal attitude determination performance, the effects of sensor measurement accuracy variations, the effects of errors in knowledge of sensor measurement accuracy, and the effects of star tracker misalignment errors. Results of a nonlinear simulation analysis of attitude determination performance are also presented. These analyses show that precision attitude determination for stellar, solar and earth pointing missions is achieved.

Murrell, J. W.↗

Stability analysis of the discrete-data Large-Space-Telescope system

This paper considers the application of the discrete describing function to the stability analysis of the Large-Space Telescope (LST) system. An analytical model of the CMG gimbal friction is derived, which is then used to arrive at a closed-form analytical expression for the discrete describing function of the nonlinearity. The analysis is used for the study of the fine-pointing stability of the LST vehicle. Simulation results corroborate the conclusions from the analytical analysis.

Kuo, B. C.↗

Sensitivity Analysis of the Integrated Medical Model for ISS Programs

Sensitivity analysis estimates the relative contribution of the uncertainty in input values to the uncertainty of model outputs. Partial Rank Correlation Coefficient (PRCC) and Standardized Rank Regression Coefficient (SRRC) are methods of conducting sensitivity analysis on nonlinear simulation models like the Integrated Medical Model (IMM). The PRCC method estimates the sensitivity using partial correlation of the ranks of the generated input values to each generated output value. The partial part is so named because adjustments are made for the linear effects of all the other input values in the calculation of correlation between a particular input and each output. In SRRC, standardized regression-based coefficients measure the sensitivity of each input, adjusted for all the other inputs, on each output. Because the relative ranking of each of the inputs and outputs is used, as opposed to the values themselves, both methods accommodate the nonlinear relationship of the underlying model. As part of the IMM v4.0 validation study, simulations are available that predict 33 person-missions on ISS and 111 person-missions on STS. These simulated data predictions feed the sensitivity analysis procedures. The inputs to the sensitivity procedures include the number occurrences of each of the one hundred IMM medical conditions generated over the simulations and the associated IMM outputs: total quality time lost (QTL), number of evacuations (EVAC), and number of loss of crew lives (LOCL). The IMM team will report the results of using PRCC and SRRC on IMM v4.0 predictions of the ISS and STS missions created as part of the external validation study. Tornado plots will assist in the visualization of the condition-related input sensitivities to each of the main outcomes. The outcomes of this sensitivity analysis will drive review focus by identifying conditions where changes in uncertainty could drive changes in overall model output uncertainty. These efforts are an integral part of the overall verification, validation, and credibility review of IMM v4.0.

medical equipment↗

Collapse analysis of a waffle plate strongback for Space Station Freedom

The structural integrity was studied of the Integrated Equipment Assembly (IEA) Strongback of the SSF for the launch environment. The strongback structure supports the electrical power system for SSF. To achieve minimum launch mass, it is essential that flight structures are designed as light as possible. A nonlinear structural analysis was conducted to determine the collapse load of the structure and the associated factor of safety against the service loads. A modeling technique is provided for simulating the load conditions and the buckling and post buckling (collapse) load is evaluated of the IEA Strongback structure, using the finite element computer code MARC. Two of four strongback panels were modeled and analyzed. The effects were examined of the following factors on the global behavior of the strongback panels: (1) load simplification and simulation; (2) type of support boundary conditions; and (3) the possibility of weight reduction of the original structure. For this purpose, several models of the two panels of the strongback were considered. The stress level and distribution in the panels for launch condition, the Eigenvalue critical buckling load and/or the collapse load were determined.

Monasa, Frank F.↗

ALESQP: An Augmented Lagrangian Equality-Constrained SQP Method for Optimization with General Constraints

Here we present a new algorithm for infinite-dimensional optimization with general constraints, called ALESQP. In short, ALESQP is an augmented Lagrangian method that penalizes inequality constraints and solves equality-constrained nonlinear optimization subproblems at every iteration. The subproblems are solved using a matrix-free trust-region sequential quadratic programming (SQP) method that takes advantage of iterative, i.e., inexact linear solvers, and is suitable for large-scale applications. A key feature of ALESQP is a constraint decomposition strategy that allows it to exploit problem-specific variable scalings and inner products. We analyze convergence of ALESQP under different assumptions. We show that strong accumulation points are stationary. Consequently, in finite dimensions ALESQP converges to a stationary point. In infinite dimensions we establish that weak accumulation points are feasible in many practical situations. Under additional assumptions we show that weak accumulation points are stationary. We present several infinite-dimensional examples where ALESQP shows remarkable discretization-independent performance in all of its iterative components, requiring a modest number of iterations to meet constraint tolerances at the level of machine precision. Also, we demonstrate a fully matrix-free solution of an infinite-dimensional problem with nonlinear inequality constraints.

97 MATHEMATICS AND COMPUTING↗

Fully nonlinear development of the most unstable goertler vortex in a three dimensional boundary layer

The nonlinear development is studied of the most unstable Gortler mode within a general 3-D boundary layer upon a suitably concave surface. The structure of this mode was first identified by Denier, Hall and Seddougui (1991) who demonstrated that the growth rate of this instability is O(G sup 3/5) where G is the Gortler number (taken to be large here), which is effectively a measure of the curvature of the surface. Previous researchers have described the fate of the most unstable mode within a 2-D boundary layer. Denier and Hall (1992) discussed the fully nonlinear development of the vortex in this case and showed that the nonlinearity causes a breakdown of the flow structure. The effect of crossflow and unsteadiness upon an infinitesimal unstable mode was elucidated by Bassom and Hall (1991). They demonstrated that crossflow tends to stabilize the most unstable Gortler mode, and for certain crossflow/frequency combinations the Gortler mode may be made neutrally stable. These vortex configurations naturally lend themselves to a weakly nonlinear stability analysis; work which is described in a previous article by the present author. Here we extend the ideas of Denier and Hall (1992) to the three-dimensional boundary layer problem. It is found that the numerical solution of the fully nonlinear equations is best conducted using a method which is essentially an adaption of that utilized by Denier and Hall (1992). The influence of crossflow and unsteadiness upon the breakdown of the flow is described.

Otto, S. R.↗

Robustness Analysis and Optimally Robust Control Design via Sum-of-Squares

A control analysis and design framework is proposed for systems subject to parametric uncertainty. The underlying strategies are based on sum-of-squares (SOS) polynomial analysis and nonlinear optimization to design an optimally robust controller. The approach determines a maximum uncertainty range for which the closed-loop system satisfies a set of stability and performance requirements. These requirements, de ned as inequality constraints on several metrics, are restricted to polynomial functions of the uncertainty. To quantify robustness, SOS analysis is used to prove that the closed-loop system complies with the requirements for a given uncertainty range. The maximum uncertainty range, calculated by assessing a sequence of increasingly larger ranges, serves as a robustness metric for the closed-loop system. To optimize the control design, nonlinear optimization is used to enlarge the maximum uncertainty range by tuning the controller gains. Hence, the resulting controller is optimally robust to parametric uncertainty. This approach balances the robustness margins corresponding to each requirement in order to maximize the aggregate system robustness. The proposed framework is applied to a simple linear short-period aircraft model with uncertain aerodynamic coefficients.

Dorobantu, Andrei↗

Inelastic Analysis of Thermomechanically Cycled Structures

Simplified inelastic analysis computer program (ANSYMP) developed for predicting stress/strain history of thermomechanically cycled structure from an elastic solution. Program uses an iterative and incremental procedure to estimate plastic strains from material stress/strain properties and simulated plasticity hardening model. Program ANSYMP developed to simplify nonlinear structural analysis using only elastic solution as input data.

Kaufman, A.↗

An integrated approach to the synthesis of geometrically non-linear structures

An integrated approach to the minimum weight design of geometrically nonlinear three-dimensional truss structures with geometric imperfections, subject to inequality constraints on static displacements, stresses, local buckling and cross sectional areas, is investigated. The integrated structural synthesis problem involves design and response quantities as independent variables and equilibrium equations, describing the finite element model, as equality constraints. The nonlinear structural analysis and the optimization are thus merged together into a single process. A computer program developed to compute the constraint values and analytical gradients is coupled with a generalized reduced gradient algorithm to solve the integrated problem. Numerical results for a geometrically nonlinear shallow dome example problem are presented for various types of imperfections. Furthermore, it is found that the algorithm is capable of detecting and guarding against system as well as element elastic instability using equilibrium information only, that is, without imposing system and local buckling inequality constraints.

Smaoui, H.↗