Research on Terahertz Band Demultiplexer Based on Metamaterials
-
摘要: 本文提出了一种基于超材料的太赫兹分波器，其结构单元由上下两层“一”字型金属线以及中间的石英介质构成。此分波器工作在太赫兹通信窗口，将中心频率为0.225 THz和0.300THz的两束太赫兹波分离，其隔离度分别为34 dB、47 dB，插入损耗分别为0.19 dB、0.04 dB，分波器的群延迟稳定。另外，也对分波器的等效参数以及结构参数变化进行了仿真分析。结果表明，所设计的太赫兹分波器具有插入损耗小、隔离度大、群延迟稳定等特点，且性能指标受参数变化的影响较小，在太赫兹通信系统中具有较好的应用前景。Abstract: This paper proposes a terahertz demultiplexer based on metamaterials; its unit cell is composed of two layers of "one" shaped metal wires and a quartz substrate in the middle. The terahertz demultiplexer works in the terahertz communication window and demultiplexes two beams of terahertz waves with center frequencies of 0.225 THz and 0.300 THz. The isolation of the terahertz demultiplexer is 34 dB and 47 dB, and its insertion loss is 0.19 dB and 0.04 dB, respectively. The group delay of the demultiplexer was relatively stable. Finally, the influence of the effective and structural parameters of the demultiplexer on its performance are discussed. This terahertz demultiplexer has low insertion loss, large isolation, and stable group delay, and its performance is less affected by parameter changes. It has significant application prospects in practical terahertz communication systems.
-
Keywords:
- terahertz /
- demultiplexer /
- metamaterial
-
0. 引言
编织复合材料是由碳纤维材料为增强体，以环氧树脂为基质，由不同方向纤维铺层叠加组成的新型复合材料，其具有高比强度、高比刚度、耐腐蚀、耐老化、耐热性的优点。编织复合材料由于制造成本低、生产效率高、机械性能突出、具有纤维互锁的网状结构[1]。编织复合材料在制造和适用过程中会碳纤维复合材料实现轻量化的同时，还能提供产品相当的安全性和可靠性，已被广泛应用于航空航天、汽车、轮船、运动器械等[2-4]。NASA喷气发动机风扇机匣、波音787飞机风扇机匣都采用的是编织复合材料。复合材料结构是由纤维铺层叠加而成，在面临碰撞和冲击时，纤维和纤维之间都是由基体传递载荷，冲击容易造成材料基体开裂、分层、脱粘等缺陷，编织层压板冲击损伤模式主要是层间分层，不会出现纤维劈裂等损伤模式[5]。在结构表面看不出来损伤缺陷，目视检测也无法检测出缺陷信息，但其内部可能因瞬时冲击能量的吸收，产生了严重的损伤，使得材料结构失效、材料性能降低，对结构承载能力造成大幅度下降[6-7]。因此，有效及时检测编织复合材料内部损伤或缺陷情况，对于保证复合材料结构件可靠性具有重要的意义。
超声红外热成像技术以超声作为热激励源，激励过程中仅界面类缺陷或闭合裂纹发热，正常区域不发热，属于选择性激励方式，更适合于不规则形状材料界面缺陷或闭合裂纹检测。因此目前超声红外热成像技术主要应用在对金属裂纹的检测[8-11]。韩梦[12]等采用超声红外热成像对铝金属的疲劳微裂纹进行检测，微裂纹最小宽度为20 μm；习小文[13]等采用超声红外热成像技术对航空发动机叶片裂纹进行检测，最小检测裂纹宽度为0.5 μm。G. Busse等[14]最早将超声红外锁相热成像检测技术引入到复合材料界面类缺陷检测。田干等[15]采用超声红外热成像技术对航空复合材料进行数值仿真研究，建立含裂纹缺陷复合材料的有限元模型；金国锋、张炜等[16-17]通过数值计算和试验研究了超声红外热成像技术对复合材料冲击损伤检测的适用性；李胤等[18]采用超声红外热像检测对复合材料层合板24 J和29 J冲击损伤进行检测研究，并与超声C扫进行对比，结果表明超声红外热像检测可检测出复合材料不同损伤形式，并且比超声C扫检测速度快、检测精度高。
目前对冲击能量损伤研究大部分基于1个或2个冲击能量对复合材料损伤的研究，为了解不同冲击能量对编织复合材料损伤及扩展情况，采用10 J、20 J、30 J、40 J、50 J冲击能量对不同编织复合材料进行冲击，并对冲击过程中的冲击能量-时间曲线进行了记录，采用超声红外热成像技术对不同冲击能量损伤试件进行检测与分析。
1. 基本原理
1.1 超声热成像技术的检测原理
超声红外热成像技术是将超声波激励和红外热成像两者结合的新型无损检测方法，属于主动式红外无损检测方法的一种。超声红外热成像技术的检测原理如图 1所示，是通过将高频超声波（20~40 kHz）耦合注入待测试件，超声波经过界面耦合在试件中传播，在超声波的作用下，试件缺陷由于塑性变形、摩擦作用、粘弹性效应等产生热量，使得裂纹、损伤区域温度升高，在热传导作用下，其裂纹或损伤位置对应的试件表面温度也相应的升高，利用红外热像仪采集试件表面的红外热序列图像，对红外序列图像进行处理便可得到其内部缺陷信息。相较于其他激励方式的红外热成像技术，例如：闪光灯激励、红外灯激励、激光激励等，这些激励方式都属于面激励，不论是缺陷区域还是无缺陷区域，试件表面被整体加热，表面吸收的能量通过热传导到试件内部，遇到隔热性缺陷时，热波反射回表面。超声红外热成像技术是对试件局部加热，其仅对缺陷的区域进行加热，无缺陷区域不进行加热，属于选择性加热方式，在超声激励下产生相对位移的缺陷区域温度都将升高，相当于是体加热方式，并且从缺陷区域通过热传导的方式传播到试件表面，仅传播了单次的距离。
1.2 超声热成像技术的检测原理
超声波在试件内界面缺陷处衰减产生热量，缺陷处则成为热源。由此可知，超声波激励产生的热流率q的传播满足热传导方程：
$$ {k_x}\frac{{{\partial ^2}T}}{{\partial {x^2}}} + {k_y}\frac{{{\partial ^2}T}}{{\partial {y^2}}} + {k_z}\frac{{{\partial ^2}T}}{{\partial {z^2}}}{\text{ + }}q{\text{ = }}\rho c\frac{{\partial T}}{{\partial t}} $$ (1) 式(1)中：T为温度，℃；kx、ky、kz分别为在三个维度上的导热率，W/(m·℃)；ρ为密度，kg/m3，c为比热容，J/(kg·℃)。
热传导方程(1)的初始条件为：
$$ T(x, y, z, t)＝T_{0}\qquad t＝0$$ (2) 式(2)中：T0为环境温度，℃。
热传导方程支配热传导及其扩散过程，试件检测过程时间比较短，表面温度变化比较小，忽视试件表面对流和热辐射，则：
$$ \left\{ \begin{array}{l} \frac{{\partial T}}{{\partial t}} - \alpha {\nabla ^2} = q(t) \hfill \\ T(0) = {T_0} \hfill \\ \end{array} \right. $$ (3) 式(3)中：α为热扩散率，m2/s；∇2为拉普拉斯算子，
$$ {{\mathit{\nabla}} ^2}{\text{ = }}\frac{{{\partial ^2}}}{{\partial {x^2}}} + \frac{{{\partial ^2}}}{{\partial {y^2}}} + \frac{{{\partial ^2}}}{{\partial {z^2}}} $$ 分别对上式进行关于空间(x, y, z)的傅里叶变换得到：
$$ \left\{ \begin{array}{l}\frac{{\rm{d}}{T}^{f}}{{\rm{d}}t}-\alpha {ϵ}^{2}{T}^{f}={q}^{f}(ϵ, t)\\ {T}^{f}(ϵ, 0)={T}_{0}{}^{f}\end{array} \right. $$ (4) 式(4)中：qf为对生热率q作温度T的傅里叶变换；$ ϵ $为傅里叶变换过程中的中间变量，没有特殊的物理含义；Tf为对温度T的傅里叶变换。
常微分方程的解析解为：
$$ {T}^{f}(ϵ, t)={T}_{0}{}^{f}{{\rm{e}}}^{-\alpha {ϵ}^{2}t}+{\displaystyle {\int }_{0}^{t}{q}^{f}}(ϵ, t){{\rm{e}}}^{-\alpha {ϵ}^{2}(t-\tau )}{\rm{d}}\tau $$ (5) 对Tf(∈, t)进行傅里叶反变换得到：
$$ \begin{array}{l}T\left(X, t\right)={T}_{0}{{\displaystyle \int }}_{{R}^{3}x}^{\text{ }}K(X-ϵ, t){\rm{d}}ϵ+\\ {{\displaystyle \int }}_{0}^{t}{\rm{d}}\tau {{\displaystyle \int }}_{{R}^{3}x}^{\text{ }}K\left(X-ϵ, t-\tau \right){q}^{f}(ϵ, \tau ){\rm{d}}ϵ\end{array} $$ (6) 式(6)为超声波在试件内界面缺陷热源公式，式子表示，热源公式是一个复杂函数表达式。
2. 试验与结果分析
2.1 碳纤维复合材料冲击损伤
试验用的碳纤维复合材料是T700编织复合材料层压板，编织技术制造网格状堆叠而成，铺层方式为－45/90/45/0，经过浸渍固化直接形成复合材料结构，试件共5件，试件尺寸均为150 mm×100 mm×4 mm，对5件试件进行冲击试验，冲击能量分别是10 J、20 J、30 J、40 J、50 J，冲击能量-时间曲线如图 2所示，冲击后试件光学图像如图 3所示，图 3(a)为试件正面，为冲击面，图 3(b)为试件背面，光学图像上很难看到试件表面有损伤。
2.2 超声热成像检测系统
超声红外热成像系统如图 4所示，系统主要包括超声激励电源、超声枪、红外热像仪、夹具、图像处理和控制系统。超声枪包括换能器、调幅器、工具杆，最大功率为1200 W，激励频率为20 kHz，超声激励时间可设置为0.2~20 s。通过超声激励电源箱控制器调节输出功率大小，预紧力大小通过机械控制，预紧力大小最大可设置1000 N。红外热像仪型号为FLIR A655sc，属于非制冷型红外热像仪，工作波段为7.5~14 μm，室温下等效噪声温差为0.04 K，图像分辨率为640 pixel×480 pixel。
试验过程中超声枪不直接作用于试件表面，试件固定在夹具上，夹具固定在平台上，超声枪顶在夹具上面，在试验过程中，需要确保超声波有效注入到夹具中，夹具和试件要接触良好，超声枪在20 kHz的激励下，振幅约为几十微米，如果超声枪与夹具之间接触不好，则超声能量无法传递到夹具上，因此设置超声枪与夹具的预紧力为200 N，预紧力是通过旋转预紧力手柄迫使预紧力系统单元向前移动，内部弹簧受到挤压，从而对超声枪头产生压力，超声枪头压紧夹具。校正好红外热像仪，超声激励时间为1.5 s，采集频率为50 Hz，采集时间为10 s。
2.3 检测结果分析
试验中共有5种不同冲击能量的试件，选取30 J能量冲击能量的试件3#进行分析，图 5所示为超声红外热成像系统检测试件3#所示的红外序列热图，从图 5中可以看出，超声激励前，看不到任何热斑信息，随着超声持续激励，图像中热斑信息逐渐清晰，0.08 s时稍微可以看到零星散点，到0.5 s时逐渐形成了较为清晰的热斑，超声激励时间为1.5 s，红外序列中时间为1.6 s时热斑最为清晰，在超声激励结束后，冲击损伤区域不能持续发热，仅能靠自身热传导把热斑区域的能量传递到低温区域，逐渐达到温度平衡，从红外序列图中表现为图像越来越模糊，对比度越来越低，到8.66 s时已不能分辨出热斑轮廓信息，随着热传导的进行，将不能检测到任何热斑信息，所检测到的红外图像与0.02 s时图像一致。
图 6为不同位置处A点、B点、C点、D点、E点的温升变化曲线，由图 6曲线可以看出，超声激励后，损伤区域温度迅速上升，且各点温度上升不尽相同，A点为试件3#的温升最高点，最高温升可达到8℃，B点和C点为其他区域温升点，这3个点温升曲线有一个特点，温升曲线都是呈现先上升再下降的过程，这3个点代表着试件内部有持续发热的热源，这些热源就是在超声激励下，损伤处由于摩擦生热等原因产生的热源。D点温升曲线只有温升的过程，没有下降的过程，从图像上看，C点位热斑区中间位置，D点处于C点上方，D点处于热斑区域的边缘，边缘位置损伤比较轻微，温升没有那么高，但依然会在持续激励过程中温度不断上升，等热激励结束以后，随着热扩散的进行温度将下降，由于其处于热斑的边缘区域，热斑处的热量会由于热扩散向低温区域传播，C点和D点的温度由于热扩散会逐渐接近，因此导致D点温度下降曲线基本是平的。E点代表着无热斑区域，即无损伤区域，温升曲线基本是平坦的，E点对应位置的内部没有发热源，与上文分析一致。
为了更形象描述不同帧之间空间温度变化曲线，分别取图 5所示前6个时间点的对应图像的空间温度变化曲线，直线所对应的位置为图 7所示红色空间直线对应的位置，左边是起始点，对应的距离为0，右边对应的距离为141.5 mm。曲线横坐标为距离，纵坐标为温度，环境温度约为22℃，t＝0.02 s对应的空间曲线是超声激励之前的空间温度曲线，温度基本相同，没有太明显变化，t＝0.08 s红外图像上可以看到零星散点，红色空间直线正好位于散点上，在空间温度曲线上正好可以看到损伤处对应位置处有凸起，随着超声激励的加载，损伤区域空间温度曲线温度逐渐升高，并且其对应的宽度稍微有所增加，但温度升高的点仅限于损伤区域，非损伤区域对应的曲线基本是不变的，重叠在一起；损伤区域t＝0.24 s对应的空间温度曲线就比t＝0.08 s对应的空间温度曲线要高，同样的损伤区域t＝0.50 s对应的空间温度曲线比t＝0.24 s对应的空间温度曲线高，到超声激励结束时间附近t＝1.6 s，空间温度曲线升高最高值，图 7对应的最高曲线。超声激励结束后，损伤区域不再主动发热，内部热源消失，温度开始下降，空间温度曲线也随之下降，t＝2.28 s损伤区域对应的曲线的值就比t＝1.6 s对应曲线值有所下降。图像中红色箭头指示的空间距离约60 mm时间为t＝0.50 s和t＝1.60 s损伤区域的边缘点，温度比其旁边两个尖峰都要低，当超声激励结束以后，随着热扩散的进行，t＝2.28 s对应的温度反而要比t＝1.60 s对应的温度要高，此时该区域对应的2个尖峰也变成了一个大而平滑的尖峰曲线，在红外序列图像中就表现为图像变模糊，与图 5所示图像相对应。
2.4 不同冲击能量试件结果与分析
图 8为不同冲击能量冲击试件的检测图像，从光学图像上都看不到表面的损伤信息；从检测图上看，冲击能量为10 J的1#试件检测图上有3个小的热斑区域，冲击能量为20 J的2#试件看出热斑区域连成了一个整体，冲击能量为30 J的3#试件可以看到热斑已延展开，冲击能量为40 J的4#试件检测图上看出形成的热斑延展成一条长的缝，冲击能量为50 J的5#试件检测热图上热斑几乎贯穿对角线，并且还有额外2条分支。试件是二维编织技术制作的二维编织复合材料，编织使得纤维之间连成了有机的整体，使得二维编织复合材料展现出了较强的整体性。在冲击能量的作用下，外表面几乎看不出损伤信息，当编织复合材料受到小的外部冲击后，不足以破坏有机的整体，损伤区域仅限于冲击点周围；当编织复合材料受到大的外部冲击后，其内部将会产生大的面积分层，损伤展现出延展性。
2.5 损伤区域识别
衡量超声红外热成像技术对编织复合材料检测能力关键是对损伤区域进行提取，采用双阈值分割，对图 5中时间为t＝0.50 s、t＝1.60 s、t＝2.28 s图像进行分割，分割效果如图 9所示，无法达到理想的分割效果。底部区域是超声激励中夹具位置，夹具在超声激励作用下温度升高，通过热传导会把试件底部的温度升高，采用双阈值分割方法无法把该区域给剔除，对比3张图分割效果可以看出，选用不同时间的图像，分割的结果不尽相同。
由于红外序列图中单帧图像由于热传导一直在发生变化，采用单帧图像分割算法会由于选取的帧不同而导致分割区域的差异，为了准确提取损伤区域，本文采用曲线分类算法，利用每一个像素点的温升曲线进行分类，定义损伤区域强弱程度，温升曲线有上升沿和下降沿的点为热斑中心区域，该区域损伤比较严重，图 6所示对应的A、B、C三个点，该像素点置红色；温升曲线有上升沿但无下降沿的点为损伤比较轻微，图 6所示对应的D点，该像素点置绿色；温升曲线基本不变化的为无损伤区域，图 6所示的E点，该像素点置黑色。不同冲击能量试件分割效果如图 10所示。对比于双阈值分割，采用曲线分类算法分割效果比较好，且图像中有两种颜色表示，红色区域基本对应于检测图像中的亮热斑区域。
![]()
图 10 不同冲击能量试件分割效果Figure 10. Effect of splitting specimens with different impact energies表 1为不同冲击能量试件损伤面积统计，从表格中可以看出，冲击能量越大，损伤面积越大，但是试件4#损伤面积反而比试件3#损伤面积小。注意到在冲击试验中，图 2所示冲击能量-时间曲线，试件4#峰值冲击能量是40 J，正常冲击结束处的能量是30 J左右，但冲击结束处的能量只有20 J，反而低于试件3#冲击结束处的能量24 J，所以造成这种对应关系不匹配的情况。
表 1 不同冲击能量试件损伤面积Table 1. Damage area of specimens with different impact energySpecimen Red area/mm2 Total area/mm2 1#-10J 630 926 2#-20J 1475 2525 3#-30J 2909 3929 4#-40J 2896 3383 5#-50J 4517 6291 由于试件4#冲击结束处能量偏低，因此把其他4个试件的冲击能量和损伤面积作为损伤面积关系图，如图 11所示，可以发现损伤面积和冲击能量成线性关系，与碳纤维层合板冲击损伤面积和冲击能量关系曲线相同[19]。
![]()
图 11 编织复合材料损伤面积与冲击能量的关系Figure 11. Damage area of braided composite versus impact energy3. 结论
编织复合材料冲击损伤主要涉及复合材料的分层损伤，这些损伤无法采用目视检测，且损伤会使得结构承载能力严重降低，从而威胁复合材料构件的安全使用。制作了编织复合材料层压板试件共5件，并对不同冲击能量下冲击损伤开展了超声红外热成像技术研究。超声激励试件时，冲击损伤区域由于摩擦等因素产生热量，相当于内部有热源的存在；无损伤区域不生热，使用红外热像仪采集红外图像序列便可实现对损伤区域检测，本文以试件3#为例对超声激励过程的温升曲线、空间温度曲线进行了详细的分析。并对10~50 J冲击能量试件进行了对比分析，最后采用曲线分类算法对损伤区域进行了定量识别，发现编织复合材料损伤面积和冲击能量成线性关系。后续的工作是采用其他技术手段（如超声C扫描无损检测技术）确认试件损伤区域面积，以及对损伤区域的评估，将对损伤定量化及曲线分类算法得到更进一步提高。
-
[1] 赵亚军, 郁光辉, 徐汉青. 6G移动通信网络: 愿景、挑战与关键技术[J]. 中国科学: 信息科学, 2019, 49(8): 963-987. https://www.cnki.com.cn/Article/CJFDTOTAL-PZKX201908003.htm ZHAO Yajun, YU Guanghui, XU Hanqing. 6G mobile communication network: vision, challenges and key technologies[J]. Science in China: Information Science, 2019, 49(8): 963-987. https://www.cnki.com.cn/Article/CJFDTOTAL-PZKX201908003.htm
[2] 陈亮, 余少华. 6G移动通信发展趋势初探[J]. 光通信研究, 2019, 45(4): 1-2. https://www.cnki.com.cn/Article/CJFDTOTAL-GTXY201904001.htm XU Liang, YU Shaohua. Preliminary study on the development trend of 6G mobile communication[J]. Optical Communication Research, 2019, 45(4): 1-2. https://www.cnki.com.cn/Article/CJFDTOTAL-GTXY201904001.htm
[3] 莫秋燕, 吴家隐. 波分复用系统的研究综述[J]. 信息通信, 2019(7): 191-192. DOI: 10.3969/j.issn.1673-1131.2019.07.087 MO Qiuyan, WU Jiayin. Review of research on wavelength division multiplexing system[J]. Information and Communication, 2019(7): 191-192. DOI: 10.3969/j.issn.1673-1131.2019.07.087
[4] 孟宪浩. 密集波分复用系统中合波器/分波器的测试[J]. 电气化铁道, 2002(1): 41-42. https://www.cnki.com.cn/Article/CJFDTOTAL-DQHD200201014.htm MENG Xianhao. Test of multiplexer/demultiplexer in dense wavelength division multiplexing system[J]. Electrified Railway, 2002(1): 41-42. https://www.cnki.com.cn/Article/CJFDTOTAL-DQHD200201014.htm
[5] Yata M, Fujita M, Nagatsuma T. Photonic-crystal diplexers for terahertz- wave applications[J]. Optics Express, 2016, 24(7): 7835-7849. DOI: 10.1364/OE.24.007835
[6] LI Jiusheng, HAN L, LE Z. Compact four-channel terahertz demultiplexer based on directional coupling photonic crystal[J]. Optics Communications, 2015, 350: 248-251. DOI: 10.1016/j.optcom.2015.04.034
[7] LI S, LIU H, SUN Q, et al. Multi-channel terahertz wavelength division demultiplexer with defects-coupled photonic crystal wave guide[J]. Journal of Modern Optics, 2016, 63(10): 955-960. DOI: 10.1080/09500340.2015.1111457
[8] Telecommunication Standardization Sector of ITU. Series G: Transmission Systems and Media, Digital Systems and Networks: ITU-T G695[S]. 2018.
[9] 丁涛, 刘占芳, 宋恺. 三维光子晶体的制备[J]. 化学进展, 2008, 20(9): 1283-1293. https://www.cnki.com.cn/Article/CJFDTOTAL-HXJZ200809005.htm DING Tao, LIU Zhanfang, SONG Kai. Fabrication of three- dimensional photonic crystals[J]. Chemical Progress, 2008, 20(9): 1283-1293. https://www.cnki.com.cn/Article/CJFDTOTAL-HXJZ200809005.htm
[10] ZHANG L, ZHANG S, LIU Y, et al. Broadband tunable frequency selective surface for steerable antenna applications[J]. IEEE Transactions on Antennas & Propagation, 2016, 64(12): 5496-5500. http://ieeexplore.ieee.org/document/7676371/
[11] Osipov A V. Minimum reflection properties of planar impedance- matched boundaries[J]. IEEE Transactions on Antennas & Propagation, 2014, 62(11): 5666-5670. http://ieeexplore.ieee.org/document/6891193/
[12] Lim J H, LIU T, Kim S S. Numerical analysis of complex impedance and microwave absorption of metamaterials composed of split cut wires on grounded dielectric substrate[J]. Applied Physics A, 2014, 117(3): 1401-1407. DOI: 10.1007/s00339-014-8562-5
[13] XIONG X, SUN W H, BAO Y J, et al. Switching the electric and magnetic responses in a metamaterial[J]. Physical Review B, 2009, 80(20): 2665-2668. http://adsabs.harvard.edu/abs/2010APS..MARB14013W
[14] Stojanović D B, Radovanović J, Milanović V. Influence of the geometry of terahertz chiral metamaterial on transmission group delays[J]. Optical and Quantum Electronics, 2016, 48(4): 272. DOI: 10.1007/s11082-016-0533-y
[15] 白胜元, 顾培夫, 刘旭, 等. 波分复用薄膜干涉窄带滤光片的设计[J]. 光学仪器, 2001(Z1): 114-119. https://www.cnki.com.cn/Article/CJFDTOTAL-GXYQ2001Z1022.htm BAI Shengyuan, GU Peifu, LIU Xu, at el. Design of WDM thin-film interference narrowband filter[J]. Optical Device, 2001(Z1): 114-119. https://www.cnki.com.cn/Article/CJFDTOTAL-GXYQ2001Z1022.htm
[16] Piran M J, Suh D Y. Learning-driven wireless communications, towards 6G[C]//2019 International Conference on Computing, Electronics & Communications Engineering(ICCECE), 2019: 219-224.
[17] Wakatsuchi H, Paul J, Greedy S, et al. Cut-wire metamaterial design based on simplified equivalent circuit models[J]. IEEE Transactions on Antennas & Propagation, 2012, 60(8): 3670-3678. http://ieeexplore.ieee.org/document/6204328/references?signout=success
[18] MIN L, HUANG L. Perspective on resonances of metamaterials[J]. Optics Express, 2015, 23(15): 19022. DOI: 10.1364/OE.23.019022
下载:
计量
- 文章访问数: 244
- HTML全文浏览量: 104
- PDF下载量: 31
