导热硅脂的老化与泵出实战:3 年热循环后台架 η 只挪了 0.7 个点,Cpk 却先破了 1.33 线——泵出的面积平方律、来料四项验收和把 10 年压进 9 周的加速试验

2026年10月23日 0 By admin

上一篇《探头贴装质量出厂验收实战》文末我留了个预告:导热硅脂的老化与泵出——贴装合格只是起点,3 年热循环后 η 掉到多少、泵出(pump-out)的现场指纹是什么、来料批次 λ 与粘度怎么进厂验收、能不能用加速试验(温度循环+压紧维持)预测 10 年后的跟随比。这篇全部答上。

事情起于今年夏天的一次 3 年期复测。2023 年 6 月出厂的那批柜子(8 台、每台 4 个温度通道),按出厂时同一台架、同一规程抽了 4 台回测:好群体均值从 0.926 掉到 0.919,σ 还是 0.031,单台最低 0.886——所有通道全部在 0.80 门内,报告可以签”合格”。但在线数据是另一幅画面:在线斜率比中位数连续 4 个季度往下挪(0.970→0.966),Cpk 从 1.36 掉到 1.28,1.33 的线已经破了。单台看都合格,合起来看过程已经漂了。

第一反应当然是”探头批次坏了”。换上 2026 年新批探头(同供方同规格),原样贴装,半年后台架复测 η=0.9192——和 3 年前那批老头子几乎一样。第二个反应”台架漂了”:参考 Pt100 送外校,偏差 0.02K 以内;用新硅脂、新探头、同一工装现贴一个 witness coupon,η=0.9295≈规范值 0.930。台架没漂,是 DUT 变了。第三个反应”在线判据变灵了”也不成立:整柜同向、全批同向,不是个别通道的事。排除了所有外部因素,剩下的解释只有一个:3 年前那层硅脂在老化。λ 在衰减,界面在微泵出,Gc 从 53.3mW/K 走到了 47.9。

结论先行,四条:

  • 3 年组合老化(λ×0.90 叠加 0.25% 空洞)台架 η 只挪 0.7 个点(0.930→0.923),10 年也挪不到 3 个点(0.903)——老化动不了 0.80 的门。能击穿门的只有工艺/来料事故:失控脂(半衰期 4 年)第 7 年破门,忘涂脂/用错胶当天破门;
  • 但 Cpk 先破线:共同模漂移下 σ 不变,均值挪 0.7 个点就把 Cpk 从 1.36 打到 1.28——老化的全部信息在趋势里,单台合格判据的灵敏度天花板决定了它天生抓不到老化;
  • 泵出有面积平方律:L² 效应下 2mm 探头天生免疫(3 年空洞 0.25%),10mm 功率界面才是主战场(3 年 8.4%)——”界面贴得越大越保险”是反的;
  • 加速试验的全部根据是 Arrhenius:125℃ 下半衰期 121 天,10 场年=1451h≈9 周;来料用 witness coupon 把 λ 规格变成 η_w≥0.92 的功能判据,防的是失控不是老化。

工具照例全离线:grease_aging_calc.py,10 个模式,56 项自检断言,文中每个数字都能复算。链路口径与《贴装出厂验收》、《NTC 采样链路故障》、《BMS 温度门限参数表》严格一致:16S 200Ah、B3950/R25 10k/上拉 10k/3.3V/12bit、Gc0=53.33mW/K、Cs=0.35J/K、Ga_bench=4mW/K、Ga_inpack=1.5mW/K。

一、三种老化机理,三种时间尺度,都打在同一个 Gc 上

先把”老化”拆开。导热硅脂在温度场里随时间变坏,机理至少三种,时间尺度差着数量级:

机理 物理过程 驱动量 时间尺度 对 Gc 的作用方式 可逆性
油分离/干涸 硅油迁移挥发、填料沉降结块 温度(Arrhenius) 月~年 λ 缓降,整体均匀 不可逆
泵出 pump-out CTE 失配,每循环把脂挤向边缘,中心长空洞 循环数×ΔT×L² 循环驱动 有效接触面积收缩 不可逆
接触应力松弛 压力下脂蠕变、压紧力衰减 时间×温度 年 接触热阻升(L2 层) 部分可逆(复压)

