第2章多体系统动力学基本理论.doc

1、第 2 章 多体系统动力学基本理论1第 2 章 多体系统动力学基本理论本章主要介绍多体系统动力学的基本理论，包括多刚体系统动力学建模、多柔体系统动力学建模、多体系统动力学方程求解及多体系统动力学中的刚性（Stiff）问题。通过本章的学习可以对多体系统动力学的基本理论有较深入的了解，为具体软件的学习打下良好的理论基础。2.1 多体系统动力学研究状况多体系统动力学的核心问题是建模和求解问题，其系统研究开始于 20 世纪 60 年代。从 60 年代到 80 年代，侧重于多刚体系统的研究，主要是研究多刚体系统的自动建模和数值求解；到了 80 年代中期，多刚体系统动力学的研究已经取得一系列成果，尤其是建

2、模理论趋于成熟，但更稳定、更有效的数值求解方法仍然是研究的热点；80 年代之后，多体系统动力学的研究更偏重于多柔体系统动力学，这个领域也正式被称为计算多体系统动力学，它至今仍然是力学研究中最有活力的分支之一，但已经远远地超过一般力学的涵义。本节将叙述多体系统动力学发展的历史和目前国内外研究的现状。2.1.1 多体系统动力学研究的发展机械系统动力学分析与仿真是随着计算机技术的发展而不断成熟的，多体系统动力学是其理论基础。计算机技术自其诞生以来，渗透到了科学计算和工程应用的几乎每一个领域。数值分析技术与传统力学的结合曾在结构力学领域取得了辉煌的成就，出现了以ANSYS、NASTRAN 等为代表的应

3、用极为广泛的结构有限元分析软件。计算机技术在机构的静力学分析、运动学分析、动力学分析以及控制系统分析上的应用，则在二十世纪八十年代形成了计算多体系统动力学，并产生了以 ADAMS 和 DADS 为代表的动力学分析软件。两者共同构成计算机辅助工程（CAE）技术的重要内容。多体系统是指由多个物体通过运动副连接的复杂机械系统。多体系统动力学的根本目的是应用计算机技术进行复杂机械系统的动力学分析与仿真。它是在经典力学基础上产生的新学科分支，在经典刚体系统动力学上的基础上，经历了多刚体系统动力学和计算多体系统动力学两个发展阶段，目前已趋于成熟。多刚体系统动力学是基于经典力学理论的，多体系统中最简单的情况

4、自由质点和一般简单的情况少数多个刚体，是经典力学的研究内容。多刚体系统动力学就是为多个刚体组成的复杂系统的运动学和动力学分析建立适宜于计算机程序求解的数学模型，并寻求高效、稳定的数值求解方法。由经典力学逐步发展形成了多刚体系统动力学，在发展过程中形成了各具特色的多个流派。早在 1687 年，牛顿就建立起牛顿方程解决了质点的运动学和动力学问题；刚体的概机械系统动力学分析及 ADAMS 应用2念最早由欧拉于 1775 年提出，他采用反作用力的概念隔离刚体以描述铰链等约束，并建立了经典力学中的牛顿-欧拉方程。 1743 年，达朗贝尔研究了约束刚体系统，区分了作用力和反作用力，达朗贝尔将约束反力称为“

5、丢失力” ，并形成了虚功原理的初步概念。1788 年，拉格朗日发表了分析力学 ，系统地研究了约束机械系统，他系统地考虑了约束，并提出了广义坐标的概念，利用变分原理考虑系统的动能和势能，得出第二类拉格朗日方程最少数量坐标的二阶常微分方程（ODE） ；并利用约束方程与牛顿定律得出带拉格朗日乘子的第一类拉格朗日方程最大数量坐标的微分代数方程（DAE） 。虚功形式的动力学普遍方程尚不能解决具有非完整约束的机械系统问题，1908 年若丹给出了若丹原理虚功率形式的动力学普遍方程，利用若丹原理可以方便地讨论碰撞问题和非完整系统的动力学问题。对于由多个刚体组成的复杂系统，理论上可以采用经典力学的方法，即以牛顿

