巖土參數(shù)反算在工程應用中的分析

王書劍

（1.長沙理工大學交通運輸工程學院，湖南 長沙 410076；2.長沙市公路管理局城郊路政執(zhí)法大隊，湖南 長沙410007）

0 引言

目前用來分析巖土工程穩(wěn)定性的方法有很多，如數(shù)值分析法、極限平衡法等。在對巖土體進行穩(wěn)定性分析時，準確的強度參數(shù)是至關重要的。工程中經常用土工試驗來獲取內摩擦角φ、黏聚力 c等強度參數(shù)。該方法得出的試驗結果可靠，且適用于大多數(shù)巖土體。但對構造復雜的土體，如碎石土等，若用土工試驗來獲取其參數(shù)，所耗時間長，且難度較大。并且有些土體強度參數(shù)的空間變異性較大，得到的試驗結果離散，很難得到準確的計算參數(shù)[1-2]。在此情況下，通過采用參數(shù)反演手段來獲取巖土體強度參數(shù)已得到廣泛的應用[3]。準確分析和評價邊坡穩(wěn)定性的關鍵是如何正確選取土工參數(shù)，進而會影響到邊坡工程的安全。而這些土工參數(shù)中，滑動面巖土體的重度γ，其選用值的誤差對邊坡穩(wěn)定性的影響遠小于粘聚力c、內摩擦角φ。因此，滑坡巖土體抗剪強度參數(shù)c、φ的取值成為影響邊坡穩(wěn)定性分析結果可靠度的關鍵，研究確定c、φ值的方法對邊坡工程的穩(wěn)定性分析與加固設計具有重大的實際意義。

1 原理分析

1.1 反算分析法

在國內外邊坡穩(wěn)定性分析通常采用極限平衡條分法，主要有瑞典法、Janbu法、簡化Bishop法和傳遞系數(shù)法等，其中傳遞系數(shù)法在滑面邊坡的應用最為廣泛。傳遞系數(shù)法不計力矩平衡條件，只考慮力的平衡，簡單易行，而且與其他極限平衡方法的計算結果相差不大，能基本滿足邊坡工程實踐的需求。反演分析采用的邊坡穩(wěn)定性分析模型不同，反算得到的強度參數(shù)值也會有所不同。

1.1.1 反分析基本前提

反演分析的可靠性主要取決于邊坡穩(wěn)定狀態(tài)和邊界條件的準確性。反演分析的基本前提：①結合實際情況選定邊坡穩(wěn)定性系數(shù)值；②勘察滑動面位置及滑動形式；③總結分析造成邊坡失穩(wěn)或坡體變形的外力因素。

1.1.2 反分析求解過程

對于確定的滑坡斷面，先假定穩(wěn)定系數(shù)Fs已知，利用傳遞系數(shù)法可以求取c、φ值的對應關系。但c、φ值是兩個不確定參數(shù)，對于單一斷面的一個極限平衡方程，c、φ值具有不確定性，必須先假定c或φ中一個，再求解另一個未知數(shù)，這樣可求出無數(shù)組c、φ值的解，但卻無法獲得反映滑坡滑動時真實狀態(tài)的唯一解。為了同時反算出c、φ值，可在滑坡體中尋找兩個或多個處于極限平衡狀態(tài)的斷面，建立不少于兩個極限平衡方程聯(lián)立求解。通常情況下，為了提高反算結果的可信度，可同時結合土工試驗、工程類比或敏感性分析等方法來確定強度參數(shù)的取值。

聯(lián)立滑坡的兩個或多個斷面進行求解，可能會出現(xiàn)兩條或多條平行直線方程的情況，以致無法求解。因此，多斷面極限平衡方程聯(lián)立反算的基本條件是選取的多個滑坡斷面必須相似。相似條件包括：①地質條件相似，滑坡的類型和形式、滑帶土的物質組成以及含水狀態(tài)要相似；②滑坡的運動過程和狀態(tài)相似。

將坡體將滑而未滑那一瞬間視為滑坡的極限平衡狀態(tài)，在未發(fā)生滑動時，其穩(wěn)定系數(shù)Fs=1。采用傳遞系數(shù)法列出滑坡斷面的極限平衡方程，求解出多個斷面的一系列不同c、φ值組合，釆用圖解法繪制多個斷面的c-φ曲線，交點即為反算的強度參數(shù)值。

圖1 c、φ值反分析圖解示意圖

1.2 基于傳遞系數(shù)法的反算分析

