丁 橋 付 強(qiáng) 陳 昆 武紅斌

（成都新欣神風(fēng)電子科技有限公司 成都 611731）

1 引言

隨著各軍事應(yīng)用平臺(tái)上電子設(shè)備的集成度越來(lái)越高，對(duì)平臺(tái)進(jìn)行系統(tǒng)級(jí)電磁兼容分析已經(jīng)成為系統(tǒng)電磁兼容設(shè)計(jì)非常重要的一個(gè)環(huán)節(jié)。當(dāng)前，在飛機(jī)、艦船乃至車載平臺(tái)上，安裝相控陣天線以提高平臺(tái)的使用性能，加強(qiáng)平臺(tái)作戰(zhàn)能力，變得更加明顯。平臺(tái)上相控陣天線作為發(fā)射天線具有功率高、增益大、輻射強(qiáng)等優(yōu)點(diǎn)，但由于其副瓣作用，也更容易引起平臺(tái)上或者周圍環(huán)境的電磁兼容問(wèn)題。一般來(lái)說(shuō)，相控陣天線作為一個(gè)復(fù)雜、昂貴的系統(tǒng)，很難有條件在其設(shè)計(jì)或者研制階段進(jìn)行實(shí)際電磁兼容試驗(yàn)測(cè)試，往往是在平臺(tái)集成過(guò)程中才能夠進(jìn)行相關(guān)試驗(yàn)；因此，為了保證相控陣安裝平臺(tái)的電磁兼容性能，預(yù)先對(duì)天線進(jìn)行電磁兼容仿真分析是有益的。相控陣天線由于其口徑面積大、工作頻率高，其安裝平臺(tái)及其附近范圍，往往處于該天線的近場(chǎng)范圍，因此對(duì)相控陣天線進(jìn)行近場(chǎng)輻射場(chǎng)計(jì)算成為相控陣天線電磁兼容分析的重要內(nèi)容。

相控陣在進(jìn)行陣列單元集成時(shí)，由于單元之間的距離較近，各單元之間存在互耦效應(yīng)，引起陣列環(huán)境中的天線單元方向圖與自由空間中的單元方向圖之間存在較大差異，這也為相控陣天線的近場(chǎng)計(jì)算帶來(lái)挑戰(zhàn)；并且受限于當(dāng)前計(jì)算機(jī)性能，對(duì)于大型陣列進(jìn)行全波的電磁仿真計(jì)算是不現(xiàn)實(shí)的。一直以來(lái)，對(duì)大型陣列天線輻射場(chǎng)的計(jì)算都是研究熱點(diǎn)，文獻(xiàn)［1～3］采用小陣外推的方法實(shí)現(xiàn)了大型陣列在考慮互耦條件下的輻射場(chǎng)和散射場(chǎng)計(jì)算，但是都是計(jì)算的陣列遠(yuǎn)場(chǎng)；文獻(xiàn)［4］給出了一種相控陣天線的電磁環(huán)境分析方法，但是該方法中采用了陣列單元在自由空間中的理想方向圖，沒有計(jì)入組陣后單元之間的互耦效應(yīng)；也沒有考慮非均勻激勵(lì)情況；文獻(xiàn)［5］提出了一種基于散射矩陣?yán)碚摰妮椛鋱?chǎng)計(jì)算公式，同樣重點(diǎn)在計(jì)算遠(yuǎn)場(chǎng)方面；文獻(xiàn)［6～7］采用純理論方法討論了近場(chǎng)區(qū)的計(jì)算問(wèn)題，但是針對(duì)口徑面天線；文獻(xiàn)［8］提出了一種計(jì)算陣列天線的近似計(jì)算方法，與本文提供的方法有所區(qū)別。

準(zhǔn)確計(jì)算陣列天線近場(chǎng)，互耦因素是必然要考慮的，文獻(xiàn)［9～11］對(duì)陣列互耦進(jìn)行了研究。本文將在現(xiàn)有研究成果基礎(chǔ)上，提出一種新的考慮互耦的大型陣列天線近場(chǎng)計(jì)算方法，并通過(guò)與全波電磁仿真軟件的計(jì)算結(jié)果進(jìn)行對(duì)比，分析了本文所給出計(jì)算方法容易引起誤差的來(lái)源，給出了本文方法的適用范圍。

