牛林峰,朱 瑞,王忠福,陳鈞龍

(1.華北水利水電大學河南省巖土力學與結(jié)構(gòu)工程重點實驗室,河南鄭州450045；2.河南省水利勘測設計研究有限公司,河南鄭州450016)

基于PIV技術(shù)的反傾層狀結(jié)構(gòu)邊坡模型試驗研究

牛林峰1,朱 瑞2,王忠福1,陳鈞龍1

(1.華北水利水電大學河南省巖土力學與結(jié)構(gòu)工程重點實驗室,河南鄭州450045；2.河南省水利勘測設計研究有限公司,河南鄭州450016)

基于PIV技術(shù),對頂部加載工況下的反傾層狀結(jié)構(gòu)邊坡的變形破壞過程進行記錄監(jiān)測,研究反傾層狀結(jié)構(gòu)邊坡的變形破壞機制。結(jié)果表明,頂部荷載較小時,邊坡模型變形較小,位移集中在模型上部和后緣,形成與巖層平行的應變集中帶；隨著荷載的加大,邊坡模型變形增大,上部位移方向反轉(zhuǎn),出現(xiàn)橫向貫通的應變集中帶,并隨時間向下推移,最終位置固定,從坡腳處向上延伸,貫通整個邊坡模型,與水平成15°～17°的夾角；邊坡模型完全破壞時,形成自上而下的3條近于平行的破裂帶,最底部的最先形成,最頂部的最后形成,變形破壞程度自上而下依次增大。

PIV技術(shù)；反傾層狀邊坡；頂部加載；變形；破裂帶

0 引 言

反傾層狀結(jié)構(gòu)邊坡的變形破壞問題廣泛存在于水利水電工程、道路交通工程、采礦工程等領域[1-2],已經(jīng)成為工程建設中亟待解決的重要工程地質(zhì)問題。對反傾層狀結(jié)構(gòu)邊坡的研究始于20世紀60年代,研究方法可分為解析分析法、數(shù)值模擬法、物理模擬法。其中,物理模擬法可以有效解決解析分析法和數(shù)值模擬法難以解決的問題,能直接記錄和觀測反傾層狀邊坡的變形破壞過程。汪小剛等[3]結(jié)合龍灘水電站左岸邊坡實例,利用離心機試驗研究了反傾巖質(zhì)邊坡的變形破壞；左保成等[4]通過室內(nèi)物理力學試驗研究了反傾層狀邊坡破壞模式和影響因素；楊國香等[5]通過大型振動臺試驗探討了反傾層狀邊坡的動力響應特征和破壞模式。

隨著計算機技術(shù)的發(fā)展,高精度圖像處理技術(shù)在試驗中得到廣泛應用。粒子圖像測速(PIV)技術(shù)被廣泛應用于工程領域的試驗研究,其測量精度高、干擾小、穩(wěn)定性強、工作量小、自動化程度高,可完好地記錄了試驗過程,加深了對研究對象的變形破壞機制的研究。利用PIV技術(shù),劉明亮等[6]研究了錨板抗拉破壞模式；姜彤等[7]對抗拔錨板的裙錨效應進行了研究；劉君等[8]將PIV技術(shù)應用于大型振動臺試驗中,并取得較好的效果；沈禮等[9]對構(gòu)造物理模型試驗進行了監(jiān)測分析。

本文采用PIV技術(shù),對頂部加載工況下的反傾層狀結(jié)構(gòu)邊坡的變形破壞過程進行分析研究,通過對采集的灰度圖像的處理,再現(xiàn)了邊坡變形破壞過程的位移場變化,進一步揭示該類邊坡的變形破壞規(guī)律,為反傾層狀結(jié)構(gòu)邊坡的治理提供參考。

1 PIV技術(shù)原理

PIV技術(shù)克服了傳統(tǒng)單點測量的限制,實現(xiàn)了動態(tài)瞬時全流場的實時測量,并且具有不干擾測量對象,測量精度高的優(yōu)點,基本原理見圖1。利用圖片采集設備采集t0時刻(變形前)和t1時刻(變形后)測量對象的圖片,將采集的t0、t1時刻的灰度圖像分割成若干均勻網(wǎng)格,利用交叉關(guān)聯(lián)函數(shù)對變形前(即t0時刻)圖像的某一網(wǎng)格在變形后(即t1時刻)圖像的指定范圍內(nèi)進行全場匹配和相關(guān)運算。交叉關(guān)聯(lián)函數(shù)R(Δx,Δy)為[10]

圖1 PIV基本原理

(1)

式中,a、b為圖像塊的尺寸；Y為t0時刻圖像中圖像塊中心坐標在(m,n)處的灰度值分布函數(shù)；X為t1時刻圖像中圖像塊中心坐標在(m+Δx,n+Δy)處的灰度值分布函數(shù)；Δx、Δy分別為x、y方向的位移增量。

