考虑系统变形的驱动桥准双曲面齿轮啮合效率计算方法
王钦1, 贺迪1, 桂良进1, 胡智宇2, 彭金3, 范子杰1    
1. 清华大学 车辆与运载学院, 汽车安全与节能国家重点实验室, 北京 100084;
2. 陕西汉德车桥有限公司, 西安 710201;
3. 汉德车桥(株洲)齿轮有限公司, 株洲 412000
摘要:准双曲面齿轮是驱动桥系统的关键零部件之一, 常在高速重载工况下工作。由于系统变形大且齿面啮合状态复杂, 重载工况下会引起齿面载荷的显著偏载。因此, 考虑系统变形对准双曲面齿轮啮合性能的影响具有重要意义。该文针对复杂工况下的驱动桥准双曲面齿轮, 提出了考虑系统变形的齿轮啮合效率计算方法。首先, 将驱动桥系统简化为多支撑轴系耦合模型, 计算不同载荷工况下系统变形引起的齿轮啮合错位量。然后, 基于齿轮机床加工参数计算实际齿面, 并采用齿轮摩擦加载接触分析(frictional loaded tooth contact analysis, FLTCA)方法求解得到实际运行工况下齿面载荷分布与啮合效率。最后, 通过驱动桥台架试验对比不同工况下齿轮副齿面载荷分布和啮合效率。试验结果表明:系统变形对准双曲面齿轮齿面载荷分布具有不可忽略的影响。同时, 该试验验证了驱动桥准双曲面齿轮啮合效率计算方法的有效性。
关键词驱动桥    准双曲面齿轮    接触分析    啮合效率    混合润滑    系统变形    
Calculation method of hypoid gear meshing efficiency of drive axles with considering system deformation
WANG Qin1, HE Di1, GUI Liangjin1, HU Zhiyu2, PENG Jin3, FAN Zijie1    
1. State Key Laboratory of Automotive Safety and Energy, School of Vehicle and Mobility, Tsinghua University, Beijing 100084, China;
2. Shaanxi Hande Axle Co., Ltd., Xi'an 710201, China;
3. Hande Axle(Zhuzhou) Gear Co., Ltd., Zhuzhou 412000, China
Abstract: [Objective] The face-hobbing hypoid gear is a crucial component of drive axles owing to its continuous indexing processing capabilities. However, in practical engineering applications, gear pairs rarely operate under ideal conditions. Instead, they frequently experience heavy-load and high-speed conditions. These conditions result in substantial system deformation and a complex meshing state of the tooth surface. Under heavy-load conditions, considerable skewing of the load on the tooth surface can occur, greatly affecting the performance of the gear pair and drive axles. Currently, the system deformation factor under actual working conditions is not sufficiently considered, making it challenging to address the aforementioned issue. Consequently, this paper proposes the analysis of the tooth surface load distribution by employing a semi-analytic loaded tooth contact analysis method to accurately predict the tooth surface load distribution. Based on a load distribution analysis, a highly accurate calculation method for gear meshing efficiency is proposed. [Methods] This paper proposes a calculation method for gear meshing efficiency under mixed-lubrication conditions for hypoid gears in drive axles operating under complex working conditions. First, a multi-support shaft system modeling method is employed to analyze the drive axles system. This method can calculate the forces acting on various components, such as gears and bearings, as well as the gear misalignment caused by system deformation under various load conditions. Second, by simulating the spatial motion process of the actual gear machining machine, the coordinates of the tool cutting point are transformed to the coordinate system of the gear blank via coordinate transformation. This process results in the correspondence of the tooth profile with the actual machining parameters. The time-varying friction coefficient distribution of the tooth surface under different working conditions is derived by combining the point contact mixed-lubrication friction coefficient model of the tooth surface with its relative motion relationship. Then, taking into account the tooth surface deformation equilibrium equation, tooth surface torque equilibrium equation, and tooth surface contact pressure equilibrium equation, a gear frictional loaded tooth contact analysis method is established. This method accurately calculates the tooth surface load distribution and mesh efficiency of gears under different working conditions through an iterative solution. Finally, the calculation results of the tooth surface load distribution under various working conditions are compared with the experimental results obtained from a loading experiment conducted on the entire drive axles. The meshing efficiency of the gear pair under various working conditions is determined by conducting a system no-load efficiency experiment and loading efficiency experiment and comparing the results with those obtained by calculations. [Results] In the gear no-load experiment, the contact patterns of the gear tooth surface were compared under forward and reverse working conditions. The experimental results were found to be in good agreement with theoretical calculations, verifying the accuracy of the tooth surface calculation method and no-load tooth contact analysis. Subsequently, a loading experiment was conducted on the drive axles system, indicating that the system deformation had a considerable impact on the load distribution of the hypoid gear tooth surface. For instance, in the experiment drive axles, under heavy-load conditions, system deformation caused the tooth surface load on the driving side of the gear pair to shift toward the outside, while the tooth surface loaded on the driven side shifts toward the inside. The meshing efficiency experiment results revealed that system deformation considerably impacted the gear meshing efficiency under heavy-load and high-speed working conditions. In addition, the vehicle speed had a considerable impact on the meshing efficiency, with an increase in speed from 10 to 80 km/h, resulting in a 1% improvement in meshing efficiency. [Conclusions] By comprehensively considering gear meshing misalignment caused by system deformation and the mixed-lubrication state of the tooth surface, the tooth surface load distribution and meshing efficiency can be accurately calculated under loaded conditions. Therefore, the proposed process enhances the accuracy of calculating the drive axles system efficiency while providing a solid foundation for gear optimization research.
Key words: drive axles    hypoid gear    contact analysis    meshing efficiency    mixed lubrication    system deformation    

