Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EQUATIONS OF MOTION”

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,117 records · Page 62

Dynamic stability of a rotor blade using finite element analysis

The aeroelastic stability of flap bending, lead-lag bending, and torsion of a helicopter rotor blade in hover is examined using a finite element formulation based on the principle of virtual work. Quasi-steady two-dimensional airfoil theory is used to obtain the aerodynamic loads. The rotor blade is discretized into beam elements, each with ten modal degrees of freedom. The resulting nonlinear equations of motion are solved for steady-state blade deflections through an iterative procedure. The flutter solution is calculated assuming blade motion to be a small perturbation about the steady solution. The normal mode method based on the coupled rotating natural modes about the steady deflections is used to reduce the number of equations in the flutter eigenanalysis. Results are presented for hingeless and articulated rotor blade configurations.

Sivaneri, N. T.↗

Numerical integration of the satellites of the outer planets

The formulation of a numerical integration program, used to construct ephemerides for satellites of the outer planets, is described. The equations of motion are derived, including (1) n massive integrated satellites, (2) m massive perturbing planets, (3) J2, J4 oblateness coefficients of the primary and each satellite, and (4) a barycentric coordinate system. Variational equations are also found: (1) 6n initial states and n masses of the integrated satellites, (2) J2, J4 of the primary planet, (3) the mass of the system, and (4) right ascension and declination of the pole. The formulation was used to construct the satellite ephemerides of Saturn for the Voyager mission.

Peters, C. F.↗

Tidal disruption and the origin of the Moon

The dynamic problem of the tidal disruption of a rocky planetismal was solved by a direct integration of the fully three-dimensional, nonlinear equations of motion. The hypothesis that any object that passes within the Roche limit is disrupted was disproven. A time dependent solution was performed numerically, treating the planetismal as a fluid with a Murnaghan equation of state in the solid regions and zero pressure otherwise. Calculations show that a rocky body which passes by the Earth on a parabolic orbit with a perigee within the Roche limit is not tidally disrupted. Objects on hyperbolic orbits would experience even less tidal disruption. The results herein do not apply to bodies with very low viscosity. It is shown, however, that tidal disruption can be ruled out as a mechanism for reducing planetismal masses. Mechanisms for forming the Moon which rely upon tidal disruption are unlikely to be correct.

Mizuno, H.↗

Modal identities for multibody elastic spacecraft: An aid to selecting modes for simulation

The question: Which set of modes furnishes a higher fidelity math model of dynamics of a multibody, deformable spacecraft (hinges-free or hinges-locked vehicle modes) is answered. Two sets of general, discretized, linear equations of motion of a spacecraft with an arbitrary number of deformable appendages, each articulated directly to the core body, are obtained using the above two families of modes. By a comparison of these equations, ten sets of modal identities are constructed which involve modal momenta coefficients and frequencies associated with both classes of modes. The sums of infinite series that appear in the identities are obtained in terms of mass, and first and second moments of inertia of the appendages, core body, and vehicle by using certain basic identities concerning appendage modes. Applying the above identities to a four-body spacecraft, the hinges-locked vehicle modes are found to yield a higher fidelity model than hinges-free modes, because the latter modes have nonconverging modal coefficients; a characteristic proved and illustrated.

Hablani, Hari B.↗

Integrated Approach to the Dynamics and Control of Maneuvering Flexible Aircraft

This work uses a fundamental approach to the problem of simulating the flight of flexible aircraft. To this end, it integrates into a single formulation the pertinent disciplines, namely, analytical dynamics, structural dynamics, aerodynamics, and controls. It considers both the rigid body motions of the aircraft, three translations (forward motion, sideslip and plunge) and three rotations (roll, pitch and yaw), and the elastic deformations of every point of the aircraft, as well as the aerodynamic, propulsion, gravity and control forces. The equations of motion are expressed in a form ideally suited for computer processing. A perturbation approach yields a flight dynamics problem for the motions of a quasi-rigid aircraft and an 'extended aeroelasticity' problem for the elastic deformations and perturbations in the rigid body motions, with the solution of the first problem entering as an input into the second problem. The control forces for the flight dynamics problem are obtained by an 'inverse' process and the feedback controls for the extended aeroservoelasticity problem are determined by the LQG theory. A numerical example presents time simulations of rigid body perturbations and elastic deformations about 1) a steady level flight and 2) a level steady turn maneuver.

Waszak, Martin R.↗

Satellite attitude motion models for capture and retrieval investigations

