Multibody Dynamics of Very Flexible Damped Systems

An efficient multibody dynamics formulation is presented for simulating the forward dynamics of open and closed loop mechanical systems comprised of rigid and flexible bodies interconnected by revolute, prismatic, free, and fixed joints. Geometrically nonlinear deformation of flexible bodies is included and the formulation does not impose restrictions on the representation of material damping within flexible bodies.The approach is based on Kane's equation without multipliers and the resulting formulation generates 2ndof+m first order ordinary differential equations directly where ndof is the smallest number of system degrees of freedom that can completely describe the system configuration and m is the number of loop closure velocity constraint equations. The equations are integrated numerically in the time domain to propagate the solution.Flexible bodies are discretized using a finite element approach. The mass and stiffness matrices for a six-degree-of-freedom planar beam element are developed including mass coupling terms, rotary inertia, centripetal and Coriolis forces, and geometric stiffening terms.The formulation is implemented in the general purpose multibody dynamics computer program flxdyn. Extensive validation of the formulation and corresponding computer program is accomplished by comparing results with analytically derived equations, alternative approximate solutions, and benchmark problems selected from the literature. The formulation is found to perform well in terms of accuracy and solution efficiency.This article develops the formulation and presents a set of validation problems including a sliding pendulum, seven link mechanism, flexible beam spin-up problem, and flexible slider crank mechanism.     

Optimal Trajectory Planning for Flexible Link Manipulators with Large Deflection Using a New Displacements Approach

The main objective of the present paper is to determine the optimal trajectory of very flexible link manipulators in point-to-point motion using a new displacement approach. A new nonlinear finite element model for the dynamic analysis is employed to describe nonlinear modeling for three-dimensional flexible link manipulators, in which both the geometric elastic nonlinearity and the foreshortening effects are considered. In comparison to other large deformation formulations, the motion equations contain constant stiffness matrix because the terms arising from geometric elastic nonlinearity are moved from elastic forces to inertial, reactive and external forces, which are originally nonlinear. This makes the formulation particularly efficient in computational terms and numerically more stable than alternative geometrically nonlinear formulations based on lower-order terms. In this investigation, the computational method to solve the trajectory planning problem is based on the indirect solution of open-loop optimal control problem. The Pontryagin’s minimum principle is used to obtain the optimality conditions, which is lead to a standard form of a two-point boundary value problem. The proposed approach has been implemented and tested on a single-link very flexible arm and optimal paths with minimum effort and minimum vibration are obtained. The results illustrate the power and efficiency of the method to overcome the high nonlinearity nature of the problem.     

Rigid-flexible coupling dynamics of three-dimensional hub-beams system

In the previous research of the coupling dynamics of a hub-beam system, coupling between the rotational motion of hub and the torsion deformation of beam is not taken into account since the system undergoes planar motion. Due to the small longitudinal deformation, coupling between the rotational motion of hub and the longitudinal deformation of beam is also neglected. In this paper, rigid-flexible coupling dynamics is extended to a hub-beams system with three-dimensional large overall motion. Not only coupling between the large overall motion and the bending deformation, but also coupling between the large overall motion and the torsional deformation are taken into account. In case of temperature increase, the longitudinal deformation caused by the thermal expansion is significant, such that coupling between the large overall motion and the longitudinal deformation is also investigated. Combining the characteristics of the hybrid coordinate formulation and the absolute nodal coordinate formulation, the system generalized coordinates include the relative nodal displacement and the slope of each beam element with respect to the body-fixed frame of the hub, and the variables related to the spatial large overall motion of the hub and beams. Based on precise strain-displacement relation, the geometric stiffening effect is taken into account, and the rigid-flexible coupling dynamic equations are derived using velocity variational principle. Finite element method is employed for discretization. Simulation of a hub-beams system is used to show the coupling effect between the large overall motion and the torsional deformation as well as the longitudinal deformation. Furthermore, conservation of energy in case of free motion is shown to verify the formulation.     

Experimental validation of rigid-flexible coupling dynamic formulation for hub–beam system

