吕程　易广威　刘子建

摘要：以圆柱形孔轴配合的结合面为研究对象，研究了新一代GPS标准体系下多公差项的孔轴配合实际误差模型。以圆柱度为例，结合GPS标准中圆柱面形状误差的评定方法，采用蒙特卡洛法模拟误差的随机性，建立了孔轴圆柱面体外拟合尺寸实际变动区间求解模型，并分析实际误差对配合性质的影响。分析了包括方向公差（或位置公差）在内的三类公差耦合作用下孔轴结合面误差的形成机理，建立了不同配合性质下的圆柱形孔轴结合面的实际误差模型，获得了结合面误差分量的实际变动区间以及实际的配合性质。实现了在设计阶段对孔轴实际装配精度的预测。以实际的孔轴结合面精度与配合性质要求为约束，孔轴加工成本为目标，进行了孔轴零件的公差优化设计。最后以一孔轴配合的装配误差分析与公差优化为例，验证了该方法的可行性与实用性。

关键词：蒙特卡洛模拟法；结合面误差；装配精度预测；配合性质；公差优化设计

中图分类号：TH115 文献标识码：A

孔轴配合是机械产品中最常见的重要配合类型之一。根据产品的功能要求确定孔轴的配合性质，进行合理的公差设计，对保证装配精度和降低制造成本具有重要的意义。孔轴的装配精度影响结合面的接触状态，从而对产品的运动精度、装配难易程度和使用寿命等产生影响。因此，进行孔轴装配精度研究具有重要的工程意义。

国内外学者对圆柱形孔轴配合精度的影响因素与建模方法进行了大量有意义的研究。徐旭松和黄芳研究了基于新一代GPS规范体系的SDT公差建模理论，探讨了装配公差分析模型。U.Roy等分析了公差的语义，建立了基于数学定义的尺寸公差、形状公差、定向与定位公差的数学模型，用于求解孔轴装配误差。Anselmetti将偏差分解为基本几何元MEDG的变动，通过建立约束几何模型和偏差矢量，计算孔轴偏差传递。周思杭通过对单个零件及零件间的偏差传递分析，建立偏差传递模型，提出了基于偏差传递模型的装配精度计算方法。刘伟东等详细阐述了装配偏差的种类、作用，偏差计算模型和评价方法，提出了孔轴配合偏差的有向图表达方法。

上述研究重点针对孔轴间的某类公差进行，从装配结合面角度对存在多种公差耦合作用的孔轴配合装配误差建模方法以及加工误差对配合性质影响的研究还比较少见。零件加工误差通过装配结合面传递、累积，因此，结合面是误差传递累积作用的关键点。本文针对圆柱形孔轴配合的装配结合面，在考虑多种误差数据随机性的基础上，提出孔轴结合面装配误差的建模方法，分析了实际误差对配合性质的影响。据此在设计阶段实现对孔轴装配精度的预测，并以此为依据进行零件的公差优化设计，保证设计要求的装配精度与功能要求的配合性质。

1 孔轴圆柱面的误差与配合性质

孔轴圆柱面的加工误差是孔轴装配结合面误差的主要来源。需首先明确圆柱表面自身的误差变动，因此，此时圆柱面为单一要素，即除尺寸公差外，只有形状公差要求。两者的关系可以为相互独立、相互影响或相互补偿，并以具体关系为依据判别加工零件是否合格。独立原则是处理两者关系的基本原则。下面主要讨论独立原则下圆柱面在尺寸公差与形状公差要求下的误差变动情况。

1.1 孔轴圆柱面误差的实际变动分析

以圆柱度作为形状公差，分析同时存在尺寸公差和形状公差要求时孔轴圆柱面的误差变动。由于形状误差的存在使圆柱面的形状发生了改变，分析这种改变对孔轴配合性质的影响，通常以孔和轴的体外拟合尺寸代替实际孔轴圆柱面尺寸进行配合误差分析。结合GPS标准体系中圆柱度误差的测量与评定方法，可知圆柱度公差作用下的孔轴圆柱面与其体外拟合尺寸的关系如图1所示。

图1中，圆柱的轴线为实际圆柱面的拟合轴线，T_F为圆柱面形状公差，D_ae为该圆柱面作为内表面（孔）时的体外拟合尺寸，d_ae为该圆柱面作为外表面（轴）时的体外拟合尺寸，d_act为任意一处圆柱面的局部实际尺寸。

