Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix equations”

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 163 records · Page 9

Development of three-dimensional spatial displays using a geographically based information system

Procedures necessary for the development of a generalized three-dimensional perspective software capability in support of graphic, topographic, and color mapping of Landsat data are reviewed. The NASA Earth Resources Laboratory developed the procedures in order to facilitate the processing and analysis of disparate, geographically oriented base maps from aircraft and satellite sensors. Perspective displays are obtained through a translation of the space-viewed object to a vantage point coordinate system, followed by a rotation through two angles for alignment along the vantage line of sight, and finally a perspective transformation to yield two-dimensional displays with no hidden lines. Matrix equations for the transformations are reviewed, including scaling, and block diagrams are provided of the data and perspective software systems. The classification data plane may be mapped onto a topographic elevation data plane.

Junkin, B. G.↗

Wave scattering from a periodic dielectric surface for a general angle of incidence

Electromagnetic waves scattered from a periodic dielectric and perfectly conducting surface are studied for a general angle of incidence. It is shown that the one-dimensional corrugated surface can be solved by using two scalar functions: the components of the electric and magnetic fields along the row direction of the surface, and appropriate boundary conditions to obtain simple matrix equations. Results are compared to the case where the incident angle wave vector is perpendicular to the row direction. Numerical results demonstrate that energy conservation and reciprocity are obeyed for scattering by sinusoidal surfaces for the general case, which checks the consistency of the formalism.

Chuang, S. L.↗

ACCESS computer program for the synthesis of large structural systems

The structural-synthesis computer code ACCESS-3, originally limited to the optimization of relatively simple problems involving truss, membrane, and shear panel elements, has been restructured to make it applicable to structures with a larger number of degrees-of-freedom, design variables, and constraints. The modifications made allow the use of six degrees-of-freedom per node and permit the out-of-core solution of large matrix equations encountered at the analysis stage. Also, finite elements which permit pure bending as well as combined bending and membrane (axial) behavior have been added to satisfy the need for a more comprehensive set of elements. Like its predecessor, the new version of the program combines the approximation concepts with primal or dual formulations to provide an efficient optimization tool.

Fleury, C.↗

Unconditionally stable implicit-explicit algorithms for coupled thermal stress waves

An unconditionally stable implicit-explicit method is proposed for the analysis of transient coupled thermal stress waves, and the computer-implementation aspects of the method are discussed. In the method proposed here, the mechanical displacement and temperature are used as independent variables. The resulting coupled finite element matrix equations are symmetric.

Liu, W. K.↗

Mixed time integration schemes for transient conduction forced-convection analysis

A partition procedure for forced-convection conduction transient problems is presented. Mixed time partitions are defined wherein coupled conduction force-matrix equations are discretized using an implicit integration method, followed by derivation of a mixed time integration technique. Explicit-implicit and explicit-explicit partitions are performed for a stability analysis for transient conditions, e.g., those found in an actively air-cooled engine and airframe structure.

Liu, W. K.↗

Finite elements and the method of conjugate gradients on a concurrent processor

An algorithm for the iterative solution of finite element problems on a concurrent processor is presented. The method of conjugate gradients is used to solve the system of matrix equations, which is distributed among the processors of a MIMD computer according to an element-based spatial decomposition. This algorithm is implemented in a two-dimensional elastostatics program on the Caltech Hypercube concurrent processor. The results of tests on up to 32 processors show nearly linear concurrent speedup, with efficiencies over 90% for sufficiently large problems.

Lyzenga, G. A.↗

Developement of an Optimum Interpolation Analysis Method for the CYBER 205

A state-of-the-art technique to assimilate the diverse observational database obtained during FGGE, and thus create initial conditions for numerical forecasts is described. The GLA optimum interpolation (OI) analysis method analyzes pressure, winds, and temperature at sea level, mixing ratio at six mandatory pressure levels up to 300 mb, and heights and winds at twelve levels up to 50 mb. Conversion to the CYBER 205 required a major re-write of the Amdahl OI code to take advantage of the CYBER vector processing capabilities. Structured programming methods were used to write the programs and this has resulted in a modular, understandable code. Among the contributors to the increased speed of the CYBER code are a vectorized covariance-calculation routine, an extremely fast matrix equation solver, and an innovative data search and sort technique.

