李 晶 朱亞杰 張嘉鑫

(1.中國石油工程建設(shè)公司華東環(huán)境巖土工程分公司，山東 青島 266071；2.青島理工大學(xué) 山東省地質(zhì)環(huán)境與效應(yīng)工程技術(shù)研究中心，山東 青島 266033)

位移矢量角在堆積層滑坡預(yù)測預(yù)報中的應(yīng)用研究

李 晶1朱亞杰1張嘉鑫2

(1.中國石油工程建設(shè)公司華東環(huán)境巖土工程分公司，山東 青島 266071；2.青島理工大學(xué) 山東省地質(zhì)環(huán)境與效應(yīng)工程技術(shù)研究中心，山東 青島 266033)

引入能反映滑坡坡體受力方向變化規(guī)律的特征參數(shù)——位移矢量角，運用有限元程序ANSYS進行數(shù)值模擬，研究堆積層滑坡在不同部位的位移矢量角變化趨勢，為滑坡預(yù)測預(yù)報提供參考依據(jù)。

堆積層滑坡，位移矢量角，數(shù)值模擬，預(yù)測

0 引言

滑坡是地殼表層巖體的一種地災(zāi)現(xiàn)象，是一種多發(fā)性的地質(zhì)災(zāi)害，給城市規(guī)劃、交通建設(shè)、水利、礦山開采等造成巨大損失，據(jù)有關(guān)數(shù)據(jù)，每年自然災(zāi)害給我國造成的直接損失高達千億元，其中滑坡等占20%左右[1]，滑坡災(zāi)害的防治問題亟待解決。

目前，國內(nèi)外一些學(xué)者對位移矢量角參數(shù)在堆積層滑坡預(yù)測預(yù)報中的應(yīng)用進行了相關(guān)的研究，并取得了一定的研究成果，其中的代表性研究方法主要有數(shù)理統(tǒng)計分析法和有限元分析法。例如，賀可強等[2]對邊坡塑性滑移階段的位移矢量角及其與滑坡穩(wěn)定性的關(guān)系進行了系統(tǒng)的分析，運用數(shù)理統(tǒng)計的趨勢位移分析原理構(gòu)建邊坡位移矢量角和位移速率參數(shù)統(tǒng)計預(yù)測預(yù)報判據(jù)；陽吉寶等[3]運用灰色理論，針對堆積層滑坡變形破壞的特點，對滑坡位移觀測數(shù)據(jù)進行了處理，得出相關(guān)灰色位移矢量角。土體擴容對邊坡穩(wěn)定性的影響方面研究，國外學(xué)者Mazarin[4]利用有限元的思想，經(jīng)過研究發(fā)現(xiàn)邊坡失穩(wěn)區(qū)域邊緣左側(cè)節(jié)點的豎直位移在體力增加到最終值時突然增大，而位于破壞區(qū)右側(cè)的節(jié)點，則表現(xiàn)出很小的變化，豎直位移突然增大說明邊坡垂直位移矢量角在邊坡失穩(wěn)時將要發(fā)生突變；國外學(xué)者Kim[5]利用有限元對邊坡進行分析中也發(fā)現(xiàn)邊坡在臨界失穩(wěn)時其位移速率矢量場發(fā)生突變，其垂直位移矢量角有明顯偏轉(zhuǎn)。

1 邊坡位移矢量角的定義及其與邊坡穩(wěn)定性的關(guān)系

邊坡位移矢量角，是邊坡位移矢量與水平面的夾角，即為垂直位移矢量角。對滑坡原始位移監(jiān)測資料中的垂直和水平位移進行二次再生處理，抽象出原始資料的規(guī)律性，便于總結(jié)邊坡運動運移規(guī)律。令x(0),y(0)分別為水平和豎直位移序列，x(1),y(1)為它們的累加序列[3,6]：

(1)

(2)

(3)

(4)

(5)

即可得到所設(shè)置觀測點的位移矢量角數(shù)據(jù)序列。

位移矢量角是一個矢量，具有方向和大小兩個數(shù)量值，兩個數(shù)量值構(gòu)成了位移矢量場。邊坡位移矢量場是邊坡運行運移狀態(tài)的直接數(shù)值反映，同時也反映了邊坡的穩(wěn)定狀態(tài)。邊坡的不同穩(wěn)定演變階段所具有的位移矢量角表現(xiàn)著不同的變化特征，因此，可以利用位移矢量角的變化規(guī)律來研究邊坡的穩(wěn)定性的演化規(guī)律[2]。