设圆柱面的直径基本尺寸为D，尺寸公差为T_D，上下偏差分别为T_DU，T_DL（T_D=T_DU-T_DL）。可知局部实际尺寸的变动区间为：

上述讨论确定了在尺寸公差与形状公差作用下，圆柱面作为孔和轴表面的体外拟合尺寸的变动区间及约束关系。但未考虑加工过程的随机性，无法确定孔轴圆柱面体外拟合尺寸的实际变动区间。

1.2 孔轴体外拟合尺寸的实际变动区间

由于圆柱面加工误差具有随机性，因此孔轴体外拟合尺寸的变动也具有随机性。本文采用蒙特卡洛法模拟孔轴体外拟合尺寸变动的随机性，进而获取孔轴体外拟合尺寸的实际变动区间。蒙特卡洛模拟法是通过对随机变量的统计实验和随机模拟来求解问题近似解的方法，具有适应性强、方便程序实现、求解误差与问题维数无关等优点。

蒙特卡洛法求解孔轴体外拟合尺寸实际变动区间的步骤如下：

1）假设孔轴体外拟合尺寸符合正态分布（有效分布范围6σ）。由式（2）得D_ae，d_ae理想分布的均值及标准差为（D+（T_DL+T_DU-T_F）/2，（T_D+T_F）/6），（D+（T_DL+T_DU+T_F）/2，（T_D+T_F）/6）。

2）进行孔轴体外拟合尺寸的模拟抽样。因约束关系式（3）的存在，d_ae，D_ae不同的变动顺序可能导致取值区间不同。因此，分别对两种变动顺序进行抽样。以d_ae→D_ae为例，抽样流程如图2所示：

图2中，k记录满足条件的样本数量，初值为0；K为每种抽样顺序要求的合格样本数，取K=10000；K1为常量，用以防止抽样过程陷入死循环。完成抽样后可得2K组孔轴体外拟合尺寸样本。

3）因d_ae，D_ae的实际分布与理想分布存在差异，因此采用x²拟合检验法对孔轴体外拟合尺寸进行假设验证。若经验证属于正态分布，则采用极大似然估计法估计分布函数的均值与方差分别为：若经验证不属于正态分布，则可通过当量正态化法将其转化为等效正态分布，在此不做详细讨论。

4）根据孔轴体外拟合尺寸的实际分布类型与标准差，查表可得孔轴体外拟合尺寸的实际变动区间带宽为：式中：G为相对分布系数，正态分布时G取1。孔轴体外拟合尺寸的实际变动区间为：

1.3 孔轴配合的实际配合性质

孔轴的配合性质与孔和轴圆柱面的尺寸和形状误差相关，直接影响运动精度、装配性能等，是保证产品功能的重要因素。由于孔轴实际表面在随机误差作用下偏离了理想的设计尺寸与形状，可能导致孔轴的配合性质超出设计范围。因此，有必要对孔轴的实际配合性质进行分析。

孔轴的配合性质与孔和轴的体外拟合尺寸直接相关。根据孔与轴各自的公差设计值，通过蒙特卡洛法分别获得孔轴体外拟合尺寸样本。以孔轴配合间隙S_c衡量孔轴配合性质。S_c≥0表示最小间隙为0的间隙配合，S_c<0表示过盈配合。可通过对样本的分析获取实际配合性质与设计要求的配合性质之间的关系。例如，通过对间隙值的分析可得实际配合性质符合设计要求的概率。

2 孔轴结合面误差建模

圆柱形孔轴结合面是指具有配合关系的圆柱形孔轴表面相互贴合形成的一对接触面。孔轴的尺寸、形状及方向位置误差伴随装配过程在结合面处耦合累积形成的结合面误差，是影响装配体误差传递、运动精度、产品性能与使用寿命等的关键因素。在产品设计阶段依据精度设计值通过合理的误差模型进行产品装配精度预测是验证精度设计合理性的重要手段。

2.1 结合面误差的形成机理

孔与轴的理想轴线的相对位姿变动体现了孔轴各类误差在结合面处综合作用的效果。下面以孔轴理想轴线的相对位姿变动讨论结合面误差，如图3所示。

图3中，孔与轴的表面为实际体外拟合圆柱面。T_P1，T_P2分别为轴与孔的方向或位置公差。a₁，a₂为轴、孔的理想轴线，a₁′，a₂′为轴、孔的实际体外拟合圆柱面轴线，l为孔轴的配合长度。