The primary purpose of this research is to provide mathematical models which may be used in the investigation of various aspects of the remote capture and retrieval of uncontrolled satellites. Emphasis has been placed on analytical models; however, to verify analytical solutions, numerical integration must be used. Also, for satellites of certain types, numerical integration may be the only practical or perhaps the only possible method of solution. First, to provide a basis for analytical and numerical work, uncontrolled satellites were categorized using criteria based on: (1) orbital motions, (2) external angular momenta, (3) internal angular momenta, (4) physical characteristics, and (5) the stability of their equilibrium states. Several analytical solutions for the attitude motions of satellite models were compiled, checked, corrected in some minor respects and their short-term prediction capabilities were investigated. Single-rigid-body, dual-spin and multi-rotor configurations are treated. To verify the analytical models and to see how the true motion of a satellite which is acted upon by environmental torques differs from its corresponding torque-free motion, a numerical simulation code was developed. This code contains a relatively general satellite model and models for gravity-gradient and aerodynamic torques. The spacecraft physical model for the code and the equations of motion are given. The two environmental torque models are described.

Cochran, John E., Jr.↗

Effects of large deflection and transverse shear on response of rectangular symmetric composite laminates subjected to acoustic excitation

Nonlinear equations of motion of symmetrically laminated anisotropic plates are derived accounting for von Karman strains. The effect of transverse shear is included in the formulation. Using a single-mode Galerkin procedure, the nonlinear modal equation is obtained. The direct equivalent linearization method is employed for solution of this equation. The response to acoustic excitation of moderately thick composite panels is studied. Further, the effects of transverse shear on large deflection vibration of laminates under random excitation are studied. Mean-square deflections and mean-square inplane stresses are obtained for symmetric graphite-epoxy laminates. Using equilibrium equations and the continuity requirements, the mean-square transverse shear stresses are calculated. The results obtained will be useful in the sonic fatigue design of composite aircraft panels. The analysis is presented in detail for simply supported plates. The analogous equations for a clamped case are also obtained.

Mei, Chuh↗

Gain in computational efficiency by vectorization in the dynamic simulation of multi-body systems

An improved technique for the identification and extraction of the exact quantities associated with the degrees of freedom at the element as well as the flexible body level is presented. It is implemented in the dynamic equations of motions based on the recursive formulation of Kane et al. (1987) and presented in a matrix form, integrating the concepts of strain energy, the finite-element approach, modal analysis, and reduction of equations. This technique eliminates the CPU intensive matrix multiplication operations in the code's hot spots for the dynamic simulation of the interconnected rigid and flexible bodies. A study of a simple robot with flexible links is presented by comparing the execution times on a scalar machine and a vector-processor with and without vector options. Performance figures demonstrating the substantial gains achieved by the technique are plotted.

Amirouche, F. M. L.↗

Magnetic field decay in isolated neutron stars

Three mechanisms that promote the loss of magnetic flux from an isolated neutron star - Ohmic decay, ambipolar diffusion, and Hall drift - are investigated. Equations of motions are solved for charged particles in the presence of a magnetic field and a fixed background of neutrons, while allowing for the creation and destruction of particles by weak interactions. Although these equations apply to normal neutrons and protons, the present interpretations of their solutions are extended to cover cases of neutron superfluidity and proton superconductivity. The equations are manipulated to prove that, in the presence of a magnetic force, the charged particles cannot be simultaneously in magnetostatic equilibrium and chemical equilibrium with the neutrons. The application of the results to real neutron stars is discussed.

Goldreich, Peter↗

Convergence Acceleration of the Navier-Stokes Equations Through Time-Derivative Preconditioning

Chorin's method of artificial compressibility is extended to both compressible and incompressible fluids by using physical arguments to define artificial fluid properties that make up a local preconditioning matrix. In particular, perturbation expansions are used to provide appropriate temporal derivatives for the equations of motion at both low speeds and low Reynolds numbers. These limiting forms are then combined into a single function that smoothly merges into the physical time derivatives at high speeds so that the equations are left unchanged at transonic, high Reynolds number conditions. The effectiveness of the resulting preconditioning procedures for the Navier-Stokes equations is demonstrated for a wide speed and Reynolds number ranges by means of stability results and computational solutions. Nevertheless, the preconditioned equations sometimes fail to provide a solution for applications for which the non-preconditioned equations converge. Often this is because the reduced dissipation in the preconditioned equations results in an unsteady solution while the more dissipative non-preconditioned equations result in a steady state. Problems of this type represent a computational challenge; it is important to distinguish between non-convergence of algorithms, and the non-existence of steady state solutions.

Merkle, Charles L.↗

Orbit, reentry, and landing attachment for globes

Navigational device, invented to aid recovery of spacecraft from any orbit, also illustrates motions of satellites relative to earth and their entry-ranging requirements. Device rapidly and accurately defines lateral range requirements for spacecraft returning to any desired site without manual or computerized calculation of orbital equations of motion.

Pritchard, E. B.↗

Vibration in Planetary Gear Systems with Unequal Planet Stiffnesses

An algorithm suitable for a minicomputer was developed for finding the natural frequencies and mode shapes of a planetary gear system which has unequal stiffnesses between the Sun/planet and planet/ring gear meshes. Mode shapes are represented in the form of graphical computer output that illustrates the lateral and rotational motion of the three coaxial gears and the planet gears. This procedure permits the analysis of gear trains utilizing nonuniform mesh conditions and user specified masses, stiffnesses, and boundary conditions. Numerical integration of the equations of motion for planetary gear systems indicates that this algorithm offers an efficient means of predicting operating speeds which may result in high dynamic tooth loads.