2 堆積層滑坡位移矢量角變化規(guī)律數(shù)值分析

1)屈服準(zhǔn)則。Drucker-Prager準(zhǔn)則[7]是目前巖土工程領(lǐng)域中常用的屈服準(zhǔn)則，同時也是巖土工程問題的解決中所應(yīng)用的大型有限元軟件所采納的，其表達式為：

(6)

2)滑坡數(shù)值模型的建立。根據(jù)地質(zhì)調(diào)查所確定的某滑坡的某典型剖面，進行平面非線性有限元分析，將坡體作為變形體，按照土的變形特性，計算土坡內(nèi)的應(yīng)力分布?；麦w后緣取至公路切坡導(dǎo)致的剪出口處，底部深度延至滑體高度的3倍左右，以滿足計算精度要求。坐標(biāo)系為：垂直向上為y軸的正向，水平指向河床為x軸的正向，原點為滑坡體后緣起點，按照平面應(yīng)變建立模型。邊界條件：在滑體前后緣的垂直邊界上為水平約束，底部約束條件為雙向固定約束。計算模型采用四邊形節(jié)點(Plane82)單元，有限元劃分有2 547個單元，節(jié)點總數(shù)為7 880，網(wǎng)格劃分情況如圖1所示。

3)滑坡位移矢量角變化規(guī)律分析。按照強度折減法的思想，對土的原始強度指標(biāo)進行一系列的折減，將折減后的c,φ值代入程序進行有限元計算，得到邊坡的位移、應(yīng)力的數(shù)值解。經(jīng)多次迭代運算，當(dāng)c=1.5×108Pa，φ=20°時，程序恰好不收斂，邊坡達到了臨界狀態(tài)。

假設(shè)此時的折減系數(shù)是1.0，根據(jù)現(xiàn)場實際情況，在滑坡的前部、中前部、中部、中后部及滑坡后部各設(shè)置了兩個觀測點；根據(jù)模擬計算結(jié)果，得出各節(jié)點位移矢量角如表1所示。

表1 滑坡各監(jiān)測點的位移數(shù)據(jù)及位移矢量角

從各節(jié)點位移矢量角來看，數(shù)據(jù)比較符合牽引式滑坡的特征，前端位移矢量角較小，說明前端位移以水平位移占優(yōu)，位于邊坡頂端的290節(jié)點甚至出現(xiàn)負角，說明該邊坡頂端滑動出現(xiàn)上揚現(xiàn)象，這是在剪出口位置才可能發(fā)生的情況，而事實上290節(jié)點正位于剪出口位置。這說明了模型參數(shù)設(shè)置合理，得到的數(shù)據(jù)比較可信，可以作為分析計算的依據(jù)。從表1可以看出，邊坡中前部位移矢量角均較小，邊坡已經(jīng)開始出現(xiàn)險情，治理迫在眉睫。

由于各個折減系數(shù)下的模擬計算是重復(fù)進行的，所以其他折減系數(shù)下的邊坡位移矢量角僅給出計算結(jié)果。根據(jù)D-P準(zhǔn)則，對c，φ值進行折減，折減系數(shù)即為折減系數(shù)F。以折減系數(shù)為1.0為基點，折減計算出其他幾個折減系數(shù)下對應(yīng)的c,φ值，見表2。

表2 在各折減系數(shù)下所對應(yīng)的c,φ值

將表2中各c,φ值代入模型進行模擬計算，給出在各折減系數(shù)下的所設(shè)各觀測點的位移矢量角(見表3)。

表3 邊坡各監(jiān)測點在折減系數(shù)下的位移矢量角

總體上看，位移矢量角是隨折減系數(shù)的變化而改變的，圖 2給出了邊坡模型前緣的位移矢量角隨邊坡穩(wěn)定性演化的曲線圖，隨著邊坡折減系數(shù)的降低，位移矢量角也在減小，在邊坡折減系數(shù)低于1.0時，位移矢量角急劇減小。當(dāng)折減系數(shù)達到0.95時，程序計算已經(jīng)不收斂了，節(jié)點290的值已經(jīng)為0，說明邊坡已經(jīng)破壞，而節(jié)點25的位移矢量角也急劇減小，并當(dāng)折減系數(shù)為0.9時也不收斂了。根據(jù)前面所述理論，邊坡位移矢量角的減小，說明土體有向臨空方向運動的趨勢，在前緣表現(xiàn)為隆起、剪出[9]。

