CN109030016A - 基于载荷谱的综合传动装置耐久性评估方法 - Google Patents

基于载荷谱的综合传动装置耐久性评估方法 Download PDF

Info

Publication number
CN109030016A
CN109030016A CN201810595308.3A CN201810595308A CN109030016A CN 109030016 A CN109030016 A CN 109030016A CN 201810595308 A CN201810595308 A CN 201810595308A CN 109030016 A CN109030016 A CN 109030016A
Authority
CN
China
Prior art keywords
load
formula
bearing
stress
under
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
Application number
CN201810595308.3A
Other languages
English (en)
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.)
China North Vehicle Research Institute
Original Assignee
China North Vehicle Research Institute
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 China North Vehicle Research Institute filed Critical China North Vehicle Research Institute
Priority to CN201810595308.3A priority Critical patent/CN109030016A/zh
Publication of CN109030016A publication Critical patent/CN109030016A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M17/00Testing of vehicles
    • G01M17/007Wheeled or endless-tracked vehicles

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明涉及一种基于载荷谱的综合传动装置耐久性评估方法,涉及履带车辆传动装置耐久性技术领域。本发明提供的基于载荷谱的综合传动装置耐久性评估计算方法,改变了以往综合传动整机耐久性计算没有统一载荷谱,各零部件耐久性计算无法关联匹配,整机耐久性评估没有方法参考的局面,从根本上解决了综合传动耐久性性定量计算的技术难题,为综合传动装置各个零部件的等寿命设计提供了技术基础。采用本发明使综合传动装置各个部件包含齿轮、轴、轴承等优化设计具有统一载荷谱,采用CAE技术分析能够得到关键部位和薄弱环节的应力变化规律,提高了计算的效率。建立综合传动装置耐久性评估数学模型,得到耐久性定量评估结果。

Description

