1. 课程导入与多体系统摩擦间隙问题全景
多体系统动力学研究的对象，是由多个刚体或柔性体通过运动副连接而成的机械系统，例如曲柄滑块机构、机械臂、车辆悬架、卫星太阳翼展开机构、航空发动机叶片榫连结构等。在理想模型中，运动副被当作无间隙、无摩擦的几何约束，关节力通过约束方程直接求解。但真实机械中，铰链孔与销轴之间存在装配间隙，接触面之间存在摩擦，润滑状态会随载荷与速度变化，构件还会发生弹性变形。这些因素叠加后，系统的动力学行为会明显偏离理想模型：出现冲击、颤振、粘滑、分岔甚至混沌，关节磨损加速，定位精度下降，噪声与振动超标。
本课程的核心任务，是建立一套能够描述"间隙—接触—摩擦—动力学响应"耦合关系的分析框架，并把它落到可计算、可判断、可验证的工程流程上。本章作为导入，先给出全局图景：学什么、为什么学、整体框架如何组织、典型工程场景长什么样。
学习目标
- 理解多体系统、运动副间隙、摩擦非线性三者的耦合机制，能说清"间隙如何引入自由度、摩擦如何引入能量耗散与状态切换"。
- 掌握含间隙运动副的建模思路：从几何约束到接触力模型，从库仑摩擦到LuGre、Stribeck等连续摩擦模型。
- 能针对给定机构（如曲柄滑块、四杆机构）建立含间隙动力学方程，并选择合适的求解器与步长策略。
- 能识别数值结果中的典型非线性特征：接触力尖峰、粘滑平台、Poincaré截面上的分岔、频谱中的边带。
- 能对分析结果做工程判断：间隙是否可接受、摩擦参数是否需修正、是否需要引入柔性体或润滑模型。
核心概念
多体系统（Multibody System）：由若干物体通过运动副连接的系统，其位形可用一组广义坐标描述。理想运动副施加完整约束，系统自由度等于广义坐标数减去独立约束数。
运动副间隙（Joint Clearance）：销轴与轴套之间的径向游隙。间隙使运动副从"约束"退化为"接触—分离"的变拓扑问题，系统自由度增加，接触状态随相对运动切换。
摩擦非线性（Friction Nonlinearity）：摩擦力随相对速度、法向力、材料与润滑状态变化，具有强非线性与不连续性。库仑摩擦在零速附近不连续，Stribeck效应使摩擦系数随速度先降后升，LuGre模型用内部状态变量描述预滑移。
接触力模型：常用Lankarani-Nikravesh模型，法向力由弹性项与滞回阻尼项组成，形式为 F_n = K δ^n + D δ̇，其中δ为穿透深度，K为接触刚度，D为阻尼系数。该模型能反映接触过程中的能量损失。
粘滑（Stick-Slip）：当驱动力不足以克服静摩擦时，接触面粘滞；当相对速度超过阈值后转为滑动。粘滑切换是摩擦诱导振动与噪声的重要来源。
整体框架
本课程按"建模—求解—分析—验证"的主线展开，整体框架可概括为下图所示的信号链与数据流关系。
图中左侧是输入与对象，中间是建模与求解，右侧是分析输出，底部是典型场景。红色虚线表示反馈：分析结果若与试验或工程经验不符，需要回到建模层修正接触刚度、摩擦参数或间隙值，而不是只调求解器。
典型应用场景与工程背景
以曲柄滑块机构为例。设曲柄长 r = 50 mm，连杆长 l = 200 mm，曲柄转速 n = 600 r/min，即角速度 ω ≈ 62.8 rad/s。若曲柄与连杆之间的铰链存在径向间隙 c = 0.1 mm，销轴半径 R = 10 mm，接触刚度 K = 1×10⁸ N/m^1.5，阻尼系数 D = 200 N·s/m，摩擦系数 μ = 0.15。在理想模型中，滑块加速度是光滑的周期函数；引入间隙后，销轴在轴套内往复碰撞，接触力出现高频尖峰，峰值可达理想关节力的3～8倍。若同时考虑摩擦，在曲柄接近死点、相对速度过零时会出现粘滑，滑块速度曲线出现平台，加速度频谱在啮合频率附近出现边带。
这类问题的工程后果很直接：关节磨损速率与接触力峰值的高次方相关，间隙会随磨损进一步增大，形成正反馈；定位精度下降会影响加工质量；噪声与振动会传递到整机。因此，分析的目的不是"算出曲线"，而是回答三个问题：当前间隙与摩擦参数下，系统是否稳定？接触力峰值是否超过许用值？若超限，应优先减小间隙、改善润滑，还是调整驱动策略？
本课程的组织方式
后续内容将按以下逻辑推进：先建立含间隙运动副的几何与接触力模型，再引入摩擦模型并讨论粘滑判别，然后给出多体系统动力学方程的组装与求解策略，接着通过曲柄滑块、四杆机构、机械臂关节等案例展示完整分析流程，最后讨论参数辨识、模型验证与工程判断准则。每一章都会给出可复现的示例参数、关键公式和结果判读方法，避免停留在概念层面。
需要强调的是，摩擦间隙非线性分析没有"一次建模、永久适用"的通用模型。接触刚度、阻尼、摩擦系数都依赖材料、表面粗糙度、润滑状态与载荷水平，必须结合试验或工程经验标定。本课程的目标，是让你具备"选模型、定参数、算响应、判结果、做修正"的完整能力，而不是记住某一个固定公式。