Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear dynamical 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 1,045 records · Page 58

A comparison of controller designs for an experimental flexible structure

Control systems design and hardware testing are addressed for an experimental structure that displays the characteristics of a typical flexible spacecraft. The results of designing and implementing various control design methodologies are described. The design methodologies under investigation include linear quadratic Gaussian control, static and dynamic dissipative controls, and H-infinity optimal control. Among the three controllers considered, it is shown, through computer simulation and laboratory experiments on the evolutionary structure, that the dynamic dissipative controller gave the best results in terms of vibration suppression and robustness with respect to modeling errors.

Lim, K. B.↗

Large-amplitude finite element flutter analysis of composite panels in hypersonic flow

Consideration is given to a finite-element approach for determining the nonlinear flutter characteristics of 3D thin laminated composite panels using the full third-order-piston transverse-loading aerodynamic theory. The unsteady hypersonic aerodynamic theory and the von Karman large-deflection-plate theory are used to formulate the aeroelasticity problem. Nonlinear flutter analyses are performed to assess the influence of the higher-order aerodynamic theory on the structure's limit-cycle amplitude and the dynamic pressure of the flow velocity. A solution procedure is presented to solve the nonlinear panel flutter and large-amplitude free-vibration finite-element equations. Nonlinear flutter analyses are performed for different boundary support-conditions and for various system parameters. Linear finite-element flutter for isotropic and composite panels and large-amplitude isotropic panel flutter results are compared with existing classical solutions. The large-amplitude panel flutter results using the full third-order piston aerodynamic theory are presented to assess the influence of the nonlinear aerodynamic theory.

Gray, Carl E., Jr.↗

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗

Control Surface Interaction Effects of the Active Aeroelastic Wing Wind Tunnel Model

This paper presents results from testing the Active Aeroelastic Wing wind tunnel model in NASA Langley s Transonic Dynamics Tunnel. The wind tunnel test provided an opportunity to study aeroelastic system behavior under combined control surface deflections, testing for control surface interaction effects. Control surface interactions were observed in both static control surface actuation testing and dynamic control surface oscillation testing. The primary method of evaluating interactions was examination of the goodness of the linear superposition assumptions. Responses produced by independently actuating single control surfaces were combined and compared with those produced by simultaneously actuating and oscillating multiple control surfaces. Adjustments to the data were required to isolate the control surface influences. Using dynamic data, the task increases, as both the amplitude and phase have to be considered in the data corrections. The goodness of static linear superposition was examined and analysis of variance was used to evaluate significant factors influencing that goodness. The dynamic data showed interaction effects in both the aerodynamic measurements and the structural measurements.

Heeg, Jennifer↗

Revolutionary Vertical Lift Technology Lift + Cruise Concept Vehicle Powertrain Dynamic Model

The Lift + Cruise (L+C) is one of several Revolutionary Vertical Lift Technology (RVLT) concept aircraft identified by NASA to investigate Urban Air Mobility (UAM) requirements. This six-passenger vehicle has a conventional tube-and-wing configuration, but with eight lift rotors (four on each wing) and a pusher rotor on the tail. This arrangement allows it to take off and land vertically as well as fly like an airplane. This report presents a novel dynamic model of the L+C powertrain built using the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS) and the Electrical Modeling and Thermal Analysis Toolbox (EMTAT), both of which are NASA-developed toolboxes, and together are capable of modeling hybrid-electric powertrain architectures. This research effort models the complex behavior of the powertrain at the mechanical timescale (15 to 20 ms). The model consists of a turboshaft engine, generator, direct current bus, energy storage management system, gearboxes, and electric propulsor systems coupled together to create a dynamic model that mimics the operation of the multidomain powertrain. Two control algorithms, designed to equalize the power consumed to the power produced, are compared to assess their practicality: (1) independent proportional-integral (PI) control and (2) linear quadratic integral regulator (LQIR) control. Simulation results highlight the key features present in the model and demonstrate its operation under the two control schemes.