傳遞系數(shù)法（即不平衡推力法）是滑坡穩(wěn)定性分析使用最廣泛的方法之一。傳遞系數(shù)法假定滑坡體為整體滑動，各條塊的條間力以集中力表示，作用在分界面中點，其作用方向與上一條塊的底面平行。在對滑坡范圍內滑動方向和速度大體一致的滑體進行計算時，可將它們視為一個計算單元，并選擇其中的典型滑坡斷面進行分析。根據不同地層性質和滑面坡度的差異，將整個滑體劃分為若干豎直條塊。由力的平衡條件可知，從第一個條塊連續(xù)向前逐條向下推求，直至最后一個條塊的推力為零，該方法通常有隱式解法和顯式解法兩種[4]。傳遞系數(shù)法有三個假定條件：①滑體不可壓縮并作整體下滑；②條塊間只傳遞推力不傳遞拉力；③取垂直滑坡主軸方向單位寬度(一般為1m)的巖土體作為計算斷面，不考慮條塊兩側的摩阻力[2]。

圖2 第i個條塊受力示意圖

取滑坡體的第i個條塊為研究對象，條塊i的受力圖如圖2所示。由力平衡條件，可得第i條塊的條間力Pi表達式：

其中，

ψi為傳遞系數(shù)，可簡化計算為：

上述式子中的c和tanφ均為滑面實際強度參數(shù)，強度不再發(fā)生變化。采用下滑力超載穩(wěn)定系數(shù)，即下滑力Ti前乘了穩(wěn)定系數(shù)Fs；然而傳遞系數(shù)ψi可簡化成一個常數(shù)，且與Fs無關。即設Fs=1，由此推出安全系數(shù)的顯式表達式為：

上式是《公路路基設計規(guī)范》中常用的公式。在滑坡處于極限平衡狀態(tài)下，傳遞系數(shù)ψi的表達式中可假設Fs=1?；趥鬟f系數(shù)法進行邊坡滑動面強度參數(shù)反算分析時，計算量龐大，可利用MATLAB數(shù)學軟件進行計算，這樣可以大大提高計算效率[5]。

2 實例計算

2.1 邊坡概況

永龍界ZK56+860～ZK56+940段邊坡巖體較破碎松散，覆蓋較厚強風化層。邊坡開挖成型后，受降雨影響，坡體有兩處小面積滑塌，滑帶巖土體主要為表層粘土（夾雜少量碎石）、強風化層軟巖、泥質頁巖等。以下為該路的2處典型滑坡示意圖：

圖3 1#滑坡體示意圖

圖4 2#滑坡體示意圖ZK56+860～ZK56+940邊坡1#、2#滑坡體示意圖

該邊坡位于斷層碎裂帶邊緣，地質狀況復雜，受地質構造及外界條件影響很大，邊坡整體穩(wěn)定性較差。巖質邊坡地處斷層碎裂帶邊緣，受地質構造影響顯著，巖體完整度較差，塊狀偏小，巖體風化較嚴重，抗剪切強度較小。巖體以節(jié)理、裂隙發(fā)育的強風化砂質頁巖和粉質粘土中夾強風化頁巖塊石兩大類為主，強風化層厚度大，巖質較軟，整體結構松散。根據現(xiàn)場滑坡的形態(tài)特征、滑帶巖土體及結構特征和勘探的地質資料，選取ZK56+860～ZK56+940段邊坡為研究對象，該處邊坡坡體中碎石土較多，土體強度參數(shù)的空間變異性較大，對室內直剪試驗結果產生較大的離散性，很難獲取準確的參數(shù)。影響邊坡穩(wěn)定性的因素較多，如巖土體的基本物理力學性質、邊坡的幾何特征、巖體的結構面參數(shù)等。它們對邊坡穩(wěn)定性都有不同程度的影響。由于這些影響因素難以量化，在進行邊坡的穩(wěn)定性分析時，很難考慮到它們的影響。而反算分析方法的思想就是把這些難以量化的因素綜合考慮在“等效參數(shù)”中，因此采用了基于傳遞系數(shù)法來反算滑動面巖土體的抗剪強度參數(shù)，為邊坡工程分析發(fā)展提出了一種新的方法。

2.2 參數(shù)分析

假定邊坡處于極限平衡狀態(tài)，采用傳遞系數(shù)法對ZK56+860～ZK56+940邊坡巖土體參數(shù)進行敏感性分析。依據勘測資料提供的c、φ取值范圍，形成不同的組合系列，利用MATLAB軟件計算對應的穩(wěn)定性系數(shù)Fs，計算結果如表1、表3所示。

