趙凌康,田振國(guó),靳利園

(燕山大學(xué)河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點(diǎn)實(shí)驗(yàn)室，河北 秦皇島 066004)

在電磁發(fā)射過(guò)程中，大功率的脈沖電流通入到軌道和電樞所形成的回路中，可以在極短的時(shí)間內(nèi)將電樞加速至很高的速度[1]。伴隨著電樞高速運(yùn)動(dòng)，軌道處于大電流、高溫、高壓和強(qiáng)磁場(chǎng)的工作環(huán)境中[2]。而由于電磁熱效應(yīng)引起的溫度變化會(huì)導(dǎo)致軌道燒蝕、變形，這不利于電磁發(fā)射軌道在大功率下的連續(xù)發(fā)射[3]。因此，了解發(fā)射過(guò)程中軌道溫度分布，以及軌道溫度隨電樞移動(dòng)的變化規(guī)律，對(duì)于軌道冷卻系統(tǒng)的設(shè)計(jì)有著重要的參考意義。

近年來(lái)，隨著電磁軌道發(fā)射技術(shù)的快速發(fā)展，其所涉及的各種問(wèn)題的研究都在不斷深入。文獻(xiàn)[4-6]關(guān)注了電磁軌道炮結(jié)構(gòu)的優(yōu)化及部件的設(shè)計(jì)，從推進(jìn)方式、結(jié)構(gòu)設(shè)計(jì)、和幾何參數(shù)等方面給出了優(yōu)化方案，文獻(xiàn)[7]分析了發(fā)射過(guò)程中炮膛內(nèi)的磁場(chǎng)分布規(guī)律，以及電樞距離與峰值磁通密度的關(guān)系，并給出了能夠有效計(jì)算軌道炮智能炮彈部位感應(yīng)電場(chǎng)和磁場(chǎng)的方法，文獻(xiàn)[8-10]關(guān)注了軌道材料性能的研究，分析了作為軌道復(fù)合層合金材料的電磁、熱及力學(xué)方面的參數(shù)特性，文獻(xiàn)[11-13]研究了軌道的強(qiáng)度、剛度、動(dòng)態(tài)響應(yīng)等問(wèn)題。尤其是電磁軌道發(fā)射時(shí)的溫度場(chǎng)問(wèn)題，得到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注，文獻(xiàn)[14]搭建了基于光纖光柵的溫度測(cè)量系統(tǒng)，測(cè)量了導(dǎo)軌多點(diǎn)溫度數(shù)據(jù)，分析導(dǎo)軌上熱量分布及傳導(dǎo)過(guò)程，文獻(xiàn)[15]利用有限元方法，對(duì)發(fā)射中軌道溫升這一瞬態(tài)過(guò)程進(jìn)行仿真計(jì)算，在電磁場(chǎng)與溫度場(chǎng)耦合的情況下，求解了軌道橫截面的電流與溫度分布，文獻(xiàn)[16]通過(guò)三維數(shù)值仿真研究了電磁軌道炮在脈沖波形電流下的電磁場(chǎng)特性，分析了電樞軌道及樞軌接觸面的溫度變化。從目前的研究成果來(lái)看，關(guān)于電磁發(fā)射軌道溫度場(chǎng)的研究主要是通過(guò)實(shí)驗(yàn)和數(shù)值仿真的方式進(jìn)行分析，理論計(jì)算結(jié)果還很少見。

本文通過(guò)建立發(fā)射軌道和電樞的熱傳導(dǎo)方程，將電流的焦耳熱看作內(nèi)熱源，應(yīng)用格林函數(shù)的方法得到溫度場(chǎng)的理論解，通過(guò)算例分析軌道和電樞的溫度分布規(guī)律，以及非傅里葉熱效應(yīng)對(duì)溫度的影響。

1 模型及研究方法

如圖1所示，電磁軌道炮由軌道、電樞、蓋板等部分組成。建立圖2所示的坐標(biāo)系，其中，d為軌道寬度，b為軌道間距離，在進(jìn)行溫度場(chǎng)分析時(shí)，將軌道沿x軸方向離散成等距的N段，其間距均為Δx，軌道沿x軸坐標(biāo)表達(dá)為(x1,x2…xi,…,xj…xN)，軌道總長(zhǎng)度為L(zhǎng)。電樞滑過(guò)每個(gè)坐標(biāo)時(shí)將電樞與對(duì)應(yīng)位置的軌道視為一個(gè)區(qū)塊(如圖2虛線所示)，整體視為多個(gè)區(qū)塊疊加而成。對(duì)于每個(gè)區(qū)塊，其溫度場(chǎng)應(yīng)由如下三部分組成：電樞加熱后對(duì)軌道的影響；電樞在區(qū)塊上作用時(shí)間內(nèi)的軌道自身電加熱；電樞離開區(qū)塊后軌道自身的電加熱。

