郭安新,王學(xué)智,劉少偉,杜翔宇

(空軍工程大學(xué) 防空反導(dǎo)學(xué)院,西安 710051)

0 引言

電磁軌道發(fā)射是利用電磁推力將負(fù)載加速至超高速的發(fā)射方式[1-4],具有響應(yīng)速度快,發(fā)射性能穩(wěn)定,推力可控等優(yōu)勢(shì)[5,6],電磁軌道炮的彈丸正從常規(guī)動(dòng)能彈擴(kuò)展到智能彈藥[7-9]。電磁軌道炮在發(fā)射過(guò)程中產(chǎn)生的磁場(chǎng)極易影響智能彈藥中元器件的正常工作,特別是軌道炮電流產(chǎn)生的強(qiáng)電磁[10],過(guò)強(qiáng)的磁場(chǎng)會(huì)嚴(yán)重?fù)p壞其電子元器件的工作性能并使其破壞。除此之外,軌道炮在工作時(shí),樞軌接觸面內(nèi)的電流集中分布會(huì)導(dǎo)致軌道局部焦耳熱過(guò)高[11-14],造成局部熱腐蝕,影響連續(xù)發(fā)射。

基于此,國(guó)內(nèi)外學(xué)者做了大量的研究,燕山大學(xué)劉暢[15]采用CircuitEditor電路作為軌道炮激勵(lì)電源,研究了接觸電阻對(duì)電磁炮電磁場(chǎng)的影響;李騰達(dá)[16]提出了銅基復(fù)合型四極電磁發(fā)射器模型并與普通型四極電磁發(fā)射器的電磁特性進(jìn)行對(duì)比分析,發(fā)現(xiàn)銅基復(fù)合型四極電磁發(fā)射器的磁場(chǎng)屏蔽效果更好;薛新鵬[17]對(duì)六極軌道發(fā)射器的推進(jìn)力與磁屏蔽效果進(jìn)行了量化研究;劉明[18]采用有限元分析法對(duì)比分析了不同軌道結(jié)構(gòu)下電樞電流分布特性;任師達(dá)[19]采用電磁-結(jié)構(gòu)耦合方法對(duì)雙曲電樞四軌電磁發(fā)射器的電接觸特性進(jìn)行了研究。

基于以上研究可以發(fā)現(xiàn):對(duì)電磁發(fā)射裝置電磁特性的分析而言,對(duì)多極電磁軌道發(fā)射器的磁場(chǎng)屏蔽效果研究較少;大多數(shù)學(xué)者都只研究對(duì)比分析相同軌道數(shù)的電磁發(fā)射器,對(duì)不同軌道數(shù)的相同性能指標(biāo)的研究基本沒(méi)有。為此,本研究中建立了相同體積的四極和六極電磁發(fā)射器,利用Simplorer軟件建立了脈沖成型網(wǎng)絡(luò),對(duì)比分析了兩者電樞的電流密度及磁場(chǎng)屏蔽效果,為電磁發(fā)射器的選取及優(yōu)化提供了一定的參考。

1 物理模型及仿真條件與方法

1.1 四極和六極電磁發(fā)射器模型

電磁發(fā)射器電樞模型如圖1所示,為了更好地對(duì)比四六極電磁發(fā)射器的電磁特性,采用控制變量法,控制兩者電樞的口徑、拋體裝載區(qū)、電樞臂以及體積大小相等,其中圖2為電樞基本結(jié)構(gòu)參數(shù)示意圖。

圖1 電磁發(fā)射器電樞模型

圖2 電樞基本結(jié)構(gòu)參數(shù)示意圖

其中,h表示電樞頭部厚度,d0表示尾翼厚度,ds表示尾翼長(zhǎng)度,dt為樞軌接觸面長(zhǎng)度,r為電樞頭部與尾翼連接處的圓角半徑。

