Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 757 records · Page 42

Influence of unbalance on the nonlinear dynamical response and stability of flexible rotor-bearing systems

In this paper, some of the effects of unbalance on the nonlinear response and stability of flexible rotor-bearing systems is presented from both a theoretical and experimental standpoint. In a linear system, operating above its stability threshold, the amplitude of motion grows exponentially with time and the orbits become unbounded. In an actual system, this is not necessarily the case. The actual amplitudes of motion may be bounded due to various nonlinear effects in the system. These nonlinear effects cause limit cycles of motion. Nonlinear effects are inherent in fluid film bearings and seals. Other contributors to nonlinear effects are shafts, couplings and foundations. In addition to affecting the threshold of stability, the nonlinear effects can cause jump phenomena to occur at not only the critical speeds, but also at stability onset or restabilization speeds.

Gunter, E. J.↗

A high-speed implementation of the random decrement algorithm

The algorithm is useful for measuring net system damping levels in stochastic processes and for the development of equivalent linearized system response models. The algorithm works by summing together all subrecords which occur after predefined threshold level is crossed. The random decrement signature is normally developed by scanning stored data and adding subrecords together. The high speed implementation of the random decrement algorithm exploits the digital character of sampled data and uses fixed record lengths of 2(n) samples to greatly speed up the process. The contributions to the random decrement signature of each data point was calculated only once and in the same sequence as the data were taken. A hardware implementation of the algorithm using random logic is diagrammed and the process is shown to be limited only by the record size and the threshold crossing frequency of the sampled data. With a hardware cycle time of 200 ns and 1024 point signature, a threshold crossing frequency of 5000 Hertz can be processed and a stably averaged signature presented in real time.

Kiraly, L. J.↗

Random field estimation approach to multibody dynamics

The use of recursive techniques similar to random field models to factor and invert the inertia matrix of a multibody system is discussed. An equivalence is established between the composite multibody system inertia matrix and the covariance of the output of a described linear system model. Conditional mean estimation and sequential estimation problems are solved along with problems of filtering and smoothing. Formulas are developed to compute the covariance of several relevant quantities. The foregoing results are used to obtain the inverse of the composite multibody system inertia in closed form.

Rodriguez, G.↗

Nonlinear input-output systems

Necessary and sufficient conditions that the nonlinear system dot-x = f(x) + ug(x) and y = h(x) be locally feedback equivalent to the controllable linear system dot-xi = A xi + bv and y = C xi having linear output are found. Only the single input and single output case is considered, however, the results generalize to multi-input and multi-output systems.

Hunt, L. R.↗

Mid-frequency Band Dynamics of Large Space Structures

High and low intensity dynamic environments experienced by a spacecraft during launch and on-orbit operations, respectively, induce structural loads and motions, which are difficult to reliably predict. Structural dynamics in low- and mid-frequency bands are sensitive to component interface uncertainty and non-linearity as evidenced in laboratory testing and flight operations. Analytical tools for prediction of linear system response are not necessarily adequate for reliable prediction of mid-frequency band dynamics and analysis of measured laboratory and flight data. A new MATLAB toolbox, designed to address the key challenges of mid-frequency band dynamics, is introduced in this paper. Finite-element models of major subassemblies are defined following rational frequency-wavelength guidelines. For computational efficiency, these subassemblies are described as linear, component mode models. The complete structural system model is composed of component mode subassemblies and linear or non-linear joint descriptions. Computation and display of structural dynamic responses are accomplished employing well-established, stable numerical methods, modern signal processing procedures and descriptive graphical tools. Parametric sensitivity and Monte-Carlo based system identification tools are used to reconcile models with experimental data and investigate the effects of uncertainties. Models and dynamic responses are exported for employment in applications, such as detailed structural integrity and mechanical-optical-control performance analyses.

dynamics↗

Application of modern control theory to the design of optimum aircraft controllers

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.

Power, L. J.↗

Guidance of Nonlinear Nonminimum-Phase Dynamic Systems

The research work has advanced the inversion-based guidance theory for: systems with non-hyperbolic internal dynamics; systems with parameter jumps; and systems where a redesign of the output trajectory is desired. A technique to achieve output tracking for nonminimum phase linear systems with non-hyperbolic and near non-hyperbolic internal dynamics was developed. This approach integrated stable inversion techniques, that achieve exact-tracking, with approximation techniques, that modify the internal dynamics to achieve desirable performance. Such modification of the internal dynamics was used (a) to remove non-hyperbolicity which is an obstruction to applying stable inversion techniques and (b) to reduce large preactuation times needed to apply stable inversion for near non-hyperbolic cases. The method was applied to an example helicopter hover control problem with near non-hyperbolic internal dynamics for illustrating the trade-off between exact tracking and reduction of preactuation time. Future work will extend these results to guidance of nonlinear non-hyperbolic systems. The exact output tracking problem for systems with parameter jumps was considered. Necessary and sufficient conditions were derived for the elimination of switching-introduced output transient. While previous works had studied this problem by developing a regulator that maintains exact tracking through parameter jumps (switches), such techniques are, however, only applicable to minimum-phase systems. In contrast, our approach is also applicable to nonminimum-phase systems and leads to bounded but possibly non-causal solutions. In addition, for the case when the reference trajectories are generated by an exosystem, we developed an exact-tracking controller which could be written in a feedback form. As in standard regulator theory, we also obtained a linear map from the states of the exosystem to the desired system state, which was defined via a matrix differential equation.

