《万方数据-数字化期刊群》全文上网期刊
CNKI《中国学术期刊(网络版)》全文收录期刊
《中文科技期刊数据库》(维普网)全文收录期刊
超星期刊域出版平台、博看网全文收录期刊
日本科学技术振兴机构数据库收录

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于性能退化的失效阈值变化速率对产品可靠性的影响

王新刚 马瑞敏 黄小凯 常苗鑫

王新刚, 马瑞敏, 黄小凯, 等. 基于性能退化的失效阈值变化速率对产品可靠性的影响[J]. 航天器环境工程, 2020, 37(2): 137-142 doi:  10.12126/see.2020.02.006
引用本文: 王新刚, 马瑞敏, 黄小凯, 等. 基于性能退化的失效阈值变化速率对产品可靠性的影响[J]. 航天器环境工程, 2020, 37(2): 137-142 doi:  10.12126/see.2020.02.006
WANG X G, MA R M, HUANG X K, et al. Effect of failure threshold change rate on the anlysis of product reliability based on performance degradation[J]. Spacecraft Environment Engineering, 2020, 37(2): 137-142 doi:  10.12126/see.2020.02.006
Citation: WANG X G, MA R M, HUANG X K, et al. Effect of failure threshold change rate on the anlysis of product reliability based on performance degradation[J]. Spacecraft Environment Engineering, 2020, 37(2): 137-142 doi:  10.12126/see.2020.02.006

基于性能退化的失效阈值变化速率对产品可靠性的影响

doi: 10.12126/see.2020.02.006
基金项目: 国家自然科学基金项目(编号:51475086);中央高校基本科研业务项目(编号:N162312001);北京卫星环境工程研究所创新基金项目(编号:CAST-BISEE 2019-019)
详细信息
    作者简介:

    王新刚(1979—),男,博士学位,研究方向为机械零部件的动态可靠性设计。E-mail: wangxingang1217@126.com

  • 中图分类号: TB114.3

