1、內容提要第4章ANSYS邊坡工程應用實例分析本章首先對邊坡工程進行了概述，然后介紹 了 ANSYS模擬邊坡穩定性分析的步驟，最后用實 例詳細介紹了 ANSYS進行邊坡穩定性分析的全 過程。本章重點號邊坡工程概述* ANSYS邊坡穩定性分析步驟ANSYS邊坡穩定性實例分析本章典型效果圖ANSYS10.0 土木工程應用實例分析4.1 邊坡工程概述4.1.1 邊坡工程邊坡指地殼表部一切具有側向臨空面的地質體，是坡面、坡頂及其下部一定深度坡體的總 稱。坡面與坡頂面下部至坡腳高程的巖體稱為坡體。傾斜的地面稱為斜坡，鐵路、公路建筑施工中，所形成的路堤斜坡稱為路堤邊坡；開挖路 塹所形成的斜坡稱為路塹邊坡；水

2、利、市政或露天煤礦等工程開挖施工所形成的斜坡也稱為邊 坡；這些對應工程就稱為邊坡工程對邊坡工程進行地質分類時，考慮了下述各點。首先，按其物質組成，即按組成邊坡的地層和巖性，可以分為巖質邊坡和土質邊坡(后者包括黃土邊坡、砂土邊坡、土石混合邊坡)。地層和巖性是決定邊坡工程地質特征的基本因素之一，也是研究區域性邊坡穩定問題的主要依據.其次，再按邊坡的結構狀況進行分類。因為在巖性相同的條件下，坡體結構是決定邊坡穩定狀 況的主要因素，它直接關系到邊坡穩定性的評價和處理方法。最后，如果邊坡已經變形，再按 其主要變形形式進行劃分。即邊坡類屬的稱謂順序是：巖性一結構一變形。邊坡工程對國民經濟建設有重要的影響：

3、在鐵路、公路與水利建設中，邊坡修建是不可避 免的，邊坡的穩定性嚴重影響到鐵路、公路與水利工程的施工安全、運營安全以及建設成本。 在路堤施工中，在路堤高度一定條件下，坡角越大，路基所占面積就越小，反之越大。在山區， 坡角越大，則路堤所需填方量越少。因此，很有必要對邊坡穩定性進行分析，4.1.2 邊坡變形破壞基本原理4.1.2.1 應力分布狀態邊坡從其形成開始，就處于各種應力作用 (自重應力、構造應力、熱應力等)之下。在邊坡 的發展變化過程中，由于邊坡形態和結構的不斷改變以及自然和人為營力的作用，邊坡的應力 狀態也隨之調整改變。根據資料及有限元法計算，應力主要發生以下變化：(1)巖體中的主應力跡線

4、發生明顯偏轉，邊坡坡面附近最大主應力方向和坡而平行，而最 小主應力方向則與坡面近于垂直，并開始出現水平方向的剪應力，其總趨勢是由內向 外增多，愈近坡腳愈高，向坡內逐漸恢復到原始應力狀態。(2)在坡腳逐漸形成明顯的應力集中帶。邊坡愈陡，應力集中愈嚴重，最大最小主應力的 差值也愈大。此外，在邊坡下邊分別形成切向應力減弱帶和水平應力緊縮帶，而在靠 近邊坡的表部所測得的應力值均大于按上覆巖體重量計算的數值。(3)邊坡坡面巖體由于側向應力近于零，實際上變為兩向受力。在較陡邊坡的坡面和頂面，出現拉應力，形成拉應力帶.拉應力帶的分布位置與邊坡的形狀和坡面的角度有關。邊坡應力的調整和拉應力帶的出現，是邊坡變形

5、破壞最初始的征兆。例如，由于坡腳應 力的集中，常是坡腳出現擠壓破碎帶的原因；由于坡面及坡頂出現拉應力帶，常是表 層巖體松動變形的原因。4.1.2.2 邊坡巖體變形破壞基本形式邊坡在復雜的內外地質營力作用下形成，又在各種因素作用下變化發展。所有邊坡都在不 斷變形過程中，通過變形逐步發展至破壞。 其基本變形破壞形式主要有：松弛張裂、滑動、崩塌、 傾倒、蠕動和流動。4.1.3 影響邊坡穩定性的因案影響邊坡穩定性的主要因素有：(1)邊坡材料力學特性參數：包括彈性模量、泊松比、摩擦角、粘結力、容重、抗剪強度等參數。(2)邊坡的幾何尺寸參數：包括邊坡高度、坡面角和邊坡邊界尺寸以及坡面后方坡體的幾何形狀，即

