某航天器火工裝置作動后殼體滯后裂紋機(jī)理研究

葉耀坤， 丁 鋒, 李曉剛， 穆慧娜

(1.北京空間飛行器總體設(shè)計部，北京 100094；2.北京理工大學(xué)，北京 100081)

0 引言

爆炸螺栓是航天器、衛(wèi)星系統(tǒng)廣泛應(yīng)用的火工爆炸分離裝置之一[1]。在某航天器的艙段分離系統(tǒng)中，使用了爆炸螺栓作為分離機(jī)構(gòu)的核心作動部件，其主要由電點(diǎn)火器、管座、活塞、密封圈、殼體等零部件組成，如圖1所示。

圖1 爆炸螺栓結(jié)構(gòu)組成圖Fig.1 Explosion bolt structure composition

上述爆炸螺栓的工作原理為：連接時，其類似于螺釘通過殼體前端螺紋將分離結(jié)構(gòu)連接在一起；分離時，依靠爆炸螺栓雷管裝藥的爆轟能量，推動活塞運(yùn)動撞擊殼體，并在爆炸螺栓殼體V形槽處將其撞斷，從而實現(xiàn)解鎖分離功能。

在上述爆炸螺栓研制試驗過程中，進(jìn)行了大、小藥量和高低溫條件下共計20發(fā)產(chǎn)品的驗證試驗。爆炸螺栓功能驗證均正常，試驗后檢查殘體并未發(fā)現(xiàn)殼體有異?，F(xiàn)象，1 h后觀察發(fā)現(xiàn)2發(fā)爆炸螺栓的殼體外圓柱面存在縱向裂紋，16 h后7發(fā)有裂紋，24 h后8發(fā)有裂紋，且部分產(chǎn)品的裂紋加大、增多。經(jīng)統(tǒng)計分析，本次試驗殼體出現(xiàn)滯后裂紋的產(chǎn)品比例為40%，均為高溫發(fā)火下的大藥量試驗件。殼體裂紋的位置均靠近活塞停止運(yùn)動部位，裂紋現(xiàn)象如圖2所示。

圖2 爆炸螺栓作動后的殼體滯后裂紋現(xiàn)象Fig.2 The hysteresis phenomenon of the shell after the explosion bolt actuated

上述現(xiàn)象表明爆炸螺栓存在一定缺陷，雖然均為完成解鎖分離功能后出現(xiàn)，但考慮爆炸螺栓作用后其內(nèi)部仍然存在高溫高壓燃?xì)?，如果殼體裂紋足夠大，可能會造成安全隱患，對航天器的安全構(gòu)型有不可估量的影響。為此，有必要針對該爆炸螺栓的殼體滯后裂紋機(jī)理開展研究，為后續(xù)的精準(zhǔn)優(yōu)化設(shè)計提供理論依據(jù)，減少試驗耗費(fèi)的成本。

為此，本文基于有限元模型的LS-DYNA數(shù)值仿真分析手段，分析爆炸螺栓密閉狀態(tài)下的瞬態(tài)爆轟現(xiàn)象，獲得產(chǎn)品作用過程中殼體或內(nèi)部結(jié)構(gòu)的應(yīng)力分布情況，為揭示爆炸螺栓作用的殼體滯后開裂機(jī)理、優(yōu)化螺栓殼體結(jié)構(gòu)設(shè)計提供理論支撐。

1 有限元模型建立

1.1 LS-DYNA建模

爆炸螺栓主要由點(diǎn)火器、雷管、活塞、墊圈、殼體等組成，爆炸工作過程為：電點(diǎn)火器通電升溫產(chǎn)生火焰，通過傳火孔引燃雷管并輸出爆轟波，活塞在爆轟波的沖擊作用下開始加速并對螺栓殼體形成機(jī)械沖擊，使殼體在預(yù)置V形槽處斷裂，實現(xiàn)產(chǎn)品的解鎖分離。將爆炸螺栓三維模型按照總—分—總的方式，先對總體模型做有限元模型簡化，再分塊建模，最后進(jìn)行模型集成。該方法有兩大優(yōu)點(diǎn)：1)通過全局考慮，對于非關(guān)鍵做功步驟所牽涉的結(jié)構(gòu)可以適當(dāng)簡化，然后進(jìn)行模型分拆，便可多人多任務(wù)同時進(jìn)行，有利于加快計算進(jìn)度；2)建模過程由易到難，有利于尋找和控制各做功過程中影響數(shù)值計算的因素。在本研究中，先按零部件進(jìn)行拆分，結(jié)合工作面切割方法，對各部分進(jìn)行六面體為主、四面體為輔的網(wǎng)格劃分，建立基本的有限元運(yùn)算模型，有限元網(wǎng)格的尺寸為0.02～0.04 cm，網(wǎng)格總數(shù)量超過10萬個，如圖3所示。

