CN106649971B - 一种弧齿锥齿轮长寿命传动疲劳可靠性的评价方法 - Google Patents

一种弧齿锥齿轮长寿命传动疲劳可靠性的评价方法 Download PDF

Info

Publication number
CN106649971B
CN106649971B CN201610937468.2A CN201610937468A CN106649971B CN 106649971 B CN106649971 B CN 106649971B CN 201610937468 A CN201610937468 A CN 201610937468A CN 106649971 B CN106649971 B CN 106649971B
Authority
CN
China
Prior art keywords
variable
grinding
sensitivity
tooth surface
value
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
Application number
CN201610937468.2A
Other languages
English (en)
Other versions
CN106649971A (zh
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.)
Shenyang Dechuang metal materials Co.,Ltd.
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN201610937468.2A priority Critical patent/CN106649971B/zh
Publication of CN106649971A publication Critical patent/CN106649971A/zh
Application granted granted Critical
Publication of CN106649971B publication Critical patent/CN106649971B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing

Landscapes

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

Abstract

本发明涉及一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,包括:(1)研究磨削工艺参数对齿面粗糙度和残余应力的灵敏度的影响;(2)研究热处理工艺参数对齿面硬度和渗碳深度的灵敏度的影响;(3)研究齿面粗糙度、残余应力、齿面硬度和渗碳深度对弧齿锥齿疲劳长寿命可靠性的灵敏度;(4)计算得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度,通过灵敏度来评价磨削和热处理。本发明可以正确的评价弧齿锥齿轮疲劳可靠性,优化加工参数;同时有效地控制了磨削残余应力,并提高了齿面质量,降低了齿面粗糙度,提高了生产效率,对于齿轮的实际生产有着重要的意义。

Description