Frater, J. L.↗

A linear quadratic tracker for Control Moment Gyro based attitude control of the Space Station

The paper discusses a design for an attitude control system for the Space Station which produces fast response, with minimal overshoot and cross-coupling with the use of Control Moment Gyros (CMG). The rigid body equations of motion are linearized and discretized and a Linear Quadratic Regulator (LQR) design and analysis study is performed. The resulting design is then modified such that integral and differential terms are added to the state equations to enhance response characteristics. Methods for reduction of computation time through channelization are discussed as well as the reduction of initial torque requirements.

Kaidy, J. T.↗

Semi-analytical theory for long-term rotational motion of uncontrolled satellites

Analytical solutions/integrals of the free rotational motion of uncontrolled satellites which contain spinning rotors are employed in conjunction with a suitable analytical/numerical averaging approach to predict long-term attitude motion. This semi-analytical method of averaging is basically similar to the one used in purely analytical 'general perturbation' theories. However, in view of the complexities involved in analytical treatment of the problem, its application is primarily numerical in character. An example is presented illustrating the use of the theory to predict the attitude motion of an axisymmetric satellite containing symmetrical constant-speed rotors under gravity-gradient torque. The use of the theory is considerably more efficient than straightforward numerical integration of the complete set of equations of motion.

Lahr, B. S.↗

Impingement of water droplets on wedges and diamond airfoils at supersonic speeds

An analytical solution has been obtained for the equations of motion of water droplets impinging on a wedge in a two-dimensional supersonic flow field with a shock wave attached to the wedge. The closed-form solution yields analytical expressions for the equation of the droplet trajectory, the local rate of impingement and the impingement velocity at any point on the wedge surface, and the total rate of impingement. The analytical expressions are utilized to determine the impingement on the forward surfaces of diamond airfoils in supersonic flow fields with attached shock waves. The results presented include the following conditions: droplet diameters from 2 to 100 microns, pressure altitudes from sea level to 30,000 feet, free-stream static temperatures from 420 degrees to 460 degrees R. Also, free-stream Mach numbers from 1.1 to 2.0, semi-apex angles for the wedge from 1.14 degrees to 7.97 degrees, thickness-to-chord ratios for the diamond airfoil from 0.02 to 0.14, chord lengths from 1 to 20 feet, and angles of attack from zero to the inverse tangent of the airfoil thickness-to-chord ratio.

SURFACE CONDITIONS - WING SECTIONS↗

Impingement of water droplets on wedges and double-wedge airfoils at supersonic speeds

An analytical solution has been obtained for the equations of motion of water droplets impinging on a wedge in a two-dimensional supersonic flow field with a shock wave attached to the wedge. The closed-form solution yields analytical expressions for the equation of the droplet trajectory, the local rate of impingement and the impingement velocity at any point on the wedge surface, and the total rate of impingement. The analytical expressions are utilized to determine the impingement on the forward surfaces of diamond airfoils in supersonic flow fields with attached shock waves. The results presented include the following conditions: droplet diameters from 2 to 100 microns, pressure altitudes from sea level to 30,000 feet, free-stream static temperatures from 420 degrees r, free stream Mach numbers from 1.1 to 2.0, semiapex angles for the wedge from 1.14 degrees to 7.97 degrees, thickness-to-chord ratios for the diamond airfoil from 0.02 to 0.14, chord lengths from 1 to 20 feet, and angles of attack from zero to the inverse tangent of the airfoil thickness-to-chord ratio.

Serafini, John S↗

Approximate optimal guidance for the advanced launch system

A real-time guidance scheme for the problem of maximizing the payload into orbit subject to the equations of motion for a rocket over a spherical, non-rotating earth is presented. An approximate optimal launch guidance law is developed based upon an asymptotic expansion of the Hamilton - Jacobi - Bellman or dynamic programming equation. The expansion is performed in terms of a small parameter, which is used to separate the dynamics of the problem into primary and perturbation dynamics. For the zeroth-order problem the small parameter is set to zero and a closed-form solution to the zeroth-order expansion term of Hamilton - Jacobi - Bellman equation is obtained. Higher-order terms of the expansion include the effects of the neglected perturbation dynamics. These higher-order terms are determined from the solution of first-order linear partial differential equations requiring only the evaluation of quadratures. This technique is preferred as a real-time, on-line guidance scheme to alternative numerical iterative optimization schemes because of the unreliable convergence properties of these iterative guidance schemes and because the quadratures needed for the approximate optimal guidance law can be performed rapidly and by parallel processing. Even if the approximate solution is not nearly optimal, when using this technique the zeroth-order solution always provides a path which satisfies the terminal constraints. Results for two-degree-of-freedom simulations are presented for the simplified problem of flight in the equatorial plane and compared to the guidance scheme generated by the shooting method which is an iterative second-order technique.

Feeley, T. S.↗