6、坡體的不連續面與開挖面的坡度及方向之間的幾何關系，它將確定坡體的各個部分是否滑動或塌落。(3)邊坡外部荷載：包括地震力、重力場、滲流場、地質構造地應力等。4.1.4 邊坡穩定性的分析方法分析邊坡穩定問題，基本上可以分為兩種方法：極限平衡方法和數值分析方法。4.1.4.1 極限平衡方法極限平衡方法的基本思想是：以摩爾一庫侖抗剪強度理論為基礎，將滑坡體劃分成若干垂直 條塊，建立作用在垂直條塊上的力的平衡方程式，求解安全系數。這種計算分析方法遵循下列基本假定 ：(1)遵循庫侖定律或由此引伸的準則。(2)將滑體作為均質剛性體考慮， 認為滑體本身不變形，且可以傳遞應力。因此只研究滑動 面上的受力大小，不

7、研究滑體及滑床內部的應力狀態。(3)將滑體的邊界條件大大簡化。如將復雜的滑體型態簡化為簡單的幾何型態；將滑面簡化為圓弧面、平面或折面；一般將立體問題簡化為平面問題，取沿滑動方向的代表性剖面， 以表征滑體的基本型態；將均布力簡化為集中力， 有時還將力的作用點簡化為通過滑體 重心。極限平衡方法包括以下幾種方法：(1)瑞典圓弧滑動法(2)簡化逼肖普法(3)簡布普通條分法(4)摩根斯坦-普賴斯法(5)不平衡推力傳遞法以上各種方法都是假定土體是理想塑性材料，把土條作為一個剛體，按照極限平衡的原則 進行力的分析，最大的不同之處在于對相鄰上條之間的內力作何種假定，也就是如何增加已知 條件使超靜定問題變成靜定

8、問題。這些假定的物理意義不一樣， 所能滿足的平衡條件也不相同，計算步驟有繁有簡，使用時必須注意他們的適用場合。極限平衡方法關鍵是對滑體的休型和滑面的形態進行分析、正確選用滑面的計算參數以及 正確引用滑體的荷載條件等。因為極限平衡方法完全不考慮土體本身的應力-應變關系，不能真 實地反映邊坡失穩時的應力場和位移場，因此而受到質疑。4.1.4.2 數值分析方法數值數值分析方法考慮土體應力應變關系，克服了極限平衡方法完全不考慮土體本身的應 力-應變關系缺點，為邊坡穩定分析提供了較為正確和深入的概念。邊坡穩定性數值分析方法主要包含以下幾種方法：(1)有限元法有限單元法是數值模擬方法在邊坡穩定評價中應用最

9、早的方法，也是目前最廣泛使用的一 種數值方法，可以用來求解彈性、彈塑性、粘彈塑性、粘塑性等問題。目前用有限元法求解邊坡穩定主要有兩種方法。a.有限元滑面搜索法：將邊坡體離散為有限單元格，按照施加的荷載及邊界條件進行有限元 計算可得到每個結點的應力張量。然后假定一個滑動面，用有限元數據給出滑動面任一點的向 正應力和剪應力，根據摩爾一庫侖準則可得該點的抗滑力，由此即能求得滑動面上每個結點的 下滑力與抗滑力，再對滑動面上下滑力與抗滑力進行積分，就可以求得每一個滑動面的安全系 數。b.有限元強度折減法：首先選取初始折減系數，將巖土體強度參數進行折減， 將折減后的參 數作為輸入，進行有限元計算，若程序收

10、斂，則巖土體仍處于穩定狀態，然后再增加折減系數，直到程序恰好不收斂，此時的折減系數即為穩定或安全系數。(2)自適應有限元法自20世紀70年代開始自適應理論被引入有限元計算，主導思想是減少前處理工作量和實現 網格離散的客觀控制。現已基本建立了一般彈性力學、流體動力學、滲流分析等領域的平面自 適應分析系統，能使計算較為快速和準確。(3)離散單元法離散單元法的突出功能是它在反映巖塊之間接觸面的滑移、分離與傾翻等大位移的同時，又能計算巖塊內部的變形與應力分布。因此，任何一種巖體材料都可引入到模型中，例如彈性、粘彈性或斷裂等均可考慮，故該法對塊狀結構、層狀破裂或一般破裂結構巖體邊坡比較合適。并且，它利用

11、顯式時間差分法(動態差分法)求解動力平衡方程，求解非線性大位移與動力問題 比較容易。離散元法在模擬過程中考慮了邊坡失穩破壞的動態過程，允許巖土體存在滑動、平移、轉 動和巖體的斷裂及松散等復雜過程，具有宏觀上的不連續性和單個巖塊休運動的隨機性，可以 較真實、動態地模擬邊坡在形成和開挖過程中應力、位移和狀態的變化，預測邊坡的穩定性， 因此在巖質高邊坡穩定性的研究中得到廣泛的應用。(4)拉格朗日元法為了克服有限元等方法不能求解大變形問題的缺陷，人們根據有限差分法的原理，提出了 FLA徵值分析方法。該方法較有限元法能更好地考慮巖土體的不連續和大變形特性，求解速度 較快。缺點是計算邊界、單元網格的劃分帶

