CN111488658A - 斜齿轮副齿宽修形优化设计方法 - Google Patents

斜齿轮副齿宽修形优化设计方法 Download PDF

Info

Publication number
CN111488658A
CN111488658A CN202010274023.7A CN202010274023A CN111488658A CN 111488658 A CN111488658 A CN 111488658A CN 202010274023 A CN202010274023 A CN 202010274023A CN 111488658 A CN111488658 A CN 111488658A
Authority
CN
China
Prior art keywords
gear pair
tooth width
tooth
modification
bevel gear
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010274023.7A
Other languages
English (en)
Inventor
刘辉
严鹏飞
张晨
吴云豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202010274023.7A priority Critical patent/CN111488658A/zh
Publication of CN111488658A publication Critical patent/CN111488658A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Gears, Cams (AREA)

Abstract

本发明提供了一种斜齿轮副齿宽修形优化设计方法,包括以下步骤:步骤A:确定斜齿轮副的设计变量;步骤B:以斜齿轮副的齿向动态载荷分布系数作为齿宽方向动态性能优化目标,并确定目标函数;步骤C:建立设计变量的边界约束;步骤D:根据设计变量、目标函数和边界条件在matlab中建立基于齿宽修形的动态性能优化数学模型,并对数学模型进行求解。本发明通过软件对齿宽进行修形,实现了斜齿轮副齿宽修形优化,提高了齿面动态载荷分布均匀度,降低了齿面动态载荷峰值。

Description