(a)全局模型

(b)1/4模型圖3 爆炸螺栓數(shù)值仿真有限元模型Fig.3 Numerical simulation finite element model of the explosion bolt

爆炸螺栓的關(guān)鍵部件殼體、活塞的有限元模型如圖4所示，通過多層次網(wǎng)格劃分后，所獲有限元單元主要以六面體為主，可有效提高仿真的效率和精度。

(a)殼體有限元模型

(b)活塞有限元模型圖4 螺栓關(guān)鍵零件的有限元模型Fig.4 The finite element model of the key parts of the bolt

1.2 算法選擇

本研究采用數(shù)值模擬軟件ANSYS/LS-DYNA進(jìn)行爆轟傳遞過程研究。使用該程序，用ANSYS建立模型，用LS-DYNA做顯式求解，然后用標(biāo)準(zhǔn)的LS-PREPOST后處理來觀看仿真分析結(jié)果。顯式有限元程序LS-DYNA和ANSYS程序強(qiáng)大的前后處理在該軟件中得到了有效結(jié)合。用LS-DYNA的顯式算法能快速求解瞬時大變形動力學(xué)、大變形和多重非線性準(zhǔn)靜態(tài)問題以及復(fù)雜的接觸碰撞問題，也可以在ANSYS和ANSYS-LS-DYNA之間傳遞幾何信息和結(jié)果信息，便于執(zhí)行連續(xù)的隱式-顯式/顯式-隱式分析。

LS-DYNA程序算法以拉格朗日(Lagrange)算法為主，兼有歐拉(Euler)和ALE(Abritrary Lagrangian Eulerian)算法。在有限元分析中，需根據(jù)模型的具體結(jié)構(gòu)、應(yīng)力應(yīng)變產(chǎn)生的主要位置、材料的類型以及計算的精度要求進(jìn)行算法的選擇。所選算法的合理性將直接影響計算結(jié)果的準(zhǔn)確度與可信度，同時還可能大幅度地改變運(yùn)算周期。由于爆炸螺栓包含火工品單元，若采用Lagrange單元，其裝藥結(jié)構(gòu)在數(shù)值模擬的過程中會產(chǎn)生大的形變，畸變的網(wǎng)格可能會減慢運(yùn)算速度，甚至?xí)?dǎo)致計算終止[2]。因此，本研究對空氣、填充物與雷管裝藥采用ALE多物質(zhì)算法，殼體等其他結(jié)構(gòu)采用Lagrange單元，并采用關(guān)鍵詞*CONSTRAINED_LAGRANGE_IN_ SOLID_TITLE完成ALE與Lagrange網(wǎng)格的耦合。

1.3 材料模型與狀態(tài)方程

LS-DYNA程序中有金屬和非金屬材料100余種可供選用，如彈性、彈塑性、超彈性、泡沫、玻璃、地質(zhì)、混凝土、土壤、復(fù)合材料、炸藥及起爆燃燒、剛性以及自定義材料[2]，在本研究中，結(jié)合不同零部件的特征選用不同的材料模型，便于兼顧仿真分析的精度與效率。

(1)主要金屬材料模型選取

F=RU

(1)

式中，R=FU-1為正交旋轉(zhuǎn)張量，U=(FTF)1/2為對稱正定張量，稱為右伸長向量。采用Green-Naghdi應(yīng)力率

(2)

式中，ω=RRT為旋轉(zhuǎn)率張量，σ為柯西應(yīng)力。

