Numerical Analysis of Flexible Multibody Systems

Flexible multibody systems featurea coupled mathematical model with ODEs or DAEs governing the gross motion and PDEs describing the elastic deformation of particular bodies.Frequently, semidiscretization of elastic members leads to a stiff mechanical system with widely different time scales. Thispaper investigates the behavior of numerical time integration methodsat flexible multibody systems and gives some recommendations.Simulation results for a slider crank mechanism and a truck modelillustrate what can go wrong and how implicit methods like RADAU5 can beeffectively applied.     

First order sensitivity analysis of flexible multibody systems using absolute nodal coordinate formulation

Design sensitivity analysis of flexible multibody systems is important in optimizing the performance of mechanical systems. The choice of coordinates to describe the motion of multibody systems has a great influence on the efficiency and accuracy of both the dynamic and sensitivity analysis. In the flexible multibody system dynamics, both the floating frame of reference formulation (FFRF) and absolute nodal coordinate formulation (ANCF) are frequently utilized to describe flexibility, however, only the former has been used in design sensitivity analysis. In this article, ANCF, which has been recently developed and focuses on modeling of beams and plates in large deformation problems, is extended into design sensitivity analysis of flexible multibody systems. The Motion equations of a constrained flexible multibody system are expressed as a set of index-3 differential algebraic equations (DAEs), in which the element elastic forces are defined using nonlinear strain-displacement relations. Both the direct differentiation method and adjoint variable method are performed to do sensitivity analysis and the related dynamic and sensitivity equations are integrated with HHT-I3 algorithm. In this paper, a new method to deduce system sensitivity equations is proposed. With this approach, the system sensitivity equations are constructed by assembling the element sensitivity equations with the help of invariant matrices, which results in the advantage that the complex symbolic differentiation of the dynamic equations is avoided when the flexible multibody system model is changed. Besides that, the dynamic and sensitivity equations formed with the proposed method can be efficiently integrated using HHT-I3 method, which makes the efficiency of the direct differentiation method comparable to that of the adjoint variable method when the number of design variables is not extremely large. All these improvements greatly enhance the application value of the direct differentiation method in the engineering optimization of the ANCF-based flexible multibody systems.     

A new version of transfer matrix method for multibody systems

In order to avoid the global dynamics equations and increase the computational efficiency for multibody system dynamics (MSD), the transfer matrix method of multibody system (MSTMM) has been developed and applied very widely in research and engineering in recent 20 years. It differs from ordinary methods in multibody system dynamics with respect to the feature that there is no need for a global dynamics equation, and it uses low-order matrices for high computational efficiency. For linear systems, MSTMM is exact even if continuous elements like beams are involved. The discrete time MSTMM, however, has to use local linearization. In order to release the method from such approximations, a new version of MSTMM is presented in this paper where translational and angular accelerations, on the one hand, and internal forces and moments, on the other hand, are used as state variables. Already linear relationships among these quantities are utilized, which results in new element transfer matrices and algorithms making the study of multibody systems as simple as the study of single bodies. The proposed approach also allows combining MSTMM with any general numerical integration procedure. Some numerical examples of MSD are given to demonstrate the proposed method.     

Analysis of Nonlinear Multibody Systems with Elastic Couplings

This paper is concerned with the dynamic analysis of nonlinear multibody systems involving elastic members made of laminated, anisotropic composite materials. The analysis methodology can be viewed as a three-step procedure. First, the sectional properties of beams made of composite materials are determined based on an asymptotic procedure that involves a two-dimensional finite element analysis of the cross-section. Second, the dynamic response of nonlinear, flexible multibody systems is simulated within the framework of energy-preserving and energy-decaying time-integration schemes that provide unconditional stability for nonlinear systems. Finally, local three-dimensional stresses in the beams are recovered, based on the stress resultants predicted in the previous step. Numerical examples are presented and focus on the behavior of multibody systems involving members with elastic couplings.     

Transfer matrix method for linear multibody system

A new method for linear hybrid multibody system dynamics is proposed in this paper. This method, named as transfer matrix method of linear multibody system (MSTMM), expands the advantages of the traditional transfer matrix method (TMM). The concepts of augmented eigenvector and equation of motion of linear hybrid multibody system are presented at first to find the orthogonality and to analyze the responses of the hybrid multibody system using modal method. If using this method, the global dynamics equation is not needed in the study of linear hybrid multibody system dynamics. The MSTMM has a small size of matrix and higher computational speed, and can be applied to linear multi-rigid-body system dynamics, linear multi-flexible-body system dynamics and linear hybrid multibody system dynamics. This method is simple, straightforward, practical, and provides a powerful tool for the study on linear hybrid multibody system dynamics. This method can be used in the following: (1) Solve the eigenvalue problem of linear hybrid multibody systems. (2) Obtain the orthogonality of eigenvectors of linear hybrid multibody systems. (3) Realize the accurate analysis of the dynamics response of linear hybrid multibody systems. (4) Find the connected parameters between bodies used in the computation of linear hybrid multibody systems. A practical engineering system is taken as an example of linear multi-rigid-flexible-body system, the dynamics model, the transfer equations and transfer matrices of various bodies and hinges; the overall transfer equation and overall transfer matrix of the system are developed. Numerical example shows that the results of the vibration characteristics and the response of the hybrid multibody system received by MSTMM and by experiment have good agreements. These validate the proposed method.     

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.     

Impacts in multibody systems: modeling and experiments

This paper outlines a novel approach to the modeling and analysis of impact involving multibody systems. This approach is based on an analysis of energy absorption and restitution during impact, using a decomposition of the kinetic energy, which decouples the parts associated with the spaces of admissible and constrained motions of the underlying unilateral constraints. Such a decomposition turns out to be useful in the analysis of energy dissipation during impact, and leads to a generalized definition of the energetic coefficient of restitution, which targets particularly collisions in multibody systems. The applicability of the approach reported is investigated by conducting an experimental study on a robotic testbed. It is shown that impact between multibody systems is considerably affected not only by the local dynamics characteristics of the interacting bodies, but also the configuration of the whole multibody system. The results reported here show that our decomposition can offer a sound characterization of impact in several problems of multibody systems.     

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.     

Nonholonomic Multibody Dynamics

The paper deals with the nonholonomic multibody system dynamics from apoint of view which is caused by some actual applications in high-tecareas like high-speed train technology or biomechanics of somedisciplines in high-performance sports. Obviously, looking at suchproblems, there are very close connections between classical analyticaldynamics, differential geometry and modern control theory. But theseconnections cannot be used to get new composed results in solvingcomplicated problems of multibody system dynamics because correspondingsoftware tools are not enough in tune with each other. This paper willgive some ideas for developing a unified basis for modeling, simulationand control of nonholonomic multibody systems.First, a derivative-free approach for generating Lagrangian motionequations of multibody systems with kinematical tree structure as wellas for constrained multibody systems is given. This has been done byusing differential-geometric concepts in a Riemannian space. Secondly,the well-known theorem of Frobenius is considered with respect to itsclassical interpretation by the so-called object of nonholonomy as wellas by its modern interpretation in the nonlinear control theory usingLie-brackets. The ideas are illustrated by the classical rollingcondition and edge condition on double-curved surfaces. Specialnumerical problems in simulation of multibody systems subject toadditional kinematic constraints are discussed. Finally threeapplications are given.