斜齿轮副齿宽修形优化设计方法
技术领域
本发明涉及机械动力学技术领域,特别是涉及一种斜齿轮副齿宽修形优化设计方法。
背景技术
斜齿轮副在传动过程中,常常面临着齿面动态载荷分布不均,且齿面动态载荷峰值较高的问题。
因此,如何实现斜齿轮副齿宽修形优化,提高齿面动态载荷分布均匀度,降低齿面动态载荷峰值为本领域技术人员亟待解决的技术问题。
发明内容
有鉴于此,本发明目的是提供一种斜齿轮副齿宽修形优化设计方法,能够实现斜齿轮副齿宽修形优化,提高齿面动态载荷分布均匀度,降低齿面动态载荷峰值。
为了达到上述目的,本发明提供如下技术方案:
一种斜齿轮副齿宽修形优化设计方法,包括以下步骤:
步骤A:确定斜齿轮副的设计变量;
步骤B:以斜齿轮副的齿向动态载荷分布系数作为齿宽方向动态性能优化目标,并确定目标函数;
步骤C:建立所述设计变量的边界约束;
步骤D:根据所述设计变量、所述目标函数和所述边界条件在matlab中建立基于齿宽修形的动态性能优化数学模型,并对所述数学模型进行求解。
在一个具体实施方案中,所述步骤A中,所述斜齿轮副中的两个齿轮的设计变量相同,所述设计变量包括:所述齿轮副中齿端两侧位置的最大修形量Δd,max1和Δd,max2,述齿轮副中齿端两侧多项式函数的弯曲指数s1和s2,所述齿轮副中齿端两侧位置的修形长度b1和b2
在另一个具体实施方案中,所述步骤B具体包括:
步骤B1:将所述齿轮副的齿轮沿齿宽方向分解为N个独立且均匀的薄片齿轮副;
步骤B2:确定目标函数fd(X)=Kβ
Figure BDA0002444148290000021
式中,(Fm,i)RMS为第i对所述薄片齿轮副动态啮合力的均方根值,[(Fm,i)RMS]max为每对所述薄片齿轮副动态啮合力的均方根值的最大值,Fn为准静态条件下齿面所受实际载荷,N为切片后处于啮合状态的所述薄片齿轮副总数。
在另一个具体实施方案中,所述步骤C中的约束条件包括边界约束和性能约束。
在另一个具体实施方案中,所述边界约束为:
Figure BDA0002444148290000022
Figure BDA0002444148290000023
为所述设计变量下限组成的列向量,
Figure BDA0002444148290000024
为所述设计变量上限组成的列向量。
在另一个具体实施方案中,所述性能约束包括齿轮弯曲疲劳强度和接触疲劳强度约束;
σH≤[σH];
σF≤[σF];
σF和σH分别为弯曲疲劳强度和接触疲劳强度,[σF]和[σH]分别为弯曲疲劳强度极限和接触疲劳强度极限。
在另一个具体实施方案中,所述最大修形量的上限取为平均转矩作用下齿面最大弹性变形量,齿宽修形长度范围为0≤bi≤d/2;
d为齿宽。
在另一个具体实施方案中,所述修形曲线的多项式函数的弯曲指数范围为1≤si≤5;
所述最大修形量的上限为33μm。
在另一个具体实施方案中,所述步骤D中的数学模型为:
min fd(X)。
根据本发明的各个实施方案可以根据需要任意组合,这些组合之后所得的实施方案也在本发明范围内,是本发明具体实施方式的一部分。
根据上述技术方案可知,本发明提供的斜齿轮副齿宽修形优化设计方法,首先,确定斜齿轮副的设计变量;接着,确定目标函数;然后,建立设计变量的边界约束;最后,在matlab中建立齿宽修形的动态性能优化数学模型,并对数学模型进行求解。本发明通过软件对齿宽进行修形,实现了斜齿轮副齿宽修形优化,提高了齿面动态载荷分布均匀度,降低了齿面动态载荷峰值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出新颖性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1是本发明提供的斜齿轮副齿宽修形优化设计方法流程图;
图2是本发明实施例二中的载荷分配因子对比图;
图3是本发明实施例二中的修形后的接触动应力分布三维图;
图4是本发明实施例二中的修形后的应力分布二维等高线图;
图5是本发明实施例二中的无修形的接触动应力分布三维图;
图6是本发明实施例二中的无修形的应力分布二维等高线图。
具体实施方式
为了使本领域的技术人员更好的理解本发明的技术方案,下面结合附图和具体实施方式对本发明作进一步的详细说明。
如图1所示,本发明公开了一种斜齿轮副齿宽修形优化设计方法,包括以下步骤:
步骤S1:确定斜齿轮副的设计变量。
具体地,斜齿轮副中的两个齿轮的设计变量相同,以任意一个齿轮的设计变量为例,设计变量包括:齿轮副中齿端两侧位置的最大修形量Δd,max1和Δd,max2,述齿轮副中齿端两侧多项式函数的弯曲指数s1和s2,齿轮副中齿端两侧位置的修形长度b1和b2
由于实际工作状态下,齿面中心两侧载荷分布有较大区别,因此,将齿轮副中齿端两侧位置的最大修形量Δd,max1和Δd,max2作为设计变量。
步骤S2:以斜齿轮副的齿向动态载荷分布系数作为齿宽方向动态性能优化目标,并确定目标函数。
具体地,目标函数为:fd(X)=Kβ。齿向动态载荷分布系数是与各薄片齿轮副齿面所受动态载荷均值及最大载荷相关参数,可以用来表达斜齿轮副动态载荷在齿宽方向分布的不均匀程度。以齿轮传动系统的齿向动态载荷分布系数作为齿宽方向动态性能优化目标。
步骤S2具体包括步骤S21和步骤S22。
步骤S21:将齿轮副的齿轮沿齿宽方向分解为N个独立且均匀的薄片齿轮副。
步骤S22:确定目标函数fd(X)=Kβ
Figure BDA0002444148290000051
式中,(Fm,i)RMS为第i对薄片齿轮副动态啮合力的均方根值,[(Fm,i)RMS]max为每对薄片齿轮副动态啮合力的均方根值的最大值,Fn为准静态条件下齿面所受实际载荷,N为切片后处于啮合状态的薄片齿轮副总数。
步骤S3:建立设计变量的边界约束。
约束条件包括边界约束和性能约束。
进一步地,本发明公开了边界约束为设计变量的边界条件:
Figure BDA0002444148290000052
其中,
Figure BDA0002444148290000053
为设计变量下限组成的列向量,
Figure BDA0002444148290000054
为设计变量上限组成的列向量。
进一步地,本发明公开了性能约束包括齿轮弯曲疲劳强度和接触疲劳强度约束。σH≤[σH],σF≤[σF],其中,σF和σH分别为弯曲疲劳强度和接触疲劳强度,[σF]和[σH]分别为弯曲疲劳强度极限和接触疲劳强度极限。
进一步地,本发明公开了最大修形量的上限取为平均转矩作用下齿面最大弹性变形量,齿宽修形长度范围为0≤bi≤d/2,d为齿宽。
具体地,本发明公开了修形曲线的多项式函数的弯曲指数范围为1≤si≤5,最大修形量的上限为33μm。
步骤S4:根据设计变量、目标函数和边界条件在matlab中建立基于齿宽修形的动态性能优化数学模型,并对数学模型进行求解。
具体地,数学模型为:min fd(X)。
本发明提供的斜齿轮副齿宽修形优化设计方法,首先,确定斜齿轮副的设计变量;接着,确定目标函数;然后,建立设计变量的边界约束;最后,在matlab中建立齿宽修形的动态性能优化数学模型,并对数学模型进行求解。本发明通过软件对齿宽进行修形,实现了斜齿轮副齿宽修形优化,提高了齿面动态载荷分布均匀度,降低了齿面动态载荷峰值。
实施例二
本实施例中,以采用多参数修形方法对斜齿轮副进行齿宽动态特性优化设计。
斜齿轮副动态性能优化前后的设计变量和目标函数分别如表1和表2所示。从结果可以看出,采用齿宽修形优化设计方法,优化后齿轮的齿宽动态载荷系数为1.58,相比初始设计变量下的动载系数2.03降低了22.17%。优化设计后使得齿宽动态载荷水平有了显著降低。
表1
Figure BDA0002444148290000071
表2
Figure BDA0002444148290000072
仿真分析了最优修形参数下斜齿轮副的动态载荷分配因子在齿宽方向的分布状态,并与无修形的标准齿轮副对比分析,结果如图2所示。从图中可以看出,相比于无修形状态,在最优齿宽修形参数下齿轮系统的动载荷在齿宽方向的变化幅值显著降低了,载荷分配因子最大值从1.39降低至1.11。此外,还可看出载荷分布主要集中在齿面的中间区域,且载荷中心位置的分布较为均匀。
仿真分析了最优修形参数下齿轮系统的齿面接触动应力,并与无修形状态下齿轮副接触动应力进行对比分析,结果如图3-6所示。从图3所示的齿面接触应力三维图可以看出,优化后的齿面应力分布主要集中在接触重合度较低的区域,且位于齿宽的中心区域。接触动应力最大峰值约为343MPa,相比于图5所示无修形状态下最大动应力419MPa,降低了约18.14%。此外,从图4所示的齿面接触应力二维等高线图可以看出,在最优修形参数下的齿轮副齿面接触动应力在超过340MPa的齿面接触位置主要集中在齿面宽度为[-17.5mm,10.5mm]的范围内,相比于图6所示无修形状态下齿面接触应力二维等高线图,有效的缓解了齿面应力中心偏离齿面中心位置和载荷中心位置应力分布不均匀的现象。
从结果分析可知,采用多参数动态性能的齿宽修形优化设计方法可以有效地降低齿面动态载荷峰值,缓解齿面动态载荷分布不均现象,并改善齿面载荷偏离齿面中心位置的状态,为高速重载、高精度的齿轮设计提供理论技术依据。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的实施例,而是要符合与本文所公开的原理和创造特点相一致的最宽的范围。

