劉彥輝

摘要： 滾石的運動特征研究是滾石災(zāi)害防治的重要前提，而坡表覆蓋層和滾石自身特性是影響其運動特征的重要因素。以斜坡覆蓋層材質(zhì)為變量，開展物理模型試驗，研究了混凝土坡面和土質(zhì)坡面下滾石試塊的運動特征;并在此基礎(chǔ)上，以在建合壁津高速公路邊坡為例，開展?jié)L石運動特征數(shù)值模擬研究，以運動速度、轉(zhuǎn)動速度、運動能量等為指標(biāo)分析了滾石形狀和尺寸對其運動過程的影響。結(jié)果表明：形狀和尺寸對滾石運動特征影響較大，滾石越接近球形、尺寸越小，其運動速度越大、轉(zhuǎn)動速度越大;滾石在下落過程中有30%以上的勢能通過坡面碰撞損耗，在滾石防治工程中應(yīng)考慮這部分損耗，而滾石轉(zhuǎn)動能量只占總動能的1/30～1/10，因此一般不單獨考慮滾石的轉(zhuǎn)動過程。研究成果可為滾石防治工程的設(shè)計及優(yōu)化提供參考。

關(guān) 鍵 詞： 滾石運動特征; 物理試驗; 離散元模擬; 合壁津高速公路

中圖法分類號： ?P642.2

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

DOI： 10.16232/j.cnki.1001-4179.2022.03.030

0 引 言

滾石運動過程是滾石災(zāi)害防治工程中的關(guān)注重點，其基本形式可劃分為滾動、斜拋和碰撞? [1-2] 。國內(nèi)外學(xué)者對斜坡滾石運動特征規(guī)律開展了廣泛的研究：黃潤秋等? [3] 通過開展?jié)L石現(xiàn)場原位試驗研究了坡角對滾石運動特征的影響;Niklaus等? [4] 通過將自研的“StoneNode”儀器植入滾石內(nèi)部獲取了滾石運動特征參數(shù)，并對其運動規(guī)律進行了分析;Giacomini等? [5] 將“立體測量影像系統(tǒng)”應(yīng)用于滾石運動軌跡分析中，通過重建的滾石三維運動軌跡實時判斷滾石的運動情況;張龍等? [6] 開展室外滾石運動試驗得出試塊形狀是影響滾石運動的最顯著因素;Spadari? [7] 通過原位試驗得到滾石恢復(fù)系數(shù)的經(jīng)驗值;楊仲元? [8] 采用數(shù)值分析手段對滾石運動恢復(fù)系數(shù)進行了討論，并提出基于衰減反彈接觸法（DROM）的滾石路徑分析方法;Liang，Clement等? [9-11] 收集了滾石運動過程中的震動信號，并通過震動信號反演得到滾石運動特征規(guī)律。此外，還有部分學(xué)者采用數(shù)值模擬手段研究了影響滾石運動過程的因素? [12-13] 。

本文依托在建合壁津高速公路危巖落石防治工程項目，采用“室內(nèi)試驗+數(shù)值分析”的手段，對滾石運動特征進行了系統(tǒng)性研究，研究成果可為該區(qū)段滾石防治工程的設(shè)計提供設(shè)計參考。

1 研究區(qū)概況

研究區(qū)域地處四川盆地東側(cè)，大巴山前側(cè)，地貌類型為山地地貌，由構(gòu)造剝蝕作用形成。區(qū)域內(nèi)地勢東南高西北低，存在較大落差，最高點在南側(cè)山脊高程約 820 m，最低點在北部溝谷地帶高程約480 m，相對高差約340 m（見圖1）。

研究區(qū)危巖落石體量和方量較大，體積從幾方至上千方不等，且多散布于邊坡上，因坡表植被和坡度影響產(chǎn)生堆積。區(qū)域內(nèi)危巖落石主要為中風(fēng)化砂巖，且多為帶菱角的不規(guī)則形狀，塊徑為2～50 cm不等，并伴有少量塊徑為80～250 cm的大塊石出現(xiàn)。落石體積集中于0.125～4.0 m? 3 之間，最大尺寸可達(dá)15 m? 3 。

2 物理模型試驗

2.1 試驗設(shè)計

