聚能切割索切割能力仿真及試驗研究

劉 寧，渠弘毅，戈慶明，王寅虎，洪 陽

聚能切割索切割能力仿真及試驗研究

劉 寧，渠弘毅，戈慶明，王寅虎，洪 陽

（北京航天長征飛行器研究所，北京，100076）

為深入研究聚能切割索的切割能力，采用數(shù)值仿真和試驗的手段研究了聚能切割索（Flexible Linear Shaped Charge，F(xiàn)LSC）切割不同材料鋁板的切割效果。仿真和試驗結果表明：材料的屈服強度對侵徹性能影響最大，隨著被切割對象屈服強度的增大，侵徹深度降低。

聚能切割索；數(shù)值仿真；屈服強度；侵徹深度

0 引 言

分離技術[1]是航天工程中的關鍵技術。線分離技術（包括聚能切割索[2]、柔性導爆索等）以其火工元件少、分離結構本身不承受大載荷、結構簡單、可靠、快速、同步性高等優(yōu)點而廣泛應用于國內(nèi)外眾多航天器的級間分離中[3]。

聚能射流的成形過程分為3個階段[4]：a）傳爆藥起爆主裝藥階段：主裝藥起爆后，爆轟波在主裝藥中高速傳播；b）藥型罩壓垮與射流成形階段：爆轟波傳至藥型罩，藥型罩內(nèi)表面金屬受爆轟波驅動，向藥型罩軸線高速壓合形成高速射流，外表面金屬形成射流杵體，在實際射流成形過程中，藥型罩底由于“角裂”現(xiàn)象的存在而向外飛出，并不參與射流成形過程；c）射流在空氣中拉伸自由運動階段：由于射流成形過程中存在速度梯度，射流速度由射流頭部沿軸線方向至射流尾部逐漸減小，因此射流在空氣中不斷拉伸，直至出現(xiàn)頸縮斷裂現(xiàn)象。

1 有限元模型

1.1 物理模型

線性聚能裝藥（Linear Shaped Charge，LSC）是聚能裝藥的一種，是利用聚能原理制成的一種帶有楔形金屬藥型罩的線性裝藥，原理是炸藥爆炸后的爆轟產(chǎn)物對藥型罩進行壓垮，形成一種薄片狀的金屬聚能刀來對目標進行切割。其基本理論是從錐形聚能裝藥理論引申和擴充而來的，用于切割金屬和其他材料。爆炸切割技術具有高速、高效、操作方便且不受環(huán)境限制等優(yōu)點，主要應用在導彈和航天飛機中的自毀系統(tǒng)和分離裝置，能夠在高溫和真空條件下對火箭、導彈和飛行器進行切割和分離。

采用Autodyn軟件建立二維有限元模型，如圖1所示。聚能射流（藥型罩、裝藥、空氣）采用Euler網(wǎng)格進行計算，鋁板采用網(wǎng)格進行計算，聚能切割索切割能力仿真采用Euler/Lagrange流固耦合算法進行求解。為了減少計算量，整個過程采用二維軸對稱計算模型。

圖1 切割索有限元模型

1.2 材料模型

a）空氣。

Euler網(wǎng)格中填充的空氣用理想氣體狀態(tài)方程描述：

b）紫銅。

紫銅具有密度大、熔點適中、聲速高、延展性好、不易斷裂等特點，是制備藥型罩的理想材料。紫銅制成的藥型罩能形成穩(wěn)定射流并保證成形過程中不發(fā)生氣化。紫銅的強度方程選用Steinberg-Guinan模型[5,6]，該模型可以描述固體在高應變率（105/s以上）下的響應，對于高應變率下的剪切模量和屈服應力有：

c）黑索金裝藥。

爆炸產(chǎn)物JWL狀態(tài)方程[7]不顯含化學反應，能夠精確描述爆炸產(chǎn)物的等熵膨脹過程，其具體形式為

式中為爆炸產(chǎn)物的相對比容；為炸藥的比內(nèi)能；，，，，，為常數(shù)。

d）鋁板。

在射流切割鋁板時，主要有3個特征：高溫、高壓和高應變率。為準確描述鋁板在此狀態(tài)下的響應規(guī)律，采用Johnson-Cook模型[8]來描述。Johnson-Cook模型常用于模擬金屬材料從低應變率到高應變率下的動態(tài)行為，該模型采用變量乘積關系描述了應變率、溫度和應變的影響，本構方程為

1.3 沖擊狀態(tài)方程

沖擊狀態(tài)方程[9]的關系式為

對于大多數(shù)材料，在沒有發(fā)生相變的情況下，沖擊波速度和質點速度的關系為

結合質量、動量和能量守恒方程可以得到：

1.4 初始及邊界條件

為避免應力波在邊界處發(fā)生反射并確保爆轟產(chǎn)物正常飛散，在Euler域邊界處設置流出邊界，模擬無限大空氣域。起爆方式選用裝藥頂端中心點起爆，建立計算模型如圖2所示。在射流運動方向上，每隔1 mm設置一個高斯點，共設置11個高斯點，記錄射流頭部速度的變化過程；在鋁板的切割方向（即豎直方向），每隔2個網(wǎng)格設置一個高斯點，記錄鋁板變形及當?shù)亓Ｗ拥募铀偾闆r；在鋁板的水平方向，由于其尺度較大，每隔5個網(wǎng)格設置一個高斯點，記錄沖擊波在鋁板內(nèi)部的傳播及反射等信息。

圖2 高斯點設置信息

數(shù)字—各監(jiān)測點的編號

2 聚能切割索切割能力仿真研究

選取了1100、6061、2024鋁3種不同強度的鋁板進行了切割能力仿真分析，計算參數(shù)見表1。

表1 模型中三種金屬材料參數(shù)

Tab.1 Three Parameters of Metal Materials in the Model

材料密度(g·cm-3)剪切模量GPa屈服應力MPa最大屈服應力MPa硬化系數(shù)β硬化指數(shù)n 1100鋁2.70727.1404804000.27 6061鋁2.70327.62906801250.1 2024鋁2.78528.62607603100.185

仿真計算結果如圖3至圖5所示。

圖3 1100鋁材料鋁板侵徹過程

圖4 6061鋁材料鋁板侵徹過程

圖5 2024鋁材料鋁板侵徹過程

材料為1100的鋁板侵徹深度最大為4 mm，材料為6061的鋁板侵徹深度為3 mm，材料為2024鋁的侵徹深度為2.6 mm。由仿真結果可知，對于6061鋁和2024鋁，射流侵徹鋁板初期，主要靠高速、高溫、高壓的射流切割鋁板，鋁板幾乎不發(fā)生變形；在侵徹過程中，射流速度不斷減小，侵徹能力不斷下降，當射流速度下降到約200 m/s時，幾乎失去侵徹能力，此時鋁板開始發(fā)生彎曲變形。射流頭部速度和平板上節(jié)點速度情況如圖6、圖7所示。由圖6、圖7可知，射流頭部速度與鋁板變形后的移動速度幾乎一致，說明此時射流已不再切割鋁板，可以認為此時銅粒子已經(jīng)與鋁板粘在儀器，保持相對靜止狀態(tài)。

圖6 射流頭部速度曲線

圖7 平板上節(jié)點速度

鋁板侵徹時的應力云圖見圖8，1100鋁板背面受到的最大拉應力約為142 MPa，遠大于其屈服應力 40 MPa，鋁板完全被“射流刀”切開；6061鋁板背面受到的最大拉應力約為249 MPa，小于其屈服應力 290 MPa，故其背面并未拉斷；2024鋁板背面受到的最大拉應力約為261 MPa，與屈服應力相當，但由于背部未切斷部位較厚，故未拉斷。

續(xù)圖8

3 試驗與結果分析

3.1 試驗方法

采用3種不同材料鋁板進行切割試驗，如圖9所示，鋁板厚度為4.2 mm。藥型罩材料為紫銅，密度為 8.93 g/cm3，高度為（3±0.1）mm，頂角壁厚大于0.15 mm，V型深度為（1.3±0.1）mm，寬度為3.5 mm，炸高約1.5 mm；裝藥為散裝黑索金，線密度為（4.5±0.3）g/m，爆速在7500 m/s以上。

采用電雷管軸向起爆的方式起爆切割索，雙邊固定鋁板。

圖9 鋁板切割試驗產(chǎn)品

3.2 結果分析

試驗后，回收鋁板殘骸，觀察切割截面并測量侵徹厚度和拉斷厚度。鋁板的切割截面劃分為：侵徹區(qū)、拉斷區(qū)和未切斷區(qū)。其中侵徹區(qū)表面光滑，有熔融痕跡；拉斷區(qū)表面不規(guī)則，有細小的尖銳角，拉斷區(qū)和侵徹區(qū)之間有明顯的分界面。

經(jīng)測量，材料為1100的鋁板侵徹厚度在 3.8~4 mm，材料為6061的鋁板侵徹深度為 2.9~3.2 mm，材料為2024鋁板侵徹深度為2.6~2.7 mm，與仿真結果吻合較好。

4 模型應用注意點及結果修正方法

4.1 模型應用注意點

由于柔性切割索本身尺寸較小，為保證計算精度，應控制藥型罩的厚度方向至少有3層網(wǎng)格，即網(wǎng)格尺寸應不大于0.05 mm。在作用過程中，炸藥及藥型罩會發(fā)生大變形，若采用Lagrange單元會產(chǎn)生網(wǎng)格畸變，導致計算終止，因此炸藥及藥型罩采用Euler網(wǎng)格建模，與被切割鋁板間采用自動流固耦合算法。此外，Lagrange網(wǎng)格尺寸應不大于Euler網(wǎng)格尺寸的2倍。

4.2 結果修正方法

金屬材料的本構模型比較成熟，需要計算校正的參數(shù)為導爆索裝藥黑索金的JWL狀態(tài)方程參數(shù)。由于柔性導爆索裝藥為散裝黑索金，無法通過圓筒試驗獲取其JWL狀態(tài)方程參數(shù)，因此，以某型炸藥的JWL狀態(tài)方程參數(shù)為基礎進行數(shù)值模擬，仿真計算結果與試驗結果進行對比，不斷調試校正以獲取散裝黑索金的JWL狀態(tài)方程參數(shù)。

5 結束語

通過對聚能切割索切割鋁板的數(shù)值仿真和試驗結果研究，聚能切割索切割鋁板主要分為3個階段：a）高速射流侵徹鋁板，鋁板在高溫、高壓、高速的射流作用下發(fā)生侵徹；b）射流速度不斷降低，最終銅粒子和鋁板粘在一起，未切斷部分部段拉伸，鋁板厚度部段減小；c）在拉應力的作用下，當拉應力大于材料屈服應力時，鋁板整體被切斷。

材料屈服強度對侵徹性能影響最大，隨著被切割鋁板屈服強度的增大，侵徹深度降低。

[1] 陳榮, 等. 爆炸分離裝置中保護罩安全性能分析[J]. 導彈與航天運載技術, 2007(4): 13-16.

Chen Rong, et al. Safety performance analysis of protective cover in explosive separation device[J]. Missiles and Space Vehicles, 2007(4): 13-16.

[2] Lan Whaley. Development of the stars II shroud separation system[R]. AIAA 2001-3769, 2001.

[3] Jiang Ruigang. Discussion on tactical missile stage linear separation[J]. Modern Defence Technology, 2001, 29(3): 31-33.

[4] 黃正祥. 聚能裝藥理論與實踐[M]. 北京: 北京理工大學出版社, 2004.

Huang Zhengxiang. Theory and practice of shaped charge[M]. Beijing: Beijing University of Technology Press, 2004.

[5] Steinberg D J, Cochran S G, Guinan M W. A constitutive model for metals

applicable at high-strain rate[J]. Journal of Applied Physics, 1980(51): 1498.

[6] Steinberg D J, Lund C M. A constitutive model for strain rates from 10-4 to10-6S-1[J]. Journal of Applied Physics, 1989(65): 1528.

[7] Lee E, Horning H, Kury J. Adiabatic expansion of high explosive detonation products[R]. Lawrence National Laboratory Report UCRI-50422, 1968.

[8] Johnson G R, Cook W H. A constitutive model and data for metals subjects to large strains, high strain rates and high temperature[C]. New York: the 7th Int. Symposiu on Ballistics, 1983.

[9] McQueen R, et al. The equation of state of solids from shock wave measurements, high-velocity impact[M]. New York: Academic Press, 1970.

Simulation and Experiment of Cutting Ability of Flexible Linear Shaped Charge

Liu Ning, Qu Hong-yi, Ge Qing-ming, Wang Yin-hu, Hong Yang

(Beijing Institute of Space Long March Vehicle, Beijing, 100076)

Material separation mechanisms are studied from the view point of cutting ability of flexible linear shaped charge. The effect of flexible linear shaped charge cutting different aluminum is analyzed by numerical simulation and experiment method. The simulation and experiment results proved that the material yield intension is the main complication for penetration performance. The penetration depth induces with the increasing of the yield strength of the cutting target.

flexible linear shaped charge; numerical simulation; yield strength; penetration depth

V416

A

1004-7182(2020)01-0038-05

10.7654/j.issn.1004-7182.20200107

2019-03-08；

2019-12-11

劉 寧（1982-），男，高級工程師，主要研究方向為飛行器設計。

渠弘毅（1986-），男，工程師，主要研究方向為飛行器設計。

戈慶明（1983-），男，高級工程師，主要研究方向為飛行器設計。

王寅虎（1986-），男，工程師，主要研究方向為飛行器火工設計。

洪 陽（1989-），男，工程師，主要研究方向為飛行器設計。