三种机理有个共同落点:它们最终都改 Gc。而 124 篇立的整套判据——跟随比 η、时间常数 τ、双投影互检、自热——全是 Gc 的投影。所以老化不需要新判据、新台架,要做的只是把时间和循环数作为新自变量加进那套体系。这是本篇能薄的原因:模型是现成的,新东西只有两个函数 λ(t) 和 f(n)。

顺手把三个投影对 Gc 的灵敏度排个序,后面反复用:

  • η=Gc/(Gc+Ga) 最钝:Ga 只有 4mW/K,Gc 掉 30%,η 才掉 2.7 个点(0.930→0.903);
  • τ=Cs/(Gc+Ga) 和自热 P/(Gc+Ga) 同比例:Gc 掉 30%,两者都涨 39%(6.10→8.46s、4.75→6.58mK);
  • 结论:想抓老化别死磕 η,自热最灵,τ 是判别器。第七节的测量方法选型直接从这个排序来。

二、λ 老化:Arrhenius 与”40℃ 半衰期 20 年”

硅脂的导热系数随时间衰减(硅油分离挥发、氧化、填料沉降),工程上用半衰期参数化最省事:

λ(t, T) = λ0 · 2^(-t / t_half(T))
t_half(T) = t_half(40℃) / AF(T)
AF(T) = exp(Ea/R · (1/T40 − 1/T))        Ea=50kJ/mol(假设)

t_half(40℃)=20 年是本篇的中心假设(对”典型填充硅脂”取的工程值,正文所有 3 年/10 年数字都挂在它上面,脚本里一行可改)。Ea 取 50kJ/mol,折算出的加速因子表:

温度 AF(相对 40℃) λ 半衰期
40℃(场均温,假设) 1.00 20.0 年
60℃(柜内热点均值) 3.17 6.3 年
85℃ 11.17 1.79 年(654 天)
100℃ 21.93 333 天
125℃ 60.34 121 天

“125℃ 下 4 个月 λ 减半”——这句话就是第八节加速试验的全部根据。顺带一个场次结论:柜内热点(60℃ 均值)的半衰期只有 6.3 年,装在电芯旁和装在柜顶的探头,老化速度差 3 倍,维护周期不能一刀切。

λ 保持率(40℃ 场均温):1 年 0.966、3 年 0.901、10 年 0.707。脂的档次敏感性(η 为台架口径):

档次 t_half(40℃) λ(3y) η(3y) Cpk(3y) λ(10y) η(10y) Cpk(10y) 破门?
好脂 40 年 0.949 0.927 1.32 0.841 0.918 1.22 否
典型 20 年 0.901 0.923 1.28 0.707 0.903 1.07 否
差脂 8 年 0.771 0.911 1.15 0.420 0.848 0.47 否(Cpk 早破)
失控(油分离) 4 年 0.595 0.888 0.91 0.177 0.700 — 第 7 年

失控脂的破门时间:0.80 门 ⇔ Gc≥16mW/K ⇔ λ≥0.60 W/(m·K),λ 从 2.0 衰到 0.60 需要 t=4·log2(2.0/0.60)=6.9 年;若同时叠加 20% 空洞提前到 5.7 年。“差脂”和”失控脂”之间没有平滑过渡——bleed 一超标,半衰期按数量级掉,这就是第九节来料验收要防的东西。

C 代码里这个模型就三行,半衰期表插值即可,别在线跑 exp:

/* λ 老化预测: 半衰期表 + 线性插值(Ea=50kJ/mol, t_half(40C)=20y 折算生成) */
#define THALF40_Y   20.0f

static const struct { int8_t t_c; float half_y; } THALF_TAB[] = {
    {  40, 20.00f }, { 60, 6.31f }, { 85, 1.79f },
    { 100, 0.91f  }, { 125, 0.33f },
};

float lam_keep(float years, float t_avg_c)
{
    float half_y = thalf_interp(THALF_TAB, t_avg_c);   /* 表间线性插值 */
    return powf(2.0f, -years / half_y);
}

三、泵出:面积是平方级地重要,探头天生免疫