物理模型試驗結(jié)合研究區(qū)背景進行相似比的考慮。通過運用相似第二定理，將應(yīng)力 σ、應(yīng)變ε、泊松比μ、彈性模量E、位移δ、接觸力F作為量綱分析參量，并選取E、δ 為基本量。

f σ，ε，μ，E，δ，v，J =0 （1）

各參數(shù)量綱可為

σ =ML? -1 T? -1?? E =ML? -1 T? -1?? δ =L? F =MLT? -1? ?（2）

由此可得到量綱關(guān)系式為

π 1= σ E? y 1 δ? y 2? = ML? -1 T? -1??? ML? -1 T? -1???? y 1 L? y 2?? π 2= ε E? y 3 δ? y 4? = ε ??ML? -1 T? -1???? y 3 L? y 4?? π 3= μ E? y 5 δ? y 6? = μ ??ML? -1 T? -1???? y 5 L? y 6?? π 4= F E? y 7 δ? y 8? = MLT? -1??? ML? -1 T? -1???? y 7 L? y 8? ??（3）

計算各式π值為

π 1= σ E ，π 2=ε，π 3=μ，π 4= F E 1δ 2 ?，

因而試驗所采用的相似準(zhǔn)則為

C ε=C μ=1 （4）

式中： C ε為應(yīng)變相似比，C μ 為泊松比相似比。

試驗以幾何相似比 C L =100為基礎(chǔ)，對其他相似比參數(shù)進行合理控制，使物理模型試驗?zāi)芘c現(xiàn)場情況較好擬合。根據(jù)研究區(qū)斜坡坡度范圍構(gòu)建物理模型試驗，如圖2所示，坡面形態(tài)與研究區(qū)基本保持一致，如圖3所示。

根據(jù)圖3可知，試驗斜坡模型高1.0 m，長0.9 m，坡度變化范圍介于26°～70°之間，分別選用混凝土和黏土作為坡表覆蓋層，用來模擬邊坡表面為裸露巖體和土體的2種工況。試塊形狀考慮對運動最有利的方形試塊，試塊尺寸為3 cm×3 cm×3 cm，在物理性質(zhì)上，試塊制作采用C30混凝土且泊松比為0.3，這與現(xiàn)場調(diào)查得到的砂巖物理性質(zhì)存在相似性。

采用如圖4所示的方法釋放試塊，將試塊置于坡頂，向臨空面緩緩?fù)频剑科渥陨碇亓ο侣?。運動全過程由攝像機側(cè)面拍攝，記錄停止位置、運動距離及偏移量。為消除試驗誤差，每組試驗重復(fù)進行60次，并選取90%的可用組數(shù)開展后續(xù)結(jié)果計算。

2.2 試驗理論計算

2.2.1 試塊運動速度

假設(shè)滾石發(fā)生運動后，于 t ?a 時刻運動到A點（x? a ，y? a ），于t? b 時刻運動到B點（x? b ，y? b ） （見圖5），由此可計算滾石豎向和橫向偏移速度：

v x= x b-x a t b-t a ，v y= y b-y a t b-t a? （5）

2.2.2 試塊恢復(fù)系數(shù)

根據(jù)塊體非彈性碰撞相關(guān)理論，滾石在運動過程中常與坡面產(chǎn)生碰撞耗能作用，導(dǎo)致其運動速度下降，速度減幅即為能量損耗的大小，計算公式為

R= u r u i? （6）

式中： R為滾石恢復(fù)系數(shù)，u? i 和u? r? 分別為滾石碰撞坡面前后的運動平均速度（見圖6）。

由于滾石在斜坡上的運動可分解為切向和法向，因此分別采用法向恢復(fù)系數(shù)和切向恢復(fù)系數(shù)表征滾石的速度變化情況：

R n= u? r，n? u? i，n? ?（7）

R t= u? r，t? u? i，t? ?（8）

2.3 試驗現(xiàn)象描述

根據(jù)對試驗現(xiàn)象的觀察，試塊下落過程可劃分為快速下落階段和緩速下落階段（見圖7）。

（1） 快速下落階段。該段為運動的前半程，斜坡坡角較大，試塊以滾動和跳躍為主要姿態(tài)，試塊碰撞耗能較低，速度在較短時間內(nèi)快速增長。

