2.3多柔体系统动力学建模2.3.1柔性体上点的位置向量、速度和加速度1．柔性体系统中的坐标系图2.5柔性体上节点的位置柔性体系统中的坐标系如图2.5所示，包括惯性坐标系（）和动坐标系（）。前者不随时间而变化，后者是建立在柔性体上，用于描述柔性体的运动。动坐标系可以相对惯性坐标系进行有限的移动和转动。动坐标系在惯性坐标系中的坐标（移动、转动）称为参考坐标。与刚体不同，柔性体是变形体，体内各点的相对位置时时刻刻都在变化，只靠动坐标系不能准确描述该柔性体在惯性坐标系中的位置，因此，引入弹性坐标来描述柔性体上各点相对动坐标系统的变形。这样柔性体上任一点的运动就是动坐标系的“刚性”运动与弹性变形的合成运动。由于柔体上各点之间有相对运动所以动坐标系的选择不是采用连体坐标系，而需要采用随着柔性体形变而变化的坐标系，即“浮动坐标系”。在研究多柔体系统时，合适的坐标系是非常重要的。在确定浮动坐标系时有两点准则：1、便于方程建立求解；2、柔性体刚体运动与变形运动的耦合尽量小。目前常见的浮动坐标系大致有如下5种，局部附着框架、中心惯性主轴框架、蒂斯拉德框架、巴克凯恩斯框架以及刚体模态框架。采用何种需因实际情况而定。2．柔性体上任意点的位置向量、速度和加速度在分析刚体平面运动的时候，把复杂的刚体平面运动分解为几种简单的运动在对柔性体的运动，尤其是在小变形的情况下，也可以采用类似的方法。如某柔性体从位置L1运动到位置L2，其间运动可以分解为：刚性移动－>刚性转动－>变形运动。对于柔性体上任意一点P，其位置向量为：（2-158）为P点在惯性坐标系中的向量；为浮动坐标系原点在惯性坐标系中的向量；根据式2.2-14，A为方向余弦矩阵；为柔性体未变形时P点在浮动坐标系中的向量;为相对变形向量，可以用不同的方法离散化，与讨论平面问题相同，对于点，该单元的变形采用模态坐标来描述，有：（2-159）式（2-159）中，为点满足里兹基向量要求的假设变形模态矩阵，为变形的广义坐标。柔性体上任一点的速度向量及加速度向量可以对式求对时间一阶导数和二阶导数得到：（2-160）（2-161）2．多柔体系统的能量（1）动能和质量矩阵考虑节点变形前后的位置、方向和模态，柔性体的广义坐标可以表示为：（2-176）速度表达式（2-160）在系统广义坐标式（2-176）的时间导数中表示为：（2-177）柔性体的动能为：（2-178）其中，和分别为节点的节点质量和节点惯性张量；，为点相对于全局坐标基的角速度在局部坐标基中的斜方阵表示。将式（2.3-20）和G()reBrB()bePsPPuP图2.6柔性体变形模型关系式代入式（2-178），得到动能的广义表达式：（2-179）上式中的质量矩阵为维的方阵，表示为：（2-180）其中下标分别表示平动、旋转和模态自由度。质量矩阵的六个独立分量分别表示为：（2-181）其中9个惯性时不变矩阵列表如下：式(2.3-24)中可以明显看出质量矩阵与模态坐标显性相关，而且由于引入转换矩阵和，质量矩阵也与系统的方向坐标显性相关。质量矩阵中的9个惯性时不变矩阵可通过计算有限元模型的个节点信息在预处理过程中一次性得到，从而简化运动微分方程的求解。节点信息包括：每个节点的质量，未变形时的位置矢量以及模态矩阵。（2）势能和刚度矩阵势能一般分为重力势能和弹性势能两部分，可用下列二次项表示：（2-182）在弹性势能中，是对应于模态坐标的结构部件的广义刚度矩阵，通常为常量。重力势能表示为：（2-183）其中表示重力加速度矢量，重力可对求导得；（2-184）（3）能量损失和阻尼矩阵阻尼力依赖于广义模态速度并可以从下列二次项中推导得出：（2-185）上式称为Rayleigh能量损耗函数。矩阵包含阻尼系数，它是常值对称阵。当引入正交模态振型时，阻尼矩阵可用对角线为模态阻尼率的对角阵来表示。对于每一个正交模态的阻尼率都可以取不同值，而且还能以该模态的临界阻尼的比值形式给出。3．多柔体动力学方程柔性体的运动方程从下列拉格朗日方程导出：（2-186）其中，为约束方程；为对应于约束方程的拉氏乘子；为如式（2.3-19）定义的广义坐标；为投影到上的广义力；为拉格朗日项，定义为，和分别表示动能和势能；表示能量损...