12、有很大的隨意性。(5)界面元法界面元法是一種基于累積單元變形于界面的界面應力元法模型，建立適用于分析不連續、 非均勻、各向異性和各類非線性問題、場問題，以及能夠完全模擬各類錨桿復雜空間布局和開 挖擾動的方法。4.2 .4.3有限元法用于邊坡穩定性分析優點有限元法考慮了介質的變形特征，真實地反應了邊坡的受力狀態。它可以模擬連續介質， 也可以模擬不連續介質；能考慮邊坡沿軟弱結構面的破壞，也能分析邊坡的整體穩定破壞。有 限元法可以模擬邊坡的圓弧滑動破壞和非圓弧滑動破壞。同時它還能適應各種邊界條件和不規 則幾何形狀，具有很廣泛的適用性。有限元法應用于邊坡工程，有其獨特的優越性。與一般解析方法相比，有限

13、元法有以下優點：(1)它考慮了巖體的應力-應變關系，求出每一單元的應力與變形，反映了巖體真實工作狀O(2)與極限平衡法相比，不需要進行條間力的簡化，巖體自始至終處于平衡狀態。(3)不需要像極限平衡法一樣事先假定邊坡的滑動面，邊坡的變形特性、塑性區形成都根據實際應力應變狀態“自然”形成。(4)若巖體的初始應力己知，可以模擬有構造應力邊坡的受力狀態。(5)不但能像極限平衡法一樣模擬邊坡的整體破壞，還能模擬邊坡的局部破壞，把邊坡的整體破壞和局部破壞納入統一的體系。(6)可以模擬邊坡的開挖過程，描述和反應巖體中存在的節理裂隙、斷層等構造面。鑒于有限元法具有如此多優點，本章借助通用有限元軟件ANSY來實

14、現對邊坡穩定性分析，用具體的邊坡工程實例詳細介紹應用 ANSY軟件分析邊坡穩定性問題。4.2 ANSYS邊坡穩定性分析步驟ANSYS邊坡穩定性分析一般分以下幾個步驟：1.1.1 理環境1.1.2 型，劃分網格，對模型的不同區域賦予特性1.1.3 條件和載荷1.1.41.1.5 （查看計算結果）4.2.1 創建物理環境在定義邊坡穩定性分析問題的物理環境時，進入ANSYS前處理器，建立這個邊坡穩定性分析的數學仿真模型。按照以下幾個步驟來建立物理環境：1、設置GUT菜單過濾如果你希望通過 GUI路徑來運行ANSYS ,當ANSYS被激活后第一件要做的事情就是選 擇菜單路徑：Main Menu>

15、Preferences,執行上述命令后，彈出一個如圖4-1所示的對話框出現后，選擇Structural。這樣ANSYS會根據你所選擇的參數來對GUI圖形界面進行過濾，選擇Structural以便在進行邊坡穩定性分析時過濾掉一些不必要的菜單及相應圖形界面。2、定義分析標題（/ TITLE ）在進行分析前，可以給你所要進行的分析起一個能夠代表所分析內容的標題，比如“ Slope stability Analysis”,以便能夠從標題上與其他相似物理幾何模型區別。用下列方法定義分析標題。命令：/TITLEGUI : Utility Menu>File>Change Title3、說明單元

16、類型及其選項（KEYOPT選項）與ANSYS的其他分析一樣，也要進行相應白單元選擇。 ANSYS軟件提供了 100種以上的 單元類型，可以用來模擬工程中的各種結構和材料，各種不同的單元組合在一起，成為具體的 物理問題的抽象模型。例如，不同材料屬性的邊坡土體用PLANE82單元來模擬。大多數單元類型都有關鍵選項 （KEYOPTS ）,這些選項用以修正單元特性。 例如，PLANE82 有如下KEYOPTS :KEYOPT（2）包含或抑制過大位移設置KEYOPT（3）平面應力、軸對稱、平面應變或考慮厚度的平面應力設置KEYOPT（9）用戶子程序初始應力設置設置單元以及其關鍵選項的方式如下：命令：ET

17、KEYOPTGUI : Main Menu> PreprocessorElement Type> Add/Edit/Delete圖4-1 GUI圖形界面過濾4定義單位結構分析只有時間單位、長度單位和質量單位三個基本單位，則所有輸入的數據都應當是這三個單位組成的表達方式。如標準國際單位制下，時間是秒（s）,長度是米（m）,質量是千克（kg）,則導出力的單位是kg? m/s2 （相當于牛頓N）,材料的彈性模量單位是 kg/m? s2 （相 當于帕Pa）。命令：/UNITS5、定義材料屬性大多數單元類型在進行程序分析時都需要指定材料特性，ANSYS程序可方便地定義各種材料的特性，如結構材

18、料屬性參數、熱性能參數、流體性能參數和電磁性能參數等。ANSYS程序可定義的材料特性有以下三種：（1）線性或非線性。（2）各向同性、正交異性或非彈性。（3）隨溫度變化或不隨溫度變化。因為分析的邊坡模型采用理想彈塑性模型（D-P模型），因此邊坡穩定性分析中需要定義邊坡中不同土體的材料屬性：容重、彈性模量、泊松比、凝聚力以及摩擦角。命令：MPGUI: Main Menu>Preprocessor>Material Props>Material Models 或 Main Menu>Solution>Load Step Opts>Other>Change M

