SEARCH · Search NASA
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.
The numerical operations involved in a currently used optimization technique are discussed and analyzed with special attention to the numerical accuracy. Alternative methods for deriving linear system transfer functions, finding the relationships between the transfer function coefficients and the design parameters, and solving a matrix equation are presented for more accurate and cost effective solutions.
For a vehicle amenable to idealization as a set of n + 1 rigid bodies interconnected by n line hinges (implying tree topology), a set of minimum-dimension matrix equations of motion is recorded in terms of variables including the relative rotations of contiguous bodies. When a substructure or portion of the vehicle is subjected to large gross deformations but is restricted to 'small' strains or local deformations, the n-body model is (for sufficiently large n) assumed to experience 'small' relative rotations of contiguous bodies within the substructure, and the equations of motion are linearized in these variables. The application of coordinate transformations and truncations to the linearized subset of the system equations of motion then produces a hybrid-coordinate formulation of the problem involving both discrete and distributed coordinates.
A wide range of eigenvalue problems encountered in practical structural engineering analyses is defined, in which the structures are assumed to be discretized by any suitable technique such as the finite-element method. A review of the usual numerical procedures for the solution of such eigenvalue problems is presented and is followed by an extensive account of recently developed eigenproblem solution procedures. Particular emphasis is placed on the new numerical algorithms and associated computer programs based on the Sturm sequence method. Eigenvalue algorithms developed for efficient solution of natural frequency and buckling problems of structures are presented, as well as some eigenvalue procedures formulated in connection with the solution of quadratic matrix equations associated with free vibration analysis of structures. A new algorithm is described for natural frequency analysis of damped structural systems.
A procedure is described for synthesis of optimal aircraft control systems by application of the concepts of optimal control theory to time-invariant linear systems with quadratic performance criteria. Essential in this synthesis procedure is the solution of the Riccati matrix equation which results in a constant linear feedback control law for an output regulator which maintains a plant in an equilibrium in the presence of impulse disturbances. An algorithm is derived for designing maneuverable output regulators with selected state variables for feedback.
Linear optimal quadratic control theory is applied to a low wing-loading STOL aircraft for ride quality and flight path following. Design criteria include minimum rms response to wind turbulence and desired transient response characteristics. Design techniques include proper choosing of design versus evaluation models, choosing appropriate performance index responses, and use of classical evaluation techniques. Results are obtained through a combination of frequency response shaping and gust observation. Effects of control rate and authority saturation are examined with a new rapid calculation of random input describing functions. Parameter sensitivity is also evaluated using a Liapunov type matrix equation.
This paper explores the concept of finite dynamic elements involving higher order dynamic correction terms in the associated stiffness and mass matrices. Such matrices are then developed for a rectangular prestressed membrane element. Next, efficient analysis techniques for the eigenproblem solution of the resulting quadratic matrix equations are described in detail. These are followed by suitable numerical examples which indicate that employment of such dynamic elements in conjunction with an efficient quadratic matric solution technique will result in a most significant economy in the free vibration analysis of structures.
A technique utilizing the finite element displacement method is developed for the static analysis of structures subjected to non-deterministic thermal loading in which the material properties, assumed isotropic, are temperature dependent. Matrix equations are developed for the first two statistical moments of the displacements using a third order series expansion for the displacements in terms of the random temperatures. Sample problems are included to demonstrate the range of applicability of the third order series solutions. These solutions are compared with results from Monte Carlo analyses and also, for some problems, with solutions obtained by numerically integrating equations for the statistical properties of the displacements. In general, it is shown that the effect of temperature dependent material properties can have a significant effect on the covariances of the displacements.
A procedure is presented for synthesizing time-invariant control logic to cause the outputs of a linear plant to track the outputs of an unforced (or randomly forced) linear dynamic system. The control logic uses feed-forward of the reference system state variables and feedback of the plant state variables. The feed-forward gains are obtained from the solution of a linear algebraic matrix equation of the Liapunov type. The feedback gains are the usual regulator gains, determined to stabilize (or augment the stability of) the plant, possibly including integral control. The method is applied here to the design of control logic for a second-order servomechanism to follow a linearly increasing (ramp) signal, an unstable third-order system with two controls to track two separate ramp signals, and a sixth-order system with two controls to track a constant signal and an exponentially decreasing signal (aircraft landing-flare or glide-slope-capture with constant velocity).
The thickness (T) of the lunar elastic lithosphere at the time (3.6 to 3.8 billion years ago) of the earliest preserved basalt flows in circular mare basins can be estimated by inverting the observed locations of extensional tectonic features in and surrounding the maria. In performing the inversion, the lithosphere is modeled as an elastic shell with a liquid interior, and the basalt load for each mare is approximated by a set of concentric cylinders. To permit solving the forward problem of placing radial limits on the positions of the rilles around a given mare, an additional parameter F, the ratio of the radial stress at the radial limits to the maximum radial stress, is introduced. T and F are chosen to give the best weighted-squares fit of the radial limits to the observations, and are used as the initial values in a linearized matrix inversion to check the resolution and estimate errors. The application of the procedure to three maria with prominent extensional features, Humorum, Orientale, and Serenitatis, gives values of T from about 40 + or - 10 to 50 + or - 10 km, and in each case the linearized matrix equation has an exact inverse.
An analytical design technique for an active flutter-suppression and gust-alleviation control system is presented. It is based on a rational approximation of the unsteady aerodynamic loads in the entire Laplace domain, which yields matrix equations of motion with constant coefficients. Some existing rational approximation schemes are reviewed, and a new technique which yields a minimal number of augmented states for a desired accuracy is presented. The state-space aeroelastic model is used to design a constant gain, partial-feedback control system, which simultaneously assures stability and optimizes any desired combination of gust response parameters throughout the entire flight envelope.
The properties of spectral methods are surveyed and their extension to solve problems in complex geometries is developed. A new iteration procedure is introduced to solve efficiently the full matrix equations resulting from spectral approximations to nonconstant coefficient boundary-value problems in complex geometries. It is shown that the work required to solve these spectral equations exceeds that of solving the lowest-order finite-difference approximation to the same problem by only O(N log N).
A method of fluid-structure coupling which provides symmetrical matrix equations of standard form solved by existing finite element computer programs is presented. The method postulates that the uncoupled vibration modes of the fluid or the structure be calculated before the coupled analysis. A numerical solution of vibration modes in an axisymmetric container demonstrated that a static approximation to higher order fluid modes can improve the accuracy of dynamic response computations using modal methods.
Evaluation of two-sided complex integrals is often required when analyzing linear systems to determine signal variances resulting from stochastic inputs and system noise bandwidths. Algebraic solutions of integrals in a closed matrix equation form, using coefficients of the numerator and denominator polynomials, are presented. The closed forms provide the possibility of obtaining some insight into parameter sensitivity in addition to greatly reducing the computational complexity required by the normal method of evaluation by residues.
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.
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.
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.
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.
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.