各向同性硬化材料塑性加載時，屈服面中心固定不變而屈服面半徑增大，隨動硬化材料塑性加載時，屈服面半徑不變而屈服中心沿塑性應(yīng)變方向平移。對于大多數(shù)實際材料，屈服硬化規(guī)律介于各向同性與隨動硬化之間，稱為混合硬化模式[2]。

由于爆炸螺栓的殼體、活塞、管座、點(diǎn)火器殼體是開裂現(xiàn)象的關(guān)鍵零部件，其材料均為30CrMnSiNi2A，并均進(jìn)行了除氫處理。為了保證仿真分析的精度，其材料模型選用彈塑性混合硬化模式。有限元計算采用cm-g-μs單位制，導(dǎo)出單位：壓力107N，應(yīng)力/壓力Mbar，速度10 km/s，密度103kg/m3。仿真分析過程中主要金屬零部件的材料參數(shù)見表1。

表1 主要金屬零部件的材料參數(shù)Tab.1 Material parameters of main metal parts

(2)流固耦合模型中的空氣域材料模型

在爆炸螺栓仿真分析模型中，建立了流固耦合模型，其中空氣域選用空白材料模型,可避免炸藥爆轟過程計算應(yīng)力、應(yīng)變[3]。

(3)爆炸螺栓裝藥材料模型與狀態(tài)方程

爆炸螺栓的點(diǎn)火器裝藥為鎂粉點(diǎn)火藥，雷管主裝藥為斯蒂芬酸鉛(150 mg)、疊氮化鉛(220 mg)、太安(310 mg)3種藥劑，其中太安作為主裝藥。通過火炸藥材料手冊可以獲得上述3種藥劑的爆速分別為3.8,4.7 ,7.2 km/s[4]。

在爆炸螺栓仿真分析模型中，所有藥劑均采用高能炸藥燃燒材料模型，輸入?yún)?shù)：爆速D，C-J壓力Pc-j與燃燒系數(shù)

Ff=max(Ff1、Ff2)

(3)

式中，te為爆轟波由起點(diǎn)傳至當(dāng)前單元中心處所需最短時間。若Ff>1，則取Ff=1。

此外，在爆炸螺栓裝藥序列的數(shù)值模型中，主發(fā)裝藥均采用*MAT_HIGH_EXPLOSIVE_BURN模型和JWL狀態(tài)方程，被激發(fā)藥采用LS-DYNA程序中為反應(yīng)速率提供的Lee-Tarver點(diǎn)火與增長反應(yīng)率模型。其中JWL狀態(tài)方程是凝聚炸藥典型的動力學(xué)狀態(tài)方程，它能描述火炸藥爆轟產(chǎn)物壓力—比容—比內(nèi)能(P-V-E)關(guān)系[5]。方程形式為

(4)

點(diǎn)火與增長狀態(tài)方程：此模型假設(shè)少量炸藥由于沖擊波加熱而被點(diǎn)火，形成微小球形燃燒區(qū)。其反應(yīng)率由壓力及表面積所控制。微小球形燃燒區(qū)的增長，使熱點(diǎn)間的炸藥在10-6s內(nèi)耗盡?；谏鲜黾僭O(shè)的反應(yīng)率方程可寫為

G1(1-λ)cλdpy+G2(1-λ)eλgpz

(5)

式(5)中右邊第一項代表點(diǎn)火項，參數(shù)a是臨界壓縮度，在一定的壓縮度情況下才允許點(diǎn)火。在點(diǎn)火項和第一增長項中的(1-λ)參數(shù)b和c等于2/3，代表向內(nèi)的球形顆粒燃燒。點(diǎn)火量作為沖擊強(qiáng)度和持續(xù)時間的函數(shù)，其大小由參數(shù)I和x控制。G1和d控制點(diǎn)火后的早期增長，G2和z確定高壓反應(yīng)速率[6]。各項參數(shù)取值如表2所示。

表2 點(diǎn)火—增長反應(yīng)速率方程主要參數(shù)Tab.2 Ignition-growth reaction rate equation main parameters

2 爆炸螺栓作用過程機(jī)理分析

爆炸螺栓的作用過程大致可以分為3個階段[7]：1)點(diǎn)火器通過管座的傳火孔起爆雷管；2)爆轟波在雷管裝藥內(nèi)部傳爆并加載于活塞端面；3)活塞在沖擊波作用下向前撞擊螺栓殼體，致使螺栓殼體在預(yù)置V形槽處斷開并實現(xiàn)分離。