Claims (9)

1.一种斜齿轮副齿宽修形优化设计方法,其特征在于,包括以下步骤:
步骤A:确定斜齿轮副的设计变量;
步骤B:以斜齿轮副的齿向动态载荷分布系数作为齿宽方向动态性能优化目标,并确定目标函数;
步骤C:建立所述设计变量的边界约束;
步骤D:根据所述设计变量、所述目标函数和所述边界条件在matlab中建立基于齿宽修形的动态性能优化数学模型,并对所述数学模型进行求解。
2.根据权利要求1所述的斜齿轮副齿宽修形优化设计方法,其特征在于,所述步骤A中,所述斜齿轮副中的两个齿轮的设计变量相同,所述设计变量均包括:所述齿轮副中齿端两侧位置的最大修形量Δd,max1和Δd,max2,述齿轮副中齿端两侧多项式函数的弯曲指数s1和s2,所述齿轮副中齿端两侧位置的修形长度b1和b2
3.根据权利要求1所述的斜齿轮副齿宽修形优化设计方法,其特征在于,所述步骤B具体包括:
步骤B1:将所述齿轮副的齿轮沿齿宽方向分解为N个独立且均匀的薄片齿轮副;
步骤B2:确定目标函数fd(X)=Kβ
Figure FDA0002444148280000011
式中,(Fm,i)RMS为第i对所述薄片齿轮副动态啮合力的均方根值,[(Fm,i)RMS]max为每对所述薄片齿轮副动态啮合力的均方根值的最大值,Fn为准静态条件下齿面所受实际载荷,N为切片后处于啮合状态的所述薄片齿轮副总数。
4.根据权利要求1所述的斜齿轮副齿宽修形优化设计方法,其特征在于,所述步骤C中的约束条件包括边界约束和性能约束。
5.根据权利要求4所述的斜齿轮副齿宽修形优化设计方法,其特征在于,所述边界约束为:
Figure FDA0002444148280000021
Figure FDA0002444148280000022
为所述设计变量下限组成的列向量,
Figure FDA0002444148280000023
为所述设计变量上限组成的列向量。
6.根据权利要求5所述的斜齿轮副齿宽修形优化设计方法,其特征在于,所述性能约束包括齿轮弯曲疲劳强度和接触疲劳强度约束;
σH≤[σH];
σF≤[σF];
σF和σH分别为弯曲疲劳强度和接触疲劳强度,[σF]和[σH]分别为弯曲疲劳强度极限和接触疲劳强度极限。
7.根据权利要求2所述的斜齿轮副齿宽修形优化设计方法,其特征在于,所述最大修形量的上限取为平均转矩作用下齿面最大弹性变形量,齿宽修形长度范围为0≤bi≤d/2;
d为齿宽。
8.根据权利要求2-7中任意一项所述的斜齿轮副齿宽修形优化设计方法,其特征在于,所述修形曲线的多项式函数的弯曲指数范围为1≤si≤5;
所述最大修形量的上限为33μm。
9.根据权利要求3所述的斜齿轮副齿宽修形优化设计方法,其特征在于,所述步骤D中的数学模型为:
min fd(X)。
CN202010274023.7A 2020-04-09 2020-04-09 斜齿轮副齿宽修形优化设计方法 Pending CN111488658A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010274023.7A CN111488658A (zh) 2020-04-09 2020-04-09 斜齿轮副齿宽修形优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010274023.7A CN111488658A (zh) 2020-04-09 2020-04-09 斜齿轮副齿宽修形优化设计方法