revolutionary vertical lift technology↗

Computational fluid dynamics on a massively parallel computer

A finite difference code was implemented for the compressible Navier-Stokes equations on the Connection Machine, a massively parallel computer. The code is based on the ARC2D/ARC3D program and uses the implicit factored algorithm of Beam and Warming. The codes uses odd-even elimination to solve linear systems. Timings and computation rates are given for the code, and a comparison is made with a Cray XMP.

Jespersen, Dennis C.↗

A methodology for formulating a minimal uncertainty model for robust control system design and analysis

In the design and analysis of robust control systems for uncertain plants, the technique of formulating what is termed an M-delta model has become widely accepted and applied in the robust control literature. The M represents the transfer function matrix M(s) of the nominal system, and delta represents an uncertainty matrix acting on M(s). The uncertainty can arise from various sources, such as structured uncertainty from parameter variations or multiple unstructured uncertainties from unmodeled dynamics and other neglected phenomena. In general, delta is a block diagonal matrix, and for real parameter variations the diagonal elements are real. As stated in the literature, this structure can always be formed for any linear interconnection of inputs, outputs, transfer functions, parameter variations, and perturbations. However, very little of the literature addresses methods for obtaining this structure, and none of this literature addresses a general methodology for obtaining a minimal M-delta model for a wide class of uncertainty. Since have a delta matrix of minimum order would improve the efficiency of structured singular value (or multivariable stability margin) computations, a method of obtaining a minimal M-delta model would be useful. A generalized method of obtaining a minimal M-delta structure for systems with real parameter variations is given.

Belcastro, Christine M.↗

How Well do Multisatellite Products Capture the Space-Time Dynamics of Precipitation? Part II: Building an Error Model Through Spectral System Identification

Satellite precipitation products, as all quantitative estimates, come with some inherent degree of uncertainty. To associate a quantitative value of the uncertainty to each individual estimate, error modeling is necessary. Most of the error models proposed so far compute the uncertainty as a function of precipitation intensity only, and only at one specific spatio-temporal scale. We propose a spectral error model which accounts for the neighboring space-time dynamics of precipitation into the uncertainty quantification. Systematic distortions of the precipitation signal and random errors are characterized distinctively in every frequency-wavenumber band in the Fourier domain, to accurately characterize error across scales. The systematic distortions are represented as a deterministic space-time linear filtering term. The random errors are represented as a non-stationary additive noise. The spectral error model is applied to the IMERG multi satellite precipitation product and its parameters are estimated empirically through a system identification approach using the GV-MRMS gauge-radar measurements as reference (“truth”) over the eastern United States. The filtering term is found to be essentially low-pass. While traditional error models attribute most of the error variance to random errors, it is found here that the systematic filtering term explains 48% of the error variance at the native resolution of IMERG. This fact confirms that, at high resolution, filtering effects in satellite precipitation products cannot be ignored, and that the error cannot be represented as a purely random additive or multiplicative term. An important consequence is that precipitation estimates derived from totally different sources shall not be expected to automatically have statistically independent errors.

Precipitation↗

A Short-Range Distance Sensor with Exceptional Linearity

A sensor has been demonstrated that can measure distance over a total range of about 300 microns to an accuracy of about 0.1 nm (resolution of about 0.01 nm). This represents an exceptionally large dynamic range of operation - over 1,000,000. The sensor is optical in nature, and requires the attachment of a mirror to the object whose distance is being measured. This work resulted from actively developing a white light interferometric system to be used to measure the depths of defects in the Space Shuttle Orbiter windows. The concept was then applied to measuring distance. The concept later expanded to include spectrometer calibration. In summary, broadband (i.e., white) light is launched into a Michelson interferometer, one mirror of which is fixed and one of which is attached to the object whose distance is to be measured. The light emerging from the interferometer has traveled one of two distances: either the distance to the fixed mirror and back, or the distance to the moving mirror and back. These two light beams mix and produce an interference pattern where some wavelengths interfere constructively and some destructively. Sending this light into a spectrometer allows this interference pattern to be analyzed, yielding the net distance difference between the two paths. The unique feature of this distance sensor is its ability to measure accurately distance over a dynamic range of more than one million, the ratio of its range (about 300 microns) to its accuracy (about 0.1 nanometer). Such a large linear operating range is rare and arises here because both amplitude and phase-matching algorithms contribute to the performance. The sensor is limited by the need to attach a mirror of some kind to the object being tracked, and by the fairly small total range, but the exceptional dynamic range should make it of interest.

