崔孟陽,王學(xué)智,丁日顯,童思遠(yuǎn),陳青榮

(1 空軍工程大學(xué)研究生院, 西安 710051； 2 空軍工程大學(xué)防空反導(dǎo)學(xué)院, 西安 710051)

0 引言

電磁軌道炮利用流經(jīng)軌道、電樞的大電流產(chǎn)生的強(qiáng)磁場,使電樞在洛倫茲力的作用下加速運(yùn)動到軌道炮出口,并拋射出去,是一種新概念動能武器。由于電樞運(yùn)動速度較快,為保證電樞和軌道始終具有較好的接觸,需要對軌道施加適當(dāng)?shù)念A(yù)緊力。在發(fā)射過程中,電樞受到電磁力并且發(fā)生變形,會對接觸的軌道面產(chǎn)生擠壓。兩側(cè)的軌道由于流過相反方向的強(qiáng)電流,也會產(chǎn)生較大的排斥力。因此,保證軌道處于較為穩(wěn)定的工作狀態(tài),并且具有較長的使用壽命,是電磁軌道發(fā)射技術(shù)的重要研究方向[1]。

對于在發(fā)射過程中軌道的變形問題,許多學(xué)者進(jìn)行了相關(guān)研究。文獻(xiàn)[2]中將電磁發(fā)射裝置的軌道簡化為移動載荷作用下彈性簡支梁來分析其受力與變形；文獻(xiàn)[3]中將方口徑電磁發(fā)射裝置的導(dǎo)軌及壁板簡化為雙層彈性基礎(chǔ)梁來分析導(dǎo)軌及壁板的動力響應(yīng)；文獻(xiàn)[4]中將簡化為彈性地基梁的導(dǎo)軌控制方程進(jìn)行Fourier變換,推導(dǎo)出發(fā)射過程中導(dǎo)軌的動態(tài)響應(yīng)方程,分析了導(dǎo)軌幾何參數(shù)、通入電流的強(qiáng)度及外包層的彈性特征等量對導(dǎo)軌動態(tài)響應(yīng)的影響。

文中利用Maxwell電磁分析軟件建立了電磁軌道炮的仿真模型,分析通入電流之后軌道炮的電磁場分布情況,得到體積力密度分布情況，并將數(shù)據(jù)導(dǎo)入Ansys Workbench的Static Structural模塊,進(jìn)行仿真得到電樞和軌道的變形情況和等效應(yīng)力分布情況。

1 軌道的受力分析模型

電磁軌道炮的結(jié)構(gòu)模型如圖1所示。電樞在兩軌道之間產(chǎn)生的強(qiáng)磁場作用下,受到較強(qiáng)的洛倫茲力作用使其沿著軌道移動。電樞在沿著軌道運(yùn)動的過程中,會對軌道造成擠壓。同時軌道還會受到相鄰軌道的排斥力,以及外包層對軌道產(chǎn)生的預(yù)緊力。軌道的受力情況如圖2所示。

圖1 軌道炮結(jié)構(gòu)模型圖

圖2 軌道的受力簡圖

在軌道的受力分析中,假設(shè)導(dǎo)軌的兩端沒有約束,外包層為軌道提供彈性支撐。軌道的動力學(xué)控制方程[7]為：

(1)

式中：w(x,t)為軌道的橫向變形;E為軌道材料的楊氏模量;Iy是軌道橫截面慣性矩;k為外包層支撐的彈性系數(shù);ρ是軌道的材料密度;A為軌道的橫截面積;q(x,t)為軌道受到的載荷。

對式(1)進(jìn)行關(guān)于x、t的傅里葉變換,得[7]:

(2)

其中,α、β是傅里葉變換中的積分變量。得：

(3)

對式(3)進(jìn)行傅里葉逆變換得：

(4)

即為在發(fā)射過程中軌道的橫向變形量。

表1 數(shù)值計(jì)算參數(shù)

2 軌道相關(guān)受力分析

軌道所受到的應(yīng)力主要有：電樞受電磁力發(fā)生變形對軌道造成擠壓應(yīng)力,相鄰軌道之間的排斥力以及電磁軌道炮外包層對軌道的預(yù)緊力。電樞與軌道的接觸應(yīng)力是相互的,分析計(jì)算電樞在運(yùn)動時所受電磁力在垂直于樞軌接觸面方向上的分力,就可以得到電樞對軌道的擠壓應(yīng)力。