以六極電樞為例,裝配后的六極電磁發(fā)射器三維模型如圖3所示,6根軌道等距離、軸對(duì)稱分布,相鄰兩根軌道對(duì)應(yīng)的圓心角都為60°,同時(shí)任意2根相鄰軌道的電流方向相反,箭頭方向?yàn)殡娏鞣较?6根導(dǎo)軌中的電流都會(huì)在空間產(chǎn)生磁場(chǎng),在電樞中心位置磁場(chǎng)會(huì)相互抵消,保證拋體免受強(qiáng)磁場(chǎng)干擾,在電樞中心以外的其他位置磁場(chǎng)得到加強(qiáng)。

圖3 六極電磁發(fā)射器三維模型

六極電磁發(fā)射器的參數(shù)如表1所示,其中a、b、L、d分別為導(dǎo)軌的截面寬度,截面長(zhǎng)度,長(zhǎng)度以及電樞口徑如表1所示。

1.2 仿真條件與方法

電磁發(fā)射裝置的發(fā)射過(guò)程是極為復(fù)雜的瞬態(tài)過(guò)程,整個(gè)過(guò)程持續(xù)時(shí)間極短,電樞需要在極短的時(shí)間內(nèi)達(dá)到超高的速度,常規(guī)電源難以滿足要求,一般使用高功率脈沖電源,因此仿真激勵(lì)采用外電路電容放電方式。

仿真的整個(gè)過(guò)程考慮鄰近效應(yīng)、趨膚效應(yīng),對(duì)電磁發(fā)射器模型內(nèi)部進(jìn)行網(wǎng)格劃分,設(shè)定網(wǎng)格單元的最大邊長(zhǎng)為模型相應(yīng)邊長(zhǎng)的1/50;為了保證仿真的準(zhǔn)確度,求解域設(shè)定為400%。

采用有限元法對(duì)四極和六極電磁發(fā)射裝置進(jìn)行仿真,對(duì)兩者電樞的推進(jìn)力、電流密度和磁場(chǎng)分布情況進(jìn)行分析比較,此方法可以實(shí)現(xiàn)對(duì)仿真的高精度的求解。

2 理論計(jì)算基礎(chǔ)

2.1 脈沖成形網(wǎng)絡(luò)

脈沖成形網(wǎng)絡(luò)(PFN)的主要功能是對(duì)初極電能進(jìn)行壓縮,調(diào)節(jié)輸出脈沖波形向負(fù)載傳輸能量。根據(jù)PFN中晶閘管所處位置的不同,可分為Ⅰ型PFN和Ⅱ型PFN,圖4為兩型PFN的電路拓?fù)浣Y(jié)構(gòu)圖。

圖4 兩型PFN的電路拓?fù)浣Y(jié)構(gòu)

其中,C為儲(chǔ)能單元;T為放電開(kāi)關(guān);RC和LC為支路雜散電阻和電感;R0和L0為負(fù)載電阻和電感;RD和LD為續(xù)流支路電阻和電感;iC、iD、iO分別為電容支路電流、續(xù)流電流、負(fù)載支路電流。

仿真電路采用Ⅰ型脈沖成形網(wǎng)絡(luò),以Ⅰ型PFN為例,其放電分為2個(gè)階段:第一階段,脈沖電容對(duì)負(fù)載放電,續(xù)流支路截止;第二階段,續(xù)流支路開(kāi)始導(dǎo)通;圖5為第一階段等效電路圖。

圖5 第一階段等效電路圖

令R1=RC+R0,L1=LC+L0,電容初始電壓為U0,根據(jù)基爾霍夫電壓定律可得電路二階微分方程:

(1)

求解二階線性齊次微分方程為

(2)

(3)

在第一階段電容電壓隨時(shí)間逐漸減小,直至uab為0,此時(shí)續(xù)流支路開(kāi)始導(dǎo)通,續(xù)流支路電壓為

(4)

令uab(t1)=0,可得續(xù)流支路開(kāi)始導(dǎo)通時(shí)刻為

(5)

(6)

令uc(t)=0,可得電容電壓降為0的時(shí)刻為

(7)

2.2 軌道炮作用力定律

圖6為軌道炮的等效電路圖,a為電樞的起始位置,經(jīng)歷dt時(shí)間后的位置為b,電感增量為dL,電壓為V,電樞位移為dx,回路電流I(假設(shè)不隨時(shí)間和位移變化)。