理想装配（零误差）时应有a₁与a₂重合，因方向或位置误差的存在使实际拟合轴线a₁′与a₂′偏离各自的理想轴线，又因配合间隙等使a₁′与a₂′不重合，形成了结合面复杂的误差形式。以轴为装配基准件，结合面的误差形成机理如图4所示。

2.2 结合面误差建模

为方便描述结合面的误差变动，采用小位移旋量（SmallDisplacementTorsor，SDT）描述结合面误差，即S=（α，β，δ，u，v，ω）。其中口，β，δ表示绕x，y，z轴转动的微小变动量，u，v，ω表示沿x，y，z轴移动的微小变动量。

根据新一代GPS标准中的恒定度概念可知，当一个几何要素受到自身形状与方位限制后，仍可沿（或绕）某一坐标轴方向平动（或转动），则其具有该方向上的恒定度，对应的SDT分量为0。以垂直度及位置度为例，设孔轴的理想轴线均为z轴方向（后文均默认轴线为z轴方向），其SDT如表1所示。

设图3中公差T_P1，T_P2均为具有圆柱形公差域的位置公差。因圆柱面具有径向同性的特点，因此以过轴线的xoz截面内的孔轴位姿为例，分析结合面的误差变动。结合面误差为理想轴线a₂相对于a₁的位姿变动（误差分量为β，ω）。结合面误差根据其形成机理，各部分误差分量的变动可统一表示如下：式中：β_iU，β_iL为第i部分误差中误差分量β变动区间的上、下边界；ω_iU，ω_iL为第i部分误差中误差分量叫变动区间的上、下边界，-l/2≤x≤l/2（后文中x取值均为该区间）。则结合面误差为：

各部分误差中，各误差分量变动区间的边界取值如下：

1）i=11时，位置公差T_P1限制了a₁到a₁′的误差变动，因此β_iU=-β_iL=T_P1/l，ω_iU=-ω_iLT_P1/2。

2）i=12时，a₁'到a₂'的误差变动与孔轴具体的配合性质直接相关，具体如下：

a）过盈配合时，因孔轴圆柱表面完全贴合，因此轴线重合，无相对位姿变动。β_i=0，ω_i=0。

b）间隙配合时，因配合间隙S_c=D₂-D₁，使孔轴产生相对位姿变动。因此有β_iU=-β_iL=S_c/l，ω_iU=-ω_iL-S_c/2。D₁，D₂为轴、孔实际体外拟合尺寸。

3）i=22时，位置公差T_P2限制了a₂'到a₂的误差变动，因此β_iU=-β_iL=T_P2/l，ω_iU=-ω_iLT_P2/2。

结合式（5）与步骤1）～3）的误差分量变动区间，通过蒙特卡洛法分别获得各β_iU与各叫ω_iU的变动区间，由式（6）可得结合面误差分量的实际变动区间。

圆柱面直径尺寸公差域和形状公差域均径向同性，若孔与轴的位置公差域也径向同性（如公差域为圆柱区域），则结合面误差的SDT为（0，β，δ，0，v，ω），且δ₁₂=β₁₂，v₁₂=ω₁₂。否则δ₁₂与v₁₂需根据各部分误差源对应的公差域特点具体分析，在此不再赘述。根据设计给定的公差值及上述结合面误差模型可实现具体孔轴配合的装配精度预测。

3 孔轴配合的公差优化设计

加工成本是进行公差设计需要考虑的主要因素。如果经误差分析与装配精度预测后，孔轴装配精度或配合性质不满足设计要求，则需改进公差设计来减小误差，提高配合质量。进行公差优化设计的目的是在满足精度与配合要求的前提下，同时满足加工的经济性。

遗传算法是源于自然遗传进化机制的搜索优化算法，已经广泛应用于公差优化研究中。本文以加工成本为目标，结合面装配精度与配合间隙为约束条件，且满足公差设计时T_D>T_P>T_S的基本原则，其中T_D为尺寸公差，T_P为位置公差，T_S为形状公差。设孔与轴的各项公差的公差域均有径向同性的特点，即结合面误差的SDT为（0，β，δ，0，V，ω），且δ=β，V=ω。按照带约束优化问题的表达形式，孔轴配合的公差优化问题可表达如下：式中：C_F（T）为加工成本函数；n为公差项数；m为约束不等式g_i（T）的数量。

