非均質(zhì)巖石單軸壓縮試驗(yàn)破壞過(guò)程細(xì)觀模擬及分形特性

李守巨,李 德,武 力,曹麗娟

(1.大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧大連 116024;2.大連交通大學(xué)機(jī)械工程學(xué)院,遼寧大連 116028;3.大連海洋大學(xué)機(jī)械與動(dòng)力工程學(xué)院,遼寧大連 116023)

非均質(zhì)巖石單軸壓縮試驗(yàn)破壞過(guò)程細(xì)觀模擬及分形特性

李守巨1,李 德1,武 力2,曹麗娟3

(1.大連理工大學(xué)工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧大連 116024;2.大連交通大學(xué)機(jī)械工程學(xué)院,遼寧大連 116028;3.大連海洋大學(xué)機(jī)械與動(dòng)力工程學(xué)院,遼寧大連 116023)

為了細(xì)觀模擬非均質(zhì)巖石試樣單軸壓縮試驗(yàn)的破壞過(guò)程,提出了一種新的有限元數(shù)值模擬方法。該方法以莫爾-庫(kù)侖破壞準(zhǔn)則為基礎(chǔ),將巖石視為非均質(zhì)和準(zhǔn)脆性材料,并假定巖石單元的抗剪強(qiáng)度參數(shù)服從正態(tài)分布。有限元細(xì)觀模擬結(jié)果表明,巖樣非均質(zhì)性對(duì)本構(gòu)關(guān)系曲線形狀有顯著影響,模擬的應(yīng)力-應(yīng)變曲線與試驗(yàn)曲線基本吻合。基于分形幾何理論,采用計(jì)盒維數(shù)方法計(jì)算單軸壓縮試驗(yàn)不同荷載條件下的巖石破壞分形維數(shù),建立了分形維數(shù)與頂部位移之間的非線性演化模型。分形分析研究表明,巖石試樣的破壞過(guò)程具有分形特性,分形維數(shù)可用來(lái)定量地表征單軸壓縮試驗(yàn)巖樣破壞過(guò)程。

非均質(zhì)巖石;單軸壓縮試驗(yàn);分形特性;有限元;莫爾-庫(kù)侖準(zhǔn)則;計(jì)盒維數(shù)

Key words:inhomogeneous rock;uniaxial compression test;fractal characteristics;finite element method;Mohr-Coulomb fracture criterion;box-counting dimension

