李煒昕, 魏維偉, 涂 建, 瞿 寧

(1.上海無線電設(shè)備研究所,上海201109;2.上海市目標(biāo)識別與環(huán)境感知工程技術(shù)研究中心,上海201109)

0 引言

靜電探測[1-3]是通過感應(yīng)目標(biāo)與探測單元之間的靜電場而獲得目標(biāo)靜電信息的一種探測方法?；诒粍邮届o電感應(yīng)機理的探測系統(tǒng)具有結(jié)構(gòu)簡單、反隱身和抗電磁干擾能力強等優(yōu)勢,在氣固兩相流參數(shù)測量、航空發(fā)動機狀態(tài)監(jiān)測和空中飛行目標(biāo)探測等領(lǐng)域具有廣泛的應(yīng)用前景[4-6]。

目前,國內(nèi)外對靜電目標(biāo)方位探測開展了大量研究工作。中科院靜電課題組[7]利用微機電系統(tǒng)(micro-electro-mechanical system,MEMS)電場敏感芯片獲取大氣電場信息,但該方法存在響應(yīng)速度慢等缺點。北京理工大學(xué)靜電探測課題組[8]采用平面圓陣被動式靜電探測方法進行目標(biāo)定位,理論上遠(yuǎn)距離定位精度可達(3～15)%,但在近距情況下方位解算誤差大,適用于地面靜電探測系統(tǒng)。德國的TRINKS等[1]利用正交靜電探測陣列電極,提出基于陣列信號相位的測向算法,但是需要利用目標(biāo)特征曲線中過零點的信息計算方位,適用于探測目標(biāo)飛行軌跡。

針對靜電探測系統(tǒng)小型化、高精度以及實時探測等特點,本文基于靜電感應(yīng)和靜電平衡機理,分析靜電場探測數(shù)學(xué)模型和電荷靈敏理論模型,建立陣列電極的靜電場矢量探測數(shù)學(xué)模型。在此基礎(chǔ)上,通過數(shù)學(xué)仿真和模擬試驗等方法,驗證靜電場矢量探測的方位解算精度。

1 靜電探測系統(tǒng)工作原理

靜電探測系統(tǒng)是基于靜電感應(yīng)和靜電平衡的原理進行設(shè)計的。處于電場中的孤立導(dǎo)體,外加電場將對導(dǎo)體中的自由電子產(chǎn)生力的作用,使電子逆著電場的方向運動,這種現(xiàn)象稱為靜電感應(yīng)。靜電感應(yīng)可以使原來不帶電的導(dǎo)體發(fā)生電荷的重新分布,并達到靜電平衡狀態(tài),該過程通常都是在極短時間內(nèi)完成的(約為10-19s)[9-11]。

靜電感應(yīng)探測原理如圖1所示。設(shè)目標(biāo)與感應(yīng)電極的垂直距離為z,感應(yīng)電極表面積為S,目標(biāo)帶電量為+Q,目標(biāo)運動時間為t。目標(biāo)以速度v從感應(yīng)電極正上方通過,可以得到感應(yīng)電流i。

圖1 靜電感應(yīng)探測原理

根據(jù)高斯定理,可得到感應(yīng)電極表面的電荷密度

式中:ε0為真空絕對介電常數(shù);εr為相對介電常數(shù);Ean為感應(yīng)電極表面垂直方向上的電場強度。

則感應(yīng)電極表面的電荷密度為

式中:x、y、z分別為目標(biāo)相對于感應(yīng)電極中心的位置坐標(biāo)。

因此,感應(yīng)電極靠近荷電目標(biāo)一側(cè)所帶的異號電荷總量為

可得到探測電極所獲取的感應(yīng)電荷量

式中:a為圓形探測電極半徑。

2 基于陣列電極的目標(biāo)方位解算模型

