多体运动方程

多体运动描述了多个互连刚体的动态行为。 刚体使用不同类型的接头耦合在一起,可以相对于彼此进行平移和旋转。 体的动态行为通过接头处的体运动方程和约束方程来描述。

两个或更多体通过一个或多个接头连接。 接头通过特定的运动学约束来定义。 因此,相对于彼此的体位置会受到运动约束条件的限制。 体连接一节中介绍了体之间的连接类型。

描述多体系统动态行为的常规方程如下:
1. EQUATION_DISPLAY
Mq¨=f
(4932)
其中,q 为广义坐标矢量,f 为广义力,M 为刚体的惯性矩阵的方块对角矩阵:
2. EQUATION_DISPLAY
M=(M100Mn)
(4933)

将运动方程与限制体的运动自由度的约束条件结合起来:

3. EQUATION_TITLE
ϕ(q,t)=0
(4934)

约束表示为体加速度的线性条件。 求解 Eqn. (4934) 的二阶导数得到:

4. EQUATION_TITLE
Jq¨=Q
(4935)

其中,Jϕ 的雅可比矩阵,q¨ 为体的加速度,Q 为不同质性。

为了强制约束的加速度条件,需要将约束力添加到系统中。 通过引入所有约束的拉格朗日乘数 λ,无效的约束力由以下方程给出:

5. EQUATION_TITLE
fc=JTλ
(4936)

需要寻找矢量 λ,以便在组合使用约束力 fc 与任何外力 fext(例如重力)时,都能够生成满足约束(由 Eqn. (4935) 给出)的运动。

根据 Eqn. (4936),可以将 Eqn. (4932) 写为:

6. EQUATION_TITLE
Mq¨=JTλ+fext
(4937)
求解 q¨
7. EQUATION_TITLE
q¨=M1JTλ+M1fext
(4938)

生成以下线性方程组:

8. EQUATION_TITLE
Aλ=b
(4939)
其中:
9. EQUATION_TITLE
A=JM1JT
(4940)
10. EQUATION_TITLE
b=JM1fext+Q
(4941)

λ 通过 Eqn. (4939) 获得。 在已知 λ 的情况下,对 Eqn. (4938) 进行两次积分,获得广义坐标矢量。

约束稳定

要获得坐标矢量 q,需要执行数值积分。 在 Simcenter STAR-CCM+ 中,将按 Eqn. (4935) 给出的加速度形式来施加约束。 对于此类方法,由于存在舍入误差和数值积分误差,因此当然会发生某些约束漂移。 约束稳定的目标是确保此类漂移不会随时间累积并保持在界限之内。

根据 Baumgarte [949],通过修改 Eqn. (4935) 应用约束稳定:

11. EQUATION_TITLE
Jq¨=Q
(4942)

且:

12. EQUATION_TITLE
Q=Q2α(Jq˙+Φt)β2Φ
(4943)

其中,αβ 为稳定法的参数。