江紅祥，杜長龍，劉送永，高魁東

?

高壓水射流沖擊破巖損傷場分析

江紅祥，杜長龍，劉送永，高魁東

(中國礦業(yè)大學(xué) 機(jī)電工程學(xué)院，江蘇 徐州，221116)

以Johnson-Holmquist modelⅡ(JH-2)脆性材料本構(gòu)模型來表征巖石的力學(xué)特性，利用光滑粒子流(SPH)和有限元(FEA)耦合算法建立高壓水射流沖擊破巖過程的數(shù)值模型，很好地模擬水射流損傷破巖過程中巖石失效、裂紋擴(kuò)展以及不同位置巖石單元的損傷程度隨時(shí)間變化的過程。此外，模擬分析射流沖擊速度、入射角以及平移速度對(duì)巖石損傷場的影響。研究結(jié)果表明：巖石的損傷破壞、裂紋擴(kuò)展是剪切和拉伸共同作用引起的，且近射流沖擊點(diǎn)的損傷破壞由剪切作用主導(dǎo)，而放射性裂紋和層狀裂紋的擴(kuò)展主要由拉伸作用主導(dǎo)；巖石單元的損傷值隨時(shí)間呈階躍性變化，且射流沖擊損傷破巖過程為微秒量級(jí)；射流沖擊損傷破巖的入射角存在一個(gè)有效范圍，當(dāng)入射角大于70o時(shí)損傷破巖效果較好；射流沖擊速度較低時(shí)以表面沖蝕巖石為主，射流沖擊速度達(dá)到一定值才能使巖石損傷破壞，出現(xiàn)放射性裂紋和層狀裂紋；射流平移速度與沖擊速度相比很小，對(duì)巖石的損傷場影響不明顯。

水射流；光滑粒子流；有限元；巖石損傷破壞；JH-2本構(gòu)模型

高壓水射流破巖技術(shù)已經(jīng)在石油鉆井、地下工程以及瓦斯抽放等領(lǐng)域得到了很好的應(yīng)用[1?4]。高壓水射流破巖的機(jī)理和實(shí)際物理過程非常復(fù)雜，其損傷破巖機(jī)理的研究一直都是熱點(diǎn)和難點(diǎn)問題[5?6]。到目前為止，在理論和實(shí)驗(yàn)研究基礎(chǔ)上主要形成了汽蝕破壞、沖擊破壞、動(dòng)壓破壞以及水楔破壞等理論學(xué)說。高壓水射流損傷破巖過程是集流體力學(xué)、巖石力學(xué)、沖擊動(dòng)力學(xué)及流固耦合等的多學(xué)科問題，其損傷破巖機(jī)理至今未形成較統(tǒng)一的學(xué)說[7]。此外，由于巖石的透明性差，難以通過實(shí)驗(yàn)方法來觀察、研究巖石的損傷破壞過程，也不能確定是拉應(yīng)力、剪應(yīng)力或壓應(yīng)力引起的損傷破壞。目前已有很多學(xué)者對(duì)高壓水射流沖擊破巖過程進(jìn)行數(shù)值模擬研究：廖華林等[8]利用標(biāo)準(zhǔn)?雙方程模型、體積控制法、各向同性彈性力學(xué)以及有限元法建立了水射流和巖石耦合數(shù)值模型，模擬研究了不同沖擊速度條件下水射流的動(dòng)力學(xué)特性以及巖石內(nèi)部的應(yīng)力分布規(guī)律；Ni等[9]利用有限元、損傷力學(xué)以及細(xì)觀力學(xué)等理論建立了水射流破巖的流固耦合模型，模擬表明射流沖擊破巖存在拉伸和卸載過程；孫清德等[10]利用動(dòng)態(tài)有限元法模擬了高壓水射流破巖過程，認(rèn)為水射流破巖的2個(gè)臨界壓力直接影響破巖機(jī)理，并分析了射流沖擊速度、橫移速度以及入射角等對(duì)破巖效率的影響規(guī)律；司鵠等[11]利用非線性動(dòng)態(tài)有限元法對(duì)水射流破巖過程中巖石應(yīng)力波效應(yīng)進(jìn)行了數(shù)值分析，表明巖石受到水射流沖擊時(shí)的不穩(wěn)定性是一個(gè)由強(qiáng)變?nèi)醯倪^程；劉佳亮等[7]利用拉格朗日?歐拉耦合算法對(duì)水射流沖擊有圍壓巖石的損傷破壞過程進(jìn)行了模擬，認(rèn)為圍壓對(duì)巖石軸向的損傷演化有較大的影響；盧義玉等[12]利用光滑粒子動(dòng)力學(xué)方法模擬了脈沖射流沖擊不同類型巖石的應(yīng)力波效應(yīng)，且應(yīng)力波導(dǎo)致不同類型巖石的破壞形式不同。以上研究在一定程度上推動(dòng)了水射流破巖機(jī)理的發(fā)展，但并未深入研究水射流損傷破巖過程以及巖石損傷場的演化和裂紋擴(kuò)展過程。為此，本文作者采用SPH/FEA耦合方法建立了高壓水射流沖擊損傷破壞巖石的過程，對(duì)巖石的損傷場演化和裂紋擴(kuò)展過程進(jìn)行模擬研究，并對(duì)影響巖石損傷場的因素進(jìn)行數(shù)值分析，以期提高對(duì)高壓水射流沖擊損傷破巖過程的認(rèn)識(shí)。