Effect of failure threshold change rate on the anlysis of product reliability based on performance degradation

  • 摘要: 性能退化分析是一种有效弥补数据不足的可靠性评估方法。采用线性变化过程描述机械产品的性能退化过程,考虑产品受极值冲击和运行冲击的情况,分析性能退化过程对突发失效的阈值影响。基于变失效阈值的性能退化过程和变失效阈值的冲击过程,建立突发失效和性能退化失效同时存在下的竞争失效可靠性模型。结合实例对可靠性模型进行分析,结果表明:若采用定失效阈值,忽略失效阈值的变化,产品的可靠度会被高估;在一定范围内,机械产品的可靠度随着失效阈值变化速率的增大而下降,随着性能退化量的变大而提升。
  • 图  1  性能退化失效过程

    Figure  1.  Failure process of performance degradation

    图  2  突发失效过程

    Figure  2.  Suddenly-occurred failure process

    图  3  运行冲击模型

    Figure  3.  Operation shock model

    图  4  极值冲击情况下不同失效阈值变化速率下的产品可靠度

    Figure  4.  Product reliability at different failure threshold degradation rates in extreme operation shock

    图  5  极值冲击情况下不同L值下的产品可靠度

    Figure  5.  Product reliability for different values of L in extreme operation shock

    图  6  运行冲击情况下不同失效阈值变化速率下的产品可靠度

    Figure  6.  Product reliability at different failure threshold degradation rates in operation shock

    图  7  运行冲击情况下不同L值下的产品可靠度

    Figure  7.  Product reliability for different values of L in operation shock

    表  1  微引擎可靠性分析参数

    Table  1.   Parameter values for micro-engine reliability analysis

    参数数据来源
    D0/GPa1.2假设
    D1/GPa1.4文献[15]
    D2/GPa1.3假设
    H0/μm30.001 25文献[15]
    L/μm37×10-4假设
    β0/μm3${\mu _{\beta 0}}$$ = 1.2 \times {10^{{\text{-9}}}}$
    ${\sigma _{\beta 0}}$$ = 1 \times {10^{{\text{-9}}}}$
    假设
    β1/μm3${\mu _{\beta 1}}$$ = 8.48 \times {10^{{\text{-9}}}}$
    ${\sigma _{\beta 1}}$$ = 6 \times {10^{{\text{-10}}}}$
    文献[15]
    β2/μm3${\mu _{\beta 2}}$$ = 1 \times {10^{{\text{-5}}}}$
    ${\sigma _{\beta 2}}$$ = 1 \times {10^{{\text{-6}}}}$
    假设
    β3/μm3${\mu _{\beta 3}}$$ = 1 \times {10^{{\text{-5}}}}$
    ${\sigma _{\beta 3}}$$ = 1 \times {10^{{\text{-6}}}}$
    假设
    λ/转5×10-5假设
    Yi/μm3${\mu _Y}$$ = 1 \times {10^{{\text{-4}}}}$
    ${\sigma _Y}$$ = 2 \times {10^{{\text{-5}}}}$
    文献[15]
    Wi/μm3${\mu _W} = 1.2$
    ${\sigma _W} = 0.2$
    假设
    P0.131 0假设
    Q0.869 0假设
    下载: 导出CSV
  • [1] 尹继东. 基于性能退化数据的产品竞争失效可靠性评估方法及理论研究[D]. 成都: 西南交通大学, 2018: 13-17
    [2] 邓爱民, 陈循, 张春华, 等. 基于性能退化数据的可靠性评估[J]. 宇航学报, 2006, 27(3): 228-234

    DENG A M, CHEN X, ZHANG C H, et al. Reliability evaluation based on performance degradation data[J]. Journal of Astronautics, 2006, 27(3): 228-234
    [3] CRK V. Reliability assessment from degradation data[C]//Reliability & Maintainability Symposium, IEEE Proceedings. Los Angeles, CA, 2000: 155-161
    [4] PARK C, PADGETT W J. Stochastic degradation models with several accelerating variables[J]. IEEE Transactions on Reliability, 2006, 55(2): 379-390 doi:  10.1109/TR.2006.874937
    [5] 袁容. 基于性能退化分析的可靠性方法研究[D]. 成都: 电子科技大学, 2015: 77-83
    [6] LI W, PHAM H. An inspection-maintenance model for systems with multiple competing processes[J]. IEEE Transactions on Reliability, 2005, 54(2): 318-327 doi:  10.1109/TR.2005.847264
    [7] WANG Z L, HUANG H Z, DU L, et al. Reliability analysis on competitive failure processes under fuzzy degradation data[J]. Applied Soft Computing, 2011, 11(3): 2964-2973 doi:  10.1016/j.asoc.2010.11.018
    [8] 刘晓娟, 王华伟, 徐璇. 考虑多退化失效和突发失效之间竞争失效的可靠性评估方法[J]. 中国机械工程, 2017, 28(1): 7-12 doi:  10.3969/j.issn.1004-132X.2017.01.002

    LIU X J, WANG H W, XU X. Reliability assessment based on competition failure considering multi-degradation and catastrophic failure[J]. China Mechanical Engineering, 2017, 28(1): 7-12 doi:  10.3969/j.issn.1004-132X.2017.01.002
    [9] AN Z W, SUN D M. Reliability modeling for systems subject to multiple dependent competing failure processes with shock loads above a certain level[J]. Reliability Engineering & System Safety, 2017, 157: 129-138
    [10] RAFIEE K, FENG Q, COIT D W. Reliability modeling for dependent competing failure processes with changing degradation rate[J]. LIE Transactions, 2014, 46(5): 483-496
    [11] 黄文平, 周经伦, 宁菊红, 等. 基于变失效阈值的竞争失效可靠性模型[J]. 系统工程与电子技术, 2017, 39(4): 941-946 doi:  10.3969/j.issn.1001-506X.2017.04.34

    HUANG W P, ZHOU J L, NING J H, et al. Reliability model for competing failure with shift thresholds[J]. Journal of Systems Engineering and Electronics, 2017, 39(4): 941-946 doi:  10.3969/j.issn.1001-506X.2017.04.34
    [12] JIANG L, FENG Q, COIT D W. Reliability and maintenance modeling for dependent competing failure processes with shifting failure thresholds[J]. IEEE Transactions on Reliability, 2012, 61(4): 932-948 doi:  10.1109/TR.2012.2221016
    [13] RAFIEE K, FENG Q, COIT D W. Reliability assessment of competing risks with generalized mixed shock models[J]. Reliability Engineering & System Safety, 2017, 159: 1-11
    [14] RAFIEE K, FENG Q, COIT D W. Reliability modeling for dependent competing failure processes with changing degradation rate[J]. IEEE Transactions on Reliability, 2014, 46(5): 483-49
    [15] TANNER D M, DUGGER M T. Wear mechanisms in a reliability methodology[J]. Proceedings of the SPIE, 2003, 4980: 22-40 doi:  10.1117/12.476345
  • [1] 种婧宜, 周昊澄, 王敏, 陈余军.  通信卫星系统潜在故障识别与应对方法 . 航天器环境工程, 2023, 40(1): 84-91. doi: 10.12126/see.2022091
    [2] 付强, 任洪鹏, 霍琳颖, 李彦恒, 朱光, 吴昊, 王然江.  新型空天货物转载支架车静强度及可靠性仿真分析 . 航天器环境工程, 2022, 39(4): 368-373. doi: 10.12126/see.2022.04.006
    [3] 黄首清, 代巍, 姚泽民, 张禹森, 黄小凯, 刘守文.  两种工程化的航天器用滚动轴承加速寿命试验方法 . 航天器环境工程, 2021, 38(4): 413-419. doi: 10.12126/see.2021.04.006
    [4] 肖波, 吕倩倩, 谷瀚天, 李培蕾, 张俊峰.  宇航用100 V厚膜DC/DC电源变换器极限评估试验研究 . 航天器环境工程, 2021, 38(2): 176-182. doi: 10.12126/see.2021.02.010
    [5] 朱爱红, 孙亮, 李明, 姜德胜, 温慧博.  地空导弹武器装备环境适应性试验 现状与未来发展 . 航天器环境工程, 2021, 38(5): 604-608. doi: 10.12126/see.2021.05.017
    [6] 刘闯斌, 许晓冬, 李昊, 庞智敏, 张金辉, 魏强.  长期存储航天产品评估的系统工程方法应用 . 航天器环境工程, 2020, 37(5): 525-530. doi: 10.12126/see.2020.05.017
    [7] 刘岐, 沈鸣杰, 董艺.  电离总剂量效应对浮栅型Flash存储器擦写耐久与数据保持特性的影响研究 . 航天器环境工程, 2019, 36(3): 264-270. doi: 10.12126/see.2019.03.011
    [8] 张义, 张红旗, 王征, 陶兆增, 范壮壮, 郭燕君, 焦美荣, 于福莉, 夏占军.  基于宇航元器件评估技术的宇航用千兆数据电缆传输性能试验评价 . 航天器环境工程, 2019, 36(2): 156-160. doi: 10.12126/see.2019.02.009
    [9] 李鹏, 刘凯, 辛敏成, 赵楠, 邹田骥, 张海涛.  NAND闪存固态硬盘空间环境效应分析及测试系统设计 . 航天器环境工程, 2018, 35(6): 581-587. doi: 10.12126/see.2018.06.012
    [10] 闫琦, 黄小凯, 张立海, 刘守文.  低温泵现场使用可靠性评价及维修保障方法 . 航天器环境工程, 2017, 34(3): 336-342. doi: 10.12126/see.2017.03.019
    [11] 夏侨丽, 侯永青, 魏传锋.  可维修航天器的可靠性度量与建模方法 . 航天器环境工程, 2016, 33(6): 649-652. doi: 10.12126/see.2016.06.013
    [12] 张华1, 宗益燕1, 韦锡峰1, 陶强1, 王芳2地球同步轨道卫星多阶段任务可靠性建模 . 航天器环境工程, 2016, 33(4): 439-445. doi: 10.12126/see.2016.04.018
    [13] 司圣平, 王珏, 刘伟, 朱维.  遥控指令锁定电路的安全裕度分析 . 航天器环境工程, 2015, 32(4): 429-433. doi: 10.12126/see.2015.04.017
    [14] 舒适1, 周志涛1, 那顺布和1, 高雪峰2, 秦冉冉1, 赵美玲1基于在轨运行状态的卫星稳定性与可靠性研究 . 航天器环境工程, 2015, 32(2): 201-205. doi: 10.12126/see.2015.02.011
    [15] 路毅, 王再成, 张彬, 刘智斌, 郑会明, 杨春生.  卫星导电铜箔接地性能退化机理研究 . 航天器环境工程, 2014, 31(4): 397-400. doi: 10.12126/see.2014.04.010
    [16] 胡洪凯, 施蕾, 董暘暘, 刘波, 叶有时.  SRAM型FPGA空间应用的抗单粒子翻转设计 . 航天器环境工程, 2014, 31(5): 510-515. doi: 10.12126/see.2014.05.010
    [17] 张华, 宗益燕.  太阳电池阵运输过程的振动环境及可靠性分析 . 航天器环境工程, 2014, 31(5): 536-542. doi: 10.12126/see.2014.05.015
    [18] 窦威, 王勇, 任国华, 黄锡宁.  柔性收集室检漏技术可靠性分析与验证 . 航天器环境工程, 2014, 31(6): 670-675. doi: 10.12126/see.2014.06.020
    [19] 金历群, 吴 东, 王金华, 康纪军.  卫星总体可靠性安全性设计审查 . 航天器环境工程, 2010, 27(4): 500-504.
    [20] 祝耀昌.  可靠性试验及其发展综述 . 航天器环境工程, 2007, 24(5): 261-269.
  • 加载中