圖1 電磁發(fā)射軌道炮結(jié)構(gòu)示意圖

因此，軌道的溫升與電樞移動(dòng)的位置相關(guān)。在tA時(shí)刻電樞沿軌道滑動(dòng)的位移按下式計(jì)算[17]

(1)

圖2 溫度場(chǎng)計(jì)算簡(jiǎn)圖

2 軌道溫度的計(jì)算

因?yàn)檐壍篮碗姌猩郎匦再|(zhì)不同，電樞從發(fā)射起就開始持續(xù)加熱，軌道則是在電樞經(jīng)過(guò)后才開始加熱，所以將樞、軌加熱時(shí)間分別表示為tA和tR。將電流發(fā)熱作為內(nèi)熱源，則軌道與電樞的熱傳導(dǎo)控制方程為[18]

(2)

其中：λ為熱導(dǎo)率；t0為對(duì)應(yīng)材料的弛豫時(shí)間；ρ為材料密度；Cp為比熱容；g(x)為電加熱熱源；J表示軌道電流密度，則電加熱作用下熱源功率為g(x)=J2/σR，考慮到軌道發(fā)射情況，則其熱源功率可表達(dá)為[19]

(3)

其中:U為加載電壓；a為電樞沿軌道長(zhǎng)度的一半；σA為電樞的電導(dǎo)率；σR為軌道電導(dǎo)率，h軌道厚度；L′為電樞滑過(guò)軌道的長(zhǎng)度；R0為外接電阻，軌道與電樞接觸的熱能傳遞方式視為熱傳導(dǎo)。

假定軌道初始溫度為f1(y)，無(wú)外部熱流通量，則其初始條件為

T|tR=0=f1(y)

(4)

(5)

考慮軌道與電樞之間進(jìn)行熱交換，即在y=0處軌道邊界受電樞熱流密度影響，而軌道與外界空氣或外包層為絕熱，則軌道的初始條件可表達(dá)為

(6)

T|y=d=0

(7)

式(6)中:q(tA)表示電樞對(duì)軌道邊界的熱流密度，其值由電樞的溫度對(duì)y求偏導(dǎo)數(shù)得到；tA表示電樞行進(jìn)的時(shí)間，即電樞自加熱時(shí)間。

式(2)為含熱源的非傅里葉導(dǎo)熱方程，若忽略熱源項(xiàng)，則得如下普適的齊次非傅里葉導(dǎo)熱方程

(8)

式(8)的可以通過(guò)分離變量法表示為：

T(y,t)=Y(y)·E(t)

(9)

將(9)代入(8)中可得如下兩個(gè)微分方程

Y″+β2Y=0

(10)

t0E″+E′+αβ2E=0

(11)

將式(10)、(11)的解代入到式(9)中，可得溫度場(chǎng)的表達(dá)式為

(12)

(13)

當(dāng)初始條件為常數(shù)時(shí)，有

(14)

(15)

(16)

則軌道溫度場(chǎng)的解為

cosβmycosβmy′f(y′)dy′

(17)

將式(12)與格林函數(shù)解的形式相對(duì)應(yīng)，可以得到本問(wèn)題的格林函數(shù)為

cosβmycosβmy′

(18)

軌道溫度場(chǎng)的解為

(19)

其中α=λ/ρCp表示熱擴(kuò)散率。解的第一項(xiàng)表示軌道初始條件的影響，第二項(xiàng)表示軌道自身電加熱的影響，第三項(xiàng)表示電樞在該區(qū)域上熱流密度的影響，體現(xiàn)為非齊次邊界條件。但該解在邊界y=b上不收斂，應(yīng)對(duì)最后一項(xiàng)做分部積分，用等價(jià)的封閉形式代替，令

(20)

對(duì)(20)進(jìn)行分部積分可得：

(21)

其中g(shù)2為電樞熱源功率，將式(20)、(21)代入到(19)中即可得出軌道溫度場(chǎng)的解。

3 電樞溫度的計(jì)算

在軌道溫度場(chǎng)的解析解中含有電樞對(duì)軌道邊界的熱流密度項(xiàng)，故需求電樞的溫度場(chǎng)，然后用電樞的溫度場(chǎng)對(duì)y求偏導(dǎo)可以得出電樞對(duì)軌道的熱流密度