19、at Props>Material Models進行邊坡穩定性分析計算時，采用強度折減法來實現。首先選取初始折減系數F,然后對邊坡土體材料強度系數進行折減，折減后凝聚力以及摩擦角分別見式4-1和式4-2。tanFtan(4-1)(4-2)C和為邊坡土體的初始凝聚力和摩擦角。C和 進行折減，輸入邊坡模型計算，若收斂，則此時邊坡是穩定的；繼續增大折減系數F,直到程序恰好不收斂，此時的折減系數即為穩定或安全系數。4.2.2 建立模型和劃分網格創建好物理環境，就可以建立模型。在進行邊坡穩定性分析時，需要建立模擬邊坡土體的 PLANE82單元。在建立好的模型各個區域內指定特性（單元類型、選項、實常

20、數和材料性質等）以后，就可以劃分有限元網格了。通過GUI為模型中的各區賦予特性：1、選擇 Main Menu> PreprocessorMeshing> Mesh AttributesPicked Areas2、點擊模型中要選定的區域。3、在對話框中為所選定的區域說明材料號、實常數號、單元類型號和單元坐標系號。4、重復以上三個步驟，直至處理完所有區域。通過命令為模型中的各區賦予特性：ASEL （選擇模型區域）MAT （說明材料號）REAL （說明實常數組號）TYPE （指定單元類型號）ESYS （說明單元坐標系號）4.2.3 施加約束和荷載在施加邊界條件和荷載時，既可以給實體模型（

21、關鍵點、線、面）也可以給有限元模型（節點和單元）施加邊界條件和荷載。在求解時，ANSYS程序會自動將加到實體模型上的邊界條件 和載荷轉遞到有限元模型上。邊坡穩定性分析中，主要是給邊坡兩側和底部施加自由度約束。命令：D施加荷載包括自重荷載以及邊坡開挖荷載。4.2.4 求解接著就可以進行求解，ANSYS程序根據現有選項的設置， 從數據庫獲取模型和載荷信息并 進行計算求解，將結果數據寫入到結果文件和數據庫中。命令：SOLVEGUI : Main Menu>Solution>Solve> Current LS4.2.5 后處理后處理的目的是以圖和表的形式描述計算結果。對于邊坡穩定性分

22、析中， 進入后處理器后,查看邊坡變形圖和節點的位移、應力和應變。隨著強度折減系數的增大，邊坡的水平位移增大, 塑性應變急劇發展，塑性區發展形成一個貫通區域時，計算不收斂，認為邊坡發生了破壞。通 過研究位移、應變和塑性區域，來綜合判斷邊坡的穩定性。命令：/POST1GUI: Main Menu> General Postproc4.2.6 補充說明邊坡的失穩破壞定義有很多種，對于采用彈塑性計算模型的邊坡，需要綜合考慮以下因素(1)把有限元計算的收斂與否作為一個重要的衡量指標，邊坡處于穩定狀態，計算收斂， 邊坡破壞時，邊坡不收斂。(2)邊坡失穩的同時還表現出位移急劇增加。(3)邊坡失穩總是伴

23、隨著塑性變形的明顯增加和塑性區的發展，塑性區的發展狀況反映了邊坡是否處于穩定狀態。此外，采用彈塑性有限元法進行計算，它具有獨特的優勢：(1)彈塑性分析假定巖體為彈塑性材料，巖體在受力初期處于彈性狀態，達到一定的屈服準 則后，處于塑性狀態。采用彈塑性模型更能反應巖體的實際工作狀態。(2)巖體所承受的荷載超過材料強度時，就會出現明顯的滑移破壞面。因此，彈塑性計算不 需要假定破壞面的形狀和位置，破壞面根據剪應力強度理論自動形成。當整個邊坡破壞時，就 會出現明顯的塑性區。(3 )能綜合考慮邊坡的局部失穩和整體失穩破壞。1.1.1.1 ANSYS邊坡穩定性實例分析4.3.1 實例描述米Jr米| 米刃米米

24、邊坡圍巖2彈塑性材料,米邊坡圍巖1彈性材料圖4-2邊坡模型邊坡實例選取國內某礦，該邊坡考慮彈性和塑性兩種材料，邊坡尺寸如圖4-2所示。分析目的是對該邊坡進行穩定性計算分析，以判斷其穩定性和計算出安全系數，該邊坡圍巖材料屬性 見表4-1。表4-1邊坡模型圍巖參數類別彈性模量/Gpa泊松比v容重/ kN / m3內聚力 /Mpa摩擦角（。）圍巖2 （彈塑性）300.2525000.942圍巖1 （彈性）310.242700-對于像邊坡這樣縱向很長的實體，計算模型可以簡化為平面應變問題。假定邊坡所承受的外力不隨Z軸變化，位移和應變都發生在自身平面內。對于邊坡變形和穩定性分析， 這種平面假設是合理的。