圖6 軌道炮等效電路圖

假設(shè)作用在電樞上的力F,則其所做的機(jī)械功為

Wm=Fdx

(8)

軌道炮的感應(yīng)磁能增量為:

Wi=dLI2/2=L′dxI2/2

(9)

根據(jù)法拉第定律:

V=d(LI)/dt=L′dxI/dt=L′Iv

(10)

傳遞給電路的功:

Wg=VIdt=L′I2vdt=L′I2dx

(11)

由能量守恒定律:

Wg=Wm+Wi

(12)

由式(8)、式(9)、式(11)得到電樞力F為

(13)

3 仿真分析

3.1 電磁發(fā)射裝置總電流及推進(jìn)力分析

對(duì)電磁發(fā)射裝置的電路等效模型各元器件設(shè)置完參數(shù)以后即可進(jìn)行求解,圖7為六極發(fā)射裝置總電流隨時(shí)間變化曲線。

圖7 六極發(fā)射裝置總電流隨時(shí)間變化曲線

從圖7中可以看出回路電流在0～0.45 ms內(nèi)迅速增大,增大到270 kA左右之后電流迅速下降,直至為0;由2.1節(jié)可知,產(chǎn)生這種現(xiàn)象的原因是:在0～0.45 ms內(nèi)電路為脈沖電容對(duì)負(fù)載放電,電容電壓由初始電壓逐漸降低到0,發(fā)射裝置電流迅速增大到峰值,此時(shí)續(xù)流支路開(kāi)始導(dǎo)通,調(diào)波電感與二極管形成回路進(jìn)行放電,發(fā)射裝置電流逐漸減小為0。

電磁發(fā)射裝置在實(shí)際發(fā)射過(guò)程中一般持續(xù)時(shí)間極短,大約為2 ms左右,結(jié)合圖7可知,在t=2.0 ms時(shí),電流為18 kA,幾乎對(duì)電樞不再產(chǎn)生推力作用;因此假設(shè)電樞在t=0～2 ms過(guò)程中受到推力作用,在t=2 ms時(shí)電樞發(fā)射出去,四極和六極發(fā)射裝置的總電流隨時(shí)間變化趨勢(shì)相同,這是因?yàn)樵谙嗤碾娐返淖饔孟?四極和六極發(fā)射器自身的電阻與電感相差不大。

發(fā)射裝置電樞推進(jìn)力F的大小決定電樞的出口速度,平滑的推進(jìn)力能減小磨損及轉(zhuǎn)捩的發(fā)生,提高射擊的準(zhǔn)確度;在同一電路下,四極和六極發(fā)射器推力隨時(shí)間的變化情況如圖8所示,其中A、B、P分別代表四極電磁發(fā)射器、六極電磁發(fā)射器、四極電樞各項(xiàng)指標(biāo)比六極電樞各項(xiàng)指標(biāo)。從圖8中可以看出四極和六極電樞的推進(jìn)力變化趨勢(shì)相同,且與電流變化大體一致,都呈現(xiàn)先上升到達(dá)最高點(diǎn)再下降的趨勢(shì);比例系數(shù)P在發(fā)射過(guò)程中始終穩(wěn)定在數(shù)值2附近,這說(shuō)明在同一電路下四極電樞推力的大小是六極電樞推力的2倍左右。

控制六極電磁發(fā)射器每個(gè)軌道通入的電流與四極電磁發(fā)射器一致,得到其電磁力變化圖如圖9所示。從圖9中可以看出,六極電樞的推力明顯大于四極電樞,最大達(dá)到了14.1 kN,比例系數(shù)P穩(wěn)定在0.85附近。

圖9 相同電流電磁力變化圖

電樞的出口速度是衡量發(fā)射器作戰(zhàn)性能的重要指標(biāo)之一,圖10為電樞速度變化圖。

圖10 電樞速度變化圖

