Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NONLINEAR PROGRAMMING”

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 667 records · Page 37

FINDS: A fault inferring nonlinear detection system. User's guide

The computer program FINDS is written in FORTRAN-77, and is intended for operation on a VAX 11-780 or 11-750 super minicomputer, using the VMS operating system. The program detects, isolates, and compensates for failures in navigation aid instruments and onboard flight control and navigation sensors of a Terminal Configured Vehicle aircraft in a Microwave Landing System environment. In addition, FINDS provides sensor fault tolerant estimates for the aircraft states which are then used by an automatic guidance and control system to land the aircraft along a prescribed path. FINDS monitors for failures by evaluating all sensor outputs simultaneously using the nonlinear analytic relationships between the various sensor outputs arising from the aircraft point mass equations of motion. Hence, FINDS is an integrated sensor failure detection and isolation system.

Lancraft, R. E.↗

Discrete-State Simulated Annealing For Traveling-Wave Tube Slow-Wave Circuit Optimization

Algorithms based on the global optimization technique of simulated annealing (SA) have proven useful in designing traveling-wave tube (TWT) slow-wave circuits for high RF power efficiency. The characteristic of SA that enables it to determine a globally optimized solution is its ability to accept non-improving moves in a controlled manner. In the initial stages of the optimization, the algorithm moves freely through configuration space, accepting most of the proposed designs. This freedom of movement allows non-intuitive designs to be explored rather than restricting the optimization to local improvement upon the initial configuration. As the optimization proceeds, the rate of acceptance of non-improving moves is gradually reduced until the algorithm converges to the optimized solution. The rate at which the freedom of movement is decreased is known as the annealing or cooling schedule of the SA algorithm. The main disadvantage of SA is that there is not a rigorous theoretical foundation for determining the parameters of the cooling schedule. The choice of these parameters is highly problem dependent and the designer needs to experiment in order to determine values that will provide a good optimization in a reasonable amount of computational time. This experimentation can absorb a large amount of time especially when the algorithm is being applied to a new type of design. In order to eliminate this disadvantage, a variation of SA known as discrete-state simulated annealing (DSSA), was recently developed. DSSA provides the theoretical foundation for a generic cooling schedule which is problem independent, Results of similar quality to SA can be obtained, but without the extra computational time required to tune the cooling parameters. Two algorithm variations based on DSSA were developed and programmed into a Microsoft Excel spreadsheet graphical user interface (GUI) to the two-dimensional nonlinear multisignal helix traveling-wave amplifier analysis program TWA3. The algorithms were used to optimize the computed RF efficiency of a TWT by determining the phase velocity profile of the slow-wave circuit. The mathematical theory and computational details of the DSSA algorithms will be presented and results will be compared to those obtained with a SA algorithm.

Wilson, Jeffrey D.↗

Evaluating the Stability of NASA's Space Launch System with Adaptive Augmenting Control

NASA's baseline Space Launch System (SLS) flight control system (FCS) design includes an adaptive augmenting control (AAC) component that modifies the attitude control system response to provide the classical gain-scheduled control architecture with additional performance and robustness. The NASA Engineering and Safety Center (NESC) teamed with the Space Launch System (SLS) Program to perform a comprehensive assessment of the stability and robustness of the FCS with AAC. This paper provides an overview of the approach, specific analysis techniques, and outcomes that were particularly relevant for the SLS Program. Multiple analysis techniques that specifically target the nonlinear AAC were commissioned as part of this assessment, which was completed outside of the Program's standard design analysis cycle. The following analyses were included, with each technique adding its own valuable insights: Lyapunov-based stability analysis, classical stability analysis with static AAC gain variations, circle criterion-based analysis of the FCS with a time-varying element, time-domain stability margin assessment, Monte Carlo simulations with expanded dispersions, and an extensive set of stressing cases. Several of the completed analyses focused on determining whether the inclusion of AAC introduced risk to the FCS, while others quantified the benefits of the adaptive augmentation.

Tannen S Vanzwieten↗

Marshall system for aerospace system simulation (MARSYAS), user's manual

The capabilities of the Marshall system for aerospace system simulation (MARSYAS) and how to use it are described. MARSYAS is a software system that allows easy setup and control of the simulation of the dynamics of large physical systems on a digital computer. The physical systems are modeled in the form of block diagrams or equations. The blocks can have multiple inputs and multiple outputs, and they can be nested to form hierarchies. The block diagrams can contain transfer functions, nonlinear and logical functions, equations, analog computer elements and FORTRAN programs. The input format of the equations can be combinations of nonlinear, time-varying differential equations and algebraic equations in their original format. MARSYAS could also serve as a storage and retrieval system for models as a basis for a model configuration control system on a central time-shared computer. The outputs of the simulation system can be not only time-responses but also other analysis data such as frequency response, power spectrum and stability parameters. The MARSYAS translator is written in FORTRAN running on the Univac 1108 computer under the EXEC 8 operating system.

