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

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

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

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

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

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

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

0 引 言

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

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

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

1 PIV技術(shù)原理

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

圖1 PIV基本原理

(1)

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

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

(2)

(3)

(4)

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

2 試驗?zāi)Ｐ图霸O(shè)備

2.1 試驗?zāi)Ｐ?/p>

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

圖2 邊坡模型

2.2 試驗設(shè)備

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

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

圖3 試驗設(shè)備組合系統(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 變形破壞各階段的位移場及應(yīng)變場

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

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

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

圖6 t=2 000 s時模型位移場、應(yīng)變場分布

圖7 t=3 400 s時模型位移場、應(yīng)變場分布

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

圖8 t=4 000 s時模型應(yīng)變場分布及模型圖像

圖9 t=5 800 s時模型應(yīng)變場

4 結(jié) 語

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

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

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

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

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

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

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

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

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

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

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

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

[9]沈禮, 賈東, 尹宏偉, 等. 構(gòu)造物理模擬和PIV有限應(yīng)變分析對構(gòu)造裂縫預(yù)測的啟示[J]. 高校地質(zhì)學(xué)報, 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．

(責(zé)任編輯楊 健)

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

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

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

TU457

：A

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