Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equation systems”

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 631 records · Page 35

Closed-form solutions for linear regulator-design of mechanical systems including optimal weighting matrix selection

This paper addresses the restriction of Linear Quadratic Regulator (LQR) solutions to the algebraic Riccati Equation to design spaces which can be implemented as passive structural members and/or dampers. A general closed-form solution to the optimal free-decay control problem is presented which is tailored for structural-mechanical systems. The solution includes, as subsets, special cases such as the Rayleigh Dissipation Function and total energy. Weighting matrix selection is a constrained choice among several parameters to obtain desired physical relationships. The closed-form solution is also applicable to active control design for systems where perfect, collocated actuator-sensor pairs exist. Some examples of simple spring mass systems are shown to illustrate key points.

Hanks, Brantley R.↗

Similar solutions for the compressible laminar boundary layer with heat transfer and pressure gradient

Stewartson's transformation is applied to the laminar compressible boundary-layer equations and the requirement of similarity is introduced, resulting in a set of ordinary nonlinear differential equations previously quoted by Stewartson, but unsolved. The requirements of the system are Prandtl number of 1.0, linear viscosity-temperature relation across the boundary layer, an isothermal surface, and the particular distributions of free-stream velocity consistent with similar solutions. This system admits axial pressure gradients of arbitrary magnitude, heat flux normal to the surface, and arbitrary Mach numbers. The system of differential equations is transformed to integral system, with the velocity ratio as the independent variable. For this system, solutions are found by digital computation for pressure gradients varying from that causing separation to the infinitely favorable gradient and for wall temperatures from absolute zero to twice the free-stream stagnation temperature. Some solutions for separated flows are also presented.

Cohen, Clarence B↗

Solving sparse triangular linear systems on parallel computers

This paper describes and compares three parallel algorithms for solving sparse triangular systems of equations. These methods involve some preprocessing overhead and are primarily of interest in solving many systems with the same coefficient matrix. The first approach is to use a fixed blocksize and form the inverse of the diagonal blocks. The second approach is to use a variable blocksize and reorder the unknowns so that the diagonal blocks are diagonal matrices. The latter technique is called level scheduling because of how it is represented in the adjacency graph, and both row-wise and jagged diagonal storage for the off-diagonal blocks are considered. These techniques are analyzed for general parallel computers and experiments are presented for the eight-processor Alliant FX/8.

Anderson, Edward↗

The dynamics and control of large flexible space structures, 2. Part A: Shape and orientation control using point actuators

The equations of planar motion for a flexible beam in orbit which includes the effects of gravity gradient torques and control torques from point actuators located along the beam was developed. Two classes of theorems are applied to the linearized form of these equations to establish necessary conditions for controlability for preselected actuator configurations. The feedback gains are selected: (1) based on the decoupling of the original coordinates and to obtain proper damping, and (2) by applying the linear regulator problem to the individual model coordinates separately. The linear control laws obtained using both techniques were evaluated by numerical integration of the nonlinear system equations. Numerical examples considering pitch and various number of modes with different combination of actuator numbers and locations are presented. The independent model control concept used earlier with a discretized model of the thin beam in orbit was reviewed for the case where the number of actuators is less than the number of modes. Results indicate that although the system is controllable it is not stable about the nominal (local vertical) orientation when the control is based on modal decoupling. An alternate control law not based on modal decoupling ensures stability of all the modes.

Bainum, P. M.↗

Linear dynamic coupling in geared rotor systems

The effects of high frequency oscillations caused by the gear mesh, on components of a geared system that can be modeled as rigid discs are analyzed using linear dynamic coupling terms. The coupled, nonlinear equations of motion for a disc attached to a rotating shaft are presented. The results of a trial problem analysis show that the inclusion of the linear dynamic coupling terms can produce significant changes in the predicted response of geared rotor systems, and that the produced sideband responses are greater than the unbalanced response. The method is useful in designing gear drives for heavy-lift helicopters, industrial speed reducers, naval propulsion systems, and heavy off-road equipment.

David, J. W.↗

Optimal inputs for system identification

Identification criteria are presented for linear dynamic systems with and without process noise. With process noise, the state equations are replaced by the Kalman filter equations. If the identification performance index is expanded in a Taylor's series with respect to the parameters to be identified, then maximizing the weighting factor of the quadratic term with respect to the inputs will insure that an identification algorithm will converge more rapidly and to a more accurate result than with nonoptimal inputs. The expectation of this weighting factor is the Fisher information matrix, and its inverse is a lower bound for the covariance of the parameters. Direct and indirect methods of calculating the information matrix are presented for systems with and without process noise.

Reid, D. B.↗

Farfield viscous effects in nonlinear noise propagation

Discussion of the method of parametric differentiation in application to predictions of farfield noise propagation in both lossless and dissipative media. It is shown that, in the lossless medium, the governing equation, transformed to parameter space, reduces to a wave equation in the farfield. In the dissipative medium, the system of nonlinear partial differential equations, transformed to parameter space, reduces to a linear partial differential equation of the propagating type which contains a third-order derivative as well as the wave operator.

Harris, W. L., Sr.↗

A decision and control multi-axis pilot model based on an urgency for action concept

