针对高承载工况下锥齿轮接触印痕偏移与性能退化问题，构建承载变形-误差耦合模型，提出基于位移场补偿的齿面拓扑修形方法。采用多目标优化算法确定修形参数，并在功率封闭式试验台上开展阶梯加载验证。结果表明，优化方案可有效降低接触中心偏移率与载荷分布离散度，抑制高载区振动与温升增长趋势。研究建立了承载变形与修形参数之间的定量映射关系，为高载齿轮接触控制提供了理论依据。

在高承载工况下，锥齿轮的啮合错位由支撑系统刚度特性及齿轮本体弹性变形共同作用。支撑系统刚度矩阵由轴承和箱体刚度构成，影响齿轮副轴线的平行度与相交度。当载荷变化时，支撑系统产生的径向位移和角度偏差导致齿轮安装误差发生动态变化。齿轮本体的弹性变形由作用扭矩引起，扭转角度随载荷增大而增大，进而影响啮合点位置，改变接触线的实际空间分布。假设齿轮副受到轴向载荷与扭矩作用，错位量可定量描述为

式中：δ为错位量，m；T为作用扭矩，N·m；b为齿宽，m；G为材料剪切模量，Pa；J为齿轮本体的极惯性矩，m⁴；K_s为支撑系统刚度矩阵；F为外载荷向量，N；e_δ为错位方向单位向量。式（1）揭示了载荷与弹性变形对啮合位置偏移的定量关系，为承载变形引发的啮合错位机制提供了基础分析框架。

重载条件下，齿面弯曲变形与接触弹性压缩叠加，啮合点随载荷增加发生连续偏移。有限元结果表明，扭矩由0.5T_n（额定扭矩）增至1.1T_n时，接触椭圆中心沿齿长方向向小端迁移0.2～0.5mm，沿齿高方向向齿顶侧偏移，长轴旋转1°～3°。齿根挠度增大致局部载荷集中，接触应力提高15%～25%，接触带呈现不对称扩展。在超额定载荷条件下出现边缘接触，印痕由椭圆向拉伸形态演化，中心偏移量与支撑系统等效刚度呈负相关，体现载荷、变形与接触分布之间的耦合特征。 

接触印痕形态反映齿面载荷分布特征，并与失效模式相对应。当印痕宽度小于齿宽的40%时，单位面积接触应力显著升高，局部应力超过材料接触疲劳极限，易产生次表层裂纹并形成点蚀。当印痕沿齿高方向上移且距齿顶边缘小于0.5mm时，接触线偏离有效啮合区，齿顶附加弯曲应力增大，诱发崩角或断齿。当印痕向小端或大端偏移超过齿长的10%时，载荷呈单侧集中，齿根应力梯度增加，促进裂纹萌生。印痕由居中椭圆向边缘拉伸形态演化，表征接触状态由均布向偏载转变，是判定潜在失效的重要几何指标。 

承载印痕优化需在可制造误差范围内建立稳定的分析基准。结合装配精度与加工能力，轴向安装距误差∆a控制在±0.03mm以内，轴线交角误差控制在±0.02°以内，轴承径向游隙折算的等效径向偏移量不超过0.02mm。齿轮加工精度按ISO 1328 6级控制，齿向误差不大于0.015mm，齿廓误差不大于0.010mm。上述误差区间构成承载印痕分析的扰动边界，用于描述实际制造与装配波动对接触位置稳定性的影响范围。 

敏感性分析选取大轮滚检转角误差∆θ与小轮安装距误差∆a为主变量，∆θ取值范围为±0.05°、步长为0.01°，∆a取值范围为±0.05mm、步长为0.01mm，构建二维误差扰动网格。通过齿面接触分析获得接触中心在齿宽方向的偏移量∆x，并考虑支撑系统等效刚度k_θ对误差放大的耦合作用，建立归一化误差响应函数为

式中：ε为印痕中心偏移率；r为节锥半径，mm；k_θ为支撑系统等效刚度，N·m·rad^-1；α₁、α₂、α₃为灵敏度系数。式（2）综合考虑几何误差与载荷变形对接触中心的耦合影响，实现误差-载荷-刚度的统一量化。当ε超过0.08时，接触位置对误差扰动呈高敏感状态，可作为后续修形参数设定的判定基准。 

承载变形预控修形以T_n下的齿面弹性位移场为设计基准，目标是在加载后形成期望接触区域。齿长方向修形需补偿齿轮副在支撑柔度与齿体弯曲共同作用下产生的挠曲分布。基于齿宽方向弯曲线可近似为二次函数，修形曲线采用连续二阶可导形式，使修形梯度与载荷引起的挠度梯度相匹配。中心鼓形量∆c₀由额定载荷下齿宽中心挠度δ_m反算确定，并考虑材料弹性模量E、齿根截面惯性矩I及受力臂长l的综合影响。

齿高方向修形用于抵消弯曲变形引起的基节误差与啮合线投影偏差。齿顶区域的法向位移δ_t与基圆半径r_b及法向模数m_n相关，压力角修形量需满足加载后啮合线位置恢复至理论设计线。综合齿长与齿高方向的耦合补偿关系，建立承载预补偿模型为

