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 1,063 records · Page 59

Effect of gimbal friction modelling technique on control stability and performance for Centaur upper stage

The powered-phase autopilot for the Centaur upper stage rocket uses an autopilot forward loop gain scheduler that decreases the proportional gain as propellant mass is depleted. Nonlinear time response simulation studies revealed that Centaur vehicles with low-gain autopilots would have large attitude error limit cycles. These limit cycles were due to the assumed presence of Coulomb friction in the engine gimbals. This situation could be corrected through the use of an harmonic dither, programmed into the on-board digital computer and added to the engine command signal. This would introduce impending motion to the engines, allowing control of the engines even under small commands. Control authority was found to be restored when dither was used. A concern arose that the Centaur could be unacceptably excited at resonances near the dither frequency, if the dither amplitude was to be chosen on the basis of friction level present, a test was conducted to measure this level. Dither characteristics were to be based on the test results. The test results showed that the gimbal friction characteristic was actually hysteretic rather than the assumed Coulomb friction. The simulation results showed that, using this new model of gimbal friction, dither would no longer be necessary.

Graham, Ronald E.↗

Parallel Computation of the Jacobian Matrix for Nonlinear Equation Solvers Using MATLAB

Demonstrating speedup for parallel code on a multicore shared memory PC can be challenging in MATLAB due to underlying parallel operations that are often opaque to the user. This can limit potential for improvement of serial code even for the so-called embarrassingly parallel applications. One such application is the computation of the Jacobian matrix inherent to most nonlinear equation solvers. Computation of this matrix represents the primary bottleneck in nonlinear solver speed such that commercial finite element (FE) and multi-body-dynamic (MBD) codes attempt to minimize computations. A timing study using MATLAB's Parallel Computing Toolbox was performed for numerical computation of the Jacobian. Several approaches for implementing parallel code were investigated while only the single program multiple data (spmd) method using composite objects provided positive results. Parallel code speedup is demonstrated but the goal of linear speedup through the addition of processors was not achieved due to PC architecture.

Rose, Geoffrey K.↗

Experimental and Analytical Studies for a Computational Materials Program

The studies supported by Grant NAG1-1780 were directed at providing physical data on polymer behavior that would form the basis for computationally modeling these types of materials. Because of ongoing work in polymer characterization this grant supported part of a larger picture in this regard. Efforts went into two combined areas of their time dependent mechanical response characteristics: Creep properties on the one hand, subject to different volumetric changes (nonlinearly viscoelastic behavior) and time or frequency dependence of dilatational material behavior. The details of these endeavors are outlined sufficiently in the two appended publications, so that no further description of the effort is necessary.

Knauss, W. G.↗

Optimal horizontal guidance law for aircraft in the terminal area.

A complete nonlinear control law is derived for guiding an aircraft in minimum time from an arbitrary initial position and heading to a prescribed terminal position and heading in the horizontal plane. The solution is obtained with the aid of the maximum principle, and is implemented by constructing a digital-computer program for the resulting switching logic.

Pecsvaradi, T.↗

Analytical and numerical investigation of structural response of compliant wall materials

Surface motion of compliant walls in drag reduction experiments was analyzed. The spectrum of surface motion indicates that membranes over deep cavities respond at low frequencies and large wavelengths. The membrane over a deep cavity is therefore found not to yield the desired reponse predicted by the postulated mechanism. The membrane over a thin air gap is found to act as a wavelength chopper, and analysis of the nonlinear response of the compliant surface indicates its possible suitability for compliant wall experiments. Periodic structures are found to lock in the desired wavelengths of motion. Laminated structures are found to be very ineffective as compliant models, except when there is no bonding between the membrane and the backing. Computer programs developed for these analyses are documented.

Goglia, G. L.↗

Structural dynamics analysis

Dynamic analysis of large space antenna systems must treat the deployment as well as vibration and control of the deployed antenna. Candidate computer programs for deployment dynamics, and issues and needs for future program developments are reviewed. Some results for mast and hoop deployment are also presented. Modeling of complex antenna geometry with conventional finite element methods and with repetitive exact elements is considered. Analytical comparisons with experimental results for a 15 meter hoop/column antenna revealed the importance of accurate structural properties including nonlinear joints. Slackening of cables in this antenna is also a consideration. The technology of designing actively damped structures through analytical optimization is discussed and results are presented.

Housner, J. M.↗

An application of high authority/low authority control and positivity