1 水射流沖擊損傷破壞巖石數(shù)值模型

1.1 SPH/FEA耦合算法

高壓水射流沖擊損傷破巖過程中，高速射流出現(xiàn)高壓和大變形問題，采用傳統(tǒng)的拉格朗日法模擬水射流容易出現(xiàn)網(wǎng)格畸變而導(dǎo)致計(jì)算終止。采用歐拉?拉格朗日耦合算法雖然可以避免網(wǎng)格的畸變，但其需要消耗更高的計(jì)算成本[13]。SPH算法是近年發(fā)展起來的一種無網(wǎng)格算法，屬于拉格朗日范疇，它是用粒子單元代替有限元網(wǎng)格，并基于空間函數(shù)和核函數(shù)將方程離散。SPH/FEA耦合算法可以很好地模擬不連續(xù)、大變形等問題，還可以克服傳統(tǒng)有限元以及歐拉法存在的網(wǎng)格畸變、計(jì)算耗時(shí)以及資源占用過多等問題[14?15]。

1.1.1 SPH算法

與有限單元法相比，SPH方法的粒子單元是離散的，采用光滑長度內(nèi)粒子代替了有限元法中的節(jié)點(diǎn)，每個(gè)粒子周圍光滑長度內(nèi)的粒子數(shù)量和分布是不確定的。對(duì)于任意連續(xù)光滑場函數(shù)，用()來近似某一點(diǎn)的場函數(shù)值，()可表述為

式中：為光滑核函數(shù)，目前最常用的光滑核函數(shù)是三次B樣條曲線函數(shù)[16]；為光滑長度；和’為空間不同位置點(diǎn)的向量。

在式(1)的基礎(chǔ)上，利用散度定理即可得到場函數(shù)的空間導(dǎo)數(shù)，進(jìn)而離散化可得到

式中：m為第個(gè)SPH粒子的質(zhì)量；為第個(gè)粒子的密度；為光滑長度范圍內(nèi)的粒子數(shù)。

通過上述過程的處理，即可用光滑核函數(shù)的場函數(shù)來近似光滑場函數(shù)的空間導(dǎo)數(shù)項(xiàng)，進(jìn)而實(shí)現(xiàn)基本方程的離散，SPH方法離散后的Navier-Stocks方程：

1.1.2 耦合邊界處理

在水射流沖擊損傷破巖過程中，固體介質(zhì)(巖石)的應(yīng)變率要比液體介質(zhì)(水)小得多。此外，由于SPH粒子單元過多會(huì)導(dǎo)致計(jì)算機(jī)所需要的內(nèi)存過大、計(jì)算時(shí)間過長，故本文的巖石介質(zhì)的應(yīng)力、應(yīng)變等狀態(tài)參數(shù)通過有限元法來求解。由于巖石內(nèi)孔隙流體與巖石耦合作用對(duì)巖石損傷破壞的影響很小[17]，可忽略不計(jì)，故本文利用節(jié)點(diǎn)?面接觸算法定義射流與巖石界面之間的耦合：SPH粒子視作節(jié)點(diǎn)單元，控制參數(shù)為節(jié)點(diǎn)編號(hào)、質(zhì)量以及空間位置，其定義為從節(jié)點(diǎn)；有限元法描述的巖石介質(zhì)部分定義為主面。射流和巖石的應(yīng)力、應(yīng)變等狀態(tài)參數(shù)同時(shí)求解，通過節(jié)點(diǎn)?面接觸算法實(shí)現(xiàn)力的傳遞，且射流與巖石之間滿足滑移條件：

式中：為射流和巖石接觸面SPH粒子速度；為射流和巖石接觸區(qū)域固體單元的速度。

1.2 材料模型

1.2.1 水射流材料模型

水射流沖擊損傷破巖過程中的狀態(tài)通過多項(xiàng)式狀態(tài)方程描述。當(dāng)水處于壓縮狀態(tài)時(shí)(≥0)，狀態(tài)方程為

當(dāng)水處于膨脹狀態(tài)時(shí)(≤0)，狀態(tài)方程為

式中：w為水壓力；w為水的密度；為水壓縮比；為水的內(nèi)能；1，2，3，1，2，0以及1為水介質(zhì)材料常數(shù)，相關(guān)參數(shù)如表1所示。

表1 水狀態(tài)方程相關(guān)參數(shù)

1.2.2 巖石材料模型

巖石材料的力學(xué)特性采用JH-2本構(gòu)模型描述，該模型可以很好的模擬材料的大變形、高應(yīng)變率以及高壓效應(yīng)。材料體積應(yīng)變與靜水壓力之間的關(guān)系用狀態(tài)方程來描述，在低速?zèng)_擊情況下靜水壓力張量較小。由于本文側(cè)重分析射流沖擊損傷破巖，而確定損傷度時(shí)需要準(zhǔn)確的高壓狀態(tài)方程，因此，該情況下巖石采用3次多項(xiàng)式狀態(tài)方程來描述。未受到?jīng)_擊損傷時(shí)，JH-2本構(gòu)模型中巖石的狀態(tài)方程[18]：

式中：r為靜水壓力；1為巖石的體積模量；2和3為巖石材料常數(shù)；為巖石的體積應(yīng)變。

JH-2強(qiáng)度模型是在引入應(yīng)變率和損傷因子基礎(chǔ)上，將材料的等效應(yīng)力表示成靜水壓力的冪函數(shù)形式，其中規(guī)范化強(qiáng)度模型為

當(dāng)巖石的損傷因子=0時(shí)，規(guī)范化等效應(yīng)力為

當(dāng)巖石的損傷因子=1時(shí)巖石完全失效，規(guī)范化等效應(yīng)力為

JH-2模型中的損傷因子可表示為：

采用改進(jìn)的最大主應(yīng)力模型[20]來描述巖石在水射流沖擊下的失效行為，當(dāng)單元的最大主應(yīng)力超過材料的抗拉強(qiáng)度或抗剪強(qiáng)度時(shí)單元失效破壞：

在水射流沖擊損傷破巖模型中，采用花崗巖作為巖石材料，因?yàn)樯鲜龃蟛糠謪?shù)都可以直接或間接獲得，而且產(chǎn)地對(duì)花崗巖力學(xué)特性影響不大，花崗巖的JH-2模型參數(shù)如表2所示[22]。

表2 花崗巖的JH-2模型參數(shù)

1.3 水射流沖擊損傷破巖幾何模型

水射流沖擊損傷破巖幾何模型如圖1所示，水射流簡化為20 mm×2 mm的矩形水束，共640個(gè)SPH粒子。

圖1 幾何模型

