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 433 records · Page 24

Toward higher order tests of the gravitational interaction

Analyses and interpretations of experiments which test post-Newtonian gravity are usually done under the assumption that gravity is a metric field phenomenon - a manifestation of space-time geometry. This, however, is unnecessary and one can start at a more primitive level - that there simply exists a phenomenological, gravitational, many-body equation of motion which must be determined by a package of observations. In fact, over the last couple decades, a diverse collection of solar system interbody tracking observations, supplemented by data from the binary pulsar system PSR 1913 + 16, has completely mapped out the first post-Newtonian order. After the fact, using empirically determined equations of motion, along with some observed properties of nongravitational clocks and rulers and conservation laws for energy, momentum and angular momentum, a post-Newtonian Lagrangian can be constructed, a geometrical space-time metric field conceptual interpretation can be developed, Lorentz invariance of the equations of motion can be shown, and the equations of motion are found to agree with the predictions of Einstein's gravitational theory, General Relativity, within experimental accuracy.

Nordtvedt, Ken↗

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗

Staggered solution procedures for multibody dynamics simulation

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Park, K. C.↗

Tracking and Orbit-Determination Program of the Jet Propulsion Laboratory

The lunar-probe tracking program at the Jet Propulsion Laboratory has two prime objectives: (1) provide real-time predictions of the direction of the probe from various observation stations; (2) establish a reliable trajectory corresponding to the actual flight path of the probe. The tracking program, although developed for use with lunar probes, can be used for interplanetary probes if certain modifications are made. The program, as developed for the IBM 704 digital computer, has two distinct phases. First, the equations of motion and the variational equations are integrated to each observation time where the elements of the equation A (sub u) equals b [linearization of the maximum likelihood equations] are computed. The second phase is concerned with the solution of a specified subset of A (sub u) equals b. Flexibility and ease of operation have been major objectives in writing the 704 program. The number of data points and tracking stations that may be used is limited only by computing time and core storage. Input formats and operating instructions are presented for utilizing the various computational options available in the program.

Carr, Russell E.↗

Pulsation theory.

Pulsation in cepheid instability strip in H-R DIAGRAM, giving equations for spherically symmetric motion

BINARY STAR↗

The theoretical tools of experimental gravitation

Theoretical frameworks for testing relativistic gravity are presented in terms of a system for analyzing theories of gravity invented as alternatives to Einstein. The parametrized post-Newtonian (PPN) formalism, based on the Dicke framework and the Eotvos-Dicke-Braginsky experiment, is discussed in detail. The metric theories of gravity, and their post-Newtonian limits are reviewed, and PPN equations of motion are derived. These equations are used to analyze specific effects and experimental tests in the solar system.

Will, C. M.↗

On the elliptic restricted three-body problem.

This article describes the most important features of the elliptic restricted three-body problem. The methods of numerical integration with recurrent power series are developed for both the equations of motion and the variational equations. The conditions for the existence of periodic orbits and families of periodic orbits are also outlined in detail, and finally illustrated with a family of symmetric periodic orbits.

Broucke, R.↗

Bifilar analysis users manual, volume 2

The digital computer program developed to study the vibration response of a coupled rotor/bifilar/airframe coupled system is described. The theoretical development of the rotor/airframe system equations of motion is provided. The fuselage and bifilar absorber equations of motion are discussed. The modular block approach used in the make-up of this computer program is described. The input data needed to run the rotor and bifilar absorber analyses is described. Sample output formats are presented and discussed. The results for four test cases, which use the major logic paths of the computer program, are presented. The overall program structure is discussed in detail. The FORTRAN subroutines are described in detail.

Cassarino, S. J.↗

Mass transfer between binary stars

The transfer of mass from one component of a binary system to another by mass ejection is analyzed through a stellar wind mechanism, using a model which integrates the equations of motion, including the energy equation, with an initial static atmosphere and various temperature fluctuations imposed at the base of the star's corona. The model is applied to several situations and the energy flow is calculated along the line of centers between the two binary components, in the rotating frame of the system, thereby incorporating the centrifugal force. It is shown that relatively small disturbances in the lower chromosphere or photosphere can produce mass loss through a stellar wind mechanism, due to the amplification of the disturbance propagating into the thinner atmosphere. Since there are many possible sources of the disturbance, the model can be used to explain many mass ejection phenomena.

Modisette, J. L.↗

Mixing of a passive scalar in isotropic and sheared homogeneous turbulence

In order to calculate the velocity and scalar fields, the three dimensional, time-dependent equations of motion and the diffusion equation were solved numerically. The following cases were treated: isotropic, homogeneous turbulence with decay of a passive scalar; and homogeneous turbulent shear flow with a passive scalar whose mean varies linearly in the spanwise direction. The solutions were obtained at relatively low Reynolds numbers so that all of the turbulent scales could be resolved without modeling. Turbulent statistics such as integral length scales, Taylor microscales, Kolmogorov length scale, one- and two-point correlations of velocity-velocity and velocity-scalar, turbulent Prandtl/Schmidt number, r.m.s. values of velocities, the scalar quantity and pressure, skewness, decay rates, and decay exponents were calculated. The results are compared with the available expermental results, and good agreement is obtained.

Shirani, E.↗