Simmons, Steven↗

A simple, mass balance model of carbon flow in a controlled ecological life support system

Internal cycling of chemical elements is a fundamental aspect of a Controlled Ecological Life Support System (CELSS). Mathematical models are useful tools for evaluating fluxes and reservoirs of elements associated with potential CELSS configurations. A simple mass balance model of carbon flow in CELSS was developed based on data from the CELSS Breadboard project at Kennedy Space Center. All carbon reservoirs and fluxes were calculated based on steady state conditions and modelled using linear, donor-controlled transfer coefficients. The linear expression of photosynthetic flux was replaced with Michaelis-Menten kinetics based on dynamical analysis of the model which found that the latter produced more adequate model output. Sensitivity analysis of the model indicated that accurate determination of the maximum rate of gross primary production is critical to the development of an accurate model of carbon flow. Atmospheric carbon dioxide was particularly sensitive to changes in photosynthetic rate. The small reservoir of CO2 relative to large CO2 fluxes increases the potential for volatility in CO2 concentration. Feedback control mechanisms regulating CO2 concentration will probably be necessary in a CELSS to reduce this system instability.

Garland, Jay L.↗

Gaussian step-pressure loading of rigid viscoplastic plates

The response of a thin, rigid viscoplastic plate subjected to a spatially axisymmetric Gaussian step pressure impulse loading was studied analytically. A Gaussian pressure distribution in excess of the collapse load was applied to the plate, held constant for a length of time, and then suddenly removed. The plate deforms with monotonically increasing deflections until the dynamic energy is completely dissipated in plastic work. The simply supported plate of uniform thickness obeys the von Mises yield criterion and a generalized constitutive equation for rigid viscoplastic materials. For the small deflection bending response of the plate, the governing system of equations is essentially nonlinear. Transverse shear stress is neglected in the yield condition and rotary inertia in the equations of dynamic equilibrium. A proportional loading technique, known to give excellent approximations of the exact solution for the uniform load case, was used to linearize the problem and to obtain the analytical solutions in the form of eigenvalue expansions. The effects of load concentration, of an order of magnitude change in the viscosity of the plate material, and of load duration were examined while holding the total impulse constant.

Hayduk, R. J.↗

Atmospheric turbulence parameters for modeling wind turbine dynamics

A model which can be used to predict the response of wind turbines to atmospheric turbulence is given. The model was developed using linearized aerodynamics for a three-bladed rotor and accounts for three turbulent velocity components as well as velocity gradients across the rotor disk. Typical response power spectral densities are shown. The system response depends critically on three wind and turbulence parameters, and models are presented to predict desired response statistics. An equation error method, which can be used to estimate the required parameters from field data, is also presented.

Holley, W. E.↗

Periodic control of the individual-blade-control helicopter rotor

This paper describes the results of an investigation into methods of controller design for linear periodic systems utilizing an extension of modern control methods. Trends present in the selection of various cost functions are outlined, and closed-loop controller results are demonstrated for two cases: first, on an analog computer simulation of the rigid out of plane flapping dynamics of a single rotor blade, and second, on a 4 ft diameter single-bladed model helicopter rotor in the MIT 5 x 7 subsonic wind tunnel, both for various high levels of advance ratio. It is shown that modal control using the IBC concept is possible over a large range of advance ratios with only a modest amount of computational power required.

Mckillip, R. M., Jr.↗