巖石模型的長與高分別為100 mm和30 mm，劃分出48 000個(gè)有限單元。為消除邊界對(duì)損傷場的影響，除巖石上表面外，其他都施加無反射邊界，并在巖石內(nèi)部定義數(shù)個(gè)測(cè)量點(diǎn)(gauges)以測(cè)量巖石內(nèi)部的應(yīng)力、損傷因子等。

2 水射流損傷破壞巖石過程演化

為研究花崗巖在射流沖擊損傷破壞過程，利用上述建立的數(shù)值模型、狀態(tài)方程以及強(qiáng)度模型等，進(jìn)行動(dòng)態(tài)模擬研究。圖2所示為花崗巖在800 m/s的射流沖擊作用下的損傷、破壞以及失效過程。在沖擊時(shí)間=2.2 μs射流沖擊巖石作用的最初階段，射流撞擊巖石表面形成動(dòng)壓力并以應(yīng)力波的形式向巖石內(nèi)部傳播，該過程中伴隨著巖石的損傷和失效。從圖2可見：射流沖撞巖石初期存在2種破壞區(qū)域，即由剪應(yīng)力導(dǎo)致的剪切破壞區(qū)以及最大拉應(yīng)力引起的破壞區(qū)。在射流連續(xù)沖擊作用下，射流作用點(diǎn)區(qū)域受到持續(xù)的壓應(yīng)力作用，使巖石局部彈性勢(shì)能不斷增加，進(jìn)而使剪切失效區(qū)域周圍萌發(fā)拉伸失效的裂紋。在=6.0~12.0 μs階段，萌生的裂紋在拉伸應(yīng)力作用下在縱向不斷擴(kuò)展形成放射性裂紋，但花崗巖內(nèi)部的層狀裂紋擴(kuò)展緩慢，幾乎不影響花崗巖的宏觀損傷破裂失效，該失效情況與文獻(xiàn)[12]中實(shí)驗(yàn)結(jié)果一致。圖2(d)所示為沖擊時(shí)間12.0 μs時(shí)巖石的損傷場。從圖2(d)可見：巖石的損傷值大致上隨距離射流入射點(diǎn)遠(yuǎn)近而變化，且距離入射點(diǎn)越遠(yuǎn)，巖石損傷度越小，這主要是應(yīng)力波損傷能力隨傳播距離逐漸降低。圖3所示為不同位置單元的損傷失效過程。由圖3可見：巖石單元的損傷失效是在瞬間內(nèi)完成的，且呈階躍性。兩單元失效時(shí)間間隔越來越大，說明射流沖擊能量在巖石內(nèi)部是急劇衰減的，且射流損傷破巖過程是微秒量級(jí)。

沖擊時(shí)間/μs: (a) 2.2; (b) 6.0; (c) 12.0; (d) 12.0

圖2 巖石損傷破壞過程

Fig. 2 Processes of rock damage and failure

圖3 巖石單元損傷過程

3 高強(qiáng)度巖石損傷過程的影響因素分析

3.1 入射角對(duì)巖石損傷場的影響

入射角是指射流軸線與巖石表面所成的夾角(≤90°)。圖4所示為射流沖擊速度為800 m/s，入射角分別為80°，70°，60°，50°，40°以及30°時(shí)，在=10 μs的巖石損傷場，其中入射角為90°時(shí)巖石損傷場如圖5(d)所示。從模擬結(jié)果可見：射流以不同的入射角沖擊損傷破巖時(shí)，其產(chǎn)生的巖石損傷場有明顯的差異，故射流入射角是影響水射流沖擊損傷破巖的1個(gè)重要參數(shù)。射流的入射角越小，由剪切作用主導(dǎo)產(chǎn)生的壓碎區(qū)越小，縱向裂紋和層向裂紋的擴(kuò)展長度也相應(yīng)減小。當(dāng)射流入射角較小時(shí)(如30°和40°)，沒有放射性裂紋產(chǎn)生，層狀裂紋也未能得到有效的擴(kuò)展，其主要是因?yàn)椋寒?dāng)射流入射角較小時(shí)，受到巖石表面液動(dòng)摩擦因數(shù)的影響，射流沖擊能量未能有效地正向撞擊損傷巖石，導(dǎo)致此時(shí)對(duì)巖石的損傷破壞效果較差。當(dāng)射流的入射角較大時(shí)(如70°和80°等)，射流沖擊損傷破巖效果較好，其主要是因?yàn)椋荷淞鹘朴诖怪睕_擊巖石時(shí)沖擊能量得到有效的利用，應(yīng)力波損傷破巖能力增強(qiáng)，巖石損傷區(qū)域增大。從以上分析可見：在射流沖擊損傷破巖過程中，射流的入射角存在1個(gè)臨界值，大于70°時(shí)沖擊損傷破巖效果較好。模擬結(jié)果與文獻(xiàn)[4]中實(shí)驗(yàn)結(jié)果變化規(guī)律類似，其中文獻(xiàn)[4]存在1個(gè)最優(yōu)值主要是實(shí)驗(yàn)中偏斜一定角度的射流容易將損傷破壞的巖石碎渣帶走引起的，而射流垂直入射則阻礙此過程。

入射角/(°): (a) 80; (b) 70; (c) 60; (d) 50; (e) 40; (f) 30

圖4 入射角對(duì)巖石損傷場的影響

Fig. 4 Rock damage field with different incidence angles

3.2 沖擊速度對(duì)巖石損傷場的影響

射流沖擊損傷破巖實(shí)際上是射流能量在巖石中傳播和衰減過程。射流速度決定射流的沖擊能量，其對(duì)巖石的損傷破壞場有直接的影響。圖5所示為花崗巖分別在200，400，600以及800 m/s沖擊速度下，=10 μs時(shí)的巖石損傷場。從數(shù)值模擬結(jié)果可見：射流沖擊速度越大，其損傷破巖能力越強(qiáng)，損傷范圍也越大。當(dāng)射流速度為200 m/s時(shí)，射流對(duì)巖石的沖擊損傷主要集中于射流入射點(diǎn)下方，且對(duì)巖石的損傷深度和范圍都比較有限，在表面形成很小的破碎坑。當(dāng)射流速度為400 m/s時(shí)，巖石損傷失效區(qū)域變大，且在失效區(qū)周圍萌生、擴(kuò)展了放射性裂紋，但此時(shí)巖石內(nèi)無層狀裂紋產(chǎn)生。當(dāng)射流速度為600 m/s時(shí)，射流損傷能力進(jìn)一步增強(qiáng)，巖石損傷范圍變大，且在巖石內(nèi)部萌生、擴(kuò)展了層狀裂紋。當(dāng)射流速度為800 m/s時(shí)，放射性裂紋擴(kuò)展更為顯著，此時(shí)射流能量可以損傷至底部的巖石單元。從數(shù)值模擬結(jié)果可見：射流速度過小難以對(duì)巖石進(jìn)行高效損傷破壞，主要以沖蝕作用破壞巖石的表面；射流速度達(dá)到一定值時(shí)才能夠使巖石出現(xiàn)層狀裂紋，產(chǎn)生體積破壞，如本文花崗巖的有效沖擊損傷破壞速度應(yīng)不小于600 m/s。

沖擊速度/(m?s?1): (a) 200; (b) 400; (c) 600; (d) 800

圖5 巖石在不同沖擊速度下的損傷場

Fig. 5 Rock damage field with different impact velocities

3.3 平移速度對(duì)對(duì)巖石損傷場的影響

在實(shí)際環(huán)境中，射流損傷破巖時(shí)通常具有一定的平移速度。圖6所示為在射流沖擊速度為600 m/s，平移速度分別為10，20，30和40m/s，=10 μs時(shí)的巖石損傷場。從圖6可見：射流平移速度對(duì)巖石的縱向損傷深度影響不大，也對(duì)放射性裂紋的萌生、擴(kuò)展影響很小，這主要由于射流平移速度相對(duì)于沖擊速度而言很小，而射流損傷破巖過程在微秒量級(jí)。平移速度對(duì)層狀裂紋的擴(kuò)展形式有一定的影響，平移速度在對(duì)平移方向上的層狀裂紋擴(kuò)展有積極的作用，這主要是射流沖擊巖石表面在平移方向上產(chǎn)生較大的張拉力引起的?？偟膩碚f，當(dāng)射流平移速度比沖擊速度小很多時(shí)，其對(duì)射流沖擊巖石損傷場影響不大。

