劉文一,焦冀光

(91550部隊,遼寧大連 116000)

?

艦載環(huán)境對固體發(fā)動機裝藥影響分析

劉文一,焦冀光

(91550部隊,遼寧大連 116000)

為了研究艦艇巡航時運動以及所經(jīng)海區(qū)溫度載荷對艦載固體發(fā)動機裝藥的影響,建立了艦艇運動模型和所經(jīng)海區(qū)的溫度變化模型;采用有限元技術(shù)計算了裝藥在這些載荷下的應(yīng)力、變形和溫度分布。結(jié)果表明:裝藥在艦艇動載荷的作用下,出現(xiàn)了較為嚴(yán)重的應(yīng)力集中和變形,溫度載荷對裝藥的影響不大;但是在艦艇動載荷和溫度載荷的耦合作用下,疊加的集中應(yīng)力有可能超過人工脫粘的許用應(yīng)力而破壞脫粘層導(dǎo)致發(fā)動機故障。

裝藥;運動;溫度;應(yīng)力

0 引言

固體發(fā)動機是艦載導(dǎo)彈的核心系統(tǒng),導(dǎo)彈被豎直裝載至艦艇執(zhí)行巡航任務(wù)時,艦艇航行時的搖擺、振動、晝夜及艦艇所經(jīng)海區(qū)的溫度變化會對貼壁澆注固體發(fā)動機裝藥產(chǎn)生影響。為了研究這些載荷對裝藥的影響,文中建立了艦艇運動模型和所經(jīng)海區(qū)的溫度變化模型;采用有限元技術(shù)計算了裝藥在這些載荷下的應(yīng)力、變形和溫度分布情況。

1 裝藥模型

發(fā)動機裝藥為等截面六星孔型,基于藥柱的幾何對稱性,縱向共有3個對稱剖面,將發(fā)動機裝藥均分成完全對稱的6等份,每份夾角為60°,取其中對稱的一份建立三維有限元計算模型,共劃分10 272個單元。其有限元模型如圖1所示。

圖1 裝藥有限元模型

推進劑具有粘彈性特性,泊松比為定值,其性能參數(shù)如表1所示.

表1 裝藥性能參數(shù)

表中應(yīng)力松弛模量E(t)用Prony級數(shù)形式表示,其級數(shù)表達式中的系數(shù)和指數(shù)需要實驗得到的松弛模量來確定。該發(fā)動機裝藥的松弛模量E(t)的Prony級數(shù)形式[1]為:

E(t)=1.799+1.431e-2.5t/aT+2.053e-2.5t/aT+

3.04e-0.25t/aT+3.886e-0.025t/aT

各溫度下的溫度-時間轉(zhuǎn)換因子aT由實驗測定。

2 艦艇運動對裝藥的影響

艦艇在水中的運動是一種具有6個自由度的空間運動,對搭載導(dǎo)彈影響最大的是橫搖、縱搖和升沉運動。發(fā)動機裝藥由于具有密度高、質(zhì)量大的特性,在這種低頻高幅載荷的作用下,發(fā)動機裝藥可能會出現(xiàn)界面脫粘。

2.1 艦艇運動模型

2.1.1 艦艇橫搖運動模型

以艦艇穩(wěn)定高度點為坐標(biāo)原點,艦尾至艦艏方向為X軸,垂直艦尾至艦艏方向為Y軸建立坐標(biāo)系。根據(jù)艦艇在海浪中的受力,列出艦艇的動平衡方程,得艦艇在海浪中的單自由度線性橫搖方程[2]:

(1)

2.1.2 艦艇縱搖運動模型

艦艇縱搖與橫搖類似,單自由度縱搖方程為:

(2)

2.1.3 艦艇升沉運動模型

艦艇在海浪區(qū)航行時,其升沉運動主要來自一階波浪力和二階波浪力的作用?？梢圆捎肏irom近似公式來計算艦艇所受波浪力[3],在Δh(t)和ψ(t)較小的情況下,艦艇所受的波浪力為:

(3)

(4)

2.2 仿真計算及結(jié)果分析

根據(jù)2.1建立的數(shù)學(xué)模型,在有限元軟件中建立合適的載荷模型,并將之加載到有限元模型上。艦艇橫搖和縱搖的角加速度、角速度和角度均為實測值,而升沉運動則為利用公式的近似值,故加載升沉運動載荷時,需在前處理軟件中建立合適的載荷場,才能使升沉運動的計算模型與實際情況吻合。

計算得到了裝藥在這些運動載荷下的應(yīng)力和變形如圖2所示。

圖2 裝藥在橫搖、縱搖和升沉運動載荷下的云圖

