邊坡動力破壞機(jī)理的振動臺試驗研究

劉振平，劉福海，趙顯波

（1.黑龍江工程學(xué)院 土木與建筑工程學(xué)院，哈爾濱 150050；

2.東北電力設(shè)計院，長春130021；3.黑龍江省水利科學(xué)研究院，哈爾濱 150080）

在土木工程和水利工程中，有眾多的自然邊坡、人工邊坡如堤壩、路基路塹等邊坡工程。中國的許多重大工程如大型水電站、南水北調(diào)西線工程、西部山區(qū)高速公路等均處于強(qiáng)震區(qū)，特別是2008年汶川“5·12”地震以來，邊坡動力破壞機(jī)理引起學(xué)者們的廣泛關(guān)注，成為巖土地震工程領(lǐng)域中重要的研究課題之一。有關(guān)邊坡地震的原型觀測資料很少，所以室內(nèi)振動臺試驗就成為研究地震作用下邊坡動力問題的重要手段之一。Lin等［1］采用大型振動臺試驗研究土坡在地震作用下的響應(yīng)，認(rèn)為砂土邊坡破壞面較淺，只是坡體表面的破壞。徐光興等［2－3］通過振動臺模型試驗研究了混合土邊坡的動力特性，并結(jié)合數(shù)值分析結(jié)果認(rèn)為在強(qiáng)震作用下單一均質(zhì)土層邊坡的破壞模式仍然是沿著某一弧形潛在滑動面失穩(wěn)。許強(qiáng)等［4］通過多組振動臺試驗認(rèn)為均質(zhì)土坡和塊狀巖質(zhì)斜坡的動力破壞模式都是坡頂拉裂，中下部剪切滑移破壞。陳新民等［5－6］采用振動臺試驗研究了粘土邊坡動力特性、動力反應(yīng)和宏觀變形，發(fā)現(xiàn)在動力作用下邊坡坡頂和坡腳出現(xiàn)裂縫。葉海林等［7］巖質(zhì)邊坡大型振動臺模型試驗，結(jié)果表明地震滑動面為上部拉裂縫和下部剪切滑移面形成貫通的破裂面。劉君等［8］把PIV技術(shù)應(yīng)用在土質(zhì)邊坡振動臺模型試驗中，獲取了邊坡破壞的完整過程。Wang等［9］在砂土邊坡振動臺試驗中采用PIV技術(shù)獲取了邊坡表面的位移場，探討了邊坡動力失穩(wěn)機(jī)理。Wartman等［10－11］通過多組粘土邊坡振動臺模型試驗研究了邊坡產(chǎn)生的永久位移，并與基于峰值和殘余強(qiáng)度的Newmark公式計算結(jié)果進(jìn)行了比較。Katz等［12］采用振動沙箱研究了砂土邊坡破壞的類型和頻率的關(guān)系。以上針對邊坡的動力特性、宏觀破壞現(xiàn)象及個別點的位移進(jìn)行的研究，而沒有涉及邊坡的位移場，劉君、Wang等雖用PIV技術(shù)得到了整個邊坡的位移場，但沒有對應(yīng)變場進(jìn)行研究。

應(yīng)變是直接對含有噪聲的離散位移進(jìn)行差分計算得到的，那么即使微小的位移測量誤差也會被急劇地放大，使計算的應(yīng)變不可靠［13］。一般的數(shù)據(jù)擬合平滑方法也不能滿足應(yīng)變計算的精度要求，潘兵等［13－14］提出了局部最小二乘擬合的應(yīng)變計算方法，其原理同Savitzky-Golay平滑方法。該方法能夠在一定程度上去除原始位移場的噪聲，缺點是由于每個位移點都要一個計算窗口平滑，計算量巨大，并且對于非均勻應(yīng)變場，應(yīng)變計算窗口的大小顯著影響應(yīng)變的計算精度，目前還沒有合適的方法來選擇最佳的窗口大小。本文采用數(shù)字圖像位移和有限元數(shù)據(jù)平滑技術(shù)，得到了整個邊坡模型的位移場和應(yīng)變場，同時探討了土質(zhì)邊坡的地震破壞機(jī)理。

1 基于有限元數(shù)據(jù)平滑的應(yīng)變計算

數(shù)字圖像位移測量技術(shù)現(xiàn)在已經(jīng)成熟，采用互相關(guān)圖像匹配技術(shù)獲取位移場后，再用有限元數(shù)據(jù)平滑方法對位移場進(jìn)行處理后，得到應(yīng)變場。