Ventre, A.↗

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.↗

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.↗

Estimation of wing nonlinear aerodynamic characteristics at supersonic speeds

A computational system for estimation of nonlinear aerodynamic characteristics of wings at supersonic speeds was developed and was incorporated in a computer program. This corrected linearized theory method accounts for nonlinearities in the variation of basic pressure loadings with local surface slopes, predicts the degree of attainment of theoretical leading edge thrust, and provides an estimate of detached leading edge vortex loadings that result when the theoretical thrust forces are not fully realized.

Carlson, H. W.↗

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.↗

A critical evaluation of the minimization techniques as applied to nonlinear structural analyses

This paper identifies the potential for unconstrained minimization algorithms of mathematical programming to be cost-effective with other conventional techniques of transient nonlinear structural analysis. With this in mind the authors have attempted to critically evaluate a few of the most commonly used algorithms for their effectiveness in solving structural problems involving geometric and material nonlinearities. The algorithms have been categorized as being zeroth order, first order and second order depending upon the order of the derivative of the function called for by the algorithm. The sensitivity of these algorithms to the accuracy of derivatives derived on the basis of finite difference operations clearly suggest using analytically derived derivatives to obtain better control on the number of minimizations required for convergence to the exact solution.

Kamat, M. P.↗

Three-dimensional vibration analysis of a uniform beam with offset inertial masses at the ends

Analysis of a flexible beam with displaced end-located inertial masses is presented. The resulting three-dimensional mode shape is shown to consist of two one-plane bending modes and one torsional mode. These three components of the mode shapes are shown to be linear combinations of trigonometric and hyperbolic sine and cosine functions. Boundary conditions are derived to obtain nonlinear algebraic equations through kinematic coupling of the general solutions of the three governing partial differential equations. A method of solution which takes these boundary conditions into account is also presented. A computer program has been written to obtain unique solutions to the resulting nonlinear algebraic equations. This program, which calculates natural frequencies and three-dimensional mode shapes for any number of modes, is presented and discussed.

Robertson, D. K.↗

Application of GRASP to nonlinear analysis of a cantilever beam

The General Rotorcraft Aeromechanical Stability Program (GRASP) was developed to analyze the steady-state and linearized dynamic behavior of rotorcraft in hovering and axial flight conditions. Because of the nature of problems GRASP was created to solve, the geometrically nonlinear behavior of beams is one area in which the program must perform well in order to be of any value. Numerical results obtained from GRASP are compared to both static and dynamic experimental data obtained for a cantilever beam undergoing large displacements and rotations caused by deformation. The correlation is excellent in all cases.

Hinnant, Howard E.↗

Nonconical Relaxation for Supersonic Potential Flow

Nonlinear, three-dimensional effects computed from full potentialflow equation. Nonconical Relaxation program, NCOREL, employs new computational technique for prediction of inviscid, nonlinear supersonic aerodynamics. Unlike conventional linear potential equations, NCOREL utilizes full potential flow equation to predict formation of supercritical crossflow regions, embedded shocks, and bow shocks. NCOREL written in FORTRAN IV for batch execution.

Siclari, M. J.↗

Formulation and validation of high-order linearized models of helicopter flight mechanics

A high-order linearized model of helicopter flight dynamics is extracted from a nonlinear time domain simulation. The model has 29 states that describe the fuselage rigid body degrees of freedom, the flap and lag dynamics in a nonrotating coordinate system, the inflow dynamics, the delayed entry of the horizontal tail into the main rotor wake, and, approximately, the blade torsion dynamics. The nonlinear simulation is obtained by extensively modifying the GENHEL computer program. The results indicate that the agreement between the linearized and the nonlinear model is good for small perturbations, and deteriorates for large amplitude maneuvers.

Kim, Frederick D.↗

Parameter Design and Optimal Control of an Open Core Flywheel Energy Storage System

In low earth orbit (LEO) satellite applications spacecraft power is provided by photovoltaic cells and batteries. To overcome battery shortcomings the University of Maryland, working in cooperation with NASA/GSFC and NASA/LeRC, has developed a magnetically suspended flywheel for energy storage applications. The system is referred to as an Open Core Composite Flywheel (OCCF) energy storage system. Successful application of flywheel energy storage requires integration of several technologies, viz. bearings, rotor design, motor/generator, power conditioning, and system control. In this paper we present a parameter design method which has been developed for analyzing the linear SISO model of the magnetic bearing controller for the OCCF. The objective of this continued research is to principally analyze the magnetic bearing system for nonlinear effects in order to increase the region of stability, as determined by high speed and large air gap control. This is achieved by four tasks: (1) physical modeling, design, prototyping, and testing of an improved magnetically suspended flywheel energy storage system, (2) identification of problems that limit performance and their corresponding solutions, (3) development of a design methodology for magnetic bearings, and (4) design of an optimal controller for future high speed applications. Both nonlinear SISO and MIMO models of the magnetic system were built to study limit cycle oscillations and power amplifier saturation phenomenon observed in experiments. The nonlinear models include the inductance of EM coils, the power amplifier saturation, and the physical limitation of the flywheel movement as discussed earlier. The control program EASY5 is used to study the nonlinear SISO and MIMO models. Our results have shown that the characteristics and frequency responses of the magnetic bearing system obtained from modeling are comparable to those obtained experimentally. Although magnetic saturation is shown in the bearings, there are good correlations between the theoretical model and experimental data. Both simulation and experiment confirm large variations of the magnetic bearing characteristics due to air gap growth. Therefore, the gap growth effect should be considered in the magnetic bearing system design. Additionally, the magnetic bearing control system will be compared to other design methods using not only parameter design but H-infinity optimal control and mu synthesis.

