, 胡全全1, 张义民2, 郭丙帅11. 东北大学 机械工程与自动化学院, 辽宁 沈阳 110819;
2. 沈阳化工大学 机械工程学院，辽宁 沈阳 110142
收稿日期：2018-01-02
基金项目：国家自然科学联合基金资助项目(U1710119);辽宁省联合基金资助项目(U1708254)。
作者简介：杨周(1979-), 女, 辽宁沈阳人, 东北大学副教授;
张义民(1958-), 男, 吉林长春人, 沈阳化工大学教授, 博士生导师。
摘要：为研究车辆零件在磨损作用影响下的可靠性问题, 将零件结构尺寸最大磨损量考虑到渐变可靠性模型中, 得到随时间变化的可靠性状态方程的数学模型, 并利用二阶矩和随机摄动法得出动态可靠性指标, 进而获得零部件的渐变可靠性及灵敏度的定量分析结果.通过算例与蒙特卡洛相结合验证了该方法的有效性, 拟合得到了零部件磨损渐变可靠性分布规律曲线和灵敏度曲线, 很好地反映了零部件几何尺寸在磨损状态下危险截面的可靠性变化情况, 为车辆零件的磨损渐变可靠性研究提供了理论参考.
关键词：车辆零件磨损渐变可靠性蒙特卡洛灵敏度
Reliability Analysis of Wear Gradient of Automobile Parts
YANG Zhou1
, HU Quan-quan1, ZHANG Yi-min2, GUO Bing-shuai11. School of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819, China;
2. College of Mechanical Engineering, Shenyang University of Chemical Technology, Shenyang 110142, China
Corresponding author: YANG Zhou, E-mail: yangzhou@mail.neu.edu.cn
Abstract: To study the reliability of automobile parts under the influence of abrasion, the biggest wear loss is taken into account in the gradual changing reliability model and the model of reliability state with the changing of time is obtained. Then, the secondary moment and the random perturbation method are adopted to give a quantitative analysis of the gradual changing reliability and the sensitivity. Finally, the effectiveness of the proposed method is confirmed by an example and the Monte Carlo. The fitting curves of the regularities of distribution of reliability and of sensitivity with abrasion are obtained, which can well reflect the change of reliability of dangerous sections of automobile parts in wear conditions and provide references for the gradual changing reliability with abrasion.
Key words: automobile partsabrasiongradual changing reliabilityMonte Carlosensitivity
随着人们对车辆质量更高层次的追求, 可靠性理论方法的研究也需要更深入的探讨.汽车可靠性主要是指在预期的使用时间内和零部件所需的工作条件下汽车零部件能够完成规定的功能的能力[1].由于车辆零部件在使用过程中不断磨损[2], 零件磨损量随着使用时间而逐渐增大, 一个零部件失效通常是由于磨损导致的.磨损的存在会使零部件几何尺寸(体积)变小, 进而导致零部件之间不能够很好地吻合, 失去原有人们期望的功能而失效.车辆零件剧烈磨损的磨损量与工作时间之间具有一定的渐变性[3-4].零部件失效以后就不能完成预定的功能, 如果继续使用甚至会产生一定的危险.因此, 零部件的可靠性是时间的函数, 并呈现出显著的渐变性[5].
本文采用随机摄动方法、二阶矩技术、可靠性理论等, 讨论了车辆零件在考虑磨损时的动态可靠性.在知道零部件尺寸、初始强度等基本随机变量的概率特性下, 能够准确地得出随时间变化的车辆磨损零件动态可靠性.这对以后车辆零件磨损的渐变可靠性的研究具有重要意义.
1 随时间变化的可靠性分析方法在静态可靠性中, 由应力-强度干涉理论可知, 其零部件临界失效状态方程为
 | (1) |
可将零部件的状态函数表示为
 | (2) |
 | (3) |
 | (4) |
但在实际状况中由于各种因素的影响使得零部件的强度与载荷不是固定不变的量, 而都是关于时间的函数.失效状态方程可表示为
 | (5) |
 | (6) |
 | (7) |
 | (8) |
 | (9) |
 | (10) |
 | (11) |
 | (12) |
 | (13) |
 | (14) |
 | (15) |
 | (16) |
2 车辆零件磨损可靠性设计实例在车辆零件可靠性设计中, 零件的磨损是影响车辆可靠性的最大因素[11].在实际使用状况中, 磨损速度主要是由汽车行驶状况、零部件所处的环境温度、润滑特性、载荷等因素所决定的.磨损速度具有分散性, 属于随机变量, 在研究零部件磨损时通常是研究其最大磨损量.若某时刻磨损最大速度vmax已知, 则经过t时间后最大磨损量为
 | (17) |