2 陣列天線輻射場(chǎng)計(jì)算理論分析

首先，陣列天線按場(chǎng)點(diǎn)離陣列中心的距離不同，可將其輻射場(chǎng)分為近區(qū)場(chǎng)和遠(yuǎn)區(qū)（Fraunhofer區(qū)）場(chǎng)，在不同場(chǎng)區(qū)范圍內(nèi)，陣列天線的輻射場(chǎng)具有不同特征。陣列天線近區(qū)和遠(yuǎn)區(qū)的范圍［4］可近似由下式進(jìn)行確定，當(dāng)場(chǎng)點(diǎn)離陣列中心距離r滿足式（1）時(shí)，則為遠(yuǎn)區(qū)，否則為近區(qū)。

其中D為口徑最大尺寸，λ為工作波長(zhǎng)。

下面，將給出能夠同時(shí)計(jì)算陣列天線近區(qū)和遠(yuǎn)區(qū)輻射場(chǎng)的計(jì)算方法。

考慮一個(gè)按矩形柵格分布的平面陣，并假定所有陣元位于xoy平面內(nèi)，并設(shè)其幾何中心為相位中心，P為場(chǎng)點(diǎn)，O為陣列中心位置，如下圖所示。

圖1 相控陣天線電磁環(huán)境分析示意圖（部分陣列單元）

對(duì)于大型平面陣，一般考慮的陣列近場(chǎng)范圍滿足陣列單元的遠(yuǎn)場(chǎng)條件。因此，平面陣的近場(chǎng)可看作是陣中每個(gè)單元的遠(yuǎn)場(chǎng)疊加。根據(jù)Maxwell方程組，陣列中第i個(gè)單元的遠(yuǎn)場(chǎng)輻射方向圖主體部分可寫成一個(gè)方向性函數(shù)fi（θ，φ）與一個(gè)球面波因子乘積，然后輔之以激勵(lì)系數(shù)和修正系數(shù)，即第i個(gè)單元的輻射場(chǎng)計(jì)算式如式（2）。

對(duì)任意一個(gè)由N個(gè)單元組成的陣列，只考慮幅值，且自由空間滿足線性、均勻條件，可使用疊加原理，因此，根據(jù)式（2），陣列在任意一點(diǎn)( )r，θ，φ的輻射場(chǎng)可寫作式（3）。

在利用式（3）進(jìn)行相控陣天線的輻射場(chǎng)計(jì)算時(shí)，需要確定以下四個(gè)未知參量：

1）修正因子K；

2）激勵(lì)電流Iiejφi；

3）源點(diǎn)與場(chǎng)點(diǎn)之間距離Ri，一旦陣列固定，利用空間幾何關(guān)系容易得到該參數(shù)；

接下來(lái)，將首先討論參數(shù)K的求解方法。

如果場(chǎng)點(diǎn)P位于輻射場(chǎng)遠(yuǎn)場(chǎng)區(qū)內(nèi)最大輻射方向，此時(shí)滿足，陣列單元?dú)w一化場(chǎng)強(qiáng)f( )0，φ=1，每個(gè)單元到場(chǎng)點(diǎn)距離滿足近似條件，可設(shè)為R，根據(jù)式（3），場(chǎng)點(diǎn)P處電場(chǎng)強(qiáng)度幅值為

同時(shí)，該點(diǎn)處的電場(chǎng)強(qiáng)度還可用式（5）計(jì)算［12］：

其中Pt和Gt分別是相控陣天線的總發(fā)射功率和發(fā)射增益。聯(lián)立式（4）和式（5），可得式（6）。

在實(shí)際計(jì)算過(guò)程中，往往很難直接確定式（6）中的總發(fā)射功率Pt和發(fā)射增益Gt，此時(shí)，我們可以用以下方法進(jìn)行估計(jì)。

