3. 常见寿命分布及其竞争失效含义
在竞争失效建模中,寿命分布不是一条“拟合曲线”那么简单,它决定了各失效机理在时间轴上的风险形状,也决定了外推寿命时的可信边界。指数分布、威布尔分布、对数正态分布、伽马分布和极值分布,是工程可靠性中最常用的五类寿命模型。理解它们的适用条件与参数含义,是正确建立竞争失效模型的前提。
指数分布:恒定风险的基准模型
指数分布的失效密度为 f(t)=λe−λt,可靠度 R(t)=e−λt,失效率 h(t)=λ 为常数。参数 λ 是单位时间的失效速率,MTTF=1/λ。它适用于失效由随机外部冲击主导、且没有明显磨损或老化趋势的阶段,例如电子元器件的偶然失效期、某些软件缺陷触发、随机过应力事件。在竞争失效中,指数分布常作为“无记忆”基准:若多个独立失效机理都近似恒定风险,则系统总失效率为各机理失效率之和,MTTF 为各机理 MTTF 的并联倒数关系。但要注意,指数分布不能描述早期失效下降或磨损上升,若误用于老化阶段,会严重高估长寿命可靠性。
威布尔分布:形状参数决定失效机理
威布尔分布是竞争失效建模中最灵活的分布。其可靠度为 R(t)=exp[−(t/η)β],其中 η 为尺度参数(特征寿命,t=η 时 R≈0.368),β 为形状参数。β<1 对应早期失效(失效率递减),β=1 退化为指数分布(恒定风险),β>1 对应磨损或老化(失效率递增)。在竞争失效中,不同失效机理往往对应不同 β:例如焊点热疲劳 β 常在 2~4,而随机过电应力 β 接近 1。若把多个机理混在一起拟合,得到的 β 可能是“平均形状”,掩盖真实机理。因此工程上常按机理分组,分别估计 β 和 η,再合成竞争风险。
对数正态分布:退化与裂纹扩展的常用模型
对数正态分布假设寿命的对数服从正态分布,其密度为 f(t)=1/(tσ√(2π))·exp[−(ln t−μ)²/(2σ²)]。参数 μ 是对数均值,σ 是对数标准差。它的失效率先升后降,适合描述疲劳裂纹扩展、某些腐蚀和磨损类退化过程,也常用于半导体器件和机械零件的寿命分析。在竞争失效中,对数正态分布适合那些“损伤累积到阈值才失效”的机理,其寿命分散性由 σ 控制:σ 越大,寿命分布越宽,早期失效与长寿命尾部同时存在。与威布尔相比,对数正态对尾部外推更敏感,需谨慎使用。
伽马分布:累积损伤过程的自然描述
伽马分布常被用作“到达失效前需累积若干次损伤”的模型。其密度含形状参数 k 和尺度参数 θ,均值 kθ,方差 kθ²。当 k 为整数时,伽马分布可理解为 k 个独立指数事件之和,因此特别适合描述多阶段退化、累积冲击损伤和维修次数驱动的失效。在竞争失效中,若某机理的失效由多次随机冲击累积触发,伽马过程是自然选择;其失效率随 k 增大而逐渐趋于恒定,介于指数与威布尔之间。参数 k 反映损伤累积的“阶段数”,θ 反映单次损伤的平均时间尺度。
极值分布:最弱环节与极端载荷
极值分布描述一组随机变量中最大值或最小值的分布。在可靠性中,最弱环节模型(串联系统)的寿命往往由最小值决定,其渐近分布为 Weibull 型极值分布;而极端载荷、最大应力导致的失效则可能服从 Gumbel 型极值分布。极值分布的价值在于:当失效由“最薄弱点”或“最大一次冲击”控制时,它比普通寿命分布更准确地刻画尾部行为。在竞争失效中,极值分布常用于解释为什么某些系统在极短或极长时间内失效概率异常高——这正是尾部竞争的直接体现。
参数解释与竞争失效适用条件
把五类分布放在竞争失效框架下,关键不是记住公式,而是判断“哪个机理对应哪个分布”。指数分布适合随机冲击型机理,参数 λ 直接给出风险水平;威布尔分布适合磨损、疲劳、腐蚀等有明确退化趋势的机理,β 是机理指纹,η 是时间尺度;对数正态分布适合裂纹扩展和损伤累积到阈值的机理,σ 控制分散性;伽马分布适合多阶段累积冲击,k 是阶段数;极值分布适合最弱环节和极端载荷,关注尾部而非中心。实际建模时,应先按失效机理分组,分别选择分布并估计参数,再通过 hsys(t)=Σhi(t) 合成系统风险。若各机理独立,系统可靠度为各机理可靠度之积;若存在共因或耦合,则需引入 copula 或共享脆弱性模型。
一个教学示例:两种机理的竞争
假设某电源模块有两种主要失效机理:随机过压冲击(近似指数,λ1=2×10−5/h)和电解电容磨损(威布尔,β=2.5,η=80000h)。在 t=20000h 时,指数机理的失效率仍为 2×10−5/h;威布尔机理的失效率为 h2(t)=(β/η)(t/η)β−1=(2.5/80000)×(20000/80000)1.5≈3.125×10−5×0.125≈3.9×10−6/h。此时总失效率约 2.39×10−5/h,指数机理仍占主导。到 t=60000h 时,威布尔失效率升至约 2.5/80000×(0.75)1.5≈3.125×10−5×0.6495≈2.03×10−5/h,与指数机理相当,竞争效应开始显现。这个示例说明:不同机理的风险随时间变化,竞争失效模型必须保留时间维度,不能只用单一 MTTF 概括。
选择分布时的检查要点
- 机理一致性:β 或 σ 是否与失效物理预期一致,例如疲劳 β 通常大于 1,随机冲击 β 接近 1。
- 数据支撑:尾部外推是否有足够失效样本,截尾数据是否已正确处理。
- 独立性与耦合:各机理是否独立,是否存在共因失效或共享载荷。
- 参数可解释性:η、λ、μ、k 等参数是否能对应到工程可理解的量,如特征寿命、失效率、阶段数。
- 竞争合成方式:串联竞争用可靠度乘积,风险叠加用失效率求和,耦合时需引入 copula 或脆弱性模型。
掌握这五类分布的风险形状和参数含义,才能在竞争失效建模中做出有依据的选择,而不是把数据丢给软件自动拟合。分布是机理的语言,选对分布,模型才具备工程解释力和外推可信度。