陳華等

摘 要： 為了尋找一種設計出特異媒質單元結構的可行方法，利用分形理論中的文法構圖法，采用不同的公理和改寫規(guī)則，設計出兩種全新的分形圖案，之后將其周期性地分別蝕刻在介質板兩面，通過電磁仿真軟件HFSS得到其散射參數(shù)的幅值和相位，利用等效參數(shù)提取軟件提取其等效參數(shù)，成功構造出特異媒質中的一類，即左手材料，且在5.5～6 GHz頻率下其等效介電常數(shù)和等效磁導率均為負。利用分形理論構造特異媒質單元結構的方法是方便可行的，且圖形有無窮可能性。

關鍵詞： 分形； 左手材料； 等效介電常數(shù)； 等效磁導率

中圖分類號： TN82?34 文獻標識碼： A 文章編號： 1004?373X（2015）08?0096?03

Left?handed metamaterials design based on fractal geometry

CHEN Hua1， JIA Jia2， JIA Yuan?yuan1， CHEN Zhe?si3

（1. Center of Disease Control and Prevention， Chengdu Military Region， Chengdu 610021， China；

2. The information Center of Comprehensive Plan Office of Air Force Equipment Department， Beijing 100843， China；

3. Department of Rehabilitation Medicine， General Hospital of Chengdu Military Region， Chengdu 610083， China）

Abstract： To find a feasible way to design metamaterial structure， by using different axioms and rewriting rules， two kinds of new fractal pattern were designed， and then periodically etched in the medium plate on both sides respectively. The amplitude and phase of scattering parameter were obtained with electromagnetic simulation software HFSS. The equivalent parameter is extracted with the equivalent parameter extraction software. The left?handed metamaterial in peculiar medium was successfully constructed， whose equivalent dielectric constant and equivalent permeability are both negative under 5.5～6 GHz frequency. The method to construct metamaterial structure by using fractal theory is convenient and feasible， and the graphics have infinite possibility.

Keywords： fractal； left?handed metamaterial； effective permittivity； effective permeability

0 引 言

“Metamaterials”是本世紀物理學領域出現(xiàn)的一個新的學術詞匯，國內稱為特異媒質，或是電磁超材料，其是指一些具有天然材料所不具備的超常物理性質的人工復合結構或復合材料[1]。其典型代表有左手材料（LHMs）、零折射率材料（ZIMs）、光子晶體（EBG）等[2?4]。左手材料是指其等效介電常數(shù)和等效磁導率同時為負值的特異媒質，由于其具有諸如如負折射現(xiàn)象[5]、逆多普勒效應、逆Cerenkov輻射[6]等不同尋常的物理性質，因此在微波平板聚焦透鏡、帶通濾波器、耦合器、天線以及隱身衣等方面具有廣泛的應用前景。未來，左手材料將會在無線通信的發(fā)展中起到不可忽略的作用。

分形（Fractal），源于拉丁詞語fractus，意思是不規(guī)則支離破碎形。其概念由法國數(shù)學家B.Mandelbrot提出，用以描述人們在研究自然和科學實驗中遇到的無法用傳統(tǒng)歐式幾何描述，卻具有某種自相似特性（局部形態(tài)與整體形態(tài)相似）的不規(guī)則結構和現(xiàn)象[7]。利用分形原理，可以設計出大量形狀各異的分形圖案。

目前，大部分特異媒質的設計為由平面單元組成的周期性結構，但在工程運用中，如何設計特異媒質卻一直是一個難點，本文正是探索利用分形圖案設計特異媒質中的左手材料。

1 兩種新的分形圖案設計

分形理論是當今十分風靡和活躍的新理論。其最基本的特點是用分形分維的數(shù)學工具來描述研究客觀事物。構造分形的方法很多：迭代函數(shù)法、文法構圖法、逃逸時間算法等。本文采用文法構圖法來構造分形，字符串替換是此法的核心思想[8]。可以理解為：若一個字符串由字母a和b組成（a和b可能出現(xiàn)多次），每一個字母對應一個改寫規(guī)則，如a→b，b→ab。假設初始字符串（稱為公理）僅由字母aab組成，那么經過兩次改寫后，字符串變成了aab→bbab→ababbab，以此類推，這種做法用來表示生長過程是很有意義的。之后，再將每一種字符賦予不同的操作，就能生成各種圖形。例如生成三階Koch曲線時，只要設置公理為F，改寫規(guī)則為F→F+F--F+F，迭代次數(shù)為3即可，之后規(guī)定字符代表的操作為：遇到“F”，向前移動一步并畫線，步長為d；遇到“+”，逆時針方向旋轉60°；遇到“-”，順時針方向旋轉60°。這樣三階Koch曲線便可以生成。如圖1所示為Koch曲線的生成過程。

圖1 Koch曲線生成過程

利用此方法，本文設計的第一種分形圖形為“箭頭”形，采用公理為M+++M+++M+++M，改寫規(guī)則為M→M++++A--A--A--A--A++++M，其中字母M代表向前畫一條長度為M的線段，字母A代表向前畫一條長度為A的線段，“+”代表逆時針旋轉30度，“-”代表順時針旋轉30°，經過一次改寫后，生成的圖形如圖2所示。

