Errors Associated with Spinning-up and Thrusting Symmetric Rigid Bodies
Errors associated with spinning-up and thrusting symmetric rigid bodies
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.
Errors associated with spinning-up and thrusting symmetric rigid bodies
Rigid body motion implicit representation in curved finite elements
In cases where the equations of motion of a structure having rigid-body freedom are cast in state-variable form, generalized state rigid-body modes may be needed. It is possible to find a linearly-independent set of generalized vectors which transform an n x n matrix into the almost-diagonal Jordan form. Attention is presently given to equations governing these generalized eigenvectors, together with illustrative examples of the damped and undamped structure cases.
The rigid body modes of the PAH-2 'Tiger' helicopter were identified using the Eigensystem Realization Algorithm (ERA). This work complements ground vibration tests performed using DLR's traditional phase resonance technique and the ISSPA (Identification of Structural System Parameters) method. Rigid body modal parameters are important for ground resonance prediction. Time-domain data for ERA were obtained by inverse Fourier transformation of frequency response functions measured with stepped-sine excitation. Mode purity (based on the Phase Resonance Criterion) was generally equal to or greater than corresponding results obtained in the ground vibration tests. All identified natural frequencies and mode shapes correlate well with corresponding ground vibration test results. The modal identification approach discussed in this report has become increasingly attractive in recent years due to the steadily declining cost and increased performance of scientific computers. As illustrated in this application, modern time-domain methods can be successfully applied to data acquired using DLR's existing test equipment. Some suggestions are made for future applications of time domain modal identification in this manner.
The present paper extends a previous analysis, in Farassat (1974), of the acoustic field of rigid bodies. A restriction on that analysis can be removed by replacing the moving body by a stationary, virtual radiating body whose shape depends on the observer time and location and on the shape and motion of the rigid body. The possible singularity of the solution is considered.
Equations of motion are derived for use in simulating a spacecraft or other complex electromechanical system amenable to idealization as a set of hinge-connected rigid bodies of tree topology, with rigid axisymmetric rotors and nonrigid appendages attached to each rigid body in the set. In conjunction with a previously published report on finite-element appendage vibration equations, this report provides a complete minimum-dimension formulation suitable for generic programming for digital computer numerical integration.
Equations of motion are derived for use in simulating a spacecraft or other complex electromechanical system amenable to idealization as a set of hinge-connected rigid bodies of tree topology, with rigid axisymmetric rotors and nonrigid appendages attached to each rigid body in the set. In conjunction with a previously published companion paper on finite-element appendage vibration equations, this paper provides a complete minimum-dimension formulation suitable for generic programming for digital computer numerical integration.
Mode coupling effect on response of rigid body to random excitation analyzed, using two degree of freedom model
Space structures, such as the Space Station solar arrays, must be extremely light-weight, flexible structures. Accurate prediction of the natural frequencies and mode shapes is essential for determining the structural adequacy of components, and designing a controls system. The tension pre-load in the 'blanket' of photovoltaic solar collectors, and the free/free boundary conditions of a structure in space, causes serious reservations on the use of standard finite element techniques of solution. In particular, a phenomenon known as 'grounding', or false stiffening, of the stiffness matrix occurs during rigid body rotation. The authors have previously shown that the grounding phenomenon is caused by a lack of rigid body rotational capability, and is typical in beam geometric stiffness matrices formulated by others, including those which contain higher order effects. The cause of the problem was identified as the force imbalance inherent in the formulations. In this paper, the authors develop a beam geometric stiffness matrix for a directed force problem, and show that the resultant global stiffness matrix contains complete rigid body mode capabilities, and performs very well in the diagonalization methodology customarily used in dynamic analysis.
Relative motion of two spheroidal homogeneous rigid bodies of isolated binary system, obtaining potential energy
Existence and stability of equilibrium orientations of rigid body with gyrostat in nonuniform gravitational field
An investigation is made of the rotational motion of a rigid body orbiting the earth, under the influence of the geogravity-gradient torque. The attitude dynamics are formulated as perturbations from a nominal case. The perturbed equations are solved asymptotically, by the multiple scales technique. The independent variable, time, is extended into a space of higher dimension by means of new scales, fast, slow, etc. Integration of the equations is carried out separately in the new variables. The rapid and slow aspects of the attitude dynamics are systematically separated, resulting in a more efficient computer implementation and enhanced physical insight. The theory is applied to predict the attitude dynamics of an asymmetric rigid body satellite. A comparison is made of the maximum errors as the step size increases as predicted by the multiple scales solution and direct numerical integration. An improvement in computational speed of an order of magnitude is demonstrated.
A procedure is described for evaluating and subtracting the contribution of the rigid body motion from the general displacement of a Large Space Structure. The shape distortions are thus exhibited and their root mean square calculated. It is well known that a rigid body motion is composed of a rotation and a translation. To evaluate the rotation matrix M, use is made of the fact that unitary matrices can be expressed as M = (I-B)/(I+B) where B is computed using the coordinates, before and after the displacement, of three rigidly attached points. The translation vector is then deduced from the same coordinates.
Leimanis (1965) has published a monograph, which contains Boedewadt's solution of Euler's equations of motion for a symmetric rigid body subject to a time-independent, self-excitement in a body-fixed direction. In addition, the monograph provided also Boedewadt's solution for the corresponding angles of rotation. It is pointed out that the solution of Euler's equations of motion provides a poor approximation to the problem of nearly symmetric rigid bodies in cases in which a high degree of accuracy (less than 1 percent error) is required. A much more accurate approximate solution has been provided by Longuski (1980). The result given by Leimanis in the case of the Eulerian angles is incorrect for arbitrary constant torques. The approximate solution given by Longuski, however, is extremely accurate and suitable for computer design work and error analysis of spacecraft performance. The present investigation is concerned with a review of the solutions and the specific regions of validity of the solution for the Eulerian angles.
Analytic error bounds are presented for the solutions of approximate models for self-excited near-symmetric rigid bodies. The error bounds are developed for analytic solutions to Euler's equations of motion. The results are applied to obtain a simplified analytic solution for Eulerian rates and angles. The results of a sample application of the range and error bound expressions for the case of the Galileo spacecraft experiencing transverse torques demonstrate the use of the bounds in analyses of rigid body spin change maneuvers.
Equations for calculating control gains with filters for rigid body motion of space vehicle
The Design Analysis Branch (NE-Ml) at the Kennedy Space Center has not had the ability to accurately couple Rigid Body Dynamics (RBD) and Computational Fluid Dynamics (CFD). OVERFLOW-D is a flow solver that has been developed by NASA to have the capability to analyze and simulate dynamic motions with up to six Degrees of Freedom (6-DOF). Two simulations were prepared over the course of the internship to demonstrate 6DOF motion of rigid bodies under aerodynamic loading. The geometries in the simulations were based on a conceptual Space Launch System (SLS). The first simulation that was prepared and computed was the motion of a Solid Rocket Booster (SRB) as it separates from its core stage. To reduce computational time during the development of the simulation, only half of the physical domain with respect to the symmetry plane was simulated. Then a full solution was prepared and computed. The second simulation was a model of the SLS as it departs from a launch pad under a 20 knot crosswind. This simulation was reduced to Two Dimensions (2D) to reduce both preparation and computation time. By allowing 2-DOF for translations and 1-DOF for rotation, the simulation predicted unrealistic rotation. The simulation was then constrained to only allow translations.
The global linearization (inversion) of rigid-body rotational dynamics is reviewed and representations in terms of quaternions and direction cosines are compared. Certain properties common to quaternions and direction cosines that make their use preferable to Euler angles and that simplify the inversion procedure are described. Applications of the inversion procedure for state estimation and attitude control are discussed. To avoid complexities caused by aerodynamics, an example of direct inversion for linear feedback control of spacecraft attitude is given.