MOSFET 寿命预测 SPICE 集成 — averaged 电热 + aging + lifetime 闭环

不仿真 PWM 开关,而是 averaged 模型 —— 电压源 / 受控源代替开关网络。损耗按公式估算:

其中 $Rds,on$ 同时包含温度依赖和aging 退化:

$\alpha$ 是温度系数(典型 0.5%/°C),$\Delta Rds,on$ 是 aging 引起的永久 increment。

因为 averaged model 视为稳态,不需要 RC ladder(transient response),只用 R 网络:

• $Rjc=0.8{}^{\circ }C/W$(junction-to-case)
• $Rca=10{}^{\circ }C/W$(case-to-ambient)

注意这与 topic-electro-thermal-simulation 的瞬态 Cauer/Foster 模型不同 —— 后者用 5 阶 RC 捕获 μs-ms 瞬态,这里只关心稳态 $Tj$(用于长期 lifetime 预测)。

Coffin-Manson 计算单个 thermal cycle 的失效循环数:

IRFP340 实验拟合参数:

• $\delta =-5.2776$
• $A1=4.9283\times 10^{13}$
• $Ea=814\mathrm{meV}$(activation energy)

Miner's rule 累积 stress:

$Q=1$ 表示寿命到达。

两种 aging 模型,paper 推荐 stress-related(更通用):

意味着 $Rds,on$ 累积上限是 0.2 × $Rinit$(即 20% increase)—— 但 IRFP340 EoL 准则更严:$\Delta Rds,on=0.05\Omega \approx 0.12Rinit$,在 $Q\approx 0.6$ 时已经到 EoL。

关键:这是自加速反馈环 — 过去的 stress 让 $Rds,on$ 升,使 $Pcond$ 升,使 $Tj$ 升,使后续 stress 进一步加速。

paper 用 IRFP340 在 4 种温度循环 stress 下做实验:

Coffin-Manson 拟合在所有 4 case 误差 $<10\%$,验证模型可信。

寿命曲线形状跟 stress 水平强相关——高 stress 后期指数爆破,低 stress 几乎线性,这直接决定在线监测的预警窗口长度:

• 高 stress(Case A):$Rds,on$ 早期平,后期指数加速爆破 — 自加速效应明显
• 低 stress(Case D):$Rds,on$ 几乎线性增长 — 在 EoL 之前都不会突变

工程含义:用 in-field $Rds,on$ 监测预警,高 stress 工况预警窗口短(几百 cycle),低 stress 工况预警窗口长(几千 cycle)。

传统 rainflow counting 需要完整波形(offline),无法在线。半 cycle counting 用每 5s 采样 $Tj$,把上升 / 下降斜率各算作 1 half cycle:

• $\Delta Tj=|Tj,\text{now}-Tj,\text{prev}|$
• $Tm=(Tj,\text{now}+Tj,\text{prev})/2$

精度比 rainflow 差 ~ 10%,但可以在线 SPICE 直接实现,这是关键工程价值。

Smaller-Step rainflow 用 LTspice 的 behavioral source 可以直接在仿真里跑,关键是用 SH() 采样 + idt() 积分凑出 Miner sum——代码不长但每行都贴 LTspice 函数语义:

; 采样 + 归一化(LTspice 限制 SH function 输入 < 10V)
A2 in clk out (sample-and-hold,每 5s 触发)
V_sampling = SH( V(T_j_mos) / 30 )

; Δ T_j 和 T_m 计算(乘回 30 抵消采样时的 /30 归一化,还原真实 ℃/K,再喂 Coffin-Manson)
V_ΔTj = 30 · abs( V(sampling) - delay(V(sampling), 5s) )
V_Tm  = 30 · ( V(sampling) + delay(V(sampling), 5s) ) / 2

; Coffin-Manson N_f
V_Nf = pow( V_ΔTj, delta ) · A_1 · exp( E_a / (k · V_Tm) )

; Miner Q 累积(half-cycle = 0.5 个 cycle)
V_Q = idt( 0.5 / V_Nf )

idt() 是 LTspice 的连续积分函数,直接实现 Miner sum;delay() 实现采样比较;SH() 是 sample-and-hold。

Aging 模型用到的 LTspice 函数集中在 6 个,理解它们的语义比 copy 代码更重要——下面是 cheat sheet:

paper Table 1 完整列出了 7 个函数 + 模型中各自用途。

仿真 120 ks(= 10 年实时压缩,无 PWM 详细)。环境温度和载阻分别用:

• $Ta=\text{random}(t,0.1)\times 5+\sin (30,30.41\mathrm{mdaily})+\sin (20,83.33\mu \mathrm{Hz}\mathrm{annual})$
• $RL$ = PWL file,每 5 s 一变(模拟 4.8 hr 实际工况变化)

仿真时间 46 秒(Intel i7-7600U)。

10 个 SPICE 节点同时可观察:

无 aging 模型(传统方法)预测:13 年到 EoL($Q=1$)。

含 aging 模型(本 paper):9.67 年到 EoL —— 提早 25%。

物理原因:aging 让 $Rds,on$ 缓慢升 → 损耗渐升 → $Tj$ 升 → 后续 thermal cycle 更狠 → aging 加速。这是正反馈失稳,被 paper 的 stress-feedback 模型捕获。

工程含义:传统不含 aging 的寿命预测 over-promise 25%,产品保修期会被 underestimate failure rate。