6、-欧拉方法为代表的矢量力学方法和以拉格朗日方程为代表的分析力学方法。这种方法对于单刚体或者少数几个刚体组成的系统是可行的，但随着刚体数目的增加，方程复杂度成倍增长，寻求其解析解往往是不可能的。后来由于计算机数值计算方法的出现，使得面向具体问题的程序数值方法成为求解复杂问题的一条可行道路，即针对具体的多刚体问题列出其数学方程，再编制数值计算程序进行求解。对于每一个具体的问题都要编制相应的程序进行求解，虽然可以得到合理的结果，但是这个过程长期的重复是让人不可忍受的，于是寻求一种适合计算机操作的程式化的建模和求解方法变得迫切需要了。在这个时候，也就是 20世纪 60 年代初期，在航天领域和机械领域，

7、分别展开了对于多刚体系统动力学的研究，并且形成了不同派别的研究方法。最具代表性的几种方法是罗伯森-维滕堡（Roberson-Wittenburg）方法、凯恩（Kane）方法、旋量方法和变分方法。罗伯森与维滕堡于 1966 年提出一种分析多刚体系统的普遍性方法，简称为 R/W 方法，这种方法的主要特点是利用图论的概念及数学工具描述多刚体系统的结构，以邻接刚体之间的相对位移作为广义坐标，导出适合于任意多刚体系统的普遍形式动力学方程，并利用增广体概念对方程的系数矩阵作出物理解释。R/W 方法以十分优美的风格处理了树结构多刚体系统，对于非树系统，通过铰切割或刚体分割方法将非树系统转变成树系统进行处理。

8、凯恩方法是在 1965 年左右形成的分析复杂系统的一种方法，其利用广义速率代替广义坐标描述系统的运动，直接利用达朗伯原理建立动力学方程，并将矢量形式的力与达朗伯惯性力直接向特定的基矢量方向投影以消除理想约束力，兼有矢量力学和分析力学的特点，既适用完整系统，也适用于非完整系统。旋量方法是一种特殊的矢量力学方法（或牛顿-欧拉方法，简称为 N/E 方法） ，其特点是将矢量与矢量矩合为一体，采用旋量的概念，利用对偶数作为数学工具，使 N/E 方程具有极其简明的表达形式，在开链和闭链空间机构的运动学和动力学分析得到广泛运用。变分方法是不同于矢量力学或分析力学的另一类分析方法，高斯最小拘束原理是变分方法的

9、基本原理，保保夫和里洛夫从这一原理出发发展了两种不同风格的计算方法。该方法有利于结合控制系统的优化进行综合分析，而且由于其不受铰的约束数目的影响，适用于带多个闭环的复杂系统。第 2 章 多体系统动力学基本理论3这几种方法构成了早期多刚体系统动力学的主要内容，借助计算机数值分析技术，可以解决由多个物体组成的复杂机械系统动力学分析问题。但是多体系统动力学在建模与求解方面的自动化程度，相对于结构有限元分析的成熟来说相差甚远。正是为了解决多体系统动力学建模与求解的自动化问题，美国 Chace 和 Haug 于 80 年代提出了适宜于计算机自动建模与求解的多刚体系统笛卡尔建模方法，这种方法不同于以罗伯森

10、-维滕堡方法为代表的拉格朗日方法，它是为以系统中每个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，再根据铰约束和动力学原理建立系统的数学模型进行求解。20 世纪 80 年代，Haug 等人确立了“计算多体系统动力学 ”这门新的学科，多体系统动力学的研究重点由多刚体系统走向侧重多柔体系统，柔性多体系统动力学成为计算多体系统动力学的重要内容。柔性多体系统动力学在 20 世纪 70 年代逐渐引起人们的注意，一些系统如高速车辆、机器人、航天器、高速机构、精密机械等其中柔性体的变形对系统的动力学行为产生很大影响