在載體表面布設(shè)靜電探測器感應(yīng)電極陣列,建立目標(biāo)方位解算模型,通過電極陣列得到空中目標(biāo)靜電場在探測器附近三維方向的感應(yīng)分量,結(jié)合測向算法進行推算,可以得到空中目標(biāo)的方位信息。擬采用靜電場矢量方位解算方法,如圖2所示。

圖2 靜電目標(biāo)方位解算算法示意圖

將靜電目標(biāo)視為oxyz空間中的一個點電荷,其攜帶電量為Q;A、B和C、D為感應(yīng)電極編號;γ為目標(biāo)方位角。假設(shè)在某一時刻,目標(biāo)處于T點,到坐標(biāo)原點的距離為R,此時由于帶電目標(biāo)靜電場的作用,兩對靜電感應(yīng)電極極板間會產(chǎn)生不同的電勢差。電場在極板對上產(chǎn)生的電勢差分別記為UAB和UCD,表達式為

式中:E為目標(biāo)電場強度矢量;LAB、LCD分別表示電極A、B和電極C、D兩對平行極板間的距離矢量。結(jié)合電場強度公式和探測系統(tǒng)的幾何關(guān)系,可獲取極板對間的電勢差分別為

式中:α為偏航角;β為俯仰角;L為平行板電極間的距離。利用靜電探測系統(tǒng)測得UAB和UCD,則可求出靜電感應(yīng)電極附近的電場強度,獲取目標(biāo)的脫靶方位信息為

因此,采用正交平行布設(shè)的兩組陣列感應(yīng)電極,復(fù)合利用各個電極感應(yīng)電荷幅值特征,就能夠?qū)崿F(xiàn)對任意時刻帶電體方位的高精度解算。

3 數(shù)值計算分析

有限元方法適用于邊界形狀或邊界條件復(fù)雜及介質(zhì)分布復(fù)雜的問題的求解。由于這是一種數(shù)值解法,不受場域邊界形狀的限制,因此特別適用于求解本研究中涉及的靜電場電荷分布問題。

3.1 數(shù)值計算模型

由于耦合模型需要進行靜態(tài)和參數(shù)化動態(tài)分析,且形狀相對復(fù)雜,簡化為二維模型較困難,故采用三維模型進行電場仿真計算。利用有限元仿真分析軟件Maxwell,建立數(shù)值計算模型。設(shè)目標(biāo)為直徑400 mm的鋁質(zhì)小球,沿z軸正方向運動;感應(yīng)電極由兩組銅板構(gòu)成,厚度為1 mm,兩組電極間距0.4 m,每組電極呈正交分布;求解域為矩形空氣場;小球激勵源為10 k V電壓源,感應(yīng)電極和導(dǎo)彈均設(shè)置為懸浮狀態(tài)。

3.2 方位解算精度分析

對本文提出的靜電目標(biāo)方位解算算法的解算精度進行數(shù)值分析,靜電目標(biāo)方位探測的仿真模型定義如圖3所示。四個感應(yīng)電極呈正交方式分布于oyz平面,感應(yīng)電極按逆時針方向順序命名,即感應(yīng)電極A、C、B、D;從彈頭方向看,象限命名遵循傳統(tǒng)數(shù)學(xué)象限規(guī)則。其中,靜電目標(biāo)位于坐標(biāo)軸z軸負(fù)半軸上,初始位置與坐標(biāo)原點相距8 m。

圖3 仿真模型定義規(guī)則

進行參數(shù)化建模計算,利用麥克斯韋(Maxwell)場計算公式獲取陣列感應(yīng)電極的感應(yīng)電荷總量,在1 000 m/s的交會速度條件下,感應(yīng)電荷變化曲線如圖4所示。從計算結(jié)果可以看出,在模擬彈目交會過程中,陣列中距離目標(biāo)較近的電極感應(yīng)電荷總量較大,在脫靶點電荷量達到最大值,感應(yīng)電流為零;y軸方向的兩對電極(電極C和電極D)由于在交會過程中與目標(biāo)的距離始終相等,計算獲取的感應(yīng)電荷量變化曲線趨于一致;由于導(dǎo)彈彈體(結(jié)構(gòu))的影響,目標(biāo)空間靜電場分布會發(fā)生改變,導(dǎo)致感應(yīng)電極電荷總量存在差異。

