第一作者劉懷亮男，博士生，1980年生

通信作者閻紹澤男，博士，教授，1964年生

航天器中爆炸切割器的爆炸斷裂及沖擊響應(yīng)分析

劉懷亮，崔德林，閻紹澤

(清華大學(xué)機(jī)械工程系摩擦學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100084)

摘要：為了研究星箭分離沖擊環(huán)境特征和爆炸沖擊對(duì)衛(wèi)星的影響，以航天器中星箭包帶式連接結(jié)構(gòu)的某型爆炸切割器為研究對(duì)象，利用LS-DYNA顯式動(dòng)力學(xué)分析方法，建立了爆炸切割器爆炸切斷過(guò)程的三維有限元模型，分析了爆炸切割過(guò)程中剪切刀的速度和加速度，討論了連接桿上典型斷裂單元上的應(yīng)力應(yīng)變變化和斷裂特性，并分析了爆炸沖擊載荷的特征和不同階段切割器的結(jié)構(gòu)響應(yīng)特性。為航天器上爆炸切割器的爆炸斷裂過(guò)程提供了一種數(shù)值仿真的方法，同時(shí),也為星箭分離過(guò)程中產(chǎn)生的爆炸沖擊環(huán)境模擬和預(yù)示提供了數(shù)值參考。

關(guān)鍵詞：爆炸切割器；沖擊載荷；斷裂；動(dòng)力學(xué)；有限元分析

基金項(xiàng)目：國(guó)家自然科學(xué)基金面上項(xiàng)目(11072123);北京市自然科學(xué)基金資助項(xiàng)目(3132030)

收稿日期：2015-03-03修改稿收到日期：2015-04-01

中圖分類號(hào):O383.3； V421.7文獻(xiàn)標(biāo)志碼：A

基金項(xiàng)目：國(guó)家自然科學(xué)基金(512777186)

Dynamics explosion fracture and shock response analysis of explosive cutters in spacecraft

LIUHuai-liang,CUIDe-lin,YANShao-ze(State Key Laboratory of Tribology, Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China)

Abstract:In order to reveal the shock characteristics and the influence of explosion on satellite during the separation process，the explosive cutters of clamp band system in spacecraft were taken as the research object. A three-dimensional finite element model for analyzing the cut-off process of the explosive cutters was established by use of the explicit dynamic software LS-DYNA. The speed and acceleration of shear knife in the process of explosion fracture were investigated. Here the change of stress and strain and fracture features of typical fracture units on the link were discussed. And the characteristics of the exploding shock loading and structural response characteristics of the cutter were analyzed in different stages. A method of numerical simulation for the fracture process of explosive cutters in spacecraft was provided, which can also be regarded as a reference to the simulation and presage of pyrotechnic shock environment in the process of satellite-rocket separation.

Key words:explosive cutter; impact load; fracture; dynamics; finite element analysis

