层压母线设计与CAE仿真
理论与物理
概述 — 什么是层压母排
老师，层压母排和普通的铜板有什么区别？就是在逆变器里看到的那些很多薄板叠在一起的东西吧？
问得好。层压母排是将P层（正极）和N层（负极）的铜板用薄绝缘膜（聚酰亚胺或PET，厚度约50~200µm）夹住并层压而成的。关键在于电流以相反方向流动。反向电流产生的磁场相互抵消，从而可以将回路电感降低到数nH级别。
数nH？那普通的母排大概是多少呢？
如果用实心铜母排布线，根据长度不同，大约在50~200nH左右。在电动汽车逆变器中，母线电压为600V，开关电流为数百A，浪涌电压由 $V_{\text{surge}} = L \cdot di/dt$ 决定。例如，若 $L = 50\,\text{nH}$，$di/dt = 5\,\text{kA/}\mu\text{s}$，则 $V_{\text{surge}} = 250\,\text{V}$。这在800V系统的SiC逆变器中叠加起来，可能会超过元件的绝缘耐压。因此，必须通过层压结构将 $L$ 抑制在5~10nH以下。
原来如此，开关速度越快，电感的影响就越大啊。但是为什么反向电流会使磁场抵消呢？
回想一下安培定律。流过电流 $I$ 的导体周围会产生 $\mathbf{H} = I/(2\pi r)$ 的磁场。如果在P层流过 $+I$，在紧邻的N层流过 $-I$，则两个磁场在导体外部几乎相互抵消。导体间距 $d$ 越小（即绝缘层越薄），抵消效果越接近完美。这就是层压母排的基本原理。从数学上讲，回路所包围的面积 $A$ 变小，而 $L \propto A$，因此电感急剧减小。
支配方程 — 部分电感与电流分布
请告诉我设计中使用的公式。电感是怎么计算的？
母排的电感计算常用部分电感法（Partial Element Equivalent Circuit, PEEC）。将导体分割成微小单元，单元 $i$ 和单元 $j$ 之间的相互部分电感 $L_{p,ij}$ 用诺依曼公式计算：
其中 $a_i, a_j$ 是各单元的截面积，$V_i, V_j$ 是各单元的体积。自部分电感（$i = j$）可以用相同公式计算，但需要适当处理 $|\mathbf{r}_i - \mathbf{r}_j| \to 0$ 的奇异性。
那和回路电感的关系是怎样的？
由P层和N层组成的母排的回路电感为：
$L_{p,\text{P}}$ 是P层的自部分电感，$L_{p,\text{N}}$ 是N层的自部分电感，$M_{p,\text{PN}}$ 是P-N层间的相互部分电感。在层压结构中，$M_{p,\text{PN}} \approx L_{p,\text{P}} \approx L_{p,\text{N}}$，因此 $L_{\text{loop}} \to 0$。这就是“磁场抵消”的定量表达。
太厉害了！就是说通过减法几乎变成零了。那决定电流分布的方程呢？
导体内的电流密度分布，从准静态麦克斯韦方程组出发。导体内的涡流方程为：
$\mathbf{A}$ 是磁矢势，$\sigma$ 是电导率，$\mathbf{J}_s$ 是源电流密度。在频域中，将 $\partial/\partial t \to j\omega$ 替换：
通过3D FEM求解此方程，即可得到开关频率下的电流密度分布 $\mathbf{J}(\mathbf{r}, \omega)$。
趋肤效应与邻近效应
趋肤效应在母排中也会成为问题吗？我以为板子薄就没关系呢。
想得太简单了。铜导体的趋肤深度 $\delta$ 为：
对于铜（$\sigma = 5.8 \times 10^7\,\text{S/m}$）：
| 频率 $f$ | 趋肤深度 $\delta$ | 备注 |
|---|---|---|
| 1 kHz | 2.09 mm | Si-IGBT基波程度 |
| 10 kHz | 0.66 mm | SiC MOSFET开关频率带 |
| 20 kHz | 0.47 mm | 典型的EV逆变器 |
| 100 kHz | 0.21 mm | GaN器件·高频DC-DC |
| 1 MHz | 0.066 mm | GaN高频变换器 |
20kHz时是0.47mm！母排的铜板厚度大多是1~3mm吧。那内部几乎没电流流过？
没错。导体厚度 $t$ 与趋肤深度 $\delta$ 的比值 $t/\delta$ 很重要，当 $t/\delta > 3$ 时，中心部分几乎没有电流。而且邻近效应（proximity effect）会雪上加霜。当P层和N层靠近时，对方导体产生的磁场会使电流进一步集中在导体表面。层压母排的结构使得邻近效应强烈作用，有时有效电阻会达到直流电阻的2~5倍。
交流电阻增大系数 $F_R$，对于矩形截面导体，近似为：
这里 $m$ 是层数（P-N-P-N时 $m=2$）。$m$ 越大，邻近效应越显著，因此不能随意增加层数。
电流分布均匀性评估
不只是电感，电流是否均匀流动也很重要吗？
非常重要。在功率模块并联配置中，如果电流不能平均分配到各模块，过大的电流会集中在部分模块上，导致热失控。使用不均匀系数 $\eta$ 作为电流均匀性的评估指标：
理想情况是 $\eta = 1.0$。实用上以 $\eta < 1.3$ 为设计目标。若 $\eta > 1.5$，则意味着特定模块流过超过50%的过电流，会对可靠性产生严重影响。
电流偏流的原因是什么？
主要有三个原因。(1) 端子位置的非对称性 — 从输入端子到各模块的阻抗不同导致电流偏流。(2) 孔、槽、凹凸引起的电流路径紊乱。(3) 趋肤效应和邻近效应导致的高频电流局部集中。通过3D FEM可视化并优化所有这些因素，是母排设计的核心。
各物理量的含义与单位
| 变量 | SI单位 | 母排设计中的典型值 |
|---|---|---|
| 回路电感 $L_{\text{loop}}$ | H | 3~50 nH（层压：3~10 nH，实心：50~200 nH） |
| 电流密度 $J$ | A/m² | DC：3~10 A/mm²，脉冲峰值：50~200 A/mm² |
| 电导率 $\sigma$（铜，20°C） | S/m | $5.8 \times 10^7$（100°C时降至约$4.3 \times 10^7$） |
| 趋肤深度 $\delta$ | m | 0.066~2.09 mm（1 kHz~1 MHz） |
| 绝缘层厚度 $d$ | m | 50~200 µm（聚酰亚胺/PET） |
| 导体厚度 $t$ | m | 0.5~3 mm |
假设条件与适用范围
- 准静态近似：忽略位移电流（$\omega\varepsilon \ll \sigma$）。当母排尺寸远小于波长时有效。100kHz、尺寸30cm时波长3000m，基本没问题。
- 线性材料：铜、铝为非磁性（$\mu_r = 1$）且线性。若包含铁制母排夹，则需要非线性B-H曲线。
- 温度依赖性：铜的电导率随温度下降约0.4%/°C。大电流时需要电磁-热耦合分析。
- PEEC法的局限：存在磁性材料或屏蔽时，仅用PEEC法不够充分。需要与FEM结合或使用混合方法。
“IGBT莫名损坏”的元凶是母排
某电动汽车制造商开发初期，原型逆变器中频繁发生IGBT模块原因不明的损坏。测量开关波形发现，关断时观测到超过800V的浪涌。对于额定1200V的IGBT，800V的浪涌本不应致命，但实际上，IGBT内部布线电感与外部母排电感之和，使得实际施加在芯片上的浪涌超过了1000V。将母排从实心结构（约80nH）改为层压结构（约8nH）后，浪涌被抑制在200V以下，问题完全解决。$V = L \cdot di/dt$ 这个简单的公式，左右了数亿日元产品可靠性的实际案例。
数值解法与实现
PEEC法电路参数提取
计算母排的电感，应该用FEM还是PEEC法？
根据用途区分使用。PEEC法只需对导体进行网格划分，不需要空气域，Ansys Q3D Extractor就采用了这种方法。在将母排的RLC参数传递给电路仿真器（如SPICE）的工作流程中，它非常有效。另一方面，3D FEM擅长包含磁性体和屏蔽影响的电流密度分布可视化。
PEEC法的离散化是怎么做的？
将导体分割为长方体单元，每个单元用电阻 $R$ 和部分电感 $L_p$ 的串联电路表示。包含单元间的相互电感 $M_p$，整体的电路方程为：
$\mathbf{R}$ 是电阻矩阵（对角），$\mathbf{L}_p$ 是部分电感矩阵（稠密矩阵），$\mathbf{I}$ 是各单元的电流向量。$\mathbf{L}_p$ 是稠密矩阵，这是PEEC法计算成本上的瓶颈，对于单元数 $N$ 需要 $O(N^2)$ 的内存。大规模问题会使用多极展开（FMM）进行加速。
3D FEM公式化
FEM的话是用边单元吧？母排也一样吗？
没错。电磁场FEM中，通常使用Nédélec边单元（棱边单元）对磁矢势 $\mathbf{A}$ 进行离散化。边单元自动保证切向分量的连续性，并能排除节点单元中出现的伪模式（非物理解）。
弱形式为：
$\mathbf{N}_i$ 是边形状函数，$\Omega_c$ 是导体区域。将其离散化得到：
$\mathbf{K}$ 是旋度-旋度刚度矩阵，$\mathbf{C}$ 是与电导率相关的质量矩阵，$\mathbf{a}$ 是边自由度向量。对于母排级别的问题规模（10万~数百万自由度），用直接法（如MUMPS）求解比较稳定。
网格策略 — 应对趋肤深度
网格方面需要注意什么？和普通的结构分析不一样吗？
最重要的是在趋肤深度内至少放入3层以上的单元。如果是20kHz的分析，$\delta = 0.47\,\text{mm}$，那么从导体表面到0.47mm的范围内，需要单元尺寸在0.15mm以下的网格。相反，导体中心部分几乎没有电流，网格可以粗一些。
| 网格要求 | 推荐值 | 理由 |
|---|---|---|
| 趋肤深度内的层数 | 3~5层 | 准确解析电流密度梯度 |
| 绝缘层的单元数 | 1~2层 | 厚度方向的电位分布（仅在需要时建模） |
| 面内方向的单元尺寸 | 导体宽度的1/20~1/50 | 捕捉端部的电流集中 |
| 空气域的范围 | 导体尺寸的3~5倍 | 充分表现磁场的衰减 |
| 推荐单元类型 | 二阶六面体边单元 | 精度与效率的平衡 |
母排的铜板又薄又宽，用六面体网格好像很容易划分。
你注意到了关键点。母排是薄板形状，因此与扫描网格（六面体单元拉伸）非常匹配。在厚度方向施加偏置网格（表面细、中央粗）会很高效。不过，螺丝孔和槽口周围常常会与四面体网格混合使用。
频域与时域的选择
应该用频域还是时域？
对于稳态交流阻抗或电感提取，频域的效率压倒性地高。每个频率只需解一次联立方程。在多个频率下扫描即可得到阻抗的频率特性。
另一方面，要直接评估开关瞬态现象（开通/关断时的 $di/dt$ 浪涌），则需要时域分析。但时间步长必须小于最高频率成分的1/20，因此计算成本会高出几个数量级。
在实际工作中，最有效的方法是两步走：首先在频域提取RLC参数，然后将其传递给SPICE等电路仿真器进行开关波形仿真。
频域与时域的区分使用 — 直观理解
频域好比“用体温计量体温” — 可以一下子了解稳态特性。时域好比“做心电图” — 能看到瞬息万变的瞬态行为，但需要长时间的记录。在母排设计中，首先用体温计（频域提取Lp）把握整体健康状况，只对可能有问题的部位做心电图（时域确认浪涌波形），是高效的推进方法。