Pang, D.↗

On local total strain redistribution using a simplified cyclic inelastic analysis based on an elastic solution

Strain redistribution corrections were developed for a simplified inelastic analysis procedure to economically calculate material cyclic response at the critical location of a structure for life prediction purposes. The method was based on the assumption that the plastic region in the structure is local and the total strain history required for input can be defined from elastic finite element analyses. Cyclic stress-strain behavior was represented by a bilinear kinematic hardening model. The simplified procedure has been found to predict stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load cycled problems. This study derived and incorporated Neuber type corrections in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was exercised on a mechanically load cycled benchmark notched plate problem. Excellent agreement was found between the predicted material response and nonlinear finite element solutions for the problem. The simplified analysis computer program used 0.3 percent of the CPU time required for a nonlinear finite element analysis.

Hwang, S. Y.↗

On local total strain redistribution using a simplified cyclic inelastic analysis based on an elastic solution

Strain redistribution corrections were developed for a simplified inelastic analysis procedure to economically calculate material cyclic response at the critical location of a structure for life prediction purposes. The method was based on the assumption that the plastic region in the structure is local and the total strain history required for input can be defined from elastic finite element analyses. Cyclic stress-strain behavior was represented by a bilinear kinematic hardening model. The simplified procedure has been found to predict stress-strain response with reasonable accuracy for thermally cycled problems but needs improvement for mechanically load cycled problems. This study derived and incorporated Neuber type corrections in the simplified procedure to account for local total strain redistribution under cyclic mechanical loading. The corrected simplified method was exercised on a mechanically load cycled benchmark notched plate problem. Excellent agreement was found between the predicted material response and nonlinear finite element solutions for the problem. The simplified analysis computer program used 0.3 percent of the CPU time required for a nonlinear finite element analysis.

Hwang, S. Y.↗

Extension of a nonlinear systems theory to general-frequency unsteady transonic aerodynamic responses

A methodology for modeling nonlinear unsteady aerodynamic responses, for subsequent use in aeroservoelastic analysis and design, using the Volterra-Wiener theory of nonlinear systems is presented. The methodology is extended to predict nonlinear unsteady aerodynamic responses of arbitrary frequency. The Volterra-Wiener theory uses multidimensional convolution integrals to predict the response of nonlinear systems to arbitrary inputs. The CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code is used to generate linear and nonlinear unit impulse responses that correspond to each of the integrals for a rectangular wing with a NACA 0012 section with pitch and plunge degrees of freedom. The computed kernels then are used to predict linear and nonlinear unsteady aerodynamic responses via convolution and compared to responses obtained using the CAP-TSD code directly. The results indicate that the approach can be used to predict linear unsteady aerodynamic responses exactly for any input amplitude or frequency at a significant cost savings. Convolution of the nonlinear terms results in nonlinear unsteady aerodynamic responses that compare reasonably well with those computed using the CAP-TSD code directly but at significant computational cost savings.

Silva, Walter A.↗

Simulation of Non-Linear Transonic Aeroelastic Behavior on the B-2

At high subsonic flight speeds, large flexible aircraft begin to encounter unsteady airloads which are not predicted by most currently available aerodynamic analysis and design methods. With increasing speed and the development of transonic flow and shocks, viscous effects quickly become very important, and flow separation can occur. The Northrop Grumman USAF B-2 Bomber encountered a nonlinear aeroelastic Residual Pitch Oscillation (RPO) under these conditions. Simulation studies were performed with the Computational Aeroelasticity Program-Transonic Small Disturbance, Viscous (CAPTSDv) computer program to evaluate its ability to predict these nonlinear aeroelastic responses. Open and closed loop simulations were performed to assess the participation of the flight control system. Control, system actuator hysteresis characteristics were modeled and found to be a significant participant in the RPO phenomenon. Simulations were also performed for varying Mach numbers and altitudes to establish the stability boundaries and compare with flight test data. These CAPTSDv simulations compared well with flight data and revealed many potential further modeling enhancements.

Residual pitch oscillation↗