孫 偉，李曉春，毛軍發(fā)

(上海交通大學(xué) 高速電子系統(tǒng)設(shè)計(jì)與電磁兼容研究教育部重點(diǎn)實(shí)驗(yàn)室，上海 200240)

隨著電路工作頻率的不斷提高和供電電壓的不斷下降，由同步開關(guān)噪聲和電壓降引起的電源完整性問題成為電源配送網(wǎng)絡(luò)設(shè)計(jì)中的重要環(huán)節(jié).電源配送網(wǎng)絡(luò)通常以一對電源地平面的形式出現(xiàn)，它為芯片供電，并為信號線提供電流返回路徑.在工作頻率很高和電源電壓很小的電路中，電源完整性分析與設(shè)計(jì)尤為重要.常用的減小電源地平面噪聲的方法有2種：一種是在板間添加去耦電容，通過降低電源地平面間的輸入阻抗提高供電電壓的穩(wěn)定性；另一種是通過將電源地平面切分來隔離噪聲，這就導(dǎo)致了電源地平面上會出現(xiàn)一些狹小的細(xì)縫.因此，準(zhǔn)確分析具有帶隙和去耦電容結(jié)構(gòu)的電源地平面成為電源完整性設(shè)計(jì)的重要環(huán)節(jié).

高速電路系統(tǒng)中的電源分配網(wǎng)絡(luò)分析方法可以分為2大類：一類是電磁場分析方法；另一類是等效電路分析方法.電磁場分析方法精度高，但效率低.文獻(xiàn)[1]中使用時(shí)域有限差分算法(FDTD)，分析電源地平面特性，該方法是一種全波方法但是算法的時(shí)間步長受限制于穩(wěn)定性約束條件.文獻(xiàn)[2]中使用非連續(xù)伽遼金時(shí)域法提取電源地平面對的散射參數(shù)，該方法同樣是一種顯式時(shí)間步的方法，并且在求解每個(gè)單元陣所需要的消耗比較大.文獻(xiàn)[3]中使用鏡像方法提取電源地平面的阻抗參數(shù)，該方法是一種解析方法.它將地平面的四周等效成理想磁壁邊界(PMC),將待求解端口鏡像到區(qū)域外，但是這種方法只適用于規(guī)則電源地平面.

等效電路方法通過提取電源地平面的等效電路參數(shù)，包括電阻(R)、電感(L)、電導(dǎo)(G)和電容(C)參數(shù)，來分析電源地平面電路的噪聲特性，具有效率高的優(yōu)點(diǎn).針對電源地平面的等效電路模型，文獻(xiàn)[4-5]中引入基于加權(quán)拉蓋爾多項(xiàng)式-時(shí)域有限差分(WLP-FDTD)方法分析電源地平面的電源完整性.WLP-FDTD算法不受FDTD穩(wěn)定性條件約束，是一種無條件穩(wěn)定的時(shí)域方法[6-7].因此，相比FDTD方法，WLP-FDTD具有效率高的優(yōu)點(diǎn).近年來，WLP-FDTD方法也被應(yīng)用于均勻傳輸線[8]、非均勻傳輸線[9].文獻(xiàn)[4]中分析了一個(gè)理想的完整電源地平面.文獻(xiàn)[5]中考慮了去耦電容的影響，但未考慮電源地平面的帶隙結(jié)構(gòu).文獻(xiàn)[10]中提出一種基于諧振腔和傳輸線理論的SPICE模型來分析電源地平面的阻抗特性.在此基礎(chǔ)上，文獻(xiàn)[11]中提出了一種場路結(jié)合的方法分析帶隙電源地平面的散射參數(shù)特性.這2種方法在求解電源地平面中均采用了電磁場數(shù)值方法，但場的求解過程中均沒有考慮金屬平面寄生電阻的影響.

