基于FLAC的路基邊坡穩定性分析

楊帆

摘 要：科學準確的邊坡穩定性分析是路基邊坡防治加固的重要決策依據，針對傳統路基邊坡穩定性分析方法存在誤差大、效率低下等弊端，本文梳理總結了提高路基邊坡穩定性分析結果精確性的技術要點，基于FLAC仿真軟件從主應力云圖、位移云圖和穩定性系數三個方面對構建的高邊坡仿真模型進行了計算與分析，據此系統評估了邊坡穩定性，為路基邊坡防治加固提供借鑒與參考。

關鍵詞：FLAC模型;路基邊坡;穩定性;仿真分析

中圖分類號：U416.14 文獻標識碼：A

0 引言

路基邊坡穩定性關系著整個公路的耐久、穩定與安全，邊坡坡度、地質巖性、地質結構、地下水發育情況以及施工等因素均對邊坡穩定性造成重要影響，而傳統路基邊坡穩定性分析方法存在誤差大、效率低下等弊端[1-3]，故有必要運用科學合理的路基邊坡穩定性分析方法，進一步確定合理的路基邊坡構型，為路基邊坡防治加固提供借鑒與參考。

1 路基邊坡穩定性分析概述

傳統路基邊坡穩定性分析方法存在誤差大、效率低下等弊端，為保證路基邊坡穩定性分析結果精確性，科學掌握斜坡的自然環境、地形和地質參數，以及合理運用匹配計算方法和參數指標尤為重要[4-6]。

1.1 參數值的合理化

路基邊坡穩定性分析需要考慮選取參數的合理性，相關分析顯示，計算方法在對邊坡穩定性的評估方面實際影響相對較小，但由于斜坡土本身存在不均勻的特性，其參數值選取的變化性必會對計算結果造成偏差，因此如何合理優化計算參數是未來研究重點。

1.2 分析方法的組合

目前尚未存在一種可以完全解決和分析所有問題的方法，將兩種或更多種方法結合起來相互融合運用是未來的發展趨勢。如對專家系統進行分析與設計時，可以利用遺傳算法進行推導和設計，并對函數選擇以及相關參數進行確定;另外，可充分組合利用有限元、離散元，數值解、解析解的組合等各自優勢解決較為復雜的梯度問題。

1.3 系統理論

路基邊坡是復雜的開放系統，故綜合考慮各個子系統之間的相互作用和聯系，并利用控制論和協同學等理論來穩定邊坡，據此研究和分析機制特性也是研究趨勢之一。

1.4 反分析方法運用

目前多依據梯度應力來研究梯度系統的不穩定性，若反向思維考慮，則可以從梯度系統的變形分析中獲得對梯度穩定性的分析，并進行穩定性評估。例如，反位移分析的方法是基于測得的位移值的反演獲得初始應力和參數，然后將其應用于工程實踐。

2 FLAC程序簡介

上述概述可知，合理運用計算方法對保證路基邊坡穩定性和安全性分析結果精確性尤為重要。FLAC軟件基于拉格朗日差分法構成，允許計算過程中材料的屈服和流變以模擬地面巖土所具有的力學屬性，故更適合解決模擬較大變形及扭曲，具有綜合性，速度性和經濟性等優點[7-8]，其相關計算流程如圖1所示：

3 公路路基邊坡穩定性分析

3.1 FLAC模型建立

某一級公路路基設計寬度為30 m，邊坡設計坡度為1：1。路基土體分為三層，基層為粘土，層厚10 m;中層為硬粘土，經壓實，層厚10 m;上層為中細沙土，層厚10 m;路基底部寬100 m。

建立幾何模型如圖2所示：

FLAC的邊界條件限定為限制節點位移，假設路基中心線左右各材料性質完全一致，中心線左右受力及位移正對稱，據此確定邊界約束條件為：地基土的下部邊界視為無位移的固定邊界，中心對稱面和側面各節點限制水平位移。