圖4 感應(yīng)電荷交會特征曲線

根據(jù)文中所建立的脫靶方位解算模型,對數(shù)值仿真獲取的感應(yīng)電荷數(shù)據(jù)進行數(shù)學(xué)計算,可以得到目標(biāo)方位解算曲線,如圖5所示。可知,方位角解算結(jié)果為90°,與彈目交會模擬過程一致。在有限元數(shù)值仿真交互過程中,由于采用了區(qū)域精細(xì)化網(wǎng)格劃分方法,通過理論計算獲取的目標(biāo)靜電場信息具有良好的一致性,因此理論計算獲取的目標(biāo)方位具有很好的探測精度。

圖5 目標(biāo)方位解算曲線

4 試驗驗證

為驗證靜電探測理論和數(shù)值計算的正確性,進行了模擬彈目交會的準(zhǔn)動態(tài)試驗。試驗系統(tǒng)由靜電探測系統(tǒng)、目標(biāo)模擬單元和數(shù)據(jù)采集系統(tǒng)構(gòu)成,其中目標(biāo)模擬單元由懸掛裝置、靜電目標(biāo)、高壓電纜、靜電高壓發(fā)生器和滑軌組成,如圖6所示。試驗?zāi)繕?biāo)為直徑400 mm的球體,試驗交會姿態(tài)為:脫靶量4.04 m,位置誤差0.18 m,偏移0.04 m,交會平均速度13 m/s。

圖6 試驗裝置

目標(biāo)相對探測系統(tǒng)的方位角為135°,交會過程中靜電探測系統(tǒng)輸出信號如圖7所示。目標(biāo)在1.35 s時到達脫靶點,交會過程中相對速度逐漸增加,探測系統(tǒng)所獲取的電壓幅值逐漸增高。

圖7 靜電探測系統(tǒng)輸出信號

為驗證探測系統(tǒng)角度解算精度,開展準(zhǔn)動態(tài)外場試驗。試驗?zāi)繕?biāo)方位角范圍為0°～360°,交會軌跡沿x軸方向。典型目標(biāo)方位解算曲線如圖8所示,解算數(shù)據(jù)如表1所示。

表1 目標(biāo)方位解算數(shù)據(jù)

圖8 目標(biāo)方位解算曲線

由試驗結(jié)果可知:

a)基于陣列電極的目標(biāo)方位解算精度具有較好的一致性,角度誤差的浮動較小;

b)在接近脫靶位置過程中所獲取的角度精度較高,方位解算精度優(yōu)于10°;

c)靜電探測系統(tǒng)輸出信號幅度越大,信干比越大,系統(tǒng)解算精度越高,誤差越小;

d)方位解算精度穩(wěn)定性有待提高,主要影響因素包括目標(biāo)位置誤差、靜電探測系統(tǒng)誤差和環(huán)境靜電場隨機誤差等。

5 結(jié)論

基于靜電感應(yīng)和靜電平衡機理,建立陣列電極的靜電場矢量探測數(shù)學(xué)模型,通過數(shù)學(xué)仿真和模擬試驗等方法,驗證了靜電場矢量探測的方位解算精度。理論分析和試驗結(jié)果表明,在接近脫靶位置的過程中系統(tǒng)所獲取的角度精度較高,方位解算精度優(yōu)于10°。后續(xù)研究中可在自適應(yīng)匹配濾波和數(shù)據(jù)處理優(yōu)化算法等方面進行深入探討,進一步提高靜電探測系統(tǒng)方位解算精度。