A computer-oriented method for reducing linearized multibody system equations by incorporating constraints

Consider a spatial multibody system with rigid and elastic bodies. The bodies are linked by rigid interconnections (e.g. revolute joints) causing constraints, as well as by flexible interconnections (e.g. springs) causing applied forces. Small motions of the system with respect to a given nominal configuration can be described by linearized dynamic equations and kinematic constraint equations. We present a computer-oriented procedure which allows to develop a minimum number of these equations. There are three problems. First: algorithmic selection of position coordinates; second: condensation of the dynamic equations; third: evaluation of the constraint forces. To demonstrate the procedure, a closed loop multibody system is used as an example.     

On the First-Order Decoupling of Dynamical Equations of Motion for Elastic Multibody Systems as Applied to a Two-Link Flexible Manipulator

An improved method for deriving elastic generalized coordinatesis considered. Then Kane's equations of motion for multibody systemsconsisting of an arbitrary number of rigid and elastic bodies ispresented. The equations are in general form and are applicable for anydesired holonomic system. Flexibility in choosing generalized speeds interms of generalized coordinate derivatives in Kane's method is used. Itis shown that proper choice of a congruency transformation betweengeneralized coordinate derivatives and generalized speeds leads toequations of motion for holonomic multibody systems consisting of anarbitrary number of rigid and elastic bodies. These equations aredecoupled in first-order terms. In order to show the use of this method,a simple system consisting of a lumped mass, a spring and a clamped-freeelastic beam is modeled. Finally, the numerical implementation ofdecoupling using congruency transformation is discussed and shown viasimulation of a two-degrees-of-freedom flexible robot.     

Dynamic analysis of flexible multi-body systems using mixed modal and tangent coordinates

This paper presents a mixed modal and tangent coordinate technique for computer aided analysis of flexible mechanical systems whose components undergo large translations and large rotations. In this model the configuration of a flexible component is identified by using two sets of generalized coordinates, namely rigid body and elastic coordinates. The rigid body coordinates define the location and orientation of a body axis, whereas the elastic coordinates define the displacement field of a component with respect to its body axis. The elastic coordinates are introduced by using finite element discretization to model flexible components with complex geometries. A modal analysis technique is used to identify the elastic mode shapes and to eliminate insignificant higher frequency modes. An orthonormalization of constraint Jacobian matrix associated with rigid body coordinates is used to identify the rigid body tangent coordinates. The resulting modal and tangent coordinates are used to develop an automated numerical integration scheme to solve the system differential and algebraic equations. Two numerical examples are considered to show the feasibility of dynamic analysis of flexible mechanical systems using this scheme.     

Quadratic and Higher-Order Constraints in Energy-Conserving Formulations of Flexible Multibody Systems

The treatment of constraints is considered here within the framework ofenergy-momentum conserving formulations for flexible multibody systems.Constraint equations of various types are an inherent component of multibodysystems, their treatment being one of the key performance features ofmathematical formulations and numerical solution schemes.Here we employ rotation-free inertial Cartesian coordinates of points tocharacterise such systems, producing a formulation which easily couples rigidbody dynamics with nonlinear finite element techniques for the flexiblebodies. This gives rise to additional internal constraints in rigid bodies topreserve distances. Constraints are enforced via a penalty method, which givesrise to a simple yet powerful formulation. Energy-momentum time integrationschemes enable robust long term simulations for highly nonlinear dynamicproblems.The main contribution of this paper focuses on the integration of constraintequations within energy-momentum conserving numerical schemes. It is shownthat the solution for constraints which may be expressed directly in terms ofquadratic invariants is fairly straightforward. Higher-order constraints mayalso be solved, however in this case for exact conservation an iterativeprocedure is needed in the integration scheme. This approach, together withsome simplified alternatives, is discussed.Representative numerical simulations are presented, comparing the performanceof various integration procedures in long-term simulations of practicalmultibody systems.     

Dynamics of Flexible Multibody Systems Using Virtual Work and Linear Graph Theory