從圖中可以看出,發(fā)動機裝藥在艦艇復(fù)合運動載荷的作用下,出現(xiàn)了應(yīng)力集中,并產(chǎn)生了一定的變形。應(yīng)力集中部位出現(xiàn)在藥柱前端人工脫粘部位,其最大值為9.69 MPa;最大變形出現(xiàn)在裝藥沿軸線中部,其最大值為1.04 mm,并且其分布范圍較廣;在裝藥沿軸線中部較大范圍內(nèi)均出現(xiàn)了較大的變形,其值在0.62～1.04 mm。

3 艦載環(huán)境溫度對裝藥的影響

遂行艦艇巡航任務(wù)的導(dǎo)彈長期立式貯存在艦艇發(fā)射筒內(nèi),所經(jīng)海區(qū)溫度變化不同,而且每天晝夜有溫差,致使藥柱內(nèi)部產(chǎn)生熱應(yīng)力。

3.1 艦載環(huán)境溫度載荷模型

裝藥溫度來自于發(fā)動機殼體所接觸的外界環(huán)境溫度,但是艦載環(huán)境溫度是隨機變化的。這主要是由于艦艇在航行中會經(jīng)過不同的海域,經(jīng)歷不同的環(huán)境溫度。并且每天天氣陰晴和晝夜的變化,都會導(dǎo)致艦載環(huán)境溫度變化。艦艇長期巡航環(huán)境溫度隨機變化的模擬,國內(nèi)外有很多方法[4],文中采用式(5)所示的方程來近似模擬長期巡航過程中環(huán)境溫度隨機變化的規(guī)律。

T=Tm-Tacos (ωat+φa)-Tdcos (ωdt+φd)

(5)

式中:Tm為年平均溫度;Ta為年均溫度變化幅值;ωa為年循環(huán)頻率;φa為年循環(huán)初始相位角;Td為日均溫度變化幅值;ωd為日循環(huán)頻率;φd為日循環(huán)初始相位角;t為時間(d)。

根據(jù)歷年實測數(shù)據(jù),建立了艦艇巡航所經(jīng)歷的兩個海域A與B的溫度載荷數(shù)學(xué)模型:

T=26-8.5cosωat-5cosωdt

(6)

T=32-10cosωat-8cosωdt

(7)

而裝藥的導(dǎo)熱微分方程為:

(8)

式中:qv、ρ、Cv、λ分別為內(nèi)熱源發(fā)射率、密度、質(zhì)量定容熱容和導(dǎo)熱系數(shù)。對于無內(nèi)熱源的第一類邊值問題,用時間差分Galerkin格式可得有限元剛度方程[5-6]:

(9)

式中:T1和T0分別為時間步長Δt前后瞬時的單元節(jié)點溫度列陣。

計算采用的空間等參數(shù)單元坐標(biāo)和溫度分別為:x=N(xi),y=N(yi),z=N(zi)。

3.2 仿真計算及結(jié)果分析

3.2.1 裝藥的溫度場分布

分別計算了裝藥在A、B兩個海區(qū)3d后的溫度分布云圖,如圖3和圖4所示。

圖3 裝藥在A海區(qū)3 d后溫度分布

圖4 裝藥在B海區(qū)3 d后溫度分布

從圖3和圖4可以看出,裝藥在A和B兩個海區(qū)內(nèi)溫度分布趨勢一致,只是數(shù)值不同;由于裝藥起始溫度均低于環(huán)境溫度,經(jīng)過相同時間的緩慢熱傳遞之后,兩者均是在靠近殼體的部分溫度較高,這是由于金屬發(fā)動機殼體直接暴露在環(huán)境溫度中,使靠近發(fā)動機殼體的裝藥升溫很快。由于裝藥的導(dǎo)熱率較低,靠近裝藥內(nèi)表面,特別是星根部位的溫度變化不大。

同時,對比圖3和圖4可以看出,裝藥在B海區(qū)時溫度變化率較大,這是由B海區(qū)的溫度變化幅度較大所致。

3.2.2 裝藥的熱應(yīng)力場分布

分別計算了裝藥在A、B兩個海區(qū)3d后溫度載荷下的熱應(yīng)力分布云圖,如圖5和圖6所示。

圖5 裝藥在A海區(qū)3 d后熱應(yīng)力

圖6 裝藥在B海區(qū)3 d后熱應(yīng)力

