(中南大學 資源與安全工程學院，湖南 長沙，410083)

工程實踐表明，天然巖體中含有大量的節(jié)理，有地殼運動所產生的地應力作用形成的節(jié)理，有風化作用形成的節(jié)理，還有人為因素(爆破、打擊等)形成的節(jié)理。節(jié)理破壞了巖體的完整性，大大降低了巖體的力學性能。因此，節(jié)理的力學性質直接影響巖體的力學性質，進而影響巖體工程的穩(wěn)定與安全。對節(jié)理、節(jié)理巖體的力學及相關性質的研究，具有極大的工程實踐意義。節(jié)理的抗剪強度遠小于巖石的抗壓強度，對節(jié)理力學性質的研究主要是對其抗剪性能的研究，眾多學者研究表明，節(jié)理的剪切強度與節(jié)理表面形貌特征有較大關系。Carr等[1?2]運用分形理論來描述復雜的節(jié)理表面形貌特征，計算了節(jié)理表面的分形參數(shù)。杜守繼等[3]研究了節(jié)理的力學性質與分形參數(shù)之間的關系。周枝華等[4?5]運用統(tǒng)計學原理分析了節(jié)理表面特征參數(shù)的統(tǒng)計分布規(guī)律。夏才初等[6?9]定義了節(jié)理表面形貌參數(shù)，定量描述節(jié)理表面的形貌特征。Hutson等[10?14]研究了剪切作用下節(jié)理表面的形貌特征變化，并定義了變化度量參數(shù)，同時建立了相應的函數(shù)式。曹平等[15]研究了節(jié)理表面形貌在剪切作用下的變化規(guī)律，發(fā)現(xiàn)節(jié)理表面輪廓角的加權平均值隨著剪切次數(shù)的增加而減小。杜守繼等[16]首次從能量的角度簡略分析了節(jié)理抗剪強度與表面損傷的關系。但是上述學者沒有分析節(jié)理抗剪強度與節(jié)理體積損失的關系。在眾多學者研究的基礎上，本文作者運用先進的三維非接觸式高精度激光形貌儀，對剪切前后節(jié)理表面進行激光掃描測試，從節(jié)理表面摩擦的角度，研究了剪切作用下節(jié)理體積損失變化規(guī)律，分析了抗剪強度與體積損失的關系。

1 節(jié)理體積計算公式

在計算體積之前，首先設定一個水平面M，將節(jié)理在水平面M以上部分稱為節(jié)理體積(Joint volume)。在計算整個節(jié)理體積時，以與最低谷點相切的平面為水平基準面M，以保證節(jié)理表面所有凸起體的體積都包含在計算范圍內，水平面M的高度可以隨著節(jié)理表面凸起體體積計算區(qū)域進行調整，節(jié)理表面某點到水平基準面M的垂直距離稱為節(jié)理高度h。為了計算節(jié)理在水平面M以上部分的體積，提出節(jié)理體積計算的積分式：

式中：Vh為節(jié)理高度為h時水平面M以上部分的節(jié)理體積；An為節(jié)理表面在水平面M上的投影面積；z(x,y)為節(jié)理表面的高度函數(shù)；dx和dy分別為An范圍內在x和y方向上的微分。

實際的節(jié)理表面是一個極為復雜的曲面，用數(shù)學函數(shù)式來描述這個曲面具有很大的困難，且節(jié)理表面特征具有多變性，不同的節(jié)理具有不同的表面形貌，需要用不同的數(shù)學函數(shù)式來表達，即z(x,y)的表達式難以求出，采用式(1)計算節(jié)理體積具有相當大的難度。因此，采用離散形式估算節(jié)理體積是一種實用的方法。

將節(jié)理表面水平投影面積An分別沿x和y方向劃分為n個邊長為ixΔ，jyΔ的矩形區(qū)域，以該矩形為底，向上延伸成一個上表面不規(guī)則的近似棱柱體。為了便于計算，將這個棱柱體的上表面簡化成一個平面，如圖1所示。為了盡可能求得接近棱柱體的真實體積，以式(2)求得棱柱體的平均高zΔ：