當關(guān)聯(lián)函數(shù)值達到峰值時,根據(jù)峰值相關(guān)系數(shù)確定該網(wǎng)格變形后的位置,進而得到該網(wǎng)格的像素位移,將像素位移按一定比例關(guān)系轉(zhuǎn)換可得到物理位移,對變形前后的所有網(wǎng)格進行類似運算就可以得到整個圖像的位移場。對得到的位移場進行計算[11-12],可得出其對應的應變場,據(jù)此可進一步分析研究對象的變形破壞機理。公式如下

(2)

(3)

(4)

式中,u為x方向的位移分量；v為y方向的位移分量；εx為沿x方向的正應變；εy為沿y方向的正應變；γxy為剪應變。

2 試驗模型及設備

2.1 試驗模型

反傾層狀結(jié)構(gòu)邊坡模型由提前預制好的薄板堆砌而成,預制板長80 cm、寬15 cm、厚1.5 cm,由石膏和經(jīng)過200目篩子篩選的黃土按石膏∶黃土∶水為1∶1∶2的比例配置而成。經(jīng)室內(nèi)試驗,其物理力學特性滿足試驗要求。由于預制板的厚度小,長寬比大,在脫模的過程要格外注意,脫模時間要適當延遲。堆砌邊坡模型的高70 cm、寬15 cm、底部長80 cm、坡頂長30 cm、坡腳距底部邊界15 cm、距模型前緣邊界18 cm、坡角60°、層面傾角75°,固定在提前定制好的模型框中,制作好的邊坡模型見圖2。

圖2 邊坡模型

2.2 試驗設備

將制作好的物理模型放置在加載裝置中,通過與加載裝置連接的壓力板對模型頂部加壓,使其變形破壞。同時,利用PIV系統(tǒng)的CCD相機獲取模型變形破壞過程的位移圖像。應用PIVview2C圖像處理軟件,分析邊坡模型破壞過程中的位移場和應變場。

加載裝置為CMT4000系列試驗機,可提供拉、壓向力,最大力30 kN,最大位移1 150 mm。通過PC機上安裝的SANS軟件,可方便控制加載力的速度和方向,并及時輸出試驗過程中力與位移的變化,力精度和位移精度分別為1/300 000 FS和0.03 μm,滿足試驗要求。CCD高速相機采用SONY科研級芯片,分辨率為2 489 pixel×2 091 pixel,曝光時間為100 μm×80 ms,采集速率可達200 fps,通過Camlink專用接口連接計算機,利用Davis8.0系列軟件儲存采集圖片。試驗設備組合系統(tǒng)見圖3。

圖3 試驗設備組合系統(tǒng)

3 試驗結(jié)果及分析

3.1 頂部加載力與時間的關(guān)系

圖4為整個試驗過程邊坡頂部受力與時間的關(guān)系。根據(jù)其變化特征可將試驗過程分為4個階段：0～1 000 s為第1階段；1 000～4 000 s為第2階段；4 000～5 000 s為第3階段； 5 000 s以后為第4階段。第1階段力與時間關(guān)系呈現(xiàn)較規(guī)則的線性關(guān)系,斜率為4.0,隨著時間的增大,模型頂部受力線性增大,最大為4 300 N；第2階段的受力在2 000～2 200 s有一次大的跌落,經(jīng)過短暫停滯,迅速回升,回升后的受力特點與跌落前的受力特點基本保持不變,大致線性變化,斜率約為4/3,最大值接近峰值；第3階段受力與時間關(guān)系曲線大致為水平直線,保持在8 500 N左右；第4階段模型受力快速下降,模型完全破壞。

圖4 力與時間的關(guān)系

3.2 變形破壞各階段的位移場及應變場

第1階段邊坡模型變形較小,變形集中在模型中上部和后緣,位移自上而下呈層狀分布,逐漸減小。應變集中帶呈傾斜條帶狀分布在模型后緣,其傾向與巖層傾向一致。圖5為t=1 000 s時模型位移場、應變場分布及模型矢量。

第2階段邊坡模型變形增大,模型下部亦產(chǎn)生較大位移,以某一貫通圓弧界面為分界線,上下位移方向相反,并在圓弧面處產(chǎn)生拉張應變集中帶。圖6為t=2 000 s時模型位移場、應變場分布。隨著時間的增大,圓弧形界面不斷向下推移,并向兩端延伸,最終在t=3 400 s時界面位置固定不再下移,由圓弧形發(fā)展為直線形,從坡腳處向上延伸貫通模型,與水平成15°～17°的夾角,邊坡在該位置形成明顯的張裂帶。圖7為t=3 400 s時模型位移場、應變場分布。圖8為t=4 000 s時模型應變場以及該時刻模型圖像。

圖5 t=1 000 s時模型位移場、應變場分布及其矢量

圖6 t=2 000 s時模型位移場、應變場分布

圖7 t=3 400 s時模型位移場、應變場分布

第3、4階段邊坡模型變形加大,破壞程度加劇,出現(xiàn)新的應變集中帶,其變化規(guī)律與之前的應變集中帶相似,最終形成3條近于平行的應變集中帶,見圖9。邊坡模型在應變集中帶處產(chǎn)生不同程度的拉張斷裂。

圖8 t=4 000 s時模型應變場分布及模型圖像

圖9 t=5 800 s時模型應變場

4 結(jié) 語

