搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

动态可调谐的频域多功能可重构极化转换超表面

黄晓俊 高焕焕 何嘉豪 栾苏珍 杨河林

黄晓俊, 高焕焕, 何嘉豪, 栾苏珍, 杨河林. 动态可调谐的频域多功能可重构极化转换超表面. 物理学报, 2022, 71(22): 224102. doi: 10.7498/aps.71.20221256
引用本文: 黄晓俊, 高焕焕, 何嘉豪, 栾苏珍, 杨河林. 动态可调谐的频域多功能可重构极化转换超表面. 物理学报, 2022, 71(22): 224102. doi: 10.7498/aps.71.20221256
Huang Xiao-Jun, Gao Huan-Huan, He Jia-Hao, Luan Su-Zhen, Yang He-Lin. Dynamically tunable frequency-domain multifunctional reconfigurable polarization conversion metasurface. Acta Phys. Sin., 2022, 71(22): 224102. doi: 10.7498/aps.71.20221256
Citation: Huang Xiao-Jun, Gao Huan-Huan, He Jia-Hao, Luan Su-Zhen, Yang He-Lin. Dynamically tunable frequency-domain multifunctional reconfigurable polarization conversion metasurface. Acta Phys. Sin., 2022, 71(22): 224102. doi: 10.7498/aps.71.20221256

动态可调谐的频域多功能可重构极化转换超表面

黄晓俊, 高焕焕, 何嘉豪, 栾苏珍, 杨河林

Dynamically tunable frequency-domain multifunctional reconfigurable polarization conversion metasurface

