CN106649971A - 一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法 - Google Patents
一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法 Download PDFInfo
- Publication number
- CN106649971A CN106649971A CN201610937468.2A CN201610937468A CN106649971A CN 106649971 A CN106649971 A CN 106649971A CN 201610937468 A CN201610937468 A CN 201610937468A CN 106649971 A CN106649971 A CN 106649971A
- Authority
- CN
- China
- Prior art keywords
- represent
- variables
- grinding
- sensitivity
- coefficient
- 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
Classifications
-
- 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
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三个参数的均值以及方差的灵敏度,
式中:βx表示y变量对x变量的灵敏度,x变量表示磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差,y变量表示表面粗糙度Ra、残余应力S11;x0为起始的x变量的值,y0为x0对应的值,x1为微小变化之后的x变量的值,y1为x1对应的值。
进一步的,所述步骤(2)中齿轮热处理时,齿轮材料、模数因素一定,渗碳时间和渗碳温度是变化的,选择这2个因素作为正交试验因素,根据热处理工艺经验推荐值,没因素分别取若干个不同的水平值,建立正交实验表。开展热处理试验,记录每组热处理试验后的齿面硬度和渗碳深度。
进一步的,所述步骤(2)中用有限差分法对齿面硬度和渗碳深度进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δa,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式。分别计算出渗碳时间和渗碳温度两个参数的均值以及方差的灵敏度,
式中:β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表示弯曲强度尺寸系数;
齿面接触应力公式为:
式中:σH表示齿面接触应力,ZM-B表示中点区域系数,ZH表示节点啮合区域系数,ZE表示弹性影响系数,ZLS表示载荷分担系数,Zβ表示螺旋角系数,ZK表示锥齿轮系数,Fmt表示切向力,KA表示使用系数,KV表示动载系数,KHβ表示齿向载荷分布系数,KHα表示端面载荷分配系数,dv1表示小齿轮的分度圆直径,lbm表示接触线长度,uv表示齿数比。
齿根弯曲应力公式为:
式中:σF表示弯曲应力,Fmt表示切向力,KA表示使用系数,KV表示动载系数,KFβ表示齿向载荷分布系数,KFα表示端面载荷分配系数,YFa表示齿形系数,YSa表示应力修正系数,Yε表示重合度系数,YK表示锥齿轮系数,YLS表示载荷分担系数,b表示工作齿宽,mmn表示小齿轮的法向模数。
表面粗糙度、齿面硬度、精度等轮齿表面质量参数直接影响上述公式中的粗糙度系数ZR、工作硬化系数ZW、齿根表面状况系数YR以及可靠性分析中某些参数的方差分布的选取,所以会对齿轮疲劳强度和可靠性造成影响。
进一步的,所述步骤(3)中研究齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数对齿轮疲劳可靠性的敏感性影响采用有限差分法对弧齿锥齿轮进行灵敏度分析。其基本做法是使设计变量有一个微小的摄动,用差分格式计算输出即齿轮可靠性对设计变量的近 似导数,采用向前差分格式,
式中βy表示R变量对y变量的灵敏度,y变量表示表面粗糙度Ra、残余应力S11两个表面参数的均值以及方差,R变量表示可靠性;y0为起始的y变量的值,R0为y0对应的值,y1为微小变化之后的y变量的值,R1为y1对应的值,
式中:βz表示R变量对z变量的灵敏度,z变量表示齿面硬度和渗碳深度两个表面参数的均值以及方差,R变量表示可靠性;z0为起始的z变量的值,R0为z0对应的值,z1为微小变化之后的z变量的值,R1为z1对应的值。
进一步的,所述步骤(4)中根据求导规则,将齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数当做中间变量消掉,最终得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度:
式中:β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三个参数的均值以及方差的灵敏度,
式中:βx表示y变量对x变量的灵敏度,x变量表示磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差,y变量表示表面粗糙度Ra、残余应力S11;x0为起始的x变量的值,y0为x0对应的值,x1为微小变化之后的x变量的值,y1为x1对应的值。
进一步的,所述步骤(2)中齿轮热处理时,齿轮材料、模数等因素一定,渗碳时间和渗碳温度是变化的,选择这2个因素作为正交试验因素,根据热处理工艺经验推荐值,没因素分别取若干个不同的水平值,建立正交实验表。开展热处理试验,记录每组热处理试验后的齿面硬度和渗碳深度。
所述步骤(2)中用有限差分法对齿面硬度和渗碳深度进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δa,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式。分别计算出渗碳时间和渗碳温度两个参数的均值以及方差的灵敏度,
式中:β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表示弯曲强度尺寸系数。
齿面接触应力公式为:
式中:σH表示齿面接触应力,ZM-B表示中点区域系数,ZH表示节点啮合区域系数,ZE表示弹性影响系数,ZLS表示载荷分担系数,Zβ表示螺旋角系数,ZK表示锥齿轮系数,Fmt表示切向力,KA表示使用系数,KV表示动载系数,KHβ表示齿向载荷分布系数,KHα表示端面载荷分配系数,dv1表示小齿轮的分度圆直径,lbm表示接触线长度,uv表示齿数比。
齿根弯曲应力公式为:
式中:σF表示弯曲应力,Fmt表示切向力,KA表示使用系数,KV表示动载系数,KFβ表示齿向载荷分布系数,KFα表示端面载荷分配系数,YFa表示齿形系数,YSa表示应力修正系数,Yε表示重合度系数,YK表示锥齿轮系数,YLS表示载荷分担系数,b表示工作齿宽,mmn表示小齿轮的法向模数。
表面粗糙度、齿面硬度、精度等轮齿表面质量参数直接影响上述公式中的粗糙度系数ZR、工作硬化系数ZW、齿根表面状况系数YR以及可靠性分析中某些参数的方差分布的选取,所以会对齿轮疲劳强度和可靠性造成影响。
所述步骤(3)中研究齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数对齿轮疲劳可靠性的敏感性影响采用有限差分法对弧齿锥齿轮进行灵敏度分析。其基本做法是使设计变量有一个微小的摄动,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,
式中:βy表示R变量对y变量的灵敏度,y变量表示表面粗糙度Ra、残余应力S11两个表面参数的均值以及方差,R变量表示可靠性;y0为起始的y变量的值,R0为y0对应的值,y1为微小变化之后的y变量的值,R1为y1对应的值。
式中:βz表示R变量对z变量的灵敏度,z变量表示齿面硬度和渗碳深度两个表面参数的均值以及方差,R变量表示可靠性;z0为起始的z变量的值,R0为z0对应的值,z1为微小变化之后的z变量的值,R1为z1对应的值。
所述步骤(4)中根据求导规则,将齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数当做中间变量消掉,最终得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度:
式中:β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磨削工艺参数与齿轮表面粗糙度数据
应用有限差分法得到磨削工艺参数对表面粗糙度和残余应力的灵敏度的结果如表2、表3所示:
表2磨削工艺参数与齿轮表面粗糙度灵敏度分析结果
从表中分析可以知道磨削深度对表面粗糙度的影响最明显,在加工过程中控制表面粗糙度应该着重的去控制磨削深度。
表3磨削工艺参数与齿轮残余灵敏度分析结果
从表中分析可以知道磨削深度对残余应力的影响最明显,在加工过程中控制残余应力应该着重的去控制磨削深度。
(2)研究热处理工艺参数对齿面硬度和渗碳深度的灵敏度的影响;
在强渗期碳势0.95~1%c、扩散期碳势0.9~0.95%c下,搜集型号齿轮的部分渗碳工艺(渗碳温度和渗碳时间)与渗碳深度数据如表4所示:
表4 0.95~1%c强渗碳势下渗碳温度、时间与渗层深度关系数据
应用有限差分法得到磨削工艺参数对表面粗糙度和残余应力的灵敏度的结果如表5所示:
表5 0.95~1%c强渗碳势下渗碳温度、时间与渗层深度的灵敏度分析结果
从表中分析可以知道渗碳时间对渗碳深度的影响最明显,在加工过程中控制渗碳深度应该着重的去控制渗碳时间。
(3)研究齿面粗糙度、残余应力、齿面硬度和渗碳深度对弧齿锥齿疲劳长寿命可靠性的灵敏度;
以研究齿面粗糙度和残余应力对弧齿锥齿疲劳长寿命可靠性的灵敏度为例,应用有限差分法得到齿面粗糙度和残余应力对弧齿锥齿疲劳长寿命可靠性的灵敏度如表6所示:
表6粗糙度和残余应力对疲劳可靠性敏感性分析结果
(4)计算得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度,通过灵敏度来评价磨削和热处理。
以研究磨削工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度为例,将齿面粗糙度、残余应力这些表面参数当做中间变量消掉,最终得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度:
表7 磨削工艺参数对疲劳可靠性敏感性分析结果
最终得到的磨削工艺参数对长寿命疲劳可靠性的评价结果。从结果中可以知道磨削深度对长寿命疲劳可靠性的影响最大,磨削速度次之,横向进给速度最小,为保证可靠性,首先得控制磨削深度。
总之,本发明针对弧齿锥齿轮长寿命传动可靠性的评价的问题,以加工工艺参数为正交试验因素,通过正交试验并对试验数据进行有限差分法处理,得到了工艺参数对表面完整性的灵敏度,和表面完整性对长寿命疲劳可靠性的灵敏度,最终计算出加工工艺参数对长寿命疲劳可靠性的灵敏度,通过灵敏度对加工工艺参数进行评价。从而为弧齿锥齿轮的寿命预测和工艺评价工作提供重要的依据。
Claims (8)
1.一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于,该方法实现步骤如下:
步骤(1)、研究磨削工艺参数对齿面粗糙度和残余应力的灵敏度的影响;
步骤(2)、研究热处理工艺参数对齿面硬度和渗碳深度的灵敏度的影响;
步骤(3)、研究齿面粗糙度、残余应力、齿面硬度和渗碳深度对弧齿锥齿疲劳长寿命可靠性的灵敏度;
步骤(4)、计算得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度,通过灵敏度来评价磨削和热处理。
2.根据权利要求1所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(1)中齿轮磨削时,齿轮材料、磨削液因素是一定的,磨削深度ap、磨削速度Vs和横向进给速度Vw是变化的,选择这3个因素作为正交试验因素,根据磨削工艺经验推荐值,每因素分别取若干个不同水平值,建立正交试验表,开展磨削试验,记录每组磨削试验后齿轮的表面粗糙度Ra、残余应力S11。
3.根据权利要求2所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(1)中用有限差分法对表面粗糙度和残余应力进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动Δx,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,分别计算出磨削深度ap、磨削速度Vs和横向进给速度Vw三个参数的均值以及方差的灵敏度,
式中:β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,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,分别计算出渗碳时间和渗碳温度两个参数的均值以及方差的灵敏度,
式中:βa表示z变量对a变量的灵敏度,a变量表示渗碳时间和渗碳温度两个参数的均值以及方差,z变量表示齿面硬度和渗碳深度;a0为起始的a变量的值,z0为a0对应的值,a1为微小变化之后的a变量的值,z1为a1对应的值。
6.根据权利要求1所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(3)中航空弧齿锥齿轮齿面接触疲劳强度公式为:
σ'Hlim=σHlimZNZLZVZRZWZx
式中:表示齿面接触疲劳强度,σHlimZN表示接触应力,ZL表示润滑剂系数,Zv表示速度系数,ZR表示粗糙度系数,ZW表示工作硬化系数,ZX表示尺寸系数;
齿根弯曲疲劳强度公式为:
σ′Flim=σFlimYNYSTYσYRYx
式中:表示齿根弯曲疲劳强度,σFlimYN表示弯曲应力,YST表示应力修正系数,Yσ表示相对齿根圆角敏感系数,YR表示相对齿根表面状况系数,Yx表示弯曲强度尺寸系数;
齿面接触应力公式为:
式中:σH表示齿面接触应力,ZM-B表示中点区域系数,ZH表示节点啮合区域系数,ZE表示弹性影响系数,ZLS表示载荷分担系数,Zβ表示螺旋角系数,ZK表示锥齿轮系数,Fmt表示切向力,KA表示使用系数,KV表示动载系数,KHβ表示齿向载荷分布系数,KHα表示端面载荷分配系数,dv1表示小齿轮的分度圆直径,lbm表示接触线长度,uv表示齿数比;
齿根弯曲应力公式为:
式中:σF表示弯曲应力,Fmt表示切向力,KA表示使用系数,KV表示动载系数,KFβ表示齿向载荷分布系数,KFα表示端面载荷分配系数,YFa表示齿形系数,YSa表示应力修正系数,Yε表示重合度系数,YK表示锥齿轮系数,YLS表示载荷分担系数,b表示工作齿宽,mmn表示小齿轮的法向模数;
表面粗糙度、齿面硬度、精度等轮齿表面质量参数直接影响上述公式中的粗糙度系数ZR、工作硬化系数ZW、齿根表面状况系数YR以及可靠性分析中某些参数的方差分布的选取,所以会对齿轮疲劳强度和可靠性造成影响。
7.根据权利要求6所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(3)中研究齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数对齿轮疲劳可靠性的敏感性影响采用有限差分法对弧齿锥齿轮进行灵敏度分析,其基本做法是使设计变量有一个微小的摄动,用差分格式计算输出即齿轮可靠性对设计变量的近似导数,采用向前差分格式,
式中:βy表示R变量对y变量的灵敏度,y变量表示表面粗糙度Ra、残余应力S11两个表面参数的均值以及方差,R变量表示可靠性;y0为起始的y变量的值,R0为y0对应的值,y1为微小变化之后的y变量的值,R1为y1对应的值,
式中:βz表示R变量对z变量的灵敏度,z变量表示齿面硬度和渗碳深度两个表面参数的均值以及方差,R变量表示可靠性;z0为起始的z变量的值,R0为z0对应的值,z1为微小变化之后的z变量的值,R1为z1对应的值。
8.根据权利要求1所述的基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法,其特征在于:所述步骤(4)中根据求导规则,将齿面粗糙度、残余应力、齿面硬度和渗碳深度这些表面参数当做中间变量消掉,最终得到磨削和热处理工艺参数对弧齿锥齿疲劳长寿命可靠性的灵敏度:
式中:β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对应的值;
最终通过灵敏度来评价磨削和热处理。
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 true CN106649971A (zh) | 2017-05-10 |
CN106649971B 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) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109271711A (zh) * | 2018-09-25 | 2019-01-25 | 重庆大学 | 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法 |
CN109472068A (zh) * | 2018-10-26 | 2019-03-15 | 北京航空航天大学 | 一种基于多失效模式的重载弧齿锥齿轮定长寿命传动疲劳工艺可靠性分析与评价方法 |
CN110851922A (zh) * | 2019-09-04 | 2020-02-28 | 湖南大学 | 一种基于齿面闪温最低的斜齿轮最坏螺旋角的确定方法 |
CN112100938A (zh) * | 2020-09-11 | 2020-12-18 | 南京航空航天大学 | 一种飞溅润滑锥齿轮温度的数值计算方法及系统 |
CN112948989A (zh) * | 2021-01-26 | 2021-06-11 | 浙江大学 | 一种数据驱动的齿面强化方法 |
CN113297703A (zh) * | 2021-06-15 | 2021-08-24 | 中国航发沈阳发动机研究所 | 一种面齿轮齿面接触疲劳强度评估方法 |
CN113742859A (zh) * | 2021-08-13 | 2021-12-03 | 湖南磐钴传动科技有限公司 | 一种基于齿面残余应力与接触应力的复合应力场构建方法 |
CN116818290A (zh) * | 2023-05-15 | 2023-09-29 | 江苏科技大学 | 一种综合考虑硬度、残余应力和粗糙度的磨削加工试件疲劳强度预测方法 |
Citations (3)
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 | 鲍珂 | 一种考虑失效模式相关性的车辆重载齿轮可靠性预计方法 |
-
2016
- 2016-10-25 CN CN201610937468.2A patent/CN106649971B/zh active Active
Patent Citations (3)
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)
Title |
---|
YANZHONG WANG 等: "Roughness model for tooth surfaces of spiral bevel gears under", 《MECHANISM AND MACHINE THEORY》 * |
王延忠 等: "基于DEFORM有限元仿真的弧齿锥齿轮热处理过程残余应力与变形分析", 《机械传动》 * |
王延忠 等: "航空弧齿锥齿轮轮齿应力谱精确获取及分布估计", 《新技术新工艺》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109271711A (zh) * | 2018-09-25 | 2019-01-25 | 重庆大学 | 一种考虑不均匀特性的渗碳硬化齿轮有限元建模方法 |
CN109472068A (zh) * | 2018-10-26 | 2019-03-15 | 北京航空航天大学 | 一种基于多失效模式的重载弧齿锥齿轮定长寿命传动疲劳工艺可靠性分析与评价方法 |
CN109472068B (zh) * | 2018-10-26 | 2022-11-18 | 北京航空航天大学 | 一种基于多失效模式的重载弧齿锥齿轮定长寿命传动疲劳工艺可靠性分析与评价方法 |
CN110851922B (zh) * | 2019-09-04 | 2022-07-19 | 湖南大学 | 一种基于齿面闪温最低的斜齿轮最坏螺旋角的确定方法 |
CN110851922A (zh) * | 2019-09-04 | 2020-02-28 | 湖南大学 | 一种基于齿面闪温最低的斜齿轮最坏螺旋角的确定方法 |
CN112100938A (zh) * | 2020-09-11 | 2020-12-18 | 南京航空航天大学 | 一种飞溅润滑锥齿轮温度的数值计算方法及系统 |
CN112100938B (zh) * | 2020-09-11 | 2021-05-14 | 南京航空航天大学 | 一种飞溅润滑锥齿轮温度的数值计算方法及系统 |
CN112948989A (zh) * | 2021-01-26 | 2021-06-11 | 浙江大学 | 一种数据驱动的齿面强化方法 |
CN113297703A (zh) * | 2021-06-15 | 2021-08-24 | 中国航发沈阳发动机研究所 | 一种面齿轮齿面接触疲劳强度评估方法 |
CN113297703B (zh) * | 2021-06-15 | 2022-08-19 | 中国航发沈阳发动机研究所 | 一种面齿轮齿面接触疲劳强度评估方法 |
CN113742859A (zh) * | 2021-08-13 | 2021-12-03 | 湖南磐钴传动科技有限公司 | 一种基于齿面残余应力与接触应力的复合应力场构建方法 |
CN113742859B (zh) * | 2021-08-13 | 2024-03-08 | 湖南磐钴传动科技有限公司 | 一种基于齿面残余应力与接触应力的复合应力场构建方法 |
CN116818290A (zh) * | 2023-05-15 | 2023-09-29 | 江苏科技大学 | 一种综合考虑硬度、残余应力和粗糙度的磨削加工试件疲劳强度预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106649971B (zh) | 2020-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106649971A (zh) | 一种基于磨削和热处理的弧齿锥齿轮长寿命传动疲劳可靠性的评价方法 | |
CN106503350B (zh) | 一种弧齿锥齿轮长寿命传动疲劳可靠性的设计方法 | |
Kishawy et al. | An energy based analysis of broaching operation: Cutting forces and resultant surface integrity | |
Li | Effects of machining errors, assembly errors and tooth modifications on loading capacity, load-sharing ratio and transmission error of a pair of spur gears | |
CN106372449A (zh) | 基于齿面接触动应力的弧齿锥齿轮多工况疲劳寿命快速预估方法 | |
Blagojevic et al. | Numerical and experimental analysis of the cycloid disc stress state | |
CN103206515A (zh) | 一种针对摆线齿锥齿轮误差加载齿面接触分析的方法 | |
Pinto et al. | Fretting fatigue under variable amplitude loading considering partial and gross slip regimes: Numerical analysis | |
Liu et al. | Influence of load spectrum on contact fatigue damage of a case carburized wind turbine gear | |
CN107451359A (zh) | 一种考虑基体裂纹影响的齿轮啮合特性有限元分析方法 | |
CN111291456A (zh) | 一种考虑强度退化和失效模式的齿轮可靠性分析方法 | |
Molaie et al. | Spiral bevel gears: nonlinear dynamic model based on accurate static stiffness evaluation | |
Chavadaki et al. | Finite element analysis of spur gear to find out the optimum root radius | |
WO2011101306A1 (en) | Method for testing the life-time of gear drives and the like | |
Wang et al. | Damage behavior due to rolling contact fatigue and bending fatigue of a gear using crystal plasticity modeling | |
CN103246759A (zh) | 一种滑动轴承的优化设计方法 | |
CN107609311A (zh) | 基于排屑力模型的深孔钻削深度优化方法 | |
Stupina et al. | Control and management by resource of rolling cutter bits in drilling rock massif | |
Cao et al. | Failure analysis of a pinion of the jacking system of a jack-up platform | |
He et al. | Shakedown analysis of a wind turbine gear considering strain-hardening and the initial residual stress | |
CN108274081A (zh) | 一种螺旋锥齿轮的变形补偿加工方法 | |
Ognjanović et al. | Assessment of Probability of Gear Tooth Side Wear of a Planetary Gearbox | |
Hanumanna et al. | Prediction of fatigue life of gear subjected to varying loads | |
CN109918618B (zh) | 齿轮可靠性设计的重要性抽样方法 | |
Lewicki | Effect of speed (centrifugal load) on gear crack propagation direction |
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 |