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.
Two of the more common procedures for analyzing the stability and forced response of equations with periodic coefficients are reviewed: the use of Floquet methods, and the use of multiblade coordinate and harmonic balance methods. The analysis procedures of these periodic coefficient systems are compared with those of the more familiar constant coefficient systems.
Two of the more common procedures for analyzing the stability and forced response of equations with periodic coefficients are reviewed: the use of Floquet methods, and the use of multiblade coordinate and harmonic balance methods. The analysis procedures of these periodic coefficient systems are compared with those of the more familiar constant coefficient systems. Previously announced in STAR as N82-23702
Lateral-stability flight tests were made over the Mach number range from 0.7 to 1.3 of models of three airplane configurations having 45deg sweptback wings. One model had a high wing; one, a low wing; and one, a high wing with cathedral. The models were otherwise identical. The lateral oscillations of the models resulting from intermittent yawing disturbances were interpreted in terms of full-scale airplane flying qualities and were further analyzed by the time-vector method to obtain values of the lateral stability derivatives. The effects of changes i n wing height on the static sideslip derivatives were fairly constant in the speed range investigated and agreed well with estimated values based on subsonic wind-tunnel tests. Effects of geometric dihedral on the rolling moment due to sideslip agreed well with theoretical and other experimental results and with a theoretical relation involving the damping in roll. The damping in roll, when compared with theoretical and other experimental results, shared good agreement at supersonic speeds but was somewhat higher at a Mach number of 1.0 and at subsonic speeds. The damping in yaw shared no large changes in the transonic region.
In another work, a method was presented, primarily conceived as an emergency backup system, that addressed the problem of a space capsule that needed to execute a safe atmospheric entry from an arbitrary initial attitude and angular rate in the absence of nominal control capability. The proposed concept permits the arrest of a tumbling motion, orientation to the heat shield forward position and the attainment of a ballistic roll rate of a rigid spacecraft with the use of control in one axis only. To show the feasibility of such concept, the technique of single input single output (SISO) feedback linearization using the Lie derivative method was employed and the problem was solved for different number of jets and for different configurations of the inertia matrix: the axisymmetric inertia matrix (I(sub xx) > I(sub yy) = I(sub zz)), a partially complete inertia matrix with I(sub xx) > I(sub yy) > I(sub zz), I(sub xz) not = 0 and a realistic complete inertia matrix with I(sub xx) > I(sub yy) > I)sub zz), I(sub ij) not= 0. The closed loop stability of the proposed non-linear control on the total angle of attack, Theta, was analyzed through the zero dynamics of the internal dynamics for the case where the inertia matrix is axisymmetric (I(sub xx) > I(sub yy) = I(sub zz)). This note focuses on the problem of the diagonal non-axisymmetric inertia matrix (I(sub xx) > I(sub yy) > I(sub zz)), which is half way between the axisymmetric and the partially complete inertia matrices. In this note, the control law for this type of inertia matrix will be determined and its closed-loop stability will be analyzed using the same methods that were used in the other work. In particular, it will be proven that the control system is stable in closed-loop when the actuators only provide a roll torque.
It is shown how the generalized characteristic locus method developed by MacFarlane (1972) can be used for the step-by-step design of multivariable closed-loop compensators. The design procedure attempts to strike a balance between four conflicting objectives: stability, integrity, interaction, and accuracy. The generalized locus method is illustrated by a jet engine example.
High-order Taylor series expansions can be used to model nonlinear dynamics at the cost of integrating a large set of variational equations to obtain high-order state-transition tensors (STTs). This paper presents an innovative technique for approximating the high-order STTs that reduces significantly the computational cost by retaining only the dominant secular terms. We propagate the low-order partial derivatives of Kepler’s equation, which only requires the integration of six additional equations to extend an n-th order approximation to order (n + 1). The approximation stems from the Lindstedt-Poincare procedure and exploits the stability properties of orbital motion. Since the method makes no dynamical assumptions, it can accommodate any source of orbital perturbations. We show how the approximation of the second-order STT significantly increases the accuracy of the linear method for uncertainty propagation with only a small computational overhead. Finally, we derive a high-order approximate extended Kalman filter that implements the proposed approximation of the STT and improves the performance of linear filters. Examples of application with different perturbation sources include the heliocentric orbit of an asteroid, an orbiter around Europa, and an Earth-orbiting satellite.
Relative comparisons of Loran-C frequency transmissions between the master station of Catanzaro (Simeri Crichi) and the X, Z slave stations of Estartit (Spain) and Lampedusa (Italy) are carrying out by the GG LORSTA monitor station of the Mediterranean Sea Loran-C chain. These comparisons are able to emphasize the relative and, under certain conditions, the absolute rate of the emitting standard frequencies of the slave stations and some relevant statistical properties of the Loran-C Method for frequency transmission and time synchronization. The stability of each Loran-C frequency standard transmission is subject to perturbations, more or less known, due to the propagation medium and other causes. Following the Allan (1966) method for data processing, the performance of the relative rate of frequency of the transmissions of the X, Z slave stations are described calculating the standard deviation of a set of N frequency measurements from its mean averaged during sampling times. This standard deviation is designated as the measure of the stability of the Loran-C frequency transmission.
The derivation of an approximate error characteristic equation describing the transient system error response is given, along with a procedure for selecting adaptive gain parameters so as to relate to the transient error response. A detailed example of the application and implementation of these methods for a space shuttle type vehicle is included. An extension of the characteristic equation technique is used to provide an estimate of the magnitude of the maximum system error and an estimate of the time of occurrence of this maximum after a plant parameter disturbance. Techniques for relaxing certain stability requirements and the conditions under which this can be done and still guarantee asymptotic stability of the system error are discussed. Such conditions are possible because the Lyapunov methods used in the stability derivation allow for overconstraining a problem in the process of insuring stability.
The finite difference formulation and method of solution is presented for a wide variety of fluid flow problems with associated heat transfer. Only a few direct results from these formulations are given as examples, since the book is intended primarily to serve a discussion of the techniques and as a starting point for further investigations; however, the formulations are sufficiently complete that a workable computer program may be written from them. In the appendixes a number of topics are discussed which are of interest with respect to the finite difference equations presented. These include a very rapid method for solving certain sets of linear algebraic equations, a discussion of numerical stability, the inherent error in flow rate for confined flow problems, and a method for obtaining high accuracy with a relatively small number of mesh points.
In this work, the feasibility of using massively parallel computation to study the response of ablative materials is investigated. Explicit and implicit finite difference methods are used on a massively parallel computer, the Thinking Machines CM-5. The governing equations are a set of nonlinear partial differential equations. The governing equations are developed for three sample problems: (1) transpiration cooling, (2) ablative composite plate, and (3) restrained thermal growth testing. The transpiration cooling problem is solved using a solution scheme based solely on the explicit finite difference method. The results are compared with available analytical steady-state through-thickness temperature and pressure distributions and good agreement between the numerical and analytical solutions is found. It is also found that a solution scheme based on the explicit finite difference method has the following advantages: incorporates complex physics easily, results in a simple algorithm, and is easily parallelizable. However, a solution scheme of this kind needs very small time steps to maintain stability. A solution scheme based on the implicit finite difference method has the advantage that it does not require very small times steps to maintain stability. However, this kind of solution scheme has the disadvantages that complex physics cannot be easily incorporated into the algorithm and that the solution scheme is difficult to parallelize. A hybrid solution scheme is then developed to combine the strengths of the explicit and implicit finite difference methods and minimize their weaknesses. This is achieved by identifying the critical time scale associated with the governing equations and applying the appropriate finite difference method according to this critical time scale. The hybrid solution scheme is then applied to the ablative composite plate and restrained thermal growth problems. The gas storage term is included in the explicit pressure calculation of both problems. Results from ablative composite plate problems are compared with previous numerical results which did not include the gas storage term. It is found that the through-thickness temperature distribution is not affected much by the gas storage term. However, the through-thickness pressure and stress distributions, and the extent of chemical reactions are different from the previous numerical results. Two types of chemical reaction models are used in the restrained thermal growth testing problem: (1) pressure-independent Arrhenius type rate equations and (2) pressure-dependent Arrhenius type rate equations. The numerical results are compared to experimental results and the pressure-dependent model is able to capture the trend better than the pressure-independent one. Finally, a performance study is done on the hybrid algorithm using the ablative composite plate problem. It is found that there is a good speedup of performance on the CM-5. For 32 CPU's, the speedup of performance is 20. The efficiency of the algorithm is found to be a function of the size and execution time of a given problem and the effective parallelization of the algorithm. It also seems that there is an optimum number of CPU's to use for a given problem.
Various tools for the aeroelastic stability and response analysis of rotor blades in hover and forward flight were developed and incorporated in a comprehensive package capable of performing aeroelastic tailoring of rotor blades in forward flight. The results indicate that substantial vibration reductions, of order 15-40%, in the vibratory hub shears can be achieved by relatively small modifications of the initial design. Furthermore the optimized blade can be up to 20% lighter than the original design. Accomplishments are reported for the following tasks: (1) finite element modeling of rotary-wing aeroelastic problems in hover and forward flight; (2) development of numerical methods for calculating the aeroelastic response and stability of rotor blades in forward fight; (3) formulation of the helicopter air resonance problem in hover with active controls; and (4) optimum design of rotor blades for vibration reduction in forward flight.
The stability and robustness of the controlled system obtained via the substructure control synthesis (SCS) method of Su et al. (1990) were examined using a six-bay truss model, and employing an LQG control design method to obtain controllers for two separate structures. It is found that the assembled controller provides a stability in this instance. A qualitative assessment of the stability robustness of the system with controller designed with the SCS method is provided by obtaining a controller using the complete truss model and comparing the robustness of the corresponding closed-loop systems.
The motivation to increase the step size with no degradation of numerical accuracy and stability has led to the discovery of particular members of the class of stiffly stable implicit linear multistep algorithms. Sufficient conditions for a consistent linear multistep method to be stiffly stable are given. These conditions involve properties of the stability mapping from the extended complex plane onto itself.
Performance degradation of multiple-input multiple-output (MIMO) control systems having limited actuators is often handled by augmenting the controller with an antiwindup mechanism, which attempts to maintain system performance when limits are encountered. The goals of this paper are: 1) to develop a method to analyze antiwindup systems to determine precisely what stability and performance degradation is incurred under limited conditions. It is shown that by reformulating limited actuator commands as resulting from multiplicative perturbations to the corresponding controller requests, mu-analysis tools can be utilized to obtain quantitative measures of stability and performance degradation. 2) To propose a linear, time invariant (LTI) criterion on which to base the antiwindup design. These analysis and design methods are illustrated through the evaluation of two competing antiwindup schemes augmenting the controller of a Short Take-Off and Vertical Landing (STOVL) aircraft in transition flight.
Performance degradation of multiple-input multiple-output (MIMO) control systems having limited actuators is often handled by augmenting the controller with an antiwindup mechanism, which attempts to maintain system performance when limits are encountered. The goals of this paper are: (1) To develop a method to analyze antiwindup systems to determine precisely what stability and performance degradation is incurred under limited conditions. It is shown that by reformulating limited actuator commands as resulting from multiplicative perturbations to the corresponding controller requests, mu-analysis tools can be utilized to obtain quantitative measures of stability and performance degradation. (2) To propose a linear, time invariant (LTI) criterion on which to base the antiwindup design. These analysis and design methods are illustrated through the evaluation of two competing antiwindup schemes augmenting the controller of a Short Take-Off and Vertical Landing (STOVL) aircraft in transition flight.
A structured singular value (mu) analysis method of computing flutter margins has robust stability of a linear aeroelastic model with uncertainty operators (Delta). Flight data is used to update the uncertainty operators to accurately account for errors in the computed model and the observed range of aircraft dynamics of the aircraft under test caused by time-varying aircraft parameters, nonlinearities, and flight anomalies, such as test nonrepeatability. This mu-based approach computes predict flutter margins that are worst case with respect to the modeling uncertainty for use in determining when the aircraft is approaching a flutter condition and defining an expanded safe flight envelope for the aircraft that is accepted with more confidence than traditional methods that do not update the analysis algorithm with flight data by introducing mu as a flutter margin parameter that presents several advantages over tracking damping trends as a measure of a tendency to instability from available flight data.
It seems desirable to have some simple method for calculating quickly and with sufficient accuracy: 1) the correct position of the center of gravity; 2) the requisite tail-group dimensions; 3) and the course of the wing and tail-group moments. In out deductions, we will first replace the biplane (disregarding the effect of stagger, decalage and induced drag) by an equivalent monoplane, whose dimensions and position in space can be approximately determined in a simple manner.
This report presents refinements in two areas of the initial design presented in the report entitled 'Conceptual Design of a Fleet of Autonomous Regolith Throwing Devices for Radiation Shielding of Lunar Habitats'. The first section presents an evaluation of the critical areas of the design and presents alternative solutions for these areas. The areas for design refinement are the traction required by the device and the stability of the device when throwing regolith. Several alternative methods are presented to solve these problems. First, the issue of required traction is covered. Next, the design is refined to provide a more stable device. The issue of stability is addressed both by presenting solutions for the configuration chosen for the computer simulation and by presenting two more device configurations. The next section presents the selected solutions. To prevent inadequate traction, the depth of dig-per-pass is reduced. A method combining a dynamic counterweight and an outrigger is chosen to provide a stable device.