泵出的物理图像:两种材料热膨胀系数失配(填充硅脂 vs 金属壳体,Δα 约 250 ppm/K,假设),每个温度循环把界面里的脂往边缘挤一次,单次边缘位移 u=Δα·ΔT·L;被挤出的体积 ∝ 周长×位移×厚度 ∝ L²·ΔT,而脂总量 ∝ L²·δ——两边一除,空洞形成时间:

n_c = n0 · (δ/δ0) · (L0/L)² · (ΔT0/ΔT)
f(n) = f_inf · (1 − exp(−n/n_c))          f_inf=0.35(中心空洞上限,假设)

标定锚(工程假设,取功率器件-散热器界面的文献口径):L0=10mm、δ0=100μm、ΔT0=100K 时 n0=800 循环出现可测中心空洞。把这个锚按两种界面的实际参数展开:

界面 n_c@场ΔT20K n_c@ΔT100K f(3y) f(10y) f(1000 循环@100K)
NTC 探头 2mm/150μm 150,000 30,000 0.25% 0.84% 1.1%
功率界面 10mm/100μm 4,000 800 8.4% 21.0% 25.0%

三个读法:

第一,探头天生免疫,别为它操心。 2mm 界面的 n_c 是 15 万次场循环,3 年只走 1095 次——空洞 0.25%,对 Gc 的贡献是 −0.25%,被 λ 衰减(−10%)完全盖住。本篇标题里”η 只挪 0.7 个点”的账,九成是 λ 老化的,泵出只出了一成。

第二,功率界面才是泵出主战场,而且大得多。 10mm 界面 3 年空洞 8.4%、10 年 21.0%——配合 λ 衰减,10 年 Gc 打到 53.3×0.707×0.79=29.8mW/K,η 0.882。对贴散热器的功率器件,这是实打实要进设计评审的量。

第三,标度假设的不确定性必须靠试验兜底。 每循环位移 ∝ ΔT 是线性标度;若文献里的平方标度成立,功率界面的 f(3y) 从 8.4% 缩到 1.9%、f(10y) 从 21.0% 缩到 5.8%(探头的两个标度下都可忽略,结论不变)。一个数量级的模型不确定度,靠嘴争不出结果——第八节试验 B 就是为它设的。

设计推论只有一条,但很硬:界面面积每翻倍,泵出敏感性 ×4(L² 律)。”把探头贴大点热阻小”在老化维度上是负优化的——面积换来的 η 是零点几个点,赔进去的是 4 倍的泵出敏感性和更大的 CTE 失配位移。探头 2mm 级、功率界面 separately 评估,两类界面不要用同一套老化假设。

/* 泵出空洞比: f(n)=finf*(1-exp(-n/nc)); nc 正比于 δ/(L^2*ΔT) */
#define FINF     0.35f
#define NC_REF   800.0f     /* cycles @ L0=10mm, δ0=100um, ΔT0=100K */

float n_cycles_half(float l_mm, float d_um, float dt_k)
{
    return NC_REF * (d_um / 100.0f)
                 * (10.0f * 10.0f / (l_mm * l_mm))
                 * (100.0f / dt_k);
}

float void_frac(uint32_t cycles, float l_mm, float d_um, float dt_k)
{
    float nc = n_cycles_half(l_mm, d_um, dt_k);
    return FINF * (1.0f - expf(-(float)cycles / nc));
}

四、组合推演:η 挪 0.7 个点,Cpk 从 1.36 到 1.28 到 1.07

两个衰减函数叠起来,就是探头界面随时间的全投影表(40℃ 场均温、日循环 365 次/年 ΔT=20K,脚本 aging 模式):

t λ 比 空洞 Gc(mW/K) η_bench τ_bench(s) η_inpack 自热(mK) 斜率比中位数 Cpk Gc/门
0y 1.000 0.00% 53.33 0.930 6.10 0.973 4.75 0.970 1.36 333%
1y 0.966 0.09% 51.47 0.928 6.31 0.972 4.91 0.969 1.33 322%
3y 0.901 0.25% 47.94 0.923 6.74 0.970 5.24 0.967 1.28 300%
10y 0.707 0.84% 37.40 0.903 8.46 0.961 6.58 0.959 1.07 234%

(Cpk 一列以 124 篇的 1.36 为锚、按 (μ−LSL) 比例折算,共同模漂移口径 σ=0.0311 不变:0y 1.36 → 1y 1.33 → 3y 1.28 → 5y 1.23 → 10y 1.07。)

