艦載環(huán)境對(duì)固體發(fā)動(dòng)機(jī)裝藥影響分析

劉文一,焦冀光

(91550部隊(duì),遼寧大連 116000)

?

艦載環(huán)境對(duì)固體發(fā)動(dòng)機(jī)裝藥影響分析

劉文一,焦冀光

(91550部隊(duì),遼寧大連 116000)

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

裝藥;運(yùn)動(dòng);溫度;應(yīng)力

0 引言

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

1 裝藥模型

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

圖1 裝藥有限元模型

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

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

表中應(yīng)力松弛模量E(t)用Prony級(jí)數(shù)形式表示,其級(jí)數(shù)表達(dá)式中的系數(shù)和指數(shù)需要實(shí)驗(yàn)得到的松弛模量來(lái)確定。該發(fā)動(dòng)機(jī)裝藥的松弛模量E(t)的Prony級(jí)數(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

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

2 艦艇運(yùn)動(dòng)對(duì)裝藥的影響

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

2.1 艦艇運(yùn)動(dòng)模型

2.1.1 艦艇橫搖運(yùn)動(dòng)模型

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

(1)

2.1.2 艦艇縱搖運(yùn)動(dòng)模型

艦艇縱搖與橫搖類(lèi)似,單自由度縱搖方程為:

(2)

2.1.3 艦艇升沉運(yùn)動(dòng)模型

艦艇在海浪區(qū)航行時(shí),其升沉運(yùn)動(dòng)主要來(lái)自一階波浪力和二階波浪力的作用。可以采用Hirom近似公式來(lái)計(jì)算艦艇所受波浪力[3],在Δh(t)和ψ(t)較小的情況下,艦艇所受的波浪力為:

(3)

(4)

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

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

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

圖2 裝藥在橫搖、縱搖和升沉運(yùn)動(dòng)載荷下的云圖

從圖中可以看出,發(fā)動(dòng)機(jī)裝藥在艦艇復(fù)合運(yùn)動(dòng)載荷的作用下,出現(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)境溫度對(duì)裝藥的影響

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

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

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

根據(jù)歷年實(shí)測(cè)數(shù)據(jù),建立了艦艇巡航所經(jīng)歷的兩個(gè)海域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ù)。對(duì)于無(wú)內(nèi)熱源的第一類(lèi)邊值問(wèn)題,用時(shí)間差分Galerkin格式可得有限元?jiǎng)偠确匠蘙5-6]:

(9)

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

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

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

3.2.1 裝藥的溫度場(chǎng)分布

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

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

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

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

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

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

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

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

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

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

4 結(jié)論

1)發(fā)動(dòng)機(jī)裝藥在艦艇橫搖、縱搖和升沉復(fù)合運(yùn)動(dòng)載荷的作用下,出現(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ā)動(dòng)機(jī)裝藥在經(jīng)過(guò)溫度變化不同的兩個(gè)海區(qū)時(shí),由于環(huán)境溫度變化而導(dǎo)致發(fā)動(dòng)機(jī)裝藥溫度變化,但在兩個(gè)海區(qū)內(nèi)發(fā)動(dòng)機(jī)裝藥的溫度分布較為一致,但發(fā)動(dòng)機(jī)裝藥在B海區(qū)時(shí)溫度變化率較大,這是由B海區(qū)的溫度變化幅度較大所致。

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

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

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

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

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

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

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

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

[6] 劉鴻雁, 黎國(guó)保. 艦載環(huán)境溫度載荷作用下固體發(fā)動(dòng)機(jī)傳熱分析 [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)分析及評(píng)定。

V435

A