图(7) / 表ll (1)
计量
  • PDF下载量:  31
  • 文章访问数:  464
  • HTML全文浏览量:  178
文章相关
  • 中图分类号:  TB114.3
  • 收稿日期:  2019-07-26
  • 修回日期:  2020-03-16
  • 网络出版日期:  2020-04-15
  • 刊出日期:  2020-04-01

基于性能退化的失效阈值变化速率对产品可靠性的影响

doi: 10.12126/see.2020.02.006
    基金项目:  国家自然科学基金项目(编号:51475086);中央高校基本科研业务项目(编号:N162312001);北京卫星环境工程研究所创新基金项目(编号:CAST-BISEE 2019-019)
    作者简介:

    王新刚(1979—),男,博士学位,研究方向为机械零部件的动态可靠性设计。E-mail: wangxingang1217@126.com

  • 中图分类号: TB114.3

摘要: 性能退化分析是一种有效弥补数据不足的可靠性评估方法。采用线性变化过程描述机械产品的性能退化过程,考虑产品受极值冲击和运行冲击的情况,分析性能退化过程对突发失效的阈值影响。基于变失效阈值的性能退化过程和变失效阈值的冲击过程,建立突发失效和性能退化失效同时存在下的竞争失效可靠性模型。结合实例对可靠性模型进行分析,结果表明:若采用定失效阈值,忽略失效阈值的变化,产品的可靠度会被高估;在一定范围内,机械产品的可靠度随着失效阈值变化速率的增大而下降,随着性能退化量的变大而提升。

