Mechanical Engineering and Technology 机械工程与技术, 2020, 9(3), 198-206 Published Online June 2020 in Hans. http://www.hanspub.org/journal/met https://doi.org/10.12677/met.2020.93021
Analysis and Optimization on Spiral Bevel Gear Pair Based on Finite Element Method
Jianhua Xue, Zhenhua Zhang, Chongyu Liu
Technology Center, Shaanxi Hande Axle Co., Ltd., Xi’an Shaanxi
thndth
Received: May 8, 2020; accepted: May 22, 2020; published: May 29, 2020
Abstract
The finite element method (FEM) is applied to the analysis and optimization of the fatigue capacity of spiral bevel gear pairs on the basis of Gleason rating method. The accurate finite element analysis model of spiral gear is established. The contact pattern, contact stress distribution, bending stress of tooth root, transmission error and the actual contact ratio are obtained based on FEM. The com-prehensive performance of spiral bevel gear pair is evaluated. Finally, the radius of cutter tip, the tool modification and the tooth thickness are optimized based on the finite FEM, and the experi-mental verification is carried out. The bench test results show that the results of the finite element method are basically the same as those of the bench test, and the FEM can effectively analyze and optimize the fatigue capacity of the spiral bevel gear pair.
Keywords
Spiral Bevel Gear Pair, FEM, Optimization of Gear Pair, LTCA
螺旋锥齿轮副有限元分析及优化
薛建华，张振华，刘宠誉
陕西汉德车桥有限公司技术中心，陕西 西安
收稿日期：2020年5月8日；录用日期：2020年5月22日；发布日期：2020年5月29日
摘 要
本文在格里森校核公式的基础上，将有限元方法应用于螺旋锥齿轮副的疲劳承载能力分析与优化。建立了精确的螺旋锥齿轮副有限元分析模型，得到了齿面啮合印迹、接触应力、齿根弯曲应力、传动误差以
文章引用: 薛建华, 张振华, 刘宠誉. 螺旋锥齿轮副有限元分析及优化[J]. 机械工程与技术, 2020, 9(3): 198-206. DOI: 10.12677/met.2020.93021
薛建华 等
及实际重合度，更全面地评价了螺旋锥齿轮副的综合性能。最后利用有限元方法对某齿轮副的刀尖圆角半径、刀具修形和齿厚进行了优化，并进行了试验。台架试验结果表明：有限元方法分析结果与台架试验基本相同，通过有限元方法能够有效地分析优化螺旋锥齿轮副的疲劳承载能力。
关键词
螺旋锥齿轮副，有限元方法，齿轮副优化，加载齿面接触分析
Copyright ? 2020 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
Open Access 1. 引言
螺旋锥齿轮副是传统车桥中最核心的承载零部件，目前加工技术已趋于成熟[1]。但是设计、校核和加工均比较复杂，既要考虑疲劳寿命，又要考虑NVH性能，因此试验成本较高，传统的设计校核方法无法准确全面评价加载情况下的啮合情况。
张华等[2]利用有限元方法研究了非零正传动弧齿锥齿轮强度分析。Mermoz等[3] [4]建立了螺旋锥齿轮副的有限元模型，计算获得了齿面接触应力，法向啮合力等参数，并研究了带壳体情况下齿轮副的优化问题。
2. 螺旋锥齿轮副格里森公式校核方法
目前螺旋锥齿轮副疲劳承载能力主要采用格里森公司给出的校核公式[5]，该方法经过格里森公司几十年的总结修改优化，安全性和可靠性较高。该公式校核方法清楚简单，适合前期齿坯设计，但是没有考虑实际齿面修形加工等因素，因此存在一定的局限性，一方面校核公式建立在齿面接触啮合印迹理想的状态下，而实际上啮合印迹常存在缺陷，比如大端小端脱出，或者齿顶边缘接触，该方法都无法评估；另一方面Gleason校核方法是基于齿轮节线进行的，而实际上锥齿轮运行过程中接触应力和弯曲应力在不断变化，因此取啮合过程中最大应力校核会更准确。因此格里森校核更适合于宏观参数优化，而不是适用于微观修形及NVH分析优化。
2.1. 齿根弯曲疲劳打齿
齿根弯曲疲劳打齿是售后市场最常见的失效形式，其主要是由于齿根弯曲应力过大、应力集中以及加工缺陷从而产生疲劳源，格里森给出的大轮与小轮齿根弯曲应力计算方法如式(1)~(2)。
σF2=2000T2KAKvKFβYXmb2d2J2 (1)
=σF12000T1KAKvKFβYXz22000T1maxKAKvKHβZXZRT1z23=σHZE·· (2) 2mb2d2J1′z1T1maxz1b2d2I式中：T1T2为小轮/大轮扭矩(Nm)，KA为使用系数，Kv为动载荷系数；KFβ为齿向载荷分布系数；YX为弯曲尺寸系数；m为大端模数，b2大轮齿宽，d2为大轮大端节圆直径，J1J2为小轮/大轮弯曲几何系数，
z1z2为小轮/大轮齿数。
DOI: 10.12677/met.2020.93021
199
机械工程与技术
薛建华 等
2.2. 齿面接触疲劳打齿
齿面接触疲劳失效主要包括疲劳点蚀和疲劳剥落打齿，是由齿面接触应力过大引起的。齿面接触应力如式(3)所示。
=σHZE2000T1maxKAKvKHβZXZR2b2d2I?3T1T1max?z2 (3) z1式中：ZE为弹性系数；T1max为小轮最大扭矩；KHβ为齿向载荷分配系数；ZX为接触尺寸系数；ZR为表面状况系数。
3. 螺旋锥齿轮有限元分析模型
目前通用的Gleason软件CAGE和克林贝格Kimos均具有加载接触分析LTCA (Load tooth contact analysis)功能，但是由于采用的是简化算法，分析结果与台架试验误差较大。有限元分析方法是一种各个行业广泛应用的分析手段，本文以某车桥准双曲面齿轮为例，通过有限元方法分析其啮合印迹、传动误差、齿根弯曲应力及实际重合度，为锥齿轮副的设计、校核以及NVH分析提供依据。
螺旋锥齿轮副有限元建模
根据螺旋锥齿轮加工原理，通过刀具参数和机床调整参数模拟加工过程可得到螺旋锥齿轮齿面和齿根点坐标[6]，进一步建立螺旋锥齿轮有限元模型，为了保证精度又节约时间，本文采用5齿有限元模型，齿长方向划分80个单元，齿高方向划分20个单元，齿根过渡曲线划分15个单元。在建模时，应该考虑螺旋锥齿轮副4个方向的错位量。正确装配的螺旋锥齿轮副有限元模型如图1所示，将大小轮的内孔耦合到其旋转轴线参考点上，在需要接触的齿面上建立接触对，设置初始干涉量，设置好载荷步。
Figure 1. FEM contact analysis model of spiral bevel gear pair 图1. 螺旋锥齿轮副有限元接触模型
4. 螺旋锥齿轮副有限元分析结果
对分析结果进行后处理得有效的啮合信息，作为齿轮副的设计和校核的依据。
4.1. 齿面接触应力和啮合印迹
分析中每隔0.05 s计算一个工况，某一啮合瞬间应力分布和位置如图2(a)所示，在啮合过程中，大轮齿面接触印迹呈椭圆状分布，该时刻有三对齿参与接触。
DOI: 10.12677/met.2020.93021
200
机械工程与技术
薛建华 等
Figure 2. Contact press and pattern of the of spiral bevel gear pair 图2. 螺旋锥齿轮副接触应力和啮合印迹
提取每一时刻接触应力最大值进行叠加，即可得到整个齿面的啮合印迹，如图2(b)所示，取中间齿接触应力如图3(a)所示，即为齿轮副的啮合印迹。在台架试验中，通过在齿面涂红丹粉，加载运行一段时间后查看啮合印迹，如图3(b)所示。可见有限元分析结果与台架试验结果吻合良好。有限元分析不仅能够得到啮合印迹区域和位置，而且能得到最大应力区域和大小，相较于台架试验，即节省了时间和经济成本，可在齿轮加工前对其模拟分析优化，有效缩短齿轮开发周期。
Figure 3. Contact pattern of spiral bevel gear pair 图3. 螺旋锥齿轮啮合印迹分布
4.2. 齿根弯曲应力结果
在螺旋锥齿轮副售后失效件中，由齿根弯曲疲劳引起的打齿失效占很大比例。有限元齿根弯曲应力结果如图4(a)所示，典型的螺旋锥齿轮齿根弯曲疲劳打齿失效件如图4(b)所示，售后旧件打齿疲劳裂源与有限元分析结果中最大应力位置弯曲吻合，均处在工作侧齿根中点略偏小端齿根圆角处。
4.3. 螺旋锥齿轮副传动误差分布
螺旋锥齿轮采用局部共轭的设计方法，传动误差TE (transmission error)是螺旋锥齿轮副设计中关键因素之一，传动误差幅值越小，齿轮传动越平稳，振动和噪音越小，但是传动误差太小则可能发生边缘接触，引起齿面接触疲劳打齿。传动误差计算如式(4)所示：
TE=φ2?φ20?()z1φ1?φ10 (4) z2()式中，φ20φ10为小轮/大轮在初始转角，φ2φ1为某一啮合瞬间的角位移，z1z2为小轮/大轮齿数。
DOI: 10.12677/met.2020.93021
201
机械工程与技术
薛建华 等
Figure 4. The FEM root bending stress and the failure part of spiral bevel gear 图4. 螺旋锥齿轮有限元齿根弯曲应力和失效件
分析过程中提取每一接触时刻的大小轮角位移，按照式(4)即可得到该工况下的齿轮副传动误差分布。有限元方法不仅能够计算空载下传动误差，也可以得到加载情况下的传动误差。有限元分析方法和传统方法计算传动误差如图5所示，两者分布趋势和幅值基本相同，说明有限元方法能够准确地计算螺旋齿轮副的传动误差，有限元相比数值计算方法优势在于可以模拟加载情况下轮齿和辐板的各种变形，更符合实际情况。
Figure 5. The comparison of transmission error between FEM and numerical calculation method 图5. 螺旋锥齿轮有限元传动误差与数值计算方法对比
4.4. 螺旋锥齿轮副法向啮合力分布和实际重合度
理论研究和实际应用均表明，重合度是影响齿轮副承载能力、振动和噪音的重要因素之一，实际重合度越大，传动越平稳，振动和噪音越小。目前CAGE和Kimos软件只能计算空载时的重合度，无法获得加载情况下实际重合度。在有限元分析中，法向啮合力反映了轮齿的接触状态，通过法向啮合力分布可以得到加载情况下的实际重合度，如图6所示。
在该情况下，重合度εγ=?tt0，t0为齿轮副啮合周期，?t为一个轮齿从啮入到啮出经历的时间，据此得到该载荷下总重合度为2.2。
DOI: 10.12677/met.2020.93021
202
机械工程与技术
螺旋锥齿轮副有限元分析及优化 - 图文
