An exact nonlinear hybrid-coordinate formulation for flexible multibody systems

The previous low-order approximate nonlinear formulations succeeded in capturing the stiffening terms, but failed in simulation of mechanical systems with large deformation due to the neglect of the high-order deformation terms. In this paper, a new hybrid-coordinate formulation is proposed, which is suitable for flexible multibody systems with large deformation. On the basis of exact strain–displacement relation, equations of motion for flexible multibody system are derived by using virtual work principle. A matrix separation method is put forward to improve the efficiency of the calculation. Agreement of the present results with those obtained by absolute nodal coordinate formulation (ANCF) verifies the correctness of the proposed formulation. Furthermore, the present results are compared with those obtained by use of the linear model and the low-order approximate nonlinear model to show the suitability of the proposed models. The project supported by the National Natural Science Foundation of China (10472066, 50475021).     

Partition method for impact dynamics of flexible multibody systems based on contact constraint

The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.     

多体系统动力学求解与动画输出的并行处理

在介绍了多体系统动力学的建模方法及数值仿真的原理之后，着重介绍如何利用局域网的技术，在多台计算机之间实现数值求解和动画计算与输出的并行处理。用一个柔性索动力学仿真算例，论证了该技术的可行性，为多体系统动力学仿真的并行处理打下了基础。     

Regularization and stability of the constraints in the dynamics of multibody systems

In the analysis of multibody dynamics, we are often required to deal with singularity problems where the constraint Jacobian matrix may become less than full rank at some instantancous configurations. This creates numerical instability which will affect the performance of the mechanical system. A modification procedure of the constraints when they vanish or become linearly dependent is proposed to regularize the dynamics of the system. A distinction between the asymptotic stability due to the representation of the constraints (at the velocity and acceleration level), and the one due to the singularity is discussed in full in this paper. It is shown that Baumgarte technique could be extended to accommodate the representation of the constraints in the neighborhood of singularity. A two link planar manipulator undergoing large motion and passing through a singular configuration is used to illustrate the proposed stability technique.     

Screw-matrix method in dynamics of multibody systems

In the present paper the concept of screw in classical mechanics is expressed in matrix form, in order to formulate the dynamical equations of the multibody systems. The mentioned method can retain the advantages of the screw theory and avoid the shortcomings of the dual number notation. Combining the screw-matrix method with the tool of graph theory in Roberson/Wittenberg formalism. We can expand the application of the screw theory to the general case of multibody systems. For a tree system, the dynamical equations for eachj-th subsystem, composed of all the outboard bodies connected byj-th joint can be formulated without the constraint reaction forces in the joints. For a nontree system, the dynamical equations of subsystems and the kinematical consistency conditions of the joints can be derived using the loop matrix. The whole process of calculation is unified in matrix form. A three-segment manipulator is discussed as an example. This work is supported by the National Natural Science Fund.     

Recursive dynamics simulator (ReDySim): A multibody dynamics solver

Recursive formulations have significantly helped in achieving real-time computations and model-based control laws. The recursive dynamics simulator (ReDySim) is a MATLAB-based recursive solver for dynamic analysis of multibody systems. ReDySim delves upon the decoupled natural orthogonal complement approach originally developed for serial-chain manipulators. In comparison to the commercially available software, dynamic analyses in ReDySim can be performed without creating solid model. The input parameters are specified in MATLAB environment. ReDySim has the capability to incorporate any control algorithm with utmost ease. In this work, the capabilities of ReDySim for solving open-loop and closed-loop systems are shown by examples of robotic gripper, KUKA KR5 industrial manipulator and four-bar mechanism. ReDySim can be downloaded for free from http://www.redysim.co.nr and can be used almost instantly.     

Fast simulation of flexible multibody dynamics with electric-hydraulic drive system

Hardware in the loop simulation (HILS) has been investigated in the field of the multibody dynamics (MBD), which combined the MBD simulation with the actual mechanical system. The fast simulation is necessary for the HILS system in order to require the real time simulation. This paper presents a fast simulation technique using the domain decomposition method with the iteration in the flexible multibody system in which flexible linkage system and electro-hydraulic drive system are coupled with each other.     

Finite element modeling of contact conditions in multibody system dynamics

