4. 齿轮传动多体动力学建模

齿轮副是多体动力学模型中最“敏感”的元件之一：它的啮合刚度、传递误差和动态激励直接决定系统振动噪声与载荷分配的仿真可信度。本节围绕直齿轮、斜齿轮、锥齿轮三类常见齿轮，说明如何把齿面接触的弹性、制造与安装偏差、时变啮合特性转化为多体动力学可用的刚度、误差与激励输入。

一、三类齿轮的啮合特性差异

直齿轮的接触线平行于轴线，啮合刚度在单齿对与双齿对交替时呈阶跃式变化，激励以“矩形波”为主，冲击明显；斜齿轮接触线倾斜，啮合过程沿齿宽逐渐进入和退出，刚度变化更平滑，但会产生轴向力；锥齿轮的接触区位于锥面上，齿面相对运动包含滚动与滑动，且节锥顶点存在偏置或相交，啮合刚度与传递误差对安装误差极为敏感。三者的差异决定了建模时刚度函数、误差投影和激励方向的构造方式不同。

二、啮合刚度的参数化建模

啮合刚度不是常数，而是随啮合位置变化的周期函数。工程上常用“单齿刚度 + 重合度”构造：先由齿根弯曲、齿面接触和剪切变形计算单齿对刚度 k_single，再根据重合度 ε 判断同时参与啮合的齿对数，得到时变刚度 k(t)。直齿轮的重合度通常在 1.2～1.8 之间，因此 k(t) 在 k_single 与 2k_single 之间跳变；斜齿轮因接触线倾斜，重合度可达 2.5 以上，k(t) 更接近多段平滑曲线。

在多体动力学软件中，常用两种方式实现：一是用样条曲线直接输入 k(t) 的周期数据；二是用“齿轮接触力”元件，由软件根据齿廓几何实时计算。前者计算快、适合参数研究；后者精度高、适合详细校核。示例参数：模数 m=3 mm，齿数 z1=20、z2=40，齿宽 B=30 mm，直齿轮单齿刚度取 k_single=1.4×10^8 N/m，重合度 ε=1.6，则 k(t) 在 1.4×10^8 与 2.8×10^8 N/m 之间交替，周期为啮合周期 T=60/(n·z1)，当转速 n=1500 r/min 时 T≈2 ms。

三、传递误差的构造与投影

传递误差 TE 定义为从动轮实际转角与理想转角之差，折算到啮合线上的等效位移。它由齿距偏差、齿形误差、安装偏心、轴系变形等组成。建模时通常把 TE 处理为沿啮合线方向的位移激励 e(t)，其幅值由实测或标准公差确定。例如某直齿轮副的齿距累积偏差 F_p=25 μm，齿形总偏差 F_α=18 μm，则 TE 的峰峰值可取 30～40 μm，按啮合频率展开为傅里叶级数后输入。

斜齿轮的 TE 还需考虑螺旋角误差和轴向窜动，误差沿接触线分布，不能简单按直齿轮方式投影；锥齿轮的 TE 则要沿节锥面分解为轴向、径向和切向分量，安装距偏差 0.05 mm 就可能使 TE 显著增大。因此，三类齿轮的误差投影矩阵不同，建模时必须明确坐标系和投影方向。

四、动态激励的合成与施加

动态啮合力由弹性力和阻尼力组成：F_dyn = k(t)·e(t) + c·v_rel，其中 c 为啮合阻尼，通常取 0.03～0.1 倍的临界阻尼，v_rel 为沿啮合线的相对速度。直齿轮的激励方向沿啮合线，无轴向分量；斜齿轮需分解为切向、轴向和径向三个分量，轴向力 F_a=F_t·tanβ；锥齿轮则产生三向力和附加弯矩，需在节点处施加力和力矩。

五、建模检查与结果判读

完成建模后，可从三个层面检查：一是刚度曲线是否与重合度一致，直齿轮应出现明显阶跃，斜齿轮应平滑过渡；二是传递误差的幅值和频率是否与齿轮精度等级匹配，例如 6 级精度齿轮的 TE 通常在 10～30 μm；三是动态啮合力的频谱是否出现啮合频率及其边频带。若直齿轮仿真中未出现阶跃刚度，或斜齿轮未出现轴向力，说明模型构造有误，应回到刚度函数和投影方向重新检查。

实际工程中，建议先用参数化刚度与误差快速跑通模型，确认激励频率和量级合理后，再切换到详细齿廓接触模型进行校核。这样既能保证多体动力学求解效率，又能兼顾齿轮啮合细节的仿真精度。