Devasia, Santosh↗

A low-rank solver for the stochastic unsteady Navier–Stokes problem

Here we study a low-rank iterative solver for the unsteady Navier–Stokes equations for incompressible flows with a stochastic viscosity. The equations are discretized using the stochastic Galerkin method, and we consider an all-at-once formulation where the algebraic systems at all the time steps are collected and solved simultaneously. The problem is linearized with Picard’s method. To efficiently solve the linear systems at each step, we use low-rank tensor representations within the Krylov subspace method, which leads to significant reductions in storage requirements and computational costs. Combined with effective mean-based preconditioners and the idea of inexact solve, we show that only a small number of linear iterations are needed at each Picard step. The proposed algorithm is tested with a model of flow in a two-dimensional symmetric step domain with different settings to demonstrate the computational efficiency.

97 MATHEMATICS AND COMPUTING↗

Decentralized control of large-scale systems: Fixed modes, sensitivity and parametric robustness

The problem of stabilization/pole placement under structural constraints of large scale linear systems is discussed. The existence of a solution to this problem is expressed in terms of fixed modes. The aim is to provide a bibliographic survey of the available results concerning the fixed modes (characterization, elimination, control structure selection to avoid them, control design in their absence) and to present the author's contribution to this problem which can be summarized by the use of the mode sensitivity concept to detect or to avoid them, the use of vibrational control to stabilize them, and the addition of parametric robustness considerations to design an optimal decentralized robust control.

Tarras, A.↗

Investigation, development and application of optimal output feedback theory. Vol. 4: Measures of eigenvalue/eigenvector sensitivity to system parameters and unmodeled dynamics

Some measures of eigenvalue and eigenvector sensitivity applicable to both continuous and discrete linear systems are developed and investigated. An infinite series representation is developed for the eigenvalues and eigenvectors of a system. The coefficients of the series are coupled, but can be obtained recursively using a nonlinear coupled vector difference equation. A new sensitivity measure is developed by considering the effects of unmodeled dynamics. It is shown that the sensitivity is high when any unmodeled eigenvalue is near a modeled eigenvalue. Using a simple example where the sensor dynamics have been neglected, it is shown that high feedback gains produce high eigenvalue/eigenvector sensitivity. The smallest singular value of the return difference is shown not to reflect eigenvalue sensitivity since it increases with the feedback gains. Using an upper bound obtained from the infinite series, a procedure to evaluate whether the sensitivity to parameter variations is within given acceptable bounds is developed and demonstrated by an example.

Halyo, Nesim↗

Identification of aerodynamic indicial functions using flight data

It is pointed out that the use of indicial function representation provides a model superior to the aerodynamic derivative model. Specific derivatives can be approximated from the indicial models. The model can also be used to compute equivalent stability and control parameters not usually available from flight data. It is shown that derivatives regarding the angle-of-attack and the side slip angle can be derived directly from the indicial functions without any identifiability problem. Attention is given to the pitch moment coefficient, linear indicial function representation, the identification problem for the pitch moment equation, the identifiability of linear systems, parametric representations of the indicial functions, an identification technique, angle-of-attack and pitch rate dynamics in the pitch plane, multivariate linear models, nonlinear aerodynamic indicial functions, measurement system accuracy, and poststall and spin-entry data from a scaled research vehicle.

Gupta, N. K.↗

A new theory for multistep discretizations of stiff ordinary differential equations: Stability with large step sizes

A large set of variable coefficient linear systems of ordinary differential equations which possess two different time scales, a slow one and a fast one is considered. A small parameter epsilon characterizes the stiffness of these systems. A system of o.d.e.s. in this set is approximated by a general class of multistep discretizations which includes both one-leg and linear multistep methods. Sufficient conditions are determined under which each solution of a multistep method is uniformly bounded, with a bound which is independent of the stiffness of the system of o.d.e.s., when the step size resolves the slow time scale, but not the fast one. This property is called stability with large step sizes. The theory presented lets one compare properties of one-leg methods and linear multistep methods when they approximate variable coefficient systems of stiff o.d.e.s. In particular, it is shown that one-leg methods have better stability properties with large step sizes than their linear multistep counter parts. The theory also allows one to relate the concept of D-stability to the usual notions of stability and stability domains and to the propagation of errors for multistep methods which use large step sizes.

Majda, G.↗

