SHPB測試中巖石試件應(yīng)力均勻性分析研究

劉 鑫， 郭連軍， 張大寧， 徐振洋， 郭 浩

(1. 遼寧科技大學(xué) a. 礦業(yè)工程學(xué)院， b. 土木工程學(xué)院, 遼寧鞍山 114051；2. 遼寧省冶金地質(zhì)勘察局四零三隊， 遼寧鞍山 114021)

?

SHPB測試中巖石試件應(yīng)力均勻性分析研究

劉 鑫1a， 郭連軍1b， 張大寧1a， 徐振洋1a， 郭 浩2

(1. 遼寧科技大學(xué) a. 礦業(yè)工程學(xué)院， b. 土木工程學(xué)院, 遼寧鞍山 114051；2. 遼寧省冶金地質(zhì)勘察局四零三隊， 遼寧鞍山 114021)

摘要：對花崗巖、千枚巖、磁鐵石英巖三種巖石進(jìn)行SHPB應(yīng)力均勻性試驗測試，研究表明：試件內(nèi)達(dá)到應(yīng)力均勻時間介于彈塑性狀態(tài)假設(shè)條件下試件內(nèi)達(dá)到應(yīng)力均勻時間之間，且試件內(nèi)達(dá)到應(yīng)力均勻時間與二者的平均值具有一定的相關(guān)性。試驗中子彈的沖擊速度在一定范圍內(nèi)，能夠使三種巖石滿足應(yīng)力均勻性的要求，且在斷裂應(yīng)變達(dá)到1%之前均有較長的恒應(yīng)變率時間段。

關(guān)鍵詞：SHPB； 應(yīng)力均勻性； 沖擊速度； 巖性； 斷裂應(yīng)變； 恒應(yīng)變率

1引 言

分離式霍普金森壓桿自Kolsky〔1〕發(fā)明以來，并在材料動態(tài)力學(xué)行為研究方面不斷改進(jìn)的基礎(chǔ)上，現(xiàn)已成為研究材料在高應(yīng)變率條件下動態(tài)特性的一種試驗裝置。SHPB試驗裝置是以兩個基本假定〔2〕為前提的：壓桿的一維應(yīng)立波假定和試件應(yīng)力均勻性假定。近年來，SHPB試驗技術(shù)被應(yīng)用于爆破等工程領(lǐng)域〔3-4〕，但由于巖石類脆性材料不同于大多數(shù)金屬類材料，巖石試件可能在加載結(jié)束之前就發(fā)生了斷裂破壞，所以評估巖石類脆性材料滿足試驗的應(yīng)力均勻性假定，是利用SHPB試驗裝置研究巖石類脆性材料動態(tài)行為的前提。

本文利用φ50mm的SHPB裝置系統(tǒng)，通過應(yīng)力均勻性控制理論計算出巖石試件達(dá)到應(yīng)力均勻性的時間，從而分析研究花崗巖、千枚巖、磁鐵石英巖三種巖石試件的動態(tài)沖擊試驗達(dá)到應(yīng)力均勻性時間的差異以及運(yùn)用彈塑性狀態(tài)假設(shè)，計算達(dá)到應(yīng)力均勻性時間的均值與達(dá)到應(yīng)力均勻性時間的關(guān)系。

2應(yīng)力均勻性控制理論

2.1彈性狀態(tài)假設(shè)

SHPB試驗中，彈性應(yīng)力波在不同波阻抗的壓桿和試件中傳播時，將會在壓桿和試件交界面上發(fā)生透、反射。若令試件兩個端面的應(yīng)力差為αk，彈性應(yīng)力波在試件與壓桿界面間的透、反射次數(shù)為k，試件和壓桿的波阻抗比β=(ρScSAS)/(ρBcBAB)，其中ρS，cS，AS和ρB，cB，AB分別為試件和壓桿材料的密度、彈性縱波波速和截面積。

在彈性狀態(tài)假設(shè)條件下〔5-6〕，應(yīng)力差αk、應(yīng)力波透、反射次數(shù)k、波阻抗比β三者可以建立如下關(guān)系：

(1)

一般可近似認(rèn)為αk≤5%時，試件達(dá)到了應(yīng)力均勻假設(shè)的要求。當(dāng)試件和壓桿的材料確定以后，波阻抗比值β就是已知的，由此可以得到試樣內(nèi)應(yīng)力均勻程度αk和應(yīng)力波透射和反射次數(shù)k的關(guān)系式，從而描繪出αk和k的曲線圖。