（2） 緩速下落階段。由于斜坡坡角逐漸減小，該段試塊運動模式以滾動為主，試塊與坡面碰撞次數(shù)增多，試塊運動增速明顯減緩。

2.4 試驗結(jié)果分析

2.4.1 土體覆蓋層對運動速度的影響

圖8為土體有無覆蓋層情況下對所有試驗組求均值后得到的速度時程曲線。試塊從坡頂開始下落后速度迅速增大，并且在這個過程中，試塊不斷旋轉(zhuǎn)、與坡體產(chǎn)生碰撞，曲線上每一個速度的突變點都表示碰撞點。

有土體覆蓋層工況下，試塊于0.7 s運動至坡底，運動速度峰值為2.7 m/s;無土體覆蓋層工況下，試塊于0.8 s運動至坡底，運動峰值速度為2.1 m/s。試塊在到達(dá)坡底后，由于碰撞耗能作用，以及試塊在平地發(fā)生的滾動摩擦耗能，速度進一步減小趨于零。

對曲線進一步分析可得到，土體覆蓋層的有無對滾石運動速度存在較大影響，以峰值運動速度為例，兩種工況差值接近于20%，試塊與坡面發(fā)生碰撞后速度的變化也十分明顯。

2.4.2? 土體覆蓋層對運動距離與偏移量的影響

圖9為滾石試塊運動距離-偏移量示意圖。根據(jù)試驗數(shù)據(jù)繪制有無覆蓋層條件下橫向偏移量和水平運動距離的散點圖如圖10所示。由圖10可知：試塊偏移量與運動距離存在強相關(guān)關(guān)系，且在多數(shù)場景下，橫向偏移量與水平運動距離成正相關(guān)關(guān)系，且散點大多位于 d =0.42l-46.3和 d =0.23l-19.7兩條曲線之間的區(qū)域。

兩種工況下，試塊的水平運動距離多集中在160～200 cm間，分別達(dá)到有土體覆蓋層工況和無土體覆蓋層工況試驗組數(shù)的75%和83%以上。但2種工況下，試塊不同水平運動距離區(qū)間存在顯著差異，在有土體覆蓋層工況下，試塊主要停留在160～180 cm之間，達(dá)到25次，占比約50%，而在無覆蓋層工況下，有17次試驗試塊停積在這個距離內(nèi)，占比約30%;在180～200 cm區(qū)間，多停積無覆蓋層工況試塊，記數(shù)為28次，超過總數(shù)的50%，無覆蓋層情況下在該區(qū)間停積試塊占比為30%，與水平運動距離相比，試塊橫向偏移也呈現(xiàn)出相似規(guī)律。據(jù)此可知，在無土體覆蓋層情況下，試塊表現(xiàn)出更強的運動特性。

此外，采用偏移比表征滾石橫向運動的能力：

α= d l? （9）

式中： l 為滾石沿最大坡降方向的運動距離， d 為滾石橫向偏移量。

圖11為試塊水平運動距離與偏移比的關(guān)系散點圖。從圖11可以看出：試塊偏移比分布于0.08～ 0.18 的區(qū)間內(nèi)，超過90%的數(shù)據(jù)分布于偏移比為0.12～0.16的區(qū)間內(nèi)。此外，無土體覆蓋層的試驗組具有較高偏移比（>0.14）的組數(shù)更多，約為69%，說明無土體覆蓋層情況下滾石具有較強的運動能力。

2.4.3 有無覆蓋層對運動恢復(fù)系數(shù)的影響

通過對每組試驗切向與法向恢復(fù)系數(shù)的典型值進行計算，繪制有無土體覆蓋層情況下恢復(fù)系數(shù)的關(guān)系如圖12所示。由圖12可知：不同覆蓋層條件下2項恢復(fù)系數(shù)有較大區(qū)別，對于有土體覆蓋層工況，法向恢復(fù)系數(shù)位于區(qū)間0.25～0.45，切向恢復(fù)系數(shù)位于區(qū)間0.36～0.75;對于無土體覆蓋層工況，法向恢復(fù)系數(shù)位于區(qū)間0.35～0.58，切向恢復(fù)系數(shù)位于區(qū)間0.63～ 0.96 。后一種工況下的恢復(fù)系數(shù)明顯大于前者，這與Asteriou文中巖質(zhì)邊坡恢復(fù)系數(shù)明顯大于土質(zhì)邊坡的結(jié)論類似。此外，無論哪種工況，恢復(fù)系數(shù)值分布較散亂，這可能與滾石與坡面碰撞時的姿態(tài)有關(guān)，有的是碰撞點，有的是碰撞面，另外由于試驗組數(shù)有限，統(tǒng)計數(shù)據(jù)較少，也會出現(xiàn)這種情況。