基于载荷谱的综合传动装置耐久性评估方法
技术领域
本发明涉及履带车辆传动装置耐久性技术领域,具体涉及一种基于载荷谱的综合传动装置耐久性评估方法。
背景技术
综合传动装置是履带车辆常用的传动系统,一般有正差速,负差速和零差速传动系统。零差速液力机械综合传动装置是当前的主流传动型式。综合传动装置虽然结构复杂,其结构件分为齿轮类、轴类、轴承类部件,这些传动部件有机的协调工作,完成综合传动直驶、转向等各种功能和性能。当前关于综合传动耐久性评估基本都是按照实车试验,公开发表的文献或者专利没有从理论角度给出可行的计算方法。对于综合传动装置各个部件耐久性的计算往往是按照额定工况或者最大扭矩工况工作一定的时间,从而计算各个部件的耐久性,各个部件之间的载荷缺乏有机的匹配和组合,无法评估整机的耐久性。因缺乏统一的载荷谱,各个部件耐久性计算缺乏可比性,从而对整机等寿命设计等带来困难。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何提供一种基于载荷谱的综合传动耐久性评估方法,用以解决现有综合传动装置整机耐久性评估理论方法缺失的问题。
(二)技术方案
为了解决上述技术问题,本发明提供了一种基于载荷谱的综合传动装置耐久性评估方法,包括以下步骤:
步骤1、根据综合传动装置的结构特点,获得零部件工作载荷测试方案;
步骤2、通过测试获得步骤1的载荷统计数据,并进行初步处理,然后按照当前行驶的各种工况和里程,分别对齿轮、轴承和轴获得当前零件损伤的工作载荷谱;
步骤3、根据步骤2的结果,基于载荷谱计算综合传动装置各个部件的损伤率和寿命;
步骤4、基于步骤3的结果进行综合传动装置整机耐久性计算模型制定。
优选地,步骤1具体为:
根据预设测点,布置特定的传感器,测试4个转矩信号、风扇传动输入扭矩T1的时间历程、左侧汇流排输出扭矩T2的时间历程、右侧汇流排输出扭矩T3的时间历程、转向机构输出扭矩T4的时间历程、3个转速信号和1个档位信号,基于测试的数据通过内部传动链传动比的分析获得各个齿轮所传递的扭矩,并基于测试的数据通过对轴系传递功率的分析获得各个轴承、轴在特定工况所受的力,在试验过程中记录对应的数据;所述4个转矩信号包括风扇传动输入扭矩T1、左侧汇流排输出扭矩T2、右侧汇流排输出扭矩T3、转向机构输出扭矩T4;所述3个转速信号包括传动装置输入转速n1,转向零轴转速n2,变速机构输出转速n3
优选地,步骤2中按照当前行驶的各种工况和里程,分别对齿轮、轴承和轴获得当前零件损伤的工作载荷谱具体包括:
步骤2.1按照步骤1中记录的测试内容,对应相应的工况和传动比,计算所有齿轮受力情况:对于轴承的载荷谱统计,首先运用CAE软件建立3D模型,在有限元软件中施加稳态工况的传递力矩T,得到啮合齿轮齿根弯曲应力极值σF和齿面接触疲劳应力极值σH,对于润滑、速度、粗糙度、齿面硬化这些工艺措施对于应力的影响,分别引入润滑系数ZL、速度系数ZV、粗糙度系数ZR、齿面硬化系数ZW,对于接触疲劳应力σH进行修正,公式如式(1)所示,对于弯曲疲劳应力引入系数相对齿根表面状况系数进行修正,公式如式(2)所示:
σH′=σHZLZVZRZW (1)
σF′=σFYRrelT (2)
根据上述分析从而建立了齿轮传递转矩和应力的对应比例关系,将齿轮工作传递转矩和对应的转速按照时域对标,将传递转矩的均值进行8级等分,将对应的转速信号提取出来,再将转速信号8级等分,从而得到任意齿根弯曲疲劳和齿面接触疲劳应力区间下任意转速区间所对应的工作时间;
齿轮每旋转一周对应的轮齿啮合一次,也就意味着应力循环一次,按照式(3)计算对应的循环次数;
Nij=60bjgtij i,j=1,2,.....8 (3)
式中,Nij表示第i级载荷、第j级转速下的循环次数,bj表示第j级齿轮转速,gtij表示在第i级载荷、第j级转速下的齿轮的工作时间;
步骤2.2将实际的轴向和径向载荷等效为一假想的当量动载荷P来处理,在当量动载荷作用下,轴承的工作寿命与在实际工作载荷下的寿命相等,根据轴系的受力特点,计算受到的等效载荷如式(4)所示:
P=XFr+YFa (4)
式中,Fr表示轴承所承受的径向载荷;Fa表示轴承所承受的轴向载荷;X,Y表示径向载荷系数和轴向载荷系数;提取对应的转速,将不同转速区间和对应的等效均值载荷进行8级分级处理,得到轴承工作的载荷谱;
轴承寿命的计算公式按照下式(5)所示:
式中,Lhij为轴承在第i级当量动载荷,第j级转速区间下以90%可靠度工作的寿命,ci为轴承第i级当量载荷;dj为轴承第j级转速;ε为轴寿命系数:ft为温度系数;C为基本额定动载荷;
步骤2.3根据步骤1所测得的转矩,结合传动链传动比获得不同以传递转矩为主的轴的载荷,传递弯矩为主的轴通过轴系的受力分析,间接计算弯矩的变化规律,通过对稳态工况在CAE软件中获得传递转矩条件下的不同部位的应力分布,获得传递载荷和对应的极值应力的对应比例关系,将轴所受的应力开展雨流计数统计,获得应力均值和幅值的8级分布,得到任意应力均值和幅值条件所对应的循环次数,按照Goodman等效方法,将不同应力幅值条件下的应力均值转化为对称循环条件下的应力等效幅值。
优选地,步骤3具体包括:
步骤3.1齿轮耐久性计算:根据步骤2.1得到的载荷谱,按照Palmgren-Miner疲劳累积等效法则,依据齿轮材料对应的弯曲疲劳SN曲线,从而计算得到任意弯曲疲劳应力区间任意转速区间的疲劳损伤率,式(6)表示任意工况下的损伤率,式(7)表示所有工况下的总损伤率:
式中,Gij第ai级载荷,第bj级转速下的损伤率,Nij表示第i级载荷,第j级转速下的循环次数,Nai表示在齿轮弯曲疲劳SN曲线上,当应力均值为ai时最大循环次数,GZ表示总损伤率,如式(7)所示:
步骤3.2轴承耐久性计算:根据步骤2.2得到的载荷谱,以及已经获得任意轴承在特定的转速和转矩区间下工作的时间和可靠工作的工作寿命,按照式(8)获得任意工况的损伤率和总损伤率:
式中,Bij表示轴承在第ci级载荷,第dj级转速下的损伤率,tij表示轴承在第ci级载荷,第dj级转速下的工作时间,Lhij表示轴承在第ci级载荷,第dj级转速下以90%可靠度工作的寿命,BZ表示该轴承总的损伤率,如式(9):
步骤3.3轴耐久性计算:根据步骤2.2得到的载荷谱,按照Palmgren-Miner疲劳累积等效法则,对于承受扭转疲劳时效的轴采用扭转疲劳SN曲线,对于弯曲疲劳失效的轴采用弯曲疲劳SN曲线,从而得到轴在对应的应力均值和幅值条件下的损伤率,计算公式如式(10)和式(11)所示,Sij表示轴在第gj级等效应力幅值下的损伤率,hij代表在相应的应力均值和幅值条件下的应力循环次数,Ngi表示轴在对称循环应力gi下的允许疲劳循环次数,SZ表示该轴的总损伤率:
优选地,步骤4具体包括:
基于步骤3的结果,首先分别统计所有齿轮、轴和轴承的损伤率;
所有齿轮的损伤率G,用式(12)表示,所有轴的损伤率S,用式(13)表示,所有轴承的损伤率B,用式(14)表示:
G=max(GZ1,GZ2,.....GZm) (12)
S=max(SZ1,SZ2,.....SZm) (13)
B=max(BZ1,BZ2,.....BZm) (14)
用Q表示整个综合传动装置耐久性评价指标,当Q等于1的车辆行驶里程就是综合传动装置的寿命;
Q=max(G,B,S) (15)。
优选地,步骤2中,进行初步处理具体为:首先对于采集的载荷统计数据采用低通滤波法去掉干扰信号,然后采用幅值门限法去掉信号中的奇异值,并去除扭矩信号中损伤小于预设值的载荷。
(三)有益效果
与现有技术相比较,本发明具备如下有益效果:
1、本发明提供的基于载荷谱的综合传动装置耐久性评估计算方法,改变了以往综合传动整机耐久性计算没有统一载荷谱,各零部件耐久性计算无法关联匹配,整机耐久性评估没有方法参考的局面,从根本上解决了综合传动耐久性性定量计算的技术难题,为综合传动装置各个零部件的等寿命设计提供了技术基础。
2、采用本发明使综合传动装置各个部件包含齿轮、轴、轴承等优化设计具有统一载荷谱,采用CAE技术分析能够得到关键部位和薄弱环节的应力变化规律,提高了计算的效率。建立综合传动装置耐久性评估数学模型,得到耐久性定量评估结果。
附图说明
图1为综合传动测点布置图;
图2为特定转矩区间的转速数据提取图;
图3为某工况雨流矩阵示意图。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
本发明提供的一种基于载荷谱的综合传动装置耐久性评估方法,包括如下步骤:
步骤1、根据综合传动装置的结构特点,获得零部件工作载荷测试方案。
根据图1所示的测点,布置特定的传感器。测试4个转矩信号,风扇传动输入扭矩T1的时间历程(该时间历程指输入扭矩T1的过程所花费的时间)、左侧汇流排输出扭矩T2的时间历程、右侧汇流排输出扭矩T3的时间历程、转向机构输出扭矩T4的时间历程,3个转速信号,传动装置输入转速n1,转向零轴转速n2,变速机构输出转速n3,1个档位信号。基于测试的数据通过内部传动链传动比的分析就可以获得各个齿轮所传递的扭矩。基于测试的数据通过对轴系传递功率的分析可以获得各个轴承,轴在特定工况所受的力。在试验过程中记录对应的数据,如表1所示。
表1数据记录表
步骤2、通过测试获得步骤1的载荷统计数据,并进行初步处理。首先对于采集的载荷统计数据采用低通滤波法去掉干扰信号,从而提高信号的信噪比,所用的采用幅值门限法去掉信号中的奇异值,去除扭矩信号中损伤较小的载荷,减少信号处理的数据量。按照当前行驶的各种工况和里程,分别对齿轮、轴承和轴获得当前零件损伤的工作载荷谱。
2.1按照表1所示的测试内容,对应相应的工况和传动比,计算所有齿轮受力情况。对于轴承的载荷谱统计,首先运用CAE软件建立3D模型,在有限元软件中施加稳态工况的传递力矩T,得到啮合齿轮齿根弯曲应力极值σF和齿面接触疲劳应力极值σH。因为CAE求解齿面接触应力和齿根弯曲应力只能考虑材料和几何结构,对于润滑、速度、粗糙度、齿面硬化等工艺措施对于应力的影响,分别引入润滑系数ZL、速度系数ZV、粗糙度系数ZR、齿面硬化系数ZW,对于接触疲劳应力σH应该进行修正,公式如式(1)所示。对于弯曲疲劳应力引入系数相对齿根表面状况系数进行修正,公式如式(2)所示。这些系数的参考取值参照ISO 6336-2008。
σH′=σHZLZVZRZW (1)
σF′=σFYRrelT (2)
根据上述分析从而建立了齿轮传递转矩和应力的对应比例关系。将齿轮工作传递转矩和对应的转速按照时域对标,将传递转矩的均值进行8级等分,将对应的转速信号提取出来,如图2所示。再将转速信号8级等分,从而得到任意齿根弯曲疲劳和齿面接触疲劳应力区间下任意转速区间所对应的工作时间,如表2和表3所示,a1~a8代表弯曲(接触)疲劳应力的8级均值分级,b1~b8代表齿轮平均转速的8级分级,gt11~gt88代表在相应的应力均值和转速均值条件下的工作时间,N11~N88代表在相应的应力均值和转速均值条件下的循环次数。
齿轮每旋转一周对应的轮齿啮合一次,也就意味着应力循环次一次。按照式(3)计算对应的循环次数。
Nij=60bjgtij i,j=1,2,.....8 (3)
式中,Nij表示第i级载荷、第j级转速下的循环次数,bj表示第j级齿轮转速(单位rpm),gtij表示在第i级载荷、第j级转速下的齿轮的工作时间(单位h)。
表2齿轮工况统计表一
表3齿轮工况统计表二
2.2对于轴承的载荷谱统计,关键是获得轴承载荷和旋转次数的关系。轴承因包含内圈、外圈、保持架和滚动体,每一个部件的受力情况都不一致,按照ISO 281标准,将轴承旋转一圈视为应力循环一次,轴承在既承受轴向载荷又承受径向载荷,为了便于研究,将实际的轴向和径向载荷等效为一假想的当量动载荷P来处理,在当量动载荷作用下,轴承的工作寿命与在实际工作载荷下的寿命相等。根据轴系的受力特点,计算受到的等效载荷如式(4)所示。
P=XFr+YFa (4)
式中Fr表示轴承所承受的径向载荷(单位N);Fa表示轴承所承受的轴向载荷(单位N);X,Y表示径向载荷系数和轴向载荷系数,通过轴承手册可以查到。按照图2所示的方法提取对应的转速。将不同转速区间和对应的等效均值载荷进行8级分级处理,得到轴承工作的载荷谱,如表4所示,c1~c8代表载荷的8级均值分级,d1~d8代表轴承平均转速的8级分级,bt11~bt88代表在相应的转矩均值和转速均值条下的工作时间。
轴承寿命的计算公式按照下式(5)所示
式中,Lhij为轴承在第i级当量动载荷,第j级转速区间下以90%可靠度工作的寿命(h),ci为轴承第i级当量载荷,单位N;dj为轴承第j级转速,rpm;ε为轴寿命系数:对于球轴承ε=3,对于滚子轴承ε=10/3;ft为温度系数;C为基本额定动载荷。轴承寿命计算表如表5所示,Lh11~Lh88代表在相应的转矩均值和转速均值条件下的寿命。
表4轴承载荷谱统计表一
表5轴承载荷谱统计表二
2.3对于轴的载荷谱统计,关键是获得轴传递的载荷和载荷循环次数。综合传动装置重要的轴,例如传动主轴、转向零轴等主要传递转矩,齿轮支撑轴大部分以传递弯矩为主。根据表1所测得转矩,结合传动链传动比可以获得不同以传递转矩为主的轴的载荷。传递弯矩为主的轴可以通过轴系的受力分析,间接计算弯矩的变化规律。通过对稳态工况在CAE软件中获得传递转矩(弯矩)条件下的不同部位的应力分布,获得传递载荷和对应的极值应力的对应比例关系。将轴所受的应力开展雨流计数统计,获得应力均值和幅值的8级分布,得到任意应力均值和幅值条件所对应的循环次数,如表6所示。按照Goodman等效方法,将不同应力幅值条件下的应力均值转化为对称循环条件下的应力等效幅值,如表7所示。e1~e8代表载荷的8级均值分级,f1~f8代表轴所受应力幅值8级分级,h11~h88代表在相应的应力均值和幅值条件下的应力循环次数,g1~g8代表当量幅值的8级分级。
表6轴的载荷统计表
表7轴的等效载荷统计表
步骤3、根据步骤2的结果,基于载荷谱计算综合传动装置各个部件的损伤率和寿命。综合传动装置因为工况复杂,不同工况下所对应的轴承、轴和齿轮的寿命都不一致,因此采用损伤率来表示当前工作历程和时间中的载荷对于零件寿命的影响情况。
3.1齿轮耐久性计算。根据表3所示的载荷谱,按照Palmgren-Miner疲劳累积等效法则,依据齿轮材料对应的弯曲(接触)疲劳SN曲线,从而计算得到任意弯曲(接触)疲劳应力区间任意转速区间的疲劳损伤率。式(6)表示任意工况下的损伤率,式(7)表示所有工况下的总损伤率,损伤率计算表如表8所示。
式中,Gij第ai级载荷,第bj级转速下的损伤率,Nij表示第i级载荷,第j级转速下的循环次数,Nai表示在齿轮弯曲(接触)疲劳SN曲线上,当应力均值为ai时最大循环次数。GZ表示总损伤率,如式(7)所示。
表8齿轮损伤率计算表
3.2轴承耐久性计算。根据表4和表5所示的载荷谱,根据已经获得任意轴承在特定的转速和转矩区间下工作的时间和可靠工作的工作寿命,按照式(8)可以获得任意工况的损伤率和总损伤率,如表9所示。
式中,Bij表示轴承在第ci级载荷,第dj级转速下的损伤率,tij表示轴承在第ci级载荷,第dj级转速下的工作时间,Lhij表示轴承在第ci级载荷,第dj级转速下以90%可靠度工作的寿命。BZ表示该轴承总的损伤率,如式(9)所示。
表9轴承损伤率计算表
3.3轴耐久性计算。根据表7所示的载荷谱,按照Palmgren-Miner疲劳累积等效法则,对于承受扭转疲劳时效的轴采用扭转疲劳SN曲线,对于弯曲疲劳失效的轴采用弯曲疲劳SN曲线,从而得到轴在对应的应力均值和幅值条件下的损伤率,如表10所示,计算公式如式(10)和式(11)所示,Sij表示轴在第gj级等效应力幅值下的损伤率,hij代表在相应的应力均值和幅值条件下的应力循环次数,Ngi表示轴在对称循环应力gi下的允许疲劳循环次数,SZ表示该轴的总损伤率。
表10轴的损伤率计算表
步骤4、综合传动装置整机耐久性计算模型制定。综合传动装置组成结构复杂,轴系中任意齿轮、轴或者轴承断裂均会产生致命故障影响综合传动整机的耐久性。基于步骤3的结果,首先分别统计所有齿轮、轴和轴承的损伤率,如表11所示,GZ1~GZm分别代表m个齿轮的损伤率,SZ1~SZn分别代表n个轴的损伤率,BZ1~BZp分别代表p个轴承的损伤率。
所有齿轮的损伤率G,可用式(12)表示,所有轴的损伤率S,可用式(13)表示,所有轴承的损伤率B,可用式(14)表示。
G=max(GZ1,GZ2,.....GZm) (12)
S=max(SZ1,SZ2,.....SZm) (13)
B=max(BZ1,BZ2,.....BZm) (14)
Q表示整个综合传动装置耐久性评价指标,当Q等于1的车辆行驶里程就是综合传动装置的寿命。
Q=max(G,B,S) (15)
表11损伤率表
项目 损伤率 项目 损伤率 项目 损伤率
齿轮1 GZ1 轴1 SZ1 轴承1 BZ1
齿轮2 GZ2 轴2 SZ2 轴承2 BZ2
....... ....... ....... ....... ....... .......
齿轮m GZm 轴n SZn 轴承p BZp
可以看出,本发明应用CAE分析技术手段得到各个部件薄弱部位的应力和变化规律,基于雨流计数法得到各个部件载荷循环的载荷谱,为耐久性计算提供了准确的边界条件。综合传动装置结构复杂任何部件的损坏都会造成整机的功能丧失,基于工作特点提出综合传动耐久性计算的模型。本发明中综合传动装置各部件耐久性计算所采用的载荷谱主要依据对车辆的道路实测和CAE分析等手段获得,计算结果置信度较高,提供了综合传动耐久性定量评估的手段。如果车辆在对应的试验跑道循环行驶,根据特定里程所获得综合传动整机损伤率,就可以获得损伤率为1时,综合传动装置当量寿命里程,可用于在线寿命预计。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。