由于爆炸螺栓流固耦合模型計算復(fù)雜，仿真分析過程中需采用先局部后總體的運(yùn)算模式，通過局部計算，調(diào)整模型網(wǎng)格的各項指標(biāo)，消除網(wǎng)格畸變、死循環(huán)等困擾，然后進(jìn)行系統(tǒng)模型的總體計算，實現(xiàn)爆炸螺栓做功全過程的仿真分析。

首先，電點(diǎn)火器產(chǎn)生的火焰通過點(diǎn)火孔引燃雷管，使得雷管內(nèi)的裝藥產(chǎn)生強(qiáng)爆轟輸出。相對而言，火焰的沖擊強(qiáng)度要小得多，且通過點(diǎn)火孔衰減，對雷管處燃腔產(chǎn)生的沖擊作用很小，因此，可弱化電點(diǎn)火器的功能作用，對其結(jié)構(gòu)進(jìn)行適當(dāng)簡化，采用等能量的高能炸藥替代，利用點(diǎn)火—增長模型進(jìn)行該過程的數(shù)值模擬[8]。點(diǎn)火器引燃后的壓力云圖分析結(jié)果如圖5所示。

圖5 點(diǎn)火器作用后的壓力云圖Fig.5 Pressure cloud map after the igniter acts

點(diǎn)火器引燃后，藥劑的化學(xué)反應(yīng)持續(xù)進(jìn)行并對外輸出能量，通過管座中的傳火孔起爆雷管組件中的斯蒂芬酸鉛，該歷程可由圖6進(jìn)行表征。

圖6 點(diǎn)火器傳火與起爆雷管第一層裝藥的壓力云圖Fig.6 Pressure cloud map of the first layer of charge of the igniter and detonator growth process

雷管的第一層裝藥被引爆后，爆轟波逐步成長并趨于穩(wěn)定后形成對下級裝藥的沖擊起爆，雷管成長過程的壓力云圖如圖7所示。

圖7 雷管第一層、第二層裝藥爆轟成長過程的壓力云圖Fig.7 Pressure cloud map of the first layer and second layer of charge initiation process

爆炸螺栓點(diǎn)火與雷管起爆過程還可以通過雷管裝藥的壓力進(jìn)行表征，在雷管的三級裝藥中，沿中軸線自上而下選取若干個有限元，繪制其壓力曲線，如圖8所示。

圖8 雷管第一層、第二層裝藥爆炸壓力增長過程曲線Fig.8 Pressure curve of the first layer and second layer of charge initiation process

上述結(jié)果表明,經(jīng)歷了約2 μs后，雷管靠近點(diǎn)火器的端面上開始出現(xiàn)壓力，說明點(diǎn)火器成功作用并開始起爆雷管，隨著時間推移，雷管裝藥上的壓力峰值不斷增大，至3.3 μs時達(dá)到最大，隨后壓力突然下降，表明第一級裝藥完全消耗。隨后，又出現(xiàn)一波逐步上升的壓力峰，表明下級裝藥被成功起爆并實現(xiàn)了能量的傳遞。各個有限元單元上的壓力值都具有相似的變化規(guī)律，先迅速上升然后快速下降，隨后出現(xiàn)一個震蕩的次級峰，其中第一個峰表示沖擊波快速通過時加載的壓力，后一個震蕩峰則是波后質(zhì)子在腔體內(nèi)形成的沖擊與壓力。

雷管中的疊氮化鉛被斯蒂芬酸鉛起爆后，輸出壓力繼續(xù)起爆主裝藥太安，輸出壓力在雷管殼輸出端蓋片、空氣等介質(zhì)中衰減后加載在活塞上，爆轟波傳遞過程如圖9所示。

圖9 主裝藥太安起爆過程壓力云圖Fig.9 Pressure cloud map of main charge Tai-an initiation process