1.1 基本原理

基本思路是引入泛函f（F）［15－19］，使該泛函取得極小值的自變函數(shù)F就是要求的光滑位移函數(shù)。

根據(jù)試驗數(shù)據(jù)點的位置分布，把區(qū)域A劃分為若干單元，在每一個單元內(nèi)用形函數(shù)逼近光滑函數(shù)F，這個過程與有限元方法是一樣的。

在第j個單元內(nèi)，把光滑函數(shù)F表示成式中：N j為第j個單元的形函數(shù)；uej為第j個單元的節(jié)點自由度。

式（9）是有限元的標(biāo)準(zhǔn)格式，其中K＝K1＋λK2，λ是平滑系數(shù)，U表示所有單元的節(jié)點自由度。平滑系數(shù)的確定方法在下文專門論述。

只要確定了平滑系數(shù)，就可從式（9）中解出U，再從式（2）中得到平滑位移函數(shù)u，其應(yīng)變的計算方法如下：

有限變形條形下的格林（Green）應(yīng)變［20］（以受壓為正）為

式中：u和v分別為上述所求的x和y方向的位移函數(shù)；εx和εy分別為x和y方向的正應(yīng)變；γxy為工程剪應(yīng)變；γmax為最大剪應(yīng)變。

小變形條件下柯西（Cauchy）應(yīng)變的計算可不考慮上述應(yīng)變表達(dá)式中的二次項。

1.2 形函數(shù)的確定

采用三角形單元，以便適應(yīng)任意形狀的數(shù)據(jù)區(qū)域。同時為提高插值函數(shù)的光滑可導(dǎo)性，采用的是五次多項式的單元插值函數(shù)［15］，其表達(dá)式為式中：X為變量項；C為系數(shù)列向量。由于沿三角形單元三條邊滿足法向?qū)?shù)連續(xù)的條件，21個系數(shù)可縮減為18個。

1.3 平滑系數(shù)的確定

為了能夠自動得到最佳的平滑系數(shù)，Golub等［21］、Craven等［22］及 Bates等［23］提 出 了 廣 義 交 互驗證GCV方法，建立了關(guān)于λ的GCV函數(shù)X是n階滿秩陣，K2是奇異陣，其秩為r。擬總體剛度陣K1是對稱正定陣，因此式（15）成立，可以證明X是滿秩陣，因此y也能夠由式（16）式唯一確定。

求取廣義交互驗證函數(shù)GCV的最小值，就可以得到最佳平滑系數(shù)λ，具體計算方法詳見文獻(xiàn)［24］。

2 試驗設(shè)備及模型設(shè)計

2.1 試驗設(shè)備

試驗在大連理工大學(xué)工程抗震試驗室的水平與豎直雙向水下振動臺上進(jìn)行，其主要性能參數(shù)如下：數(shù)字式控制方式，臺面尺寸3 m×4 m，最大載重100 k N，工作頻率0.1～50 Hz；滿載工況下最大加速度±1.0g（水 平 向）、±0.7g（豎 向）、最 大 速 度±50 cm／s（水平向）、±35 cm／s（豎向）、最大位移±75 mm（水平向）、±50 mm（豎向）。

2.2 試驗?zāi)Ｐ?/h3>

在鋼制模型箱內(nèi)堆制模型進(jìn)行試驗，為便于觀察和圖像采集，一側(cè)為有機(jī)玻璃。模型壩壩高1.3 m，邊坡坡率1∶1.6，其模型如圖1所示。模型填土基本性能參數(shù)見表1。

表1 試驗填土基本性能參數(shù)

圖1 邊坡振動臺模型示意圖（單位：cm）

試驗加載的是頻率10 Hz的正弦增幅波（喇叭波），如圖2所示。從0開始逐漸增加，直至邊坡完全失穩(wěn)。

采用的圖像采集設(shè)備是CANON EOS 450D高清數(shù)碼相機(jī)，最大分辨率4 272×2 848像素，配備10～22 mm廣角鏡頭，可以在距模型2 m處拍攝整個模型的圖片。存儲為JPG格式的圖片文件，最大連拍速度約3.3幀／s。振動過程中連續(xù)的采集圖像，用于數(shù)字圖像位移和應(yīng)變識別。

圖2 輸入的加速度時程

3 試驗結(jié)果分析

3.1 模型的位移應(yīng)變計算流程