11、。二十多年来柔性多体系统动力学一直是研究热点，这期间产生了许多新的概念和方法，有浮动标架法、运动-弹性动力学方法、有限段方法以及最新提出的绝对节点坐标法等，其中浮动标架法最早是在航天领域研究中提出来的。计算多体系统动力学是指用计算机数值手段来研究复杂机械系统的静力学分析、运动学分析、动力学分析以及控制系统分析的理论和方法。相比于多刚体系统，对于柔性体和多体与控制混合问题的考虑是其重要特征。其具体任务为：1建立复杂机械系统运动学和动力学程式化的数学模型，开发实现这个数学模型的软件系统，用户只需输入描述系统的最基本数据，借助计算机就能自动进行程式化处理。2开发和实现有效的处理数学模型的计算方法与数

12、值积分方法，自动得到运动学规律和动力学响应。3实现有效的数据后处理，采用动画显示、图表或其他方式提供数据处理结果。计算多体系统动力学的产生极大地改变了传统机构动力学分析的面貌，使工程师从传统的手工计算中解放了出来，只需根据实际情况建立合适的模型，就可由计算机自动求解，并可提供丰富的结果分析和利用手段；对于原来不可能求解或求解极为困难的大型复杂问题，现可利用计算机的强大计算功能顺利求解；而且现在的动力学分析软件提供了与其它工程辅助设计或分析软件的强大接口功能，它与其它工程辅助设计和分析软件一起提供了完整的计算机辅助工程（CAE）技术。2.1.2 多体系统动力学研究活动自 20 世纪 60 年代以

13、来，国内外在多体系统动力学方面多次召开了深具意义的会议。国际范围内的会议主要有：1977 年国际理论和应用力学学会（International Union of Theoretical and Applied Mechanics - IUTAM）发起在德国慕尼黑由 Magnus 主持召开第一次多刚体系统动力学讨论会；。1983 年北大西洋公约组织与美国国家科学基金委等（NATO-NSF-ARD）联合组织机械系统动力学分析及 ADAMS 应用4在美国爱阿华由 Haug 主持召开“机械系统动力学计算机辅助分析与优化高级研讨会” ；。1985 年第八届国际车辆动力学协会（International

14、Association of Vehicle System Dynamics - IAVSD）会议， Kortum 和 Schiehlen 发表了用于车辆动力学仿真的多体软件。1985 年 IUTAM 与国际机器及机构理论联合会（IFTOMM）联合在意大利 Udine 由Bianchi 和 Schiehlen 主持举行了第二届国际多体系统动力学讨论会，这次会议总结了该领域的进展，标志多刚体系统动力学已趋于成熟。1989 年由德国斯图加特大学主持对当时比较先进的大型软件进行测试，编辑出版了“多体系统手册” ；以后几乎每年都有国际的多体系统动力学的会议，并出现了多体系统动力学的专门的刊物。在国内召

15、开的关于多体系统动力学方面的重要会议主要有：1986 年由中国力学学会一般力学专业委员会在北京主持召开“多刚体系统动力学”研讨会。1988 年在长春召开“柔性多体系统动力学研讨会” 。1992 年在上海召开“全国多体系统动力学理论、计算方法与应用学术会议” 。1996 年由中国力学学会一般力学专业委员会与中国空间学会空间机械委员会联合在山东长岛召开“全国多体系统动力学与控制学术会议” 。除了国内外重要的会议之外，国内外还出版了多种关于多体系统动力学的优秀教材或专著。国外的优秀教材和专著主要有：1977 年，Witenburg 出版了 “Dynamics of Systems of Rigid

16、Bodies”，这是第一本关于多体系统的教材，其奠定了多刚体系统动力学拉朗日方法的基础，并将图论方法引入到多体系统动力学中。1986 年，Schiehlen 出版了 “Technische Dynamik”（德语） ，其将多体系统作为与有限元系统及连续系统相当的系统来统一考虑。1988 年，Nikravesh 出版了“Computer-Aided Analysis of Mechanical Systems”。1988 年，Roberson 和 Schwertassek 出版了“Dynamics of Multibody System” 。1989 年，Haug 出版了“Computer-Ai

