楊 超

(武漢船舶設(shè)計(jì)研究院有限公司，湖北 武漢 430064)

0 引 言

潛艇圍殼區(qū)域空間有限，天線布置緊湊。天線強(qiáng)輻射近場(chǎng)不僅會(huì)對(duì)人員造成電磁安全性危害，也會(huì)使其他天線、設(shè)備受到電磁能量阻塞甚至毀傷；在載體的影響下，天線輻射方向圖會(huì)發(fā)生畸變。在潛艇總體論證和設(shè)計(jì)初期，電磁仿真預(yù)測(cè)作為獲取電磁安全性數(shù)據(jù)、預(yù)報(bào)天線布局方案可行性以及指導(dǎo)優(yōu)化設(shè)計(jì)的重要手段，對(duì)于有針對(duì)性地確定效費(fèi)比最佳的圍殼天線布局和優(yōu)化設(shè)計(jì)方案具有重要的意義。

20世紀(jì)末，美國(guó)洛氏國(guó)際海試科學(xué)中心提出了潛艇電磁兼容系統(tǒng)工程方案，在這個(gè)方案中首先提出了采用計(jì)算機(jī)技術(shù)對(duì)潛艇進(jìn)行電磁設(shè)計(jì)。美國(guó)海上系統(tǒng)司令部干舷設(shè)計(jì)小組也為美海軍潛艇綜合電磁設(shè)計(jì)開(kāi)發(fā)了一些軟件系統(tǒng)，建立了許多潛艇電磁環(huán)境分析的計(jì)算機(jī)模型。美海軍海上系統(tǒng)司令部同海軍指揮、控制和海洋監(jiān)測(cè)中心試制部以及洛氏國(guó)際公司一起研制開(kāi)發(fā)了二維和三維科學(xué)數(shù)據(jù)顯象工具，利用這些計(jì)算機(jī)圖形顯象技術(shù)幫助潛艇設(shè)計(jì)者進(jìn)行直觀的優(yōu)化設(shè)計(jì)。

1 建模及計(jì)算方法

本文參考Suboff潛艇進(jìn)行幾何建模[1-5]，建立的潛艇幾何模型如圖1所示，包含主體和指揮臺(tái)圍殼。Suboff項(xiàng)目由美國(guó)國(guó)防高等研究計(jì)劃署（DARPA）提出，其潛艇模型由軸對(duì)稱體、圍殼、尾翼等組合而成。

Suboff模型總長(zhǎng)4.356 m，最大直徑0.508 m。指揮臺(tái)圍殼長(zhǎng)0.368 m，高0.206 m。計(jì)算時(shí)將其放大20倍。

矩量法（MoM）適合求解線面結(jié)合問(wèn)題，其最早用來(lái)求解線問(wèn)題，最后推廣用于求解三維問(wèn)題。這類方法的一個(gè)優(yōu)點(diǎn)是采用面剖分的方法，求解物體表面的感應(yīng)電流；不需要設(shè)置邊界條件，不需要對(duì)全部計(jì)算空間進(jìn)行體剖分；不存在數(shù)值色散誤差；適合處理弱非均勻性目標(biāo)。1982年，S.M.Rao，D.R.Wilton和A.W.Glisson提出的RWG基函數(shù)極大地推動(dòng)了矩量法（MoM）的發(fā)展。RWG基函數(shù)能夠很好地描述表面電流的物理特性，而且其基于平面三角形網(wǎng)格，可以對(duì)任意形狀目標(biāo)進(jìn)行剖分，是目前應(yīng)用最多的一種基函數(shù)。為了加快算法的計(jì)算效率，提出了高階基函數(shù)和曲面基函數(shù)等。Lu，Song和Chew提出了MLFMA，使計(jì)算效率大大提高，計(jì)算機(jī)內(nèi)存需求和計(jì)算時(shí)間降至O(NlogN)。在此基礎(chǔ)上，Chew與Cui引入快速遠(yuǎn)場(chǎng)近似和射線傳播技術(shù)，使MLFMA的計(jì)算效率進(jìn)一步提高。

圖1 建立的潛艇幾何模型Fig.1 Geometric model of submarine

根據(jù)導(dǎo)體表面的邊界條件，由式（1）得電場(chǎng)積分方程：

將式（3）代入式（2），并運(yùn)用伽略金法，得到下面的矩陣方程：

2 仿真

計(jì)算時(shí)網(wǎng)格剖分模型如圖2所示。在圍殼上布置一根單極短波天線，該單極天線發(fā)射時(shí)，艇體表面電流分布仿真結(jié)果如圖3所示。

圖2 網(wǎng)格剖分模型Fig.2 Mesh splitting model

圖3 艇體表面電流分布仿真結(jié)果Fig.3 Simulation results of current distribution on submarine surface

設(shè)置無(wú)窮大導(dǎo)電平面模擬海水，半潛狀態(tài)下輻射方向圖仿真結(jié)果如圖4所示。

圖4 半潛時(shí)輻射方向圖仿真結(jié)果Fig.4 Simulation results of radiation pattern at half latency

半潛狀態(tài)下圍殼上單極短波天線發(fā)射時(shí)，輻射電場(chǎng)分布云圖如圖5所示。

圖5 半潛時(shí)輻射電場(chǎng)分布云圖Fig.5 Radiation electric field distribution nephogram at half latency

艇體潛入水下3/4，圍殼上單極短波天線發(fā)射時(shí)輻射方向圖如圖6所示。

圖6 潛入水下3/4時(shí)輻射方向圖仿真結(jié)果Fig.6 Simulation results of radiation pattern at 3/4

艇體潛入水下3/4，圍殼上單極短波天線發(fā)射時(shí)輻射場(chǎng)強(qiáng)分布云圖如圖7所示。

圖7 潛入水下3/4時(shí)輻射電場(chǎng)分布云圖Fig.7 Radiation electric field distribution nephogram at 3/4

3 仿真分析

在上述輻射特性仿真的基礎(chǔ)上，模擬了潛艇不同下潛狀態(tài)下在垂直剖面和水平剖面上的輻射電場(chǎng)分布云圖，將短波天線周圍的遠(yuǎn)場(chǎng)和近場(chǎng)數(shù)值以更加直觀的顏色區(qū)分加以體現(xiàn)，可為潛艇電磁兼容性設(shè)計(jì)提供支撐。

建立如圖8所示坐標(biāo)系，設(shè)艇首指向艇尾為x軸正向，潛艇右舷為y軸正向，潛艇左舷為y軸負(fù)向，短波天線位置為原點(diǎn)。

圖8 潛艇模型坐標(biāo)軸模型Fig.8 Submarine model coordinate axis model

依據(jù)潛艇模型不同潛態(tài)的仿真結(jié)果[7]，可得電場(chǎng)強(qiáng)度分布曲線圖如圖9所示。

圖9 短波天線輻射x軸電場(chǎng)強(qiáng)度Fig.9 Electric field strength of x-axis radiated by short-wave antenna

圖10 短波天線輻射y軸電場(chǎng)強(qiáng)度Fig.10 Electric field strength of y-axis radiated by short-wave antenna

由圖9和圖10可看出，輻射電場(chǎng)X分量以艇體y=0的縱剖面為對(duì)稱平面成類似對(duì)稱分布，天線正上方為正峰值，依次向首尾方向衰減；輻射電場(chǎng)Y分量以艇體x=0的橫剖面為對(duì)稱平面成反對(duì)稱分布，在天線正上方為一中間均值，在兩舷側(cè)對(duì)稱各有一個(gè)反向峰值。仿真結(jié)果表明輻射感應(yīng)電場(chǎng)具有明顯的空間分布特性，且艇體周圍區(qū)域輻射電場(chǎng)具有一定的規(guī)律性，本文僅試驗(yàn)了單支天線的輻射電場(chǎng)分布特性，在多支天線同時(shí)存在時(shí)互相增益或衰減特性則有所不同，艦載雷達(dá)天線間的散射、繞射、諧振和耦合現(xiàn)象也會(huì)影響最終的電場(chǎng)分布結(jié)果。本文仿真分析結(jié)果與其他文獻(xiàn)提供的實(shí)測(cè)結(jié)果[8]相比較，具有較好的一致性。

4 結(jié) 語(yǔ)

在潛艇論證設(shè)計(jì)階段，電磁兼容優(yōu)化布局需要將圍殼上的天線系統(tǒng)進(jìn)行綜合考慮，進(jìn)行合理的布置，以避免相互影響。借助仿真預(yù)測(cè)技術(shù)，對(duì)各種天線布置方案的性能優(yōu)劣進(jìn)行評(píng)估，可大大減小工作量，最后確定最佳的方案。在仿真預(yù)測(cè)中，還有許多關(guān)鍵技術(shù)有待突破，如圍殼復(fù)合材料的建模與仿真，海面及海浪的模擬，預(yù)測(cè)精度的試驗(yàn)驗(yàn)證等。