(22)

電樞在電磁作用下的熱源功率為

(23)

應(yīng)用式(8)來(lái)求解電樞溫度場(chǎng)，假定軌道電樞溫度為f2(y)，其初始條件為

T|tA=0=f2(y)

(24)

(25)

若不考慮電樞與外界的熱交換，則電樞邊界條件為

(26)

T|y=-b=0

(27)

電樞溫度場(chǎng)的解為

(28)

將式(28)代入到式(22)中，即可得到軌道邊界熱流密度，進(jìn)而得到軌道溫度場(chǎng)的完整表達(dá)式。

為了對(duì)比非傅里葉和傅里葉效應(yīng)計(jì)算結(jié)果的區(qū)別，按照同樣的條件給出軌道在僅考慮傅里葉效應(yīng)時(shí)的溫度場(chǎng)計(jì)算結(jié)果，傅里葉傳熱的控制方程如下

(29)

邊界條件同式(6)、(7)，(24)、(25)。用格林函數(shù)解法可以得到如下解析解

(30)

4 算例分析

取軌道總長(zhǎng)L=3 m，軌道上每段區(qū)域長(zhǎng)Δx=0.02 m，電樞的尺寸為0.2 m×0.2 m×0.2 m，導(dǎo)軌間距b=0.2 m，導(dǎo)軌厚度h=0.2 m。電樞材質(zhì)為鋁，電導(dǎo)率為σA=3.6×107(Ω·m)-1，密度ρA=2 700 kg/m3,傳熱系數(shù)λA=238 W/(m·℃)。軌道材質(zhì)為銅，電導(dǎo)率σR=5.88×107(Ω·m)-1，密度ρR=8 900 kg/m3，比熱容c=386 J/(kg·K)。傳熱系數(shù)λR=377 W/(m·℃)。軌道電感梯度Li′=0.451 24 μH/m，電樞初始溫度設(shè)為30 ℃，軌道初始溫度設(shè)為20 ℃，材料的弛豫時(shí)間t0=1×10-5s，系統(tǒng)加載電壓U=8 000 V，不計(jì)摩擦。

圖3表示用非傅里葉和傅里葉的計(jì)算公式，分別計(jì)算得到的電樞內(nèi)y=0-m處和y=-0.05 m處的溫度隨電樞在軌道內(nèi)滑動(dòng)時(shí)的變化曲線，可以看出，在同一時(shí)刻越靠近邊界即坐標(biāo)軸，溫度越高。并且，對(duì)于電樞溫度，非傅里葉和傅里葉計(jì)算的結(jié)果差距并不明顯。

圖3 l=0- m與l=-0.05 m處電樞溫度隨時(shí)間變化曲線

圖4為電樞到達(dá)xi=1 m時(shí)，對(duì)應(yīng)位置的軌道沿寬度方向的溫度分布曲線，可以看出電樞在同一時(shí)刻同一位置時(shí)，按照非傅里葉效應(yīng)計(jì)算軌道的整體溫度場(chǎng)結(jié)果比傅里葉計(jì)算的結(jié)果值略大，這是因?yàn)榭紤]到材料的弛豫效應(yīng)，造成溫度傳遞的延遲，體現(xiàn)為非傅里葉熱傳導(dǎo)時(shí)軌道的局部溫度相對(duì)較高。

圖4 電樞在xi=1 m處軌道溫度沿y分布曲線

圖5為電樞從xi=1 m分別移動(dòng)到xi=1.2 m，xi=1.4 m時(shí)，軌道上xi=1 m處的溫度分布變化曲線?？梢钥闯?，在電樞離開后，原位置繼續(xù)電加熱會(huì)引起該區(qū)塊內(nèi)溫度繼續(xù)升高。

圖5 電樞在xi=1 m處繼續(xù)移動(dòng)軌道溫度沿y分布曲線

圖6為電樞分別移動(dòng)到xi=1.2 m，xi=1.4 m時(shí)，軌道上xi=1 m處軌道電加熱的溫度分布曲線，可以看出，對(duì)于非傅里葉和傅里葉兩種算法，電加熱影響項(xiàng)計(jì)算的結(jié)果幾乎相等。造成這種現(xiàn)象的原因是，電樞極快滑過(guò)這兩個(gè)位置，在一個(gè)對(duì)應(yīng)軌道區(qū)塊內(nèi)，對(duì)于軌道整體溫度場(chǎng)加熱時(shí)間極短，邊界y=0+m處熱擾動(dòng)還未來(lái)得及傳遞。