點(diǎn)火器開始作用約3 μs后，雷管被引爆，約5 μs后雷管的三級裝藥均作用完畢，爆轟波到達(dá)活塞端面，活塞在沖擊波的作用下向前推進(jìn)，并高速撞擊螺栓殼體，殼體發(fā)生拉伸作用。48 μs左右，預(yù)制V形溝槽處開始出現(xiàn)頸縮，隨后開始出現(xiàn)環(huán)形裂紋，并于86 μs左右完全斷裂，爆炸螺栓被一分為二，實現(xiàn)了解鎖分離功能，如圖10所示。

圖10 爆炸螺栓殼體V形槽斷裂過程仿真結(jié)果Fig.10 Simulation result of shell-shaped fracture process of explosion bolt

上述爆炸螺栓的作用全過程仿真分析表明，爆炸螺栓的現(xiàn)有設(shè)計可以滿足解鎖分離功能要求，同時有效揭示了爆炸螺栓的內(nèi)部作用過程，可解決作用過程不可測難題，同時也為爆炸螺栓的殼體滯后開裂研究建立了可行的研究手段。

3 殼體滯后裂紋機(jī)理分析

采用第2章中的訓(xùn)練模型進(jìn)行數(shù)值仿真，獲得爆炸螺栓的殼體依次在2.599,5.599,15.99 μs等不同時刻的應(yīng)力云圖，如圖11所示。

圖11 爆炸螺栓的殼體應(yīng)力云圖Fig.11 Shell stress cloud diagram of explosion bolt

針對螺栓殼體受沖擊程度的不同，可將主要受壓區(qū)域按不同的壁厚分為a,b，c3段，其對應(yīng)的壁厚值依次為2.75,3.75,3 mm，此外a段與管座殼體重合，c段與活塞重合，b段部分處于空白域，如圖12所示。

圖12 爆炸螺栓的殼體壁厚分段Fig.12 Shell wall thickness section of the explosion bolt

從仿真分析結(jié)果中提取殼體不同段的有限元單元應(yīng)力曲線，如圖13所示。

圖13 爆炸螺栓殼體應(yīng)力曲線Fig.13 Stress curve of explosion bolt shell

圖13曲線表明，殼體腔內(nèi)的應(yīng)力峰值在3 000～4 000 MPa之間，應(yīng)力較大區(qū)域為b段空白區(qū)域，原因是雷管爆炸后爆轟波直接作用于殼體內(nèi)壁，使之往外膨脹，所受沖擊較強(qiáng)，因而產(chǎn)生的應(yīng)力也較大。

借鑒薄壁筒應(yīng)力計算方法，殼體理論可承受的應(yīng)力按以下公式進(jìn)行估算

σ=(πD2/4-πd2/4)σb

(6)

式中，D和d分別表示殼體圓筒的內(nèi)外徑，σb表示螺栓殼體熱處理后的強(qiáng)度，此處按照實際熱處理強(qiáng)度值取1 390 MPa。

依據(jù)式(6)可以計算出各段的拉伸應(yīng)力，結(jié)果如表3所示。

表3 爆炸螺栓殼體應(yīng)力計算結(jié)果Tab.3 Calculation results of shell stress of explosion bolts

從表中數(shù)據(jù)可以看出，除a段應(yīng)力2 972 MPa小于3 000 MPa以外，b，c兩段的拉伸應(yīng)力均大于或正處于仿真計算的結(jié)果范圍內(nèi)，因此b，c兩段通常不會率先開裂或失效。對于a段單從計算結(jié)果與仿真計算的結(jié)果來看，在強(qiáng)沖擊過后應(yīng)該會出現(xiàn)裂紋或失效，但螺栓殼體的a段處于管座的外圍，管座的壁厚約5 mm，所能承載的拉伸或剪切力為4 200 MPa左右，兩者相加為7 172 MPa，遠(yuǎn)遠(yuǎn)高于仿真結(jié)果應(yīng)力值，因而不會產(chǎn)生裂紋失效。c段由于活塞與殼體存在重合段，可以避免直接的爆炸沖擊損傷。理論上，爆炸螺栓殼體可以承受4 000 MPa應(yīng)力，最薄弱環(huán)節(jié)為殼體的b段，略小于實際應(yīng)力水平，存在強(qiáng)度臨界狀態(tài)。