巖石是由多種礦物組成的非均質(zhì)體,含有大量的節(jié)理和裂隙,從微觀和細(xì)觀的角度來(lái)衡量,巖石具有非均質(zhì)性、各向異性和不連續(xù)性等特性。研究表明,在單軸壓縮情況下,巖石試樣的破壞過(guò)程是隨分布強(qiáng)度參數(shù)的巖石顆粒的剛度逐漸劣化的結(jié)果。從連續(xù)介質(zhì)力學(xué)觀點(diǎn)分析,巖石試樣的破壞方式可以歸結(jié)為壓縮破壞、拉伸破壞和剪切破壞等方式。關(guān)于巖石試樣和巖體結(jié)構(gòu)的破壞過(guò)程與機(jī)理已經(jīng)進(jìn)行了廣泛研究,取得了許多研究成果。Zhao[1]采用試驗(yàn)方法研究了巖石試樣裂紋的萌生和擴(kuò)展過(guò)程,確定了裂紋分形維數(shù)與應(yīng)力之間的關(guān)系,建立了巖石分形損傷本構(gòu)模型。Innaurato等[2]采用實(shí)驗(yàn)和數(shù)值模擬相結(jié)合的方法研究了TBM滾刀的破巖過(guò)程,建立了滾刀貫入深度與貫入力之間的關(guān)系,研究了滾刀作用下巖石的剪切破壞問(wèn)題。Li等[3]建立了表征巖石宏觀變形與細(xì)觀斷裂之間關(guān)系的本構(gòu)模型,并使用該模型研究了巖石微裂紋擴(kuò)展、非線性變形、加載-卸載過(guò)程的滯變現(xiàn)象以及彈性模量和泊松比的變化等問(wèn)題。夏毅敏等[4]利用虛擬樣機(jī)技術(shù),采用Pro/E軟件三維參數(shù)化建模,ADAMS軟件可視化仿真,建立了一個(gè)鈷結(jié)殼螺旋采集式采礦頭的虛擬樣機(jī)。結(jié)合虛擬樣機(jī)和相似理論建立了一個(gè)具有雙螺旋線、可安裝36個(gè)截齒的采礦頭試驗(yàn)臺(tái),并通過(guò)大量試驗(yàn)驗(yàn)證了所建立虛擬樣機(jī)的正確性。Jiang等[5]指出了含有節(jié)理巖體的地下結(jié)構(gòu)開挖的變形和失穩(wěn)機(jī)理是由這些結(jié)構(gòu)面的幾何分布特性所控制的,并采用計(jì)盒維數(shù)方法定量化評(píng)估了結(jié)構(gòu)面分布的分形特性。Harthong等[6]采用離散元和離散裂隙網(wǎng)絡(luò)耦合模型數(shù)值模擬了巖石試樣三維壓縮試驗(yàn)的破壞過(guò)程,研究了巖體的強(qiáng)度特性并指出裂隙巖體的變形機(jī)理主要取決于裂隙強(qiáng)度。朱萬(wàn)成等[7]采用RFPA軟件、并假定巖石的彈性模量和強(qiáng)度指標(biāo)服從Weibull分布,數(shù)值模擬了地下結(jié)構(gòu)圍巖裂隙的萌生、擴(kuò)展和貫通過(guò)程。Liu等[8]數(shù)值模擬了單壓頭和雙壓頭兩種工況下的破巖過(guò)程,提出了通過(guò)優(yōu)化壓頭間距形成較大巖塊的可能性,從而控制裂紋的擴(kuò)展方向,實(shí)現(xiàn)降低巖石破碎比功的目的。黃達(dá)等[9]基于不同卸荷速率和初始圍壓條件下三軸高應(yīng)力大理巖卸圍壓試驗(yàn),結(jié)合分形理論和能量原理,研究了高應(yīng)力卸荷條件下巖石破裂塊度分布規(guī)律及其與能量耗散和釋放的相關(guān)性。分形理論已經(jīng)應(yīng)用于采動(dòng)裂隙網(wǎng)絡(luò)實(shí)時(shí)演化特性表征、砂巖細(xì)觀力學(xué)效應(yīng)分析、三維空間巖石裂紋分形表征和煤巖沖擊破碎塊度特征分析等巖石力學(xué)領(lǐng)域[10-14]。

盡管如此,還有許多問(wèn)題沒(méi)有得到很好解決。本文的目的在于提出一種基于有限元理論的巖石破壞過(guò)程細(xì)觀數(shù)值模擬方法,數(shù)值模擬非均質(zhì)巖石試樣單軸壓縮試驗(yàn)破壞過(guò)程,研究巖石破壞過(guò)程的分形表征方法,建立荷載與分形維數(shù)之間的非線性映射關(guān)系模型。

1 非均質(zhì)巖石試樣單軸壓縮試驗(yàn)的有限元細(xì)觀模擬

巖石試樣單軸壓縮試驗(yàn)如圖1所示,H為巖樣高度,W為巖樣寬度,P為頂部垂直荷載。

圖1 巖石試樣單軸壓縮試驗(yàn)Fig.1 Rock specimen under uniaxial compression test

將巖石試樣離散為有限元單元模型。在加載過(guò)程中,每個(gè)單元的應(yīng)力狀態(tài)如圖2所示。

圖2 單元可能破壞面的法向和切向應(yīng)力Fig.2 Normal and tangent stresses in some plane for any element