參照文獻(xiàn)〔7-10〕 的方法，利用αk與k關(guān)系曲線確定k值大小，則在彈性狀態(tài)假設(shè)條件下試件內(nèi)達(dá)到應(yīng)力均勻性要求的時間te為：

te=kLS/c

(2)

式中：LS為試件厚度;c是彈性波波速。

2.2塑性狀態(tài)假設(shè)

一般在波傳播過程中彈性先驅(qū)波后面跟著塑性波，若假設(shè)試件是在塑性狀態(tài)下達(dá)到應(yīng)力均勻性要求，試件達(dá)到應(yīng)力均勻性要求時塑性波需要在試件內(nèi)傳播π次〔11〕，則在塑性狀態(tài)假設(shè)條件下試件內(nèi)達(dá)到應(yīng)力均勻性要求的時間tp為：

(3)

2.3應(yīng)力均勻時間

由霍普金森試驗裝置系統(tǒng)的一維應(yīng)力波假設(shè)可知，試件左側(cè)端面的壓力σL(t)和試件右端面的壓力σR(t)為：

(4)

(5)

式中：EB是壓桿的彈性模量；εi(t)是實測的入射波應(yīng)變；εr(t)實測的反射波應(yīng)變；εt(t)實測的透射波應(yīng)變。

Ravichandran等〔7〕定義試件兩個端面應(yīng)力的差值與其平均值的比值α(t)為應(yīng)力均勻度來度量試件的應(yīng)力均勻性，則：

(6)

一般α(t)-t曲線上存在較明顯的轉(zhuǎn)折點〔7〕，轉(zhuǎn)折點之后的曲線上應(yīng)力不均勻度數(shù)值小于轉(zhuǎn)折點的應(yīng)力不均勻度數(shù)值，并認(rèn)為轉(zhuǎn)折點的應(yīng)力不均勻度數(shù)值所對應(yīng)的時間t是滿足應(yīng)力均勻性要求的時間。

3試驗分析

3.1試驗應(yīng)力均勻性分析

SHPB沖擊試驗時間應(yīng)變率一般為102s-1～104s-1，巖石類脆性材料總體可以認(rèn)為經(jīng)歷了4個階段，彈性變形階段、非線性變形階段、很小的塑性變形階段(或者不存在)和破壞階段。試件內(nèi)經(jīng)歷了復(fù)雜的應(yīng)力波傳播過程，按彈性狀態(tài)假設(shè)計算的應(yīng)力均勻時間是以理想的矩形波為基礎(chǔ)，所以按彈性狀態(tài)假設(shè)計算的應(yīng)力均勻時間要比實際試驗系統(tǒng)中試件達(dá)到應(yīng)力均勻時間的值小。反之，按塑性狀態(tài)假設(shè)計算的應(yīng)力均勻時間要比實際試驗系統(tǒng)中試件達(dá)到應(yīng)力均勻時間的值大。那么，以彈塑性波傳播為基礎(chǔ)，根據(jù)上述假設(shè)條件，試件達(dá)到應(yīng)力均勻性要求的時間t是介于te和tp之間，并與試件開始破壞時間tk比較分析。

3.2實際試驗系統(tǒng)存在的問題及偏差

SHPB應(yīng)力均勻性試驗分析中忽略試件與壓桿波阻抗比值變化的影響；考慮試件的加工難易程度和最佳的長徑比，選取試件的直徑為50mm，厚度為25mm；入射脈沖波形整形技術(shù)一般采用厚度為1mm、直徑為10mm的黃銅波形整形器。所以，文中忽略了試件與壓桿的波阻抗比值、試件與壓桿的面積比、試件厚度和波形整形器尺寸對試驗應(yīng)力均勻性的影響。

4 SHPB試驗

4.1沖擊壓縮試驗及數(shù)據(jù)

選用φ50mm的SHPB裝置系統(tǒng)，試驗用的花崗巖、千枚巖、磁鐵石英巖三種巖石材料取自齊大山鐵礦，密度分別是2.56g/cm3、2.72g/cm3、3.63g/cm3，對按要求加工好的三大組試件分別進(jìn)行不同速度的沖擊試驗，試驗數(shù)據(jù)見表1，試驗裝置示意圖如圖1所示。

4.2試驗數(shù)據(jù)分析

試件達(dá)到應(yīng)力均勻性要求的時間t與按彈性狀態(tài)假設(shè)條件下和塑性狀態(tài)假設(shè)條件下計算的應(yīng)力均勻時間的平均值(te+ tp)/2的關(guān)系如圖2所示，擬合的5組曲線函數(shù)方程見表2。

表1　三種巖石的試驗數(shù)據(jù)

