王剛?cè)A，謝 龍，趙海龍，闞明先，肖 波，何 勇，宋盛義

（中國工程物理研究院流體物理研究所，四川 綿陽 621999）

現(xiàn)代先進兵器對彈丸速度的要求越來越高，彈丸超高速發(fā)射是增大射程、縮短飛行時間的重要前提，是未來戰(zhàn)場取得優(yōu)勢的保證。然而，傳統(tǒng)火炮受限于火藥氣體的膨脹聲速，發(fā)射速度很難超過2 km/s[1]。電磁發(fā)射技術(shù)在發(fā)射潛力上能滿足超高速發(fā)射的需求，在原理上不存在聲速限制，是各軍事強國重點關(guān)注的未來武器技術(shù)路線之一[2-6]。為了提升彈丸初速和發(fā)射效率，近年來，針對電樞與軌道的作用機理開展了大量研究，如電樞性能、電樞與導(dǎo)軌的接觸問題和匹配性能等[7-12]。串聯(lián)增強型軌道設(shè)計是在固有儲能條件下提高彈丸初速和發(fā)射效率的一種有效途徑。這種設(shè)計方案通過電路串聯(lián)電流提高了電樞上的磁場強度，從而提高了發(fā)射能力，后來Lü等將這一技術(shù)發(fā)展到了多匝串聯(lián)并列軌道炮[13-14]。

模擬分析電磁發(fā)射過程中電樞/軌道上的電磁場演化對于優(yōu)化和改進軌道、電樞設(shè)計極為重要，是控制軌道、電樞溫升、電樞轉(zhuǎn)捩的主要依據(jù)。很多學者通過靜磁場或時諧分析獲得了許多有價值的結(jié)論[15-17]，但是更好地模擬磁場演化需要瞬態(tài)分析。因為電磁軌道炮問題的三維有限元模擬是一個復(fù)雜的多場耦合的大型問題，它對于計算機資源（內(nèi)存、CPU 時間）的耗費很高。電磁軌道炮的發(fā)射過程中，軌道長度在米的量級，而電樞尺寸上的細微結(jié)構(gòu)往往不到1 mm，而且，電磁場在這些地方的分布非常重要，這就使得計算網(wǎng)格規(guī)模十分龐大，即使使用自適應(yīng)網(wǎng)格技術(shù)也很難完成全過程電磁場演化模擬。對于串聯(lián)增強型或多匝串聯(lián)并列軌道炮而言，構(gòu)型更復(fù)雜，計算量也更大。

為了更有效率地使用計算機資源，從而可以在有限資源情況下盡可能提高程序所能達到的最大解題規(guī)模，本文中將有限元方法和邊界元方法相結(jié)合[18-20]，從而免去對空氣域劃分網(wǎng)格的需要。此外，本文中所關(guān)心的主要物理現(xiàn)象僅發(fā)生在電樞-軌道接觸面附近，而在遠離軌道-電樞接觸面的地方，場量分布具有漸進平移不變的特性。利用該特點，可采用所謂“移動窗口法”[21]來減小實際計算域。其要點是：在不考慮末端效應(yīng)的情況下，假設(shè)軌道為無限長，并且，在每個時間步，僅將電樞-軌道接觸面附近的一小段軌道納入計算域，從而大幅度減少了未知量。采用該方法后，計算網(wǎng)格始終保持不變，因此邊界元方法中的迦遼金表面積分只需要計算一次。此后，即使電樞位置發(fā)生改變，也不需要再重新計算上述積分，這就大大減少了每個時間步的機時。實際計算表明，該方法有效提高了計算效率，通過該方法研制的Railgun3D程序可以模擬電磁軌道炮的全過程。在此基礎(chǔ)上，本文中將Railgun3D程序拓展到串聯(lián)增強型軌道的計算，并對其電磁場演化特征進行細致分析，觀察電流渦結(jié)構(gòu)演化和電流速度趨膚效應(yīng)等。

1 Railgun3 D程序簡介

