雙排抗滑樁在滑坡治理中的支護優化探討

甘怡，楊波

（1第三軍醫大學校務部營房處，重慶 400041；2重慶第十建設有限公司，重慶 400041）

0 前言

埋入式抗滑樁是指樁頂在地下一定深度的抗滑樁，是一種常見的抗滑結構形式，常用于滑坡的治理。埋入式抗滑樁的優點是能節省大量材料，有很好的經濟性。使用埋入式抗滑樁來治理滑坡工程時，比較重要的就是確定合理的支護方案。通常所說的埋入式抗滑樁的合理支護方案是指在滿足一定要求（達到一定安全系數）下的，確定抗滑樁的埋深和在滑坡體中所處的水平位置。對于單排埋入式抗滑樁的支護，唐曉松等人利用遺傳算法嵌套支持向量機已經做了研究[1]。但是關于雙排埋入式抗滑樁的合理支護的研究還較少，其難點是需要確定雙排樁各自的埋深和水平位置。本文利用遺傳算法對雙排埋入式抗滑樁的合理支護方案進行了優化探討。

1 雙排埋入式抗滑樁優化思路

1.1 雙排埋入式抗滑樁優化的基本原則

雙排埋入式抗滑樁的優化涉及的因素較多，不僅受到滑坡參數及外力作用的影響，還與樁本身的位置與材料等有關，比較復雜。通過研究，我們提出雙排埋入式抗滑樁合理支護優化的兩個基本原則。第一個基本原則是雙排埋入式抗滑樁各排樁推力和最小。因為抗滑樁的優化就是降低工程造價，節約材料。而抗滑樁的配筋多少、截面大小又決定于滑坡推力的大小，所以采用各排樁的推力作為重要的參考因素。第二個基本原則是治理后的安全系數最大化。因為滑坡治理的首要任務是安全可靠，坡體的穩定安全系數必須要大于滑坡治理中設計要求的安全系數，而在設計中，為了確保安全，一般按照穩定安全系數最大化原則。

1.2 雙排埋入式抗滑樁優化的步驟

雙排埋入式抗滑樁的優化就是要確定雙排樁各自的埋深和水平位置，涉及4個參數。通過研究我們發現，雙排埋入式抗滑樁中的前排樁對安全系數的貢獻很大，前排樁如果設置成埋入式，常常會導致滑坡安全系數很小，達不到設計要求，所以雙排埋入式抗滑樁的前排樁應設計為全長式抗滑樁，后排樁設計為埋入式抗滑樁。所謂全長式抗滑樁，是指樁頂在地面的抗滑樁。因此，雙排埋入式抗滑樁的優化就是確定3個參數：雙排樁各自的水平位置和后排樁的埋深（前排樁為全長式，埋深為0）。

雙排埋入式抗滑樁的優化一般可以分成兩步來進行。第一步是雙排埋入式抗滑樁各自水平位置的優化。先不考慮埋深，考慮雙排樁都為全長式抗滑樁時的情況。計算出雙排全長式抗滑樁的一組樣本數據，然后取定安全系數最大時的位置為雙排樁的最優水平位置。第二步是雙排抗滑樁后排樁埋深的優化，取第一步中計算出的雙排樁的最優水平位置，變化后排樁的埋深，應用遺傳算法嵌套支持向量機的方法對此位置的后排樁進行埋深優化。

1.3 雙排埋入式后排樁埋深的優化方法

通過MATLAB支持向量機（SVM）工具箱LS_SVMlab形成設計參數與滑坡安全系數及埋入式抗滑樁的推力間的映射關系模型[2-3]。根據這一數值模型，對于任一組給定的支護參數（后排樁的埋深），均可以用支持向量機泛化推廣的功能得到相應的安全系數和樁所受的推力。然后，把支持向量機嵌套入遺傳算法中，以支護加固后滑坡體需要達到的最小安全系數等作為約束條件，以樁所受推力作為評價指標，給定搜索的范圍，即可對支護參數進行尋優。

2 映射關系和樣本的建立

要得到安全系數與支護參數以及抗滑樁推力之間的映射關系模型，首先需構建一定數量的基礎樣本，而且要保證樣本的準確性。傳統的計算方法（如不平衡推力法等）由于不能考慮樁-土之間的各種共同作用，除可以計算滑坡體的推力以外，其余有用的設計參數大都無法確定。近年來，有限元強度折減法得到了廣泛地應用，鄭穎人及其學生率先將有限元法應用于抗滑樁的設計，取得了良好效果[4-6]，這為抗滑樁的推力及對應的安全系數樣本的準確性提供了保證。因此，本文樁體采用實體單元、線彈性材料模擬，土體采用D-P彈塑性模型模擬，利用ANSYS有限元程序對雙排埋入式抗滑樁進行穩定性分析。

3 工程算例

3.1 計算模型與參數

該滑坡原型為陽平關滑坡，地處陜西省寶成線鳳州-陽平關段，滑坡編號DK271+255。該條線路處于秦巴山區和陜北黃土高原溝壑地區，若遇暴雨、連陰雨有可能發生滑坡、崩塌、泥石流等災害，歷史上多次發生滑坡傷人事件。該類滑坡順層面或結構面滑動，主滑段為直線型，前緣抗滑段較小、阻滑特征不明顯。計算模型為此滑坡的典型剖面，滑坡長170m，滑坡高105m，前排樁長25m，寬4m。后排樁長45m，寬4m。計算參數如表1所示，計算模型如圖1所示。

表1 巖土體和樁的物理力學參數

表2 不同排距下兩排樁的推力表（單寬）

圖1 簡化計算模型

3.2 模型的優化