平移速度/(m?s?1): (a) 10; (b) 20; (c) 30; (d) 40

圖6 巖石在不同平移速度下的損傷場

Fig. 6 Rock damage field at different translation velocities

4 結(jié)論

1) 引入Johnson-Holmquist modelⅡ(JH-2)表征巖石力學(xué)特性，利用SPH/FEA耦合算法建立了高壓水射流沖擊破巖模型，很好地重現(xiàn)了水射流損傷破巖過程中巖石失效、裂紋擴(kuò)展以及單元損傷破壞過程，對(duì)水射流沖擊破巖機(jī)理的研究有一定的參考意義。

2) 高壓水射流沖擊損傷破壞巖石過程中存在2種失效破壞區(qū)域，即剪切破壞區(qū)和拉伸破壞區(qū)：剪切破壞區(qū)主要集中在射流入射點(diǎn)下方，形成巖石的沖擊粉碎區(qū)；拉伸破壞主要出現(xiàn)在遠(yuǎn)離射流入射點(diǎn)區(qū)域，其主導(dǎo)放射性裂紋和層狀裂紋的萌生、擴(kuò)展。此外，巖石單元的損傷破壞呈階躍性，損傷過程為微秒量級(jí)，且距射流入射點(diǎn)越近的巖石損傷程度越大。

3) 射流的入射角和沖擊速度對(duì)巖石損傷場影響較大，而射流平移速度對(duì)巖石損傷場影響不大：射流沖擊速度較小時(shí)以沖蝕破壞巖石為主，射流達(dá)到一定值才能夠有效損傷破壞巖石，產(chǎn)生放射性裂紋和層狀裂紋，本文射流的有效沖擊速度為600 m/s；射流沖擊損傷破壞巖石過程中，射流入射角存在一個(gè)最佳區(qū)域，大于70°時(shí)損傷破壞巖石效果較好；由于射流平移速度較沖擊速度小得多，其對(duì)射流沖擊巖石損傷場影響不明顯。

[1] 廖華林, 李根生, 熊偉. 超高壓水射流輔助破巖鉆孔研究進(jìn)展[J]. 巖石力學(xué)與工程學(xué)報(bào), 2002, 21(2): 2583?2587.
LIAO Hualin, LI Gengsheng, XIONG Wei.[J]. Chinese Journal of Rock Mechanics and Engineering, 2002, 21(2): 2583?2587.

[2] 李根生, 廖華林. 超高壓水射流沖蝕切割巖石斷口微觀斷裂機(jī)理實(shí)驗(yàn)研究[J]. 高壓物理學(xué)報(bào), 2005, 19(4): 337?343.
LI Gensheng, LIAO Hualin. Micro-Failure mechanism analysis and test study for rock failure surface under water jet impact[J]. Chinese Journal of High Pressure Physics, 2005, 19(4): 337?343.

[3] Guo C W, Dong L. Energy consumption in communication of mica with cavitation abrasive water jet[J]. Journal of China University of Mining & Technology, 2007, 17(2): 251?254.

[4] 馮云春, 徐依吉, 趙付國. 超高壓淹沒射流破巖規(guī)律實(shí)驗(yàn)研究[J]. 高壓物理學(xué)報(bào), 2005, 19(1): 66?71.
FENG Yunchun, XU Yiji, ZHAO Fuguo. Experimental study on the breaking of rock by ultra-high pressure water jetting[J]. Chinese Journal of High Pressure Physics, 2005, 19(1): 66?71.

