抛物面反射镜天线的CAE分析
抛物面反射镜天线CAE的理论基础
概述 — 为什么抛物天线分析很困难
抛物天线分析有什么难的？不就是电波在抛物面上反射吗？
好问题。电波被转换成平面波的原理本身很简单。但在实际分析中，对于在电气上很大的天线——即口径是波长的100倍以上的那种——全波FEM分析会导致内存完全不足。比如Ku波段（12 GHz，波长25 mm）的1.2 m直径卫星广播天线，$D/\lambda = 48$ 左右。用3D FEM试图求解这个问题，自由度会达到数千万。
哎，直径仅仅1米多的天线会这么重吗？
因为电磁场网格需要约波长的1/10～1/6的要素尺寸。$D/\lambda = 48$ 的话，单方向约300个单元，三维下是几百万单元。直径超过10 m的射电望远镜时，$D/\lambda$ 达到数千，全波法就完全不现实了。所以需要物理光学(PO)或GTD(几何光学衍射理论)等渐近法与混合方法结合。副反射镜的遮挡和馈源位相中心偏差也会直接影响增益，这些细节也很关键。
抛物面的几何和焦点特性
先从基础讲起。"抛物面"为什么能集中电波？
抛物面可以用圆柱坐标 $(r, \phi, z)$ 表示为：
其中 $f$ 是焦距。这个抛物面在数学上的重要性质是，从焦点出发的球面波在反射面上的所有光程都相等。也就是说，反射后的波前是完全平面波。
我明白了，位相一致所以通过干涉相加。焦点距和口径的比(f/D比)怎么确定？
$f/D$ 比是设计中最基本的参数。总结一般范围和特点如下：
| $f/D$ 范围 | 半开角 $\theta_0$ | 特点 | 应用例 |
|---|---|---|---|
| 0.25～0.35 | 90°～70° | 深碟形。馈源指向性要求松散 | 卫星通信地面站 |
| 0.35～0.50 | 70°～53° | 最常见。溅射和照射效率平衡良好 | 雷达、VSAT |
| 0.50～0.80 | 53°～35° | 浅碟形。交叉极化特性良好 | 射电望远镜 |
半开角 $\theta_0$ 和 $f/D$ 的关系为：
开口效率和子效率的分解
"开口效率"经常听到，具体是什么概念？
抛物天线的增益用下式表示：
其中 $\eta_{ap}$ 是开口效率(aperture efficiency)，表示相对理想均匀照射的实际增益比率。这个 $\eta_{ap}$ 可以分解为多个子效率的乘积：
| 记号 | 名称 | 物理意义 | 典型值 |
|---|---|---|---|
| $\eta_{ill}$ | 照射效率 | 开口面上电场分布的均匀性 | 0.75～0.85 |
| $\eta_{sp}$ | 溅射效率 | 馈源功率从反射镜外漏的损耗 | 0.85～0.95 |
| $\eta_{ph}$ | 相位效率 | 开口面上相位误差导致的损耗 | 0.95～0.99 |
| $\eta_{pol}$ | 极化效率 | 交叉极化分量导致的损耗 | 0.95～0.99 |
| $\eta_{bl}$ | 遮挡效率 | 副反射镜和支承结构的阻挡 | 0.90～0.98 |
| $\eta_{sf}$ | 表面精度效率 | 镜面制造误差(Ruze公式) | 0.85～0.98 |
表面精度效率那里的"Ruze公式"是什么公式？
Ruze的公式是对反射镜面RMS误差 $\epsilon$ 对增益的影响进行量化的著名公式：
例如，$\epsilon/\lambda = 1/20$ 时，$\eta_{sf} \approx 0.67$（降低30%）。当 $\epsilon/\lambda = 1/50$ 时，$\eta_{sf} \approx 0.94$。射电望远镜ALMA的反射镜要求RMS 25 μm精度的原因正是这个，最短观测波长0.3 mm时需要维持 $\epsilon/\lambda \approx 1/12$ 。
支配方程 — 开口面积分和增益
从反射镜开口面上的等效电流，用开口面积分(aperture integration)计算远场。设开口面上的等效电场分布为 $\mathbf{E}_{ap}(x',y')$，则远场由下面的傅里叶变换型积分给出：
这里 $k = 2\pi/\lambda$ 是波数，$A$ 是开口面积。对于轴对称的馈源图样，这个二维积分可以简化为关于 $r$ 的一维积分(贝塞尔函数变换)。增益可表示为：
明白了，只要知道开口面上的电场分布就能计算远场。这个"开口面电场"怎么求，就是分析方法的区别吧？
正是这样。用PO法近似求解，还是用FEM精确求解，就由"电气尺寸"决定。
抛物面的焦点收敛——为什么不用"抛物线"就无法追踪通信卫星
当把馈源元件放在抛物面焦点时，镜面反射的电波作为平行光（平面波）辐射。这是由抛物线的定义"焦点距离与准线距离相等"引出的性质，从焦点出射的球面波能以完全相同的光程到达镜面，所以位相一致。椭圆面或双曲面的这个性质单独不成立，因此为了以高增益接收远处点源（卫星），抛物面是基本配置。不过卡塞格伦或格里高利结构用双曲面和椭圆面副反射镜与主镜搭配，可以改变等效焦距。
离散化方法
这个方程在电脑上怎么实际求解？
用有限元法(FEM)进行空间离散。组装单元刚度矩阵，构建整体刚度方程。
进行弱形式(变分形式)变换，用试验函数和形状函数利用Galerkin法进行定式化。单元类型的选择(低阶单元 vs. 高阶单元、完全积分 vs. 低减积分)直接关系到解的精度和计算成本的折衷。
矩阵解法算法
矩阵解法算法具体是什么意思？
用直接法(LU分解、Cholesky分解)或迭代法(CG法、GMRES法)求解联立方程。对大规模问题，前处理结合的迭代法很有效。
| 求解方法 | 分类 | 内存使用 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky分解 | 直接法(对称正定) | O(n²) | 小～中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG前处理 | 前处理 | O(n) | 超大规模 |
换句话说，在有限元法这块偷工减料会后来吃苦头吧。铭记于心！
商用工具中的实现
要做抛物反射镜天线得用什么软件？
| 工具名称 | 开发方/现在 | 主要文件格式 |
|---|---|---|
| Ansys HFSS | Ansys Inc. | .aedt, .hfss |
| CST Studio Suite | Dassault Systèmes SIMULIA | .cst |
| COMSOL Multiphysics | COMSOL AB | .mph |
供应商系谱和产品整合经过
各个软件的由来，有没有什么戏剧化的故事？
Ansys HFSS
接下来讲Ansys HFSS的事，内容是什么？
Ansoft公司开发的3D高频电磁场模拟器。2008年Ansys收购了Ansoft。
现在所属: 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
文件格式和互操作性
不同软件间传递数据时有什么注意事项吗？
| 格式 | 扩展名 | 种类 | 概述 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | ISO 10303准的3D CAD数据交换格式。支持几何+PMI。 |
| IGES | .igs/.iges | 中立CAD | 初期的CAD数据交换规格。曲面数据互操作有问题。STEP向转移。 |
| STL | .stl | 网格 | 用三角形小面近似3D形状。3D打印标准。CAE网格不太适合。 |
不同求解器间转换模型时，要注意单元类型对应关系、材料模型兼容性、荷载和边界条件的表示差异。特别是高阶单元或特殊单元(内聚单元、用户定义单元等)往往不能在求解器间直接转换。
所以文件格式表面简单，实际上藏着深情啊。
实务上的注意点
教科书没有的"现场智慧"之类的东西有吗？
网格收敛性确认、边界条件合理性验证、材料参数灵敏度分析非常重要。
呃，抛物反射镜天线真是博大精深啊… 但多亏老师的讲解，总算整理清楚了！
嗯，很好！实际动手试验是最好的学习。有不明白的地方随时来问。
抛物面的焦点收敛——为什么不用"抛物线"就无法追踪通信卫星
当把馈源元件放在抛物面焦点时，镜面反射的电波作为平行光（平面波）辐射。这是由抛物线的定义"焦点距离与准线距离相等"引出的性质，从焦点出射的球面波能以完全相同的光程到达镜面，所以位相一致。椭圆面或双曲面单独不具这个性质，所以为以高增益接收遥远点源（卫星），抛物面成为基本配置。不过在卡塞格伦或格里高利构型中，用双曲面和椭圆面副反射镜与主镜搭配，可改变等效焦距。
抛物面反射镜天线CAE的数值计算方法
数值方法的详细
具体怎么用什么算法求解抛物反射镜天线？
哦，原来这样！抛物反射镜天线是这样的机制啊。
离散化的定式化
用形状函数 $N_i$ 近似未知量:
用公式表示是这样的。
基础方程式的离散形
用公式表示是这样的。
嗯，光看公式有点摸不着头脑… 这是在表示什么？
连续体的支配方程离散化后，得到下面的代数方程组:
这里 $[K]$ 是整体刚度矩阵(或等效的系统矩阵)，$\{u\}$ 是未知节点变量向量，$\{F\}$ 是外力向量。
哦，原来如此！连续体的支配方程这样求解啊。
单元技术
听说过"单元技术"，但可能没有完全理解…
| 单元类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1次 | 线性 | 4 | 低(剪切锁定) | 低 |
| 四面体2次 | 二次 | 10 | 高 | 中 |
| 六面体1次 | 线性 | 8 | 中 | 中 |
| 六面体2次 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中 |
积分方案
积分方案具体是什么意思？
听到这儿，终于明白为啥单元类型这么重要了！
收敛性和稳定性
不收敛了的话，首先检查什么？
收敛速度: 二次单元时误差以 $O(h^2)$ 阶减小(平滑解的情况)
所以，细化网格看似简单，实际上也很深奥啊。
求解器设置的推荐
具体用什么算法求解抛物反射镜天线？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法的收敛判定 | $10^{-6}$ | 残差范数基准 |
| 前处理手法 | ILU(0) or AMG | 按问题规模 |
| 最大迭代次数 | 1000 | 不收敛时需重新设置 |
| 内存模式 | In-core | 尽可能使用 |
边单元(Nedelec单元)
为电磁场解析设计的单元。自动保证切向分量的连续性，消除虚假模式。3D高频分析的标准。
节点单元
用节点值插值标量场。用于电位、温度分析。对矢量场会有虚假解问题。
FEM vs BEM(边界单元法)
FEM: 对非线性材料和非均质媒质适用。BEM: 无限领域(开领域问题)能自然处理。混合FEM-BEM也有效。
非线性收敛(磁饱和)
B-H曲线的非线性用牛顿-拉夫逊法处理。残差基准: $||R||/||R_0|| < 10^{-4}$ 一般惯例。
频域分析
时间谐波假设简化为稳态问题。需复数运算，但宽带特性用时域分析取得。
时间域的时间步长
最高频率分量的1/20以下的时间步长必要。隐式时间积分虽然步长可大些，但精度需注意。
抛物面反射镜天线CAE的实务应用
实践指南
老师，请教我"实践指南"！
讲解抛物反射镜天线的实务分析流程和注意点。
哦，原来如此！抛物反射镜天线是这样的机制啊。
分析流程
从第一步开始教起！首先做什么？
1. 前处理(Pre-processing)
- CAD数据导入和形状简化
- 材料特性定义
- 网格生成(单元类型和尺寸决定)
- 边界条件和荷载条件设置
2. 求解(Solving)
- 求解器设置(解法、收敛基准、输出控制)
- 作业投入和计算执行
- 收敛监测
3. 后处理(Post-processing)
- 结果可视化(位移、应力等物理量)
- 结果验证和合理性确认
- 报告编制
网格生成的最佳实践
怎样判断网格的好坏？
单元品质指标
请教我"单元品质指标"！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 长宽比 | 1.0 | < 5.0 | 精度下降 |
| 雅可比比 | 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公司开发的3D高频电磁场模拟器。2008年Ansys收购了Ansoft。
现在所属: Ansys Inc.
CST Studio Suite
CST Studio具体是怎么回事？
Computer Simulation Technology(德国)开发。2016年Dassault Systèmes收购并整合入SIMULIA。
现在所属: Dassault Systèmes SIMULIA
老师的说明好清楚！工具名称的疑惑晴了。
常见失败和对策
初学者容易犯什么错？想事先知道！
| 现象 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、边界条件不适当 | 改善网格、修改约束条件 |
| 应力异常大 | 应力奇异点、网格依赖 | 避开奇异点、局部网格细分 |
| 位移不现实 | 材料常数误误、单位系混杂 | 确认输入数据 |
| 计算时间过长 | 不必要的细分、低效求解 | 网格优化、并行计算 |
品质保证检查清单
教科书没有的"现场智慧"之类的东西有吗？
呃，抛物反射镜天线真是博大精深啊… 但多亏老师的讲解，总算整理清楚了！
嗯，很好！实际动手试验是最好的学习。有不明白的地方随时来问。
电波望远镜ALMA的反射镜精度——0.025mm的误差映射宇宙
在智利阿塔卡玛沙漠设置的电波望远镜ALMA(阿尔玛)的抛物反射镜，直径12 m的镜面表面精度管理在RMS 25μm(0.025mm)以下。这相当于观测波长(最短0.3 mm)的约1/12。反射镜的变形分析采用有限元法的结构分析与电磁场分析联成的手法，气温变化和重力方向变化导致的反射镜弯曲通过主动表面机制实时补正。"电磁场CAE与结构CAE的联成"实际决定着太空观测精度这一事实，让设计者腰板挺得笔直。
抛物面反射镜天线CAE的软件比较
商用工具比较
有这么多软件呀？各自的特点教我！
详述对应抛物反射镜天线的主要商用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公司开发的3D高频电磁场模拟器。2008年Ansys收购了Ansoft。
现在所属: 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开发的OSS求解器 |
选择的指针
到底该选哪个，判断基准请教我？
抛物反射镜天线工具选择时要考虑以下因素:
呃，抛物反射镜天线真是博大精深啊… 但多亏老师的讲解，总算整理清楚了！
嗯，很好！实际动手试验是最好的学习。有不明白的地方随时来问。
JPL的格里高利反射镜——深空通信的仿真规模
NASA喷气推进研究所(JPL)的深空通信设施运作着直径70 m的巨大抛物天线。这样超大型反射镜的电磁场分析对常规FDTD和MoM来说计算量太大，所以采用射线追踪法(Ray Tracing)和物理光学法(PO)的混合解法。探测器旅行者1号目前离地球240多亿公里，却能用这个70 m天线仅用几W发射功率通信。"能把信号传这么远"是用仿真优化天线设计的成果。
抛物面反射镜天线CAE的前沿研究
前沿课题与研究动向
抛物反射镜天线这个领域以后怎么发展？
看看抛物反射镜天线中最新研究动向和先进手法。
哦，原来如此！抛物反射镜天线是这样的机制啊。
最新的数值方法
接下来讲最新数值方法的事，内容是什么？
嗯，光看公式有点摸不着头脑… 这是在表示什么？
高性能计算(HPC)的对应
| 并行化手法 | 概述 | 适用求解器 |
|---|---|---|
| MPI(领域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU(CUDA/OpenCL) | GPGPU活用。陆解法特别有效 | LS-DYNA, Fluent等 |
| 混合MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
抛物面反射镜天线CAE的故障排除
故障排除
哦，原来如此！抛物反射镜天线是这样的机制啊。
常见错误和对策
老师也在抛物反射镜天线里通宵debug过吗？(笑)
1. 收敛失败
收敛失败具体是什么意思？
现象: 求解器未在指定迭代回数内收敛异常结束
可能原因:
- 网格品质不足(过度扭曲单元)
- 材料参数设置不当
- 不恰当的初始条件
- 非线性太强(荷载步不足)
对策:
- 执行网格品质检查(长宽比、雅可比)
- 确认材料参数的单位系
- 将荷载分割为多个步(增加子步数)
- 放松收敛判定基准(但要注意精度)
换句话说，在收敛失败这块偷工减料会后来吃苦头吧。铭记于心！
2. 非物理的结果
接下来讲非物理的结果的事，内容是什么？
现象: 应力/位移/温度等物理上不现实
可能原因:
- 边界条件误设
- 单位系混杂(SI单位与工学单位混淆)
- 不适当的单元类型选择
- 应力奇异点存在
对策:
- 确认反力合计(力的平衡)
- 确认单位系的一致性
- 重新考虑单元类型的适当性
- 去除奇异点或进行子建模
前辈说"收敛失败一定要好好做"的含义明白了。
3. 计算时间超过
计算时间超过具体是什么意思？
现象: 计算花费想定时间的数倍
对策:
- 优化网格的粗密分布
- 活用对称性(1/2、1/4模型)
- 优化求解器设置(迭代法、前处理的选择)
- 活用并行计算
4. 内存不足
请教我关于"内存不足"！
现象: Out of Memory 错误
前辈说"收敛失败一定要好好做"的含义明白了。
对策:
- 使用核外求解法
- 削减网格规模
- 确认64bit版求解器
- 增加内存配置
哇，收敛失败的话题，相当有意思！请再多讲讲。
Nastran代表错误
代表错误具体是什么意思？
Abaqus代表错误
请教我关于"代表错误"！
那么，把工具名称弄好了的话就大致没问题吧？
"解析不匹配"时的做法
- 首先深呼吸——不要慌张地随便改设置，问题会更复杂
- 制作最小再现情况——把抛物反射镜天线的问题简化到最简单形式。"减法调试"最有效
- 一次改一个——同时进行多个改变会搞不清哪个有效果。像科学实验一样做"对照实验"原则
- 回到物理——如果计算结果像"违反重力的东西浮起来"这样非物理，怀疑输入数据根本性错误
价值
更多
错误