式中，z1，z2，z3，z4分別表示所劃分計算區(qū)域的近似棱柱體4個角點的高度。

圖1 計算區(qū)域示意圖Fig.1 Schematic diagram of computational domain

Δxi和Δyj分別為計算區(qū)域在x和y方向上的長度，這也是節(jié)理體積的最小計算單位；i，j分別表示x，y方向上第i，j個計算區(qū)域，那么計算區(qū)域棱柱體體積Δvi j的計算公式如下：

在An范圍內，對圖1所示的所有計算區(qū)域棱柱體的體積求和得

式中：Ni為在An范圍內x方向劃分小矩形個數(shù)總和；Nj為在An范圍內y方向上劃分小矩形個數(shù)總和。

在TalyMap5.0中，將ixΔ和jyΔ設定為單位長度。因此，上文所劃分的計算區(qū)域為單位區(qū)域，單位區(qū)域棱柱體體積與其高度直接相關，式(4)簡化為：

2 舉例分析

2.1 掃描儀器介紹

Talysurf CLI 2000三維非接觸式高精度激光形貌儀由英國Taylor Hobson公司研制生產。該儀器的主要功能是測試物件表面形態(tài)，其外觀尺寸(長×寬×高)為800 mm×600 mm×970 mm，質量為280 kg，主要由激光掃描測量裝置、控制單元、數(shù)據(jù)采集軟件Talysurf CLI、計算機、數(shù)據(jù)分析軟件 TalyMap5.0組成[15]。

2.2 用TalyMap5.0計算節(jié)理體積

圖2所示為節(jié)理表面某個區(qū)域經過放大處理的三維圖，節(jié)理的最大高度為1.2 mm。從圖2可見：整個節(jié)理表面由微小凸起體組成。節(jié)理表面受到摩擦時，即這些凸起體的磨損破壞。因此，用TalyMap5.0計算節(jié)理體積具有很高精度。

圖2 節(jié)理表面三維圖Fig.2 3D image of joint surface

節(jié)理體積計算示意圖如圖3所示。在TalyMap5.0中，右邊工具條中的箭頭可以隨鼠標上下移動，移動范圍為0～hmax。箭頭的上下移動代表示水平面M的上下移動，當箭頭停止在某個高度時，水平面以上部分的體積和最大高度等參數(shù)能夠計算出來。同時，可以在節(jié)理表面最大高度范圍內輸入數(shù)值來精確確定水平面M的高度，從而計算節(jié)理在該高度以上部分的體積。從圖3可知：節(jié)理表面在中間范圍較低，是凹谷部分，左右兩端較高，分布著大量凸起體。同時也是容易受到磨損，提供抗剪強度的部分。高度越大，相對應的體積投影面積就越少。利用該功能可以計算和分析節(jié)理在不同高度范圍內的體積變化以及凸起體磨損情況。

2.3 實例分析

在節(jié)理剪切試驗過程中發(fā)現(xiàn)，如果節(jié)理上下盤耦合較差，且法向應力較低時，節(jié)理表面磨損較少，且節(jié)理的抗剪強度較小，節(jié)理的抗剪強度曲線沒有明顯的峰值抗剪強度。當節(jié)理上下盤耦合良好時，節(jié)理表面有明顯的摩擦痕跡，且抗剪強度曲線有明顯的峰值抗剪強度。本文選取節(jié)理上下盤耦合良好的情況予以分析。

圖 4所示為示例節(jié)理的三維圖。該節(jié)理在x，y方向的長度和寬度分別為72.0 mm和51.5 mm，在z方向的最大高度隨著剪切作用發(fā)生了變化(如表 1所示)。表2和3所示為在經過2次不同法向應力的剪切作用后，且每次剪切前后，用Talysurf CLI 2000掃描測試，并使用 TalyMap5.0中的體積分析功能計算的結果。

圖3 節(jié)理體積計算示意圖Fig.3 Illustration of computational joint volume on horizontal plane

表1 節(jié)理面最大高度值Table 1 Maximum height of joint surface mm

圖4 計算節(jié)理的三維圖Fig.4 3D image of computational joint