Claims (6)

1.一种基于载荷谱的综合传动装置耐久性评估方法,其特征在于,包括以下步骤:
步骤1、根据综合传动装置的结构特点,获得零部件工作载荷测试方案;
步骤2、通过测试获得步骤1的载荷统计数据,并进行初步处理,然后按照当前行驶的各种工况和里程,分别对齿轮、轴承和轴获得当前零件损伤的工作载荷谱;
步骤3、根据步骤2的结果,基于载荷谱计算综合传动装置各个部件的损伤率和寿命;
步骤4、基于步骤3的结果进行综合传动装置整机耐久性计算模型制定。
2.如权利要求1所述的方法,其特征在于,步骤1具体为:
根据预设测点,布置特定的传感器,测试4个转矩信号、风扇传动输入扭矩T1的时间历程、左侧汇流排输出扭矩T2的时间历程、右侧汇流排输出扭矩T3的时间历程、转向机构输出扭矩T4的时间历程、3个转速信号和1个档位信号,基于测试的数据通过内部传动链传动比的分析获得各个齿轮所传递的扭矩,并基于测试的数据通过对轴系传递功率的分析获得各个轴承、轴在特定工况所受的力,在试验过程中记录对应的数据;所述4个转矩信号包括风扇传动输入扭矩T1、左侧汇流排输出扭矩T2、右侧汇流排输出扭矩T3、转向机构输出扭矩T4;所述3个转速信号包括传动装置输入转速n1,转向零轴转速n2,变速机构输出转速n3
3.如权利要求2所述的方法,其特征在于,步骤2中按照当前行驶的各种工况和里程,分别对齿轮、轴承和轴获得当前零件损伤的工作载荷谱具体包括:
步骤2.1按照步骤1中记录的测试内容,对应相应的工况和传动比,计算所有齿轮受力情况:对于轴承的载荷谱统计,首先运用CAE软件建立3D模型,在有限元软件中施加稳态工况的传递力矩T,得到啮合齿轮齿根弯曲应力极值σF和齿面接触疲劳应力极值σH,对于润滑、速度、粗糙度、齿面硬化这些工艺措施对于应力的影响,分别引入润滑系数ZL、速度系数ZV、粗糙度系数ZR、齿面硬化系数ZW,对于接触疲劳应力σH进行修正,公式如式(1)所示,对于弯曲疲劳应力引入系数相对齿根表面状况系数进行修正,公式如式(2)所示:
σH′=σHZLZVZRZW (1)
σF′=σFYRrelT (2)
根据上述分析从而建立了齿轮传递转矩和应力的对应比例关系,将齿轮工作传递转矩和对应的转速按照时域对标,将传递转矩的均值进行8级等分,将对应的转速信号提取出来,再将转速信号8级等分,从而得到任意齿根弯曲疲劳和齿面接触疲劳应力区间下任意转速区间所对应的工作时间;
齿轮每旋转一周对应的轮齿啮合一次,也就意味着应力循环一次,按照式(3)计算对应的循环次数;
Nij=60bjgtij i,j=1,2,.....8 (3)
式中,Nij表示第i级载荷、第j级转速下的循环次数,bj表示第j级齿轮转速,gtij表示在第i级载荷、第j级转速下的齿轮的工作时间;
步骤2.2将实际的轴向和径向载荷等效为一假想的当量动载荷P来处理,在当量动载荷作用下,轴承的工作寿命与在实际工作载荷下的寿命相等,根据轴系的受力特点,计算受到的等效载荷如式(4)所示:
P=XFr+YFa (4)
式中,Fr表示轴承所承受的径向载荷;Fa表示轴承所承受的轴向载荷;X,Y表示径向载荷系数和轴向载荷系数;提取对应的转速,将不同转速区间和对应的等效均值载荷进行8级分级处理,得到轴承工作的载荷谱;
轴承寿命的计算公式按照下式(5)所示:
式中,Lhij为轴承在第i级当量动载荷,第j级转速区间下以90%可靠度工作的寿命,ci为轴承第i级当量载荷;dj为轴承第j级转速;ε为轴寿命系数:ft为温度系数;C为基本额定动载荷;
步骤2.3根据步骤1所测得的转矩,结合传动链传动比获得不同以传递转矩为主的轴的载荷,传递弯矩为主的轴通过轴系的受力分析,间接计算弯矩的变化规律,通过对稳态工况在CAE软件中获得传递转矩条件下的不同部位的应力分布,获得传递载荷和对应的极值应力的对应比例关系,将轴所受的应力开展雨流计数统计,获得应力均值和幅值的8级分布,得到任意应力均值和幅值条件所对应的循环次数,按照Goodman等效方法,将不同应力幅值条件下的应力均值转化为对称循环条件下的应力等效幅值。
4.如权利要求3所述的方法,其特征在于,步骤3具体包括:
步骤3.1齿轮耐久性计算:根据步骤2.1得到的载荷谱,按照Palmgren-Miner疲劳累积等效法则,依据齿轮材料对应的弯曲疲劳SN曲线,从而计算得到任意弯曲疲劳应力区间任意转速区间的疲劳损伤率,式(6)表示任意工况下的损伤率,式(7)表示所有工况下的总损伤率:
式中,Gij第ai级载荷,第bj级转速下的损伤率,Nij表示第i级载荷,第j级转速下的循环次数,Nai表示在齿轮弯曲疲劳SN曲线上,当应力均值为ai时最大循环次数,GZ表示总损伤率,如式(7)所示:
步骤3.2轴承耐久性计算:根据步骤2.2得到的载荷谱,以及已经获得任意轴承在特定的转速和转矩区间下工作的时间和可靠工作的工作寿命,按照式(8)获得任意工况的损伤率和总损伤率:
式中,Bij表示轴承在第ci级载荷,第dj级转速下的损伤率,tij表示轴承在第ci级载荷,第dj级转速下的工作时间,Lhij表示轴承在第ci级载荷,第dj级转速下以90%可靠度工作的寿命,BZ表示该轴承总的损伤率,如式(9):
步骤3.3轴耐久性计算:根据步骤2.2得到的载荷谱,按照Palmgren-Miner疲劳累积等效法则,对于承受扭转疲劳时效的轴采用扭转疲劳SN曲线,对于弯曲疲劳失效的轴采用弯曲疲劳SN曲线,从而得到轴在对应的应力均值和幅值条件下的损伤率,计算公式如式(10)和式(11)所示,Sij表示轴在第gj级等效应力幅值下的损伤率,hij代表在相应的应力均值和幅值条件下的应力循环次数,Ngi表示轴在对称循环应力gi下的允许疲劳循环次数,SZ表示该轴的总损伤率:
5.如权利要求6所述的方法,其特征在于,步骤4具体包括:
基于步骤3的结果,首先分别统计所有齿轮、轴和轴承的损伤率;
所有齿轮的损伤率G,用式(12)表示,所有轴的损伤率S,用式(13)表示,所有轴承的损伤率B,用式(14)表示:
G=max(GZ1,GZ2,.....GZm) (12)
S=max(SZ1,SZ2,.....SZm) (13)
B=max(BZ1,BZ2,.....BZm) (14)
用Q表示整个综合传动装置耐久性评价指标,当Q等于1的车辆行驶里程就是综合传动装置的寿命;
Q=max(G,B,S) (15)。
6.如权利要求2所述的方法,其特征在于,步骤2中,进行初步处理具体为:首先对于采集的载荷统计数据采用低通滤波法去掉干扰信号,然后采用幅值门限法去掉信号中的奇异值,并去除扭矩信号中损伤小于预设值的载荷。
CN201810595308.3A 2018-06-11 2018-06-11 基于载荷谱的综合传动装置耐久性评估方法 Pending CN109030016A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810595308.3A CN109030016A (zh) 2018-06-11 2018-06-11 基于载荷谱的综合传动装置耐久性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810595308.3A CN109030016A (zh) 2018-06-11 2018-06-11 基于载荷谱的综合传动装置耐久性评估方法

