滾石運(yùn)動(dòng)特性的物理模型試驗(yàn)研究

劉彥輝

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

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

中圖法分類號(hào)： ?P642.2

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

DOI： 10.16232/j.cnki.1001-4179.2022.03.030

0 引 言

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

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

1 研究區(qū)概況

研究區(qū)域地處四川盆地東側(cè)，大巴山前側(cè)，地貌類型為山地地貌，由構(gòu)造剝蝕作用形成。區(qū)域內(nèi)地勢(shì)東南高西北低，存在較大落差，最高點(diǎn)在南側(cè)山脊高程約 820 m，最低點(diǎn)在北部溝谷地帶高程約480 m，相對(duì)高差約340 m（見(jiàn)圖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 物理模型試驗(yàn)

2.1 試驗(yàn)設(shè)計(jì)

物理模型試驗(yàn)結(jié)合研究區(qū)背景進(jìn)行相似比的考慮。通過(guò)運(yùn)用相似第二定理，將應(yīng)力 σ、應(yīng)變?chǔ)?、泊松比μ、彈性模量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）

計(jì)算各式π值為

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

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

C ε=C μ=1 （4）

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

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

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

采用如圖4所示的方法釋放試塊，將試塊置于坡頂，向臨空面緩緩?fù)频?，靠其自身重力下落。運(yùn)動(dòng)全過(guò)程由攝像機(jī)側(cè)面拍攝，記錄停止位置、運(yùn)動(dòng)距離及偏移量。為消除試驗(yàn)誤差，每組試驗(yàn)重復(fù)進(jìn)行60次，并選取90%的可用組數(shù)開(kāi)展后續(xù)結(jié)果計(jì)算。

2.2 試驗(yàn)理論計(jì)算

2.2.1 試塊運(yùn)動(dòng)速度

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

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)理論，滾石在運(yùn)動(dòng)過(guò)程中常與坡面產(chǎn)生碰撞耗能作用，導(dǎo)致其運(yùn)動(dòng)速度下降，速度減幅即為能量損耗的大小，計(jì)算公式為

R= u r u i? （6）

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

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

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

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

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

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

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

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

2.4 試驗(yàn)結(jié)果分析

2.4.1 土體覆蓋層對(duì)運(yùn)動(dòng)速度的影響

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

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

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

2.4.2? 土體覆蓋層對(duì)運(yùn)動(dòng)距離與偏移量的影響

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

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

此外，采用偏移比表征滾石橫向運(yùn)動(dòng)的能力：

α= d l? （9）

式中： l 為滾石沿最大坡降方向的運(yùn)動(dòng)距離， d 為滾石橫向偏移量。

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

2.4.3 有無(wú)覆蓋層對(duì)運(yùn)動(dòng)恢復(fù)系數(shù)的影響

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

3 滾石運(yùn)動(dòng)過(guò)程數(shù)值模擬研究

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

3.1 計(jì)算模型校正

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

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

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

3.2 模擬計(jì)算工況

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

3.3 計(jì)算結(jié)果分析

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

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

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

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

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

4 結(jié) 論

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

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

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

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

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

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

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

參考文獻(xiàn)：

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

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

[3] 黃潤(rùn)秋，劉衛(wèi)華，周江平，等.滾石運(yùn)動(dòng)特征試驗(yàn)研究[J].巖土工程學(xué)報(bào)，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á).滾石運(yùn)動(dòng)特征影響因素分析[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] 楊仲元.道路邊坡危巖落石運(yùn)動(dòng)路徑研究[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é)報(bào)，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