17、ded Kinematics and Dynamics of Mechanical Systems”，提出了多刚体系统建模的笛卡尔方法。1989 年，Shabana 出版了“Dynamics of Multibody Systems”，这是一本关于多柔体动力学的专著。1990 年，Huston 出版了“ Multibody Dynamics”，其讨论了多体系统动力学、数值方法以及柔性多体系统动力学。1992 年，Bremer 和 Pfeiffer 出版了“Elastische Mehrkorpersystem”（德语） ，其主要讨论了柔性体，并给出了诸多工程实例。1992 年，Amirouch

18、e 出版了“Computational Methods for Multibody Dynamics”，重点讨论了矩阵方法。1994 年，Garcia 和 Bayo 出版了 “Kinematic and Dynamic Simulation of Multibody 第 2 章 多体系统动力学基本理论5Systems”，提出了完全笛卡尔坐标方法，并给出了一种求解效率高的计算方法用于实时仿真的需要。1994 年，Shabana 出版了“Computational Dynamics” ，主要是有关多刚体系统。1996 年，Stejskal 和 Valasek 出版了“Kinematics and

19、Dynamics of Machinery”，从空间机构的 CAD 设计入手，讨论了自由刚体问题，给出了高副和低副运动学约束的描述，并讨论了动力学分析以及数值计算方面的问题。国内有影响的教材和专著主要有： 刘延柱，洪嘉振，杨海兴. 多刚体系统动力学. 北京：高等教育出版社，1989 黄文虎 邵成勋 等. 多柔体系统动力学. 北京：科学出版社，1996 陆佑方. 柔性多体系统动力学. 北京：高等教育出版社，1996 洪嘉振. 计算多体系统动力学. 北京：高等教育出版社，19992.1.3 多体系统动力学研究现状计算多体系统动力学中所研究的多体系统，根据系统中物体的力学特性可分为多刚体系统、柔性多

20、体系统和刚柔混合多体系统。多刚体系统是指可以忽略系统中物体的弹性变形而将其当作刚体来处理的系统，该类系统常处于低速运动状态；柔性多体系统是指系统在运动过程中会出现物体的大范围运动与物体的弹性变形的耦合，从而必须把物体当作柔性体处理的系统，大型、轻质而高速运动的机械系统常属此类；如果柔性多体系统中有部分物体可以当作刚体来处理，那么该系统就是刚柔混合多体系统，这是多体系统中最一般的模型。12.1.3.1 多体系统建模理论对于多刚体系统，从 20 世纪 60 年代到 80 年代，在航天和机械两个领域形成了两类不同的数学建模方法，分别称为拉格朗日方法和笛卡尔方法；20 世纪 90 年代，在笛卡尔方法的