Publications (1)

Publication Number Publication Date
CN109030016A true CN109030016A (zh) 2018-12-18

Family

ID=64612510

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810595308.3A Pending CN109030016A (zh) 2018-06-11 2018-06-11 基于载荷谱的综合传动装置耐久性评估方法

Country Status (1)

Country Link
CN (1) CN109030016A (zh)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109977459A (zh) * 2019-02-11 2019-07-05 中国第一汽车股份有限公司 一种应用cae技术提高强化载荷谱精度的方法
CN110823561A (zh) * 2019-11-13 2020-02-21 南京工业大学 一种基于材料模型库的非稳态载荷下轧机传动轴系统关键件疲劳寿命实时监控的方法
CN111896255A (zh) * 2020-08-12 2020-11-06 上海理工大学 基于多轴随机道路载荷下轮毂轴承服役寿命快速评估方法
CN111914370A (zh) * 2020-08-07 2020-11-10 上海理工大学 一种电驱动系统高速轴承可靠性试验载荷谱的编制方法
CN112307630A (zh) * 2020-11-02 2021-02-02 中国航发沈阳发动机研究所 一种航空发动机主轴承寿命分析载荷谱编制方法
CN112613198A (zh) * 2021-03-08 2021-04-06 中国空气动力研究与发展中心低速空气动力研究所 一种去除风洞风扇干扰的数据处理方法
CN113704880A (zh) * 2021-08-07 2021-11-26 中国航空工业集团公司沈阳飞机设计研究所 一种静振联合加载的疲劳谱编制方法
CN114199431A (zh) * 2021-12-07 2022-03-18 上海中联重科桩工机械有限公司 卷扬减速机的载荷谱获取方法及系统和旋挖钻机
CN114925468A (zh) * 2022-05-06 2022-08-19 中国北方车辆研究所 一种基于历史数据的修正系数可靠性分配方法
CN115270583A (zh) * 2022-09-19 2022-11-01 麦格纳动力总成(江西)有限公司 汽车差速器扭转冲击疲劳分析方法及系统
CN115659763A (zh) * 2022-12-08 2023-01-31 麦格纳动力总成(江西)有限公司 一种轴承和轴承座耐久分析方法
CN116337484A (zh) * 2023-05-30 2023-06-27 麦格纳动力总成(江西)有限公司 驱动桥耐久试验方法、系统、计算机及可读存储介质
CN116735196A (zh) * 2023-08-16 2023-09-12 盛瑞传动股份有限公司 针对轴承损伤的变速箱强化载荷谱试验方法
CN117932987A (zh) * 2024-03-25 2024-04-26 华南农业大学 一种电动汽车齿轮箱疲劳载荷计算方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1777647A1 (en) * 2005-10-13 2007-04-25 Ford Global Technologies, LLC Transmission life prediction method
CN104101548A (zh) * 2013-04-09 2014-10-15 中国人民解放军第二炮兵工程大学 一种适用于低成本的无人机机体结构定寿方法
CN106644464A (zh) * 2016-11-18 2017-05-10 南京工业大学 一种基于载荷谱分析的轧机传动系统关键零部件的疲劳寿命预警方法
CN106886638A (zh) * 2017-01-22 2017-06-23 北京理工大学 一种基于核密度估计的履带车辆传动轴载荷谱编制方法
CN107590311A (zh) * 2017-08-08 2018-01-16 明阳智慧能源集团股份公司 一种行星传动齿轮箱强度计算集成系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1777647A1 (en) * 2005-10-13 2007-04-25 Ford Global Technologies, LLC Transmission life prediction method
CN104101548A (zh) * 2013-04-09 2014-10-15 中国人民解放军第二炮兵工程大学 一种适用于低成本的无人机机体结构定寿方法
CN106644464A (zh) * 2016-11-18 2017-05-10 南京工业大学 一种基于载荷谱分析的轧机传动系统关键零部件的疲劳寿命预警方法
CN106886638A (zh) * 2017-01-22 2017-06-23 北京理工大学 一种基于核密度估计的履带车辆传动轴载荷谱编制方法
CN107590311A (zh) * 2017-08-08 2018-01-16 明阳智慧能源集团股份公司 一种行星传动齿轮箱强度计算集成系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘新猛: "汽车变速器关键部件疲劳寿命预测与软件开发", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
刘永臣 等: "基于实测载荷谱的副车架疲劳寿命估算方法", 《汽车技术》 *

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109977459A (zh) * 2019-02-11 2019-07-05 中国第一汽车股份有限公司 一种应用cae技术提高强化载荷谱精度的方法
CN110823561A (zh) * 2019-11-13 2020-02-21 南京工业大学 一种基于材料模型库的非稳态载荷下轧机传动轴系统关键件疲劳寿命实时监控的方法
CN110823561B (zh) * 2019-11-13 2021-03-26 南京工业大学 一种基于材料模型库的非稳态载荷下轧机传动轴系统关键件疲劳寿命实时监控的方法
CN111914370A (zh) * 2020-08-07 2020-11-10 上海理工大学 一种电驱动系统高速轴承可靠性试验载荷谱的编制方法
CN111914370B (zh) * 2020-08-07 2021-03-02 上海理工大学 一种电驱动系统高速轴承可靠性试验载荷谱的编制方法
US11835419B2 (en) 2020-08-12 2023-12-05 University Of Shanghai For Science And Technology Method for rapid assessment of service life of hub bearing based on multiaxial random road loading
CN111896255A (zh) * 2020-08-12 2020-11-06 上海理工大学 基于多轴随机道路载荷下轮毂轴承服役寿命快速评估方法
CN112307630A (zh) * 2020-11-02 2021-02-02 中国航发沈阳发动机研究所 一种航空发动机主轴承寿命分析载荷谱编制方法
CN112613198A (zh) * 2021-03-08 2021-04-06 中国空气动力研究与发展中心低速空气动力研究所 一种去除风洞风扇干扰的数据处理方法
CN113704880A (zh) * 2021-08-07 2021-11-26 中国航空工业集团公司沈阳飞机设计研究所 一种静振联合加载的疲劳谱编制方法
CN113704880B (zh) * 2021-08-07 2023-12-15 中国航空工业集团公司沈阳飞机设计研究所 一种静振联合加载的疲劳谱编制方法
CN114199431A (zh) * 2021-12-07 2022-03-18 上海中联重科桩工机械有限公司 卷扬减速机的载荷谱获取方法及系统和旋挖钻机
CN114925468A (zh) * 2022-05-06 2022-08-19 中国北方车辆研究所 一种基于历史数据的修正系数可靠性分配方法
CN115270583B (zh) * 2022-09-19 2022-12-27 麦格纳动力总成(江西)有限公司 汽车差速器扭转冲击疲劳分析方法及系统
CN115270583A (zh) * 2022-09-19 2022-11-01 麦格纳动力总成(江西)有限公司 汽车差速器扭转冲击疲劳分析方法及系统
CN115659763A (zh) * 2022-12-08 2023-01-31 麦格纳动力总成(江西)有限公司 一种轴承和轴承座耐久分析方法
CN115659763B (zh) * 2022-12-08 2023-03-10 麦格纳动力总成(江西)有限公司 一种轴承和轴承座耐久分析方法
CN116337484A (zh) * 2023-05-30 2023-06-27 麦格纳动力总成(江西)有限公司 驱动桥耐久试验方法、系统、计算机及可读存储介质
CN116337484B (zh) * 2023-05-30 2023-09-26 麦格纳动力总成(江西)有限公司 驱动桥耐久试验方法、系统、计算机及可读存储介质
CN116735196A (zh) * 2023-08-16 2023-09-12 盛瑞传动股份有限公司 针对轴承损伤的变速箱强化载荷谱试验方法
CN116735196B (zh) * 2023-08-16 2023-10-31 盛瑞传动股份有限公司 针对轴承损伤的变速箱强化载荷谱试验方法
CN117932987A (zh) * 2024-03-25 2024-04-26 华南农业大学 一种电动汽车齿轮箱疲劳载荷计算方法
CN117932987B (zh) * 2024-03-25 2024-06-28 华南农业大学 一种电动汽车齿轮箱疲劳载荷计算方法