English Abstract

王新刚, 马瑞敏, 黄小凯, 等. 基于性能退化的失效阈值变化速率对产品可靠性的影响[J]. 航天器环境工程, 2020, 37(2): 137-142 doi:  10.12126/see.2020.02.006
引用本文: 王新刚, 马瑞敏, 黄小凯, 等. 基于性能退化的失效阈值变化速率对产品可靠性的影响[J]. 航天器环境工程, 2020, 37(2): 137-142 doi:  10.12126/see.2020.02.006
WANG X G, MA R M, HUANG X K, et al. Effect of failure threshold change rate on the anlysis of product reliability based on performance degradation[J]. Spacecraft Environment Engineering, 2020, 37(2): 137-142 doi:  10.12126/see.2020.02.006
Citation: WANG X G, MA R M, HUANG X K, et al. Effect of failure threshold change rate on the anlysis of product reliability based on performance degradation[J]. Spacecraft Environment Engineering, 2020, 37(2): 137-142 doi:  10.12126/see.2020.02.006
    • 科学技术的进步和生产技术的发展使产品的可靠性不断提高,在允许的时间内,很难获得产品性能研究所需的失效数据[1]。性能退化分析作为一种有效弥补机械产品失效数据不足的方法,可以间接预测产品的寿命[2]

      基于机械产品性能退化的可靠度分析,关键在于产品退化模型的建立以及失效标准即失效阈值的确定。通常采用退化轨迹和随机冲击过程这2种方法来描述机械产品的性能退化过程。退化轨迹模型是指采用一些应用较为成熟的数学模型来拟合产品性能退化量和时间的关系[3-4]。一般情况下产品的失效往往是由多种失效模式导致的,当突发失效或性能退化失效之一达到其失效阈值都会导致产品的失效[5]。Li等[6]分析了产品同时包含随机冲击过程与性能退化过程且两者相互独立下的可靠性。Wang等[7]考虑了外部冲击载荷对退化过程的影响,并在模糊退化数据下进行了退化分析,建立了基于模糊退化数据的可靠性模型。刘晓娟等[8]在外部冲击过程和内部性能退化过程同时存在的情况下,采用时变copula函数建立多个退化过程和冲击过程下的竞争失效可靠性模型,研究表明产品在工作过程中,其自身的性能退化过程和所受到的外部冲击过程相互影响。An等[9]分析了冲击载荷引起的性能退化增量与冲击载荷大小之间的关系,研究表明系统有抵抗微小载荷的能力,只有超过一定范围的冲击载荷才会对系统的性能退化量产生影响。Rafiee等[10]在假定性能退化过程为线性的前提下,研究了外界冲击分别为极值冲击和运行冲击等情况下对线性退化速率的影响。黄文平等[11]认为突发失效的失效阈值是变化的,在极值冲击过程和线性退化过程同时存在以及突发失效的阈值受到性能退化过程影响的前提下,建立了微引擎在变失效阈值下的竞争失效可靠度模型。

      在上述分析的基础上,考虑到产品本身的复杂性和工作环境的不稳定性,失效阈值为固定值不符合实际情况。随着性能退化量的逐渐增大,产品抵抗外载荷冲击的能力会逐渐减弱,即突发失效的失效阈值为随时间推移而降低的变量。温度等因素也会对材料的物理性能产生影响,因此性能退化过程的失效阈值亦为随时间推移而降低的变量。本文提出了在冲击载荷分别为极值冲击和运行冲击时,多个失效模式下的时变失效阈值可靠度模型,分析了性能退化模型和突发失效模型,建立了极值冲击和运行冲击下产品的可靠度模型,最后结合实例对所提出的模型进行了验证。

    • 图1所示,失效阈值随着时间推移逐渐降低,在工作应力的作用下,机械产品的性能退化量会随着工作的进行逐渐增加,在机械产品遭受外界冲击的时候(t1t2t3分别表示外界冲击载荷作用的时刻),如果外界冲击不能使得产品立即失效,就会造成产品性能退化量突然增加[12]。外界随机冲击载荷造成的产品性能退化量的突然增加和产品自身的性能退化量之和构成了产品总的性能退化量[13]。当总的性能退化量超过失效阈值H的时候,产品便会退化失效。

      图  1  性能退化失效过程

      Figure 1.  Failure process of performance degradation

      为了研究方便,令产品的失效阈值变化过程和性能退化过程都服从线性轨迹模型:

      $$X(t) = \phi + {\beta _1}t\text{,}{\rm{ }}{\beta _1}\sim N({\mu _{{\beta _1}}},{\delta _{{\beta _1}}})\text{;}$$ (1)
      $$H(t) = {H_0} - {\beta _0}t\text{,}{\rm{ }}{\beta _0}\sim N({\mu _{{\beta _0}}},{\delta _{{\beta _0}}})\text{。}$$ (2)

      式中:X(t)为t时刻产品自身的性能退化量;$\phi $为产品的性能退化量在0时刻的数值;β1为产品的性能退化速率,并且服从正态分布;H(t)为t时刻产品的失效阈值;H0为0时刻产品的失效阈值;β0为失效阈值的变化速率,且服从正态分布。

      假设产品所受的冲击载荷发生概率服从参数为λ的齐次泊松分布。令Yi表示在时间t内第i次冲击导致的性能退化量,其中Yi相互独立并且服从同一分布函数,即${Y_i}\sim N({\mu _Y},\delta _Y^2)$。则,产品在时间t内由外界载荷导致的性能退化增量为

      $$S(t) = \left( \begin{aligned} & 0, \;\;N(t) = 0\\ & \sum\limits_{i = 1}^{N(t)} {{Y_i}} ,\;N(t) \text{>} 0 \end{aligned} \right.\text{,}$$ (3)

      总的性能退化量为

      $${X_S}(t) = X(t) + S(t)\text{;}$$ (4)

      产品在时间t内不发生失效的概率为

      $$\begin{split} & {F_X}(x,t) = \sum\limits_{i = 0}^\infty {P(X(t) + S(t)\text{<} H(t)|N(t) = i)}P(N(t) = i) = \\ & \sum\limits_{i = 0}^\infty {\varPhi \left( {\frac{{H - {\mu _{{\beta _0}}}t - ({\mu _{{\beta _1}}}t + \varphi + i{\mu _y})}}{{\sqrt {{\delta _{{\beta _0}}}^2{t^2}{\rm{ + }}{\delta _{{\beta _1}}}^2{t^2} + i\delta _y^2} }}} \right)} \frac{{\exp ({\text{-}}\lambda t){{(\lambda t)}^i}}}{{i!}}\text{。}\\[-25pt] \end{split}$$ (5)
    • 图2所示,用wi表示外界冲击载荷的幅值,D(t)表示突发失效阈值,外界冲击载荷的幅值服从正态分布,发生率服从泊松分布。

      图  2  突发失效过程

      Figure 2.  Suddenly-occurred failure process

      产品的性能退化过程会影响突发失效阈值的变化过程:随着性能退化量的不断增加,产品抵抗外界冲击载荷的能力减弱,即产品的突发失效阈值逐渐降低,当外界冲击载荷超过失效阈值时,产品发生突发失效。突发失效阈值的变化情况同样适用于其他模型。为了研究方便,假设突发失效的阈值变化服从二阶段的线性变化过程,即在产品性能退化量小于L时以β2的速率递减,在产品性能退化量大于L时以β3的速率递减,β2β3均为服从正态分布函数的变量。

      在时间t内,突发失效的阈值变化情况如图2所示。设在u时刻性能退化量正好到达L,则在时间t内突发失效的阈值分布函数为

      $$D\left( t \right) = {D_1} - {\beta _2}u - {\beta _3}\left( {t - u} \right)\text{,}$$ (6)

      式中:${\beta _2}\sim N({\mu _{{\beta _2}}},\delta _{{\beta _2}}^2)$${\beta _3}\sim N({\mu _{{\beta _3}}},\delta _{{\beta _3}}^2)$

    • 在时间t内,考虑产品的性能退化失效和突发失效2种模式同时存在的条件下,产品的可靠度可以表示为

      $$R(t) = \sum\limits_{i = 0}^\infty {R(t|N(t) = i)P(N(t) = i)\text{,}} $$ (7)

      式中:R(t|N(t) = i)为t时刻发生i次冲击下产品的可靠度;P(N(t) = i)为时间t内外界冲击载荷发生i次的概率。

      R(t)的计算较为复杂,为了分析方便分以下3种情况讨论:

      1)在时间t内,产品只发生了性能退化,则其可靠度为

      $$\begin{split} &{R_1}(t) = P(X(t)\text{<}H(t)|N(t) = 0)P(N(t) = 0) = \\ & \varPhi \left( {\frac{{H - {\mu _{{\beta _0}}}t - ({\mu _{{\beta _1}}}t + \varphi )}}{{\sqrt {{\delta _{{\beta _0}}}^2{t^2}{\rm{ + }}{\delta _{{\beta _1}}}^2{t^2}} }}} \right)\exp ({\text{-}}\lambda t)\text{。} \end{split}$$ (8)

      2)在时间t内,总的性能退化量小于L且有外界冲击载荷作用于产品,则产品可靠度为

      $$\begin{split} & {R_2}(t) = \sum\limits_{i = 1}^\infty {P(X(t) + S(t)\text{<}L(t)|N(t) = i)} \cdot \\ & P(N(t) = i) \cdot P\left(\bigcap\limits_{i = 1}^\infty {\left\{ {{W_i} \text{<}D(t)} \right\}} \right) = \\ & \sum\limits_{i = 1}^\infty {\varPhi \left( {\frac{{L - {\mu _{{\beta _0}}}t - ({\mu _{{\beta _1}}}t + \varphi + i{\mu _Y})}}{{\sqrt {{\delta _{{\beta _0}}}^2{t^2}{\rm{ + }}{\delta _{{\beta _1}}}^2{t^2} + i\delta _Y^2} }}} \right)} \cdot \\ & \varPhi {\left( {\frac{{{D_0} - {\rm{(}}{\mu _{{\beta _2}}}t{\rm{ + }}{\mu _w}{\rm{)}}}}{{\sqrt {{\delta _w}^2 + \delta _{{\beta _2}}^2{t^2}} }}} \right)^i}\frac{{\exp ({\text{-}}\lambda t){{(\lambda t)}^i}}}{{i!}}\text{。} \end{split}$$ (9)

      3)在时间t内,总的性能退化量大于L且小于H,且有外界冲击载荷作用于产品,则产品可靠度为

      $$ \begin{split} {R_3}(t) = & \sum\limits_{i = 1}^\infty {\sum\limits_{k = 0}^i {P({X_{S1}}(u)\text{<}L(u)|N(u) = k)} }\cdot \\ & P(L(t - u)\text{≤} {X_{S2}}(t - u)\text{<}H(t - u)|N(t - u) = i - k)\cdot \\ & P(N(u) = k) \cdot P(N(t - u) = i - k) \cdot \\ & P\left(\bigcap\limits_{i = 1}^\infty {\left\{ {{W_i}\text{<}D(t)} \right\}} \right)\text{。}\\[-20pt] \end{split} $$ (10)

      为了分析方便需要假定几个参数,把性能退化量刚好到达L的时刻记为u,时间t内的冲击次数记为i,在时间u内冲击的次数记为k,那么[u, t]区间内的冲击次数为(i k)。在时间u内,把性能退化量小于L且外界冲击次数为k次条件下产品不发生性能退化失效的概率函数记为FX(k, x, u),则

      $$\begin{split} {F_X}(k,x,u) = & \sum\limits_{{\rm{i = }}0}^\infty {P{\rm{(}}X(u) + S(u)\text{<}L(u){\rm{|}}N(u) = k{\rm{) = }}} \\ & \sum\limits_{i = 0}^\infty {\varPhi \left( {\frac{{L - {\mu _{{\beta _0}}}u - ({\mu _{{\beta _1}}}u + \varphi + k{\mu _y})}}{{\sqrt {{\delta _{{\beta _0}}}^2{u^2}{\rm{ + }}{\delta _{{\beta _1}}}^2{u^2} + k\delta _y^2} }}} \right)}\text{,} \end{split}$$ (11)

      产品在时间t内的可靠度为

      $$\begin{split} {R_3}(t) = & {\sum\limits_{i = 1}^\infty {\sum\limits_{k = 0}^i {\int_0^t {\varPhi \left( {\frac{{{D_0} - {\rm{[}}{\mu _{\beta 2}}u{\rm{ + }}{\mu _{\beta 3}}{\rm{(}}t - u{\rm{) + }}{\mu _w}]}}{{\sqrt {{\delta _w}^2 + \delta _{\beta 2}^2{u^2} + \delta _{\beta 3}^2{{(t - u)}^2}} }}} \right)} } } ^i} \cdot \\ & \frac{{\exp [{\text{-}}\lambda (t - u)]{{[\lambda (t - u)]}^{i - k}}}}{{(i - k)!}} \cdot \\ & \varPhi \left( {\frac{{H - L - {\mu _{{\beta _0}}}(t - u) - [{\mu _{{\beta _1}}}(t - u) + \varphi + (i - k){\mu _Y}]}}{{\sqrt {{\delta _{{\beta _1}}}^2{{(t - u)}^2} + \delta _{{\beta _0}}^2{t^2} + (i - k)\delta _{_Y}^2} }}} \right) \cdot \\ & \exp ({\text{-}}\lambda u)\frac{{{{(\lambda u)}^k}}}{{k!}}{f_X}(k,x,u){\rm{d}}u\text{。}\\[-15pt] \end{split}$$ (12)

      机械产品在外界冲击载荷为极值冲击的前提下,在工作时间t内其工作状态是上述3种情况之一,因此时间t内产品的可靠度为

      $$R(t) = {R_1}(t) + {R_2}(t) + {R_3}(t)\text{。}$$ (13)
    • 图3所示,在运行冲击模型中,如果有连续m个外载荷大于某个固定阈值D2或者大于突发失效阈值D(t)的时候,就会发生运行冲击失效[14]

      图  3  运行冲击模型

      Figure 3.  Operation shock model

      为了分析方便首先定义Ui为在i次外界冲击载荷中不会有连续m个冲击载荷大于某个固定阈值D2的概率,Ui应该满足

      $$ \begin{split}{U_i} = & Q{U_{i - 1}} + PQ{U_{i - 2}} + {P^2}Q{U_{i - 3}} + \cdots + \\ & {P^{m - 1}}Q{U_{i - m}}\text{。}\end{split}$$ (14)

      式中:P为冲击载荷大于固定阈值D2并且小于失效阈值D1的概率;Q为冲击载荷小于固定阈值D2的概率。PQ的表达式分别为:

      $$\begin{split} P = & P({W_i}\text{>}{D_2}|{W_i}\text{<}{D_1}) = \\ & \frac{{P({D_2}\text{<}{W_i}\text{<}{D_1})}}{{P({W_i} \text{<}{D_1})}} = \frac{{{F_W}({D_1}) - {F_W}({D_2})}}{{{F_W}({D_1})}}\text{;}\\ & \!\!\!\!\!\!\!\!\!\!\!\!\!Q = 1 - P = P({W_i}\text{<}{D_2}|{W_i}\text{<}{D_1}) = \\ & \frac{{P({W_i}\text{<}{D_2})}}{{P({W_i}\text{<}{D_1})}} = \frac{{{F_W}({D_2})}}{{{F_W}({D_1})}}\text{。} \end{split}$$ (15)

      分3种情况考虑产品在时间t内的可靠度:

      1)在时间t内,产品只发生了性能退化,则其可靠度为

      $$\begin{split} {R_1}(t) = & P(X(t)\text{<}H|N(t) = 0)P(N(t) = 0) = \\ & \varPhi \left( {\frac{{H - {\mu _\beta }_{_0}t - ({\mu _\beta }_{_1}t + \varphi )}}{{{\delta _{{\beta _1}}}t}}} \right)\exp ({\text{-}}\lambda t)\text{。} \end{split}$$ (16)

      2)在时间t内,总的性能退化量小于L且有外界冲击载荷作用于产品,则产品可靠度为

      $$\begin{split} {R_2}(t) = & \sum\limits_{i = 1}^\infty {P{\rm{(}}X(t) + S(t)\text{<} H(t){\rm{|}}N(t) = i{\rm{)}}} \cdot \\ & P{\rm{(}}N(t) = i{\rm{)}}P\left(\bigcap\limits_{i = 1}^\infty {\left\{ {{W_i}\text{<}D(t)} \right\}} \right){U_i}{\rm{ = }}\\ & \sum\limits_{i = 1}^\infty {\varPhi \left( {\frac{{L - {\mu _{{\beta _0}}}t - ({\mu _{{\beta _1}}}t + \varphi + i{\mu _Y})}}{{\sqrt {{\delta _{{\beta _1}}}^2{t^2} + i\delta _Y^2} }}} \right)} \cdot \\ & \varPhi {\left( {\frac{{{D_1} - {\rm{(}}{\mu _{{\beta _2}}}t{\rm{ + }}{\mu _w}{\rm{)}}}}{{\sqrt {{\delta _w}^2 + \delta _{{\beta _2}}^2{t^2}} }}} \right)^i}\frac{{\exp ({\text{-}}\lambda t){{(\lambda t)}^i}}}{{i!}}{U_i}\text{。} \end{split}$$ (17)

      3)在时间t内,总的性能退化量大于L且小于H,且有外界冲击载荷作用于产品,则产品可靠度为

      $$\begin{split} {R_3}(t) = & {\sum\limits_{i = 1}^\infty {\sum\limits_{k = 0}^i {\int_0^t {\varPhi \left( {\frac{{{D_1} - {\rm{[}}{\mu _{{\beta _2}}}u{\rm{ + }}{\mu _{{\beta _3}}}{\rm{(}}t - u{\rm{) + }}{\mu _w}]}}{{\sqrt {{\delta _w}^2 + \delta _{{\beta _2}}^2{u^2} + \delta _{{\beta _3}}^2{{(t - u)}^2}} }}} \right)} } } ^i} \cdot \\ & \frac{{\exp [{\text{-}}\lambda (t - u)]{{[\lambda (t - u)]}^{i - k}}}}{{(i - k)!}} \cdot \\ & \varPhi \left( {\frac{{H - L - {\mu _{{\beta _0}}}(t - u) - [{\mu _{{\beta _1}}}(t - u) + \varphi + (i - k){\mu _Y}]}}{{\sqrt {{\delta _{{\beta _1}}}^2{{(t - u)}^2} + \delta _{{\beta _0}}^2{{(t - u)}^2} + (i - k)\delta _Y^2} }}} \right) \cdot \\ & \exp ({\text{-}}\lambda u)\frac{{{{(\lambda u)}^k}}}{{k!}}{f_X}(k,x,u){U_i}{\rm{d}}u\text{。} \\[-15pt] \end{split}$$ (18)

      综合以上3种情况,产品在时刻t的可靠度为

      $$R(t) = {R_1}(t) + {R_2}(t) + {R_3}(t)\text{。}$$ (19)
    • 根据美国Sandia国家实验室进行的一项实验研究发现,一种包括正交梳状驱动执行机构和相互机械连接的旋转齿的微型发动机,其梳齿传动的线性位移通过销传递到齿轮上。微引擎参数如表1所示,随着工作过程的进行,齿轮与销连接的表面会有明显的磨损。产品正常工作时的载荷和外载荷的冲击导致磨损,外部冲击会导致微型发动机的零部件突然断裂,这2个失效模式都可能导致微型发动机失效。

      表 1  微引擎可靠性分析参数

      Table 1.  Parameter values for micro-engine reliability analysis

      参数数据来源
      D0/GPa1.2假设
      D1/GPa1.4文献[15]
      D2/GPa1.3假设
      H0/μm30.001 25文献[15]
      L/μm37×10-4假设
      β0/μm3${\mu _{\beta 0}}$$ = 1.2 \times {10^{{\text{-9}}}}$
      ${\sigma _{\beta 0}}$$ = 1 \times {10^{{\text{-9}}}}$
      假设
      β1/μm3${\mu _{\beta 1}}$$ = 8.48 \times {10^{{\text{-9}}}}$
      ${\sigma _{\beta 1}}$$ = 6 \times {10^{{\text{-10}}}}$
      文献[15]
      β2/μm3${\mu _{\beta 2}}$$ = 1 \times {10^{{\text{-5}}}}$
      ${\sigma _{\beta 2}}$$ = 1 \times {10^{{\text{-6}}}}$
      假设
      β3/μm3${\mu _{\beta 3}}$$ = 1 \times {10^{{\text{-5}}}}$
      ${\sigma _{\beta 3}}$$ = 1 \times {10^{{\text{-6}}}}$
      假设
      λ/转5×10-5假设
      Yi/μm3${\mu _Y}$$ = 1 \times {10^{{\text{-4}}}}$
      ${\sigma _Y}$$ = 2 \times {10^{{\text{-5}}}}$
      文献[15]
      Wi/μm3${\mu _W} = 1.2$
      ${\sigma _W} = 0.2$
      假设
      P0.131 0假设
      Q0.869 0假设
    • 产品在极值冲击的工作环境下,将表1中的数据代入式(13),得到微引擎的可靠度函数,如图4图5所示。

      图  4  极值冲击情况下不同失效阈值变化速率下的产品可靠度

      Figure 4.  Product reliability at different failure threshold degradation rates in extreme operation shock

      根据式(13)可知,突发失效阈值的变化速率${\mu _{\beta 2}}$${\mu _{\beta 3}}$都对产品可靠度有影响。为了便于分析,分别给出了3组数据。通过图4可以看到:不考虑突发失效阈值的变化,产品的可靠度会被高估;当失效阈值变化速率变大时,可靠度曲线向左偏移,说明在一定范围内随着突发失效阈值变化速率的增大,产品的可靠度降低。

      根据上文分析可以得知,当性能退化量达到L时,产品抵抗冲击能力减弱,即突发失效阈值变化速率发生改变。在其他参数和工作时间不变的情况下,通过改变L的大小来研究其对产品可靠度的影响。图5分别给出了3组不同的L值,可以看到:在一定范围内随着L的逐渐变大,产品的可靠度也随之提高。这是因为L越小,突发失效阈值变化速率增大的条件越容易被满足。

      图  5  极值冲击情况下不同L值下的产品可靠度

      Figure 5.  Product reliability for different values of L in extreme operation shock

    • 在外界冲击为运行冲击的情况下,将表1中的数据代入式(19),得到微引擎的可靠度函数,如图6图7所示。

      图  6  运行冲击情况下不同失效阈值变化速率下的产品可靠度

      Figure 6.  Product reliability at different failure threshold degradation rates in operation shock

      失效阈值变化速率对运行冲击模型可靠度的影响与对极值冲击模型可靠度的影响相似。在其他参数不变的条件下,给出了3种失效阈值的变化速率。通过图6可以看到,在一定范围内,随着失效阈值变化速率的增加,产品的可靠度降低。

      数值L对运行冲击模型可靠度的影响与对极值冲击模型可靠度的影响相似。在其他参数不变的条件下,给出了3组不同的L值。通过图7可以看到,在一定范围内随着L变大,产品的可靠度提高。

      图  7  运行冲击情况下不同L值下的产品可靠度

      Figure 7.  Product reliability for different values of L in operation shock

    • 本文研究了在外界冲击为极值冲击和运行冲击的情况下,突发失效与性能退化同时存在下机械产品失效阈值变化速率对产品可靠度的影响。

      结果表明:1)忽略机械产品突发失效阈值的变化,机械产品的可靠度会被高估;2)突发失效阈值的变化速率会影响机械产品的可靠度,在一定范围内,随着变化速率的逐渐增大,机械产品的可靠度逐渐下降;3)决定失效阈值变化速率的性能退化量L对产品可靠度有影响。在一定范围内,随着性能退化量L变大,产品的可靠度提高。

参考文献 (15)

目录

    /

    返回文章
    返回