The aim of this paper is to compare the accuracy of the absolute nodal coordinate formulation and the floating frame of reference formulation for the rigid-flexible coupling dynamics of a three-dimensional Euler–Bernoulli beam by numerical and experimental validation. In the absolute nodal coordinate formulation, based on geometrically exact beam theory and considering the torsion effect, the material curvature of the beam is derived, and then variational equations of motion of a three-dimensional beam are obtained, which consist of three position coordinates, two slope coordinates, and one rotational coordinate. In the floating frame of reference formulation, the displacement of an arbitrary point on the beam is described by the rigid-body motion and a small superimposed deformation displacement. Based on linear elastic theory, the quadratic terms of the axial strain are neglected, and the curvatures are simplified to the first order. Considering both the linear damping and the quadratic air resistance damping, the equations of motion of the multibody system composed of air-bearing test bed and a cantilevered three-dimensional beam are derived based on the principle of virtual work. In order to verify the results of the computer simulation, two experiments are carried out: an experiment of hub–beam system with large deformation and a dynamic stiffening experiment. The comparison of the simulation and experiment results shows that in case of large deformation, the frequency result obtained by the floating frame of reference formulation is lower than that obtained by the experiment. On the contrary, the result obtained by the absolute nodal coordinate formulation agrees well with that obtained by the experiment. It is also shown that the floating frame of reference formulation based on linear elastic theory cannot reveal the dynamic stiffening effect. Finally, the applicability of the floating frame of reference formulation is clarified.     

Modeling Spatial Motion of 3D Deformable Multibody Systems with Nonlinearities

A computational strategy for modeling spatial motion of systems of flexible spatial bodies is presented. A new integral formulation of constraints is used in the context of the floating frame of reference approach. We discuss techniques to linearize the equations of motion both with respect to the kinematical coupling between the deformation and rigid body degrees of freedom and with respect to the geometrical nonlinearities (inclusion of stiffening terms). The plastic behavior of bodies is treated by means of plastic multipliers found as the result of fixed-point type iterations within a time step. The time integration is based on implicit Runge Kutta schemes with arbitrary order and of the RadauIIA type. The numerical results show efficiency of the developed techniques.     

A distributed parameter model for the dynamics of flexible-link robots

In this paper a model is developed for kinematic and dynamic analysis of flexible robots undergoing general three-dimensional motion. For modeling robotic links, distributed mass and flexibility are considered without discretization. Some modeling issues are discussed, and parameters characterizing the real design of a robot are introduced into the analysis. The concept of a fictitious rigid link is presented to consider the rigid body motion of a link separately, and to account for possibly complex link shapes. Based on Jourdain's principle, an alternative formulation is proposed to derive the dynamic equations of flexible robots. The equations of motion are developed and analyzed in detail. The vibrations of links are described by linear, inhomogeneous partial differential equations, with homogeneous, nonlinear, time-dependent boundary conditions. © 1998 John Wiley & Sons, Inc.     

Modeling of flexible-link manipulators with prismatic joints

The axially translating flexible link in flexible manipulators with a prismatic joint can be modeled using the Euler-Bernoulli beam equation together with the convective terms. In general, the method of separation of variables cannot be applied to solve this partial differential equation. In this paper, we present a nondimensional form of the Euler-Bernoulli beam equation using the concept of group velocity and present conditions under which separation of variables and assumed modes method can be used. The use of clamped-mass boundary conditions lead to a time-dependent frequency equation for the translating flexible beam. We present a novel method to solve this time-dependent frequency equation by using a differential form of the frequency equation. We then present a systematic modeling procedure for spatial multi-link flexible manipulators having both revolute and prismatic joints. The assumed mode/Lagrangian formulation of dynamics is employed to derive closed form equations of motion. We show, using a model-based control law, that the closed-loop dynamic response of modal variables become unstable during retraction of a flexible link, compared to the stable dynamic response during extension of the link. Numerical simulation results are presented for a flexible spatial RRP configuration robot arm. We show that the numerical results compare favorably with those obtained by using a finite element-based model.     

Nonlinear generalized equations of motion for multi-link inverted pendulum systems

In this paper, the equations of motion for a general multi-link inverted pendulum system are derived. Assumptions previously employed to simplify such formulation are removed. The pendulum system is more general and includes nonlinear friction terms to suit various engineering applications. The generalized equations are first developed in the absolute coordinate system using Lagrange's technique, then a simple linear transformation is proposed to obtain the set of nonlinear equations in the DevanitHartenberg coordinate system. The equations of motion for double and triple link inverted pendulum systems are given as examples for such dynamics equations.     

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

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