By combining linear graph theory with the principle of virtualwork, a dynamic formulation is obtained that extends graph-theoreticmodelling methods to the analysis of flexible multibody systems. Thesystem is represented by a linear graph, in which nodes representreference frames on rigid and flexible bodies, and edges representcomponents that connect these frames. By selecting a spanning tree forthe graph, the analyst can choose the set of coordinates appearing inthe final system of equations. This set can include absolute, joint, orelastic coordinates, or some combination thereof. If desired, allnon-working constraint forces and torques can be automaticallyeliminated from the dynamic equations by exploiting the properties ofvirtual work. The formulation has been implemented in a computerprogram, DynaFlex, that generates the equations of motion in symbolicform. Three examples are presented to demonstrate the application of theformulation, and to validate the symbolic computer implementation.     

Hybrid control of flexible multi-body systems

In constrained systems of rigid and flexible bodies, the gross rigid body motion and elastic deformation cannot be controlled independently because of the coupling between these two motions. A hybrid control method for suppressing the vibration of a geometrically nonlinear flexible multi-body system is proposed in this paper. This method utilizes both the passive and active control concepts. In the passive control strategy, flexible components in the system are manufactured from fiber-reinforced composite laminates which have high strength-to-weight and stiffness-to-weight ratios. On the other hand, the active control scheme used in this paper utilizes measurable velocity and acceleration signals to produce the command signals required to activate the actuator forces. A small number of sensors and controllers with constant gain factors are used in order to obtain a low-cost and simple control system. The generalized active control forces associated with the system generalized coordinates are developed using the virtual work and are written in terms of the coupled set of reference and elastic coordinates. The system differential equations of motion are developed using Lagrange's equation and the Jacobian matrix of the nonlinear algebraic constraint equations describing mechanical joints in the system is used to identify a set of independent generalized coordinates. The associated independent differential equations are identified and are written in the state space formulation. The characteristics of the proposed hybrid control are evaluated through computer simulations of a seven-body flexible vehicle. The performance characteristics of the hybrid control are also compared to the performance characteristics of the passive and active controls.     

Pulse control of flexible multibody systems

An active pulse control method is developed to reduce the vibrations of multibody systems resulting from impact loadings. The pulse, which is a function of system generalized coordinates and velocities, is determined analytically using energy and momentum balance equations of the impacting bodies. Elastic components in the multibody system are discretized using the finite element method. The system equations of motions and nonlinear algebraic constraint equations describing mechanical joints between different components are written in the Lagrangian formulation using a finite set of coupled reference position and local elastic generalized coordinates. A set of independent differential equations are identified by the generalized coordinate partitioning of the constraint Jacobian matrix. These equations are written in the state space formulation and integrated forward in time using a direct numerical integration method. Dependent coordinates are then determined using the constraint kinematic relations. Points in time at which impact occurs are monitored by an impact predictor function, which controls the integration algorithms and forces for the solution of the momentum relation, to define the jump discontinuities in the composite velocity vector as well as the system reaction forces. The effectiveness of the active pulse control in reducing the vibration of flexible multibody aircraft during the touchdown impact is investigated and numerical results are presented.     

Chain Vibration and Dynamic Stress in Three-Dimensional Multibody Tracked Vehicles

A three-dimensional computational finite element procedure for the vibration and dynamic stress analysis of the track link chains of off-road vehicles is presented in this paper. The numerical procedure developed in this investigation integrates classical constrained multibody dynamics methods with finite element capabilities. The nonlinear equations of motion of the three-dimensional tracked vehicle model in which the track link s are considered flexible bodies, are obtained using the floating frame of reference formulation. Three-dimensional contact force models are used to describe the interaction of the track chain links with the vehicle components and the ground. The dynamic equations of motion are first presented in terms of a coupled set of reference and elastic coordinates of the track links. Assuming that the structural flexibility of the track links does not have a significant effect on their overall rigid body motion as well as the vehicle dynamics, a partially linearized set of differential equations of motion of the track links is obtained. The equations associated with the rigid body motion are used to predict the generalized contact, inertia, and constraint forces associated with the deformation degrees of freedom of the track links. These forces are introduced to the track link flexibility equations which are used to calculate the deformations of the links resulting from the vehicle motion. A detailed three-dimensional finite element model of the track link is developed and utilized to predict the natural frequencies and mode shapes. The terms that represent the rigid body inertia, centrifugal and Coriolis forces in the equations of motion associated with the elastic coordinates of the track link are described in detail. A computational procedure for determining the generalized constraint forces associated with the elastic coordinates of the deformable chain links is presented. The finite element model is then used to determine the deformations of the track links resulting from the contact, inertia, and constraint forces. The results of the dynamic stress analysis of the track links are presented and the differences between these results and the results obtained by using the static stress analysis are demonstrated.     