圖2 “箭頭”形分形圖案

設計的第二種分形圖形為“十”字形，采用公理為H-H-H-H，規(guī)則為H→H-H--H-H-H--H-H。字母H代表向前畫一條長度為F+K，寬度為K 的矩形條，“+”代表逆時針旋轉90°，“-”代表順時針旋轉90°。經過一次改寫后，生成的圖形如圖3所示。

2 等效參數(shù)提取軟件的正確性驗證

等效參數(shù)提取軟件是一款基于Smith方法而自行開發(fā)的軟件，其原理是通過散射參數(shù)的幅度和相位計算出折射率和阻抗，之后再獲得等效介電常數(shù)和等效磁導率值的方法[9]。為了保證后期計算出來數(shù)據(jù)的可靠性，首先對等效電磁參數(shù)提取軟件的正確性進行驗證。

圖3 “十”字形分形圖案

驗證采用的模型是著名的Srr?Rod結構，將開口諧振環(huán)（Srr）和金屬細線（Rod）分別印刷在微波基板兩面，具體結構如圖4（a）所示，根據(jù)文獻[10]中的數(shù)據(jù)，單元結構尺寸分別為：矩形環(huán)外環(huán)邊長w=3 mm，矩形環(huán)條帶寬度c=0.25 mm，內外環(huán)間距d=0.5 mm，矩形環(huán)開口寬度g=0.5 mm，金屬細線寬度g=0.5 mm，金屬細線長度a=5 mm，單元間距4.8 mm，微波基板的介電常數(shù)為3.84，損耗角正切為0.018，長度和寬度均為5 mm，厚度0.25 mm。圖4（b）為仿真模型俯視圖，波沿x方向傳播。

圖4 Srr?Rod單元結構示意圖和Srr?Rod模型俯視圖

經過仿真計算之后得到的散射參數(shù)幅值和相位如圖5所示。用等效參數(shù)提取軟件提取出等效介電常數(shù)和等效磁導率的實部和虛部，之后分別與參考文獻中得出的等效參數(shù)結果對比，如圖6所示，可以發(fā)現(xiàn)兩者結果十分吻合。

3 設計模型的參數(shù)提取

在介電常數(shù)2.65，邊長L=12 mm，厚度1 mm的微波基板兩側分別蝕刻如圖2、圖3所示的“箭頭”形和“十”字形兩種分形圖形。其中，參數(shù)M=4.3 mm，A=2.4 mm，F(xiàn)=1.1 mm，K=1.2 mm，電磁波沿z軸垂直介質板入射。之后利用該軟件對此基于分形設計出來的模型進行參數(shù)提取。

圖5 S參數(shù)幅值和相位

圖6 等效介電常數(shù)實部和虛部

根據(jù)仿真計算的散射參數(shù)幅值和相位結果，得到了此材料的等效介電常數(shù)、等效磁導率以及折射率的實部和虛部具體值。結果表明，在5.5～6 GHz頻率范圍內，此模型的等效介電常數(shù)、等效磁導率及折射率均為負，即呈現(xiàn)左手材料性質，如表1所示。

表1 模型在5.5～6.0 GHz處的等效電磁參數(shù)實部和虛部數(shù)據(jù)

圖7為該左手材料的等效介電常數(shù)、等效磁導率以及折射率的實部和虛部曲線圖。

圖7 左手材料實部和虛部

4 結 語

特異媒質的單元結構設計一直是個難點。但在本文中，利用分形理論中文法構圖法，一種新型左手材料被成功設計出來，表明此方法是有效可行的。更多特異媒質結構，如零折射率材料，EBG結構等都可嘗試利用此法進行設計。

參考文獻

[1] 姜濤.Metamaterial的若干電磁性質及其應用研究[D].杭州：浙江大學，2011.

[2] NADER Engheta， RICHARD Z W. Metamaterials： physics and engineering explorations [M]. [S.l.]： IEEE Press， 2006.

[3] ROTMAN W. Plasma simulation by artifcialdielectrics and parallelolate media [J]. IRE Transactions on Antennas Propagat， 1962， 4： 82?84.

[4] AL? Andrea. Epsilon?near?zero metamaterials and electromagnetic source：Tailoring the radiation phase pattern [J]. Physical Review B， 2007， 6： 102?105.

[5] PARAZZOLI C G， GREEGO R B， LI K， et al. Experimental verification and simulation of negative index of refraction using snell′s law [J]. Physcis Review Letters， 2003， 12： 107?110.

[6] KONG J A. Electromagnetics [M]. New York： Wiley， 1975.

[7] 宋小弟，馮恩信.分形天線工程及其新進展[J].無線通信技術，2008（2）：48?54.

[8] 朱華，姬翠翠.分形理論及其應用[M].北京：科學出版社，2011.

[9] SCHURIG D， MOCK J J， JUSTICE B J， et al. Metamaterial electromagnetic cloak at microwave frequencies [J]. Science， 2006， 314： 977?980.

[10] LUBKOWSKI G， SCHUHMANN R， WEILAND T. Extraction of effective metamaterial parameters by parameter fitting of dispersive models [J]. Microwave and Optical Technology Letters， 2007， 49： 23?26.