爆炸切割器是一種航天器上常見(jiàn)的火工分離裝置?；鸸し蛛x裝置因其具有體積小、結(jié)構(gòu)緊湊、能量密度高等優(yōu)點(diǎn)，在航天工程領(lǐng)域得到了廣泛的應(yīng)用。航天器級(jí)間分離、整流罩分離、星箭分離、太陽(yáng)能帆板展開(kāi)等動(dòng)作，都是依靠火工分離裝置來(lái)完成的?；鸸し蛛x裝置隨著航天技術(shù)的發(fā)展而衍生出多種類型，按作用的能量來(lái)源主要分為兩類[1]：點(diǎn)源裝置和線源裝置。點(diǎn)源裝置主要包括爆炸螺栓、分離螺母、螺栓拉拔器、爆炸切割器、電爆閥等，線源裝置主要包括柔性爆炸索、線性切割器等。Brauer[2]對(duì)各種火工裝置的作用原理、結(jié)構(gòu)組成、材料特性、作用形式等進(jìn)行了系統(tǒng)的總結(jié)，NASA-STD-7003A和美軍標(biāo)MIL-HDBK-83578 也對(duì)火工裝置進(jìn)行了詳細(xì)的闡述。Shmuel等[3-5]利用LS-DYNA3D對(duì)爆炸切割器、電爆閥等火工裝置進(jìn)行了仿真，并對(duì)此類裝置作用過(guò)程的影響參數(shù)進(jìn)行了分析。星箭包帶式連接結(jié)構(gòu)通過(guò)點(diǎn)源如爆炸螺栓連接，在星箭分離時(shí)包帶的應(yīng)變能釋放引起的沖擊環(huán)境又類似于線源，因此可以將包帶連接結(jié)構(gòu)看作點(diǎn)源和線源的混合類型[6]。關(guān)于包帶連接結(jié)構(gòu)解鎖分離引起的沖擊特性研究的文獻(xiàn)較少。譚雪峰等[7-8]對(duì)星箭包帶式連接結(jié)構(gòu)進(jìn)行了動(dòng)力學(xué)理論建模和計(jì)算，得到了包帶結(jié)構(gòu)沖擊響應(yīng)特性。秦朝燁等[9]研究了星箭包帶式連接結(jié)構(gòu)的整體耦合動(dòng)力學(xué)特性。武新峰等[10]利用LS-DYNA對(duì)星箭包帶式連接裝置的分離過(guò)程進(jìn)行了沖擊響應(yīng)分析。但是文獻(xiàn)[7-10]中，都將爆炸螺栓引起的沖擊簡(jiǎn)化為一個(gè)沖擊力和沖量，沒(méi)有爆炸螺栓引起的沖擊過(guò)程進(jìn)行詳細(xì)研究。初文華等[11]利用SPH-FEM方法和有限元方法建立了爆炸螺栓的三維數(shù)值模型，模擬分析了爆炸螺栓的解鎖分離過(guò)程及結(jié)構(gòu)在爆炸沖擊激勵(lì)下的動(dòng)響應(yīng)特性，得到螺栓體的加速度響應(yīng)時(shí)程曲線。

本文側(cè)重于研究爆炸切割器的爆炸動(dòng)態(tài)斷裂過(guò)程和沖擊響應(yīng)分析。以長(zhǎng)征三號(hào)甲型運(yùn)載火箭上包帶式連接結(jié)構(gòu)中常見(jiàn)的爆炸切割器[12]為例，利用LS-DYNA非線性有限元軟件，建立了爆炸切割器的三維有限元模型，對(duì)爆炸切割器的斷裂的整個(gè)作用過(guò)程進(jìn)行了數(shù)值模擬，得到了切割器近場(chǎng)環(huán)境下的爆炸沖擊載荷特征與結(jié)構(gòu)速度加速度響應(yīng)特征。本文的仿真結(jié)果不僅可以加深對(duì)爆炸切割器作用機(jī)理和作用過(guò)程的認(rèn)識(shí)，而且，通過(guò)改變裝藥量的多少，可以得到連接桿不同程度地?cái)嗔?，為爆炸切割器裝藥量的優(yōu)化設(shè)計(jì)提供參考。如果裝藥量過(guò)少，連接桿的斷裂不完全，會(huì)引起包帶解鎖的失敗；如裝藥量過(guò)多，在連接桿完全斷裂后會(huì)引起較大的沖擊環(huán)境，會(huì)對(duì)衛(wèi)星造成危害。因此，如何時(shí)選擇合適的裝藥量，保證包帶爆炸解鎖的可靠性，顯得十分必要。同時(shí)，仿真結(jié)果還可以為航天器分離沖擊環(huán)境預(yù)示提供數(shù)值上的參考。

1爆炸切割器的結(jié)構(gòu)及工作原理