根據(jù)力的平衡原理,單元法向應(yīng)力和切向應(yīng)力與主應(yīng)力之間的關(guān)系為

式中,σN為單元可能破壞面的法向應(yīng)力;σT為單元可能破壞面的切向應(yīng)力;σ1為有限元計(jì)算得到的單元第3主應(yīng)力(取負(fù)號(hào)),即巖石力學(xué)中的大主應(yīng)力; σ3為有限元計(jì)算得到的單元第1主應(yīng)力(取負(fù)號(hào)),即巖石力學(xué)中的小主應(yīng)力;α為每個(gè)單元最可能破壞面外法線與垂直軸的夾角,即破壞面與水平軸的夾角,其大小與單元的內(nèi)摩擦角有關(guān)。

為表征單軸壓縮試驗(yàn)巖樣的剪切破壞特性,采用Mohr-Coulomb準(zhǔn)則來(lái)評(píng)價(jià)巖石的破壞與否,巖石單元的Mohr-Coulomb破壞準(zhǔn)則為

式中,c為巖石材角。

非均質(zhì)巖石試樣單軸壓縮試驗(yàn)?zāi)M分析是在ANSYS有限元軟件平臺(tái)上完成的,編寫了APDL參數(shù)化語(yǔ)言程序。巖石為中晶大理巖,粒徑為 1～4 mm,平均粒徑為1.5 mm,實(shí)際巖石試樣的尺寸為高100 mm、直徑50 mm[15]?？紤]到軸對(duì)稱性,有限元模型簡(jiǎn)化為二維模型,其高100 mm、寬50 mm,如圖1所示。有限元模型網(wǎng)格尺寸為1 mm,與巖石顆粒的粒徑相近,共劃分為5 000個(gè)單元。

有限元細(xì)觀模擬巖石試樣的基本力學(xué)特性參數(shù)引自文獻(xiàn)[15],詳見(jiàn)表1,表中抗剪強(qiáng)度指標(biāo)參數(shù)為平均值。為表征巖石材料的非均質(zhì)性、各向異性、不連續(xù)性等特性,假定巖石單元抗剪強(qiáng)度指標(biāo)c,φ滿足正態(tài)分布[16],并利用其平均值與標(biāo)準(zhǔn)差隨機(jī)生成5 000個(gè)不同c和φ值,并賦給5 000個(gè)單元,進(jìn)而實(shí)現(xiàn)對(duì)每個(gè)單元賦予不同的抗剪強(qiáng)度指標(biāo)以模擬真實(shí)巖石材料的非均質(zhì)性等特性。為模擬巖石試樣的加載過(guò)程,在有限元模型的頂部施加垂直位移荷載,位移加載步長(zhǎng)為0.025 mm,共計(jì)20個(gè)加載步,累計(jì)垂直位移為0.50 mm。

表1 巖石試樣力學(xué)性能參數(shù)[15]Table 1 Basic mechanical properties for rock specimen[15]

根據(jù)Mohr-Coulomb剪切破壞準(zhǔn)則,在每個(gè)荷載步分別判斷有限元模型中每個(gè)單元的當(dāng)前應(yīng)力狀態(tài),當(dāng)滿足Mohr-Coulomb剪切破壞準(zhǔn)則即式(4)時(shí),認(rèn)為該單元破壞,進(jìn)而修改破壞單元的材料屬性,對(duì)破壞單元的彈性模量進(jìn)行剛度弱化處理。根據(jù)損傷力學(xué)理論,破壞后單元的彈性模量為