開發(fā)了一套數(shù)字圖像測量分析程序，可以批處理多種格式的圖片，可以得到一系列的位移場和應(yīng)變場，還能以等值線圖、云圖和矢量圖等形式輸出各種位移場、應(yīng)變場和速度場。其計算流程和步驟如下：

位移計算流程：圖像像素塊劃分→計算總位移→計算模型箱位移→計算相對位移（總位移減去模型箱剛體位移）→刪除計算的位移壞點→輸出位移場等后處理。

得到位移場后，用有限元平滑位移場再微分計算應(yīng)變場。其計算流程如下：生成單元和節(jié)點信息→計算單元剛度陣和單元節(jié)點荷載列陣→集成總體剛度陣和總體節(jié)點荷載列陣→計算平滑系數(shù)→求解節(jié)點未知參數(shù)→平滑位移場微分求應(yīng)變→輸出應(yīng)變場等后處理。

3.2 邊坡模型的破壞過程描述

邊坡模型破壞過程基本可分3個階段，第1階段是整體變形階段（圖3（a）），當(dāng)輸入加速度較小時，整個模型以均勻沉降變形為主，同時伴有相對于模型箱的水平往復(fù)運動；第2階段是滑移變形階段（圖3（b）），當(dāng)輸入加速度大約在0.1g時，模型開始沿邊坡一側(cè)向下滑移，隨著輸入加速度的增大，滑移量也逐漸增大，在輸入加速度在0.45g時，坡頂出現(xiàn)張拉裂縫，坡體內(nèi)部也出現(xiàn)剪切裂隙，形成明顯的滑動帶，邊坡達(dá)到極限破壞狀態(tài)；第3階段是破壞階段，滑動帶上的裂隙使得滑動體與邊坡相對分離，滑動體加速下滑，最后整個坡體坍塌，邊坡完全破壞。

從圖3（c）可以看出，邊坡模型動力失穩(wěn)存在深層滑動帶，滑動面大概呈圓弧狀。

3.3 位移分析

圖4、圖5和圖6分別是8.42 s（對應(yīng)輸入加速度為0.45g）時即在極限破壞狀態(tài)時的水平位移、豎向位移和總位移的等值線圖，其水平位移出現(xiàn)在坡面中部，豎向位移在坡頂附近，水平位移和總位移等值線基本平行于坡面，均符合一般的邊坡破壞時的位移分布規(guī)律。

圖3 邊坡模型位移矢量

圖4 8.42 s時（輸入加速度0.45 g）時水平位移場（單位：mm）

圖5 8.42 s時（輸入加速度0.45 g）時豎向位移場（單位：mm）

圖6 8.42 s時（輸入加速度0.45 g）時總位移場（單位：mm）

3.4 應(yīng)變分析

從圖7和圖8最大剪應(yīng)變和水平方向應(yīng)變分布上可以看出，坡體中部到坡腳是剪切破壞，坡頂一定深度是拉剪破壞。其變形破壞模式與汶川地震觀察到的邊坡破壞現(xiàn)象相吻合［25］：土坡主要是坡頂向下一定深度內(nèi)的拉破壞，坡腳向上延伸形成剪切滑移帶，最終二者連通形成貫通的破裂面。這也與許強(qiáng)等［4］的振動臺試驗結(jié)果以及鄭穎人等［26－27］用 FLAC程序做的數(shù)值分析結(jié)果基本一致。

圖7 8.42 s時（輸入加速度0.45 g）時水平應(yīng)變場（單位：%）

8 8.42 s時（輸入加速度0.45 g）時最大剪應(yīng)變場（單位：%）

3.5 位移應(yīng)變時程分析

A點位于滑動體上，B點位于滑動體外（圖1），從圖9的位移時程可以看出，位移隨時間是逐漸增加的，沒有明顯的突變點，說明邊坡在地震作用下是漸進(jìn)式破壞的。

圖9 位移時程曲線

由于位移時程曲線沒有明顯的突變拐點，如何判斷模型破壞的時刻就成為一個難題。地震作用下邊坡形成滑移帶后，滑動體與坡體相對分離，由于滑移帶的阻隔，滑動體以外的壩體振動對其影響將大大減小，滑動體將比較平穩(wěn)的運動，而不會隨輸入地震波而振動。所以可以把位移時程曲線的曲率做為判斷壩坡失穩(wěn)滑移的一個物理量，這里所指的曲率是有正負(fù)之分的，正號表示是凹曲線，負(fù)號表示是凸曲線，可定義為廣義曲率。