包帶式連接結(jié)構(gòu)中的爆炸切割器的典型結(jié)構(gòu)示意圖見(jiàn)圖1，它由切割器殼體、剪切刀、連接桿和火炸藥四部分組成。切割器通過(guò)螺母與包帶的端頭連接，包帶的端頭除了切割器連接桿連接外，還通過(guò)連接彈簧連接，連接彈簧的作用是在切割器的連接桿斷裂后對(duì)包帶進(jìn)行連接限位約束。爆炸切割器的工作原理是：在預(yù)定時(shí)刻，電引爆器通過(guò)引爆火炸藥，炸藥爆炸產(chǎn)生大量的高溫高壓氣體驅(qū)動(dòng)剪切刀運(yùn)動(dòng)，剪切刀運(yùn)動(dòng)一段距離后達(dá)到一定的速度，依靠剪切刀的動(dòng)能與連接桿發(fā)生碰撞侵徹，從而沖斷連接桿，連接桿斷裂后，包帶瞬間失去預(yù)緊力的約束，迅速散開(kāi)，從而實(shí)現(xiàn)星箭連接結(jié)構(gòu)的解鎖。

圖1　爆炸切割器結(jié)構(gòu)示意圖 Fig.1 Structure diagram of explosive cutter

2爆炸切割器的有限元建模

2.1有限元模型

本文采用ANSYS/LS-DYNA來(lái)模擬爆炸切割器的剪切斷裂問(wèn)題。在建模過(guò)程中，將切割器簡(jiǎn)化為切割器殼體、主裝藥、剪切刀、連接桿四部分，省略去起爆引線、密封圈、連接彈簧及包帶連接件等。由于切割器結(jié)構(gòu)是關(guān)于軸線完全對(duì)稱，為了簡(jiǎn)化計(jì)算，取爆炸切割器簡(jiǎn)化模型的1/4模型來(lái)研究(見(jiàn)圖2)。圖2中，切割體外徑粗處5.6 cm, 細(xì)處4.4 cm, 內(nèi)徑粗處3.2 cm,細(xì)處2.0 cm, 剪切刀總長(zhǎng)6.5 cm, 剪切刀端頭與切割體緊密相貼，細(xì)處真徑為1.8 cm, 尖頭處真徑0.9 cm, 切割器與剪切刀之間組成藥室，藥室半徑1.6 cm, 藥室寬0.8 cm, 裝藥為扁平圓柱形裝藥，裝藥半徑1.2 cm,裝藥長(zhǎng)0.6 cm, 連接桿直徑1.2 cm, 半連接桿長(zhǎng)8 cm。

圖2　爆炸切割器簡(jiǎn)化模型示意圖(單位：cm) Fig.2 Simplified model of explosive cutter(unit:cm)

在ANSYS/LS-DYNA中，將圖2所示的切割器1/4簡(jiǎn)化模型，建立其三維有限元模型(見(jiàn)圖3)。其中，切割器殼體、剪切刀、連接桿、主裝藥，均采用SOLID164實(shí)體單元，拉格朗日網(wǎng)格。網(wǎng)格大部分采用掃掠規(guī)則劃分，切割器上被連接桿穿過(guò)的部分處采用局部細(xì)分，模型共計(jì)生成24 219個(gè)節(jié)點(diǎn)，20 014個(gè)單元。

圖3　爆炸切割器簡(jiǎn)化模型有限元三維模型 Fig.3 Simplified FEM of explosive cutter

2.2材料模型的選擇

材料模型包括切割器殼體、剪切刀、連接桿和炸藥。炸藥選用的是常用的黑索金(RDX)，其相關(guān)參數(shù)(見(jiàn)表1)。

表1　RDX炸藥的基本參數(shù)

在LS-DYNA中，炸藥材料采用內(nèi)置的高性能炸藥描述，壓力用JWL(Jones-Wilkins-Lee)狀態(tài)方程描述[13]。JWL狀態(tài)方程的標(biāo)準(zhǔn)形式為：

式中：p和V分別表示爆轟產(chǎn)物的壓力和相對(duì)比容，E0為單位體積爆轟產(chǎn)物的初始內(nèi)能，式中的右端三項(xiàng)分別在高、中、低壓力區(qū)起主要作用。A、B、R1、R2、ω為表征炸藥爆炸特性的參數(shù)，一般通過(guò)爆炸實(shí)驗(yàn)數(shù)據(jù)擬合得到；對(duì)于不同成分的炸藥，炸藥參數(shù)的取值不同，這些參數(shù)常常通過(guò)圓筒實(shí)驗(yàn)來(lái)確定。本文中炸藥采用常見(jiàn)的黑索金炸藥RX-04-DS，參數(shù)來(lái)自于autodyn 13.0 材料庫(kù)中的數(shù)據(jù)。