21、基础上又形成了完全笛卡尔方法。这几种建模方法的主要区别在于对刚体位形描述的不同。航天领域形成的拉格朗日方法，是一种相对坐标方法，以 Roberson-Wittenburg 方法为代表，是以系统每个铰的一对邻接刚体为单元，以一个刚体为参考物，另一个刚体相对该刚体的位置由铰的广义坐标（又称拉格朗日坐标）来描述，广义坐标通常为邻接刚体之间的相对转角或位移。这样开环系统的位置完全可由所有铰的拉格朗日坐标阵 所确定。q其动力学方程的形式为拉格朗日坐标阵的二阶微分方程组，即 ),(),(tqBtA(2.1-1)这种形式首先在解决拓扑为树的航天器问题时推出。其优点是方程个数最少，树系统的坐标数等于系统自由度

22、，而且动力学方程易转化为常微分方程组（ODEs - Ordinary Differential Equations） 。但方程呈严重非线性，为使方程具有程式化与通用性，在矩阵与 中常常包含描述系统拓扑的信息，其形式相当复杂，而且在选择广义坐标时需人AB为干预，不利于计算机自动建模。不过目前对于多体系统动力学的研究比较深入，现在有几种应用软件采用拉格朗日的方法也取得了较好的效果。机械系统动力学分析及 ADAMS 应用6对于非树系统，拉格朗日方法要采用切割铰的方法以消除闭环，这引入了额外的约束，使得产生的动力学方程为微分代数方程，不能直接采用常微分方程算法去求解，需要专门的求解技术。机械领域形成的

23、笛卡尔方法是一种绝对坐标方法，即 Chace 和 Haug 提出的方法，以系统中每一个物体为单元，建立固结在刚体上的坐标系，刚体的位置相对于一个公共参考基进行定义，其位置坐标（也可称为广义坐标）统一为刚体坐标系基点的笛卡尔坐标与坐标系的方位坐标，方位坐标可以选用欧拉角或欧拉参数。单个物体位置坐标在二维系统中为 3 个，三维系统中为 6 个（如果采用欧拉参数为 7 个） 。对于由 N 个刚体组成的系统，位置坐标阵 中的坐标个数为 3N（二维）或 6N（或 7N） （三维） ，由于铰约束的存在，q这些位置坐标不独立。系统动力学模型的一般形式可表示为. 0),(tBATq(2.1-2)式中 为位置坐

24、标阵 的约束方程， 为约束方程的雅可比矩阵， 为拉格朗日乘子。q这类数学模型就是微分-代数方程组（ DAEs - Differential Algebraic Equations） ，也称为欧拉-拉格朗日方程组（Euler-Lagrange Equations） ，其方程个数较多，但系数矩阵呈稀疏状，适宜于计算机自动建立统一的模型进行处理。笛卡尔方法对于多刚体系统的处理不区分开环与闭环（即树系统与非树系统） ，统一处理。目前国际上最著名的两个动力学分析商业软件 ADAMS 和 DADS 都是采用这种建模方法。完全笛卡尔坐标方法，由 Garcia 和 Bayo 于 1994 年提出，是另一种形式

25、的绝对坐标方法。这种方法的特点是避免使用一般笛卡尔方法中的欧拉角或欧拉参数，而是利用与刚体固结的若干参考点和参考矢量的笛卡尔坐标描述刚体的空间位置与姿态。参考点选择在铰的中心，参考矢量沿铰的转轴或滑移轴，通常可由多个刚体共享而使未知变量减少。完全笛卡尔坐标所形成的动力学方程与一般笛卡尔方法本质相同，只是其雅可比矩阵为坐标线性函数，便于计算。至于柔性多体系统，从计算多体系统动力学角度看，柔性多体系统动力学的数学模型首先应该和多刚体系统与结构动力学有定的兼容性。当系统中的柔性体变形可以不计时，即退化为多刚体系统。当部件间的大范围运动不存在时，即退化为结构动力学问题。柔性多体系统不存在连体基，通常选

26、定一浮动坐标系描述物体的大范围运动，物体的弹性变形将相对该坐标系定义。弹性体相对于浮动坐标系的离散将采用有限单元法与现代模态综合分析方法。在用集中质量有限单元法或一致质量有限单元法处理弹性体时，用结点坐标来描述弹性变形。在用正则模态或动态子结构等模态分析方法处理弹性体时用模态坐标描述弹性变形。这就是莱肯斯首先提出的描述柔性多体系统的混合坐标方法。即用坐标阵 描述系统的位形，其中 q 为浮动坐标系的位形坐标，a 为变形坐标。考Taqp)(虑到多刚体系统的两种流派，在柔性多体系统动力学中也相应提出两种混合坐标，即浮动坐标系的拉格朗日坐标加弹性坐标与浮动坐标系的笛卡尔坐标加弹性坐标。根据动力学基本原

27、理推导的柔性多体系统动力学方程，形式同式（2.1-1）和（2.1-2 ） ，第 2 章 多体系统动力学基本理论7只是将 q 用 p 代替。即，柔性多体系统具有与多刚体系统类同的动力学数学模型。2.1.3.22 .多体系统动力学数值求解多刚体系统拉格朗日方法产生的形如式（2.1-1）的动力学数学模型，是形式复杂的二阶常微分方程组（ODEs） ，系数矩阵包含描述系统拓扑的信息。对于该类问题的求解，通常采用符号-数值相结合的方法或者全数值的方法。符号 -数值方法是先采用基于计算代数的符号计算方法，进行符号推导，得到多刚体系统拉格朗日模型系数矩阵简化的数学模型，再用数值方法求解 ODE 问题。鉴于计算

28、机技术的发展，目前全数值方法也较为流行，就是将多刚体系统拉格朗日数学模型当作一般 ODE 问题进行求解，这方面的技术已经较为成熟。多刚体系统笛卡尔方法产生的形如式（2.1-2）的动力学数学模型，是著名的微分- 代数方程组（DAEs） 。DAE 问题是计算多体系统动力学领域的热点问题。柔性多体系统的动力学数学模型，其形式与多刚体系统相同，可以借鉴多刚体系统数学模型的求解方法。只是混合坐标中描述浮动坐标系运动的刚体坐标 q 通常是慢变大幅值的变量，而描述相对于浮动坐标系弹性变形的坐标 a 却为快变微幅的变量，两类变量出现在严重非线性与时变的耦合动力学方程中，其数值计算呈病态，将出现多刚体系统中见不

29、到的数值计算困难。综上所述，多体系统动力学问题的求解集中于微分-代数方程组的求解，下面将简要地介绍一下 DAE 问题的求解方法，更具体的介绍，将在（2.5）节进行。1、 （1）微分-代数方程组的特性多刚体系统采用笛卡尔方法建模生成的微分-代数方程组为： 0),(),(),( tqQtqtMT(2.1-3). 0(2.1-4)其中， 、 、 分别是系统位置、速度、加速度向量， 是拉格朗日乘nRmR子， 是时间， 为机械系统惯性矩阵， 为约束雅可比矩阵，t nmq为外力向量， 为位置约束方程。nRQm将式(1.4)对时间求一阶和二阶导数，得到速度和加速度约束方程：. 0),(),(),( tqtt

30、qq(2.1-5). (2.1-6)其中， 称为速度右项， 称为加速度右项。),(t tqtq2)(给定方程组初始条件：. 0)(q机械系统动力学分析及 ADAMS 应用8(2.1-7)微分-代数方程组的特性和需要注意的问题有：微分- 代数方程问题不是常微分方程（ ODE）问题；由式(2.1-3)和(2 .1-4)组成的微分- 代数方程组是指标 3 问题，通过对约束方程求导化为由式(2 .1-3)-( 2.1-6)组成的微分- 代数方程组后，其指标降为 1；微分- 代数方程数值求解的关键在于避免积分过程中代数方程的违约现象；初值式(2.1-7)与位置约束式(2.1-4)及速度约束式(2 .1-

31、5)的相容性；微分- 代数方程组的刚性问题。（2）2、微分-代数方程组积分技术自二十世纪七十年代以来，国际上对微分-代数方程问题作了大量的研究，时至如今，新的算法仍不断涌现。根据对位置坐标阵和拉格朗日乘子处理技术的不同，可以将微分-代数方程组问题的处理方法分为增广法和缩并法。传统的增广法是把广义坐标加速度 和拉格朗日乘子 作为未知量同时求解，再对加q速度 进行积分求出广义坐标速度 及广义坐标位置 ，包括直接积分法和约束稳定法。q近十年来，在传统增广法的基础上又发展形成了超定微分- 代数方程组（ODAEs）方法等新的一类算法。直接积分法：将式(2.1-3) 和(2.1-6)联立在一起，同时求出

32、与 ，然后对 积分得qq和 。该方法未考虑式(2.1 -4)和(2.1-5)的坐标和速度违约问题，积分过程中误差积累严重，很易发散。在实际的数值计算过程中，并不直接采用直接积分法，但在直接积分法的基础上发展了一系列控制违约现象的数值方法。约束稳定法：将控制反馈理论引入微分-代数方程组的数值积分过程以控制违约现象。通过把式(2 .1-6)右边量替换为含位置约束和速度约束的参数式，保证位置约束和速度约束在式(2.1 -3)和(2.1-6) 联立求解时恒满足。该方法稳定性好，响应快，但如何选择参数式中速度项和位置项适当的系数是一个问题。超定微分-代数方程组（ODAEs）法：将系统速度作为变量引入微分

33、-代数方程组，从而将原来的二阶 DAE 化为超定的一阶 DAE，再为所得方程组引入未知参数，根据模型的相容性消除系统的超定性，如此可使数值计算的稳定性明显改变。或者将系统位置、速度、加速度向量和拉格朗日乘子向量联立作为系统广义坐标，再将由式(2.1-3) 、(2.1-4)、(2 .1-5)和(2.1-6)组成的微分-代数方程组及速度与位置、加速度与速度的微分关系式作为约束，化二阶 DAE 为超定的一阶 DAE，再根据系统相容性引入二个未知参数，消除超定性，这样所得的最终约化模型更为简单，但方程组要多 n 个。在 ODAE 方法的基础上产生了一系列新的更为有效的算法。解耦 ODAE 法：在 OD

34、AE 方法的基础上，发展形成了一类解耦思想，就是在 ODAEs基础上，对常用的隐式 ODE 方法采用预估式，再按加速度、速度和位置的顺序进行求解。后来进一步发展形成了无需对隐式 ODE 方法利用预估式的解耦思想，更一步地提高了效率。缩并法就是通过各种矩阵分解方法将描述系统的 n 个广义坐标用 p 个独立坐标表达，从而将微分-代数方程组从数值上化为与式 (2.1-1)类似的数学模型，如此易于用 ODE 方法第 2 章 多体系统动力学基本理论9进行求解。传统的缩并法包括 LU 分解法、QR 分解法、SVD 分解法以及零空间方法等，后来在传统缩并法的基础上产生了局部参数化缩并方法等新的算法。缩并法中

35、的这些具体方法，分别对应着约束雅可比矩阵的不同分解。LU 分解法：又称为广义坐标分块法。把广义位置坐标 用相关坐标 和独立坐标qu分块表示，再将约束雅可比矩阵 用 LU 分解法分块，得到广义坐标速度 、加速度vqq用独立坐标速度 、加速度 表达的式子。将这两个表达式代入式(2.1-3) ，就可得到形qv如式(2.1 -1)的关于独立坐标加速度 的二阶微分方程。该算法可靠、精确，并可控制误差，v但效率稍低。QR 分解法：通过对约束雅可比矩阵 正交分解的结果作微分流型分析，得到可选q作受约束系统独立速度的 ，并将微分-代数方程组化作形如式(2.1-1)的关于 的二阶微分z z方程，如此可保证在小时

36、间间隔内由 积分引起的广义坐标的变化不会导致大的约束违约。zSVD 分解法：把约束雅可比矩阵 作奇异值分解所得结果分别用于式(2.1-3) 和(2.1-6) ，q得到缩并后的系统动力学方程。在该方法推导过程中没有用到式(2.1-4) 和(2.1-5)，所以也存在位置和速度违约问题，可用约束稳定法改善其数值性态。可微零空间法：通过 Gram-Schmidt 正交化过程自动产生约束雅可比矩阵 的可微、q唯一的零空间基，来对系统方程降阶。具体做法是对由 和任意矩阵nmqR构造的矩阵 采用 Gram-Schmidt 正交化过程，将 化为正交非奇异nmRB)( nRPP矩阵 。再引入新的速度矢量 ，使满

37、足 ，将新速度矢量 和加速度矢量VzVzTz按正交化结果分块，得到新的独立速度矢量 和加速度矢量 。如此可将微分-代数方z II程组化为关于新的独立加速度矢量 的动力学方程。I局部参数化缩并方法：先将式(2.1-3 )-( )( 2.1-6)改写为等价的一阶形式，再用微分流形理论的切空间局部参数化方法将等价的欧拉-拉格朗日方程降为参数空间上的常微分方程。总的说来，微分-代数方程组数值求解的方法都可归为增广法或缩并法，除了上面所介绍的这些增广法和缩并法所运用的增广和缩并技术外，近几年来还出现了不少独具特色的处理算法，它们或者是在数值求解算法中独具匠心，或者针对某些具体情况作了专门研究。（3、 ）

38、相容性问题和刚性问题初值相容性问题：在微分-代数方程组的数值求解过程中，给定的位置和速度初始条件与微分- 代数方程组中的位置和速度约束的相容性是值得注意的一个问题。相容性是微分-代数方程组有解的必要条件。刚性问题：由于现代机械系统的复杂性，会由于系统的耦合而使所得到的微分-代数方程组呈现刚性特性。对于刚性问题的求解，目前最常用的方法是隐式方法，隐式方法不仅用于求解刚性问题，而且相比于显式方法具有更好的稳定性和计算精度。近几年来，无论是在 LU 分解法基础上发展起来的新缩并法，还是基于 ODAE 方法的增广法，或是基于多体系统正则方程的解法，应用的无不是隐式方法。关于刚性问题的进一步介绍见 2.

39、5 节。机械系统动力学分析及 ADAMS 应用102.2 多刚体系统动力学建模计算多体系统动力学分析，首先在于提供一个友好方便的界面以利于建立多体系统的力学模型，并在系统内部由多体系统力学模型得到动力学数学模型；再者需要有一个优良的求解器对数学模型进行求解，求解器要求效率高、稳定性好，并具有广泛的适应性；最后还需要对求解结果提供丰富的显示查询手段。这其中的关键技术就是自动建模技术和求解器设计，所谓自动建模就是由多体系统力学模型自动生成其动力学数学模型，求解器的设计则必须结合系统的建模，以特定的动力学算法对模型进行求解。本节首先给出多体系统动力学的一些基本概念，再介绍计算多体系统动力学建模与求解

40、的一般过程，然后作为重点，按豪格（Haug）的笛卡尔方法给出多体系统的运动学和动力学数学模型，这是考虑到目前最负盛名的两个多体系统软件都是采用笛卡尔方法，最后简要叙述自动建模技术。2.2.1 多体系统动力学基本概念物理模型：这里也称力学模型，由物体、铰、力元和外力等要素组成并具有一定拓扑构型的系统。拓扑构型：多体系统中各物体的联系方式称为系统的拓扑构型，简称拓扑。根据系统拓扑中是否存在回路，可将多体系统分为树系统与非树系统。系统中任意两个物体之间的通路唯一，不存在回路的，称为树系统；系统中存在回路的称为非树系统。物体：多体系统中的构件定义为物体。在计算多体系统动力学中，物体区分为刚性体（刚体）

41、和柔性体（柔体） 。刚体和柔体是对机构零件的模型化，刚体定义为质点间距离保持不变的质点系，柔体定义为考虑质点间距离变化的质点系。约束：对系统中某构件的运动或构件之间的相对运动所施加的限制称为约束。约束分为运动学约束和驱动约束，运动学约束一般是系统中运动副约束的代数形式，而驱动约束则是施加于构件上或构件之间的附加驱动运动条件。铰：也称为运动副，在多体系统中将物体间的运动学约束定义为铰。铰约束是运动学约束的一种物理形式。力元：在多体系统中物体间的相互作用定义为力元，也称为内力。力元是对系统中弹簧、阻尼器、致动器的抽象，理想的力元可抽象为统一形式的移动弹簧-阻尼器- 致动器（TSDA ） ，或扭转弹簧-阻尼器- 致动器（RSDA） 。外力（偶）：多体系统外的物体对系统中物体的作用定义为外力（偶） 。数学模型：分为静力学数学模型、运动学数学模型和动力学数学模型，是指在相应条件下对系统物理模型（力学模型）的数学描述。机构：装配在一起并允许作相对运动的若干个刚体的组合。运动学：研究组成机构的相互联接的构件系统的位置、速度和加速度，其与产生运动