Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear programming problem”

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 307 records · Page 17

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) for matching single-input single-output (SISO) transfer function models to empirical frequency response data is described. The program was written in MATLAB and has a graphical user interface (GUI). It is interactive in that the user manually builds the transfer function model using ``elementary factors'' (gain, delay, differentiators and integrators, and first- and second-order poles and zeros) and adjusts their values with sliders or entry fields. A nonlinear optimization can also be used to determine maximum-likelihood estimates of the transfer function parameters and their associated uncertainties. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and other system identification problems. FRFit is demonstrated using example problems, including the identification of aircraft transfer functions and rational function approximations of Theodorsen's function.

Frequency response↗

Deterministic Design Optimization of Structures in OpenMDAO Framework

Nonlinear programming algorithms play an important role in structural design optimization. Several such algorithms have been implemented in OpenMDAO framework developed at NASA Glenn Research Center (GRC). OpenMDAO is an open source engineering analysis framework, written in Python, for analyzing and solving Multi-Disciplinary Analysis and Optimization (MDAO) problems. It provides a number of solvers and optimizers, referred to as components and drivers, which users can leverage to build new tools and processes quickly and efficiently. Users may download, use, modify, and distribute the OpenMDAO software at no cost. This paper summarizes the process involved in analyzing and optimizing structural components by utilizing the framework s structural solvers and several gradient based optimizers along with a multi-objective genetic algorithm. For comparison purposes, the same structural components were analyzed and optimized using CometBoards, a NASA GRC developed code. The reliability and efficiency of the OpenMDAO framework was compared and reported in this report.

Coroneos, Rula M.↗

Deterministic Design Optimization of Structures in OpenMDAO Framework

Nonlinear programming algorithms play an important role in structural design optimization. Several such algorithms have been implemented in OpenMDAO framework developed at NASA Glenn Research Center (GRC). OpenMDAO is an open source engineering analysis framework, written in Python, for analyzing and solving Multi-Disciplinary Analysis and Optimization (MDAO) problems. It provides a number of solvers and optimizers, referred to as components and drivers, which users can leverage to build new tools and processes quickly and efficiently. Users may download, use, modify, and distribute the OpenMDAO software at no cost. This paper summarizes the process involved in analyzing and optimizing structural components by utilizing the framework s structural solvers and several gradient based optimizers along with a multi-objective genetic algorithm. For comparison purposes, the same structural components were analyzed and optimized using CometBoards, a NASA GRC developed code. The reliability and efficiency of the OpenMDAO framework was compared and reported in this report.

Coroneos, Rula M.↗

A direct method for synthesizing low-order optimal feedback control laws with application to flutter suppression

A direct method of synthesizing a low-order optimal feedback control law for a high order system is presented. A nonlinear programming algorithm is employed to search for the control law design variables that minimize a performance index defined by a weighted sum of mean square steady state responses and control inputs. The controller is shown to be equivalent to a partial state estimator. The method is applied to the problem of active flutter suppression. Numerical results are presented for a 20th order system representing an aeroelastic wind-tunnel wing model. Low-order controllers (fourth and sixth order) are compared with a full order (20th order) optimal controller and found to provide near optimal performance with adequate stability margins.

Mukhopadhyay, V.↗

Parallel processors and nonlinear structural dynamics algorithms and software

Techniques are discussed for the implementation and improvement of vectorization and concurrency in nonlinear explicit structural finite element codes. In explicit integration methods, the computation of the element internal force vector consumes the bulk of the computer time. The program can be efficiently vectorized by subdividing the elements into blocks and executing all computations in vector mode. The structuring of elements into blocks also provides a convenient way to implement concurrency by creating tasks which can be assigned to available processors for evaluation. The techniques were implemented in a 3-D nonlinear program with one-point quadrature shell elements. Concurrency and vectorization were first implemented in a single time step version of the program. Techniques were developed to minimize processor idle time and to select the optimal vector length. A comparison of run times between the program executed in scalar, serial mode and the fully vectorized code executed concurrently using eight processors shows speed-ups of over 25. Conjugate gradient methods for solving nonlinear algebraic equations are also readily adapted to a parallel environment. A new technique for improving convergence properties of conjugate gradients in nonlinear problems is developed in conjunction with other techniques such as diagonal scaling. A significant reduction in the number of iterations required for convergence is shown for a statically loaded rigid bar suspended by three equally spaced springs.