1-子彈Bullet;2-入射桿The incident bar;3-透射桿The transmission bar;4-吸收桿Aborbing bar;5-緩沖桿Bumper bar;6-平行光源Source of parallel light; 7-波形整形器Waveform shaper;8-應(yīng)變片Strain gauge;9-試件Specimen;10-激光測速儀Laser velocimetry;11-超動態(tài)應(yīng)變儀High dynamic strain indicator;12-波形存儲器Wave memorizer圖1　試驗裝置示意圖Fig.1　Schematic diagram of experimental apparatus

圖2　t與(te+ tp)/2的關(guān)系Fig.2　The relationship between t and (te+ tp)/2

對比R2可知，t與(te+ tp)/2是多項式關(guān)系，具有一定的相關(guān)性，其表達(dá)式為：

t=0.0028[(te+ tp)/2]2+0.5062[(te+ tp)/2]+9.0244

(7)

表2　(te+ tp)/2與t的擬合曲線方程

注：令tep=(te+tp)/2

圖3表明三種巖石的SHPB試驗中，子彈分別對花崗巖、千枚巖、磁鐵石英巖的沖擊速度在11.5m/s～14.0m/s、10.9m/s～14.3m/s、9.6m/s～14.1m/s范圍內(nèi)，SHPB試驗中巖石試件在破裂應(yīng)變達(dá)到1%之前，應(yīng)變時間曲線有較長一段時間是線性的，試件處于恒應(yīng)變率狀態(tài)，即達(dá)到應(yīng)力均勻性的要求。

圖3　三種巖石滿足應(yīng)力均勻性要求的應(yīng)變時間曲線Fig.3　Strain-time curves of rocks meeting the stress uniformity requirements

花崗巖、千枚巖、磁鐵石英巖在沖擊速度分別是8.7m/s、6.7m/s、7.3m/s時都沒有達(dá)到應(yīng)力均勻性要求，原因可能是沖擊速度太小，使試件產(chǎn)生較大的振蕩，導(dǎo)致試件破裂或有明顯的宏觀裂紋，隨后應(yīng)力波繼續(xù)作用在試件上使試件破碎成更多的碎塊。如圖4所示，隨著時間的增加，應(yīng)力時間曲線先是出現(xiàn)一個較小的峰值，隨后曲線迅速下降，緊接著上升直到出現(xiàn)更大的應(yīng)力峰值?；◢弾r、千枚巖、磁鐵石英巖在沖擊速度分別是17.6m/s、16.8m/s、16.9m/s時也沒有達(dá)到應(yīng)力均勻性要求，如圖5所示，SHPB試驗中三種巖石試件在破裂應(yīng)變達(dá)到1%之前，不存在較長時間的恒應(yīng)變率狀態(tài)，因為子彈對三種試件的沖擊速度太大，體現(xiàn)了脆性材料在高應(yīng)變率條件下易斷裂的特點。

圖4　較小沖擊速度三種巖石應(yīng)力時間曲線Fig.4　Stress-time curves of rocks at smaller impact velocity in three kinds of rocks

圖5　較大沖擊速度三種巖石應(yīng)變時間曲線Fig.5　Strain-time curves of rocks at larger impact velocity in three kinds of rocks

5結(jié) 論

(1)SHPB試驗中子彈的沖擊速度在一定范圍內(nèi)，能夠使三種巖石滿足應(yīng)力均勻性的要求，且在斷裂應(yīng)變達(dá)到1%之前均有較長的恒應(yīng)變率時間段。

(2)試件內(nèi)達(dá)到應(yīng)力均勻時間要比彈性狀態(tài)假設(shè)條件下試件內(nèi)達(dá)到應(yīng)力均勻時間長，且試件內(nèi)達(dá)到應(yīng)力均勻時間與彈性狀態(tài)假設(shè)條件下和塑性狀態(tài)假設(shè)條件下試件內(nèi)達(dá)到應(yīng)力均勻時間的平均值具有一定的相關(guān)性，其表達(dá)式為t=0.0028[(te+ tp)/2]2+0.5062[(te+ tp)/2]+9.0244。

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

〔1〕 KOLSKY H. An investigation of the mechanical properties of materials at very high rates of loading[J].Proceedings of the Physical Society B, 1949, 62(11): 676-700.

〔2〕 王禮立.應(yīng)力波基礎(chǔ)[M]. 2版. 北京:國防工業(yè)出版社,2010:5-64.

WANG Li-li. Foundations of stress waves[M]. 2nd ed. Beijing: National Defense Industry Press, 2010:5-64.

