王叢飛，邢力超，許 光，滿 滿

(北京宇航系統(tǒng)工程研究所，北京 100076)

引言

在航天領(lǐng)域，破裂膜片廣泛應(yīng)用于運載、武器等液體型號的自生增壓管路中，達到隔離與打開的目的，破裂膜片能否正常打開直接關(guān)系到自生增壓系統(tǒng)是否正常工作，具有十分重要的作用[1-2]?，F(xiàn)役型號廣泛采用傳統(tǒng)鋁制平板形破裂膜片，純鋁材料耐腐蝕性差，同時，焊接易出現(xiàn)氣泡等缺陷。針對以上問題，擬采用新型不銹鋼破裂膜片對破裂膜片產(chǎn)品進行更新?lián)Q代，該膜片采用不銹鋼材料，產(chǎn)品的耐腐蝕性明顯改善，且具有較好的焊接性能，產(chǎn)品的示意圖如圖1所示。結(jié)構(gòu)為反拱形式，并在表面預(yù)制一定形狀的刻痕，在膜片打開的過程引導(dǎo)產(chǎn)品規(guī)則打開，可以有效防止多余物的產(chǎn)生。

不同于平板形和正拱形破裂膜片單純依靠強度破壞實現(xiàn)打開[3-5]，反拱形破裂膜片依靠結(jié)構(gòu)失穩(wěn)實現(xiàn)打開，影響結(jié)構(gòu)穩(wěn)定性的主要因素包括：原材料力學(xué)性能、膜片厚度、拱高、刻痕處剩余厚度等[6]。本研究采用有限元的方法對破裂膜片的打開壓力進行仿真分析研究，同時，通過試驗對產(chǎn)品的打開壓力性能進行了驗證。

1 產(chǎn)品幾何參數(shù)與性能指標

1.1 幾何參數(shù)

反拱帶槽形破裂膜片結(jié)構(gòu)參數(shù)圖如圖2所示?？毯畚挥诠皟?nèi)側(cè)，為“十”字形，減弱槽采用激光刻蝕加工，為矩形凹槽，刻痕尺寸規(guī)整，有利于進行參數(shù)測量。

圖1 反拱帶槽形破裂膜片結(jié)構(gòu)示意圖

圖2 膜片主要參數(shù)定義

膜片的主要參數(shù)：膜片外徑D1；起拱外徑D0；成形拱高H；膜片厚度t1；刻痕處剩余厚度t2；刻痕寬度B；成拱半徑R。

1.2 性能指標

火箭增壓管路用破裂膜片指標要求：通徑90 mm；正向打開壓力(0.4±0.05)MPa；反向承壓能力不小于0.5 MPa。

2 方案仿真分析

非線性穩(wěn)定性分析方法是結(jié)構(gòu)領(lǐng)域計算結(jié)構(gòu)最大承載能力的一種重要手段[7]。同時，羅代明[8]、閆照峰[9]均采用ANSYS有限元軟件非線性的方法對反拱形膜片的打開壓力進行過研究分析。

本研究針對反拱帶槽形破裂膜片采用ABAQUS軟件進行了穩(wěn)定性分析。實際結(jié)構(gòu)在工程上不可能是完善的，存在一定的缺陷，非線性穩(wěn)定性分析一般需引入初始缺陷，初始缺陷一般以屈曲模態(tài)為基礎(chǔ)，選取一定的缺陷比例因子。因此，首先對破裂膜片進行了模態(tài)分析，計算了破裂膜片產(chǎn)品的前幾階模態(tài)；在此基礎(chǔ)上，對破裂膜片進行非線性穩(wěn)定性分析，計算破裂膜片的最大承載能力。

2.1 有限元模型

通過經(jīng)驗分析，膜片厚度在0.2 mm左右，產(chǎn)品較薄，因此，選用了殼shell單元進行計算分析。模型為軸對稱結(jié)構(gòu)，為了簡化計算量，選取1/4模型進行分析，設(shè)置對稱邊界約束。破裂膜片的有限元模型如圖3所示。

圖3 破裂膜片有限元模型

其中，起拱外徑D0為90 mm，刻痕直徑D2為86 mm，膜片外徑D1為139 mm，該參數(shù)主要與上下?lián)跞懊芊饨Y(jié)構(gòu)相關(guān)，對破裂膜片打開壓力性能無影響，破裂膜片厚度、拱高、刻痕剩余厚度為設(shè)計參數(shù)。

