3. 动力学建模的力学基础

在上一节中，我们已经把曲柄连杆机构的几何运动关系梳理清楚：曲柄转角、连杆摆角、活塞位移之间存在确定的运动学约束。但运动学只回答“怎么动”，动力学还要回答“为什么这样动”——驱动力矩如何被气体力、惯性力、摩擦阻力共同消耗，各铰链处到底承受多大的力。要回答这些问题，必须先把真实机构抽象成可计算的力学模型，而这一抽象过程依赖四个核心概念：牛顿-欧拉方程、达朗贝尔原理、约束、自由度。

一、牛顿-欧拉方程：刚体动力学的两条基本方程

曲柄连杆机构中的曲轴、连杆、活塞在低速工况下可近似视为刚体。刚体动力学由两条方程共同描述：

• 牛顿方程（平动）：合外力等于质量乘以质心加速度，ΣF = m·a_c。它决定质心的平动规律。
• 欧拉方程（转动）：合外力矩等于转动惯量矩阵乘以角加速度，ΣM_c = I_c·α + ω × (I_c·ω)。它决定刚体绕质心的转动规律。

对平面机构，欧拉方程可简化为 ΣM_c = I_c·α，其中 I_c 是绕质心且垂直于运动平面的转动惯量。以某单缸机为例，设曲柄质量 m_cr = 2.4 kg、绕质心转动惯量 I_cr = 0.012 kg·m²，连杆质量 m_rod = 1.1 kg、I_rod = 0.008 kg·m²，活塞组质量 m_pis = 0.65 kg。这些参数是后续所有动力学计算的输入，必须与三维模型或实测值一致。

需要特别强调：牛顿-欧拉方程是矢量方程，使用前必须先定义坐标系。本课程统一约定：x 轴沿气缸轴线、指向曲轴中心为正，y 轴水平，转角 θ 从 x 轴正方向逆时针度量。全文与图示均遵循此约定，避免正负号混乱。

二、达朗贝尔原理：把动力学问题“降级”为静力学问题

直接对每个刚体列牛顿-欧拉方程，会得到一组含铰链内力的耦合方程，求解繁琐。达朗贝尔原理提供了一条捷径：把惯性项 -m·a_c 和 -I_c·α 当作“惯性力”和“惯性力矩”加到刚体上，则动力学方程在形式上变成静力平衡方程：

ΣF + (-m·a_c) = 0
ΣM_c + (-I_c·α) = 0

这样，整个机构在任意瞬时都可以按“受力平衡”的方式处理。对曲柄连杆机构，最常用的做法是把连杆质量按质心位置等效分配到曲柄端和活塞端，得到“二质量替代系统”：

• 往复质量 m_rec = m_pis + m_rod·(l_b/l)，随活塞做往复运动；
• 旋转质量 m_rot = m_cr + m_rod·(l_a/l)，随曲柄做旋转运动。

其中 l_a、l_b 分别为连杆质心到小头、大头中心的距离，l = l_a + l_b。以 l = 0.16 m、l_a = 0.055 m 为例，则 m_rec = 0.65 + 1.1×0.105/0.16 ≈ 1.372 kg，m_rot = 2.4 + 1.1×0.055/0.16 ≈ 2.778 kg。往复惯性力 F_rec = -m_rec·a_pis，旋转惯性力 F_rot = m_rot·r·ω²，方向沿曲柄向外。这两项是曲柄连杆机构动力学分析中最主要的惯性载荷来源。

三、约束与自由度：机构能怎么动，由约束决定

自由度描述系统独立运动参数的个数。一个在平面内自由运动的刚体有 3 个自由度（x、y 平动加绕 z 转动）。曲柄连杆机构由曲轴、连杆、活塞三个刚体组成，若完全自由应有 9 个自由度。但实际机构被一系列运动副约束：

• 曲轴与机架之间为转动副，约束 2 个平动自由度，保留 1 个转动自由度；
• 曲柄与连杆之间为转动副，约束 2 个相对平动自由度；
• 连杆与活塞之间为转动副，约束 2 个相对平动自由度；
• 活塞与气缸壁之间为移动副，约束 2 个自由度（1 个平动、1 个转动），保留沿气缸轴线的移动自由度。

按平面机构自由度公式 F = 3n - 2p_L - p_H，n = 3、低副数 p_L = 4、高副数 p_H = 0，得 F = 3×3 - 2×4 = 1。这说明整个机构只有一个独立运动参数，通常取曲柄转角 θ。只要给定 θ 和曲轴角速度 ω，全部构件的位移、速度、加速度都能唯一确定——这正是运动学分析能先行、动力学分析能简化的根本原因。

四、从真实机构到多刚体系统：抽象步骤

把曲柄连杆机构抽象为多刚体系统，可按以下顺序进行：

1. 确定构件与连接关系：曲轴—连杆—活塞—机架，共 4 个物体、4 个运动副。
2. 建立坐标系与广义坐标：取 θ 为广义坐标，其余坐标由约束方程消去。
3. 计算质量与惯量参数：各构件质量、质心位置、转动惯量，必要时做二质量替代。
4. 施加外力与力矩：气缸气体压力 F_g = p·A_pis、往复惯性力、旋转惯性力、负载力矩、摩擦力。
5. 列写动力学方程：对每个刚体用牛顿-欧拉方程，或整体用达朗贝尔原理列平衡方程。
6. 求解与校核：解出铰链反力、曲轴输出力矩，检查量纲、数量级和边界工况是否合理。

以某工况为例：缸径 D = 80 mm，活塞面积 A_pis ≈ 5.027×10⁻³ m²，缸内压力 p = 4.5 MPa，则气体力 F_g ≈ 22.6 kN。若曲柄半径 r = 0.04 m、转速 n = 3000 r/min，则 ω ≈ 314 rad/s，旋转惯性力 F_rot = 2.778×0.04×314² ≈ 10.96 kN。可见惯性力与气体力处于同一量级，绝不能忽略。

五、抽象结果的检验与常见误区

完成抽象后，应做三项基本检验。第一，自由度检验：若算得 F 不等于 1，说明约束数或构件数有误，常见错误是把活塞与气缸壁之间的移动副漏算，或把曲轴与机架的支承重复计入。第二，量纲与数量级检验：气体力约 10¹~10² kN，惯性力随转速平方增长，若算出的铰链反力比气体力小两个数量级，多半是质量参数或加速度符号出错。第三，边界工况检验：在 θ = 0°（上止点）时活塞加速度最大、往复惯性力最大；在 θ = 180°（下止点）时气体力通常较小。若仿真结果与此规律矛盾，应回到约束方程和坐标系定义检查。