隐身战机以超音速掠过雷达网,机身蒙皮却像一块精密的"智能滤网",友方频段的信号畅通无阻地穿透,而敌方雷达波则被悉数反射或吸收,让这一切成为可能的,正是频率选择表面(FSS)。
FSS本质这一特性,它正成为天线罩带通滤波、电磁隐身、反射阵天线等高端应用中不可替代的关键技术路径。然而当工作频率不断升高,FSS 的微小单元需要集成到飞机、舰船等大型平台上时,传统仿真手段便遭遇了严峻挑战:全模型网格规模呈指数级爆炸,再强的计算资源也难以承受,工程师们往往面临"算不出、算不快、算不准"的困境。 如何破局?
01 案例背景
本文以 INTESIM 高频电磁仿真模块为利器,系统展示了一套从FSS单元到有限大阵列的完整仿真解决方案:先用 Floquet 端口+周期边界精准提取单元特性,再通过 SBR 弹跳射线法结合 Fresnel 边界,将复杂阵列等效为均匀阻抗表面。这一方案不仅精度媲美国际主流商软,更让电大尺寸平台上的 FSS 隐身性能评估从"耗时难收敛"变为"高效可量化",为装备隐身设计提供了可靠的工程依据。
图1 FSS结构的应用场景
02 案例功能特点
所属物理场模块:高频电磁分析
主要验证功能:周期边界、Floquet端口、色散材料以及SBR的Fresnel边界
03 技术实现
在电磁仿真中,无限大周期阵列的严格求解依赖于对基本单元的精确建模,而电大尺寸有限阵列的系统级响应则需兼顾计算效率。本技术方案采用“微观单元精确提取到宏观阵列等效替代”的分层求解策略:
通过Floquet周期边界与端口技术,从单个FSS单元全波仿真中提取出包含入射角、频率及极化的全角域S参数矩阵;在电大平台仿真中,将S参数矩阵以Fresnel边界的形式赋予待分析的阵列区域,从而将复杂的周期结构等效为具有相同电磁响应的均匀阻抗表面。
在保证阵列边缘绕射和多次反射精度的前提下,大幅降低了网格数量与求解自由度,为电大尺寸FSS阵列的RCS快速分析提供了可行途径,具体实现过程如下。
(3.1)FSS单元案例
(3.1.1)几何模型
INTESIM软件自带几何建模功能,同时也支持外部几何模型和网格模型导入。本案例FSS单元几何模型采用外部导入几何的方式,导入的FSS模型如下图所示:
图2 FSS单元几何模型
频选表面模型由圆环金属贴片、介质基板构成,通过Floquet端口施加激励,完成高频电磁场分析。
(3.1.2)材料参数
本案例用到的材料为介质、空气和色散材料,金属部分均采用PEC边界来代替。如下表所示:
材料属性表1
(3.1.3)边界及激励
(1) 周期边界
以单元为中心建立空气盒子,在x、y两个水平方向施加主从边界条件(Master/Slave),使电磁场满足周期性拓展,模拟无限大阵列。
沿z向在空气盒子顶面和底面分别设置Floquet端口。端口自动将入射平面波展开为TE和TM的Floquet模式序列,求解后可直接获得S参数。单元的周期边界和Floquet端口设置如下图所示:
本案例用到的材料为介质、空气和色散材料,金属部分均采用PEC边界来代替。如下表所示:
图3 周期边界和Floquet端口激励设置
图4 编辑源界面
(3.1.4)网格部分
周期单元模型中,对金属图形边缘、介质交界面等场强突变区域进行网格加密,周期边界中主从边界上的网格节点一一对应。
图5 单元周期网格
(3.1.5)仿真结果
仿真得出的周期边界FSS单元S参数如下图所示:
图6 S参数随频率变化曲线结果与商软对比
FSS结构的S参数在2GHz附近S11曲线低点接近-29 dB。INTESIM的计算结果与商软结果高度吻合,两条曲线的形状、通带位置基本一致,表明求解器具有高精度。
电磁波入射到FSS结构时的电场强度分布:
图7 金属面电场Mag_E分布
电场幅值(Mag_E)分布可以看出,电磁波在FSS层内部及金属贴片边缘产生了明显的能量聚集,金属贴片边缘电场强度达到14504 V/m。
(3.2)FSS阵列案例
(3.2.1)几何模型
几何模型采用外部导入方式构建,其结构由频率选择表面(FSS)单元分别排列为10×10与100×100两种阵列。首先以10×10单元阵列及等效平板模型为对象,基于S参数反演法获取均匀平板的等效电磁参数,对等效前后性能对比验证等效方法的有效性与精度。再将提取得到的S参数矩阵作为Fresnel边界条件,赋给100×100电大尺寸阵列。在阵列上施加平面波激励,评估散射性能。
10×10 FSS阵列和等效平板模型:
图8 10×10单元阵列和等效平板模型
100×100 FSS阵列等效为Fresnel边界平板模型:
图9 100×100单元阵列等效为Fresnel边界
(3.2.2)材料参数
10×10FSS阵列材料参数见3.1.2详细说明。等效后的材料参数通过S参数反演法得到等效电磁参数,具体参数如下:
材料属性表2
针对100×100阵列将复杂的周期结构等效为具有相同电磁响应的均匀阻抗表面,将S参数矩阵以Fresnel边界的形式赋予待分析的阵列区域,无材料属性。
(3.2.3)边界及激励
(1)SBR中得Fresnel边界
在系统级电大平台仿真中,将局部区域设定为Fresnel边界。通过单元仿真计算得到的全角域S参数矩阵(随频率、入射角、极化变化的反射/透射系数),SBR算法采用Fresnel边界等效替代真实的周期物理模型,极大提高计算效率。
图10 Fresnel边界
(2)平面波激励
图11 平面波激励
(3.2.4)网格部分
SBR算法仿真中使用三角面元划分网格:
图12 三角面元网格
(3.2.5)仿真结果
10×10 FSS模型平面波垂直入射时的单站RCS结果:
图13 10×10FSS模型单站RCS随频率变化曲线与等效平板对比
INTESIM 计算得到10×10阵列FSS模型的单站RCS结果与等效平板模型基本一致,在2GHz处RCS值低至-28dB,在此频点FSS结构的隐身效果好,两条曲线的一致性较好也验证了基于S参数矩阵等效方法的有效性与高精度。
100×100 FSS模型平面波垂直入射下设置Fresnel边界的双站RCS结果:
图14 100×100平板模型双站RCS随角度变化曲线与商软对比
INTESIM 计算得到的电大尺寸100×100阵列FSS模型双站RCS结果曲线在俯仰角Theta=180°时RCS值大,表明在此角度FSS结构的隐身效果相对较差。INTESIM计算得到的双站RCS结果曲线与对标商软高度一致,验证了Fresnel边界计算的准确性和可靠性。
04 总结展望
案例展示了INTESIM高频电磁仿真中核心功能在FSS设计中的协同应用:周期边界与Floquet端口实现了无限大阵列的等效高效模拟,大幅缩短单元优化周期;SBR算法加Fresnel边界打通了从单元到电大平台的设计链路,让FSS的隐身、天线罩等系统级性能得到快速计算。随着5G/6G通信及先进隐身技术的不断发展,INTESIM将继续深耕高频电磁领域助力高端装备研发降本增效。