首先，每個(gè)陣列單元由已知電流Ii激勵(lì)，且在計(jì)算過(guò)程中認(rèn)為各個(gè)輻射單元是匹配的，其輸入阻抗ri為已知的，不考慮陣列單元的熱損耗，每個(gè)陣列單元的輸入功率全部轉(zhuǎn)化為單元的輻射功率（天線輻射效率為1），再利用疊加原理，可得陣列的總輻射功率如式（7）所示。

其次，陣列天線的方向性主要由陣列口徑形狀和陣列口徑照射確定，對(duì)有均勻幅相照射的口徑，平面口徑具有最大方向性，其中 A為陣列口徑面積，λ為工作波長(zhǎng)，當(dāng)一個(gè)平面口徑由大量單元組成（如陣列天線口徑）時(shí)，各個(gè)單元的激勵(lì)系數(shù)不同，這時(shí)往往需要定義一個(gè)稱之為口徑效率εA的參數(shù)來(lái)聯(lián)系Dmax和實(shí)際陣列方向性，根據(jù)文獻(xiàn)［13］，口徑效率的計(jì)算公式可寫作如式（9）所示，相應(yīng)的陣列方向性如式（8）所示，該式也就是陣列天線增益的計(jì)算公式。式（9）中，AFuniform_max和AFmax分別表示均勻激勵(lì)和非均勻加權(quán)激勵(lì)下方向圖函數(shù)的最大值，Ii_uniform和Ii分別均勻和非均勻加權(quán)的激勵(lì)系數(shù)。

至此，文中已給出了修正因子K；而激勵(lì)電流系數(shù)由待求陣列的輸入條件確定。在下一節(jié)將討論方向性函數(shù)fi( )θi，φi的求解方法。

3 中心陣元方向圖法

在自由空間中，天線單元的方向性函數(shù)是確定的；當(dāng)天線處于陣列環(huán)境中時(shí)，由于單元之間互耦效應(yīng)的影響，陣列中的天線單元方向性函數(shù)與自由空間中天線單元的方向性函數(shù)不同；同時(shí)，還會(huì)改變各陣列單元的激勵(lì)電流。

對(duì)大型陣列進(jìn)行精確的互耦分析是一個(gè)復(fù)雜、耗時(shí)的過(guò)程，在保證計(jì)算精度的同時(shí)，為了簡(jiǎn)化互耦分析，本節(jié)將采用中心陣元方向性函數(shù)方法得到在陣列環(huán)境中各單元的方向性函數(shù)。

在計(jì)算陣列天線輻射場(chǎng)時(shí)，可將陣列單元按所處的位置劃分為邊緣單元和非邊緣單元兩種。一方面，非邊緣位置的各陣列單元可近似看作處于相同的陣列環(huán)境中，認(rèn)為其單元方向性函數(shù)近似相等，可統(tǒng)一用中心陣列單元的方向性函數(shù)代替；另一方面，隨著陣列規(guī)模的擴(kuò)大，陣中屬于非邊緣單元的比例會(huì)逐漸增加，如7×7元陣列中非邊緣陣元數(shù)只占全陣單元總數(shù)中的51%，21×21元陣列中非邊緣陣元數(shù)占全陣單元總數(shù)的82%，而60×60元陣列中非邊緣陣元數(shù)已經(jīng)占到了93%。說(shuō)明在大型陣列中，非邊緣陣列單元處于主導(dǎo)地位，主要影響陣列的輻射場(chǎng)特性。因此，在計(jì)算大型陣列的輻射場(chǎng)過(guò)程中，可統(tǒng)一用中心陣列單元的方向性函數(shù)代替全陣任意位置的單元方向性函數(shù)同時(shí)保持不考慮互耦條件下的激勵(lì)電流不變。

中心陣元的方向性函數(shù)可以利用其有源單元方向圖（Active Element Pattern，AEP）代替。所謂AEP是指在陣列環(huán)境內(nèi)，僅該單元以單位源激勵(lì)、其余所有單元接匹配負(fù)載條件下的陣列總輻射場(chǎng)。而研究表明［1］，對(duì)于陣列環(huán)境中的指定單元而言，只有其鄰近的幾個(gè)單元對(duì)其方向性函數(shù)影響較大，因而，當(dāng)陣列足夠大時(shí)，中心單元的AEP已經(jīng)可以近似等于任意更大規(guī)模的陣列中心單元的AEP。因此，求解中心陣元的AEP，可采用一小規(guī)模的陣列進(jìn)行計(jì)算。

通過(guò)以上近似方法，將陣列復(fù)雜的互耦效應(yīng)全部用中心單元的方向性函數(shù)來(lái)代替，同時(shí)保持了計(jì)入互耦條件下與不計(jì)入互耦條件下的陣列激勵(lì)電流相同，簡(jiǎn)化了大型陣列天線輻射場(chǎng)的求解過(guò)程。

4 數(shù)值驗(yàn)證

4.1 陣列單元

選取微帶貼片天線作為大型陣列天線的組成單元，該天線的FEKO模型如圖2（a）所示，其駐波比曲線如圖2（b）所示，其E面和H面方向圖如圖3所示。選擇3.1GHz為天線的中心頻率。

圖2天線單元輻射性能示意圖

圖3 天線單元方向圖（E面和H面）

4.2 輻射場(chǎng)求解模型

通過(guò)對(duì)由圖2（a）所示單元構(gòu)成的小陣面分析，當(dāng)陣列規(guī)模分別為7×7和9×9時(shí)，兩者的中心陣元AEP差別已經(jīng)不大。因此構(gòu)建如圖4所示的7×7典型小陣面，圖中①號(hào)單元屬于中心陣元，②、③、④號(hào)單元均屬于邊緣陣元，計(jì)算①號(hào)陣元AEP，由該陣元AEP替代大型陣列天線任意位置單元的方向性函數(shù)fi( )θ，φ，其E面和H面的方向圖如圖5所示。比較圖3和圖5，單元在陣列環(huán)境中與自由空間環(huán)境中的單元方向圖不同，主要差別在于陣列環(huán)境中的單元方向圖隨角度的變化不再呈單調(diào)趨勢(shì)，而是有較小的波動(dòng)現(xiàn)象，這主要是由于其余單元對(duì)該單元的互耦引起的。

圖4 7×7典型小陣面模型

圖5 中心陣元AEP（H面和E面）

下面，構(gòu)建一個(gè)15×15的平面陣，陣列FEKO模型如下：

1）陣列布置形式為矩形柵格口徑，如圖3（a）所示，口徑面大小為0.7307×0.7462 m2。

2）工作頻率3.1GHz，陣列被均勻激勵(lì)，輻射總功率1W。

采用文中模型計(jì)算的相關(guān)參數(shù)如下：

1）根據(jù)功率關(guān)系式（7）得到均勻激勵(lì)電流約為9.4×10-3A。

2）根據(jù)式（8），方向性約為28.6dB。

3）修正因子K=113.0762。

4）方向性函數(shù)為圖4中①號(hào)單元AEP，如圖5。

根據(jù)式（1），其遠(yuǎn)場(chǎng)距離約在離陣列中心22.556m處。其輻射場(chǎng)求解模型如圖6所示，分別采用FEKO和基于Matlab的文中方法計(jì)算：

1）主軸：x=0，y=0，z=0.05m～50m。

2）H 面：x=0.5m，y=0，z=0.05m～50m；x=2m，y=0，z=0.05m～50m。

3）E面：x=0，y=0.5m，z=0.05m～50m；x=0，y=2m，z=0.05m～50m。

圖6 輻射場(chǎng)求解模型

為了研究陣列天線的近場(chǎng)輻射特性，首先定義陣列偏軸這一概念：所謂偏軸，即是在三維直角坐標(biāo)系中，以陣列所在平面為xoy面，保持x、y坐標(biāo)不變，場(chǎng)點(diǎn)沿z軸變化所構(gòu)成的一條直線。當(dāng)x、y、z位于陣列主波束照射范圍時(shí)，則稱為主軸。當(dāng)x、y、z相對(duì)于陣列的遠(yuǎn)場(chǎng)距離比較小時(shí)，陣列的偏軸輻射場(chǎng)很好地反映了陣列的近場(chǎng)輻射特性。采用本文方法和FEKO計(jì)算方法得到的仿真結(jié)果如圖7～11所示。

圖7 主軸輻射場(chǎng)計(jì)算結(jié)果對(duì)比

4.3 15×15陣列模型驗(yàn)證

從圖7～圖11分析可見：

1）采用中心陣元方向圖法，計(jì)算得到的主軸上的輻射場(chǎng)精度較高，與矩量法計(jì)算結(jié)果相比誤差小于0.5dB（不考慮零陷處值的誤差）。

2）采用中心陣元方向圖法計(jì)算偏軸輻射場(chǎng)，計(jì)算結(jié)果顯示精度較高，與矩量法計(jì)算結(jié)果相比，其峰值場(chǎng)強(qiáng)分布誤差小于1dB，零陷處場(chǎng)強(qiáng)之間的誤差稍大一些；但兩者的計(jì)算結(jié)果的變化規(guī)律呈現(xiàn)出較好的一致性，這有利于我們采用文中方法對(duì)輻射場(chǎng)的變化規(guī)律進(jìn)行分析。

圖8 均勻激勵(lì)情況下，偏軸0.5m處電場(chǎng)分布曲線（H面）

3）圖中分別計(jì)算了偏離主軸0.5m和2m時(shí)的輻射場(chǎng)分布，結(jié)果顯示，隨著偏軸位置遠(yuǎn)離主軸，離陣列較近場(chǎng)點(diǎn)處的計(jì)算值與基于矩量法的結(jié)果相比呈現(xiàn)零陷偏移的現(xiàn)象，而隨著z值的增加，場(chǎng)點(diǎn)逐漸靠近陣列天線主瓣范圍，兩者誤差減小。

5 誤差分析

如圖10所示，計(jì)算了z=5m時(shí)，x距離不同偏軸位置處，基于FEKO和本文方法的計(jì)算結(jié)果，圖中證實(shí)了隨著偏軸距離增加，零陷偏移增大。不難得出，當(dāng)y和z值固定，隨著x值增加，場(chǎng)點(diǎn)位置的θ角逐漸增加，這時(shí)陣列邊緣單元對(duì)輻射場(chǎng)的貢獻(xiàn)增加，采用中心陣元的方向圖計(jì)算會(huì)帶來(lái)較大誤差。

再者，隨著場(chǎng)點(diǎn)位置的θ角增大，天線單元方向圖變化率增大，如圖5所示，影響了計(jì)算過(guò)程中的對(duì)方向圖的取值，同樣會(huì)給計(jì)算帶來(lái)誤差。

圖9 均勻激勵(lì)情況下，偏軸2m處電場(chǎng)分布曲線

但是，從圖7～10整體不難看出，無(wú)論待求場(chǎng)點(diǎn)位于空間何處，各峰值之間的誤差依然在6dB范圍內(nèi)，滿足工程應(yīng)用的要求。

綜上所述，當(dāng)陣列規(guī)模達(dá)到15×15時(shí)，采用中心陣元方向圖法計(jì)算的結(jié)果已經(jīng)具有較高的精度。對(duì)大規(guī)模陣列，可以采用該方法進(jìn)行輻射場(chǎng)計(jì)算，以尋找輻射場(chǎng)的變化規(guī)律，為復(fù)雜平臺(tái)上陣列天線的電磁兼容分析提供理論依據(jù)。

6 結(jié)語(yǔ)

大型陣列天線的近場(chǎng)計(jì)算對(duì)于實(shí)現(xiàn)陣列天線電磁兼容分析具有重要意義，同時(shí)也是電磁兼容計(jì)算過(guò)程中的難點(diǎn)之一，想要進(jìn)行精確的計(jì)算是比較困難的。本文從工程應(yīng)用的角度出發(fā)，探討了一種考慮互耦、快速的大型陣列天線近場(chǎng)計(jì)算方法，并進(jìn)行驗(yàn)證。驗(yàn)證結(jié)果表明，該方法能夠有效計(jì)算場(chǎng)區(qū)范圍內(nèi)的電場(chǎng)峰值，可以為大型陣列天線的電磁兼容分析提供數(shù)據(jù)支撐。同時(shí)，該方法還可拓展應(yīng)用到三角柵格任意口徑形式的情況下。