爆炸切割器的殼體部分選擇高強(qiáng)度合金鋼42CrMo，剪切刀選擇硬度和強(qiáng)度大的鉭鎢合金(TANT10%W)，連接桿選擇合金鋼4340，三都均采用彈性—理想塑性帶應(yīng)變失效的本構(gòu)模型[14]，模型參數(shù)按照autodyn 13.0材料庫(kù)中的數(shù)據(jù)及《機(jī)械設(shè)計(jì)手冊(cè)》中的參數(shù)，材料的物理力學(xué)參數(shù)見(jiàn)表2。

表2　材料的基本參數(shù)

2.3接觸的定義

炸藥與殼體和剪切刀之間定義了僅滑動(dòng)面面接觸，剪切刀與連接桿之間為侵蝕接觸，連接桿與殼體之間為自動(dòng)面面接觸。

2.4其它參數(shù)設(shè)置

如圖3所示的有限元模型中，共生成四個(gè)部件(PART)。由于1/4模型的對(duì)稱性，顯然，模型是關(guān)于XOZ平面循環(huán)對(duì)稱的，將XOZ平面上的所有節(jié)點(diǎn)定義為對(duì)稱邊界條件；同理將YOZ平面上的所有節(jié)點(diǎn)定義為對(duì)稱邊界條件。為了抑制零能模式，使用EDBIVS命令對(duì)沙漏進(jìn)行控制，采用全局調(diào)整體積粘性，設(shè)置線性參數(shù)和二次參數(shù)，設(shè)為默認(rèn)值1.5和0.06。同時(shí)，使用EDHGLS命令來(lái)增加沙漏系數(shù)，

對(duì)全局增加彈性剛度來(lái)控制沙漏，設(shè)置沙漏系數(shù)為0.1。為了防止仿真中出現(xiàn)負(fù)體積，本文的爆炸仿真中將時(shí)間步長(zhǎng)縮放系數(shù)設(shè)置為0.6。

3計(jì)算結(jié)果與分析

3.1爆炸切割器爆炸斷裂過(guò)程中裝藥量與藥室平均壓強(qiáng)之間的關(guān)系

通過(guò)改變主裝藥藥量的大小，可以仿真得到不同裝藥量下的爆炸切割器爆炸斷裂情況。當(dāng)裝藥量減少時(shí)，爆轟氣體產(chǎn)物的做功能力不足，剪切刀沒(méi)有足夠的速度和動(dòng)能來(lái)切斷連接桿；當(dāng)裝藥量增加時(shí)，爆轟氣體做功能力增加，爆炸載荷增大，剪切刀的速度加速度均增加，在行程不變的情況下，剪切刀與連接桿開(kāi)始接觸的時(shí)間與完全切斷連接桿的時(shí)間均減小，而在連接桿完全斷裂后，切割器殼體和剪切刀上的速度加速度響應(yīng)均增大。

本文中選擇的藥量是在連接桿能夠完全斷裂的情況下選取的。根據(jù)文中有限元模型中各部件體積和材料密度計(jì)算得到各部件的質(zhì)量，其中，切割器外殼為 333.95 g，剪切刀為 76.15 g，連接桿為 35.40 g，主裝藥為 1.26 g。

而在工程上根據(jù)經(jīng)典內(nèi)彈道理論，常用定容條件下密閉爆發(fā)器內(nèi)火藥燃燒最大壓強(qiáng)的經(jīng)驗(yàn)公式，來(lái)估算燃?xì)獾淖畲髩簭?qiáng)。引入火藥力f和裝填密度Δ，經(jīng)驗(yàn)公式為[15]：

(2)

式中：α為氣體余容，ω為火藥藥量，V0為藥室初始容積。

黑索金等容爆熱6 025 J/g，等效火藥力1 510 J/g，藥室初始容積V0為 1.61 cm3，氣體余容取 0.7 cm3/g，黑索金裝藥量ω為1.26 g，代入式(2)，得到Pm為2.62 GPa。根據(jù)式(1)，爆轟產(chǎn)物充滿藥室時(shí)，相對(duì)比容為2.37，得到壓強(qiáng)為2.51 GPa，兩者相比較，相差5%。

隨著剪切刀的運(yùn)動(dòng)，藥室的體積是增加的，當(dāng)剪切刀完全切斷連接桿時(shí)，剪切刀向前移動(dòng)了1.3 cm，根據(jù)式(2)，此時(shí)壓強(qiáng)為0.57 GPa。此時(shí)相對(duì)比容為6.22，根據(jù)式(1)，得到此時(shí)的壓強(qiáng)為0.66 GPa；兩者相比較，相差11%。

3.2爆炸切割器剪切斷裂過(guò)程的應(yīng)力分析

將切割器的有限元模型生成的K文件，加入炸藥參數(shù)編輯修改后，利用ANSYS/LS-DYNA求解，將結(jié)果文件導(dǎo)入軟件LS-Prepost中進(jìn)行后處理。炸藥從起爆時(shí)刻t=0到連接桿完全斷裂時(shí)刻t=220 μs期間的爆炸切割器的等效應(yīng)力變化情況如圖4(a)～圖4(d)所示。

圖4(a)中，在t=7.3 μs 時(shí)刻，炸藥的爆轟產(chǎn)物氣體導(dǎo)致充滿藥室，高壓氣體開(kāi)始推動(dòng)剪切刀運(yùn)動(dòng)，剪切刀與氣體接觸的端頭部分的等效應(yīng)力最大，爆炸氣體產(chǎn)生的應(yīng)力波沿切割器和剪切刀傳播。

(a) t=7.3μs炸藥氣體膨脹充滿藥室,推動(dòng)剪切刀前進(jìn)(b) t=35.2μs剪切刀開(kāi)始與連接桿相撞

(c) t=110.2μs剪切刀與連接桿碰撞斷裂期間,接觸區(qū)應(yīng)力較大(d) t=212.2μs連接桿發(fā)生沖塞斷裂,完全斷裂圖4 爆炸切割器爆炸斷裂過(guò)程中等效應(yīng)力變化云圖Fig.4Contoursofeffectivestressintheprocessofexplosivecutter’sfracture

圖4(b)中，在t=35.2 μs 時(shí)刻，剪切刀運(yùn)動(dòng)一段距離后，開(kāi)始與連接桿相撞。此時(shí)，切割器殼體藥室由于爆轟氣體的高壓作用，發(fā)生一定的“鼓包”膨脹，剪切刀與切割器內(nèi)壁摩擦接觸處有應(yīng)力集中，等效應(yīng)力最大。

圖4(c)中，在t=110.2 μs 時(shí)刻，剪切刀與連接桿相撞，剪切刀侵徹到連接桿中一定深度，剪切刀頭與連接桿接觸處的等效應(yīng)力最大，彈性應(yīng)力波在切割器殼體、剪切刀和連接桿之間來(lái)回反射、振蕩。

圖4(d)中，在t=212.2 μs 時(shí)刻，連接桿發(fā)生“沖塞”斷裂，連接桿斷裂為三部分。等效應(yīng)力最大處仍發(fā)生在剪切刀與殼體內(nèi)壁接觸摩擦部分。

3.3爆炸切割器連接桿斷裂處的應(yīng)力應(yīng)變分析

為了更好地研究連接桿的斷裂過(guò)程，選取連接桿上A、B、C、D四個(gè)典型單元，分析單元在斷裂過(guò)程中的應(yīng)力應(yīng)變變化，如圖5(a)所示。

圖5　爆炸切割器中連接桿上典型單元斷裂過(guò)程分析 Fig.5 Fracture analysis of marked elements in the link of the explosive cutter

圖5(b)中，連接桿上單元A在 45 μs達(dá)到失效應(yīng)變，單元B在 125 μs達(dá)到失效應(yīng)變，單元D在 200 μs達(dá)到失效應(yīng)變。

圖5c中，單元A在等效應(yīng)力達(dá)到屈服極限之后應(yīng)變隨即達(dá)到失效應(yīng)變，單元失效；單元B在等效應(yīng)力達(dá)到屈服極限之后, 等效應(yīng)力經(jīng)歷一個(gè)先減小后增加的振蕩的過(guò)程，應(yīng)變逐漸增加到失效應(yīng)變，單元失效；單元D在等效應(yīng)力達(dá)到屈服極限之后, 經(jīng)歷多次先減小后增加的振蕩的過(guò)程，而應(yīng)變幾乎保持為零，在失效時(shí)刻，應(yīng)變突然躍升達(dá)到失效應(yīng)變；單元C在等效應(yīng)力達(dá)到屈服極限之后, 經(jīng)歷多次先減小后增加的振蕩的過(guò)程之后，而應(yīng)變緩慢上升，最大值為 0.124 ，并未失效，單元C保持一個(gè)殘余應(yīng)力和殘余速度。

3.4爆炸切割器斷裂過(guò)程中各部件的整體速度分析

為了更好的分析剪切刀和連接桿的碰撞斷裂過(guò)程，以切割器殼體、剪切刀、連接桿為研究對(duì)象，研究三部件(PART)分別作為整體時(shí)質(zhì)心的速度變化規(guī)律，根據(jù)仿真結(jié)果文件，得到三部件質(zhì)心在三個(gè)方向上的整體速度隨時(shí)間變化曲線見(jiàn)圖6。

圖6(a)中，剪切刀在炸藥的瞬時(shí)爆轟作用下，軸向速度(即Z向)迅速增加，在 14 μs 時(shí)刻即達(dá)到 45 m/s, 爆轟氣體推動(dòng)切刀前進(jìn)一段距離后與連接桿相碰，在 39 μs 時(shí)刻切刀達(dá)到最大速度 45.8 m/s, 隨著切刀和連接桿的碰撞，連接桿開(kāi)始變形、斷裂，而切刀的速度逐漸減小，約在 210 μs 時(shí)刻連接桿斷裂完成，剪切刀的速度減小到 26.0 m/s。

圖6(b)和圖6(c)中，約在7.3 μs時(shí)刻，炸藥爆轟膨脹充滿藥室；此時(shí)，切割器殼體的藥室因?yàn)楸Z氣體的作用而發(fā)生“鼓包”膨脹現(xiàn)象。此時(shí)，剪切刀的速度最大值22.6 m/s，殼體的速度最大值16.5 m/s。然后，隨著炸藥氣體與殼體和剪切刀的耦合作用，殼體與切刀的速度在Y方向上下振蕩，逐漸衰減到零。連接桿的速度在斷裂后，在Y方向上存在殘余速度約為4.5 m/s，在X向上的始終為零。

(a) 殼體、剪切刀、連接桿沿軸向(Z軸)速度曲線(b) 殼體、剪切刀、連接桿沿徑向(Y軸)速度曲線(c) 殼體、剪切刀、連接桿沿切向(X軸)速度曲線圖6 爆炸切割器中殼體、剪切刀、連接桿的三向速度曲線Fig.6Velocitycurvesofcase,knifeandlinkoftheexplosivecutter

(a) 殼體、剪切刀、連接桿沿軸向(Z軸)加速度曲線(b) 殼體、剪切刀、連接桿沿軸向(Z軸)加速度曲線

(c) 殼體、剪切刀、連接桿沿徑向(Y軸)加速度曲線(d) 殼體、剪切刀、連接桿沿切向(Y軸)加速度曲線圖7 爆炸切割器中殼體、剪切刀、連接桿的三向加速度曲線Fig.7Accelerationcurvesofcase,knifeandlinkoftheexplosivecutter

3.5爆炸切割器斷裂過(guò)程中各部件的整體加速度分析

選取切割器殼體、剪切刀、連接桿為研究對(duì)象，研究三者的整體質(zhì)心加速度變化規(guī)律，根據(jù)仿真結(jié)果文件，得到切割器三部件質(zhì)心在三個(gè)方向上的整體加速度隨時(shí)間變化曲線(見(jiàn)圖7)。

從圖7(a)可知，從0 μs時(shí)刻炸藥起爆開(kāi)始，剪切刀在爆轟氣體作用于1.75 μs 時(shí)刻達(dá)到最大Z向加速度，約為1 600 000 g, 與此同時(shí)，切割體殼體也同時(shí)達(dá)到反向Z向最大加速度，約為360 000 g。連接桿的Z向最大加速度發(fā)生在30 μs 時(shí)刻，約為110 000 g。在40 μs 時(shí)刻，剪切刀運(yùn)動(dòng)一段時(shí)程后與連接桿相遇，在剪切刀和連接桿碰撞過(guò)程中，由圖7(b)可知，連接桿在Z向的加速度保持在±50 000 g上下振蕩衰減。在圖7(c)中，在1.75 μs 時(shí)刻，切割器殼體的Y向最大加速度約為600 000 g；約在4.7 μs 時(shí)刻，剪切刀的Y向最大加速度約為650 000 g。由圖7(d)可知，在剪斷連接桿的過(guò)程中，剪切刀和切割器殼體在Y向上的加速度在±50 000 g上下振蕩衰減。由于對(duì)稱性，切割器殼體和剪切刀的X向加速度特性大小與Y向幾乎一樣。

4結(jié)論

本文以星箭連接結(jié)構(gòu)中典型的爆炸切割器為例，建立了爆炸切割器的三維有限元模型，分析了切割器爆炸沖擊斷裂過(guò)程中連接桿失效單元的應(yīng)力應(yīng)變變化規(guī)律，通過(guò)數(shù)值模擬得到了剪切刀和連接桿在炸藥爆炸載荷下的速度和加速度響應(yīng)，得到幾點(diǎn)結(jié)論：

(1)炸藥爆炸瞬間完成，在10 μs之內(nèi)炸藥能量釋放完畢，通過(guò)爆轟氣體與相鄰結(jié)構(gòu)的作用，通過(guò)應(yīng)力波傳播，轉(zhuǎn)化為相鄰結(jié)構(gòu)的應(yīng)變能和動(dòng)能，在炸藥附近的結(jié)構(gòu)會(huì)發(fā)生較大的塑性變形。

(2)切割器中連接桿的斷裂是依靠剪切刀的動(dòng)能侵徹和爆轟氣體余壓的共同作用完成的，連接桿斷裂處發(fā)生了彈塑性斷裂。

(3)連接桿上的單元在達(dá)到屈服極限時(shí)并不立即斷裂，而是進(jìn)入硬化強(qiáng)化階段，塑性應(yīng)變?cè)黾?，在?yīng)變積累到一定程度時(shí)，斷裂發(fā)生。連接桿的斷裂是由最大等效塑性應(yīng)變控制的。

(4)在其它條件不變的情況下，改變裝藥量的多少，可以得到連接桿斷裂的不同程度。仿真結(jié)果可以應(yīng)用到此類爆炸切割器裝置裝藥量的優(yōu)化設(shè)計(jì)中。

(5)與內(nèi)彈道理論中的經(jīng)驗(yàn)公式相比，仿真中狀態(tài)方程的壓強(qiáng)與經(jīng)驗(yàn)公式得到的壓強(qiáng)相差在5%～11%之間，誤差范圍相差不大。

本文提出的仿真方法，對(duì)于爆炸切割器爆炸斷裂動(dòng)態(tài)全過(guò)程中的物理力學(xué)參數(shù)能夠進(jìn)行全面的描述，詳細(xì)揭示了爆炸動(dòng)態(tài)斷裂過(guò)程中的動(dòng)力學(xué)行為，對(duì)于爆炸切割器的斷裂機(jī)理的認(rèn)識(shí)，具有一定的參考價(jià)值。

