4. 齿轮传动误差与啮合冲击建模

齿轮传动系统在理想情况下应保持恒定的瞬时传动比，但制造误差、安装偏差、轮齿受载变形以及齿面磨损等因素，会使实际啮合状态偏离理想位置，从而产生传递误差与啮合冲击。这些误差与冲击是齿轮系统振动噪声的主要激励源，也是动力学建模中必须准确描述的关键输入。本节围绕静态传递误差、动态传递误差、齿距偏差、齿形偏差、啮入啮出冲击以及误差激励的数学描述展开，建立从误差物理量到动力学激励的完整映射关系。

一、静态传递误差与动态传递误差

传递误差（Transmission Error, TE）定义为从动轮实际角位置与理想角位置之差，沿啮合线方向折算为线位移量，单位为微米（μm）。其数学表达为：

TE = r_b2 · θ_2 - r_b1 · θ_1 / i_12

其中 r_b1、r_b2 分别为主、从动轮基圆半径，θ_1、θ_2 为对应转角，i_12 为理想传动比。静态传递误差（STE）是在准静态、轻载或空载条件下测得的传递误差，主要反映齿轮副的几何制造精度，包括齿距偏差、齿形偏差、基节偏差等综合作用。动态传递误差（DTE）则是在实际载荷与转速下，考虑轮齿弹性变形、轴承刚度、轴系扭振等因素后的传递误差，其频谱中除包含静态误差的基频与谐波外，还会出现由系统固有频率调制产生的边频成分。

在建模实践中，STE 通常作为系统的位移激励输入，而 DTE 是系统响应的一部分。二者关系可理解为：STE 是激励源，DTE 是激励经过齿轮-轴-轴承系统传递后的响应表现。若将齿轮副简化为单自由度扭转模型，则 DTE 的动力学方程可写为：

m_eq · ẍ + c_m · ẋ + k_m(t) · x = F_ext + k_m(t) · e_ste(t)

式中 m_eq 为等效质量，c_m 为啮合阻尼，k_m(t) 为时变啮合刚度，e_ste(t) 为静态传递误差激励，x 为动态传递误差。该式表明，STE 通过时变啮合刚度调制后进入系统，成为动态响应的直接激励。

二、齿距偏差与齿形偏差

齿距偏差（f_pt）是指实际齿距与理论齿距之差，分为单齿距偏差和累积齿距偏差。单齿距偏差影响单对齿啮合时的瞬时传动比波动，累积齿距偏差则影响齿轮旋转一周内的低频误差分量。在建模中，齿距偏差通常以转角位置的函数形式引入，表现为啮合线方向的周期性位移扰动。对于齿数为 z 的齿轮，第 k 个齿的齿距偏差可表示为 f_pt(k)，其离散傅里叶变换后得到各阶谐波幅值，其中一阶谐波对应齿轮旋转频率，高阶谐波对应齿频及其倍频。

齿形偏差（f_fα）是指实际齿廓偏离理论渐开线的量，包括齿廓总偏差、齿廓形状偏差和齿廓倾斜偏差。齿形偏差直接影响啮合刚度的变化规律和啮入啮出时刻的接触状态。在动力学建模中，齿形偏差通常转化为沿啮合线方向的等效位移激励，与齿距偏差共同构成静态传递误差的几何分量。实际工程中，齿形偏差常通过齿廓修形来主动控制，例如齿顶修缘和齿根挖根，以改善啮入啮出冲击。

三、啮入啮出冲击

啮入冲击发生在轮齿进入啮合的瞬间。由于齿距偏差和受载变形，主动轮齿顶与从动轮齿根之间可能存在间隙或过盈，导致接触瞬间产生速度突变，形成冲击力。啮出冲击则发生在轮齿脱离啮合时，由于载荷突然转移至下一对齿，造成接触力突变。冲击力的大小与啮合线速度、齿距偏差、齿形偏差及啮合刚度突变程度相关，可近似表示为：

F_impact = k_m · Δx + c_m · Δv

其中 Δx 为啮入时刻的间隙或过盈量，Δv 为接触瞬间的相对速度差。在时域仿真中，啮入啮出冲击表现为啮合力曲线上的尖峰，其频谱覆盖较宽频带，容易激发齿轮箱体高频模态。

四、误差激励的数学描述

将上述误差源统一纳入动力学方程，静态传递误差 e_ste(t) 可展开为傅里叶级数形式：

e_ste(t) = e_0 + Σ [A_n · cos(n·ω_m·t + φ_n)]

式中 e_0 为平均误差分量，ω_m 为啮合角频率，A_n 和 φ_n 分别为第 n 阶谐波的幅值和相位。对于齿距偏差和齿形偏差，其谐波阶次与齿数相关，主要能量集中在啮合频率及其倍频处。啮入啮出冲击则作为瞬态激励叠加在稳态误差激励之上，在时域中表现为周期性脉冲序列，在频域中表现为宽带成分。

在仿真实现中，通常将 e_ste(t) 以位移激励形式施加于啮合线方向，通过时变啮合刚度 k_m(t) 调制后进入系统。k_m(t) 本身也随啮合位置周期性变化，其变化频率与啮合频率一致。因此，误差激励与刚度激励在频域上存在耦合，会产生丰富的边频成分。建模时需注意采样频率应至少覆盖 5 倍以上最高关注频率，以保证冲击脉冲的时域分辨率。

五、教学示例：单级直齿轮副误差激励设置

以模数 m=3 mm、齿数 z1=25、z2=40、压力角 α=20°、齿宽 b=20 mm 的单级直齿轮副为例。假设主动轮转速 n1=1500 r/min，则啮合频率 f_m = z1·n1/60 = 625 Hz。设定静态传递误差一阶幅值 A1=2.5 μm，二阶幅值 A2=1.2 μm，三阶幅值 A3=0.6 μm，相位均为 0。啮合阻尼比取 ζ=0.08，时变啮合刚度均值 k_avg=2.5×10^8 N/m，刚度波动幅值取均值的 15%。

在仿真中，将上述参数代入动力学方程，采用四阶龙格-库塔法求解，时间步长取 1×10^-6 s，仿真时长 0.1 s。观察动态传递误差的时域波形和频谱，可发现：时域波形中除基频周期成分外，在啮入时刻存在明显尖峰；频谱中除 625 Hz 及其倍频外，在 625±f_n 处出现边频带，其中 f_n 为系统固有频率。通过对比不同齿距偏差和齿形偏差组合下的仿真结果，可定量评估各误差源对动态传递误差的贡献占比，为齿轮精度等级选择和修形方案设计提供依据。