雙排埋入式抗滑樁的支護方案優化可以分成兩步來進行。第一步是雙排埋入式抗滑樁水平位置的優化，第二步是后排抗滑樁埋深的優化。

3.2.1 雙排埋入式抗滑樁水平位置的優化

先用有限元強度折減法建模計算得到排距變化時的雙排樁的推力，如表2所示。安全系數最大值處在間距90m的地方，在此處安全系數為1.54，雙排埋入式抗滑樁的水平位置確定。在這里采用安全系數最大化準則有如下兩點考慮：一是安全系數最大，在一定程度上來說更加保險和安全；二是為后排樁的埋入留下余地。因為隨后排樁埋入深度的增加，安全系數會逐漸變小，直到降低到設計安全系數。而且經過研究我們發現，埋入式抗滑樁樁上的推力與抗滑樁在滑體中的有效長度成正比。而且樁在滑坡體中的有效長度對推力的影響比抗滑樁的水平位置對推力的影響要大，說明增大埋深對減小兩排樁上的推力和效果更好。

3.2.2 雙排埋入式抗滑樁后排樁埋深的優化

首先應用有限元強度折減法進行建模，改變雙排埋入式抗滑樁后排樁的埋深，計算出對應的樁所受推力和滑坡體的安全系數，共得到了9組樣本，完成樣本的構建，如表3所示。 然后調用initlssvm函數對模型進行初始化，通過trainlssvm對模型進行訓練，得到后排埋入式抗滑樁的埋深與滑坡的安全系數及雙排樁所受推力的映射關系。最后調用simlssvm函數通過回歸來驗證所建立的映射關系的準確性，得到樣本的輸入值和目標值如表4所示。

表3 不同排距下兩排樁的推力表（單寬）

表4 兩排樁樁頂中心處的豎直坐標

表5 支持向量機的目標值與預測值

為了驗證所建立的映射關系的正確性，本文利用支持向量機泛化推廣的功能，根據文獻[7]中的方法對雙排抗滑樁的推力和數值進行了預測，并和文獻方法中的數值模擬結果進行了對比，如表5所示。從表5中可以看出，預測值和數值模擬得到的目標值十分接近（誤差小于4.00%），由此說明通過支持向量機建立的映射函數關系是符合精度要求的。

為進一步驗證所建立的映射關系的準確性，這里任意取一組坐標（x=11，y=0）作為輸入值，通過有限元強度折減法建立AYSYS模型，計算得到的雙排埋入式抗滑樁的推力之和的實際值為9400kN。而利用支持向量機泛化推廣的功能得到的雙排埋入式抗滑樁的推力之和的預測值為9430kN，兩者非常接近，如表6所示。

在建立了后排埋入式樁的埋深與安全系數及目標函數推力間的映射關系后，就可以根據已經建立的雙排埋入式抗滑樁的數學模型采用遺傳算法進行尋優了。由于遺傳算法只能搜索矩形區域，因此需要把擬搜索的區域近似為若干個矩形。目前對于埋入式抗滑樁的合理支護的研究尚處于欠成熟階段，所以保守地取設計安全系數為1.4。

通過遺傳算法進行尋優，當樁頂中心處水平坐標x=10.3，y=0時，目標函數（雙排埋入式抗滑樁的推力和）最小，數值為9452kN，此時安全系數為1.4。采用有限元強度折減法用AYSYS建立同樣的滑坡模型得到的滑坡體的安全系數等于1.41，計算得到雙排埋入式抗滑樁所承受的推力和等于9300kN，和遺傳算法全局尋優得到的最優解比較接近，此時滑坡體潛在滑移面如圖2所示。通過兩歩的優化，得到埋入式抗滑樁的最優位置為如圖2所示，后排樁距離前排樁90m，后排樁的埋深10.3m。此時前排樁的推力為4538kN，后排樁的推力為4828kN。

圖2 雙排埋入式抗滑樁最優時的塑性應變云圖

表6 預測值和實際值的對比分析

4 結論

應用遺傳算法嵌套支持向量機的方法，通過在一定區域內尋優，對雙排埋入式抗滑樁的合理的支護方案進行了優化研究。建議雙排埋入式抗滑樁合理支護優化的二個基本原則，一是雙排埋入式抗滑樁各排樁推力和最小，二是治理后的安全系數最大化。并通過具體滑坡的優化實例分析，給出了確定雙排埋入式抗滑樁合理支護方案的優化步驟，即先確定雙排埋入式抗滑樁的最優水平位置，再保持前排樁全長式，優化后排樁的埋深，得到最終的優化方案。為雙排埋入式抗滑樁在滑坡工程中的合理支護做了新的理論嘗試。

[1]唐曉松，鄭穎人，段永勝，等.人工智能在埋入式抗滑樁中的應用[J].地下空間與工報，2010，15（4）:120-126.

[2]許建華，張學工，李衍達.支持向量機的新進展[J].控制與決策，2004，19（5）:481-493.

[3]馬文濤.支持向量機在膨脹土分類中的應用[J].巖土力學，2005，26（11）:1790-1792.

[4]鄭穎人，趙尚毅，鄧楚鍵.有限元極限分析法發展及其在巖土工程中的應用研究[J].中國工程科學，2006，8（12）:35-61.

[5]鄭穎人，趙尚毅.用有限元強度折剪法求滑（邊）坡支擋結構的內力[J].巖石力學與工程學報，2004，23（20）:3552-3558.

[6]鄭穎人，趙尚毅.巖土工程極限分析有限元法及其應用[J].土木工程學報，2005，38（1）:91-99.

[7]雷文杰，鄭穎人，馮夏庭.滑坡治理中抗滑樁樁位分析[J].巖土力學，2006，6（27）:950-954.