上述问题的约束不等式多为复杂的隐式函数，因此优化过程也非常复杂。为兼顾工程经济性与计算效率，本文将以一合理的计算代数作为优化过程的结束条件，因此需保证优化过程中的每组迭代点均处于约束的可行域内。而内点惩罚函数法是在处理不等式约束优化问题中满足上述需求的一种十分有效的方法。因此，本文选取内点惩罚函数法将上述问题转化为非约束问题。根据约束g_i（T）≤0，可取惩罚项为：式中：c为下降系数，c<1，一般可取c=0.5～0.7。初始惩罚因子r^（0）的选取应适当，r^（0）选取过小时，可能导致函数形态变坏，r^（0）选取过大将使收敛速度下降，计算效率降低，但可保证计算过程的稳定性。通常，初始点的选择应使惩罚项在适应度函数中不起到主导作用。遗传算法公差优化流程如图5所示。

图5中G为优化过程的计算代数。Φ_min为计算过程中的优化结果，其初始值由当前公差设计值根据式（7）计算。

4 实例分析

图6为孔轴配合零件的公差设计图。要求孔轴结合面误差分量β₁₂，δ₁₂位于区间[-0.005，0.005]内，v₁₂，ω₁₂位于区间[-0.1200，0.1200]内，且满足设计图给定配合间隙的概率不小于98%。

由图5可知，设计的配合间隙为[0，0.25]。设计给出的轴的尺寸公差、位置公差、形状公差分别为T_D1=0.1（T_OU1=0，T_DL1=0.1），T_P1=0.08，T_P1=0.04，孔的尺寸公差、位置公差、形状公差分别为T_D2=0.15（T_DU2=0.15，T_DL2=0），T_P2=0.1，T_P2=0.04。孔与轴的各项公差的公差域均径向同性。根据前文所述方法，通过MATLAB编程实现结合面误差求解、配合性质分析及公差优化设计。过程如下：

1）通过蒙特卡洛模拟分别求得孔轴实际体外作用尺寸的实际变动区间，通过对模拟结果的样本分析获得孔轴实际配合间隙的分布，如图7所示。

图7中，满足配合要求的概率为98.9%，满足设计要求，即虽然产生了超出配合要求的过盈配合，但概率很小，可忽略其影响。

2）由于孔与轴的位置公差域均径向同性，可知结合面误差的SDT为（0，β，δ，1，v，ω），且δ₁₂=β₁₂，v₁₂=ω₁₂，因此，可仅针对β₁₂与叫ω₁₂进行求解。结合面误差分量的实际变动区间为：

另外，为验证结合面误差求解方法的精确效果，通过极值法求得结合面误差分量的变动区间如下：对比两种方法所求的结合面误差，如表2所示。

由表2可知，本文所述的方法较极值法求得的结合面误差的变动区间宽度平均下降了18.18%。但实际的精度仍然未达到精度设计要求。在此，进行公差的优化设计以提高精度，满足装配精度要求。

3）孔轴表面的公差优化设计。采用文献给出的加工成本模型，也是加工成本最常用的模型：

a）轴的尺寸与形状公差成本函数：b）孔的尺寸及形状公差成本函数：

代人计算得总公差加工成本为15.21。

根据图5所示流程，通过MATLAB编程实现公差的优化设计。遗传算法的基本参数为：种群规模M-100，交叉概率R-0.8，变异概率Pm-0.05，计算代数G-200。优化结果（保留两位小数）如表3所示。

根据公差优化值可得优化后的配合间隙满足配合要求的概率为98.10%，零件公差加工成本为11.72。可见，通过公差的优化设计，在满足精度要求的前提下大幅降低了加工成本。

本文所述的结合面误差求解方法考虑了多种公差的耦合作用，可更准确地预测装配精度和实际配合性质，并提供了公差优化设计方法，从而满足孔轴配合的精度要求，并降低了加工成本。

5 结论

本文采用蒙特卡洛模拟法模拟圆柱面实际加工误差的随机变动，求得了孔与轴的体外拟合尺寸的实际变动区间，分析了实际的配合性质。以孔轴配合的结合面误差为研究对象，研究了包括尺寸公差、形状公差、位置（或方向）公差在内的多种公差耦合作用的圆柱形孔轴结合面误差的形成机理，建立了孔轴结合面的实际误差模型，为在设计阶段较准确地预测孔轴实际配合性质与装配精度提供了有效方法。同时给出了基于遗传算法的孔轴各项公差的优化设计方法。论文通过圆柱形孔轴结合面误差分析与公差优化设计实例，验证了本文研究方法的可行性、实用性与经济性。