從表1可知：經過2次不同法向應力的剪切試驗，節(jié)理最大高度減小，由剪切前的9.740 mm變?yōu)?.544 mm，這是節(jié)理體積變化的直接特征。在實際節(jié)理剪切過程中，由于剪脹作用的存在，受到磨損的主要是節(jié)理面上較高的凸起體，而較低的凸起體不那么容易受到摩擦。從表2可知：剪切作用下節(jié)理體積因磨損而減少，由剪切前的20 343.0 mm3減少到剪切1次后的17 240.0 mm3，共減少約3 103.0 mm3，平均高度減少約0.858 mm，最大高度減少1.196 mm，與試驗現(xiàn)象一致。第2次剪切后，節(jié)理表面最大高度只有8.544 mm，因此，在9 mm以上的節(jié)理體積為0 mm3。

表2 水平面在不同高度上的節(jié)理體積Table 2 Joint volume on horizontal plane at different heights mm3

表3 水平面在不同高度范圍內的節(jié)理體積Table 3 Joint volume between two horizontal planes at different heights mm3

從表3可知：剪切前，在0～1 mm，1～2 mm和2～3 mm高度范圍內，節(jié)理體積相差不大，說明這段高度范圍內節(jié)理表面的凸起體較少，所以體積變化不大。2次剪切后，在0～3 mm高度范圍內，節(jié)理體積變化很少，共減少140.5 mm3，再次說明0～3 mm高度范圍內節(jié)理表面凸起體較少。節(jié)理體積損失主要集中在高度為3～7 mm范圍內，由于最高凸起體的體積較小，損失不多，說明節(jié)理在耦合良好的情況下，主要是中間高度范圍內的凸起體受到摩擦。

計算水平面M在不同高度情況下節(jié)理體積的差，可以求得某一指定范圍內節(jié)理體積，同時可以分析該高度范圍內節(jié)理面受到的磨損情況，這對于材料摩擦特性的研究具有較大意義。圖5所示為水平面位于8 mm處的節(jié)理體積計算示意圖。由圖5可知：剪切作用下，同一節(jié)理在相同高度處的投影面積減少，這說明剪切作用使得節(jié)理表面的凸起體被剪斷，節(jié)理表面形貌發(fā)生了變化。從表2可得，剪切前，8 mm以上節(jié)理體積有 127.60 mm3，但是剪切后，只剩下 6.90 mm3，共損失約120.70 mm3，損失的體積導致8 mm以上節(jié)理面高度和投影面積的減小，水平面在 8.544 mm以上的節(jié)理體積完全被磨損掉。

圖5 水平面在8 mm處節(jié)理體積計算示意圖Fig.5 Illustration of computing joint volume on horizontal plane at height of 8 mm

3 結論

(1) 使用Talysurf CLI 2000掃描節(jié)理表面具有很高精度，用TalyMap5.0能夠準確快速計算指定高度的節(jié)理體積。

(2) 在剪切作用下，節(jié)理表面高度減小，體積也相應減小，不同法向應力的剪切作用下，相同高度上節(jié)理表面體積變化不一致。

(3) 節(jié)理上下盤耦合良好時，較高凸起體更易受到磨損，節(jié)理體積損失主要集中在中間高度范圍內。

[1] Carr J R, Warriner J B. Relationship between the fractal dimension and joint roughness coefficient[J]. Bulletion the Association of Engineering Geologist, 1989, 1151(2): 253?263.

[2] XIE He-ping, WANG Jin-an, XIE Wei-hong. Fractal effect of surface roughness on the mechanical behavior of rock joints[J].Chaos, Solitons & Fractals, 1997, 8(2): 221?252.

[3] 杜守繼, 孫鈞, 江崎哲郎, 等. 巖體不連續(xù)面的幾何特性與剪切強度的關系[J]. 同濟大學學報, 1999, 27(4): 379?383.DU Shou-ji, SUN Jun, Tetsuro E, et al. Relation of surface roughness and shear strength of rock joints[J]. Journal of Tongji University, 1999, 27(4): 379?383.

