加工公差波动对齿轮NVH性能的鲁棒性影响机制与量化分析

1 引言  

在新能源电驱齿轮箱的NVH开发流程中，加工公差的不确定性是影响量产NVH一致性的关键因素。工程实践表明，若在设计前期未对加工公差开展系统的鲁棒性分析，仅依赖经验化的公差分配，会导致齿轮的NVH在批量生产中呈现显著离散性。尤其是在高速、紧凑化设计趋势下，齿轮微观修形、轴承座位置度、装配预紧量等参数对NVH的敏感度进一步提升，若缺乏前期的鲁棒性评估，量产阶段的EOL下线检测中，产品NVH指标的标准差往往会明显偏大，表现为批次间振动指标不稳定、噪声水平差异大、EOL指标难制定等工程问题。因此，在设计阶段开展齿轮NVH鲁棒性分析，提前识别关键公差对EOL一致性的影响规律，具有重要的工程必要性与现实价值。

齿轮NVH鲁棒性分析旨在评估加工与装配公差波动对系统振动噪声性能的影响程度，常用方法主要包括极值法、正交试验法、响应面法及蒙特卡洛模拟法。极值法操作简单，但仅考虑边界组合，难以反映真实统计分布；正交试验与响应面法可减少试验次数，但对非线性较强的齿轮啮合问题拟合精度有限；蒙特卡洛方法基于随机抽样与统计分析，能够直观反映多公差耦合下的性能分布规律，适用性强、结果可靠，是工程上进行公差鲁棒性分析的主流方法。本文聚焦对齿轮啮合影响最显著的关键加工公差，固定单一典型工况，采用蒙特卡洛仿真这一高效的统计分析方法，利用我们公司自研的齿轮NVH工艺鲁棒性分析软件对齿轮NVH离散性和公差敏感度进行分析，重点阐述分析思路、模型构建、仿真实施及鲁棒性评价的完整流程，相关的分析方法可向全公差带与全工况拓展，为同类齿轮NVH鲁棒性分析提供方法参考。

 

图1-1 公司自研齿轮NVH工艺鲁棒性分析软件

2 研究对象与关键加工公差选取

2.1 研究对象

本文以两级平行轴齿轮传动系统为研究对象，如图2.1所示。该系统广泛应用于电驱总成、工业减速机等领域，其NVH性能直接影响整机舒适性与可靠性。研究聚焦一级、二级齿轮副的啮合过程，以传递误差（TE）作为NVH鲁棒性的核心评价指标，通过分析公差波动对TE的影响，量化齿轮NVH鲁棒性水平。

 

图2-1 分析对象：两级平行轴减速器

2.2 关键加工公差选取

以典型的两级减速平行轴电驱齿轮箱为例，其与NVH性能强相关的加工公差变量至少包含三大类，总计53个核心参数，具体如下：

2.2.1 微观修形参数：40个变量

齿轮NVH性能的优劣高度依赖齿面微观修形设计，以消除边缘接触、改善载荷分布并降低啮合冲击。两级平行轴结构通常包含8个齿面（一级主动、一级从动、二级主动、二级从动，各含两个齿面），每个齿面需定义至少5个修形参数，构成微观修形参数体系。具体包括：

齿向鼓形量（Cβ）：用于补偿齿向弯曲与接触变形，实现齿面载荷均匀分布；

齿向斜度（fHβ）：表征齿向线性斜率误差，控制齿向修形梯度；

渐开线鼓形量（Cα）：优化齿形接触分布，减小啮合冲击与传递误差波动；

渐开线斜度（fHα）：表征渐开线线性斜率误差，定义鼓形变化速率；

渐开线齿顶修缘量（Cαa，行业俗称fKo）：在齿顶进行微量切削，避免齿顶干涉与啮合冲击；

 

上述5个参数×8个齿面，共计40个微观修形参数。

 

图2-2 微观修形参数对NVH的影响

2.2.2 轴承座位置度参数：12个变量