一种弧齿锥齿轮长寿命传动疲劳可靠性的评价方法
技术领域
本发明涉及一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,适用于复杂工况下的弧齿锥齿轮长寿命传动可靠性的的评价以及对磨削和热处理工艺的优化。
背景技术
航空弧齿锥齿轮是直升机主减传动系统最重要、最复杂的、最薄弱的动力元件,不同于其它常用的直齿轮等。常处于高速重载高温、有体积重量限制的条件下工作,其工况恶劣,载荷谱变化复杂,疲劳应力环境和失效模式复杂。目前航空弧齿锥齿轮长寿命可靠性设计与评估技术已经成为制约航空传动系统向高可靠、长寿命发展瓶颈,也成为规范齿轮抗疲劳设计、优化加工工艺评估的关键。针对国内外在发动机高载荷质量比弧齿锥齿轮疲劳可靠性设计与评估技术方面的差距,在现有疲劳可靠性理论基础上,将可靠性概率理论与实际寿命计算结合起来,将疲劳可靠性能通过可靠性模型与齿轮的设计和加工参数联系起来,建立合理有效的疲劳寿命可靠性分析评估方法,指导具有疲劳可靠性的要求的齿轮的设计和工艺安排,从而为航空弧齿锥齿轮轻量化长寿命可靠性设计制造提供保障。
发明内容
本发明要解决的技术问题为:针对给定复杂工况下的弧齿锥齿轮长寿命疲劳可靠性的设计,以齿面粗糙度、残余应力、齿面硬度和渗碳深度等性能的为中间变量,综合考虑弧齿锥齿轮在实际加工过程中的磨削和热处理工艺参数对齿面粗糙度、残余应力、齿面硬度和渗碳深度等性能的灵敏度,以及齿面粗糙度、残余应力、齿面硬度和渗碳深度等性能的对可靠性的灵敏度,为齿轮传动系统的寿命设计和评估提供了一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,有效降低了长寿命传动疲劳可靠性的评价所带来的难度和成本。
本发明采用的技术方案是:一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其方法流程如下:
步骤(1)、研究磨削工艺参数对齿面粗糙度和残余应力的灵敏度的影响;
步骤(2)、研究热处理工艺参数对齿面硬度和渗碳深度的灵敏度的影响;
步骤(3)、研究齿面粗糙度、残余应力、齿面硬度和渗碳深度对弧齿锥齿疲劳长寿命可靠性的灵敏度;
步骤(4)、计算得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度,通过灵敏度来评价磨削和热处理。
进一步的,所述步骤(1)中齿轮磨削时,齿轮材料、磨削液因素是一定的,磨削深度ap、磨削速度Vs和横向进给速度Vw是变化的,选择这3个因素作为正交试验因素,根据磨削工艺经验推荐值,每因素分别取若干个不同水平值,建立正交试验表。开展磨削试验,记录每组磨削试验后齿轮的表面粗糙度Ra、残余应力S11。
进一步的,所述步骤(1)中用有限差分法对表面粗糙度和残余应力进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δx,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式。分别计算出磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差的灵敏度,
Figure GDA0002199192450000021
式中:βx表示y变量对x变量的灵敏度,x变量表示磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差,y变量表示表面粗糙度Ra、残余应力S11;x0为起始的x变量的值,y0为x0对应的值,x1为微小变化之后的x变量的值,y1为x1对应的值。
进一步的,所述步骤(2)中齿轮热处理时,齿轮材料、模数因素一定,渗碳时间和渗碳温度是变化的,选择这2个因素作为正交试验因素,根据热处理工艺经验推荐值,没因素分别取若干个不同的水平值,建立正交实验表。开展热处理试验,记录每组热处理试验后的齿面硬度和渗碳深度。
进一步的,所述步骤(2)中用有限差分法对齿面硬度和渗碳深度进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δa,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式。分别计算出渗碳时间和渗碳温度两个参数的均值以及方差的灵敏度,
Figure GDA0002199192450000022
式中:βa表示z变量对a变量的灵敏度,a变量表示渗碳时间和渗碳温度两个参数的均值以及方差,z变量表示齿面硬度和渗碳深度;a0为起始的a变量的值,z0为a0对应的值,a1为微小变化之后的a变量的值,z1为a1对应的值。
进一步的,所述步骤(3)中航空弧齿锥齿轮齿面接触疲劳强度公式为:
σ′Hlim=σHlimZNZLZVZRZWZx
式中:σ`Hlim表示齿面接触疲劳强度,σHlimZN表示接触应力,ZL表示润滑剂系数,Zv表示速度系数,ZR表示粗糙度系数,ZW表示工作硬化系数,ZX表示尺寸系数。
齿根弯曲疲劳强度公式为:
σ′Flim=σFlimYNYSTYσYRYx
式中:σ`Flim表示齿根弯曲疲劳强度,σFlimYN表示弯曲应力,YST表示应力修正系数,Yσ表示相对齿根圆角敏感系数,YR表示相对齿根表面状况系数,Yx表示弯曲强度尺寸系数;
齿面接触应力公式为:
Figure GDA0002199192450000031
式中:σH表示齿面接触应力,ZM-B表示中点区域系数,ZH表示节点啮合区域系数,ZE表示弹性影响系数,ZLS表示载荷分担系数,Zβ表示螺旋角系数,ZK表示锥齿轮系数,Fmt表示切向力,KA表示使用系数,KV表示动载系数,K表示齿向载荷分布系数,K表示端面载荷分配系数,dv1表示小齿轮的分度圆直径,lbm表示接触线长度,uv表示齿数比。
齿根弯曲应力公式为:
Figure GDA0002199192450000032
式中:σF表示弯曲应力,Fmt表示切向力,KA表示使用系数,KV表示动载系数,K表示齿向载荷分布系数,K表示端面载荷分配系数,YFa表示齿形系数,YSa表示应力修正系数,Yε表示重合度系数,YK表示锥齿轮系数,YLS表示载荷分担系数,b表示工作齿宽,mmn表示小齿轮的法向模数。
表面粗糙度、齿面硬度、精度等轮齿表面质量参数直接影响上述公式中的粗糙度系数ZR、工作硬化系数ZW、齿根表面状况系数YR以及可靠性分析中某些参数的方差分布的选取,所以会对齿轮疲劳强度和可靠性造成影响。
进一步的,所述步骤(3)中研究齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数对齿轮疲劳可靠性的敏感性影响采用有限差分法对弧齿锥齿轮进行灵敏度分析。其基本做法是使设计变量有一个微小的摄动,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,
Figure GDA0002199192450000041
式中βy表示R变量对y变量的灵敏度,y变量表示表面粗糙度Ra、残余应力S11两个表面参数的均值以及方差,R变量表示可靠性;y0为起始的y变量的值,R0为y0对应的值,y1为微小变化之后的y变量的值,R1为y1对应的值,
Figure GDA0002199192450000042
式中:βz表示R变量对z变量的灵敏度,z变量表示齿面硬度和渗碳深度两个表面参数的均值以及方差,R变量表示可靠性;z0为起始的z变量的值,R0为z0对应的值,z1为微小变化之后的z变量的值,R1为z1对应的值。
进一步的,所述步骤(4)中根据求导规则,将齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数当做中间变量消掉,最终得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度:
Figure GDA0002199192450000043
Figure GDA0002199192450000044
式中:βx(R)表示R变量对x变量的灵敏度,x变量表示磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差,y变量表示表面粗糙度Ra、残余应力S11;x0为起始的x变量的值,y0为x0对应的值,x1为微小变化之后的x变量的值,y1为x1对应的值;
式中:βa(R)表示R变量对a变量的灵敏度,a变量表示渗碳时间和渗碳温度两个参数的均值以及方差,z变量表示齿面硬度和渗碳深度;a0为起始的a变量的值,z0为a0对应的值,a1为微小变化之后的a变量的值,z1为a1对应的值,
式中参数同以上所述一样。
最终通过灵敏度来评价磨削和热处理。
本发明的原理:基于正交试验和使设计变量有一个微小的摄动,用差分格式计算输出,得到不同工艺参数对于可靠性的灵敏度的影响权重,对可靠性进行评价,并对工艺参数进行优化。
本发明与现有技术相比的有益效果是:首先,应用本发明可以知道不同工艺参数对长寿命疲劳可靠性的影响权重,为优化工艺提供了方向;其次,应用本专利还能了解不同工艺参数对齿面粗糙度和残余应力等表面质量参数的影响权重,为提高弧齿锥齿表面完整性有很大帮助;再次,目前尚未有较成熟的由加工工艺参数评价齿轮长寿命可靠性的评价方法,本发明可操作性强、精确度高,只需经过简单的试验和相应的分析计算,即可得到弧齿锥齿轮在真实工况下的长寿命疲劳可靠性灵敏度,并对磨削和热处理工艺参数进行评价。
附图说明
图1为本发明的方法流程图;
具体实施方式
下面结合附图以及具体实施例进一步说明本发明。
本发明一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其方法流程如下:
步骤(1)、研究磨削工艺参数对齿面粗糙度和残余应力的灵敏度的影响;
步骤(2)、研究热处理工艺参数对齿面硬度和渗碳深度的灵敏度的影响;
步骤(3)、研究齿面粗糙度、残余应力、齿面硬度和渗碳深度对弧齿锥齿疲劳长寿命可靠性的灵敏度;
步骤(4)、计算得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度,通过灵敏度来评价磨削和热处理。
所述步骤(1)中齿轮磨削时,齿轮材料、磨削液等因素是一定的,磨削深度ap、磨削速度Vs和横向进给速度Vw是变化的,选择这3个因素作为正交试验因素,根据磨削工艺经验推荐值,每因素分别取若干个不同水平值,建立正交试验表。开展磨削试验,记录每组磨削试验后齿轮的表面粗糙度Ra、残余应力S11。
所述步骤(1)中用有限差分法对表面粗糙度和残余应力进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δx,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式。分别计算出磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差的灵敏度,
Figure GDA0002199192450000051
式中:βx表示y变量对x变量的灵敏度,x变量表示磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差,y变量表示表面粗糙度Ra、残余应力S11;x0为起始的x变量的值,y0为x0对应的值,x1为微小变化之后的x变量的值,y1为x1对应的值。
进一步的,所述步骤(2)中齿轮热处理时,齿轮材料、模数等因素一定,渗碳时间和渗碳温度是变化的,选择这2个因素作为正交试验因素,根据热处理工艺经验推荐值,没因素分别取若干个不同的水平值,建立正交实验表。开展热处理试验,记录每组热处理试验后的齿面硬度和渗碳深度。
所述步骤(2)中用有限差分法对齿面硬度和渗碳深度进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δa,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式。分别计算出渗碳时间和渗碳温度两个参数的均值以及方差的灵敏度,
Figure GDA0002199192450000061
式中:βa表示z变量对a变量的灵敏度,a变量表示渗碳时间和渗碳温度两个参数的均值以及方差,z变量表示齿面硬度和渗碳深度;a0为起始的a变量的值,z0为a0对应的值,a1为微小变化之后的a变量的值,z1为a1对应的值。
所述步骤(3)中航空弧齿锥齿轮齿面接触疲劳强度公式为:
σ′Hlim=σHlimZNZLZVZRZWZx
式中:σ`Hlim表示齿面接触疲劳强度,σHlimZN表示接触应力,ZL表示润滑剂系数,Zv表示速度系数,ZR表示粗糙度系数,ZW表示工作硬化系数,ZX表示尺寸系数。
齿根弯曲疲劳强度公式为:
σ′Flim=σFlimYNYSTYσYRYx
式中:σ`Flim表示齿根弯曲疲劳强度,σFlimYN表示弯曲应力,YST表示应力修正系数,Yσ表示相对齿根圆角敏感系数,YR表示相对齿根表面状况系数,Yx表示弯曲强度尺寸系数。
齿面接触应力公式为:
Figure GDA0002199192450000062
式中:σH表示齿面接触应力,ZM-B表示中点区域系数,ZH表示节点啮合区域系数,ZE表示弹性影响系数,ZLS表示载荷分担系数,Zβ表示螺旋角系数,ZK表示锥齿轮系数,Fmt表示切向力,KA表示使用系数,KV表示动载系数,K表示齿向载荷分布系数,K表示端面载荷分配系数,dv1表示小齿轮的分度圆直径,lbm表示接触线长度,uv表示齿数比。
齿根弯曲应力公式为:
Figure GDA0002199192450000071
式中:σF表示弯曲应力,Fmt表示切向力,KA表示使用系数,KV表示动载系数,K表示齿向载荷分布系数,K表示端面载荷分配系数,YFa表示齿形系数,YSa表示应力修正系数,Yε表示重合度系数,YK表示锥齿轮系数,YLS表示载荷分担系数,b表示工作齿宽,mmn表示小齿轮的法向模数。
表面粗糙度、齿面硬度、精度等轮齿表面质量参数直接影响上述公式中的粗糙度系数ZR、工作硬化系数ZW、齿根表面状况系数YR以及可靠性分析中某些参数的方差分布的选取,所以会对齿轮疲劳强度和可靠性造成影响。
所述步骤(3)中研究齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数对齿轮疲劳可靠性的敏感性影响采用有限差分法对弧齿锥齿轮进行灵敏度分析。其基本做法是使设计变量有一个微小的摄动,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,
Figure GDA0002199192450000072
式中:βy表示R变量对y变量的灵敏度,y变量表示表面粗糙度Ra、残余应力S11两个表面参数的均值以及方差,R变量表示可靠性;y0为起始的y变量的值,R0为y0对应的值,y1为微小变化之后的y变量的值,R1为y1对应的值。
Figure GDA0002199192450000073
式中:βz表示R变量对z变量的灵敏度,z变量表示齿面硬度和渗碳深度两个表面参数的均值以及方差,R变量表示可靠性;z0为起始的z变量的值,R0为z0对应的值,z1为微小变化之后的z变量的值,R1为z1对应的值。
所述步骤(4)中根据求导规则,将齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数当做中间变量消掉,最终得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度:
Figure GDA0002199192450000074
Figure GDA0002199192450000081
式中:βx(R)表示R变量对x变量的灵敏度,x变量表示磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差,y变量表示表面粗糙度Ra、残余应力S11;x0为起始的x变量的值,y0为x0对应的值,x1为微小变化之后的x变量的值,y1为x1对应的值;
式中:βa(R)表示R变量对a变量的灵敏度,a变量表示渗碳时间和渗碳温度两个参数的均值以及方差,z变量表示齿面硬度和渗碳深度;a0为起始的a变量的值,z0为a0对应的值,a1为微小变化之后的a变量的值,z1为a1对应的值。
最终通过灵敏度来评价磨削和热处理。
具体的,本发明的流程图如图1所示。下面以如下所示的某弧齿锥齿轮传动系统为例,具体说明本发明方法,但本发明的保护范围不限于下述实例:
小轮齿数26,大论齿数31,大端模数8.654mm,齿宽57mm,中点螺旋角35°,压力角20°,传递功率1051KW,输入转速2200r/min,使用寿命500h,精度等级6级,小轮表面粗糙度0.8,大轮表面粗糙度0.8,40℃时的润滑油黏度177.2mm2/s,齿轮材料属性为渗碳淬火的渗碳钢,材料硬度HRC59,拉伸极限1180MPa,接触疲劳极限1500MPa,弯曲疲劳极限480Mpa,材料密度7.88E-6kg/mm3,弹性模量2.07E+5MPa,泊松比0.3。
(1)研究磨削工艺参数对齿面粗糙度和残余应力的灵敏度的影响;
通过对上述齿轮磨削加工过程一组工艺参数(磨削速度、磨削深度和横向进给速度)与齿轮粗糙度结果如表1所示:
表1磨削工艺参数与齿轮表面粗糙度数据
Figure GDA0002199192450000082
Figure GDA0002199192450000091
应用有限差分法得到磨削工艺参数对表面粗糙度和残余应力的灵敏度的结果如表2、表3所示:
表2磨削工艺参数与齿轮表面粗糙度灵敏度分析结果
Figure GDA0002199192450000092
从表中分析可以知道磨削深度对表面粗糙度的影响最明显,在加工过程中控制表面粗糙度应该着重的去控制磨削深度。
表3磨削工艺参数与齿轮残余灵敏度分析结果
Figure GDA0002199192450000093
从表中分析可以知道磨削深度对残余应力的影响最明显,在加工过程中控制残余应力应该着重的去控制磨削深度。
(2)研究热处理工艺参数对齿面硬度和渗碳深度的灵敏度的影响;
在强渗期碳势0.95~1%c、扩散期碳势0.9~0.95%c下,搜集型号齿轮的部分渗碳工艺(渗碳温度和渗碳时间)与渗碳深度数据如表4所示:
表4 0.95~1%c强渗碳势下渗碳温度、时间与渗层深度关系数据
Figure GDA0002199192450000094
Figure GDA0002199192450000101
应用有限差分法得到磨削工艺参数对表面粗糙度和残余应力的灵敏度的结果如表5所示:
表5 0.95~1%c强渗碳势下渗碳温度、时间与渗层深度的灵敏度分析结果
Figure GDA0002199192450000102
从表中分析可以知道渗碳时间对渗碳深度的影响最明显,在加工过程中控制渗碳深度应该着重的去控制渗碳时间。
(3)研究齿面粗糙度、残余应力、齿面硬度和渗碳深度对弧齿锥齿疲劳长寿命可靠性的灵敏度;
以研究齿面粗糙度和残余应力对弧齿锥齿疲劳长寿命可靠性的灵敏度为例,应用有限差分法得到齿面粗糙度和残余应力对弧齿锥齿疲劳长寿命可靠性的灵敏度如表6所示:
表6粗糙度和残余应力对疲劳可靠性敏感性分析结果
Figure GDA0002199192450000103
(4)计算得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度,通过灵敏度来评价磨削和热处理。
以研究磨削工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度为例,将齿面粗糙度、残余应力这些表面参数当做中间变量消掉,最终得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度:
表7磨削工艺参数对疲劳可靠性敏感性分析结果
Figure GDA0002199192450000111
最终得到的磨削工艺参数对长寿命疲劳可靠性的评价结果。从结果中可以知道磨削深度对长寿命疲劳可靠性的影响最大,磨削速度次之,横向进给速度最小,为保证可靠性,首先得控制磨削深度。
总之,本发明针对弧齿锥齿轮长寿命传动可靠性的评价的问题,以加工工艺参数为正交试验因素,通过正交试验并对试验数据进行有限差分法处理,得到了工艺参数对表面完整性的灵敏度,和表面完整性对长寿命疲劳可靠性的灵敏度,最终计算出加工工艺参数对长寿命疲劳可靠性的灵敏度,通过灵敏度对加工工艺参数进行评价。从而为弧齿锥齿轮的寿命预测和工艺评价工作提供重要的依据。

