CN113704907A - 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 - Google Patents
一种基于齿面应力场的齿轮接触疲劳寿命预测方法 Download PDFInfo
- Publication number
- CN113704907A CN113704907A CN202110921947.6A CN202110921947A CN113704907A CN 113704907 A CN113704907 A CN 113704907A CN 202110921947 A CN202110921947 A CN 202110921947A CN 113704907 A CN113704907 A CN 113704907A
- Authority
- CN
- China
- Prior art keywords
- tooth surface
- gear
- stress
- node
- stress field
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 239000002131 composite material Substances 0.000 claims abstract description 22
- 239000000463 material Substances 0.000 claims abstract description 22
- 238000012360 testing method Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000005452 bending Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000009661 fatigue test Methods 0.000 claims description 2
- 230000002706 hydrostatic effect Effects 0.000 claims description 2
- 230000006872 improvement Effects 0.000 description 4
- 229920001971 elastomer Polymers 0.000 description 2
- 239000000806 elastomer Substances 0.000 description 2
- 238000002441 X-ray diffraction Methods 0.000 description 1
- 244000309464 bull Species 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/68—Analysis of geometric attributes of symmetry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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/04—Ageing analysis or optimisation against ageing
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
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)
- Computer Vision & Pattern Recognition (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于齿面应力场的齿轮接触疲劳寿命预测方法,包括如下步骤:步骤一、将齿面离散成多个齿面节点并且构建齿面复合应力场;步骤二、基于Papadopoulos多轴疲劳准则,结合齿轮材料疲劳特性参数,计算应力场下齿面各齿面节点的多轴等效应力:步骤三、结合齿轮材料S‑N曲线,根据各齿面节点的多轴等效应力预测齿轮接触疲劳寿命。本发明基于Papadopoulos多轴疲劳准则,可更为准确地预测齿面应力场下的接触疲劳寿命,为评估与提升齿轮使用寿命提供指导,且本发明考虑了齿面残余应力与接触应力的复合应力场,使得其更接近于齿轮的实际应力状态,从而对于齿轮的疲劳寿命预测更加准确。
Description
技术领域
本发明属于齿轮评估技术领域,尤其涉及一种基于齿面应力场的齿轮接触疲劳寿命预测方法。
背景技术
目前对于齿轮的寿命预测,通常是基于ISO与AGMA等齿轮标准中寿命预测方法,其均基于累积损伤准则与大量试验,成本较高且对复杂载荷下齿面接触疲劳寿命预测精度较低。
此外,现有的齿轮进行寿命预测前需要分析其应力场。但是目前ISO与AGMA的应力场分析通常是基于Hertz接触理论,而Hertz接触理论通过把两弹性体的接触问题等效为一个当量弹性体和一个刚性平面的接触,对复杂载荷下齿面接触应力计算精度较低,工程实际中没有对齿轮表面与亚表面的残余应力场进行考虑,不能准确反映齿面残余应力与接触应力的复合应力状态。因此得到的应力场与实际情况差别较大。这也导致对齿轮寿命的预测不够精确。
Papadopoulos多轴疲劳准则见参考文献:Papadopoulos IV.Long life fatigueunder multiaxial loading.International Journal of Fatigue.2001;23。
发明内容
为解决上述问题,本发明公开一种基于齿面应力场的齿轮接触疲劳寿命预测方法,本发明基于Papadopoulos多轴疲劳准则,可更为准确地预测齿面应力场下的接触疲劳寿命,为评估与提升齿轮使用寿命提供指导,且本发明考虑了齿面残余应力与接触应力的复合应力场,使得其更接近于齿轮的实际应力状态,从而对于齿轮的疲劳寿命预测更加准确。
为实现上述目的,本发明的技术方案如下:
一种基于齿面应力场的齿轮接触疲劳寿命预测方法,包括以下步骤:
步骤一、将齿面离散成多个齿面节点并且构建齿面复合应力场;
步骤二、基于Papadopoulos多轴疲劳准则,结合齿轮材料疲劳特性参数,计算复合应力场下各齿面节点的多轴等效应力:
式中,σn为第n个齿面节点的多轴等效应力,τ-1为齿轮材料的弯曲疲劳极限,σ-1为齿轮材料的扭转疲劳极限,Tn为第n个齿面节点的广义剪切应力幅值,σn1与σn2分别第n个齿面节点为静水压力的幅值与均值;
步骤三、结合齿轮材料S-N曲线,根据各齿面节点的多轴等效应力预测齿轮接触疲劳寿命:
a·(2Nn )b=σn (1.2)
其中,Nn为第n个齿面节点的疲劳寿命,σn为第n个齿面节点的多轴等效应力,a为疲劳强度系数,b为疲劳强度指数;a与b由齿轮材料S-N曲线得到;其中最小的Nn即为齿轮的使用寿命。
进一步的改进,所述齿轮材料S-N曲线得到方法如下:通过多个不同应力水平下的疲劳试验,以测得的疲劳寿命为横坐标,以应力为纵坐标,拟合各数据点得到齿轮材料S-N曲线。
进一步的改进,所述齿轮材料疲劳特性参数包括齿轮材料的弯曲疲劳极限τ-1和扭转疲劳极限σ-1。
进一步的改进,所述齿面复合应力场为基于齿面残余应力与接触应力的复合应力场。
进一步的改进,所述复合应力场的建立方法如下:
步骤1.1建立齿轮三维几何模型;然后输入齿轮材料参数,导入齿轮三维几何模型,对齿面接触区域进行网格细化,构建齿轮三维接触分析模型;
步骤1.2、测试得到齿轮试件沿齿厚方向各测量点的残余应力,拟合得到残余应力分布曲线;基于齿轮三维几何模型,将齿面离散成若干齿面节点,根据残余应力分布曲线,得到在齿面坐标系下各齿面节点的残余应力;
步骤1.3、结合各齿面节点的方向余弦矩阵,将齿面坐标系下的残余应力转换到全局坐标系中,得到在全局坐标系下的齿面节点的残余应力σ′n0:
σ′n0=[βn]-1[σn0]([βn]T)-1 (1.3)
其中,σ′n0为全局坐标系下第n个齿面节点的残余应力,σn0为齿面坐标系下第n个齿面节点的残余应力,βn为齿面节点的方向余弦,n表示节点编号;
步骤1.4、将全局坐标系下的各吃点节点的残余应力数据导入ABAQUS平台,耦合齿面接触应力场与残余应力场,构建齿面复合应力场。
本发明的优点:
本发明基于Papadopoulos多轴疲劳准则,可更为准确地预测齿面应力场下的接触疲劳寿命,为评估与提升齿轮使用寿命提供指导,且本发明考虑了齿面残余应力与接触应力的复合应力场,使得其更接近于齿轮的实际应力状态,从而对于齿轮的疲劳寿命预测更加准确。
附图说明
图1齿轮三维有限元接触模型;
图2齿面残余应力分布曲线;
图3齿面复合应力场;
图4距齿面不同深度处沿齿宽方向多轴等效应力;
图5距齿面不同深度处沿齿宽方向疲劳寿命。
具体实施方式
下面结合附图以及实施例对本发明进行更为详细的说明。
步骤1.1、输入齿轮几何参数见表1,采用CATIA平台建立齿槽曲面,并通过旋转阵列及布尔操作构建齿轮三维几何模型。
表1齿轮几何参数
将几何模型导入ABAQUS平台,定义齿轮材料属性见表2。对不同齿轮区域分别进行结构化网格划分,网格单元类型选择线性缩减积分单元C3D8R。为提高计算精度,对齿面接触表面与次表面区域进行细化网格。
表2齿轮材料属性
设置转速为2000r/min,转矩为1200N·m,且将转速与转矩设置在不同的分析步中。首先固定小轮所有自由度,在大齿轮轴线自由度施加较小的转动位移,使小轮齿面与大轮齿面接触。然后释放大轮绕其轴线的旋转自由度,并施加转矩;最后在小轮轴线上施加一定的转速,构建齿轮三维有限元接触模型见图1。
步骤1.2、采用X射线衍射残余应力测试与电解腐蚀剥层,测量经渗碳淬火-磨削-喷丸工艺齿轮试件表面与次表面的残余应力,得到距齿面不同深度的残余应力分布曲线如图2所示。将齿面离散成多个齿面节点,根据测试数据插值拟合各节点的残余应力,构建齿面残余应力场。
步骤1.3、结合各齿面节点的方向余弦矩阵,通过式(1.3)将齿面坐标系下的残余应力转换到全局坐标系中,得在全局坐标系下的齿面上各节点的残余应力。
将全局坐标系下的齿面残余应力数据导入ABAQUS平台,耦合齿面接触应力场与残余应力场,构建齿面复合应力场如图3所示。输出齿面各节点复合应力数据,可得各齿面节点同时存在正应力和切应力,且正应力与切应力峰值不同时到达,为多轴非比例加载。
步骤二、基于Papadopoulos多轴疲劳准则,结合齿轮材料疲劳特性参数,通过式计算复合应力场下齿面各节点的多轴等效应力,得到距齿面不同深度处沿齿宽方向多轴等效应力分布曲线见图4。
步骤三、结合齿面节点的多轴等效应力与齿轮材料S-N曲线,通过式计算齿轮接触疲劳寿命,得到距齿面不同深度处沿齿宽方向疲劳寿命分布曲线见图5。
上述仅为本发明的一个具体导向实施方式,但本发明的设计构思并不局限于此,凡利用此构思对本发明进行非实质性的改动,均应属于侵犯本发明的保护范围的行为。
Claims (5)
1.一种基于齿面应力场的齿轮接触疲劳寿命预测方法,其特征在于,包括以下步骤:
步骤一、将齿面离散成多个齿面节点并且构建齿面复合应力场;
步骤二、基于Papadopoulos多轴疲劳准则,结合齿轮材料疲劳特性参数,计算复合应力场下各齿面节点的多轴等效应力:
式中,σn为第n个齿面节点的多轴等效应力,τ-1为齿轮材料的弯曲疲劳极限,σ-1为齿轮材料的扭转疲劳极限,Tn为第n个齿面节点的广义剪切应力幅值,σn1与σn2分别第n个齿面节点为静水压力的幅值与均值;
步骤三、结合齿轮材料S-N曲线,根据各齿面节点的多轴等效应力预测齿轮接触疲劳寿命:
a·(2Nn)b=σn (1.2)
其中,Nn为第n个齿面节点的疲劳寿命,σn为第n个齿面节点的多轴等效应力,a为疲劳强度系数,b为疲劳强度指数;a与b由齿轮材料S-N曲线得到;其中最小的Nn即为齿轮的使用寿命。
2.如权利要求1所述的基于齿面应力场的齿轮接触疲劳寿命预测方法,其特征在于,所述齿轮材料S-N曲线得到方法如下:通过多个不同应力水平下的疲劳试验,以测得的疲劳寿命为横坐标,以应力为纵坐标,拟合各数据点得到齿轮材料S-N曲线。
3.如权利要求1所述的基于齿面应力场的齿轮接触疲劳寿命预测方法,其特征在于,所述齿轮材料疲劳特性参数包括齿轮材料的弯曲疲劳极限τ-1和扭转疲劳极限σ-1。
4.如权利要求1所述的基于齿面应力场的齿轮接触疲劳寿命预测方法,其特征在于,所述齿面复合应力场为基于齿面残余应力与接触应力的复合应力场。
5.如权利要求4所述的基于齿面应力场的齿轮接触疲劳寿命预测方法,其特征在于,所述齿面复合应力场的建立方法如下:
步骤1.1建立齿轮三维几何模型;然后输入齿轮材料参数,导入齿轮三维几何模型,对齿面接触区域进行网格细化,构建齿轮三维接触分析模型;
步骤1.2、测试得到齿轮试件沿齿厚方向各测量点的残余应力,拟合得到残余应力分布曲线;基于齿轮三维几何模型,将齿面离散成若干齿面节点,根据残余应力分布曲线,得到在齿面坐标系下各齿面节点的残余应力;
步骤1.3、结合各齿面节点的方向余弦矩阵,将齿面坐标系下的残余应力转换到全局坐标系中,得到在全局坐标系下的齿面节点的残余应力σ′n0:
σ′n0=[βn]-1[σn0]([βn]T)-1 (1.1)
其中,σ′n0为全局坐标系下第n个齿面节点的残余应力,σn0为齿面坐标系下第n个齿面节点的残余应力,βn为齿面节点的方向余弦,n表示节点编号;
步骤1.4、将全局坐标系下的各吃点节点的残余应力数据导入ABAQUS平台,耦合齿面接触应力场与残余应力场,构建齿面复合应力场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110921947.6A CN113704907A (zh) | 2021-08-13 | 2021-08-13 | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110921947.6A CN113704907A (zh) | 2021-08-13 | 2021-08-13 | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113704907A true CN113704907A (zh) | 2021-11-26 |
Family
ID=78652391
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110921947.6A Pending CN113704907A (zh) | 2021-08-13 | 2021-08-13 | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113704907A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114754913A (zh) * | 2022-04-15 | 2022-07-15 | 浙江大学 | 一种基于弯扭信息融合的寿命预测方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101416151A (zh) * | 2006-03-08 | 2009-04-22 | 动力学专家有限公司 | 可靠性仿真方法和系统 |
KR101576007B1 (ko) * | 2015-09-10 | 2015-12-11 | 충남대학교산학협력단 | 트랙터 pto축 어댑터를 이용한 pto기어 잔여수명 예측장치 및 그 방법 |
CN106202647A (zh) * | 2016-06-29 | 2016-12-07 | 北京科技大学 | 电主轴的多轴疲劳寿命预测方法及疲劳寿命可靠性评估方法 |
CN106886663A (zh) * | 2017-03-29 | 2017-06-23 | 北京理工大学 | 齿轮弯曲疲劳寿命预测方法及装置 |
CN110147624A (zh) * | 2019-05-24 | 2019-08-20 | 重庆大学 | 一种基于载荷谱的齿轮接触疲劳寿命预测方法 |
-
2021
- 2021-08-13 CN CN202110921947.6A patent/CN113704907A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101416151A (zh) * | 2006-03-08 | 2009-04-22 | 动力学专家有限公司 | 可靠性仿真方法和系统 |
KR101576007B1 (ko) * | 2015-09-10 | 2015-12-11 | 충남대학교산학협력단 | 트랙터 pto축 어댑터를 이용한 pto기어 잔여수명 예측장치 및 그 방법 |
CN106202647A (zh) * | 2016-06-29 | 2016-12-07 | 北京科技大学 | 电主轴的多轴疲劳寿命预测方法及疲劳寿命可靠性评估方法 |
CN106886663A (zh) * | 2017-03-29 | 2017-06-23 | 北京理工大学 | 齿轮弯曲疲劳寿命预测方法及装置 |
CN110147624A (zh) * | 2019-05-24 | 2019-08-20 | 重庆大学 | 一种基于载荷谱的齿轮接触疲劳寿命预测方法 |
Non-Patent Citations (2)
Title |
---|
张 帆等: "渐开线直齿轮弹流润滑条件下的多轴疲劳寿命预估", 摩 擦 学 学 报, vol. 37, no. 2, pages 263 - 269 * |
郑晓猛 等: "基于Matlab 计算滚动轴承滚滑接触内部应力分布", 润滑与密封, vol. 45, no. 7, pages 52 - 61 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114754913A (zh) * | 2022-04-15 | 2022-07-15 | 浙江大学 | 一种基于弯扭信息融合的寿命预测方法及装置 |
CN114754913B (zh) * | 2022-04-15 | 2023-02-28 | 浙江大学 | 一种基于弯扭信息融合的寿命预测方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109165425B (zh) | 一种齿轮接触疲劳可靠性分析方法 | |
CN106096073B (zh) | 一种基于损伤力学非概率区间分析模型的金属疲劳裂纹全寿命预估方法 | |
CN112765748B (zh) | 机理-数据异构信息融合的旋转机械数字孪生建模方法 | |
Wei | Stresses and deformations in involute spur gears by finite element method | |
CN105956368B (zh) | 一种复合材料层合结构的非概率动力可靠性评估方法 | |
CN107506535A (zh) | 一种基于临界应变损伤参量的多轴疲劳寿命预测方法 | |
CN107977516B (zh) | 一种考虑多轴载荷非比例度的缺口件局部应力应变确定方法 | |
Podrug et al. | Influence of different load models on gear crack path shapes and fatigue lives | |
CN109271711B (zh) | 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法 | |
Wang et al. | Determination of the minimum sample size for the transmission load of a wheel loader based on multi-criteria decision-making technology | |
CN113704907A (zh) | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 | |
Li et al. | Compilation of load spectrum of machining center spindle and application in fatigue life prediction | |
CN106907534A (zh) | 一种多层非粘结柔性管抗拉性能快速评价方法 | |
Gao et al. | Reliability and availability models of belt drive systems considering failure dependence | |
CN108491640B (zh) | 一种多轴疲劳寿命预测模型 | |
CN114254533B (zh) | 考核疲劳振动对产品组部件固定角度影响和预测的方法 | |
CN106909709A (zh) | 一种多层非粘结柔性管抗弯性能快速评价方法 | |
CN115238394B (zh) | 复合材料结构多源不确定性混合可靠性数字孪生建模方法 | |
CN108444697B (zh) | 渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法 | |
Svensson | Fatigue prediction models of Dynamic Power Cables by laboratory testing and FE analysis | |
CN113063589B (zh) | 一种基于神经网络的齿轮微观误差对振动的预测方法 | |
Iorga et al. | Engineering design from a safety perspective | |
CN115070774A (zh) | 一种基于中值神经网络的工业机器人故障模式识别方法 | |
Wang et al. | Establishment method of a mixture model and its practical application for transmission gears in an engineering vehicle | |
CN112651153A (zh) | 一种确定晶体塑性有限元模型材料参数的方法 |
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 |