第3讲_第三讲 运动副精确建模

f d n
－接触处滑动摩擦力； μ d －动摩擦系数； v ( q , q , t ) －接触点处的相对滑移速度。
Ff
当间隙铰构件碰撞接触后在接触区的切平面内运动状态为接触静止模式时， Fc 切向接触力 小于最大静摩擦力。 存在静摩擦力： Δ Fs μ s F n Fc 0 v (q , q , t ) 0
经典碰撞模式的动力学模型
在间隙铰的销轴与孔体接触前，系统动力学方程可表达为
自由运动
(1)
接触
T M q Kq Φ q λ Q Φ q , t 0
碰撞 过渡
Φ 式中，q为广义坐标列阵；M、K、 q 和分别为多体系统的广义 质量阵、刚度阵、约束方程的雅可比阵和广义速度二次项及广 义力阵； λ 为Lagrange乘子列阵。 约束方程:
S (q ( t ), t ) ε
则认为时刻 t 为拓扑结构切换点。
接触变形模式的动力学模型
－－运动模式识别方程
考察球铰或平面旋转铰（假定轴线不偏斜），设Pi和Pj分别为轴销i 和孔体j的轴心点，其相对位置矢量 d ij 为：
d
r
i
ij
r A d
i i
i
r
j
A d
j
j
和 r j 分别为体i和体j的体坐标系在系统惯性坐标系下的位置矢量； j i i j d 和 d 分别为点 P 和 P 在体坐标系下的位置矢量； i A 和A j 分别为轴销i和孔体j的局部坐标系与系统惯性坐标系的转换矩阵 矢量 d ij 的方向确定了法向接触力的作用方向 间隙铰k的运动模式识别方程：
Fn K
n
－－接触点处法向接触力
－－等效接触刚度
－－与有关的阻尼因子
e 1
e －－幂指数
δ －－接触点处法向穿透深度
δ －－接触点处法向相对速度
可通过对 K n 、 C (δ ) 和 e 的取值 得到不同类型的间隙模型
接触变形模式的动力学模型
－－碰撞接触力的处理
两体切向运动分为两种情况：相对滑动和粘滞。 当间隙铰构件碰撞接触后在接触区的切平面内运动状态为接触滑移模 式时，摩擦模型可选择库仑摩擦定律，即滑动摩擦力与正压力成正比， 其方向与相对滑移速度方向相反， 滑动摩擦力: F μ F sgn( v ) v (q , q , t ) 0
设 S ( q ( t Δ t ), t Δ t ) 和 S ( q ( t ), t ) 分别为在 t Δ t 和 t 时刻的间隙铰销轴 与孔体欲接触对应点的相对位置列阵，若满足
S ( q ( t ), t ) S ( q ( t Δ t ), t Δ t ) 0
T
则在该时间间隔内 t Δ t , t 至少存在一个切换点。 数值解算时，给定一个小的正数 ε ，当满足
的综合作用。对于接触面较小的碰撞，此种模型比较适用。
分离状态
接触状态
两状态模型
轴孔 分离状态 接触状态
销轴
三、含间隙机械系统动力学模型
经典碰撞模式 接触变形模式
分离状态
自由运动
接触
接触状态
轴孔 分离状态 接触状态
销轴
碰撞
过渡
1. 经典碰撞模式的动力学模型：着重考察系统动力学的整体行为，不 能预测碰撞过程中铰关节局部接触力的变化历程。 2. 接触变形模式的多体系统动力学模型：考虑了碰撞体间的法向和切 向接触力，以及阻力矩，该模型能够求解系统运动过程中间隙铰轴 销与孔体的碰撞力。
t
i
P
λ
P－－广义力的广义冲量列阵
H S / q
Lim λ dt Δ t 0 t
t
i
i
F(t)－－间隙铰碰撞力
S (q , t )
P Lim F ( t )dt Δ t 0 t
i
－－碰撞点间的距离列阵
经典碰撞模式的动力学模型
对上述方程数值求解，可求出广义速度增量和碰 撞冲量，得到碰撞后的初始状态。
经典碰撞模型
假定：