Railgun3D程序是中國工程物理研究院流體物理研究所研制的一套電磁軌道炮數(shù)值模擬程序，采用了基于移動窗口法的有限元/邊界元耦合模擬技術(shù)，可以模擬電磁軌道炮發(fā)射過程中電磁場演化、電樞運動、溫度場演化等過程，能模擬發(fā)射中的焦耳加熱效應(yīng)、磁擴散效應(yīng)等。

Railgun3D程序中電磁場演化求解如下方程：

式中：m為電樞的質(zhì)量，Vs為彈丸體積，J為電流密度，μf和σyy分別為軌道/電樞接觸面S5上的滑動摩擦系數(shù)和法向應(yīng)力，V為電樞所包含的體積。

該程序與目前常用的商用軟件如Ansys、Comsol 等相比，穩(wěn)態(tài)和瞬態(tài)電磁場計算均具有同等精度[12]，瞬態(tài)分析計算在速度上具有明顯的優(yōu)勢。

2 串聯(lián)增強型軌道炮發(fā)射過程電磁場演化模擬分析

串聯(lián)增強型軌道炮的軌道和電樞設(shè)計示意圖如圖1所示，該設(shè)計中外軌為平面軌道，內(nèi)軌道中間使用了弧形軌道，發(fā)射過程中，電樞下表面采用弧形設(shè)計，安裝過程中通過預(yù)加應(yīng)力，與內(nèi)軌道弧面緊密接觸。為了便于分析發(fā)射過程中的電磁場演化過程，計算中使用了如圖2所示的理想梯形電流波形。由于模型的幾何對稱性，計算中使用通過上下對稱面和左右對稱面切割獲得的四分之一模型，以減小計算量。

圖1 串聯(lián)增強型軌道炮的軌道和電樞設(shè)計示意圖（四分之一模型）Fig.1 Schematic diagram of rail and armature design for a series enhanced railgun (1/4 model)

圖2 加載電流波形Fig.2 Loading current waveform

圖3～6給出了電流上升段中間時刻、剛到達電流最大值時刻、電流平臺中間時刻、電流下降段中間時刻等幾個典型時刻軌道/電樞上的磁場和電流密度分布。從圖中可以看出，內(nèi)軌道弧形部分與軌道邊緣處磁場和電流分布較強，與一般的單軌結(jié)構(gòu)不同，由于增強軌道的存在，驅(qū)動電流在增強軌道上產(chǎn)生了較大的磁場，該磁場遍布整條外軌道，由于電磁感應(yīng)，在內(nèi)軌道上會產(chǎn)生相應(yīng)的感應(yīng)電流，不過，在內(nèi)軌道的左端，其上的電流主要由流經(jīng)軌道的驅(qū)動電流決定，內(nèi)軌道的右端則主要是感應(yīng)電流，而對于一般的單軌軌道炮而言，其右端電流基本為零。增強軌道的引入使得軌道電樞上的磁場和電流分布變得更復(fù)雜，在進行電樞設(shè)計與分析時必須慎重考慮。

圖3 48 ns時刻磁場與電流密度分布Fig.3 Distributions of magnetic field and current density at 48 ns

圖4 108 ns時刻磁場與電流密度分布Fig.4 Distributions of magnetic field and current density at 108 ns

圖5 204 ns時刻磁場與電流密度分布Fig.5 Distributions of magnetic field and current density at 204 ns

對比圖3～6，內(nèi)軌道右端感應(yīng)電流在驅(qū)動電流上升段較大，在驅(qū)動電流平臺時間段內(nèi)，內(nèi)軌道右端感應(yīng)電流變小，在驅(qū)動電流下降段，該感應(yīng)電流再次變大，這一變化過程正好反映了計算獲得的感應(yīng)電流與該處磁通量變化相關(guān)，是比較合理的物理圖像。

圖7給出了電樞附近48、204 ns時刻的電流方向和電流密度分布，圖中帶箭頭線長度一致，只代表電流方向。從圖7可以清晰地看出前面所述的軌道上各部分的電流走向。在電樞上，該結(jié)構(gòu)設(shè)計的電樞，其頭部外側(cè)形成了電流渦結(jié)構(gòu)，并隨著電流增大，驅(qū)動電流從電樞后表面向頭部擴散，該電流渦結(jié)構(gòu)向頭部漂移。該渦結(jié)構(gòu)的形成與內(nèi)軌道上產(chǎn)生的感應(yīng)電流有關(guān)，是軌道上電流進一步在電樞上產(chǎn)生的感應(yīng)電流。不過總的看來，產(chǎn)生渦流處電流密度較小。

