2. 沈阳化工大学 机械工程学院,辽宁 沈阳 110142
2. College of Mechanical Engineering, Shenyang University of Chemical Technology, Shenyang 110142, China
随着人们对车辆质量更高层次的追求, 可靠性理论方法的研究也需要更深入的探讨.汽车可靠性主要是指在预期的使用时间内和零部件所需的工作条件下汽车零部件能够完成规定的功能的能力[1].由于车辆零部件在使用过程中不断磨损[2], 零件磨损量随着使用时间而逐渐增大, 一个零部件失效通常是由于磨损导致的.磨损的存在会使零部件几何尺寸(体积)变小, 进而导致零部件之间不能够很好地吻合, 失去原有人们期望的功能而失效.车辆零件剧烈磨损的磨损量与工作时间之间具有一定的渐变性[3-4].零部件失效以后就不能完成预定的功能, 如果继续使用甚至会产生一定的危险.因此, 零部件的可靠性是时间的函数, 并呈现出显著的渐变性[5].
本文采用随机摄动方法、二阶矩技术、可靠性理论等, 讨论了车辆零件在考虑磨损时的动态可靠性.在知道零部件尺寸、初始强度等基本随机变量的概率特性下, 能够准确地得出随时间变化的车辆磨损零件动态可靠性.这对以后车辆零件磨损的渐变可靠性的研究具有重要意义.
1 随时间变化的可靠性分析方法在静态可靠性中, 由应力-强度干涉理论可知, 其零部件临界失效状态方程为
(1) |
式中:r为材料强度; δ为附加载荷; Y为参数向量, Y=(Y1, Y2, …, Yn)T.
可将零部件的状态函数表示为
(2) |
根据可靠性设计的随机摄动法, 可求出静态可靠性指标β[6-8]:
(3) |
若Y中所有随机参数都是正态分布时, 可得静态一阶估计量的可靠性,
(4) |
但在实际状况中由于各种因素的影响使得零部件的强度与载荷不是固定不变的量, 而都是关于时间的函数.失效状态方程可表示为
(5) |
g(Y, t)为状态函数, 可表示车辆零件的两种状态:
(6) |
通过式(3), 式(4)可知, 零部件的时变可靠性指标与时变可靠性计算式分别为
(7) |
(8) |
时变失效率表达式为
(9) |
根据应力-强度干涉理论, 可靠性对初始参数均值μx和标准差σx的灵敏度为
(10) |
(11) |
式中:
(12) |
(13) |
(14) |
(15) |
(16) |
再代入式(10)和式(11), 就计算出基础变量的可靠性灵敏度∂R/∂XT和∂R/∂Var(X).
2 车辆零件磨损可靠性设计实例在车辆零件可靠性设计中, 零件的磨损是影响车辆可靠性的最大因素[11].在实际使用状况中, 磨损速度主要是由汽车行驶状况、零部件所处的环境温度、润滑特性、载荷等因素所决定的.磨损速度具有分散性, 属于随机变量, 在研究零部件磨损时通常是研究其最大磨损量.若某时刻磨损最大速度vmax已知, 则经过t时间后最大磨损量为
(17) |
当其最大磨损量超过许用磨损量Umax(mm)时零部件就可能无法完成人们预先设计的功能, 即产生失效现象.因此磨损是零部件可靠性的影响因素.
人们为了实现某一特定的预期功能而设计出某一零件, 但由于理论、材料、工艺等方面的因素, 使得机械零件在应用中会产生各式各样的不可靠的现象.由此, 必须大力进行总结开拓现代的机械可靠性设计理论与方法.现以车辆拉杆和扭杆为例, 研究车辆零部件基于磨损的动态可靠性及相关变量的灵敏度问题.
2.1 考虑磨损失效的车辆拉杆渐变可靠性设计模型图 1是截面为管形的某车辆拉杆结构, 且承受拉压载荷作用.
管形直拉杆的拉应力为
(18) |
式中: Q为载荷; d为管形截面的内径; D为管形截面的外径.
应用应力极限状态表示出车辆拉杆的状态方程为
(19) |
式中:r为零件强度值;σmax为变量Q, D和d的函数.
2.2 考虑磨损失效的车辆扭杆渐变可靠性设计模型图 2是扭杆弹簧结构, 连接于车架与车轮之间的一种弹性元件.在汽车行驶中当遇到障碍物汽车向上跳动时, 车轮就会通过摆臂使扭杆发生弹性变形, 并存储能力, 使得车轮在脱离障碍物时能够迅速回位.按其截面可将扭杆分为:片状扭杆、圆形扭杆、管形扭杆及组合式扭杆等.其中圆形扭杆最为常见, 而因管形扭杆能够节省材料降低成本与质量, 应用也较为广泛.
通过图 2对圆形和管形扭杆进行扭力分析, 其扭杆所受应力为
(20) |
式中:T为附加给扭杆的扭矩; D, d分别为管形扭杆横截面的外径和内径(圆形扭杆内径为零).
应用应力极限状态表示出扭杆的状态方程为
(21) |
式中:r为扭杆的材料强度; τmax为变量T, D和d的函数.
3 数值算例 3.1 拉杆在一般实际情况中, 机械零部件几何尺寸变量和所使用材料的强度变量都是服从正态分布的随机变量.某车辆拉杆的基本随机参数如表 1所示, 且变量间相互独立.同时考虑车辆拉杆的磨损速度为0.001 156 8 mm/h.
通过运用上述动态可靠性方法得出随时间变化的拉杆可靠性和灵敏度的变化曲线, 并分别通过Monte Carlo(虚线)进行验证.
图 3表示的是运用时变磨损可靠性模型(实线)计算得到的车辆拉杆动态可靠性曲线与运用Monte Carlo(虚线)仿真实验计算得到的车辆拉杆动态可靠性曲线的拟合结果.从图中可以看出两者的变化趋势相一致, 吻合较好.从而说明运用时变磨损可靠性模型计算得到的可靠性结果是正确的.
图 4和图 5分别表示的是均值灵敏度∂μxi曲线和标准差灵敏度∂σxi曲线与Monte Carlo(虚线)仿真实验计算得到的灵敏度曲线的拟合结果.从图中可以看出两者的变化趋势相一致, 吻合较好.这说明运用时变磨损可靠性模型计算得到的车辆拉杆动态灵敏度结果是正确的.
以某汽车的管型扭杆为例, 其基本随机参数如表 2所示, 扭杆所受的扭矩和材料的疲劳极限均为服从正态分布的随机变量.扭杆的磨损速度为0.006 035 24 mm/h.
通过运用上述动态可靠性方法得到此扭杆随时间变化的动态可靠性曲线和灵敏度曲线, 并分别通过Monte Carlo(虚线)进行验证.
图 6中表示的是运用时变磨损可靠性模型(实线)计算得到的车辆扭杆动态可靠性曲线与运用Monte Carlo(虚线)仿真实验计算得到的车辆扭杆动态可靠性曲线的拟合结果, 两者的变化趋势相一致, 吻合较好.从而说明运用时变磨损可靠性模型计算得到的汽车扭杆动态可靠性结果是正确的.
图 7和图 8分别表示的是运用时变磨损可靠性模型计算得到的扭杆均值灵敏度∂μxi曲线和标准差灵敏度∂σxi曲线与Monte Carlo(虚线)仿真实验计算得到的灵敏度曲线的拟合结果, 从图中可以看出两者的变化趋势相一致, 吻合较好.从而说明运用时变磨损可靠性模型计算得到的车辆扭杆动态灵敏度结果是正确的.
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. ) |