式中：κ為廣義曲率；y′是位移的一階導(dǎo)數(shù)；y″是位移的二階導(dǎo)數(shù)。

從A點的位移時程曲線曲率（圖10）可以看出，在8.42 s之前，曲率以零線為中心上下波動，說明在此之前沒有發(fā)生連續(xù)的滑動，而在此之后，曲率迅速趨于零，不再發(fā)生波動，該點此時已經(jīng)失穩(wěn)破壞，曲率發(fā)生突變表明滑移帶已經(jīng)貫通形成，邊坡整體失穩(wěn)，達(dá)到了極限抗震狀態(tài)。位于滑動體以外的B點，曲率自始自終都處于振動狀態(tài)，沒有發(fā)生失穩(wěn)破壞。由此可知，用位移時程曲線曲率做為判斷邊坡動力破壞的物理量是可行的。

圖10 位移時程曲線的曲率

4 結(jié) 論

將有限元數(shù)據(jù)平滑方法引入到數(shù)字圖像位移應(yīng)變測量中是可行的，得到了邊坡振動臺模型整個試驗過程的位移場和應(yīng)變場。初步表明土質(zhì)邊坡的變形是漸進(jìn)式的，坡體中部到坡腳是剪切破壞，坡頂一定深度是拉剪破壞，破壞時有深層的圓弧狀滑動面，用位移時程的廣義曲率做為判斷邊坡動力破壞的物理量是可行的。

［1］Lin M L，Wang K L.Seismic slope behavior in a largescale shaking table model test ［J］.Engineering Geology，2006，86（2／3）：118-133.

［2］徐光興，姚令侃，高召寧，等.邊坡動力特性與動力響應(yīng)的大型振動臺模型試驗研究［J］.巖石力學(xué)與工程學(xué)報，2008，27（3）：624-631.

Xu G X，Yao L K，Gao Z N，et al.Large-scale shaking table model test study on dynamic characteristics and dynamic responses of slope［J］.Chinese Journal of Rock Mechanics and Engineering，2008，27（3）：624-631.

［3］徐光興，姚令侃，李朝紅，等.邊坡地震動力響應(yīng)規(guī)律及地震動參數(shù)研究［J］.巖土工程學(xué)報，2008，30（6）：918-923.

Xu G X，Yao L K，Li C G，et al.Dynamic response of slopes under earthquakes and influence of ground motion parameters［J］.Chinese Journal of Geotechnical Engineering，2008，30（6）：918-923.

［4］許強(qiáng)，陳建君，馮文凱，等.斜坡地震響應(yīng)的物理模擬試驗研究［J］.四川大學(xué)學(xué)報：工程科學(xué)版，2009，41（3）：266-272.

Xu J，Chen J J，F(xiàn)eng W K，et al.Study of the seismic response of slopes by physical modeling ［J］.Journal of Sichuan University：Engineering Science Edition，2009，41（3）：266-272.

［5］陳新民，沈建，魏平，等.下蜀土邊坡地震穩(wěn)定性的大型振動臺試驗研究（I）-模型試驗設(shè)計［J］.防災(zāi)減災(zāi)工程學(xué)報，2010，30（5）：492-502.

Chen X M，Shen J，Wei P，et al.Large-scale shaking table test of seismic stability of Xiashu loess slopes（I）：design of model test［J］.Journal of Disaster Prevention and Mitigation Engineering，2010，30（5）：492-502.

［6］陳新民，沈建，魏平，等.下蜀土邊坡地震穩(wěn)定性的大型振動臺試驗研究（II）-試驗結(jié)果及分析［J］.防災(zāi)減災(zāi)工程學(xué)報，2010，30（6）：587-594.

Chen X M，Shen J，Wei P，et al.Large-scale shaking table test of seismic stability of Xiashu loess slopes（II）：analysis of test results ［J］.Journal of Disaster Prevention and Mitigation Engineering，2010，30（6）：587-594.

［7］葉海林，鄭穎人，杜修力，等.邊坡動力破壞特征的振動臺模型試驗與數(shù)值分析［J］.土木工程學(xué)報，2012，45（9）：128-135.

Ye H L，Zheng Y R，Du X L，et al.Shaking table model test and numerical analysis on dynamicfailure characteristics of slope［J］.China Civil Engineering Journal，2012，45（9）：128-135.

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

Liu J，Liu F H，Kong X J，et al.Application of PIV in largescale shaking table model tests ［J］.Chinese Journal of Geotechnical Engineering，2010，32（3）：368-374.

