冀維林 劉仲亞 李 莉

(1.國(guó)家無(wú)線電監(jiān)測(cè)中心，北京 100037；2.北京郵電大學(xué)電子工程學(xué)院，北京 100876)

引 言

在高速大規(guī)模集成電路和多芯片組件中，多導(dǎo)體傳輸線常被用作信號(hào)連接線.通常，在直流或低頻情況下，信號(hào)連接線可以看作是簡(jiǎn)單的金屬導(dǎo)體，僅僅起著電連通的作用.但是，隨著半導(dǎo)體材料科學(xué)和電子信息技術(shù)的發(fā)展，信號(hào)脈沖的上升時(shí)間和寬度已達(dá)到皮秒量級(jí)，對(duì)應(yīng)的頻譜已進(jìn)入微波、毫米波波段，高頻分量的波長(zhǎng)和傳輸線的尺寸已處于同一量級(jí)；由于色散、終端失配及線間的分布耦合，還將引起信號(hào)的畸變和線間耦合，影響到高速電路系統(tǒng)的總體性能指標(biāo).因此，必須對(duì)這些現(xiàn)象加以分析和研究，而多導(dǎo)體傳輸線時(shí)域響應(yīng)分析便可較為準(zhǔn)確地分析信號(hào)連接線上各點(diǎn)電壓和電流[1-4].同時(shí)，在實(shí)際電子設(shè)備和系統(tǒng)中，由于各種傳輸線系統(tǒng)結(jié)構(gòu)的復(fù)雜性，通信線與電力線存在許多非均勻的，或者是端接非線性負(fù)載結(jié)構(gòu)的傳輸線.因此，必須對(duì)廣義傳輸線進(jìn)行分析.

隨著計(jì)算機(jī)技術(shù)的發(fā)展，數(shù)值方法在研究和分析電磁場(chǎng)問(wèn)題中發(fā)揮了越來(lái)越重要的作用.例如頻域內(nèi)的有限元法、矩量法和單矩法等；時(shí)域內(nèi)的時(shí)域有限差分法、傳輸線矩陣法和時(shí)域積分法等.其中，時(shí)域有限差分法(FDTD)是由有限差分法發(fā)展而來(lái)的數(shù)值分析方法，直接由麥克斯韋方程對(duì)電磁場(chǎng)進(jìn)行計(jì)算機(jī)模擬.近年來(lái)，時(shí)域有限差分法在生物電磁劑量學(xué)、電磁熱療系統(tǒng)的計(jì)算機(jī)模擬、天線輻射特性的計(jì)算問(wèn)題以及設(shè)備和系統(tǒng)間的電磁兼容特性分析和預(yù)測(cè)等方面得到廣泛的應(yīng)用和迅速的發(fā)展[5-6].論文主要應(yīng)用FDTD法對(duì)不均勻多導(dǎo)體傳輸線的瞬態(tài)響應(yīng)、波形特點(diǎn)以及終端失配情況進(jìn)行仿真分析.文中給出了FDTD算法處理該問(wèn)題的迭代公式，通過(guò)建模、仿真和分析，驗(yàn)證了FDTD算法處理不均勻多導(dǎo)體傳輸線瞬態(tài)響應(yīng)問(wèn)題的有效性和正確性.

1 不均多導(dǎo)體傳輸線方程及差分形式

在不考慮外界電磁場(chǎng)作用的情況下，多導(dǎo)體傳輸線上的電流和電壓滿足如下電報(bào)方程[7-9]

(1)

式中：V和I分別為多導(dǎo)體傳輸線上的電壓和電流矩陣；R和G分別為多導(dǎo)體傳輸線上的電阻和電導(dǎo)矩陣；L和C分別為多導(dǎo)體傳輸線上的電感和電容矩陣.對(duì)于無(wú)耗多導(dǎo)體傳輸線(R=0，G=0)，應(yīng)用FDTD法對(duì)多導(dǎo)體傳輸線方程進(jìn)行差分，并結(jié)合傳輸線始端和終端邊界條件可得到傳輸線上各點(diǎn)的迭代公式：