Claims (5)

1.一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于,该方法实现步骤如下:
步骤(1)、研究磨削工艺参数对齿面粗糙度和残余应力的灵敏度的影响;
步骤(2)、研究热处理工艺参数对齿面硬度和渗碳深度的灵敏度的影响;
步骤(3)、研究齿面粗糙度、残余应力、齿面硬度和渗碳深度对弧齿锥齿疲劳长寿命可靠性的灵敏度;
所述步骤(3)中研究齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数对齿轮疲劳可靠性的敏感性影响采用有限差分法对弧齿锥齿轮进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,
Figure FDA0002361317650000011
式中:βy表示R变量对y变量的灵敏度,y变量表示齿面粗糙度Ra、残余应力S11两个表面参数的均值以及方差,R变量表示可靠性;y0为起始的y变量的值,R0为y0对应的值,y1为微小变化之后的y变量的值,R1为y1对应的值,
Figure FDA0002361317650000012
式中:βz表示R变量对z变量的灵敏度,z变量表示齿面硬度和渗碳深度两个表面参数的均值以及方差,R变量表示可靠性;z0为起始的z变量的值,R0为z0对应的值,z1为微小变化之后的z变量的值,R1为z1对应的值;
步骤(4)、计算得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度,通过灵敏度来评价磨削和热处理;
所述步骤(4)中根据求导规则,将齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数当做中间变量消掉,最终得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度:
Figure FDA0002361317650000013
Figure FDA0002361317650000014
式中:βx(R)表示R变量对x变量的灵敏度,βx表示y变量对x变量的灵敏度,x变量表示磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差,y变量表示齿面粗糙度Ra、残余应力S11;x0为起始的x变量的值,y0为x0对应的值,x1为微小变化之后的x变量的值,y1为x1对应的值;
式中:βa(R)表示R变量对a变量的灵敏度,βa表示z变量对a变量的灵敏度,a变量表示渗碳时间和渗碳温度两个参数的均值以及方差,z变量表示齿面硬度和渗碳深度;a0为起始的a变量的值,z0为a0对应的值,a1为微小变化之后的a变量的值,z1为a1对应的值;
最终通过灵敏度来评价磨削和热处理。
2.根据权利要求1所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(1)中齿轮磨削时,齿轮材料、磨削液因素是一定的,磨削深度ap、磨削速度Vs和横向进给速度Vw是变化的,选择这3个因素作为正交试验因素,根据磨削工艺经验推荐值,每因素分别取若干个不同水平值,建立正交试验表,开展磨削试验,记录每组磨削试验后齿轮的齿面粗糙度Ra、残余应力S11。
3.根据权利要求2所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(1)中用有限差分法对齿面粗糙度和残余应力进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δx,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,分别计算出磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差的灵敏度,
Figure FDA0002361317650000021
式中:βx表示y变量对x变量的灵敏度,x变量表示磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差,y变量表示齿面粗糙度Ra、残余应力S11;x0为起始的x变量的值,y0为x0对应的值,x1为微小变化之后的x变量的值,y1为x1对应的值。
4.根据权利要求1所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(2)中齿轮热处理时,齿轮材料、模数因素一定,渗碳时间和渗碳温度是变化的,选择这2个因素作为正交试验因素,根据热处理工艺经验推荐值,每因素分别取若干个不同的水平值,建立正交实验表,开展热处理试验,记录每组热处理试验后的齿面硬度和渗碳深度。
5.根据权利要求4所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(2)中用有限差分法对齿面硬度和渗碳深度进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δa,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,分别计算出渗碳时间和渗碳温度两个参数的均值以及方差的灵敏度,
Figure FDA0002361317650000031
式中:βa表示z变量对a变量的灵敏度,a变量表示渗碳时间和渗碳温度两个参数的均值以及方差,z变量表示齿面硬度和渗碳深度;a0为起始的a变量的值,z0为a0对应的值,a1为微小变化之后的a变量的值,z1为a1对应的值。
CN201610937468.2A 2016-10-25 2016-10-25 一种弧齿锥齿轮长寿命传动疲劳可靠性的评价方法 Active CN106649971B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610937468.2A CN106649971B (zh) 2016-10-25 2016-10-25 一种弧齿锥齿轮长寿命传动疲劳可靠性的评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610937468.2A CN106649971B (zh) 2016-10-25 2016-10-25 一种弧齿锥齿轮长寿命传动疲劳可靠性的评价方法