In this paper a new method is developed for the dynamic analysis of contact conditions in flexible multibody systems undergoing a rolling type of motion. The relative motion between the two contacting bodies is treated as a constraint condition describing their kinematic and geometric relations. Equations of motion of the system are presented in a matrix form making use of Kane's equations and finite element method. The method developed has been implemented in a general purpose program called DARS and applied to the simulation and analysis of a rotating wheel on a track. Both the bodies are assumed flexible and discretized using a three dimensional 8-noded isoparametric elements. The time variant constraint conditions are imposed on the nodal points located at the peripheral surfaces of the bodies under consideration. The simulation is carried out under two different boundary conditions describing the support of the track. The subsequent constraint forces associated with the generalized coordinates of the system are computed and plotted. The effects of friction are also discussed.     

Relationship between multibody dynamics computation and nonlinear vibration theory

This paper presents the ground-work of implementing the multibody dynamics codes to analyzing nonlinear coupled oscillators. The recent developments of the multibody dynamics have resulted in several computer codes that can handle large systems of differential and algebraic equations (DAE). However, these codes cannot be used in their current format without appropriate modifications. According to multibody dynamics theory, the differential equations of motion are linear in the acceleration, and the constraints are appended into the equations of motion through Lagrange's multipliers. This formulation should be able to predict the nonlinear phenomena established by the nonlinear vibration theory. This can be achieved only if the constraint algebraic equations are modified to include all the system kinematic nonlinearities. This modification is accomplished by considering secondary nonlinear displacements which are ignored in all current codes. The resulting set of DAE are solved by the Gear stiff integrator. The study also introduced the concept of constrained flexibility and uses an instantaneous energy checking function to improve integration accuracy in the numerical scheme. The general energy balance is a single scalar equation containing all the energy component contributions. The DAE solution is then compared with the solution predicted by the nonlinear vibration theory. It also establishes new foundation for the use of multibody dynamics codes in nonlinear vibration problems. It is found that the simulation CPU time is much longer than the simulation of the original equations of the system.     

Impact dynamics of multibody systems with frictional contact using joint coordinates and canonical equations of motion

This paper presents a methodology in computational dynamics for the analysis of mechanical systems that undergo intermittent motion. A canonical form of the equations of motion is derived with a minimal set of coordinates. These equations are used in a procedure for balancing the momenta of the system over the period of impact, calculating the jump in the body momentum, velocity discontinuities and rebounds. The effect of dry friction is discussed and a contact law is proposed. The present formulation is extended to open and closed-loop mechanical systems where the jumps in the constraints' momenta are also solved. The application of this methodology is illustrated with the study of impact of open-loop and closed-loop examples.     

A nonlinear finite-element approach for kineto-static analysis of multibody systems

Methods that treat rigid/flexible multibody systems undergoing large motion as well as deformations are often accompanied with inefficiencies and instabilities in the numerical solution due to the large number of state variables, differences in the magnitudes of the rigid and flexible body coordinates, and the time dependencies of the mass and stiffness matrices. The kineto-static methodology of this paper treats a multibody mechanical system to consist of two collections of bulky (rigid) bodies and relatively flexible ones. A mixed boundary condition nonlinear finite element problem is then formulated at each time step whose known quantities are the displacements of the nodes at the boundary of rigid and flexible bodies and its unknowns are the deformed shape of the entire structure and the loads (forces and moments) at the boundary. Partitioning techniques are used to solve the systems of equations for the unknowns, and the numerical solution of the rigid multibody system governing equations of motion is carried out. The methodology is very much suitable in modelling and predicting the impact responses of multibody system since both nonlinear and large gross motion as well as deformations are encountered. Therefore, it has been adopted for the studies of the dynamic responses of ground vehicle or aircraft occupants in different crash scenarios. The kineto-static methodology is used to determine the large motion of the rigid segments of the occupant such as the limbs and the small deformations of the flexible bodies such as the spinal column. One of the most dangerous modes of injury is the amount of compressive load that the spine experiences. Based on the developed method, a mathematical model of the occupant with a nonlinear finite element model of the lumbar spine is developed for a Hybrid II (Part 572) anthropomorphic test dummy. The lumbar spine model is then incorporated into a gross motion occupant model. The analytical results are correlated with the experimental results from the impact sled test of the dummy/seat/restraint system. With this extended occupant model containing the lumbar spine, the gross motion of occupant segments, including displacements, velocities and accelerations as well as spinal axial loads, bending moments, shear forces, internal forces, nodal forces, and deformation time histories are evaluated. This detailed information helps in assessing the level of spinal injury, determining mechanisms of spinal injury, and designing better occupant safety devices.     

约束多体系统的基于离散零空间的隐式龙格库塔法