本文采用Mohr-Coulumb模型對邊坡穩定性進行模擬計算，此模型計算效率較高，適用于松散狀和黏結狀的粒狀散體材料，如土體、巖石、混凝土等。主要參數指標有：巖土體彈性模量E、泊松比v、剪脹角ψ、體積模量B、剪切模量S等均對穩定性分析結果造成影響。每個土壤層的物理和機械參數如下表1所示：

3.2 運行結果及分析

運行FLAC軟件自動完成仿真，獲取主應力云圖、位移云圖和穩定性系數，可依次對邊坡進行應力分析，位移分析和變形分析，據此評估邊坡的穩定性。

3.2.1 應力分析

輸出的主應力云圖如下圖3、4所示。

圖3和圖4是FLAC中的模擬應力云，分析上圖可知，顯示y方向上的垂直應力是壓應力，隨著土壤厚度的增加而增加，并且與場應力的變化一致。斜率周期性地從外到內逐漸增大，其應力值達到最小值或零，從其頂部到底部逐漸增大，并且應力達到最小值或零。x方向的側向應力有規律地分布在傾斜體上，并且隨著傾斜軸承的增加，最大壓應力值為150 kPa。拉應力也出現在傾斜的頂面下方的某些深度（下部中心）。路基斜坡的特征是在橫向方向上受拉和受壓，主要是受壓，但在垂直方向受壓。

3.2.2 邊坡位移

輸出的位移云圖如下圖5、6所示。

圖5和圖6分別是斜率的垂直位移圖和橫向位移圖，分析上圖可知，較低坡度在y方向上的最大位移達60 mm，發生在較低坡度的中心表面附近。隨著土壤層厚度的增加，位移值逐漸減小，直到可以忽略不計或為零。與原始路基設計斜坡相比，斜坡上沒有形成最大位移區域，這使斜坡變慢，從而有助于分散應力。x方向上的最大位移發生在傾斜軸承內部，最大位移值為7.5 mm。與原始路基設計坡度相比，最大位移區域從坡度附近的坡度移動到坡度內部以創建拐角。傾斜不會導致傾斜體的整體減小。與原始線性斜率相比，y和x方向上的階躍斜率的位移值分別減小了65.7%和94%。

3.2.3 穩定系數

FLAC軟件運行完成之后，由數值模擬的結果分析可知，穩定系數數據結果為1.58～1.59，邊坡穩定。

綜上，基于FLAC仿真軟件構建土壤土質高邊坡模型進行計算與分析，斜坡形狀從直線線性改變為折線線性，即設置臺階后y和x方向上的階躍斜率的位移值分別減小了65.7%和94%，且穩定系數數據結果為1.58～1.59，邊坡穩定。

4 結語

基于主應力云圖、位移云圖和穩定性系數三個方面對構建的高邊坡仿真模型進行了計算與分析，斜坡形狀從直線線性改變為折線線性，即設置臺階后y和x方向上的階躍斜率的位移值分別減小了65.7%和94%，且穩定系數數據結果為1.58～1.59，邊坡穩定。為路基邊坡防治加固提供借鑒與參考。

參考文獻：

[1]張寶龍，范文.土-巖混合邊坡的破壞模式研究[J].公路交通科技，2018（4）：33-39.

[2]吳忠廣，申瑞君，萬福茂，等.巖質高邊坡運營安全風險源辨識方法[J].公路交通科技，2018（3）：8-15+27.

[3]高連生，趙煉恒，李亮，等.帶臺階的多級邊坡穩定性上限分析[J].公路交通科技，2014（6）：1-10+22.

[4]張圣東.高邊坡整體穩定性分析研究[J].黑龍江交通科技，2020（3）：39-41.

[5]袁冬艷，張太樂，彭大江，等.山區路基邊坡穩定性分析及生態防護[J].住宅與房地產，2019（36）：183.

[6]楊敏，方詩圣，張云峰.期望成本最優的土質高邊坡坡角設計研究[J].工程與建設，2016（1）：5-8.

[7]陳海.基于FLAC3D的陡坡路基不均勻沉降控制措施分析[J].山東交通科技，2020（1）：76-79.

[8]李漢鋒.基于FLAC3D的土質路基邊坡動力響應規律研究[J].公路交通科技（應用技術版），2019（7）：87-89.