3 滾石運動過程數(shù)值模擬研究

為彌補物理實驗組數(shù)較少造成的結(jié)果局限性，使結(jié)果更為真實可靠，因此采用基于離散元方法的PFC軟件選取研究區(qū)內(nèi)邊坡研究滾石的運動過程。

3.1 計算模型校正

以室內(nèi)試驗為藍(lán)本，1 ∶1構(gòu)建離散元模型，模型由運動坡面和停積平臺兩大部分組成，坡體尺寸如圖13所示。參照室內(nèi)試驗將坡體分為7段，停積平臺根據(jù)需要取1.0～1.5 m。

模擬過程中，不考慮滾石的變形和破裂，不考慮坡體表面的形變，采用線性接觸模型，離散元參數(shù)取值參考蔣景彩等? [14] 所提出的滾石離散元數(shù)值模擬參數(shù)標(biāo)定方法進行標(biāo)定，計算取值如表1所列。

圖14為土質(zhì)覆蓋層和混凝土覆蓋層2種工況下的物理模型試驗與數(shù)值計算速度趨勢對比分析圖。由圖14可知：數(shù)值計算結(jié)果一定程度上還原了物理模型試驗，絕大部分?jǐn)?shù)據(jù)點的速度誤差值控制在13%以內(nèi)，因此可認(rèn)為該離散元的模擬分析方法和所選取的計算參數(shù)適宜，可用于滾石運動過程的分析。

3.2 模擬計算工況

選取擬定研究區(qū)邊坡為研究對象，在此基礎(chǔ)上開展數(shù)值計算工作，以擬定研究區(qū)1 ∶500地形圖為參照構(gòu)建1 ∶1邊坡三維計算模型。對所選邊坡開展野外調(diào)查，將滾石大致劃分為“條狀”和“塊狀”，并結(jié)合調(diào)查結(jié)果加入較不利的“似球狀”滾石，每種形狀的試塊分別采用50，100，150 cm 3種尺寸。采用clump軟件生成滿足要求的試塊模型（見圖15）。所選定的斜坡坡度較緩、坡表覆蓋層主要為強風(fēng)化砂巖、松軟土體，可參照表1中土質(zhì)坡面的離散元參數(shù)取值。

3.3 計算結(jié)果分析

（1） 運動速度。圖16為條狀、塊狀、似球狀滾石的運動速度時程曲線。從圖16可以看出：各塊徑下滾石試塊釋放后其速度迅速的增大，后迅速減小至零。小尺寸滾石峰值運動速度往往最大。條狀、塊狀和似球狀小尺寸滾石運動速度分別為13.7，13.8 m/s和17.8 m/s，最大塊徑的峰值速度分別為8.2，8.6 m/s和11.3 m/s，這可能是受慣性影響，較小質(zhì)量滾石易改變運動狀態(tài)，更容易獲得較大速度。此外，形狀對速度也有影響，越是接近球形，滾石運動速度越大;同尺寸的滾石，似球狀峰值速度分別比條狀和塊狀滾石速度大約20%～37%。

（2） 轉(zhuǎn)動速度。滾石在下落過程中轉(zhuǎn)動和平動交替或同時存在。圖17為不同形狀下滾石轉(zhuǎn)動速度模擬曲線圖。對于條狀、塊狀、似球狀塊徑50 cm的滾石，其峰值轉(zhuǎn)動速度分別為50，55 rad/s和58 rad/s;而對于塊徑250 cm的不同形狀滾石，其峰值轉(zhuǎn)動速度僅為5.2，6.3 rad/s和7.1 rad/s，具有顯著的差異。這也表明對于下落的滾石，其質(zhì)量較大時速度的變化多由平動速度影響，而質(zhì)量較小時轉(zhuǎn)動速度對合速度的影響很大。似球狀滾石轉(zhuǎn)動速度-時間曲線不同于其他兩種工況，其曲線較平滑，無明顯或很少存在“階梯”，表明這種形狀更易發(fā)生運動狀態(tài)，特別是轉(zhuǎn)動速度的變化，穩(wěn)定性較差。