〔3〕 郭連軍,楊躍輝,華悅含,等. 沖擊荷載作用下花崗巖動力特性試驗分析[J]. 工程爆破, 2014, 20(1): 1-4.

GUO Lian-jun, YANG Yue-hui, HUA Yue-han, et al. Experimental study of dynamic characteristics of granite under impact loading[J]. Engineering Blasting, 2014, 20(1): 1-4.

〔4〕 梁書鋒,吳帥峰,李勝林,等. 巖石材料SHPB實驗試件尺寸確定的研究[J]. 工程爆破, 2015, 21(5):1-5.

LIANG Shu-feng, WU Shuai-feng, LI Sheng-lin, et al. Study on the determination of specimen size in SHPB experiments of rock materials[J]. Engineering Blasting, 2015, 21(5):1-5.

〔5〕 WANG Li-li. Influences of stress propagation upon studying dynamic response of materials at high strain rates[J]. Journal of Beijing Institute of Technology, 2004, 13(3): 225-235.

〔6〕 王禮立,王永剛. 應(yīng)力波在用SHPB研究材料動態(tài)本構(gòu)特性中的重要作用[J]. 爆炸與沖擊,2005,25(1):17-25.

WANG Li-li, WANG Yong-gang.The important role of stress waves in the study on dynamic constitutive behavior of materials by SHPB[J]. Explosion and Shock Waves，2005,25(1):17-25.

〔7〕 RAVICHANDRAN G, SUBHASH G. Critical apprasisal of limiting strain rates for compression testing of ceramics in a split Hopkinson bar[J]. Journal of American Ceramic Society, 1994, 77(1): 263-267.

〔8〕 毛勇建,李玉龍. SHPB試驗中試件的軸向應(yīng)力均勻性[J]. 爆炸與沖擊,2008,28(9):448-454.

MAO Yong-jian, LI Yu-long. Axial stress uniformity in specimens of SHPB tests[J]. Explosion and Shock Waves, 2008,28(9):448-454.

〔9〕 宋力，胡時勝. SHPB測試中的均勻性問題及恒應(yīng)變率[J]. 爆炸與沖擊，2005,25(3):207-216.

SONG Li, HU Shi-sheng. Stress uniformity and constant strain rate in SHPB test[J]. Explosion and Shock Waves，2005,25(3):207-216.

〔10〕 劉瑞堂,果春煥,張智峰. 分離式Hopkinson壓桿系統(tǒng)沖擊壓縮過程中試樣的應(yīng)力均勻化過程[J]. 機(jī)械工程材料,2009,33(2):25-27.

LIU Rui-tang, GUO Chun-huan, ZHANG Zhi-feng.Stress uniformization process in specimens during compression process with split Hopkinson pressure bar system[J]. Materials for Mechanical Engineering, 2009,33(2):25-27.

〔11〕 DAVIES E D H, HUNTER S C. The dynamic compression testing of solids by the method of the split Hopkinson pressure bar[J]. J. Mech. Phys. Solids, 1963(11):115-179.

Experimental study for stress uniformity of rock specimen in SHPB test

LIU Xin1a, GUO Lian-jun1b, ZHANG Da-ning1a, XU Zhen-yang1a, GUO Hao2

(1. a. School of Mine Engineering, b. School of Civil Engineering,University of Science and Technology Liaoning, Anshan 114051, Liaoning， China；2. The 403 Team of Metallurgical Geology Exploration Bureau of Liaoning, Anshan 114021, Liaoning， China)

ABSTRACT：SHPB device was used to test stress uniformity of granite, phyllite and magnetite. The results indicated that the specimen stress uniformity time was between the time of the specimens reaching stress uniformity under elastic-plastic state assumption, and the specimen stress uniformity time had some relevance with the average of them. In a certain range of the bullet impact velocity, the SHPB test could satisfy the requirement of stress uniformity of the three kinds of rock specimens, and there was a long time of constant strain rate before the fracture strain reached 1%.

KEY WORDS:SHPB; Stress uniformity; Impact velocity; Rock properties; Fracture strain; Constant strain rate

文章編號：1006-7051(2016)02-0024-04

收稿日期：2015-12-20

基金項目：國家自然科學(xué)基金面上項目(51474123)

作者簡介：劉 鑫(1991-)，男，碩士，主要研究爆破優(yōu)化與設(shè)計。E-mail：2232199160@qq.com通訊作者：郭連軍(1963-)，男，博士、教授，從事爆破工程和巖土工程研究與教學(xué)。E-mail：glj0412@126.com

中圖分類號：TD235； TU45

文獻(xiàn)標(biāo)識碼：A

doi：10.3969/j.issn.1006-7051.2016.02.005