本文使用無條件穩(wěn)定的WLP-FDTD方法分析帶隙電源地平面的噪聲傳播問題.WLP-FDTD方法以加權(quán)拉蓋爾多項(xiàng)式為基函數(shù)，通過伽遼金法消除時(shí)間變量，從而使算法不再受到穩(wěn)定性條件約束.與傳統(tǒng)FDTD算法相比，該方法可以選擇任意的時(shí)間步長，解決了由于計(jì)算單元尺寸變小帶來的計(jì)算效率低問題.相比于文獻(xiàn)[6-7]中采用WLP-FDTD方法分析規(guī)則的電源地平面結(jié)構(gòu)，本文考慮了電源地平面的帶隙結(jié)構(gòu)和去耦電容，在電源地平面的電路模型中引入帶隙的等效耦合電路模型，并且將去耦電容應(yīng)用到等效電路中，拓展了WLP-FDTD在電源地平面的應(yīng)用，實(shí)現(xiàn)了基于WLP-FDTD方法的帶隙和去耦電容電源地平面的噪聲傳播分析.計(jì)算結(jié)果表明，本文提出的基于WLP-FDTD的帶隙電源地平面噪聲分析方法與FDTD方法具有相同精度，但計(jì)算效率有顯著提高.

1 帶隙和去耦電容電源地平面的建模

圖1所示為一典型具有帶隙和去耦電容的電源地平面結(jié)構(gòu)，上下表面由大片金屬分別組成電源平面和地平面，中間為填充介質(zhì)，去耦電容Cm連接在兩塊金屬板之間，帶隙位于電源平面上.圖中：w和l分別為電源地平面的寬度和長度；h為金屬層的厚度；d為介質(zhì)層的厚度；s為縫隙的寬度；εr為填充介質(zhì)的相對介電常數(shù)；σc為金屬的電導(dǎo)率，P為輸出點(diǎn).帶隙電源地平面結(jié)構(gòu)可以看作由兩部分組成，即平面部分和帶隙部分，分別對這兩部分進(jìn)行建模.

圖1 具有帶隙和去耦電容電源地平面三維結(jié)構(gòu)Fig.1 The 3-D geometry of power/ground planes with slots and decoupling capacitors

將無隙電源地平面等效為二維RLGC傳輸線模型，即分成M×N個(gè)單元，板內(nèi)介質(zhì)損耗可以忽略，即G=0，因此每個(gè)單元由RLC電路組成.圖2所示為電源地平面單元模型.圖中：Inoise(i,j,t)為外部噪聲源(i=1,2,…,M；j=1,2,…,N).單元中x和y方向的電容、電感、電阻分別為[7]

C(i,j)=ε0εrmn/d

(1)

Lx(i,j)=μ0dm/n

(2)

Ly(i,j)=μ0dn/m

(3)

(4)

(5)

式中：m和n分別為單元(i,j)在x和y方向的長度；ε0和μ0分別為自由空間的介電常數(shù)和磁導(dǎo)率.

圖2 電源地平面單元模型Fig.2 The circuit model of power/ground planes

根據(jù)基爾霍夫電壓電流定律，點(diǎn)(i,j)處沿x和y方向的電流和該點(diǎn)的電壓之間的關(guān)系[5]為

U(i,j,t)-U(i+1,j,t)=

(6)

U(i,j,t)-U(i,j+1,t)=

(7)

Ix(i,j,t)-Ix(i+1,j,t)+Iy(i,j,t)-

Iy(i,j+1,t)+Inoise(i,j,t)=

(8)

Inoise(i,j,t)表示該位置的電流噪聲.

由于d通常很小，遠(yuǎn)小于電路中工作的波長，平面邊緣可被視為開路，邊界被等效為理想磁壁邊界.因此，電源地平面的邊界條件可以表示為[5]

Ix(1,j)=0

(9)

Ix(M+1,j)=0

(10)

Ix(i,1)=0

(11)

Iy(i,N+1)=0

(12)

電源地平面工作在高頻段，帶隙之間的耦合效應(yīng)不能忽略，因此帶隙部分也需要精確建模.帶隙周圍會激發(fā)電場和磁場，電場和磁場的耦合效應(yīng)分別可以用電容和電感表示[10-11].因此，可將帶隙部分等效成耦合電路，如圖3所示.耦合電路包含2K個(gè)節(jié)點(diǎn)，K的大小與帶隙長度有關(guān)，每個(gè)單元的自電容Cs、耦合電容Cc、自電感Ls、耦合電感Lc[10-11]為