[4] 周枝華, 杜守繼. 巖石節(jié)理表面幾何特性的三維統(tǒng)計分析[J].巖土力學, 2005, 26(8): 1227?1232.ZHOU Zhi-hua, DU Shou-ji. 3D statistic analysis of geometrical properties of a rock joint[J]. Rock and Soil Mechanics, 2005,26(8): 1227?1232.

[5] 李慶松, 杜守繼. 巖石節(jié)理表面三維粗糙特性參數(shù)的定量研究[C]//中國巖石力學與工程學會第七次學術大會論文集. 北京: 科學技術出版社, 2002: 91?93.LI Qing-song, DU Shou-ji. Quantitative research of threedimensional parameters for rock joint surface roughness[C]//The 7th Academic Congress of Chinese Society for Rock Mechanics and Engineering. Beijing: China Science and Technology Press,2002: 91?93.

[6] 夏才初. 巖石結構面表面形態(tài)特征研究[J]. 工程地質學報,1996, 4(3): 71?78.XIA Cai-chu. A study on the surface morphological features of rock structure faces[J]. Journal of Engineering Geology, 1996,4(3): 71?78.

[7] 夏才初, 孫宗頎, 潘長良. 若干理想表面的形貌參數(shù)特征[J].中南礦冶學院學報, 1993, 24(6): 731?737.XIA Cai-chu, SUN Zong-qi, PAN Chang-liang. The parameters characteristics of topography of several ideal surfaces[J]. Journal of Central South Institute of Mining Metallurgy, 1993, 24(6):731?737.

[8] 夏才初, 孫宗頎. 工程巖體節(jié)理力學[M]. 上海: 同濟大學出版社, 2002: 100?201.XIA Cai-chu, SUN Zong-qi. Joint mechanics of engineering rock mass[M]. Shanghai: Tongji University Press, 2002: 100?201.

[9] Belem T, Homand F, Souley M. Quantitative parameters for rock joint surface roughness[J]. Rock Mechanics and Rock Engineering, 2000, 33(4): 217?242.

[10] Huston R W, Dowding C H. Joint asperity degradation during cyclic shear[J]. International Journal of Rock Mechanics and Mining Science and Geomechanics, 1990, 27(2): 109?119.

[11] Homand F, Belem T, Souley M. Friction and degradation of rock joint surfaces under shear loads[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2001,25(10): 973?999.

[12] Riss J, Gentier S, Sirieix C, et al. Degradation characterization of sheared joint wall surface morphology[C]//American Rock Mechanics Association, 2nd North American Rock Mechanics Symposium Rock Mechanics. Aubertin, Hassani, Mitri(eds).Balkema: Rotterdam, 1996: 1343?1349.

[13] Fox D J, Kana D D, Hsiung S M. Influence roughness on dynamic shear behaviour in jointed rock[J]. International Journal of Rock Mechanics and Mining Sciences, 1998, 35(7): 923?940.

[14] Lee H S, Park Y J, Cho T F, et al. Influence of asperity degradation on the mechanical behavior of rough rock joints under cyclic shear loading[J]. International Journal of Rock Mechanics and Mining Sciences, 2001, 38(7): 967?980.

[15] 曹平, 范祥, 蒲成志, 等. 節(jié)理剪切試驗及其表面形貌特征變化分析[J]. 巖石力學與工程學報, 2011, 30(3): 480?485.CAO Ping, FAN Xiang, PU Cheng-zhi, et al. Shear test of joint and analysis of morphology characteristic evolution of joint surface[J]. Chinese Journal of Rock Mechanics and Engineering,2011, 30(3): 480?485.

[16] 杜守繼, 江崎哲郎, 蔣宇靜, 等. 巖石節(jié)理的剪切破壞特性的定量分析[C]//中國巖石力學與工程學會第五次學術大會論文集. 北京: 中國科學技術出版社, 1998: 239?244.DU Shou-ji, Tetsuro E, JIANG Yu-jing, et al. The quantitative evaluation of surface damage of rock joint during shear processes[C]//The 7th Academic Congress of Chinese Society for Rock Mechanics and Engineering. Beijing: China Science and Technology Press, 1998: 239?244.