2.1 軌道間的相互作用力

因?yàn)檐壍琅谠诎l(fā)射的過程中,存在明顯的電流趨膚效應(yīng)。因此作出假設(shè)：1)電流在軌道內(nèi)表面存在一定的趨膚深度；2)電流集中在趨膚深度的中心線位置。

實(shí)際應(yīng)用中電樞尾翼的電流段相對軌道上的電流段較短,因此忽略其對軌道的排斥作用。

圖3 軌道炮的原理圖

在t時刻,DE段的電流在GF段的內(nèi)點(diǎn)P(x1,-(hb+δr)/2)處生成的磁場強(qiáng)度[5]可通過以下式子計(jì)算：

(5)

式中：δr為軌道的電流趨膚深度；hb為兩側(cè)導(dǎo)軌的間距；x1為點(diǎn)P距軌道尾部端面的距離；s為電樞的位移。在軌道通入電流之后,點(diǎn)P處單位長度的軌道受到另一軌道的電磁排斥力為:

Fr=BDEI

(6)

2.2 電樞與軌道接觸應(yīng)力

電流經(jīng)軌道流向電樞,并形成閉合回路。根據(jù)通入電流的大小,可以計(jì)算出發(fā)射場內(nèi)任一點(diǎn)的磁感應(yīng)強(qiáng)度。對電樞所在位置的磁場強(qiáng)度進(jìn)行計(jì)算,并根據(jù)流經(jīng)電樞的電流大小計(jì)算出電樞受到的洛倫茲力。

將載流長方體軌道看成由許多載流直導(dǎo)線組成,當(dāng)軌道和電樞中通入電流時,左側(cè)軌道電流源點(diǎn)坐標(biāo)為m(x′,y′,z′),取電樞上中心點(diǎn)坐標(biāo)為q(x(t),0,0),則m點(diǎn)的電流元Idxi在q點(diǎn)產(chǎn)生的磁感應(yīng)強(qiáng)度[7]為:

(7)

對式(7)進(jìn)行積分,可得左側(cè)軌道在q點(diǎn)產(chǎn)生的磁場強(qiáng)度。

(-z′j+y′k)dx

(8)

同理可得右側(cè)軌道在q點(diǎn)產(chǎn)生的磁場強(qiáng)度為：

(9)

以上式中i、j、k分別為X、Y、Z軸方向上的單位向量；μ0為真空磁導(dǎo)率。

根據(jù)磁場疊加原理,可得軌道電流在q點(diǎn)產(chǎn)生的磁場強(qiáng)度為:

B=B1+B2

(10)

電樞在磁場中的受力為：

(11)

根據(jù)牛頓第三運(yùn)動定律,電樞與軌道的擠壓應(yīng)力是相互的。電樞對軌道的擠壓應(yīng)力[4]為：

Fc(t)=βF1

(12)

式中：β為電樞沿垂直于軌道接觸面方向上的電磁力分力轉(zhuǎn)化為擠壓應(yīng)力的轉(zhuǎn)換系數(shù)。

2.3 軌道所受應(yīng)力計(jì)算

軌道所受應(yīng)力的計(jì)算公式[10]為：

F(t)=Fc(t)+Fr(t)+F0

(13)

式中：F0是軌道炮外包層對軌道外部的機(jī)械預(yù)緊力。

3 仿真與分析

在Ansys Maxwell電磁分析軟件中建立電磁軌道炮的仿真模型,軌道材料為銅合金,電樞材料為鋁合金。軌道炮的其他參數(shù)如表2所示。

表2 軌道炮基本參數(shù)

3.1 建模仿真

在Maxwell中進(jìn)行電磁仿真,可以得到軌道和電樞上的電流分布情況、膛內(nèi)的磁場分布情況以及軌道炮各部位的體積力密度。其中體積力密度的矢量圖如圖4所示,可以看到電樞和軌道不同部位的受力情況：兩側(cè)軌道均受到向外的排斥力；電樞尾翼受到的電磁力也朝向與其接觸的軌道一側(cè),該電磁力會使電樞發(fā)生位移或變形擠壓軌道。