三个读法,一个比一个反直觉:

第一,η 是最钝的投影,而门限恰恰设在 η 上。 3 年 η 从 0.930 挪到 0.923,0.7 个点;同期 τ 涨 10%、自热涨 10%。0.80 的台架门对应 Gc≥16mW/K,只是规范值的 30%——10 年老化走到 Gc 37.4,连门的一半路程都没走完,门冗余还剩 2.3 倍。这就是为什么”每台都合格”和”Cpk 破线”同时成立:门是给事故设的(忘涂脂当天 Gc 掉 154 倍),不是给老化设的。

第二,Cpk 比任何单台判据都先报警。 共同模漂移下 σ 不变,均值 −0.7 个点,Cpk 就从 1.36 到 1.28——破 1.33 只用 3 年。而单台视角里,0.886 的最低值和 0.926 的均值差 1.3σ,混在正常批次波动里永远看不出来。老化的信息在分布的位置里,不在任何单次测量里,这是它和 123 篇那些故障(事件驱动、单通道、当场判得出)的本质区别。

第三,在线斜率比中位数是唯一零成本的老化探针,但灵敏度也就那样。 在役好群体中位数 0.970,10 年挪到 0.959——一季挪 0.003 个点。作为趋势量合格(连续 4 个季度同向就能定性),作为判据永远不够(噪声都比它大)。它的价值在回看:出事时按循环史把中位数曲线翻出来,能回答”这批柜子是不是从第 14 个季度开始整体漂的”。

最后是 124 篇那条标定传递链的老化终点:台架好群体分布(η=0.926±0.0311)→ 台架门 0.80 → 现场判据线 0.70/0.30。老化不改这条链的结构,只是把整条分布沿时间轴平移——链上的每个值都要按”出厂 + N 年”重新理解,这也是维护周期的算法依据(第十节)。

五、现场指纹五条与一个判别器:τ 动不动

老化在监控上长什么样?五条指纹,够用了:

  1. 整柜同向:全柜 16 个通道一起降,降幅几乎一样(±0.005 以内)——同批脂、同循环史、同环境。反例是 123 篇的单通道故障:随机单点、无关联;
  2. 双投影同向:η 降 + τ 涨,且 Gc(η)=ηGa/(1−η) 与 Gc(τ)=Cs/τ−Ga 两个反解同步下降、相对差仍在 25% 内。只有 η 动而 τ 不动,就不是 Gc 的事;
  3. 自热抬升:同电流下探头自热 4.75→5.24mK(3 年),双偏置法下指示温差 169.7→187.3mK(第七节)——最灵敏的投影;
  4. 与累计循环积分相关,而非日历:η(t) 跟 Σ|ΔT|(温差绝对值积分)走,夏季末最差;入冬不恢复(老化不可逆),但在线测量值有季节调制(环境温度改变信号幅度,不改老化趋势本身);
  5. 换批台阶:供方/批号变化 → 整线台阶式变化,新批上台架复测恢复——我这次就是靠指纹五排除了”探头批次坏了”。

判别器只有一个,但一刀切开三类病:τ 动不动。

症状组合 τ 自热 结论 去哪查
η 整柜齐降 齐涨 齐涨 Gc 真降(老化/来料) 循环史 + 来料批次
η 整柜齐降 不动 不动 公共路径(上拉/基准/软件系数) 123 篇口径:全柜等量平移
η 齐降 + 双投影差>25% 矛盾 矛盾 分布式接触/内部缺陷 拆开分析,不是老化
η 单通道降 单通道涨 单通道涨 123 篇故障树 按 123 篇 triage

公共路径为什么 τ 不动:上拉电阻 ±1%、基准漂移、软件系数错,改的是”raw→温度”的换算,Gc 一点没动——η(相对电芯的增益)和 τ(热时间常数)都该原地不动,动的只有读数整体平移。把公共路径的等量平移当成老化,就会去换一批硅脂,然后发现新脂上去照旧——这是本篇最贵的一条指纹教训,我拿 4 台柜子和半年时间换的。

六、自热是最灵敏的投影:双偏置法把 4.7mK 的自热放大成 169.7mK 判别信号