其中,EF為破壞后單元的彈性模量;R為破壞比[7];E為單元未破壞時(shí)的彈性模量;Df為損傷指標(biāo),Df=0,表示單元處于無(wú)損傷狀態(tài),Df=1,表示單元完全喪失抵抗變形的能力,Df的大小與巖石的脆性性質(zhì)、延性性質(zhì)以及荷載大小相關(guān)。在三軸壓縮條件下,損傷指標(biāo)還與圍壓相關(guān)。在實(shí)際問(wèn)題中,由于相鄰單元的約束作用,巖樣局部破壞后,單元的彈性模量不可能立即降為0,也就是說(shuō)滿足Mohr-Coulomb剪切破壞準(zhǔn)則的單元仍具有部分抵抗變形的能力,其損傷指標(biāo)不可能馬上變?yōu)?,根據(jù)相關(guān)文獻(xiàn),本文取破壞比R= 1[7]。

在巖石內(nèi)聚力標(biāo)準(zhǔn)差為5.0 MPa,內(nèi)摩擦角標(biāo)準(zhǔn)差為5.0°的情形下有限元模擬非均質(zhì)巖石試樣破壞的過(guò)程如圖3所示,其中,顏色較深部分的單元為破壞單元。有限元細(xì)觀模擬結(jié)果表明,巖石試樣單軸壓縮破壞形成了一系列剪切破壞帶,剪切帶外法線與垂直軸的夾角約為45°,而實(shí)驗(yàn)值為60°左右[15],數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果基本一致。

圖3 不同荷載條件下巖石試樣破壞區(qū)域Fig.3 Distributions of fracturing region of rock specimen in different load case

不同剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差下數(shù)值模擬的頂部荷載與頂部位移關(guān)系曲線與試驗(yàn)曲線對(duì)比如圖4(a)所示,模擬的應(yīng)力-應(yīng)變曲線與試驗(yàn)曲線對(duì)比如圖4(b)所示。其中,模擬曲線1代表的巖樣剪切強(qiáng)度參數(shù)內(nèi)聚力標(biāo)準(zhǔn)差 Sc=5.00 MPa,內(nèi)摩擦角標(biāo)準(zhǔn)差 Sφ= 5.00°;模擬曲線2代表的巖樣剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差Sc=2.50 MPa,Sφ=2.50°;模擬曲線3代表的巖樣剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差Sc=1.25 MPa,Sφ=1.25°。由模擬曲線3計(jì)算巖石抗壓強(qiáng)度為122 MPa,而實(shí)驗(yàn)值為114 MPa,相對(duì)誤差為7%,兩者比較接近。

圖4 模擬曲線與試驗(yàn)曲線對(duì)比Fig.4 Comparisons between simulated curves and test curves

由圖4(b)可以看出,巖石試樣的非均質(zhì)性對(duì)模擬曲線的結(jié)果有明顯的影響,隨著剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差數(shù)值的降低,模擬曲線的峰值強(qiáng)度逐漸增大,其峰后曲線越來(lái)越陡,應(yīng)力跌落程度加劇。巖樣破壞具有明顯的脆性破壞特征,這與有限元模擬將巖樣視為準(zhǔn)脆性材料假設(shè)相符,試驗(yàn)曲線則具有一定的延性破壞特征。模擬曲線與試驗(yàn)曲線在彈性階段十分吻合,在峰后殘余強(qiáng)度階段體現(xiàn)出不同的破壞特性。峰后殘余強(qiáng)度階段,當(dāng)巖樣接近破壞時(shí),模擬曲線會(huì)有一個(gè)緩慢上升的趨勢(shì),這與實(shí)驗(yàn)值差異較大。其原因在于當(dāng)巖樣接近破壞時(shí),其破壞比R可能會(huì)急劇下降,R隨著軸向應(yīng)變的增加而減小,但在數(shù)值模擬模型中始終假定為一常數(shù)。

為了進(jìn)一步驗(yàn)證所提出模擬方法的有效性,數(shù)值模擬了一個(gè)算例,其中大理巖樣應(yīng)力-應(yīng)變?cè)囼?yàn)曲線數(shù)據(jù)引自文獻(xiàn)[17],相關(guān)計(jì)算參數(shù)見(jiàn)表2,取破壞比R=0.1。