從圖5和圖6可以看出,裝藥在A、B兩個海區(qū)的交變溫度載荷下,均出現(xiàn)了應(yīng)力集中現(xiàn)象,應(yīng)力集中位置均在裝藥的前后端人工脫粘部位。裝藥在A海區(qū)時最大熱應(yīng)力為0.166MPa,裝藥在B海區(qū)時最大熱應(yīng)力為0.260MPa,裝藥在溫度變化幅度大的B海區(qū)出現(xiàn)的熱應(yīng)力要比溫度變化幅度小的A海區(qū)出現(xiàn)的熱應(yīng)力大。

4 結(jié)論

1)發(fā)動機裝藥在艦艇橫搖、縱搖和升沉復(fù)合運動載荷的作用下,出現(xiàn)了應(yīng)力集中和變形。應(yīng)力集中部位出現(xiàn)在藥柱前端的人工脫粘部位,其最大值為9.69MPa;

最大變形出現(xiàn)在裝藥沿軸線中部,其最大值為1.04 mm,并且其分布范圍較廣;在裝藥沿軸線中部較大范圍內(nèi),均出現(xiàn)了較大的變形,其值在0.62～1.04 mm。

2)發(fā)動機裝藥在經(jīng)過溫度變化不同的兩個海區(qū)時,由于環(huán)境溫度變化而導(dǎo)致發(fā)動機裝藥溫度變化,但在兩個海區(qū)內(nèi)發(fā)動機裝藥的溫度分布較為一致,但發(fā)動機裝藥在B海區(qū)時溫度變化率較大,這是由B海區(qū)的溫度變化幅度較大所致。

3)發(fā)動機裝藥在A、B兩個海區(qū)的溫度載荷下,均出現(xiàn)了應(yīng)力集中現(xiàn)象,應(yīng)力集中位置均在裝藥的前后端人工脫粘部位。裝藥在A海區(qū)時最大熱應(yīng)力為0.166 MPa,裝藥在B海區(qū)時最大熱應(yīng)力為0.260 MPa,裝藥在溫度變化幅度大的B海區(qū)出現(xiàn)的熱應(yīng)力要比溫度變化幅度小的A海區(qū)出現(xiàn)的熱應(yīng)力大,但這兩個應(yīng)力值較小,不會對裝藥產(chǎn)生嚴(yán)重影響。

4)綜合1)和3)分析,發(fā)動機應(yīng)力集中均出現(xiàn)在

藥柱前段的人工脫粘部位,并且在艦艇運動載荷和溫度載荷的耦合作用下,疊加的集中應(yīng)力有可能超過人工脫粘的許用應(yīng)力而破壞脫粘層導(dǎo)致發(fā)動機故障。

[1] 李金飛, 黃衛(wèi)東, 李瑞亮. 基于實測艦載環(huán)境溫度的固體發(fā)動機藥柱累積損傷分析 [J]. 四川兵工學(xué)報, 2012, 33(10): 7-9.

[2] 戴余良. 潛艇在隨機海浪中搖蕩運動的仿真研究 [J]. 計算機仿真, 2001, 18(5): 42-45.

[3] 劉勝,荊兆壽. 高海情下船舶減搖鰭控制系統(tǒng)的仿真 [J]. 船舶工程, 1995(2): 37-41.

[4] 周紅梅, 高吉力, 齊強, 等. 長期貯存條件下固體發(fā)動機藥柱應(yīng)力分析 [J]. 海軍航空工程學(xué)院學(xué)報, 2010, 25(1): 54-56.

[5] 朱智春, 蔡峨. 固體火箭發(fā)動機藥柱三維溫度場有限元分析 [J]. 推進技術(shù), 1997, 18(2): 21-26.

[6] 劉鴻雁, 黎國保. 艦載環(huán)境溫度載荷作用下固體發(fā)動機傳熱分析 [J]. 戰(zhàn)術(shù)導(dǎo)彈技術(shù), 2011(1): 44-48.

Analysis of Solid Rocket Motor Propellant in Marine Environment

LIU Wenyi,JIAO Jiguang

(No.91550 Unit, Liaoning Dalian 116000, China)

When a naval vessel cruises, its motion and temperature load of the sea would affect the missile’s solid rocket motor propellant. The naval vessel cruise’s kinematic model was established to research this effect. Finite-element method was used to analyze stress, deformation and temperature distribution of propellant under those loads. It proves that severe stress concentration and deformation occurred on the propellant in motion of the naval vessel, and the temperature load slightly affects it, but under coupling of the ship’s movement and temperature loads, superimposed stress concentration would be likely to exceed allowable artificial debonding stress, the artificial debonding layer would be damaged, thus, the missile’s solid rocket motor would fail.

propellant; motion; temperature; stress

2014-10-24

劉文一(1982-),男,甘肅甘谷人,工程師,碩士,研究方向:飛行器結(jié)構(gòu)分析及評定。

V435

A