η 钝、τ 中、自热灵——但常态自热只有 4.75mK(0.165mA 激励、25℃),读数噪声(±0.5mK 级)占比 10%,3 年 10% 的抬升直接淹掉。解法是经典的双偏置法:同一探头分别在两个激励电流下读数,自热造成的指示温差与 (I²−I1²) 成正比,总热导一把反解:

G_tot = (I2² − I1²)·R25 / ΔT_ind

取 I1=0.165mA(常态 3.3V/20k)、I2=1mA(诊断模式),功耗差 9.73mW,指示温差被放大到百 mK 级:

t Gc(mW/K) 自热@0.165mA(mK) 自热@1mA(mK) 双偏置指示温差(mK) 提取 G_tot(mW/K)
0y 53.33 4.75 174.4 169.7 57.3
1y 51.47 4.91 180.3 175.4 55.5
3y 47.94 5.24 192.5 187.3 51.9
10y 37.40 6.58 241.6 235.0 41.4

判别力:3 年 vs 新鲜差 17.6mK、10 年 vs 新鲜差 65.3mK——0.5mK 的读数分辨率下分别是 35σ 和 130σ。3 年的老化,一次双偏置测量就能可信判定,比任何台架复测都便宜(台架要把柜子拆了拉回去,双偏置在柜子上就能跑)。

三个使用限制说清楚,不然又要被当成银弹:

  • 1mA 只能诊断模式用。常态开着 1mA,读数本身就被自热抬高 174mK,温度门限整体偏 0.17K——0℃ 禁充门能偏出 0.2 度去。测完必须回 0.165mA;
  • 必须同环境比。提取出的 G_tot=Gc+Ga,Ga 由引线路径决定(台架 4mW/K、在役 1.5mW/K),台架值和柜内值差 2.5mW/K,跨环境比绝对值会得出错误结论——比趋势、比同批同位置;
  • 判不出”对不对”。和 124 篇的结论一样:双偏置测”贴得好不好、老没老”,测不了”贴得对不对”(位置错误),后者靠定位工装。
/* 双偏置自热法: G_tot = (I2^2-I1^2)*R25 / ΔT_ind */
#define R25       10000.0f
#define I1_MA     0.165f      /* 常态激励 3.3V/20k */
#define I2_MA     1.0f        /* 诊断激励 */

float g_tot_from_dual_bias(float t_ind1, float t_ind2)
{
    float dp = (I2_MA * I2_MA - I1_MA * I1_MA) * 1e-6f * R25;  /* W, ≈9.73mW */
    float dt = t_ind2 - t_ind1;                                /* K */
    if (dt < 1e-4f) return 0.0f;          /* 判别力不足, 不判 */
    return dp / dt;                       /* W/K */
}

七、加速试验:把 10 年压进 9 周

老化模型最大的敌人是”等不起”:10 场年 = 87,600 小时 = 4.2 年的连续观测,比产品生命周期还长。加速试验的全部根据就是第二节那张 Arrhenius 表——温度换时间,循环数换循环数。三个试验,各管一种机理:

试验 A:λ 老化(烘箱)。witness coupon 与产线同几何同工艺(4mm²、0.15mm、同工装压力、同批表面处理),进烘箱,定点取出上台架测 η(124 篇同一规程):

条件 等效场年 预测 λ 比 时长 用途
500h @ 125℃ 3.4 0.888 3 周 快速筛批次
1000h @ 100℃ 2.5 0.917 6 周 AF 验证点(对预测)
1000h @ 125℃ 6.9 0.788 6 周 强应力档
1451h @ 125℃ 10.0 0.707 9 周 10 年全预测点

判据:1451h@125℃ 后 η ≥ 0.85(= 门 0.80 以上留 5 个点余量)。这一条能把三档脂分开:典型脂 0.903 过、差脂 0.848 压线毙、失控脂 0.700 直接毙。1000h@100℃ 那个点不判合格,用来验证 AF 本身——实测 λ 比若偏离预测 0.917 超过 ±10%,Arrhenius 假设在该脂上不成立,外推作废重标。