Nestler, M. S.↗

Direct computation of optimal control of forced linear system

It is known that the optimal control of a forced linear system may be reduced to that of tracking the system without forces. The solution of the tracking problem is available via the costate variables method. This procedure is computationally expensive for large order systems. It requires solution of matrix Riccati equation and two final value problems. An alternate approach is outlined for the direct computation of the optimal control. Instead of Riccati equation, a matrix Volterra integral must be solved. For this purpose two computational schemes are described, and an illustrative example is given. The results compare favorably with the classical solution. This alternative approach may be especially useful for the control of large space structure where large order models are required.

Utku, S.↗

The self-consistent dynamic pole tide in global oceans

The dynamic pole tide is characterized in a self-consistent manner by means of introducing a single nondifferential matrix equation compatible with the Liouville equation, modelling the ocean as global and of uniform depth. The deviations of the theory from the realistic ocean, associated with the nonglobality of the latter, are also given consideration, with an inference that in realistic oceans long-period modes of resonances would be increasingly likely to exist. The analysis of the nature of the pole tide and its effects on the Chandler wobble indicate that departures of the pole tide from the equilibrium may indeed be minimal.

Dickman, S. R.↗

PIFCGT: A PIF autopilot design program for general aviation aircraft

This report documents the PIFCGT computer program. In FORTRAN, PIFCGT is a computer design aid for determing Proportional-Integral-Filter (PIF) control laws for aircraft autopilots implemented with a Command Generator Tracker (CGT). The program uses Linear-Quadratic-Regulator synthesis algorithms to determine feedback gains, and includes software to solve the feedforward matrix equation which is useful in determining the command generator tracker feedforward gains. The program accepts aerodynamic stability derivatives and computes the corresponding aerodynamic linear model. The nine autopilot modes that can be designed include four maneuver modes (ROLL SEL, PITCH SEL, HDG SEL, ALT SEL), four final approach models (APR GS, APR LOCI, APR LOCR, APR LOCP), and a BETA HOLD mode. The program has been compiled and executed on a CDC computer.

Broussard, J. R.↗

Finite elements and the method of conjugate gradients on a concurrent processor

An algorithm for the iterative solution of finite element problems on a concurrent processor is presented. The method of conjugate gradients is used to solve the system of matrix equations, which is distributed among the processors of a MIMD computer according to an element-based spatial decomposition. This algorithm is implemented in a two-dimensional elastostatics program on the Caltech Hypercube concurrent processor. The results of tests on up to 32 processors show nearly linear concurrent speedup, with efficiencies over 90 percent for sufficiently large problems.

Lyzenga, G. A.↗

Polynomial compensation, inversion, and approximation of discrete time linear systems

The least-squares transformation of a discrete-time multivariable linear system into a desired one by convolving the first with a polynomial system yields optimal polynomial solutions to the problems of system compensation, inversion, and approximation. The polynomial coefficients are obtained from the solution to a so-called normal linear matrix equation, whose coefficients are shown to be the weighting patterns of certain linear systems. These, in turn, can be used in the recursive solution of the normal equation.

Baram, Yoram↗

Sparse matrix methods research using the CSM testbed software system

Research is described on sparse matrix techniques for the Computational Structural Mechanics (CSM) Testbed. The primary objective was to compare the performance of state-of-the-art techniques for solving sparse systems with those that are currently available in the CSM Testbed. Thus, one of the first tasks was to become familiar with the structure of the testbed, and to install some or all of the SPARSPAK package in the testbed. A suite of subroutines to extract from the data base the relevant structural and numerical information about the matrix equations was written, and all the demonstration problems distributed with the testbed were successfully solved. These codes were documented, and performance studies comparing the SPARSPAK technology to the methods currently in the testbed were completed. In addition, some preliminary studies were done comparing some recently developed out-of-core techniques with the performance of the testbed processor INV.