Erratum to: New efficient method for dynamic modeling and simulation of flexible multibody systems moving in plane

Efficient, precise dynamic analysis for general flexible multibody systems has become a research focus in the field of flexible multibody dynamics. In this paper, the finite element method and component mode synthesis are introduced to describe the deformations of the flexible components, and the dynamic equations of flexible bodies moving in plane are deduced. By combining the discrete time transfer matrix method of multibody system with these dynamic equations of flexible component, the transfer equations and transfer matrices of flexible bodies moving in plane are developed. Finally, a high-efficient dynamic modeling method and its algorithm are presented for high-speed computation of general flexible multibody dynamics. Compared with the ordinary dynamics methods, the proposed method combines the strengths of the transfer matrix method and finite element method. It does not need the global dynamic equations of system and has the low order of system matrix and high computational efficiency. This method can be applied to solve the dynamics problems of flexible multibody systems containing irregularly shaped flexible components. It has advantages for dynamic design of complex flexible multibody systems. Formulations as well as a numerical example of a multi-rigid-flexible-body system containing irregularly shaped flexible components are given to validate the method.     

Transient analysis of flexible multi-body systems. Part II: Application to aircraft landing

This investigation presents a method for transient analysis of a large-scale multi-body aircraft consisting of interconnected rigid and flexible bodies that undergo large angular rotations. Elastic components of the aircraft are discretized using the finite element method. The system equations of motion and nonlinear algebraic constraint equations describing joints between different components are written in the Lagrangian formulation using a finite set of coupled reference and modal coordinates. The system differential equations of motion and algebraic constraint equations are computer-generated and integrated forward in time using an explicit-implicit direct numerical integration algorithm coupled with a Newton-Raphson type iteration in order to check on constraint violations. Impact and intermittent motion events are accounted for by using a generalized momentum balance that predict jump discontinuities in the generalized velocities as well as jump discontinuities in the system reaction forces. The formulation presented and the computer program developed are used to simulate the impact between the landing gear and the runway. The method is also used to predict the dynamic behavior of the aircraft during the traverse of an abrupt elevation change in the runway.     

New efficient method for dynamic modeling and simulation of flexible multibody systems moving in plane

Efficient, precise dynamic analysis for general flexible multibody systems has become a research focus in the field of flexible multibody dynamics. In this paper, the finite element method and component mode synthesis are introduced to describe the deformations of the flexible components, and the dynamic equations of flexible bodies moving in plane are deduced. By combining the discrete time transfer matrix method of multibody system with these dynamic equations of flexible component, the transfer equations and transfer matrices of flexible bodies moving in plane are developed. Finally, a high-efficient dynamic modeling method and its algorithm are presented for high-speed computation of general flexible multibody dynamics. Compared with the ordinary dynamics methods, the proposed method combines the strengths of the transfer matrix method and finite element method. It does not need the global dynamic equations of system and has the low order of system matrix and high computational efficiency. This method can be applied to solve the dynamics problems of flexible multibody systems containing irregularly shaped flexible components. It has advantages for dynamic design of complex flexible multibody systems. Formulations as well as a numerical example of a multi-rigid-flexible-body system containing irregularly shaped flexible components are given to validate the method.     

Spatial joint constraints for the absolute nodal coordinate formulation using the non-generalized intermediate coordinates