2.2 材料屬性

破裂膜片材料選用316不銹鋼材料，彈性模量195 GPa，泊松比0.3，材料的塑性參數(shù)見表1所示。

表1 316不銹鋼材料塑性參數(shù)

2.3 分析步設(shè)置

第一步計算破裂膜片的屈曲模態(tài)，選擇lanczos分析步，提取破裂膜片的前幾階模態(tài)。膜片起拱的部位加載一定的壓力載荷，計算分析得到破裂膜片前三階屈曲模態(tài)。

第二步進行非線性穩(wěn)定性分析，選取Risk分析方法，邊界條件及載荷設(shè)置與buckling分析步設(shè)置相同，并通過編輯關(guān)鍵字設(shè)置缺陷因子，以第一階模態(tài)作為初始變形的參考依據(jù)。

缺陷因子對打開壓力有著重要影響，因此，正確設(shè)置缺陷因子對最終打開壓力有顯著影響[10]。羅代明[8]在計算反拱帶刀形破裂膜片的時候，通過對比分析數(shù)據(jù)與試驗數(shù)據(jù)，根據(jù)不同結(jié)構(gòu)選取的缺陷比例因子有0.07,0.08，且結(jié)構(gòu)不同選取的缺陷因子有一定區(qū)別。閆照峰[9]在計算反拱帶刀形破裂膜片的時候，選擇缺陷系數(shù)0.006，但未對缺陷因子對打開壓力的影響開展過研究。本研究的結(jié)構(gòu)為反拱帶槽形結(jié)構(gòu)，通過減弱槽實現(xiàn)破裂打開，與反拱帶刀形破裂膜片有一定區(qū)別，缺陷因子在選取上有一定的不同，通過研究缺陷因子對打開壓力的影響，并結(jié)合試驗情況，確定合適的缺陷因子。

2.4 仿真分析結(jié)果

首先對破裂膜片進行了屈曲分析，提取了產(chǎn)品的前三階屈曲模態(tài)，如圖4所示。將一階屈曲模態(tài)的變形位移乘以一定的缺陷因子作為初始缺陷進行非線性分析。

圖4 破裂膜片的前三階屈曲模態(tài)

通過非線性穩(wěn)定性分析，得到了破裂膜片在失穩(wěn)前結(jié)構(gòu)的變形分布，如圖5所示。該變形分布與一階模態(tài)的變形形式類似，說明結(jié)構(gòu)是沿著一階模態(tài)的形式進行擴展變形。

圖5 破裂膜片翻轉(zhuǎn)前的變形分布

通過變形分布發(fā)展，破裂膜片在壓力持續(xù)增加的過程中，在4個角部逐漸塌陷，在結(jié)構(gòu)的平衡轉(zhuǎn)換點，結(jié)構(gòu)突然失穩(wěn)，在失穩(wěn)塌陷的過程，結(jié)構(gòu)在刻槽位置破裂打開，實現(xiàn)翻轉(zhuǎn)爆破打開的目的。

分別研究了初始缺陷因子k、原材料厚度e1、膜片拱高h、刻痕剩余厚度e2等4個參數(shù)對破裂膜片打開壓力的影響。

圖6給出了缺陷因子對破裂膜片打開壓力的影響，結(jié)果顯示，在一定的范圍內(nèi)，缺陷因子越大，打開壓力越小，且缺陷因子對打開壓力分析結(jié)果有較為顯著的影響，在進行破裂膜片分析設(shè)計過程中，必須選取合適的缺陷因子。

圖6 缺陷因子對打開壓力的影響

圖7給出了反拱形破裂膜片拱高對打開壓力的影響，可以看出，在一定的范圍內(nèi)，打開壓力與拱高基本程線性關(guān)系，打開壓力隨起拱高度的增加而提高。

圖7 打開壓力隨拱高的變化

圖8給出了破裂膜片打開壓力隨著膜片厚度變化的趨勢，從數(shù)據(jù)可以得出，膜片打開壓力隨著厚度的增大而增大，之間的關(guān)系近似線性關(guān)系。

圖8 打開隨膜片厚度的變化趨勢

圖9給出了破裂膜片打開壓力隨刻痕剩余厚度的關(guān)系，結(jié)果顯示，在一定范圍內(nèi)，刻痕處剩余厚度越大，打開壓力越高。

