接触碰撞动力学的多变量选取方法

在柔性多体的接触碰撞动力学问题中,多变量方法基于附加约束的接触模型,将柔性体的变形用不同变量来描述：接触局部区域的变形用有限元节点坐标描述,非接触局部区域的变形用模态坐标描述,兼顾了计算精度和效率. 将该方法推广到三维空间碰撞问题,对两柔性杆纵向碰撞过程进行动力学仿真,数值结果与实验结果吻合良好,验证了该方法的有效性. 针对柔性体各自区域的变量如何选取的问题,研究了节点取法、模态阶数以及材料参数对计算结果精度的影响,寻找到合理的多变量选取方法,保证精度的同时使自由度得到最大程度的缩减.     

Investigation on the choice of boundary conditions and shape functions for flexible multi-body system

The objective of this investigation is to examine the correctness and efficiency of the choice of boundary conditions when using assumed mode approach to simulate flexible multi-body systems. The displacement field due to deformation is approximated by the Rayleigh-Ritz assumed modes in floating frame of reference (FFR) formulation. The deformations obtained by the absolute nodal coordinate (ANC) formulation which are transformed by two sets of reference coordinates are introduced as a criterion to verify the accuracy of the simulation results by using the FFR formulation. The relationship between the deformations obtained from different boundary conditions is revealed. Numerical simulation examples demonstrate that the assumed modes with cantilevered-free, simply-supported and freefree boundary conditions without inclusion of rigid body modes are suitable for simulation of flexible multi-body system with large over all motion, and the same physical deformation can be obtained using those mode functions, differ only by a coordinate transformation. It is also shown that when using mode shapes with statically indeterminate boundary conditions, significant error may occur. Furthermore, the slider crank mechanism with rigid crank is accurate enough for investigating boundary condition problem of flexible multi-body system, which cost significant less simulating time.     

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.     

柔性多体碰撞问题的多变量方法

针对柔性多体系统碰撞问题的特点, 提出了柔性多体系统碰撞问题的多变量方法. 在未碰撞阶段和分离后的阶段, 多变量方法考虑了物体大范围运动与变形的耦合.在碰撞阶段, 考虑非碰撞区域的变形, 对非碰撞区域用浮动参考系方法求解, 对发生碰撞的局部区域用非线性有限元方法求解. 非碰撞区域采用浮动基位形坐标和柔性体相对浮动基的变形模态坐标描述, 大大减缩了变量的维数,提高了计算效率. 碰撞局部区域采用非线性有限元节点坐标描述, 可以得到碰撞局部区域高精度的应力、应变响应, 而且可以反映碰撞局部区域的大变形、塑性等非线性响应. 该方法既能体现碰撞对系统大范围运动和长期动力学仿真的影响, 又能精确反映出碰撞发生的微小时间段内物体碰撞局部区域的应力、应变响应, 表现出碰撞过程中物体接触区域的演化历程.介绍了两杆正碰撞与铅垂面内双摆撞击自由圆盘的算例, 设计了两柔性杆正碰撞实验,多变量方法的数值仿真结果与实验测量结果很好地吻合.      

多体系统接触碰撞问题的牛顿-欧拉线性互补方法

基于非光滑动力学方法的多体系统接触碰撞分析是目前多体系统动力学的研究热点.本文采用牛顿-欧拉方法建立多体系统接触、碰撞问题的动力学模型,给出一种牛顿-欧拉型线性互补公式.该建模方法与目前一般采用的拉格朗日建模方法的不同之处是约束条件中除了库仑摩擦、单边约束之外还含有光滑等式约束.在建立系统动力学模型时,首先解除摩擦约束和单边约束得到原系统对应的基本系统.牛顿-欧拉方法采用最大数目坐标建立基本系统的动力学方程,由于坐标不相互独立,因此基本系统中带有等式约束,其数学模型为一组微分代数方程.借助约束雅可比矩阵,在基本系统微分代数方程中添加摩擦接触和单边约束对应的拉氏乘子,就可以得到系统全局运动的具有变拓扑结构特征的动力学方程,再结合非光滑约束互补条件便可构成完备的系统动力学模型.完备的动力学模型由动力学微分方程以及等式约束和不等式约束组成.线性互补公式采用分块矩阵形式进行推导,简化了推导过程.数值计算采用基于线性互补的时间步进算法.时间步进算法是目前流行的非光滑数值算法,其突出特点是可以免去数值积分中繁琐的事件检测过程,而数值积分过程中通过对线性互补问题的求解可以确定系统的触-离状态.通过对典型的曲柄滑块间隙机构进行数值分析,验证本文方法的有效性.     

多体系统动力学方程违约修正的数值计算方法

