班力壬,候宇航,杜偉升,俞 縉,戚承志, ,單仁亮
(1.北京建筑大學(xué) 土木與交通工程學(xué)院,北京 100044;2.煤炭科學(xué)研究總院有限公司 深地科學(xué)院,北京 100013;3.華僑大學(xué) 巖土工程研究所,福建 廈門 361021;4.中國礦業(yè)大學(xué)(北京)力學(xué)與建筑工程學(xué)院,北京 100083)
節(jié)理裂隙是控制巖體的變形破壞不可忽視的因素,節(jié)理的剪切力學(xué)特性比完整巖體的剪切力學(xué)特性更應(yīng)引起工程關(guān)注[1-2]。一般而言,影響巖石節(jié)理剪切力學(xué)特性的因素包括節(jié)理的表面形貌特征、法向應(yīng)力、節(jié)理面的接觸狀態(tài)以及充填物質(zhì)等。對于吻合程度較好、不含充填物質(zhì)的巖石節(jié)理而言,節(jié)理表面的形貌粗糙特征和法向應(yīng)力則是影響其剪切力學(xué)行為的主要因素[3-4]。
節(jié)理粗糙度作為影響剪切強度的關(guān)鍵因素,其表征方法一直都是巖石力學(xué)領(lǐng)域的熱點和難點[5]。粗糙度的定量表征方法包括統(tǒng)計參數(shù)法和分形法。常見的統(tǒng)計參數(shù)有:一階導(dǎo)數(shù)均方根Z2、結(jié)構(gòu)函數(shù)SF、跡線長度Rp、均方根RMS、均方值VMS、中心線平均值VCLA等[6]。許多統(tǒng)計參數(shù)雖然表征了粗糙度的各向異性特性,然而大部分參數(shù)均基于二維輪廓線研究節(jié)理表面的起伏度。采用分形理論表征的粗糙度公式很多,其優(yōu)勢在于粗糙度的表征可以獨立于測量尺度[7-8]。但有學(xué)者認為分形理論獲得的粗糙度的公式間通用性較差,分形法沒有體現(xiàn)節(jié)理粗糙度的三維特性和各向異性特征[9]。這限制了傳統(tǒng)統(tǒng)計參數(shù)和分形理論在節(jié)理剪切理論中的應(yīng)用和發(fā)展。
基于粗糙度的表征方法建立的節(jié)理面峰值抗剪強度模型可大致分為3 類:Barton 系模型、Grasselli系模型、其他類型模型。Barton 系模型主要以Barton 提出的JRC-JCS 模型為基礎(chǔ),其研究的重點在于如何定量、精確地確定節(jié)理的粗糙度參數(shù)JRC[10]。盡管近年來數(shù)字化處理技術(shù)[11]、小波變換理論[12]等方法在對節(jié)理粗糙度的表征和JRC 的獲取方面做了新突破,但獲取JRC 最常用的方法為:一是與10 條標(biāo)準(zhǔn)輪廓線對比,以主觀比較的方式確定JRC;二是通過直剪試驗得到的峰值抗剪強度反算JRC。第1 種方法評估粗糙度方式過于主觀,第2 種方法則顛倒了工程應(yīng)用的主次,即研究的目的在于預(yù)測峰值剪強度而非通過強度反分析粗糙度。當(dāng)法向應(yīng)力趨于0 時,Barton 系模型失去了數(shù)學(xué)意義。此外,JRC-JCS 模型最突出的問題在于只能截取若干二維輪廓線估算節(jié)理面的二維形貌,這種平均化思想過度簡化了節(jié)理的粗糙度信息,無法充分反映節(jié)理形貌的三維特征。
Grasselli 系的節(jié)理峰值抗剪強度模型主要基于其提出的三維粗糙度指標(biāo)體系。GRASSELLI 等[13]通過激光掃描技術(shù)獲得了節(jié)理表面的三維點云,將節(jié)理面粗糙不平的起伏近似為三角形的微凸體,運用三角形算法擬合了節(jié)理面的有效視傾角分布函數(shù),并提出最大接觸面積比A0、最大有效視傾角,粗糙度分布參數(shù)C等3 個函數(shù)作為表面形貌參數(shù)來預(yù)測和估算節(jié)理面的峰值抗剪強度。YANG 模型[14]、XIA 模型[15]、TATONE 模型[16]均采用了GRASSELLI 提出的3 個表面形貌參數(shù),TIAN 模型[17]將GRASSELLI提出的粗糙度分布參數(shù)C修正為C′,上述模型均可視為GRASSELLI 系模型。
除了BARTON 系模型和GRASSELLI 系模型外,學(xué)者們還提出了其他類型模型[18-20]。葛云峰[18]提出了一個節(jié)理粗糙度的經(jīng)驗?zāi)P?,其中的粗糙度參?shù)BAP(光亮面積百分比)與BARTON 提出的JRC 的相關(guān)性較差。ZHANG 等[19]提出的模型能較好反映節(jié)理在正向和反向剪切方向上的各向異性特性,且峰值剪脹角符合節(jié)理剪切的邊界條件,然其中的2 個粗糙度參數(shù)修正的一階導(dǎo)數(shù)均方根Z2與退化參數(shù)Cm是2個獨立參數(shù),其形式復(fù)雜且不易獲取,預(yù)測精度未得到廣泛證實。GHAZVINIAN 模型[20]只能通過直剪試驗反算初始剪脹角i0,因而不能預(yù)測節(jié)理的峰值剪切強度,在實際中應(yīng)用中還存在不小局限性。而由GRASSELLI 等[21]提出及后續(xù)發(fā)展的系列模型,由一組關(guān)鍵的形貌參數(shù)表征節(jié)理表面的粗糙度信息,兼顧對三維形態(tài)特征和節(jié)理各向異性特性的考慮,且預(yù)測的峰值抗剪強度精度較高,是目前較為合適的節(jié)理峰值抗剪強度模型。
節(jié)理在剪切過程中的損傷退化機理也在逐步被揭示。在建立節(jié)理峰值抗剪強度模型時,一般需要考慮節(jié)理在剪切方向上的各向異性特性,包括在相反方向上的各向異性(在某一方向和該方向旋轉(zhuǎn)180°對應(yīng)的方向)。此外,傳統(tǒng)的JRC-JCS 模型以JCS(壁面抗壓強度,對于未風(fēng)化新鮮巖石等于巖石的單軸抗壓強度)作為材料強度解釋直剪試驗過程中微凸體的受壓破壞尚存在不足。近年來XIA 等[15]、ZHANG 等[19]通過試驗證實,節(jié)理面表面微凸體多發(fā)生拉伸破壞而非受壓破壞,因此峰值剪切強度模型中的材料強度項應(yīng)為節(jié)理面壁抗拉強度。GRASSELLI 等[13]指出,在粗糙節(jié)理剪切過程中,節(jié)理的滑移、啃斷等只是一小部分接觸區(qū)域引起的,接觸面上面向剪切方向、最陡的區(qū)域與節(jié)理的損傷破壞區(qū)域及峰值剪切強度有密切聯(lián)系。BAN 等[22]基于GRASSELLI 形貌分布理論與赫茲接觸理論,提出考慮實際接觸節(jié)理粗糙度的峰值剪脹角模型。這些機理的揭示促進了節(jié)理峰值抗剪強度模型的發(fā)展。
然而上述研究僅把實際接觸微凸體角度的平均值視為節(jié)理剪脹角,不同形貌的微元對節(jié)理剪切強度的貢獻并未考慮。為合理預(yù)測節(jié)理峰值剪脹角,筆者重點開展GRASSELLI 系列模型的研究,通過考慮實際接觸節(jié)理微凸體對剪切強度不同貢獻比例,推導(dǎo)出了整個節(jié)理微凸體等效實際接觸節(jié)理平均角。將節(jié)理微凸體等效實際接觸節(jié)理平均角等同于節(jié)理峰值剪脹角,提出新的節(jié)理峰值剪脹角模型。通過現(xiàn)有的研究成果論證了筆者所提峰值剪脹角的合理性,進一步揭示了節(jié)理抗剪行為的機理,為相關(guān)研究提供參考。
GRASSELLI[9]研究發(fā)現(xiàn)節(jié)理面的真實形貌可用三角形網(wǎng)格微元近似表示,三角形微元可通過節(jié)理面測量的數(shù)據(jù)點生成,而三角形微元的傾角θ、有效視傾角 θ?、剪切方向與微元傾向的夾角α之間存在一定的幾何關(guān)系,如圖1 所示。
圖1 剪切面上三角形微元受剪示意[9]Fig.1 Shear state of the triangle asperity on the shear plane[9]
通過統(tǒng)計數(shù)據(jù)擬合,發(fā)現(xiàn)三角形微元有效視傾角大于 θ?的微元面積總和與節(jié)理面面積比Aθ?與 θ?的關(guān)系見式(1)[9]:
其中,A0為最大接觸面積比,是所有節(jié)理微元等效傾角大于0°時的面積總和與節(jié)理表面面積總和之比;C為粗糙度分布參數(shù),描述視傾角的分布情況。當(dāng)C>1 時,式(1)中的函數(shù)關(guān)系可用圖2 表示。
圖2 GRASSELLI 視傾角分布函數(shù)中 Aθ? 與 θ?的關(guān)系Fig.2 Relationship between Aθ? and θ? in GRASSELLI’s distribution function of apparent dip angle
由圖2 可知,隨著有效視傾角 θ?的增大,有效視傾角大于 θ?的所有微元面積比Aθ?在減小,最終減小到0。
式中,τp為節(jié)理峰值抗剪強度;σn為法向應(yīng)力;φb為基本摩擦角;φr為殘余摩擦角;β為節(jié)理面的傾角(一般β=0);σt為單軸抗拉強度。
TATONE 等[23]計算了圖2 曲線下方的閉合面積,對式(1)進行積分,得
式中,l為采樣間距,mm。
XIA 等[15]認為GRASSELLI 模型與TATONE 模型不滿足摩爾-庫倫形式,因此提出式(5)來預(yù)測節(jié)理峰值抗剪強度:
YANG 等[14]提出節(jié)理峰值抗剪強度模型,見式(6):
式中,σc為單軸抗壓強度。
BAN 等[22]基于GRASSELLI 形貌分布理論與赫茲接觸理論,提出考慮實際接觸節(jié)理粗糙度的峰值剪脹角模型,見式(7):
式中,ip為峰值剪脹角。
BAN 等[22]提出峰值剪脹角即為實際接觸節(jié)理微元的平均角度,然而不同形貌的微元對節(jié)理剪切強度的貢獻并未考慮。
節(jié)理的峰值抗剪強度符合摩爾-庫倫準(zhǔn)則[25],見式(8):
將式(8)與既有的GRASSELLI 系列模型聯(lián)立可反推出各模型峰值剪脹角。
陳曦等[25]研究發(fā)現(xiàn)節(jié)理在剪切時受到的阻力只是一小部分較陡的微凸體產(chǎn)生的,式(1)中C的增大會導(dǎo)致較陡的視傾角面積比含量降低,節(jié)理粗糙度也會降低。當(dāng)C增大,節(jié)理中較陡視傾角的減少才是引起節(jié)理粗糙度降低的根本原因,而非較平緩視傾角的增大。因此,對節(jié)理粗糙度的研究應(yīng)當(dāng)重點關(guān)注較陡的視傾角而非較平緩的視傾角。
對于巖石節(jié)理在法向應(yīng)力下的剪切問題,實際接觸面積Ac可以近似認為是作用在節(jié)理面上法向荷載N除以巖石的單軸抗壓強度σc[22],見式(9):
式中,A為節(jié)理面名義面積。
式中,Aθ?cr/A0為面向剪切方向接觸面積與面向剪切方向所有節(jié)理微凸體面積的比,約等于Ac/A[22]。
基于以上假設(shè),結(jié)合式(9)、(10),得到所有接觸微凸體中的最小傾角,有
由式(11)可確定實際接觸節(jié)理的最小角度。剪切過程中實際接觸的微凸體才對剪切強度有貢獻,那么對剪切強度有貢獻的部分為視傾角在的微凸體。對于視傾角在(0,)的微凸體在研究節(jié)理剪切強度時應(yīng)該予以剔除。
KWON 等[26]研究了長方體微凸體的破壞模式,如圖3 所示,長方體微凸體高度為h,長度為b,幾何參數(shù)m=h/b,法向荷載為N,切向荷載為T。當(dāng)m
圖3 長方體微凸體破壞模式Fig.3 Failure modes of rectangular-shaped asperity
臨界幾何參量mc如式(12)[26]所示:
式中,c、φf分別為完整巖石的黏聚力與峰值摩擦角。
單個長方體微凸體的剪切強度[26]如式(13)所示。由式(13)可知,對于確定材料的巖石節(jié)理微凸體,當(dāng)m小于臨界幾何參數(shù)mc時,剪切強度隨著m增加線性增加,達到臨界幾何參數(shù)后保持不變。
上述是將節(jié)理微凸體等效為長方形微凸體進行分析,對于微小的單元等效為長方形單元與三角形單元是相似的。因此將節(jié)理等效為三角形微凸體也應(yīng)具有以上類似性質(zhì)。假設(shè)微凸體的臨界形狀參數(shù)mc對應(yīng)的傾角為臨界角度,則對于確定材料的巖石節(jié)理三角形微凸體,當(dāng)?shù)刃A角 θ?小于臨界角度時,剪切強度隨著 θ?增加線性增加,即剪切強度正比于 θ?;等效傾角 θ?大于臨界角度時,剪切強度不變,正比于。
綜合2.1 節(jié),在巖石節(jié)理面上抵抗剪切的阻力是由面向剪切方向上所接觸微凸體產(chǎn)生的,所接觸微凸體最小角度為。對于 θ?在的微凸體,其破壞模式為剪脹破壞,對剪切強度的貢獻正比于 θ?;對于 θ?在的微凸體,其破壞模式為剪斷破壞,對剪切強度的貢獻正比于。
陳曦等[25]研究表明GRASSELLI 分布函數(shù)中的C與能夠描述有效視傾角的分布情況和“含量”情況,類似于土壤顆粒粒度分析過程。將GRASSELLI分布函數(shù)進行變式,將A0移到等式左側(cè)得
圖4 與θ?的關(guān)系Fig.4 Relationship between and θ?
故所有剪脹破壞微凸體的剪切強度τ1,可由式(16)求得
對于 θ?在間的微凸體,其破壞模式為剪斷破壞,對剪切強度的貢獻正比于。則所有剪斷破壞的三角形微凸體強度正比于與所有等效傾角為θ?的三角形單元含量乘積之和。所有剪斷破壞微凸體的剪切強度τ2,可由式(17)求得
則所有面向剪切方向微凸體的抗剪強度τ3,如式(18)所示:
BAN 等[22]發(fā)現(xiàn)實際接觸節(jié)理微凸體的平均角度為節(jié)理峰值剪脹角。然而上述研究僅從實際接觸節(jié)理形貌分布角度來提出峰值剪脹角,節(jié)理峰值剪脹角是為預(yù)測節(jié)理剪切強度服務(wù)的,峰值剪脹角應(yīng)該從物理意義上與節(jié)理剪切強度聯(lián)系起來。本文考慮實際接觸節(jié)理微凸對剪切強度的不同貢獻,提出了一個可與剪切強度相聯(lián)系的等效實際接觸節(jié)理平均角;實際接觸節(jié)理平均傾角即為剪脹角,所提等效實際接觸節(jié)理平均角即為新的峰值剪脹角模型,如式(19)所示。
表1 各學(xué)者試驗數(shù)據(jù)以及模型計算結(jié)果Table 1 Test data and calculation results of different scholars
為對各個模型進行評價,選取相對誤差平均值來定量描述各模型的精度。
相對誤差平均值δ的計算,如式(20)所示:
經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為14%,BAN 模型平均誤差為14%,GRASSELLI 模型估算的平均誤差為17%,XIA 模型估算的相對平均誤差為13%,YANG 模型估算的相對平均誤差為21%,TATONE 模型估算的相對平均誤差為14%。
鑒于節(jié)理試樣的巖性對剪切行為的影響,下面用不同學(xué)者們?yōu)椴煌瑤r性的節(jié)理試樣實施剪切試驗的結(jié)果來驗證本模型。
根據(jù)YANG 等[14]發(fā)表的10 組花崗巖試驗結(jié)果,采用最小二乘法擬合得到=28°。表1 展示了YANG 的試驗結(jié)果與上述6 種模型的預(yù)測結(jié)果對比。經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為3%,BAN 模型平均誤差為8%,GRASSELLI模型估算的平均誤差為13%,XIA 模型估算的相對平均誤差為13%,YANG 模型估算的相對平均誤差為25%,TATONE 模型估算的相對平均誤差為12%。
根據(jù)YANG 等[14]發(fā)表的10 組砂巖試驗結(jié)果,采用最小二乘法擬合得到=30°。表1 展示了YANG的試驗結(jié)果與上述6 種模型的預(yù)測結(jié)果對比。經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為6%,BAN 模型平均誤差為11%,GRASSELLI 模型估算的平均誤差為24%,XIA 模型估算的相對平均誤差為22%,YANG 模型估算的相對平均誤差為8%,TATONE 模型估算的相對平均誤差為26%。
根據(jù)XIA 等[15]發(fā)表的10 組砂巖試驗結(jié)果,采用最小二乘法擬合得到=40°。表1 展示了XIA 的試驗結(jié)果與上述6 種模型的預(yù)測結(jié)果對比。經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為14%,BAN 模型平均誤差為14%,GRASSELLI 模型估算的平均誤差為19%,XIA 模型估算的相對平均誤差為13%,YANG 模型估算的相對平均誤差為16%,TATONE 模型估算的相對平均誤差為19%。根據(jù)TATONE 等[23]發(fā)表的6 組砂巖試驗結(jié)果,采用最小二乘法擬合得到=40°。表1 展示了TATONE 的試驗結(jié)果與上述6 種模型的預(yù)測結(jié)果對比。經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為4%,BAN 模型平均誤差為4%,GRASSELLI 模型估算的平均誤差為8%,XIA 模型估算的相對平均誤差為11%,YANG 模型估算的相對平均誤差為22%,TATONE 模型估算的相對平均誤差為6%。
根據(jù)GRASSELLI[9]發(fā)表的7 組石灰?guī)r節(jié)理試驗結(jié)果,采用最小二乘法擬合得到=28°。表1 展示了GRASSELLI 試驗的試驗結(jié)果與上述6 種模型的預(yù)測結(jié)果對比。經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為4%,BAN 模型平均誤差為15%,GRASSELLI 模型估算的平均誤差為5%,XIA模型估算的相對平均誤差為9%,YANG 模型估算的相對平均誤差為6%,TATONE 模型估算的相對平均誤差為5%。
根據(jù)GRASSELLI[9]發(fā)表的7 組花崗巖節(jié)理試驗結(jié)果,采用最小二乘法擬合得到=35°。表1 展示了GRASSELLI 的試驗結(jié)果與上述6 種模型的預(yù)測結(jié)果對比。經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為5%,BAN 模型平均誤差為23%,GRASSELLI 模型估算的平均誤差為8%,XIA 模型估算的相對平均誤差為5%,YANG 模型估算的相對平均誤差為8%,TATONE 模型估算的相對平均誤差為6%。
根據(jù)Grasselli[9]發(fā)表的11 組大理巖節(jié)理試驗結(jié)果,采用最小二乘法擬合得到=28°。表1 展示了GRASSELLI 的試驗結(jié)果與上述6 種模型的預(yù)測結(jié)果對比。經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為12%,BAN 模型平均誤差為14%,GRASSELLI 模型估算的平均誤差為16%,XIA 模型估算的相對平均誤差為11%,YANG 模型估算的相對平均誤差為11%,TATONE 模型估算的相對平均誤差為16%。
根據(jù)GRASSELLI[9]發(fā)表的3 組節(jié)理試驗結(jié)果,采用最小二乘法擬合得到=28°。表1 展示了Grasselli的試驗結(jié)果與上述6 種模型的預(yù)測結(jié)果對比。經(jīng)計算本文提出的峰值剪脹角新模型估算的相對平均誤差為21%,BAN 模型平均誤差為23%,GRASSELLI模型估算的平均誤差為66%,XIA 模型估算的相對平均誤差為40%,YANG 模型估算的相對平均誤差為22%,TATONE 模型估算的相對平均誤差為9%。
將BAN 等[22]、GRASSELLI[9]、XIA 等[15]、YANG等[14]、TATONE 等[23]的89 組試驗數(shù)據(jù)匯總分析平均誤差,如圖5 所示(圖5 與圖6 中斜直線為數(shù)據(jù)擬合直線,直線越接近y=x,表明模型精度越高)。本文提出的峰值剪脹角新模型估算的相對平均誤差為10%,BAN 模型估計的誤差為13%,GRASSELLI 模型估算的平均誤差為17%,XIA 模型估算的相對平均誤差為13%,YANG 模型估算的相對平均誤差為13%,TATONE 模型估算的相對平均誤差為14%。由此可知在數(shù)據(jù)比較多時本文模型精度屬于最好的。
圖5 所得節(jié)理峰值剪脹角試驗結(jié)果與本模型計算結(jié)果對比Fig.5 Comparison of test results and calculation results of the peak shear strength model
圖6 =32°時所得節(jié)理峰值剪脹角試驗結(jié)果與本模型計算結(jié)果對比Fig.6 Comparison of test results and calculation results of the peak shear strength model with =32°
根據(jù)第4 節(jié)驗證結(jié)果,對于部分試驗數(shù)據(jù),BAN模型與新模型預(yù)測精度相差不大;而對于另外的試驗結(jié)果,新模型預(yù)測結(jié)果明顯優(yōu)于BAN 模型。為解釋上述情況,筆者探討不同C、對的影響。當(dāng)=70°,σn=1 MPa,A0=0.5,σc=27.5 MPa,取C=3、5、7 時,分析隨著變化,的變化,如圖7 所示。
圖7 與預(yù)測值關(guān)系Fig.7 Relationship between and
圖7 中實線為BAN 模型預(yù)測值,點劃線為新模型預(yù)測值。由圖7 可知,當(dāng)C一定時,隨著增大,以類似拋物線的形式逐漸變大,并且趨于BAN 模型預(yù)測值;隨著C的增大,BAN 模型與新模型的預(yù)測值都在降低。其中BAN 模型由39°減小為27°、21°。當(dāng)=32°時,新模型由32°減小為26°、20°。這是因為C的增大會導(dǎo)致較陡的視傾角的面積比含量降低,節(jié)理粗糙度也會降低,進而節(jié)理峰值剪脹角減小。同時由圖7 可知,隨著C的增大,2 個模型的差值也會降低。當(dāng)=32°的情況下,C=3 時,兩模型的差值為8°;C=5 時,兩模型的差值為4°;而C=7 時,兩模型的差值為1°。這表明,在C較大時,兩模型的預(yù)測結(jié)果較為接近;而C較小時,新模型的預(yù)測精度要高于BAN 模型。因此是否考慮對于不同形貌的節(jié)理有著不同的影響。
由圖7 可知,當(dāng)C較小時,對模型的預(yù)測結(jié)果影響較大。通過89 組實驗數(shù)據(jù)的擬合結(jié)果,筆者確定了約為32°。根據(jù)共89 組試驗數(shù)據(jù)匯總分析平均誤差僅為11%,與其他模型比較所得結(jié)果精度最高。上述僅從試驗結(jié)果擬合確定了,并未從物理機理上去解釋。因此有必要對取值在理論上的解釋進行一些闡述。
對于長方體微凸體,由式(12)可知臨界幾何參量與巖石材料性質(zhì)、法向應(yīng)力有關(guān)。在低法向應(yīng)力下[27],式(12)可進行簡化,即
由式(21)可知,mc僅與峰值內(nèi)摩擦角關(guān)系較大。式(12)是節(jié)理形貌由連續(xù)長方體微凸體劃分所得的理論計算結(jié)果。如采用三角形劃分節(jié)理形貌,將微凸體由長方體轉(zhuǎn)變?yōu)槿切挝⑼贵w,則可將微凸體最高點等效為三角形微凸體頂點,劃分間距為三角形的最底部長度。圖8 表示微凸體為剪斷破壞模式時長方體與三角形微凸體的劃分情況,綠色實體部分為劃分的三角形微凸體。
圖8 采用長方形與三角形微凸體劃分網(wǎng)格時的示意Fig.8 Schematic diagram of grid generation with rectangular and triangular asperities
對于長方形微凸體,其臨界幾何參數(shù)為mc,破壞面沿微凸體的根部。由圖8 可知,將節(jié)理面等效為連續(xù)三角形微凸體時,三角形臨界角度如式(22)所示。
常見巖石內(nèi)摩擦角通常為21°~40°,可取30°[28]。代入式(22)可得為25°~35°,通常為30°左右。因此,由89 組試驗結(jié)果擬合獲得的=32°有一定合理性。
新模型的應(yīng)用需要獲取一些節(jié)理微觀參數(shù),對于一些工程中這些參數(shù)的獲取較為困難。以上針對上述情形對模型進行了簡化,取試驗擬合的=32°來統(tǒng)一確定新的峰值剪脹角模型,上述對模型的簡化的研究可以拓展新模型的適用范圍。
(1)在巖石節(jié)理上抵抗剪切的阻力是由面向剪切方向的所接觸微凸體產(chǎn)生的。對于 θ?在(,)的微凸體,其破壞模式為剪脹破壞,對剪切強度的貢獻正比于 θ?。而 θ?在的微凸體的破壞模式為剪斷破壞,對剪切強度的貢獻正比于。
(2)根據(jù)實際接觸節(jié)理微凸對剪切強度不同貢獻比例,推導(dǎo)出了等效實際接觸節(jié)理平均角。將等效實際接觸節(jié)理平均角等同于節(jié)理峰值剪脹角,提出新的節(jié)理峰值剪脹角模型。該指標(biāo)物理意義明確,TATONE與BAN 所提粗糙度指標(biāo)僅僅是本文所提新指標(biāo)的一種特殊情況。
(3)用以往學(xué)者的89 組試驗數(shù)據(jù)對分析了模型的預(yù)測能力。本文提出的峰值剪脹角新模型估算的相對平均誤差為10%,BAN 模型估計的誤差為13%,GRASSELLI 模型估算的平均誤差為17%,XIA 模型估算的相對平均誤差為13%,YANG 模型估算的相對平均誤差為13%,TATONE 模型估算的相對平均誤差為14%。在數(shù)據(jù)較多時,本文模型精度較好。