The significance of error dynamics in model-following for flight control design

The role of the system error dynamics in model-following control systems is discussed, along with the use of handling quality specifications, actuation bandwidth constraints, stability, and closed-loop performance requirements in flight control design. The model-following problem is formulated using both direct state-space and linear, quadratic, optimization techniques. The results are then demonstrated using several examples involving a generic forward-swept-wing vehicle and a conventional flight vehicle with large parameter uncertainty in order to illustrate the trade-off in closed-loop performance and control law complexity.

Schmidt, David K.↗

Low order dynamical models for mixing layers

The goal of this project was to investigate the feasibility of developing useful low order dynamical systems to allow one to determine cheaply the effect of various initial conditions on the evolution of mixing layers. A mathematically rigorous approach for a time evolving base flow as found in a mixing layer seems to be the singular value decomposition of the linear operator transforming some initial perturbations into the solution at a given later time. Unfortunately, this is very costly, if at all feasible, and only deals with the linear evolution of a disturbance. A different approach based on weakly nonlinear analysis was attempted, but numerous difficulties were encountered. Several formulations were proposed which may be ill-posed. The reasons for this possible failure of the attempted approach is related to the fundamental limitations of weakly nonlinear analysis. The clarification of these issues is attempted in this report.

Riahi, D. N.↗

Adaptive attitude control and momentum management for large-angle spacecraft maneuvers

The fully coupled equations of motion are systematically linearized around an equilibrium point of a gravity gradient stabilized spacecraft, controlled by momentum exchange devices. These equations are then used for attitude control system design of an early Space Station Freedom flight configuration, demonstrating the errors caused by the improper approximation of the spacecraft dynamics. A full state feedback controller, incorporating gain-scheduled adaptation of the attitude gains, is developed for use during spacecraft on-orbit assembly or operations characterized by significant mass properties variations. The feasibility of the gain adaptation is demonstrated via a Space Station Freedom assembly sequence case study. The attitude controller stability robustness and transient performance during gain adaptation appear satisfactory.

Parlos, Alexander G.↗

Reply to the comment by Y. Ricard, R. Sabadini, and G. Spada

We welcome this opportunity which Drs. Ricard, Sabadini and Spada (RSS) have offered us to elaborate on some points concerning the connection between convective circulations and rotational dynamics. The relationship of the linear viscoelastic normal-mode theory to other approaches, such as viscous and the non-linear viscoelastic treatments, can be found in Moser et al. (1993). We will address the issues raised by RSS by first going over the reference system. This will be followed by a discussion concerning long-term polar wander and relative angular momentum. We will also go over the issues concerning the off-diagonal terms in the total moment of inertia and polar wander.

Moser, Jiri↗

Analog Techniques for Measuring the Frequency Response of Linear Physical Systems Excited by Frequency-Sweep Inputs

Data-reduction methods using general-purpose analog computer equipment and compatible testing techniques for determining the frequency response of linear physical systems are examined. The techniques considered may be classed as steady state or transient depending on the method of excitation. The relative merits of periodic, slow sweep, and transient (rapid sweep) forcing functions are discussed and applications are given that relate to dynamic-response tests of aeroelastic systems. Two frequency-sweep-input methods are considered in detail. one case the sweep rate is sufficiently slow that the response is approximately the same as that for steady-state conditions. With this input the frequency response can be evaluated and displayed in real time while the test is in progress. Errors due to treating sweep data as steady state can be eliminated, when desired, by reanalyzing tape-recorded time histories of the input and output as transient rather than as periodic data. In the second method the frequency-response function is deter- mined from the system's transient response to a very rapid sweep input. The purpose of frequency sweep in this case is to provide sufficient harmonic content in the input to overcome noise while keeping the test time as short as possible. tests and limited flight-test data presented herein, it appears that a transient-type rapid-sweep forcing function offers a considerable saving in test time while preserving the accuracy possible with steady-state sinusoidal inputs.

Wilmer H. Reed III↗