圖6 348 ns時刻磁場與電流密度分布Fig.6 Distributions of magnetic field and current density at 348 ns

圖7 兩個不同時刻的電流密度分布與電流方向Fig.7 Current density distributions and current directionsat two different moments

在電流下降段，負的電流變化率使得電樞后表面電流出現(xiàn)反向，緊挨電樞頭部處軌道電流密度較大，電樞中段側(cè)面出現(xiàn)電流渦結(jié)構(gòu)，如圖8所示。電樞上后表面電流方向的改變使得電樞受力情況變得復(fù)雜，如圖9所示，盡管電樞總體受力方向沒有改變，但后表面上洛倫茲力轉(zhuǎn)向炮膛中央，這有可能導(dǎo)致電樞與軌道接觸應(yīng)力不足，甚至出現(xiàn)電樞轉(zhuǎn)捩，而電樞轉(zhuǎn)捩可能是對軌道造成較大的破壞性損傷、炮管震動、電樞翻轉(zhuǎn)等的重要誘因，進而會影響軌道炮的重復(fù)發(fā)射性能、發(fā)射精度等。Wang 等[22]也觀察到了這種在電流下降段誘發(fā)的渦流現(xiàn)象，并用下坡轉(zhuǎn)捩理論進行了解釋。

圖8 348 ns 時刻電流密度分布與電流方向Fig.8 Current density distribution and current direction at 348 ns

圖9 348 ns時刻電流密度分布與洛倫茲力方向Fig.9 Current density distribution and Lorentz force direction at 348 ns

電樞加速過程中，軌道/電樞上的電流分布受擴散和速度趨膚效應(yīng)兩種機制控制。為了便于觀察，在軌道/電樞對稱面上作電流密度分布二維視圖，見圖10。在48 ns時，電樞運動速度低，電流擴散效應(yīng)顯著，是決定電流分布的主導(dǎo)因素。至108 ns，電樞速度有所提高，電流進一步擴散的同時，可以觀察到速度趨膚效應(yīng)的影響，該時刻電流密度整體增大是由于電流上升。至204 ns，電樞速度已經(jīng)較高，速度趨膚效應(yīng)比較明顯，成為電流相對大小分布的主導(dǎo)因素，但由于電流依然向軌道和電樞內(nèi)部擴散，使得該觀察窗口內(nèi)電流密度峰值有所下降。隨著電樞速度進一步提高，速度趨膚效應(yīng)愈加明顯。在電流下降段的348 ns時刻，電流峰值下降導(dǎo)致電流密度下降，軌道/電樞部分區(qū)域電流反向，只是電流分布更復(fù)雜。

圖10 4 個不同時刻電流密度的二維分布Fig.10 Two-dimensional distributions of current density at four different moments

3 結(jié) 論

利用流體物理研究所研發(fā)的Railgun3D電磁軌道炮數(shù)值模擬程序?qū)Υ?lián)增強型軌道炮進行了模擬，詳細分析了一復(fù)雜構(gòu)型的軌道/電樞在梯形驅(qū)動電流加載下電磁場的演化過程。由于增強軌道的存在，驅(qū)動電流在增強軌道上產(chǎn)生了較大的磁場，由于電磁感應(yīng)，在內(nèi)軌道上會產(chǎn)生相應(yīng)的感應(yīng)電流，即在炮口一端的軌道上有顯著的磁場和電流分布，感應(yīng)電流的大小依賴于驅(qū)動電流的變化率。

給出了多個時刻電樞附近的電流方向分布，觀察到了電流渦結(jié)構(gòu)的演化過程，并在電流下降段給出了電樞后表面上電流反向結(jié)果，指出該效應(yīng)可能是導(dǎo)致電樞與軌道接觸應(yīng)力不足、甚至出現(xiàn)電樞轉(zhuǎn)捩的重要因素。

中心對稱面上電流密度云圖顯示出磁擴散與速度趨膚效應(yīng)在整個過程中的競爭機制。