[5] 李根生, 廖華林, 黃中偉, 等. 超高壓水射流作用下巖石損傷破碎機(jī)理[J]. 機(jī)械工程學(xué)報(bào), 2009, 45(10): 285?293.
LI Gensheng, LIAO Hualin, HUANG Zhongwei, et al. Rock damage mechanisms under ultra high pressure water jet impact[J]. Journal of Mechanical Engineering, 2009, 45(10): 285?293.

[6] Momber A W. Deformation and fracture of rock due to high-speed liquid impingement[J]. International Journal of Fracture, 2004, 130: 683?704.

[7] 劉佳亮, 司鵠. 高壓水射流破碎高圍壓巖石損傷場的數(shù)值模擬[J]. 重慶大學(xué)學(xué)報(bào), 2011, 34(4): 40?46.
LIU Jialiang, SI Hao. Numerical simulation on damage field of high pressure water jet breaking rock under high ambient pressure[J]. Journal of Chongqing University, 2011, 34(4): 40?46.

[8] 廖華林, 李根生. 超高壓水射流沖擊巖石的流固耦合分析[J]. 水動(dòng)力學(xué)研究與進(jìn)展, 2004, 19(4): 453?457.
LIAO Hualin, LI Gensheng. Fluid-structure interaction of high pressure jets impinging on rock[J]. Journal of Hydrodynamics, 2004, 19(4): 453?457.

[9] Ni H J, Wang R H, Zhang Y Q. Numerical simulation study on rock breaking mechanism and process under high pressure water jet[J]. Applied Mathematics and Mechanics (English Edition), 2005, 26(12): 1445?1452.

[10] 孫清德, 汪志明, 于軍泉, 等. 高壓水射流破巖規(guī)律的數(shù)值模擬研究[J]. 巖土力學(xué), 2005, 26(6): 978?982.
SUN Qingde, WANG Zhiming, YU Junquan, et al. A disquisition on breaking mechanism of high pressure jet impacting on rock[J]. Rock and Soil Mechanics, 2005, 26(6): 978?982.

[11] 司鵠, 王丹丹, 李曉紅. 高壓水射流破巖應(yīng)力波效應(yīng)的數(shù)值模擬[J]. 重慶大學(xué)學(xué)報(bào), 2008, 31(8): 942?945.
SI Hao, WANG Dandan, Li Xiaohong. Stress wave effect in numerical simulation on rock breaking under high pressure water jet[J]. Journal of Chongqing University, 2008, 31(8): 942?945.

[12] 盧義玉, 張賽, 劉勇, 等. 脈沖水射流破巖過程中的應(yīng)力波效應(yīng)分析[J]. 重慶大學(xué)學(xué)報(bào), 2012, 35(1): 117?124.
LU Yiyu Y, ZHANG Sai, LIU Yong, et al. Analysis on stress wave effect during the process of rock breaking by pulsed jet[J]. Journal of Chongqing University, 2012, 35(1): 117?124.

[13] Wang J M, Gao N, Gong W J. Abrasive water-jet machining simulation by coupling smoothed particle hydrodynamics/finite element method[J]. Chinese Journal of Mechanical Engineering, 2010, 23(1): 1?6.

[14] 馬利. 無網(wǎng)格法及液體射流高速碰撞與侵徹模擬[D]. 杭州: 浙江大學(xué)航空航天學(xué)院, 2007: 61?62, 76.
MA Li. Mesh-free method and its application on high velocity fluid jet impact and penetration simulation[D]. Hangzhou: Zhejiang University. School of Aeronautics and Astronautics, 2007: 61?62, 76.

[15] Ma L, Bao R, Guo Y. Waterjet penetration simulation by hybrid code of SPH and FEA[J]. International Journal of Impact Engineering, 2008, 35(9): 1035?1042.

[16] Limido J, Espinosa C, Salaun M, et al.[J]. International Journal of Mechanical Sciences, 2007, 49: 898?908.