表2 第2個(gè)巖石試樣力學(xué)性能參數(shù)Table 2 Mechanical properties of the second rock specimen

有限元細(xì)觀模擬曲線與試驗(yàn)曲線對(duì)比如圖4(c)所示。其中,巖樣剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差為 Sc= 1.25 MPa,Sφ=1.25°。由圖4(c)可看出,模擬曲線與試驗(yàn)曲線在彈性階段十分吻合,在峰后殘余強(qiáng)度階段有所差別。由模擬曲線計(jì)算巖石抗壓強(qiáng)度為134 MPa,而文獻(xiàn)[17]所給巖石抗壓強(qiáng)度對(duì)應(yīng)為130 MPa,相對(duì)誤差為3%,兩者十分相近。有限元模擬結(jié)果與算例結(jié)果吻合地較好,進(jìn)一步驗(yàn)證了數(shù)值模擬方法的有效性。

2 非均質(zhì)巖石試樣破壞過(guò)程的分形特性

式中,Ds為破壞區(qū)域的自相似分形維數(shù);r為正方形盒子邊長(zhǎng);N(r)為用邊長(zhǎng)為r的正方形盒子覆蓋整個(gè)圖形中破壞區(qū)域所需盒子的數(shù)目。

非均質(zhì)巖石試樣單軸壓縮試驗(yàn)破壞分形維數(shù)的計(jì)算是在Matlab數(shù)值計(jì)算軟件平臺(tái)上完成的,編寫了Matlab語(yǔ)言程序。算法主要流程如下:①前處理。讀入原始圖像,圖像灰度化預(yù)處理,圖像二值化預(yù)處理,存儲(chǔ)相關(guān)數(shù)據(jù);②求解。以邊長(zhǎng)r大小的正方形盒子覆蓋二值圖,統(tǒng)計(jì)巖樣破壞區(qū)域正方形盒子數(shù)N(r),保存相關(guān)數(shù)據(jù),其中,取r=0.14,0.28,0.56, 1.11,2.22,4.44,8.89,17.78 mm;③ 后處理。對(duì)r和N(r)取雙對(duì)數(shù)坐標(biāo)進(jìn)行回歸分析,得到回歸方程

式中,a,b為常數(shù),b值即為巖石破壞分形維Ds,其與方程log N(r)=log a-blog r是等價(jià)的,只需對(duì)式(7)兩邊求對(duì)數(shù)即可。

由圖3可以看出,單軸壓縮試驗(yàn)巖石試樣破壞區(qū)域的發(fā)展過(guò)程具有統(tǒng)計(jì)自相似性,近似符合分形特性,可用分形維數(shù)表征。頂部位移為0.30 mm時(shí)的巖石破壞區(qū)域的分形特性擬合曲線如圖5所示,相關(guān)系數(shù)R2=0.997 2,說(shuō)明分形維數(shù)具有較高的可信度,分形維數(shù)Ds約為1.867 1。

分形理論可有效描述自然界中許多外表不規(guī)律事物的內(nèi)在規(guī)律性,在生物學(xué)、物理、化學(xué)、天文學(xué)、經(jīng)濟(jì)學(xué)、巖土工程、計(jì)算機(jī)圖形學(xué)等諸多領(lǐng)域展現(xiàn)出其獨(dú)特的應(yīng)用前景[18-19]。

分形維數(shù)是分形理論的核心內(nèi)容之一,其使用分?jǐn)?shù)維代替?zhèn)鹘y(tǒng)歐式空間的整數(shù)維,能夠定量地表征具有非整數(shù)值維數(shù)的事物,如Cantor集、Koch曲線等經(jīng)典分形。分形維數(shù)具有多種定義方式,本文采用自相似分形維數(shù),也稱為盒維數(shù)或計(jì)盒維數(shù)(box-counting dimension),計(jì)算公式如下:

圖5 頂部位移為0.30 mm時(shí)的破壞區(qū)域的分形特性Fig.5 Fractal characteristics of fracturing region with a top displacement of 0.30 mm