運用PIV技術(shù),記錄了反傾層狀邊坡模型在頂部加載工況下的變形破壞過程,通過后期對采集圖像的處理分析,研究邊坡模型的位移場和應變場的變化,得出以下幾點結(jié)論：

(1)頂部荷載較小時,邊坡模型變形較小,位移集中在模型中上部和后緣,沒有橫向貫通的應變集中帶,產(chǎn)生的應變集中帶呈傾斜條帶狀分布在邊坡后緣,其傾向與巖層傾向一致。

(2)隨著頂部荷載的加大,邊坡變形增大,在頂部位移方向反轉(zhuǎn),出現(xiàn)橫向貫通拉張應變集中帶,并隨時間向下推移,最終位置固定,貫通模型從坡腳向上延伸,與水平成15°～17°的夾角,邊坡模型在該處出現(xiàn)明顯拉張破壞。

(3)邊坡模型完全破壞時,形成自上而下的3條近于平行的破裂帶,最底部的破裂帶最先形成,最頂部的最后形成,其變形破壞程度自上而下依次增大。

[1]王思敬. 金川露天礦邊坡變形機制及過程[J]. 巖土工程學報, 1982, 4(1)： 76- 83．

[2]程東幸, 劉大安, 丁恩寶, 等. 層狀反傾巖質(zhì)邊坡影響因素及反傾條件分析[J]. 巖土工程學報, 2005, 27(11)： 1362- 1366．

[3]汪小剛, 張建紅, 趙毓芝, 等. 用離心模型研究巖石邊坡的傾倒破壞[J]. 巖土工程學報, 1996, 18(5)： 14- 21．

[4]左保成, 陳從新, 劉小巍, 等. 反傾巖質(zhì)邊坡破壞機理物理模型研究[J]. 巖石力學與工程學報, 2005, 24(19)： 3505- 3511．

[5]楊國香, 葉林海, 伍法權(quán), 等. 反傾層狀結(jié)構(gòu)巖質(zhì)邊坡動力響應特性及破壞機制振動臺模型試驗研究[J]. 巖石力學與工程學報, 2012, 31(11)： 2214- 2221．

[6]劉明亮, 朱珍德, 劉金元. 錨板抗拉破壞機制試驗研究[J]. 巖土力學, 2011, 32(3)： 697- 702．

[7]姜彤, 任淼, 張昕. 抗拔錨板群錨基礎模型試驗研究[J]. 地下空間與工程學報, 2016, 12(4)： 968- 974．

[8]劉君, 劉福海, 孔憲京, 等. PIV技術(shù)在大型振動臺模型試驗中的應用[J]. 巖土工程學報, 2010, 32(3)： 368- 373．

[9]沈禮, 賈東, 尹宏偉, 等. 構(gòu)造物理模擬和PIV有限應變分析對構(gòu)造裂縫預測的啟示[J]. 高校地質(zhì)學報, 2016, 22(1)： 171- 182．

[10]THOMAS S H, JAGDISHKUMAR K A. Image Se-quence Analysis [M]. New York: Springer-Verlag, 1981．

[11]ROSCOE K H, ARTHUR J R F, JAMES R G. The determination of strains in soils by X-ray method[J]. Civil Engineering & Public Works Review, 1963, 58: 873- 876, 1009- 1012．

[12]YAMAMOTO K, KUSUDA K. Failure mechanisms and bearing capacities of reinforced foundations[J]. Geotextiles and Geomembranes, 2001, 19(3): 127- 162．

(責任編輯楊 健)

ExperimentalStudyofAnti-dipLayeredSlopeModelBasedonPIVTechnology

NIU Linfeng1, ZHU Rui2, WANG Zhongfu1, CHEN Junlong1
(1. Henan Key Laboratory of Geomechanics and Structural Engineering, North China University of Water Resources and Electric Power, Zhengzhou 450045, Henan, China; 2. Henan Water & Power Engineering Consulting Co., Ltd., Zhengzhou 450016, Henan, China)

The process of deformation and failure of anti-dip layered slope under pressure load at slop top are measured based on Particle Image Velocity (PIV) to analyze deformation and failure mechanism. The results show that, (a) when the pressure load is relatively small, there is little deformation, the displacements appear at the upper part and the back of model slope, and the strain concentration is formed in parallel with rock formation; (b) with the increasing of load, the deformation is also increasing, the displacement sense of upper part of slope becomes opposite, and a horizontal run-through strain concentration appears which gradually moves down until it runs through the slope from slope toe to slope back with a angle of 15-17 degree to horizontal line; and (c) when the slope is at complete failure stage, there are three parallel fracture zones in model slope, the one at the bottom of slope appears firstly and the one at the uppermost part of slope appears at the latest, and the failure degree gets increasing from top to bottom of slope．

PIV technology; anti-dip layered slope; tope pressure load; deformation; fracture zone

2017- 02- 21

華北水利水電大學研究生教育創(chuàng)新計劃(YK2016- 04)

牛林峰(1992—),男,河南南陽人,碩士研究生,主要從事邊坡穩(wěn)定性分析研究．

TU457

：A

：0559- 9342(2017)06- 0100- 05