Belytschko, Ted↗

Parallel solver for trajectory optimization search directions

A key algorithmic element of a real-time trajectory optimization hardware/software implementation is presented, the search step solver. This is one piece of an algorithm whose overall goal is to make nonlinear trajectory optimization fast enough to provide real-time commands during guidance of a vehicle such as an aeromaneuvering orbiter or the National Aerospace Plane. Many methods of nonlinear programming require the solution of a quadratic program (QP) at each iteration to determine the search step. In the trajectory optimization case, the QP has a special dynamic programming structure. The algorithm exploits this special structure with a divide- and conquer type of parallel implementation. The algorithm solves a (p.N)-stage problem on N processors in O(p + log2 N) operations. The algorithm yields a factor of 8 speed-up over the fastest known serial algorithm when solving a 1024-stage test problem on 32 processors.

Psiaki, M. L.↗

Application of symbolic/numeric matrix solution techniques to the NASTRAN program

The matrix solving algorithm of any finite element algorithm is extremely important since solution of the matrix equations requires a large amount of elapse time due to null calculations and excessive input/output operations. An alternate method of solving the matrix equations is presented. A symbolic processing step followed by numeric solution yields the solution very rapidly and is especially useful for nonlinear problems.

Buturla, E. M.↗

Elevated temperature crack growth

The objective of the Elevated Temperature Crack Growth Program is to evaluate proposed nonlinear fracture mechanics methods for application to hot section components of aircraft gas turbine engines. Progress during the past year included linear-elastic fracture mechanics data reduction on nonlinear crack growth rate data on Alloy 718. The bulk of the analytical work centered on thermal gradient problems and proposed fracture mechanics parameters. Good correlation of thermal gradient experimental displacement data and finite element prediction was obtained.

Malik, S. N.↗

The pEst version 2.1 user's manual

This report is a user's manual for version 2.1 of pEst, a FORTRAN 77 computer program for interactive parameter estimation in nonlinear dynamic systems. The pEst program allows the user complete generality in definig the nonlinear equations of motion used in the analysis. The equations of motion are specified by a set of FORTRAN subroutines; a set of routines for a general aircraft model is supplied with the program and is described in the report. The report also briefly discusses the scope of the parameter estimation problem the program addresses. The report gives detailed explanations of the purpose and usage of all available program commands and a description of the computational algorithms used in the program.

Murray, James E.↗

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 quasi-Newton approach to optimization problems with probability density constraints

A quasi-Newton method is presented for minimizing a nonlinear function while constraining the variables to be nonnegative and sum to one. The nonnegativity constraints were eliminated by working with the squares of the variables and the resulting problem was solved using Tapia's general theory of quasi-Newton methods for constrained optimization. A user's guide for a computer program implementing this algorithm is provided.

Tapia, R. A.↗

Three dimensional inelastic finite element analysis of laminated composites

Formulations of the inelastic response of laminated composites to thermal and mechanical loading are used as the basis for development of the computer NALCOM (Nonlinear Analysis of Laminated Composites) computer program which uses a fully three dimensional isoparametric finite element with 24 nodes and 72 degrees of freedom. An incremental solution is performed with nonlinearities introduced as pseudoloads computed for initial strains. Equilibrium iteration may be performed at every step. Elastic and elastic-plastic response of boron/epoxy and graphite/epoxy graphite/epoxy and problems of curing 0/90 sub s Gr/Ep laminates with and without circular holes are analyzed. Mechanical loading of + or - 45sub s Gr/Ep laminates is modeled and symmetry conditions which exist in angle-ply laminates are discussed. Results are compared to experiments and other analytical models when possible. All models are seen to agree reasonably well with experimetnal results for off-axis tensile coupons. The laminate analyses show the three dimensional effects which are present near holes and free corners.

Griffin, O. H., Jr.↗

Using Hough harmonics to validate and assess nonlinear shallow-water models

