微带线路的电磁场分析
理论与物理
什么是微带线？
微带线说白了就是印刷电路板上的布线图案吧？为什么需要电磁分析呢？
问得好。对于低频数字电路，布线只是简单的“连接”而已。但到了GHz频段，波长与基板尺寸相当，布线会表现为传输线。如果阻抗不匹配，信号会反射导致波形失真。
具体有多严格呢？
例如USB 3.2或PCIe Gen5要求90 $\Omega$ 差分阻抗控制在 $\pm 5\%$ 以内。超出此范围会导致眼图闭合、比特误码率增加。5G基站射频前端有时需要将50 $\Omega$ 控制在 $\pm 2\%$ 以内。因此电磁场仿真不可或缺。
微带线是一种平面型传输线结构，在介质基板的一面布置导体带，另一面布置接地平面。它广泛应用于PCB（印刷电路板）和MMIC（单片微波集成电路）。
由于其结构是开放的（顶部暴露在空气中），无法支持纯粹的TEM模式，而是传播准TEM模式（quasi-TEM）。低频时TEM近似成立，但随着频率升高，TE/TM分量不可忽略，会产生色散。
特性阻抗的理论公式
那么，那个阻抗是怎么计算的呢？需要直接使用CAD工具吗？
首先有可用于手算的近似公式。Hammerstad-Jensen（1980）公式在实际工作中最常用，根据 $w/h > 1$ 和 $w/h \leq 1$ 分情况讨论。在 $\varepsilon_r \leq 16$、$0.05 \leq w/h \leq 20$ 范围内，精度约为1%以内。
窄线宽（$w/h \leq 1$）的情况：
宽线宽（$w/h > 1$）的情况：
此外，包含导体厚度 $t$ 的IPC-2141标准简易公式也被广泛使用：
在FR4基板上要实现50 $\Omega$，具体需要多宽的线宽呢？
对于FR4（$\varepsilon_r \approx 4.3$），基板厚度 $h = 0.2\,\text{mm}$，铜箔厚度 $t = 35\,\mu\text{m}$ 时，$w \approx 0.36\,\text{mm}$ 即可得到50 $\Omega$。因为 $w/h \approx 1.8$，所以使用宽线宽公式。实际工作中，通常用2D截面求解器（例如SIwave）进行验证。
有效介电常数与色散
公式中出现的 $\varepsilon_{eff}$ 是什么？和基板的介电常数 $\varepsilon_r$ 不一样吗？
因为微带线的电场跨越了介质和空气，电磁波“感受到”的介电常数是 $\varepsilon_r$ 和1（空气）之间的中间值。这就是有效介电常数 $\varepsilon_{eff}$。
例如，对于FR4（$\varepsilon_r = 4.3$），当 $w/h = 2$ 时，$\varepsilon_{eff} \approx 3.3$。因此传播速度降低为 $v = c/\sqrt{\varepsilon_{eff}} \approx 1.65 \times 10^8\,\text{m/s}$，大约是真空中的55%。
然而，随着频率升高，电场更集中于介质中，因此 $\varepsilon_{eff}$ 随频率逐渐趋近于 $\varepsilon_r$。描述这种色散（dispersion）的Kirschning频率相关公式为：
其中 $f_p = Z_0 / (2\mu_0 h)$ 是频率参数，$G$ 是与几何形状相关的常数。这种色散会导致脉冲信号的高频分量延迟，引起边缘钝化，产生ISI（码间干扰）。
损耗机制
在实际基板中信号会衰减吧？有哪些类型的损耗呢？
主要分为三种：导体损耗、介质损耗和辐射损耗。在低GHz频段，导体损耗占主导；超过10 GHz，介质损耗的影响会变得显著。
导体损耗（$\alpha_c$）由于趋肤效应，随频率的平方根增加而增大：
其中 $R_s$ 是表面电阻，$\sigma$ 是电导率（铜：$5.8 \times 10^7\,\text{S/m}$）。此外，实际铜箔存在表面粗糙度，需要使用Hammerstad-Bekkadal修正将 $R_s$ 增大1.4至2倍。
介质损耗（$\alpha_d$）由损耗角正切 $\tan\delta$ 描述：
| 基板材料 | $\varepsilon_r$ | $\tan\delta$ | 用途 |
|---|---|---|---|
| FR4 | 4.2〜4.5 | 0.02 | 一般数字电路（〜3 GHz） |
| Megtron 6 | 3.7 | 0.002 | 高速数字电路（〜25 GHz） |
| Rogers RO4003C | 3.55 | 0.0027 | 射频/微波 |
| Rogers RT/duroid 5880 | 2.20 | 0.0009 | 毫米波（77 GHz车载雷达） |
| 氧化铝（$\text{Al}_2\text{O}_3$） | 9.8 | 0.0001 | MMIC基板 |
支配方程
CAE中求解的“方程”，归根结底是麦克斯韦方程组吗？
没错。所有电磁场分析的出发点都是麦克斯韦方程组。假设时谐场（$e^{j\omega t}$ 依赖），则得到频域形式。
将其关于 $\mathbf{E}$ 消去，得到矢量波动方程：
其中 $k_0 = \omega\sqrt{\mu_0\varepsilon_0} = 2\pi f / c$ 是自由空间波数。使用FEM、FDTD、MoM等数值方法离散化并求解此方程，就是微带线的电磁场仿真。
“准TEM模式”这一妥协的产物
微带线不支持真正的TEM波（纵向分量为零）。只要存在基板（介质）和空气两种媒质，就会传播混有TM分量的“准TEM模式”。特性阻抗可以用Hammerstad-Jensen近似公式（1980）在1%精度内求得，但在高频下色散会变得显著。实际上，在1970年代带状线的电气特性更优，但微带线便于元件安装，最终在“物理妥协”与“制造便利性”的权衡中，微带线取得了胜利。在CAE中，使用FEM、FDTD、MoM进行严格分析，用于“色散特性”的实际评估，并精确估算设计频率下的有效介电常数和损耗。
各变量的物理意义
- $w$（线宽）：导体带的宽度。是对 $Z_0$ 影响最大的设计参数。加宽线宽会降低阻抗（因为电容增大）。
- $h$（基板厚度）：介质基板的厚度。增加厚度会提高阻抗，同时辐射损耗也会增加。制造公差 $\pm 10\%$ 会导致 $Z_0$ 产生 $\pm 5\%$ 的变动。
- $\varepsilon_r$（相对介电常数）：基板材料的固有值。越大越能缩小线路尺寸，但损耗和色散会增加。
- $t$（导体厚度）：铜箔厚度（通常18〜35 $\mu$m）。当 $w/h$ 较小时，对 $Z_0$ 的影响较大。
- $\tan\delta$（损耗角正切）：介质损耗指标。FR4较大，为0.02；Rogers系列较小，为0.001〜0.003。
近似公式的适用限制
- 准TEM近似仅在 $f \ll f_{TE1}$（最低次TE模式的截止频率）时有效
- 超过 $f_{TE1} \approx c / (2w\sqrt{\varepsilon_r - 1})$ 会激励高次模，无法使用近似公式
- 线间耦合（串扰）无法用单根线路的公式评估。需要进行耦合线分析
- 过孔、弯曲、锥形等不连续部分必须进行3D电磁场分析
- 忽略基板材料的频率依赖性（因果性Dk/Df）会导致宽带特性误差增大
数值解法与实现
基于FEM的频域分析
用于微带线分析的数值方法，该如何区分使用呢？
主要有三种。FEM（有限元法）在频域中逐个频率求解。网格自适应（Adaptive Mesh Refinement）会自动运行，因此即使对于复杂的3D结构也易于使用。它是Ansys HFSS的主力引擎。
在FEM中，将矢量波动方程转换为弱形式，并使用边单元（Nedelec单元）进行离散化。由于是频域，在每个频率点求解
这样的线性方程组。其中 $[S]$ 是刚度矩阵（源于 $\nabla \times$ 算子），$[T]$ 是质量矩阵（源于 $\varepsilon_r$），$\{b\}$ 是端口激励向量。
Ansys HFSS的自动网格自适应：直到S参数的变化量 $|\Delta S_{ij}|$ 低于设定阈值（默认0.02），自动细化电场梯度较大区域的网格。通常需要3〜6次迭代（网格更新次数）收敛。
基于FDTD的时域分析
FDTD和FEM有什么区别呢？
FDTD在时域中直接对麦克斯韦方程组进行时间推进。只需一次注入高斯脉冲，经过FFT后即可一次性获得宽带的频率特性。CST Studio的时域求解器就是基于此方法。
那么FDTD效率更高吗？
视情况而定。FDTD使用正交网格（Yee网格），因此倾斜结构的表示精度依赖于网格。此外，对于高Q结构（如谐振器），振铃衰减时间很长，计算可能无法结束。另一方面，一次性获得宽带特性是其强大的优势。
Yee差分格式（1966）将 $\mathbf{E}$ 和 $\mathbf{H}$ 在空间和时间上交错半个网格放置：
稳定性条件（CFL条件）: $\Delta t \leq \frac{1}{c\sqrt{1/\Delta x^2 + 1/\Delta y^2 + 1/\Delta z^2}}$
| 方法 | 域 | 擅长结构 | 不擅长结构 | 代表工具 |
|---|---|---|---|---|
| FEM | 频域 | 复杂3D、曲面、非均匀材料 | 超宽带（多频率点） | Ansys HFSS |
| FDTD | 时域 | 宽带、瞬态响应 | 高Q谐振器、精细结构 | CST Studio |
| MoM | 频域 | 平面结构、开放问题 | 3D体积结构 | Sonnet, Momentum |
| 2.5D Planar | 频域 | 多层基板、高速设计 | 3D不连续部分 | Ansys SIwave |
MoM与Planar 3D
也经常听到MoM，它与微带线的匹配性如何？
MoM（矩量法）对平面结构非常强大。它使用格林函数直接求解导体表面的电流分布，因此不需要空气区域的网格，效率很高。Sonnet Suite和Keysight Momentum采用此方法。对于PCB设计者，被称为“Planar 3D”的2.5D方法很受欢迎，它将多层基板的每一层作为平面问题处理，并通过过孔连接层间。SIwave和Cadence Clarity 3D属于此类。
边单元与伪模
刚才提到了“边单元”，和普通的有限元有什么区别呢？
如果用通常的节点单元（将标量变量分配给节点）求解电磁场，会产生大量物理上不存在的虚假解（伪模）。边单元（Nedelec单元）通过将矢量场的切向分量分配给边，自动满足 $\nabla \cdot \mathbf{B} = 0$，从而排除伪模。在高频FEM中，边单元是事实上的标准。
频域 vs 时域的区分使用
频域分析类似于“将收音机调到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率分量的瞬态现象，但计算成本高。如果只需要微带线S参数的10个频率点，FEM更高效；如果需要DC〜40 GHz的全频带，FDTD更高效。
实践指南
分析流程
实际分析微带线时，应该从何开始呢？
基本流程有5个步骤。要点是首先通过手算得出目标值。如果直接跳到3D分析，将无法验证结果的正确性。
- 明确规格：确定目标阻抗、频率范围、插入损耗上限、基板叠层结构
- 手算/2D求解器进行初始设计：使用Hammerstad公式或SIwave确定线宽 $w$
- 3D建模：创建包含过孔、连接器、锥形等不连续部分的模型
- 执行电磁场仿真：获取S参数、电场分布、TDR阻抗剖面
- 后处理与设计迭代：眼图仿真、灵敏度分析、最终验证
网格设计要点
网格需要细化到什么程度呢？
微带线有特定的网格要求。导体边缘电场会奇异集中，因此边缘附近需要细化到最小波长 $\lambda_{min}$ 的 $1/20$ 以下。此外，基板厚度方向至少需要4〜6层网格。别忘了使用比趋肤深度 $\delta$ 更细的网格来评估导体损耗。
| 区域 | 推荐网格尺寸 | 理由 |
|---|---|---|
| 导体边缘附近 | $\lambda_{min}/20$ 以下 | 准确评估电场集中 |
| 基板厚度方向 | 4〜6层 | 确保 $\varepsilon_{eff}$ 精度 |
| 导体内部（评估损耗时） | 趋肤深度 $\delta$ 以下 | 表面电流密度精度 |
なった
詳しく
報告