(2)

(3)

(4)

(5)

式中：L和C根據(jù)不均勻傳輸線特性分段為L(zhǎng)k和Ck； Δz、Δt分別為空間和時(shí)間步長(zhǎng)；k、n分別表示空間和時(shí)間的劃分?jǐn)?shù).通常，為了保證算法的穩(wěn)定，必須滿足Δt≤Δz/v.當(dāng)取等號(hào)時(shí)，Δt定義為最佳時(shí)間步長(zhǎng)，其中v為電磁波傳播的最大模式速度，在多導(dǎo)體傳輸線情況下，會(huì)有多個(gè)模式速度存在，可通過(guò)計(jì)算LC的特征值求得.

2 模型的建立與仿真分析

圖1 不均勻多導(dǎo)體傳輸線模型

2.1 情況一：線性負(fù)載

如圖1所示的不均勻多導(dǎo)體傳輸線模型(銅導(dǎo)帶寬度為50×(1+k(z)) μm，厚度為5 μm).其中，多導(dǎo)體傳輸線長(zhǎng)l=0.5 m，始端和終端負(fù)載分別為RS1=50 Ω、RS2=50 Ω、RL1=50 Ω、RL2=50 Ω.信號(hào)源為上升/下降沿為0.5 ns，脈沖寬度為5 ns的方波.非均勻傳輸線的分布參數(shù)為：

式中：l(z)=387/(1+K(z))；lm(z)=K(z)L(z)；c(z)=104.13/(1-K(z))；cm(z)=-K(z)C(z)；K(z)=0.25(1+0.6sin(πz+π/4)).

應(yīng)用FDTD法處理電磁學(xué)問(wèn)題時(shí)，由于在計(jì)算機(jī)的計(jì)算空間中電磁波只能跳躍式地傳播，不可能連續(xù)傳播，因此，對(duì)于不同波長(zhǎng)的電磁波，在特定的網(wǎng)格劃分情況下，計(jì)算機(jī)中電磁波傳播比起實(shí)際物理空間中的電磁波會(huì)出現(xiàn)誤差，包括速度和方向.針對(duì)數(shù)值色散問(wèn)題，可以通過(guò)減小空間和時(shí)間步長(zhǎng)，減小數(shù)值色散的影響.一般當(dāng)空間步長(zhǎng)滿足如下條件時(shí)，差分近似帶來(lái)的數(shù)值色散將非常小[5]

Δz≤λmin/12，

(6)

(a)

(b)圖2 不均勻多導(dǎo)體傳輸線各端口瞬態(tài)電壓響應(yīng)(光速)

(a)

(b)圖3 不均勻多導(dǎo)體傳輸線端口瞬態(tài)電壓響應(yīng)(最大模式速度)

圖4 不均勻多導(dǎo)體傳輸線特性阻抗曲線