式中：C(x,y)为齿面任意点(x,y)的修形量，mm；δ(x,y,T_n)为额定扭矩T_n作用下该点的弹性位移，mm；δ_ref为期望接触基准位移，mm。式（3）将齿面空间位移场直接映射为修形量，实现基于承载变形的场函数级补偿设计，使修形参数由标量经验值提升为分布函数形式。

修形方案解算以齿面空间修形函数参数向量p=[∆c₀, k_c,∆α,k_α]为设计变量，其中∆c₀为齿长方向中心鼓形量，k_c为齿长分布系数，∆α为齿高方向压力角修形量（mm），k_α为齿高梯度系数。变量取值范围依据制造能力与误差敏感区间设定：∆c₀∈[0,0.030]mm，∆α∈[0,0.020]mm。优化目标围绕承载接触状态的稳定性与均匀性构建，采用接触中心偏移率、载荷分布离散度及传动误差波动量作为表征量。基于齿面接触分析模型获得额定载荷下的接触应力场与位移场，建立统一标度的多目标泛函为

式中：Φ(p)为多目标优化函数；A_c为实际接触区域；q(x,y)为局部接触载荷密度，N·mm^-2；为平均载荷密度；∆x为接触中心偏移量，mm；τ(x,y)为局部传动误差，μm；τ₀为基准误差；λ₁、λ₂为权重系数。约束条件包括最大接触应力不超过材料允许值、修形曲面连续可导且曲率受限、印痕长度占比处于设定区间。求解采用基于齿面接触分析耦合有限元的迭代优化框架，结合遗传算法进行全局搜索，再通过梯度法进行局部收敛修正，实现修形参数在载荷-误差耦合场中的最优匹配。 

试验采用功率封闭式锥齿轮加载试验台，变频电机驱动，额定功率为45W，最高转速为3000r·min^-1。液压加载系统最大输出扭矩为1500N·m，扭矩测量精度为±0.5%。齿轮副安装于刚性箱体内，轴承预紧力按设计值设定，轴向安装距控制在±0.01mm以内。按阶梯方式加载，依次施加0.25T_n、0.50T_n、0.75T_n、1.00T_n和1.10T_n的载荷，每级稳态运行10min后采集数据。接触印痕采用红丹粉显影，高清摄影的空间分辨率不高于20μm（每个像素对应的实际尺寸不超过20μm），通过图像处理提取印痕长度、宽度及中心偏移量。试验参数与优化模型一致，用于形成承载工况下的实测对比数据。 

依据第3.1节的加载程序，对原始设计（轻载修形）与优化设计（承载修形）在不同扭矩等级下的接触印痕进行量化分析。齿宽b为40mm，印痕长度占比η₁与宽度占比η_b由图像识别算法提取，中心偏移量∆x由印痕质心与齿宽中心线之间的距离计算，归一化得到偏移率ε=∆x/b。不同载荷等级下的关键参数统计结果见表1。

表1 不同载荷下印痕参数对比

由表1可见，低载区间两种设计的印痕参数差异有限，中心偏移率均低于0.01。载荷升至0.75T_n后，原始设计的∆x明显增大，偏移率接近0.04，接触带出现单侧扩展趋势。达到1.00T_n及以上时，原始设计偏移率超过0.08，印痕由中部向小端及齿顶迁移并伴随边缘接触。优化设计在相同载荷区间内，偏移率始终控制在0.05以内，接触区域扩展形态保持对称。载荷提高至1.10T_n时，两种设计的印痕差异进一步扩大，原始设计表现出明显的接触集中与边界逼近特征，优化设计仍维持相对稳定的接触分布状态。 

在印痕形态对比基础上，对两种设计在额定及超额定载荷下的传动效率、温升、振动加速度及跑合后齿面状态进行综合评价。效率通过输入输出扭矩差计算，温升由齿轮箱稳定运行30min后的油温与环境温度差值获得；振动信号采用加速度传感器测量，采样频率为10kHz，提取啮合频率及其2倍频处的幅值。测试结果见表2。

表2 综合性能指标对比

由表2可见，原始设计在额定及超额定载荷下表现出效率下降与温升升高趋势，啮合频率及倍频处振幅值明显增大，跑合后齿面局部出现非均匀磨痕。优化设计在相同工况下各项指标变化幅度较小，振动频谱能量分布更为集中，齿面状态保持一致性。综合结果表明，接触印痕的空间分布特征与效率损失、振动响应及表面磨损形态存在一致性关联，高载下印痕偏移与边缘接触程度直接反映在动力学与热特性指标上。

本文建立了承载变形诱发啮合错位与印痕迁移的分析模型，明确了载荷、误差与接触位置之间的耦合关系。构建基于位移场补偿的齿面拓扑修形函数，将额定载荷下的弹性位移映射为修形分布，实现几何参数的定量预控。通过多目标优化求解修形参数组合，在约束接触应力与印痕占比条件下获得稳定解。试验结果表明，优化设计在高载区间保持较低的接触中心偏移率与振动幅值，效率与温升变化趋势趋于平稳。研究形成了承载印痕控制的系统方法框架。未来可进一步结合热弹耦合与磨损演化模型，拓展长期服役条件下的修形适应性研究。

END