[17] 倪紅堅(jiān), 王瑞和, 葛洪魁. 高壓水射流破巖的數(shù)值模擬分析[J]. 巖石力學(xué)與工程學(xué)報(bào), 2004, 23(4): 550? 554.
NI Hongjian, WANG Ruihe, GE Hongkui. Numerical simulation on rock breaking under high pressure waterjet[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(4): 550?554.

[18] Century Dynamics Inc. Autodyn theory manual[Z]. Concord, CA: Century Dynamics Inc, 2003: 177?178.

[19] Holmquist T J, Templeton D W, Bishnoi K D. Constitutive modeling of aluminum nitride for large strain, high strain rate, and high-pressure application[J]. International Journal of Impact Engineering, 2001, 25(3): 211?231.

[20] 朱哲明, 王超, 王蒙, 等. 爆炸載荷作用下缺陷巖體破壞特征的數(shù)值模擬研究[J]. 四川大學(xué)學(xué)報(bào)(工程科學(xué)版), 2011, 43(2): 1?8.
ZHU Zheming, WANG Chao, WANG Meng, et al. Numerical predication on the fracture characteristics of rock with defects under blasting loads[J]. Journal of Sichuan University (Engineering Science Edition), 2011, 43(2): 1?8.

[21] 朱哲明, 李元鑫, 周志榮, 等. 爆炸載荷下缺陷巖體的動(dòng)態(tài)響應(yīng)[J]. 巖石力學(xué)與工程學(xué)報(bào), 2011, 30(6): 1157?1167.
ZHU Zheming, LI Yuanxin, ZHOU Zhirong, et al. Dynamic response of defected rock under blasting load[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(6): 1157?1167.

[22] Dehghan Banadaki M M, Mohanty B. Numerical simulation of stress wave induced fractures in rock[J]. International Journal of Impact Engineering, 2012, 40(2): 16?25.

Numerical analysis on damage field of rock fragmentation with water jet

JIANG Hongxiang, DU Changlong, LIU Songyong, GAO Kuidong

(School of Mechanical and Electrical Engineering, China University of Mining & Technology, Xuzhou 221116, China)

A constitutive model, Johnson-Holmquist modelⅡ(JH-2) for brittle material was adopted to characterize the mechanical properties of rock. And the numerical model of rock breaking with high-pressure water jet impacting was established with coupling Smoothed Particle Hydrodynamics (SPH) method and finite element analysis (FEA) method. It well simulated the rock failure, cracks propagation and the variations of rock element damage in different positions. Moreover, the effects of impact velocity, incidence angle and translation velocity of water jet on rock damage field were investigated. The results show that rock damage, failure and cracks propagation are caused by shear and tensile action combinations. The damage and failure zone of rock near the jet impacting point is mainly due to shear action; however, the radial and spall cracks propagation are dominated by tensile action. The damage of rock element represents a step changing with time, and the response of rock damage and failure process under water jet impacting is within microseconds. The incidence angle of water jet impacting rock has an effective range, and better damage effect can be obtained when the incidence angle is more than 70o. The surface erosion of rock is primary with low impact velocity of water jet, actual damage and failure of rock as radial and spall cracks occurrence can be achieved only when impact velocity of water jet reaches a certain value. The effect of water jet translation velocity on rock damage field is not remarkable because the translation velocity of water jet is much smaller than its impact velocity.

water jet; SPH; FEA; rock damage and failure; JH-2 constitutive model

TD825

A

1672?7207(2015)01?0287?08

2014?02?13；

2014?04?20

國家高科技研究發(fā)展計(jì)劃(863計(jì)劃)項(xiàng)目(2012AA062104)；國家自然科學(xué)基金面上資助項(xiàng)目(51375478)；江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程資助項(xiàng)目(SZBF2011-6-B35) (Project(2012AA062104) supported by the National High Technology Research and Development Program (863 Program); Project(51375478) supported by the National Natural Science Foundation of China; Project(SZBF2011-6-B35) supported by the Priority Academic Program Development of Jiangsu High Education Institute of China)

杜長龍，教授，博士生導(dǎo)師，從事高壓水射流破巖技術(shù)研究；E-mail: jdjxx3@cumt.edu.cn

10.11817/j.issn.1672?7207.2015.01.039

(編輯 楊幼平)