蔡尚曄， 姜春蘭， 毛亮， 王在成， 胡榕， 葉勝

(北京理工大學(xué) 爆炸科學(xué)與技術(shù)國家重點(diǎn)實(shí)驗(yàn)室， 北京 100081)

0 引言

聚四氟乙烯/鋁(PTFE/Al)復(fù)合材料在一般環(huán)境中呈惰性狀態(tài)[1]，但在撞擊加載條件下，沖擊能量會促使活性金屬Al與氟聚物PTFE發(fā)生劇烈的化學(xué)反應(yīng)，形成類爆轟效應(yīng)。富鋁含量PTFE/Al含能材料由于其材料密度和聲速較高，用作彈藥戰(zhàn)斗部毀傷元材料具備更加廣闊的應(yīng)用前景。

近年來，國內(nèi)外研究學(xué)者針對PTFE基含能材料的組分配比、顆粒度、致密度、材料反應(yīng)特性以及細(xì)觀模擬等方面開展了大量的研究工作。Shen[2]采用泰勒桿撞擊實(shí)驗(yàn)，研究了材料配比以及桿撞擊速度對PTFE/Ti含能材料撞擊引發(fā)變形及點(diǎn)火時間的規(guī)律。Wu等[3]采用落錘試驗(yàn)表征PTFE/Al的反應(yīng)特性，得出特征落差高度和點(diǎn)火能量隨著Al粒徑的增加而增加。Mcgregor等[4]研究了致密度為40%的PTFE/Al含能材料撞擊下的沖擊反應(yīng)閾值。Dolgoborodov等[5]發(fā)現(xiàn)疏松PTFE/Al(密度0.4～0.5 g/cm3)含能材料在1 GPa的沖擊壓力下可以形成穩(wěn)定爆轟。Koch等[6]通過監(jiān)測金屬氟碳基反應(yīng)材料反應(yīng)過程，得出金屬氟碳基體系中材料能否反應(yīng)與溫度有較大關(guān)聯(lián)。Feng等[7]開展了PTFE/Al反應(yīng)材料的準(zhǔn)靜態(tài)壓縮實(shí)驗(yàn)，根據(jù)實(shí)驗(yàn)現(xiàn)象提出了PTFE/Al裂紋誘導(dǎo)的引發(fā)機(jī)制。麥麥提圖爾蓀等[8]建立了包含顆粒和孔洞的隨機(jī)分布細(xì)觀力學(xué)有限元模型，對該材料的準(zhǔn)靜態(tài)壓縮力學(xué)行為進(jìn)行模擬。

孔隙作為材料制備中重要的考慮因素，對材料的沖擊溫升會產(chǎn)生重要的影響，從而會影響材料的沖擊反應(yīng)閾值、反應(yīng)時間等。鑒于前人在孔隙度對PTFE/Al含能材料沖擊溫升影響方面的研究甚少，尤其是針對富鋁含量的PTFE/Al(質(zhì)量配比50/50)含能材料。本文將采用含孔隙PTFE/Al沖擊溫升效應(yīng)模型進(jìn)行理論分析，并建立含孔隙PTFE/Al細(xì)觀模型進(jìn)行數(shù)值模擬，隨后采用分離式霍普金森壓桿(SHPB)實(shí)驗(yàn)驗(yàn)證的方法對富鋁含量PTFE/Al含能材料的溫升規(guī)律進(jìn)行研究，以期對富鋁含量PTFE/Al含能材料的工程化應(yīng)用提供參考。

1 理論分析

1.1 含孔隙PTFE/Al沖擊溫升效應(yīng)模型

含孔隙PTFE/Al含能材料屬于多孔洞材料，在沖擊加載下材料內(nèi)部的孔洞會發(fā)生塌陷、閉合，材料的不均勻性引起內(nèi)能的增加集中于孔洞周邊，形成“熱點(diǎn)”，最終“熱點(diǎn)”導(dǎo)致發(fā)生點(diǎn)火反應(yīng)[9]。由于孔洞形狀、數(shù)量和位置等變化多樣，在受到?jīng)_擊加載時研究孔洞的塌陷、溫升等問題也十分復(fù)雜。本文在Whitworth[10]所提出的粘塑性孔洞塌縮溫升模型基礎(chǔ)上，引入孔洞周邊介質(zhì)的融化效應(yīng)[11]，由此建立了含孔隙PTFE/Al的溫升效應(yīng)模型，如圖1所示，其中p為胞元外邊界周期性壓力，a、b分別為t時刻孔洞、胞元的半徑(μm)。首先，作如下假設(shè)：

1)孔洞附近的流場是球?qū)ΨQ的；

2)粘塑性區(qū)域以及融化區(qū)域?yàn)槊軐?shí)態(tài)，且不可壓縮；

3)以胞元外邊界周期性壓力p近似沖擊加載過程；

4)PTFE/Al復(fù)合材料在PTFE熔點(diǎn)處發(fā)生融化，且孔洞附近的溫升與半徑有關(guān)。

圖1 考慮了熔化效應(yīng)的一維粘塑性孔洞塌縮模型示意圖Fig 1 One-dimensional viscoplastic cavity collapse model considering melting effect

胞元的初始孔隙度α為

(1)

式中：ρt為材料的理論密度(g/cm3)；ρ為材料的實(shí)測密度(g/cm3).

胞元的孔隙塌縮動力學(xué)方程由(2)式[11]描述：

(2)

(3)

式中：pg0為孔隙中的初始?xì)怏w壓力，取標(biāo)準(zhǔn)大氣壓pg0=101 kPa=1.0110-6GPa；γ為絕熱指數(shù)，取γ=1.4. 由于要考慮熔化效應(yīng)，Y(pT,T)和μ(pT,T)不能取作常數(shù)，而必須與pT=(r,t)和T=(r,t)有關(guān)。因此，它們又可寫成Y=(r,t)和μ=(r,t)。取熔點(diǎn)溫度Tm(對于PTFE/Al混合物熔點(diǎn)溫度為600 K),屈服應(yīng)力Y=(r,t)和粘度μ=(r,t)的函數(shù)形式如下：

(4)

(5)

式中：Y0為初始屈服應(yīng)力；μ0為初始粘度；Ta為粘性流體的激活溫度(K)。為使問題簡化，在(6)式中忽略了粘度對壓力的依賴關(guān)系。

孔洞內(nèi)徑塌縮運(yùn)動而產(chǎn)生的溫升[10]可表示為

(6)

對(2)式和(6)式聯(lián)立進(jìn)行求解，便可獲得含孔隙PTFE/Al含能材料在受到外界壓力作用時，內(nèi)部孔洞的孔徑變化以及所引起的溫升變化。

1.2 含孔隙PTFE/Al沖擊溫升計算結(jié)果

基于上述沖擊溫升分析模型，針對本文所研究的孔隙度為10%、20%、30%的PTFE/Al含能材料進(jìn)行計算，其初始計算參數(shù)如表1[11-13]所示。

表1 PTFE/Al材料參數(shù)

圖2給出了經(jīng)過數(shù)據(jù)處理后所得到的在周期性壓力加載下孔隙度為10%、20%、30%的PTFE/Al含能材料內(nèi)部孔洞的孔徑和溫度隨時間的變化曲線。Wang等[14]和李尉等[15]證實(shí)了在SHPB加載未反應(yīng)的PTFE/Al為周期性加載，這種方式的加載會在材料反應(yīng)之前撞擊材料約2～3次。參考文獻(xiàn)[16]，這種周期性加載會在計算孔徑和溫度隨時間變化時存在兩個延時平臺。

圖2 不同孔隙度下PTFE/Al內(nèi)部孔洞孔徑、溫度隨時間的理論計算變化曲線Fig.2 Diameter-time and temperature-time curves of PTFE/Al materials with different porosities

從圖2(a)中可以看出，孔隙度為10%、20%及30%的PTFE/Al試件其孔洞尺寸隨時間逐漸減小，在孔隙塌縮的最后階段孔隙附近介質(zhì)被加速至較高的速度，孔洞半徑迅速減小?？紫抖容^高的試件孔洞變化率的值首先達(dá)到最大，宏觀上表現(xiàn)為材料容易被壓縮。

為進(jìn)一步計算沖擊溫升與孔隙度的關(guān)系，根據(jù)圖2中的10%孔隙度左右溫升較高的結(jié)果，在1.2%～20%孔隙度之間取若干個孔隙度數(shù)值在800 μs時刻進(jìn)行計算。表2和圖3給出了在800 μs時刻材料的溫度- 孔隙度計算結(jié)果和曲線，除1.2%、10%、20%和30%孔隙度的材料外，其余孔隙度的屈服強(qiáng)度采用線性插值法計算得出。從圖3可以看出：孔隙度在1.2%～14%的區(qū)域內(nèi)，到達(dá)800 μs時的溫度呈現(xiàn)出隨著孔隙度的上升而上升；在14%～20%的區(qū)域內(nèi)，到達(dá)800 μs時的溫度呈現(xiàn)出隨著孔隙度的上升而下降。因此，分析認(rèn)為能夠使溫升速率最高的孔隙度在10%～20%之間。

表2 800 μs時PTFE/Al孔隙度與沖擊溫度

圖3 800 μs時的溫度- 孔隙度曲線Fig.3 Temperature-porosity curve at 800 μs

2 數(shù)值模擬

2.1 含孔隙PTFE/Al有限元模型建立

為進(jìn)一步得出化學(xué)反應(yīng)被激發(fā)之前不同孔隙度PTFE/Al材料沖擊點(diǎn)火規(guī)律，需要對材料內(nèi)部溫度與時間的關(guān)系進(jìn)行分析。受到測試技術(shù)的制約，動態(tài)加載下材料溫度隨時間變化規(guī)律的系統(tǒng)性數(shù)據(jù)很難獲得，因此本文基于AUTODYN-2D的仿真平臺，采用任意拉格朗日- 歐拉算法對SHPB加載下的富鋁含量PTFE/Al含能試件進(jìn)行數(shù)值仿真，材料本構(gòu)模型選取Johnson-Cook模型[17]。該Johnson-Cook模型的材料準(zhǔn)靜態(tài)力學(xué)參數(shù)取自文獻(xiàn)[18]配比為50/50的PTFE/Al，結(jié)合實(shí)驗(yàn)所得動態(tài)力學(xué)參數(shù)擬合得到的材料的特定本構(gòu)方程參數(shù)[19]。模擬的試件尺寸為φ8 mm×5 mm. 參考理論分析結(jié)果，模擬試件的孔隙度為10%、20%和20%孔隙度。作如下假設(shè)：

1) 根據(jù)SHPB結(jié)構(gòu)的對稱性，為提高計算效率，本文建立如圖4的1/2仿真模型；

2) 材料孔洞采用的是理想的球形孔，孔洞大小以及孔洞位置在一定范圍內(nèi)隨機(jī)得出；

3) 所建立的模型中不同孔隙度的表現(xiàn)為材料孔洞占整個材料面積之比；

4) 入射桿設(shè)定的是周期性恒定壓力加載，而材料中的沖擊溫升基于文獻(xiàn)[16]實(shí)驗(yàn)中的入射桿波形周期性來分析。

圖4 細(xì)觀仿真模型Fig.4 Mesoscopic simulation model

本文帶不同孔隙材料的建模方法為：利用文獻(xiàn)[20]的方法得出m個隨機(jī)點(diǎn)(xm，ym，Rm)作為帶孔材料的孔洞半徑和圓心，孔洞半徑分布于集合[Rmin,Rmax]之間；并將所得出的m個隨機(jī)點(diǎn)(xm，ym，Rm)以命令流的方式導(dǎo)入到網(wǎng)格劃分軟件Truegrid中，對模型進(jìn)行合理的網(wǎng)格劃分后，生成zon格式文件。

2.2 仿真計算結(jié)果與分析

圖5～圖8為5 150 s-1應(yīng)變率下不同孔隙度PTFE/Al材料的沖擊溫升圖。從圖5～圖8中可以看出：在0～200 μs內(nèi)，材料在沖擊加載下開始發(fā)生變形且材料內(nèi)部孔洞主要被壓合的方向?yàn)檠貤U的運(yùn)動方向，材料內(nèi)部溫度隨著入射桿周期性的加載逐漸升高，這種材料升溫主要集中在孔隙邊緣；在200～500 μs內(nèi)，隨著時間的增加材料內(nèi)部溫度進(jìn)一步升高，孔洞在加載方向被壓合；在500～800 μs內(nèi)，除材料內(nèi)部溫度升高外，由于入射桿與透射桿尺寸的限制，材料擠出桿的邊緣，此時邊緣部分材料呈不規(guī)則狀。

圖5 5 150 s-1應(yīng)變率下不同孔隙度PTFE/Al材料在0 μs時刻的沖擊溫度云圖Fig.5 Impact temperatures of PTFE/Al materials with different porosities at 0 μs under 5 150 s-1 strain rate

圖6 5 150 s-1應(yīng)變率下不同孔隙度PTFE/Al材料在200 μs時刻的沖擊溫度云圖Fig.6 Impact temperatures of PTFE/Al materials with different porosities at 200 μs under 5 150 s-1 strain rate

圖7 5 150 s-1應(yīng)變率下不同孔隙度PTFE/Al材料在500 μs時刻的沖擊溫度云圖Fig.7 Impact temperatures of PTFE/Al materials with different porosities at 500 μs under 5 150 s-1 strain rate

圖8 5 150 s-1應(yīng)變率下不同孔隙度PTFE/Al材料在800 μs時刻的沖擊溫度云圖Fig.8 Impact temperatures of PTFE/Al materials with different porosities at 800 μs under 5 150 s-1 strain rate

圖9給出了不同孔隙度材料沖擊壓縮下溫度- 時間曲線，材料的溫度值取沖擊溫度云圖中最大溫度值。從圖9可以看出：在0～500 μs時間段內(nèi)，在入射桿周期性加載下，材料內(nèi)部溫度隨孔洞的塌縮呈現(xiàn)出上升趨勢；在500～800 μs時間段內(nèi)，此時材料受到單向載荷的壓縮孔洞被壓扁，此時材料內(nèi)部為密實(shí)狀態(tài)，在SHPB的加載下材料內(nèi)部溫度繼續(xù)升高但溫升幅度較小[21]。相比于20%及30%孔隙度的PTFE/Al含能材料，10%孔隙度材料在0～800 μs時間段內(nèi)溫度升高最高。

圖9 數(shù)值模擬所得不同孔隙度PTFE/Al溫度- 時間曲線Fig.9 Simulated time-temperature curves of PTFE/Al materials with different porosities

數(shù)值模擬所得數(shù)據(jù)與理論計算對比如圖10所示。從圖10中可以看出：在0～500 μs區(qū)間內(nèi)仿真數(shù)據(jù)所得溫度- 時間與理論計算所得溫度- 時間數(shù)據(jù)吻合較好；在700～900μs時間段內(nèi)，仿真所得數(shù)據(jù)與理論計算有些許差別，仿真值在時間上稍微“滯后”于理論計算值，這是由于孔洞塌縮溫升模型進(jìn)行了孔隙附近的流場是球?qū)ΨQ以及胞元外邊界壓力為恒定p的假設(shè)，當(dāng)材料受到單向載荷的壓縮時，孔洞會呈現(xiàn)出橢圓形的壓縮方式(見圖6)，這樣的壓縮方式會使得材料孔洞壓縮的面積略小于圓形孔洞塌縮的面積，從而在材料內(nèi)部溫升達(dá)到相同值時仿真的溫升時間較為“滯后”。

圖10 材料沖擊溫度- 時間理論計算與數(shù)值模擬結(jié)果對比Fig.10 Comparison of theoretical calculated and simulated impact temperature-time data

2.3 實(shí)驗(yàn)驗(yàn)證

為驗(yàn)證數(shù)值模擬的準(zhǔn)確性，本文以10%、20%以及30%孔隙度的含能材料作為研究對象(質(zhì)量配比50/50，試件尺寸φ8 mm×5 mm)，利用20 000幀/s的高速相機(jī)記錄PTFE/Al含能材料在SHPB撞擊加載下的沖擊引發(fā)化學(xué)反應(yīng)過程。

圖11 孔隙度為10%的PTFE/Al試樣在5 121 s-1應(yīng)變率加載下的高速攝影圖像Fig.11 High-speed photographic image of a 10% porosity PTFE/Al specimen under 5 121 s-1 strain rate

圖11～圖13為5 121～5 194 s-1應(yīng)變率下0～1 000 μs高速攝影圖像。假定入射桿初始撞擊試件的時刻為0 μs，從高速攝影圖像中可以看出：在0～500 μs時，試件在撞擊后受到強(qiáng)烈壓縮，發(fā)生連續(xù)變形，此時試樣開始失效，材料準(zhǔn)備飛濺碎屑；在500～800 μs時，試樣被進(jìn)一步壓縮，造成破裂并產(chǎn)生一定的碎屑，這與圖4～圖7數(shù)值模擬結(jié)果具有相同趨勢；在850 μs時，可以觀察到10%孔隙度PTFE/Al試樣激發(fā)發(fā)生化學(xué)反應(yīng)，并產(chǎn)生火光；在1 000 μs時，火光仍未結(jié)束，而20%、30%孔隙度PTFE/Al試樣未觀測到火光的產(chǎn)生，即此時未發(fā)生化學(xué)反應(yīng)。

圖12 孔隙度為20%的PTFE/Al試樣在5 154 s-1應(yīng)變率加載下的高速攝影圖像Fig.12 High-speed photographic image of a 20% porosity PTFE/Al specimen under 5 154 s-1 strain rate

圖13 孔隙度為30%的PTFE/Al試樣在5 194 s-1應(yīng)變率加載下的高速攝影圖像Fig.13 High-speed photographic image of a 30% porosity PTFE/Al specimen under 5 194 s-1 strain rate

文中PTFE/Al材料設(shè)計的為含有孔隙的材料，符合沖擊加載下能夠產(chǎn)生“熱點(diǎn)”并激發(fā)反應(yīng)的材料[22]且類似實(shí)驗(yàn)中也能夠看到火光產(chǎn)生[23]。材料在壓縮過程中會發(fā)生沖擊溫升且材料在這個沖擊加載下是等熵絕熱的[21]。因此，材料沖擊溫升為點(diǎn)火起爆過程的主要因素。這種10%孔隙度PTFE/Al試樣激發(fā)發(fā)生化學(xué)反應(yīng)，而20%、30%孔隙度PTFE/Al試樣未發(fā)生化學(xué)反應(yīng)的實(shí)驗(yàn)現(xiàn)象驗(yàn)證了數(shù)值模擬中相同應(yīng)變率下10%孔隙度PTFE/Al材料能夠產(chǎn)生較高溫升的結(jié)果。

3 結(jié)論

本文采用考慮熔化效應(yīng)的一維粘塑性孔洞塌縮模型對含孔隙富鋁PTFE/Al含能材料的沖擊溫升進(jìn)行了理論分析，并借助AUTODYN非線性動力有限元軟件進(jìn)行了沖擊加載的細(xì)觀數(shù)值模擬，最后通過實(shí)驗(yàn)進(jìn)行了驗(yàn)證。得到的主要結(jié)論如下：

2) 基于AUTODYN-2D數(shù)值模擬，含孔隙PTFE/Al含能材料在壓縮過程中孔隙首先被壓合，隨后材料再有較大變形；相比于孔隙度為20%和30%的富鋁含量PTFE/Al含能材料，孔隙度為10%的富鋁PTFE/Al含能材料(質(zhì)量配比50/50，試件尺寸φ8 mm×5 mm)在材料壓縮過程中溫度升高最高。

3) 含孔隙富鋁PTFE/Al含能材料在SHPB撞擊加載下的數(shù)值模擬與理論分析計算結(jié)果表明，材料內(nèi)部溫度隨入射桿的周期性加載總體呈現(xiàn)出間歇性升高的現(xiàn)象，且與實(shí)驗(yàn)結(jié)果吻合較好。

參考文獻(xiàn)(References)

[1] KOCH E C. Metal-fluorocarbon based energetic materials [M]. Weinheim, Germany: Wiley-VCH Verlag GmbH & Co. KGaA, 2012: 6-17.

[2] SHEN Y B. The chemical and mechanical behaviors of polymer/reactive metal systems under high strain rates [D].Atlanta, GA， US: Georgia Institute of Technology, 2012: 9-34.

[3] WU J X, FANG X, GAO Z R, et al. Investigation on mechanical properties and reaction characteristics of Al-PTFE composites with different Al particle size[J]. Advances in Materials Science and Engineering, 2018: ID 2767563.

[4] MCGREGOR N M, SUTHERLAND G T. Plate impact experiments on a porous teflon aluminum mixture[J]∥AIP Conference Proceedings， 2004, 706(1): 10.1063/1.1780406.

[5] DOLGOBORODOV A Y, MAKHOV M N, KOLBANEV I V, et al. Detonation in an aluminum-Teflon mixture[J]. Journal of Experimental & Theoretical Physics Letters, 2005, 81(7):311-314.

[6] KOCH E C，WEBB R，WEISER V. Review on thermo chemical codes [M]. Brussels, Belgium: NATO-Munitions Safety Information Analysis Center, 2011: 7-31.

[7] FENG B, LI Y C, WU S Z, et al. A crack-induced initiation mechanism of Al-PTFE under quasi-static compression and the investigation of influencing factors[J]. Materials & Design, 2016, 108:411-417.

[8] 烏布力艾散·麥麥提圖爾蓀, 葛超, 董永香, 等. 基于Al/PTFE真實(shí)細(xì)觀特性統(tǒng)計模型的宏觀力學(xué)性能模擬[J]. 復(fù)合材料學(xué)報, 2016, 33(11):2528-2536.

MAIMAITITUERSUN W，GE C，DONG Y X, et al. Simulation on mechanical properties of Al/PTFE based on mesoscopic statistical model [J]. Acta Materiae Compositae Sinica, 2016, 33(11):2528-2536. (in Chinese)

[9] DUAN Z P, WEN L J, LIU Y, et al. A pore collapse model for hot-spot ignition in shocked multi-component explosives [J]. International Journal of Nonlinear Sciences and Numerical Simulation, 2010, 11(Supplement):19-23.

[10] WHITWORTH N J. Mathematical and numerical modeling of shock initiation in heterogeneous solid explosives [D]. Bedfordshire, UK: Cranfield University, 2008.

[11] 張振宇, 浣石, 盧芳云. 多孔TNT炸藥中熱點(diǎn)形成的粘塑性塌縮機(jī)理[J]. 含能材料, 1994, 2(2): 36-48.

ZHANG Z Y, HUAN S, LU F Y. Viscoplastic collapse mechanism of hot-spot formation in porous TNT explosives[J]. Chinese Journal of Energetic materials, 1994, 2(2): 36-48. (in Chinese)

[12] MEYERS M A. Mechanical behavior of materials [J]. Upper Saddle River, NJ, US: Prentice Hall, 1998: 62-97.

[13] FREY R B. The initiation of explosive charges by rapid shear [C]∥Proceedings of the 7th Symposium on Detonation. Annapolis, MD, US: Narel Surface Weapon Center, 1981.

[14] WANG L, LIU J X, LI S K, et al. Investigation on reaction energy, mechanical behavior and impact insensitivity of W-PTFE-Al composites with different W percentage[J]. Materials & Design, 2016, 92:397-404.

[15] 李尉,任會蘭,寧建國,等 . Al/PTFE 活性材料的動態(tài)力學(xué)行為和撞擊點(diǎn)火特性[J]. 含能材料,2020,28(1):38-45.

LI W, REN H L, NING J G, et al. Dynamic mechanical beha-vior and impact ignition characteristics of Al/PTFE reactive materials[J]. Chinese Journal of Energetic Materials, 2020, 28(1): 38-45. (in Chinese)

[16] JIANG C L, CAI S Y, MAO L， et al. Effect of porosity on dynamic mechanical properties and impact response characteristics of high aluminum content PTFE/Al energetic materials [J]. Materials, 2020, 13(1): 140-150.

[17] 門建兵，蔣建偉，王樹有. 爆炸沖擊數(shù)值模擬技術(shù)基礎(chǔ)[M]. 北京: 北京理工大學(xué)出版社，2015: 97-146.

MEN J B, JIANG J W, WANG S Y. Fundamentals of numerical simulation for explosion and shock problems [M]. Beijing: Beijing Institute of Technology Press, 2015: 97-146. (in Chinese)

[18] 徐松林. Al/PTFE含能反應(yīng)材料力學(xué)性能研究[D].長沙: 國防科學(xué)技術(shù)大學(xué), 2010: 109-120.

XU S L. Study on the mechanical performance of polytetrafluoroethylene/Al energetic reactive materials[D]. Changsha: National University of Defense Technology, 2010: 109-120. (in Chinese)

[19] 李建光，施琪，曹結(jié)東. Johnson-Cook本構(gòu)方程的參數(shù)標(biāo)定[J]. 蘭州理工大學(xué)學(xué)報，2012，38(2)：164-167.

LI J G, SHI Q, CAO J D. Parameters calibration for Johnson-Cook constitutive equation [J]. Journal of Lanzhou University of Technology, 2012, 38(2): 164-167. (in Chinese)

[20] 王芳, 蔣建偉, 門建兵. 鎢銅射流成形的細(xì)觀數(shù)值模擬分析[J]. 兵工學(xué)報, 2018, 39(2):245-253.

WANG F, JIANG J W, MEN J B. Mesoscopic numerical simulation on the formation of tunsgsten-copper shaped charge jet [J]. Acta Armamentarii, 2018, 39(2):245-253. (in Chinese)

[21] 喬良. 多功能含能結(jié)構(gòu)材料沖擊反應(yīng)與細(xì)觀特性關(guān)聯(lián)機(jī)制研究[D]. 南京：南京理工大學(xué), 2013:13-34.

QIAO L. Study on the effects of mesoscale characteristic on shock reaction behavior of multifunctional energetic structural materials (MESMs) [D]. Nanjing: Nanjing University of Science and Technology, 2013: 13-34. (in Chinese)

[22] FYE P M . Initiation and growth of explosion in liquids and solids[J]. Journal of Chemical Education, 1953, 30(6):326.

[23] LI Y, WANG Z C, JIANG C L, et al. Experimental study on impact-induced reaction characteristics of PTFE/Ti composites enhanced by W particles[J]. Materials, 2017, 10(2): 175.