基于MASTA的电动拖拉机行星齿轮箱优化仿真
关键词：MASTA；行星齿轮；传动误差；宏几何优化；多目标优化
一、文章简要介绍
电动拖拉机没有发动机噪声做遮掩，齿轮啸叫会直接传进驾驶室，而农用工况负载又重又乱，齿轮箱既要扛得住也要转得静。韩国忠北大学团队在Scientific Reports上发表的工作，针对一台双电机电动拖拉机样机的复合行星齿轮箱做宏几何优化：在齿数比和箱体尺寸都锁死的前提下，用MASTA传动分析软件建出齿轮-轴-轴承耦合模型，拿NSGA-II算法同时压低峰峰值静态传动误差、抬高啮合效率，强度按ISO 6336校核兜底，最后用CRITIC权重法从Pareto前沿里挑出最终方案。优化后全部齿轮副的传动误差降低14%到16%，啮合效率不降反升（DOI：10.1038/s41598-026-43864-3）。
图1：电动拖拉机样机与动力总成。牵引电机与PTO电机经复合行星齿轮箱实现功率分流合成，齿轮箱尺寸与传动比受样机箱体约束。
样机采用牵引电机加PTO电机的双电机架构，通过复合行星齿轮箱实现功率分流与合成。行星齿轮箱由PGS1和PGS2两级组成，传动比固定为7.8，太阳轮15齿、行星轮51齿、齿圈117齿，材料SCM420H，润滑油ISO VG 220。
图2：复合行星齿轮箱CPGT构型。PGS1与PGS2两级行星排串联，太阳轮15齿、行星轮51齿、齿圈117齿，传动比固定7.8。
二、MASTA建模与仿真步骤
步骤一：建模。在MASTA（v13.0.1）里建CPGT整体模型，软件内部实现ISO 6336强度计算与齿面接触分析，行星轮系做耦合求解，太阳轮轴、行星轮轴、齿圈轴上的轴承按刚性约束固定到箱体。
图3：MASTA中建立的CPGT仿真模型。含齿轮、轴与轴承的耦合模型，内置ISO 6336强度计算与齿面接触分析，行星轮系耦合求解。
步骤二：定义目标函数。峰峰值静态传动误差PPSTE按啮合线上的角位移几何定义计算，它是齿轮啸叫的主要激励源；啮合效率GME由耦合行星轮系求解器直接输出。
图4：传动误差几何定义。沿啮合线度量啮合角位移，峰峰值静态传动误差PPSTE是齿轮啸叫的主要激励源。
步骤三：设置工况与约束。载荷来自样机实测的载荷持续分布数据，犁耕、旋耕、行走按4比4比2的比例拆成24个载荷工况，按年作业342小时、寿命8年共2736小时加权。约束条件：总重合度1.2到2.0，宽径比1.0到2.0，ISO 6336弯曲安全系数不低于1.2、接触安全系数不低于1.0。
步骤四：NSGA-II优化。设计变量为两级行星轮系的法向模数、压力角、螺旋角和齿宽，种群与代数各设100，交叉与变异概率0.5，Python脚本自动调MASTA完成每次个体评价。
图5：NSGA-II宏几何优化流程。设计变量为模数、压力角、螺旋角与齿宽，个体评价由Python脚本自动调用MASTA完成。
三、关键结果与发现
PGS1的Pareto前沿显示，太阳轮-行星轮副PPSTE在2.16到4.75微米之间、效率最高到99.39%，且传动误差与啮合效率呈明显正相关——想压低误差就得在效率上让一点步，回归拟合优度达到0.963。
图6：PGS1的Pareto前沿。太阳轮-行星轮副与行星轮-齿圈副的传动误差和效率呈正相关，回归优度0.963，最高效率99.39%。
PGS2的规律更稳健：行星轮-齿圈副平均效率99.24%、波动只有0.25个百分点，对制造公差和变载荷工况更不敏感。CRITIC法加权排序后，PGS2第一名方案把两对齿轮副的传动误差分别降了15.1%和14.6%，效率还提高了0.16个百分点；最终方案确定PGS1模数3毫米、压力角19.2度、螺旋角15.1度，PGS2模数3.5毫米、压力角24.1度、螺旋角15.4度。
图7：PGS2的Pareto前沿。行星轮-齿圈副平均效率99.24%且波动最小，CRITIC排序第一名方案误差降14.6%、效率升0.16个百分点。
