CN108444697B - 渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法 - Google Patents
渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法 Download PDFInfo
- Publication number
- CN108444697B CN108444697B CN201810193031.1A CN201810193031A CN108444697B CN 108444697 B CN108444697 B CN 108444697B CN 201810193031 A CN201810193031 A CN 201810193031A CN 108444697 B CN108444697 B CN 108444697B
- Authority
- CN
- China
- Prior art keywords
- gear
- fatigue
- contact
- risk
- stress
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/02—Gearings; Transmission mechanisms
- G01M13/021—Gearings
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Gears, Cams (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法,它包括以下步骤:1、根据齿轮材料的硬度曲线以及硬度与屈服强度的线性关系,得出齿轮局部屈服强度;2、依据齿轮副的几何参数,将啮合位置的接触状态简化为刚性半圆与柔性体二维接触模型,同时基于ABAQUS平台建立该二维接触模型;3、基于Python编程语言和材料局部屈服强度,为二维接触模型添加材料属性;4、结合齿轮工作工况,计算齿面以下的Dang Van等效应力;5、以Dang Van等效应力除以齿轮局部材料屈服强度,得出关键接触区域任一点处的疲劳失效风险值。本发明的技术效果是:预测硬化齿轮弹塑性接触疲劳点蚀与疲劳剥落的风险,对硬化齿轮的设计、使用和制造提供指导作用。
Description
技术领域
本发明属于机械零部件疲劳失效风险的预测方法,具体涉及一种预测以渗碳淬火为代表的表面硬化齿轮副弹塑性接触疲劳点蚀与剥落失效风险的方法。
背景技术
在高速、重载机械中,渗碳硬化齿轮的接触疲劳失效问题已成为限制装备可靠性的重要瓶颈之一,硬化齿轮的接触疲劳失效形式常见的有点蚀和剥落两种,这两种失效问题会带来极大地安全隐患和经济损失,而硬化齿轮的接触疲劳失效影响因素众多,如工况、材料力学特征的梯度、残余应力等,此外齿轮接触过程中伴随时变多轴应力状态,使得目前对于渗碳重载齿轮的接触疲劳机理的认识与控制不够深入,工程实际中对于失效风险的预测存在很大困难。
发明内容
本发明所要解决的技术问题是提供一种预测表面硬化齿轮弹塑性接触疲劳点蚀与剥落失效风险的方法,它能预测齿轮在弹塑性接触状态下疲劳点蚀和疲劳剥落的失效风险,预测结果在工程实际中具有指导作用,避免因齿轮失效而导致的突发事故、因装置设备停机造成生产效益的损失。
本发明所要解决的技术问题是通过这样的技术方案实现的,它包括以下步骤:
步骤1、根据齿轮材料的硬度曲线和硬度与屈服强度的线性关系,计算得到齿轮材料局部屈服强度;
步骤2、依据齿轮副的几何参数,将啮合位置的接触状态简化为二维接触模型,同时基于ABAQUS平台建立该模型;
步骤3、基于Python编程语言和齿轮材料局部屈服强度,为二维接触模型添加材料属性;
步骤4、根据齿轮的材料和Dang Van多轴应力准则,结合齿轮的具体工况,计算齿轮齿面下的Dang Van等效应力;
步骤5、以Dang Van等效应力除以齿轮材料局部屈服强度,得出关键接触区域任一点处的疲劳失效风险值;根据失效风险值在不同区域相对大小和最大值所在的深度位置预测点蚀或剥落失效的风险。
本发明的技术效果是:
解决了机械行业长期存在的技术难题,能在考虑齿轮硬化后材料力学性能梯度变化的条件下,分析齿轮在弹塑性状态下的点蚀和剥落的风险,为齿轮的设计、制造、使用提供依据。
附图说明
本发明的附图说明如下:
图1为渗碳硬化齿轮失效风险区分布图;
图2为表面硬化齿轮硬度随浅层深度变化的曲线图;
图3为实施例中齿轮材料局部屈服强度沿深度分布曲线;
图4为齿轮啮合接触状态的简化示意图;
图5为实施例中有限元模型示意图;
图6为实施例中某兆瓦级风电齿轮箱的传动图;
图7为实施例中Dang Van等效应力沿深度分布的曲线;
图8为实施例的齿轮的接触疲劳失效风险值曲线。
具体实施方式
下面结合附图和实施例对本发明作进一步说明:
如图1所示,根据工程经验,在距离齿面几百微米处为疲劳点蚀的高发区,而距离齿面距离为1~3mm为疲劳剥落的高发区。图示中的曲线代表疲劳失效风险值沿深度方向的变化曲线,根据不同区域失效风险的大小以及最大失效风险值的深度位置可以用于预测齿轮发生疲劳点蚀和疲劳剥落失效的风险,若最大值发生在点蚀高发区,则齿轮在给定工况下的点蚀风险较高,若最大失效风险值在剥落高发区发生,则齿轮剥落失效风险较高。
如图2所示,曲线为使用维氏硬度测试法得到的硬度曲线。纵坐标为维氏硬度值,横坐标为距离齿面的深度值,单位为毫米(mm),总测量深度为5mm,测量值为一系列数据点,根据测量数据点拟合出一条曲线,拟合函数为:
HV(y)=a·y4+b·y3+c·y2+d·y+e (1)
式(1)中,a=-1.733,b=25.37,c=-118,d=144.8,e=589.5,y为深度值。
图2所示的硬度值对于不同的齿轮需要重新测量,因为该值与齿轮的材料和热处理工艺等有关,而测量硬度值的方法众多,这里只提供一个参考的实施例,其他方法所获得的硬度值也可以用于本专利提供的方法进行失效风险预测,这里不再赘述。
本发明包括以下步骤:
步骤1、根据齿轮材料的硬度曲线和硬度与屈服强度的线性关系,计算得到齿轮材料局部屈服强度;
根据“Correlation of Yield Strength and Tensile Strength with Hardnessfor Steels”,E.J.Pavlina and C.J.V.Tyne,Journal of Materials Engineering&Performance,vol.17,pp.888-893,2008.(“钢的抗拉强度与屈服强度之间的关系”,E.J.Pavlina and C.J.V.Tyne,材料工程与性能,第888-893页,2008年)记载了维氏硬度与屈服极限的线性关系为:
σYS(y)=-90.7+2.876HV(y) (2)
式中,y代表齿轮深度,σYS为齿轮材料局部屈服强度,HV为沿深度y分布的维氏硬度值。
由公式(2)计算得到的齿轮材料屈服强度沿深度的变化如图3所示。
步骤2、依据齿轮副的几何参数,将啮合位置的接触状态简化为二维接触模型,同时基于ABAQUS平台建立该模型;
齿轮副啮合接触状态简化为二维接触模型的过程如图4所示,左图代表要分析的啮合齿轮副,根据两齿轮啮合点处的曲率半径将分析对象转化为中间图所示的两圆接触,该模型又可进一步简化为一刚性半圆与柔性体的接触(如右图)。此时的刚性半圆的综合曲率半径和柔性体的等效弹性模型的计算方法为:
R=R1R2/(R1+R2) (3)
上式中,R1,R2为两齿轮接触位置的曲率半径,R为综合曲率半径,E1,E2为两齿轮的弹性模量,E为等效弹性模量,υ1,υ2为两齿轮的泊松比。“AGMA information sheet 908-B89,1989.“Geometry factors for determining the pitting resistance and bendingstrength of spur,helical and herringbone gear teeth”.”(美国齿轮标准1989年的“判定直齿轮,斜齿轮,人字形齿轮轮齿的抗点蚀能力和弯曲强度的几何影响因子”)在第5-7页记载了R1,R2的计算方法。
步骤3、基于Python编程语言和步骤1得到的齿轮材料局部屈服强度,为二维接触模型添加材料属性;
基于Python语言对ABAQUS有限元软件进行二次开发,编写对柔性体进行分层的程序。
如图5所示,左图表示使用Python语言对柔性体上部分进行分层后的局部放大图,单层厚度为0.01mm,中间图代表在ABAQUS平台中建立的二维接触模型,右图为该模型添加齿轮材料局部屈服强度后的示意图,右图中的曲线代表从齿面到芯部的屈服强度分布。
在该步骤中,可以根据实际工程需要,对采用不同的单层厚度,均属于本专利的范畴,在此不再赘述。
步骤4、依据步骤3中建立的模型,计算齿面下的Dang Van等效应力;
"On a new multiaxial fatigue limit criterion:theory and application,"K.Dang Van,B.Griveau,and O.Message,Biaxial and Multiaxial Fatigue,M.W.Brownand K.J.Miller,Eds.,ed:Mechanical Engineering Publications,London,1989,pp.479-496.(“一种新的多轴疲劳极限准则:理论与应用”,K.Dang Van,B.Griveau,andO.Message,两轴与多轴疲劳,M.W.Brown and K.J.Miller,编辑,机械工程出版,伦敦,1989年,479-496页)中提出了Dang Van等效应力的计算公式:
τmax(θ,t)+ασh(t)=σequi,DangVan(y) (5)
上述式子中,σequi,DangVan(y)是Dang Van多轴疲劳准则的等效应力,τmax(θ,t)为最大剪应力幅值,σh(t)为水静应力,α是材料参数;τ-1是扭转疲劳极限,σ-1是弯曲疲劳极限,σ1、σ2和σ3是第一主应力、第二主应力和第三主应力。
"Prediction of contact fatigue for the rough surfaceelastohydrodynamic lubrication line contact problem under rolling and slidingconditions",Q.Hua,PhD Thesis,Cardiff University,UK,2005.(“滑动和滚动条件下的粗糙表面弹流润滑线接触问题的接触疲劳预测”,Q.Hua,博士论文,卡迪夫大学,英国,2005年)中指出最大剪应力幅值的求法为:
τ(θ,t)=τxy(cos2θ-sin2θ)+(σy-σx)sinθcosθ (8)
τmax(θ,t)=max{τa(θ,t)=|τ(θ,t)-τm|} (9)
式(8)、式(9)中,τ(θ,t)是材料点瞬时剪应力,τxy是x-y平面剪应力,σx是x轴向正应力,σy是y轴向正应力,θ代表剪应力幅值达到最大值的平面与滚动方向之间的夹角;τmax(θ,t)是最大剪应力幅值,τa(θ,t)是剪应力幅值,τm是剪应力均值。
步骤5、以Dang Van等效应力除以齿轮局部材料屈服强度,得出关键接触区域任一点处的疲劳失效风险值。疲劳失效风险值越大,疲劳失效的可能性越大,根据风险值在不同区域的大小,预测点蚀或剥落的失效风险。疲劳失效风险值的计算式如下:
式中Aff(y)为疲劳失效风险值。
实施例
如6图所示,样品齿轮来自于某兆瓦级风力发电机齿轮箱中间及齿轮副,它在工程实际应用时,失效的概率明显高于其他齿轮。
齿轮副的主要参数如下:
图2为该样本所对应的拟合硬度曲线,根据实施例齿轮参数,以齿轮啮合线上节点处啮合瞬时状态为例,逐步计算出齿轮副接触的疲劳时效风险值:
步骤1、根据齿轮硬度曲线和公式(2),可得齿轮材料局部屈服强度,如图3所示。
步骤2、采用图4所示的简化过程将啮合齿轮简化为二维接触模型。根据公式(3)~(4),由R1=236.28mm,R2=44.18mm,计算得到齿轮副综合曲率半径R=37.22mm;等效弹性模量E=1.15×1011Pa;同时在ABAQUS平台建立该二维接触模型。
步骤3、使用Python语言和步骤1得到的齿轮材料局部屈服强度,为二维接触模型添加材料参数。Python语言程序示例如下:
for i in range(1,500):
s.Line(point1=(-20.0,10-0.01*i),point2=(20.0,10-0.01*i))
p=mdb.models['Model-1'].parts['deformable']
f=p.faces
pickedFaces=f.findAt(((6.666667,-6.666667,0.0),))
e1,d2=p.edges,p.datums
p.PartitionFaceBySketch(faces=pickedFaces,sketch=s)
s.unsetPrimaryObject()
del mdb.models['Model-1'].sketches['__profile__']
步骤4、在ABAQUS平台使用公式(5)~(9)计算Dang Van等效应力,得到Dang Van等效应力沿深度分布如图7所示。Dang Van等效应力在表面以下达到最大值,最大值深度位置约为500um,整条曲线先上升后下降。
步骤5、根据公式(10)计算出材料的局部疲劳失效风险值,如图8所示,失效风险值在齿面以下达到最大值,最大值深度位置约为500um,最大失效风险值约为0.3,失效风险值较低。整条曲线先上升后下降,最大失效风险值在疲劳点蚀高发区。
“Tooth Flank Fracture–Basic Principles and Calculation Model for aSub-Surface-Initiated Fatigue Failure Mode of Case-Hardened Gears”,I.Boiadjiev,J.Witzig,T.Tobie and K.Stahl,Gear Technology,pp.59-64,2015(“齿面断裂-表面硬化齿轮的次表面疲劳失效的基本准则与计算模型”,I.Boiadjiev,J.Witzig,T.Tobie和K.Stahl,齿轮技术,第59-64页,2015年)记载了考虑齿面硬化的失效风险值的实验,实验得到的沿深度分布的失效风险值曲线趋势与本发明得到的趋势基本一致,由此验证了本方法发明的实用性和可靠性。
Claims (6)
1.渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法,其特征是,包括以下步骤:
步骤1、根据齿轮材料的硬度曲线和硬度与屈服强度的线性关系,计算得到齿轮材料局部屈服强度;
步骤2、依据齿轮副的几何参数,将啮合位置的接触状态简化为二维接触模型,同时基于ABAQUS平台建立该模型;
步骤3、基于Python编程语言和齿轮材料局部屈服强度,为二维接触模型添加材料属性;
步骤4、根据齿轮的材料和Dang Van多轴应力准则,结合齿轮的具体工况,计算齿轮齿面下的Dang Van等效应力;
步骤5、以Dang Van等效应力除以齿轮材料局部屈服强度,得出关键接触区域任一点处的疲劳失效风险值;根据疲劳失效风险值在不同区域相对大小和最大值所在的深度位置预测点蚀或剥落失效的风险。
2.根据权利要求1所述的渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法,其特征是,在步骤1中,所述的齿轮材料局部屈服强度为:
σYS(y)=-90.7+2.876HV(y)
式中,y代表沿齿面向下的深度,σYS(y)为齿轮材料局部屈服强度,HV(y)为沿深度y分布的维氏硬度值。
4.根据权利要求3所述的渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法,其特征是,在步骤3中,所述的为二维接触模型添加材料属性的方法为:
将二维接触模型中的柔性体,使用Python编程语言进行分层,为每一层添加不同的屈服强度。
5.根据权利要求4所述的渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法,其特征是,在步骤4中,所述的Dang Van等效应力计算方法为:
τmax(θ,t)+ασh(t)=σequi,DangVan(y)
式中σequi,DangVan(y)是Dang Van等效应力,τmax(θ,t)为最大剪应力幅值,σh(t)为水静应力,α是材料参数,τ-1是扭转疲劳极限,σ-1是弯曲疲劳极限,σ1、σ2和σ3是第一主应力、第二主应力和第三主应力;
τ(θ,t)=τxy(cos2θ-sin2θ)+(σy-σx)sinθcosθ
τmax(θ,t)=max{τa(θ,t)=|τ(θ,t)-τm|}
式中,τ(θ,t)是材料点瞬时剪应力,τxy是x-y平面剪应力,σx是x轴向正应力,σy是y轴向正应力,θ代表剪应力幅值达到最大值的平面与滚动方向之间的夹角;τa(θ,t)是剪应力幅值,τm是剪应力均值。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711254805 | 2017-12-04 | ||
CN2017112548059 | 2017-12-04 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108444697A CN108444697A (zh) | 2018-08-24 |
CN108444697B true CN108444697B (zh) | 2020-04-07 |
Family
ID=63193908
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810193031.1A Active CN108444697B (zh) | 2017-12-04 | 2018-03-09 | 渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108444697B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109271711B (zh) * | 2018-09-25 | 2023-03-28 | 重庆大学 | 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法 |
CN109299559B (zh) * | 2018-10-08 | 2023-05-30 | 重庆大学 | 一种表面硬化齿轮磨损及疲劳失效竞争机制分析方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006326634A (ja) * | 2005-05-25 | 2006-12-07 | Daido Steel Co Ltd | 金型の損傷を予測する方法 |
CN104699976B (zh) * | 2015-03-20 | 2017-07-18 | 北京航空航天大学 | 一种包含平均应力影响的金属材料多轴高周疲劳失效预测方法 |
CN105466772B (zh) * | 2015-11-15 | 2018-04-27 | 北京工业大学 | 一种基于临界面法的多轴短裂纹扩展寿命预测方法 |
CN106126856B (zh) * | 2016-07-06 | 2019-03-01 | 西南交通大学 | 适用于负应力比的基于低周疲劳性能参数的裂纹扩展速率预测方法 |
CN106886663B (zh) * | 2017-03-29 | 2020-04-03 | 北京理工大学 | 齿轮弯曲疲劳寿命预测方法及装置 |
CN106979861B (zh) * | 2017-03-30 | 2019-04-23 | 北京理工大学 | 齿轮接触疲劳全寿命评估方法及装置 |
-
2018
- 2018-03-09 CN CN201810193031.1A patent/CN108444697B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108444697A (zh) | 2018-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Evans et al. | Analysis of micro-elastohydrodynamic lubrication and prediction of surface fatigue damage in micropitting tests on helical gears | |
Wang et al. | Evaluation of rolling contact fatigue of a carburized wind turbine gear considering the residual stress and hardness gradient | |
CN109299559B (zh) | 一种表面硬化齿轮磨损及疲劳失效竞争机制分析方法 | |
Li et al. | Micro-pitting fatigue lives of lubricated point contacts: Experiments and model validation | |
Zhang et al. | Numerical simulation of competing mechanism between pitting and micro-pitting of a wind turbine gear considering surface roughness | |
Liu et al. | Study on contact fatigue of a wind turbine gear pair considering surface roughness | |
Li et al. | Prediction of surface wear of involute gears based on a modified fractal method | |
Bonaiti et al. | Gear root bending strength: A comparison between Single Tooth Bending Fatigue Tests and meshing gears | |
CN108318244A (zh) | 考虑残余应力的渗碳硬化齿轮接触疲劳风险评估方法 | |
CN108444697B (zh) | 渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法 | |
Deng et al. | Finite element analysis of contact fatigue and bending fatigue of a theoretical assembling straight bevel gear pair | |
Abbasi et al. | Effect of out-of-phase loading on fretting fatigue response of Al7075-T6 under cyclic normal loading using a new testing apparatus | |
Bonaiti et al. | Estimation of gear SN curve for tooth root bending fatigue by means of maximum likelihood method and statistic of extremes | |
Clarke et al. | Effects of profile errors on lubrication performance of helical gears | |
Zhao et al. | Application study of fractal theory in mechanical transmission | |
CN108256241A (zh) | 一种重载齿轮次表面裂纹萌生的预测方法 | |
CN112651562B (zh) | 基于信噪比和改进灰色关联度激光熔覆工艺优化方法 | |
Xie et al. | Investigations on transitions of lubrication states for water lubricated bearing. Part I: Determination of friction coefficients and film thickness ratios | |
Wang et al. | Evaluation of contact fatigue risk of a carburized gear considering gradients of mechanical properties | |
Hu et al. | Investigation on load distribution among rollers of planetary roller screw mechanism considering machining errors: Analytical calculation and machine learning approach | |
Ahmadi et al. | A three-dimensional finite element damage mechanics model to simulate fretting wear of hertzian line and circular contacts in partial slip regime | |
Honkalas et al. | A review on design and efficiency improvement of worm and worm wheel of a gear motor | |
Cortabitarte et al. | Application of the theory of critical distances for fatigue life assessment of spur gears | |
Fabre et al. | Microgeometrical influences on micropitting fatigue damage: multi-scale analysis | |
Zhang et al. | Interaction between gear surface wear and contact fatigue |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |