任立剛， 朱武兵， 高中峰，周茜

（1．海軍92925部隊(duì)，山西長治046011; 2．海軍工程大學(xué)電氣工程學(xué)院，武漢 430033;3．海軍91446部隊(duì), 河北諑州 072750）

0 引 言

隨著現(xiàn)代信息處理技術(shù)、微電子技術(shù)和傳感器技術(shù)的發(fā)展，對艦船目標(biāo)特性的研究進(jìn)一步深入，艦船軸頻電磁場已經(jīng)成為重要的水下信號。文獻(xiàn)[1]通過研究表明深海環(huán)境下艦船軸頻電磁場可以通過運(yùn)動的垂直時諧偶極子模擬，但是只研究了艦船是單個螺旋槳的情形下的建模與仿真?，F(xiàn)役的中型、大型艦艇都可能采用多槳的推進(jìn)方式，因此只研究單槳艦艇的軸頻電磁場顯然是不夠的。對于四槳的大型艦船，在采用多槳的推進(jìn)方式時，四個槳經(jīng)常是不同步的，這樣就會導(dǎo)致四個螺旋槳旋轉(zhuǎn)時各自產(chǎn)生的軸頻電磁場具有不同的相位、幅值、頻率。因此研究四槳運(yùn)動艦船產(chǎn)生的軸頻電磁場具有十分重要的實(shí)際意義。對軸頻電磁場的實(shí)際應(yīng)用需要對這種四槳的艦船的電磁場進(jìn)行建模。并為未來軸頻電磁場用于艦艇的遠(yuǎn)程探測提供了基礎(chǔ)。本文采用四個運(yùn)動垂直時諧電偶極子在固定場點(diǎn)產(chǎn)生的電磁場對深海中運(yùn)動四槳艦船的電磁場進(jìn)行模擬。

1 深海中運(yùn)動垂直時諧偶極子電磁場表達(dá)式

根據(jù)文獻(xiàn)[1]推導(dǎo)出深海中運(yùn)動垂直時諧電偶極子電磁場的表達(dá)式：

式（1）～式（8）中的相關(guān)參數(shù)及物理意義與文獻(xiàn)[1]一致。運(yùn)用基于漢克爾變換的 FFT算法[9]得到了電磁場三分量的頻域曲線和時域曲線。那么就可以得到任意時刻t電磁場的頻域值和時域值記t時刻的電磁場三分量的頻域值和時域值為：

2 四個運(yùn)動垂直時諧電偶極子電磁場表達(dá)式

這里假設(shè)采用4個處于同一位置的運(yùn)動垂直時諧偶極子對軸頻電磁場進(jìn)行建模。時諧電偶極子與單個電偶的極子表示方法類似?？赏茖?dǎo)出四個運(yùn)動垂直時諧電偶極子在固定場點(diǎn)產(chǎn)生的電磁場的表達(dá)式如下：對于電場x分量：

對于磁場x分量：

以上給出固定點(diǎn)電場x分量、磁場x分量的電磁場的表達(dá)式，其他電磁場分量的公式與之類似。在這里為了節(jié)省篇幅，不作贅述。

3 相位不同步的四槳艦船的軸頻電磁場仿真

對于4槳的運(yùn)動艦船在航行的過程中，假設(shè)四個槳的轉(zhuǎn)速相同，但是四個槳的相位不同，艦船的運(yùn)行速度為5 m/s。這里可以用四個頻率相同的運(yùn)動垂直時諧偶極子對艦船的軸頻電磁場進(jìn)行建模。建立文獻(xiàn)[1]的洛倫茲模型，空氣、海水、海床的介質(zhì)參數(shù)同文獻(xiàn)[1]。各時諧運(yùn)動的垂直電偶極子的頻率為 1 Hz，各電偶極子的初始坐標(biāo)（x=-1000 m，y=0，z=10 m），終點(diǎn)坐標(biāo)（x=-1000 m，y=0，z=10 m），強(qiáng)度均為1 A*m。 本文主要給出x分量的仿真圖（一方面已經(jīng)可以說明問題、另一方面節(jié)省篇幅）。

通過分析仿真結(jié)果圖，可以知道：

1）分析圖 1(a)、(c)可知，相角的變化并不影響頻域的幅值。