圖9 打開壓力隨刻痕剩余厚度的關(guān)系

但刻痕剩余厚度0.06 mm處出現(xiàn)了不同趨勢，對該參數(shù)下破裂膜片的模態(tài)進行了分析研究，發(fā)現(xiàn)在該參數(shù)下的一階屈曲模態(tài)如圖10所示，與圖6失穩(wěn)前的變形分布有著顯著不同，說明該參數(shù)下，結(jié)構(gòu)的失穩(wěn)形式發(fā)生了變化，從而影響到最終的失穩(wěn)壓力。

通過以上分析，初步確定了破裂膜片的結(jié)構(gòu)參數(shù)，如表2所示。

圖10 刻痕剩余厚度0.06 mm時膜片的一階屈曲模態(tài)

表2 破裂膜片結(jié)構(gòu)參數(shù) mm

3 試驗研究

以仿真結(jié)果為基礎(chǔ)，開展了破裂膜片的加工及試驗驗證工作，為了保證產(chǎn)品的加工進度及一致性，采用激光刻蝕的方法。增壓速率對破裂膜片打開壓力有一定的影響，為了保證試驗的可靠性，選擇緩慢增壓的方式，選用自動化壓力試驗系統(tǒng)[11-13]。試驗原理如圖11所示，包括氣源、壓力控制及記錄設(shè)備、破裂膜片安裝法蘭等。

圖11 破裂膜片打開前后狀態(tài)

對該形破裂膜片開展了常溫狀態(tài)下的打開壓力試驗。破裂膜片在加壓過程，起拱側(cè)壓力逐漸增加，在增加到一定程度時，結(jié)構(gòu)發(fā)生屈曲失穩(wěn)，在失穩(wěn)的同時，膜片迅速翻轉(zhuǎn)，在翻轉(zhuǎn)過程，破裂膜片在刻痕位置扯開，形成泄漏通道，實現(xiàn)導(dǎo)通的目的。圖12給出了打開前后破裂膜片的狀態(tài)，表3給出了試驗結(jié)果。

圖12 破裂膜片試驗裝置示意圖

表3 破裂膜片打開壓力試驗結(jié)果 MPa

試驗結(jié)果顯示，90 mm通徑破裂膜片正向打開壓力在0.40 MPa左右，反向極限壓力均大于2.0 MPa，具有良好的反向承壓能力。

根據(jù)試驗結(jié)果，結(jié)合圖6缺陷因子對打開壓力的影響趨勢，認為該結(jié)構(gòu)在選取缺陷因子0.05時，仿真結(jié)果與試驗結(jié)果最為接近。根據(jù)分析及試驗情況，同類型反拱帶槽形破裂膜片在進行非線性穩(wěn)定性分析時，推薦缺陷因子0.05。

將缺陷因子在0.05時的分析結(jié)果與試驗結(jié)果進行了對比分析，對比數(shù)據(jù)見表4所示。

表4 仿真結(jié)果與試驗結(jié)果數(shù)據(jù)對比

以上對比數(shù)據(jù)顯示，仿真結(jié)果與試驗結(jié)果的誤差在5%左右，結(jié)果具有很高的可信度。非線性穩(wěn)定性的計算方法可以較好預(yù)示反拱形破裂膜片的打開壓力，可以作為產(chǎn)品設(shè)計、分析的依據(jù)。

4 結(jié)論

對反拱形破裂膜片開展了仿真分析工作，研究了缺陷因子、主要結(jié)構(gòu)參數(shù)對打開壓力的影響。同時開展了試驗驗證工作，主要結(jié)論如下：

(1) 采用非線性穩(wěn)定性分析的方法，可實現(xiàn)對反拱形破裂膜片打開壓力仿真分析，誤差不超過5%。對同類型、相似規(guī)格的破裂膜片產(chǎn)品在進行非線性溫度性分析時，推薦缺陷因子0.05;

(2) 破裂膜片的打開壓力主要與拱高、膜片厚度、刻痕剩余厚度相關(guān)，在一定范圍內(nèi)，隨著拱高、膜片厚度的增加，打開壓力線性增加，在一定范圍內(nèi)，刻痕剩余厚度越大，打開壓力越高，但剩余厚度超過該范圍時，會導(dǎo)致破裂膜片的模態(tài)發(fā)生變化，因此變形方式方法變化，將導(dǎo)致打開壓力的突變。