(13)

(14)

(15)

(16)

式中：Ce和Co分別為偶模等效電容和奇模等效電容；上標(biāo)a和d分別表示空氣和介質(zhì)下的電容，文獻(xiàn)[11]中給出了2種模式下電容的計(jì)算公式.

圖3 帶隙部分的等效電路Fig.3 The equivalent circuit of slots

通過帶隙結(jié)構(gòu)的耦合等效電路，可以將帶隙和平面級聯(lián).假設(shè)節(jié)點(diǎn)(i,j)位于帶隙等效電路的某個(gè)節(jié)點(diǎn)處，令相應(yīng)的鏡像節(jié)點(diǎn)位于(i+1,j)處，這樣就建立起兩個(gè)相鄰平板間的聯(lián)系，帶隙部分的電壓和電流關(guān)系為

(17)

U(i,j+1,t)=

(18)

Ix(i,j,t)-Ix(i+1,j,t)+Iy(i,j,t)-

(19)

基于等效電路方法可以很方便地處理電源地平面中具有去耦電容的情況.網(wǎng)格剖分的時(shí)候，選擇將節(jié)點(diǎn)應(yīng)用到去耦電容處.對于存在耦合電容的節(jié)點(diǎn)處，式(8)可以替換為

Ix(i,j,t)-Ix(i+1,j,t)+Iy(i,j,t)-

Iy(i,j+1,t)+Inoise(i,j,t)=

(20)

2 應(yīng)用WLP-FDTD求解

不同于傳統(tǒng)的FDTD算法，本文采用WLP-FDTD方法分析帶隙電源地平面噪聲的瞬態(tài)響應(yīng).WLP-FDTD方法引入拉蓋爾多項(xiàng)式，并利用多項(xiàng)式加權(quán)正交性構(gòu)造一組正交多項(xiàng)式組成基函數(shù)[6]

φp(st)=e-st/2Lp(st)

(21)

式中：s為時(shí)間尺度因子.基函數(shù)具有正交性：

(22)

利用這組基函數(shù)，電壓電流可以展開為

(23)

(24)

(25)

以式(23)為例，對時(shí)間的一階微分可以表示為

(26)

式中：W(t)為平面任意位置的電壓或電流.

將式(26)代入式(6)、(7)和(17)～(20)，可得

(27)

(28)

(29)

(30)

(31)

(32)

在式(27)～(32)的兩邊乘以φq(st)并在st=[0,∞)積分，可得

(33)

(34)

0.5s(C(i,j)+Cc)Uq(i,j)+

(35)

(36)

Uq(i,j+1)=Uq(i,j)+(Ls-Lc)×

(37)

(38)

(39)

無窮積分上界可以用噪聲持續(xù)時(shí)間Td代替.邊界條件式(9)～(12)可以表示為

(40)

(41)

(42)

(43)

電源地平面中的電壓電流關(guān)系為式(33)～(35)，帶隙部分的電壓電流關(guān)系為式(36)～(38)，邊界條件為式(40)～(43).根據(jù)式(33)～(38)，將q-1階項(xiàng)和噪聲項(xiàng)移到等式右邊，整理可以得到電源地平面的電壓電流矩陣形式:

AWq=Inoise+bq-1

(44)

(45)

可以看出，式(44)右邊的階數(shù)低于左邊，因此從q=0 開始求解，采用遞歸方法求解式(44),可以計(jì)算出電源地平面上的電壓電流值.此外，A與階數(shù)q無關(guān)，因此整個(gè)迭代過程中只需要在計(jì)算開始時(shí)對矩陣A求逆一次.

式(33)～(38)的求和上界q通常用有限值Q來代替.Q值的大小對計(jì)算精度和計(jì)算時(shí)間都有影響，雖然Q值越大，計(jì)算精度越高，但消耗的計(jì)算時(shí)間也越多，因此需要選擇合適的Q值.假設(shè)電路中工作的信號帶寬有限，截止頻率為fB，則有[5]

Qmin=2fBTf+1

(46)

式中：Tf為信號的時(shí)間周期.

3 數(shù)值結(jié)果