然而，比較爆炸螺栓殼體的外徑與厚度，也可將其視為厚壁筒，對于受內(nèi)壓作用的厚壁筒，塑性區(qū)由內(nèi)壁開始向外擴(kuò)展，形式為內(nèi)層為塑性區(qū)，外層為彈性區(qū)。由于外層彈性區(qū)的約束，內(nèi)層塑性區(qū)的變形仍與彈性變形為同一量級，一旦全截面均進(jìn)入塑性狀態(tài)，無限制的塑性流動便成為可能。因此，如果爆炸螺栓在作用過程中未達(dá)到全截面塑性狀態(tài)，而內(nèi)部塑性區(qū)厚度較大，在存放過程中由于殼體容腔內(nèi)的爆炸壓力、溫度等環(huán)境條件逐步變化，引起殼體的內(nèi)部應(yīng)力水平緩慢發(fā)生波動變化，并隨著殼體溫度應(yīng)力的完全釋放以及內(nèi)部爆炸壓力的緩慢卸載，導(dǎo)致塑性區(qū)域的加大，從而出現(xiàn)裂紋，最終殼體應(yīng)力水平達(dá)到均勻平衡狀態(tài)。

此外，從材料疲勞裂紋的萌生與擴(kuò)展機(jī)理可知，金屬材料疲勞裂紋起源于應(yīng)變集中的局部顯微區(qū)域，即所謂疲勞源區(qū)，如圖14所示，爆炸螺栓在完成做功后，殼體的內(nèi)壁出現(xiàn)了幾個應(yīng)變集中局部顯微區(qū)域，主要分布在a,b,c段的交界處。

圖14 爆炸螺栓殼體的應(yīng)力應(yīng)變集中區(qū)域Fig.14 Stress and strain concentration area of explosion bolt shell

疲勞條紋的形成過程可以用Laird模型來說明，當(dāng)交變應(yīng)力為零時，循環(huán)開始時裂紋處于閉合狀態(tài)，拉應(yīng)力增大后，裂紋張開，且頂端沿最大切應(yīng)力方向產(chǎn)生滑移，當(dāng)應(yīng)力繼續(xù)增大至最大值，裂紋張開最大，相應(yīng)的塑性變形范圍也隨之?dāng)U大[9]。由于塑性變形，裂紋頂端鈍化，應(yīng)力集中減少，當(dāng)應(yīng)力反向時，滑移方向也改變，裂紋表面被壓攏。到壓應(yīng)力最大時，裂紋完全閉合，并恢復(fù)到這一周次的開始狀態(tài)，但裂紋卻擴(kuò)展了一個相當(dāng)裂紋擴(kuò)展速率數(shù)值的增量。爆炸螺栓在完成做功后，相當(dāng)長一段時間內(nèi)都處于如圖15所示的交變應(yīng)力震蕩狀態(tài)下，且隨著時間推移，殼體內(nèi)部爆熱、爆壓下降導(dǎo)致殼體發(fā)生縮變，其中爆熱下降持續(xù)約2 h后趨于平衡，爆壓下降持續(xù)約30 h后趨于平衡，在兩者綜合作用下，使得殼體b段的滯后開裂成為可能，這也與爆炸螺栓試驗過程b段出現(xiàn)滯后裂紋現(xiàn)象吻合。

圖15 螺栓殼體受力情況Fig.15 Stress condition of bolt shell

結(jié)合上述機(jī)理，為有效解決爆炸螺栓的殼體b段滯后開裂問題，一方面可以通過減少主裝藥量，降低爆炸產(chǎn)生的應(yīng)力水平；另一方面可以通過增加殼體b段的壁厚，提升殼體承受的極限應(yīng)力能力[10]。

4 結(jié)論

本文針對某航天器爆炸螺栓發(fā)火試驗后出現(xiàn)殼體滯后裂紋問題，采用了LS-DYNA有限元仿真手段，建立了爆炸螺栓作用全過程仿真模型，完整地模擬了爆炸螺栓從點(diǎn)火到輸出的全動態(tài)作用過程，得到了爆炸螺栓作用過程中各階段作動時序以及結(jié)構(gòu)應(yīng)力分布，揭示了殼體薄弱段的滯后裂紋機(jī)理，并提出了優(yōu)化改進(jìn)措施，為工程研制提供了理論支撐。