组成间隙铰的销轴与孔体是刚性的 碰撞瞬时完成，在碰撞前后碰撞体有速度“阶跃”和能量损失 采用动量定理和恢复系数计算碰撞分离后的动量分配。
这种间隙铰模型不考虑碰撞过程细节，是一种刚性碰撞模型，不能计 算实际的碰撞载荷大小，只能采用反作用力冲量衡量冲击造成的严重程度。
含间隙转动副 运动的三个状 态
sk d
ij
d
ij
δk 0
2
δk
－铰k处的间隙大小
四、含间隙机构运动方程的数值解法
运动方程特点： 1. 强非线性、病态、刚性方程组； 2. 拓扑结构时变性，易于发生突变，使得计算效率低、不稳 定，严重影响数值算法的可靠性； 3. 构件柔性导致系统产生振动冲击，引起宽频段响应，为模 态的选择带来困难； 4. 方程对误差和扰动非常敏感，易于产生误差和误差累积， 使结果失去可信度。 在含间隙的机构动力学中，离散方法及数值计算的精 度、效率以及稳定性等问题还需系统研究。
自由运动
接触
碰撞
过渡
接触变形模型
接触变形模型考虑精细的碰撞过程，采用无惯性的弹簧和阻尼器描述碰 撞体的变形和能量损失，目前一般采用Dubowsky等提出的“碰撞铰”模型。
此类模型需要构造接触体力－变形位移关系，由Hertz提出的静态接触 问题经典描述被多数研究者所采用，能量耗散一般考虑变形位移和变形速度
M
M
(i) v
(i) v
μ v Fn
(i)
－接触点i处的等效滚动摩擦力矩
M
(i) ω
旋转摩擦力矩：
M
(i) ω
μ ω ( Fn )
(i )Βιβλιοθήκη nωnω－接触点i处的等效转动摩擦力矩
－幂指数 －转动摩擦系数
接触变形模式的动力学模型
－－运动模式识别方程
对间隙铰轴销与孔体的“分离－接触”过程的切换点不能提前 预知，可由间隙铰两接触体对应点的相对位置监测系统运行状态。
一、含间隙机构动力学研究背景
桁架式展开机构、太阳能帆板、对接机构、锁定机构等
一、含间隙机构动力学研究背景
由于连接铰间隙导致的工程问题举例
1.
2.
3.
1997年5月，我国发射的东三通信卫星由于铰关节间隙的 存在，出现颤振 1999年2月，俄国宇航员在太空中未展开一面巨大的太空 镜子，使“人造月亮计划”失败 1990年，美国Hubble太空望远镜入轨后，因太阳能帆板连 接处存在间隙再加之昼夜温差的影响而导致帆板的“卡死 －滑动”周期性运动 ，使其不能精确定位
空间机构动力学影响因素 1.变拓扑结构（自由度变化） 2.间隙单边约束 3.空间环境(无根、失重、光压、真空、温差等） 4.碰撞、阻尼、摩擦等干扰
一、含间隙机构动力学研究背景
1.
2.
3.
连接间隙导致两种状态：带有摩擦的滑动和空转下 的内撞击。 运动副间隙的接触模型及其碰撞特性是最基本的问 题之一。 含铰接间隙的机构动力学，从力学本质上看是一种 非定常的、含碰撞的、变结构的、非线性的动力学 过程，连接铰间隙是结构非线性和无源阻尼的主要 来源。
接触变形模式的动力学模型
由于间隙铰元素间一般为低速内碰撞，可将间隙铰销轴与孔体的接触 －变形－恢复－脱离接触的接触碰撞过程归结为“自由运动－接触变形” 两种状态。当间隙铰两元素碰撞接触时，在碰撞体间引入等效接触力。
自由运行阶段：
T M q Kq Φ q λ Q Φ q , t 0
经典碰撞模式的动力学模型的缺陷
1. 2.
3.
4.
经典碰撞理论是在球状刚体碰撞实验中建立的，其理论还有待进一步 完善 从碰撞前后的速度、冲量、能量损失等出发，对恢复系数有不同的理 解和定义形式； 对碰撞过程中含摩擦的多体动力学建模将更为复杂，Kane在将这一理 论推广到处理一般性物体的碰撞问题时发现某些系统碰撞后的动能大 于碰撞前的动能，导致系统能量的增加，即所谓Kane动力学之谜； 经典碰撞理论是在一系列简化假定基础上建立的，将其原封不动地引 入含间隙多体系统动力学模型是有缺陷的，特别是对于非线性较强的 柔性系统，经典假设“碰撞过程瞬间完成，碰撞物体均作为刚体”不 能满足。
C2
C2 C1 B1
C1 B1 B2
B2
A2 A1 D1 D2
A2
A1
D1
D2
带铰间隙四杆机构
机构等效模型
有限元模型
间隙有限元模型需建立碰撞体接触后的准确动边界条件， 采用弹塑性有限元方法全面求解两物体的碰撞－接触变形过程。 这种建模方法可以得到碰撞体间作用力随时间和随空间的分布 规律。 但这种模型复杂，计算量较大，加之其处理刚体大位移的 局限性，目前还难以用于对复杂多体机械系统动力学数值模拟。
Δ F s －在接触处的静摩擦力与滑移力之差； μ s
－静摩擦系数。
接触变形模式的动力学模型
－－碰撞接触力的处理
考虑两接触体的切向变形，切向接触力：
F c K t ε (t )
K
t
－等效切向接触刚度
ε (t )
－接触点的切向局部弹性变形， 当有相对滑动时，为常数
μ v －滚动摩擦系数
滚动摩擦力矩：
分离状态
接触状态
接触变形阶段：
M q Kq Φ q λ Q F g Φ q , t 0
T
轴孔 分离状态 接触状态
销轴
F g －－接触力
F 相对于广义坐标 q 的广义力列阵
接触变形模式的动力学模型
－－碰撞接触力的处理
1.
2.
3.
4.
两构件接触变形，产生一个接触区，存在法向接触力、切 向接触力和阻力矩。 法向接触力可以采用接触力学理论确定，由于其模型复杂， 常简化为弹簧阻尼力表征； 切向接触力包括摩擦力（切向变形产生的接触力）和切向 阻尼力。 对于两接触体之间有滚动、旋转运动，还存在阻力矩，常 常表现为摩擦力矩。