Search NASASearch

SEARCH · Search NASA

Results for “multibody”

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 127 records · Page 7

Connection forces in deformable multibody dynamics

In the dynamic formulation of holonomic and nonholonomic systems based on D'Alembert-Lagrange equation, the forces of constraints are maintained in the dynamic equations by introducing auxiliary variables, called Lagrange multipliers. This approach introduces a set of generalized reaction forces associated with the system generalized coordinates. Different sets of variables can be used as generalized coordinates and accordingly, the generalized reactions associated with these generalized coordinates may not be the actual reaction forces at the joints. In rigid body dynamics, the generalized reaction forces and the actual reaction forces at the joints represent equipollent systems of forces since they produce the same total forces and moments at and about any point on the rigid body. This is not, however, the case in deformable body analyses wherein the generalized reaction forces depend on the system generalized reference and elastic coordinates. In this paper, a method for determining the actual reaction forces at the joints from the generalized reaction forces in deformable multibody systems is presented.

Shabana, A. A.

A computational procedure for multibody systems including flexible beam dynamics

A computational procedure suitable for the solution of equations of motions for flexible multibody systems has been developed. A fully nonlinear continuum approach capable of accounting for both finite rotations and large deformations has been used to model a flexible beam component. The beam kinematics are referred directly to an inertial reference frame such that the degrees of freedom embody both the rigid and flexible deformation motions. As such, the beam inertia expression is identical to that of rigid body dynamics. The nonlinear coupling between gross body motion and elastic deformation is contained in the internal force expression. Numerical solution procedures for the integration of spatial kinematic systems can be directily applied to the generalized coordinates of both the rigid and flexible components. An accurate computation of the internal force term which is invariant to rigid motions is incorporated into the general solution procedure.

Downer, J. D.

Multibody model reduction by component mode synthesis and component cost analysis

The classical assumed-modes method is widely used in modeling the dynamics of flexible multibody systems. According to the method, the elastic deformation of each component in the system is expanded in a series of spatial and temporal functions known as modes and modal coordinates, respectively. This paper focuses on the selection of component modes used in the assumed-modes expansion. A two-stage component modal reduction method is proposed combining Component Mode Synthesis (CMS) with Component Cost Analysis (CCA). First, each component model is truncated such that the contribution of the high frequency subsystem to the static response is preserved. Second, a new CMS procedure is employed to assemble the system model and CCA is used to further truncate component modes in accordance with their contribution to a quadratic cost function of the system output. The proposed method is demonstrated with a simple example of a flexible two-body system.

Spanos, J. T.

Some fundamentals regarding kinematics and generalized forces for multibody dynamics

In order to illustrate the various forms in which generalized forces can arise from diverse subsystem analyses in multibody dynamics, intrinsic dynamical equations for the rotational dynamics of a rigid body are derived from Hamilton's principle. Two types of generalized forces are derived: (1) those associated with the virtual rotation vector in some orthogonal basis, and (2) those associated with varying generalized coordinates. As one physical or kinematical result (such as a frequency or a specific direction cosine) cannot rely on this selection, a 'blind' coupling of two models in which generalized forces are calculated in different ways would be wrong. Both types should use the same rotational coordinates and should denote the virtual rotation on a similar basis according to method 1, or in terms of common rotational coordinates and their diversifications as in method 2. Alternatively, the generalized forces and coordinates of one model may be transformed to those of the other.

Hodges, Dewey H.

Projection and assembly method for multibody component model reduction

The problem addressed is that of obtaining reduced-order component models for use in simulating the dynamics of a multibody system. In certain cases, nonlinear system models may be constructed using linear dynamic models for each component, but allowing large angle motion between components. Without some form of model reduction, system models constructed in this manner may be too large for use in control system design and simulation trades. This paper analyzes one method of component model reduction that allows systems level requirements (e.g., capturing the effect of body 1 reaction wheel noise on body 2 camera pointing) to aid in the selection of the reduced-order component models. Briefly stated, important modes are selected at the system level and projected onto the components, and reduced-order components are then assembled into a reduced-order system model that retains the projected modes.

Bernard, Douglas E.

Selection of component modes for the simulation of flexible multibody spacecraft