圖6 電樞在xi=1.2 m和xi=1.4 m軌道自身電加熱影響曲線

圖7為當(dāng)電樞移動(dòng)到xi=2 m時(shí)，對(duì)應(yīng)位置的軌道溫度場(chǎng)分布曲線。與圖4相比，在軌道內(nèi)側(cè)y=0+m位置上按照非傅里葉效應(yīng)計(jì)算軌道的溫度場(chǎng)結(jié)果比傅里葉計(jì)算的結(jié)果差值明顯增加，而遠(yuǎn)離邊界位置的計(jì)算結(jié)果趨于一致。造成這個(gè)結(jié)果的原因是，在軌道末端，電樞運(yùn)動(dòng)速度極快，對(duì)應(yīng)的這個(gè)區(qū)塊的軌道自身電加熱時(shí)間極短，軌道內(nèi)部溫度較低，軌道外表面為絕熱條件，與電樞接觸表面受到電樞的快速、高溫差的熱沖擊。

圖7 電樞在xi=2 m處軌道溫度沿y分布圖

圖8為軌道與電樞接觸面上不同位置處的溫度隨時(shí)間變化曲線?？梢钥闯鲭姌袕某跏嘉恢瞄_始移動(dòng)到軌道終點(diǎn)的過(guò)程中，引起xi處軌道內(nèi)表面y=0+m位置的溫度變化。其中，線1、2分別表示在xi=0.5 m處非傅里葉效應(yīng)和傅里葉效應(yīng)計(jì)算對(duì)比的結(jié)果，其含義是電樞在未到達(dá)指定區(qū)域時(shí)，軌道一直保持初始溫度，電樞到達(dá)后邊界溫度瞬間有一個(gè)階躍，可以看出，兩種計(jì)算方法在電樞到達(dá)該區(qū)塊瞬間，軌道內(nèi)表面溫度存在階躍差值。線3、4、5、6、7、8分別表示在xi=1 m，xi=2 m及xi=2.5 m位置時(shí)兩種計(jì)算方法的溫度變化，可以看出，隨著電樞運(yùn)動(dòng)時(shí)間的增加，兩種方法計(jì)算的軌道內(nèi)表面溫度的階躍差值不斷增加。造成這個(gè)結(jié)果的原因是，隨著電樞移動(dòng)的距離增加，電樞表面y=0-m的溫度不斷增加，而電樞還未滑過(guò)的部分軌道還未形成通路，始終保持初始溫度，當(dāng)電樞溫度與軌道溫度差越大時(shí)，對(duì)軌道內(nèi)表面引起的溫度階躍差值越大。

圖8 軌道內(nèi)表面y=0+m處溫度隨時(shí)間變化曲線

圖9為電樞運(yùn)動(dòng)過(guò)程中對(duì)軌道xi=1 m和xi=2 m兩個(gè)位置由熱流密度產(chǎn)生的影響曲線，可以看出，對(duì)于各個(gè)相等的(x=0.02 m長(zhǎng)度的區(qū)塊，隨著電樞溫度的增加，軌道邊界的熱流密度同時(shí)增大，由非傅里葉效應(yīng)和傅里葉效應(yīng)計(jì)算的對(duì)溫度場(chǎng)影響的差值也增大。

圖9 電樞在xi=1 m和xi=2 m處熱流密度影響曲線

5 結(jié)論

1)靠近發(fā)射端的軌道加熱時(shí)間長(zhǎng)，電樞對(duì)溫度影響較小；而出口端的軌道加熱時(shí)間短，電樞溫度相對(duì)較高，對(duì)軌道的影響顯著。因而，靠近軌道發(fā)射端的溫度高，沿寬度分布均勻，靠近出口端的溫度值相對(duì)較低，沿寬度分布不均勻，靠近與電樞交界面處的相對(duì)較高。

2)由非傅里葉熱效應(yīng)計(jì)算的電磁炮發(fā)射時(shí)的溫度場(chǎng)與傅里葉定律計(jì)算的溫度場(chǎng)存在差值。電樞的非傅里葉熱效應(yīng)并不明顯，而軌道的非傅里葉熱效應(yīng)則較為顯著，尤其是在靠近出口一端，這是因?yàn)槌掷m(xù)加熱的電樞溫度較高，速度較快，與軌道進(jìn)行熱交換的時(shí)間短、溫差大。