從圖10中可以看出,2種電樞的出口速度都呈現(xiàn)先迅速增大再緩慢上升的變化趨勢(shì),這與電樞所受電磁力有關(guān);最終四極電樞的出口速度可達(dá)1.67 km/s,而六極電樞僅為0.86 km/s;在發(fā)射全程中,四極電樞的速度都為六極電樞的2倍左右,由此可見(jiàn)在同一電路下,四極發(fā)射器的能量利用率高于六極發(fā)射器。

3.2 電磁發(fā)射裝置電流密度分析

電流密度分布是發(fā)射器重要的電磁性能指標(biāo)。發(fā)射器電流密度分布的好壞影響焦耳熱的大小,電流密度分布不均勻、產(chǎn)生集中現(xiàn)象會(huì)對(duì)軌道產(chǎn)生嚴(yán)重的燒蝕和熱損傷,進(jìn)而影響軌道炮的連續(xù)發(fā)射;由圖8可知,在t=0.45 ms時(shí),電路RLC放電結(jié)束,電路電流達(dá)到峰值,此時(shí)發(fā)射器電流密度達(dá)到最大;以t=0.45 ms時(shí)刻的電流密度為代表,圖11和圖12分別為四極和六極軌道電流密度分布圖。

圖11 四極軌道電流密度分布圖

圖12 六極軌道電流密度分布圖

從圖11、圖12中可以看出,對(duì)四極和六極軌道而言,由于電流趨膚效應(yīng)的作用,電流主要分布在軌道的外表面附近,軌道中間的電流分布較少;受電流鄰近效應(yīng)的影響,相鄰軌道的兩邊電流更加集中,即軌道內(nèi)側(cè)棱邊的電流密度大于外側(cè)棱邊的電流密度,最終在軌道的4個(gè)棱角處產(chǎn)生了電流的集中。

但是經(jīng)過(guò)對(duì)比發(fā)現(xiàn),四極軌道的電流密度比六極軌道的電流密度分布更加集中,六極軌道的電流密度為2.70×109A/m2,而四極可達(dá)4.48×109A/m2,是六極的1.5倍左右;這是因?yàn)樵诳傠娏鱅相同的情況下,六極單個(gè)軌道電流大小為1/3I,而四極為1/2I,是六極的1.5倍。

為更好地反映四極和六極電磁發(fā)射器對(duì)電樞電流分布的影響,選取電樞及樞軌接觸面仿真,結(jié)果如圖13和圖14所示。

圖14 六極發(fā)射器電樞電流密度分布

從圖中可以看出,受電流最短路徑的影響,四極和六極電樞的電流都集中分布在引流弧及樞軌接觸面的內(nèi)側(cè)棱邊上,且變化趨勢(shì)相同;從其中可以看出,四極電樞的電流密度可達(dá)到3.50×109A/m2,而六極電樞僅為3.38×109A/m2,說(shuō)明六極電樞整體的電流密度小比四極電樞小,但六極的分布更加均勻,可以減小磨損及轉(zhuǎn)捩的發(fā)生。

3.3 磁場(chǎng)數(shù)值屏蔽分析

在脈沖電流峰值時(shí)刻,發(fā)射器電樞前后端面的磁場(chǎng)強(qiáng)度云圖如圖15所示。

圖15 發(fā)射器電樞端面磁場(chǎng)強(qiáng)度

由圖15(a)、圖15(b)可知,2種電樞前端面的磁場(chǎng)強(qiáng)度均在引流弧處出現(xiàn)了磁場(chǎng)集中現(xiàn)象,但是六極電樞沒(méi)有四極電樞的明顯;四極電樞磁場(chǎng)強(qiáng)度的最大值為5.45 T,而六極電樞為4.05 T,且六極電樞中的磁場(chǎng)屏蔽區(qū)域明顯大于四極電樞,屏蔽效果更好。由圖15(c)、圖15(d)可知,四極電樞后端面磁場(chǎng)強(qiáng)度為9.90 T,而六極電樞為8.55 T,四極電樞在引流弧及電樞臂頭部產(chǎn)生了嚴(yán)重的磁場(chǎng)集中現(xiàn)象,而六極電樞的磁場(chǎng)分布更加均勻,有利于電樞更加精準(zhǔn)穩(wěn)定的發(fā)射。