试验 B:泵出循环(温箱)。DUT 带实际贴装进温箱,−40↔125℃ 或 0↔100℃(ΔT=100~165K),dwell 10min,斜坡 ≤10K/min,恒力工装保压(0.5MPa,假设),每 100 循环取出、回 25℃ 稳态后测 η。线性标度下 1 次 ΔT=100K 循环 ≡ 5 次场循环:

循环数@ΔT100K 场等效 f(探头) f(大界面 10mm)
100 1.4 年 0.12% 4.1%
500 6.8 年 0.58% 16.3%
1000 13.7 年 1.15% 25.0%

判据:探头界面 Δη≤0.02;功率界面 ΔRth≤15%(文献口径 10~20% 的下沿取)。试验 B 是第三节那个”标度不确定”的粉碎机:线性标度下 1000 循环 = 13.7 场年,平方标度下 = 68 场年,实测落在哪,模型就用哪个——一次试验定标,后面全线的预测都有了根。

试验 C:应力松弛(保压)。85℃ 恒位移保压,记录压力衰减(p(t)=p0·exp(−t/τc),τc=2000h 假设):500h 保压 78%、1000h 61%、2000h 37%。判据:1000h 保压 ≥60%。这一项只管压力依赖型结构(纯压紧的功率界面);RTV/卡扣承担机械的探头结构可以不办,但办了能提前抓住”工装夹持力设计不足”这类隐性缺陷。

样本与对照:每条件 n≥8(同 124 篇的 Cpk 口径),每轮必须带金批参考样(存续的合格批),用来把”试验系统漂移”和”样品真退化”分开——没有参考样,试验 A 的第 1000 小时数据谁也不敢签。成本账:A 全预测 9 周、B 4~5 周、C 2 周,可并行;一次完整跑下来管一个供方 10 年,比售后换探头便宜两个数量级。

八、来料验收:四项指标加一个 witness coupon

来料验收要防的不是老化(老化防不住,只能预测),是失控——半衰期从 20 年掉到 4 年的那种批次事故,它在第 7 年破门,之前毫无症状。四项指标:

项 方法 判据 防的失效模式
λ 导热系数 型式:ASTM D5470 热板/激光闪射(供方报告+年度第三方);批次:witness coupon 台架阶跃 λ≥1.725 W/(m·K) ⇔ η_w≥0.920 劣质/干涸批次
粘度/触变 旋转流变仪 25℃ 对单 供方窗口内(典型 300~800 Pa·s@1s⁻¹,假设) 泵出快/涂不匀
油分离 bleed 100℃×24h 烘箱或离心 2000rpm×24h 质量损失 ≤2% 干涸快(半衰期骤降)
CTE TMA/膨胀计 ≤250 ppm/K(假设) 泵出 n_c(第三节)

外加两项管理项:垂直涂片 80℃×24h 不流淌(立式安装防流失);批号追溯+先进先出,保质期 12~24 个月(假设)。

witness coupon 是这套体系里最便宜的一块:把 λ 这个”材料参数”变成 η_w 这个”功能参数”,用产线上已有的台架阶测就能放行,不需要热板仪。几何必须和产线严格一致(4mm²、0.15mm、同压力、同表面处理),否则 η_w 和 λ 的映射关系要重标:

λ (W/mK) Gc_w(mW/K) η_w 判读
2.00(规范) 53.3 0.930 基准
1.80(−10%) 48.0 0.923 合格
1.725(−14%) 46.0 0.920 批次门
1.60(−20%) 42.7 0.914 低于门,毙
1.00(劣质) 26.7 0.870 witness 直接毙

判据的来历:η_w=0.92 反解 Gc_w=46mW/K,对应 λ=1.725,即规范值 −15%。为什么门不开在 −20%:λ 掉 20% 时 η_w 只掉 1.6 个点(0.930→0.914),台架分辨率 ±0.004 完全判得出,但现场在役 10 年 λ 还要再走 −29%,来料就把冗余用掉一半,10 年终点就贴门了。每批 witness coupon 的 η_w 落盘,它同时是供方趋势数据和事故追溯数据——”这批 bleed 超标”这句话,要有 η_w 曲线作证。

九、设计侧缓解:五条,第一条反直觉

  1. 界面面积别贪大(L² 律):探头 bead 2mm 级天然免疫泵出;功率界面面积每翻倍泵出敏感性 ×4。”贴大点热阻小”在老化维度是负优化,两类界面分开评估;
  2. 脂层 0.1~0.2mm 甜点:n_c ∝ δ(厚层耐泵出)但 Gc ∝ 1/δ(薄层热阻小),权衡表:
脂厚 δ Gc(mW/K) η_bench n_c@ΔT100K 评价
0.10mm 80.0 0.952 20,000 热阻最小,泵出敏感
0.15mm 53.3 0.930 30,000 规范值,甜点
0.20mm 40.0 0.909 40,000 硅脂量 +25% 档
0.30mm 26.7 0.870 60,000 过量涂覆,η 已贴门

探头应用里泵出冗余有 15 万次场循环,η 该赢这场权衡(走 0.1~0.15mm);功率界面 n_c 只有几千次,反过来该偏厚(0.2mm 起)。同一支硅脂,两种界面,两个最优厚度——这就是”数值不可互搬”在厚度维度的翻版;

  1. 机械与热分工:RTV/卡扣承担机械应力,硅脂只管热界面。压力依赖型结构最怕应力松弛(试验 C),探头上这是个便宜的结构性改进;
  2. 别在固件里补偿老化:把 Gc(t) 折进采样链路做温度补偿,等于把不可观测的漂移变成隐性错误——补偿系数本身也会错,而且错了没有任何判据能发现。老化要监测,不要掩盖;
  3. 趋势量落盘:Cpk 按批、在线斜率比中位数按季、双偏置 G_tot 按维护周期,全部按累计循环史归档。老化的全部信息在趋势里,不落盘等于没有。

十、排障 triage:看到”齐降”先别换探头

现场/产线速查(脚本 triage 模式):

症状 τ 自热 分布特征 结论 动作
η 整柜齐降 齐涨 齐涨 同批脂+同循环史 真老化(λ+泵出) 算循环史,落盘趋势,排维护周期
η 整柜齐降 不动 不动 等量平移 公共路径(上拉/基准/系数) 查电路与软件系数(123 篇)
η 台阶式齐降 齐涨 齐涨 换供方/换批号 来料批次事故 冻结批次,查 bleed/λ 报告
η 单通道降 单通道涨 单通道涨 随机单点 123 篇故障树 按 123 篇 triage
η 降+双投影差>25% 矛盾 矛盾 随机 分布式接触/内部缺陷 拆开分析,不是老化
台架 η<0.80 或 τ>60s >60s >10mK 整批双峰 工艺/来料事故 查挤出机/介质,100% 重扫
在线斜率比中位数季漂 — — 全线缓慢 老化趋势量 Cpk 趋势+循环史回看

第一行和最后一行的动作都包含”不换件”——门内的老化不是缺陷,是寿命在按模型走路;该做的是把预测终点和门限的余量算出来,写进维护周期(10 年预测 η=0.903,门内,到期轮换;若某供方预测 0.85 以下,第 7 年轮换)。把趋势当时刻、把寿命当故障,是比不监控更贵的错误。

十一、跟前几篇的口径对齐与边界

四篇的时间轴分工,一句话各表:

篇 管什么 驱动量 判据性质
122 门限参数表 参数消费侧 事件(越限) 规则
123 NTC 故障判据 在役数据质量 事件(故障) 在线三层
124 贴装出厂验收 出生质量 工艺(贴装) 台架绝对
本篇(125) 时间维度 时间/循环 趋势量+加速试验

继承锚点(全部与 124 篇一致,未重算):Gc0=53.33mW/K、Cs=0.35J/K、Ga_bench=4mW/K、Ga_inpack=1.5mW/K、台架门 0.80/60s、好群体 0.926±0.0311、Cpk 1.36、斜率比中位数 0.97、电链 B3950/R25 10k/3.3V/12bit、P_ntc(25℃)=0.272mW。

本篇新增假设(全部集中在脚本头部,一行可改,正文均已标注):t_half(40℃)=20y、Ea=50kJ/mol、场均温 40℃、日循环 365 次 ΔT=20K、泵出 n0=800@10mm/100μm/100K、f_inf=0.35、Δα=250ppm/K、τc(85℃)=2000h、λ 批次门 1.725、bleed≤2%、CTE≤250ppm/K、保压 0.5MPa。