Transient analysis of flexible multi-body systems. Part I: Dynamics of flexible bodies

The formulation for the dynamic analysis of undamped linear structural systems using the finite element method results in two element matrices; the mass and stiffness matrices, that describe the element inertia and stiffness properties. However, these matrices are not sufficient to describe the dynamics of structures that undergo large rigid-body motion. Other element matrices, in addition to the mass and stiffness matrices, are required to account for the inertia coupling between gross motion and elastic deformation. These matrices are time-invariant and can be generated and assembled in the same manner as the mass and stiffness matrices are assembled in linear structural dynamics. An inherent relation between these matrices and the deformable body mean axes exists. This paper is the first of two parts. It presents the two-dimensional and three-dimensional formulation of the system equations of motion of inertia-variant flexible bodies. In particular, Euler parameters are employed to describe the rotations of the body reference in the spatial analysis. In Part II [13], this formulation is applied to the impact analysis of a large-scale constrained flexible aircraft which are modeled as a multi-body system consisting of interconnected rigid and flexible components.     

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.     

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.     

Dynamics modeling for a rigid-flexible coupling system with nonlinear deformation field

In this paper, a moving flexible beam, which incorporates the effect of the geometrically nonlinear kinematics of deformation, is investigated. Considering the second-order coupling terms of deformation in the longitudinal and transverse deflections, the exact nonlinear strain-displacement relations for a beam element are described. The shear strains formulated by the present modeling method in this paper are zero, so it is reasonable to use geometrically nonlinear deformation fields to demonstrate and simplify a flexible beam undergoing large overall motions. Then, considering the coupling terms of deformation in two dimensions, finite element shape functions of a beam element and Lagrange’s equations are employed for deriving the coupling dynamical formulations. The complete expression of the stiffness matrix and all coupling terms are included in the formulations. A model consisting of a rotating planar flexible beam is presented. Then the frequency and dynamical response are studied, and the differences among the zero-order model, first-order coupling model and the new present model are discussed. Numerical examples demonstrate that a ‘stiffening beam’ can be obtained, when more coupling terms of deformation are added to the longitudinal and transverse deformation field. It is shown that the traditional zero-order and first-order coupling models may not provide an exact dynamic model in some cases.     

Flexible Multibody Systems Models Using Composite Materials Components

The use of a multibody methodology to describe the large motion of complex systems that experience structural deformations enables to represent the complete system motion, the relative kinematics between the components involved, the deformation of the structural members and the inertia coupling between the large rigid body motion and the system elastodynamics. In this work, the flexible multibody dynamics formulations of complex models are extended to include elastic components made of composite materials, which may be laminated and anisotropic. The deformation of any structural member must be elastic and linear, when described in a coordinate frame fixed to one or more material points of its domain, regardless of the complexity of its geometry. To achieve the proposed flexible multibody formulation, a finite element model for each flexible body is used. For the beam composite material elements, the sections properties are found using an asymptotic procedure that involves a two-dimensional finite element analysis of their cross-section. The equations of motion of the flexible multibody system are solved using an augmented Lagrangian formulation and the accelerations and velocities are integrated in time using a multi-step multi-order integration algorithm based on the Gear method.     

Solution schemes for the system equations of flexible robots

In this work, a method for generating the dynamic equations of flexible robots with open-chain linkage mechanisms is developed. A general transformation matrix associated with the elastic deformation is introduced. In determining the elastic response, a method of separation of variables and the Galerkin's approach are suggested for the boundary-value problem with time-dependent boundary conditions. Besides the formulation scheme, the present work also studies the difficulty of dealing with the inverse kinematic problem, in which the unknowns involve the rigid-body displacements and the elastic deflections. Finally, the ideas presented here have been implemented in a computer simulation, and the formulation of the boundary-value problem has been employed to obtain the equations of motion of a flexible robot. Simulation results are presented.     

空间柔性机械臂的真-伪混合坐标形式的运动方程