（3） ?能量變化。滾石下落過程中，隨著運動速度的變化，能量也隨之改變。圖18為不同形狀、不同塊徑滾石運動過程中機械能的變化。運動初始階段，滾石的全部機械能為重力勢能，隨著沿斜坡向下運動，勢能轉(zhuǎn)化為動能，后者根據(jù)其運動形式又可分為平動動能和轉(zhuǎn)動動能。轉(zhuǎn)動動能在總動能中只占極小一部分約1/30～1/10。因此在滾石防治設(shè)計中往往忽略這部分能量。

（4） 運動軌跡。圖19表示了不同形狀滾石運動軌跡圖，其中 條狀、塊狀滾石受其形狀因素影響明顯，

呈現(xiàn)出不同的運動軌跡，運動距離也存在較大的偏差。而不同塊徑似球狀滾石的斜坡運動軌跡相似，但到達(dá)坡底時，大尺寸滾石因接觸面積大會受到更大的阻力，而小尺寸滾石受阻力較小，因而運動距離較遠(yuǎn)。表2統(tǒng)計了不同形狀滾石運動過程參數(shù)，在僅考慮運動速度、轉(zhuǎn)動速度、運動距離情況下，小尺寸滾石運動能力、活動性更強。

4 結(jié) 論

為探究坡表覆蓋層及滾石自身特性對滾石運動特征的影響，本文采用室內(nèi)試驗及離散元模擬的方法，系統(tǒng)研究了坡表不同覆蓋層和滾石塊徑、形狀對滾石運動特性的影響，結(jié)論如下。

（1） 坡表覆蓋層性質(zhì)對滾石運動速度有顯著影響，混凝土坡面較土質(zhì)坡面滾石運動速度提高約20%。下落過程中，斜坡坡度對速度的變化影響較大，較陡時（試驗條件下大于40°）滾石呈現(xiàn)出顯著加速狀態(tài)，小于這個角度則加速度明顯減小。

（2） 滾石運動距離與其運動速度有較好的相關(guān)性，大多數(shù)情況下，混凝土坡面滾石的運動距離更大，平均比土質(zhì)坡面大約20 cm，同時其橫向偏移量也更大約是另一種工況下的1.2倍。試驗中，約90%滾石的運動偏移比在0.08～0.16區(qū)間，同時近70%的混凝土坡面運動試塊的偏移比大于土質(zhì)坡面，表明在混凝土坡面上的滾石具有較強的運動能力。

（3） 對邊坡滾石運動的數(shù)值模擬表明，滾石運動速度受形狀和塊徑影響，相同塊徑下滾石越接近球形，其運動速度越快，模擬中似球狀滾石可比同塊徑下的條狀滾石運動速度大約30%;相同形狀下的滾石，較小質(zhì)量的可比大質(zhì)量速度大近60%。

（4） 滾石下落運動過程中，運動狀態(tài)包括平動和轉(zhuǎn)動，和速度有類似的趨勢，小質(zhì)量滾石由于慣性小，易改變運動狀態(tài)，其在運動過程中轉(zhuǎn)動速度往往較大，相反，大質(zhì)量滾石一般具有較小的轉(zhuǎn)動速度，但這種變化與滾石形狀關(guān)系不大。

（5） 滾石下落過程中涉及復(fù)雜的能量轉(zhuǎn)化，重力勢能轉(zhuǎn)化為動能和與坡面碰撞損耗的內(nèi)能，動能則包括平動動能和轉(zhuǎn)動動能，后者一般是前者的1/10～ 1/30 ，因此在滾石災(zāi)害防治過程中，通常不考慮轉(zhuǎn)動能量的影響。

（6） 綜上所述，在實際工程中，可根據(jù)本文研究得到的不同覆蓋層性質(zhì)、坡面類型、滾石形狀和塊徑對滾石的運動距離、運動速度、偏移比的影響，為滾石防治結(jié)構(gòu)的安置位置、修筑高度和能級等提供合理化建議。

