CN104615878A - 一种计算带有边缘接触的齿间载荷分配方法 - Google Patents

Abstract

本发明涉及一种计算带有边缘接触的齿间载荷分配方法，首先通过不考虑边缘接触的轮齿加载接触分析(LTCA)计算，得到了每个瞬时齿轮副的齿面接触情况、目标面的接触椭圆长短轴；然后通过齿面接触分析(TCA)判断在轮齿变形条件下的齿面接触情况，同时可以判断出各个齿面上发生边缘接触的情况；通过定义接触区域实部和接触区域虚部的概念，以及定义接触区域实载荷与接触区域虚载荷的概念，可以将齿轮的有效计算齿宽(或有效齿高、亦或有效齿宽和有效齿高)进行延长，使边缘接触的区域位于延长后的虚拟齿面上，则可继续按不考虑边缘接触的轮齿加载接触分析计算，得到考虑边缘接触的齿间载荷分配结果。

Description

一种计算带有边缘接触的齿间载荷分配方法

技术领域

本发明属于齿轮几何学与力学分析交叉领域，具体涉及一种计算带有边缘接触的齿间载荷分配方法，特别是获得齿轮在拟真实工况下发生边缘接触的受力情况。

背景技术

进行齿轮副的准静力学分析，其中一个关键问题是如何求解各个瞬态的齿间载荷分配，当轮齿发生边缘接触时，这种求解计算将变得非常复杂，也难于求解。

目前，齿轮的受力计算分析主要采用ANSYS等有限元软件，通过有限元的方法可以计算齿轮在任何接触瞬时的齿面接触情况、齿轮应力应变情况、齿面接触压力等等，同时也能较好的处理齿轮副发生边缘接触的情况。但是采用ANSYS等有限元软件计算齿轮加载接触，通常非常麻烦，在齿轮有限元模型不精确时，得到的计算结果十分不可靠，要想得到较为精确的结果通常需要将有限元模型划分得十分细密，而导致的代价是计算所需要的机时十分长，同时对计算机的性能要求较高才能处理。可以说，采用有限元方法进行计算，其计算结果的好坏程度上很大一部分取决于有限元模型的建立，具有计算的不确定性。

方宗德学者在文献“齿轮轮齿承载接触分析(LTCA)的模型和方法，机械传动，1998，22(2)，p1-3”中基于柔度系数法和变形协调条件提出了一种轮齿加载接触分析方法(LTCA)，虽然在计算上相对于有限元方法具有更好的确定性，同时计算复杂度相对较小，但该方法忽略了接触柔度系数的非线性影响，以及将接触区域载荷近似处理为线接触载荷来替代实际的区域载荷，同时也没有说明对边缘接触的处理方式，因此该方法还不够完善和准确，其求解结果也不可信。

发明内容

要解决的技术问题

为了避免现有技术的不足之处，本发明提出一种计算带有边缘接触的齿间载荷分配方法。

技术方案

一种计算带有边缘接触的齿间载荷分配方法，其特征在于步骤如下：

步骤1、进行不考虑边缘接触的齿轮加载接触分析：

通过齿面接触分析方法得到：未施加载荷的每个啮合瞬时大小齿轮副的齿面接触时参与啮合齿面的接触位置、静态传动误差、接触点位置主曲率、主方向和诱导法曲率；

根据接触位置、静态传动误差、接触点位置主曲率、主方向和诱导法曲率建立各个接触位置上接触椭圆大小、齿面压缩量与载荷三者的关系式，并得到接触椭圆的方向；

根据多对相邻齿面静态传动误差的小齿轮转角区间，得到当小轮转角为给定任意角度φ₁时，所对应的小大齿轮的接触齿面对数m；

对各接触齿面对上施加的载荷分别为P₁,P₂,……P_m，根据变形协调条件和总载荷不变得到方程组：

