2026/9/7 11:22:43
FEATURED · 精选文章 
简介一份关于相位梯度超表面结构操控电磁波束的学术研究文献面向电磁超材料、微波通信与雷达技术方向的研究人员、工程师及高校师生。内容从高阻抗表面单元尺寸与反射相位的关系入手设计了一种相位梯度超表面结构并通过数值仿真验证其在10 GHz频段可实现约22°反射波束偏转体现了厚度薄、重量轻且可灵活操控波前相位的特点。文献还系统梳理了相位梯度超表面的代表性进展如广义斯涅尔定律、涡旋波束生成、异常反射和表面波耦合等为读者理解超表面波束调控原理和分析方法提供了完整脉络。该论文刊载于《浙江科技学院学报》2016年第28卷第3期作者单位涵盖同济大学与浙江科技学院结构包含摘要、关键词、理论分析和数值仿真等内容。资源为1个PDF文件大小约2.04MB已有219人浏览学习适合作为相关课题的参考文献或专业学习材料。 很多人刚接触超表面时第一反应都是先跑仿真在 HFSS 或者 CST 里画个单元扫几个参数然后急着排阵列。但折腾几轮之后你会慢慢意识到相位梯度超表面对电磁波束的操控本质上是一个“相位工程”问题。你前期把相位算得越明白后续的仿真、加工、实测就越顺。这篇文章我结合自己做相位梯度超表面的实际项目经历把从原理到设计、再到仿真和实测的完整链路拆开讲。内容涉及广义斯涅尔定律的工程理解、单元库搭建、波束偏折/聚焦/涡旋的设计方法以及加工和测试里那些文献不会写、但非常影响结果的坑。适合刚入门超表面的研究生也适合想把固定结构往可调方向延伸的工程师参考。1. 从相位不连续说起超表面为什么能实现波束操控1.1 广义斯涅尔定律超表面的底层逻辑传统的透镜、棱镜操控电磁波靠的是电磁波在介质中传播时积累的连续相位。波穿过一块透镜中心路径和边缘路径走过的光学长度不同出射面上就形成一个相位分布波前因此被弯曲。这套逻辑很成熟缺点是器件体积大、重量重、剖面高。超表面的思路完全不一样。它不需要波在体内走很长的距离而是在一个极薄的界面上让每个单元对电磁波施加一个突变相位。这个相位突变沿界面的变化率就是常说的“相位梯度”。2011年Capasso组的广义斯涅尔定律给出了定量描述Bsinθt - nisinθi (λ0 / 2π) · dΦ/dx这个公式看着长其实物理图像很清晰入射波和透射波在界面切线方向的波矢分量之差由界面的相位梯度来补偿。右边那一项 dΦ/dx 就是相位沿界面的变化率它的存在意味着即使垂直入射θi0透射波也可以被“推”向一个非零的出射角。换句话说相位梯度相当于给透射波一个额外的横向动量。我在项目里习惯把这条定律反过来用先根据目标偏折角算需要的相位梯度再反推每个单元应该落在这个梯度曲线的哪个位置。这个过程很像给波前画刻度尺每格分多少度相位排多少个单元计算清楚了再动手画结构。1.2 与相控阵的类比一块“超薄的移相阵面”理解超表面最直观的方式其实是把它想成一块无源的相控阵天线。相控阵靠每个通道的移相器给单元馈电通过控制通道间的相位差把主波束指向某个方向。超表面做的事情在数学上完全一致只是它的“移相器”变成了每一个亚波长单元结构而控制手段不是给电信号而是通过改变单元的几何尺寸、形状或者朝向来实现不同的相位响应。这个类比在工程上很有用。阵列天线里的副瓣电平、栅瓣抑制、单元间距问题都可以平移过来指导超表面设计。比如单元周期必须小于工作波长目的就是避免出现栅瓣相位量化级数太少会产生寄生波束这和相控阵里数字移相器的量化误差是同一个物理问题。我做第一个超表面样件时因为当时手里有一块标准的8×8微带阵列所以我对“相位差决定波束指向”这件事特别熟悉。后来把这种思维带到超表面上后很多概念一下子对上了。可以这么说超表面就是一张把移相网络直接“长”在辐射面上的阵面只是每个单元的“馈电”是自由空间里的入射波。2. 相位怎么“铺”到表面上单元结构与相位覆盖设计2.1 单元怎么选从H形到方环的取舍超表面单元的选择决定了你能获得的相位范围、幅度稳定性以及加工容差。目前最常用的单元可以归成几类贴片型矩形贴片、圆贴片、缝隙型H形缝隙、十字缝隙、环型方环、圆环、双开口环以及介质型介质柱等。项目里最常用的还是金属贴片/缝隙类因为它们可以直接用PCB工艺做出来。这些单元各有各的脾气。矩形贴片加工最简单、设计变量少但单个贴片的相位覆盖范围有限想凑出完整0-2π往往要依赖多个变量的组合H形缝隙可以通过改变两个分支的长度来实现较大的相位变化但对极化方向敏感排布时要注意入射波的极化方向一致性方环单元的相位曲线比较平滑而且通过外环尺寸和槽宽两个自由度能较好地兼顾相位范围和幅度稳定。我自己在10GHz频段做透射型超表面时最常用的是方环和十字缝隙的组合。介质基板选的是F4B相对介电常数2.65损耗角正切0.001左右厚度1mm。这个组合的好处是相位变化范围大而且F4B这种材料在C波段到X波段的介电常数比较稳定价格也适中不像Rogers板材那么贵。下面这个表格是我在实际项目中对比过的几种单元方案方便你根据需求快速做取舍单元类型相位覆盖能力透射/反射效率加工难度带宽表现适用场景矩形贴片中等透射一般低窄快速验证H形缝隙较大透射较高中中等线极化透射阵方环大较高中中等宽带偏折/聚焦方环互补结构很大高中高较宽需要大相位梯度的场景介质柱大高高宽高频/毫米波2.2 参数扫描的核心既要覆盖2π也要幅度均匀选定单元类型后下一步就是对单元参数做扫描。扫描的目标有两个第一找出能覆盖0-2π的相位变化范围第二所有工作点的透射/反射幅度尽量高且均匀。很多新手做单元仿真时就盯着相位曲线看到相位扫出360°就觉得万事大吉。但实际排完阵列后一做仿真方向图增益很低一查原因幅度响应在单元之间差了5dB甚至更多。幅度不均匀意味着阵列的口径效率直接受损因为远场叠加时各个单元“出力”不一致。我在做单元库时一般按以下步骤进行在CST或HFSS里建立单个周期单元模型边界条件使用主从边界unit cell端口设置为Floquet端口。固定其他参数不变扫描一个关键几何参数。比如方环结构固定环宽w不变扫描外环边长L从4.0mm到7.2mm步长0.2mm。提取每个样本的S21幅度和相位保存成工程文件。筛选出幅度大于阈值我一般取0.85以上的样本点再检查这些点对应的相位是否覆盖了完整的2π。如果相位覆盖不够增加第二个扫描变量。比如同时调槽宽或者加入介质层的厚度变化用二维扫描来拓宽相位范围。这里有个关键细节单一参数扫描的相位曲线往往是S形或者单调但范围不足。比如方环的谐振会影响相位在谐振点附近相位变化非常剧烈但幅度掉得也厉害。太过靠近谐振点的工作状态虽然能拿到必要的相位值但损耗会成为问题。所以最终单元库里每个入选的结构都要尽量避开谐振点太近的位置。2.3 透射还是反射两种工作模式的分水岭决定超表面是透射型还是反射型是整个项目最开始就要确定的事因为它影响后续所有的单元设计、馈电方式和测试方案。透射型的优势是馈电方便喇叭天线从法向照射透射波在另一侧直接形成需要的波前测试时收发天线分居样品两侧比较好布置而且可以做成类似传统透镜的器件。劣势也很明显透射型单元要在“匹配”和“相位变化”之间做折衷很难在获得大相位调节范围的同时保持很高的透射率。尤其是介质板两侧都有金属结构时谐振带来的插损往往有1-2dB。反射型的优势是效率高通常能拿到90%以上的反射率而且相位覆盖范围大劣势是需要考虑馈源遮挡双站测试时收发喇叭要分处同侧而且反射型超表面往往需要加金属背板剖面会略高一些。做偏折器或者聚焦透镜两边都能做但如果效率是首要指标我会选反射型如果项目要求空间馈电方便、需要直接观测透射波前那就只能做透射型。好在绝大多数发论文的场景透射型已经足够说明问题所以我后面的举例也以透射型为主。3. 从单元到阵列偏折、聚焦、涡旋三种典型的波束操控3.1 波束偏折一个可手工验证的例子波束偏折是相位梯度超表面最经典的演示。设计目标很简单让垂直入射的平面波经过超表面后以给定的角度θt出射。根据广义斯涅尔定律需要的相位梯度是dΦ/dx k0 · sinθt其中k0是自由空间波数。工程上我们不会让相位连续变化而是按单元周期离散化。假设工作频率f10GHz自由空间波长λ030mm单元周期p7.5mmλ0/4希望偏折角θt30°那么相邻单元之间的相位差应该是ΔΦ k0 · p · sinθt (2π / 30) × 7.5 × 0.5 π/4 45°也就是说8个单元组成一个完整的相位周期0°、45°、90°、135°、180°、225°、270°、315°之后按周期重复。这个计算非常简单我在项目里通常先用一行MATLAB代码算完再把这些相位值映射到单元库里对应的结构上。值得提醒的是这个公式是理想情况下的结果。单元间耦合、介质损耗、相位量化误差都会让实际偏折角偏离理论值但偏差一般不大实测偏折角通常在理论值±2°以内。如果做完方向图测试发现偏折角偏差超过5°那大概率是相位分布和实际单元相位响应不一致了后面第4章会详细分析。3.2 聚焦透镜双曲相位分布的推导与落位波束偏折只需要一维的线性相位梯度而聚焦透镜需要二维的双曲相位分布。设计目标是把入射平面波聚焦到距离超表面为F的焦点上。为达到这个效果超表面上每个位置(x, y)需要提供的相位补偿是φ(x, y) k0 · (√(x² y² F²) − F)这个公式的来源是光程差从超表面上某个位置到焦点的路径长度比从透镜中心到焦点的路径长度要多一段为了等光程就需要相位补偿把这段多出来的路程补回来。我在做聚焦样件时F取的是100mm约3.3个波长。把8×8的单元阵列铺开后每个单元的相位都由上面公式算出。相位分布导出来画成图可以看到中央相位滞后、边缘相位领先的同心环状分布就像一枚摊平的球面透镜。这个相对相位差在阵列边缘能达到很大值。例如阵列半张角对应的相位补偿在F100mm、阵列边长60mm的条件下边缘最大需要的相位补偿约2.5个周期。这意味着单元库里的相位覆盖必须足够富裕不能只刚好卡在0-2π最好留一点冗余给这些边缘单元。3.3 涡旋波束给波束“拧”上轨道角动量涡旋波束是这三类里最有意思的也是相位梯度超表面展示“相位工程”威力的最佳舞台。涡旋波束的波前呈螺旋状携带轨道角动量。它的相位分布是φ(x, y) l · arctan(y / x)其中l称为拓扑荷对应波束的轨道角动量阶数。l1时相位绕阵列中心走一圈恰好变化2π中心处相位不确定所以中心场强为零远场方向图上表现为一个空心环。设计涡旋超表面时我的做法是先用MATLAB把每个单元的(x, y)坐标对应的相位值算出来再量化到单元库。这里最容易犯的错是忘记考虑拓扑荷的旋转方向。l的符号决定涡旋的旋转方向设计前需想清楚是l1还是l-1实测时可以用相位探针或者干涉法验证。如果方向反了只需要把相位分布做一次镜像翻转重新映射单元库即可。另一个注意点涡旋波束的质量对相位量化误差很敏感。如果只用4个相位等级做了量化涡旋纯度会明显下降表现为方向图空心环中间出现亮斑。想要获得干净的涡旋波束我建议至少用8个相位等级也就是相位步进45°。3.4 阵列建模仿真时的边界条件设置单元仿真和阵列仿真之间的衔接是整个全波仿真流程中最容易出问题的地方。单元仿真用的是主从边界加Floquet端口这个设置模拟的是无限周期阵列相当于每单位都处于完全相同的电磁环境中。可实际阵列只有8×8或者16×16边缘单元会“看到”不同的环境因此边缘单元的相位响应会与单元仿真结果有偏差。大阵列比如20×20以上这个偏差可以忽略但小阵列尤其要注意。全波建阵列模型时如果直接把单元库里的结构原样搬到阵列里你会发现Simulation结果和Unit Cell仿真结果有出入。最直接的解决办法是先跑一次完整阵列仿真提取中心单元的实际相位响应与单个单元仿真结果对比确认偏差大小。如果偏差超过10°就要在映射相位分布时加入修正量。另外CST里用Unit Cell边界做单元仿真时端口偏振方向必须和实际入射波的极化方向一致。TE和TM两种极化的响应在高频下会明显不同如果你做的是线极化设计却用圆极化端口扫描参数相位曲线会完全不对这种低级错误很常见排查起来却费时间。4. 仿真、加工、实测三连跳必须踩的坑都在这里4.1 仿真阶段最容易被忽视的单元耦合效应单元之间除了边界失配问题还存在直接的近场耦合。超表面单元间距只有λ/4到λ/2相邻单元的贴片边缘可能只有1-2mm的距离。在这种间距下单元之间的互耦不容忽视。尤其是相位相差较大的相邻单元感受的耦合环境差异非常大。我在一个偏折器项目里就遇到过这种情况单元库扫描时一切正常每个点的相位误差都在3°以内但把16个单元按45°逐级递增排成一行后再用场监视器观察出射波前在相邻单元边界处出现了不规则的跳跃最后用近场扫描仪测出来实际相位分布和设计的偏差最大处达到了17°。解决这个问题有两个思路。第一种是直接在完整阵列模型里提取每个单元的实际相位然后反过来微调单元几何修正这种做法费时但最可靠。第二种是用“耦合感知设计”在单元仿真时就把相邻单元的相位差考虑进去用两个3×3或5×5子阵替代单单元进行扫描虽然实现难度大但结果更贴近真实工作状态。4.2 介质板选型和加工精度对频响的影响仿真模型永远是完美的真实的介质基板却充满不确定性。我在项目里用最多的是F4B板材名义介电常数是2.65但批次之间的实际值可能在2.60到2.70之间浮动。不要小看这个浮动相位梯度的工作频点会跟着介质常数漂移。定量估算一下假如谐振片的设计基于εr2.65实际板材εr2.70那么谐振频率会偏移约1%在10GHz频段就是100MHz。而相位曲线的斜率在谐振点附近非常陡100MHz的偏移可能导致相位误差达到20-30°。这个误差直接会造成波束指向偏移、聚焦效率下降。所以稳妥的做法是多买两批次板材先做介电常数实测再建模型或者故意把工作频点选在板材标称值的中心留出足够裕量。加工精度同样要考虑。普通PCB线宽/槽宽的精度约为±0.05mm到±0.1mm最小线宽建议不要小于0.15mm否则良率会掉得很厉害。如果加工的是10GHz的单元结构0.1mm的误差对谐振频率的影响一般能控制在1%-2%以内是可以接受的但如果做的是Ka波段30GHz的毫米波超表面0.1mm的误差占比就太大了建议考虑带状线工艺或者更高精度的光刻工艺。4.3 实测中绕射、空气隙和参考面的处理实测超表面比仿真麻烦得多。先说样品测试环境透射型样品尺寸有限电磁波照射到样品边缘会产生绕射绕射波叠加到透射波上会污染测量结果。实测时我习惯在样品周围加一圈微波吸收材料比如吸波海绵只留出测试口径这样边缘绕射的影响会小很多。然后说样品夹具。透射型超表面样品垂直放置时一般需要一个泡沫支架。泡沫的介电常数接近1影响很小但如果是用亚克力板之类的材料当支架影响就大了。亚克力的介电常数在2.7左右放到透射波通路里会显著改变相位读数。所以我建议支架一律用高密度泡沫而且样品边缘至少留出2cm的支撑区让金属结构远离支架。相位测量时要特别小心参考面。我早期用矢量网络分析仪做自由空间法测试两个喇叭天线之间隔着样品测S21问题在于如果不做参考面校准测到的相位会包含空气路径的一部分相位这个相位会随喇叭间距、样品厚度变化直接导致测到的相位分布整体偏斜。解决办法是用“穿过测试”校准先把样品拿走测一次“空测”作为参考再放上样品测S21两次相位的差值就是样品本身的相位响应。这种方法测不出绝对相位但测相位分布和相对偏差是足够的。要注意样品两次放置的角度一致夹具要设计好定位槽。5. 从固定到可调动态超表面的方向与我的实战心得5.1 可调超表面的基本思路固定设计的相位梯度超表面做出来性能就定型了波束指向、焦距、拓扑荷都不能变。如果需要同一个器件实现多种功能就要加入可调元件。最常见的做法是给单元结构加载变容二极管或者PIN管。PIN管只有两个状态对应相位0或者180°所以这种方案通常做成1bit可编程超表面通过FPGA给每个单元加高低电平实现数字化的波束扫描。变容管可以在一定偏压范围内连续改变结电容相位响应能连续调谐但它的非线性、损耗一致性都是麻烦事。我在一个动态波束偏折项目中用的方案是单元上集成一对反向并联的变容管通过改变反向偏压改变谐振频率。控制电压从0V到18V变化时相位能实现约300°的调谐范围。但损耗还是比被动结构高很多插入损耗在2.5dB到4dB之间浮动。所以需要根据应用目的在相位调节范围和效率之间取舍。实际项目一般是“静态设计定性能上限动态调节补灵活性”。5.2 经验总结项目中最值得坚持的几个习惯这几年的超表面项目中我逐渐形成了一套固定的工作流程这里把最有价值的几条分享出来第一先手写相位分布再碰全波仿真。我用MATLAB或者Python先搭一个相位计算脚本把目标功能偏折、聚焦、涡旋转换成每个单元的相位目标值用图形画出来检查合理性再进全波仿真软件操作。这一步看起来慢但能避免大量盲目的仿真迭代。第二单元库文件做好版本管理。单元库里的每个结构都对应一个相位值和幅度值我习惯用一张表格记录结构类型、几何参数、S21幅度、S21相位、是否入选工作点。表格里还要留备注列记录这个点是偏谐振点还是非谐振点。这样后面阵列设计遇到问题时回溯定位很快。第三阵列排布时尽量保留相位分布的“平滑性”。也就是说不要为了凑相位而让相邻单元的物理结构跳变太快。如果两个相邻单元需要45°相位差那就选择两个几何参数接近的单元来实现而不是用一个方环和一个十字缝隙紧挨着排。结构剧烈跳变会在局部增加耦合导致实际相位响应偏离期望值。第四加工前至少做一次“加工可行性检查”。最小线宽过小小于0.15mm、槽宽太窄、铜皮占比太高等问题都要在发板前和板厂工程师确认清楚。省下的时间往往比仿真省的时间多得多。最后再提一个小经验测试整板前先用网络分析仪单独测几个“陪测样”也就是和超表面同批加工但没有排相位的小块单元阵列。通过陪测样的透射幅度和相位可以在正式测试前就把介质板实际参数校准进模型里。这个步骤能省下至少两轮“仿真-加工-测试”的重复迭代尤其是当你第一次使用的板材或工艺时。本文还有配套的精品资源点击获取
相位梯度超表面设计实战：从广义斯涅尔定律到波束操控
简介一份关于相位梯度超表面结构操控电磁波束的学术研究文献面向电磁超材料、微波通信与雷达技术方向的研究人员、工程师及高校师生。内容从高阻抗表面单元尺寸与反射相位的关系入手设计了一种相位梯度超表面结构并通过数值仿真验证其在10 GHz频段可实现约22°反射波束偏转体现了厚度薄、重量轻且可灵活操控波前相位的特点。文献还系统梳理了相位梯度超表面的代表性进展如广义斯涅尔定律、涡旋波束生成、异常反射和表面波耦合等为读者理解超表面波束调控原理和分析方法提供了完整脉络。该论文刊载于《浙江科技学院学报》2016年第28卷第3期作者单位涵盖同济大学与浙江科技学院结构包含摘要、关键词、理论分析和数值仿真等内容。资源为1个PDF文件大小约2.04MB已有219人浏览学习适合作为相关课题的参考文献或专业学习材料。 很多人刚接触超表面时第一反应都是先跑仿真在 HFSS 或者 CST 里画个单元扫几个参数然后急着排阵列。但折腾几轮之后你会慢慢意识到相位梯度超表面对电磁波束的操控本质上是一个“相位工程”问题。你前期把相位算得越明白后续的仿真、加工、实测就越顺。这篇文章我结合自己做相位梯度超表面的实际项目经历把从原理到设计、再到仿真和实测的完整链路拆开讲。内容涉及广义斯涅尔定律的工程理解、单元库搭建、波束偏折/聚焦/涡旋的设计方法以及加工和测试里那些文献不会写、但非常影响结果的坑。适合刚入门超表面的研究生也适合想把固定结构往可调方向延伸的工程师参考。1. 从相位不连续说起超表面为什么能实现波束操控1.1 广义斯涅尔定律超表面的底层逻辑传统的透镜、棱镜操控电磁波靠的是电磁波在介质中传播时积累的连续相位。波穿过一块透镜中心路径和边缘路径走过的光学长度不同出射面上就形成一个相位分布波前因此被弯曲。这套逻辑很成熟缺点是器件体积大、重量重、剖面高。超表面的思路完全不一样。它不需要波在体内走很长的距离而是在一个极薄的界面上让每个单元对电磁波施加一个突变相位。这个相位突变沿界面的变化率就是常说的“相位梯度”。2011年Capasso组的广义斯涅尔定律给出了定量描述Bsinθt - nisinθi (λ0 / 2π) · dΦ/dx这个公式看着长其实物理图像很清晰入射波和透射波在界面切线方向的波矢分量之差由界面的相位梯度来补偿。右边那一项 dΦ/dx 就是相位沿界面的变化率它的存在意味着即使垂直入射θi0透射波也可以被“推”向一个非零的出射角。换句话说相位梯度相当于给透射波一个额外的横向动量。我在项目里习惯把这条定律反过来用先根据目标偏折角算需要的相位梯度再反推每个单元应该落在这个梯度曲线的哪个位置。这个过程很像给波前画刻度尺每格分多少度相位排多少个单元计算清楚了再动手画结构。1.2 与相控阵的类比一块“超薄的移相阵面”理解超表面最直观的方式其实是把它想成一块无源的相控阵天线。相控阵靠每个通道的移相器给单元馈电通过控制通道间的相位差把主波束指向某个方向。超表面做的事情在数学上完全一致只是它的“移相器”变成了每一个亚波长单元结构而控制手段不是给电信号而是通过改变单元的几何尺寸、形状或者朝向来实现不同的相位响应。这个类比在工程上很有用。阵列天线里的副瓣电平、栅瓣抑制、单元间距问题都可以平移过来指导超表面设计。比如单元周期必须小于工作波长目的就是避免出现栅瓣相位量化级数太少会产生寄生波束这和相控阵里数字移相器的量化误差是同一个物理问题。我做第一个超表面样件时因为当时手里有一块标准的8×8微带阵列所以我对“相位差决定波束指向”这件事特别熟悉。后来把这种思维带到超表面上后很多概念一下子对上了。可以这么说超表面就是一张把移相网络直接“长”在辐射面上的阵面只是每个单元的“馈电”是自由空间里的入射波。2. 相位怎么“铺”到表面上单元结构与相位覆盖设计2.1 单元怎么选从H形到方环的取舍超表面单元的选择决定了你能获得的相位范围、幅度稳定性以及加工容差。目前最常用的单元可以归成几类贴片型矩形贴片、圆贴片、缝隙型H形缝隙、十字缝隙、环型方环、圆环、双开口环以及介质型介质柱等。项目里最常用的还是金属贴片/缝隙类因为它们可以直接用PCB工艺做出来。这些单元各有各的脾气。矩形贴片加工最简单、设计变量少但单个贴片的相位覆盖范围有限想凑出完整0-2π往往要依赖多个变量的组合H形缝隙可以通过改变两个分支的长度来实现较大的相位变化但对极化方向敏感排布时要注意入射波的极化方向一致性方环单元的相位曲线比较平滑而且通过外环尺寸和槽宽两个自由度能较好地兼顾相位范围和幅度稳定。我自己在10GHz频段做透射型超表面时最常用的是方环和十字缝隙的组合。介质基板选的是F4B相对介电常数2.65损耗角正切0.001左右厚度1mm。这个组合的好处是相位变化范围大而且F4B这种材料在C波段到X波段的介电常数比较稳定价格也适中不像Rogers板材那么贵。下面这个表格是我在实际项目中对比过的几种单元方案方便你根据需求快速做取舍单元类型相位覆盖能力透射/反射效率加工难度带宽表现适用场景矩形贴片中等透射一般低窄快速验证H形缝隙较大透射较高中中等线极化透射阵方环大较高中中等宽带偏折/聚焦方环互补结构很大高中高较宽需要大相位梯度的场景介质柱大高高宽高频/毫米波2.2 参数扫描的核心既要覆盖2π也要幅度均匀选定单元类型后下一步就是对单元参数做扫描。扫描的目标有两个第一找出能覆盖0-2π的相位变化范围第二所有工作点的透射/反射幅度尽量高且均匀。很多新手做单元仿真时就盯着相位曲线看到相位扫出360°就觉得万事大吉。但实际排完阵列后一做仿真方向图增益很低一查原因幅度响应在单元之间差了5dB甚至更多。幅度不均匀意味着阵列的口径效率直接受损因为远场叠加时各个单元“出力”不一致。我在做单元库时一般按以下步骤进行在CST或HFSS里建立单个周期单元模型边界条件使用主从边界unit cell端口设置为Floquet端口。固定其他参数不变扫描一个关键几何参数。比如方环结构固定环宽w不变扫描外环边长L从4.0mm到7.2mm步长0.2mm。提取每个样本的S21幅度和相位保存成工程文件。筛选出幅度大于阈值我一般取0.85以上的样本点再检查这些点对应的相位是否覆盖了完整的2π。如果相位覆盖不够增加第二个扫描变量。比如同时调槽宽或者加入介质层的厚度变化用二维扫描来拓宽相位范围。这里有个关键细节单一参数扫描的相位曲线往往是S形或者单调但范围不足。比如方环的谐振会影响相位在谐振点附近相位变化非常剧烈但幅度掉得也厉害。太过靠近谐振点的工作状态虽然能拿到必要的相位值但损耗会成为问题。所以最终单元库里每个入选的结构都要尽量避开谐振点太近的位置。2.3 透射还是反射两种工作模式的分水岭决定超表面是透射型还是反射型是整个项目最开始就要确定的事因为它影响后续所有的单元设计、馈电方式和测试方案。透射型的优势是馈电方便喇叭天线从法向照射透射波在另一侧直接形成需要的波前测试时收发天线分居样品两侧比较好布置而且可以做成类似传统透镜的器件。劣势也很明显透射型单元要在“匹配”和“相位变化”之间做折衷很难在获得大相位调节范围的同时保持很高的透射率。尤其是介质板两侧都有金属结构时谐振带来的插损往往有1-2dB。反射型的优势是效率高通常能拿到90%以上的反射率而且相位覆盖范围大劣势是需要考虑馈源遮挡双站测试时收发喇叭要分处同侧而且反射型超表面往往需要加金属背板剖面会略高一些。做偏折器或者聚焦透镜两边都能做但如果效率是首要指标我会选反射型如果项目要求空间馈电方便、需要直接观测透射波前那就只能做透射型。好在绝大多数发论文的场景透射型已经足够说明问题所以我后面的举例也以透射型为主。3. 从单元到阵列偏折、聚焦、涡旋三种典型的波束操控3.1 波束偏折一个可手工验证的例子波束偏折是相位梯度超表面最经典的演示。设计目标很简单让垂直入射的平面波经过超表面后以给定的角度θt出射。根据广义斯涅尔定律需要的相位梯度是dΦ/dx k0 · sinθt其中k0是自由空间波数。工程上我们不会让相位连续变化而是按单元周期离散化。假设工作频率f10GHz自由空间波长λ030mm单元周期p7.5mmλ0/4希望偏折角θt30°那么相邻单元之间的相位差应该是ΔΦ k0 · p · sinθt (2π / 30) × 7.5 × 0.5 π/4 45°也就是说8个单元组成一个完整的相位周期0°、45°、90°、135°、180°、225°、270°、315°之后按周期重复。这个计算非常简单我在项目里通常先用一行MATLAB代码算完再把这些相位值映射到单元库里对应的结构上。值得提醒的是这个公式是理想情况下的结果。单元间耦合、介质损耗、相位量化误差都会让实际偏折角偏离理论值但偏差一般不大实测偏折角通常在理论值±2°以内。如果做完方向图测试发现偏折角偏差超过5°那大概率是相位分布和实际单元相位响应不一致了后面第4章会详细分析。3.2 聚焦透镜双曲相位分布的推导与落位波束偏折只需要一维的线性相位梯度而聚焦透镜需要二维的双曲相位分布。设计目标是把入射平面波聚焦到距离超表面为F的焦点上。为达到这个效果超表面上每个位置(x, y)需要提供的相位补偿是φ(x, y) k0 · (√(x² y² F²) − F)这个公式的来源是光程差从超表面上某个位置到焦点的路径长度比从透镜中心到焦点的路径长度要多一段为了等光程就需要相位补偿把这段多出来的路程补回来。我在做聚焦样件时F取的是100mm约3.3个波长。把8×8的单元阵列铺开后每个单元的相位都由上面公式算出。相位分布导出来画成图可以看到中央相位滞后、边缘相位领先的同心环状分布就像一枚摊平的球面透镜。这个相对相位差在阵列边缘能达到很大值。例如阵列半张角对应的相位补偿在F100mm、阵列边长60mm的条件下边缘最大需要的相位补偿约2.5个周期。这意味着单元库里的相位覆盖必须足够富裕不能只刚好卡在0-2π最好留一点冗余给这些边缘单元。3.3 涡旋波束给波束“拧”上轨道角动量涡旋波束是这三类里最有意思的也是相位梯度超表面展示“相位工程”威力的最佳舞台。涡旋波束的波前呈螺旋状携带轨道角动量。它的相位分布是φ(x, y) l · arctan(y / x)其中l称为拓扑荷对应波束的轨道角动量阶数。l1时相位绕阵列中心走一圈恰好变化2π中心处相位不确定所以中心场强为零远场方向图上表现为一个空心环。设计涡旋超表面时我的做法是先用MATLAB把每个单元的(x, y)坐标对应的相位值算出来再量化到单元库。这里最容易犯的错是忘记考虑拓扑荷的旋转方向。l的符号决定涡旋的旋转方向设计前需想清楚是l1还是l-1实测时可以用相位探针或者干涉法验证。如果方向反了只需要把相位分布做一次镜像翻转重新映射单元库即可。另一个注意点涡旋波束的质量对相位量化误差很敏感。如果只用4个相位等级做了量化涡旋纯度会明显下降表现为方向图空心环中间出现亮斑。想要获得干净的涡旋波束我建议至少用8个相位等级也就是相位步进45°。3.4 阵列建模仿真时的边界条件设置单元仿真和阵列仿真之间的衔接是整个全波仿真流程中最容易出问题的地方。单元仿真用的是主从边界加Floquet端口这个设置模拟的是无限周期阵列相当于每单位都处于完全相同的电磁环境中。可实际阵列只有8×8或者16×16边缘单元会“看到”不同的环境因此边缘单元的相位响应会与单元仿真结果有偏差。大阵列比如20×20以上这个偏差可以忽略但小阵列尤其要注意。全波建阵列模型时如果直接把单元库里的结构原样搬到阵列里你会发现Simulation结果和Unit Cell仿真结果有出入。最直接的解决办法是先跑一次完整阵列仿真提取中心单元的实际相位响应与单个单元仿真结果对比确认偏差大小。如果偏差超过10°就要在映射相位分布时加入修正量。另外CST里用Unit Cell边界做单元仿真时端口偏振方向必须和实际入射波的极化方向一致。TE和TM两种极化的响应在高频下会明显不同如果你做的是线极化设计却用圆极化端口扫描参数相位曲线会完全不对这种低级错误很常见排查起来却费时间。4. 仿真、加工、实测三连跳必须踩的坑都在这里4.1 仿真阶段最容易被忽视的单元耦合效应单元之间除了边界失配问题还存在直接的近场耦合。超表面单元间距只有λ/4到λ/2相邻单元的贴片边缘可能只有1-2mm的距离。在这种间距下单元之间的互耦不容忽视。尤其是相位相差较大的相邻单元感受的耦合环境差异非常大。我在一个偏折器项目里就遇到过这种情况单元库扫描时一切正常每个点的相位误差都在3°以内但把16个单元按45°逐级递增排成一行后再用场监视器观察出射波前在相邻单元边界处出现了不规则的跳跃最后用近场扫描仪测出来实际相位分布和设计的偏差最大处达到了17°。解决这个问题有两个思路。第一种是直接在完整阵列模型里提取每个单元的实际相位然后反过来微调单元几何修正这种做法费时但最可靠。第二种是用“耦合感知设计”在单元仿真时就把相邻单元的相位差考虑进去用两个3×3或5×5子阵替代单单元进行扫描虽然实现难度大但结果更贴近真实工作状态。4.2 介质板选型和加工精度对频响的影响仿真模型永远是完美的真实的介质基板却充满不确定性。我在项目里用最多的是F4B板材名义介电常数是2.65但批次之间的实际值可能在2.60到2.70之间浮动。不要小看这个浮动相位梯度的工作频点会跟着介质常数漂移。定量估算一下假如谐振片的设计基于εr2.65实际板材εr2.70那么谐振频率会偏移约1%在10GHz频段就是100MHz。而相位曲线的斜率在谐振点附近非常陡100MHz的偏移可能导致相位误差达到20-30°。这个误差直接会造成波束指向偏移、聚焦效率下降。所以稳妥的做法是多买两批次板材先做介电常数实测再建模型或者故意把工作频点选在板材标称值的中心留出足够裕量。加工精度同样要考虑。普通PCB线宽/槽宽的精度约为±0.05mm到±0.1mm最小线宽建议不要小于0.15mm否则良率会掉得很厉害。如果加工的是10GHz的单元结构0.1mm的误差对谐振频率的影响一般能控制在1%-2%以内是可以接受的但如果做的是Ka波段30GHz的毫米波超表面0.1mm的误差占比就太大了建议考虑带状线工艺或者更高精度的光刻工艺。4.3 实测中绕射、空气隙和参考面的处理实测超表面比仿真麻烦得多。先说样品测试环境透射型样品尺寸有限电磁波照射到样品边缘会产生绕射绕射波叠加到透射波上会污染测量结果。实测时我习惯在样品周围加一圈微波吸收材料比如吸波海绵只留出测试口径这样边缘绕射的影响会小很多。然后说样品夹具。透射型超表面样品垂直放置时一般需要一个泡沫支架。泡沫的介电常数接近1影响很小但如果是用亚克力板之类的材料当支架影响就大了。亚克力的介电常数在2.7左右放到透射波通路里会显著改变相位读数。所以我建议支架一律用高密度泡沫而且样品边缘至少留出2cm的支撑区让金属结构远离支架。相位测量时要特别小心参考面。我早期用矢量网络分析仪做自由空间法测试两个喇叭天线之间隔着样品测S21问题在于如果不做参考面校准测到的相位会包含空气路径的一部分相位这个相位会随喇叭间距、样品厚度变化直接导致测到的相位分布整体偏斜。解决办法是用“穿过测试”校准先把样品拿走测一次“空测”作为参考再放上样品测S21两次相位的差值就是样品本身的相位响应。这种方法测不出绝对相位但测相位分布和相对偏差是足够的。要注意样品两次放置的角度一致夹具要设计好定位槽。5. 从固定到可调动态超表面的方向与我的实战心得5.1 可调超表面的基本思路固定设计的相位梯度超表面做出来性能就定型了波束指向、焦距、拓扑荷都不能变。如果需要同一个器件实现多种功能就要加入可调元件。最常见的做法是给单元结构加载变容二极管或者PIN管。PIN管只有两个状态对应相位0或者180°所以这种方案通常做成1bit可编程超表面通过FPGA给每个单元加高低电平实现数字化的波束扫描。变容管可以在一定偏压范围内连续改变结电容相位响应能连续调谐但它的非线性、损耗一致性都是麻烦事。我在一个动态波束偏折项目中用的方案是单元上集成一对反向并联的变容管通过改变反向偏压改变谐振频率。控制电压从0V到18V变化时相位能实现约300°的调谐范围。但损耗还是比被动结构高很多插入损耗在2.5dB到4dB之间浮动。所以需要根据应用目的在相位调节范围和效率之间取舍。实际项目一般是“静态设计定性能上限动态调节补灵活性”。5.2 经验总结项目中最值得坚持的几个习惯这几年的超表面项目中我逐渐形成了一套固定的工作流程这里把最有价值的几条分享出来第一先手写相位分布再碰全波仿真。我用MATLAB或者Python先搭一个相位计算脚本把目标功能偏折、聚焦、涡旋转换成每个单元的相位目标值用图形画出来检查合理性再进全波仿真软件操作。这一步看起来慢但能避免大量盲目的仿真迭代。第二单元库文件做好版本管理。单元库里的每个结构都对应一个相位值和幅度值我习惯用一张表格记录结构类型、几何参数、S21幅度、S21相位、是否入选工作点。表格里还要留备注列记录这个点是偏谐振点还是非谐振点。这样后面阵列设计遇到问题时回溯定位很快。第三阵列排布时尽量保留相位分布的“平滑性”。也就是说不要为了凑相位而让相邻单元的物理结构跳变太快。如果两个相邻单元需要45°相位差那就选择两个几何参数接近的单元来实现而不是用一个方环和一个十字缝隙紧挨着排。结构剧烈跳变会在局部增加耦合导致实际相位响应偏离期望值。第四加工前至少做一次“加工可行性检查”。最小线宽过小小于0.15mm、槽宽太窄、铜皮占比太高等问题都要在发板前和板厂工程师确认清楚。省下的时间往往比仿真省的时间多得多。最后再提一个小经验测试整板前先用网络分析仪单独测几个“陪测样”也就是和超表面同批加工但没有排相位的小块单元阵列。通过陪测样的透射幅度和相位可以在正式测试前就把介质板实际参数校准进模型里。这个步骤能省下至少两轮“仿真-加工-测试”的重复迭代尤其是当你第一次使用的板材或工艺时。本文还有配套的精品资源点击获取 