參考文獻(xiàn)：

[1]? 胡厚田.崩塌與落石[M].北京：中國鐵道出版社，1989：1-19.

[2] 趙麗娜，周科平，高峰，等.露天礦邊坡滾石運動特征及控制[J].災(zāi)害學(xué)，2008（3）：76-79.

[3] 黃潤秋，劉衛(wèi)華，周江平，等.滾石運動特征試驗研究[J].巖土工程學(xué)報，2007（9）：1296-1302.

[4] NIKLAUS P，BIRCHLERT，AEBI T，ET A L.StoneNode：a low-power sensor device for induced rockfall experiments[C]∥2017 IEEE Sensors Applications Symposium（SAS）.IEEE，2017.

[5] GIACOMINI A，THOENI K，LAMBERT C，et al.Experimental study on rockfall drapery systems for open pit highwalls[J].International Journal of Rock Mechanics & Mining ences，2017（56）：171-181.

[6] 張龍，雍睿，倪衛(wèi)達(dá).滾石運動特征影響因素分析[J].安全與環(huán)境工程，2011，18（3）：17-21.

[7] SPADARI M，GIACOMINI A，BUZZI O，et al.In situ rockfall testing in New South Wales，Australia[J].International Journal of Rock Mechanics & Mining Sciences，2012，49（2）：84-93.

[8] 楊仲元.道路邊坡危巖落石運動路徑研究[J].公路交通科技，2010，27（1）：34-38.

[9] LIANG F，VERONICA P，EMANUELE I，et al.Rockfall seismic features analysis based on in situ tests：frequency，amplitude，and duration[J].Journal of Mountain Science，2019，16（5）：955-970.

[10]? CLMENT H，ANNE MA，GILLES G，et al.Spatio-temporal evolution of rockfall activity from 2007 to 2011 at the Piton de la Fournaise volcano inferred from seismic data[J].Journal of Volcanology and Geothermal Research，2017，333-334：36-52.

[11] CLMENT H，JEAN-PHILIPPE M，F(xiàn)RANCK Br，et al.Single-block rockfall dynamics inferred from seismic signal analysis，2017（5）：283-292.

[12] LEINE R I，SCHWEIZER A，CHRISTEN M，et al.Simulation of rockfall trajectories with consideration[J].Multibody Syst Dyn，2014，32：241-271.

[13] HU J，LI S，LI L P，et al.Field，experimental，and numerical investigation of a rockfall above a tunnel portal in southwestern China[J].Bull Eng Geol Environ，2017，77：1365-1382.

[14] 蔣景彩，能野一美，山上拓男.滾石離散元數(shù)值模擬的參數(shù)反演（英文）[J].巖石力學(xué)與巖土工程學(xué)報，2008，27（12）：2418-2430.

（編輯：劉 媛）

Physical modelling experiment on rockfall motion characteristics

LIU Yanhui

（ China Railway First Survey and Design Institute Group Co.，Ltd.，Xian 710043，China ）

Abstract：

The research on the motion characteristics of rockfall is an essential vital premise to preventing and controlling the disaster of rockfall.Furthermore，both the covering layer of the slope surface and the rockfall characteristics are essential factors that affect the motion characteristics.Taking the material of the slope covering layer as a variable，a physical test was conducted to study the motion characteristics of rockfall on the concrete slope and soil slope.On this basis，taking the slope of the Hechuan-Bishan-Jiangjin Expressway as an example，the numerical simulation study of the motion characteristics of rockfall was carried out.The influence of the shape and size of the rockfall on the movement process was analyzed by considering the velocity，rotation velocity，and kinetic energy as indicators.The results showed that the shape and size significantly influence the motion characteristics of the rockfall.The closer to a ball-likely shape and the smaller size，the greater the movement speed and rotation speed.More than 30% of the potential energy of the rockfall was lost during the collision，which should be considered in the rockfall stone control project.The rotational kinetic energy of the rockfall only accounted for 130～110 of the total kinetic energy，so the rotating process of rockfall was generally not considered separately.The research results can provide reference for the design and optimization of the rockfall prevention project.

Key words：

rockfall motion characteristics;physical modelling experiment;discrete element simulation;Hechuan-Bishan-Jiangjin Expressway