口径差声明:λ 衰减是”半衰期+Arrhenius”的经验参数化,不是机理模型;泵出 n_c 的标定来自功率器件文献口径的工程假设;场循环参数(365 次/年、ΔT=20K)是柜内 BMS 舱的典型值。换产品、换供方、换柜体结构,这套数要重标——重标的方法就是第七节三个试验,别直接搬。

十二、结语与预告

收个尾。硅脂老化的数字化,本质是承认一件事:贴装界面不是一个常数,是一条随时间走的曲线。这条曲线可测(η/τ/自热全是探针)、可预测(Arrhenius+循环计数)、可验收(来料四项+加速试验三件套):

  • 一个落点:三种老化机理都改 Gc,判据体系全部复用,只加 λ(t) 和 f(n) 两个函数;
  • 一个排序:η 最钝守门、τ 是判别器(动不动)、自热最灵(双偏置法放大 36 倍);
  • 一条铁律:老化是趋势量不是事故量——η 挪 0.7 个点而 σ 有 3.1 个点,单台合格判据天生抓不到它;Cpk、在线中位数、双偏置 G_tot 才是它的探针;
  • 一个根据:加速试验的全部底气是 Arrhenius——125℃ 半衰期 121 天,10 场年=1451h≈9 周;来料验收防的不是老化,是半衰期 4 年的失控脂(第 7 年破门)。

工具全离线:grease_aging_calc.py,10 个模式(base/lam/pump/aging/cpk/accel/witness/selfheat/triage/table),56 项自检断言,含 124 篇基线锚点(Gc0=53.33、η=0.9302、τ=6.105s、自热 4.75mK)、Arrhenius 锚点(AF 100℃=21.93、125℃=60.34)、泵出对称性锚点(1000 循环@100K 对探头的损伤 ≡ 5000 场循环)和失控脂破门锚点(6.9y/5.7y)。参数改成本家产品的口径(介质 λ、半衰期、循环史、界面尺寸)就能直接用。

下期预告:多通道 NTC 巡采的 MUX 时序与同步采样——123/124 两篇都把它挂在备选位,这次转正:多路复用一个 ADC 时,MUX 的电荷注入与串扰、采样保持窗口、通道间 skew、巡检vs同步采样的判据表,以及”第 5 通道总是偏低 0.3 度”这类现场病的定位。备选仍是那桩悬案:09-26《ESD/浪涌防护实战》的数值错误回修(12.5kV 应为 12.5V,连挂多期,编辑非新写)。

自查清单(来料与老化管理,12 项):

  • ☐ 来料 λ 有型式报告 + 批次 witness coupon 双重把关
  • ☐ witness coupon 与产线同几何同工艺,台架阶跃提取 η_w,门 η_w≥0.920
  • ☐ 油分离 bleed≤2%,每批验,曲线落盘
  • ☐ 粘度/流淌在供方窗口内,批批对单
  • ☐ CTE 有数据(泵出模型输入),未填充硅脂不上功率大界面
  • ☐ 加速试验 A/B/C 各一样:A 验 λ 老化、B 验泵出、C 验应力松弛
  • ☐ 每轮试验带金批参考样,n≥8
  • ☐ 试验 A 的 1000h@100℃ 点用来验证 AF,偏离 ±10% 重标
  • ☐ Cpk 按批落盘,在线斜率比中位数按季落盘,按循环史可回看
  • ☐ 维护周期按 10 年预测 η 定,不按”感觉”定
  • ☐ 看到整柜齐降先查 τ 动不动,别先换探头
  • ☐ 固件不做老化补偿,只做监测;1mA 双偏置仅诊断模式用,测完回 0.165mA

以上判据与数值基于文中给定的模型口径(16S 200Ah、Gc0=53.33mW/K、Cs=0.35J/K、Ga_bench=4mW/K、Ga_inpack=1.5mW/K、接触面积 4mm²、t_half(40℃)=20y、Ea=50kJ/mol、泵出 n0=800@10mm/100μm/100K、f_inf=0.35、场循环 365 次/年 ΔT=20K)与典型工程假设推算,不同厂家探头的封装热阻、硅脂物性批次、柜体循环史会改变具体数值;落地前请以本方产品的台架复测、加速试验与在役标定结果为准。