In this investigation, a systematic procedure that can be used for modeling joint constraints for the absolute nodal coordinate formulation is developed. To this end, the non-generalized intermediate coordinates are introduced to derive a mapping between the generalized gradient coordinates and the non-generalized rotation parameters. With this mapping, a wide variety of joint constraints can be defined for the absolute nodal coordinate formulation in terms of the non-generalized reference coordinates and, therefore, existing well-developed constraint libraries formulated for the rigid body reference coordinates can be directly employed without significant modifications in existing codes. Furthermore, in order to define a rigid surface at the joint definition point, a set of orthonormality conditions is imposed on the gradient coordinates. This leads to not only accurate modeling of interface to mechanical joint, but also a simpler definition of the joint coordinate system obtained by the orthonormal gradient vectors. For this reason, a simpler form of constraint Jacobian and quadratic velocity vectors can be obtained as compared to those of the existing approach which requires the use of highly nonlinear joint coordinate system. A systematic procedure for eliminating the non-generalized coordinates and the dependent Lagrange multipliers associated with the coordinate mapping equations from the equations of motion is presented. As a result, a standard augmented form of the equations of motion can be obtained in terms of the generalized coordinates only. Several numerical examples are presented in order to demonstrate the use of the joint constraint formulation developed in this investigation.     

A recursive, numerically stable, and efficient simulation algorithm for serial robots with flexible links

A methodology for the formulation of dynamic equations of motion of a serial flexible-link manipulator using the decoupled natural orthogonal complement (DeNOC) matrices, introduced elsewhere for rigid bodies, is presented in this paper. First, the Euler Lagrange (EL) equations of motion of the system are written. Then using the equivalence of EL and Newton–Euler (NE) equations, and the DeNOC matrices associated with the velocity constraints of the connecting bodies, the analytical and recursive expressions for the matrices and vectors appearing in the independent dynamic equations of motion are obtained. The analytical expressions allow one to obtain a recursive forward dynamics algorithm not only for rigid body manipulators, as reported earlier, but also for the flexible body manipulators. The proposed simulation algorithm for the flexible link robots is shown to be computationally more efficient and numerically more stable than other algorithms present in the literature. Simulations, using the proposed algorithm, for a two link arm with each link flexible and a Space Shuttle Remote Manipulator System (SSRMS) are presented. Numerical stability aspects of the algorithms are investigated using various criteria, namely, the zero eigenvalue phenomenon, energy drift method, etc. Numerical example of a SSRMS is taken up to show the efficiency and stability of the proposed algorithm. Physical interpretations of many terms associated with dynamic equations of flexible links, namely, the mass matrix of a composite flexible body, inertia wrench of a flexible link, etc. are also presented. The work has been carried out in the Dept. of Mechanical Engineering, Indian Institute of Technology Delhi, New Delhi 110016, India.     

考虑刚-柔-热耦合的板结构多体系统的动力学建模

对热载荷作用下中心刚体与大变形薄板多体系统的动力学建模问题进行研究.基于Kirchhoff假设,从格林应变和曲率与绝对位移的非线性关系式出发,推导了非线性广义弹性力阵,用绝对节点坐标法建立了大变形矩形薄板的有限元离散的动力学变分方程.为了考虑刚体姿态运动、弹性变形和温度变化的相互耦合作用,推导了热流密度与绝对节点坐标之间的关系式.引入系统的运动学约束方程,建立了中心刚体-矩形板多体系统的考虑刚-柔-热耦合的热传导方程和带拉格朗日乘子的第一类拉格朗日动力学方程.为了有效地提高计算效率,将改进的中心差分法和广义-α法相结合,求解热传导方程和动力学方程,差分后的方程通过牛顿迭代法耦合求解.对刚-柔耦合和刚-柔-热三者耦合两种模型的仿真结果进行比较表明,刚体运动对温度梯度和热变形的影响显著.此外,本文建模方法考虑了几何非线性项,因此也考虑了热膨胀引起的轴向变形对横向变形的影响.     

A systematic formulation of dynamic equations for robot manipulators with elastic links