25、實測經驗表明，邊坡的影響范圍在2倍坡高范圍，因此本文計算區域為邊坡體橫向 延伸2倍坡高，縱向延伸3倍坡高。兩側邊界水平位移為零，下側邊界豎向位移為零。彈性有限 元的計算模型如圖4-2所示。采用雙層模型，模型上部為理想彈塑性材料，下部為彈性材料，左右邊界水平位移為 零，下邊界豎向位移為零。雙層模型考慮土體的彈塑性變形，其塑性區的發展，應力的分布更符合實際情況。考慮雙層模型，塑性區下部的單元可以產生一定的垂直變形和水平變形，基本消除了 由于邊界效應在邊坡下部出現的塑性區，更好地模擬了邊坡的變形和塑性區的發展。4.3.2 GUI操作方法1.1.1.2 創建物理環境1）在【開始】菜單中依次選取【所有程

26、序】/ ANSYS10.0 /【ANSYS Product Launcher,得到“10.0ANSYS Product Launcher”對話框。2）選中File Management】，在"Working Directory ”欄輸入工作目錄 "D:ansysexample4-1”, 在"Job Name”欄輸入文件名“Slope”。3）單擊“ RUN”按鈕，進入 ANSYS10.0的GUI操作界面。4）過濾圖形界面：Main Menu> Preferences,彈出"Preferences for GUI Filtering”對話框，選中“St

27、ructural”來對后面的分析進行菜單及相應的圖形界面過濾。5）定義工作標題：Utility Menu> File> Change Title ,在彈出的對話框中輸入"Slope stabilityAnalysis”,單擊“OK”,如圖 4-3。圖4-3定義工作標題6）定義單元類型：a.定義 PLANE82 單元：Main Menu> PreprocessorElement Type> Add/Edit/Delete ,彈出一 個單元類型對話框，單擊“ Add”按鈕。彈出如圖4-4所示對話框。在該對話框左面滾動欄中選 擇“Solid”，在右邊的滾動欄中選擇“

28、 Quad 8node 82',單擊“Apply"，就定義了 “ PLANE82 ” 單元。圖4-4定義PLANE82單元對話框b.設定 PLANE82 單元選項：Main Menu> PreprocessorElement Type> Add/Edit/Delete ,彈 出一個單元類型對話框， 選中“Type 2 PLANE82 "，單擊“Options”按鈕，彈出一個"PLANE82 element Type options”對話框，如圖 4-5所示。在"Element behavior K3 ”欄后面的下拉菜單中 選取“ Pl

29、ane strain”，其它欄后面的下拉菜單采用ANSYS默認設置就可以，單擊“ OK”按鈕。PLANE82 element type 0Pti0rlsOptions fcr PL4NE日2, Element Type Ret No. 1Element behaviorlane strainExtr曰 element outputExtra surface outputOK| No outputK6CancelNo extra output圖4-5 PLANE82單元庫類型選項對話框通過設置PLANE82單元選項“ K3”為“Plane strain”來設定本實例分析采取平面應 變模型進行分析

30、。因為邊坡是縱向很長的實體，故計算模型可以簡化為平面應變問題。8節點PLANE8邪元每個節點有 UX和UY兩個自由度，比4節點PLANE4卻元具有更 高的精確性，對不規則網格適應性更強。7）定義材料屬性a.定義邊坡圍巖 1 材料屬性：Main MenuPreprocessorMaterial Props> Material Models ,彈 出“Define Material Model Behavior ”對話框，如圖 4-6所示。圖4-6定義材料本構模型對話框在圖4-6中右邊欄中連續雙擊"StructuralLinearElastic>Isotropic”后，又彈出

31、如圖 4-7 所示 a Linear Isotropic Properties for Material Number 1 ”對話框，在該對話框中"EX"后面的輸 入欄輸入“ 3E10”，在“PRXY”后面的輸入欄輸入“ 0.25”，單擊“OK”。再在選中“ Density” 并雙擊，彈出如圖4-8所示“Density for Material Number 1 ”對話框，在"DENS"后面的欄中輸 入邊坡土體材料的密度“ 2500”，單擊“ OK”按鈕。再次在圖 4-6 中右邊的欄中連續雙擊“Structural>Nonlinear> In

32、elastic>Non-metalplasticity>drucker-prager”后，又彈出一個如圖4-9所示對話框。在"Cohesion”欄添入邊坡圍巖材料1的內聚力“ 0.9E6”，在“ Fric Angle”欄添入邊坡內摩擦角“ 42”，單擊“ OK”按鈕。Ri imp/ RiiTnpicPRRerUes for MafterlBl FW Density for Material Number 1Lreor LsDincpc Material 4叩crtim for Materid Ibmtcr 1Density fcr Macer al Vimbe* 1EXF