An all digital multiple axis, multiple input, multiple output pilot model flies five different tasks in aircraft of different classes while encountering turbulence represented by the Dryden spectral model. Six degree of freedom linearized aircraft equations of motion are used together with a generalized stability augmentation system. The pilot model consists of two parts: (1) a decision maker, and (2) a control action implementer. This decision maker selects the critical variable and predicts the error at action implementation. The control action implementer decides the magnitude of control to be applied and applies it. This system has been successfully compared with hybrid/analog man in the loop simulations of the F-5, A-7, T-33, and 707 aircraft, thus validating the all digital simulation and the decision and control pilot model.

Pollard, J. J.↗

A block iterative LU solver for weakly coupled linear systems

A hybrid technique, called the block iterative LU solver, is proposed for solving the linear equations resulting from a finite element numerical analysis of certain fluid dynamics problems where the equations are weakly coupled between distinct sets of variables. Either the block Jacobi iterative method or the block Gauss-Seidel iterative solver is combined with LU decomposition.

Cooke, C. H.↗

Algorithms for linear-systems control

Control-theory design package aids design of linear-quadratic-Gaussian (LOG) controllers and optimal filters. It is applicable to systems that can be modeled by linear time-invariant differential or difference equations.

Armstrong, E. S.↗

Stochastic optimal attitude control of spacecraft with movable appendages

This paper deals with the application of linear optimal control and filtering theory to control a spinning spacecraft with movable telescoping appendages. The equations of motion are linearized about the desired final state. A feedback control system is designed to maintain this final state, with plant noise and measurement noise present in the system. Analytic results are obtained for special cases and numerical results are presented for the general case.

Sellappan, R.↗

Reanalysis and design in structural dynamics

A unified efficient formulation for the static and dynamic reanalysis of locally modified structures is presented. The reanalysis problem is to find the structural responses when some of the element properties are adjusted as a result of design modifications or when additional structures are appended to the original system. The reanalysis problem is formulated as a problem of much lower order than the original system. This is achieved by utilizing the linearity property of the structure using the pseudoload concept together with the solution of the original system. The modifications to the structure are treated as displacement-dependent pseudoloads of the system. By expressing the modified system response as linear combinations of the response of the original system and a term depending on the pseudoload, a reduced set of response equations can be obtained. In the static and sinusoidal steady state analysis, this leads to a set of linear algebraic equations. For a free vibration analysis, this results in a reduced eigenvalue problem. Using these reduced equations, a design problem can be solved; that is, the magnitude of a specified local modification for given dynamic response characteristics can be calculated. Several examples are presented to illustrate the general formulations.

Wang, B. P.↗

Solution of 3-dimensional time-dependent viscous flows. Part 1: Investigation of candidate algoriths

There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow to aid in the design of helicopter rotors. Numerical algorithms are examined to determine their overall suitability for the efficient and routine solution of an appropriate system of partial differential equations. It is concluded that a consistently split time linearized block implicit scheme using either quintic B spline collocation or the generalized operator compact implicit approach to generate a fourth order accurate algorithm is particularly well suited for use on the present problem. High cell Reynolds number behavior leads to favoring the generalized operator compact implicit approach over the quintic B spline collocation method.

Weinberg, B. C.↗

Optimal slewing maneuvers for flexible spacecraft using a closed form solution for the linear tracking problem

The problem of maneuvering a flexible spacecraft through a large angle while requiring the feedback control system to track a desired output state is considered. The desired output state is assumed to be provided by the corresponding open-loop maneuvering solution for the linear system model. Closed form solutions are provided for the Riccati and prefilter equations defining the optimal control. Example maneuvers are presented where control-rate penalties have been included in the performance index, in order to smooth both the open- and closed-loop control commands.

Turner, J. D.↗

On the relation between 'mixing length' and 'direct interaction approximation' theories of turbulence

The capabilities of direct interaction approximations (DIA) and weak-coupling approximations (WCA) in modeling turbulence in plasma physics studies are examined. DIA and WCA formalisms are reviewed and a DIA expression for homogeneous incompressible hydrodynamics is simplified to derive an energy spectrum in an arbitrary wavenumber space. A solution is defined for the expression when applied to long wavelength phenomena generated by an instability. The equation is further generalized to yield a transport equation for modeling the spectrum generated by homogeneous turbulence in a linearly unstable system undergoing nonlinear interaction. Finally, the technique is applied to two-dimensional turbulence in a low pressure, weakly ionized plasma embedded in a homogeneous magnetic field and experiencing instabilities caused by electron density and electric potential gradients, such as found in the ionospheric E region.

Sudan, R. N.↗

A computational method for the determination of the response of a linear system

A number of issues are discussed relating to the analysis of a linear damped oscillator equation with a forcing term where the right side of the equation, g(t), is known only at fixed, equal time intervals. (A particular class of such equations has application to civil earthquake engineering.) It is noted that the equation is not a differential equation because the condition of uniqueness does not obtain. A general computational method is presented for calculating x(t), based on an unconditionally stable finite difference technique. The work generalizes the results given in the recent paper of Ly (1984).

Mickens, R. E.↗