不同頂部位移情形下巖樣破壞分形維數(shù)演化曲線如圖6所示。由圖6可看出,巖石試樣加載過(guò)程中破壞分形維演化曲線受巖樣非均質(zhì)性影響顯著。其中,分形曲線1代表的巖樣剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差Sc= 5.00 MPa,Sφ=5.00°;分形曲線2代表的巖樣剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差Sc=2.50 MPa,Sφ=2.50°;分形曲線3代表的巖樣剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差Sc=1.25 MPa,Sφ= 1.25°。隨著巖樣剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差的降低,巖石試樣的力學(xué)特性趨于更均勻,破壞單元開始出現(xiàn)時(shí)的頂部位移增大,在3組不同標(biāo)準(zhǔn)差條件下分別為0.125,0.175,0.200 mm。同時(shí),隨著巖樣剪切強(qiáng)度參數(shù)標(biāo)準(zhǔn)差的降低,破壞分形維演化曲線變得越來(lái)越陡。分形維數(shù)曲線最后趨于穩(wěn)定,隨著頂部位移的增加分形維數(shù)基本保持不變,表明巖樣已經(jīng)完全破壞。對(duì)應(yīng)3組不同的巖樣剪切強(qiáng)度標(biāo)準(zhǔn)差,巖樣完全破壞時(shí)的分形維數(shù)分別為1.840,1.850和1.924。說(shuō)明隨著巖樣非均質(zhì)性降低,巖樣完全破壞分形維數(shù)增大,分形維數(shù)有效地表征了單軸壓縮試驗(yàn)巖樣破壞的特性。

圖6 巖石試樣加載過(guò)程中破壞分形維數(shù)演化曲線Fig.6 Evolutional curves of fractal dimensions forfracturing process of rock specimen

巖石試樣加載過(guò)程中破壞分形維數(shù)隨荷載演化過(guò)程可采用Sigmoid型演化模型來(lái)表征[19],表達(dá)式為

式中,Dmax,Dmin,α,dt,dt0為分形維數(shù)演化模型中的參數(shù),采用回歸分析方法確定,見(jiàn)表3。

表3 分形維數(shù)Sigmoid演化模型中的參數(shù)Table 3 Model parameters of evolutional process of fractal dimension

巖樣破壞分形維數(shù)演化模型預(yù)測(cè)值與理論值對(duì)比如圖7所示。由圖7可看出,巖樣破壞分形維數(shù)演化模型預(yù)測(cè)值與理論值基本吻合,驗(yàn)證了該演化模型的有效性。

圖7 巖樣破壞分形維數(shù)演化模型預(yù)測(cè)值與理論值對(duì)比Fig.7 Comparison between predicted values and theoretical ones for fracturing process of rock specimen

3 結(jié) 論

(1)提出了非均質(zhì)巖石壓縮試驗(yàn)細(xì)觀模擬方法, 2個(gè)數(shù)值算例模擬的巖石極限抗壓強(qiáng)度與實(shí)驗(yàn)值的相對(duì)誤差分別為3%和7%,驗(yàn)證了數(shù)值模擬方法的有效性。

(2)研究結(jié)果表明,非勻質(zhì)巖樣單軸壓縮破壞具有統(tǒng)計(jì)自相似性,巖石試樣的非均質(zhì)性對(duì)破壞分形維數(shù)演化模型參數(shù)有較大影響,分形維數(shù)可有效地表征巖石試樣的破壞特性。

(3)建立了巖石破壞過(guò)程的分形維數(shù)Sigmoid演化模型,采用非線性回歸分析方法確定了模型中的參數(shù),該模型能夠計(jì)算分形維數(shù)隨著巖石試樣破壞過(guò)程的變化規(guī)律。

(4)進(jìn)一步的研究工作包括進(jìn)行三軸壓縮條件下巖石破壞過(guò)程分形特性分析和建立破壞分形維數(shù)演化模型參數(shù)與巖石力學(xué)特性參數(shù)之間的映射關(guān)系。