采用2個(gè)帶隙電源地平面結(jié)構(gòu)作為算例，驗(yàn)證所提出的基于WLP-FDTD算法的計(jì)算精度與效率.

算例1不帶去耦電容，平面尺寸為50 mm×50 mm,帶隙位于電源平面的中心，帶隙尺寸為 0.5 mm×20 mm,中間介質(zhì)厚度為 0.05 mm,相對介電常數(shù)為 4.2，金屬層厚度 0.018 mm.平面左側(cè)存在噪聲電流，其形式為正弦調(diào)制的高斯脈沖

Inoise(t)=

(47)

式中：Td=fc/2；Tc=3Td；fc=1 GHz.高斯脈沖的周期Tf=11.71 ns，截止頻率fB=5 GHz,代入式(46)得Qmin，本文中Qmin值選取為120，K值選取為20.通過考察如圖1所示右側(cè)邊緣中心P點(diǎn)的電壓來分析噪聲對平面性能的影響.當(dāng)采用FDTD方法分析該問題時(shí)，由于FDTD受穩(wěn)定性條件約束[6]，時(shí)間空間步長必須滿足

(48)

式中：vmax為介質(zhì)中最大相速；Δx和Δy分別為x和y方向的網(wǎng)格尺寸.當(dāng)網(wǎng)格尺寸選擇1 mm，F(xiàn)DTD的時(shí)間步長必須小于 2.357 ps.圖4給出了SPICE仿真軟件結(jié)果、FDTD和WLP-FDTD算法的計(jì)算結(jié)果.可以看出，三者幾乎完全吻合，驗(yàn)證了WLP-FDTD方法的準(zhǔn)確性.針對算例1，表1給出了FDTD和WLP-FDTD的計(jì)算時(shí)間對比，仿真高斯脈沖的周期數(shù)為20個(gè).由于WLP-FDTD沒有時(shí)間步長影響，因此計(jì)算效率遠(yuǎn)高于FDTD.

圖4 未加去耦電容的電源地平面中P點(diǎn)電壓Fig.4 The output voltage of the power ground planes without coupling capacitor at the point P

表1 算例1計(jì)算時(shí)間對比Tab.1 Comparison of computation time of the 1st example

算例2平面尺寸與算例1相同，增加了2對去耦電容.去耦電容之間的間距為2 mm,分別位于帶隙兩側(cè)和輸出點(diǎn)P的兩側(cè).去耦電容為10 pF.由于每對去耦電容的間距較小，網(wǎng)格劃分比較密，網(wǎng)格尺寸最大為1 mm×1 mm，此時(shí)FDTD的時(shí)間步長必須小于 2.357 ps.圖5給出了SPICE軟件、FDTD和WLP-FDTD對算例2的仿真計(jì)算結(jié)果.可見，三者具有極高的吻合度.表2給出了相同空間步長情況下FDTD和WLP-FDTD對算例2的計(jì)算時(shí)間，仿真高斯脈沖的周期數(shù)為20個(gè).可見，WLP-FDTD算法的計(jì)算效率遠(yuǎn)高于FDTD算法，且其優(yōu)勢隨網(wǎng)格劃分密度的增加而更加顯著.

圖5 添加去耦電容的電源地平面中P點(diǎn)電壓Fig.5 The output voltage of the power ground planes with coupling capacitors at the point P

表2 算例2計(jì)算時(shí)間對比Tab.2 Comparison of computation time of the 2nd example

4 結(jié)語

本文基于一種無條件穩(wěn)定的WLP-FDTD時(shí)域有限差分算法，分析帶隙電源地平面中噪聲的傳播問題.首先，基于電源地平面的等效電路模型和帶隙的等效電路模型建立電路方程；然后采用WLP-FDTD方法求解方程，從而得到復(fù)雜電源地平面結(jié)構(gòu)噪聲傳播的瞬態(tài)響應(yīng).理論與實(shí)驗(yàn)結(jié)果證明，與傳統(tǒng)FDTD方法相比，WLP-FDTD方法能夠?qū)崿F(xiàn)無條件穩(wěn)定，因此在分析復(fù)雜結(jié)構(gòu)電源地平面噪聲時(shí)可以顯著提高計(jì)算效率.