參考文獻(xiàn)

[1]NASA-STD-7003A. Pyroshock Testcriteria[S].NASA Technical Standard. Washington, DC, USA, 2011.

[2]Brauer K O. Handbook of pyrotechnics[S].New York: Chemical Publishing Co. INC, USA,1974.

[3]Shmuel B S， Goldstein S. Numerical simulation of a bolt cutter [C]//34th AIAA/ASME/SAE/ASEE Joint Propulsion Conference, Cleveland, OH, July 1998..

[4]Shmuel B S. Sizing Bolt cutter assemblies for specific applications by hydrocode analysis [C]//36th AIAA/ ASME/SAE/ASEE Joint Propulsion Conference,Huntsville, AL, July 2000.

[5]ShmuelB S. Performance analysis of a normally closed pyrovalve [C]//38th AIAA/ASME/SAE/ASEE Joint Propulsion Conferenceand Exhibit, Indianapolis, IN, July, 2002.

[6]張建華. 航天產(chǎn)品的爆炸沖擊環(huán)境技術(shù)綜述[J]. 導(dǎo)彈與航天運(yùn)載技術(shù)，2005,3:30-36.

ZHANG Jian-hua. Pyroshock environment of missiles and launch vehicles[J]. Missles and Space Vehicles, 2005,3:30-36.

[7]譚雪峰，閻紹澤. 星箭包帶式連接結(jié)構(gòu)動(dòng)力學(xué)研究進(jìn)展[J]. 導(dǎo)彈與航天運(yùn)載技術(shù)，2010,1:1-6.

TAN Xue-feng，YAN Shao-ze. Evolution of research on dynamics of clamp band system between satellite and launch vehicle[J]. Missles and Space Vehicles, 2010,1:1-6.

[8]Li Jun-lan, Yan Shao-ze, Tan Xue-feng. Dynamic envelope analysis of clamp-band joint considering pyroshock of satellite separation. [J]. Journal of Spacecraft and Rockets, 2014, 51(5):1390-1400.

[9]Qin Zhao-ye, Yan Shzo-ze,Chu Fu-lei. Dynamic characteristics of launch vehicle and spacecraft connected by clamp band[J]. Journal of Sound and Vibration, 2011,330(10):2161-2173.

[10]武新峰，劉觀日，雷勇軍,等. 基于LS-DYNA的包帶式星箭連接裝置分離過(guò)程和沖擊響應(yīng)分析[J]. 振動(dòng)與沖擊，2013，32(24)：174-179.

WU Xin-feng, LIU Guan-ri, LEI Yong-jun, et al.Separation process and shock response analysis of clamp band systems based on LS-DYNA[J].Journal of Vibration and Shock, 2013，32(24)：174-179.

[11]初文華，張阿漫，時(shí)付紅，等. SPH-FEM耦合算法在爆炸螺栓解鎖分離過(guò)程中的應(yīng)用[J].振動(dòng)與沖擊，2012，31(23)：197-202.

CHU Wen-hua, ZHANG A-man, MING Fu-ren, et al. Application of three-dimensional SPH-FEM coupling method in unlocking process of an explosion bolt[J].Journal of Vibration and Shock, 2012，31(23)：197-202.

[12]Cen Zheng, Gao Ruo-fei, Shen Lei.LM-3A series launch vechicle user’s manual[M]. Beijing: CALT, CGWIC, CLTC, 2011.

[13]尚曉江，蘇建宇，王化鋒,等. ANSYS/LS-DYNA動(dòng)力分析方法與工程實(shí)例[M].北京：中國(guó)水利出版社，2008.

[14]郝好山，胡仁喜，康士廷,等. ANSYS12.0 LS-DYNA非線性有限元分析從入門到精通[M].北京：機(jī)械工業(yè)出版社，2010.

[15]金志明. 槍炮內(nèi)彈道學(xué)[M].北京：北京理工大學(xué)出版社，2004.