Although a variety of formulation schemes for the dynamic equations of robot manipulators with rigid links can be found in the literature, an efficient method of the formulation for robot manipulators with elastic links is not well known. Accordingly, this work presents the derivation of the equations of motion for application to mechanical manipulators with elastic links. The formulation is conducted analytically using Hamilton's principle. The resultant equations consist of the terms of inertial, Coriolis, centrifugal, gravitational, and exerted forces. They are expressed in terms of a set of independent generalized coordinates. In contrast to conventional variational approaches, the present method provides an efficient and systematic way for obtaining the compact symbolic equations of flexible manipulator systems. Two examples are presented to illustrate the proposed methodology. Firstly, a three-link flexible manipulator with three revolute joints is studied. A flexible manipulator consisting of a prismatic joint and a discrete mass is the second model.     

Absolute nodal coordinate plane beam formulation for multibody systems dynamics

A new plane beam dynamic formulation for constrained multibody system dynamics is developed. Flexible multibody system dynamics includes rigid body dynamics and superimposed vibratory motions. The complexity of mechanical system dynamics originates from rotational kinematics, but the natural coordinate formulation does not use rotational coordinates, so that simple dynamic formulation is possible. These methods use only translational coordinates and simple algebraic constraints. A new formulation for plane flexible multibody systems are developed utilizing the curvature of a beam and point masses. Using absolute nodal coordinates, a constant mass matrix is obtained and the elastic force becomes a nonlinear function of the nodal coordinates. In this formulation, no infinitesimal or finite rotation assumptions are used and no assumption on the magnitude of the element rotations is made. The distributed body mass and applied forces are lumped to the point masses. Closed loop mechanical systems consisting of elastic beams can be modeled without constraints since the loop closure constraints can be substituted as beam longitudinal elasticity. A curved beam is modeled automatically. Several numerical examples are presented to show the effectiveness of this method.     

Maximum allowable load of flexible manipulators for given dynamic trajectory

This paper presents the formulation and numerical solution of the dynamic load carrying capacity (DLCC) problem of flexible manipulators. For manipulators under the rigid body assumption, the major limiting factor in determining the maximum allowable load (load mass and load moment of inertia) for a prescribed dynamic trajectory (positions, velocities and accelerations) is the joint actuator capacity. But for a flexible robot, an additional constraint on allowable deformation at the end effector must be imposed because either lighter-weight links or operating at a higher speed could cause unacceptable fluctuations when moving along a trajectory. A Lagrangian assumed mode method was used to model the manipulator and load dynamics, including both joint and deflection motions. The deflection equations are then coupled with robot kinematics to solve for the generalized coordinates. A strategy to determine the DLCC subject to both constraints mentioned above is formulated where the end effector deflection is specified in terms of a series of spherical bounds with a radius equal to the allowable deformation. A general computational procedure for the multiple-link case given arbitrary trajectories is described in detail. Symbolic derivation and simulation by using a PC-based symbolic language MATHEMATICA® was carried out for a two-link planer robot. The results confirmed the necessity of the dual constraints and showed that which constraint is more critical for a given robot and trajectory depends on the required accuracy.     

Coupling characteristics of rigid body motion and elastic deformation of a 3-PRR parallel manipulator with flexible links

Modeling of multibody dynamics with flexible links is a challenging task, which not only involves the effect of rigid body motion on elastic deformations, but also includes the influence of elastic deformations on rigid body motion. This paper presents coupling characteristics of rigid body motions and elastic motions of a 3-PRR parallel manipulator with three flexible intermediate links. The intermediate links are modeled as Euler–Bernoulli beams with pinned-pinned boundary conditions based on the assumed mode method (AMM). Using Lagrange multipliers, the fully coupled equations of motions of the flexible parallel manipulator are developed by incorporating the rigid body motions with elastic motions. The mutual dependence of elastic deformations and rigid body motions are investigated from the analysis of the derived equations of motion. Open-loop simulation without joint motion controls and closed-loop simulation with joint motion controls are performed to illustrate the effect of elastic motion on rigid body motions and the coupling effect amongst flexible links. These analyses and results provide valuable insight to the design and control of the parallel manipulator with flexible intermediate links.