本文从柔性体的真-伪混合坐标形式的拉格朗日方程出发,导出了空间双柔性连杆机械臂系统的真-伪混合坐标形式的动力学方程.由于采用伪坐标和矩阵运算,避免了繁冗的对方位参数的求导过程,使方程结构紧凑,便于建立递推形式,并且在选择不同参数描述整体刚性运动情况下,方程形式保持不变,因而具有一般性.     

Dynamic Modeling and Analysis of a Robot Manipulator Intercepting and Capturing a Moving Object,with the Consideration of Structural Flexibility

In this paper we investigate the dynamics of robotic interception and capture of a moving object. This problem, i.e., the interception and capture of a moving object by a robot, is called dynamic mass capture. The effects of structural flexibility of the robot is taken into consideration. In terms of time, the analysis is divided into three phases: before capture (finite motion), at the vicinity of interception and capture (impulsive motion), and after capture (finite motion). Special attention is paid to the modeling of the second phase when the robot captures the target and it becomes part of the end effector, thus, the system's degrees of freedom suddenly change. To decribe this event, a novel approach is proposed. This is based on the use of a class of impulsive constraints, the so-called inert constraints. Jourdain's principle is employed to derive the dynamic equations for both finite and impulsive motions. Simulation results are presented for two examples: a single flexible link and a two-link manipulator capturing moving objects. In the example of the single link, the results are compared with the observations of an experiment, and good agreement is found between experimental and simulation results.     

Modeling and simulation of structural components in recursive closed-loop multibody systems

Passenger cars, transit buses, railroad vehicles, off-highway trucks, earth moving equipment and construction machinery contain structural and light-fabrications (SALF) components that are prone to excessive vibration due to rough terrains and work-cycle loads’ excitations. SALF components are typically modeled as flexible components in the multibody system allowing the analysts to predict elastic deformation and hence the stress levels under different loading conditions. Including SALF component in the multibody system typically generates closed-kinematic loops. This paper presents an approach for integrating SALF modeling capabilities as a flexible body in a general-purpose multibody dynamics solver that is based on joint-coordinates formulation with the ability to handle closed-kinematic loops. The spatial algebra notation is employed in deriving the spatial multibody dynamics equations of motion. The system kinematic topology matrix is used to project the Cartesian quantities into the joint subspace, leading to a condensed set of nonlinear equations with minimum number of generalized coordinates. The proposed flexible body formulation utilizes the component mode synthesis approach to reduce the large number of finite element degrees of freedom to a small set of generalized modal coordinates. The resulting reduced flexible body model has two main characteristics: the stiffness matrix is constant while the mass matrix depends on the elastic modal coordinates. A consistent set of pre-computed inertia shape integrals are identified and used to update the modal mass matrix at each time step. The implementation of the component mode synthesis approach in a closed-loop recursive multibody formulation is presented. The kinematic equations are modified to include the effect of the flexible body modal elastic coordinates. Also, modified constraint equations that include the effect of flexibility at the joint connections and the necessary details of the Jacobian matrix are presented. Baumgarte stabilization approach is used to stabilize the constraint equations without using iterative schemes. A sample results for flexible body impeded in a closed system will be presented to demonstrate the above mentioned approach.     

Finite element formulation for dynamics of planar flexible multi-beam system

In some previous geometric nonlinear finite element formulations, due to the use of axial displacement, the contribution of all the elements lying between the reference node of zero axial displacement and the element to the foreshortening effect should be taken into account. In this paper, a finite element formulation is proposed based on geometric nonlinear elastic theory and finite element technique. The coupling deformation terms of an arbitrary point only relate to the nodal coordinates of the element at which the point is located. Based on Hamilton principle, dynamic equations of elastic beams undergoing large overall motions are derived. To investigate the effect of coupling deformation terms on system dynamic characters and reduce the dynamic equations, a complete dynamic model and three reduced models of hub-beam are prospected. When the Cartesian deformation coordinates are adopted, the results indicate that the terms related to the coupling deformation in the inertia forces of dynamic equations have small effect on system dynamic behavior and may be neglected, whereas the terms related to coupling deformation in the elastic forces are important for system dynamic behavior and should be considered in dynamic equation. Numerical examples of the rotating beam and flexible beam system are carried out to demonstrate the accuracy and validity of this dynamic model. Furthermore, it is shown that a small number of finite elements are needed to obtain a stable solution using the present coupling finite element formulation.     

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.