將在Maxwell建立好的電磁軌道炮模型和分析結(jié)果導(dǎo)入Ansys Workbench,建立電磁-應(yīng)力耦合仿真分析模型,在Static Structural模塊進(jìn)行應(yīng)力變形有限元分析。在只對軌道的兩端作固定約束的情況下,模擬軌道和電樞的變形情況。

圖4 體積力密度分布圖

當(dāng)電樞運(yùn)動距離分別為800 mm和1 700 mm時,電樞和軌道的等效應(yīng)力云圖和變形情況如圖5所示。隨著電樞的位移,膛內(nèi)的磁場會發(fā)生變化,軌道之間的電磁排斥力逐漸增大；軌道的外部兩側(cè)應(yīng)力較大,中間部位應(yīng)力較??；電樞的尾翼因?yàn)槭艿綇?qiáng)電磁力發(fā)生形變對軌道擠壓,是樞軌接觸壓力增大的重要原因。

3.2 樞軌接觸面的應(yīng)力分析

在進(jìn)行電磁-應(yīng)力耦合仿真時,對兩個軌道的外側(cè)施加固定約束,通過Sweep方法對軌道和電樞進(jìn)行網(wǎng)格劃分。當(dāng)電樞位移距離為1 100 mm時,電樞與軌道接觸面的應(yīng)力分布情況如圖6和圖7所示。由圖6可得：電樞表面的應(yīng)力分布由末端到肩部是逐漸增大的,電樞的肩部存在應(yīng)力集中,因此在實(shí)際應(yīng)用中也需要對電樞的肩部進(jìn)行特殊的處理。

圖5 軌道炮的等效應(yīng)力與變形圖

軌道上的應(yīng)力分布如圖7所示,與電樞接觸部位的應(yīng)力分布情況基本和電樞相同：較大應(yīng)力集中在與電樞肩部接觸的部位；軌道與電樞尾部的接觸面外側(cè)應(yīng)力要大于中間部位的；電樞運(yùn)動過的軌道部位應(yīng)力值較小,為軌道之間的排斥力；未通入電流的軌道段幾乎無應(yīng)力分布。

當(dāng)施加的激勵電流幅值分別為150 kA和300 kA時,軌道內(nèi)側(cè)表面中心線位置的應(yīng)力與位置關(guān)系如圖8所示。以通入300 kA電流為例,電樞運(yùn)動過的軌道部位排斥力約為1.49 MPa,軌道上所受的最大應(yīng)力約為4.46 MPa。電樞未滑過的軌道部位壓力逐漸減小,這與前面的理論分析保持一致。隨著電樞位置的移動,樞軌接觸面上的最大應(yīng)力分布如圖9所示,電樞運(yùn)動距離的增加使最大接觸應(yīng)力有所增大,但增幅不大；激勵電流幅值為300 kA時的最大接觸應(yīng)力約為150 kA的4倍。

圖6 電樞的等效應(yīng)力圖

圖8 軌道上的應(yīng)力與位置的關(guān)系

4 結(jié)論

利用Ansys Workbench有限元仿真軟件建立電磁-應(yīng)力多物理場耦合模型,通過Maxwell仿真得到體積力密度、電流和磁場分布并導(dǎo)入到Static Structural模塊進(jìn)行仿真分析,得到發(fā)射過程中軌道受到的應(yīng)力分布與變形情況。

圖9 電樞在不同位置接觸面上的最大應(yīng)力

1)電樞滑過區(qū)域的軌道,因?yàn)殡姌械臄D壓和軌道之間的排斥力,變形撓度較大；未滑過區(qū)域的軌道變形相對較小。

2)在發(fā)射過程中,電樞的尾翼受到較強(qiáng)電磁力會往外側(cè)張開,電樞和軌道的材料在滿足耐磨損的條件下,可以保證電樞和軌道有良好的接觸。

3)電樞對軌道的擠壓、軌道之間的排斥力都會導(dǎo)致軌道變形。

為了保證發(fā)射精度,需要控制軌道炮的口徑變化。增強(qiáng)軌道的抗彎剛度,減小軌道的撓度幅值可以改善軌道炮工作的穩(wěn)定性,有利于實(shí)現(xiàn)電磁軌道炮的安全連續(xù)發(fā)射。