Publications (1)

Publication Number Publication Date
CN111488658A true CN111488658A (zh) 2020-08-04

Family

ID=71792677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010274023.7A Pending CN111488658A (zh) 2020-04-09 2020-04-09 斜齿轮副齿宽修形优化设计方法

Country Status (1)

Country Link
CN (1) CN111488658A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112989514A (zh) * 2021-03-05 2021-06-18 明阳智慧能源集团股份公司 一种风力发电机传动齿轮的不均载系数校核工具

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6080199A (en) * 1996-10-02 2000-06-27 Toyota Jidosha Kabushiki Kaisha Gear designing method, gear manufacturing method, and gear manufactured by the same method
US20180003283A1 (en) * 2016-05-31 2018-01-04 Barry James Planetary Gearset
CN108856910A (zh) * 2018-08-21 2018-11-23 中车大连机车车辆有限公司 一种齿轮修缘的方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6080199A (en) * 1996-10-02 2000-06-27 Toyota Jidosha Kabushiki Kaisha Gear designing method, gear manufacturing method, and gear manufactured by the same method
US20180003283A1 (en) * 2016-05-31 2018-01-04 Barry James Planetary Gearset
CN108856910A (zh) * 2018-08-21 2018-11-23 中车大连机车车辆有限公司 一种齿轮修缘的方法及装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
杨硕文 等: "基于KISSsoft软件的人字齿轮传动修形设计研究", 《机械传动》 *
王成: "基于齿轮传动系统横-扭-摆耦合", 《振动与冲击》 *
王春燕: "基于KISSsoft斜齿轮齿面的优化修形", 《机械传动》 *
王海 等: "数字样机技术在升船机主提升减速器设计中的应用", 《中国重型装备》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112989514A (zh) * 2021-03-05 2021-06-18 明阳智慧能源集团股份公司 一种风力发电机传动齿轮的不均载系数校核工具
CN112989514B (zh) * 2021-03-05 2023-03-14 明阳智慧能源集团股份公司 一种风力发电机传动齿轮的不均载系数校核工具

