電磁波伝搬
理论与物理
麦克斯韦方程与波动方程
老师，请告诉我电磁波的支配方程。
从麦克斯韦方程导出的波动方程：
相速度$v_p = 1/\sqrt{\mu\varepsilon}$。真空中光速$c = 3 \times 10^8$ m/s。
在电介质中速度会变慢呢。
在相对介电常数$\varepsilon_r$的介质中，$v_p = c/\sqrt{\varepsilon_r}$。FR-4基板（$\varepsilon_r \approx 4.4$）中约为$0.48c$。波长也会变短，因此在高频电路中尺寸会成为问题。
总结
- 波动方程 — 从麦克斯韦方程导出
- 相速度$v_p = c/\sqrt{\varepsilon_r}$ — 在介质中减速
- 有损介质 — 将$\varepsilon$扩展为复数$\varepsilon' - j\varepsilon''$
麦克斯韦在1865年预言的电磁波——直到赫兹证明花了23年
麦克斯韦从方程预言电磁波的存在是在1865年。然而，许多物理学家对“公式上的存在”持怀疑态度，实际上海因里希·赫兹在实验室实现电磁波的发送和接收验证是在1888年，整整23年之后。赫兹使用的是火花放电产生的数百MHz电波。如今高频CAE处理的数GHz～数百GHz电磁波理论，正是建立在这两位巨人的工作之上。每当凝视麦克斯韦方程时，请回想一下这23年的时间重量。
各项的物理意义
- 电场项 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：法拉第电磁感应定律。随时间变化的磁通密度产生电动势。【日常例子】自行车发电机（发电机）通过旋转磁铁使附近的线圈产生电压——这是磁场随时间变化会感应出电场这一定律的直接应用。IH电磁炉也基于相同原理，高频磁场的变化在锅底感应出涡流，通过焦耳热加热。
- 磁场项 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：安培-麦克斯韦定律。电流和位移电流产生磁场。【日常例子】电线通电时周围产生磁场——这就是安培定律。电磁铁根据此原理工作，通过线圈通电产生强磁场。智能手机的扬声器也应用了此定律：电流→磁场→振膜的力。在高频（GHz频段天线等）下，位移电流 $\partial D/\partial t$ 变得不可忽视，它描述了电磁波的辐射。
- 高斯定律 $\nabla \cdot \mathbf{D} = \rho_v$：表明电荷是电通量的发散源。【日常例子】用垫板摩擦头发会产生静电使头发竖起——带电的垫板（电荷）放射状地发出电力线，对轻的头发施加力。电容器（电容）设计中，用此定律计算电极间的电场分布。ESD（静电放电）对策也基于高斯定律的电场分析。
- 磁通守恒 $\nabla \cdot \mathbf{B} = 0$：表示不存在磁单极子。【日常例子】将条形磁铁切成两半，也无法得到只有N极或只有S极的磁铁——N极和S极总是成对存在。这意味着磁力线描绘的是“没有起点也没有终点的闭合回路”。在数值分析中，为了满足此条件，采用矢量势 $\mathbf{B} = \nabla \times \mathbf{A}$ 的公式化方法，自动保证磁通守恒。
假设条件与适用范围
- 线性材料假设：磁导率·介电常数不依赖于磁场·电场强度（在饱和区域需要非线性B-H曲线）
- 准静态近似（低频）：可忽略位移电流项（$\omega \varepsilon \ll \sigma$）。涡流分析中常用
- 2D假设（截面分析）：电流方向均匀且可忽略边缘效应时有效
- 各向同性假设：对于各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的案例：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
数值解法与实现
数值解法
用什么数值方法来求解电磁波？
| 方法 | 公式化 | 擅长的问题 |
|---|---|---|
| FEM | 频域 | 谐振腔、波导、复杂形状 |
| FDTD | 时域 | 宽带、瞬态响应、大规模 |
| MoM | 积分方程 | 开放空间、天线 |
| FIT | 积分形式麦克斯韦 | CST Studio Suite的基础 |
HFSS基于FEM，CST Studio Suite基于FIT/FDTD，FEKO（Altair）基于MoM。
网格尺寸如何决定？
以波长$\lambda$的1/10以下为基准。FEM使用二阶单元时，$\lambda/5$也能获得良好精度。FDTD需要满足CFL条件$\Delta t \leq \Delta x/(c\sqrt{3})$。
总结
- FEM (HFSS) — 复杂形状的频域分析
- FDTD (CST) — 宽带的时域分析
- $\lambda/10$ 规则 — 网格尺寸的基准
射线追踪法与移动通信基站的无线电波设计
在城市地区的移动电话基站设计中，考虑建筑物·地形引起的电磁波反射·绕射·散射的无线电波传播模拟不可或缺。全波FDTD精度虽高，但对于数百米见方的城市模型计算成本大到不切实际。此时发挥作用的是射线追踪（Ray Tracing）法，它通过光线追踪的类比，高速计算各路径的传播损耗·延迟。目前仍用于决定5G基站最优布局的“区域设计”，建筑物模型的精度极大地影响模拟精度。
边单元（Nedelec单元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料·非均匀介质。BEM: 自然处理无限区域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿·拉夫森法处理B-H曲线的非线性。残差基准: $||R||/||R_0|| < 10^{-4}$ 为一般标准。
频域分析
通过时间谐波假设归结为稳态问题。需要进行复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的区分使用
频域分析类似于“将收音机调到特定频率”——能高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——能再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实务中的应用
天线设计、波导设计、雷达散射截面（RCS）计算、无线电波传播预测是典型应用。
实务检查清单
传播损耗的测量——弗里斯传输方程与现实的差距
自由空间中天线间的接收功率可以用弗里斯（Friis）传输方程计算，但在实际现场，损耗通常比这个理论值大10〜40 dB。原因在于多径干扰·障碍物造成的遮蔽·地面反射等。“现场接收功率 vs 模拟值”的校准实测不可避免，测量点的选取·天线设置高度·周边环境的建模左右着精度。实践中，通常首先将弗里斯计算值作为上限来把握，再用环境系数（路径损耗指数）进行修正。
分析流程的比喻
电机的电磁场分析类似于“给吉他调音”。调整琴弦粗细（线圈匝数）和琴桥位置（磁铁配置），以引出最美妙的音色（高效的扭矩特性）。改变一个参数，整体的平衡就会改变——因此参数化研究很重要。
初学者容易陷入的陷阱
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“因为磁力线也会扩散到铁芯之外”。如果将分析区域紧贴铁芯，无处可去的磁通会“撞上”边界并反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球不断撞墙反弹的状态。
边界条件的思考方式
远方的边界条件看似不起眼却至关重要。需要在数值上表现“从这里开始是无限广阔的空间”。如果设置错误，磁通就会像撞上“看不见的墙”一样被反射回来。
软件比较
工具
| 工具 | 特点 |
|---|---|
| Ansys HFSS | 基于FEM。自适应网格。擅长复杂形状 |
| CST Studio Suite | FIT/FDTD。宽带。S参数提取 |
| COMSOL RF | 多物理场耦合。热-电磁耦合 |
| Altair FEKO | MoM为基础。 |
なった
詳しく
報告