槽隙天线的电磁场分析
理论与物理
巴比涅原理与互补性
缝隙天线真的只要在金属板上开个口就能发射电波吗？总觉得有点反直觉……
很好的问题。关键在于巴比涅原理。这原本是光学中的原理，说的是“与屏幕开口形状相同的障碍物会产生相同的衍射图案”。将其扩展到电磁波后，会发现有趣的现象。
光学原理和天线有关？怎么联系起来的？
假设在无限大导体平板上开一个长度为 $l$、宽度为 $w$ 的缝隙。布克在1946年指出，这个缝隙与相同尺寸的偶极子天线形成“互补”关系。具体来说，缝隙的电场分布对应于偶极子的磁场分布，反之亦然。
E场和H场互换……？那辐射方向图会怎样？
回想一下半波偶极子的辐射方向图。E面（电场方向）是8字形图案，H面是全向性的对吧。对于缝隙天线，E面和H面会完全互换。也就是说，偶极子的E面图案会成为缝隙的H面图案。极化方向也正交——如果偶极子是垂直极化，缝隙就是水平极化。
原来如此，图案形状相同，但电场和磁场的角色互换了。非常对称……
巴比涅原理的电磁波版本由布克（1946年）和索末菲公式化。在无限大理想导体平板上开的缝隙与同形状的金属片（偶极子）存在互补关系，并满足以下条件：
巴比涅原理（电磁波版）
$$\mathbf{E}_\text{slot}(\theta,\phi) = \pm\,\frac{1}{\eta}\,\mathbf{H}_\text{dipole}(\theta,\phi)$$其中 $\eta = \sqrt{\mu_0/\varepsilon_0} \approx 376.73\;\Omega$ 是自由空间阻抗。缝隙的电场与偶极子的磁场成正比，极化方向旋转90°。
缝隙阻抗的推导
方向图互换我明白了。那阻抗会怎样？看起来会和偶极子的73Ω有很大不同吧。
问得好。从巴比涅原理推导出的最重要关系式是这个：
布克的互补阻抗关系
$$Z_\text{slot} \cdot Z_\text{dipole} = \frac{\eta^2}{4}$$$Z_\text{dipole}$ 是同尺寸偶极子的输入阻抗，$Z_\text{slot}$ 是缝隙的输入阻抗。半波偶极子的谐振阻抗约为 $73 + j42\;\Omega$，因此在谐振时（忽略电抗部分，仅用实部计算）：
也就是说，半波缝隙的输入阻抗约为486Ω。比偶极子的73Ω高得多。这意味着与50Ω系统匹配很困难，实际应用中需要通过偏置馈电或调整谐振长度来处理。
486Ω好高啊！直接接到50Ω同轴电缆上大部分都会反射吧……怎么进行匹配呢？
实际应用中主要有三种方法：
- 偏置馈电： 将馈电点从缝隙中心偏移。越靠近电压分布的波节，阻抗越低。波导缝隙阵列中就是通过将缝隙从波导宽壁中心偏移来利用这一点。
- 调整缝隙宽度： 宽度 $w$ 增大，阻抗会变化。加宽宽度也能增加带宽，但存在交叉极化恶化的权衡。
- 背腔结构： 在缝隙背面设置腔体，实现单面辐射的同时，通过腔体尺寸调整阻抗。
谐振长度与带宽设计
缝隙的谐振长度和偶极子一样是半波长吗？
基本是这样。理想无限大导体平板上薄缝隙的谐振长度是 $l \approx \lambda/2$。但实际中会安装在基板上或切在波导壁上，因此会受到有效介电常数 $\varepsilon_\text{eff}$ 的影响：
缝隙的谐振条件
$$l_\text{res} \approx \frac{\lambda_0}{2\sqrt{\varepsilon_\text{eff}}} = \frac{c}{2f_r\sqrt{\varepsilon_\text{eff}}}$$其中 $c$ 是光速，$f_r$ 是谐振频率，$\varepsilon_\text{eff}$ 是有效介电常数。
原来如此。例如，在10GHz、$\varepsilon_\text{eff} = 2.2$（特氟龙基板）的情况下，谐振长度是多少？
我们来算一下。$\lambda_0 = c/f = 3\times10^8 / 10\times10^9 = 30\;\text{mm}$，所以：
$$l_\text{res} = \frac{30}{2\sqrt{2.2}} \approx \frac{30}{2 \times 1.483} \approx 10.1\;\text{mm}$$比自由空间的半波长15mm缩短了约33%。如果忽略这种缩短效应进行设计，谐振频率会产生很大偏移，必须注意。
带宽是怎么决定的？
缝隙的带宽主要由缝隙宽度 $w$ 与长度 $l$ 的比值 $w/l$ 决定。大致来说：
- 窄缝隙（$w/l < 0.05$）： 带宽为数百分比。谐振尖锐，用于需要滤波特性的场合。
- 宽缝隙（$w/l > 0.1$）： 带宽可达10%~20%。适用于宽带通信。
- 领结型缝隙： 中央加宽的形状，可实现40%以上的超宽带。
对于背腔结构，腔体深度（通常为 $\lambda/4$）也会影响带宽。深度越深带宽越宽，但存在外形变厚的权衡。
辐射方向图特性
无限大导体平板上的缝隙是向两面辐射的吧，但实际应用中单面辐射的情况更多，对吗？
没错。无限大平板上的缝隙向正反两面等量辐射（增益约 $5\;\text{dBi}$，高于半波偶极子的 $2.15\;\text{dBi}$）。但实际上，绝大多数场合需要单面辐射。因此采用以下方法：
- 背腔： 用腔体覆盖背面，只向前方辐射。雷达和相控阵中的标准做法。
- 波导壁缝隙： 波导本身起到腔体作用，缝隙向外辐射。
- 带接地板的基板缝隙： 在微带基板的接地板上开缝隙的方式。广泛用于WiFi和IoT天线。
波导缝隙阵列是通过排列缝隙来聚束波束的吧。一般排列多少个？
这取决于用途。船舶雷达可能是32~64个缝隙，军用相控阵阵列有时会达到数百甚至数千个缝隙。通过单独调整每个缝隙的偏移量，施加泰勒分布或切比雪夫分布的幅度锥削来抑制旁瓣，这是设计的关键。这正是CAE的用武之地，必须进行包含所有缝隙相互耦合的全波分析。
飞机的“隐形天线”——齐平安装的秘密
仔细观察喷气式客机的机身，会发现表面光滑没有突起，但ATC应答机和气象雷达却能工作。其秘密就在于缝隙天线。在金属机身上切割的细缝，基于巴比涅原理，能作为与偶极子等效的辐射器工作，无需外接天线。这可以避免增加空气阻力，因此被称为“齐平安装”，是超音速飞机和隐形飞机的必备技术。设计难点在于机身曲率会影响阻抗，导致与平板的理论值偏差很大，因此必须对实际形状进行CAE分析。
各项的物理意义
- 缝隙阻抗 $Z_\text{slot}$：缝隙馈电点处的电压与电流之比。无限大平板上的半波缝隙约为486Ω（纯电阻）。受有限尺寸平板、基板介质、邻近结构物的影响，此值会变化。
- 自由空间阻抗 $\eta = \sqrt{\mu_0/\varepsilon_0}$：真空中电场与磁场之比，约376.73Ω。是布克公式中的核心常数。
- 有效介电常数 $\varepsilon_\text{eff}$：对于基板上的缝隙，电场会穿过基板和空气，因此取基板介电常数 $\varepsilon_r$ 和空气介电常数1.0的中间值。通常使用修正后的微带线经验公式。
- 谐振长度 $l_\text{res} \approx \lambda/(2\sqrt{\varepsilon_\text{eff}})$：缝隙上形成半波长驻波的条件。由于边缘效应，有效长度比物理长度略长，因此实际谐振长度通常比计算值短2%~5%。
假设条件与适用范围
- 无限大导体平板假设：巴比涅原理以无限大平板为前提。对于有限尺寸平板，衍射效应会导致方向图变形。平板尺寸小于 $2\lambda$ 时尤为明显。
- 理想导体假设：忽略金属的有限电导率。在毫米波及以上频段，趋肤效应会导致增益下降，因此需要考虑有效电导率。
- 薄金属板假设：板厚 $t \ll \lambda$ 时布克公式有效。板厚变厚时，缝隙本身会表现出波导特性，出现截止效应。
- 单一缝隙假设：阵列排列时，单元间的相互耦合不可忽略，仅靠单独的巴比涅分析是不够的。
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 阻抗 $Z$ | Ω（欧姆） | 复数 $Z = R + jX$。$R$ 是辐射电阻，$X$ 是电抗。 |
| 自由空间阻抗 $\eta$ | Ω | $\eta_0 = 376.73\;\Omega$（≈120π Ω） |
| 波长 $\lambda$ | m | $\lambda = c/f$。10GHz时为30mm，5GHz时为60mm。 |
| 电场 $E$ | V/m | 缝隙上的电场分布呈正弦波状。 |
| 磁场 $H$ | A/m | 与缝隙边缘的电流密度直接相关。 |
| 增益 $G$ | dBi | 以各向同性天线为基准。半波缝隙（无限大平板）≈5 dBi。 |
数值解法与实现
矩量法（MoM）公式化
用计算机分析缝隙天线时，最经典的方法是什么？
从历史和实用性来看，首先是矩量法。缝隙天线本质上是口径问题——求解金属面开口部（缝隙）的电场分布作为未知函数。MoM在缝隙上假设等效磁流 $\mathbf{M}_s$，建立积分方程。
等效磁流是什么？物理上不存在磁流吧？
这是数学上的方便。根据缝隙上的切向电场 $\mathbf{E}_t$ 定义等效磁流：
$$\mathbf{M}_s = -2\hat{n} \times \mathbf{E}_t$$其中 $\hat{n}$ 是金属面的法向矢量。将这个 $\mathbf{M}_s$ 用基函数（通常是分段正弦波或RWG函数）展开，并用伽辽金法加权，就得到线性方程组 $[Z]\{M\} = \{V\}$。这里 $[Z]$ 是阻抗矩阵（稠密矩阵），$\{V\}$ 是激励向量。
稠密矩阵的话，单元数增加计算量不会爆炸性增长吗？
这正是MoM的弱点，内存需求为 $O(N^2)$，直接解法计算量为 $O(N^3)$。但对于缝隙天线，只需要对缝隙开口面和金属面表面进行网格划分，因此未知数远少于需要3D体网格的FEM。单个缝隙的话，数百到数千个未知数就足以获得足够的精度。对于大规模阵列，也有使用快速多极子法压缩到 $O(N\log N)$ 的方法。
FEM进行3D全波分析
像Ansys HFSS等使用的FEM在什么场合下会用到？
FEM在缝隙周围存在复杂3D结构时能发挥威力。例如：
- 背腔缝隙——包含腔体内部谐振模式的分析。
- 介质基板上的缝隙——精确建模基板的非均匀性或多层结构。
- 飞机机身上的安装——曲面导体上缝隙的辐射特性。
FEM将分析区域用四面体网格分割，并用边单元（Nedelec单元）对电场矢量进行插值。边界上配置吸收边界条件或完美匹配层来模拟无限空间。
边单元和FEM的节点单元有什么不同？
通常的节点单元适合处理标量（温度、电位），但像电场、磁场这样的矢量量如果用节点单元插值，会产生称为“伪模”的非物理解。边单元在单元的边上定义电场的切向分量，因此能自动满足 $\nabla \cdot \mathbf{E} = 0$，排除伪解。在高频电磁场分析中，边单元是事实上的标准。
用FEM求解频域麦克斯韦方程组时，弱形式如下：
FEM弱形式（频域）
$$\int_\Omega \left(\frac{1}{\mu_r}\nabla \times \mathbf{N}_i\right) \cdot \left(\nabla \times \mathbf{E}\right)\,d\Omega - k_0^2 \int_\Omega \varepsilon_r \mathbf{N}_i \cdot \mathbf{E}\,d\Omega = -j\omega\mu_0 \int_\Omega \mathbf{N}_i \cdot \mathbf{J}_s\,d\Omega$$其中 $\mathbf{N}_i$ 是边单元的基函数，$k_0 = \omega\sqrt{\mu_0\varepsilon_0}$ 是自由空间波数。
FDTD进行宽带分析
CST Studio使用的FDTD，和FEM怎么区分使用？
FDTD（时域有限差分法）直接在时域求解麦克斯韦方程组。输入高斯脉冲这样的宽带信号，追踪时间演变，然后对结果进行FFT，就能一次性得到宽频率范围的S参数和辐射方向图。对于缝隙阵列的频率扫描非常高效。
那全部用FDTD不就好了吗？
FDTD使用正交网格（Yee网格），因此对于斜面或曲面不得不采用阶梯近似。如果缝隙宽度在波长的1/100以下这样非常窄，就需要在那里使用极细的网格，导致计算成本飙升。相反，FEM使用非结构网格，可以只在缝隙周围局部细化。区分使用的标准是：
- FEM： 单一缝隙的精密分析、复杂结构内的缝隙、窄带谐振分析。
- FDTD： 宽带特性、阵列的快速筛选、瞬态现象。
- MoM： 大型
なった
詳しく
報告