多体系统动力学方程为微分代数方程,一般将其转化成常微分方程组进行数值计算,在数值积分的过程中约束方程的违约会逐渐增大.本文对具有完整、定常约束的多体系统,在修改的带乘子Lagrange正则形式的方程的基础上,根据Baumgarte提出的违约修正的方法,给出了一种多体系统微分代数方程违约修正法和系统的动力学方程的矩阵表达式.通过对曲柄-滑块机构的数值仿真,计算结果表明本文给出的方法在计算精度和计算效率上好于Baumgarte提出的两种违约修正的方法.     

Numerical method for dynamics of multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints

Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints is presented. In particular, a wheeled multi-body system is considered. Here,the state transition of stick-slip between wheel and ground is transformed into a nonlinear complementarity problem(NCP). An iterative algorithm for solving the NCP is then presented using an event-driven method. Dynamical equations of the multi-body system with holonomic and nonholonomic constraints are given using Routh equations and a constraint stabilization method. Finally, an example is used to test the proposed numerical method. The results show some dynamical behaviors of the wheeled multi-body system and its constraint stabilization effects.     

Partition method and experimental validation for impact dynamics of flexible multibody system

The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one of the main difficulties in many engineering applications. The numerical approaches being used widely in impact analysis are mainly from two fields: multibody system dynamics (MBS) and computational solid mechanics (CSM). Approaches based on MBS provide a more efficient yet less accurate analysis of the contact/impact problems, while approaches based on CSM are well suited for particularly high accuracy needs, yet require very high computational effort. To bridge the gap between accuracy and efficiency in the dynamic simulation of a flexible multibody system with contacts/impacts, a partition method is presented considering that the contact body is divided into two parts, an impact region and a non-impact region. The impact region is modeled using the finite element method to guarantee the local accuracy, while the non-impact region is modeled using the modal reduction approach to raise the global efficiency. A three-dimensional rod-plate impact experiment is designed and performed to validate the numerical results. The principle for how to partition the contact bodies is proposed: the maximum radius of the impact region can be estimated by an analytical method, and the modal truncation orders of the non-impact region can be estimated by the highest frequency of the signal measured. The simulation results using the presented method are in good agreement with the experimental results. It shows that this method is an effective formulation considering both accuracy and efficiency. Moreover, a more complicated multibody impact problem of a crank slider mechanism is investigated to strengthen this conclusion.     

An Augmented Formulation for Mechanical Systems with Non-Generalized Coordinates: Application to Rigid Body Contact Problems

In computational multibody algorithms, the kinematic constraintequations that describe mechanical joints and specified motiontrajectories must be satisfied at the position, velocity andacceleration levels. For most commonly used constraint equations, onlyfirst and second partial derivatives of position vectors with respect tothe generalized coordinates are required in order to define theconstraint Jacobian matrix and the first and second derivatives of theconstraints with respect to time. When the kinematic and dynamicequations of the multibody systems are formulated in terms of a mixedset of generalized and non-generalized coordinates, higher partialderivatives with respect to these non-generalized coordinates arerequired, and the neglect of these derivatives can lead to significanterrors. In this paper, the implementation of a contact model in generalmultibody algorithms is presented as an example of mechanical systemswith non-generalized coordinates. The kinematic equations that describethe contact between two surfaces of two bodies in the multibody systemare formulated in terms of the system generalized coordinates and thesurface parameters. Each contact surface is defined using twoindependent parameters that completely define the tangent and normalvectors at an arbitrary point on the body surface. In the contact modeldeveloped in this study, the points of contact are searched for on lineduring the dynamic simulation by solving the nonlinear differential andalgebraic equations of the constrained multibody system. It isdemonstrated in this paper that in the case of a point contact andregular surfaces, there is only one independent generalized contactconstraint force despite the fact that five constraint equations areused to enforce the contact conditions.     

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.     

A chassis-based kinematic formulation for flexible multi-body vehicle dynamics

This research proposed an efficient implicit integration method for the real-time simulation of flexible multi-body vehicle dynamics models. The equations of motion for the flexible bodies were formulated with respect to the moving chassis-body reference frame instead of the fixed inertial reference frame. The proposed approach does not require evaluation of system Jacobian and its LU-decomposition in time loof of simulation. This is one of the key aspects that enable high computational efficiency of the proposed method. The numerical simulation results of the proposed approach were matched up with those of the conventional approach but the computation time can be reduced by applying the proposed method. The joint constraint and generalized force equations are the same as the equations for rigid vehicle dynamics models because the joints and forces between flexible bodies are connected by the RBE (rigid body element). On the various driving conditions, the numerical simulation results show that the proposed approach yields almost exact solutions compared to the conventional approach. And the proposed approach spends only 22.9% of conventional approach on computation time under CPU 3.2 GHz personal computer.     

A general method for impact dynamic analysis of a planar multi-body system with a rolling ball bearing joint

