ANSYS 齿轮仿真：直齿圆柱齿轮齿根弯曲与赫兹接触应力随转矩变化

项目过程

本项目为某减速机厂的汽车变速箱齿轮副强度校核项目，需求是评估直齿圆柱齿轮副在 0~300 N·m 转矩范围内的齿根弯曲疲劳与齿面接触疲劳裕度，并输出最大接触压力分布以评估点蚀风险。客户给定：齿轮模数 m=3 mm、齿数 z=23、面宽 b=28 mm、材料 16MnCr5（渗碳淬火、表面硬度 60 HRC、心部 30 HRC）、压力角 20°、齿顶高系数 1.0、变位系数 0。

本项目采用 ANSYS Workbench Mechanical + ANSYS Fatigue 静力 + 疲劳叠加。齿轮副单齿啮合区用三维接触模型：主动轮刚性、从动轮柔性（50 万四面体单元 Solid187），接触 *CONTACT_TARGE170 + *CONTACT_CONTA174，法向接触刚度因子 0.1，摩擦系数 μ=0.05。载荷通过 *LOAD_BOLT_PRE_TENSION 在从动轮内孔施加转矩 T，扫描 0~300 N·m，每 50 N·m 一个工况，输出齿根最大弯曲应力 σf（ISO 6336 标准路径）与赫兹接触应力 σH。

案例分析

两条应力曲线随转矩 T 增长呈现典型非线性特征：齿根弯曲应力 σf 从 0 上升到 280 MPa（T=300 N·m），曲线接近线性（斜率 0.93 MPa/N·m），这是因为弯曲应力在弹性范围内与载荷成正比；赫兹接触应力 σH 从 0 单调上升到 620 MPa（T=300 N·m），但曲线呈现”幂次 < 1″的凹形——这是 Hertz 接触理论中 σH∝(F/L)^0.5 的体现（F 是法向力、L 是接触线长度）。相同载荷下 σH 比 σf 大约 2.2 倍，与 16MnCr5 渗碳淬火齿轮的工程经验吻合。

接触压力沿啮合线分布呈现完美 Hertz 半椭圆：T=300 N·m 时最大接触压力 p_max=620 MPa，接触半宽 a=0.42 mm（齿轮法向），沿啮合线方向呈左右对称分布，载荷作用中心（p_max 处）位于啮合线中点，载荷边缘压力快速线性下降到零。最大接触压力 620 MPa 已接近 60 HRC 钢的接触疲劳极限（约 1500 MPa），安全裕度 2.4。

总结教训

本项目最容易踩的坑是接触刚度因子 K 的设定。客户最初用 ANSYS 默认 K=1.0（无罚因子），结果接触穿透达 0.05 mm，应力偏低 25%。本项目改用 K=0.1 + 接触刚度 *FKN=200 后与 Hertz 解析解误差 <3%。

另一坑是齿根应力提取的”应力奇异点”——单齿啮合区末端齿根圆角应力集中最严重，但 FEA 在圆角根部会出现伪峰值（奇异点）。本项目按 ISO 6336 在齿根 30°切线方向提取 σf，避开奇异点。

项目经验

经验一：σf/σH≈2.2 是直齿圆柱齿轮的典型比值；如果换成斜齿（β=15°），σH 会下降 10%、σf 下降 5%，比值变为 2.0；如果换成锥齿轮，比值进一步降到 1.8。

经验二：T=300 N·m 工况下 σf=280 MPa 仍 < 材料允许弯曲应力 600 MPa（60 HRC），FOS≈2.1，可承受短期 2 倍过载；如果客户长期在 350 N·m 运行，需升级为 17CrNiMo6。

经验三：接触压力半宽 a=0.42 mm 是判断齿面润滑的关键尺度——润滑油膜厚度需 ≥ 0.42 μm 才能形成弹流润滑，常规 68# 工业齿轮油膜厚 0.6~1.0 μm，满足。

有限元仿真