選取兩倍的脈沖信號(hào)轉(zhuǎn)折頻率(1 GHz)作為信號(hào)最高頻率求得最小波長(zhǎng)，代入式(6)，近似得到空間步長(zhǎng)為0.01 m，可得到仿真網(wǎng)格數(shù)為50.通常，仿真選取的電磁波傳播速度只要大于電磁波在傳輸線上的最大模式速度，仿真就可以順利進(jìn)行.首先比較當(dāng)仿真時(shí)選取的電磁波的傳播速度不同時(shí)對(duì)仿真結(jié)果的影響.電磁波的傳播速度分別選取光速c和最大模式速度cm.運(yùn)用MATLAB軟件編程，傳輸線兩端瞬態(tài)電壓響應(yīng)如圖2和3所示.從仿真結(jié)果可以看出，在圖2中，信號(hào)和串?dāng)_波形都具有比較明顯的抖動(dòng)紋波，而圖3則很平滑.可見(jiàn)，采用最大模式速度作為求解相關(guān)參數(shù)(Δt)并進(jìn)行仿真，可以得到最優(yōu)的結(jié)果.同時(shí)，信號(hào)經(jīng)過(guò)3.285 ns后到達(dá)傳輸線的終端，這說(shuō)明不均勻傳輸線對(duì)信號(hào)的傳播速度沒(méi)有影響.然而，與均勻多導(dǎo)體傳輸線各端口的電壓響應(yīng)不同，由于傳輸線上各個(gè)點(diǎn)的特征阻抗存在一定的差異，所以信號(hào)在傳播的過(guò)程中會(huì)有反射現(xiàn)象，需要經(jīng)過(guò)較長(zhǎng)的時(shí)間才能趨于穩(wěn)定.為了更好地說(shuō)明不均勻多導(dǎo)體傳輸線各個(gè)端口電壓響應(yīng)波形，我們畫(huà)出了傳輸線上各點(diǎn)的特征阻抗值和(Cm/C)+(Lm/Lw)值的曲線，如圖4和5所示.從圖中看出，傳輸線特性阻抗值在傳輸線長(zhǎng)度方向上先減小后變大，變化速率也是先變小后變大.信號(hào)經(jīng)過(guò)傳輸線時(shí)，由于傳輸線匹配特性遵循由差變好再變差的規(guī)律，所以信號(hào)線始端信號(hào)波形呈現(xiàn)由大變小再變大的特點(diǎn)，可見(jiàn)信號(hào)線始端信號(hào)與特性阻抗的分布有關(guān)；另外，受擾線近端串?dāng)_與(Cm/C)+(Lm/Lw)的值成正比，所以近端串?dāng)_具有先增大后減小的變化趨勢(shì).

圖5 (Cm/C)+(Lm/Lw)曲線

圖6 端接非線性負(fù)載不均勻多導(dǎo)體傳輸線模型

2.2 情況二：非線性負(fù)載

端接非線性負(fù)載的不均勻多導(dǎo)體傳輸線的模型如圖6所示(作為參考地面的導(dǎo)體未畫(huà)出).在公式(2)、(3)、(4)、(5)的基礎(chǔ)上，只需將N+1處的電壓迭代公式修改為式(7)，便可得到端接非線性負(fù)載的多導(dǎo)體傳輸線差分公式為.

(a)

(b)圖7 端接非線性負(fù)載多導(dǎo)體傳輸線各端口瞬態(tài)電壓響應(yīng)

(a)

(b)圖8 各端口瞬態(tài)電壓響應(yīng)(RS2=RL2=0 Ω)

(7)

非線性負(fù)載的伏安特性關(guān)系為：

式中:IA=10 nA；VT=25 mV.

多導(dǎo)體傳輸線始端和終端負(fù)載分別為RS1=50 Ω、RS2=50 Ω、RL1=50 Ω、RL2=50 Ω，信號(hào)源為上升/下降沿為0.5 ns，脈沖寬度為5 ns的方波.運(yùn)用MATLAB軟件編程，仿真結(jié)果如圖7.可以看出，非線性負(fù)載的接入，使得多導(dǎo)體傳輸線終端匹配條件變差，信號(hào)的反射增大，信號(hào)需要經(jīng)過(guò)很長(zhǎng)的時(shí)間才能趨于穩(wěn)定.同時(shí)，信號(hào)線輸出端信號(hào)較之未接非線性負(fù)載時(shí)增大0.16 V.

當(dāng)受擾線兩端阻抗RS2=RL2=0 Ω的情況，仿真結(jié)果如圖8所示.未接非線性負(fù)載時(shí)，受擾線可以看作電位為零的參考地面，受擾線上是不會(huì)有串?dāng)_信號(hào)的.接入非線性負(fù)載后，受擾線始端全反射，信號(hào)為零，信號(hào)線的始端電壓可以逐漸趨于零.終端處，無(wú)論是信號(hào)線還是受擾線，信號(hào)呈現(xiàn)振蕩的狀態(tài)， 脈沖每經(jīng)過(guò)兩個(gè)傳輸線傳播時(shí)間反相一次.因此，為了防止出現(xiàn)信號(hào)在傳輸線上振蕩的現(xiàn)象，在實(shí)際中應(yīng)盡量避免受擾線僅有非線性負(fù)載連接的情況.