邊坡中部、后部觀測點的位移矢量角隨折減系數(shù)變化曲線圖如圖3,圖4所示。

從總體趨勢來看，隨著時間的推移，c,φ值也逐漸降低，位移矢量角也在逐漸降低，即位移矢量角是隨折減系數(shù)的降低而降低的，與邊坡前緣相比，邊坡前中前緣、中部以及中后緣的位移矢量角的變化突發(fā)性依次減小。從數(shù)值上看，當(dāng)c,φ值降低到某個程度時，即折減系數(shù)低于1.0時，位移矢量角開始急劇減小，當(dāng)折減系數(shù)調(diào)整為0.93時，位移矢量角數(shù)值急轉(zhuǎn)直下，當(dāng)折減系數(shù)調(diào)整為0.90時，程序計算時絕大數(shù)節(jié)點不收斂。

3 結(jié)語

通過上述分析研究，可以得出以下結(jié)論：1)基于邊坡位移矢量場的原理，結(jié)合數(shù)學(xué)建模進行數(shù)值分析，通過對邊坡位移矢量場動態(tài)變化的分析和研究，能夠?qū)Χ逊e層滑坡的穩(wěn)定性進行評價。2)在彈塑性力學(xué)分析方法基礎(chǔ)上，運用邊坡的位移信息分析研究堆積層滑坡不同部位的位移矢量角的變化趨勢，以位移矢量角的形式反映邊坡穩(wěn)定性的變化規(guī)律，在對滑坡治理方案進行優(yōu)化設(shè)計上有數(shù)理支持價值。

[1] 黃潤秋，許 強.工程地質(zhì)廣義系統(tǒng)科學(xué)分析原理及其應(yīng)用[M].北京：地質(zhì)出版社，1997.

[2] 賀可強，陽吉寶，王思敬.堆積層邊坡表層位移矢量角及其在穩(wěn)定性預(yù)測中的作用與意義[J].巖石力學(xué)與工程學(xué)報,2003,22(12):1976-1983.

[3] 陽吉寶，鐘正雄.位移矢量角在堆積層滑坡時間預(yù)報中的應(yīng)用[J].山地研究，1995，13(1):49-54.

[4] Mazarin M T,Nour M A.Significance of soil dilatancy in slope stability analysis[J].Journal of Geotechnical and Geo-environmental Engineering,ASCE,2000,126(1):76-79.

[5] Kim J.Limit analysis of soil slope subjected to pore water pressures[J].Journal of Geotechnical and Geo-environmental Engineering,ASCE,1999,125(1):49-57.

[6] 賀可強，陽吉寶.堆積層邊坡位移矢量角的形成作用機制及其與穩(wěn)定性演化關(guān)系的研究[J].巖石力學(xué)與工程學(xué)報，2002，21(2)：185-192.

[7] 王勖成，邵 敏.有限單元法基本原理與數(shù)值方法[M].第2版.北京：清華大學(xué)出版社，2002.

[8] 龔曉南.土工計算機分析[M].北京：中國建筑工業(yè)出版社，2000.

On application of displacement vector angle in forecasting of accumulation landslide

Li Jing1Zhu Yajie1Zhang Jiaxin2

(1.EastChinaEnvironmentalGeotechnicalBranch,ChinaPetroleumEngineering&ConstructionCorp,Qingdao266071,China;2.ShandongGeologicalEnvironmentandEffectEngineeringTechnicalResearchCenter,QingdaoUniversityofTechnology,Qingdao266033,China)

The paper leads the displacement vector angle, the characteristic parameter, which can reflect the law of changes for stressed direction of landsclides, adopts the finite element program, ANSYS, to undertake the numeric simulation, and researches the change tendency for the displacement vector angles of accumulation landslides, so as to provide some reference for forecasting the landslides.

accumulation slope, displacement vector angle, numerical simulation, forecasting

2015-01-08

李 晶(1980- )，男，碩士，工程師； 朱亞杰(1981- )，男，工程師； 張嘉鑫(1989- )，男，在讀碩士

1009-6825(2015)08-0090-03

P642.22

A