耦合问题的降阶模型(ROM)
理论:用少数时间常数表示响应
概述
降阶模型会降低精度,为什么还用于耦合分析?
因为耦合中“速度”与精度同样重要。例如逆变器电路仿真为了追踪开关动作,要用细小时间步计算数十万步,每一步都需要求半导体温度,并把随温度变化的损耗返回电路。直接放入详细的有限元热模型,计算根本算不完。所以要建只用少数变量再现热模型输入(发热)与输出(温度)关系的降阶模型。传热是线性的,响应通常由少数时间常数决定,做得好的话几个变量就能得到与原模型几乎相同的结果。关键是要确认它在什么范围内正确。
用福斯特网络降阶
原模型是在1 cm²面积上叠加硅0.3 mm、焊料0.1 mm、铜3 mm、导热脂0.1 mm、铝20 mm,底面以换热系数2,000 W/m²K冷却的一维热模型(200自由度,热阻6.45 K/W,最慢时间常数32.1秒)。对其阶跃响应拟合3项。
与随温度增加的损耗耦合
P₀ = 10 W,环境25℃。温度每升高1 K损耗增加0.6%(假设导通电阻的温度依赖)。
福斯特与考尔
用电阻、电容等组合的电路表示给定响应的方法,诞生于1920年代的电话电路研究。福斯特于1924年、考尔于1926年分别提出了不同形式的网络。福斯特网络把响应表示为指数函数之和,拟合简单;考尔网络为梯形,容易与层的排列对应。百年前的电话理论,如今每天都作为功率半导体热降阶模型使用。
算例
例1:拟合3级福斯特网络
| 级 | 热阻 | 时间常数 | 主要表示 |
|---|---|---|---|
| 1 | 0.061 K/W | 0.0029秒 | 芯片与焊料 |
| 2 | 0.446 K/W | 0.68秒 | 铜板与导热脂 |
| 3 | 5.943 K/W | 32.1秒 | 散热器与冷却 |
| 时间 | 原模型 | 降阶模型 |
|---|---|---|
| 0.01秒 | 0.060 K/W | 0.067 K/W |
| 0.1秒 | 0.146 K/W | 0.140 K/W |
| 1秒 | 0.580 K/W | 0.587 K/W |
| 10秒 | 2.101 K/W | 2.098 K/W |
| 100秒 | 6.185 K/W | 6.186 K/W |
整体拟合误差最大0.008 K/W(热阻的0.1%),但在0.01秒时相对偏差超过一成。
例2:电-热耦合计算
| 工况 | 原模型(200自由度) | 降阶模型(3状态) |
|---|---|---|
| 通电5秒·停5秒共300秒,最高温度 | 71.79℃ | 71.76℃ |
| 连续通电稳态温度(有耦合) | 130.2℃ | 130.2℃ |
| 连续通电稳态温度(损耗恒定) | 89.5℃ | 89.5℃ |
计入耦合后稳态温度从89.5℃升到130.2℃。损耗每K只增0.6%,为什么差这么多?
因为存在正反馈:温度升高损耗增加,损耗增加温度进一步升高。损耗恒定时温升为热阻×损耗=64.5 K,计入耦合则为64.5/(1 − 0.006×64.5) = 105 K。分母中的0.39越接近1温度上升越急,达到1时温度不再稳定持续上升(热失控)。所以电-热耦合中只用损耗恒定的热计算会得到偏危险的估计。降阶模型正是在这种场合发挥作用,在电路仿真中每一步交换温度和损耗,计算时间几乎不增加。但降阶模型以线性导热为前提,材料热导率随温度变化大或冷却条件改变时,要用原模型确认,或按条件重新建立降阶模型。
建模流程
- 用详细热模型计算各发热源的阶跃响应。
- 用指数函数之和拟合,确认时间常数范围覆盖使用的时间尺度。
- 有多个发热源时,也计入相互影响(互热阻)。
- 嵌入电路仿真,交换温度和损耗。
- 在代表性条件下确认与原模型的差别。
“在短脉冲上失准的预测”
某电机驱动装置中,降阶模型预测的芯片温度在电机启动时几毫秒的大电流下比实测低近10 K。降阶模型是按秒级响应建立的,最短时间常数比几毫秒长,表示不了芯片内的快速升温。增加2个短时间常数项重建后,启动时温度也与实测一致了。
常见错误
错误与对策
| 错误 | 后果 | 对策 |
|---|---|---|
| 只用长时间拟合 | 短脉冲温度出错 | 在使用时间范围内拟合 |
| 损耗取常数 | 温度低估 | 与温度耦合 |
| 只对一个发热源降阶 | 漏掉相邻元件影响 | 计入互热阻 |
| 冷却条件改变后继续沿用 | 时间常数不符 | 按条件重建 |
| 使用福斯特节点温度 | 没有物理意义 | 只用输入输出关系 |
我想了解相关内容。
相关页面参考模型降阶方法(POD)、数字孪生、耦合分析的收敛、数据驱动的多物理场分析和联合仿真耦合。