龙伯透镜天线这个题目,我在做多波束天线方案预研的时候反复翻出来过好几回。它最吸引人的地方很直接:把馈源沿着球面挪个位置,波束就跟着转,中间不需要移相器、不需要机械转动、也不需要有源器件,纯靠一块球形的空间变化介电材料把球面波“掰”成平面波。听起来像魔法,本质上就是梯度折射率介质在做光程补偿。真要把它落到CST 电磁仿真里,麻烦就来了——仿真工具默认一个实体配一种均匀材料,而龙伯透镜要求介电常数从球心的 2 连续降到球面的 1,这种“空间变量”怎么建模、怎么离散、离散多少层才算够、网格怎么划、结果怎么判断可信,全是坑。这篇内容就是我围绕这套流程整理的一份实操记录,把分层离散、空间相关材料定义、网格策略、结果判读、3D 打印落地这几块串起来讲透。刚接触梯度折射率天线的同学可以照着走一遍流程,做过类似结构的同行也可以对照看看我踩过的那些坑有没有重合。
1. 项目概述与核心需求拆解
1.1 龙伯透镜真正要解决的问题
多波束、宽角扫描这两个需求,在天线圈里一直是“又贵又难”的典型。相控阵能打,但通道数一多,T/R 组件的成本就压不住;反射面天线便宜,可惜扫描得靠机械转动,响应速度上不去,转动机构本身也占体积。龙伯透镜走的是第三条路:一块球对称的梯度折射率介质球,球面上任意一点既是“焦点”也是“馈电点”。把馈源贴在球面某处,球对面的口径面就会辐射出一个笔形波束;馈源沿球面滑动,波束就在空间里扫描,中间没有任何有源器件。
这个特性来自它的折射率分布:球心折射率最高,越往外越低,到球面正好等于自由空间的折射率 1。球面波从馈源出发,在球内被连续弯折,等相位面在出射口被拉平,等效成一个均匀照射的口径。由于球面处折射率与空气匹配,理论上没有反射,这也是它比介质棒透镜、介质球透镜更“干净”的根本原因。工程上常见的做法是把球切成两半,用半球加一块金属地板,或者整球加一个馈源支撑结构,前者省材料、后者对称性更好。
1.2 空间变化介电材料:难点到底在哪
如果这是一块均匀介质,事情十分钟就做完了:画个球、赋个材料、加个波导端口、跑。但“空间变化”这四个字把问题彻底改了性质。材料的相对介电常数 εr 是位置的函数,球内每一点的电磁参数都不同,这意味着:
- 几何上,它不再是一个“实体+材料”的组合,而是一个参数场;
- 网格上,材料梯度方向必须被足够细的网格捕捉,否则阶梯误差会直接吃掉增益;
- 材料定义上,多数仿真平台的默认材料模型是常数,要表达连续变化,得绕道走。
三条主流路线摆在这里:分层离散,把球切成若干个同心球壳,每层赋一个常数 εr;空间相关材料,直接用随坐标变化的表达式或外部数据表定义 εr;体素/网格导入,把外部生成的三维材料分布整体灌进仿真。每一条都有明确的适用场景和代价,选错了后面全是返工。
1.3 仿真目标与验收指标
在动手建模之前,我习惯先把验收指标写死,不然跑完一堆远场图根本不知道好坏。这个案例我设的中心频率 10 GHz,球半径 100 mm,指标如下表。
| 指标项 | 目标值 | 说明 |
|---|---|---|
| 中心频率 | 10 GHz | 对应自由空间波长 30 mm |
| 口径直径 | 200 mm | 约 6.7 个波长 |
| 峰值增益 | 不低于 24 dBi | 理论方向性约 26.4 dBi,口径效率按 60% 估 |
| 半功率波束宽度 | 9 度到 12 度 | 粗略估算 65 乘以波长除以口径直径 |
| 输入回波损耗 | 优于 -10 dB | 馈源与透镜之间的匹配 |
| 副瓣电平 | 低于 -15 dB | 与口径场锥削和边缘衍射相关 |
| 正负 60 度扫描增益跌落 | 小于 2 dB | 验证球对称性是否被破坏 |
| 口径效率 | 高于 55% | 综合判断分层与网格是否到位 |
提示:指标里最容易忽略的是“扫描增益跌落”。很多人只盯着正对时的增益,结果一转到 60 度发现掉了 5 dB,回头查半天才发现是分层太粗导致斜入射时阶梯散射加剧。这一项要一开始就写进验收表。
2. 原理与建模思路:为什么这样选型
2.1 折射率分布公式的来历与用法
龙伯透镜的经典分布是 n(r) 等于根号下 2 减去 (r/R) 的平方,其中 r 是到球心的距离,R 是球半径。这个式子不是拍脑袋来的,它来自费马原理加球对称假设:要求从球面上一点发出的所有射线,汇聚到球面对径点时走过的光程完全相等,解这个积分方程就得到上面的形式。换成电磁参数,磁导率取 1,相对介电常数就是 n 的平方,也就是 2 减去 (r/R) 的平方。
几个关键数值先记住:球心 εr 等于 2,球面 εr 等于 1,中间单调递减。球面处等于 1 这一点极其重要,它保证了波从介质进入空气时没有阻抗突变,反射系数理论上为零。这也是为什么龙伯透镜的端口回波可以做得很低,前提是馈源本身匹配好。
有一点容易搞混:龙伯透镜的焦点在球面上,不在球心,也不在球外。所以馈源要贴着球面放,相位中心对准球面点。我见过有人把喇叭口放在离球面半波长的地方,结果口径相位弯曲,副瓣直接抬高 6 dB,还以为是分层数据错了。
2.2 连续渐变、分层离散、空间相关材料:三条路怎么选
这三条路线我在不同阶段都试过,结论是各有各的活路,但优先级明确。
| 路线 | 建模难度 | 物理保真度 | 网格负担 | 适用阶段 |
|---|---|---|---|---|
| 同心球壳分层 | 低,宏批量生成 | 中,阶梯近似 | 中,层数越多越大 | 方案验证、参数扫描 |
| 空间相关材料 | 中,需定义表达式或数据表 | 高,连续分布 | 高,梯度区必须加密 | 最终确认、精细分析 |
| 体素或外部模型导入 | 高,需要预处理 | 最高,能含缺陷 | 很高,文件大 | 实物对标、3D 打印件复核 |
分层的核心优势是简单直接,每一层就是一个普通实体,任何求解器都吃得下,布尔运算出错的概率也低。缺点是阶梯近似会引入寄生散射,层与层之间的界面相当于一圈圈小的阻抗台阶,斜入射时尤其明显。
空间相关材料是“物理上更对”的做法,材料的 εr 随坐标连续变化,理论上没有台阶。代价是网格:梯度材料要求网格在材料变化方向上足够细,否则求解器采样不到梯度,效果和分层没什么区别,还白白多花了建模时间。
体素导入最适合实物复现。3D 打印件一定有填充率误差、层间缝隙、支撑痕迹,把这些数字化后灌进仿真,能提前预判实物和仿真的差距。但文件动辄几百兆,几何清理非常耗时间,一般放在项目后期做。
我的建议:先用分层跑通全流程、把指标摸清楚,再用空间相关材料做一次交叉验证,两者增益差在 0.3 dB 以内就说明分层够用了,可以放心往下走。
2.3 分层数量到底取多少:给一个能算的数
层数不是拍脑袋定的,有一个够用的工程估算。设总层数 N,每层厚度 t 等于 R 除以 N。相邻层之间的折射率差,在最外层最大,约为 t 除以 R,也就是 1 除以 N。单层引入的相位误差近似等于波数乘以折射率差再乘以层厚,整理后得到:
最大相位误差 ≈ k0 乘以 R 除以 N 的平方
代入本例:k0 等于 2π 除以 30 mm,约 0.2094 每毫米;R 等于 100 mm,所以 k0 乘以 R 约等于 20.94 弧度。取常用判据,把最大相位误差控制在 π/8(约 0.393 弧度)以内,解得 N 大于等于 7.3,也就是 8 层起步。
| 层数 N | 最大相位误差(弧度) | 折算角度 | 评价 |
|---|---|---|---|
| 5 | 0.838 | 48 度 | 完全不够,波束会散 |
| 8 | 0.327 | 18.7 度 | 勉强,副瓣偏高 |
| 10 | 0.209 | 12.0 度 | 可用,适合快速验证 |
| 16 | 0.082 | 4.7 度 | 比较稳 |
| 20 | 0.052 | 3.0 度 | 推荐,增益接近收敛 |
| 30 | 0.023 | 1.3 度 | 收益递减,网格压力大 |
这个估算是“最坏情况上界”,实际射线在球内是弯曲的,穿过每层时的入射角都不同,误差不会简单线性累加,所以 10 到 20 层在实际项目里已经够用。但如果你要做 60 度斜扫描,建议往上取,我一般用 20 层。
2.4 层内介电常数怎么赋值:取中点还是做积分平均
这是个小细节,但影响不小。最朴素的做法是取每层中心半径处的 εr。层数多时没问题,层数少时会有系统性偏差。更稳的做法是对每层做径向积分平均:
ε_层 = 2 减去 (r1 平方 + r1 乘 r2 + r2 平方) 除以 (3 乘 R 平方)
其中 r1、r2 是该层的内外半径。这个式子其实就是把 2 减去 (r/R) 平方沿层厚积分再除以层厚,结果是一阶精确的,能抵消掉取中点带来的偏差。层数少的时候(比如 8 层),两种做法的增益差能到 0.2 dB 量级。既然公式就一行,没必要省这个事。
2.5 坐标系与馈源位置必须先约定清楚
我吃过这个亏:模型建完跑出远场图,波束指向和预期差了 90 度,排查半天发现是馈源放在 +x 轴、而远场图默认的观察坐标系是另一套。后来养成习惯,建模第一步就把约定写进参数表:
- 球心固定在原点;
- 馈源放在 +x 轴上,相位中心落在 (R, 0, 0);
- 主波束沿 -x 方向出射;
- 远场监视器的 theta 从 0 到 180 度、phi 从 0 到 360 度;
- 扫描验证时,馈源绕 z 轴旋转,等效波束在 xoy 面内扫描。
先把这套东西写死,后面看方向图切片、比较不同扫描角的时候才不会乱。
3. 建模实操:从几何到材料赋值
3.1 参数化变量表怎么搭
模型一定要参数化,尤其层数 N,后面要做收敛性扫描,改一个数字整个模型重建,非常省事。我用的参数清单大致如下。
| 参数名 | 含义 | 取值 |
|---|---|---|
| R | 球半径 | 100 mm |
| N | 分层数 | 20 |
| f0 | 中心频率 | 10 GHz |
| eps_max | 球心介电常数 | 2.0 |
| eps_min | 球面介电常数 | 1.0 |
| tand | 损耗角正切 | 0.005 |
| gap | 馈源与球面间距 | 0 mm |
| pad | 背景空间 padding | 0.25 倍波长 |
把这些放进参数列表之后,球壳生成、材料赋值、边界设置全部引用变量,改层数只需要改 N 一个数字。
3.2 同心球壳的批量生成
二十层手动做布尔运算,做三遍就烦了,肯定要写宏。思路很简单:外层球减去内层球,得到一层壳。但更省事的办法是用抽壳命令,先生成一个半径为 r2 的实心球,抽掉厚度为 r2 减 r1 的壳,直接得到球壳,然后删掉辅助体。
下面这段是思路示例,语法细节以你手上的版本为准:
Dim R As Double, N As Integer, i As Integer Dim rIn As Double, rOut As Double, t As Double Dim epsL As Double, nm As String R = 100 N = 20 For i = 0 To N - 1 rOut = R * (1 - i / N) rIn = R * (1 - (i + 1) / N) t = rOut - rIn ' 用积分平均公式计算该层介电常数 epsL = 2 - (rIn^2 + rIn*rOut + rOut^2) / (3 * R^2) nm = "Shell_" & CStr(i) ' 生成实心球 Sphere.Reset Sphere.Name "tmp_" & CStr(i) Sphere.Radius rOut Sphere.Create ' 抽壳成球壳 Solid.Shell "tmp_" & CStr(i), nm, t ' 赋材料 ' 这里按 epsL 新建或调用预定义材料 ' 清理辅助体 Next i有几点实操经验值得单独说。第一,循环顺序要从外往内还是从内往外,取决于你用的命令是否会产生重叠实体,抽壳法两种都行,布尔差法建议从外往内。第二,材料命名要有规律,比如用介质加层号,方便后面批量改参数。第三,生成完记得做一次相交检查,二十层里只要有一层的法向反了,远场就会出现莫名其妙的裂瓣,查起来很痛苦。
注意:抽壳厚度不能小于网格步长的量级。如果你的最外层厚度只有 0.5 mm,而全局网格步长是 2 mm,这一层等于被网格“吃掉”了,等于没做。层厚和网格的比值至少要 1.5 以上。
3.3 用空间相关材料替代分层
如果版本支持空间相关材料,路径大致是:新建材料,类型选择随空间变化的那一档,然后给相对介电常数、损耗角正切各写一个关于 x、y、z 的表达式,或者挂接外部数据表。表达式直接写 2 减去 (x 平方 + y 平方 + z 平方) 除以 R 平方,外面套一个条件,半径超过 R 的区域置 1。
这个做法在频域求解器里效果不错,因为频域求解器允许网格和材料梯度共存。时域求解器也能跑,但要求材料变化在一个网格单元内不能太剧烈,否则会引入数值反射。实际经验是:梯度区的网格步长要控制在十分之一波长以内,否则精度和分层差不多。
外部数据表的方式更通用,格式就是每行一个点,列是坐标和电磁参数:
# x(mm) y(mm) z(mm) eps_r tan_d 0.000 0.000 0.000 2.000 0.005 1.000 0.000 0.000 1.999 0.005 2.000 0.000 0.000 1.998 0.005数据量要提前算:半径 100 mm、间距 1 mm 的笛卡尔网格,包住球的体积大概是 200 的三次方,八百万个点,纯文本几百兆。实际用的时候没必要这么密,间距 2 mm 就够,或者干脆用柱坐标/球坐标采样,点数能少一个量级。
3.4 损耗和色散怎么处理
真实材料不可能是零损耗。3D 打印用的光敏树脂,损耗角正切一般在 0.005 到 0.02 之间,工程塑料略高一些。这个量级在 10 GHz 上会让增益掉多少?粗略估一下:单程光程约 200 mm,衰减常数近似正比于频率乘以损耗角正切再乘以根号下介电常数,折算下来单程衰减在 0.3 到 0.8 dB 之间。也就是说,一个仿真里增益 26 dBi 的理想透镜,换成有损耗材料可能只剩 25.2 dBi。
所以材料定义里一定要带损耗,别图省事填零。仿真和实测对不上,十有八九是损耗这块低估了。
色散方面,大部分打印材料在 8 到 12 GHz 内介电常数变化小于 2%,第一轮仿真可以按常数处理。如果要做宽带分析,最好用波导法或者谐振腔法先测几个频点的 εr,再拟合成一条曲线挂进材料模型。
4. 求解设置与网格策略
4.1 时域还是频域:看你要干什么
这两个求解器我一般是配合着用。时域求解器一次激励能拿到整个频段的 S 参数和远场,适合快速扫 8 到 12 GHz 看带宽趋势,代价是网格必须覆盖整个频段,最苛刻的是最高频,网格量偏大。频域求解器一次只算一个频点,但单点精度高,尤其适合验证口径相位这种细节。
我的流程是:先用频域求解器在 10 GHz 单点跑一遍,确认几何、材料、馈源都没错,看近场相位平不平;确认没问题之后,切到时域求解器做宽带扫描。这样能省掉大量“跑了两小时发现模型建错”的时间。
还有一点:分层模型的层间界面是真实存在的几何面,网格必须与这些面共形,否则阶梯散射会被放大。时域求解器的六面体网格在这种情况下会自动加密,网格量会比均匀球大不少,内存要提前留够。
4.2 边界条件、背景和对称性
透镜天线是有限尺寸辐射结构,边界用开放边界加背景空间,padding 至少四分之一波长,我一般给到 0.3 倍波长,避免近场被截断。远场监视器要提前加上,位置就放在开放边界内表面上。
对称性这块要小心。球本身是球对称的,但馈源一放上去,对称性就只剩一个面了。如果馈源是标准的矩形波导口径,且极化沿 z 方向,那么 xoy 平面是一个磁对称面,可以用磁壁切一半,网格量直接减半。但如果馈源带支撑结构或者偏置,对称性被破坏,老老实实算全模型。
提示:用了对称面之后,远场图的 phi 覆盖范围会变成 180 度,别以为是求解出错,这是正常的。
4.3 网格:这个案例最容易翻车的地方
先算一个数:球内最中心的介质 εr 等于 2,对应介质波长约 21 mm。按每波长 10 到 12 条网格线算,网格步长应该在 1.8 到 2.1 mm。而层厚是 R 除以 N,20 层时是 5 mm,一层里面只有两三个网格单元,勉强够。如果只做 10 层,层厚 10 mm,一层里五个网格单元,看着够,但阶梯误差本身就大了。
麻烦在于网格量。球的体积是 4/3 乘以 π 乘以 R 的三次方,约 4.19 乘以 10 的 6 次方立方毫米。按 2 mm 步长,每个单元体积 8 立方毫米,理论上五万多个单元,加上海量曲面网格和层间界面加密,实际单元数会到几十万到上百万。这个量级在普通工作站上还能跑,但再加密一档就吃力了。
| 网格步长 | 每波长网格数(介质内) | 估算单元量级 | 内存需求 | 建议用途 |
|---|---|---|---|---|
| 5 mm | 4 | 十万级 | 8 GB 以内 | 快速试算,只看趋势 |
| 3 mm | 7 | 数十万 | 16 GB | 常规验证 |
| 2 mm | 10 | 百万级 | 32 GB | 精细分析 |
| 1.5 mm | 14 | 两百万以上 | 64 GB | 收敛性检查 |
实操上我一般做两级网格:全局 3 mm,透镜内部加一个局部网格区,加密到 2 mm。这样既保证球内介质梯度被捕捉,又不至于让整个背景空间也跟着加密。
曲面逼近精度也要调。球面在六面体网格下是阶梯化的,默认精度下球面会有一圈圈锯齿,这些锯齿本身就是散射源。把曲面逼近精度调高一到两档,球面会平滑很多,远场副瓣能降 1 到 2 dB。
4.4 激励设置与监视器布局
馈源我用的是标准矩形波导口径,工作在 10 GHz,端口用波导端口激励。馈源口径面贴着球面,相位中心尽量对准球面上的焦点位置。如果馈源有喇叭段,喇叭口面到球面的距离要给准,这个距离对口径相位影响很大。
监视器布局如下:
- 远场监视器设在 10 GHz,时域求解时可以直接拿到该频点的方向图;
- 近场监视器设一个平面,位置在球外 x 等于 R 加 10 mm 处,用来观察口径场相位;
- 另一个近场监视器设在 xoy 平面,观察球内射线弯曲情况;
- 端口监视器记录 S11。
近场相位这个监视器特别有用。把口径面上的电场相位画出来,如果相位起伏在正负 15 度以内,说明透镜工作正常;如果看到明显的同心圆环状相位起伏,那基本就是分层太粗,或者层内介电常数赋值有系统性偏差。
5. 结果判读与验证
5.1 先看近场,再看远场
很多人一上来就看远场方向图,其实近场信息量大得多。我的判读顺序是:先看球内 xoy 平面的电场幅度分布,正常情况应该能看到射线从馈源出发、在球内逐渐弯折、最后在出射面附近形成近似平行的等相位面。如果射线在球内就走成了直线,说明介电常数根本没起作用,八成是材料没赋上或者赋成了 1。
再看口径面的相位分布。理想情况下,出射面的相位应该接近常数,允许有轻微的锥削。实测中常见的两种异常:相位呈中心对称的抛物面形状,说明等效焦距不对,多半是折射率分布公式里的半径用错了;相位呈环状起伏,那是分层阶梯的痕迹,层数加一加就能改善。
5.2 远场指标怎么和理论对
理论方向性用 4π 乘以口径面积除以波长平方,口径面积是 π 乘以 0.1 的平方,约 0.0314 平方米,波长平方是 0.0009,算下来约 438,换算成 dBi 是 26.4。乘上口径效率,60% 对应 24.2 dBi。
仿真跑出来如果只有 22 dBi,差了 2 dB 多,就要找原因了。常见的原因按概率排序:分层太少、网格太粗、馈源遮挡、馈源与球面之间有间隙、材料损耗设太高。逐个排查,一般能找回来。
波束宽度对照 65 乘以波长除以口径直径,约 9.8 度。仿真值在 9 到 12 度之间都算正常,偏宽说明口径场锥削太厉害,偏窄说明边缘照射过强,副瓣会高。
副瓣这块要特别注意。理想龙伯透镜的口径场接近均匀分布,均匀分布的第一副瓣理论值是 -13.2 dB。实际能做到 -15 dB 以下,靠的是馈源方向图本身带来的自然锥削。如果你的仿真副瓣只有 -10 dB,那不是透镜的问题,是馈源照射太宽或者边缘衍射太强。
5.3 扫描特性怎么验证
扫描验证是这个案例的重点。做法是把馈源绕球心旋转,分别在 0 度、30 度、60 度三个位置各跑一次,看主波束指向是否跟着转、增益跌落多少。
理想情况下,球对称结构在任意扫描角下增益应该完全一致。实际会有跌落,来源有三个:馈源本身在不同指向上的方向图不对称、馈源支撑结构的遮挡、分层阶梯在斜入射时的散射加剧。前两个是结构问题,第三个是建模精度问题。
| 扫描角 | 仿真增益 | 相对 0 度跌落 | 判断 |
|---|---|---|---|
| 0 度 | 24.6 dBi | 基准 | 正常 |
| 30 度 | 24.4 dBi | 0.2 dB | 优秀 |
| 60 度 | 23.9 dBi | 0.7 dB | 可接受 |
| 75 度 | 22.8 dBi | 1.8 dB | 边缘遮挡开始显现 |
如果 60 度就掉了 2 dB 以上,先查分层数够不够,再查馈源在斜位置时有没有被球的其它部分挡住。有些设计里馈源支架会伸到球外,扫描到一定角度就挡在波束路径上,这个在仿真里很容易看出来。
5.4 收敛性和可信度检查
仿真跑完不算完,得证明这个结果可信。我固定做三件事。
第一件,网格收敛。把全局网格从 3 mm 加密到 2.2 mm,增益变化小于 0.2 dB、副瓣变化小于 1 dB,说明网格够了。如果加密后增益涨了 1 dB,那之前的结果就是不可信的。
第二件,分层收敛。把层数从 20 加到 30,看增益和副瓣的变化。变化小于 0.15 dB,说明 20 层已经收敛。
第三件,能量平衡。看端口吸收加上辐射功率是否等于输入功率,差值应该在百分之几以内。如果能量对不上,多半是有源区没包住或者材料设置有问题。
这三件事做完,结果才能拿去写报告。我见过太多人拿一个没做收敛性检查的结果去指导加工,最后实测差了 3 dB,返工成本比仿真时间高几十倍。
6. 常见问题与排查技巧实录
6.1 问题速查表
| 现象 | 最可能原因 | 排查动作 |
|---|---|---|
| 增益比理论低 2 dB 以上 | 分层太少或网格太粗 | 层数加倍、网格加密各跑一次 |
| 副瓣异常高,且有规律裂瓣 | 层间界面阶梯散射 | 层数加倍,或改空间相关材料 |
| 主波束指向偏离预期 | 馈源位置或坐标系约定错 | 检查馈源坐标和远场 theta 定义 |
| 端口回波很差 | 馈源与球面间隙不当 | 微调间隙,观察 S11 变化趋势 |
| 口径相位呈环状起伏 | 层内介电常数赋值偏差 | 改用积分平均公式赋值 |
| 远场图完全对称但增益极低 | 材料没赋上,等效于空气球 | 检查材料列表和实体材料绑定 |
| 求解器报网格过多 | 全局网格过密 | 改局部网格,背景空间放粗 |
| 扫描时增益突然掉很多 | 馈源支撑结构遮挡 | 隐藏支撑件单独跑一次对比 |
6.2 几个反直觉的坑
第一个坑:层数加多不一定更好。有一次我把层数从 20 加到 40,增益反而掉了 0.3 dB。原因是层厚变成 2.5 mm,而网格步长是 3 mm,每层不到一个网格单元,求解器根本分辨不出层间界面,等效成了一团模糊介质。这个教训是:层厚必须显著大于网格步长,两者要配套调,不能单方面加层。
第二个坑:球面曲面精度比网格密度更影响副瓣。同样是 3 mm 网格,曲面逼近精度调高一档,副瓣能改善 1.5 dB,而单纯加密网格只改善 0.5 dB。原因是球面锯齿本身就是散射源,把面拟合好,比把体分细更划算。
第三个坑:馈源口径不要贴着球面画。理论上馈源相位中心在球面上,但馈源本身有物理尺寸,口径面如果和球面重合,几何上会重叠,布尔运算容易出错,网格也会在交界处畸变。我的做法是把馈源口径设在球面外 0.5 mm,让相位中心仍落在球面上,几何上留一点缝,误差可以忽略,建模却干净很多。
第四个坑:损耗角正切别拍脑袋。我一开始图省事填了 0.02,结果增益比实测低了 1.2 dB。后来测了实际材料,是 0.006,重新跑一遍就对上了。材料参数这件事,能测就测,测不了就查供应商数据,别猜。
还有一个容易被忽略的点:分层球壳之间的“缝”。如果宏生成的球壳之间留有微米级间隙,仿真里会被当成空气薄层,形成额外的阻抗台阶。生成完做一次间隙检查,把相邻层的半径严格对齐,能省掉后面一堆莫名其妙的反射。
7. 从仿真到实物:3D 打印的等效介质换算
7.1 用填充率做出连续介电常数
仿真里可以随手写一个连续变化的 εr,实物做不出来。3D 打印能控制的只有填充率,也就是材料实体占总体积的比例。要把填充率换算成等效介电常数,常用的是混合介质等效模型,形式大致是:
(εeff 减 1) 除以 (εeff 加 2) 等于 填充率 乘以 (εp 减 1) 除以 (εp 加 2)
其中 εp 是打印材料本身的介电常数。反解出填充率,就能知道每一层该打多密。以 εp 等于 2.7 的树脂为例,反算结果如下表。
| 目标 εeff | 需要的填充率 | 实现难度 |
|---|---|---|
| 2.00 | 69% | 中,接近实心 |
| 1.80 | 55% | 中 |
| 1.50 | 40% | 容易 |
| 1.20 | 17% | 容易,但结构强度差 |
| 1.00 | 0% | 纯空气,需支撑 |
可以看到,球心附近要求 69% 的填充率,接近实心;球面附近接近纯空气。这个梯度对打印工艺是个挑战,尤其是外层填充率低的地方,结构强度会很差,容易塌陷。实践中常见两种处理:一是外面加一层薄的实心蒙皮做支撑,二是把球切成上下两半分别打印再粘合。
外加固蒙皮会引入一个 εr 偏高的薄壳,对电性能有影响。仿真里可以把这层壳也建出来,厚度 1 mm 左右,εr 按实心材料设。这样仿真和实物的一致性会好很多。
7.2 打印误差怎么折算成电性能偏差
打印精度一般在 0.2 mm 量级,层厚 0.1 mm。对 10 GHz 来说,0.2 mm 是波长的百分之一不到,单看几何误差不大。但填充率误差会被放大:填充率偏差 5%,等效 εr 可能偏差 3% 到 4%,累积到整个口径上,相位误差就不容忽视了。
粗略估一下:εr 偏差 3%,折射率偏差约 1.5%,穿过 200 mm 光程,相位偏差约 0.03 乘以 200 再乘以波数,折算下来十几度。听起来不大,但如果偏差是系统性的(比如整体填充率都偏低),相位面会整体倾斜,波束指向就会偏。
所以实物做完一定要测。最简单的测法是远场方向图,看增益和波束指向;讲究一点可以切开一个截面,用近场扫描看口径相位。测出来的数据和仿真对不上,先别急着改设计,把实测 εr 代回仿真重跑一遍,通常就吻合了。
7.3 装配环节的几个细节
半球拼装的时候,粘接剂的介电常数往往比空气高不少,胶层厚度控制不好会在赤道面上形成一圈高介电常数带,等效于在口径中间插了一道相位台阶。我的做法是用低介电常数的胶,涂薄薄一层,压紧后把溢出的刮干净。仿真里可以在赤道面加一个 0.2 mm 厚、εr 为 3 的薄圆环,看看影响有多大,心里有个数。
馈源安装也要注意。馈源要固定在球面上,支架材料最好是低介电常数、低损耗的,金属支架会直接挡住一部分口径。有些设计把馈源做成可滑动的小车,沿着球面导轨移动来实现扫描,导轨本身的散射要在仿真里评估一下。
8. 一点个人体会
这套流程我从头到尾跑过好几遍,最大的感受是:龙伯透镜的仿真难点不在电磁理论,而在“连续”和“离散”之间的转换。理论公式就那一行,谁都会写;真正决定结果好坏的,是分层取多少、层内怎么赋值、网格和层厚怎么配合、曲面精度调多高。这几个参数之间是耦合的,单独调一个往往没用。
我现在固定的做法是:先用 10 层、3 mm 网格快速跑一版,看整体趋势和指标量级;确认没问题后升到 20 层、2 mm 网格做正式计算;最后用空间相关材料跑一版做交叉验证。三版结果的增益差控制在 0.3 dB 以内,我才认为这个模型是可信的。这套流程看起来啰嗦,但比“一版跑到底、结果不对再回头查”要快得多。
另外就是别迷信仿真。分层模型再细,也做不出打印件的填充率误差和装配缝隙。仿真给出的 24.6 dBi,实物能做到 23.5 dBi 就已经很不错了。心里有这条线,验收的时候才不会自己跟自己较劲。