轴承座的位置精度直接决定轴系的对中状态，是影响齿轮啮合几何关系的关键。对于两级平行轴结构，其壳体上存在输入轴、中间轴、输出轴共三组轴承孔，每个轴承孔的位置度均需独立控制，总计12个位置度参数。这些参数的波动会直接导致中心距偏差、轴系平行度偏差及同轴度偏差，进而显著改变齿轮的传递误差（TE），如图2-3所示。

 

图2-3 副轴、差轴位置度对二级TE的影响

2.2.3 锥轴承预紧量参数：1个变量

为保证高速旋转轴系的刚性与稳定性，圆锥滚子轴承通常采用预紧设计。锥轴承的轴向预紧量作为关键装配公差，其波动会影响轴承接触角、接触刚度，从而间接改变齿轮啮合状态与振动水平，如图2-3所示，该模型中锥轴承预紧量对二级TE的影响较小，但作为轴系刚度与啮合状态的重要装配参数，其波动仍会间接改变齿轮啮合条件。为保证鲁棒性分析的全面性与工程完整性，本文仍将锥轴承预紧量纳入公差变量体系一并分析。

 

图2-3 锥轴承预紧量对NVH的影响

结合工程实践，选取对齿轮啮合状态、传递误差影响最显著的4类关键加工公差，避开次要公差以简化分析、聚焦核心方法，共计25个加公共差参数，具体如下：

轴承座位置度12个；

齿向鼓形4个；

齿向斜度4个；

渐开线鼓形4个；

锥轴承预紧量1个。

2.3 工况设定

为聚焦分析方法本身，避免多工况带来的分析冗余，本文固定单一典型工况：选取齿轮传动系统最常用的30%POT作为本次仿真工况，传递误差的波动特征能最直观地反映公差对NVH鲁棒性的影响，工况参数根据常规两级平行轴齿轮系统设计标准设定，确保分析的工程代表性。

2.4 目标设定

根据仿真验证及工程经验，通过对比不同传递误差（TE）水平对应的噪声表现，可以发现TE与齿轮系统噪声具有较强的相关性。若直接以噪声作为NVH鲁棒性分析目标，其建模与计算流程复杂、工作量极大，工程实现效率较低。因此在齿轮NVH鲁棒性分析中，通常采用传递误差TE作为替代评价指标。

3 齿轮NVH工艺鲁棒性影响分析

蒙特卡洛仿真是一种基于随机抽样的统计分析方法，通过模拟各公差参数的随机波动，重复进行大量仿真计算，统计分析输出参数（传递误差）的分布特征，进而量化公差波动对NVH鲁棒性的影响，其核心优势是能高效处理多参数随机波动问题，贴合实际加工过程中公差的随机性，且无需复杂的解析推导，工程可操作性强。本文通过蒙特卡洛方法生成加公共差范围内的随机样本，运行带有公差的系统级DOE仿真模型，选择传递误差作为评估目标，最后通过自研齿轮NVH工艺鲁棒性分析软件对数据进行分析，获取该加公共差的鲁棒性结果，具体实施流程如图3-1所示。

 

图3-1 技术路线图

3.1 齿轮传动系统仿真模型

以典型两级平行轴电驱总成为研究对象，采用Romax软件建立系统刚柔耦合动力学模型。齿轮宏观设计参数如表3-1所示，齿轮微观修形参数及圆锥轴承预紧量方案如表3-2所示，获取理想设计方案下的一级、二级TE如图3-2所示。

表3-1 齿轮宏观微观参数

 

表3-2 驱动面齿轮微观参数

 

图3-2 理想设计方案下齿轮的传递误差

3.2 齿轮NVH工艺鲁棒性仿真实施

基于搭建的仿真模型和公差分布设定，利用蒙特卡洛仿真方法开展鲁棒性分析，关键设定如下：

3.2.1 样本数量设定

