偶极子天线的理论、数值解析与实践指南
偶极子天线的理论基础
概要 ― 什么是偶极子天线
偶极子天线是教科书中最先出现的天线，但它真的被广泛应用吗？
被广泛应用。Wi-Fi路由器上的棒状天线就是偶极子——更准确地说是基于接地平面的单极子，但原理相同。半波长偶极子是FM广播和业余无线电的标准天线。而且它还是衡量其他所有天线增益的参考基准（0dBd = 2.15dBi）。偶极子天线可以说是天线工程中的"标准米尺"。
那在结构上是怎样的呢？
两根导体棒（单元）排成一条直线，从中央给电。全长为半波长（$\lambda/2$）的称为半波长偶极子。比如2.4GHz的Wi-Fi，波长是125mm，所以每根棒约31mm。结构很简单，但要真正理解为什么它的输入阻抗是 $73 + j42.5\,\Omega$ 就需要扎实的电磁学基础。这往往成为检验电磁理论掌握程度的试金石。
微小偶极子的辐射
直接讲半波长太难了，能从更简单的模型开始吗？
好主意。我们先从微小偶极子（Hertzian dipole）讲起。当导体相对于波长很短（$l \ll \lambda$）时，假设流过恒定电流 $I_0$，远场电场为：
其中 $\eta_0 = 120\pi\,\Omega$ 是自由空间的阻抗，$k = 2\pi/\lambda$ 是波数。关键在于 $\sin\theta$ 项：在 $\theta = 90°$（垂直于天线轴）时最大，在 $\theta = 0°, 180°$（沿天线轴）时为零。这就是甜甜圈形的辐射方向图。
辐射阻抗怎么计算呢？
从远场的Poynting矢量 $\mathbf{S} = \frac{1}{2\eta_0}|E_\theta|^2 \hat{r}$ 开始，在全立体角上积分得到辐射功率 $P_{rad}$，然后用 $P_{rad} = \frac{1}{2}R_{rad}|I_0|^2$ 反推辐射阻抗：
比如当 $l = \lambda/10$ 时，$R_{rad} \approx 7.9\,\Omega$。相对波长越短，辐射阻抗越小，意味着辐射效率越差。这就是为什么人们常说"电气上短的天线效率不好"。手机天线设计中的很多困难都源于此。
半波长偶极子的电流分布与辐射方向图
微小偶极子过于理想化了。真正的半波长偶极子电流分布是怎样的？
说得好。半波长偶极子（全长 $L = \lambda/2$）的电流不均匀，而是正弦波分布：
在中点（给电点）达到最大值 $I_0$，在两个端点为零。从这个电流分布计算远场：
与微小偶极子的 $\sin\theta$ 相比，半波长偶极子的方向图 $\cos(\frac{\pi}{2}\cos\theta)/\sin\theta$ 稍微"收窄"了一些。最大辐射方向（$\theta = 90°$）相同，但沿轴向的零点更深。这个差异正是指向性差异的来源。
两端的电流为零这个直观上容易理解——开路端没有电流。但两端的电压（电场）却最大，对吗？
完全正确。这跟传输线路的开路端边界条件一样。电流为零，电压最大。偶极子天线的本质就是"开放式传输线路"，其中的电压和电流形成驻波。这个观点很关键。
辐射阻抗73Ω与输入阻抗的推导
73Ω这个数字看起来很奇怪。为什么是这个值？
从远场 $E_\theta$ 计算Poynting矢量 $\mathbf{S} = \frac{1}{2\eta_0}|E_\theta|^2 \hat{r}$，对全球面积分得到辐射功率：
这个积分没有闭式解，需要数值积分得到约1.219。由 $P_{rad} = \frac{1}{2}R_{rad}|I_0|^2$ 可得：
电磁教科书中"约73Ω"就是这样来的。有趣的是，75Ω同轴电缆的特征阻抗是当时为了与半波长偶极子的辐射阻抗相匹配而设定的，体现了天线和传输线的紧密联系。
输入阻抗中的虚数部分 $j42.5\,\Omega$ 从何而来？
天线近场中存储着电磁能量的无功成分。在恰好半波长（$L = \lambda/2$）的情况下，储存的电场能略大于磁场能，所以输入阻抗呈感性（正的电抗 $j42.5\,\Omega$）。
实务中为了消除这个电抗、使输入阻抗变成纯电阻，通常将天线长度缩短至约 $0.47\lambda$～$0.48\lambda$（约半波长的95%～96%）。这样电抗就趋于零，$Z_{in} \approx 73\,\Omega$（纯电阻）。但这个"共振长度"会因为线径和周围环境而微妙变化，最终需要通过仿真来精确调优。
指向性1.64（2.15dBi）的物理意义
天线规格书里经常看到2.15dBi，这到底是什么？
指向性 $D$ 定义为"最大辐射方向的辐射强度"除以"假如能量均匀向全方向辐射时的辐射强度"：
半波长偶极子的指向性为 $D = 1.643$。用dB表示：
意思就是"相比各向同性天线，在最强辐射方向强64%"。相同功率输入下，偶极子会把能量集中在赤道方向（甜甜圈的"赤带"上），因此该方向强度更高。电视天线（八木天线）增益能到10～12dBi，就是因为它进一步在偶极子基础上聚焦了波束。
我看过dBi和dBd两种记号。有什么区别吗？
dBd是以半波长偶极子为参考的增益。所以 $0\,\mathrm{dBd} = 2.15\,\mathrm{dBi}$。业余无线电圈子喜欢用dBd，但学术论文和CAE工具用dBi作标准。转换公式是 $G_{\mathrm{dBi}} = G_{\mathrm{dBd}} + 2.15$。
Hertz的实验（1888年）——最早的"天线"就是偶极子
1888年，海因里希·Hertz用两根金属棒（正是偶极子）首次实现了电磁波的发送和接收。他用火花放电激励高频振荡，在几米外的另一个接收偶极子上观察到火花——这是首次实验证实了Maxwell理论。波长约66厘米，相当于现在的450MHz频段。Hertz当时并不看好这项技术的实用价值，但仅十年后Marconi就用它改变了世界。
偶极子天线的数值计算手法
数值手法的详细说明
具体怎样用算法求解偶极子天线呢？
现在我理解为什么前辈说"偶极子天线一定要好好做"了。
离散化的定式
用形状函数 $N_i$ 来近似未知量：
用数式表示如下。
基础方程的离散形式
用数式表示如下。
仅仅看公式还是不太明白…这表示什么呢？
连续体的控制方程离散化后，得到如下代数方程组：
这里 $[K]$ 是总体刚度矩阵（或等效的系统矩阵），$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。
啊，我明白了！连续体的控制方程就是这样离散化的呀。
单元技术
听说过"单元技术"，但还没完全理解…
| 单元类型 | 次数 | 节点数(3D) | 精度 | 计算费用 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中等 |
| 六面体1阶 | 线性 | 8 | 中等 | 中等 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中等 |
积分方案
积分方案具体是怎么回事？
到这里为止，单元类型为什么重要，总算理解了！
收敛性和稳定性
不收敛了，首先要查什么？
收敛速度：二次单元误差以 $O(h^2)$ 的量级递减（光滑解的情况下）
啊！网格细分虽然看起来简单，其实很深奥啊。
求解器设置建议
具体怎样用算法求解偶极子天线呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数标准 |
| 预处理方法 | ILU(0) 或 AMG | 依问题规模 |
| 最大迭代次数 | 1000 | 不收敛时需重新审视设置 |
| 内存模式 | In-core | 尽可能采用 |
边单元（Nedelec单元）
为电磁场解析专门设计的单元。自动保证切向分量连续性，排除虚假模式。3D高频解析的标准。
节点单元
用于标量势定式化。在静磁场的标量势法或静电场解析中有效。
FEM vs BEM（边界单元法）
FEM：适应非线性材料和非均质介质。BEM：能自然处理无限域（开放区域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用Newton-Raphson法处理B-H曲线的非线性性。残差标准通常为 $||R||/||R_0|| < 10^{-4}$。
频域解析
在时间谐波假设下退化为稳态问题。需要复数运算，但宽带特性由时间域解析获取。
时间域的时间步长
需要时间步长 $\Delta t$ 小于最高频率分量周期的1/20。隐式时间积分允许更大步长，但要注意精度。
偶极子天线的实务应用
实践指南
先生，请讲讲"实践指南"！
讲解偶极子天线的实务解析流程和注意事项。
现在我理解为什么前辈说"偶极子天线一定要好好做"了。
解析流程
从第一步开始教我！我该从哪里开始？
1. 前处理 (Pre-processing)
- CAD数据导入和形状简化
- 材料特性定义
- 网格生成（单元类型、尺寸决定）
- 边界条件和荷载条件设定
2. 求解 (Solving)
- 求解器设置（求解法、收敛基准、输出控制）
- 作业投入和计算执行
- 收敛监控
3. 后处理 (Post-processing)
- 结果可视化（位移、应力、其他物理量）
- 结果验证和合理性确认
- 报告编制
网格生成的最佳实践
怎样判断网格好坏？
单元品质指标
请讲讲"单元品质指标"！
| 指标 | 理想值 | 许可范围 | 影响 |
|---|---|---|---|
| 长宽比 | 1.0 | < 5.0 | 精度下降 |
| Jacobian比 | 1.0 | > 0.3 | 单元退化 |
| 翘曲 | 0° | < 15° | 精度下降 |
| 斜度 | 0° | < 45° | 收敛性恶化 |
| 锥度比 | 0 | < 0.5 | 精度下降 |
网格密度的决定
网格密度的决定具体怎么操作？
边界条件的设置指南
听说边界条件搞错了全盘皆输…
啊！过约束注意——我总算明白这个机制了。
商用工具别的实现步骤
有那么多软件，各自的特点是什么？
| 工具名 | 开发方/现所属 | 主要文件格式 |
|---|---|---|
| Ansys HFSS | ANSYS Inc. | .aedt, .hfss |
| CST Studio Suite | Dassault Systèmes SIMULIA | .cst |
| COMSOL Multiphysics | COMSOL AB | .mph |
Ansys HFSS
接下来是Ansys HFSS的内容吧。讲讲看？
由Ansoft Corporation开发的3D高频电磁场仿真器。2008年被ANSYS收购。
现归属：ANSYS Inc.
CST Studio Suite
CST Studio具体是怎样的？
由Computer Simulation Technology（德国）开发。2016年被Dassault Systèmes收购，纳入SIMULIA。
现归属：Dassault Systèmes SIMULIA
先生讲得清楚！工具名字的迷茫消散了。
COMSOL Multiphysics
请讲讲"COMSOL Multiphysics"！
1986年在瑞典成立。作为与MATLAB相互配合的FEMLAB开始，后改名为COMSOL。多物理特别强。
现归属：COMSOL AB
功能比较矩阵
预算和时间都有限，成本效益最高的是哪个？
| 功能 | HFSS | CST | COMSOL |
|---|---|---|---|
| 基本功能 | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ |
| GPU支持 | △ | △ | ○ |
转换的风险
转换时的风险具体是什么？
啊！不同工具间的模型转移就是这样工作的呀。
许可形式
听说过"许可形式"，但还没完全理解…
| 工具 | 许可 | 特点 |
|---|---|---|
| 商业FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持收费 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的开源求解器 |
选定指南
最后，怎样判断该选哪个呢？
选择偶极子天线工具时要考虑：
偶极子天线的全貌我掌握了！明天开始在实务中用上。
好！实际动手是最好的学习。有什么不懂的随时问我。
偶极子与巴伦——"不装巴伦，同轴电缆就成了天线"
用同轴电缆给偶极子天线供电时，如果省去巴伦（平衡-不平衡变换器），同轴外屏会产生电流，外屏本身也成了辐射单元。这叫"外屏电流"，会导致辐射方向图畸变，还可能超过设计值的电磁露出。实务中常见"仿真OK但实测方向图偏"的情况，罪魁祸首往往就是这个。在仿真中是否包含巴伦，对解析精度有很大影响。
偶极子天线的软件工具比较
商用工具比较
有那么多软件，各自的特点是什么？
对偶极子天线适配的主要商业CAE工具的功能比较和各产品历史背景做详细说明。
现在我理解为什么前辈说"偶极子天线一定要好好做"了。
适配工具清单
那做偶极子天线解析用哪些软件呢？
| 工具名 | 开发方/现所属 | 主要文件格式 |
|---|---|---|
| Ansys HFSS | ANSYS Inc. | .aedt, .hfss |
| CST Studio Suite | Dassault Systèmes SIMULIA | .cst |
| COMSOL Multiphysics | COMSOL AB | .mph |
Ansys HFSS
接下来是Ansys HFSS的内容吧。讲讲看？
由Ansoft Corporation开发的3D高频电磁场仿真器。2008年被ANSYS收购。
现归属：ANSYS Inc.
CST Studio Suite
CST Studio具体是怎样的？
由Computer Simulation Technology（德国）开发。2016年被Dassault Systèmes收购，纳入SIMULIA。
现归属：Dassault Systèmes SIMULIA
COMSOL Multiphysics
请讲讲"COMSOL Multiphysics"！
1986年在瑞典成立。作为与MATLAB相互配合的FEMLAB开始，后改名为COMSOL。多物理特别强。
现归属：COMSOL AB
功能比较矩阵
预算和时间都有限，成本效益最高的是哪个？
| 功能 | HFSS | CST | COMSOL |
|---|---|---|---|
| 基本功能 | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ |
| GPU支持 | △ | △ | ○ |
转换的风险
转换时的风险具体是什么？
啊！不同工具间的模型转移就是这样工作的呀。
许可形式
听说过"许可形式"，但还没完全理解…
| 工具 | 许可 | 特点 |
|---|---|---|
| 商业FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持收费 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的开源求解器 |
选定指南
最后，怎样判断该选哪个呢？
选择偶极子天线工具时要考虑：
偶极子天线的全貌我掌握了！明天开始在实务中用上。
好！实际动手是最好的学习。有什么不懂的随时问我。
NEC（Numerical Electromagnetics Code）——免费且易上手的老牌求解器
偶极子天线解析的经典工具之一是美国Lawrence Livermore国家实验室在1970年代开发的NEC（Numerical Electromagnetics Code）。基于矩量法，能以高精度计算线状导体天线的电流分布和辐射方向图。至今仍免费发布，配合GUI前端XNec2c或NEC2c，从业余无线电到卫星通信天线设计都能用。与商业工具相比精度毫不逊色，"想用低成本做正经解析"的工程师的强大后盾。
偶极子天线的先端研究
先端主题与研究动向
偶极子天线领域今后会怎样发展？
看看偶极子天线领域的最新研究动向和先进手法。
现在我理解为什么前辈说"偶极子天线一定要好好做"了。
最新的数值手法
接下来是最新数值手法的内容吧。讲讲看？
仅仅看公式还是不太明白…这表示什么呢？
高性能计算 (HPC) 的应对
| 并行化手法 | 概述 | 适用求解器 |
|---|---|---|
| MPI (区域分割) | 分布式内存型。大规模问题标配 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU利用。特别对显式法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
偶极子天线的故障排查
故障排查
现在我理解为什么前辈说"偶极子天线一定要好好做"了。
常见错误与对策
先生，你有过偶极子天线问题导致的通宵调试吗？（笑）
1. 不收敛
不收敛具体是怎样的情况？
症状：求解器在指定迭代次数内未收敛，异常退出
可能的原因：
- 网格品质不足（过度扭曲的单元）
- 材料参数设置不恰当
- 初始条件不合理
- 非线性太强（荷载步骤不足）
对策：
- 进行网格品质检查（长宽比、Jacobian）
- 确认材料参数的单位
- 分段加载（增加子步数）
- 放宽收敛判定基准（但要注意精度）
也就是在不收敛的地方加工夫，否则后面要吃苦头。记住了！
2. 非物理结果
接下来是非物理结果的内容吧。讲讲看？
症状：应力/位移/温度等出现物理上非现实的值
可能的原因：
- 边界条件设置错误
- 单位制混用（SI单位与工程单位混淆）
- 单元类型选择不适当
- 应力奇异点存在
对策：
- 验证反力合计（力的平衡）
- 检查单位制一致性
- 重新审视单元类型选择
- 排除奇异点或子模型化
前辈说"不收敛一定要好好做"的意思现在明白了。
3. 计算时间超长
计算时间超长具体是什么情况？
症状：计算耗时数倍于预期
对策：
- 优化网格粗密分布
- 利用对称性（1/2、1/4模型）
- 优化求解器设置（迭代法、预处理选择）
- 利用并行计算
4. 内存不足
请讲讲"内存不足"！
症状：Out of Memory 错误
前辈说"不收敛一定要好好做"的意思现在明白了。
对策：
- 使用出核求解法
- 缩减网格规模
- 确认64位版求解器
- 增加内存分配
啊～，关于不收敛的讨论，超级有意思！请多讲些。
Nastran代表错误
代表错误具体是什么样的？
Abaqus代表错误
请讲讲"代表错误"！
那如果记住工具名，基本就没问题吧？
"解析不合"时应该怎么办
- 先深呼吸——焦躁中随意改设置会把问题复杂化
- 制作最小重现案例——用最简形式重现偶极子天线问题。"减法式调试"最高效
- 一次只改一个——同时改多个会分不清哪个有效。学科学实验，做"对照实验"
- 回到物理本质——结果若"重力下物体反向浮起"这样违背物理，要怀疑输入数据根本错误
更多
错误