為了更好地分析四極和六極電樞的磁場(chǎng)屏蔽效果,選取如圖16所示的path1、path2,分析其峰值時(shí)刻的磁場(chǎng)強(qiáng)度變化,如圖17所示。

圖16 電樞路徑圖

圖17 path1、2路徑磁場(chǎng)分布圖

由圖17(a)可知,在path1路徑上,2種電樞的磁場(chǎng)強(qiáng)度都為拋物線式變化,即由拋體裝載區(qū)外沿向內(nèi)逐漸減小,且為對(duì)稱分布,四極電樞磁場(chǎng)強(qiáng)度的最小值為1.75 mT,六極電樞僅為0.53 mT,幾乎為0;觀察比例系數(shù)P,P的大小均在1.46以上,在裝載區(qū)中心附近變化更加明顯,最大達(dá)到了12.40。由圖17(b)可以看出,在path2路徑上,四極電樞的磁場(chǎng)強(qiáng)度呈現(xiàn)波浪起伏式上升,變化特別明顯,而六極電樞始終穩(wěn)定在1 mT左右。因此可以看出,六極發(fā)射器的磁場(chǎng)屏蔽效果遠(yuǎn)遠(yuǎn)好于四極發(fā)射器,可有效減小對(duì)拋體裝載區(qū)中電子元器件中的干擾,有利于智能彈藥的精準(zhǔn)發(fā)射。

為進(jìn)一步探究電樞的磁場(chǎng)屏蔽情況,分別選取距離電樞中軸線1 mm和2 mm的兩點(diǎn),即圖16的point1和point2,觀察2種電樞發(fā)射過(guò)程的磁場(chǎng)變化,如圖18所示。在電磁發(fā)射器的發(fā)射過(guò)程中,2種電樞在考察點(diǎn)point1和point2的磁場(chǎng)強(qiáng)度變化趨勢(shì)相同,即先增大再減小,這與電流的變化趨勢(shì)相同,但是觀察圖7知,磁場(chǎng)的變化滯后于電流的變化;六極電樞在考察點(diǎn)point1的磁場(chǎng)強(qiáng)度最大為37 mT,而四極電樞為147 mT,六極電樞的磁場(chǎng)強(qiáng)度僅為四極電樞的25.2%,六極電樞在考察點(diǎn)point2的磁場(chǎng)強(qiáng)度最大為56 mT,而四極電樞為218 mT,六極電樞的磁場(chǎng)強(qiáng)度僅為四極電樞的25.6%;四極電樞考察點(diǎn)的磁場(chǎng)強(qiáng)度變化速度明顯大于六極電樞的變化速度,這極易在拋體裝載區(qū)產(chǎn)生感應(yīng)電流,影響發(fā)射裝置電子元器件的正常工作,降低發(fā)射準(zhǔn)確度。

4 結(jié)論

通過(guò)對(duì)六極電磁發(fā)射器進(jìn)行電磁仿真并與四極發(fā)射器進(jìn)行對(duì)比分析,可得以下結(jié)論:

1) 在同一激勵(lì)電路下,四極發(fā)射器的推進(jìn)力及出口速度是六極發(fā)射器的兩倍,具有較高的能量利用率。

2) 在通入單個(gè)軌道電流相同的情況下,六極電樞產(chǎn)生的推力大于四極電樞,大約是四極電樞的1.125倍左右。

3) 六極電磁發(fā)射器的結(jié)構(gòu)設(shè)計(jì)有效降低了樞軌接觸面及引流弧電流密度的集中現(xiàn)象,降低了軌道的熱損傷,延長(zhǎng)了軌道的使用壽命,有利于發(fā)射裝置的連續(xù)發(fā)射。

4) 六極電磁發(fā)射器能夠達(dá)到良好的磁場(chǎng)屏蔽效果,相較于四極電磁發(fā)射器而言,其磁場(chǎng)屏蔽范圍更大,可有效減小對(duì)拋體裝載區(qū)中電子元器件中的干擾,有利于智能彈藥的精準(zhǔn)發(fā)射。