样本数量的合理设定是保证蒙特卡洛仿真精度与计算效率的关键，本文采用“工程经验结合后续验证”的设定方法。结合齿轮箱NVH鲁棒性分析的工程常规经验，初步确定抽样样本数量为10000次，该数量可初步兼顾仿真精度与工程计算效率，确保能充分覆盖各公差参数的波动范围，准确反映传递误差的统计分布规律。关于样本数量对鲁棒性分析结果的具体影响，将在后续章节中专门开展对比分析，验证本章节样本数量设定的合理性与可靠性，并明确不同样本数量下的仿真精度与效率差异。

3.2.2加工公差设定

加工过程中，公差波动服从随机分布，结合机械加工工艺水平，确定各关键公差参数的分布类型、公差范围，确保符合实际生产场景，具体设定如下：

分布类型：所有公差参数均服从正态分布（机械加工中，微观修形、轴承孔位置度、锥轴承预紧量的波动多符合正态分布，均值为设计名义值，标准差根据公差等级确定）；

 

图3-3 加工共差参数设计

位置度（中间轴减速器侧X）	锥轴承预紧量
驱动一级小齿齿向鼓形	驱动二级大齿齿向斜度

图3-4 加工共差参数正态分布图

后续操作需依次设定分析工况、明确分析指标，随后启动运行计算，并对分析过程中的相关数据进行保存，为后续统计分析工作提供完整的数据支撑，如图3-5所示。

 

图3-5 二级驱动齿轮齿向鼓形参数样本分布

3.3齿轮NVH工艺鲁棒性仿真结果分析

3.3.1 鲁棒性分析指标定义

为精准量化齿轮箱NVH鲁棒性水平，结合蒙特卡洛仿真获取的传递误差（TE）数据，从波动鲁棒性与目标贴合性鲁棒性两个维度，制定明确的鲁棒性评价标准，具体如下：

（1）波动鲁棒性

波动鲁棒性用于评价设计参数受加工公差扰动时，TE的离散程度，核心反映TE在公差波动范围内的稳定性。其判定依据如下：TE的标准差处于较低水平，且TE的3σ合格率接近理论计算值，表明加工公差扰动对TE的离散影响较小，波动鲁棒性优良。

（2）目标贴合性鲁棒性

目标贴合性鲁棒性用于评价加工公差扰动下，传递误差（TE）贴近其理论设计值的能力，核心反映TE与设计目标的契合程度。其判定依据如下：

1）偏移率满足≤±10%；

 

式中：

2）TE理论设计值落在其3σ分布区间的中心附近，确保公差扰动下TE始终围绕设计值波动，贴合设计目标。

3.3.2齿轮NVH工艺鲁棒性分析软件

为高效处理蒙特卡洛仿真所得的大批量公差与传递误差数据，本文应用公司自主研发的齿轮NVH工艺鲁棒性分析软件，可实现输入参数与响应指标的快速统计分析，自动完成敏感度计算、标准差分析、3σ区间判定等鲁棒性评价工作，大幅提升了齿轮NVH鲁棒性分析的数据处理效率与结果可靠性。如图3-6所示。

图3-6 公司自研齿轮NVH工艺鲁棒性软件分析结果

3.3.3 样本数量影响分析

为确保工艺鲁棒性仿真结果的准确性与可靠性，本文分别开展5000、10000、100000次样本数仿真试验，验证样本数量对鲁棒性分析结果的影响。如图3-7和表3-3所示，结果表明：

样本量100000：一级TE均值偏移率-1.75%，标准差0.1μm，合格率99.991%，几乎100%，二级TE均值偏移率11.86%，标准差0.167μm，合格率99.317%，虽统计精度最高，但仿真耗时75小时，远高于5000次样本量的1h，时间成本过高，工程实用性较低；

样本量5000：一级TE均值偏移率-1.49%，标准差0.099μm，合格率100%，二级TE均值偏移率12.57%，标准差0.165μm，合格率99.42%，计算时间最短，但与样本量100000的二级TE结果误差较大，仅适用于前期试算；

样本量10000：一级TE均值偏移率-1.41%，标准差0.1μm，合格率100%，二级TE均值偏移率11.64%，标准差0.166μm，合格率99.34%，数据稳定、与样本量100000的结果误差小，能精准量化加工公差对NVH的影响，兼顾精度与效率，是最优选择。

