含摩擦滑移铰及驱动约束多刚体系统数值算法

研究了具有驱动约束及非光滑滑移铰多体系统动力学方程的建模与数值计算方法．将驱动约束视为非定常约束，非光滑滑移铰视为双边定常约束，滑移铰的摩擦模型采用库仑摩擦模型；应用第一类Lagrange方程建立系统的动力学方程，应用距离函数建立滑移铰的约束方程；将线性互补方法和Baumgarte约束稳定化方法引入，以解决滑移铰法向约束力的计算以及约束方程违约问题．最后应用曲柄摇杆机构作为算例，说明该方法的有效性．     

含干摩擦多体系统Lagrange方程的数值算法

利用第一类Lagrange方程建立了固定约束面含干摩擦的多体系统动力学方程,将摩擦力的广义力用矩阵形式描述．利用增广法。将微分-代数方程转换成常微分方程,并用矩阵形式给出,提高了计算效率．最后用算例说明该方法的有效性．     

基于微分/代数方程的多体系统动力学设计灵敏度分析的伴随变量方法

基于多体系统动力学微分/代数方程数学模型和通用积分形式的目标函数,建立了多体系统动力学设计灵敏度分析的伴随变量方法,避免了复杂的设计灵敏度计算,对于设计变量较多的多体系统灵敏度分析具有较高的计算效率.文中给出了通用公式以及具体的计算过程和验证方法,并将目标函数及其导数积分形式的计算转化为微分方程的初值问题,进一步提高了计算效率和精度.文末通过一曲柄-滑块机构算例对算法的有效性进行了验证.     

A finite element penalty method for the linearized viscoelastic Oldroyd fluid motion equations

In this paper, a fully discrete finite element penalty method is considered for the two-dimensional linearized viscoelastic fluid motion equations, arising from the Oldroyd model for the non-Newton fluid flows. With the finite element method for the spatial discretization and the backward Euler scheme for the temporal discretization, the velocity and pressure are decoupled in this method, which leads to a large reduction of the computational scale. Under some realistic assumptions, the unconditional stability of the fully discrete scheme is proved. Moreover, the optimal error estimates are obtained, which are better than the existing results. Finally, some numerical results are given to verify the theoretical analysis. The difference between the motion of the Newton and non-Newton fluid is also observed.     

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.     

具有奇异位置的多体系统动力学方程的改进算法

多体系统进行数值仿真时,很多选择了微分代数混合方程作为多体系统动力学数学模型.本文在现有的约束稳定化理论基础上,提出了针对具有奇异位置的多体系统动力学方程的改进算法.算法通过修正速度违约和控制稳定项,讨论了具有奇异位置的微分代数混合方程的数值仿真问题并给出了稳定项中相关系数的建议值,从而有效克服了求解混合方程时因为构型奇异给计算造成的困难.算例分别采用改进算法与ADAMS软件进行仿真,计算结果的比较表明了改进算法的有效性.本文给出的基于能量守恒的能量差曲线也证明了改进算法的有效性.     

A parallel simulator for multibody systems based on group equations

The Journal of Supercomputing - Multibody systems consist of a set of components connected through some joints, where the movement of the system is determined by those of its components. Their...     

Motion safety and constraints compatibility for multibody robots

In this paper we propose a methodology to ensure safe behaviors of multibody robots in reactive control frameworks. The permanent satisfaction of constraints being insufficient to ensure safety, this approach focuses on the constraints expression: the compatibility between these constraints is studied, and safe alternatives are ensured when compatibility cannot be established. Case studies involving obstacles, joint position, velocity and acceleration limits illustrates the approach. A particular method is developed to take full advantage of a smooth state of the art avoidance techniques (Faverjon and Tournassoud in Proceedings of the 1987 IEEE international conference on robotics and automation, pp. 1152–1159, 1987) while maintaining safety. Experiments involving a 6-DOF manipulator operating in a cluttered environment illustrate the reliability of the approach and validate the expected performances.     

基于隐式微分/代数方程的多体系统动力学设计灵敏度分析方法

基于一般性的积分型目标函数、隐式相容初始条件及终止时刻表达式,系统建立了含设计参数的用隐式微分／代数方程表达的多体系统动力学设计灵敏度分析的直接微分方法和伴随变量方法．为降低目标函数及其对设计变量导数的计算复杂性。将其积分形式的计算转化为微分形式．所得到的结果可方便地应用于高效的间接最优化设计方法．最后通过采用绝对坐标建模的平面两连杆机械臂模型对该方法进行了验证．     

一种降低OFDM系统峰均比新方法

文章简要分析了正交频分复用(OFDM)系统中峰均比(PAPR)的成因,并提出了一种新的基于加扰技术的方法.仿真结果表明,该方法能够较好的抑制OFDM系统的峰均比,且具有复杂度低,易实现等优点.     

基于祖冲之类方法的多体动力学方程保能量保约束积分

针对一类多体动力学问题导出的微分 代数方程,提出一种保能量、保约束的算法.该算法基于祖冲之类方法和欧拉中点保辛差分,利用祖冲之类方法保证在时间格点上精确满足约束方程,避免约束违约问题;并进一步证明该算法在时间格点上可以精确保能量.数值算例进一步验证该算法的可靠性.     

基于祖冲之类方法的多体动力学方程保能量保约束积分

针对一类多体动力学问题导出的微分-代数方程,提出一种保能量、保约束的算法.该算法基于祖冲之类方法和欧拉中点保辛差分,利用祖冲之类方法保证在时间格点上精确满足约束方程,避免约束违约问题;并进一步证明该算法在时间格点上可以精确保能量.数值算例进一步验证该算法的可靠性.     

A collaborative benchmarking framework for multibody system dynamics

Despite the importance given to the computational efficiency of multibody system (MBS) simulation tools, there is a lack of standard benchmarks to measure the performance of these kinds of numerical simulations. This works proposes a collaborative benchmarking framework to measure and compare the performance of different MBS simulation methods. The framework is made up of two main components: (a) an on-line repository of test problems with reference solutions and standardized procedures to measure computational efficiency and (b) a prototype implementation of a collaborative web-based application to collect, organize and share information about performance results in an intuitive and graphical form. The proposed benchmarking framework has been tested to evaluate the performance of a commercial MBS simulation software, and it proved to be an effective tool to collect and analyze information about the numerous factors which affect the computational efficiency of dynamic simulations of multibody systems.     

A method for the numerical integration of mechanical systems with unilateral constraints: study of impact in multibody systems

Study of mechanical systems with unilateral constraints is associated with forming two systems of equations, a system of differential equations and a system of algebraical equations. Differential equations are used to describe the motion until the moment of impact, i.e. until activation of unilateral constraints. Algebraic equations are used to describe the impact. During numerical integration, transition from one system to another occurs at the points of impact. Even in simple problems, forming algebraic equations represents a complex task.This paper presents a method, the so-called Reduction Method, which provides for the analysis of these systems without forming the algebraic equations. They are substituted by a new system which is easily derived from equations of motion. Compared to methods based on the classical impact theory, using Reduction Method,velocities after the impact are easily computed regardless of the degrees of freedom.     

非完整约束多体系统时间离散变分积分法

基于连续Galerkin方法,给出非完整约束下多体系统时间离散的变分数值积分方法.首先对非完整多体系统Hamilton正则方程的弱形式进行时间离散,得到变分积分公式,然后讨论该积分方法对能量及约束的保持,最后以蛇板为例对该方法进行数值验证和比较.