Similar Documents

Publication Publication Date Title
CN107256300B (zh) 基于齿面应力边棱作用和齿轮歪斜变形的圆柱直齿轮齿向修形方法
EP1803974B1 (en) Gear wheel with chamfered portions
CN109630652B (zh) 一种三圆弧谐波齿轮插齿刀及其齿廓设计方法
CN111488682B (zh) 渐开线斜齿轮副齿宽修形动力学模型建立方法
JP5894173B2 (ja) かさ歯車伝動装置のかさ歯車における最適化されたクラウニング
CN104889505B (zh) 非对称滚刀及其设计方法
US10584783B2 (en) High-strength gear
CN111488658A (zh) 斜齿轮副齿宽修形优化设计方法
US6230578B1 (en) Variable face width gearing
CN106369139A (zh) 一种满足高次传动误差的准双曲面齿轮加工参数获取方法
Jia et al. Design and analysis of double-crowned high-contact-ratio cylindrical gears considering the load sharing of the multi-pair contact
CN110188508A (zh) 基于实际载荷的齿面接触斑点齿轮的设计方法
CN111488660B (zh) 渐开线斜齿轮螺旋角修形优化设计方法
CN111219473B (zh) 一种大重合度弧齿锥齿轮设计方法
CN110083939B (zh) 纯电动车差速器齿轮修形优化方法
US20020170373A1 (en) Variable face width gearing
CN115809527A (zh) 基于切片法的斜齿轮时变啮合刚度及传递误差计算方法
CN111079233B (zh) 一种弹性变形优化剃齿刀设计方法
CN114769741B (zh) 一种齿顶全圆弧滚刀及其设计方法
CN109931367B (zh) 一种直线型伞齿轮的啮合结构、减速机及该齿轮的制造方法
CN111120624A (zh) 一种基于接触应力均化的摆线轮齿廓修形方法
Walton et al. A note on tip relief and backlash allowances in non-metallic gears
US20240167556A1 (en) Gear pair
CN113059238A (zh) 一种降低轮边桥等高螺旋锥齿轮传动误差的方法
Kumar et al. Material optimization of spur gear tooth

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200804