摆线准双曲面齿轮因其连续分度加工的特点,成为驱动桥系统的核心零部件之一[1]。传动效率是驱动桥的关键性能指标之一,高负载工况下齿轮啮合功耗是驱动桥系统总功耗的主要组成部分[2]。因此,准确计算齿轮啮合效率对提高驱动桥整桥品质有重要的意义。

目前,有关锥齿轮啮合效率的研究已得到广泛的关注。Xu等[3]基于回归分析方法得到了流体润滑状态下齿面摩擦因数模型,并提出了准双曲面齿轮啮合效率计算方法。在此基础上,Kolivand等[4]提出了基于混合弹流润滑理论的准双曲面齿轮啮合效率计算方法。Simon[5]分析了混合润滑状态下准双曲面齿轮加工参数对齿轮啮合效率的影响,并给出了改进建议。Ding等[6]基于热弹流润滑理论,提出了弧齿锥齿轮啮合效率的计算方法。谷建功等[7]推导了弧齿锥齿轮混合弹流润滑状态下的时变摩擦功率损失和啮合效率计算公式。曹伟等[8]针对弧齿锥齿轮,分析了一个啮合周期内齿轮时变摩擦因数和啮合效率的变化情况。蒋进科等[9]采用混合弹流润滑摩擦因数,基于齿面拓扑修形方法获得了提高齿轮啮合效率的最佳修形面。魏冰阳等[10]基于齿面修形设计并结合混合弹流润滑理论给出了齿轮时变摩擦功耗计算方法。唐效贵等[11]考虑直齿锥齿轮齿面复杂润滑状态,建立了瞬时功率损耗和啮合效率的计算公式。针对驱动桥系统效率分析,Kakavas等[12]考虑了驱动桥准双曲面齿轮摩擦、拖曳与轴承等因素的影响,提出了系统效率计算方法,并与试验结果进行了对比。王旺平[13]以微型汽车传动系统为研究对象,分析了影响系统传动效率的主要因素,并提出了系统效率经验计算方法。

齿轮啮合效率研究成果丰富,在理论层面上提供了不同的计算方法,然而现有研究仍存在不完善的地方。研究人员多数是以齿轮个体为研究对象,对实际使用工况下系统变形这一因素的考虑尚不充分。由于试验条件限制,准双曲面齿轮啮合效率难以直接测量,导致理论计算与试验结果难以直接进行对比。为解决上述问题,本文提出了一种考虑系统变形的驱动桥准双曲面齿轮啮合效率计算方法。首先,综合考虑齿轮、轴承、轴和驱动桥壳体等因素的影响,建立多支撑轴系耦合分析模型,计算得到系统各个部件受力和齿轮啮合错位量。然后,基于齿轮副的实际几何参数和机床加工参数,通过模拟机床运动得到实际齿面,并采用FLTCA方法计算实际工况下齿面载荷分布与啮合效率。最后,采用驱动桥加载试验台架,对不同工况下的齿面载荷分布与啮合效率进行了测试。通过与试验结果进行对比,验证了本文计算方法的有效性。

1 齿轮摩擦加载接触分析

齿轮摩擦加载接触分析(frictional loaded tooth contact analysis, FLTCA)是计算齿轮副啮合性能的有效方法,对此已有广泛研究。在计算齿轮啮合效率前,需精确计算每个啮合时刻齿轮副的齿面载荷分布和齿面摩擦因数。针对驱动桥实际结构,本节提出了一种考虑系统变形的齿轮摩擦加载接触分析方法。

1.1 驱动桥系统变形与齿轮啮合错位量

驱动桥系统结构复杂,包括轴、轴承、差速器壳体、桥壳和齿轮等多个部件[14]。工程实践表明:在受载状态下,驱动桥系统产生变形,导致齿轮副无法在理想的啮合状态下正常工作,齿轮副之间发生错位。错位将导致齿轮副产生偏载、脱出等不良工况,极端情况下甚至导致齿轮断裂,进而造成驱动桥的损坏。

多支撑轴系耦合分析模型[15]是计算驱动桥系统变形的有效方法。图 1为某型号的商用车后驱动桥结构示意图。其中,主动齿轮由前轴承和中轴承支承,被动齿轮固定于差速器壳体上,差速器齿轮位于差速器壳体中,并通过花键与左、右输出半轴相连。基于多支撑轴系耦合分析模型,建立系统刚度方程:

$ \boldsymbol{K} \cdot \boldsymbol{\delta}=\boldsymbol{F}. $ (1)
图 1 某商用车后驱动桥系统结构

其中:K表示由轴、轴承、齿轮和壳体等零部件的刚度矩阵组集得到的系统刚度矩阵;δ表示系统节点位移;F表示系统所受外力。求解上述方程可得到系统在任意节点处的位移和节点力。

图 2为一对准双曲面齿轮啮合错位示意图。由式(1)计算得到主动齿轮中心节点处位移δp和从动齿轮中心节点处位移δg

$ \boldsymbol{\delta}_{\mathrm{p}}=\left(\delta_{\mathrm{p} x}, \delta_{\mathrm{p} y}, \delta_{\mathrm{p} z}, \theta_{\mathrm{p} x}, \theta_{\mathrm{p} y}, \theta_{\mathrm{p} z}\right), $ (2)
$ \boldsymbol{\delta}_{\mathrm{g}}=\left(\delta_{\mathrm{g} x}, \delta_{\mathrm{g} y}, \delta_{\mathrm{g} z}, \theta_{\mathrm{g} x}, \theta_{\mathrm{g} y}, \theta_{\mathrm{g} z}\right) . $ (3)
图 2 准双曲面齿轮啮合错位量