3 結(jié) 論

隨著傳輸線系統(tǒng)結(jié)構(gòu)的不斷復(fù)雜，信號(hào)頻率的不斷提高，終端失配、色散及線間的分布耦合引起的信號(hào)畸變和線間耦合問(wèn)題越來(lái)越嚴(yán)重.針對(duì)此問(wèn)題，論文提出了通過(guò)建立不均勻多導(dǎo)體傳輸線FDTD法差分模型，運(yùn)用MATLAB軟件對(duì)不均勻多導(dǎo)體傳輸線瞬態(tài)響應(yīng)進(jìn)行仿真的算法.在此基礎(chǔ)上進(jìn)一步從理論上對(duì)非均勻多導(dǎo)體傳輸線各端口瞬態(tài)響應(yīng)的波形特點(diǎn)進(jìn)行了研究.另外，對(duì)端接非線性負(fù)載的不均勻多導(dǎo)體傳輸線端口出現(xiàn)失配的情況進(jìn)行了仿真，為實(shí)際工程設(shè)計(jì)和電磁兼容預(yù)測(cè)提供了理論參考.仿真結(jié)果驗(yàn)證了算法的正確性，在今后的工作中，可以結(jié)合工程實(shí)際，不斷拓展FDTD算法處理電磁兼容問(wèn)題的范圍.

[1] 姜仁波.多導(dǎo)體傳輸線電磁兼容的分析與研究[D].北京:華北電力大學(xué)，2004． JIANG Renbo.The Analysis and Research of EMC on Multiconductor Transmission Lines[D].Beijing:North China Electric Power University，2004.(in Chinese)

[2] 高本慶.時(shí)域有限差分法[M].北京:國(guó)防工業(yè)出版社，1995．

[3] 毛軍發(fā).高速集成電路中多導(dǎo)體傳輸線時(shí)域響應(yīng)分析[D].上海:上海交通大學(xué)，1992．

MAO Junfa.Analysis of the Time-Domain Response of Multiconductor Transmission Lines in High-Speed Integrated Circuits[D].Shanghai:Shanghai Jiao Tong University，1992.(in Chinese)

[4] 李 莉，冀維林.端接非線性負(fù)載的不等長(zhǎng)傳輸線瞬態(tài)分析[J].電波科學(xué)學(xué)報(bào)，2009，24(3):529-532．

LI Li，JI Weilin.Transient analysis of unequal length multiconductor transmission lines loaded by nonlinear devices[J].Chinese Journal of Radio Science，2009，24(4):529-532.(in Chinese)

[5] 葛德彪，閆玉波.電磁波時(shí)域有限差分方法[M].西安:西安電子科技大學(xué)出版社，2002．

[6] 齊 磊，盧鐵兵，崔 翔.端接非線性負(fù)載的非均勻傳輸線瞬態(tài)分析[J].電波科學(xué)學(xué)報(bào)，2003，18(2):153-156．

QI Lei，LU Tiebing，CUI Xiang，Transient analysis for nonuniform transmission lines with nonlinear loads[J].Chinese Journal of Radio Scienc，2003，18(2):153-156.(in Chinese)

[7] PAUL C R.Analysis of Multiconductor Transmission Lines[M].New York:John Wiley & Sons Press，1994．

[8] ORIANDI A，PAUL C R.FDTD analysis of lossy，multiconductor transmission lines terminated in arbitrary loads[J].IEEE Trans on Electromagnetic Compatibility，1996，38(3):388-399．

[9] TAEB A，ABDIPOUR A，MOHHAMADI A.FDTD analysis of the lossy coupled transmission lines loaded by nonlinear devices[C]// Asia-Pacific Microwave Conference Proceedings，December 4-7，2005，5．