Chu, Eleanor↗

Model reduction for discrete bilinear systems

A model reduction method for discrete bilinear systems is developed which matches q sets of Volterra and covariance parameters. These parameters are shown to represent both deterministic and stochastic attributes of the discrete bilinear system. A reduced order model which matches these q sets of parameters is defined to be a q-Volterra covariance equivalent realization (q-Volterra COVER). An algorithm is presented which constructs a class of q-Volterra COVERs parameterized by solutions to a Hermitian, quadratic, matrix equation. The algorithm is applied to a bilinear model of a robot manipulator.

King, A. M.↗

Modal techniques for analyzing airplane dynamics

A series of techniques are presented for analyzing airplane control and output characteristics. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled, modal coordinates. Techniques for analyzing the impacts of system inputs on the fundamental modes of motion and the appearance of these modes in the system outputs are explained. Scaled modal matrices are used to compute steady-state control inputs which optimize the steady-state response of selected system outputs. Graphics which promote quick understanding of the analysis are presented to display the resulting vectors and matrices. The defined analysis techniques are applied to an example airplane model, illustrating the insight which can be acquired using the described modal techniques.

Arbuckle, P. Douglas↗

Characteristic-Wave Approach Complements Modal Analysis

Aspects of estimation of unmodeled dynamics discussed. Report discusses solution of nonhomogeneous governing matrix equation for dynamics of short vibrational pulses propagating as characteristic waves in large structure. Applied to analyze response, to repeated pulses, of beam clamped at one end and free at other. Shows all qualitative characteristics occuring under arbitrary periodic excitations of beam and those of quasi-periodic excitations, in as much as such excitations obtained by linear superpositions of periodic excitations.

Zak, Michail↗

H2, fixed architecture, control design for large scale systems

The H2, fixed architecture, control problem is a classic linear quadratic Gaussian (LQG) problem whose solution is constrained to be a linear time invariant compensator with a decentralized processing structure. The compensator can be made of p independent subcontrollers, each of which has a fixed order and connects selected sensors to selected actuators. The H2, fixed architecture, control problem allows the design of simplified feedback systems needed to control large scale systems. Its solution becomes more complicated, however, as more constraints are introduced. This work derives the necessary conditions for optimality for the problem and studies their properties. It is found that the filter and control problems couple when the architecture constraints are introduced, and that the different subcontrollers must be coordinated in order to achieve global system performance. The problem requires the simultaneous solution of highly coupled matrix equations. The use of homotopy is investigated as a numerical tool, and its convergence properties studied. It is found that the general constrained problem may have multiple stabilizing solutions, and that these solutions may be local minima or saddle points for the quadratic cost. The nature of the solution is not invariant when the parameters of the system are changed. Bifurcations occur, and a solution may continuously transform into a nonstabilizing compensator. Using a modified homotopy procedure, fixed architecture compensators are derived for models of large flexible structures to help understand the properties of the constrained solutions and compare them to the corresponding unconstrained ones.

Mercadal, Mathieu↗

A controls engineering approach for analyzing airplane input-output characteristics

An engineering approach for analyzing airplane control and output characteristics is presented. State-space matrix equations describing the linear perturbation dynamics are transformed from physical coordinates into scaled coordinates. The scaling is accomplished by applying various transformations to the system to employ prior engineering knowledge of the airplane physics. Two different analysis techniques are then explained. Modal analysis techniques calculate the influence of each system input on each fundamental mode of motion and the distribution of each mode among the system outputs. The optimal steady state response technique computes the blending of steady state control inputs that optimize the steady state response of selected system outputs. Analysis of an example airplane model is presented to demonstrate the described engineering approach.

Arbuckle, P. Douglas↗