综合工程精度与效率需求，本文确定10000次样本数作为后续鲁棒性分析的基准样本数量，确保仿真结果可靠且符合工程实际。

表3-3 基于样本量的NVH鲁棒性对比

 

 

 

一级TE公差鲁棒性（样本数量-5000）

 

二级TE公差鲁棒性（样本数量-5000）

 

一级TE公差鲁棒性（样本数量-100000）

 

二级TE公差鲁棒性（样本数量-100000）

 

一级级TE公差鲁棒性（样本数量-10000）

 

二级TE公差鲁棒性（样本数量-10000）

图3-7 基于样本数量的NVH鲁棒性分析对比

3.3.4 各公差参数对TE的影响程度分析

在通过10000样本量的蒙特卡洛公差仿真，获取各加工公差参数对应的齿轮TE鲁棒性数据（包括均值、偏移率、标准差及3σ合格率）后，为明确各公差参数对TE的影响程度、识别核心敏感参数，需开展针对性的敏感性分析。

传统单一变量分析法难以反映多误差耦合与实际加工工况，本文采用基于蒙特卡洛仿真的多参数TE敏感性分析，考虑了各加工公差的随机耦合作用，更贴合真实场景，结果更具工程合理性。计算各公差与TE的线性相关程度，评估其对NVH鲁棒性的影响权重，为后续公差优化及工艺可行性平衡提供数据支撑与理论依据。

采用公司自研的齿轮NVH工艺鲁棒性分析软件，获得各加工公差的相关系数如图3-8所示。

 

加工公差对一级TE敏感性分析

 

加工公差对一级TE敏感性分析

图3-8 齿轮NVH鲁棒性敏感性分析

由图3-8敏感性分析结果可见，在位置度、轴承预紧量、齿向鼓形、齿向斜度及渐开线鼓形等加工参数中，各参数对一、二级传递误差TE的影响程度差异显著。

对一级TE而言，大齿轮齿向鼓形（-0.579）、小齿轮齿向鼓形（-0.509）、大齿轮渐开线鼓形（-0.435）、小齿轮渐开线鼓形（-0.413）相关系数绝对值明显更高，为核心敏感参数，由于一级TE的3σ合格率接近100%，针对该齿轮设计参数而言，鲁棒性表现优异，因此无需对其对应公差进行调整。

对二级TE而言，敏感性排序依次为大齿轮齿向鼓形（-0.585）、小齿轮齿向鼓形（-0.532）、大齿轮渐开线鼓形（-0.307）、小齿轮渐开线鼓形（-0.295）。整体上，齿向鼓形与渐开线鼓形对TE的影响显著高于位置度、轴承预紧量及齿向斜度，且均呈负相关关系，表明适当增大修形鼓形量可有效降低TE，改善电驱齿轮NVH性能。同时，由于上述参数敏感性较高，加工波动易引起TE较大偏离（偏移率达到11.64%），因此仍需对其公差带进行收紧，以提升批量生产的鲁棒性。

将齿向鼓形公差由±3μm收紧至±2μm，渐开线鼓形公差由±3μm收紧至 ±1.5μm，并在此基础上重新开展蒙特卡洛鲁棒性仿真分析，验证公差优化效果,如图3-9所示。

 

图3-9 鼓形收严的NVH鲁棒性分析

公差收紧后，加工公差影响的二级TE均值偏移率显著下降，标准差降低，合格率提升，详细分析在后文的鲁棒性评价中介绍。

3.3.5 齿轮NVH鲁棒性评价

为精准量化齿轮箱NVH鲁棒性水平，结合工艺参数仿真获取的传递误差（TE）数据，从波动鲁棒性与目标贴合性鲁棒性两个维度，制定明确的鲁棒性评价标准。根据仿真所得数据可知，加工公差对一级、二级齿轮NVH鲁棒性的影响情况如表3-3所示。

表3-3 两级齿轮加工参数的NVH鲁棒性对比

 

（1）一级齿轮TE鲁棒性评估