［9］Wang K L，Lin M L.Initiation and displacement of landslide induced by earthquake-a study of shakingtable model slope test［J］.Engineering Geology，2011，122：106-114.

［10］Wartman J.Physical model studies of seismically induced deformations in slopes ［D］.California：University of California-Berkeley，1999.

［11］Wartman J，Seed R B，Bray J D.Shaking table modeling of seismically induced deformations in slopes ［J］.Journal of Geotechnical and Geoenvironmental Engineering，2005，131（5）：610-622.

［12］Katz O，Aharonov E.Andslides in vibrating sand box：what controls types of slope failure and frequencymagnitude relations？［J］.Earth and Planetary Science Letters，2006，247（3）：280-294.

［13］潘兵，謝惠民.數(shù)字圖像相關(guān)中基于位移場局部最小二乘擬合的全場應(yīng)變測量［J］.光學(xué)學(xué)報，2007，27（11）：1980-1986.

Pan B，Xie H M.Full-field strain measurement based on least-square fitting of local displacement for digital image correlation method［J］.Acta Optica Sinica，2007，27（11）：1980-1986.

［14］Pan B，Xie H M，Guo Z Q，et al.Full-field strain measurement using a two-dimensional Savitzky-Golay digital differentiator in digital image correlation ［J］.Optical Engineering，2007，46（3）：1-10.

［15］Segalman D J，Woyak D B，Rowlands R E.Smooth spline-like finite-element differentiation of full-field experimental data over arbitrary geometry ［J］.Experimental Mechanics，1979，19（12）：429-437.

［16］Engelstad M J，Chambless D A，Swinson W F，et al.Hybrid stress analysis of vibrating plates using holographic interferometry and finite elements ［J］.Experimental Mechanics，1987，27（1）：23-30.

［17］Feng Z，Rowlands R E.Continuous full-field representation and differentiation of three-dimensional experimental vector data［J］.Computers and Structures，1987，26（6）：979-990.

［18］Sutton M A，Turner J L，Burck H A，et al.Full-field repersentation of discretely smapled surface deformational for displacement and strain analysis ［J］.Experimental Mechnaics，1991，31（2）：168-177.

［19］Freese C E，Gee L.A multilevel treatment of moiré fringe data using finite elements ［J］.Experimental Mechanics，1999，39（4）：304-310.

［20］陸明萬，羅學(xué)富.彈性理論基礎(chǔ)［M］.北京：清華大學(xué)出版社，2001：38-43.

［21］Golub G H，Heath M，Wahba G.Generalized cross validation as a method for choosing a good ridge parameter［J］.Technometrics，1979，21（2）：215-223.

［22］Craven P，Wahba G.Smoothing noisy data with spline functions：estimating the correct degree of smoothing by the method of generalized cross-validation ［J］.Numerische Mathematik，1979，31：377-403.

［23］Bates D M，Lindstrom M J，Wahba G，et al.Gcvpak-routines for generalized cross validation ［J］.Communications in Statistics-Simulation and Computation，1987，16（1）：263-297.

［24］Kent J T，Mohammadzadeh M.Global optimization of the generalized cross-validation criterion ［J］.Statistics and Computing，2000，10（3）：231-236.

［25］許強(qiáng)，董秀軍.汶川地震大型滑坡成因模式［J］.中國地質(zhì)大學(xué)學(xué)報：地球科學(xué)，2011，36（6）：1134-1142.

Xu J，Dong X J.Genetic types of large-scale landslides induced by Wenchuan Earthquake［J］.Journal of China University of Geosciences：Earth Science，2011，36（6）：1134-1142.

［26］鄭穎人，葉海林，黃潤秋.地震邊坡破壞機(jī)制及其破裂面的分析探討［J］.巖石力學(xué)與工程學(xué)報，2009，28（8）：1714-1723.

Zheng Y R，Ye H L，Huang R Q.Analysis and discussion of failure mechanism and fracture surface of slope under earthquake ［J］.Chinese Journal of Rock Mechanics and Engineering，2009，28（8）：1714-1723.

［27］鄭穎人，葉海林，黃潤秋，等.邊坡地震穩(wěn)定性分析探討［J］.地震工程與工程振動，2010，30（2）：173-180.

Zheng Y R，Ye H L，Huang R Q，et al.Study on the seismic stability analysis of a slope ［J］.Journal of Earthquake Engineering and Engineering Vibration，2010，30（2）：173-180.

（編輯王秀玲）