$\left\{\begin{array}{c}\mathrm{\Δ}{\mathrm{\θ}}_1\left(P_1\right)=\mathrm{\Δ}{\mathrm{\θ}}_2\left(P_2\right)=...=\mathrm{\Δ}{\mathrm{\θ}}_i\left(P_i\right)...=\mathrm{\Δ}{\mathrm{\θ}}_m\left(P_m\right)\\ R_{G1}{P}_1\cos {\mathrm{\α}}_{t1}+R_{G2}{P}_2\cos {\mathrm{\α}}_{t2}+......+R_{\mathrm{Gm}}{P}_m\cos {\mathrm{\α}}_{\mathrm{tm}}=T\end{array}\right.---\left(1\right)$

式中，R_Gi为第i对齿轮接触点距转动中心的转动半径；α_ti表示载荷与切平面之间的夹角；Δθ_i为加载传动误差，T为小大齿轮传递总扭矩；

所述加载传动误差Δθ_i为：

Δθ_i(P_i)＝Δθ_io(P_i)+Δθ_ic(P_i)+Δθ_ib(P_i)      (2)

式中，Δθ_io(P_i)为无加载传动误差，Δθ_ic(P_i)为接触变形引起的加载传动误差，Δθ_ib(P_i)为轮齿弯曲变形引起的加载传动误差；

对三种传动误差分量进行求解：

①求解Δθ_io(P_i)

$\mathrm{\Δ}{\mathrm{\θ}}_{\mathrm{io}}\left(P_i\right)=({\mathrm{\φ}}_{i2}-{\mathrm{\φ}}_{i2}^{*})-({\mathrm{\φ}}_1-{\mathrm{\φ}}_1^{*})\frac{z_1}{z_2}---\left(3\right)$

式中，z₁,z₂为小、大齿轮的齿数；φ₁为小轮上基准齿面的转角；φ_i2为大轮上第i对面的转角；为小轮上基准齿面的初始转角；为大轮上第i对面的初始转角；

②求解Δθ_ic(P_i)

按Hertz接触理论建立起两轮齿压缩量δ_ic与载荷P_i的关系，δ_ic＝f(P_i)，则有：

$\mathrm{\Δ}{\mathrm{\θ}}_{\mathrm{ic}}\left(P_i\right)=\frac{2{\mathrm{\δ}}_{\mathrm{ic}}}{R_{\mathrm{Gi}}}---\left(4\right)$

①求解Δθ_ib(P_i)

由弯曲变形引起的传动误差分量为：

$\mathrm{\Δ}{\mathrm{\θ}}_{\mathrm{bi}}\left(P_i\right)=\frac{{\mathrm{\δ}}_{\mathrm{bi}1}+{\mathrm{\δ}}_{\mathrm{bi}2}}{R_{\mathrm{Gi}}}---\left(5\right)$

δ_bi1,δ_bi2分别代表小大齿轮一对齿轮副在接触点位置的轮齿弯曲变形量，通过有限元方法可以方便准确的求得；

联立式(1)～(5)进行求解，得到一对齿轮副某瞬时的齿间载荷分布P₁,P₂,……P_m；

根据一对齿轮副某瞬时的齿间载荷分布P₁,P₂,……P_m，通过Hertz理论公式，求得基准接触齿面上接触椭圆的大小；

将计算得到的齿间载荷分布P₁,P₂,……P_m施加到各接触齿面，计算弯曲变形和接触变形，通过齿面接触分析方法得到施加载荷的每个啮合瞬时大小齿轮副的齿面接触时的传动误差；

根据多对相邻齿面传动误差的小齿轮转角区间，得到当小轮转角为给定任意角度φ₁时，所对应的小大齿轮的接触齿面对数k；若k>m，则将k赋给m，重新本步骤进行运算，若k＝m，则进行下一环节的计算；

步骤2、建立带有边缘接触的轮齿加载接触分析计算模型：

提取一个基准齿面的接触椭圆大小及方向，重新以每个接触点为椭圆中心，再结合所求的基准齿面椭圆；

当接触齿面出现边缘接触时，定义接触椭圆位于有效齿宽或高内的区域为接触区域实部，接触椭圆位于有效齿宽或高外的区域为接触区域虚部，接触区域实部Dr内的压力集合为接触区域实载荷，接触区域虚部Dv内的压力集合为接触区域虚载荷，得到接触区域实载荷F_r为：

$F_r=\underset{D_r}{\∫\∫}{p}_0{\{1-{(x/a)}^2-{(y/b)}^2\}}^{1/2}\mathrm{dxdy}---\left(6\right)$

a，b为接触椭圆的长短轴半径；p₀为接触椭圆中心压力；

得到接触区域虚载荷F_v如下所示：

$F_v=\underset{D_v}{\∫\∫}{p}_0{\{1-{(x/a)}^2-{(y/b)}^2\}}^{1/2}\mathrm{dxdy}---\left(7\right)$

将有效计算齿宽进行延长，使边缘接触的整个椭圆区域全部位于虚拟齿宽内；

当引入虚拟齿宽的边缘接触情况时，在计算接触变形量引起的加载传动误差时，以P_i'＝xsm_iF_ri取代发生边缘接触的齿间载荷分布P₁,P₂,……P_m；其中：xsm为载荷占有比，F_ri为齿间载荷在边缘接触齿上的边缘接触处分量P_i，i表示第i个接触齿面对；

步骤3、齿间载荷分配计算：按照未考虑边缘接触的齿轮加载接触分析计算方法计算上述修正后的带有边缘接触的计算模型，得到带有边缘接触的齿轮副齿间载荷分配结果。

有益效果

本发明提出的一种计算带有边缘接触的齿间载荷分配方法，首先通过不考虑边缘接触的轮齿加载接触分析(LTCA)计算，得到了每个瞬时齿轮副的齿面接触情况、目标面的接触椭圆长短轴；然后通过齿面接触分析(TCA)判断在轮齿变形条件下的齿面接触情况，同时可以判断出各个齿面上发生边缘接触的情况；通过定义接触区域实部和接触区域虚部的概念，以及定义接触区域实载荷与接触区域虚载荷的概念，可以将齿轮的有效计算齿宽(或有效齿高、亦或有效齿宽和有效齿高)进行延长，使边缘接触的区域位于延长后的虚拟齿面上，则可继续按不考虑边缘接触的轮齿加载接触分析计算，得到考虑边缘接触的齿间载荷分配结果。

本发明的有益效果是：结合带有轮齿变形的齿面接触分析计算(TCA)，较好的解决了对带有边缘接触的齿面接触情况判断；然后，通过定义接触区域实部和接触区域虚部的概念，以及定义接触区域实载荷与接触区域虚载荷的概念，建立了带有边缘接触的轮齿加载接触分析计算模型，该方法使发生边缘接触的计算和没有发生边缘接触的计算趋于一致，较好的解决了齿轮在发生边缘接触时所带来的计算复杂性问题。本发明为实现齿轮的准静力学分析打下了良好的基础，使现有基于ANSYS等有限元软件的齿轮受力分析变得更加简便和具体，相对于现有的计算手段，该方法的计算机时将大大缩减。

附图说明

图1(a)接触椭圆中心在有效齿宽(高)外的边缘接触；

图1(b)接触椭圆中心在有效齿宽(高)内的边缘接触；

图2接触椭圆区域局部压力的计算；

图3(a)接触椭圆中心在有效齿宽(高)外的虚拟齿宽(高)延长；

图3(b)接触椭圆中心在有效齿宽(高)外的虚拟齿宽(高)延长。

具体实施方式

现结合实施例、附图对本发明作进一步描述：

(1)进行不考虑边缘接触的齿轮加载接触分析

通过齿面接触分析(TCA)可以得到，在未施加载荷的每个啮合瞬时大小齿轮副的齿面接触情况，即此时参与啮合齿面的接触位置、静态传动误差、接触点位置主曲率、主方向、诱导法曲率，从而可以建立起各个接触位置上接触椭圆大小、齿面压缩量与载荷的关系，同时，也可以得到接触椭圆的方向。

根据上述轮齿无载荷接触分析(TCA)得到的多对齿面的传动误差，可以判断小齿轮在转动某转角φ₁时，小大齿轮的接触齿面对数m。

假设各接触齿面对上分担的载荷分别为P₁,P₂,……P_m，则根据变形协调条件和总载荷不变得到方程组：

$\left\{\begin{array}{c}\mathrm{\Δ}{\mathrm{\θ}}_1\left(P_1\right)=\mathrm{\Δ}{\mathrm{\θ}}_2\left(P_2\right)=...=\mathrm{\Δ}{\mathrm{\θ}}_i\left(P_i\right)...=\mathrm{\Δ}{\mathrm{\θ}}_m\left(P_m\right)\\ R_{G1}{P}_1\cos {\mathrm{\α}}_{t1}+R_{G2}{P}_2\cos {\mathrm{\α}}_{t2}+......+R_{\mathrm{Gm}}{P}_m\cos {\mathrm{\α}}_{\mathrm{tm}}=T\end{array}\right.---\left(1\right)$

式中，R_Gi为第i对齿轮接触点距转动中心的转动半径；α_ti表示载荷与切平面之间的夹角；Δθ_i为加载传动误差，T为小大齿轮传递总扭矩。

已知，加载传动误差Δθ_i由三部分组成，表示为:

Δθ_i(P_i)＝Δθ_io(P_i)+Δθ_ic(P_i)+Δθ_ib(P_i)    (2)

式中，Δθ_io(P_i)为无加载传动误差，Δθ_ic(P_i)为接触变形引起的加载传动误差，Δθ_ib(P_i)为轮齿弯曲变形引起的加载传动误差。下面对这三种传动误差分量进行求解。

①求解Δθ_io(P_i)

根据上述无载荷轮齿接触分析，可以得到:

$\mathrm{\Δ}{\mathrm{\θ}}_{\mathrm{io}}\left(P_i\right)=({\mathrm{\φ}}_{i2}-{\mathrm{\φ}}_{i2}^{*})-({\mathrm{\φ}}_1-{\mathrm{\φ}}_1^{*})\frac{z_1}{z_2}---\left(3\right)$

式中，z₁,z₂为小、大齿轮的齿数；φ₁为小轮上基准齿面的转角；φ_i2为大轮上第i对面的转角；为小轮上基准齿面的初始转角；为大轮上第i对面的初始转角。

②求解Δθ_ic(P_i)

按Hertz接触理论建立起两轮齿压缩量δ_ic与载荷P_i的关系，即δ_ic＝f(P_i)，则有：

$\mathrm{\Δ}{\mathrm{\θ}}_{\mathrm{ic}}\left(P_i\right)=\frac{2{\mathrm{\δ}}_{\mathrm{ic}}}{R_{\mathrm{Gi}}}---\left(4\right)$

①求解Δθ_ib(P_i)

由弯曲变形引起的传动误差分量为：

$\mathrm{\Δ}{\mathrm{\θ}}_{\mathrm{bi}}\left(P_i\right)=\frac{{\mathrm{\δ}}_{\mathrm{bi}1}+{\mathrm{\δ}}_{\mathrm{bi}2}}{R_{\mathrm{Gi}}}---\left(5\right)$

δ_bi1,δ_bi2分别代表小大齿轮一对齿轮副在接触点位置的轮齿弯曲变形量，通过有限元方法可以方便准确的求得。

联立式(1)～(5)进行求解，则可以得到一对齿轮副某瞬时的齿间载荷分布P₁,P₂,……P_m。

通过Hertz理论公式，求得基准接触齿面上接触椭圆的大小，由此进一步判断加载后两齿轮的接触情况。

重新计算计入弯曲变形和接触变形后的齿轮接触情况，判断小齿轮在转动某转角φ1时，小大齿轮的接触点对数k，若k>m，则将k赋给m，重新进行整节的运算，若k＝m，则进行下一环节的计算。

(2)建立带有边缘接触的轮齿加载接触分析计算模型

通过上一节不考虑边缘接触的轮齿加载接触分析计算，提取一个基准齿面的接触椭圆大小及方向，下面以重新以每个接触点为椭圆中心，再结合所求的基准齿面椭圆大小，进行区域分析计算判断，观察哪些接触齿面出现边缘接触。

如图1(a)所示，为接触椭圆中心在有效齿宽(高)外的边缘接触，称为第一类边缘接触，此时，小大齿轮的接触点对数m＝m+1；如图1(b)所示，为接触椭圆中心在有效齿宽(高)内的边缘接触，称为第二类边缘接触，此时，小大齿轮的接触点对数m＝m。

为了求得边缘接触时的加载传动误差，现提出接触区域实部和接触区域虚部的概念。如图1(a)和图1(b)所示，接触区域实部即发生边缘接触时计算的接触椭圆位于有效齿宽(高)内的区域；接触区域虚部即发生边缘接触时计算的接触椭圆位于有效齿宽(高)外的区域。

同时，提出接触区域实载荷与接触区域虚载荷的概念，接触区域实载荷即接触区域实部Dr内的压力集合；接触区域虚载荷即接触区域虚部Dv内的压力集合。如图2所示，可以得到接触区域实载荷F_r如下表示：

$F_r=\underset{D_r}{\∫\∫}{p}_0{\{1-{(x/a)}^2-{(y/b)}^2\}}^{1/2}\mathrm{dxdy}---\left(6\right)$

a，b为接触椭圆的长短轴半径；p₀为接触椭圆中心压力。

同理，如图2所示，接触区域虚载荷F_v如下所示：

$F_v=\underset{D_v}{\∫\∫}{p}_0{\{1-{(x/a)}^2-{(y/b)}^2\}}^{1/2}\mathrm{dxdy}---\left(7\right)$

下面，在引入接触区域实部D_r、接触区域虚部D_v、接触区域实载荷F_r、接触区域虚载荷F_v等概念后，如图3(a)和图3(b)所示，将有效计算齿宽(有些需要沿齿高方向，或者同时沿齿宽和齿高方向，下面以齿宽为例进行说明)进行延长，使边缘接触的整个椭圆区域全部位于虚拟齿宽类，则可以按不考虑边缘接触的情况进行计算，只是计算时略有差别。

对于引入虚拟齿宽的边缘接触情况，在计算弯曲变形量和接触变形量引起的加载传动误差时，需要如下处理。

已知接触椭圆中心接触变形量是由整个椭圆部分的压力引起，而椭圆区域中心的弯曲变形量则仅由接触区域实载荷F_ri产生，接触区域实载荷F_ri即为齿间载荷在边缘接触齿上的边缘接触处分量P_i，i表示第i个接触齿面对。

下面，引入一个系数——载荷占有比xsm，其表示为：

$\mathrm{xsm}=\frac{F_r+F_v}{F_r}---\left(8\right)$

则边缘接触椭圆的混合载荷P_i’可以表达为

P_i'＝xsm_iF_ri          (9)

下面只需对上述计算两轮齿接触压缩量δ_ic的载荷P_i进行修正即可，即用混合载荷P_i’代替原有的P_i进行接触变形量计算即可。

至此，带有边缘接触的轮齿加载接触分析计算模型已经建立起来。

(3)齿间载荷分配计算

按照未考虑边缘接触的齿轮加载接触分析计算方法计算上述修正后的带有边缘接触的计算模型，即可得到带有边缘接触的齿轮副齿间载荷分配结果。

本发明所述的计算带有边缘接触的齿间载荷分配方法，很好的解决了轮齿加载接触分析技术中对边缘接触处理的难题。通过有限元方法的验证，已经初步证实了该方法的正确性。该方法的实现，意味着可以比较完整的实现轮齿加载接触分析技术，从而奠定了齿轮准静力学分析的基础。

Claims (1)

1.一种计算带有边缘接触的齿间载荷分配方法，其特征在于步骤如下：

步骤1、进行不考虑边缘接触的齿轮加载接触分析：

通过齿面接触分析方法得到：未施加载荷的每个啮合瞬时大小齿轮副的齿面接触时参与啮合齿面的接触位置、静态传动误差、接触点位置主曲率、主方向和诱导法曲率；

根据接触位置、静态传动误差、接触点位置主曲率、主方向和诱导法曲率建立各个接触位置上接触椭圆大小、齿面压缩量与载荷三者的关系式，并得到接触椭圆的方向；

根据多对相邻齿面静态传动误差的小齿轮转角区间，得到当小轮转角为给定任意角度φ₁时，所对应的小大齿轮的接触齿面对数m；

对各接触齿面对上施加的载荷分别为P₁,P₂,……P_m，根据变形协调条件和总载荷不变得到方程组：

$\left\{\begin{array}{c}{\mathrm{\Δ\θ}}_1\left(P_1\right)=\mathrm{\Δ}{\mathrm{\θ}}_2\left(P_2\right)=\·\·\·={\mathrm{\Δ\θ}}_i\left(P_i\right)\·\·\·=\mathrm{\Δ}{\mathrm{\θ}}_m\left(P_m\right)\\ R_{G1}{P}_1\cos {\mathrm{\α}}_{t1}+R_{G2}{P}_2\cos {\mathrm{\α}}_{t2}+\·\·\·\·\·\·+R_{\mathrm{Gm}}{P}_m\cos {\mathrm{\α}}_{\mathrm{tm}}=T\end{array}\right.---\left(1\right)$

式中，R_Gi为第i对齿轮接触点距转动中心的转动半径；α_ti表示载荷与切平面之间的夹角；Δθ_i为加载传动误差，T为小大齿轮传递总扭矩；

所述加载传动误差Δθ_i为：

Δθ_i(P_i)＝Δθ_io(P_i)+Δθ_ic(P_i)+Δθ_ib(P_i)          (2)

式中，Δθ_io(P_i)为无加载传动误差，Δθ_ic(P_i)为接触变形引起的加载传动误差，Δθ_ib(P_i)为轮齿弯曲变形引起的加载传动误差；

对三种传动误差分量进行求解：

①求解Δθ_io(P_i)

${\mathrm{\Δ\θ}}_{\mathrm{io}}\left(P_i\right)=({\mathrm{\φ}}_{i2}-{\mathrm{\φ}}_{i2}^{*})-({\mathrm{\φ}}_1-{\mathrm{\φ}}_1^{*})\frac{z_1}{z_2}---\left(3\right)$

式中，z₁,z₂为小、大齿轮的齿数；φ₁为小轮上基准齿面的转角；φ_i2为大轮上第i对面的转角；为小轮上基准齿面的初始转角；为大轮上第i对面的初始转角；

②求解Δθ_ic(P_i)

按Hertz接触理论建立起两轮齿压缩量δ_ic与载荷P_i的关系，δ_ic＝f(P_i)，则有：

${\mathrm{\Δ\θ}}_{\mathrm{ic}}\left(P_i\right)=\frac{{2\mathrm{\δ}}_{\mathrm{ic}}}{R_{\mathrm{Gi}}}---\left(4\right)$

①求解Δθ_ib(P_i)

由弯曲变形引起的传动误差分量为：

${\mathrm{\Δ\θ}}_{\mathrm{bi}}\left(P_i\right)=\frac{{\mathrm{\δ}}_{\mathrm{bi}1}+{\mathrm{\δ}}_{\mathrm{bi}2}}{R_{\mathrm{Gi}}}---\left(5\right)$

δ_bi1,δ_bi2分别代表小大齿轮一对齿轮副在接触点位置的轮齿弯曲变形量，通过有限元方法可以方便准确的求得；

联立式(1)～(5)进行求解，得到一对齿轮副某瞬时的齿间载荷分布P₁,P₂,……P_m；

根据一对齿轮副某瞬时的齿间载荷分布P₁,P₂,……P_m，通过Hertz理论公式，求得基准接触齿面上接触椭圆的大小；

将计算得到的齿间载荷分布P₁,P₂,……P_m施加到各接触齿面，计算弯曲变形和接触变形，通过齿面接触分析方法得到施加载荷的每个啮合瞬时大小齿轮副的齿面接触时的传动误差；

根据多对相邻齿面传动误差的小齿轮转角区间，得到当小轮转角为给定任意角度φ₁时，所对应的小大齿轮的接触齿面对数k；若k>m，则将k赋给m，重新本步骤进行运算，若k＝m，则进行下一环节的计算；

步骤2、建立带有边缘接触的轮齿加载接触分析计算模型：

提取一个基准齿面的接触椭圆大小及方向，重新以每个接触点为椭圆中心，再结合所求的基准齿面椭圆；

当接触齿面出现边缘接触时，定义接触椭圆位于有效齿宽或高内的区域为接触区域实部，接触椭圆位于有效齿宽或高外的区域为接触区域虚部，接触区域实部Dr内的压力集合为接触区域实载荷，接触区域虚部Dv内的压力集合为接触区域虚载荷，得到接触区域实载荷F_r为：

$F_v=\underset{D_v}{\∫\∫}{p}_0{\{1-{(x/a)}^2-{(y/b)}^2\}}^{1/2}\mathrm{dxdy}---\left(6\right)$

a，b为接触椭圆的长短轴半径；p₀为接触椭圆中心压力；

得到接触区域虚载荷F_v如下所示：

$F_v=\underset{D_v}{\∫\∫}{p}_0{\{1-{(x/a)}^2-{(y/b)}^2\}}^{1/2}\mathrm{dxdy}---\left(7\right)$

将有效计算齿宽进行延长，使边缘接触的整个椭圆区域全部位于虚拟齿宽内；

当引入虚拟齿宽的边缘接触情况时，在计算接触变形量引起的加载传动误差时，以P_i'＝xsm_iF_ri取代发生边缘接触的齿间载荷分布P₁,P₂,……P_m；其中：xsm为载荷占有比，F_ri为齿间载荷在边缘接触齿上的边缘接触处分量P_i，i表示第i个接触齿面对；

步骤3、齿间载荷分配计算：按照未考虑边缘接触的齿轮加载接触分析计算方法计算上述修正后的带有边缘接触的计算模型，得到带有边缘接触的齿轮副齿间载荷分配结果。