一级齿轮TE作为NVH性能的核心评价指标之一，其鲁棒性表现优良。

波动鲁棒性：对于本次分析的齿轮齿轴参数模型，从加工公差扰动影响来看，该级齿轮TE的标准差为0.09927μm，离散程度较低，说明加工公差波动对其影响较小，数据分布集中，符合波动鲁棒性的合格要求。

目标贴合性：一级TE样本平均值为0.24769μm，与理论设计值的偏移率远低于10%的合格阈值，且TE理论设计值处于3σ分布区间的中心位置，贴合度较高。综合来看，一级齿轮TE在加工公差扰动下，性能稳定性强，波动风险低，鲁棒性表现优秀，无需额外调整加工公差，可满足量产一致性要求。

（2）二级齿轮TE鲁棒性评估

二级齿轮TE的鲁棒性整体达标，略逊于一级齿轮。

对比公差收紧前后的仿真结果可知：对于本次分析的齿轮齿轴参数模型，优化前二级 TE均值偏移率为11.64%，标准差0.16633μm，3σ合格率99.34%；将齿向鼓形公差由 ±3μm收紧至±2μm、渐开线鼓形公差由±3μm收紧至±1.5μm后，二级TE均值偏移率降至7.34%，标准差减小至0.12349μm，合格率提升至99.42%。

波动鲁棒性：标准差相比优化前下降了25.76个百分点，3σ合格率提高了0.08个百分点，参数波动对TE的影响得到有效抑制，批量生产的一致性与稳定性显著改善。

目标贴合性：TE均值偏移率较优化前下降了4.3个百分点，优化后更接近设计目标值，指标控制精度明显提升。

4 结论

本文提出一种基于蒙特卡洛仿真的齿轮NVH工艺鲁棒性分析方法，以两级平行轴齿轮传动系统为研究载体，针对轴承座位置度、齿轮微观修形参数等关键加工公差变量，通过量化不同公差波动对一、二级齿轮NVH性能的影响程度，明确了关键公差参数的敏感性层级，为后续公差优化设计提供了量化依据，主要结论如下：

（1）蒙特卡洛仿真可高效模拟加工公差的随机波动，精准反映加工公差与齿轮NVH 性能的关联，是齿轮NVH鲁棒性分析的有效手段，工程可操作性强，结合工程实际需求，建议选取10000次或更多样本量开展加工公差的工艺鲁棒性仿真分析；

（2）不同加工公差对齿轮TE的影响程度存在差异，对于本次分析的齿轮齿轴参数分析模型，齿向鼓形、渐开线鼓形是影响两级平行轴齿轮NVH鲁棒性的关键公差参数；

（3）对于本次分析的齿轮齿轴参数分析模型，一级TE的3σ合格率接近100%，鲁棒性表现优秀；二级TE标准差0.1663μm，偏移率约+11.64%，3σ合格率99.34%，鲁棒性整体合格但略逊于一级，公差设置参数具备进一步优化空间；

（4）根据灵敏度相关性分析结果对二级传动关键公差进行收紧优化，将齿向鼓形公差由±3μm收紧至±2μm，渐开线鼓形公差由±3μm收紧至±1.5μm。优化后二级TE均值偏移率降低4.3个百分点，标准差波动幅度降低25.76%，3σ合格率提升0.08个百分点，有效减小了二级TE的偏离程度与波动离散性，齿轮NVH鲁棒性显著改善。

通过齿轮NVH工艺参数的仿真分析，本文系统研究了关键加工公差对齿轮NVH性能的影响规律。‌结合公司自主研发的齿轮NVH工艺鲁棒性分析工具，构建了一套标准化的“多体动力学建模-公差区间定义-蒙特卡洛随机抽样-NVH统计分析”闭环分析流程。该方法可有效识别加工误差与NVH性能的映射关系，提供可复用的齿轮产品离散性控制策略，为电驱齿轮系统NVH的离散性控制提供了科学分析手段，在提升电驱系统齿轮设计效率和批量生产NVH一致性方面具有显著的工程应用价值。