CN105224744A - 一种基于啮合刚度的剥落齿轮啮合模型 - Google Patents
一种基于啮合刚度的剥落齿轮啮合模型 Download PDFInfo
- Publication number
- CN105224744A CN105224744A CN201510634304.8A CN201510634304A CN105224744A CN 105224744 A CN105224744 A CN 105224744A CN 201510634304 A CN201510634304 A CN 201510634304A CN 105224744 A CN105224744 A CN 105224744A
- Authority
- CN
- China
- Prior art keywords
- integral
- engagement
- rsqb
- lsqb
- 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.)
- Granted
Links
Landscapes
- Gears, Cams (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
一种基于啮合刚度的剥落齿轮啮合模型,考虑剥落轮齿在齿轮啮合时的受力情况,将啮合过程分为了三个阶段,并分别计算了三个啮合阶段下的剥落轮齿啮合刚度,克服了现有剥落齿轮模型仅考虑剥落啮合区与非剥落啮合区两种情况的弊端;同时考虑轮齿齿根的影响,将轮齿等效为齿根圆上的悬臂梁,准确计算出了剥落轮齿啮合刚度,提高了其精度,建立了更为准确的剥落齿轮啮合模型,本发明真实准确地反映了剥落齿轮啮合过程,能有效地应用于剥落齿轮的动力学建模研究。
Description
技术领域
本发明涉及齿轮***动力学建模领域,具体涉及一种基于啮合刚度的剥落齿轮啮合模型。
背景技术
齿轮***是现代机械工业中的常用机构,建立齿轮***动力学模型,研究齿轮***的振动特性非常必要。由于齿轮时变啮合刚度是齿轮振动的主要激励源,因此准确计算齿轮时变啮合刚度是齿轮动力学建模中的关键技术。
剥落故障是齿轮运行中的常见故障,准确计算剥落故障下的齿轮啮合刚度并建立剥落齿轮啮合模型是研究剥落齿轮振动响应的重要基础。目前大多数剥落齿轮啮合模型中使用的齿轮啮合刚度都把轮齿简化为基圆上的悬臂梁,并未考虑轮齿齿根对啮合刚度的影响。并且现有的啮合模型都只将剥落轮齿的计算简单地分为剥落区域啮合与非剥落区域啮合两种情况,并未考虑剥落齿轮的实际啮合过程。
发明内容
为了克服上述现有技术存在的缺点,本发明的目的在于提供一种基于啮合刚度的剥落齿轮啮合模型,提高啮合刚度精度,建立准确的剥落齿轮啮合模型。
为了达到上述目的,本发明采取的技术方案为:
一种基于啮合刚度的剥落齿轮啮合模型,包括以下步骤:
1)假设剥落出现在齿轮的节线位置,剥落区域为一个B×S的矩形区域,剥落厚度为t;
2)将承受啮合力的剥落轮齿简化为齿根圆上的悬臂梁;
3)根据剥落轮齿的啮合过程,将其分为三个啮合阶段,啮合阶段1为齿根啮合区,啮合阶段2为剥落区域啮合区,啮合阶段3为齿尖啮合区;
4)在啮合阶段1时,剥落轮齿的齿根部位受力,轮齿剥落区域对齿轮啮合没有影响,分别计算啮合阶段1时的弯曲势能Ub_1、剪切势能Us_1、轴向压力势能Ua_1及赫兹刚度kh_1:
其中,F为轮齿间啮合力,可以分解为径向力Fa与切向力Fb。d为啮合点到基圆的距离,x为轮齿上任意点到基圆的距离,x1为齿根上任意点到齿根圆的距离,h为啮合点到轮齿中线的距离,E为弹性模量,G为剪切模量,Rb为基圆半径,Rr为齿根圆半径,Ix_1为轮齿上任意点处的转动惯量,Ix为齿根上任意点处的转动惯量,Ax_1为轮齿上任意点处的横截面积,Ax为齿根上任意点处的横截面积,L为齿宽,υ为泊松比,
求得势能后,即求得相应的弯曲刚度kb、剪切刚度ks、轴向压缩刚度ka:
5)经过啮合阶段1后,剥落故障齿轮进一步转动,齿轮剥落故障部分进入啮合,此时为啮合阶段2,齿轮剥落区域进入啮合时,该区域并不受力,此时的接触线长度由齿宽L变为L-B,啮合阶段2时的弯曲势能Ub_2、剪切势能Us_2、轴向压缩势能Ua_2及赫兹刚度kh_2为:
6)经过啮合阶段2后,齿轮剥落故障部分脱出啮合进入啮合阶段3,此时齿轮剥落故障区域虽然不直接参与啮合,但剥落区域属于悬臂梁的一部分,因而在啮合阶段3时齿轮剥落故障区域同样会影响时变啮合刚度,在计算时应减去剥落区域的刚度,啮合阶段3时的弯曲势能Ub_3、剪切势能Us_3、轴向压缩势能Ua_3及赫兹刚度kh_3为:
其中,ds1为剥落起始位置到基圆的距离,ds2为剥落终止位置到基圆的距离,Isx_3为剥落故障区域某点的转动惯量,Asx_3为剥落故障区域某点的横截面积;
7)计算完各个啮合阶段时的刚度后,需要知道各个啮合阶段的啮合时间,假设剥落轮齿刚进入双齿啮合时为初始时刻,剥落轮齿对应角度为起始角度,此时θ=0;齿轮继续转动,当转动到剥落故障区域刚进入啮合时,此时θ=αs1;当剥落故障区域刚退出啮合时,此时θ=αs2,因此,0<θ<αs1时为剥落轮齿啮合阶段1,αs1<θ<αs2时为啮合阶段2,θ>αs2时为啮合阶段3,在不同的啮合阶段应选用不同的公式计算啮合刚度,最后即得到齿轮一个周期内的综合啮合刚度:
式中j=1表示第一对轮齿啮合,j=n表示第n对轮齿啮合,
8)基于步骤7)中得到的剥落齿轮综合啮合刚度,按照齿轮动力学建模方法,即可得到基于啮合刚度的剥落齿轮啮合模型。
本发明充分考虑剥落轮齿在齿轮啮合时的受力情况,将啮合过程分为了三个阶段,并分别计算了三个啮合阶段下的剥落轮齿啮合刚度,克服了现有剥落齿轮啮合模型仅考虑剥落啮合区与非剥落啮合区两种情况的弊端;同时考虑轮齿齿根的影响,将轮齿等效为齿根圆上的悬臂梁,准确计算出了剥落轮齿啮合刚度,提高了其精度,建立了更为准确的剥落齿轮啮合模型。
附图说明
图1为轮齿剥落示意图。
图2为剥落轮齿受力图。
图3为剥落齿轮啮合阶段1示意图。
图4为剥落齿轮啮合阶段2示意图。
图5为剥落齿轮啮合阶段3示意图。
图6为不同剥落宽度下的齿轮啮合刚度曲线。
图7为不同剥落长度下的齿轮啮合刚度曲线。
图8为不同剥落深度下的齿轮啮合刚度曲线。
具体实施方式
下面结合附图对本发明做进一步详细描述:
一种基于啮合刚度的剥落齿轮啮合模型,包括以下步骤:
1)由于轮齿剥落一般出现在节线附近靠近齿根部分的表面上,因此假设剥落出现在齿轮的节线位置,剥落区域为一个B×S的矩形区域,剥落厚度为t,剥落示意图如图1所示;
2)考虑齿根对啮合刚度的影响,将承受啮合力的剥落轮齿简化为齿根圆上的悬臂梁,对剥落轮齿进行受力分析,受力图如图2所示,由图2可知,将基圆与齿根圆之间的过渡圆弧作为悬臂梁的一部分,从而考虑了齿根对啮合刚度的影响;
3)根据剥落轮齿的啮合过程,将其分为三个啮合阶段,啮合阶段1为齿根啮合区,啮合阶段2为剥落区域啮合区,啮合阶段3为齿尖啮合区;
4)在啮合阶段1时,剥落轮齿的齿根部位受力,轮齿剥落区域对齿轮啮合没有影响,受力示意图如图3所示,左图中,齿轮O1为剥落齿轮,为主动轮;齿轮O2为正常齿轮,为被动轮,右图为从A视角观察剥落齿面看到的视图,图中阴影部分为剥落区域,AB为啮合接触线长度,当剥落齿轮刚进入啮合,即d<ds1时,齿轮中的剥落部位没有受到啮合力F,
啮合阶段1时的弯曲势能Ub_1、剪切势能Us_1、轴向压力势能Ua_1及赫兹刚度kh_1:
其中,F为轮齿间啮合力,可以分解为径向力Fa与切向力Fb。d为啮合点到基圆的距离,x为轮齿上任意点到基圆的距离,x1为齿根上任意点到齿根圆的距离,h为啮合点到轮齿中线的距离,E为弹性模量,G为剪切模量,Rb为基圆半径,Rr为齿根圆半径,Ix_1为轮齿上任意点处的转动惯量,Ix为齿根上任意点处的转动惯量,Ax_1为轮齿上任意点处的横截面积,Ax为齿根上任意点处的横截面积,L为齿宽,υ为泊松比,
求得势能后,即求得相应的弯曲刚度kb、剪切刚度ks、轴向压缩刚度ka:
5)经过啮合阶段1后,剥落故障齿轮进一步转动,齿轮剥落故障部分进入啮合,此时为啮合阶段2,其受力示意图如图4,从图4可以看出,此时啮合力F的作用点在剥落区域,即ds1<d<ds2,此时剥落区域并不受力,接触线长度由齿宽L变为L-B,啮合阶段2时的弯曲势能Ub_2、剪切势能Us_2、轴向压缩势能Ua_2及赫兹刚度kh_2为:
6)经过啮合阶段2后,剥落齿轮进入啮合阶段3后,其受力示意图如图5,图中可以看出此时齿轮剥落故障区域虽然不直接参与啮合,合力F的作用点在齿尖区域,即d>ds2,但剥落区域属于悬臂梁的一部分,因而在啮合阶段3时齿轮剥落故障区域同样会影响时变啮合刚度,在计算时应减去剥落区域的刚度。啮合阶段3时的弯曲势能Ub_3、剪切势能Us_3、轴向压缩势能Ua_3及赫兹刚度kh_3为:
其中,ds1为剥落起始位置到基圆的距离,ds2为剥落终止位置到基圆的距离,Isx_3为剥落故障区域某点的转动惯量,Asx_3为剥落故障区域某点的横截面积;
7)计算完各个啮合阶段时的刚度后,需要知道各个啮合阶段的啮合时间,假设剥落轮齿刚进入双齿啮合时为初始时刻,剥落轮齿对应角度为起始角度,此时θ=0;齿轮继续转动,当转动到剥落故障区域刚进入啮合时,此时θ=αs1;当剥落故障区域刚退出啮合时,此时θ=αs2。αs1与αs2可以根据图2所示的几何关系求得,因此,0<θ<αs1时为剥落轮齿啮合阶段1,αs1<θ<αs2时为啮合阶段2,θ>αs2时为啮合阶段3,在不同的啮合阶段应选用不同的公式计算啮合刚度,最后即可得到齿轮一个周期内的综合啮合刚度:
式中j=1表示第一对轮齿啮合,j=n(n=1~2)表示第n对轮齿啮合,n=1时表示单齿啮合,n=2时表示双齿啮合,
8)基于步骤7)中得到的剥落齿轮综合啮合刚度,按照齿轮动力学建模方法,即可建立基于啮合刚度的剥落齿轮啮合模型。
将上述发明应用于下表所示的齿轮副中,其中主动轮为剥落齿轮。不同剥落尺寸下的齿轮啮合刚度分别如图6~图8所示。
由图6~图8可以看出,当轮齿剥落区域进入与离开啮合瞬间,齿轮的啮合刚度都将产生突变,从而产生冲击激励。图6表明,剥落长度的增加会导致剥落起始角度的减小与剥落终止角度的增大,从而增大剥落区域作用的时间。从图中可以看出,剥落长度越大,剥落区域对应的作用时间越长。图7表明,剥落宽度沿齿宽方向增大时,虽然剥落区域进入和离开啮合的时间保持不变,但是由于啮合线长度发生变化,啮合刚度会随着剥落宽度的增加而降低。由图8可看出,剥落深度增大对剥落区域啮合刚度无影响,而影响区域后的啮合刚度,剥落深度越大,相应啮合刚度越小,但总体来说剥落深度影响不大。
以上内容是结合具体实例对本发明所作的进一步详细说明,不能认定本发明的具体实施方式仅限于此,还可用于不同的齿轮***模型的建立。实施者只需对本模型相应参数进行适当调整,以适应不同产品的应用需求。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单的推演或替换,都应当视为属于本发明由所提交的权利要求书确定专利保护范围。
Claims (1)
1.一种基于啮合刚度的剥落齿轮啮合模型,其特征在于,包括以下步骤:
1)假设剥落出现在齿轮的节线位置,剥落区域为一个B×S的矩形区域,剥落厚度为t;
2)将承受啮合力的剥落轮齿简化为齿根圆上的悬臂梁;
3)根据剥落轮齿的啮合过程,将其分为三个啮合阶段,啮合阶段1为齿根啮合区,啮合阶段2为剥落区域啮合区,啮合阶段3为齿尖啮合区;
4)在啮合阶段1时,剥落轮齿的齿根部位受力,轮齿剥落区域对齿轮啮合没有影响,分别计算啮合阶段1时的弯曲势能Ub_1、剪切势能Us_1、轴向压力势能Ua_1及赫兹刚度kh_1:
其中,F为轮齿间啮合力,可以分解为径向力Fa与切向力Fb。d为啮合点到基圆的距离,x为轮齿上任意点到基圆的距离,x1为齿根上任意点到齿根圆的距离,h为啮合点到轮齿中线的距离,E为弹性模量,G为剪切模量,Rb为基圆半径,Rr为齿根圆半径,Ix_1为轮齿上任意点处的转动惯量,Ix为齿根上任意点处的转动惯量,Ax_1为轮齿上任意点处的横截面积,Ax为齿根上任意点处的横截面积,L为齿宽,υ为泊松比,
求得势能后,即求得相应的弯曲刚度kb、剪切刚度ks、轴向压缩刚度ka:
5)经过啮合阶段1后,剥落故障齿轮进一步转动,齿轮剥落故障部分进入啮合,此时为啮合阶段2,齿轮剥落区域进入啮合时,该区域并不受力,此时的接触线长度由齿宽L变为L-B,啮合阶段2时的弯曲势能Ub_2、剪切势能Us_2、轴向压缩势能Ua_2及赫兹刚度kh_2为:
6)经过啮合阶段2后,齿轮剥落故障部分脱出啮合进入啮合阶段3,此时齿轮剥落故障区域虽然不直接参与啮合,但剥落区域属于悬臂梁的一部分,因而在啮合阶段3时齿轮剥落故障区域同样会影响时变啮合刚度,在计算时应减去剥落区域的刚度,啮合阶段3时的弯曲势能Ub_3、剪切势能Us_3、轴向压缩势能Ua_3及赫兹刚度kh_3为:
其中,ds1为剥落起始位置到基圆的距离,ds2为剥落终止位置到基圆的距离,Isx_3为剥落故障区域某点的转动惯量,Asx_3为剥落故障区域某点的横截面积;
7)计算完各个啮合阶段时的刚度后,需要知道各个啮合阶段的啮合时间,假设剥落轮齿刚进入双齿啮合时为初始时刻,剥落轮齿对应角度为起始角度,此时θ=0;齿轮继续转动,当转动到剥落故障区域刚进入啮合时,此时θ=αs1;当剥落故障区域刚退出啮合时,此时θ=αs2,因此,0<θ<αs1时为剥落轮齿啮合阶段1,αs1<θ<αs2时为啮合阶段2,θ>αs2时为啮合阶段3,在不同的啮合阶段应选用不同的公式计算啮合刚度,最后即得到齿轮一个周期内的综合啮合刚度:
式中j=1表示第一对轮齿啮合,j=n表示第n对轮齿啮合,
8)基于步骤7)中得到的剥落齿轮综合啮合刚度,按照齿轮动力学建模方法,即可得到基于啮合刚度的剥落齿轮啮合模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510634304.8A CN105224744B (zh) | 2015-09-29 | 2015-09-29 | 一种基于啮合刚度的剥落齿轮啮合建模的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510634304.8A CN105224744B (zh) | 2015-09-29 | 2015-09-29 | 一种基于啮合刚度的剥落齿轮啮合建模的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105224744A true CN105224744A (zh) | 2016-01-06 |
CN105224744B CN105224744B (zh) | 2018-03-16 |
Family
ID=54993710
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510634304.8A Active CN105224744B (zh) | 2015-09-29 | 2015-09-29 | 一种基于啮合刚度的剥落齿轮啮合建模的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105224744B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844818A (zh) * | 2016-11-09 | 2017-06-13 | 北京工业大学 | 基于粗糙表面的直齿轮三维接触刚度计算方法 |
CN107314845A (zh) * | 2017-07-03 | 2017-11-03 | 重庆理工大学 | 齿轮动态啮合力的测试方法 |
CN107729626A (zh) * | 2017-09-26 | 2018-02-23 | 西安交通大学 | 一种基于概率分布的齿轮点蚀模型建模方法 |
CN107798200A (zh) * | 2017-11-10 | 2018-03-13 | 西安电子科技大学 | 一种考虑轴向变形的斜齿圆柱齿轮时变啮合刚度计算方法 |
US10088386B2 (en) | 2016-11-09 | 2018-10-02 | Beijing University Of Technology | Device and method for measuring three-dimensional contact stiffness of spur gear based on rough surface |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8061229B2 (en) * | 2005-07-05 | 2011-11-22 | Roman Vasilyevich Novikov | Gear drive |
CN103577687A (zh) * | 2013-09-23 | 2014-02-12 | 北京工业大学 | 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法 |
CN104820756A (zh) * | 2015-05-18 | 2015-08-05 | 东北大学 | 一种考虑延长啮合的裂纹齿轮转子***动力参数确定方法 |
-
2015
- 2015-09-29 CN CN201510634304.8A patent/CN105224744B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8061229B2 (en) * | 2005-07-05 | 2011-11-22 | Roman Vasilyevich Novikov | Gear drive |
CN103577687A (zh) * | 2013-09-23 | 2014-02-12 | 北京工业大学 | 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法 |
CN104820756A (zh) * | 2015-05-18 | 2015-08-05 | 东北大学 | 一种考虑延长啮合的裂纹齿轮转子***动力参数确定方法 |
Non-Patent Citations (1)
Title |
---|
赵树滨等: "势能法计算剥落故障齿轮时变啮合刚度", 《机械设计与制造》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106844818A (zh) * | 2016-11-09 | 2017-06-13 | 北京工业大学 | 基于粗糙表面的直齿轮三维接触刚度计算方法 |
CN106844818B (zh) * | 2016-11-09 | 2018-03-30 | 北京工业大学 | 基于粗糙表面的直齿轮三维接触刚度计算方法 |
WO2018086160A1 (zh) * | 2016-11-09 | 2018-05-17 | 北京工业大学 | 基于粗糙表面的直齿轮三维接触刚度计算方法 |
US10088386B2 (en) | 2016-11-09 | 2018-10-02 | Beijing University Of Technology | Device and method for measuring three-dimensional contact stiffness of spur gear based on rough surface |
CN107314845A (zh) * | 2017-07-03 | 2017-11-03 | 重庆理工大学 | 齿轮动态啮合力的测试方法 |
CN107729626A (zh) * | 2017-09-26 | 2018-02-23 | 西安交通大学 | 一种基于概率分布的齿轮点蚀模型建模方法 |
CN107729626B (zh) * | 2017-09-26 | 2020-03-17 | 西安交通大学 | 一种基于概率分布的齿轮点蚀模型建模方法 |
CN107798200A (zh) * | 2017-11-10 | 2018-03-13 | 西安电子科技大学 | 一种考虑轴向变形的斜齿圆柱齿轮时变啮合刚度计算方法 |
CN107798200B (zh) * | 2017-11-10 | 2019-12-24 | 西安电子科技大学 | 一种考虑轴向变形的斜齿圆柱齿轮时变啮合刚度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105224744B (zh) | 2018-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105224744A (zh) | 一种基于啮合刚度的剥落齿轮啮合模型 | |
Ma et al. | An improved analytical method for mesh stiffness calculation of spur gears with tip relief | |
CN103577687B (zh) | 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法 | |
CN103698124B (zh) | 一种测量齿轮传动误差的方法 | |
CN110375054B (zh) | 一种基于齿廓倾斜偏差的非对称齿轮设计方法 | |
Han et al. | Influences of tooth spalling or local breakage on time-varying mesh stiffness of helical gears | |
CN105181327A (zh) | 一种裂纹轮齿啮合刚度计算的方法 | |
CN104568083A (zh) | 机械计数器及使用其测量燃气表机电转换耐久性的方法 | |
CN104565289B (zh) | 大齿向修形量齿轮的设计方法及大齿向修形量齿轮 | |
CN203035890U (zh) | 一种新型齿轮消间隙机构 | |
CN103174814A (zh) | 传动*** | |
CN109101705A (zh) | 一种基于齿廓普遍方程的行星齿轮时变啮合刚度计算方法 | |
CN107131282A (zh) | 融合齿背接触机理的高速圆柱齿轮动态啮合刚度计算方法 | |
CN204921855U (zh) | 一种低噪声的柴油机齿轮传动结构 | |
CN101015872B (zh) | 一种可以解决z<17的渐开线齿轮“根切”问题的滚刀齿形 | |
CN204704376U (zh) | 一种钢丝加载的双片齿轮消隙机构 | |
CN203962916U (zh) | 具有抗弯曲和耐冲击性能的直齿圆柱齿轮 | |
CN201965421U (zh) | 一种擒纵机构 | |
CN205715469U (zh) | 一种防轴裂与齿裂的行星减速机 | |
CN205013633U (zh) | 适用于装载机变矩器的改良结构 | |
CN202790338U (zh) | 一种共轴直升机用减速齿轮组 | |
Ohshima et al. | Study on High Reduction Face Gears(1 st Report, Meshing Analysis) | |
CN203809606U (zh) | 一种升速机 | |
CN203585218U (zh) | 多点啮合外齿轮减速机 | |
CN212055726U (zh) | 一种用于康复器械传动的齿轮 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |