基于內(nèi)聚裂紋模型的高聚物粘結(jié)炸藥模擬材料動態(tài)斷裂行為研究

崔云霄， 陳鵬萬， 郭保橋， David A. Cendón， 周忠彬

(1.北京理工大學(xué) 爆炸科學(xué)與技術(shù)國家重點實驗室， 北京 100081； 2.西北核技術(shù)研究所， 陜西 西安 710024； 3.馬德里理工大學(xué)， 西班牙 馬德里 28040)

基于內(nèi)聚裂紋模型的高聚物粘結(jié)炸藥模擬材料動態(tài)斷裂行為研究

崔云霄1，2， 陳鵬萬1， 郭保橋1， David A. Cendón3， 周忠彬1

(1.北京理工大學(xué) 爆炸科學(xué)與技術(shù)國家重點實驗室， 北京 100081； 2.西北核技術(shù)研究所， 陜西 西安 710024； 3.馬德里理工大學(xué)， 西班牙 馬德里 28040)

對某高聚物粘結(jié)炸藥(PBX)模擬材料的動態(tài)拉伸斷裂行為進(jìn)行研究。針對該材料開展了基于霍普金森壓桿的動態(tài)帶預(yù)制裂紋半圓盤彎曲實驗，并結(jié)合高速攝像與數(shù)字圖像相關(guān)方法，得到了試樣動態(tài)破壞過程中的位移場和應(yīng)變場?；趦?nèi)聚裂紋模型，對其動態(tài)拉伸破壞過程進(jìn)行了數(shù)值模擬。數(shù)值模擬與實驗結(jié)果進(jìn)行對比后發(fā)現(xiàn)，拉伸應(yīng)力曲線、試樣破壞前后變形場等結(jié)果符合較好。根據(jù)數(shù)值模擬結(jié)果，分析了PBX試樣在動態(tài)預(yù)制裂紋半圓盤彎曲實驗過程中的裂紋擴(kuò)展演化規(guī)律，得到裂紋寬度比實驗結(jié)果偏小約15%的結(jié)論。

爆炸力學(xué)； 內(nèi)聚裂紋模型； 半圓盤彎曲實驗； 高聚物粘結(jié)炸藥； 數(shù)字圖像相關(guān)方法

0 引言

高聚物粘結(jié)炸藥(PBX)是戰(zhàn)斗部的關(guān)鍵組成部分，一般由奧克托今或黑索今炸藥顆粒、聚合物粘結(jié)劑和添加劑組成。準(zhǔn)確預(yù)測PBX的動態(tài)力學(xué)行為，特別是其內(nèi)部動態(tài)損傷斷裂，對于裝藥安全性評估具有重要意義[1]。PBX的拉伸強(qiáng)度較低，測試其動態(tài)拉伸力學(xué)性能多采用動態(tài)巴西實驗方法，但由于PBX材料的非均勻性，試樣在加載位置容易產(chǎn)生局部破壞，影響結(jié)果的可重復(fù)性[2]。采用帶預(yù)制裂紋半圓盤試樣進(jìn)行動態(tài)拉伸實驗，具有加工簡便、對應(yīng)力波傳播影響小的優(yōu)點，國際巖石力學(xué)學(xué)會建議采用該方法獲取巖石的動態(tài)斷裂韌性[3]。Zhou等[4]和Chen等[5]基于分離式霍普金森壓桿(SHPB)實驗系統(tǒng)，對PBX的半圓盤試樣進(jìn)行了動態(tài)彎曲實驗，結(jié)合高速攝像結(jié)果和數(shù)字圖像相關(guān)(DIC)方法，獲取了該材料的斷裂韌性。

PBX材料的力學(xué)行為復(fù)雜，通過數(shù)值模擬再現(xiàn)其破壞過程、揭示其內(nèi)在機(jī)理具有較大難度。Dienes等[6]基于耦合力- 熱- 化學(xué)的統(tǒng)計斷裂力學(xué)模型研究了多重沖擊實驗，認(rèn)為裂紋方向?qū)嚎s加載條件下PBX的脆性響應(yīng)有明顯影響。但是，該模型需要標(biāo)定多達(dá)200個參數(shù)，限制了其實際應(yīng)用。傅華等[7]建立了PBX的細(xì)觀計算模型，利用離散元方法模擬了PBX的動態(tài)巴西實驗，給出了試樣中損傷演化發(fā)展的過程，對于研究PBX的斷裂破壞有一定的啟示。王鵬飛等[8]將巖土材料研究常用的線性Durcker-Prager模型引入PBX的彈塑性變形研究中，模擬得到的單軸壓縮和雙軸壓縮下PBX的應(yīng)力和應(yīng)變關(guān)系與實驗結(jié)果較為符合。上述研究中對PBX的動態(tài)破壞行為研究的比較少，因此有必要開展更多的研究，以揭示其破壞機(jī)理及行為特征。

PBX與混凝土等準(zhǔn)脆性材料具有類似的宏觀力學(xué)性質(zhì)[9]，借鑒并改進(jìn)混凝土類材料的本構(gòu)模型可為研究PBX的動態(tài)斷裂行為提供新的途徑。本文基于內(nèi)聚裂紋模型對PBX動態(tài)半圓盤彎曲實驗進(jìn)行數(shù)值模擬，通過對比實驗測得的應(yīng)變信號，綜合高速攝像和數(shù)字圖像相關(guān)(DIC)方法獲得的位移場和應(yīng)變場驗證了模型的有效性，進(jìn)而根據(jù)模擬結(jié)果分析了PBX的裂紋擴(kuò)展過程。

1 內(nèi)聚裂紋模型的建立

內(nèi)聚裂紋模型最早提出是用于研究混凝土等準(zhǔn)脆性材料的拉伸斷裂問題[10-11]，經(jīng)過不斷發(fā)展和改進(jìn)，已用于PBX斷裂行為的模擬[12]。該模型可以描述材料任意位置拉伸裂紋的起裂或擴(kuò)展，以及裂紋擴(kuò)展導(dǎo)致的破壞行為。其基本思想是通過擴(kuò)展單元形函數(shù)來描述裂紋行為，一旦滿足起裂判據(jù)，即在單元內(nèi)部引入“虛擬”裂紋，并將內(nèi)聚力模型直接作用于“虛擬”裂紋兩側(cè)，避免了處理裂紋尖端奇異的困難。內(nèi)聚裂紋模型描述斷裂過程的示意圖如圖1所示。圖1中，t為裂紋兩側(cè)的拉應(yīng)力矢量，w為裂紋張開導(dǎo)致的位移間斷。

如果單元內(nèi)部沒有產(chǎn)生裂紋，則假設(shè)材料行為是彈性的。裂紋的起裂判據(jù)采用最大主應(yīng)力準(zhǔn)則，即一旦單元的最大主應(yīng)力超過給定拉伸強(qiáng)度閾值，就在單元內(nèi)部引入垂直于最大主應(yīng)力方向的“虛擬”裂紋。此時，單元處于完整材料與完全開裂材料之間的斷裂過程區(qū)。裂紋產(chǎn)生后材料開始軟化，隨著裂紋張開位移的增加，裂紋兩側(cè)的應(yīng)力開始下降。裂紋兩側(cè)的應(yīng)力t與裂紋張開位移w之間的關(guān)系如下：

(1)

下面簡述內(nèi)聚裂紋模型對單元內(nèi)部出現(xiàn)裂紋后的計算過程。

首先，為處理裂紋間斷，對單元形函數(shù)進(jìn)行擴(kuò)展，將單元的位移場u(x)描述分解為連續(xù)部分和非連續(xù)部分。連續(xù)部分表示節(jié)點的傳統(tǒng)形函數(shù)部分，非連續(xù)部分表示裂紋附加的自由度，即

(2)

式中：α為單元的節(jié)點編號；A為含裂紋單元，被“虛擬”裂紋分為A+和A-兩個子域；Nα(x)為節(jié)點形函數(shù)；uα為節(jié)點位移；H(x)為Heaviside函數(shù)；N+(x)定義為

(3)

其次，假設(shè)單元內(nèi)部的應(yīng)力等于“虛擬”裂紋兩側(cè)的應(yīng)力，得到內(nèi)聚裂紋模型的基本方程：

(4)

式中：E為彈性模量張量；n為裂紋張開方向的單位法向矢量；εα為根據(jù)節(jié)點位移得到的應(yīng)變張量；上標(biāo)s表示取張量的對稱部分；

(5)

通過迭代求解(4)式，可以得到單元應(yīng)變和裂紋張開位移，進(jìn)而得到單元應(yīng)力。利用LS-DYNA的二次開發(fā)接口，將內(nèi)聚裂紋模型嵌入程序中，通過輸出歷史變量可以給出計算模型的斷裂演化情況。

2 動態(tài)半圓盤點彎曲實驗

動態(tài)半圓盤彎曲實驗在SHPB實驗裝置上進(jìn)行，實驗裝置如圖2所示。子彈、入射桿和透射桿的材料均為LC4鋁合金，直徑為20 mm，子彈長200 mm，入射桿長2 000 mm，透射桿長1 000 mm. 加工圓柱支座安裝于透射桿端部，兩個圓柱支座直徑2 mm，跨度16 mm. 為了使試樣斷裂后可以自由轉(zhuǎn)動，在支座中心位置還預(yù)制了直徑9 mm的半圓凹槽。試樣為PBX模擬材料，加工成含預(yù)制裂紋的φ20 mm×10 mm的半圓盤，厚度10 mm，預(yù)制裂紋寬0.3 mm、長4 mm. 將含預(yù)制裂紋的半圓盤試樣放置于兩圓柱支座和入射桿之間，通過入射桿在試樣的頂面施加載荷，實現(xiàn)對試樣的動態(tài)三點彎曲。

為了保證試樣兩端受力平衡，在入射桿端部加入整形器，延緩加載波的上升前沿。由入射桿上的應(yīng)變片記錄的入射波和反射波信號以及透射桿應(yīng)變片記錄的透射波信號來檢驗試樣在加載過程中應(yīng)力是否平衡。實驗時，采用高速相機(jī)記錄試樣的變形破壞過程，相機(jī)的采集速率為43 200幀/s，圖像分辨率為384像素×144像素。根據(jù)高速攝像結(jié)果，通過DIC方法計算得到半圓盤試樣不同時刻的變形場。實驗后，采用DIC方法獲取試樣破壞前表面的位移場、應(yīng)變場以及位移矢量場。

3 動態(tài)半圓盤彎曲實驗的數(shù)值模擬

3.1 計算模型建立及材料參數(shù)

根據(jù)實驗裝置建立1/2對稱的三維模型，網(wǎng)格尺寸約為0.5 mm，模型示意如圖3所示。由于實驗中采用了波形整形技術(shù)，將實驗得到的應(yīng)變波形換算為軸向應(yīng)力后，作為輸入載荷施加于入射桿端面。將兩個圓柱支座簡化為兩個半圓支座，半圓支座與透射桿之間采用固連接觸。試樣和霍普金森桿為自動面面接觸?；羝战鹕瓧U采用彈性本構(gòu)模型，主要材料參數(shù)為：密度2.77 g/cm3，彈性模量71 GPa，泊松比0.29. PBX模擬材料試樣采用嵌入的內(nèi)聚裂紋模型，主要材料參數(shù)為：密度1.829 g/cm3，體積模量8.71 GPa，剪切模量2.69 GPa，拉伸強(qiáng)度σb為6 MPa，斷裂能Gf為120 J/m2，采用線性軟化曲線，其函數(shù)形式為

(6)

3.2 試樣兩端受力

通過數(shù)值模擬得到的入射桿和透射桿中部的軸向應(yīng)變波形曲線如圖4所示。由圖4可見，入射波與反射波相加后，與透射波基本一致，說明試樣兩端應(yīng)力平衡。通過數(shù)值模擬得到透射桿中部軸向應(yīng)變與實驗結(jié)果的對比如圖5所示。由圖5可見，數(shù)值模擬得到的拉伸應(yīng)變時程與實驗結(jié)果比較吻合，峰值誤差小于5%，持續(xù)時間基本相同，模擬得到拉伸應(yīng)變峰值對應(yīng)時間比實驗結(jié)果提前約9 μs.

3.3 試樣表面的變形場

通過數(shù)值模擬得到開裂前后試樣表面的位移場與實驗結(jié)果的對比如圖6和圖7所示，其中：圖6(a)和圖7(a)、圖6(b)和圖7(b)為數(shù)值模擬結(jié)果；圖6(c)和圖7(c)、圖6(d)和圖7(d)為DIC方法給出的實驗結(jié)果。

由圖6和圖7可以看出：數(shù)值模擬得到的位移場分布與實驗得到的位移場分布較為一致，開裂前后位移場具有較好的對稱性。開裂后，數(shù)值模擬得到的x方向位移場峰值為0.10 mm，集中于試樣預(yù)制裂紋所在區(qū)域，比DIC方法得到的峰值0.08 mm偏大約25%，預(yù)制裂紋兩側(cè)的x方向位移約0.03 mm，與實驗得到的0.04 mm較接近；y方向位移場峰值為0.053 mm，位于半圓盤的底部邊緣，比DIC方法得到的峰值0.080 mm偏小約33%. 注意到，實驗過程中，通過DIC方法得到的x方向位移場對稱性并不理想，其原因可能是在實驗過程中，試樣在兩個圓柱支座位置受力不均勻?？紤]到PBX的非均質(zhì)特性，在壓制過程中產(chǎn)生有大量微觀缺陷，動載下這些缺陷會影響試樣的宏觀響應(yīng)及裂紋的擴(kuò)展，因此二者存在差異是可以理解的。

數(shù)值模擬得到開裂前后試樣表面的拉應(yīng)變場分布如圖8所示，其中：圖8(a)和圖8(b)為數(shù)值模擬結(jié)果；圖8(c)和圖8(d)為DIC方法給出的實驗結(jié)果。

由圖8可以看出：開裂前，在試樣的預(yù)制裂紋尖端沿加載方向出現(xiàn)明顯的局部拉應(yīng)變集中，與y方向位移場的中間分界區(qū)域相重疊，表面應(yīng)變比內(nèi)部應(yīng)變小。開裂前，試樣表面拉應(yīng)變峰值約為3.1%，與實驗得到的拉應(yīng)變峰值3.0%比較吻合；開裂后，試樣表面拉應(yīng)變峰值約為15%，與實驗得到的拉應(yīng)變峰值15%非常吻合。此時，局部拉應(yīng)變集中區(qū)域發(fā)生演化，使得預(yù)制裂紋尖端出現(xiàn)向圓弧頂部擴(kuò)展的裂紋。

3.4 試樣中的裂紋演化過程

由于內(nèi)聚裂紋模型可以得到PBX試樣內(nèi)部裂紋擴(kuò)展隨時間變化的信息，可以進(jìn)一步分析試樣中裂紋的演化過程。數(shù)值模擬得到的試樣中裂紋的無量綱寬度演化如圖9所示，以寬度值0.04 mm對其歸一化。

由圖9可以看出：在0.45 ms時，試樣預(yù)制裂紋尖端出現(xiàn)一條微裂紋，最大裂紋寬度約0.002 9 mm，由預(yù)制裂紋尖端沿半徑向頂部加載方向發(fā)展，即數(shù)值模擬得到的起裂時間約為0.45 ms；在0.02 ms之后即0.47 ms時，微裂紋向入射桿方向擴(kuò)展2.7 mm，最大裂紋寬度為0.038 0 mm；隨著傳入應(yīng)力波的進(jìn)一步加載，到0.70 ms時裂紋已貫穿試樣，試樣的張開角度約17.0°；在0.88 ms時試樣張開角度約26.3°；根據(jù)高速攝像，試樣沿加載方向出現(xiàn)肉眼可見裂紋的時間為0.47 ms，在0.70 ms時試樣張開角度約16.9°，在0.89 ms時試樣張開角度約27.2°. 上述不同時刻PBX試樣的張開角度如圖10所示，表明計算結(jié)果與實驗結(jié)果比較一致。

圖11給出了模擬得到的預(yù)制裂紋尖端位置單元的y方向應(yīng)力和裂紋寬度變化。

由圖11可以看出，在0.44 ms時，應(yīng)力波傳播到試樣預(yù)制裂紋尖端，該位置單元的y方向應(yīng)力達(dá)到峰值6 MPa，與給定的拉伸強(qiáng)度一致，隨后試樣起裂。開裂后，裂紋寬度快速增加、豎直應(yīng)力快速下降。在該下降階段，數(shù)值模擬得到的裂紋寬度比通過DIC方法得到的結(jié)果偏小約15%，當(dāng)應(yīng)力降為0后，數(shù)值模擬結(jié)果與實驗結(jié)果的偏差逐漸增大。

綜合來看，內(nèi)聚裂紋模型可以較好地模擬PBX在動態(tài)半圓盤彎曲實驗過程中的裂紋擴(kuò)展過程，表明應(yīng)用該模型分析PBX的拉伸破壞行為是可行的。在動態(tài)加載下，PBX試樣呈現(xiàn)出與準(zhǔn)靜態(tài)實驗類似的劈裂樣式，內(nèi)聚裂紋模型較好地再現(xiàn)了PBX試樣的主要失效機(jī)制，即沿加載直徑的拉伸破壞，最終裂紋沿直徑方向擴(kuò)展貫穿試樣。

4 結(jié)論

本文基于內(nèi)聚裂紋模型，對PBX模擬材料制成的帶預(yù)制裂紋半圓盤試樣在沖擊作用下的動態(tài)拉伸變形和破壞過程進(jìn)行了數(shù)值模擬，通過對比實驗測量得到的應(yīng)變結(jié)果、試樣表面位移場、應(yīng)變場，以及不同時刻試樣張開角度的高速攝影結(jié)果，驗證了內(nèi)聚裂紋模型的可靠性，表明：1)內(nèi)聚裂紋模型適用于模擬分析PBX材料的動態(tài)拉伸破壞行為；2)動態(tài)加載下，PBX試樣的主要失效機(jī)制為與加載方向相關(guān)的拉伸破壞，主要失效模式為內(nèi)部缺陷位置起裂，主裂紋貫穿破壞。

)

[1] 陳鵬萬，黃風(fēng)雷. 含能材料損傷理論及應(yīng)用[M]. 北京: 北京理工大學(xué)出版社, 2006.

CHEN Peng-wan, HUANG Feng-lei. Damage theory and application of energetic materials[M]. Beijing: Beijing Institute of Technology Press, 2006. (in Chinese)

[2] Grantham S G, Siviour C R, Proud W G, et al. High-strain rate Brazilian testing of an explosive simulant using speckle metrology[J]. Measurement Science and Technology, 2004, 15(9): 1867-1870.

[3] Zhou Y X, Xia K, Li X B, et al. Suggested methods for determining the dynamic strength parameters and mode-I fracture toughness of rock materials[J]. International Journal of Rock Mechanics & Mining Sciences, 2012, 49:105-112.

[4] Zhou Z B, Chen P W, Duan Z B. Comparative study of the fracture toughness determination of a polymer-bonded explosive stimulant [J]. Engineering Fracture Mechanics, 2011,78(17): 2991-2997.

[5] Chen R, Cheng L, Lin Y, et al. Studies on the dynamic fracture properties and failure modes of a PBX[J]. International Journal of Applied Mechanics， 2014, 6(4):1450039.

[6] Dienes J K, Zuo Q H, Kershner J D. Impact initiation of explosives and propellants via statistical crack mechanics[J]. Journal of the Mechanics and Physics of Solids, 2006, 54(6):1237-1275.

[7] 傅華，李俊玲，譚多望. PBX炸藥動態(tài)Brazilian試驗及數(shù)值模擬研究[J]. 高壓物理學(xué)報, 2012, 26(2): 148-154.

FU Hua, LI Jun-ling, TAN Duo-wang. Dynamic Brazilian test and simulation of plastic-bonded explosives [J]. Chinese Journal of High Pressure Physics, 2012, 26(2): 148-154. (in Chinese)

[8] 王鵬飛，黃西成，何穎波. 基于線性Drucker-Prager模型的PBX準(zhǔn)靜態(tài)彈塑性變形分析[J]. 含能材料, 2015, 23(5):464-471.

WANG Peng-fei, HUANG Xi-cheng, HE Ying-bo. Quasi-static elastoplastic deformation analysis of PBX based on linear Drucker-Prager model[J]. Chinese Journal of Energetic Materials, 2015, 23(5):464-471.(in Chinese)

[9] 謝鳳英，張丘，劉維，等. 高聚物粘結(jié)炸藥銑削時的邊緣崩塊形成機(jī)理分析[J]. 兵工學(xué)報, 2016, 37(5): 823-828.

XIE Feng-ying, ZHANG Qiu, LIU Wei, et al. Preliminary research on formation mechanisms of edge spalling in PBX milling [J]. Acta Armamentarii, 2016, 37(5): 823-828.(in Chinese)

[10] Planas J, Elices M, Guinea G V, et al. Generalizations and specializations of cohesive crack models [J]. Engineering Fracture Mechanics, 2003, 70(14): 1759-1776.

[11] Reyes E, Gálvez J C, Casati M J, et al. An embedded cohesive crack model for finite element analysis of brickwork masonry fracture [J]. Engineering Fracture Mechanics, 2009, 76(12): 1930-1944.

[12] 崔云霄，陳鵬萬，David A. Cendón，等.基于內(nèi)聚裂紋模型的高聚物粘結(jié)炸藥模擬材料動態(tài)巴西實驗的數(shù)值模擬[J]. 兵工學(xué)報, 2016,37 (9): 1639-1645.

CUI Yun-xiao, CHEN Peng-wan, Cendón D A, et al. Numerical simulation of dynamic Brazilian test of polymer bonded explosive simulant based on cohesive crack model[J]. Acta Armamentarii, 2016, 37 (9): 1639-1645. (in Chinese)

ResearchonDynamicFractureBehaviorofPolymerBondedExplosiveSimulantBasedonCohesiveCrackModel

CUI Yun-xiao1, 2, CHEN Peng-wan1, GUO Bao-qiao1, David A. Cendón3, ZHOU Zhong-bin1

(1.State Key Laboratory of Explosion Science and Technology，Beijing Institute of Technology, Beijing 100081, China; 2.Northwest Institute of Nuclear Technology, Xi’an 710024, Shaanxi，China; 3.Technical University of Madrid，Madrid 28040, Spain)

The dynamic tensile fracture behavior of polymer bonded explosive simulant was investigated. The experimental test was performed on notched semi-circular bending (NSCB) samples of the material by using the split Hopkinson pressure bar(SHPB) apparatus. With the aide of a high-speed camera, the dynamic displacement and strain field during the dynamic fracture process were obtained by the digital image correlation (DIC) technique. Based on cohesive crack model (CCM), the numerical simulation of the same conditions was performed. Compared with the results obtained by simulation and experiment, the good agreement was found for the time history of tensile stress and the displacement and strain fields during dynamic fracture. The behavior of crack propagation of PBX during NSCB test was analyzed by the CCM. Furthermore, the width of crack in NSCB sample is about 15% smaller than the experimental results.

explosion mechanics; cohesive crack model; dynamic NSCB experiment; polymer bonded explosive; technique of digital image correlation

O347.3; TQ560.71

A

1000-1093(2017)12-2379-07

10.3969/j.issn.1000-1093.2017.12.011

2017-01-10

國家自然科學(xué)基金項目(11202027、11521062、U1330202)

崔云霄(1980—)，男，博士研究生。E-mail: yunxiaocui@163.com

陳鵬萬(1971—)，男，教授，博士生導(dǎo)師。E-mail: pwchen@bit.edu.cn