Control Dynamics Company (CDy), in conjunction with NASA Marshall Space Flight Center (MSFC), has supported the U.S. Air Force Wright Aeronautical Laboratory (AFWAL) in conducting an investigation of the implementation of several DOD controls techniques. These techniques are to provide vibration suppression and precise attitude control for flexible space structures. AFWAL issued a contract to Control Dynamics to perform this work under the Active Control Technique Evaluation for Spacecraft (ACES) Program. The High Authority Control/Low Authority Control (HAC/LAC) and Positivity controls techniques, which were cultivated under the DARPA Active Control of Space Structures (ACOSS) Program, were applied to a structural model of the NASA/MSFC Ground Test Facility ACES configuration. The control systems design were accomplished and linear post-analyses of the closed-loop systems are provided. The control system designs take into account effects of sampling and delay in the control computer. Nonlinear simulation runs were used to verify the control system designs and implementations in the facility control computers. Finally, test results are given to verify operations of the control systems in the test facility.

Seltzer, S. M.↗

A survey of propfan noise propagation at cruise

Acoustic measurements collected during the Propfan Test Assessment (PTA) program in 1989 had demonstrated that cruising Propfans have undesirably high far-field noise level. This study attempts to investigate the source of these noise disturbances generated by a single-rotating propeller with supersonic tip speed. Various atmospheric and propagation effects responsible for shaping the far-field acoustics will be examined. A methodology for analyzing Propfan propagation noise is proposed based on the theory of geometrical acoustics. This approach allows acoustic rays to be traced into the far-field via Snell's law. Amplitude variation along the rays are evaluated by applying the Blokhintzev energy invariant principle. Additional terms are appended to the energy invariant principle to account for propagation and atmospheric effects. A FORTRAN program is written to analyze the rays emitted from a helicoidal source path traced by supersonic helical speed blade tip. Earlier papers by Myers and Farassat have demonstrated the presence of Mach envelopes that lead to formation of caustics. Linear geometric acoustics is used to ascertain the existence of these nonlinear caustics in the far-field. The importance of caustics to Propfan noise generation is then studied.

Sim, Ben WEL-C.↗

Traj_opt User's Guide

Trajectory optimization program Traj_opt was developed at Ames Research Center to help assess the potential benefits of ultrahigh temperature ceramic materials applied to reusable space vehicles with sharp noses and wing leading edges. Traj_opt loosely couples the Ames three-degrees-of-freedom trajectory package Traj (see NASA-TM-2004-212847) with the SNOPT optimization package (Stanford University Technical Report SOL 98-1). Traj_opt version January 22, 2003 is covered by this user guide. The program has been applied extensively to entry and ascent abort trajectory calculations for sharp and blunt crew transfer vehicles. The main optimization variables are control points for the angle of attack and bank angle time histories. No propulsion options are provided, but numerous objective functions may be specified and the nonlinear constraints implemented include a distributed surface heating constraint capability. Aero-capture calculations are also treated with an option to minimize orbital eccentricity at apoapsis. Traj_opt runs efficiently on a single processor, using forward or central differences for the gradient calculations. Results may be displayed conveniently with Gnuplot scripts. Control files recommended for five standard reentry and ascent abort trajectories are included along with detailed descriptions of the inputs and outputs.

Saunders, David A.↗

The prediction of pressure distributions on an arrow-wing configuration including the effect of camber, twist, and a wing fin

Wind tunnel tests of an arrow wing body configuration consisting of flat, twisted, and cambered twisted wings were conducted at Mach numbers from 0.40 to 2.50 to provide an experimental data base for comparison with theoretical methods. A variety of leading and trailing edge control surface deflections were included in these tests, and in addition, the cambered twisted wing was tested with an outboard vertical fin to determine its effect on wing and control surface loads. Theory experiment comparisons show that current state of the art linear and nonlinear attached flow methods were adequate at small angles of attack typical of cruise conditions. The incremental effects of outboard fin, wing twist, and wing camber are most accurately predicted by the advanced panel method PANAIR. Results of the advanced panel separated flow method, obtained with an early version of the program, show promise that accurate detailed pressure predictions may soon be possible for an aeroelasticity deformed wing at high angles of attack.

Bobbitt, P. J.↗

A rapid method for the computation of equilibrium chemical composition of air to 15000 K

A rapid computational method has been developed to determine the chemical composition of equilibrium air to 15000 K. Eleven chemically reacting species, i.e., O2, N2, O, NO, N, NO+, e-, N+, O+, Ar, and Ar+ are included. The method involves combining algebraically seven nonlinear equilibrium equations and four linear elemental mass balance and charge neutrality equations. Computational speeds for determining the equilibrium chemical composition are significantly faster than the often used free energy minimization procedure. Data are also included from which the thermodynamic properties of air can be computed. A listing of the computer program together with a set of sample results are included.