人们为了实现某一特定的预期功能而设计出某一零件, 但由于理论、材料、工艺等方面的因素, 使得机械零件在应用中会产生各式各样的不可靠的现象.由此, 必须大力进行总结开拓现代的机械可靠性设计理论与方法.现以车辆拉杆和扭杆为例, 研究车辆零部件基于磨损的动态可靠性及相关变量的灵敏度问题.
2.1 考虑磨损失效的车辆拉杆渐变可靠性设计模型图 1是截面为管形的某车辆拉杆结构, 且承受拉压载荷作用.
图 1(Fig. 1)
 | 图 1 拉杆结构Fig.1 Structure of tension rod |
管形直拉杆的拉应力为
 | (18) |
应用应力极限状态表示出车辆拉杆的状态方程为
 | (19) |
2.2 考虑磨损失效的车辆扭杆渐变可靠性设计模型图 2是扭杆弹簧结构, 连接于车架与车轮之间的一种弹性元件.在汽车行驶中当遇到障碍物汽车向上跳动时, 车轮就会通过摆臂使扭杆发生弹性变形, 并存储能力, 使得车轮在脱离障碍物时能够迅速回位.按其截面可将扭杆分为:片状扭杆、圆形扭杆、管形扭杆及组合式扭杆等.其中圆形扭杆最为常见, 而因管形扭杆能够节省材料降低成本与质量, 应用也较为广泛.
图 2(Fig. 2)
 | 图 2 扭杆弹簧Fig.2 Torsion bar spring (a)—扭杆结构；(b)—圆形扭杆；(c)—管形扭杆. |
通过图 2对圆形和管形扭杆进行扭力分析, 其扭杆所受应力为
 | (20) |
应用应力极限状态表示出扭杆的状态方程为
 | (21) |
3 数值算例3.1 拉杆在一般实际情况中, 机械零部件几何尺寸变量和所使用材料的强度变量都是服从正态分布的随机变量.某车辆拉杆的基本随机参数如表 1所示, 且变量间相互独立.同时考虑车辆拉杆的磨损速度为0.001 156 8 mm/h.
表 1(Table 1)
 表 1 某型车辆拉杆危险截面基本随机参数均值和标准差值Table 1 Mean and standard deviation of the basic random parameters of dangerous cross-section of a tension rod
| 表 1 某型车辆拉杆危险截面基本随机参数均值和标准差值 Table 1 Mean and standard deviation of the basic random parameters of dangerous cross-section of a tension rod |
通过运用上述动态可靠性方法得出随时间变化的拉杆可靠性和灵敏度的变化曲线, 并分别通过Monte Carlo(虚线)进行验证.
图 3表示的是运用时变磨损可靠性模型(实线)计算得到的车辆拉杆动态可靠性曲线与运用Monte Carlo(虚线)仿真实验计算得到的车辆拉杆动态可靠性曲线的拟合结果.从图中可以看出两者的变化趋势相一致, 吻合较好.从而说明运用时变磨损可靠性模型计算得到的可靠性结果是正确的.
图 3(Fig. 3)
 | 图 3 车辆拉杆可靠性变化曲线Fig.3 Reliability changing curve of vehicle tension rod |
图 4和图 5分别表示的是均值灵敏度?μxi曲线和标准差灵敏度?σxi曲线与Monte Carlo(虚线)仿真实验计算得到的灵敏度曲线的拟合结果.从图中可以看出两者的变化趋势相一致, 吻合较好.这说明运用时变磨损可靠性模型计算得到的车辆拉杆动态灵敏度结果是正确的.
图 4(Fig. 4)
 | 图 4 车辆拉杆均值灵敏度变化曲线Fig.4 Mean sensitivity changing curve of vehicle tension rod |
图 5(Fig. 5)
 | 图 5 车辆拉杆标准差灵敏度变化曲线Fig.5 Standard deviation sensitivity changing curve of vehicle tension rod |
3.2 扭杆以某汽车的管型扭杆为例, 其基本随机参数如表 2所示, 扭杆所受的扭矩和材料的疲劳极限均为服从正态分布的随机变量.扭杆的磨损速度为0.006 035 24 mm/h.
表 2(Table 2)
表 2 某型汽车的扭杆危险截面基本随机参数均值和标准差Table 2 Mean and standard deviation of the basic random parameters of dangerous cross-section of a torsion bar
| 表 2 某型汽车的扭杆危险截面基本随机参数均值和标准差 Table 2 Mean and standard deviation of the basic random parameters of dangerous cross-section of a torsion bar |
通过运用上述动态可靠性方法得到此扭杆随时间变化的动态可靠性曲线和灵敏度曲线, 并分别通过Monte Carlo(虚线)进行验证.
图 6中表示的是运用时变磨损可靠性模型(实线)计算得到的车辆扭杆动态可靠性曲线与运用Monte Carlo(虚线)仿真实验计算得到的车辆扭杆动态可靠性曲线的拟合结果, 两者的变化趋势相一致, 吻合较好.从而说明运用时变磨损可靠性模型计算得到的汽车扭杆动态可靠性结果是正确的.
图 6(Fig. 6)
 | 图 6 扭杆可靠性变化曲线Fig.6 Reliability changing curve of vehicle torsion bar |