Publications (2)

Publication Number Publication Date
CN106649971A CN106649971A (zh) 2017-05-10
CN106649971B true CN106649971B (zh) 2020-05-26

Family

ID=58821112

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610937468.2A Active CN106649971B (zh) 2016-10-25 2016-10-25 一种弧齿锥齿轮长寿命传动疲劳可靠性的评价方法

Country Status (1)

Country Link
CN (1) CN106649971B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109271711B (zh) * 2018-09-25 2023-03-28 重庆大学 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法
CN109472068B (zh) * 2018-10-26 2022-11-18 北京航空航天大学 一种基于多失效模式的重载弧齿锥齿轮定长寿命传动疲劳工艺可靠性分析与评价方法
CN110851922B (zh) * 2019-09-04 2022-07-19 湖南大学 一种基于齿面闪温最低的斜齿轮最坏螺旋角的确定方法
CN112100938B (zh) * 2020-09-11 2021-05-14 南京航空航天大学 一种飞溅润滑锥齿轮温度的数值计算方法及系统
CN112948989A (zh) * 2021-01-26 2021-06-11 浙江大学 一种数据驱动的齿面强化方法
CN113297703B (zh) * 2021-06-15 2022-08-19 中国航发沈阳发动机研究所 一种面齿轮齿面接触疲劳强度评估方法
CN113742859B (zh) * 2021-08-13 2024-03-08 湖南磐钴传动科技有限公司 一种基于齿面残余应力与接触应力的复合应力场构建方法
CN116818290A (zh) * 2023-05-15 2023-09-29 江苏科技大学 一种综合考虑硬度、残余应力和粗糙度的磨削加工试件疲劳强度预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11230312A (ja) * 1998-02-13 1999-08-27 Nippon Steel Corp 高接触疲労寿命歯車及びその製造方法
CN105710454A (zh) * 2016-03-31 2016-06-29 大连交通大学 小型弧齿锥齿轮数控蜗盘砂轮研齿机及其加工步骤
CN105718606A (zh) * 2014-08-18 2016-06-29 鲍珂 一种考虑失效模式相关性的车辆重载齿轮可靠性预计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11230312A (ja) * 1998-02-13 1999-08-27 Nippon Steel Corp 高接触疲労寿命歯車及びその製造方法
CN105718606A (zh) * 2014-08-18 2016-06-29 鲍珂 一种考虑失效模式相关性的车辆重载齿轮可靠性预计方法
CN105710454A (zh) * 2016-03-31 2016-06-29 大连交通大学 小型弧齿锥齿轮数控蜗盘砂轮研齿机及其加工步骤

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Roughness model for tooth surfaces of spiral bevel gears under;YanZhong Wang 等;《Mechanism and Machine Theory》;20160604;第104卷;第17-30页 *
基于DEFORM有限元仿真的弧齿锥齿轮热处理过程残余应力与变形分析;王延忠 等;《机械传动》;20160115;第40卷(第1期);第148-152页 *
航空弧齿锥齿轮轮齿应力谱精确获取及分布估计;王延忠 等;《新技术新工艺》;20160925(第9期);第16-21页 *