The implementation of a technique for locating programming errors in shallow-water codes, establishing the correctness of the code, and assessing the performance of the numerical model under various flow conditions is described. The right-hand side of the differential equations is modified in such a way that the exact solution of the nonlinear initial-value problem is known, so that the truncation errors of the numerical scheme can be studied in detail. The exact solution is prescribed to be any linear combination of Hough harmonics which propagate in time according to their natural frequencies.

Dee, Dick P.↗

TROJID: A portable software package for upper-stage trajectory optimization

Performance optimization for upper-stage exoatmospheric vehicles often is performed within the framework of a full capability trajectory simulation package requiring either a large mainframe computer or powerful work-station. Since these software packages tend to include capabilities providing for high-fidelity boost and reentry simulations, the programs usually are quite large and not very portable. The program TROJID is an attempt to provide an environment for the optimization of upper-stage trajectories within a small package capable of being run on a standard desktop microcomputer. Utilizing a state-of-the-art nonlinear programming algorithm and a trajectory simulator implementing impulsive burns and an analytic coast phase propagator, TROJID is capable of producing trajectories for optimal multi-burn upper-stage orbit transfers. The package has been designed to allow full generality in definition of both the trajectory simulator and the parameter optimization problem.

Hammes, Steven M.↗

Communications, Signal Processing, and Telemetering Research Program Review

A Communications, Signal Processing, and Telemetering Research Program Review was held on February 23, 1999. Research conducted under the grant was presented and reviewed, for progress, and for possible technology transfers. The research reviewed was in the following areas: (1) Bandwidth-efficient Modulation and nonlinear equalization; (2) Investigation of an architecture for parallel signal processing applicable to communications problems; (3)Coded partial response over satellites; (4) synchronization at Low SNR; (5) Serial concatenated convolutional codes and some implementation issues on high rate turbo codes; (6) Flight experiments; (7) Real time doppler tracking; (8) Space protocol testing; (9) Lightweight optical communications without carrying a laser in space. The presentations are given by the graduate students who performed the research.

Source record↗

A Comparison of Trajectory Optimization Methods for the Impulsive Minimum Fuel Rendezvous Problem

In this paper we present a comparison of optimization approaches to the minimum fuel rendezvous problem. Both indirect and direct methods are compared for a variety of test cases. The indirect approach is based on primer vector theory. The direct approaches are implemented numerically and include Sequential Quadratic Programming (SQP), Quasi-Newton, Simplex, Genetic Algorithms, and Simulated Annealing. Each method is applied to a variety of test cases including, circular to circular coplanar orbits, LEO to GEO, and orbit phasing in highly elliptic orbits. We also compare different constrained optimization routines on complex orbit rendezvous problems with complicated, highly nonlinear constraints.

Hughes, Steven P.↗

Transmission of sound through nonuniform circular ducts with compressible mean flows

An acoustic theory is developed to determine the sound transmission and attenuation through an infinite, hard-walled or lined, circular duct carrying compressible, sheared, mean flows and having a variable cross section. The theory is applicable to large as well as small axial variations, as long as the mean flow does not separate. Although the theory is described for circular ducts, it is applicable to other duct configurations - annular, two dimensional, and rectangular. The theory is described for the linear problem, but the technique is general and has the advantage of being applicable to the nonlinear case as well as the linear case. The technique is based on solving for the envelopes of the quasi-parallel acoustic modes that exist in the duct instead of solving for the actual wave. A computer program was developed. The mean flow model consists of a one dimensional flow in the core and a quarter-sine profile in the boundary layer. Results are presented for the reflection and transmission coefficients in ducts with varying slopes and carrying different mean flows.

Nayfeh, A. H.↗

Use of global functions for improvement in efficiency of nonlinear analysis

A method for improving the efficiency of nonlinear structural analysis by the use of global displacement functions is presented. The computer programs include options to define the global functions as input or let the program automatically select and update these functions. The program was applied to a number of structures: (1) 'pear-shaped cylinder' in compression, (2) bending of a long cylinder, (3) spherical shell subjected to point force, (4) panel with initial imperfections, (5) cylinder with cutouts. The sample cases indicate the usefulness of the procedure in the solution of nonlinear structural shell problems by the finite element method. It is concluded that the use of global functions for extrapolation will lead to savings in computer time.

Almroth, B. O.↗