将离散零空间理论应用于多体系统动力学方程的数值计算,可降低多体系统动力学方程的维数。通过给出离散零空间理论与IRK法相结合的一般数学框架,提出了多体系统动力学的基于离散零空间理论的IRK法。数值算例表明：该算法可获得较满意的数值结果,约束违约程度很小,三种积分算法算例的范数均在10^-16之内。     

Unilateral problems of dynamics

Summary Contact processes may be described by local discretizations, by rigid representation or by mixed methods incorporating both ideas. A rigid body approach is proposed for the dynamics of mechanical systems, achieving good results also for multiple-contact problems. Contacts in multibody systems are mainly considered, with the corresponding contact constraints varying with time, thus generating structure-variant systems. The equations of motion for dynamical systems with such unilateral behavior are discussed, solution methods and applications are presented. Received 3 March 1999; accepted for publication 5 May 1999     

平面柔性多体系统正碰撞动力学建模理论研究

针对目前柔性多体系统碰撞动力学建模方法存在的不足,对影响碰撞动力学仿真的主要因素如柔性体建模和碰撞初始条件进行分析,建立起基于变约束的柔性体碰撞动力学方程。首先,为了解决子结构法在处理碰撞界面搜索时面临的难题,引入多体系统柔性体有限元描述方法,推导出凸形柔性体接触点间法向位移约束的二阶导数形式。其次,从碰撞引起的接触界面速度不连续机理出发,结合连续介质力学间断面理论,给出碰撞瞬时由物体本身物理性质决定的接触位置处速度跳跃公式。最后对两弹性圆盘低速碰撞问题进行数值仿真。结果表明本文提出的改进方法符合力学基本原理,仿真结果满足收敛性要求。     

An automatic constraint violation stabilization method for differential/ algebraic equations of motion in multibody system dynamics

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｄｙｎａｍｉｃｅｑｕａｔｉｏｎｓｏｆｍｏｔｉｏｎｏｆｍｕｌｔｉｂｏｄｙｓｙｓｔｅｍｓｗｉｔｈｃｏｎｓｔｒａｉｎｔｓａｒｅｔｈｅｆｏｌｌｏｗｉｎｇｄｉｆｆｅｒｅｎｔｉａｌ／ａｌｇｅｂｒａｉｃｅｑｕａｔｉｏｎｓ,ｉ．ｅ．,Ｅ...     

刚体单元及其在多体系统动力学中的应用

多体系统动力学分析软件要求人工输入形状复杂物体的质量、质心位置和转动惯量,而实际上这些参量并不容易获得。本文探索了一种以组成物体的刚体单元为基本要素的新方法,并结合实际需要具体构造了刚性四面体和刚性梁单元。以刚体单元为基础并内嵌网格剖分模块的分析软件能够自动获得这些参数,从而具备处理任何复杂系统的能力。仿真结果的对比分...     

Ride comfort evaluation for road vehicle based on rigid-flexible coupling multibody dynamics

In the present research two different whole vehicle multibody models are established respectively, including rigid and rigid-flexible coupling multibody vehicle models. The former is all composed by rigid bodies while in the later model, the flexible rear suspension is built based on the finite element method (FEM) and mode superposition method, in which the deformations of the components are considered. The ride simulations with different speeds are carried out on a 3D digitalized road, and the weighted root mean square (RMS) of accelerations on the seat surface, backrest and at the feet are calculated. The comparison between the responses of the rigid and rigid-flexible coupling multibody models shows that the flexibility of the vehicle parts significantly affects the accelerations at each position, and it is necessary to take the flexibility effects into account for the assessment of ride comfort.     

Large deformations of nonlinear viscoelastic and multi-responsive beams

This study presents analyses of deformations in nonlinear viscoelastic beams that experience large displacements and rotations due to mechanical, thermal, and electrical stimuli. The studied beams are relatively thin so that the effect of the transverse shear deformation is neglected, and the stretch along the transverse axis of the beams is also ignored. It is assumed that the plane that is perpendicular to the longitudinal axis of the undeformed beam remains plane during the deformations. The nonlinear kinematics of the finite strain beam theory presented by Reissner [27] is adopted, and a nonlinear viscoelastic constitutive relation based on a quasi-linear viscoelastic (QLV) model is considered for the beams. Deformation in beams due to mechanical, thermal, and electric field inputs are incorporated through the use of time integral functions, by separating the time-dependent function and nonlinear measures of field variables. The nonlinear measures are formulated by including higher order terms of the field variables, i.e. strain, temperature, and electric field. Responses of beams under mechanical, thermal, and electrical stimuli are illustrated and the effects of nonlinear constitutive relations on the overall deformations of the beams are highlighted.