表1 c、φ與Fs的影響關系

圖5 Fs-c-φ關系曲線圖

由上圖可通過擬合曲線得到y(tǒng)=k1x+b公式，得出不同內摩擦角的情況下對邊坡穩(wěn)定性影響的敏感度k1，如下表2所示：

表2 φ與k1的影響關系

表3 c、φ與Fs的影響關系

圖6 Fs-c-φ關系曲線圖

由上圖可通過擬合曲線得到y(tǒng)=k2x+b公式，得出不同粘聚力的情況下對邊坡穩(wěn)定性影響的敏感度k2，如下表4所示：

表4 c與k2的影響關系

由圖5、圖6可知，滑動面巖土體的c和φ值對穩(wěn)定性系數(shù)Fs的影響均較大，比較表2與表4可知，斜率k2＞k1，表明對于ZK56+860～ZK56+940邊坡滑坡體，穩(wěn)定性系數(shù)Fs對φ值的影響更敏感，且隨內摩擦角φ的變大，穩(wěn)定性系數(shù)Fs增加速率更快。因此，在進行反演分析時可先假定c值，然后再通過MATLAB軟件的極限平衡穩(wěn)定分析程序來反算φ值。

2.3 反算參數(shù)的確定

根據邊坡穩(wěn)定性傳遞系數(shù)法的穩(wěn)定系數(shù)計算公式，利用MATLAB軟件編寫了邊坡穩(wěn)定性分析電算程序，考慮 “自重+水”的工況，對ZK56+860～ZK56+940邊坡進行滑動面抗剪強度參數(shù)c和φ值的反演計算。由以上的邊坡參數(shù)敏感性分析可知c值的敏感度較低，根據室內試驗數(shù)據可確定c值的范圍，將c值進行區(qū)間分段，其中選取c值分別為 13、14、15、16、17、18、19、20kPa， 在滑坡的極限平衡狀態(tài)及邊界條件下，運用傳遞系數(shù)法反算φ值。當穩(wěn)定系數(shù)Fs=1.0時，分別反算出ZK56+860～ZK56+940邊坡1#和2#滑坡體在不同c值條件下的φ值。計算結果如表5、6所示。

表5 1#滑坡體傳遞系數(shù)法反算結果

表6 2#滑坡體傳遞系數(shù)法反算結果

圖7 1#和2#滑坡體c-φ曲線反演分析

聯(lián)立ZK56+860～ZK56+940邊坡1#和2#兩個滑坡體極限平衡方程，求解該邊坡滑動面抗剪強度參數(shù)反演計算值。如圖7中兩條c-φ曲線的交點坐標c=15kPa，φ=24.48°即為該邊坡滑動面抗剪強度參數(shù)反演計算值。

3 結論

本文介紹了獲取滑動面巖土體抗剪強度參數(shù)的反算分析法，以ZK56+860～ZK56+940邊坡兩處滑坡體為研究對象，基于傳遞系數(shù)法對該滑坡體的抗剪強度參數(shù)進行了反演計算，得到合理參數(shù)。通過分析影響邊坡穩(wěn)定性參數(shù)內摩擦角φ和黏聚力 c對邊坡穩(wěn)定性的影響得出，內摩擦角φ較黏聚力c對邊坡的影響更為敏感。在反演計算中，假定穩(wěn)定系數(shù)Fs=1.0時，分別反算不同c值條件下的φ值，聯(lián)立兩處滑坡體極限平衡狀態(tài)方程，得出抗剪強度參數(shù)為c=15kPa，φ=24.48°，為龍永高速后續(xù)邊坡開挖施工行為及防護加固時機選擇提供指導性建議。

［1］沈永飛.邊坡位移反分析及其工程應用研究[D].重慶大學,2010.

［2］童志怡.巖質邊坡滑動面力學參數(shù)的取值理論和方法研究[D].中國科學院研究生院(武漢巖土力學研究所),2009.

［3］孫志彬.Mohr-Coulomb準則下基于滑動面深度的邊坡參數(shù)反分析方法[M].巖土力學,2014,4.

［4］賈恒國.抗滑樁設計計算的程序實現(xiàn)[D].太原理工大學,2010,5.

［5］李清.基于極限平衡理論的滑坡抗剪強度參數(shù)反演數(shù)據表法研究[M].路基工程,2015,4.