辐射方向图解析
辐射方向图的理论基础
概要
我经常看到辐射方向图的3D绘图,那到底是什么意思?
那是天线电界强度角度分布的展示。通俗讲,就是"电磁波在哪个方向以多大强度辐射出去"的立体图。一眼就能看出主瓣的方向和宽度、旁瓣电平、零点位置。
这在实际工作中有那么重要吗? 似乎只要知道"朝向哪里"就够了。
想得太简单了。比如卫星通信的抛物面天线,旁瓣电平必须控制在$-30\,\text{dB}$以下。旁瓣太高的话,会对相邻卫星造成干扰。军用雷达更严格,如果旁瓣泄漏信号,隐身性就完全没有了。FEM分析的精度直接影响设计品质。
那用什么物理量来评估呢?
主要有4个指标。远场电界$\mathbf{E}(\theta,\phi)$、指向性$D$、增益$G$、半功率波束宽度HPBW。我们逐个讲解。
远场电界表达式
什么叫"远场"?有明确的距离定义吗?
有的。天线最大尺寸为$D_{\max}$,波长为$\lambda$,那么远场的边界是
在这个距离以外,电磁波接近球面波,$\mathbf{E}$和$\mathbf{H}$都与$r$成反比,相互垂直,波阻抗均匀约为$\eta_0 \approx 377\,\Omega$的横向均匀波。
远场电界怎么数学表达?
从天线表面电流分布$\mathbf{J}(\mathbf{r}')$推导远场电界,得到下式:
用方向图函数$\mathbf{F}(\theta,\phi)$简化表示为
关键是,$\mathbf{F}(\theta,\phi)$相当于电流分布的傅里叶变换。因此阵列天线的方向图合成可以看作"空间频率滤波"。
傅里叶变换! 跟信号处理一样的思路。那均匀矩形电流分布的话,方向图是不是像sinc函数?
正是如此。等幅等相均匀分布主瓣最窄,但旁瓣只能到$-13.2\,\text{dB}$。用泰勒分布或切比雪夫分布给电流加渐进,主瓣稍微宽一点,旁瓣可以降到$-20$~$-40\,\text{dB}$。
辐射强度与指向性
有了方向图函数,接下来要定量评估"有多强"吧?
对的。首先定义辐射强度$U(\theta,\phi)$,就是单位立体角的辐射功率:
第二个等号很重要,说明在远场$U$不依赖距离,只是纯粹的角度函数。指向性$D$就是与全向辐射体相比的"峰值聚集度":
$P_{\text{rad}}$是向全方向辐射的总功率吧? 数值怎么算?
在环绕天线的球面$S$上对波印廷矢量进行数值积分。商用求解器比如HFSS,先把近场变换到远场(NFTF),再用辛普森法则或高斯求积在$(\theta,\phi)$网格上积分。通常$\theta$和$\phi$都取1°步长就足够精确了。
增益与HPBW
指向性$D$明白了。但数据表上通常写的是"增益$G$"吧? 有什么不同?
增益$G$包含了天线的损耗(导体损、介质损、表面波损等),是实际性能指标。用辐射效率$\eta$($0 < \eta \leq 1$)表示为
比如指向性$D = 30\,\text{dBi}$的抛物面天线,如果开口效率(照明效率+溅射损+表面精度损等)$\eta = 0.55$($-2.6\,\text{dB}$),实际增益就约$27.4\,\text{dBi}$了。
HPBW(半功率波束宽度)怎么求?
就是主瓣顶部功率下降$-3\,\text{dB}$处的两点的角度差:
比如半波长振子的E面HPBW约$78°$。而20波长口径的天线,$\text{HPBW} \approx 51\lambda/(D_{\max}) \approx 2.55°$,波束就很尖锐了。波束越尖锐,指向性越高,对远处目标的能量集中度就越好。
旁瓣与零点
旁瓣电平(SLL)越低越好对吧? 为啥要这么严格?
有几个原因。首先是干扰回避。卫星通信的ITU规范(比如ITU-R S.580)规定了旁瓣包络,$G(\theta) \leq 32 - 25\log_{10}\theta$ [dBi]必须满足,否则认证通不过。其次是杂波抑制。气象雷达的旁瓣如果拾取地面杂波,降雨量推估精度就会大打折扣。
零点(null)怎么用?
零点是辐射功率几乎为零的方向。用零点指向技术,可以把零点故意指向干扰源,从而消除干扰。5G基地站的Massive MIMO中,还有把零点对准人体来降低SAR(比吸收率)的技术已经商用了。
为什么用"极坐标图"而不是直角坐标
初次接触辐射方向图的人,常问为什么不用直角坐标而要用极坐标。直角坐标确实更容易读取数值,但极坐标的妙处在于"以天线为中心,直观看出各方向有多强"。旁瓣、反向瓣、零点的位置一眼明白,设计评审时"这个方向图触犯了ITU规范"能瞬间看出来。实务做法是用直角坐标比较精确数值,用极坐标检查整体形状。
辐射方向图的数值计算方法
解析方法的分类
方向图计算除了FEM,还有其他方法吧?
大概有3种主流方法,各有所长:
| 方法 | 原理 | 适用场景 | 计算复杂度 |
|---|---|---|---|
| FEM(有限元法) | 体积分割、边元素 | 复杂形状、非均质介质 | $O(N^{1.5})$~$O(N^2)$ |
| MoM(矩量法) | 表面积分方程 | 金属构造体、导线 | $O(N^2)$~$O(N^3)$ |
| FDTD(时域差分法) | Yee网格时间演进 | 宽带特性、过渡过程 | $O(N)$(单步) |
怎么选呢?
实务做法是这样:微带贴片天线、导波管槽阵等有密闭区域的结构用FEM。导线天线、金属板散射这类开放区域只需离散化金属表面的用MoM。超宽带或全频段特性一次计算的话用FDTD。现在混合方法(FEM-BIE、MLFMM加速的MoM)也很主流。
FEM定式化
FEM求方向图时,要解什么方程? 跟结构分析的$[K]\{u\}=\{F\}$一样吗?
思路差不多,但要解的是矢量亥姆霍兹方程:
其中$k_0 = \omega/c$是自由空间波数,$Z_0 = \sqrt{\mu_0/\varepsilon_0}$是自由空间波阻抗。FEM用边元素(Nédélec元素)离散化弱形式。边元素自动保证电界切线分量的连续性,消除了节点元素会出现的虚假模式(非物理的寄生解)。
边元素是指自由度在边上的吧。离散化后的矩阵方程是什么样?
用加权余量法离散后,最终得到:
$[S]$是旋度-旋度矩阵,$[T]$是质量矩阵,$\{b\}$是励振源项。形式跟结构分析$[K]\{u\}=\{F\}$一样,但矩阵是复数且非正定。因此不能用CG法,要用GMRES或MINRES这样的迭代求解器。
MoM定式化
矩量法跟FEM有什么本质区别?
FEM是把整个体积都分割网格,MoM只把金属表面上的电流分布作为未知量。电场积分方程(EFIE)是
$G(\mathbf{r},\mathbf{r}') = e^{-jk_0|\mathbf{r}-\mathbf{r}'|}/(4\pi|\mathbf{r}-\mathbf{r}'|)$是自由空间格林函数。用RWG(Rao-Wilton-Glisson)基函数离散化。MoM的强项是不用给空气区域网格。导线天线、偶极子这种开放结构效率高得多。缺点是得到密矩阵,大规模问题要用MLFMM(多层快速多极子法)加速到$O(N\log N)$。
近场到远场的变换
FEM算出来的结果是有限计算域内的,怎么得到远场?
这很关键。FEM/FDTD的计算域是有限的(用ABC/PML边界封闭),直接得不到远场。因此要用近场-远场变换(NFTF)。在环绕天线的虚拟闭曲面$S'$上取出$\mathbf{E}$和$\mathbf{H}$,转换为等效电流、磁流,再做傅里叶积分。
其中$\mathbf{J}_s = \hat{n}\times\mathbf{H}$、$\mathbf{M}_s = -\hat{n}\times\mathbf{E}$是等效电流、磁流。HFSS的这个变换是自动化的,用户只需设置"Radiation Boundary"就行。
网格与收敛
方向图计算的网格跟结构分析感觉不一样?
完全不一样。结构分析在应力集中处细分就行,但电磁场单元大小由波长决定。经验法则:
- 至少$\lambda/10$以下的单元(二阶元素)
- $\lambda/20$以下最安心(特别是介质边界附近)
- 曲面网格曲率半径的$1/5$以下
- 馈电点附近做局部加密(约$\lambda/50$)
HFSS的自适应网格(AMR)功能会自动加密单元,直到$\Delta S$(S参数变化量)低于阈值(通常0.02)为止。
收敛怎么判定?
方向图分析必须同时检查S参数和增益。S参数收敛了,增益(特别是旁瓣)还在动的情况常见。网格加倍后增益变化超过$0.3\,\text{dB}$就要警惕。最后要把变化控制在$0.1\,\text{dB}$以内才算安全。
辐射方向图的实务应用
解析流程
实际做方向图分析,要怎么一步步来?
大概5个阶段:
- 几何建模:天线本体+馈电结构+周边构件(罩、外壳)用CAD建立。有对称面就充分利用降低计算量
- 材料定义:导体(PEC或有限电导率)、介质基板($\varepsilon_r$、$\tan\delta$)、铁氧体等。注意$\tan\delta$的频率依赖性
- 励振设置:波导口(导波管、微带)或集总口(同轴、分立)。指定口的模式次数
- 边界+网格:辐射边界(PML或吸收边界)、对称面、初始网格设定
- 计算+后处理:S参数、方向图(E/H面切线、3D)、增益、效率全部检查
边界条件设置
辐射边界放在哪里? 紧挨天线可以吗?
千万别。近场的无源成分会反射回来搞坏结果。指导原则:
- PML(完全匹配层):离天线表面至少$\lambda/4$。PML本身厚度也要$\lambda/4$~$\lambda/2$
- 吸收边界条件(ABC):离天线$\lambda/4$~$\lambda/2$。但大入射角反射率会升高,需注意
- 对称面:PEC对称($\hat{n}\times\mathbf{E}=0$)和PMC对称($\hat{n}\times\mathbf{H}=0$)要按偏波方向区分用
PEC和PMC搞反了会怎样?
偏波分量会整个反转。比如振子的E面图如果PEC/PMC用反了,$E_\theta$和$E_\phi$互换,峰值方向就偏了90°。这是菜鸟最常犯的错误。正确做法是先拿简单的振子模型验证一遍,然后再上复杂模型。
网格设计要点
网格品质用什么指标评估?
| 指标 | 理想值 | 允许范围 | 电磁场特有注意 |
|---|---|---|---|
| 单元大小/λ | < 1/10 | < 1/6 | 介质中用λ_d = λ/√ε_r |
| 长宽比 | 1.0 | < 5.0 | PML层内可例外 |
| 雅可比比 | 1.0 | > 0.3 | 单元退化对边元素有害 |
| 曲率分割数 | ≥ 8 | ≥ 4 | 影响圆柱贴片的曲面精度 |
关键一点:介质内波长变短。$\varepsilon_r = 10$的基板中波长只有自由空间的$1/\sqrt{10} \approx 0.316$倍,网格也要同比细分。忽略这点的话,微带贴片的共振频率会偏几个百分点。
验证与有效性确认
分析结果怎么验证? 测方向图又很费时…
分阶段来:
- 理论值对比:半波振子($D = 2.15\,\text{dBi}$、$R_{\text{rad}} = 73.1\,\Omega$)或矩形开口(sinc方向图)验证求解器配置是否正确
- 交叉工具验证:同一模型用FEM和MoM都算一遍,增益和旁瓣在$0.5\,\text{dB}$以内就基本对
- 网格收敛:至少3个网格密度等级,增益变化$0.1\,\text{dB}$以内收敛
- 与测量对比:消声室或紧凑靶场测量。测量系统误差通常$\pm 0.5$~$\pm 1.0\,\text{dBi}$
天线方向图测量室的秘密
仿真值和测量值差异大多数来自测量环境。理想消声室误差在±0.5dBi以内,但实际的EMC试验室会因吸音材料老化和地面反射导致±1~2dBi偏差。拥有大型紧凑靶场的测试设备国内不超50套,价格以亿元计。"仿真正确、测量有问题"的情况在天线业界是公开秘密。
辐射方向图的软件比较
主要工具比较
方向图分析用什么商用工具比较好?
主流工具对比如下:
| 工具 | 开发方 | 主求解器 | 方向图支持 | 特点 |
|---|---|---|---|---|
| Ansys HFSS | Ansys Inc. | FEM(频域) | 3D/2D切线,NFTF自动 | 自适应网格强大。业界标准 |
| CST Studio Suite | Dassault (SIMULIA) | FDTD、FEM、MoM | 远场监视器自动 | 宽带分析强。EMC集成 |
| FEKO | Altair | MoM、MLFMM、PO/UTD | 3D、大规模阵列 | 混合方法强 |
| COMSOL (RF模块) | COMSOL AB | FEM | 远场计算节点 | 多物理场耦合容易 |
| NEC (4nec2) | 开源 | MoM | 3D方向图 | 免费。仅限导线天线 |
HFSS和CST的结果能对得上吗?
主瓣增益通常在$0.2$~$0.5\,\text{dB}$以内。旁瓣可能差$1$~$2\,\text{dB}$。特别是$-30\,\text{dB}$以下的低旁瓣,网格质量和边界条件微妙差别就会显著影响。用多个工具交叉验证是最保险的。
选择指南
按用途怎么选?
- 微带贴片、集成天线:HFSS首选。自适应网格收敛稳定
- 超大阵列(100素子以上):FEKO的MLFMM+区域分解。MoM开放域自然处理
- UWB、超宽带脉冲:CST时域求解器(一次计算全频段)
- 天线+筐体热耦合:COMSOL的RF+Heat Transfer模块
- 学习、导线天线概念:NEC(免费)足够
辐射方向图的先进研究
大规模MIMO与波束成形
5G的Massive MIMO跟方向图分析怎么关联?
Massive MIMO是64~256素子阵列,动态成形波束。方向图分析面临3大挑战:
- 计算规模:256素×数GHz自由度达数千万。需用区域分解或等效原理子域法
- 互耦合:素子间相互阻抗改变各素子的方向图(有源素子方向图)。单体方向图×AF不够
- 波束权重量化:数字移相器是6~8比特量化,零点方向会偏离设计值。仿真常用连续值,与实机脱节
256素子全波全场仿真能做吗?
实际不行。业界用的是混合法。中心几素子用全波求解(HFSS/CST)得有源素子方向图,再复制到阵列各位置,与阵列因子合成。周边素子的边缘效应用子阵补正。这样计算时间削减100倍以上。
机器学习融合
最近听说用机器学习预测方向图…
有两个热点方向:
- 代理模型:从参数扫描数据训练DNN快速预测方向图。优化中成本函数计算快1000倍。但外推能力弱
- PINN(物理信息神经网络):将麦克斯韦方程编入损失函数保证物理一致。无需网格,大规模问题有优势,但精度还达不到FEM
实际应用中,两者都是"用FEM/MoM替代",而是"快速缩小初期设计探索空间"的工具。
5G基站"绕过人体"的零点控制实际
5G NR的大规模MIMO已实现让方向图的零点指向人体方向来降低SAR的技术。基站通过GPS和惯导确定用户端末位置,实时更新波束权重。但实际问题是波束权重的量子化误差让零点容易"堵上"。精确的仿真与测量对比是必须的。这是个对"仿真≠现实"的深刻洞察。
辐射方向图的故障排除
方向图计算的典型问题
方向图仿真常见的"坑"有哪些?
| 现象 | 原因 | 对策 |
|---|---|---|
| 增益比理论值高3dB以上 | 对称面设置错(全空间功率未正确除) | 检查PEC/PMC,去掉对称重新算 |
| 旁瓣左右不对称 | 网格非对称或馈电点偏移 | 网格对称化、确认几何精度 |
| 反向瓣异常大 | 辐射边界距离不足或PML参数不当 | 离天线λ/4以上,检查PML设置 |
| 共振频率偏差 | 介质网格不足或材料参数错 | λ_d/10细分,确认ε_r、tanδ |
| 交叉偏波理论上为零却有输出 | 网格非对称或数值误差 | 对称条件严格设置,细分网格 |
对称面造成3dB差别,那用对称计算真的能省计算量吗?
能省一半。但前提是$P_{\text{rad}}$的积分要自动扩展到全空间。HFSS有自动补正,但COMSOL需手工乘2。这是新手最常掉的坑。
工具别错误对策
HFSS和CST特有的方向图错误有哪些?
常遇到的:
Ansys HFSS:
- "自适应网格未收敛":初始网格太粗或辐射边界太近。手动设λseeding
- "波导口模式计算失败":波导口尺寸不当(过大或过小)。调整口宽到基板的5~10倍
- "负辐射功率":辐射边界(Radiation Boundary)与PML干扰。检查距离
CST Studio Suite:
- "远场监视无效数据":远场频率超出计算范围或时间信号衰减不足(能量衰减<0.01%)
- "馈电点网格太小":自动生成网格在口处坍缩。Local Mesh指定最小单元大小
归根结底"仿真合不上"的万能法则是什么?
错误