2. 微带线特性阻抗全波提取与解析公式对比
为什么微带线的特性阻抗需要“全波”提取
微带线是准TEM传输线：低频时电场主要集中在导体带与地平面之间，可用准静态解析公式估算特性阻抗 Z0 与有效介电常数 εeff。但当频率升高，纵向场分量不可忽略，能量在介质与空气之间重新分配，εeff 随频率上升而增大，Z0 随之下降，这就是色散效应。Hammerstad 公式属于准静态闭式解，不含频率变量，因此高频段必然出现偏差。
本例以一条 50 Ω 量级的微带线为对象，用全波本征模求解提取传播常数 γ=α+jβ，再由 β 反推 εeff 与 Z0，并与 Hammerstad 解析值逐频点对比，量化色散偏差。
建模对象与参数设定
取 FR4 基板，相对介电常数 εr=4.4，损耗角正切 tanδ=0.02，基板厚度 h=1.6 mm，导体厚度 t=35 μm（1 oz 铜）。目标特性阻抗约 50 Ω，由 Hammerstad 反推得线宽 W≈3.0 mm。仿真频段 0.1–10 GHz，覆盖从准静态区到明显色散区。
| 参数 | 符号 | 取值 |
|---|---|---|
| 相对介电常数 | εr | 4.4 |
| 基板厚度 | h | 1.6 mm |
| 导体厚度 | t | 35 μm |
| 线宽 | W | 3.0 mm |
| 仿真频段 | f | 0.1–10 GHz |
全波提取流程与解析对比
全波提取的核心是求解给定截面上的本征模传播常数。对每个频点，设置本征模求解器，在端口截面上求解 Maxwell 方程，得到基模的复传播常数 γ=α+jβ。由 β 反推有效介电常数：
εeff(f) = (β·c / (2πf))²
Z0(f) = 1 / Re(γ / (jω·C')) → 工程上常用 Z0 = V²/(2P) 或 Z0 = 1/Y' 提取实际操作中，全波求解器直接输出 Z0 与 εeff，无需手工反推。关键设置：端口用波端口（wave port），截面尺寸足够大以避免高次模污染；网格在导体边缘加密，因为边缘场对 Z0 影响显著。
Hammerstad 解析公式（准静态）给出与频率无关的 Z0_ana 和 εeff_ana。将全波结果与解析值逐频点对比，即可得到色散偏差曲线。
关键判读：低频端（<1 GHz）全波与解析应吻合良好，偏差通常 <2%；频率升高后 εeff 上升、Z0 下降，偏差随频率单调增大。若低频端就出现大偏差，优先检查端口设置与网格，而非色散本身。
结果判读与偏差分析
将全波提取的 Z0(f) 与 εeff(f) 同 Hammerstad 解析值逐点比较，重点关注三个区间：
- 准静态区（0.1–1 GHz）：全波与解析偏差应小于 2%。若此处偏差偏大，说明端口或网格设置有问题，需先修正模型再讨论色散。
- 过渡区（1–6 GHz）：色散开始显现，
εeff缓慢上升，Z0缓慢下降，偏差逐渐增大到 3%–5%。 - 强色散区（6–10 GHz）：偏差可达 5%–10%，解析公式不再可靠，必须使用全波结果进行阻抗匹配设计。
工程建议：当工作频率高于 f_disp ≈ 0.3·c/(2π·h·√εr) 时，准静态公式的误差开始不可忽略。对 h=1.6 mm、εr=4.4 的 FR4 微带线，该阈值约在 3–4 GHz 附近，与仿真观察一致。
可复现的提取与对比清单
| 步骤 | 操作 | 检查点 |
|---|---|---|
| 1 | 建立截面模型，设置 εr、h、t、W | 尺寸单位统一为 mm |
| 2 | 设置波端口，截面延伸至空气层足够高 | 端口边缘远离导体，避免高次模 |
| 3 | 导体边缘网格加密，介质内至少 3 层网格 | 网格收敛性：加密后 Z0 变化 <0.5% |
| 4 | 逐频点求解本征模，记录 β 与 Z0 | 确认基模为准 TEM 模 |
| 5 | 计算 εeff = (βc/2πf)²，与解析值对比 | 低频端偏差 <2% |
| 6 | 绘制 ΔZ0 与 Δεeff 随频率曲线 | 偏差随频率单调增大 |