2） 通過分析四個槳相位均為0時的仿真結(jié)果圖與單槳相位為0時的仿真結(jié)果圖1(a)、(b)、(e)、(f)?？梢园l(fā)現(xiàn)四個相同的電偶極子的強(qiáng)度產(chǎn)生的電磁場的時域值是單個同參數(shù)的電偶極子的4倍。由此可以說明這四個1 A*m電偶極子可以等效為強(qiáng)度為4 A*m的電偶極子。

3） 對比兩組不同相位的仿真結(jié)果，可以發(fā)現(xiàn)當(dāng)四個電偶極子的相位均為0時產(chǎn)生的電磁場最大，圖1(c)、(d)所示的情況產(chǎn)生的電磁場的數(shù)量級非常小，以現(xiàn)在的電磁場探測水平可以把這種情況產(chǎn)生的電磁場視為零?；旧峡梢哉J(rèn)為實(shí)現(xiàn)了軸頻電磁場的抵消。說明了軸頻電磁場由于相位不同而存在抵消效應(yīng)。

4 轉(zhuǎn)速不同的四槳艦船軸頻電磁場仿真

通過分析仿真結(jié)果，可以知道：

1）分析仿真結(jié)果的頻域圖2(b)（頻率越高，電場值越大），發(fā)現(xiàn)在電場z方向由于頻率的不同導(dǎo)致有較大的幅值差異。在頻域圖 2(a)中發(fā)現(xiàn)有的沒有明顯變化，這是因?yàn)轭l率對軸頻電磁場的影向并不明顯。

2）從圖2(c)軸頻電磁場時域圖可以發(fā)現(xiàn)，電磁場的包絡(luò)線沒有頻率相同時光滑，這說明頻率對軸頻電磁場的時域值有較大的影響。

3）分析仿真結(jié)果圖1(a)、(b)（相角為0，頻率均為1 Hz），圖2(a)、(b)，提出當(dāng)頻率不同時的軸頻電磁場的建模思路：首先求出兩個槳產(chǎn)生的不同的軸頻電磁場頻率，然后采用具有這兩個頻率不同的兩個虛擬時諧偶極子對軸頻電磁場進(jìn)行建模。

5 結(jié)論

本文通過仿真運(yùn)動垂直電偶極子產(chǎn)生的電磁場對軸頻電磁場進(jìn)行了模擬。模擬得到的結(jié)果不同相位和不同轉(zhuǎn)速時四槳運(yùn)動艦船產(chǎn)生的軸頻電磁場，并對產(chǎn)生的軸頻電磁場進(jìn)行了分析。分析得到了四槳艦船產(chǎn)生的軸頻電磁場特性以及為軸頻電磁場的建模提供了思路。值得說明的是，可能實(shí)際艦船不會出現(xiàn)完全反相得情況，或者頻率不會相差這么大。但是在理想情況進(jìn)行的仿真，對艦船產(chǎn)生的軸頻電磁場有了以上這些有價(jià)值的結(jié)論。

[1]毛偉，周萌，余刃. 兩層介質(zhì)中運(yùn)動垂直時諧偶極子產(chǎn)生的電磁場[J]. 武漢理工大學(xué)學(xué)報(bào)，2011，35(5):1081-1085.

[2]Weaver J T. The quasi-static field of an electric dipole embedded in a two-layer conducting hal-f-space. Can.J. Phys., 1967, 45: 1981-2002.

[3]羅耀煌, 林公源, 包德修. 諧變水平電偶極子在四分層媒質(zhì)中的電磁波[J]. 云南民族學(xué)院學(xué)報(bào)(自然科學(xué)版), 1997, 6(2): 18 - 27.

[4]雷銀照. 時諧電磁場解析方法[M]. 北京: 科學(xué)出版社, 2000.

[5]柳超, 梁高權(quán). 用氣－水－土 3層模型分析海水中水平電偶極子源的電磁場分布.

[6]胡俊, 聶在平. 索末菲爾德積分新方法—快速漢克爾變換. 電子學(xué)報(bào), 1998, 26(3): 126-128.

[7]劉文寶. 基于電偶極子模型的艦船靜電場換算研究[D]．武漢：海軍工程大學(xué)，2010．