3. 齿轮传动动力学基础
齿轮传动在理想情况下应像一对纯滚动节圆一样平稳传递运动，但真实齿轮系统在运转中始终伴随振动与噪声。理解这些振动的来源、传递路径和响应特征，是建立多体动力学模型、设置激励与判断仿真结果是否合理的基础。本节围绕四个相互关联的核心概念展开：振动来源、时变啮合刚度、传递误差、动态啮合力，并说明它们如何共同决定系统的响应特征。
一、齿轮系统振动的来源
齿轮系统的振动并非单一原因造成，而是内部激励与外部激励叠加的结果。内部激励源于齿轮副本身的结构与啮合特性，主要包括：
- 时变啮合刚度：啮合齿对数随转角周期变化，单齿对与双齿对交替承载，使综合啮合刚度呈周期性波动，直接产生参数激励。
- 传递误差：由于齿形误差、基节偏差、受载变形等，从动轮实际转角与理想转角之间存在偏差，形成位移激励。
- 啮入啮出冲击：轮齿进入和退出啮合时载荷突变，产生高频冲击成分。
- 齿面摩擦与阻尼：滑动速度方向在节点处反向，摩擦力方向随之改变，也会激发振动。
外部激励则来自原动机的扭矩波动、负载变化、轴承与箱体的支承特性，以及轴系不对中、不平衡等。在多体动力学建模中，通常把时变啮合刚度和传递误差作为内部激励输入，把外部扭矩和转速作为边界条件，二者共同作用后由系统动力学方程求解响应。
二、时变啮合刚度
啮合刚度定义为使一对轮齿在啮合线方向产生单位弹性变形所需的力，单位为 N/m。它随啮合位置变化：在单齿啮合区，仅一对齿承载，刚度较低；进入双齿啮合区后，两对齿并联承载，综合刚度升高。因此，综合啮合刚度在一个啮合周期内呈近似矩形波或梯形波变化。
工程上常用石川公式、韦伯公式或有限元法计算单齿刚度，再按啮合齿对数叠加得到综合刚度。以某模数 m=3 mm、齿数 z₁=20、z₂=40 的直齿轮副为例，假设单齿对平均刚度为 1.4×10⁸ N/m，双齿区两对齿并联后综合刚度约为 2.6×10⁸ N/m，则综合刚度波动幅值约为 1.2×10⁸ N/m。这个波动就是参数激励的来源，其基频为啮合频率 f_m = n·z/60（n 为转速，z 为齿数）。
在多体动力学中，时变啮合刚度通常以随转角变化的函数或查表形式施加到啮合力元上。若忽略其时变性而取平均值，仿真结果会丢失啮合频率及其倍频处的振动成分，导致响应谱过于干净，与实际不符。
三、传递误差
传递误差（Transmission Error, TE）定义为从动轮实际转角与理想转角之差，折算到啮合线方向后以位移量表示，单位为 μm。它由两部分组成：一是制造与安装误差，如齿距偏差、齿形误差、偏心；二是受载变形，包括轮齿弯曲、接触变形和轴系变形。传递误差相当于在啮合线上施加了一个位移激励，即使外载恒定，也会激起系统振动。
传递误差与动态啮合力通过啮合刚度耦合：动态啮合力近似等于综合啮合刚度乘以传递误差。因此，减小传递误差或降低其波动，是降低齿轮振动最直接的途径之一。实际中常通过齿廓修形（如齿顶修缘、齿根修形）使受载后的传递误差趋于平缓，从而抑制啮合冲击。
四、动态啮合力与系统响应特征
动态啮合力是齿轮副在运转过程中实际传递的力，等于静态载荷与动态分量之和。其动态分量由时变啮合刚度和传递误差共同激励产生，并受系统惯量、阻尼和支承刚度调制。动态啮合力过大，会加剧齿面疲劳、点蚀和噪声。
系统响应特征主要体现在以下几个方面：
- 时域响应：啮合力或振动加速度呈现明显的周期性波动，波动周期对应啮合周期，包络可能受轴频调制。
- 频域响应：频谱中啮合频率 f_m 及其谐波 2f_m、3f_m 处出现峰值，两侧可能出现以轴频为间隔的边频带，边频反映调制与故障。
- 共振：当啮合频率或其谐波接近系统固有频率时，响应幅值显著放大，需通过调整转速或修改支承刚度避开。
- 参数稳定性：时变啮合刚度可能导致参数振动，在特定转速区间出现不稳定，表现为响应突然增大。
下面用一张信号链与参数关系图，把振动来源、时变啮合刚度、传递误差、动态啮合力和系统响应之间的因果关系集中呈现出来。
五、建模与结果判读要点
在建立齿轮传动多体动力学模型时，建议按以下顺序处理本节涉及的要素。首先确定齿轮副的基本参数：模数、齿数、压力角、齿宽、转速和负载扭矩。其次计算或导入时变啮合刚度曲线，确保其基频与啮合频率一致。然后根据齿形误差和受载变形构造传递误差激励，可简化为正弦叠加或实测谱。最后将二者耦合到啮合力元中，设置合理的阻尼比（齿轮副通常取 0.03～0.1）和支承刚度。
仿真完成后，判读结果时重点关注三点：一是时域啮合力是否呈现与啮合周期一致的波动；二是频谱中 f_m 及其谐波幅值是否合理，边频带是否与轴频对应；三是扫频或变转速工况下是否出现共振峰或参数不稳定区。若频谱过于干净，应检查是否忽略了时变刚度或传递误差；若响应幅值异常大，应检查阻尼和支承参数是否合理。
综上，齿轮系统振动是内部激励与外部激励共同作用的结果，时变啮合刚度和传递误差是两条核心内部激励路径，动态啮合力是二者耦合后的直接表现，而系统响应特征则通过时域波动、频域啮合频率及边频带、共振与参数稳定性来体现。掌握这四者之间的因果关系，才能在多体动力学仿真中正确设置激励、合理解读结果，并为后续的减振降噪设计提供依据。