Computer program for aerodynamic and blading design of multistage axial-flow compressors

A code for computing the aerodynamic design of a multistage axial-flow compressor and, if desired, the associated blading geometry input for internal flow analysis codes is presented. Compressible flow, which is assumed to be steady and axisymmetric, is the basis for a two-dimensional solution in the meridional plane with viscous effects modeled by pressure loss coefficients and boundary layer blockage. The radial equation of motion and the continuity equation are solved with the streamline curvature method on calculation stations outside the blade rows. The annulus profile, mass flow, pressure ratio, and rotative speed are input. A number of other input parameters specify and control the blade row aerodynamics and geometry. In particular, blade element centerlines and thicknesses can be specified with fourth degree polynomials for two segments. The output includes a detailed aerodynamic solution and, if desired, blading coordinates that can be used for internal flow analysis codes.

Crouse, J. E.↗

Application of modal control to wing-flutter suppression

A discrete modal control design approach that is applied to a single control surface, unswept aircraft wing subject to bending torsion flutter is described. The modal approach is a mathematical method to decouple the equations of motion into isolated differential equations. In this paper, a pole-placement approach is then applied to determine stability gains in the discrete plane using only the two complex-conjugate flutter-mode equations. A fixed gain Kalman filter is used to estimate the model amplitudes using three measurements. Results are presented for a full-state estimator (36 states) and two reduced state estimators using two different closed-loop pole locations. The control law is designed for a dynamic pressure that is 50 percent greater than the uncontrolled-flutter dynamic pressure. With constant control-law gains, the closed-loop system remains stable over the dynamic-pressure range from flutter onset to approximately an 80-percent increase in pressure.

Ostroff, A. J.↗

Post-instability behavior of solids

The necessity of model reformulation in elasticity results from the failure of hyperbolicity of the governing equations of motion for classical models. The reformulation is based upon the introduction of additional kinematical microstructures in the form of multivalued displacement and velocity field (or fractal functions) which are generated by the mechanism of the instability. The small scale motions describing this microstructure interact with the original large scale motion and restore the hyperbolicity of new governing equations of motion. The applications of the reformulated models to the problem of vibrational control and impact energy absorption are discussed.

Zak, M.↗

Dynamics of articulated structures

An articulated structure is defined as an assembly of flexible bodies that may be coupled by kinematic connections and force elements that permit large relative displacement and rotation. Kinematics of such systems is defined using one reference frame for each body in the system and deformation modal coordinates that define displacement fields within flexible bodies. Deformation kinematics are defined by both elastic vibration and static correction deformation modes. Linear elastic deformation is presumed; i.e., a linear stress-strain relation is valid and relative displacements within each elastic component are small enough so that the theory of linear elasticity applies. Coupling of reference and modal coordinates leads to a system of nonlinear equations of motion. Methods of automatically generating and solving these equations of motion are outlined.

Haug, Edward J.↗

Optimal trajectories for the aeroassisted flight experiment, 1988-89

Research is summarized on optimal trajectories for the aeroassisted flight experiment, performed by the Aero-Astronautics Group of Rice University during the period 1988 through 1989. This research includes the following topics: (1) equations of motion in an Earth-fixed system; (2) equations of motion in an inertial system; (3) formultion of the optimal trajectory problem; (4) results on the optimal trajectory problem; and (5) guidance implications.

Miele, A.↗

Implementation of 3-D isoparametric finite elements on supercomputer for the formulation of recursive dynamical equations of multi-body systems

A computational algorithmic procedure is developed and implemented for the dynamic analysis of a multibody system with rigid/flexible interconnected bodies. The algorithm takes into consideration the large rotation/translation and small elastic deformations associated with the rigid-body degrees of freedom and the flexibility of the bodies in the system respectively. Versatile three-dimensional isoparametric brick elements are employed for the modeling of the geometric configurations of the bodies. The formulation of the recursive dynamical equations of motion is based on the recursive Kane's equations, strain energy concepts, and the techniques of component mode synthesis. In order to minimize CPU-intensive matrix multiplication operations and speed up the execution process, the concepts of indexed arrays is utilized in the formulation of the equations of motion. A spin-up maneuver of a space robot with three flexible links carrying a solar panel is used as an illustrative example.

Shareef, N. H.↗

Modeling the Benchmark Active Control Technology Wind-Tunnel Model for Application to Flutter Suppression

This paper describes the formulation of a model of the dynamic behavior of the Benchmark Active Controls Technology (BACT) wind-tunnel model for application to design and analysis of flutter suppression controllers. The model is formed by combining the equations of motion for the BACT wind-tunnel model with actuator models and a model of wind-tunnel turbulence. The primary focus of this paper is the development of the equations of motion from first principles using Lagrange's equations and the principle of virtual work. A numerical form of the model is generated using values for parameters obtained from both experiment and analysis. A unique aspect of the BACT wind-tunnel model is that it has upper- and lower-surface spoilers for active control. Comparisons with experimental frequency responses and other data show excellent agreement and suggest that simple coefficient-based aerodynamics are sufficient to accurately characterize the aeroelastic response of the BACT wind-tunnel model. The equations of motion developed herein have been used to assist the design and analysis of a number of flutter suppression controllers that have been successfully implemented.

Waszak, Martin R.↗