4.3 有限元計算結果分析
   
    4.3.1 滑坡體應力計算結果分析
   
    每個計算剖面的表面拉應力區和關鍵點拉應力值如表2所示，表2中所給出的拉應力區是指從起點到終點連續出現拉應力的區域。
 

   
    滑坡體的應力計算可通過滑坡土體容重變化模擬天然狀態和飽水狀態。由于兩種狀態下土體容重相差不大且僅考慮材料的線彈性特性，因此，滑體內部應力在兩種狀態下也相差不大。從各滑坡體拉應力計算結果來看，各滑坡體關鍵剖面的表面不同高程處均會出現連續拉應力區或者局部拉應力點，正是由于滑體表面存在這些連續拉應力區才會使得滑坡體出現拉應力裂縫，最終導致滑體產生滑動。
   
    滑坡體的拉應力區基本上分布在滑坡體的中后緣部位，表面最大拉應力出現在滑體的后緣，最大拉應力處是滑體在自重作用下最容易出現拉裂縫的地方。若遇降水，滑坡土體處于一定程度的飽水狀態，除土體自重增加以外，土體的抗剪強度也要下降，滑坡體的穩定性降低，上述拉應力區的任何部位都有出現拉裂縫的可能性。因此，在進行滑坡治理時應對各滑坡體出現拉應力的區域給予重視。
   
    4.3.2 滑坡體穩定安全系數計算分析
   
    根據有限元計算結果可以得到各個計算剖面的滑動面上受力情況，表3給出了主滑動面在不同工況下的具體受力計算結果。
   
    根據選擇的各計算剖面上抗滑力和下滑力大小可以采用式(3)計算滑坡體各計算剖面的穩定安全系數，式 (3) 中的Li、Fni和Fri二由有限元計算程序計算直接得到，各計算剖面的穩定安全系數計算結果如表4所示。同時為了驗證計算結果的合理性，對主滑動面穩定安全系數還進行了推力系數法計算，計算結果也列于表4。
  

  
    從上述計算結果可以看到，4＃滑坡體主滑動面平均傾角為32.06，滑床地勢較陡，表4中列出的計算剖面穩定系數均較小，該滑體是一個不穩定體，飽水狀態下穩定安全儲備更小。5個剖面的拉應力區分布較大，從人工陡坎向上到滑體后緣都是拉應力區，各個滑動剖面的最大拉應力較大且都處于滑體的最后緣，因此該滑體有兩種滑動的可能:一是整個滑坡體在滑體后緣形成拉裂縫整體下滑;二是人工陡坎以上坡體首先以陡坎為剪出口產生滑動，滑動后的堆積體堆積在施工便道以下的滑體上使其自重增加，從而可能產生推移式滑動。
   
    該滑坡體主滑動面穩定安全系數的有限元法計算結果和推力系數法計算結果比較接近，表明將有限元法應用于滑坡穩定性分析中是可行有效的。同時，本文所建立的滑坡穩定性分析有限元方法中考慮了滑動面上應力釋放與應力重分配的迭代計算，使得有限元方法計算的主滑動面穩定安全系數小于推力系數法計算結果，表明本文的有限元方法計算結果更接近滑坡體實際受力狀態，滑坡體穩定安全系數更符合實際情況。
   
    5 結論
   
    通過滑坡穩定性有限元分析可以得到如下結論:
   
    (1)二次開發程序可以大大減小滑坡穩定性分析有限元法建模工作量。
    (2)有限元分析方法能夠有效地進行滑坡穩定性分析，可以計算滑坡體內部應力，并直接求解滑動面上的法向應力和切向應力。
    (3)采用硬彈簧和軟彈簧描述滑動面接觸摩擦模型能夠有效模擬滑坡體滑動過程。
    (4)考慮滑動面上出現裂紋后應力釋放與應力重分配的迭代計算，使得有限元方法計算的穩定安全系數更接近滑坡體實際受力狀態，滑坡體穩定安全系數更符合實際情況。
    (5)通過有限元計算能夠得到滑坡體內部拉應力區分布范圍，有效指導滑坡治理。