[1] Zhao Yonghong.Crack pattern evolution and a fractal damage constitutive model for rock[J].International Journal of Rock Mechanics and Mining Sciences,1998,35(3):349-366.

[2] Innaurato N,Oggeri C,Oreste P P,et al.Experimental and numerical studies on rock breaking with TBM tools under high stress confinement[J].Rock Mechanics and Rock Engineering,2007,40(5): 429-451.

[3] Li Chunlin,Richard Prikry,Erling Nordlunda.The stress-strain behavior of rock material related to fracture under compression[J].Engineering Geology,1998,49(3-4):293-302.

[4] 夏毅敏,卜英勇,袁富貴.深海鈷結(jié)殼螺旋采集式采礦頭仿真研究[J].煤炭學(xué)報(bào),2004,29(2):246-248.

Xia Yimin,Bu Yingyong,Yuan Fugui.Simulation study on the deep sea bed screw-collection mining head for cobalt crust[J].Journal of China Coal Society,2004,29(2):246-248.

[5] Jiang Yujing,Tanabashi Yosihiko,Li Bo,et al.Influence of geometrical distribution of rock joints on deformational behavior of underground opening[J].Tunneling and Underground Space Technology, 2006,21(5):485-491.

[6] Barthélémy Harthong,Luc Scholtès,Frédéric-Victor Donzé.Strength characterization of rock masses,using a coupled DEM-DFN model [J].Geophysical Journal International,2012,191(2):467-480.

[7] Zhu W C,Liu J,Tang C A,et al.Simulation of progressive fracturing processes around underground excavations under biaxial compression [J].Tunnelling and Underground Space Technology,2005,20(3): 231-247.

[8] Liu H Y,Kou S Q,Lindqvist P-A,et al.Numerical simulation of the rock fragmentation process induced by indenters[J].International Journal of Rock Mechanics and Mining Sciences,2002,39(4):491-505.

[9] 黃 達(dá),譚 清,黃潤(rùn)秋.高圍壓卸荷條件下大理巖破碎塊度分形特征及其與能量相關(guān)性研究[J].巖石力學(xué)與工程學(xué)報(bào), 2012,31(7):1379-1389.

Huang Da,Tan Qing,Huang Runqiu.Fractal characteristics of fragmentation and correlation with energy of marble under unloading with high confining pressure[J].Chinese Journal of Rock Mechanics and Engineering,2012,31(7):1379-1389.

[10] 高明忠,金文城,鄭長(zhǎng)江.采動(dòng)裂隙網(wǎng)絡(luò)實(shí)時(shí)演化及連通性特征[J].煤炭學(xué)報(bào),2012,37(9):1535-1540.

Gao Mingzhong,Jin Wencheng,Zheng Changjiang.Real-time evolution and connectivity of mined crack network[J].Journal of China Coal Society,2012,37(9):1535-1540.

[11] 姜立春,溫 勇,吳愛(ài)祥.AMD蝕化砂巖細(xì)觀力學(xué)效應(yīng)[J].煤炭學(xué)報(bào),2012,37(6):931-935.

Jiang Lichun,Wen Yong,Wu Aixiang.Mesoscopic mechanical effect on sandstone corroded by acid mine drainage[J].Journal of China Coal Society,2012,37(6):931-935.

[12] 周宏偉,張 濤,薛東杰.長(zhǎng)壁工作面覆巖采動(dòng)裂隙網(wǎng)絡(luò)演化特征[J].煤炭學(xué)報(bào),2011,36(12):1957-1962.

Zhou Hongwei,Zhang Tao,Xue Dongjie.Evolution of mining-induced crack network in overburden strata of longwall face[J].Journal of China Coal Society,2011,36(12):1957-1962.

