3. Duane模型的数学形式
Duane模型的核心思想可以用一句话概括:在可靠性增长试验中,累积失效率随累积试验时间的增加而按幂律下降。这个看似简单的假设,构成了后续所有参数估计、增长评估和趋势判断的数学基础。理解Duane模型,关键不在于记住公式,而在于理解公式中每一个量在试验现场对应什么、幂指数为什么能刻画学习效应、以及为什么取对数后模型会变成一条直线。
一、累积失效率的定义与幂律关系
设累积试验时间为 T(通常以小时、循环次数或等效任务时间计),到时刻 T 为止累计发生的失效数为 N(T)。Duane模型定义累积失效率为累计失效数与累积试验时间之比:
C(T) = N(T) / T
注意这里的 C(T) 是“累积平均”意义上的失效率,而不是瞬时失效率。Duane通过观察多类产品的可靠性增长试验数据发现,C(T) 与 T 在双对数坐标下近似呈直线关系,即:
C(T) = a · T−m, 0 < m < 1,a > 0
这就是Duane模型的幂律形式。它表明:随着累积试验时间增加,累积失效率以幂函数形式递减。指数为负,说明时间越长、累积失效率越低;指数绝对值小于1,说明下降速度逐渐放缓,而不是无限加速下降。
二、模型参数 a 与 m 的含义
参数 a 是尺度参数。当 T = 1(单位时间)时,C(1) = a,因此 a 可理解为“单位累积试验时间处的累积失效率水平”。它反映了增长过程的起始高度或整体失效率量级。在工程上,a 越大,说明早期失效密度越高,或产品初始可靠性水平越低。
参数 m 是形状参数,也称增长率或学习因子。它决定累积失效率随时间的下降速率。对幂律式两边取对数可得:
ln C(T) = ln a − m ln T
因此 m 就是双对数坐标下直线的斜率绝对值。若 m = 0,则 C(T) 恒等于 a,表示没有可靠性增长;若 m 接近1,表示累积失效率下降很快,增长效果显著。实际工程中,m 通常落在0.2~0.7之间。例如某型机载电子设备在增长试验中拟合得到 m = 0.45,意味着累积试验时间每增加一个数量级,累积失效率约降至原来的 10−0.45 ≈ 0.355,即下降约64.5%。
三、对数线性化表达与参数估计
Duane模型最重要的实用价值之一,是可以通过对数变换转化为线性回归问题。令:
y = ln C(T), x = ln T, b = ln a
则幂律关系变为:
y = b − m x
这是一条截距为 b、斜率为 −m 的直线。给定一组试验数据 (Ti, Ni),先计算 Ci = Ni / Ti,再取对数得到 (xi, yi),然后用最小二乘法估计斜率和截距,最后反解出 m 和 a。具体地:
m = − [ nΣ(x_i y_i) − Σx_i Σy_i ] / [ nΣ(x_i²) − (Σx_i)² ]
b = ( Σy_i + m Σx_i ) / n
a = exp(b)
下面给出一个教学示例。假设某产品在累积试验时间 100、200、400、800、1600 小时处,累计失效数分别为 12、18、25、33、42。计算累积失效率并取自然对数,得到下表(数值为示例假设,用于说明计算流程):
| 累积时间 T | 累计失效 N | C=N/T | x=lnT | y=lnC |
|---|---|---|---|---|
| 100 | 12 | 0.1200 | 4.605 | −2.120 |
| 200 | 18 | 0.0900 | 5.298 | −2.408 |
| 400 | 25 | 0.0625 | 5.991 | −2.773 |
| 800 | 33 | 0.0413 | 6.685 | −3.187 |
| 1600 | 42 | 0.0263 | 7.378 | −3.639 |
用上述数据做最小二乘拟合,可得到斜率约为 −0.50,即 m ≈ 0.50;截距 b ≈ 0.18,故 a ≈ 1.20。于是模型为 C(T) ≈ 1.20 T−0.50。若要在图上判断拟合好坏,应观察双对数散点是否大致沿直线分布;若明显弯曲,说明幂律假设可能不成立,或试验过程中存在重大设计变更、试验条件变化等分段效应。
四、从模型到工程判断
在实际可靠性增长管理中,Duane模型的对数线性形式不仅用于拟合历史数据,还用于外推和决策。例如,若当前累积试验时间为 T0,希望累积失效率降到目标值 Ctarget,则可由 Ctarget = a T−m 反解所需累积时间:
T = ( a / Ctarget )1/m
仍以上述示例 a ≈ 1.20、m ≈ 0.50 为例,若目标累积失效率为 0.01,则所需累积试验时间约为 (1.20/0.01)2 = 14400 小时。这个结果说明:在 m = 0.5 的增长速率下,把累积失效率从约0.026降到0.01,需要再投入大量试验时间。工程上可据此评估试验周期、资源需求和增长计划是否可行。
需要特别强调的是,Duane模型描述的是累积平均失效率,而不是瞬时失效率。瞬时失效率 λ(T) 与累积失效率的关系为 λ(T) = dN/dT。在幂律假设下,可以推导出 λ(T) = (1−m) C(T),即瞬时失效率始终低于累积失效率,且两者相差一个常数因子 (1−m)。这一点在把模型结果用于可靠性预计或寿命评估时必须注意,不能把 C(T) 直接当作当前时刻的瞬时失效率使用。
总结来说,Duane模型的数学形式包含三层:第一层是幂律关系 C(T) = a T−m,描述累积失效率随累积试验时间下降的总体规律;第二层是参数 a 与 m 的工程含义,分别对应初始水平和增长速率;第三层是对数线性化 y = b − m x,使参数估计转化为线性回归问题。掌握这三层,就掌握了Duane模型从数据到判断的基本链路。