Also Published As

Publication number Publication date
CN106649971A (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
CN106649971B (zh) 一种弧齿锥齿轮长寿命传动疲劳可靠性的评价方法
CN106503350B (zh) 一种弧齿锥齿轮长寿命传动疲劳可靠性的设计方法
GB2564919B (en) Gearsets
Wang et al. Evaluation of rolling contact fatigue of a carburized wind turbine gear considering the residual stress and hardness gradient
CN109766670B (zh) 基于Copula函数的工业机器人用谐波减速器可靠性分析方法
Molaie et al. Spiral bevel gears: nonlinear dynamic model based on accurate static stiffness evaluation
Zhou et al. Investigation on stress microcycles and mild wear mechanism in gear contact fatigue
Errichello et al. Wind turbine tribology seminar-a recap
Mao et al. Numerical calculation method of meshing stiffness for the beveloid gear considering the effect of surface topography
Honkalas et al. A review on design and efficiency improvement of worm and worm wheel of a gear motor
Mo et al. Elastohydrodynamic lubrication characteristics of spiral bevel gear subjected to shot peening treatment
CN102029442B (zh) 一种提高螺旋锥齿轮润滑性能的方法
Krantz et al. Pitting and bending fatigue evaluations of a new case-carburized gear steel
Glodež et al. Determination of service life of sintered Powder metallurgy gears in regard to tooth bending fatigue
CN108274081B (zh) 一种螺旋锥齿轮的变形补偿加工方法
Patel et al. Topology Optimization of Automotive Gear using Fea
CN112287554B (zh) 一种预测齿轮表面接触疲劳损伤的方法
Cui et al. Fatigue life analysis of spur gears with precise tooth profile surfaces
de Sousa Neves et al. Analyzing micro-geometric deviations of a helical gear profile as a function of material and machining parameters in the hobbing, shaving and carbonitriding processes
Franco et al. Experimental uncertainty analysis of gear fatigue life theoretical prediction
CN112084600A (zh) 一种重载复杂工况螺旋锥齿轮的喷丸工艺参数优选方法
Huang et al. Study on Chaotic Characteristics of the Friction Process between High Hardness Alloy Steel and Cemented Carbide under C60 Nanoparticle Fluid Lubrication.
Velex International Gear Conference 2014: 26th-28th August 2014, Lyon
Kopiláková et al. Comparison HRC and CC gearing for damage to pitting
Al Tubi Effects of Variable Load and Rotational Speed Conditions on Gear Micropitting in Wind Turbine Gearboxes

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210817

Address after: Room 917, No. 58-1, Beizhan Second Road, Shenhe District, Shenyang, Liaoning 110013

Patentee after: Shenyang Dechuang metal materials Co.,Ltd.

Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road

Patentee before: BEIHANG University