[13] 薛東杰,周宏偉,易海洋.三維空間巖石裂紋分形曲線的分維估算[J].煤炭學(xué)報(bào),2012,37(9):1522-1527.

Xue Dongjie,Zhou Hongwei,Yi Haiyang.Fractal dimension estimate of rock crack curve in three-dimensional space[J].Journal of China Coal Society,2012,37(9):1522-1527.

[14] 劉曉輝,張 茹,劉建鋒.不同應(yīng)變率下煤巖沖擊動(dòng)力試驗(yàn)研究[J].煤炭學(xué)報(bào),2012,37(9):1528-1534.

Liu Xiaohui,Zhang Ru,Liu Jianfeng.Dynamic test study of coal rock under different strain rates[J].Journal of China Coal Society, 2012,37(9):1528-1534.

[15] 楊圣奇.裂隙巖石力學(xué)特性研究及時(shí)間效應(yīng)分析[M].北京:科學(xué)出版社,2011:86-128.

[16] Evert Hoek.Reliability of Hoek-Brown estimates of rock mass properties and their impact on design[J].International Journal of Rock Mechanics and Mining Sciences,1998,35(1):63-68.

[17] Fang Z,Harrison J P.A mechanical degradation index for rock[J].International Journal of Rock Mechanics and Mining Sciences, 2001,38(8):1193-1199.

[18] Behr V de,Hornick J L,Cabaraux J F.Growth patterns of Belgian Blue replacement heifers and growing males in commercial farms [J].Livestock Production Science,2001,71:121-130.

[19] Narushin V G,Takma C.Sigmoid model for the evaluation of growth and production curves in laying hens[J].Biosystems Engineering, 2003,84(3):343-348.

Meso-simulation and fractal characteristics for uniaxial compression test of inhomogeneous rock

LI Shou-ju1,LI De1,WU Li2,CAO Li-juan3
(1.State Key Laboratory of Structural Analysis for Industrial Equipment,Dalian University of Technology,Dalian 116024,China;2.Institute of Mechanical Engineering,Dalian Jiaotong University,Dalian 116028,China;3.School of Mechanical and Power Engineering,Dalian Ocean University,Dalian 116023, China)

In order to simulate fracture process of inhomogeneous rock specimen in meso-scale under uniaxial compression test,a new numerical simulation procedure based on finite element method was proposed.The simulation procedure was based on Mohr-Coulomb fracture criterion.Taking rock specimen as inhomogeneous and quasi-brittle material,the shear strength parameters of rock element were subjected to normal distribution.The FEM analysis investigation shows that inhomogeneity of rock specimen has a significant influence on the simulated stress-strain curve and the simulated curve agrees well with the experimental curve.Based on fractal theory,the fractal dimension of rock fracture for different displacement was calculated by using box-counting dimension method.The nonlinear relationship between fractal dimension and displacement at the top was proposed.The fractal analysis investigation show that the fracture of rock specimen has fractal characteristics.Fractal dimension can be used to quantitatively characterize the fracture process of rock specimen under uniaxial compression test.

TD315

A

0253-9993(2014)05-0849-06

李守巨,李 德,武 力,等.非均質(zhì)巖石單軸壓縮試驗(yàn)破壞過(guò)程細(xì)觀模擬及分形特性[J].煤炭學(xué)報(bào),2014,39(5):849-854.

10.13225/j.cnki.jccs.2013.0673

Li Shouju,Li De,Wu Li,et al.Meso-simulation and fractal characteristics for uniaxial compression test of inhomogeneous rock[J].Journal of China Coal Society,2014,39(5):849-854.doi:10.13225/j.cnki.jccs.2013.0673

2013-05-20 責(zé)任編輯:常 琛

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973)資助項(xiàng)目(2013CB035402);國(guó)家自然科學(xué)基金資助項(xiàng)目(51105048);中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(DUT13LK14)

李守巨(1960—),男,遼寧沈陽(yáng)人,教授,博士。E-mail:lishouju@dlut.edu.cn