用热阻 $Rth$ 监测老化是直觉上最合理的选择,但实测中 $Rth$ 随 aging 漂移只有 5-10%,而 $Rds,on$ 漂移 12-20%,敏感度高 2-3×。$Rds,on$ 还有额外优势:常规电路就能在线算 $Vds/Id$,无需打断运行;而 $Rth$ 需要独立热测试台,难以连续监测。

switched model 与 averaged model 在短时(< 8 ks)仿真精度差 < 5%,但长期(10 ks,约 10 年)仿真 switched 慢 5600×——同一任务 switched 跑 13 分钟,averaged 只需 0.15 秒。switched 仿真速度随频率线性恶化:50 kHz 比 10 kHz 再慢 5×;而 averaged model 完全频率无关。10 年 mission profile 在普通笔记本上可行的唯一路径是 averaged model。

Coffin-Manson 3 个参数($\delta$、$A1$、$Ea$)是对特定器件 + 特定封装 + 特定失效机理的经验拟合。IRFP340 的 $\delta =-5.2776$ 对另一型号可能偏差一个数量级,跨型号直接挪用没有任何保证。正确做法:对目标器件做至少 4 个不同 stress level($\Delta Tj$、$Tm$ 组合)的加速老化实验再拟合——少于 4 组实验时参数欠定,误差无界。

Miner Q = 1 是统计均值,不是确定性门限。器件间制造离散直接传导到寿命上:同批 IRFP340 EoL 分布宽度 $\pm 30\%$。工程含义:如果目标寿命是 $L$ 年,$Q=1$ 只够均值器件——最坏 case 器件($-1\sigma$)可能在 $0.7L$ 就到 EoL。多器件系统(如 6-phase 逆变器)的第一个失效服从最坏 case,正式产品应将 Coffin-Manson 结果送入 Weibull 分布,给出 B1/B10 而非单点预测。

去掉 aging 反馈(传统方法)让 $Rds,on$ 保持初值不升,导致 Miner sum 用偏低的 stress 积分,最终寿命预测 over-promise 25%:同一器件含 aging 预测 9.67 年,无 aging 预测 13 年。对于 8-15 年设计寿命的汽车 / 工业产品,这个误差足以让现场失效率在保修期内超过预算——特别是 aging 正反馈在后期指数加速的高 stress 工况。

Coffin-Manson + Miner 只捕获渐进退化失效——wire bond 疲劳、$Rds,on$ 漂移、die attach 裂纹。瞬发失效(gate oxide breakdown、雪崩击穿、ESD、latch-up)完全不在模型覆盖范围内。实际产品可靠性评估是两条腿:lifetime model 算退化,FMEA + FIT 率算灾难性失效,两者独立、不可互替。详见 topic-fit-fmeda-calculation。

IRFP340 的 Coffin-Manson 参数来自 Si 器件的 wire bond 铝线疲劳。SiC MOSFET 主要退化机理是 BPD(basal plane dislocation)迁移——在正向电流下 BPD 从 SiC 外延层扩展,引起 body diode $Vf$ 升,而非 $Rds,on$ 主导劣化。物理机理不同导致 $\delta$ 和 $Ea$ 差异可达一倍以上;直接套用 Si 参数会大幅高估 SiC 寿命。SiC 器件需要 SiC 专属老化实验和参数集。

Coffin-Manson 的 $Nf\propto \Delta T{j}^{\delta }$,$\delta =-5.28$ 意味着 $\Delta Tj$ 每缩小 2×,寿命增长 $2^{5.28}\approx 39\times$。对比参照:Case D($\Delta Tj=80{}^{\circ }\mathrm{C}$,$Nf=35,200$);若 $\Delta Tj$ 降至 40°C,同 $Tm$ 下:

年循环 20,000 次下对应寿命约 68 年,远超器件物理寿命。结论:$\Delta Tj<40{}^{\circ }\mathrm{C}$ 工况下 wire bond 疲劳不是失效限制因素;FMEA 中 catastrophic failure 的 FIT 率占主导,lifetime model 退为次要工具。

paper §IV 明确器件间制造离散导致 EoL scatter $\pm 30\%$。多器件系统的工程含义超过单器件:6-phase 逆变器 6 个 MOSFET 中的第一个失效时间服从最坏 case 分布,而非均值——即使均值寿命满足要求,最坏 case 器件仍可能提前 30% 失效。正式产品应将 Coffin-Manson 的单点 $Nf$ 送入 Weibull 分布,用 shape factor $\beta \approx 3-5$(低周疲劳典型值)建模 spread,给出 B1/B10 寿命指标而非单点预测。

本页参数体系来自 IRFP340 Si MOSFET 的 wire bond 铝线疲劳——热机械应力驱动铝线根部裂纹扩展,$Rds,on$ 是最敏感的 health precursor。SiC MOSFET 主导失效机理是 BPD(basal plane dislocation)迁移:正向电流下 BPD 从 SiC 外延层向内扩展,引起 body diode $Vf$ 升高,而 $Rds,on$ 变化远小于 Si。两个直接推论:① SiC 的 $\delta$ 和 $Ea$ 与 Si 相差一倍以上,激活能不可互用;② SiC 的 health precursor 应是 body diode $Vf$ 而非 $Rds,on$。将 IRFP340 参数直接套 SiC 高压逆变器是高频陷阱,会系统性地高估 SiC 寿命。