Impact affects the dynamic characteristics of mechanical multi-body systems and damages those rotating parts, such as the joint rolling element bearings, which are high-precision, defect intolerant components. Based on multi-body dynamic theory, Hertzian contact theory, and a continuous contact model, this study proposed a modelling method that can describe the dynamic behaviour of planar mechanical multi-body systems containing a rolling ball bearing joint under impact. In this method, the rigid bodies and bearing joint were connected according to their joint force constraints; the impact constraint between the multi-body system and the target rigid body was constructed using a continuous contact force model. Based on this method, the reflection relationship between the external impacts of the mechanical multi-body system and the variation law governing the dynamic load on the rolling bearing joint were revealed. Subsequently, an impact multi-body system, which was composed of a sliding–crank mechanism containing a rolling ball bearing joint and the target rigid body with an elastic support, was analysed to explore the dynamic response of such a complex discontinuous dynamic system andthe relevant relationship governing the dynamic load on the rolling bearing joint. In addition, a multi-body dynamic simulation software was used to build a virtual prototype of the impact slider–crank system. Compared with the theoretical model, the prototype had an additional deep groove ball bearing. That is to say, the prototype model took account of the specific geometric structural characteristics and the complex contact relationship of the inner and outer races, rolling balls, and bearing cage. Finally, the effectiveness of the theoretical method proposed in this study was verified by comparative analysis of the results. The results suggested that the external impact of a mechanical multi-body system was prone to induce sudden changes in the equivalent reaction force on its bearing joint and the dynamic load carried on its rolling balls. This study provided an effective method for exploring the distribution characteristics of dynamic loads on rolling ball bearing joints under working impact load conditions. Moreover, it offered support for the parameter optimisation of geometric structure, performance evaluation, and dynamic design of the rolling ball bearings.     

平面柔性多体碰撞阶段附加接触约束方法

含碰撞现象的柔性多体动力学涵盖碰撞前、碰撞中和碰撞后3个不同拓扑的过程, 针对这个特点, 采用附加接触约束方法可将上述3个过程的动力学方程形式统一. 引入碰撞检测机制, 实现拓扑模型的自动变化. 与传统的罚函数方法不同, 以物体间互不嵌入的思想建立附加接触约束, 计算碰撞期间的动力学过程. 通过对发生多次碰撞问题的仿真,从兼顾计算精度和计算效率的角度, 对提出的附加接触约束方法和LS-DYNA罚函数方法进行了对比, 改进了LS-DYNA罚函数方法存在的不足.      

作大范围回转运动柔性梁斜碰撞动力学研究

为正确估计由于碰撞引起的多柔体系统动力学特性的变化,针对作大范围回转运动的柔性梁与一固定斜面发生斜碰撞的情况,在考虑刚柔耦合效应的多柔体系统动力学建模理论的基础上,利用假设模态法建立起重力场作用下的柔性梁一致线性化动力法向碰撞过程中系统的动力行为。基于Ｈｅｒｔｚ接触理论和非线性阻尼项建立法向碰撞接触模型,基于线性切向接触刚度建立柔性梁切向碰撞接触模型,提出的数值算法保证了计算结果的合理性,给出的仿     

多柔体系统碰撞动力学研究综述

多柔体系统碰撞动力学研究具有重要的研究价值和工程实际意义,本文针对多柔体系统碰撞动力学研究中的几个基本问题进行了全面的分析和评述,其中包括多柔体系统动力学方程的描述、碰撞模型的建立、铰接间隙引起的碰撞问题、数值算法、实验研究、控制等几个方面,并根据目前的发展现状和研究中存在的问题,指出了今后多柔体系统动力学碰撞研究中的发展方向      

基于LuGre摩擦模型的接触约束法旋转柔性梁斜碰撞研究

基于接触约束法和LuGre摩擦模型对在重力场作用下作大范围旋转运动的柔性梁系统和斜坡发生含摩擦斜碰撞的动力学问题进行研究.首先运用刚柔耦合的多体系统动力学理论对大范围运动的柔性梁进行离散化和动力学建模,在碰撞时采用冲量动量法求出跳跃速度,其次在法向上引入接触约束求解出碰撞力,在切向上采用LuGre摩擦模型分两种方式求解...     

Dynamic Bifurcation of Multibody Systems

In this paper, the process of loss of stability of multibody systems and structures is analyzed. A novel approach is presented and applied to the statically loaded spatial systems for the analysis of a dynamic response of systems imposed on impact, high velocity compulsive motion, or percussive forces. The analysis is based on the solution of the dynamic equations and eigenvalue problem of systems, and of the resultant motion simulation. The flexible systems are discretized using the finite element method. The dynamic equations are derived with respect to the relative coordinates of the finite elements. Large flexible deflections due to a loss of stability are simulated. The initial forms of the possible deformations are defined by the computed eigenvectors solving the eigenvalue problem for the system stiffness matrix. The critical forces and system deflections are then analyzed. Examples of bifurcation of beam and beam structure imposed on compulsive motion, percussive forces, and impact are presented.