Huang Xiao-Jun, Gao Huan-Huan, He Jia-Hao, Luan Su-Zhen, Yang He-Lin
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 设计任意调控极化的电磁器件是一个研究热点, 其中多功能可重构电磁器件在雷达、卫星通信等领域有着广泛的应用. 本文设计了一种基于正本征负极(PIN)二极管可调谐的多功能可重构极化转换超表面, 可以实现不同频段内的线极化波转换、线极化波-圆极化波转换和全反射功能的切换, 在斜入射角小于30°时, 多功能转换器能保持高效的宽带极化转换特性. 这种转换和重构特性主要是由于结构本身的各向异性和PIN管不同状态时耦合模式的改变. 此外, 利用表面电流解释了偏振转换的物理机理, 电谐振和磁谐振的共同作用导致了偏振转换. 最后, 对该结构样品进行实验验证, 其结果与仿真吻合较好. 该器件在极化调控、频率控制、智能反射面设计和天线设计等方面具有潜在的应用价值.
    The design of electromagnetic device with arbitrary polarization manipulation is the hot spot of the current research. Multifunctional reconfigurable electromagnetic devices have been put into wide applications in radar, satellite communication and other fields. In this work designed is a multifunctional reconfigurable polarization conversion metasurface based on two PIN diodes, which can realize linear polarization conversion, linear-circular polarization conversion and total reflection switching in the different frequency bands, and the multi-function converter can still maintain the high-efficient broadband conversion characteristics when the oblique incidence angle is less than 30 degrees. The polarization conversion and reconfigurability are mainly due to the anisotropy of the structure and the changes of coupling mode in different states of PIN diodes. In addition, the physical mechanism of polarization conversion is explained by surface current. The combined action of electric resonance and magnetic resonance results in polarization conversion. Finally, the experimental results of the sample are in good agreement with the simulation results. The device has the potential application values in polarization manipulation, frequency control, intelligent reflecting surface design, and antenna design and so on.
      PACS:
      41.20.Jb(Electromagnetic wave propagation; radiowave propagation)
      78.67.Pt(Multilayers; superlattices; photonic structures; metamaterials)
      42.25.Ja(Polarization)
      通信作者: 黄晓俊, hxj@xust.edu.cn ; 杨河林, emyang@mail.ccnu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 61974119)和陕西省自然科学基金(批准号: 2021JM-395)资助的课题.
      Corresponding author: Huang Xiao-Jun, hxj@xust.edu.cn ; Yang He-Lin, emyang@mail.ccnu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 61974119), Natural Science Foundation of Shaanxi Province (Grant No. 2021JM-395).

    极化是雷达通信、雷达目标识别、光通信和卫星通信系统中一个重要的物理属性. 目前, 设计任意调控极化功能的电磁器件成为电磁领域的一个热点课题, 传统方法是利用自然材料的双折射效应, 通过增大厚度来实现相位的累积, 但会导致器件的体积过大, 难以小型化; 另外, 自然材料的损耗较大, 会导致极化调控的效率较低. 电磁超表面作为超材料的二维平面结构, 利用金属结构单元的电磁谐振实现对电磁波的振幅、相位、频率和极化的调控, 为灵活高效地操纵电磁波提供了一种有效的途径, 并在诸多领域得到了广泛的应用. 目前, 在极化转换[1-7]、电磁吸收[8-10]、波束调控[11-16]、天线设计[17,18]和传感器[19]等领域已经取得了相当多的成果. Wang等[20]设计了一种可拆卸的编码超表面, 通过手动重新配置编码元件的阵列布局, 实现了对振幅、相位和偏振状态的独立控制, 但存在重构速度慢, 集成度差的问题. Tian等[21]提出了一种可重构的超宽带反射型极化转换超表面, 实现了线极化转换和线-圆极化转换功能. He等[22]设计了一种基于PIN二极管的线极化可重构超表面, 实现了线极化全反射与线-圆极化转换的双功能. 此外, 更高集成度的多功能电磁器件也相继被设计与验证, 如Li等[3]提出了一种宽带多功能可切换极化转换超表面, 实现了包括极化分离与极化转换在内的多种极化功能的可重构功能. Liu等[23]提出了一种在频域可操纵电磁波的可重构超表面, 实现了不同频率下的线-圆极化转换功能. Wang等[24]提出了一种电磁波全空间控制的可重构多功能超表面, 实现了电磁波的反射、透射和吸收功能, 特别是在全反射状态下验证了可重构多功能超表面在隐身天线罩中的潜在应用. 上述成果都是利用PIN二极管的通断实现了多种电磁波极化可重构这一功能. 随着大规模通信设备的发展和部署, 电子通信设备之间的交互变得越来越复杂, 这就要求电子设备系统中必须集成多种电磁调控器件, 然而集成多种功能的器件必然会出现设备难以小型化这一问题. 因此, 设计具有多功能集成于一体的电磁调控器件对于实现未来通信设备小型化具有重要意义. 另外, 同时实现极化和频率重构的多功能电磁器件的研究仍然很少, 在这方面还存在诸多挑战.

    本文设计了一种基于PIN二极管的可调谐多功能可重构极化转换超表面, 实现了不同频带内的极化可重构. 在3.70—4.16 GHz内实现了相对带宽为12.1%、转换率大于0.9的线极化转换; 在4.25—4.95 GHz实现了线极化波到圆极化波的转换, 轴比小于3 dB的相对带宽达到15.1%; 在6.02—7.16 GHz内实现了全反射功能, 相对带宽达到17.5%, 极化转换率和吸收率均小于0.1; 在7.30—8.15 GHz实现了相对带宽为11.0%、转换率大于0.9的线极化转换. 这种具有频域极化选择的超表面不仅具有频率可重构特性, 还实现了不同极化波的转换. 该器件在极化调控、频率控制、智能反射面设计和天线设计等方面具有一定的应用价值.

    在频域中, 电磁波可以调制为不同频率分量的散射波. 根据频域调控电磁波的相关理论, 超表面的反射电场和入射电场之间可用散射参数S来表征, 散射参数S主要与频率相关. 具体的表达式如下:

    [ErxEry]=S[EixEiy],
    (1)

    其中, i和r分别表示电磁场的入射与反射, xy分别表示电磁波的极化状态. 对于第n个结构单元的散射矩阵可表示为

    Sn=[Sfn(f)00Sfn(f)],
    (2)

    其中, Sfn(f)表示第n个单元的散射参数. 由于散射矩阵S满足Jordan矩阵, 所以在各个频率范围内的电磁响应将互不干扰. 因此反射波的总电场可表示为

    E=n=1Ern=(n=1Sn)Ein.
    (3)

    Sfn(f)可根据泰勒级数进行展开, 其定义如下:

    Sfn(f)=P=0(a1mpP!(ff1)p+a2mpP!(ff2)p+a3mpP!(ff3)p+a4mpP!(ff4)p+L),
    (4)

    式中, f1, f2, f3, f4Sfn(f)不同的谐振频率, a1mp, a2mp, a3mp, a4mp是不同谐振频率下的谐振系数. 经以上分析, 固定的超表面结构很难实现不同频带内极化功能的切换, 因此, 本文通过在无源超表面结构中嵌入PIN二极管来实现极化和频率同时可重构.

    本文设计了一种多缝隙的超表面结构, 图1(a)是该结构的立体图, 图1(b)为主视图, 图1(c)x-z平面的侧视图. 基本单元由两层介质和三层金属层组成, 其中金属层包括有源切换贴片层、金属地层和偏置馈线层. 金属层的材料均为铜, 厚度为0.035 mm, 电导率为5.8 × 107 S/m, 介质基板为F4B材料, 介电常数为2.65, 损耗正切值为0.001. 具体几何参数如下: P = 15.0 mm, l0 = 14.2 mm, l1 = 4 mm, l2 = 14 mm, w = 0.5 mm, t1 = 3.0 mm, t2 = 0.5 mm, g = 1.5 mm, via = 0.8 mm, n1 = 7.2 mm, n2 = 4.9 mm, n3 = 1.9 mm, n4 = 10.2 mm, m1 = 10.1 mm, m2 = 13.1 mm, m3 = 17.3 mm, m4 = 14.3 mm. 有源切换层是由正方形贴片开槽构成的, 其通过中间的导通孔与金属地层相连接, 两个提供偏置电压的导通孔距离中心导通孔的距离为5∶4, 两个二极管位于两个开槽的中心位置. 值得注意的是, 有源切换层上的3个贴片的大小保证了超表面实现频域可重构的功能, 且不同频段的谐振频率为f1, f2, f3, f4; 多条缝隙型结构会导致更多的谐振点的产生, 而这些谐振点的叠加有效拓宽了不同极化功能下的频带带宽. 此外, 合理设计偏置馈线有利于降低整个结构的剖面, 也降低了更大的损耗.

    图 1 单元结构示意图 (a) 立体图; (b) 主视图; (c) 侧视图\r\nFig. 1. Schematics of unit cell: (a) Solid Shape; (b) main view; (c) side view.
    图 1  单元结构示意图 (a) 立体图; (b) 主视图; (c) 侧视图
    Fig. 1.  Schematics of unit cell: (a) Solid Shape; (b) main view; (c) side view.

    使用CST Microwave studio 对单元结构进行全波仿真计算, x, y方向分别设置为unit cell, –z方向设置为电边界, +z方向设置为open and space边界. 结构中加载的二极管型号为SMP2019-079 LF, 在仿真中, 二极管导通时, 定义为“1”状态, 将其等效为电感Lp = 0.7 nH与电阻RS = 0.5 Ω的串联电路; 二极管断开时, 定义为“0”状态, 可等效为电感Lp = 0.5 nH与电阻Ct = 0.24 pF的串联电路[25]. 基本单元结构上两个二极管的状态可用2 bit编码来表示(“10”, “01”, “00”, “11”), 例如“10”表示PIN1导通PIN2断开. 图2所示为二极管位于不同状态下单元结构的等效电路, 以TE波入射为例, 3个缝隙贴片等效为电阻-电感-电容(RLC)串联电路, 图2(a)可得, 周期多缝隙型结构平行于电场方向表现为电感, 分别等效为L1, L2, L3; C1, C2, C3分别为周期多缝隙贴片之间的等效电容; 两个相邻缝隙型金属贴片之间的耦合呈现为电容特性, 则可等效为Cg1, Cg2; 长度为h的传输线表示介质基板; Z0Z1分别表示自由空间和介质基板的等效阻抗; 图2(b)蓝色区域展示了二极管处于导通与截止状态下的等效电路, 当二极管导通时, 可等效为电感Lp与电阻RS的串联; 二极管截止时, 等效为电感Lp与电阻Ct的串联.

    图 2 等效电路 (a) 等效电路分析; (b) 等效电路模型\r\nFig. 2. Equivalent circuit: (a) Equivalent circuit analysis; (b) equivalent circuit model.
    图 2  等效电路 (a) 等效电路分析; (b) 等效电路模型
    Fig. 2.  Equivalent circuit: (a) Equivalent circuit analysis; (b) equivalent circuit model.

    在频域中, 可重构超表面可等效为一个多功能极化转换器, 通过计算轴比(AR)、极化转换率(PCR)和吸收率(A)来描述超表面在不同频率范围内的极化转换能力, 具体定义如下:

    AR=|ryx|2+|rxx|2+a|ryx|2+|rxx|2a,
    (5)

    其中:

    a=|ryx|4+|rxx|4+2|ryx|2|rxx|2cos(2Δφyx),
    (6)
    Δφyx=φyxφxx,
    (7)
    PCRx=|ryx|2|ryx|2+|rxx|2,
    (8)
    A=1|ryx|2|rxx|2,
    (9)
    A=1|ryx|2|rxx|2.
    (10)

    图3x极化波和y极化波入射时, 两个二极管处于不同状态时反射系数的仿真结果. 当PIN1导通PIN2断开时, 仿真结果如图3(a)所示, 从仿真结果可以清晰地看出, 交叉极化的反射系数ryxrxy在3.70—4.16 GHz接近0 dB, 而共极化反射系数rxxryy小于–10 dB. 由此说明, 在二极管状态为“10”时, 入射的线极化波转化成为交叉极化波, 超表面作为线极化转换器. 当PIN1断开PIN2导通时, 仿真结果如图3(b)所示, 在4.35—4.90 GHz时, 交叉极化ryxrxy和共极化rxxryy的反射系数接近–3 dB, 图3(c)中可以看出相应的相位差保持在–π/2. 由此说明, 在二极管状态为“01”时, 入射的线极化波转换为左旋圆极化波, 超表面是一个线-圆极化转换器. 当PIN1和PIN2同时断开时, 仿真结果如图3(d)所示, 交叉极化的反射系数ryxrxy在6.02—7.16 GHz小于–10 dB, 而共极化反射系数rxxryy接近0 dB. 由此表明, 在二极管状态为“00”时, 入射的线极化波以共极化的形式全反射, 超表面等效为一个完美电导体. 当PIN1和PIN2同时导通时, 仿真结果如图3(e)所示, 共极化反射系数rxxryy在7.30—8.05 GHz小于–10 dB, 而交叉极化反射系数ryxrxy接近0 dB. 由此说明, 在二极管状态为“11”时, 入射的线极化波转换为交叉极化波, 超表面等效为一个线极化转换器. 因此, 可以看出极化转换器的这种多功能性在频域上表现出可重构特性.

    图 3 二极管不同状态下的仿真曲线 (a) “10”-反射系数; (b) “01”-反射系数; (c) “01”-反射相位; (d) “10”-反射系数; (e) “11”-反射系数\r\nFig. 3. Simulation curve of diodes in different states: (a) “10”-reflectance; (b) “01”-reflectance; (c) “01”-reflection phase; (d) “10”- reflectance; (e) “11”- reflectance.
    图 3  二极管不同状态下的仿真曲线 (a) “10”-反射系数; (b) “01”-反射系数; (c) “01”-反射相位; (d) “10”-反射系数; (e) “11”-反射系数
    Fig. 3.  Simulation curve of diodes in different states: (a) “10”-reflectance; (b) “01”-reflectance; (c) “01”-reflection phase; (d) “10”- reflectance; (e) “11”- reflectance.

    为了更好地描述多功能频域可重构超表面的性能, 根据反射系数计算了相应的极化转换率、吸收率和轴比, 如图4所示. 图4(a)显示了PIN二极管在“10”和“01”两种状态下的转换率和轴比, 可以看出当二极管处于“10”状态时, PCR在3.70—4.16 GHz大于0.9, 这意味着将入射的x线极化波调控为与其垂直的y极化波, 且相对带宽达到12.1%. 当二极管处于“01”状态时, AR在4.25—4.95 GHz小于3 dB, 相对带宽达到了15.1%, 这意味着入射的x线极化波调控为左旋圆极化波. 同样, 图4(b)计算了PIN二极管在“00”和“11”两种状态下的转换率和吸收率, 当二极管处于“00”状态时, PCR和A在6.02—7.16 GHz小于0.1, 相对带宽达到了17.5%, 这表明入射的线极化波在当前二极管状态下可以实现完美的全反射; 当二极管处于“11”状态时, PCR在7.30—8.15 GHz大于0.9, 这意味着将入射的x线极化波调控为y极化波, 实现了线极化转换.

    图 4 二极管不同状态下的极化转换率、吸收率和轴比 (a) “10”与“01”; (b) “00”与“11”\r\nFig. 4. Polarization conversion rate, absorptance and axial ratio of diodes in different states: (a) “10” and “01”, (b) “00” and “11”.
    图 4  二极管不同状态下的极化转换率、吸收率和轴比 (a) “10”与“01”; (b) “00”与“11”
    Fig. 4.  Polarization conversion rate, absorptance and axial ratio of diodes in different states: (a) “10” and “01”, (b) “00” and “11”.

    为了验证本设计的合理性和正确性, 图5(a)为一个15 × 15结构单元的样品模型, 图5(b)为该模型的馈线. 采用PCB工艺对样品进行加工, 同样包含了15 × 15个基本单元以及底部馈电层, 如图5(c)(d)所示, 每个单元结构中嵌入两个相同的PIN二极管, 导通孔与馈电层相连接, 从而便于为样品模型上的二极管提供偏置电压, 两个二极管的负极采用共地的方式连接到金属地层. 在实际测试中, 将样品固定于暗室的吸波材料上, 测试环境如图6所示, 将一组工作在2—18 GHz频段内的喇叭天线固定在三角支架上, 通过同轴馈线连接到罗德施瓦茨矢量网络分析仪上, 天线距离可重构超表面样品为450 mm, 喇叭天线与样品应处于同一水平位置上, 天线之间的夹角为5°左右, 为了保证二极管的导通状态, 标准直流电源箱为PIN二极管提供+3 V的偏置电压, 从而实现可重构极化功能的切换.

    图 5 可重构超表面样品结构与实际测试样品 (a) 顶层结构; (b) 底部馈线; (c) 实验样品顶层; (d) 实验样品底部馈线\r\nFig. 5. Reconfigurable metasurface sample structure with actual test sample: (a) Top structure; (b) bottom feeder; (c) top layer of experimental sample; (d) bottom feeder of experimental sample.
    图 5  可重构超表面样品结构与实际测试样品 (a) 顶层结构; (b) 底部馈线; (c) 实验样品顶层; (d) 实验样品底部馈线
    Fig. 5.  Reconfigurable metasurface sample structure with actual test sample: (a) Top structure; (b) bottom feeder; (c) top layer of experimental sample; (d) bottom feeder of experimental sample.
    图 6 实验测试环境\r\nFig. 6. Experimental test environment.
    图 6  实验测试环境
    Fig. 6.  Experimental test environment.

    入射为x极化波的实验测试结果如图7所示. 二极管状态为“10”时, 在频段3.80—4.22 GHz内极化转换率超过0.9, 且相对带宽为10.9%, 如图7(a)所示; 二极管状态为“01”时, 在4.25—5.0 GHz范围内实现线-圆极化转换, 且轴比小于3 dB的相对带宽为16.2%, 如图7(b)所示; 二极管状态为“00”时, 在5.8—7.3 GHz实现全反射功能, 且相对带宽为22.9%, 极化转换率小于0.1, 如图7(c)所示; 二极管状态为“11”时, 在7.35—8.10 GHz范围内的极化转换率超过0.9, 且相对带宽为9.1%, 如图7(d)所示. 相比4个功能的仿真结果, 可以清楚地看出实验测试的频谱发生蓝移, 幅值略有偏差, 但二者结果总体上吻合较好. 可认为造成二者误差的主要原因是加工误差, 一方面来自于介质基板的介电常数与正切损耗的误差; 另一方面来自于焊接二极管时引入的寄生电容与电感效应.

    图 7 二极管不同状态下实验与仿真的性能对比 (a) “10”; (b) “01”; (c) “00”; (d) “11”\r\nFig. 7. Comparison of experimental and simulated performance of diodes in different states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.
    图 7  二极管不同状态下实验与仿真的性能对比 (a) “10”; (b) “01”; (c) “00”; (d) “11”
    Fig. 7.  Comparison of experimental and simulated performance of diodes in different states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.

    一般情况下, 电磁波的极化转换特性采用本征模进行分析. 假设将x极化波入射到超表面, 此时, 入射电场Ei可以分解为u, v方向上相互正交的电场分量Eiu, Eiv; 同样, 反射电场Er分别在u, v方向上存在两个相互正交的电场分量EruErv, 仿真结果如图8所示. 从图8(a)可以看出, 共极化反射系数ruurvv在3.7—4.1 GHz内几乎相等且接近于1, 相应的反射相位差近似于π. 图8(d)中也显示了共极化反射系数ruurvv在7.3—8.0 GHz几乎相等且为1, 相应的相位差也近似为π. 因此, 这两个频段下将入射的x极化波转化为与其垂直的y极化波. 于是反射波的电场可表示为

    图 8 不同二极管状态下在u、v方向上的反射系数和相位差 (a) “10”; (b) “01”; (c) “00”; (d) “11”\r\nFig. 8. Reflectances and phase differences in u and v directions for different diode states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.
    图 8  不同二极管状态下在uv方向上的反射系数和相位差 (a) “10”; (b) “01”; (c) “00”; (d) “11”
    Fig. 8.  Reflectances and phase differences in u and v directions for different diode states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.
    Er=euruuEiuej(φ+φuu)+evrvvEivej(φ+φuuπ )=(euev)E0ej(φ+φuu)=yEoej(φ+φuu).
    (11)

    图8(b)可以看出, 共极化ruurvv的反射系数在4.25—5.0 GHz内保持相等且为1, 相位差大约为–π/2, 这意味着入射的线极化波转换成左旋圆极化波, 于是反射电场可表示为

    Er=euruuEiuej(φ+φuu)+evrvvEivej(φ+φuuπ/2)=(euEiu+evjEiv)ej(φ+φuu).
    (12)

    图8(c)可以看出, 共极化的反射系数在6.0—7.1 GHz几乎相等且接近1, 相应的相位差近似为0, 这意味着该频段下的超表面结构可以等效为完美电导体, 将入射的x极化波以相同的极化方式全反射. 于是反射波的电场可表示为:

    Er=euruuEiuej(φ+φuu)+evrvvEivej(φ+φvv)=(eu+ev)E0ej(φ+φuu)=xEoej(φ+φuu).
    (13)

    在实际的电磁环境中, 电磁波除了垂直入射到物体的表面, 还会存在任意角度入射的情况. 因此, 研究电磁波的斜入射特性对于超表面具有重要意义. 以TE波为例, 不同入射角对不同二极管状态下的极化特性影响结果如图9所示, 从图9(a)可以看出, 二极管处于“10”状态, 在3.7—4.16 GHz内斜入射角小于30°时, 极化转换效率仍能保持在0.83以上, 随着斜入射角的增大, 转换效率也逐渐降低, 但带宽略有增大. 图9(b)可以看出, 二极管处于“01”状态, 斜入射角小于45°且频率在4.25—4.95 GHz范围内, 轴比均小于3 dB, 这表明线-圆极化转换具有良好的广角特性, 这可归因于斜入射电场在水平和垂直方向上基本一致的分量和几乎没有变化的相位差. 图9(c)可以看出, 在6.02—7.16 GHz频带下二极管处于“00”状态, 斜入射角小于30°时转换率和吸收率均小于0.1, 这表明入射的线极化波也具有良好的广角特性, 随着角度的增大, 极化转化率和吸收率也在逐渐增大. 从图9(d)可以看出, 二极管处于“11”状态且斜入射角小于30°, 在7.30—8.15 GHz频带范围时极化转化率保持在0.83以上, 随着斜入射角的增大, 转换率急剧下降且带宽逐渐变窄, 这主要是因为有源切换层与金属地层之间磁耦合效应的减弱, 使得金属层之间的谐振逐渐降低.

    图 9 不同入射角对不同二极管状态下可重构极化调控性能的影响 (a) “10”; (b) “01”; (c) “00”; (d) “11”\r\nFig. 9. Effect of different incidence angles on the reconfigurable polarization modulation performance in different diode states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.
    图 9  不同入射角对不同二极管状态下可重构极化调控性能的影响 (a) “10”; (b) “01”; (c) “00”; (d) “11”
    Fig. 9.  Effect of different incidence angles on the reconfigurable polarization modulation performance in different diode states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.

    考虑到实际中电磁波入射时的方位角可能沿任意方向, 因此, 图10展示了不同二极管状态下时在其他方位角下的斜入射特性. 图10(a), (b)(d)可以看出, 二极管处于“10”, “01”与“11”状态, 方位角ϕ为15°时, 该结构的斜入射响应较为良好; 图10(c)展示了二极管处于“00”状态, 方位角ϕ为10°时, 该结构的斜入射响应较为良好. 由此也可知, 对于电磁波偏振沿其他方向入射时, 对该器件的性能有一定的影响.

    图 10 其他方位角的斜入射响应 (a) “10”状态; (b) “01”状态 (c) “00”状态; (d) “11”状态\r\nFig. 10. Oblique incident response at other azimuths: (a) “10”; (b) “01”; (c) “00”; (d) “11”.
    图 10  其他方位角的斜入射响应 (a) “10”状态; (b) “01”状态 (c) “00”状态; (d) “11”状态
    Fig. 10.  Oblique incident response at other azimuths: (a) “10”; (b) “01”; (c) “00”; (d) “11”.

    为了解释极化转换的物理机理, 仿真了不同二极管状态下各功能谐振点处的表面电流分布. 当二极管处于“10”状态时, 监测两个谐振点f1 = 3.85 GHz, f2 = 4.13 GHz的表面电流分布, 如图11(a)(b)所示. 当入射的电磁波为x极化时, 有源切换层的表面电流与金属地层表面电流方向相反, 使得感生磁场分量方向与原生磁场方向产生90°的偏转, 这也是入射的x极化波调控为y极化波的原因. 当二极管处于“01”状态时, 两个谐振点f1 = 4.50 GHz和f2 = 4.75 GHz的表面电流分布如图11(c)(d)所示. 有源切换层的表面电流和金属地层的表面电流分别等效为I1I2, 由图可知, 表面电流I1I2相互垂直, 此时输入阻抗呈高电感特性, 有源切换层与金属地层之间谐振引入了90°的相位差, 导致入射的线极化波转换为圆极化波, 同时, 多个谐振点的叠加拓宽了该功能的有效带宽. 当二极管处于“00”状态时, 图11(e)(f)显示了f1 = 6.30 GHz, f2 = 7.07 GHz时的表面电流分布, 有源切换层上的表面电流与金属地层的表面电流方向一致, 构成了电偶极子, 感生电场的方向与原生电场的方向保持一致, 这表明入射电场被反射后并没有发生变化, 可重构超表面则表现为全反射功能. 当二极管处于“11”状态时, 图11(g)(h)展示了f1 = 7.36 GHz, f2 = 7.89 GHz的表面电流分布, 这与二极管处于“10”状态时的物理机理一致, 也就是将入射的电场调控到与其相互垂直的方向上, 此时入射波的极化状态发生了变化. 综上所述, 不同二极管状态下有源切换层与金属地层表面电流分布的不同会导致各谐振点处耦合响应的不同, 主要是将入射波的电场调控到不同方向上, 从而表现为电磁波极化状态的转换.

    图 11 不同二极管状态下各极化转换功能在谐振点处的表面电流分布 (a) (b) “10”; (c) (d) “01”; (e) (f) “00”; (g) (h) “11”\r\nFig. 11. Surface current distribution of each polarization conversion function at the resonance point under different diode states: (a) (b) “10”; (c) (d) “01”; (e) (f) “00”; (g) (h) “11”.
    图 11  不同二极管状态下各极化转换功能在谐振点处的表面电流分布 (a) (b) “10”; (c) (d) “01”; (e) (f) “00”; (g) (h) “11”
    Fig. 11.  Surface current distribution of each polarization conversion function at the resonance point under different diode states: (a) (b) “10”; (c) (d) “01”; (e) (f) “00”; (g) (h) “11”.

    本文基于频域调控理论设计了一种极化与频率同时重构的多缝隙型超表面, 通过加载PIN二极管实现了不同频段内的不同极化转换功能, 3.80—4.22 GHz和7.35—8.10 GHz实现了极化转换, 极化转换效率均大于0.9; 4.25—5.0 GHz频带内实现了线-圆极化转换; 5.8—7.3 GHz实现了全反射功能. 对加工后的样品进行测试, 其结果与仿真结果吻合较好. 本文所提出的结构可进一步应用于空间波的极化调控和实现波束的控制, 同时, 可应用于新型数字通信系统.

    [1]

    Gao X, Han X, Cao W P, Li H O, Ma H F, Cui T J 2015 IEEE Trans. Antennas Propag. 63 3522Google Scholar

    [2]

    Huang X J, Yang H L, Zhang D H, Luo Y 2019 IEEE Trans. Antennas Propag. 67 4636Google Scholar

    [3]

    Li Y, Cao Q, Wang Y 2018 IEEE Antennas Wirel. Propag. Lett. 17 1314Google Scholar

    [4]

    Nama L, Nilotpal, Bhattacharyya S, Jain P K 2021 IEEE Antennas Propag. Mag. 63 100Google Scholar

    [5]

    Zheng Q, Guo C, Ding J 2018 IEEE Antennas Wirel. Propag. Lett. 17 1459Google Scholar

    [6]

    李海鹏, 吴潇, 丁海洋, 辛可为, 王光明, 徐进, 李荣强, 蒋小平, 王身云, 韩天成 2021 物理学报 70 027803Google Scholar

    Li H P, Wu X, Ding H Y, Xin K W, Wang G M, Xu J, Li R Q, Jiang X P, Wang S Y, Han T C 2021 Acta Phys. Sin. 70 027803Google Scholar

    [7]

    徐进, 李荣强, 蒋小平, 王身云, 韩天成 2019 物理学报 68 117081Google Scholar

    Xu J, Li R Q, Jiang X P, Wang S Y, Han T C 2019 Acta Phys. Sin. 68 117081Google Scholar

    [8]

    Li Y, Lin J, Guo H, Sun W, Xiao S, Zhou L 2020 Adv. Opt. Mater. 8 1901548Google Scholar

    [9]

    Shen Z Y, Huang X J, Yang H L, Xiang T Y, Wang C W, Yu Z T, Wu J 2018 J. Appl. Phys. 123 225106Google Scholar

    [10]

    Negm A, Bakr M, Howlader M, Ali S 2021 Smart Mater. Struct. 30 075011Google Scholar

    [11]

    高喜, 唐李光 2021 物理学报 70 038101Google Scholar

    Gao X, Tang L G 2021 Acta Phys. Sin. 70 038101Google Scholar

    [12]

    李国强, 施宏宇, 刘康, 李博林, 衣建甲, 张安学, 徐卓 2021 物理学报 70 188701Google Scholar

    Li G Q, Shi H N, Liu K, Li B L, Yi J J, Zhang A X, Xu Z 2021 Acta Phys. Sin. 70 188701Google Scholar

    [13]

    李小兵, 陆卫兵, 刘震国, 陈昊 2018 物理学报 67 184101Google Scholar

    Li X B, Lu W B, Liu Z G, Chen H 2018 Acta Phys. Sin. 67 184101Google Scholar

    [14]

    张娜, 赵健民, 陈克, 赵俊明, 姜田, 冯一军 2021 物理学报 70 178101Google Scholar

    Zhang N, Zhao J M, Chen K, Zhao J M, Jiang T, Feng Y J 2021 Acta Phys. Sin. 70 178101Google Scholar

    [15]

    Castaldi G, Zhang L, Moccia M, Hathaway A Y, Tang W X, Cui T J, Galdi V 2021 Adv. Funct. Mater. 31 2007620Google Scholar

    [16]

    马婧, 刘冬冬, 王继成, 冯延 2018 物理学报 67 094102Google Scholar

    Ma J, Liu D D, Wang J C, Feng Y 2018 Acta Phys. Sin. 67 094102Google Scholar

    [17]

    Sun W J, Yang W W, Guo L, Qin W, Chen J X 2020 IEEE Antennas Wirel. Propag. Lett. 19 1088Google Scholar

    [18]

    Xia Z X, Leung K W, Yang N, Lu K 2021 IEEE Trans. Antennas Propag. 69 2031Google Scholar

    [19]

    庞慧中, 王鑫, 王俊林, 王宗利, 刘苏雅拉图, 田虎强 2021 物理学报 70 168101Google Scholar

    Pang H Z, Wang X, Wang J L, Wang Z L, Liu S Y L T, Tian H Q 2021 Acta Phys. Sin. 70 168101Google Scholar

    [20]

    Wang H, Sui S, Li Y, Chen H, Wang J, Zhang J, Qu S 2020 Smart Mater. Struct. 29 015029Google Scholar

    [21]

    Tian J, Cao X, Gao J, Yang H, Han J, Yu H, Wang S, Jin R, Li T 2019 J. Appl. Phys. 125 135105Google Scholar

    [22]

    He J H, Wang S Q, Li X W, Fan J D, Guo L Y, Guo L, Huang X J 2021 Phys. Scr. 96 125846Google Scholar

    [23]

    Liu G Y, Li L, Han J Q, Liu H X, Gao X H, Shi Y, Cui T J 2020 ACS Appl. Mater. Interfaces 12 23554Google Scholar

    [24]

    Wang H L, Ma H F, Chen M, Sun S, Cui T J 2021 Adv. Funct. Mater. 31 2100275Google Scholar

    [25]

    Zhang N, Chen K, Zheng Y, Hu Q, Qu K, Zhao J, Wang J, Feng Y 2020 IEEE J. Em. Sel. Top. C 10 20

    期刊类型引用(4)

    1. 张胜,曹为为,解鑫,李大妮,于梦杰,程德强. 基于超表面的交叉极化转换器设计. 电子元件与材料. 2024(02): 223-229 . 百度学术
    2. 许丽佳,郑宇晗,郭迎辉,蒲明博,李雄,马晓亮,罗先刚. 超构表面偏振调控最新研究进展(特邀). 光学学报. 2024(10): 210-221 . 百度学术
    3. 曹为为,解鑫,李大妮,于梦杰,张胜. 基于超表面的多频带双极化转换器设计. 光电子·激光. 2024(07): 683-690 . 百度学术
    4. 于瑞涛,符道临,梁雄帆,王佳毅,陈展野,熊伟,陈珲. 基于超表面的宽波瓣阵列天线. 海军航空大学学报. 2024(06): 765-772 . 百度学术

    其他类型引用(1)

  • 图 1  单元结构示意图 (a) 立体图; (b) 主视图; (c) 侧视图

    Fig. 1.  Schematics of unit cell: (a) Solid Shape; (b) main view; (c) side view.

    图 2  等效电路 (a) 等效电路分析; (b) 等效电路模型

    Fig. 2.  Equivalent circuit: (a) Equivalent circuit analysis; (b) equivalent circuit model.

    图 3  二极管不同状态下的仿真曲线 (a) “10”-反射系数; (b) “01”-反射系数; (c) “01”-反射相位; (d) “10”-反射系数; (e) “11”-反射系数

    Fig. 3.  Simulation curve of diodes in different states: (a) “10”-reflectance; (b) “01”-reflectance; (c) “01”-reflection phase; (d) “10”- reflectance; (e) “11”- reflectance.

    图 4  二极管不同状态下的极化转换率、吸收率和轴比 (a) “10”与“01”; (b) “00”与“11”

    Fig. 4.  Polarization conversion rate, absorptance and axial ratio of diodes in different states: (a) “10” and “01”, (b) “00” and “11”.

    图 5  可重构超表面样品结构与实际测试样品 (a) 顶层结构; (b) 底部馈线; (c) 实验样品顶层; (d) 实验样品底部馈线

    Fig. 5.  Reconfigurable metasurface sample structure with actual test sample: (a) Top structure; (b) bottom feeder; (c) top layer of experimental sample; (d) bottom feeder of experimental sample.

    图 6  实验测试环境

    Fig. 6.  Experimental test environment.

    图 7  二极管不同状态下实验与仿真的性能对比 (a) “10”; (b) “01”; (c) “00”; (d) “11”

    Fig. 7.  Comparison of experimental and simulated performance of diodes in different states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.

    图 8  不同二极管状态下在uv方向上的反射系数和相位差 (a) “10”; (b) “01”; (c) “00”; (d) “11”

    Fig. 8.  Reflectances and phase differences in u and v directions for different diode states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.

    图 9  不同入射角对不同二极管状态下可重构极化调控性能的影响 (a) “10”; (b) “01”; (c) “00”; (d) “11”

    Fig. 9.  Effect of different incidence angles on the reconfigurable polarization modulation performance in different diode states: (a) “10”; (b) “01”; (c) “00”; (d) “11”.

    图 10  其他方位角的斜入射响应 (a) “10”状态; (b) “01”状态 (c) “00”状态; (d) “11”状态

    Fig. 10.  Oblique incident response at other azimuths: (a) “10”; (b) “01”; (c) “00”; (d) “11”.

    图 11  不同二极管状态下各极化转换功能在谐振点处的表面电流分布 (a) (b) “10”; (c) (d) “01”; (e) (f) “00”; (g) (h) “11”

    Fig. 11.  Surface current distribution of each polarization conversion function at the resonance point under different diode states: (a) (b) “10”; (c) (d) “01”; (e) (f) “00”; (g) (h) “11”.

  • [1]

    Gao X, Han X, Cao W P, Li H O, Ma H F, Cui T J 2015 IEEE Trans. Antennas Propag. 63 3522Google Scholar

    [2]

    Huang X J, Yang H L, Zhang D H, Luo Y 2019 IEEE Trans. Antennas Propag. 67 4636Google Scholar

    [3]

    Li Y, Cao Q, Wang Y 2018 IEEE Antennas Wirel. Propag. Lett. 17 1314Google Scholar

    [4]

    Nama L, Nilotpal, Bhattacharyya S, Jain P K 2021 IEEE Antennas Propag. Mag. 63 100Google Scholar

    [5]

    Zheng Q, Guo C, Ding J 2018 IEEE Antennas Wirel. Propag. Lett. 17 1459Google Scholar

    [6]

    李海鹏, 吴潇, 丁海洋, 辛可为, 王光明, 徐进, 李荣强, 蒋小平, 王身云, 韩天成 2021 物理学报 70 027803Google Scholar

    Li H P, Wu X, Ding H Y, Xin K W, Wang G M, Xu J, Li R Q, Jiang X P, Wang S Y, Han T C 2021 Acta Phys. Sin. 70 027803Google Scholar

    [7]

    徐进, 李荣强, 蒋小平, 王身云, 韩天成 2019 物理学报 68 117081Google Scholar

    Xu J, Li R Q, Jiang X P, Wang S Y, Han T C 2019 Acta Phys. Sin. 68 117081Google Scholar

    [8]

    Li Y, Lin J, Guo H, Sun W, Xiao S, Zhou L 2020 Adv. Opt. Mater. 8 1901548Google Scholar

    [9]

    Shen Z Y, Huang X J, Yang H L, Xiang T Y, Wang C W, Yu Z T, Wu J 2018 J. Appl. Phys. 123 225106Google Scholar

    [10]

    Negm A, Bakr M, Howlader M, Ali S 2021 Smart Mater. Struct. 30 075011Google Scholar

    [11]

    高喜, 唐李光 2021 物理学报 70 038101Google Scholar

    Gao X, Tang L G 2021 Acta Phys. Sin. 70 038101Google Scholar

    [12]

    李国强, 施宏宇, 刘康, 李博林, 衣建甲, 张安学, 徐卓 2021 物理学报 70 188701Google Scholar

    Li G Q, Shi H N, Liu K, Li B L, Yi J J, Zhang A X, Xu Z 2021 Acta Phys. Sin. 70 188701Google Scholar

    [13]

    李小兵, 陆卫兵, 刘震国, 陈昊 2018 物理学报 67 184101Google Scholar

    Li X B, Lu W B, Liu Z G, Chen H 2018 Acta Phys. Sin. 67 184101Google Scholar

    [14]

    张娜, 赵健民, 陈克, 赵俊明, 姜田, 冯一军 2021 物理学报 70 178101Google Scholar

    Zhang N, Zhao J M, Chen K, Zhao J M, Jiang T, Feng Y J 2021 Acta Phys. Sin. 70 178101Google Scholar

    [15]

    Castaldi G, Zhang L, Moccia M, Hathaway A Y, Tang W X, Cui T J, Galdi V 2021 Adv. Funct. Mater. 31 2007620Google Scholar

    [16]

    马婧, 刘冬冬, 王继成, 冯延 2018 物理学报 67 094102Google Scholar

    Ma J, Liu D D, Wang J C, Feng Y 2018 Acta Phys. Sin. 67 094102Google Scholar

    [17]

    Sun W J, Yang W W, Guo L, Qin W, Chen J X 2020 IEEE Antennas Wirel. Propag. Lett. 19 1088Google Scholar

    [18]

    Xia Z X, Leung K W, Yang N, Lu K 2021 IEEE Trans. Antennas Propag. 69 2031Google Scholar

    [19]

    庞慧中, 王鑫, 王俊林, 王宗利, 刘苏雅拉图, 田虎强 2021 物理学报 70 168101Google Scholar

    Pang H Z, Wang X, Wang J L, Wang Z L, Liu S Y L T, Tian H Q 2021 Acta Phys. Sin. 70 168101Google Scholar

    [20]

    Wang H, Sui S, Li Y, Chen H, Wang J, Zhang J, Qu S 2020 Smart Mater. Struct. 29 015029Google Scholar

    [21]

    Tian J, Cao X, Gao J, Yang H, Han J, Yu H, Wang S, Jin R, Li T 2019 J. Appl. Phys. 125 135105Google Scholar

    [22]

    He J H, Wang S Q, Li X W, Fan J D, Guo L Y, Guo L, Huang X J 2021 Phys. Scr. 96 125846Google Scholar

    [23]

    Liu G Y, Li L, Han J Q, Liu H X, Gao X H, Shi Y, Cui T J 2020 ACS Appl. Mater. Interfaces 12 23554Google Scholar

    [24]

    Wang H L, Ma H F, Chen M, Sun S, Cui T J 2021 Adv. Funct. Mater. 31 2100275Google Scholar

    [25]

    Zhang N, Chen K, Zheng Y, Hu Q, Qu K, Zhao J, Wang J, Feng Y 2020 IEEE J. Em. Sel. Top. C 10 20

  • [1] 杨国, 施鸿强, 黎小聪, 肖如奇, 齐世山, 吴文. 大频率比的毫米波频率可重构滤波天线. 物理学报, 2025, 74(1): 018401. doi: 10.7498/aps.74.20241494
    [2] 王哲飞, 吴杰, 万发雨, 曾庆生, 侯建强, 傅佳辉, 吴群, 宋明歆, TayebA. Denidni. 基于类电磁诱导透明效应的极化转换滤波器. 物理学报, 2024, 73(18): 188101. doi: 10.7498/aps.73.20240632
    [3] 夏兆生, 刘宇行, 包正, 王丽华, 吴博, 王刚, 王辉, 任信钢, 黄志祥. 基于准连续域束缚态的强圆二色性超表面. 物理学报, 2024, 73(17): 178102. doi: 10.7498/aps.73.20240834
    [4] 张向, 王玥, 张婉莹, 张晓菊, 罗帆, 宋博晨, 张狂, 施卫. 单壁碳纳米管太赫兹超表面窄带吸收及其传感特性. 物理学报, 2024, 73(2): 026102. doi: 10.7498/aps.73.20231357
    [5] 王丹, 李九生, 郭风雷. 宽带吸收与极化转换可切换的太赫兹超表面. 物理学报, 2024, 73(14): 148701. doi: 10.7498/aps.73.20240525
    [6] 范辉颖, 罗杰. 非厄密电磁超表面研究进展. 物理学报, 2022, 71(24): 247802. doi: 10.7498/aps.71.20221706
    [7] 龙洁, 李九生. 相变材料与超表面复合结构太赫兹移相器. 物理学报, 2021, 70(7): 074201. doi: 10.7498/aps.70.20201495
    [8] 王明照, 王少杰, 许河秀. 基于剪纸方法的一种可重构线极化转换空间序构超表面. 物理学报, 2021, 70(15): 154101. doi: 10.7498/aps.70.20210188
    [9] 郭泽旭, 曹祥玉, 高军, 李思佳, 杨欢欢, 郝彪. 一种复合型极化转换表面及其在天线辐射散射调控中的应用. 物理学报, 2020, 69(23): 234102. doi: 10.7498/aps.69.20200797
    [10] 谢智强, 贺炎亮, 王佩佩, 苏明样, 陈学钰, 杨博, 刘俊敏, 周新星, 李瑛, 陈书青, 范滇元. 基于Pancharatnam-Berry相位超表面的二维光学边缘检测. 物理学报, 2020, 69(1): 014101. doi: 10.7498/aps.69.20191181
    [11] 李晓楠, 周璐, 赵国忠. 基于反射超表面产生太赫兹涡旋波束. 物理学报, 2019, 68(23): 238101. doi: 10.7498/aps.68.20191055
    [12] 徐进, 李荣强, 蒋小平, 王身云, 韩天成. 基于方形开口环的超宽带线性极化转换器. 物理学报, 2019, 68(11): 117801. doi: 10.7498/aps.68.20190267
    [13] 于惠存, 曹祥玉, 高军, 杨欢欢, 韩江枫, 朱学文, 李桐. 一种宽带可重构反射型极化旋转表面. 物理学报, 2018, 67(22): 224101. doi: 10.7498/aps.67.20181041
    [14] 李唐景, 梁建刚, 李海鹏, 牛雪彬, 刘亚峤. 基于单层线-圆极化转换聚焦超表面的宽带高增益圆极化天线设计. 物理学报, 2017, 66(6): 064102. doi: 10.7498/aps.66.064102
    [15] 李唐景, 梁建刚, 李海鹏. 基于单层反射超表面的宽带圆极化高增益天线设计. 物理学报, 2016, 65(10): 104101. doi: 10.7498/aps.65.104101
    [16] 郭文龙, 王光明, 李海鹏, 侯海生. 单层超薄高效圆极化超表面透镜. 物理学报, 2016, 65(7): 074101. doi: 10.7498/aps.65.074101
    [17] 范亚, 屈绍波, 王甲富, 张介秋, 冯明德, 张安学. 基于交叉极化旋转相位梯度超表面的宽带异常反射. 物理学报, 2015, 64(18): 184101. doi: 10.7498/aps.64.184101
    [18] 余积宝, 马华, 王甲富, 冯明德, 李勇峰, 屈绍波. 基于开口椭圆环的高效超宽带极化旋转超表面. 物理学报, 2015, 64(17): 178101. doi: 10.7498/aps.64.178101
    [19] 李勇峰, 张介秋, 屈绍波, 王甲富, 吴翔, 徐卓, 张安学. 圆极化波反射聚焦超表面. 物理学报, 2015, 64(12): 124102. doi: 10.7498/aps.64.124102
    [20] 陈谦, 江建军, 别少伟, 王鹏, 刘鹏, 徐欣欣. 含有源频率选择表面可调复合吸波体. 物理学报, 2011, 60(7): 074202. doi: 10.7498/aps.60.074202
  • 期刊类型引用(4)

    1. 张胜,曹为为,解鑫,李大妮,于梦杰,程德强. 基于超表面的交叉极化转换器设计. 电子元件与材料. 2024(02): 223-229 . 百度学术
    2. 许丽佳,郑宇晗,郭迎辉,蒲明博,李雄,马晓亮,罗先刚. 超构表面偏振调控最新研究进展(特邀). 光学学报. 2024(10): 210-221 . 百度学术
    3. 曹为为,解鑫,李大妮,于梦杰,张胜. 基于超表面的多频带双极化转换器设计. 光电子·激光. 2024(07): 683-690 . 百度学术
    4. 于瑞涛,符道临,梁雄帆,王佳毅,陈展野,熊伟,陈珲. 基于超表面的宽波瓣阵列天线. 海军航空大学学报. 2024(06): 765-772 . 百度学术

    其他类型引用(1)

计量
  • 文章访问数:  6641
  • PDF下载量:  243
  • 被引次数: 5
出版历程
  • 收稿日期:  2022-06-27
  • 修回日期:  2022-07-27
  • 上网日期:  2022-11-11
  • 刊出日期:  2022-11-20

/

返回文章
返回