This paper describes a procedure for the selection of component modes employed in discretization of component deformation in a flexible multibody spacecraft. The emphasis is placed on the selection of modes which adequately represent the interaction of the various on-board control systems with the vehicle structural flexibility. The method combines the component mode synthesis approaches of Craig-Bampton (1968), MacNeal (1971) and Rubin (1975), and Benfield-Hruda (1971) with the modal balancing method of Moore (1981) and Gregory (1984). The procedure, which is applicable to both articulating and nonarticulating systems, was used to develop a low-order model of the three-body articulating Galileo spacecraft.

Spanos, John T.

Order (n) DISCOS for multibody systems with gear reduction

Recent developments in O(n) algorithms (where n is the number of bodies in the system) and parallel processing have drastically reduced the computer time needed to simulate systems involving many bodies. This paper presents a gear-reduction model for the O(n) version of DISCOS, a standard software package for simulation and analysis of flexible multibody systems. The gear-reduction model allows the accurate modeling of harmonic drives, which are commonly used in robot joints. The formulation has been implemented and validated with known results. The gear model can also be used for gear-train, rack-and-pinion, and screw joints.

Chun, Hon M.

Constraint treatment techniques and parallel algorithms for multibody dynamic analysis

Computational procedures for kinematic and dynamic analysis of three-dimensional multibody dynamic (MBD) systems are developed from the differential-algebraic equations (DAE's) viewpoint. Constraint violations during the time integration process are minimized and penalty constraint stabilization techniques and partitioning schemes are developed. The governing equations of motion, a two-stage staggered explicit-implicit numerical algorithm, are treated which takes advantage of a partitioned solution procedure. A robust and parallelizable integration algorithm is developed. This algorithm uses a two-stage staggered central difference algorithm to integrate the translational coordinates and the angular velocities. The angular orientations of bodies in MBD systems are then obtained by using an implicit algorithm via the kinematic relationship between Euler parameters and angular velocities. It is shown that the combination of the present solution procedures yields a computationally more accurate solution. To speed up the computational procedures, parallel implementation of the present constraint treatment techniques, the two-stage staggered explicit-implicit numerical algorithm was efficiently carried out. The DAE's and the constraint treatment techniques were transformed into arrowhead matrices to which Schur complement form was derived. By fully exploiting the sparse matrix structural analysis techniques, a parallel preconditioned conjugate gradient numerical algorithm is used to solve the systems equations written in Schur complement form. A software testbed was designed and implemented in both sequential and parallel computers. This testbed was used to demonstrate the robustness and efficiency of the constraint treatment techniques, the accuracy of the two-stage staggered explicit-implicit numerical algorithm, and the speed up of the Schur-complement-based parallel preconditioned conjugate gradient algorithm on a parallel computer.

Chiou, Jin-Chern

Spatial operator approach to flexible multibody system dynamics and control

The inverse and forward dynamics problems for flexible multibody systems were solved using the techniques of spatially recursive Kalman filtering and smoothing. These algorithms are easily developed using a set of identities associated with mass matrix factorization and inversion. These identities are easily derived using the spatial operator algebra developed by the author. Current work is aimed at computational experiments with the described algorithms and at modelling for control design of limber manipulator systems. It is also aimed at handling and manipulation of flexible objects.

Rodriguez, G.

Nonlinear dynamic analysis of flexible multibody systems

Two approaches are developed to analyze the dynamic behavior of flexible multibody systems. In the first approach each body is modeled with a modal methodology in a local non-inertial frame of reference, whereas in the second approach, each body is modeled with a finite element methodology in the inertial frame. In both cases, the interaction among the various elastic bodies is represented by constraint equations. The two approaches were compared for accuracy and efficiency: the first approach is preferable when the nonlinearities are not too strong but it becomes cumbersome and expensive to use when many modes must be used. The second approach is more general and easier to implement but could result in high computation costs for a large system. The constraints should be enforced in a time derivative fashion for better accuracy and stability.

Bauchau, Olivier A.

A natural partitioning scheme for parallel simulation of multibody systems

A parallel partitioning scheme based on physical-coordinate variables is presented to systematically eliminate system constraint forces and yield the equations of motion of multibody dynamics systems in terms of their independent coordinates. Key features of the present scheme include an explicit determination of the independent coordinates, a parallel construction of the null space matrix of the constraint Jacobian matrix, an easy incorporation of the previously developed two-stage staggered solution procedure, and Schur complement based parallel preconditioned conjugate gradient numerical algorithm.

Chiou, J. C.

Bilinear modeling and state estimation of slew-induced deformations of a flexible multibody system

A model of the nonlinear dynamics of a deformable maneuvering multibody system is described, whereby elastic deformations are modeled by restoring forces and dissipative forces at point mass appendages. This model is brought into bilinear form. Estimation of deformations occasioned by rapid slewing maneuvers is carried out by a filter based on a globally equivalent linear model of the bilinear dynamics, and is shown to be an improvement over the extended Kalman filter. To further alleviate the computational burden, the estimated deformation state is propagated between observations by a low dimensional operator spline interpolator of bilinear system Volterra series, which is easily implemented.

Dwyer, Thomas A. W., III

Technique to eliminate computational instability in multibody simulations employing the Lagrange multiplier

A programming technique to eliminate computational instability in multibody simulations that use the Lagrange multiplier is presented. The computational instability occurs when the attached bodies drift apart and violate the constraints. The programming technique uses the constraint equation, instead of integration, to determine the coordinates that are not independent. Although the equations of motion are unchanged, a complete derivation of the incorporation of the Lagrange multiplier into the equation of motion for two bodies is presented. A listing of a digital computer program which uses the programming technique to eliminate computational instability is also presented. The computer program simulates a solid rocket booster and parachute connected by a frictionless swivel.

Watts, G.

Multibody interference at transonic Mach numbers

An experimental and computational study has been made of the forces and moments on multibody configurations at Mach numbers from 0.6 to 1.2. The interference forces and moments that occur on symmetric double and triple body configurations at zero angle of attack are examined in detail. In all cases, the normal interference force is always in a direction to pull the bodies together and the moment about the 0.5L point is always to rotate the tips of the bodies away from the center of the configuration. The axial interference force (pressure drag at zero angle of attack) increases the drag per store as additional stores are added to the configuration. The effect of angle of attack on the interference is relatively small.

Chapman, Gary T.

Finite element approach for transient analysis of multibody systems

A three-dimensional, finite element based formulation for the transient dynamics of constrained multibody systems with trusslike configurations is presented. A convected coordinate system is used to define the rigid-body motion of individual elements in the system. Deformation of each element is defined relative to its convected coordinate system. The formulation is oriented toward joint-dominated structures. Through a series of sequential transformations, the joint degree of freedom is built into the equations of motion of the element to reduce geometric constraints. Based on the derivation, a general-purpose code has been developed. Two examples are presented to illustrate the application of the code.

Wu, Shih-Chin

Sub-domain decomposition methods and computational controls for multibody dynamical systems

This paper presents a concurrent methodology to simulate the dynamics of flexible multibody systems with a large number of degrees of freedom. A general class of open-loop structures is treated and a redundant coordinate formulation is adopted. A range space method is used in which the constraint forces are calculated using a preconditioned conjugate gradient method. By using a preconditioner motivated by the regular ordering of the directed graph of the structures, it is shown that the method is order N in the total number of coordinates of the system. The overall formulation has the advantage that it permits fine parallelization and does not rely on system topology to induce concurrency. It can be efficiently implemented on the present generation of parallel computers with a large number of processors. Validation of the method is presented via numerical simulations of space structures incorporating large number of flexible degrees of freedom.

Menon, R. G.

Multibody-plasma interactions - Charging in the wake

Multibody-plasma interactions refers to two or more charged bodies simultaneously interacting with the surrounding plasma as well as each other. This paper concerns a basic type of such interactions in space: the charging of a free flyer in the wake of a large structure. The conditions for severe charging are discussed quantitatively. Computer particle simulations are carried out to obtain the charging potentials in the wake of both floating and biased plate. It is shown that a severe charging zone exists in the near wake of a floating plate within which a free flyer is charged to the KV range under the sun shadow/auroral electron condition. Whether the plate is biased or floating, a large potential difference always exists between the plate and a docking free flyer in the wake. The effects of wake charging on spacecraft docking operations are discussed.

Wang, J.

Dynamic modeling and control of multibody mechanical systems which are structurally flexible

This paper concerns the dynamic modeling and control of multibody, open-chained, structurally flexible, mechanical systems where the bodies are connected by revolute joints. The equations of motion are formulated based on a matrix form of Lagrange's equations for inertial quasi-coordinates. Each body is treated as a substructure of the system. For the purposes of simulation and control, the equations of motion are separated into two sets of equations using a perturbation approach: one to describe large rigid-body motions of the articulated system and the other to describe small linear motions of the bodies about the large motions. A biologically natural control strategy is used for vibration suppression and tracking the prescribed motion.

Quinn, R. D.