Eigenvalue assignment strategies in rotor systems

The work done to establish the control and direction of effective eigenvalue excursions of lightly damped, speed dependent rotor systems using passive control is discussed. Both second order and sixth order bi-axis, quasi-linear, speed dependent generic models were investigated. In every case a single, bi-directional control bearing was used in a passive feedback stabilization loop to resist modal destabilization above the rotor critical speed. Assuming incomplete state measurement, sub-optimal control strategies were used to define the preferred location of the control bearing, the most effective measurement locations, and the best set of control gains to extend the speed range of stable operation. Speed dependent control gains were found by Powell's method to maximize the minimum modal damping ratio for the speed dependent linear model. An increase of 300 percent in stable speed operation was obtained for the sixth order linear system using passive control. Simulations were run to examine the effectiveness of the linear control law on nonlinear rotor models with bearing deadband. The maximum level of control effort (force) required by the control bearing to stabilize the rotor at speeds above the critical was determined for the models with bearing deadband.

Youngblood, J. N.↗

Juqbox.jl

The Juqbox.jl package implements functionality for solving the quantum optimal control problem for realizing logical gates in closed quantum systems. The dynamics of the quantum system is modeled by Schroedinger's equation, which takes to form of a linear system of ordinary differential equations (ODE). Juqbox.jl solves this ODE by numerical time stepping and applies a gradient-based optimization technique to determine control pulses for driving an initial state to a final state, according to the desired logical gate transformation. To evaluate the gradient, Juqbox.jl applies the ``first discretize, then optimize'' approach based on a discrete adjoint time stepping technique. The actual optimization is performed by the open source Ipopt libaray. Juqbox.jl is written in the Julia programming language which, among many other features, provides a convenient interface to the Ipopt library.

PETERSSON, NILSA.↗

Aeroelastic stability of complete rotors with application to a teetering rotor in forward flight

The derivation of a set of nonlinear coupled flap-lag-torsion equations of motion for moderately large deflections of an elastic, two bladed teetering helicopter rotor in forward flight is concisely outlined. The following degrees of freedom are included in the mathematical model: rigid body flapping, rigid body lead lag, elastic bending in flap and lead-lag blade root torsion and shaft torsion. Quasi-steady aerodynamic loads are considered and the effects of reversed flow are included. The aeroelastic stability of the complete rotor is investigated using a linearized system of equations of motion. The equilibrium position about which the equations are linearized is obtained by considering the trim state of the helicopter, in true or simulated forward flight conditions. The sensitivity of the aeroelastic stability boundaries to interblade structural and mechanical coupling is illustrated by comparing the complete rotor stability boundaries with those obtained from a single blade analysis for a number of hover and forward flight cases.

Shamie, J.↗

Aeroelastic stability of complete rotors with application to a teetering rotor in forward flight

The derivation of a set of non-linear coupled flap-lag-torsion equations of motion for moderately large deflections of an elastic, two-bladed teetering helicopter rotor in forward flight is concisely outlined. The following degrees of freedom are included in the mathematical model: rigid body flapping, rigid body lead-lag, elastic bending in flap and lead-lag, blade root torsion, and shaft torsion. Quasi-steady aerodynamic loads are considered and the effects of reversed flow are included. The aeroelastic stability of the complete rotor is investigated by using a linearized system of equations of motion. The equilibrium position about which the equations are linearized is obtained by considering the trim state of the helicopter, in true or simulated forward flight conditions. The sensitivity of the aeroelastic stability boundaries to interblade structural and mechanical coupling is illustrated by comparing the complete rotor stability boundaries with those obtained from a single blade analysis for a number of hover and forward flight cases.

Shamie, J.↗

Control of nonlinear time-varying systems

Necessary and sufficient conditions are given for a time-varying nonlinear system of specific form to be transformed into a time-invariant controllable linear system. Since the present work will be in a neighborhood of the origin, it is unnecessary to name specific sets and it is assumed that all assumptions, conditions and results hold in an open set in the appropriate Euclidean space that contains the origin. This theory can be combined with the global inverse function theorems to produce global results.

Hunt, L. R.↗

Development and evaluation of a generic active helicopter vibration controller

A computerized generic active controller is developed, which alleviates helicopter vibration by closed-loop implementation of higher harmonic control (HHC). In the system, the higher harmonic blade pitch is input through a standard helicopter swashplate; for a four-blade helicopter rotor the 4/rev vibration in the rotorcraft is minimized by inducing cyclic pitch motions at 3, 4, and 5/rev in the rotating system. The controller employs the deterministic, cautious, and dual control approaches and two linear system models (local and global), as well as several methods of limiting control. Based on model testing, performed at moderate to high values of forward velocity and rotor thrust, reductions in the rotor test apparatus vibration from 75 to 95 percent are predicted, with HHC pitch amplitudes of less than one degree. Good performance is also noted for short-duration maneuvers.

Davis, M. W.↗