其中:δpxδpyδpzθpxθpyθpz分别表示主动轮中心节点在xyz方向的平动位移和角位移;δgxδgyδgzθgxθgyθgz分别表示从动轮中心节点在xyz方向的平动位移和角位移。由上述节点位移可求得主动轮和从动轮错位量:

$ \left\{\begin{array}{l} {\Delta} P_{1}=\delta_{\mathrm{p} x}, \\ {\Delta} W_{1}=-\delta_{\mathrm{p} y}+\overline{\mathrm{WP}} \cos \alpha \sin \theta_{\mathrm{p} z}, \\ {\Delta} Y_{1}=-\delta_{\mathrm{p} z}-\overline{\mathrm{WP}} \cos \alpha \sin \theta_{\mathrm{p} y}, \\ {\Delta} \mathit{\Sigma}_{1}=-\theta_{\mathrm{p} z}, \end{array}\right. $ (4)
$ \left\{\begin{array}{l} {\Delta} P_{2}=-\delta_{\mathrm{g} x}-\overline{\mathrm{GP}} \cos \beta \sin \theta_{\mathrm{g} z}, \\ {\Delta} W_{2}=\delta_{\mathrm{g} y}, \\ {\Delta} Y_{2}=\delta_{\mathrm{g} z}-\overline{\mathrm{GP}} \cos \beta \sin \delta_{\mathrm{g} x}, \\ {\Delta} \mathit{\Sigma}_{2}=-\theta_{\mathrm{g} z}. \end{array}\right. $ (5)

其中:ΔP1和ΔP2分别表示主动齿轮和从动齿轮沿主动齿轮轴线方向的位移;ΔW1和ΔW2分别表示主动齿轮和从动齿轮沿从动齿轮轴线方向的位移;ΔY1和ΔY2分别表示主动齿轮和从动齿轮沿齿轮偏置方向的位移;ΔΣ1和ΔΣ2分别表示主动齿轮和从动齿轮沿轴交角方向的角位移;αβ分别表示主动齿轮和从动齿轮的偏置角;$ \overline{\mathrm{WP}}$$ \overline{\mathrm{GP}}$分别表示齿轮啮合点到主动齿轮和从动齿轮中心的距离[16]

进一步计算齿轮副总错位量。主、从动齿轮沿主动齿轮轴线方向的相对位移ΔP、沿从动齿轮轴线方向的相对位移ΔW、沿齿轮偏置方向的相对位移ΔY、沿齿轮轴交角方向的相对角位移ΔΣ分别为

$ \begin{aligned} {\Delta} P & ={\Delta} P_{1}+{\Delta} P_{2}, \\ {\Delta} W & ={\Delta} W_{1}+{\Delta} W_{2}, \\ {\Delta} Y & ={\Delta} Y_{1}+{\Delta} Y_{2}, \\ {\Delta} \mathit{\Sigma} & ={\Delta} \mathit{\Sigma}_{1}+{\Delta} \mathit{\Sigma}_{2}. \end{aligned} $ (6)
1.2 准双曲面齿轮齿面

利用机床加工参数,可以模拟刀具在加工过程中的运动,从而获得摆线准双曲面齿轮齿面点的空间坐标[17],这些坐标通过坐标系之间的变换,将刀盘上刀具刃面点坐标变换至轮坯坐标系中。图 3为以主动齿轮为例建立加工机床坐标系Om-imjmkm、刀盘坐标系Oc-icjckc和轮坯坐标系Olp-ilpjlpklp之间空间位置关系[18]Ob-ibjbkb表示轮坯轴线ilp和机床坐标下jm-km平面的交点所建立的坐标系。点A表示齿面参考点,X1表示水平轮位,Lcz表示垂直轮位,XB1表示床位,δM1表示机床安装根锥角,s1表示径向刀位,q1表示角向刀位。

图 3 机床加工示意图

准双曲面齿轮齿面几何形状由机床加工参数确定,轮坯坐标下齿面点径向和法向坐标可表示为机床加工参数的相关函数[19]

$ \left\{\begin{array}{l} \boldsymbol{r}=f_{r}\left(\theta_{\mathrm{d}}, \varphi_{\mathrm{d}}, \boldsymbol{\xi}_{\mathrm{c}}\right) ;\\ \boldsymbol{n}=f_{n}\left(\theta_{\mathrm{d}}, \varphi_{\mathrm{d}}, \boldsymbol{\xi}_{\mathrm{c}}\right). \end{array}\right. $ (7)

其中:θd表示刀盘绕刀盘中心转角,φd表示刀盘中心绕摇台转角, ξc表示机床加工参数向量。根据式(7)得到离散的齿面点,采用Ferguson曲线三次样条函数进行曲面拟合,拟合后的齿面用于空载和加载接触分析。

1.3 混合润滑摩擦因数

高速重载驱动桥准双曲面齿轮工况恶劣,齿面啮合点处于点接触混合润滑状态,本文采用基于加权函数的摩擦因数对啮合齿面摩擦因数μML进行描述[9, 20]

$ \mu_{\mathrm{ML}}=\mu_{\mathrm{FL}} f_{\mathrm{K}}^{1.2}+\mu_{\mathrm{DC}}\left(1-f_{\mathrm{K}}\right), $ (8)
$ f_{\mathrm{K}}=\frac{1.21 \lambda^{0.64}}{1+0.37 \lambda^{1.26}}, $ (9)
$ \lambda=\frac{h_{0}}{S}, $ (10)
$ S=\sqrt{S_{\mathrm{r} 1}^{2}+S_{\mathrm{r} 2}^{2}} . $ (11)

其中:μFL示流体润滑摩擦因数;μDC表示边界润滑摩擦因数;fK表示接触载荷承载因数;λ表示油膜比厚;Sr1Sr2分别表示主动齿轮和从动齿轮的齿面粗糙度;S表示等效粗糙度;h0表示啮合处中心润滑油油膜厚度[21]h0可表示如下:

$ h_{0}=2.69 W^{-0.067} U^{0.67} G^{0.53}\left(1-0.61 \mathrm{e}^{-0.73 l}\right), $ (12)
$ \begin{gathered} W=\frac{F_{\mathrm{m}}}{E_{\mathrm{eq}} R_{x}^{2}}, \\ U=\frac{V_{\mathrm{e}} u_{0}}{E_{\mathrm{eq}} R_{x}}, \\ G=\alpha E_{\mathrm{eq}} ,\\ E_{\mathrm{eq}}=2\left[\frac{1-v_{1}^{2}}{E_{1}}+\frac{1-v_{2}^{2}}{E_{2}}\right]^{-1}, \\ l=1.03\left(\frac{R_{x}}{R_{y}}\right)^{0.64} . \end{gathered} $ (13)

其中:W表示无量纲载荷因数;U表示无量纲速度因数;G表示无量纲材料因数;Eeq表示材料等效弹性模量;l表示接触点椭圆率;RxRy分别表示沿接触椭圆长轴和短轴方向的当量曲率半径,m;v1v2分别表示主动齿轮和从动齿轮材料Poisson比;E1E2分别表示主动齿轮和从动齿轮材料弹性模量, Pa;α表示润滑油黏压系数,Pa-1u0表示润滑油黏度, Pa·s;Ve表示卷吸速度, m/s;Fm表示润滑膜承载量, N。

流体润滑摩擦因数为[22]

$ \begin{array}{c} \mu_{\mathrm{FL}}=\mathrm{e}^{f\left(\mathrm{SR}, \mathrm{Ph}, u_{0}, S\right)}(\mathrm{Ph} / \mathrm{MPa})^{b_{2}}|\mathrm{SR}|^{b_{3}}\left(V_{\mathrm{e}} /\right. \\ \left.\left(\mathrm{m} \cdot \mathrm{s}^{-1}\right)\right)^{b_{6}} u_{0}^{b_{7}} R^{b_{8}}, \end{array} $ (14)
$ \begin{array}{c} f\left(\mathrm{SR}, \mathrm{Ph}, u_{0}, \mathrm{~S}\right)= \\ b_{1}+b_{4}|\mathrm{SR}|(\mathrm{Ph} / \mathrm{MPa}) \lg \left(u_{0} /(\mathrm{Pa} \cdot \mathrm{s})\right)+ \\ b_{5} \mathrm{e}^{-|\mathrm{SR}|(\mathrm{Ph} / \mathrm{MPa}) \lg \left(u_{0} /(\mathrm{Pa} \cdot \mathrm{s})\right)}+b_{9} \mathrm{e}^{{S}}, \end{array} $ (15)
$ \mathrm{SR}=2 V_{\mathrm{s}} / V_{\mathrm{r}} . $ (16)

其中:b1~b9表示方程回归系数;Ph表示Hertz接触压力;SR表示齿面速度滑滚比;VsVr分别表示齿面啮合点处的滑动速度和滚动速度。摩擦因数模型回归系数的值如表 1所示。

表 1 摩擦因数模型回归系数
系数 数值 系数 数值 系数 数值
b1 -8.916 465 b4 -0.354 068 b7 0.752 755
b2 1.033 03 b5 2.812 084 b8 -0.390 958
b3 1.036 077 b6 -0.100 601 b9 0.620 305

1.4 驱动桥系统齿轮FLTCA方法

图 4为驱动桥系统准双曲面齿轮FLTCA计算方法流程示意图。具体步骤如下。

图 4 驱动桥系统齿轮FLTCA方法流程图

步骤1  建立并求解驱动桥系统整桥静力学模型,考虑齿轮、轴承、轴和壳体等零部件的影响,得到系统内齿轮啮合错位量和轴承受力。

步骤2  根据齿轮基本几何参数和齿轮副加工参数计算主动齿轮和从动齿轮的实际齿面。

步骤3  确定润滑油基本参数和齿轮齿面粗糙度。

步骤4  考虑不同工况下系统变形,对齿轮副进行空载接触分析(tooth contact analysis, TCA),得到理论状态下的齿面接触点和空载接触印迹。

步骤5  计算齿轮受载后的弯曲变形δb[23]、齿轮剪切变形δs和接触变形δc[24],并判断是否满足变形协调方程:

$ \left\|\boldsymbol{\delta}_{\mathrm{b}}+\boldsymbol{\delta}_{\mathrm{s}}+\boldsymbol{\delta}_{\mathrm{c}}-\left(\boldsymbol{Z}-\boldsymbol{d}_{0}\right)\right\|<\varepsilon_{1} . $ (17)

其中:Z表示齿轮刚体转动位移,d0表示齿面接触点初始间距。

步骤6  计算每个啮合时刻的齿面接触力,并判断是否满足载荷平衡方程:

$ \left|\sum\limits_{i=1}^{m} \sum\limits_{j=1}^{n}\left(\boldsymbol{F}_{\mathrm{fij}}+\boldsymbol{F}_{\mathrm{N} i j}\right) \boldsymbol{r}_{i j} \boldsymbol{p}-\boldsymbol{T}_{\mathrm{load}}\right|<\varepsilon_{2}. $ (18)

其中:FN表示齿面接触法向力; Ff表示齿面摩擦力; m表示同一啮合时刻下同时参与啮合的齿数,n表示将接触椭圆沿长轴方向的离散点数; p表示从动轮转动方向矢量; Tload表示齿轮副负载扭矩; ε1表示变形协调方程收敛容差;ε2表示载荷平衡方程收敛容差。

步骤7  由Hertz接触公式计算齿面接触压力分布,并判断接触压力是否收敛。若不收敛,则重新计算接触压力分布值和齿面摩擦因数,并返回步骤5。

步骤8  计算完成,得到齿面接触力、齿面接触印迹和齿面摩擦因数等结果。

2 啮合效率

本文采用驱动桥传动效率试验台架进行试验,控制逻辑为:输入端控制转速,输出端控制扭矩。因此,齿轮啮合效率ηg, fric可表示为:

$ \eta_{\mathrm{g}, \text { fric }}=\frac{p_{\text {out }}}{p_{\text {out }}+p_{\mathrm{g}, \text { fric }}}=\frac{\left|\boldsymbol{T}_{\text {out }}\right| \omega_{\text {out }}}{\left|\boldsymbol{T}_{\text {out }}\right| \omega_{\text {out }}+p_{\mathrm{g}, \text { fric }}}, $ (19)
$ p_{\mathrm{g}, \text { fric }}=\frac{1}{\theta_{1}-\theta_{2}} \int_{\theta_{1}}^{\theta_{2}} p_{\text {inst }} \mathrm{d} \theta . $ (20)

其中:pout表示输出功率;Tout表示输出扭矩;ωout表示输出转速;pg,fric为齿轮副总啮合功率损失; θ1θ2分别表示齿轮副一个啮合周期的啮入转角和啮出转角; pinst表示每个啮合时刻瞬时功率损失。pinst可表示为:

$ p_{\mathrm{inst}}=\sum\limits_{i=1}^{m} \sum\limits_{j=1}^{n}\left|\boldsymbol{F}_{\mathrm{f} i j}\right|\left|\boldsymbol{V}_{\mathrm{s} i j}\right|. $ (21)

其中Vsij表示齿面相对滑动速度。

3 试验验证

本文对驱动桥准双曲面齿轮的接触特性和啮合效率进行了试验验证。图 5展示了用于试验的驱动桥整桥加载试验台结构。

图 5 驱动桥整桥加载试验台结构

表 25分别给出了试验驱动桥准双曲面齿轮基本参数、机床加工主动轮加工参数、机床加工从动轮加工参数和轴承参数。

表 2 齿轮基本参数
基本参数 主动轮 从动轮
齿数 10 39
轴交角/(°) 90.00 90.00
平均压力角/(°) 22.50 22.50
偏置距/mm 30.00 30.00
齿顶高因数 0.975 0.975
齿根高因数 1.225 1.225
齿宽/mm 55.80 50.50
节圆直径/mm 77.00 300.00
中点螺旋角/(°) 44.00 30.80
旋向 左旋 右旋
刀盘名义半径/mm 125.00 125.00
刀组数 17 17
齿面平均粗糙度/μm 0.80 0.80
齿轮侧隙/mm 0.30 0.30

表 3 主动轮加工参数
加工参数 凹面 凸面
角向刀位/(°) 56.03 56.03
径向刀位/mm 162.08 162.08
垂直轮位/mm 26.17 26.17
水平轮位/mm -1.34 -1.34
机床安装根锥角/(°) -2.04 -2.04
滚比 3.73 3.73
床位/mm 24.34 24.34
刀转角/(°) 335.06 335.06
刀倾总角/(°) 23.62 23.62
刀具齿形角/(°) 26.26 18.80
参考点半径/mm 115.78 116.12
参考点高度/mm 6.46 6.46
圆弧刀具半径/mm 5 714.02 3 132.78
刀具修缘角度/(°) 3.50 3.50
刀具修缘高度/mm 3.50 3.50

表 4 从动轮加工参数
加工参数 凸面 凹面
水平刀位/mm 113.74 113.74
垂直刀位/mm 123.56 123.56
水平轮位/mm 3.63 3.63
床位/mm 0 0
机床安装根锥角/(°) 70.27 70.27
刀具齿形角/(°) 20.40 24.60
参考点半径/mm 115.45 116.03
参考点高度/mm 6.78 6.78
圆弧刀具半径/mm 5 000.00 1 500.00
刀具修缘角度/(°) 3.50 3.50
刀具修缘高度/mm 3.50 3.50

表 5 轴承参数
轴承参数 前轴承 中轴承 左轴承 右轴承
内径/mm 50.00 60.00 85.00 70.00
外径/mm 105.00 130.00 130.00 110.00
轴承宽/mm 34.50 43.00 29.00 25.00
接触角/(°) 30.00 30.00 16.42 16.17

3.1 系统加载试验

首先进行驱动桥加载试验,验证理论计算接触印迹和实际齿轮接触印迹的一致性。试验包括正车试验和反车试验。正车试验工况包括空载和20%、40%、60%、80%、100%满载;反车试验工况包括空载和20%、40%、60%满载。计算不同满载工况下的齿轮副错位量,结果如表 6所示。

表 6 不同满载工况下齿轮副错位量
工况 ΔP/μm ΔW/μm ΔY/μm ΔΣ/μrad
正车20%满载 101.6 -17.8 -128.6 0.22
正车40%满载 191.7 -31.7 -236.6 0.38
正车60%满载 279.8 -44.9 -341.3 0.51
正车80%满载 366.8 -57.7 -444.1 0.63
正车100%满载 453.2 -70.3 -545.9 0.75
反车20%满载 -173.9 231.1 173.7 0.21
反车40%满载 -326.6 411.4 312.2 0.49
反车60%满载 -475.7 583.5 445.9 0.75

为了直观展示驱动桥加载过程中的接触印迹,在加载试验时去除部分驱动桥桥包,并使部分齿轮裸露在外面。图 67分别对比了驱动桥在空载状态下正车和反车的齿面接触印痕。试验与理论计算结果较吻合,验证了齿面计算和空载接触计算的准确性。

图 6 正车空载齿面接触印痕

图 7 反车空载齿面接触印痕

为保障齿轮的寿命和驱动桥的使用安全,工程实践中通常要求在正车满载和反车60%满载状态下,确保驱动桥准双曲面齿轮的齿轮副接触印迹不在小端脱出。图 89分别对比了上述2个工况下接触印迹的计算和试验结果。在试验时,驱动桥主锥采用双轴承的形式进行布局,系统错位量大,满载工况下齿面接触印迹向齿顶靠近。可以看出,在不考虑错位量时,齿面接触印迹和试验结果差异较大;在考虑错位量后,正、反车工况下加载接触印迹分布与试验结果较一致。

图 8 正车满载接触印迹

图 9 反车60%满载接触印迹

试验结果表明:本文采用考虑系统变形的齿轮接触计算方法,能够准确地计算齿轮副在实际使用工况下的接触状态。

3.2 啮合效率试验

对驱动桥系统不同工况下的系统效率进行试验,试验车速范围10~80 km/h,试验负载功率范围10~100 kW,试验润滑油牌号85W90,试验温度通过台架温控系统控制在(80±5)℃。在该温度下,润滑油动力黏度为27 Pa·s,黏压系数为2.2×1010 Pa-1,边界润滑摩擦因数取0.115[25],润滑油液面高度为从动轮中心下方26 mm。为确保试验结果的可靠性,选取2组相同型号的驱动桥进行重复试验,并对试验结果取平均值。

图 10中,驱动桥整桥功率损耗ploss由负载相关功率损耗ploss, fric和负载无关功率损耗ploss, drag构成[12]。其中,负载相关功率损耗包括齿轮摩擦功耗pg, fric和轴承摩擦功耗pb, fric。负载无关功耗主要由驱动桥系统内齿轮搅油功耗pg, drag和轴承搅油功耗pb, drag等组成。负载无关功耗与部件尺寸、转速、润滑油属性等因素相关,通常不随载荷变化。因此,驱动桥系统总效率η和齿轮啮合效率ηg, firc分别表示如下[26]

$ \begin{gathered} \eta=1-\frac{p_{\mathrm{g}, \text { drag }}+p_{\mathrm{b}, \text { drag }}+p_{\mathrm{g}, \text { fric }}+p_{\mathrm{b}, \text { fric }}}{p_{\text {in }}}= \\ \eta_{\text {drag }}-\left(1-\eta_{\mathrm{g}, \text { fric }}\right)-\frac{p_{\mathrm{b}, \text { fric }}}{p_{\text {in }}}, \end{gathered} $ (22)
$ \eta_{\mathrm{g}, \text { fric }}=1+\eta-\eta_{\text {drag }}+\frac{p_{\mathrm{b}, \text { fric }}}{p_{\text {in }}} . $ (23)
图 10 驱动桥整桥功率损耗构成

其中:ηdrag表示驱动桥负载无关功耗对应传动效率, pin表示驱动桥输入功率。采用SKF方法计算pb,fric[27]

$ p_{\mathrm{b}, \text { fric }}=M_{\mathrm{b}} \omega_{\mathrm{b}} , $ (24)
$ M_{\mathrm{b}}=\phi_{\text {ish }} \phi_{\mathrm{rs}} M_{\mathrm{rr}}+M_{\mathrm{sl}}+M_{\text {seal }} . $ (25)

其中:Mb表示轴承摩擦力矩;ωb表示轴承转速;ϕish表示切入发热减少系数;ϕrs表示运动贫油减少系数;Mrr表示滚动摩擦力矩;Msl表示滑动摩擦力矩;Mseal表示密封件摩擦力矩。

齿轮啮合效率的试验步骤如下。

步骤1  试验负载功率为零,测试不同试验车速条件下驱动桥输入功率,该功率即不同车速下负载无关功率损耗试验值,包括pg,dragpb,drag

步骤2  在试验车速和试验负载功率范围内,依次测量驱动桥输入功率,得到不同工况下驱动桥系统的ηηdrag

步骤3  由式(23)计算得到不同工况下ηg, fric

图 11展示了在不同车速下驱动桥系统在负载功率80 kW时各项损耗的构成,包括齿轮摩擦损耗、轴承摩擦损耗和负载无关损耗。结果表明:在低速工况下,齿轮摩擦损耗是系统主要的功率损耗构成部分。随着车速的增加,齿轮摩擦损耗减小。然而,轴承摩擦损耗和负载无关损耗随车速的增加而增大。

图 11 负载功率80 kW时不同车速系统功率损耗构成

图 12对比了负载功率80 kW时,齿轮啮合效率随车速的变化趋势。可知随着车速的增加,啮合效率呈非线性上升趋势。当车速从10 km/h增加至80 km/h时,齿轮啮合效率增加约1%。图 13对比了负载力矩为1 354 N·m时,齿轮啮合效率随车速的变化趋势。可知啮合效率变化趋势和恒功率工况的规律类似,随着车速的增加,啮合效率非线性增加。当车速从10 km/h增加至80 km/h时,齿轮啮合效率增加约1.5%。图 14对比了车速60 km/h时,齿轮啮合效率随负载功率的变化趋势,可知啮合效率随负载功率变化不显著。上述结果表明车速和负载均影响齿面润滑状态。图 15对比了负载功率80 kW时一个啮合周期内不同车速下齿面油膜厚度h0和齿面摩擦因数μML的计算结果。随着车速的增加,齿面载荷减小,齿面相对滑动速度升高,齿面润滑油油膜厚度增加,齿面摩擦因数下降。

图 12 负载功率80 kW时不同车速啮合效率

图 13 负载力矩1 354 N·m时不同车速啮合效率

图 14 车速60 km/h不同负载功率啮合效率

图 15 负载功率80 kW时平均油膜厚度与平均摩擦因数

图 1214可知,齿轮啮合过程的错位影响啮合效率。在负载恒定的情况下,相比高速工况,低速工况下齿面相对滑动速度变化较小,考虑与不考虑齿轮啮合错位量计算的结果相差不大。随着车速的增加,齿面相对滑动速度增大,考虑齿轮啮合错位量后的计算结果与试验结果较接近。在车速恒定的情况下,随着负载的增加,系统错位量增大,齿面载荷分布和齿面相对滑动速度均发生变化,考虑齿轮啮合错位量后的计算结果和试验结果较接近。因此,对于高速重载驱动桥准双曲面齿轮齿面载荷分布和啮合效率的计算,考虑系统变形后的计算结果更准确。

综上,本文提出的计算方法与试验结果吻合度较高,可有效计算驱动桥齿准双曲面齿轮啮合效率。

4 结论

本文基于齿轮FLTCA方法,计算了实际工况下驱动桥准双曲面齿轮的齿面载荷分布,并提出了齿轮啮合效率计算方法。通过考虑系统变形引起的齿轮啮合错位,可以准确地计算真实状态下的齿面载荷分布。因此,本文采用摆线准双曲面齿轮机床运动过程模拟实际齿面,考虑驱动桥系统变形对实际工况下齿轮啮合状态的影响并进行接触分析。同时,采用系统加载试验验证了齿面载荷分布计算方法的准确性。在此基础上,综合考虑齿面摩擦因数分布和啮合点相对运动关系,对齿轮啮合效率进行了计算。驱动桥齿轮啮合效率的试验结果验证了本文提出的计算方法的有效性,为下一步齿轮设计优化提供了基础。

参考文献
[1]
杜进辅. 摆线齿准双曲面齿轮齿面啮合性能分析及设计技术研究[D]. 西安: 西北工业大学, 2015.
DU J F. Tooth contact pattern analysis and design technology research of cycloid hypoid gears[D]. Xi'an: Northwestern Polytechnical University, 2015. (in Chinese)
[2]
TALBOT D, KAHRAMAN A, LI S, et al. Development and validation of an automotive axle power loss model[J]. Tribology Transactions, 2016, 59(4): 707-719. DOI:10.1080/10402004.2015.1110862
[3]
XU H, KAHRAMAN A. Prediction of friction-related power losses of hypoid gear pairs[J]. Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-Body Dynamics, 2007, 221(3): 387-400. DOI:10.1243/14644193JMBD48
[4]
KOLIVAND M, LI S, KAHRAMAN A. Prediction of mechanical gear mesh efficiency of hypoid gear pairs[J]. Mechanism and Machine Theory, 2010, 45(11): 1568-1582. DOI:10.1016/j.mechmachtheory.2010.06.015
[5]
SIMON V V. Improvements in the mixed elastohy-drodynamic lubrication and in the efficiency of hypoid gears[J]. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2020, 234(6): 795-810. DOI:10.1177/1350650119866027
[6]
DING H, LI H P, CHEN S Y, et al. Energy loss and mechanical efficiency forecasting model for aero-engine bevel gear power transmission[J]. International Journal of Mechanical Sciences, 2022, 231: 107569. DOI:10.1016/j.ijmecsci.2022.107569
[7]
谷建功, 方宗德, 苏进展, 等. 混合弹流润滑下弧齿锥齿轮传动啮合效率计算方法[J]. 农业机械学报, 2010, 41(5): 188-192.
GU J G, FANG Z D, SU J Z, et al. Calculation of meshing efficiency for spiral bevel gears under the condition of mixed elastohydrodynamic lubrication[J]. Transactions of the Chinese Society for Agricultural Machinery, 2010, 41(5): 188-192. DOI:10.3969/j.issn.1000-1298.2010.05.038 (in Chinese)
[8]
曹伟, 蒲伟, 王家序, 等. 弧齿锥齿轮摩擦系数与啮合效率研究[J]. 摩擦学学报, 2018, 38(3): 247-255.
CAO W, PU W, WANG J X, et al. Friction coefficient and mesh efficiency in spiral bevel gears[J]. Tribology, 2018, 38(3): 247-255. DOI:10.16078/j.tribology.2018.03.001 (in Chinese)
[9]
蒋进科, 方宗德, 刘钊. Ease-off拓扑修形准双曲面齿轮齿面多目标优化设计方法[J]. 西安交通大学学报, 2019, 53(6): 44-53.
JIANG J K, FANG Z D, LIU Z. Design of multi-objective tooth optimization for hypoid gear with ease-off topological modification[J]. Journal of Xi'an Jiaotong University, 2019, 53(6): 44-53. (in Chinese)
[10]
魏冰阳, 王振, 杨建军, 等. Ease-off拓扑修形齿面拟Hertz接触与摩擦特性分析[J]. 机械工程学报, 2021, 57(1): 61-67.
WEI B Y, WANG Z, YANG J J, et al. Analysis on quasi Hertzian contact and friction characteristics of tooth surface modified by ease-off topology[J]. Journal of Mechanical Engineering, 2021, 57(1): 61-67. (in Chinese)
[11]
唐效贵, 王斌, 李佩鸿, 等. 基于复杂润滑状态的直齿锥齿轮啮合效率分析[J]. 机械传动, 2017, 41(7): 33-41.
TANG X G, WANG B, LI P H, et al. Analysis of meshing efficiency of straigh bevel gear based on integrated lubrication[J]. Journal of Mechanical Transmission, 2017, 41(7): 33-41. (in Chinese)
[12]
KAKAVAS I, OLVER A V, DINI D. Hypoid gear vehicle axle efficiency[J]. Tribology International, 2016, 101: 314-323. DOI:10.1016/j.triboint.2016.04.030
[13]
王旺平. 微型汽车传动系统效率建模及试验研究[D]. 武汉: 武汉理工大学, 2015.
WANG W P. Efficiency modeling and experimental research of mini-car transmission system[D]. Wuhan: Wuhan University of Technology, 2015. (in Chinese)
[14]
ZHOU C, WANG Q, DING W Q, et al. Numerical simulation of drive axles considering the nonlinear couplings between gears and bearings[J]. Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, 2019, 233(4): 851-861. DOI:10.1177/0954407018755598
[15]
田程, 周驰, 桂良进, 等. 考虑轴承刚度耦合性和非线性的多支撑轴系有限元分析方法[J]. 机械工程学报, 2015, 51(17): 90-95.
TIAN C, ZHOU C, GUI L J, et al. Method of finite element analysis for multi-support shaft system based on the coupling and nonlinearity of bearing stiffness[J]. Journal of Mechanical Engineering, 2015, 51(17): 90-95. (in Chinese)
[16]
田程, 周驰, 丁炜琦, 等. 轴系的热膨胀对于锥齿轮错位量的影响[J]. 清华大学学报(自然科学版), 2016, 56(6): 565-571.
TIAN C, ZHOU C, DING W Q, et al. Influence of the thermal expansion of a shaft on the misalignment of bevel gears[J]. Journal of Tsinghua University (Science and Technology), 2016, 56(6): 565-571. (in Chinese)
[17]
DING H, TANG J Y, ZHONG J, et al. A hybrid modification approach of machine-tool setting considering high tooth contact performance in spiral bevel and hypoid gears[J]. Journal of Manufacturing Systems, 2016, 41: 228-238.
[18]
WANG Q, ZHOU C, GUI L J, et al. Optimization of the loaded contact pattern of spiral bevel and hypoid gears based on a kriging model[J]. Mechanism and Machine Theory, 2018, 122: 432-449.
[19]
田程, 丁炜琦, 桂良进, 等. 基于回归分析的准双曲面齿轮齿面误差修正[J]. 清华大学学报(自然科学版), 2017, 57(2): 141-146.
TIAN C, DING W Q, GUI L J, et al. Flank error correction of hypoid gears based on regression analyses[J]. Journal of Tsinghua University (Science and Technology), 2017, 57(2): 141-146. (in Chinese)
[20]
CASTRO J, SEABRA J. Coefficient of friction in mixed film lubrication: Gears versus twin-discs[J]. Proceedings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2007, 221(3): 399-411.
[21]
HAMROCK B J, SCHMID S R, JACOBSON B O. Fundamentals of fluid film lubrication[M]. New York: Marcel Dekker, 2004.
[22]
XU H. Development of a generalized mechanical efficiency prediction methodology for gear pairs[D]. Columbus: The Ohio State University, 2005.
[23]
王琪, 周驰, 桂良进, 等. 准双曲面齿轮LTCA弯曲变形约束方法[J]. 航空动力学报, 2017, 32(11): 2800-2807.
WANG Q, ZHOU C, GUI L J, et al. Constraint method of calculating tooth bending deformation of hypoid gears in LTCA[J]. Journal of Aerospace Power, 2017, 32(11): 2800-2807. (in Chinese)
[24]
WILCOX L E, CHIMNER T D, NOWELL G C. Improved finite element model for calculating stresses in bevel and hypoid gear teeth[R]. Alexandria: AGMA, 1997.
[25]
HAMMAMI M, MARTINS R, ABBES M S, et al. Axle gear oils: Friction behaviour under mixed and boundary lubrication regimes[J]. Tribology International, 2017, 116: 47-57.
[26]
WANG C, WANG S R, WANG G Q. A method for calculating gear meshing efficiency by measured data from gear test machine[J]. Measurement, 2018, 119: 97-101.
[27]
Bearings SKF. The SKF Model for Calculating the Frictional Moment[J/OL]. (2012-12-15)[2023-02-07]. https://www.skf.com/binaries/pub12/Images/0901d1968065e9e7-The-SKF-model-for-calculating-the-frictional-moment_tcm_12-299767.