图 7和图 8分别表示的是运用时变磨损可靠性模型计算得到的扭杆均值灵敏度?μxi曲线和标准差灵敏度?σxi曲线与Monte Carlo(虚线)仿真实验计算得到的灵敏度曲线的拟合结果, 从图中可以看出两者的变化趋势相一致, 吻合较好.从而说明运用时变磨损可靠性模型计算得到的车辆扭杆动态灵敏度结果是正确的.
图 7(Fig. 7)
 | 图 7 车辆扭杆均值灵敏度变化曲线Fig.7 Mean sensitivity variation curve of vehicle torsion bar |
图 8(Fig. 8)
 | 图 8 车辆扭杆标准差灵敏度变化曲线Fig.8 Standard deviation sensitivity changing curve of vehicle torsion bar |
4 结论1) 本文在车辆零件可靠性研究的基础上, 将影响车辆零件可靠性的磨损因素考虑到渐变可靠性数学模型中, 使之进一步符合实际情况.
2) 在考虑零件的磨损基础上建立了渐变功能数学模型, 结合随机摄动法和二阶矩得出相应的动态可靠性指标和随时间变化的渐变可靠性曲线和渐变灵敏度曲线, 并用Monte Carlo方法验证其方法的可行性.
3) 在磨损的条件下, 几何尺寸对灵敏度的影响最大, 载荷最小.即几何尺寸对可靠性影响最大, 这与静态情况下一致.
参考文献
| [1] | 何在庚. 机械动态与渐变可靠性理论与技术评述[J].科技展望, 2016(6): 37. ( He Zai-geng. Dynamic mechanical reliability theory and technique with gradients review[J].Science and Technology, 2016(6): 37.) |
| [2] | 张义民. 机械动态与渐变可靠性理论与技术评述[J].机械工程学报, 2013, 49(20): 101–114. ( Zhang Yi-min. Dynamic mechanical reliability theory and technique with gradients review[J].Journal of Mechanical Engineering, 2013, 49(20): 101–114.) |
| [3] | Zhang Y M, Liu Q L. Practical reliability-based analysis of coil tube-spring[J].Journal of Mechanical Engineering Science, 2002, 216(C2): 179–182. |
| [4] | Zhang Y M, Liu Q L. Reliability-based design of automobile components[J].Journal of Automobile Engineering, 2002, 216(D6): 455–471. |
| [5] | Yang Z, Zhang Y M, Zhang X F, et al. Reliability sensitivity-based correlation coefficient calculation in structural reliability analysis[J].Chinese Journal of Mechanical Engineering, 2012, 25(3): 608–614.DOI:10.3901/CJME.2012.03.608 |
| [6] | Jun K J, Park T W, Lee S H, et al. Prediction of fatigue life and estimation of its reliability on the parts of an air suspension system[J].International Journal Automotion Technology, 2008, 9(6): 741–747.DOI:10.1007/s12239-008-0088-4 |
| [7] | Lee S B, Baik S, Yim H J. Optimal reliability design for thin-walled beam of vehicle structure considering vibration[J].International Journal Automotion Technology, 2003, 4(3): 135–140. |
| [8] | Melchers R E, Ahammed M. A fast approximate method for parameter sensitivity estimation in Monte Carlo structural reliability[J].Computers & Structures, 2004, 82(1): 55–61. |
| [9] | Zhang Y M, He X D, Liu Q L, et al. Reliability sensitivity of automobile components with arbitrary distribution parameters[J].Journal of Automobile Engineering, 2005, 219(2): 165–182.DOI:10.1243/095440705X5894 |
| [10] | Zhang Y M, Liu Q L, Wen B C. Reliability sensitivity of multi-degree-of-freedom uncertain nonlinear systems with independence failure modes[J].Journal of Mechanical Science and Technology, 2007, 21(6): 908–912.DOI:10.1007/BF03027067 |
| [11] | 陆克久.车用发动机金属零件磨损状态的判别研究[C]//2014中国汽车工程学会年会论文集.上海, 2014: 154-156. ( Lu Ke-jiu.Research on metal part wear states discrimination of vehicle engine[C]//2014 China Automotive Engineering Society Annual Meeting Proceedings.Shanghai, 2014: 154-156.) |