Similar Documents

Publication Publication Date Title
CN109030016A (zh) 基于载荷谱的综合传动装置耐久性评估方法
Zhao et al. Deep convolutional neural network based planet bearing fault classification
CN109000922B (zh) 基于道路测试的综合传动装置台架耐久性试验方法
CN100573085C (zh) 滚动轴承的剩余寿命诊断方法及剩余寿命诊断装置
CN111896255A (zh) 基于多轴随机道路载荷下轮毂轴承服役寿命快速评估方法
CN104297148B (zh) 一种润滑摩擦特征测量分析系统及其工作方法
Barbieri et al. Analysis of automotive gearbox faults using vibration signal
CN112326236B (zh) 一种齿轮箱运行状态在线监测方法、系统及存储介质
CN105701348A (zh) 一种谐波减速器可靠寿命预测方法
CN112364706A (zh) 一种基于类不平衡的小样本轴承故障诊断方法
CN107101827A (zh) 一种低速重载齿轮裂纹故障在线检测方法
Cao et al. Deterioration state diagnosis and wear evolution evaluation of planetary gearbox using vibration and wear debris analysis
CN110056640A (zh) 基于加速度信号和边缘计算的减速机无线故障诊断方法
CN105571994B (zh) 一种动力传动系统滑油金属屑参数检测方法
CN110646224A (zh) 车辆驱动轮瞬间输出功率的测量方法
CN207066742U (zh) 一种汽车变速器轴承试验台
CN109883692B (zh) 基于内置编码器信息的广义差分滤波方法
CN115586023A (zh) 一种轨道车辆传动系统故障诊断方法及系统
CN110196165B (zh) 一种基于K-means聚类和证据推理的旋转机械滚珠轴承故障诊断方法
CN115879245A (zh) 一种中碳钢齿轮长寿命传动疲劳可靠性评价方法
CN106053275A (zh) 一种摩擦副磨损过程摩擦振动检测方法
CN116007930B (zh) 一种汽车传动性能的测试方法及系统
CN110160781B (zh) 一种旋转机械故障分类的测试集重构和预测方法
CN102592172A (zh) 一种基于敏感测点数据融合的行星轮故障识别方法
CN116499742A (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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20181218