Prabhu, Ramadas K.↗

Symbolic modeling of flexible manipulators

An algorithm is presented to symbolically derive the full nonlinear dynamic equations of motion of multilink flexible manipulators. Lagrange's assumed modes method is the basis of new algorithm and adapted in a way suitable for symbolic manipulation by digital computers. It is applied to model a two-link flexible arm by means of a commercially available symbolic manipulation program. The advantages of the algorithm and simulation results are discussed.

Cetinkunt, Sabri↗

On 3-D inelastic analysis methods for hot section components

The objective of this program is to produce a series of new computer codes that permit more accurate and efficient three-dimensional inelastic structural analysis of combustor liners, turbine blades, and turbine vanes. Each code embodies a progression of mathematical models for increasingly comprehensive representation of the geometrical features, loading conditions, and forms of nonlinear material response that distinguish these three groups of hot section components.

Todd, E. S.↗

Three-dimensional inelastic analysis methods for hot section components

The objective of this program is to produce a series of new computer codes that permit more accurate and efficient three-dimensional inelastic structural analysis of combustor liners, turbine blades, and turbine vanes. Each code embodies a progression of mathematical models for increasingly comprehensive representation of the geometrical features, loading conditions, and forms of nonlinear material response that distinguish these three groups of hot section components.

Todd, E. S.↗

Development of an integrated BEM for hot fluid-structure interaction

A boundary element approach is developed for problems of hot fluid-structure interaction. The structure is idealized as a thermoelastic solid, and a time-domain, boundary-only integral formulation is described. Meanwhile, several approaches are discussed for the hot fluid. Integrals equations are developed for both compressible and incompressible thermoviscous flow. Due to the presence of nonlinear convective terms in the governing differential equations, domain discretization is generally required. However, with the introduction of reference velocities, volume modeling often can be confined to regions near obstacles and walls. All formulations are implemented for two-dimensional problems in GPO-BEST, a general purpose boundary element computer program. An overview of this numerical implementation is provided, along with several illustrative examples. The present fluid formulations are appropriate in the low to medium Reynolds number ranges; however, some enhancements required for higher speed simulations are noted.

Dargush, G. F.↗

A computer analysis program for interfacing thermal and structural codes

A software package has been developed to transfer three-dimensional transient thermal information accurately, efficiently, and automatically from a heat transfer analysis code to a structural analysis code. The code is called three-dimensional TRansfer ANalysis Code to Interface Thermal and Structural codes, or 3D TRANCITS. TRANCITS has the capability to couple finite difference and finite element heat transfer analysis codes to linear and nonlinear finite element structural analysis codes. TRANCITS currently supports the output of SINDA and MARC heat transfer codes directly. It will also format the thermal data output directly so that it is compatible with the input requirements of the NASTRAN and MARC structural analysis codes. Other thermal and structural codes can be interfaced using the transfer module with the neutral heat transfer input file and the neutral temperature output file. The transfer module can handle different elemental mesh densities for the heat transfer analysis and the structural analysis.

Thompson, R. L.↗

A computer analysis program for interfacing thermal and structural codes

A software package has been developed to transfer three-dimensional transient thermal information accurately, efficiently, and automatically from a heat transfer analysis code to a structural analysis code. The code is called three-dimensional TRansfer ANalysis Code to Interface Thermal and Structural codes, or 3D TRANCITS. TRANCITS has the capability to couple finite difference and finite element heat transfer analysis codes to linear and nonlinear finite element structural analysis codes. TRANCITS currently supports the output of SINDA and MARC heat transfer codes directly. It will also format the thermal data output directly so that it is compatible with the input requirements of the NASTRAN and MARC structural analysis codes. Other thermal and structural codes can be interfaced using the transfer module with the neutral heat transfer input file and the neutral temperature output file. The transfer module can handle different elemental mesh densities for the heat transfer analysis and the structural analysis.

Thompson, R. L.↗

Solution of steady-state, two-dimensional conservation laws by mathematical programming

A truly two-dimensional algorithm is created for solving the steady-state two-dimensional conservation-law problem. An overdetermined system of algebraic equations is obtained through discretization by finite-volume formulas. These equations are perturbed nonsingularly and are solved by an efficient geometrically oriented l(1) procedure. The basic algorithm and the theory for the linear case f(u) = u are presented, and computational results for the nonlinear case f(u) = sq u are also analyzed. It is noted that the l(1) procedure captures boundary shocks as well as oblige and zigzag interior shocks in bands that are one cell wide, and the solution values are accurate up to the edge of the shock.

Lavery, John E.↗