33、T1nTErrperatLreECENS feODAdd 'epeiire behteCncel I-卸除.Het圖4-7線彈性材料模型對話框圖4-8材料密度輸入對話框圖4-9定義邊坡材料1DP模型對話框b.定義邊坡圍巖2材料屬性：在圖4-6對話框中，單擊“ MaterialNew Model”，彈出一 個“Define Material ID ”對話框，在“ID”欄后面輸入材料編號“ 2”，單擊“ OK”按鈕。彈出 一個定義材料模型對話框對話框，選中" Material Model Number 2 ",和定義邊坡圍巖 1材料一 樣，在右邊的欄中連續雙擊 “Str

34、ucturalLinearElastic>Isotropic ”后，又彈出一個"Linea門sotropic Properties for Material Number 2”對話框，在該對話框中“ EX”后面的輸入欄輸入“ 3.2E10”， 在“PRXY”后面的輸入欄輸入“0.24”，單擊“OK”。再選中“Density”并雙擊，彈出一個"Density for Material Number 2”對話框，在“ DENS”后面的欄中輸入隧道圍巖材料的密度“2700”，再單擊“OK”按鈕，彈出一個定義材料模型對話框。c.復制邊坡圍巖1材料性質：在圖4-6對話框中，用鼠

35、標點擊“Edit>copy .”,彈出一個“Copy Material Model”對話框，如圖4-10所示。在"from Material number”欄后面的下拉菜單中選取 "1",在"TO Material number”欄后面輸入“3",單擊"Apply”按鈕。又彈出如土 4-10所示 對話框，然后依次在 “TO Material number”欄后面輸入“4”、“5”、“6”、“7”、“8”“9”、“10”，“11”、“12”、“13”，每輸入一個數，就單擊“ Apply”按鈕一次。圖4-10復制本構模型對話框最后彳

36、#到10個復制圍巖1的邊坡材料本構模型，如圖4-11所示。圖4-11定義強度折減后材料模型對話框Drucker Prager for Material Number 3Drucker Prager Table for Material Number 3Cohesion Frit Angle Flow Angle37,7GraphOK | Sncel | Hdp |圖4-12定義強度折減系數F=1.2時圍巖材料對話框d.定義10個強度折減后材料本構模型：首先定義強度折減系數 F=1.2后邊坡圍巖材料模型， 在圖4-11對話框中，在鼠標依次雙擊“ Material Model Number 3/D

37、rucper-Prager"。彈出一個“Drucker- Prager Material Number 3",如圖 4-12 所示，在"Cohesion”欄添入強度折減系數 F=1.2 后邊坡圍巖材料1的內聚力“ 0.75E6”，在“Fric Angle”欄添入折減后邊坡內摩擦角“ 37.7”， 單擊“ OK”按鈕。用相同方法定義強度折減系數分別為：F=1.4、F=1.6、F=1.8、F=2.0、F=2.2、F=2.4、F=2.6、F=2.8、F=3.0的邊坡圍巖材料本構模型。定義強度折減后本構模型目的是為了分析邊坡穩定性。強度折減就是降低內聚力和摩擦角，根據式

38、4-1和式4-2進行折減。1.1.1.3 建立模型和劃分網格1）創建邊坡線模型a.輸入關鍵點：Main Menu> PreprocessorModelingCreate>Keypoints>In Active CS ,彈出“Creae Keypoints in Active Cooedinate System”對話框，如圖 4-13所示。在"NPT keypoint number” 欄后面輸入 “1”，在 "X, Y, Z Location in active CS ” 欄后面輸入 “(0, 0, 0)”，單擊 “Apply” 按鈕，這樣就創建了關鍵點1。

39、再依次重復在“NPT keypoint number”欄后面輸入“2、3、4、5、6、7、8、9”，在對應 “X, Y, Z Location in active CS” 欄后面輸入 “(-800, 0, 0)、(-800, -800, 0)、(-800, -1200, 0)、(1200, -1200, 0)、(1200, -800, 0)、( 1200, 0, 0)、(1200, 378, 0)、(430, 378, 0),最后單擊 “ OK” 按鈕。圖4-13在當前坐標系創建關鍵點對話框b.創建邊坡線模型：Main Menu> PreprocessorModelingCreate&g

40、t;Lines>Straight line ,彈出“Creae straight lines'對話框，用鼠標依次點擊關鍵點1、2,單擊“Apply”按鈕，這樣就創建了直接 L1,同樣分別連接關鍵點“2、3”，“3、4”，“4、5”，“5、6”，“6、7”，“7、4”，“7、8”,8、9”，“9、1”，“9、2”，最后單擊“ OK”按鈕，就得到邊坡線模型，如圖 4-14所示。圖4-14 邊坡線模型3）創建邊坡面模型a.打開面編號顯示： Utility Menu> PlotCtrlsNumbering ,彈出"Plot Numbering Controls” 對話框，

41、如圖4-15所示。選中"Aares Numbers”選項，后面白文字由"of"變為"on”,單擊“OK”關閉窗口。圖4-15打開面編號對話框b.創建邊坡面模型：Main Menu> PreprocessorModelingCreate>Areas> Arbitraryby line,彈出一個"Create Area by lines”對話框,在圖形中選取線L4、L5、L3和L11,點擊“Apply按鈕，就生成了邊坡彈性材料區域面積A1 ;再依次用鼠標在圖形中選取線L1、L2、L6、L10、和L11,點擊“Apply ”按鈕，就

42、生成了邊坡塑性材料區域面積A2;再依次用鼠標在圖形中選取線L7、L8、L10和L9,點擊“ OK ”按鈕，就生成了邊坡開挖掉區域面積A3。最后得到邊坡模型的面模型，如圖4-16所示。4）劃分邊坡圍巖2單元網格a.給邊坡圍巖 2 賦予材料特性：Main Menu> PreprocessorMeshing> MeshTool,彈出“MeshTool”對話框，如圖4-17所示。在"Element Attributes”后面的下拉式選擇欄中選擇 "Areas", 按"Set"按鈕，彈出一個"Areas Attributes”面拾取

43、框，在圖形界面上拾取邊坡圍巖2區域，單擊拾取框上的"OK"按鈕，又彈出一個如圖4-18所示的“Areas Attributes”對話框，在"Material number"后面的下拉式選擇欄中選取"2",在"Element type number " 后面的下拉式選擇欄 中選取 “ 2 PLANE82 "，單擊 “ Apply ”。圖4-17網格劃分工具欄圖4-18定義單元屬性對話框b.設置網格劃分份數：在圖 4-17工具欄中“ Size Control”欄，用鼠標點擊“lines”后面的 "S

44、et",彈出一個選擇對話框，用鼠標在圖形選擇線 L3和L5。彈出一個"Element Sizes on PickedLines",對話框，如圖4-19所示，在"No of element division”欄后面輸入"5",單擊"Apply" 按鈕，再選才i線L4和L11,又彈出圖4-19對話框，在"No of element division ”欄后面輸入“26”, 單擊“OK”按鈕。c.劃分單元網格：在圖4-17網格劃分工具欄中單擊“ Mesh”按鈕，彈出一個拾取面積對話 框，拾取面積A1,單擊拾取框

45、上的“ OK”按鈕，生成邊坡圍巖 2單元網格。圖4-19設置網格份數對話框5）劃分邊坡圍巖1單元網格Hot Lever 8門troIsk廣 UjtpAuh廣 EonL Polyann 廣 ca.teleL Ln。2C umIf ft¥ZVl1TI¥l=210 1HiiiiwAAU.JLQin. Nu-r*5 Li«t a t 工 ft. 他L Hiey Hn芫,Tixe圖4-20選取線對話框圖4-21設置網格分數對話框a.設置網格份數：Main Menu>Preprocessor>Meshing>Size Cntrls>ManualSize

46、>Layers>PickedLines,彈出一個"Set Layer Controls”對話框，如圖4-20所示，用鼠標選取線 L1、L1和L6單 擊"OK”按鈕。彈出一個"Area Layer Mesh Control on Picked lines”對話框，如圖 4-21 所示， 在"No of line division ”欄后面輸入"10",單擊"OK"按鈕。相同方法設置線L8和L10分割份數為16;設置線L7和L9線的分割份數為12。b.給邊坡圍巖 1 賦予材料特性：Main Menu>

47、 PreprocessorMeshing> MeshTool,彈出“MeshTool”對話框，如圖4-17所示。在"Element Attributes”后面的下拉式選擇欄中選擇 "Areas", 按“Set”按鈕，彈出一個" Areas Attributes "面拾取框，在圖形界面上拾取面A2和A3 ,單擊拾取框上的"OK"按鈕，又彈出一個 "Areas Attributes"對話框，在"Material number"后面 的下拉式選擇欄中選取"1",在&

48、quot;Element type number " 后面的下拉式選擇欄中選取“2PLANE82"，單擊 “Apply”。c.劃分單元網格：在圖3-97網格劃分工具欄中單擊“ Mesh”按鈕，彈出一個拾取面積對話框，拾取圍巖，單擊拾取框上的“ OK”按鈕，生成邊坡圍巖1單元網格。最后得到邊坡模型單元網格，如圖 4-22所示。圖4-22邊坡模型單元網格4.3.2.3施加約束和荷載1）給邊坡模型施加約束a.給邊坡模型兩邊施加約束：執行Main Menu>Solution>Define Loads>Apply>Structural>Displacem

49、ent>on Nodes,彈出在節點上施加位移約束對話框，用鼠標選取隧道模型兩側邊界上 所有節點，單擊“OK”按鈕。彈出“Apply U,ROT on Nodes ”對話框，如圖4-23所示，在“DOFS to be constrained'欄后面中選取"UX",在"Apply as”欄后面的下拉菜單中選取 "Constant value", 在"Displacement value”欄后面輸入“0"值，然后單擊"OK"按鈕。圖4-23給模型兩側施加位移約束對話框圖4-24給模型底部施加位移

50、約束對話框b.給模型底部施加約束：執行 Main Menu>Solution>Define LoadsApply>Structural> Displacement>on Nodes,彈出在節點上施加位移約束對話框，用鼠標選取隧道模型底部邊界上 所有節點，單擊"OK”按鈕。彈出圖4-24所示對話框，在"DOFS to be constrained”欄后面中 選取"UX、UY",在"Apply as”欄后面的下拉菜單中選取"Constant value”,在 a Displacement value”欄后面輸

51、入“ 0”值，然后單擊“ OK”按鈕。節點選擇，可以先選擇節點上線，再選擇附在線上的節點。2)施加重力加速度： Main Menu>Solution>Define Loads>Apply>Structural>Inertia>Gravity , 彈出"Apply(Gravitational)Acceleration ”對話框，如圖 4-25 所示。只需在“Global Cartesian Y-comp 欄后面輸入重力加速度值“ 9.8”就可以，單擊“ OK”按鈕，就完成了重力加速度的施加。A/Xpplv (G日¥it目忖。nal) Acc

52、elerationACEL Apply (Gnanitathxal)心匚匚cgraticinACELX Global Cartesian X-conpAC ELY Global Cartesian Y-coimpg,日ACEL G lobal Cartes ian Z-cornp口口長 |Canizel |H.lp |圖4-25施加重力加速度對話框這時就可以得到施加約束和重力加速度后隧道有限元模型，如圖4-25所示。圖4-25施加約束和重力荷載后邊坡模型4.3.2.4 求解1）求解設置a.指定求解類型： Main Menu>Solution>Analysis Type>New

53、 Analysis ,彈出一個如圖 4-25 所 示對話框，在"Type of analysis”欄后面選中"Static",單擊"OK"按鈕。MNsw ArkalysisMITYPF Typp of wn曲擅Stade ；C fVbcUIr MarmcnicC TransientU SfiectruniC Eigen Buckling SdbsvuctLirng<LMSCK ICanoe I IHelp I圖4-25指定求解類型對話框b.設置載荷步：Main Menu>Preprocessor>Loads>Analys

54、is Type>Sol'n Controls,彈出一個Solution Controls”對話框，用鼠標單擊"Basic”選項，如圖 4-26 所示，在"Number of Substeps”欄后面輸入"5",在"Max no. of substeps”欄后面輸入“100",在"Min no. of substeps”欄后 面輸入“1”，單擊“ OK”按鈕。圖4-26設置載荷步對話框c.設置線性搜索： Main Menu>Solution>Analysis Type>Sol'n Con

55、trols,彈出一個"Solution Controls”對話框，用鼠標單擊"Nonlinear”選項，如圖4-27所示，在"Line search"欄后面下 拉菜單選中“ ON圖4-27設置線性搜索對話框d.設定牛頓-拉普森選項：Main Menu>Solution>Analysis Type>Analysis Options,彈出一個 “Static or Steady-static Analysis”對話框，如圖 4-28 所示，在"New-Raphson option”欄后面下 拉菜單選中“ Full N-R”，單擊“

56、 OK”按鈕。'' i'"k c； .il-Ai 11 rr'>isR 繇髓鬻Mmirear 卬口 ms忖L里DIM .a-qe defa-m erf=ct=off一W0PT1 wflor-Rsftscri opricn-MsutrvEdKsM回7小生,二Lirwr Options：LJvTTvIl 3白 Limped HliCE 即口MXyr NnEQELV Equattir scUver|prcgnm dTjaenTolgrarciiAEl-xald 怕。all lEwcept Fro-ital and 5pa*se § "

57、;hjefs版i 中 liar-iion hf for Fi esnciU口n CGFRECISICfJ' S fiylt R'wl teiCri *off-raid an 勺 rir FrEciniitinri CGMSAE lYemorj Save -廠 off-問 rl on ly mr Freconcflinn 06nya-ECKPimtiCh«ck-vaidory r?r Frortai u方汨 aE 匚。u匚"白令rTTlpsiRF1-.匚;*K1 *e而吟鼻力NJtu; Ifr-ICEOM.CT-J hjii &JL BSTtF OM7

58、OFFGT TeiYperaiufe ctffirencfl-csixal-tfiMBI JlXQ Ute Hr 口 一 d : CI J Lfx：u ISfflp SCdU3圖4-28設定牛頓-拉普森選項e.打開大位移求解：Main Menu>Preprocessor>Loads>Analysis Type>Sol'n Controls>Basic ,得 到圖4-26所示對話框，在"Analysis Options”欄后面下拉菜單選中大位移"Large Displacement static”，單擊 “ OK” 按鈕。f.設置收斂條件

59、： Main Menu>Preprocessor>Loads>Load Step Opts>Nonlinear>Convergence Crit,彈出一個"Default Nonlinear Convergence Criteria”,如圖 4-29 所示。圖中顯示 ANSYS 默 認的收斂條件：分別設置了力和力矩的收斂條件。圖4-29 ANSYS默認收斂條件為了使求解順利進行和得到較好解，可以修改默認收斂設置，可以分別設置力、力矩和位 移收斂條件。單擊圖 4-29中的“Replace",彈出一個"Nonlinear Convergence Criteria”對話框， 如圖4