3. 齿轮啮合刚度激励建模
齿轮传动之所以会产生振动和噪声，最根本的内部激励源之一就是时变啮合刚度。当一对轮齿从进入啮合到退出啮合，参与承载的齿对数、接触位置和齿面弹性变形量都在不断变化，导致系统在恒定外载荷下也会产生周期性的动态响应。理解并准确建模这一刚度变化规律，是后续动力学仿真能否反映真实振动特征的关键前提。
3.1 时变啮合刚度的物理来源
直齿轮副在啮合过程中，重合度通常在1到2之间。这意味着在大部分时间里只有一对齿承载，而在交替区间内会出现两对齿同时接触。单对齿承载时，轮齿的弯曲变形、剪切变形、接触变形和齿根基础变形共同决定了该啮合位置的等效刚度；当第二对齿进入接触后，载荷被两对齿分担，系统整体刚度显著上升。因此，刚度曲线呈现出以啮合周期为单位的周期性波动，其波形与重合度、齿廓修形和载荷分布密切相关。
斜齿轮的情况更为复杂。由于齿宽方向上的螺旋角，接触线从齿顶逐渐扫向齿根，啮合刚度沿轴向呈连续变化，其波动幅度通常小于直齿轮，但刚度曲线的形状不再是简单的阶跃式变化，而是带有明显的过渡段。
3.2 势能法：从变形能反推刚度
势能法是一种半解析方法，其核心思想是将轮齿视为变截面悬臂梁，分别计算弯曲势能、剪切势能、轴向压缩势能和接触势能，再通过能量守恒关系得到单齿刚度。对于一对啮合齿，总势能包括两齿各自的变形能和赫兹接触能：
U_total = U_b1 + U_s1 + U_a1 + U_b2 + U_s2 + U_a2 + U_h其中下标 b、s、a、h 分别代表弯曲、剪切、轴向压缩和赫兹接触。对总势能关于法向力求二阶导数，即可得到该啮合位置的等效刚度。势能法的优势在于计算速度快、物理意义明确，适合在参数化研究中快速评估齿宽、模数、压力角等参数对刚度的影响。但其精度依赖于梁模型的假设，对于大齿宽或复杂齿廓的齿轮，误差会有所增大。
3.3 有限元法：高精度但代价较高
有限元法通过直接离散轮齿实体，施加单位法向力后提取接触点的法向位移，再取倒数得到刚度。实际操作中，通常建立单齿或三齿模型，在齿面接触点施加单位载荷，通过求解静力学方程获得变形量。有限元法的优势在于能够考虑齿根过渡曲线、轮缘厚度、腹板结构等细节，精度最高。但每计算一个啮合位置就需要一次完整的静力求解，若要在动力学仿真中实时调用，计算成本难以接受。因此，有限元法更多用于校验解析模型或生成刚度曲线的基准数据。
3.4 解析近似公式与工程修正
工程中广泛使用的是基于ISO 6336或AGMA标准的解析近似公式。这类公式将单齿刚度表示为齿数、模数、齿宽、压力角和变位系数的函数，并通过经验系数修正。例如，ISO 6336中给出的单齿刚度近似为：
c' = c_th * C_M * C_R * C_B * cos(β)其中 c_th 为理论刚度，C_M 为材料修正系数，C_R 为轮缘修正系数，C_B 为齿宽修正系数，β 为螺旋角。这类公式的优点是计算极快，适合嵌入到动力学方程中作为时变参数调用。但需要注意，不同标准给出的系数适用范围不同，超出适用范围时误差可能达到20%以上。
3.5 刚度曲线的周期特性与参数影响
时变啮合刚度的周期特性可以用啮合周期和重合度两个核心参数来描述。对于直齿轮，刚度曲线在一个啮合周期内呈现“单齿—双齿—单齿”的交替模式，其基频为啮合频率，即轴频乘以齿数。斜齿轮由于接触线连续变化，刚度波动幅值较小，但频谱中会出现以轴向重合度为特征的边频成分。
参数影响方面，齿宽增加会提高平均刚度但也会增大波动幅值；重合度增大则使双齿承载区间变长，刚度曲线更平滑；齿廓修形可以主动改变刚度曲线的形状，在高速轻载工况下常用来降低啮入冲击。此外，载荷大小也会影响实际接触齿对数，轻载时可能只有部分齿面接触，导致刚度曲线偏离理论值。
上图对比了直齿轮与斜齿轮的典型刚度曲线形态。直齿轮的刚度曲线呈明显的阶跃特征，在单齿区和双齿区之间跳变；斜齿轮由于螺旋角的存在，接触线连续变化，刚度曲线呈现平滑过渡。图中还标注了重合度对应的双齿承载区间，以及齿宽、重合度和齿廓修形对曲线形态的影响方向。
3.6 建模方法的选择与组合策略
在实际动力学仿真中，选择哪种刚度建模方法取决于分析目的和计算资源。如果目标是快速评估系统固有频率和主共振区间，势能法或解析近似公式已经足够；如果需要精确复现齿根应力或接触疲劳，有限元法更为合适；如果要在时域积分中实时更新刚度，则通常采用“有限元生成基准曲线 + 解析公式插值”的组合策略。
一个实用的做法是：先用有限元法计算一个啮合周期内若干离散位置的刚度值，拟合出刚度随啮合位置变化的傅里叶级数，再将前几阶谐波系数代入动力学方程。这样既保留了有限元法的精度，又大幅降低了在线计算量。对于斜齿轮，还需要注意轴向重合度带来的相位差，不同齿对之间的刚度曲线存在时间延迟，不能简单叠加。
总结来说，时变啮合刚度建模的核心在于：明确物理来源，选择合适的计算方法，并正确描述刚度曲线的周期特性与参数敏感性。只有把刚度激励建模准确，后续的动力学仿真才能真实反映齿轮系统的振动行为。