CN106979861B - 齿轮接触疲劳全寿命评估方法及装置 - Google Patents
齿轮接触疲劳全寿命评估方法及装置 Download PDFInfo
- Publication number
- CN106979861B CN106979861B CN201710200764.9A CN201710200764A CN106979861B CN 106979861 B CN106979861 B CN 106979861B CN 201710200764 A CN201710200764 A CN 201710200764A CN 106979861 B CN106979861 B CN 106979861B
- Authority
- CN
- China
- Prior art keywords
- gear
- gear contact
- model
- stress
- fatigue life
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000003068 static effect Effects 0.000 claims abstract description 40
- 230000000977 initiatory effect Effects 0.000 claims description 23
- 208000037656 Respiratory Sounds Diseases 0.000 claims description 22
- 239000000463 material Substances 0.000 claims description 19
- 238000009826 distribution Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 6
- 125000004122 cyclic group Chemical group 0.000 claims description 6
- 238000010008 shearing Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000002929 anti-fatigue Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009661 fatigue test Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000001050 lubricating effect Effects 0.000 description 1
- 239000010687 lubricating oil Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000005482 strain hardening Methods 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提供一种基于裂纹萌生‑扩展的齿轮接触疲劳全寿命评估方法及装置,该方法包括:根据最大接触应力模型计算齿轮接触面上的最大接触应力;基于数值计算理论及等效边界条件,分别构建齿轮的二维静态模型、三维静态模型、二维动态模型及三维动态模型四个模型,并基于四个模型分别获得对应的齿轮接触面上的最大接触应力;将基于四个模型获得的对应最大接触应力分别与基于最大接触应力模型计算得到的最大接触应力进行比较,确定最佳数值计算模型;构建齿轮接触疲劳萌生寿命评估模型;确定有效应力强度因子值范围;构建齿轮接触疲劳扩展寿命评估模型;根据齿轮接触疲劳萌生寿命评估模型及齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估。
Description
技术领域
本发明是关于齿轮接触疲劳全寿命评估技术,特别是关于一种基于裂纹萌生扩展的齿轮接触疲劳全寿命评估方法及装置。
背景技术
齿轮作为动力传动的主体,齿面的接触疲劳破坏是最普遍的一种疲劳失效模式,明确齿轮真实接触应力分布,评估齿轮接触疲劳寿命,已成为齿轮抗疲劳设计的重要依据。
传统评估齿轮接触疲劳寿命的方法是依据大量齿轮接触疲劳试验获得齿轮接触S-N曲线,进而评估齿轮的接触疲劳寿命。但这种传统方法很难揭示齿轮接触疲劳失效的规律及机理,尤其是疲劳裂纹的萌生-扩展举止。而且,也极少考虑表面残余应力、齿轮实际承受的多轴应力状态等因素的影响。此外,从试验的角度,传统方法主要建立在大量试验的基础上,成本较高且周期较长。其方法的可靠性和适用性密切相关于特定尺寸齿轮的接触疲劳试验数据。当齿轮材料、模数及齿数等参数发生变化时,试验必须重新开展,这就需花费大量的时间和成本。
因此,亟待一种从实际齿轮应力状态出发,虑及残余应力的影响,减小对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性,能比较稳定与准确地评估齿轮接触疲劳全寿命的方法。
发明内容
本发明实施例提供了一种基于裂纹萌生扩展的齿轮接触疲劳全寿命评估方法及装置,以比较稳定与准确地评估齿轮接触疲劳全寿命,减小对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性。
为了实现上述目的,本发明实施例提供了一种基于裂纹萌生扩展的齿轮接触疲劳全寿命评估方法,该齿轮接触疲劳全寿命评估方法包括:
根据最大接触应力模型计算齿轮接触面上的最大接触应力;
基于数值计算理论及等效边界条件,分别构建齿轮的二维静态模型、三维静态模型、二维动态模型及三维动态模型,并基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型分别获得对应的齿轮接触面上的最大接触应力;
将基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型获得的对应最大接触应力分别与基于最大接触应力模型计算得到的最大接触应力进行比较,确定最佳数值计算模型;
基于位错-能量法及裂纹尺寸构建齿轮接触疲劳萌生寿命评估模型;
在半球空间内,基于裂纹扩展中的应力分布及滑开型裂纹确定有效应力强度因子值范围;
基于Paris方程、裂纹扩展角度、裂纹增量尺寸、齿轮材料硬度及所述有效应力强度因子值范围,构建齿轮接触疲劳扩展寿命评估模型;
根据所述齿轮接触疲劳萌生寿命评估模型及所述齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估。
一实施例中,所述齿轮接触疲劳萌生寿命评估模型为:
其中,Ni为接触疲劳萌生寿命,l为椭圆滑移带的半长轴,G为剪切模量,Δτ为剪切应力范围,k为位错摩擦应力,c为裂纹长度,d为晶粒尺寸。
一实施例中,该齿轮接触疲劳全寿命评估方法还包括:基于所述最佳数值计算模型对所述齿轮接触疲劳萌生寿命评估模型的拟合系数及椭圆滑移带的半长轴进行参数修正。
一实施例中,该齿轮接触疲劳全寿命评估方法还包括:
基于渗碳层中残余应力对齿轮接触疲劳萌生寿命的影响,对参数修正后的所述齿轮接触疲劳萌生寿命评估模型进行修正;
根据所述齿轮接触疲劳萌生寿命评估模型及所述齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估,包括:
根据所述齿轮接触疲劳扩展寿命评估模型及修正后的所述齿轮接触疲劳萌生寿命评估模型进行齿轮接触疲劳全寿命评估。
一实施例中,该齿轮接触疲劳全寿命评估方法还包括:
基于刚塑性滑移线理论,结合每次循环载荷作用后最大剪切应力,确定裂纹扩展角度;
基于裂纹尖端塑性区域确定所述裂纹增量尺寸。
一实施例中,所述齿轮接触疲劳扩展寿命评估模型为:
其中,Np为齿轮接触疲劳扩展寿命,C和m分别为材料参数,L为裂纹尺寸,xL为裂纹尖端在X轴方向上的坐标,ξ为裂纹扩展增量,τmax为齿轮次表面中最大剪切应力,Hb及HL分别为齿轮整体硬度及局部硬度,a0为裂纹初始长度。
为了实现上述目的,本发明实施例提供了一种基于裂纹萌生扩展的齿轮接触疲劳全寿命评估装置,包括:
第一最大接触应力计算单元,用于根据最大接触应力模型计算齿轮接触面上的最大接触应力;
第二最大接触应力计算单元,用于基于数值计算理论及等效边界条件,分别构建齿轮的二维静态模型、三维静态模型、二维动态模型及三维动态模型,并基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型分别获得对应的齿轮接触面上的最大接触应力;
最佳数值计算模型确定单元,用于将基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型获得的对应最大接触应力分别与基于最大接触应力模型计算得到的最大接触应力进行比较,确定最佳数值计算模型;
齿轮接触疲劳萌生寿命评估模型构建单元,用于基于位错-能量法及裂纹尺寸构建齿轮接触疲劳萌生寿命评估模型;
有效应力强度因子值范围确定单元,用于在半球空间内,基于裂纹扩展中的应力分布及滑开型裂纹确定有效应力强度因子值范围;
齿轮接触疲劳扩展寿命评估模型构建单元,用于基于Paris方程、裂纹扩展角度、裂纹增量尺寸、齿轮材料硬度及所述有效应力强度因子值范围,构建齿轮接触疲劳扩展寿命评估模型;
齿轮接触疲劳全寿命评估单元,用于根据所述齿轮接触疲劳萌生寿命评估模型及所述齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估。
一实施例中,还包括:参数修正单元,用于基于所述最佳数值计算模型对所述齿轮接触疲劳萌生寿命评估模型的拟合系数及椭圆滑移带的半长轴进行参数修正。
一实施例中,该齿轮接触疲劳全寿命评估装置还包括:
模型修正单元,用于基于渗碳层中残余应力对齿轮接触疲劳萌生寿命的影响,对参数修正后的所述齿轮接触疲劳萌生寿命评估模型进行修正;
所述齿轮接触疲劳全寿命评估单元具体用于:
根据所述齿轮接触疲劳扩展寿命评估模型及修正后的所述齿轮接触疲劳萌生寿命评估模型进行齿轮接触疲劳全寿命评估。
一实施例中,该齿轮接触疲劳全寿命评估装置还包括:
裂纹扩展角度确定单元,用于基于刚塑性滑移线理论,结合每次循环载荷作用后最大剪切应力,确定裂纹扩展角度;
裂纹增量尺寸确定单元,用于基于裂纹尖端塑性区域确定所述裂纹增量尺寸。
利用本发明,可以比较稳定与准确地评估齿轮接触疲劳全寿命,减小对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例基于裂纹萌生扩展的齿轮接触疲劳全寿命评估方法流程图;
图2本发明实施例裂纹扩展过程中相关参数定义示意图;
图3为本发明实施例基于裂纹萌生扩展的齿轮接触疲劳全寿命评估装置的结构示意图一;
图4为本发明实施例基于裂纹萌生扩展的齿轮接触疲劳全寿命评估装置的结构示意图二;
图5为本发明实施例基于裂纹萌生扩展的齿轮接触疲劳全寿命评估装置的结构示意图三。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例基于裂纹萌生扩展的齿轮接触疲劳全寿命评估方法流程图,如图1所示,包括:
S101:根据最大接触应力模型计算齿轮接触面上的最大接触应力。
根据国标GB/T 14229-1003《齿轮接触疲劳强度试样方法》记载的最大接触应力模型,可以计算啮合齿轮在接触面上的最大接触应力σH,最大接触应力模型为如下公式:
上述公式(1)为经验公式,公式(1)中,Ft为端面内节圆周上的名义切向力,单位为N;ZH为节点区域系数;ZE为弹性系数;Zε为接触强度计算的重合度系数;Zβ为接触强度计算的螺旋角系数;ZV为速度系数;ZL为润滑油系数;ZR为粗糙度系数;ZW为齿面工作硬化系数;ZX为接触强度计算的尺寸系数;u为齿数比;KA为使用系数;KV为动载系数;KHα为接触强度计算的齿间载荷分布系数;KHβ为接触强度计算的齿向载荷分布系数;d1为试验齿轮小轮节圆直径;b为工作齿宽。
S102:基于数值计算理论及等效边界条件,分别构建齿轮的二维静态模型、三维静态模型、二维动态模型及三维动态模型,并基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型分别获得对应的齿轮接触面上的最大接触应力。
具体实施时,可以基于ABAQUS中Standard和Explicit两种求解分析模块,结合等效边界条件及载荷,分别构建齿轮二维静态模型、三维静态模型、二维动态模型及三维动态模型四个模型,并依次得到四个模型分别对应的齿轮接触面上的最大接触应力。其中,Standard模块主要用于求解一些静态问题,不考虑质量和冲击载荷的影响。而Explicit模块则是用来模拟高速冲击问题,且在求解中考虑惯性的影响。
S103:将基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型获得的对应最大接触应力分别与基于最大接触应力模型计算得到的最大接触应力进行比较,确定最佳数值计算模型。
将基于S102中的四种模型的最大接触应力计算结果与步骤1中对最大接触应力结果进行比较,确定最佳数值计算模型。比较分析可知,由于三维动态模型中虑及了动态冲击载荷及润滑状态的影响,因此静态分析模型的计算精度不如动态分析模型的计算精度高,三维动态模型具有最好的预测精度。
S104:基于位错-能量法及裂纹尺寸构建齿轮接触疲劳萌生寿命评估模型。
具体实施时,可以基于位错滑移带内部分位错偶极子用于疲劳萌生裂纹形成的假设,考虑裂纹尺寸(长度),构建齿轮接触疲劳萌生寿命评估模型,该齿轮接触疲劳萌生寿命评估模型为:
公式(2)中,Ni为接触疲劳萌生寿命,l为椭圆滑移带的半长轴,G为剪切模量,Δτ为剪切应力范围,k为位错摩擦应力,c为裂纹长度,d为晶粒尺寸。
一实施例中,在得到上述齿轮接触疲劳萌生寿命评估模型后,还可以基于S103得到的最佳数值计算模型对所述齿轮接触疲劳萌生寿命评估模型的拟合系数及椭圆滑移带的半长轴进行参数修正。
具体地,基于Von Misses八面体剪应力屈服准则,在单轴拉压状态时,八面体剪应力τoct和材料屈服强度σy的关系,可以表达为:
基于公式(3),可以得到:
其中,σw为疲劳极限,τf为临界剪应力,σMises为等效应力。由于k也可以认为是临界剪应力τf,因此公式(2)可改写为:
公式(6)开方可得:
为提高预测寿命与试验数据之间的拟合精度,可以用拟合系数α(0<α<1)对萌生寿命模型进行修正,即公式(7)可变形为:
因此,基于位错能量法的接触疲劳萌生寿命模型可表示为:
因此,通过最佳数值计算模型可以计算得到等效应力σMises,然后根据试验数据点以及公式(9)或(10),可以反求椭圆滑移带的半长轴l和拟合系数α。
一实施例中,可以基于渗碳层中残余应力对齿轮接触疲劳萌生寿命的影响,对参数修正后的齿轮接触疲劳萌生寿命评估模型进行修正。修正后的齿轮接触疲劳萌生寿命评估模型如下:
式中,σrs为残余应力,σb为抗拉强度。
S105:在半球空间内,基于裂纹扩展中的应力分布及滑开型裂纹确定有效应力强度因子值范围。
在半球空间内,基于裂纹扩展中的应力分布,针对滑开型裂纹形态,应力强度因子值KL的计算公式可被表达式为:
计算过程中相关参数的定义可以参考图2所示。图2中,L为裂纹长度,xL为裂纹尖端在X轴方向上的坐标,xi为裂纹初始点在X轴方向上的坐标,ξ为裂纹扩展增量,PH(max)为压应力,θ为裂纹扩展角度。
此外,在不考虑齿轮啮合过程中摩擦力影响时,齿轮次表面中最大剪切应力τmax的计算公式可表示为:
式中,σx为齿轮次表面内切向应力,σy为齿轮次表面内法向应力,表达式分别为:
式中,yi为裂纹扩展尖端距齿轮表面的距离,r为齿轮接触面半宽长度。当考虑由渗碳工艺引起的残余应力影响时,公式(12)可改写为:
因此,裂纹尖端应力强度因子范围可表示为:
ΔK=Kmax-Kmin (16)
式中,Kmax及Kmin分别为任意裂纹长度下裂纹尖端最大、最小应力强度因子值。
S106:基于Paris方程、裂纹扩展角度、裂纹增量尺寸、齿轮材料硬度及所述有效应力强度因子值范围,构建齿轮接触疲劳扩展寿命评估模型。
首先,需要基于最大剪应力准则确定裂纹扩展角度及裂纹增量尺寸。
基于刚塑性滑移线理论,结合每次循环载荷作用后最大剪切应力,裂纹扩展角度可表示为:
此外,基于裂纹尖端塑性区域,裂纹扩展增量尺寸Δa可表示为:
Δa=0.1*rp (18)
其中,rp为裂纹尖端塑性区域半径,其表达式为:
式中,χ为与应力状态相关的系数,对应平面应力状态时,其值为2。
基于Prais公式,结合齿轮材料硬度,接触疲劳裂纹扩展寿命可表示为:
公式(20)中,Hb及HL分别为齿轮整体硬度及局部硬度。基于公式(11),在积分区间[a0,L]内对公式(20)进行积分,可得接触疲劳裂纹扩展寿命模型,其表达式为:
公式(21)中,C和m分别为材料参数。a0为裂纹初始长度,表达式可表示为:
式中,ΔKth为应力强度因子门槛值范围,可由经验公式得到:
ΔKth=2.45+3.41×10-3HV (23)
S107:根据齿轮接触疲劳萌生寿命评估模型及齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估。
S107中的齿轮接触疲劳萌生寿命评估模型可以为公式(2)中的模型,也可以为公式(9)或(10)中模型。
一实施例中,进行齿轮接触疲劳全寿命评估时,可以根据公式(10)的齿轮接触疲劳萌生寿命评估模型及公式(21)的齿轮接触疲劳扩展寿命评估模型得到齿轮接触疲劳全寿命评估模型:
基于上述齿轮接触疲劳全寿命评估模型,可以进行齿轮接触疲劳全寿命评估。
利用本发明,可以比较稳定与准确地评估齿轮接触疲劳全寿命,减小对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性。由于本发明虑及裂纹萌生-扩展及残余应力,可以提高计算准确度,降低试验成本低。
基于与上述齿轮接触疲劳全寿命评估方法相同的发明构思,本申请提供一种齿轮接触疲劳全寿命评估装置,如下面实施例所述。由于该齿轮接触疲劳全寿命评估装置解决问题的原理与齿轮接触疲劳全寿命评估方法相似,因此该齿轮接触疲劳全寿命评估装置的实施可以参见齿轮接触疲劳全寿命评估方法的实施,重复之处不再赘述。
图3为本发明实施例基于裂纹萌生扩展的齿轮接触疲劳全寿命评估装置的结构示意图,如图3所示,该齿轮接触疲劳全寿命评估装置包括:
第一最大接触应力计算单元301,用于根据最大接触应力模型计算齿轮接触面上的最大接触应力;
第二最大接触应力计算单元302,用于基于数值计算理论及等效边界条件,分别构建齿轮的二维静态模型、三维静态模型、二维动态模型及三维动态模型,并基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型分别获得对应的齿轮接触面上的最大接触应力;
最佳数值计算模型确定单元303,用于将基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型获得的对应最大接触应力分别与基于最大接触应力模型计算得到的最大接触应力进行比较,确定最佳数值计算模型;
齿轮接触疲劳萌生寿命评估模型构建单元304,用于基于位错-能量法及裂纹尺寸构建齿轮接触疲劳萌生寿命评估模型;
有效应力强度因子值范围确定单元305,用于在半球空间内,基于裂纹扩展中的应力分布及滑开型裂纹确定有效应力强度因子值范围;
齿轮接触疲劳扩展寿命评估模型构建单元306,用于基于Paris方程、裂纹扩展角度、裂纹增量尺寸、齿轮材料硬度及所述有效应力强度因子值范围,构建齿轮接触疲劳扩展寿命评估模型;
齿轮接触疲劳全寿命评估单元307,用于根据所述齿轮接触疲劳萌生寿命评估模型及所述齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估。
一实施例中,所述齿轮接触疲劳萌生寿命评估模型为:
其中,Ni为接触疲劳萌生寿命,l为椭圆滑移带的半长轴,G为剪切模量,Δτ为剪切应力范围,k为位错摩擦应力,c为裂纹长度,d为晶粒尺寸。
一实施例中,如图4所示,该齿轮接触疲劳全寿命评估装置还包括:参数修正单元401,用于基于所述最佳数值计算模型对所述齿轮接触疲劳萌生寿命评估模型的拟合系数及椭圆滑移带的半长轴进行参数修正。
一实施例中,该齿轮接触疲劳全寿命评估装置还包括:模型修正单元402,用于基于渗碳层中残余应力对齿轮接触疲劳萌生寿命的影响,对参数修正后的所述齿轮接触疲劳萌生寿命评估模型进行修正。该实施例中,齿轮接触疲劳全寿命评估单元307具体用于:根据所述齿轮接触疲劳扩展寿命评估模型及修正后的所述齿轮接触疲劳萌生寿命评估模型进行齿轮接触疲劳全寿命评估。
一实施例中,如图5所示,该齿轮接触疲劳全寿命评估装置还包括:
裂纹扩展角度确定单元501,用于基于刚塑性滑移线理论,结合每次循环载荷作用后最大剪切应力,确定裂纹扩展角度;
裂纹增量尺寸确定单元502,用于基于裂纹尖端塑性区域确定所述裂纹增量尺寸。
所述齿轮接触疲劳扩展寿命评估模型为:
其中,Np为齿轮接触疲劳扩展寿命,C和m分别为材料参数,L为裂纹尺寸,xL为裂纹尖端在X轴方向上的坐标,ξ为裂纹扩展增量,τmax为齿轮次表面中最大剪切应力,Hb及HL分别为齿轮整体硬度及局部硬度,a0为裂纹初始长度。
利用本发明,可以比较稳定与准确地评估齿轮接触疲劳全寿命,减小对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性。由于本发明虑及裂纹萌生-扩展及残余应力,可以提高计算准确度,降低试验成本低。
本领域内的技术人员应明白,本发明的实施例可提供为方法、***、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种齿轮接触疲劳全寿命评估方法,其特征在于,包括:
根据最大接触应力模型计算齿轮接触面上的最大接触应力;
基于数值计算理论及等效边界条件,分别构建齿轮的二维静态模型、三维静态模型、二维动态模型及三维动态模型,并基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型分别获得对应的齿轮接触面上的最大接触应力;
将基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型获得的对应最大接触应力分别与基于最大接触应力模型计算得到的最大接触应力进行比较,确定最佳数值计算模型;
基于位错-能量法及裂纹尺寸构建齿轮接触疲劳萌生寿命评估模型;
在半球空间内,基于裂纹扩展中的应力分布及滑开型裂纹确定有效应力强度因子值范围;
基于Paris方程、裂纹扩展角度、裂纹增量尺寸、齿轮材料硬度及所述有效应力强度因子值范围,构建齿轮接触疲劳扩展寿命评估模型;
根据所述齿轮接触疲劳萌生寿命评估模型及所述齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估。
2.根据权利要求1所述的齿轮接触疲劳全寿命评估方法,其特征在于,所述齿轮接触疲劳萌生寿命评估模型为:
其中,Ni为接触疲劳萌生寿命,l为椭圆滑移带的半长轴,G为剪切模量,Δτ为剪切应力范围,k为位错摩擦应力,c为裂纹长度,d为晶粒尺寸。
3.根据权利要求2所述的齿轮接触疲劳全寿命评估方法,其特征在于,还包括:基于所述最佳数值计算模型对所述齿轮接触疲劳萌生寿命评估模型的拟合系数及椭圆滑移带的半长轴进行参数修正。
4.根据权利要求3所述的齿轮接触疲劳全寿命评估方法,其特征在于,还包括:
基于渗碳层中残余应力对齿轮接触疲劳萌生寿命的影响,对参数修正后的所述齿轮接触疲劳萌生寿命评估模型进行修正;
根据所述齿轮接触疲劳萌生寿命评估模型及所述齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估,包括:
根据所述齿轮接触疲劳扩展寿命评估模型及修正后的所述齿轮接触疲劳萌生寿命评估模型进行齿轮接触疲劳全寿命评估。
5.根据权利要求4所述的齿轮接触疲劳全寿命评估方法,其特征在于,还包括:
基于刚塑性滑移线理论,结合每次循环载荷作用后最大剪切应力,确定裂纹扩展角度;
基于裂纹尖端塑性区域确定所述裂纹增量尺寸。
6.根据权利要求2所述的齿轮接触疲劳全寿命评估方法,其特征在于,所述齿轮接触疲劳扩展寿命评估模型为:
其中,Np为齿轮接触疲劳扩展寿命,C和m分别为材料参数,L为裂纹尺寸,xL为裂纹尖端在X轴方向上的坐标,ξ为裂纹扩展增量,τmax为齿轮次表面中最大剪切应力,Hb及HL分别为齿轮整体硬度及局部硬度,a0为裂纹初始长度。
7.一种齿轮接触疲劳全寿命评估装置,其特征在于,包括:
第一最大接触应力计算单元,用于根据最大接触应力模型计算齿轮接触面上的最大接触应力;
第二最大接触应力计算单元,用于基于数值计算理论及等效边界条件,分别构建齿轮的二维静态模型、三维静态模型、二维动态模型及三维动态模型,并基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型分别获得对应的齿轮接触面上的最大接触应力;
最佳数值计算模型确定单元,用于将基于所述二维静态模型、三维静态模型、二维动态模型及三维动态模型获得的对应最大接触应力分别与基于最大接触应力模型计算得到的最大接触应力进行比较,确定最佳数值计算模型;
齿轮接触疲劳萌生寿命评估模型构建单元,用于基于位错-能量法及裂纹尺寸构建齿轮接触疲劳萌生寿命评估模型;
有效应力强度因子值范围确定单元,用于在半球空间内,基于裂纹扩展中的应力分布及滑开型裂纹确定有效应力强度因子值范围;
齿轮接触疲劳扩展寿命评估模型构建单元,用于基于Paris方程、裂纹扩展角度、裂纹增量尺寸、齿轮材料硬度及所述有效应力强度因子值范围,构建齿轮接触疲劳扩展寿命评估模型;
齿轮接触疲劳全寿命评估单元,用于根据所述齿轮接触疲劳萌生寿命评估模型及所述齿轮接触疲劳扩展寿命评估模型进行齿轮接触疲劳全寿命评估。
8.根据权利要求7所述的齿轮接触疲劳全寿命评估装置,其特征在于,还包括:参数修正单元,用于基于所述最佳数值计算模型对所述齿轮接触疲劳萌生寿命评估模型的拟合系数及椭圆滑移带的半长轴进行参数修正。
9.根据权利要求8所述的齿轮接触疲劳全寿命评估装置,其特征在于,还包括:
模型修正单元,用于基于渗碳层中残余应力对齿轮接触疲劳萌生寿命的影响,对参数修正后的所述齿轮接触疲劳萌生寿命评估模型进行修正;
所述齿轮接触疲劳全寿命评估单元具体用于:
根据所述齿轮接触疲劳扩展寿命评估模型及修正后的所述齿轮接触疲劳萌生寿命评估模型进行齿轮接触疲劳全寿命评估。
10.根据权利要求9所述的齿轮接触疲劳全寿命评估装置,其特征在于,还包括:
裂纹扩展角度确定单元,用于基于刚塑性滑移线理论,结合每次循环载荷作用后最大剪切应力,确定裂纹扩展角度;
裂纹增量尺寸确定单元,用于基于裂纹尖端塑性区域确定所述裂纹增量尺寸。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710200764.9A CN106979861B (zh) | 2017-03-30 | 2017-03-30 | 齿轮接触疲劳全寿命评估方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710200764.9A CN106979861B (zh) | 2017-03-30 | 2017-03-30 | 齿轮接触疲劳全寿命评估方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106979861A CN106979861A (zh) | 2017-07-25 |
CN106979861B true CN106979861B (zh) | 2019-04-23 |
Family
ID=59339072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710200764.9A Expired - Fee Related CN106979861B (zh) | 2017-03-30 | 2017-03-30 | 齿轮接触疲劳全寿命评估方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106979861B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107885907A (zh) * | 2017-10-17 | 2018-04-06 | 重庆大学 | 一种评估表面硬化齿轮接触疲劳失效风险的方法 |
CN108444697B (zh) * | 2017-12-04 | 2020-04-07 | 重庆大学 | 渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法 |
CN108256241A (zh) * | 2018-01-23 | 2018-07-06 | 重庆大学 | 一种重载齿轮次表面裂纹萌生的预测方法 |
CN108318244B (zh) * | 2018-01-23 | 2019-12-31 | 重庆大学 | 考虑残余应力的渗碳硬化齿轮接触疲劳风险评估方法 |
CN109165425B (zh) * | 2018-08-03 | 2022-04-12 | 湖南大学 | 一种齿轮接触疲劳可靠性分析方法 |
CN109271711B (zh) * | 2018-09-25 | 2023-03-28 | 重庆大学 | 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法 |
CN110147624B (zh) * | 2019-05-24 | 2023-05-09 | 重庆大学 | 一种基于载荷谱的齿轮接触疲劳寿命预测方法 |
CN111122358B (zh) * | 2020-01-13 | 2022-05-31 | 上海工程技术大学 | 一种考虑滞弹性能的镁合金疲劳寿命的确定方法 |
CN111751199B (zh) * | 2020-06-09 | 2021-07-09 | 西北工业大学 | 基于eifs分布的疲劳寿命预测方法 |
CN113514351B (zh) * | 2020-12-23 | 2024-04-19 | 中国特种设备检测研究院 | 一种考虑预应力重分布的疲劳裂纹扩展行为预测方法 |
CN113094957B (zh) * | 2021-04-15 | 2023-03-31 | 中南大学 | 一种评估rv减速器曲柄轴接触疲劳寿命的方法 |
CN113591268B (zh) * | 2021-06-24 | 2023-06-23 | 内蒙古工业大学 | 变幅载荷下齿轮接触疲劳寿命可靠性评估方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1955043A1 (en) * | 2005-12-02 | 2008-08-13 | United Technologies Corporation | Systems and methods for modeling surface properties of a mechanical component |
CN101839904A (zh) * | 2009-03-12 | 2010-09-22 | 通用汽车环球科技运作公司 | 预测铝合金在多轴加载下的疲劳寿命的***和方法 |
CN103616179A (zh) * | 2013-12-05 | 2014-03-05 | 广西大学 | 一种基于缺陷建模的传动齿轮疲劳寿命评估方法 |
CN105466772A (zh) * | 2015-11-15 | 2016-04-06 | 北京工业大学 | 一种基于临界面法的多轴短裂纹扩展寿命预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6612181B2 (en) * | 2001-09-04 | 2003-09-02 | Jalees Ahmad | Method and system for determining crack nucleation of a part subject to fretting fatigue |
-
2017
- 2017-03-30 CN CN201710200764.9A patent/CN106979861B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1955043A1 (en) * | 2005-12-02 | 2008-08-13 | United Technologies Corporation | Systems and methods for modeling surface properties of a mechanical component |
CN101839904A (zh) * | 2009-03-12 | 2010-09-22 | 通用汽车环球科技运作公司 | 预测铝合金在多轴加载下的疲劳寿命的***和方法 |
CN103616179A (zh) * | 2013-12-05 | 2014-03-05 | 广西大学 | 一种基于缺陷建模的传动齿轮疲劳寿命评估方法 |
CN105466772A (zh) * | 2015-11-15 | 2016-04-06 | 北京工业大学 | 一种基于临界面法的多轴短裂纹扩展寿命预测方法 |
Non-Patent Citations (1)
Title |
---|
油润滑渐开线斜齿轮摩擦动力学特性及疲劳寿命预估;董辉立;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20140615;第C029-16页 |
Also Published As
Publication number | Publication date |
---|---|
CN106979861A (zh) | 2017-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106979861B (zh) | 齿轮接触疲劳全寿命评估方法及装置 | |
CN109165425B (zh) | 一种齿轮接触疲劳可靠性分析方法 | |
Mohammed et al. | Improving mesh stiffness calculation of cracked gears for the purpose of vibration-based fault analysis | |
US20190197211A1 (en) | Method for estimating stress intensity factors and method for calculating associated service life | |
Lewicki et al. | Consideration of moving tooth load in gear crack propagation predictions | |
CN106202647A (zh) | 电主轴的多轴疲劳寿命预测方法及疲劳寿命可靠性评估方法 | |
Liu et al. | Crack modeling of rotating blades with cracked hexahedral finite element method | |
CN105008887A (zh) | 使用非接触测量和动态响应重构技术的涡轮机叶片疲劳寿命分析 | |
Raptis et al. | Rating of spur gear strength using photoelasticity and the finite element method | |
US9081921B2 (en) | Method for simulating rubber compound | |
CN105260574A (zh) | 一种基于临界面法疲劳破坏准则的高周多轴疲劳寿命预测方法 | |
Rama et al. | Investigation of the effect of load distribution along the face width and load sharing between the pairs in contact on the fracture parameters of the spur gear tooth with root crack | |
CN107180123B (zh) | 一种高强度钢潜水器耐压球壳极限承载力估算方法 | |
Podrug et al. | Numerical modelling of crack growth in a gear tooth root | |
Yin et al. | Lifetime prediction for the subsurface crack propagation using three-dimensional dynamic FEA model | |
Hiung et al. | Implementation of XFEM in the study of gear crack propagation behaviour using the SIF on different moments | |
Zhao | Increasing bending strength in spur gears using shape optimisation of cutting tool profile | |
CN107066728B (zh) | 一种钛合金潜水器耐压球壳极限承载力估算方法 | |
Behzad et al. | A new continuous model for flexural vibration analysis of a cracked beam | |
Roberts et al. | SGBEM modeling of fatigue crack growth in particulate composites | |
Guagliano et al. | Mode I stress intensity factors for curved cracks in gears by a weight functions method | |
Nazari et al. | Double cracks identification in functionally graded beams using artificial neural network | |
Dotti et al. | Vibration of open thin-walled beams with a breathing crack | |
Yu | On the fatigue crack growth analysis of spliced aircraft wing panels under sequential axial and shear loads | |
Stimpfle et al. | A project oriented approach to determine membrane properties |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190423 |
|
CF01 | Termination of patent right due to non-payment of annual fee |