CN113569332B - 托带轮疲劳可靠度计算系统、方法、履带车辆 - Google Patents
托带轮疲劳可靠度计算系统、方法、履带车辆 Download PDFInfo
- Publication number
- CN113569332B CN113569332B CN202110673262.4A CN202110673262A CN113569332B CN 113569332 B CN113569332 B CN 113569332B CN 202110673262 A CN202110673262 A CN 202110673262A CN 113569332 B CN113569332 B CN 113569332B
- Authority
- CN
- China
- Prior art keywords
- working condition
- unit
- mileage
- damage amount
- fatigue
- 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
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 230000006378 damage Effects 0.000 claims description 88
- 238000009826 distribution Methods 0.000 claims description 40
- 230000004044 response Effects 0.000 claims description 19
- 238000001228 spectrum Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000008458 response to injury Effects 0.000 abstract description 2
- 238000003466 welding Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft 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
-
- 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明公开了一种托带轮疲劳可靠度计算系统、方法、履带车辆,应用于履带车辆技术领域,基于工况分散性,通过多体动力学仿真模型、损伤响应面模型和MATLAB程序计算履带车辆托带轮的疲劳可靠度,解决了传统车辆结构疲劳可靠度预计方法无法基于工况强随机性和载荷强非线性的问题,该方法计算简便、易懂、运算量小,成本低,具有较强工程适用性。
Description
技术领域
本发明涉及履带车辆技术领域,具体涉及一种履带车辆托带轮疲劳可靠度计算系统、方法、履带车辆。
背景技术
履带车辆托带轮属于履带车辆行动系统,其功能是通过托带轮支架安装在车体侧甲板上,用于托举履带上支段,以保持车辆行驶过程中履带的总体形态。履带车辆行驶过程中,履带受路面不平度激励,对托带轮施加垂直于地面方向的交变载荷。在交变载荷的作用下,由于表面粗糙度、托带轮本身的缺陷和设计机构等原因,托带轮支架与车体侧甲板连接位置存在应力集中的部位将产生损伤并逐渐积累导致裂纹的产生、裂纹扩展直到造成托带轮断裂从而发生疲劳损伤失效,直接影响车体可靠性。
疲劳损伤失效对零部件的使用条件和外部环境非常敏感。在履带车辆行驶过程中,复合路面不平度与车辆结构载荷具有一致的随机性,它是车辆结构产生疲劳损伤的重要原因。但随着车速的提高和路面等级的降低,路面激励引起的路面动位移逐渐增加,而增加趋势是非线性的。对此,现有技术中存在一种结合粒子群算法和模糊数学理论对系统的可靠性进行研究的方法,或者采用多体动力学仿真与有限元分析相结合的方法对履带车辆结构进行疲劳分析,但是其仍然无法解决工况强随机性和载荷强非线性的问题,而且计算方式复杂,运算量大。因此,需要提出一种解决上述问题的托带轮支架疲劳可靠度计算方法及系统。
发明内容
有鉴于此,本发明提供了一种托带轮疲劳可靠度计算系统、方法、履带车辆,能够解决传统车辆结构疲劳可靠度预计方法无法基于工况强随机性和载荷强非线性的问题。
本发明采用的技术方案如下:
一种履带车辆托带轮疲劳可靠度计算系统,包括第一确定单元、第二确定单元、第三确定单元、模型建立单元、统计分布确定单元和可靠度计算单元,其中:
所述第一确定单元用于检测及确定给定任务剖面涉及的不同工况类型及不同工况类型下的里程占比分布;
所述第二确定单元用于基于履带车辆托带轮结构疲劳机理分析及工况参数的显著性分析,确定履带车辆托带轮的薄弱单元及对应的工况敏感参数;
所述第三确定单元用于在基于工况分散性的前提下,利用履带车辆托带轮的多体动力学仿真模型,获取履带车辆托带轮载荷谱并确定履带车辆托带轮的薄弱单元单位里程疲劳损伤量;
所述模型建立单元用于基于所述履带车辆托带轮的薄弱单元单位里程疲劳损伤量与所述工况敏感参数确定二者之间的定量关系,建立所述不同工况类型下的响应面模型;
所述统计分布确定单元用于对所述工况敏感参数进行离散化处理,确定所述不同工况类型下的工况敏感参数对应的统计分布;
所述可靠度计算单元用于从所述工况敏感参数对应的统计分布中随机抽取n个样本,代入响应面模型中确定履带车辆托带轮的单位里程疲劳损伤量,基于所述单位里程疲劳损伤量计算所述任务剖面下的累积疲劳损伤量,将所述累积疲劳损伤量与损伤阈值进行比较,然后经过n次抽取和n次比较,统计未失效样本数量,计算所述任务剖面下的履带车辆托带轮疲劳可靠度。
进一步地,给定总里程数,根据所述里程占比分布确定不同工况类型下的里程数;所述工况类型为路面等级。
进一步地,所述工况敏感参数为在对所述薄弱单元的疲劳损伤影响最大的工况影响因素。
进一步地,所述第三确定单元根据所述载荷谱进行雨流计数法,确定所述履带车辆托带轮的薄弱单元单位里程疲劳损伤量。
进一步地,所述响应面模型为变量为不同路面等级下的履带车辆托带轮单位里程疲劳损伤量与行驶速度的多项式响应面模型;所述不同路面等级包括4种,分别为C级路面、D级路面、E级路面和F级路面,其中,
Qc=1.991×10-7×Vc-5.12×10-6
式中,Qc表示C级路面下单位里程疲劳损伤量,Vc表示C级路面下履带车辆行驶速度;
Qd=3.256×10-8×Vd 2-2.641×10-6×Vd+5.675×10-5
式中,Qd表示D级路面下单位里程疲劳损伤量,Vd表示D级路面下履带车辆行驶速度;
Qe=-5.539×10-10×Ve 3+4.723×10-8×Ve 2-3.079×10-5
式中,Qe表示E级路面下单位里程疲劳损伤量,Ve表示E级路面下履带车辆行驶速度;
Qf=-4.727×10-9×Vf 3+6.993×10-7×Vf 2-3.313×10-5×Vf+0.0005293
式中,Qf表示F级路面下单位里程疲劳损伤量,Vf表示F级路面下履带车辆行驶速度,其中,行驶速度单位为km/h。
进一步地,所述任务剖面下的累积疲劳损伤量为Qmi,
Qmi=Yci×Qci+Ydi×Qdi+Yei×Qei+Yfi×Qfi
式中,Yci表示第i个样本在C级路面的里程,Qci表示第i个样本在C级路面的单位里程损伤量;Tdi表示第i个样本在D级路面的里程,Qdi表示第i个样本在D级路面的单位里程损伤量;Yei表示第i个样本在E级路面的里程,Qei表示第i个样本在E级路面的单位里程损伤量;Yfi表示第i个样本在F级路面的里程,Qfi表示第i个样本在F级路面的单位里程损伤量。
进一步地,所述履带车辆托带轮疲劳可靠度为R,
R=1-m/M
其中,m为所述履带车辆托带轮疲劳失效样本数,M为根据所述工况敏感参数统计分布进行抽样的次数。
一种基于所述的履带车辆托带轮疲劳可靠度计算系统的计算方法,该方法具体步骤如下:
步骤S101:检测及确定给定任务剖面涉及的不同工况类型及不同工况类型下的里程占比分布;
步骤S102:基于履带车辆托带轮结构疲劳机理分析及工况参数的显著性分析,确定履带车辆托带轮的薄弱单元及对应的工况敏感参数;
步骤S103:在基于工况分散性的前提下,利用履带车辆托带轮的多体动力学仿真模型,获取履带车辆托带轮载荷谱并确定履带车辆托带轮的薄弱单元单位里程疲劳损伤量;
步骤S104:基于所述履带车辆托带轮的薄弱单元单位里程疲劳损伤量与所述工况敏感参数确定二者之间的定量关系,建立不同工况类型下的响应面模型;
步骤S105:对所述工况敏感参数进行离散化处理,确定不同工况类型下工况敏感参数对应的统计分布;
步骤S106:通过MATLAB程序从所述工况敏感参数统计分布中随机抽取n个样本,代入响应面模型中确定履带车辆托带轮的单位里程疲劳损伤量,基于所述单位里程疲劳损伤量计算任务剖面下的累积疲劳损伤量,将所述累积疲劳损伤量与损伤阈值进行比较,然后经过n次抽取和n次比较,统计未失效样本数量,计算所述任务剖面下的履带车辆托带轮疲劳可靠度。
一种履带车辆,所述履带车辆包括托带轮及如上述的履带车辆托带轮疲劳可靠度计算系统。
有益效果:
1、本发明基于工况分散性,通过多体动力学仿真模型、损伤响应面模型和MATLAB程序计算履带车辆托带轮的疲劳可靠度,解决了传统车辆结构疲劳可靠度预计方法无法基于工况强随机性和载荷强非线性的问题,该方法计算简便、易懂、运算量小,成本低,具有较强工程适用性。
2、本发明提供了托带轮类产品任务疲劳可靠度计算模型,为同类产品任务疲劳可靠度计算提供了方法依据,该方法科学合理、可操作性强,可直接用于工程实践。
附图说明
图1为履带车辆托带轮结构组成示意图。
图2是本发明疲劳可靠度计算系统组成框图。
图3是本发明疲劳可靠度计算方法流程图。
图4是本发明任务可靠度抽样仿真流程图。
其中,1-托带轮总成、2-托带轮悬臂、3-托带轮支架、4-车体侧甲板、5-焊缝、101-第一确定单元、102-第二确定单元、103-第三确定单元、104-模型建立单元、105-统计分布确定单元、106-可靠度计算单元。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明的托带轮疲劳可靠度计算系统以履带车辆托带轮为对象,如图1所示,履带车辆托带轮结构包括托带轮总成1、托带轮悬臂2、托带轮支架3、车体侧甲板4和托带轮支架3与车体侧甲板4之间的焊缝5。
一种履带车辆托带轮疲劳可靠度计算系统,基于工况分散性,应用在履带车辆托带轮控制装置中。履带车辆包括托带轮及履带车辆托带轮疲劳可靠度计算系统。本实施例还提供了一种计算机程序产品,包括用于执行该履带车辆托带轮疲劳可靠度计算系统的计算机程序代码。
具体地,如图2所示,托带轮疲劳可靠度计算系统包括第一确定单元101、第二确定单元102、第三确定单元103、模型建立单元104、统计分布确定单元105、可靠度计算单元106,其中:
第一确定单元101,用于检测及确定给定任务剖面涉及的不同工况类型及不同工况类型下的里程占比分布。
第二确定单元102,用于基于履带车辆托带轮结构疲劳机理分析及工况参数的显著性分析,确定履带车辆托带轮的薄弱单元及对应的工况敏感参数。
第三确定单元103,用于在基于工况分散性的前提下,利用履带车辆托带轮的多体动力学仿真模型,获取履带车辆托带轮载荷谱并确定履带车辆托带轮的薄弱单元单位里程疲劳损伤量。
模型建立单元104,用于基于履带车辆托带轮的薄弱单元单位里程疲劳损伤量与工况敏感参数确定二者之间的定量关系,建立不同工况类型下的响应面模型。
统计分布确定单元105,用于对工况敏感参数进行离散化处理,确定不同工况类型下工况敏感参数对应的统计分布。
可靠度计算单元106,用于通过MATLAB程序从工况敏感参数统计分布中随机抽取n个样本,代入响应面模型中确定履带车辆托带轮的单位里程疲劳损伤量,基于单位里程疲劳损伤量计算任务剖面下的累积疲劳损伤量,将累积疲劳损伤量与损伤阈值进行比较,然后经过n次抽取和n次比较,统计未失效样本数量,计算该任务剖面下的履带车辆托带轮疲劳可靠度。
此外,本发明还提供了一种基于工况分散性的履带车辆托带轮疲劳可靠度计算方法,如图3所示,该方法具体步骤如下:
步骤S101:任务剖面分析:检测及确定给定任务剖面涉及的不同工况类型及不同工况类型下的里程占比分布,如图3所示,某履带车辆托带轮任务剖面的分析结果如表1所示,不同工况类型包括4种路面等级,分别为C级路面、D级路面、E级路面和F级路面。四种路面等级按面层材料的组成、结构强度、路面所能承担的交通任务和使用品质进行划分,可以为高级路面、次高级路面、中级路面和低级路面,其中C-F级路面与划分的路面等级相对应。给定总里程数,根据里程占比分布确定不同工况类型下的里程数,本实施例中,任务剖面总里程为10000km,C级路面里程数为1000km,D级路面里程数为3000km,E级路面里程数为3000km,F级路面里程数为3000km。通过不同工况类型下的里程数考虑工况分散性,避免因工况强随机性影响疲劳可靠度计算。
表1某履带车辆托带轮任务剖面分析表
步骤S102:工况敏感参数分析:在步骤S101的基础上,结合履带车辆托带轮的性能参数和结构原理,托带轮焊缝5处在不同路况下均易出现局部应力集中,且断裂频率高,即确定托带轮焊缝5处为薄弱单元,易发生疲劳断裂。基于薄弱单元出发计算托带轮的疲劳可靠度,保证托带轮整体结构的可靠性,避免因部分结构失效影响整体结构可靠性。在基础上,进行工况显著性分析,分析结果显示,路面等级和行驶速度对托带轮焊缝5处疲劳损伤影响最为剧烈,因此,将路面等级和行驶速度作为工况敏感参数。
步骤S103:疲劳损伤计算:在步骤S102的基础上,建立履带车辆托带轮的多体动力学模型,获取履带车辆托带轮载荷谱并确定履带车辆托带轮的薄弱单元单位里程疲劳损伤量。以总里程10000km为前提,开展全因子仿真试验,并计算各试验条件下履带车辆托带轮焊缝5处的单位里程疲劳损伤量,结果如表2所示。
表2履带车辆车体托带轮焊缝处疲劳损伤量仿真结果
步骤S104:响应面模型构建:在步骤S103的基础上,确定单位里程疲劳损伤量与工况敏感参数之间的定量关系,构建对应工况类型下的响应面模型,具体地即路面等级下的响应面模型。
根据表2仿真结果进行多项式拟合,进一步获得各路面等级下的变量为履带车辆托带轮单位里程疲劳损伤量与行驶速度的多项式响应面模型。
Qc=1.991×10-7×Vc-5.12×10-6 (1)
式中,Qc表示C级路面下单位里程疲劳损伤量,Vc表示C级路面下履带车辆行驶速度,单位为km/h。
Qd=3.256×10-8×Vd 2-2.641×10-6×Vd+5.675×10-5 (2)
式中,Qd表示D级路面下单位里程疲劳损伤量,Vd表示D级路面下履带车辆行驶速度,单位为km/h。
Qe=-5.539×10-10×Ve 3+4.723×10-8×Ve 2-3.079×10-5 (3)
式中,Qe表示E级路面下单位里程疲劳损伤量,Ve表示E级路面下履带车辆行驶速度,单位为km/h。
Qf=-4.727×10-9×Vf 3+6.993×10-7×Vf 2-3.313×10-5×Vf+0.0005293 (4)
式中,Qf表示F级路面下单位里程疲劳损伤量,Vf表示F级路面下履带车辆行驶速度,单位为km/h。
步骤S105:工况随机分布确定:在步骤S104的基础上,将工况敏感参数离散化处理,选取不同工况类型下的对应的统计分布。即对一个给定的参数值,赋予它一个离散系数,形成包含此参数值的一个参数域,进行离散化处理后确定一个随机分布的相关参数值。下面假设履带车辆在该地域行驶时,在不同路面级别下行驶速度均服从正态分布,给出各路面等级的统计分布的示意如下:
C级路面下,行驶速度Vc服从均值为45km/h,标准差为10km/h的正态分布,即Vc~Nc(45,10);
D级路面下,行驶速度Vd服从均值为40km/h,标准差为15km/h的正态分布,即Vd~Nd(40,15);
E级路面下,行驶速度Ve服从均值为35km/h,标准差为15km/h的正态分布,即Ve~Ne(35,15);
F级路面下,行驶速度Vf服从均值为25km/h,标准差为20km/h的正态分布,即Vf~Nf(25,20)。
步骤S106:考虑分散性的履带车辆托带轮疲劳可靠度计算:在步骤S104和S105的基础上,通过MATLAB程序从工况敏感参数统计分布中随机抽取n个样本,代入不同工况类型下的多项式响应面模型中,得到不同工况类型下的单位里程损伤量。
根据线性疲劳损伤累积理论,在一定里程内,若托带轮焊缝结构损伤累积大于1,则认为结构发生失效。假定在任务总里程Y内,各路面级别对应的里程分别为Yckm/h、Ydkm/h、Yekm/h和Yfkm/h。根据行驶速度的正态分布进行M次抽样,在第i次抽样时,其托带轮焊缝结构累积损伤Qmi为
Qmi=Yci×Qci+Ydi×Qdi+Yei×Qei+Yfi×Qfi (5)
式中,Yci表示第i个样本在C级路面的里程,Qci表示第i个样本在C级路面的单位里程疲劳损伤量;Ydi表示第i个样本在D级路面的里程,Qdi表示第i个样本在D级路面的单位里程疲劳损伤量;Yei表示第i个样本在E级路面的里程,Qei表示第i个样本在E级路面的单位里程疲劳损伤量;Yfi表示第i个样本在F级路面的里程,Qfi表示第i个样本在F级路面的单位里程疲劳损伤量。
将计算获得的托带轮焊缝结构在任务剖面下的累积损伤量与Miner理论规定的损伤阈值进行比较,然后通过n次抽取和n次比较,统计未失效样本数量,确定该任务剖面下疲劳可靠度,如图4所示。
当Qmi>1时,说明此样本发生失效,累计失效数加1;当Qmi≤1时,说明此样本完好,未发生失效。直到所有样本抽取完毕,最后统计托带轮疲劳失效样本数量为m,从而得到托带轮疲劳任务可靠度R为
R=1-m/M
其中,M为根据行驶速度的正态分布进行抽样的次数。
在给定任务剖面中,即总里程Y是10000km时,抽样1000000次,共发生189790次失效,此时,托带轮疲劳任务可靠度是0.8102。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种履带车辆托带轮疲劳可靠度计算系统,其特征在于,包括第一确定单元、第二确定单元、第三确定单元、模型建立单元、统计分布确定单元和可靠度计算单元,其中:
给定总里程数,根据里程占比分布确定不同工况类型下的里程数;所述工况类型为路面等级;
所述第一确定单元用于检测及确定给定任务剖面涉及的不同工况类型及不同工况类型下的里程占比分布;
所述第二确定单元用于基于履带车辆托带轮结构疲劳机理分析及工况参数的显著性分析,确定履带车辆托带轮的薄弱单元及对应的工况敏感参数;
所述第三确定单元用于在基于工况分散性的前提下,利用履带车辆托带轮的多体动力学仿真模型,获取履带车辆托带轮载荷谱并确定履带车辆托带轮的薄弱单元单位里程疲劳损伤量;
所述模型建立单元用于基于所述履带车辆托带轮的薄弱单元单位里程疲劳损伤量与所述工况敏感参数确定二者之间的定量关系,建立所述不同工况类型下的响应面模型;所述响应面模型为变量为不同路面等级下的履带车辆托带轮单位里程疲劳损伤量与行驶速度的多项式响应面模型;所述不同路面等级包括4种,分别为C级路面、D级路面、E级路面和F级路面,其中,
Qc=1.991×10-7×Vc-5.12×10-6
式中,Qc表示C级路面下单位里程疲劳损伤量,Vc表示C级路面下履带车辆行驶速度;
Qd=3.256×10-8×Vd 2-2.641×10-6×Vd+5.675×10-5
式中,Qd表示D级路面下单位里程疲劳损伤量,Vd表示D级路面下履带车辆行驶速度;
Qe=-5.539×10-10×Ve 3+4.723×10-8×Ve 2-3.079×10-5
式中,Qe表示E级路面下单位里程疲劳损伤量,Ve表示E级路面下履带车辆行驶速度;
Qf=-4.727×10-9×Vf 3+6.993×10-7×Vf 2-3.313×10-5×Vf+0.0005293
式中,Qf表示F级路面下单位里程疲劳损伤量,Vf表示F级路面下履带车辆行驶速度,其中,行驶速度单位为km/h;
所述统计分布确定单元用于对所述工况敏感参数进行离散化处理,确定所述不同工况类型下的工况敏感参数对应的统计分布;
所述可靠度计算单元用于从所述工况敏感参数对应的统计分布中随机抽取n个样本,代入响应面模型中确定履带车辆托带轮的单位里程疲劳损伤量,基于所述单位里程疲劳损伤量计算任务剖面下的累积疲劳损伤量,所述任务剖面下的累积疲劳损伤量为Qmi,
Qmi=Yci×Qci+Ydi×Qdi+Yei×Qei+Yfi×Qfi
式中,Yci表示第i个样本在C级路面的里程,Qci表示第i个样本在C级路面的单位里程损伤量;Ydi表示第i个样本在D级路面的里程,Qdi表示第i个样本在D级路面的单位里程损伤量;Yei表示第i个样本在E级路面的里程,Qei表示第i个样本在E级路面的单位里程损伤量;Yfi表示第i个样本在F级路面的里程,Qfi表示第i个样本在F级路面的单位里程损伤量;将所述累积疲劳损伤量与损伤阈值进行比较,然后经过n次抽取和n次比较,统计未失效样本数量,计算所述任务剖面下的履带车辆托带轮疲劳可靠度;所述履带车辆托带轮疲劳可靠度为R,
R=1-m/M
其中,m为所述履带车辆托带轮疲劳失效样本数,M为根据所述工况敏感参数统计分布进行抽样的次数。
2.如权利要求1所述的履带车辆托带轮疲劳可靠度计算系统,其特征在于,所述工况敏感参数为在对所述薄弱单元的疲劳损伤影响最大的工况影响因素。
3.根据权利要求1所述的履带车辆托带轮疲劳可靠度计算系统,其特征在于,所述第三确定单元根据所述载荷谱进行雨流计数法,确定所述履带车辆托带轮的薄弱单元单位里程疲劳损伤量。
4.一种基于权利要求1-3中任一项所述的履带车辆托带轮疲劳可靠度计算系统的计算方法,其特征在于,该方法具体步骤如下:
步骤S101:检测及确定给定任务剖面涉及的不同工况类型及不同工况类型下的里程占比分布;
步骤S102:基于履带车辆托带轮结构疲劳机理分析及工况参数的显著性分析,确定履带车辆托带轮的薄弱单元及对应的工况敏感参数;
步骤S103:在基于工况分散性的前提下,利用履带车辆托带轮的多体动力学仿真模型,获取履带车辆托带轮载荷谱并确定履带车辆托带轮的薄弱单元单位里程疲劳损伤量;
步骤S104:基于所述履带车辆托带轮的薄弱单元单位里程疲劳损伤量与所述工况敏感参数确定二者之间的定量关系,建立不同工况类型下的响应面模型;
步骤S105:对所述工况敏感参数进行离散化处理,确定不同工况类型下工况敏感参数对应的统计分布;
步骤S106:通过MATLAB程序从所述工况敏感参数统计分布中随机抽取n个样本,代入响应面模型中确定履带车辆托带轮的单位里程疲劳损伤量,基于所述单位里程疲劳损伤量计算任务剖面下的累积疲劳损伤量,将所述累积疲劳损伤量与损伤阈值进行比较,然后经过n次抽取和n次比较,统计未失效样本数量,计算所述任务剖面下的履带车辆托带轮疲劳可靠度。
5.一种履带车辆,其特征在于,所述履带车辆包括托带轮及如权利要求1-3中任一项所述的履带车辆托带轮疲劳可靠度计算系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110673262.4A CN113569332B (zh) | 2021-06-17 | 2021-06-17 | 托带轮疲劳可靠度计算系统、方法、履带车辆 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110673262.4A CN113569332B (zh) | 2021-06-17 | 2021-06-17 | 托带轮疲劳可靠度计算系统、方法、履带车辆 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113569332A CN113569332A (zh) | 2021-10-29 |
CN113569332B true CN113569332B (zh) | 2024-03-19 |
Family
ID=78162225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110673262.4A Active CN113569332B (zh) | 2021-06-17 | 2021-06-17 | 托带轮疲劳可靠度计算系统、方法、履带车辆 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113569332B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115311759B (zh) * | 2022-07-08 | 2023-09-05 | 东风汽车集团股份有限公司 | 一种车辆耐久目标获取方法、装置、设备及存储介质 |
CN116148104B (zh) * | 2023-04-18 | 2023-08-01 | 中汽研汽车检验中心(天津)有限公司 | 基于实际工况的车轮双轴疲劳试验载荷谱确定方法及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4583791A (en) * | 1984-03-13 | 1986-04-22 | Kabushiki Kaisha Komatsu Seisakusho | Crawler vehicle undercarriage |
CN106768997A (zh) * | 2015-11-19 | 2017-05-31 | 黑龙江恒能自控科技有限公司 | 一种履带式装甲底盘侧减速器齿轮接触疲劳分析方法 |
JP2019096041A (ja) * | 2017-11-22 | 2019-06-20 | 横浜ゴム株式会社 | シミュレーション方法、その装置およびプログラム |
CN112069599A (zh) * | 2020-08-28 | 2020-12-11 | 北京科技大学 | 一种履带车辆行驶性能可靠性预测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE60300521T2 (de) * | 2002-06-13 | 2006-02-23 | Société de Technologie Michelin | Stabilitätsregelungssystem eines Fahrzeuges unter Verwendung einer für alle Reifen charakteristischen Invarianzfunktion |
US7469200B2 (en) * | 2003-05-14 | 2008-12-23 | Ford Global Technologies, Llc | Method and apparatus for predicting belt separation failure in aging tires by computer simulation |
-
2021
- 2021-06-17 CN CN202110673262.4A patent/CN113569332B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4583791A (en) * | 1984-03-13 | 1986-04-22 | Kabushiki Kaisha Komatsu Seisakusho | Crawler vehicle undercarriage |
CN106768997A (zh) * | 2015-11-19 | 2017-05-31 | 黑龙江恒能自控科技有限公司 | 一种履带式装甲底盘侧减速器齿轮接触疲劳分析方法 |
JP2019096041A (ja) * | 2017-11-22 | 2019-06-20 | 横浜ゴム株式会社 | シミュレーション方法、その装置およびプログラム |
CN112069599A (zh) * | 2020-08-28 | 2020-12-11 | 北京科技大学 | 一种履带车辆行驶性能可靠性预测方法 |
Non-Patent Citations (3)
Title |
---|
Rigid-Flexible Coupling Simulation and Test Research of Thin-Skinned Body of Tracked Vehicle;Hongchang Xu, Yang Zhang,Ming Ma,Chunhui Zhao,Xiaoce Sun;2020 IEEE 3rd International Conference on Information Systems and Computer Aided Education (ICISCAE);114-117 * |
坦克装甲车辆对试验测试技术的需求;毛明,李振平,邓高寿,范凤明;计算机测量与控制;第24卷(第3期);8-10 * |
履带式装甲车辆车体疲劳寿命预测仿真研究;尚其刚;王红岩;杨涛;马伟标;;计算机仿真(第06期);19-23 * |
Also Published As
Publication number | Publication date |
---|---|
CN113569332A (zh) | 2021-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113569332B (zh) | 托带轮疲劳可靠度计算系统、方法、履带车辆 | |
CN106441530B (zh) | 一种基于长标距光纤光栅传感技术的桥梁动态称重方法及动态称重系统 | |
US20220415166A1 (en) | Power distribution and vehicle self-learning-based truck overload identification method | |
Bevan et al. | Development and validation of a wheel wear and rolling contact fatigue damage model | |
CN104990709B (zh) | 用于检测机车轴承故障的方法 | |
CN110222437B (zh) | 列车车列健康状态的评估方法、装置及存储介质 | |
CN105241660A (zh) | 基于健康监测数据的高铁大型桥梁性能评定方法 | |
CN108664715B (zh) | 一种路面积水车辙三维评价与行车安全分析方法 | |
CN110222436B (zh) | 列车零部件健康状态的评估方法、装置及存储介质 | |
CN108515984B (zh) | 一种车轮伤损检测方法及装置 | |
CN104309435B (zh) | 一种路面不平度在线辨识方法 | |
CN110210161B (zh) | 列车车辆健康状态的评估方法、装置及存储介质 | |
CN105923014B (zh) | 一种基于证据推理规则的轨道高低不平顺幅值估计方法 | |
CN105067209A (zh) | 基于桥梁健康监测变形数据判定桥梁结构刚度变化的方法 | |
CN103543020B (zh) | 一种基于适时技术检测数据评估二手汽车成新率的方法 | |
CN106018557A (zh) | 一种基于小波分析的轨道表面凹陷检测方法 | |
CN102890143B (zh) | 一种融合车速信息和前后轴箱加速度信息的轨道局部缺陷车载检测方法 | |
CN111942397A (zh) | 危险驾驶行为监控方法、装置、存储介质 | |
CN110203249A (zh) | 列车修程的处理方法、装置及存储介质 | |
CN113158141B (zh) | 基于大数据的新能源汽车超载检测方法 | |
CN104850678A (zh) | 基于行车走行性的公路桥梁伸缩装置走行服役性能评定方法 | |
CN104991986A (zh) | 公路桥梁支座及伸缩装置的竖向抗冲击服役性能评定方法 | |
CN114179858B (zh) | 基于车轮健康状态的车轮镟修方法及装置 | |
CN114275003B (zh) | 在建超长tbm隧道临时运输轨道不平顺性检测评价方法 | |
CN115081784A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |