ANSYS邊坡工程應(yīng)用實(shí)例分析

1、第4章 ANSYS邊坡工程應(yīng)用實(shí)例分析本章重點(diǎn) 邊坡工程概述 ANSYS邊坡穩(wěn)定性分析步驟 ANSYS邊坡穩(wěn)定性實(shí)例分析本章典型效果圖4.1 邊坡工程概述 邊坡工程邊坡指地殼表部一切具有側(cè)向臨空面的地質(zhì)體，是坡面、坡頂及其下部一定深度坡體的總稱。坡面與坡頂面下部至坡腳高程的巖體稱為坡體。傾斜的地面稱為斜坡，鐵路、公路建筑施工中，所形成的路堤斜坡稱為路堤邊坡；開挖路塹所形成的斜坡稱為路塹邊坡；水利、市政或露天煤礦等工程開挖施工所形成的斜坡也稱為邊坡；這些對(duì)應(yīng)工程就稱為邊坡工程對(duì)邊坡工程進(jìn)行地質(zhì)分類時(shí)，考慮了下述各點(diǎn)。首先，按其物質(zhì)組成，即按組成邊坡的地層和巖性，可以分為巖質(zhì)邊坡和土質(zhì)邊坡(后者包

2、括黃土邊坡、砂土邊坡、土石混合邊坡)。地層和巖性是決定邊坡工程地質(zhì)特征的基本因素之一，也是研究區(qū)域性邊坡穩(wěn)定問題的主要依據(jù).其次，再按邊坡的結(jié)構(gòu)狀況進(jìn)行分類。因?yàn)樵趲r性相同的條件下，坡體結(jié)構(gòu)是決定邊坡穩(wěn)定狀況的主要因素，它直接關(guān)系到邊坡穩(wěn)定性的評(píng)價(jià)和處理方法。最后，如果邊坡已經(jīng)變形，再按其主要變形形式進(jìn)行劃分。即邊坡類屬的稱謂順序是:巖性 結(jié)構(gòu)變形。邊坡工程對(duì)國(guó)民經(jīng)濟(jì)建設(shè)有重要的影響：在鐵路、公路與水利建設(shè)中，邊坡修建是不可避免的，邊坡的穩(wěn)定性嚴(yán)重影響到鐵路、公路與水利工程的施工安全、運(yùn)營(yíng)安全以及建設(shè)成本。在路堤施工中，在路堤高度一定條件下，坡角越大，路基所占面積就越小，反之越大。在山區(qū)，坡角

3、越大，則路堤所需填方量越少。因此，很有必要對(duì)邊坡穩(wěn)定性進(jìn)行分析，4.1.2 邊坡變形破壞基本原理 應(yīng)力分布狀態(tài)邊坡從其形成開始，就處于各種應(yīng)力作用(自重應(yīng)力、構(gòu)造應(yīng)力、熱應(yīng)力等)之下。在邊坡的發(fā)展變化過程中，由于邊坡形態(tài)和結(jié)構(gòu)的不斷改變以及自然和人為營(yíng)力的作用，邊坡的應(yīng)力狀態(tài)也隨之調(diào)整改變。根據(jù)資料及有限元法計(jì)算，應(yīng)力主要發(fā)生以下變化：（1)巖體中的主應(yīng)力跡線發(fā)生明顯偏轉(zhuǎn)，邊坡坡面附近最大主應(yīng)力方向和坡而平行，而最小主應(yīng)力方向則與坡面近于垂直，并開始出現(xiàn)水平方向的剪應(yīng)力，其總趨勢(shì)是由內(nèi)向外增多，愈近坡腳愈高，向坡內(nèi)逐漸恢復(fù)到原始應(yīng)力狀態(tài)。(2) 在坡腳逐漸形成明顯的應(yīng)力集中帶。邊坡愈陡，應(yīng)力集

4、中愈嚴(yán)重，最大最小主應(yīng)力的差值也愈大。此外，在邊坡下邊分別形成切向應(yīng)力減弱帶和水平應(yīng)力緊縮帶，而在靠近邊坡的表部所測(cè)得的應(yīng)力值均大于按上覆巖體重量計(jì)算的數(shù)值。(3) 邊坡坡面巖體由于側(cè)向應(yīng)力近于零，實(shí)際上變?yōu)閮上蚴芰?。在較陡邊坡的坡面和頂面，出現(xiàn)拉應(yīng)力，形成拉應(yīng)力帶.拉應(yīng)力帶的分布位置與邊坡的形狀和坡面的角度有關(guān)。邊坡 應(yīng)力的調(diào)整和拉應(yīng)力帶的出現(xiàn)，是邊坡變形破壞最初始的征兆。例如，由于坡腳應(yīng)力的集中，常是坡腳出現(xiàn)擠壓破碎帶的原因；由于坡面及坡頂出現(xiàn)拉應(yīng)力帶，常是表層巖體松動(dòng)變形的原因。 邊坡巖體變形破壞基本形式邊坡在復(fù)雜的內(nèi)外地質(zhì)營(yíng)力作用下形成，又在各種因素作用下變化發(fā)展。所有邊坡都在不斷變形

5、過程中，通過變形逐步發(fā)展至破壞。其基本變形破壞形式主要有:松弛張裂、滑動(dòng)、崩塌、傾倒、蠕動(dòng)和流動(dòng)。4.1.3 影響邊坡穩(wěn)定性的因案影響邊坡穩(wěn)定性的主要因素有:(1)邊坡材料力學(xué)特性參數(shù)：包括彈性模量、泊松比、摩擦角、粘結(jié)力、容重、抗剪強(qiáng)度等參數(shù)。(2)邊坡的幾何尺寸參數(shù)：包括邊坡高度、坡面角和邊坡邊界尺寸以及坡面后方坡體的幾何形狀，即坡體的不連續(xù)面與開挖面的坡度及方向之間的幾何關(guān)系，它將確定坡體的各個(gè)部分是否滑動(dòng)或塌落。(3)邊坡外部荷載： 包括地震力、重力場(chǎng)、滲流場(chǎng)、地質(zhì)構(gòu)造地應(yīng)力等。4.1.4 邊坡穩(wěn)定性的分析方法分析邊坡穩(wěn)定問題，基本上可以分為兩種方法：極限平衡方法和數(shù)值分析方法。 極限

6、平衡方法極限平衡方法的基本思想是:以摩爾一庫侖抗剪強(qiáng)度理論為基礎(chǔ)，將滑坡體劃分成若干垂直條塊，建立作用在垂直條塊上的力的平衡方程式，求解安全系數(shù)。這種計(jì)算分析方法遵循下列基本假定: (1)遵循庫侖定律或由此引伸的準(zhǔn)則。(2)將滑體作為均質(zhì)剛性體考慮，認(rèn)為滑體本身不變形，且可以傳遞應(yīng)力。因此只研究滑動(dòng)面上的受力大小，不研究滑體及滑床內(nèi)部的應(yīng)力狀態(tài)。(3)將滑體的邊界條件大大簡(jiǎn)化。如將復(fù)雜的滑體型態(tài)簡(jiǎn)化為簡(jiǎn)單的幾何型態(tài)；將滑面簡(jiǎn)化為圓弧面、平面或折面；一般將立體問題簡(jiǎn)化為平面問題，取沿滑動(dòng)方向的代表性剖面，以表征滑體的基本型態(tài)；將均布力簡(jiǎn)化為集中力，有時(shí)還將力的作用點(diǎn)簡(jiǎn)化為通過滑體重心。極限平衡方

7、法包括以下幾種方法：(1)瑞典圓弧滑動(dòng)法(2)簡(jiǎn)化逼肖普法(3)簡(jiǎn)布普通條分法(4)摩根斯坦-普賴斯法(5)不平衡推力傳遞法以上各種方法都是假定土體是理想塑性材料，把土條作為一個(gè)剛體，按照極限平衡的原則進(jìn)行力的分析，最大的不同之處在于對(duì)相鄰上條之間的內(nèi)力作何種假定，也就是如何增加已知條件使超靜定問題變成靜定問題。這些假定的物理意義不一樣，所能滿足的平衡條件也不相同，計(jì)算步驟有繁有簡(jiǎn)，使用時(shí)必須注意他們的適用場(chǎng)合。極限平衡方法關(guān)鍵是對(duì)滑體的休型和滑面的形態(tài)進(jìn)行分析、正確選用滑面的計(jì)算參數(shù)以及正確引用滑體的荷載條件等。因?yàn)闃O限平衡方法完全不考慮土體本身的應(yīng)力-應(yīng)變關(guān)系，不能真實(shí)地反映邊坡失穩(wěn)時(shí)的應(yīng)

8、力場(chǎng)和位移場(chǎng)，因此而受到質(zhì)疑。4.1.4.2 數(shù)值分析方法數(shù)值數(shù)值分析方法考慮土體應(yīng)力應(yīng)變關(guān)系，克服了極限平衡方法完全不考慮土體本身的應(yīng)力-應(yīng)變關(guān)系缺點(diǎn)，為邊坡穩(wěn)定分析提供了較為正確和深入的概念。邊坡穩(wěn)定性數(shù)值分析方法主要包含以下幾種方法：(1)有限元法有限單元法是數(shù)值模擬方法在邊坡穩(wěn)定評(píng)價(jià)中應(yīng)用最早的方法，也是目前最廣泛使用的一種數(shù)值方法，可以用來求解彈性、彈塑性、粘彈塑性、粘塑性等問題。目前 用有限元法求解邊坡穩(wěn)定主要有兩種方法。a.有限元滑面搜索法:將邊坡體離散為有限單元格，按照施加的荷載及邊界條件進(jìn)行有限元計(jì)算可得到每個(gè)結(jié)點(diǎn)的應(yīng)力張量。然后假定一個(gè)滑動(dòng)面，用有限元數(shù)據(jù)給出滑動(dòng)面任一點(diǎn)的

9、向正應(yīng)力和剪應(yīng)力，根據(jù)摩爾一庫侖準(zhǔn)則可得該點(diǎn)的抗滑力，由此即能求得滑動(dòng)面上每個(gè)結(jié)點(diǎn)的下滑力與抗滑力，再對(duì)滑動(dòng)面上下滑力與抗滑力進(jìn)行積分，就可以求得每一個(gè)滑動(dòng)面的安全系數(shù)。b.有限元強(qiáng)度折減法:首先選取初始折減系數(shù)，將巖土體強(qiáng)度參數(shù)進(jìn)行折減，將折減后的參數(shù)作為輸入，進(jìn)行有限元計(jì)算，若程序收斂，則巖土體仍處于穩(wěn)定狀態(tài)，然后再增加折減系數(shù)，直到程序恰好不收斂，此時(shí)的折減系數(shù)即為穩(wěn)定或安全系數(shù)。(2)自適應(yīng)有限元法自20世紀(jì)70年代開始自適應(yīng)理論被引入有限元計(jì)算，主導(dǎo)思想是減少前處理工作量和實(shí)現(xiàn)網(wǎng)格離散的客觀控制。現(xiàn)已基本建立了一般彈性力學(xué)、流體動(dòng)力學(xué)、滲流分析等領(lǐng)域的平面自適應(yīng)分析系統(tǒng)，能使計(jì)算較為

10、快速和準(zhǔn)確。(3)離散單元法離散單元法的突出功能是它在反映巖塊之間接觸面的滑移、分離與傾翻等大位移的同時(shí)，又能計(jì)算巖塊內(nèi)部的變形與應(yīng)力分布。因此，任何一種巖體材料都可引入到模型中，例如彈性、粘彈性或斷裂等均可考慮，故該法對(duì)塊狀結(jié)構(gòu)、層狀破裂或一般破裂結(jié)構(gòu)巖體邊坡比較合適。并且，它利用顯式時(shí)間差分法(動(dòng)態(tài)差分法)求解動(dòng)力平衡方程，求解非線性大位移與動(dòng)力問題比較容易。離散元法在模擬過程中考慮了邊坡失穩(wěn)破壞的動(dòng)態(tài)過程，允許巖土體存在滑動(dòng)、平移、轉(zhuǎn)動(dòng)和巖體的斷裂及松散等復(fù)雜過程，具有宏觀上的不連續(xù)性和單個(gè)巖塊休運(yùn)動(dòng)的隨機(jī)性，可以較真實(shí)、動(dòng)態(tài)地模擬邊坡在形成和開挖過程中應(yīng)力、位移和狀態(tài)的變化，預(yù)測(cè)邊坡的

11、穩(wěn)定性，因此在巖質(zhì)高邊坡穩(wěn)定性的研究中得到廣泛的應(yīng)用。(4)拉格朗日元法為了克服有限元等方法不能求解大變形問題的缺陷，人們根據(jù)有限差分法的原理，提出了FLAC數(shù)值分析方法。該方法較有限元法能更好地考慮巖土體的不連續(xù)和大變形特性，求解速度較快。缺點(diǎn)是計(jì)算邊界、單元網(wǎng)格的劃分帶有很大的隨意性。(5)界面元法界面元法是一種基于累積單元變形于界面的界面應(yīng)力元法模型，建立適用于分析不連續(xù)、非均勻、各向異性和各類非線性問題、場(chǎng)問題，以及能夠完全模擬各類錨桿復(fù)雜空間布局和開挖擾動(dòng)的方法。4.1.4.3 有限元法用于邊坡穩(wěn)定性分析優(yōu)點(diǎn) 有限元法考慮了介質(zhì)的變形特征，真實(shí)地反應(yīng)了邊坡的受力狀態(tài)。它可以模擬連續(xù)介

12、質(zhì)，也可以模擬不連續(xù)介質(zhì)；能考慮邊坡沿軟弱結(jié)構(gòu)面的破壞，也能分析邊坡的整體穩(wěn)定破壞。有限元法可以模擬邊坡的圓弧滑動(dòng)破壞和非圓弧滑動(dòng)破壞。同時(shí)它還能適應(yīng)各種邊界條件和不規(guī)則幾何形狀，具有很廣泛的適用性。有限元法應(yīng)用于邊坡工程，有其獨(dú)特的優(yōu)越性。與一般解析方法相比，有限元法有以下優(yōu)點(diǎn)：(1) 它考慮了巖體的應(yīng)力-應(yīng)變關(guān)系，求出每一單元的應(yīng)力與變形，反映了巖體真實(shí)工作狀態(tài)。(2)與極限平衡法相比，不需要進(jìn)行條間力的簡(jiǎn)化，巖體自始至終處于平衡狀態(tài)。(3)不需要像極限平衡法一樣事先假定邊坡的滑動(dòng)面，邊坡的變形特性、塑性區(qū)形成都根據(jù)實(shí)際應(yīng)力應(yīng)變狀態(tài)“自然”形成。(4)若巖體的初始應(yīng)力己知，可以模擬有構(gòu)造應(yīng)

13、力邊坡的受力狀態(tài)。(5)不但能像極限平衡法一樣模擬邊坡的整體破壞，還能模擬邊坡的局部破壞，把邊坡的整體破壞和局部破壞納入統(tǒng)一的體系。(6)可以模擬邊坡的開挖過程，描述和反應(yīng)巖體中存在的節(jié)理裂隙、斷層等構(gòu)造面。鑒于有限元法具有如此多優(yōu)點(diǎn)，本章借助通用有限元軟件ANSYS來實(shí)現(xiàn)對(duì)邊坡穩(wěn)定性分析，用具體的邊坡工程實(shí)例詳細(xì)介紹應(yīng)用ANSYS軟件分析邊坡穩(wěn)定性問題。4.2 ANSYS邊坡穩(wěn)定性分析步驟ANSYS邊坡穩(wěn)定性分析一般分以下幾個(gè)步驟：1、創(chuàng)建物理環(huán)境2、建立模型，劃分網(wǎng)格，對(duì)模型的不同區(qū)域賦予特性3、加邊界條件和載荷4、求解5、后處理（查看計(jì)算結(jié)果）4.2.1創(chuàng)建物理環(huán)境在定義邊坡穩(wěn)定性分析問

14、題的物理環(huán)境時(shí)，進(jìn)入ANSYS前處理器，建立這個(gè)邊坡穩(wěn)定性分析的數(shù)學(xué)仿真模型。按照以下幾個(gè)步驟來建立物理環(huán)境：1、 設(shè)置GUT菜單過濾如果你希望通過GUI路徑來運(yùn)行ANSYS，當(dāng)ANSYS被激活后第一件要做的事情就是選擇菜單路徑：Main Menu>Preferences，執(zhí)行上述命令后，彈出一個(gè)如圖4-1所示的對(duì)話框出現(xiàn)后，選擇Structural。這樣ANSYS會(huì)根據(jù)你所選擇的參數(shù)來對(duì)GUI圖形界面進(jìn)行過濾，選擇Structural以便在進(jìn)行邊坡穩(wěn)定性分析時(shí)過濾掉一些不必要的菜單及相應(yīng)圖形界面。2、 定義分析標(biāo)題（TITLE）在進(jìn)行分析前，可以給你所要進(jìn)行的分析起一個(gè)能夠代表所分析內(nèi)

15、容的標(biāo)題，比如“Slope stability Analysis”，以便能夠從標(biāo)題上與其他相似物理幾何模型區(qū)別。用下列方法定義分析標(biāo)題。命令：TITLEGUI：Utility Menu>File>Change Title3、 說明單元類型及其選項(xiàng)（KEYOPT選項(xiàng)）與ANSYS的其他分析一樣，也要進(jìn)行相應(yīng)的單元選擇。ANSYS軟件提供了100種以上的單元類型，可以用來模擬工程中的各種結(jié)構(gòu)和材料，各種不同的單元組合在一起，成為具體的物理問題的抽象模型。例如，不同材料屬性的邊坡土體用PLANE82單元來模擬。大多數(shù)單元類型都有關(guān)鍵選項(xiàng)（KEYOPTS），這些選項(xiàng)用以修正單元特性。例如，

16、PLANE82有如下KEYOPTS：KEYOPT(2) 包含或抑制過大位移設(shè)置KEYOPT(3) 平面應(yīng)力、軸對(duì)稱、平面應(yīng)變或考慮厚度的平面應(yīng)力設(shè)置KEYOPT(9) 用戶子程序初始應(yīng)力設(shè)置設(shè)置單元以及其關(guān)鍵選項(xiàng)的方式如下：命令：ET KEYOPTGUI：Main Menu> Preprocessor> Element Type> Add/Edit/Delete圖 4-1 GUI圖形界面過濾4定義單位結(jié)構(gòu)分析只有時(shí)間單位、長(zhǎng)度單位和質(zhì)量單位三個(gè)基本單位，則所有輸入的數(shù)據(jù)都應(yīng)當(dāng)是這三個(gè)單位組成的表達(dá)方式。如標(biāo)準(zhǔn)國(guó)際單位制下，時(shí)間是秒（s），長(zhǎng)度是米（m），質(zhì)量是千克（kg），則

17、導(dǎo)出力的單位是kgm/s2（相當(dāng)于牛頓N），材料的彈性模量單位是kg/ms2（相當(dāng)于帕Pa）。命令：/UNITS5、定義材料屬性大多數(shù)單元類型在進(jìn)行程序分析時(shí)都需要指定材料特性，ANSYS程序可方便地定義各種材料的特性，如結(jié)構(gòu)材料屬性參數(shù)、熱性能參數(shù)、流體性能參數(shù)和電磁性能參數(shù)等。ANSYS程序可定義的材料特性有以下三種：（1）線性或非線性。（2）各向同性、正交異性或非彈性。（3）隨溫度變化或不隨溫度變化。因?yàn)榉治龅倪吰履Ｐ筒捎美硐霃椝苄阅Ｐ停―-P模型），因此邊坡穩(wěn)定性分析中需要定義邊坡中不同土體的材料屬性：容重、彈性模量、泊松比、凝聚力以及摩擦角。命令：MPGUI：Main Menu>

18、;Preprocessor>Material Props>Material Models 或 Main Menu>Solution>Load Step Opts>Other>Change Mat Props>Material Models進(jìn)行邊坡穩(wěn)定性分析計(jì)算時(shí)，采用強(qiáng)度折減法來實(shí)現(xiàn)。首先選取初始折減系數(shù)F，然后對(duì)邊坡土體材料強(qiáng)度系數(shù)進(jìn)行折減，折減后凝聚力以及摩擦角分別見式4-1和式4-2。 (4-1) (4-2)u 和為邊坡土體的初始凝聚力和摩擦角。u 對(duì)和進(jìn)行折減，輸入邊坡模型計(jì)算，若收斂，則此時(shí)邊坡是穩(wěn)定的；繼續(xù)增大折減系數(shù)F，直到程序恰好不收斂

19、，此時(shí)的折減系數(shù)即為穩(wěn)定或安全系數(shù)。4.2.2 建立模型和劃分網(wǎng)格 創(chuàng)建好物理環(huán)境，就可以建立模型。在進(jìn)行邊坡穩(wěn)定性分析時(shí)，需要建立模擬邊坡土體的PLANE82單元。在建立好的模型各個(gè)區(qū)域內(nèi)指定特性（單元類型、選項(xiàng)、實(shí)常數(shù)和材料性質(zhì)等）以后，就可以劃分有限元網(wǎng)格了。通過GUI為模型中的各區(qū)賦予特性：1、選擇Main Menu> Preprocessor> Meshing> Mesh Attributes> Picked Areas2、點(diǎn)擊模型中要選定的區(qū)域。3、在對(duì)話框中為所選定的區(qū)域說明材料號(hào)、實(shí)常數(shù)號(hào)、單元類型號(hào)和單元坐標(biāo)系號(hào)。4、重復(fù)以上三個(gè)步驟，直至處

20、理完所有區(qū)域。通過命令為模型中的各區(qū)賦予特性：ASEL（選擇模型區(qū)域）MAT （說明材料號(hào)）REAL（說明實(shí)常數(shù)組號(hào)）TYPE（指定單元類型號(hào)）ESYS（說明單元坐標(biāo)系號(hào)） 4.2.3 施加約束和荷載在施加邊界條件和荷載時(shí)，既可以給實(shí)體模型（關(guān)鍵點(diǎn)、線、面）也可以給有限元模型（節(jié)點(diǎn)和單元）施加邊界條件和荷載。在求解時(shí)，ANSYS程序會(huì)自動(dòng)將加到實(shí)體模型上的邊界條件和載荷轉(zhuǎn)遞到有限元模型上。邊坡穩(wěn)定性分析中，主要是給邊坡兩側(cè)和底部施加自由度約束。命令：D施加荷載包括自重荷載以及邊坡開挖荷載。 求解 接著就可以進(jìn)行求解，ANSYS程序根據(jù)現(xiàn)有選項(xiàng)的設(shè)置，從數(shù)據(jù)庫獲取模型和載荷信息并進(jìn)行計(jì)算求解，將

21、結(jié)果數(shù)據(jù)寫入到結(jié)果文件和數(shù)據(jù)庫中。命令：SOLVEGUI：Main Menu>Solution>Solve> Current LS4.2.5 后處理 后處理的目的是以圖和表的形式描述計(jì)算結(jié)果。對(duì)于邊坡穩(wěn)定性分析中，進(jìn)入后處理器后，查看邊坡變形圖和節(jié)點(diǎn)的位移、應(yīng)力和應(yīng)變。隨著強(qiáng)度折減系數(shù)的增大，邊坡的水平位移增大，塑性應(yīng)變急劇發(fā)展，塑性區(qū)發(fā)展形成一個(gè)貫通區(qū)域時(shí)，計(jì)算不收斂，認(rèn)為邊坡發(fā)生了破壞。通過研究位移、應(yīng)變和塑性區(qū)域，來綜合判斷邊坡的穩(wěn)定性。命令：POST1GUI： Main Menu> General Postproc4.2.6 補(bǔ)充說明邊坡的失穩(wěn)破壞定義有很多種，

22、對(duì)于采用彈塑性計(jì)算模型的邊坡，需要綜合考慮以下因素:（1）把有限元計(jì)算的收斂與否作為一個(gè)重要的衡量指標(biāo)，邊坡處于穩(wěn)定狀態(tài)，計(jì)算收斂，邊坡破壞時(shí)，邊坡不收斂。（2）邊坡失穩(wěn)的同時(shí)還表現(xiàn)出位移急劇增加。 （3）邊坡失穩(wěn)總是伴隨著塑性變形的明顯增加和塑性區(qū)的發(fā)展，塑性區(qū)的發(fā)展?fàn)顩r反映了邊坡是否處于穩(wěn)定狀態(tài)。此外，采用彈塑性有限元法進(jìn)行計(jì)算，它具有獨(dú)特的優(yōu)勢(shì)：(1)彈塑性分析假定巖體為彈塑性材料，巖體在受力初期處于彈性狀態(tài)，達(dá)到一定的屈服準(zhǔn)則后，處于塑性狀態(tài)。采用彈塑性模型更能反應(yīng)巖體的實(shí)際工作狀態(tài)。(2)巖體所承受的荷載超過材料強(qiáng)度時(shí)，就會(huì)出現(xiàn)明顯的滑移破壞面。因此，彈塑性計(jì)算不需要假定破壞面的形狀

23、和位置，破壞面根據(jù)剪應(yīng)力強(qiáng)度理論自動(dòng)形成。當(dāng)整個(gè)邊坡破壞時(shí)，就會(huì)出現(xiàn)明顯的塑性區(qū)。(3 )能綜合考慮邊坡的局部失穩(wěn)和整體失穩(wěn)破壞。4.3 ANSYS邊坡穩(wěn)定性實(shí)例分析4.3.1 實(shí)例描述圖4-2 邊坡模型邊坡實(shí)例選取國(guó)內(nèi)某礦，該邊坡考慮彈性和塑性兩種材料，邊坡尺寸如圖4-2所示。分析目的是對(duì)該邊坡進(jìn)行穩(wěn)定性計(jì)算分析，以判斷其穩(wěn)定性和計(jì)算出安全系數(shù)，該邊坡圍巖材料屬性見表4-1。表4-1 邊坡模型圍巖參數(shù)類別彈性模量/Gpa泊松比容重/內(nèi)聚力/Mpa摩擦角（。）圍巖2（彈塑性）300.2525000.942圍巖1（彈性）310.242700-對(duì)于像邊坡這樣縱向很長(zhǎng)的實(shí)體，計(jì)算模型可以簡(jiǎn)化為平面應(yīng)

24、變問題。假定邊坡所承受的外力不隨Z軸變化，位移和應(yīng)變都發(fā)生在自身平面內(nèi)。對(duì)于邊坡變形和穩(wěn)定性分析，這種平面假設(shè)是合理的。實(shí)測(cè)經(jīng)驗(yàn)表明，邊坡的影響范圍在2倍坡高范圍，因此本文計(jì)算區(qū)域?yàn)檫吰麦w橫向延伸2倍坡高，縱向延伸3倍坡高。兩側(cè)邊界水平位移為零，下側(cè)邊界豎向位移為零。彈性有限元的計(jì)算模型如圖4-2所示。采用雙層模型，模型上部為理想彈塑性材料，下部為彈性材料，左右邊界水平位移為零，下邊界豎向位移為零。u 雙層模型考慮土體的彈塑性變形，其塑性區(qū)的發(fā)展，應(yīng)力的分布更符合實(shí)際情況。u 考慮雙層模型，塑性區(qū)下部的單元可以產(chǎn)生一定的垂直變形和水平變形，基本消除了由于邊界效應(yīng)在邊坡下部出現(xiàn)的塑性區(qū)，更好地模

25、擬了邊坡的變形和塑性區(qū)的發(fā)展。4.3.2 GUI操作方法4.3.2.1 創(chuàng)建物理環(huán)境1) 在【開始】菜單中依次選取【所有程序】/【ANSYS10.0】/【ANSYS Product Launcher】，得到“10.0ANSYS Product Launcher”對(duì)話框。2）選中【File Management】，在“Working Directory”欄輸入工作目錄“D:ansysexample4-1”，在“Job Name”欄輸入文件名“Slope”。3）單擊“RUN”按鈕，進(jìn)入ANSYS10.0的GUI操作界面。4）過濾圖形界面：Main Menu> Preferences，彈出“P

26、references for GUI Filtering”對(duì)話框，選中“Structural”來對(duì)后面的分析進(jìn)行菜單及相應(yīng)的圖形界面過濾。5）定義工作標(biāo)題：Utility Menu> File> Change Title，在彈出的對(duì)話框中輸入“Slope stability Analysis”，單擊“OK”，如圖4-3。圖4-3 定義工作標(biāo)題6）定義單元類型：a.定義PLANE82單元：Main Menu> Preprocessor> Element Type> Add/Edit/Delete，彈出一個(gè)單元類型對(duì)話框，單擊“Add”按鈕。彈出如圖4-4所示對(duì)話框。

27、在該對(duì)話框左面滾動(dòng)欄中選擇“Solid”，在右邊的滾動(dòng)欄中選擇“Quad 8node 82”，單擊“Apply”，就定義了“PLANE82”單元。圖4-4 定義PLANE82單元對(duì)話框 b.設(shè)定PLANE82單元選項(xiàng)：Main Menu> Preprocessor> Element Type> Add/Edit/Delete，彈出一個(gè)單元類型對(duì)話框，選中“Type 2 PLANE82”，單擊“Options”按鈕，彈出一個(gè)“PLANE82 element Type options”對(duì)話框，如圖4-5所示。在“Element behavior K3”欄后面的下拉菜單中選取“Pl

28、ane strain”，其它欄后面的下拉菜單采用ANSYS默認(rèn)設(shè)置就可以，單擊“OK”按鈕。圖4-5 PLANE82單元庫類型選項(xiàng)對(duì)話框u 通過設(shè)置PLANE82單元選項(xiàng)“K3”為“Plane strain”來設(shè)定本實(shí)例分析采取平面應(yīng)變模型進(jìn)行分析。因?yàn)檫吰率强v向很長(zhǎng)的實(shí)體，故計(jì)算模型可以簡(jiǎn)化為平面應(yīng)變問題。u 8節(jié)點(diǎn)PLANE82單元每個(gè)節(jié)點(diǎn)有UX和UY兩個(gè)自由度，比4節(jié)點(diǎn)PLANE42單元具有更高的精確性，對(duì)不規(guī)則網(wǎng)格適應(yīng)性更強(qiáng)。7）定義材料屬性a.定義邊坡圍巖1材料屬性：Main Menu> Preprocessor> Material Props> Material

29、Models，彈出“Define Material Model Behavior”對(duì)話框，如圖4-6所示。圖 4-6 定義材料本構(gòu)模型對(duì)話框 在圖4-6中右邊欄中連續(xù)雙擊“Structural> Linear> Elastic>Isotropic”后，又彈出如圖4-7所示“Linear Isotropic Properties for Material Number 1”對(duì)話框，在該對(duì)話框中“EX”后面的輸入欄輸入“3E10”，在“PRXY” 后面的輸入欄輸入“0.25”，單擊“OK”。再在選中“Density”并雙擊，彈出如圖4-8所示“Density for Materi

30、al Number 1”對(duì)話框，在“DENS”后面的欄中輸入邊坡土體材料的密度“2500”，單擊“OK”按鈕。再次在圖4-6中右邊的欄中連續(xù)雙擊“Structural>Nonlinear> Inelastic>Non-metal plasticity>drucker-prager”后，又彈出一個(gè)如圖4-9所示對(duì)話框。在“Cohesion”欄添入邊坡圍巖材料1的內(nèi)聚力“0.9E6”，在“Fric Angle” 欄添入邊坡內(nèi)摩擦角“42”，單擊“OK”按鈕。圖 4-7 線彈性材料模型對(duì)話框 圖 4-8材料密度輸入對(duì)話框圖4-9 定義邊坡材料1DP模型對(duì)話框 b.定義邊坡圍巖

31、2材料屬性：在圖4-6對(duì)話框中，單擊“Material> New Model”， 彈出一個(gè)“Define Material ID”對(duì)話框，在“ID”欄后面輸入材料編號(hào)“2”，單擊“OK”按鈕。彈出一個(gè)定義材料模型對(duì)話框?qū)υ捒?，選中“Material Model Number 2”，和定義邊坡圍巖1材料一樣，在右邊的欄中連續(xù)雙擊“Structural> Linear> Elastic>Isotropic”后，又彈出一個(gè)“Linear Isotropic Properties for Material Number 2”對(duì)話框，在該對(duì)話框中“EX”后面的輸入欄輸入“3.2E

32、10”，在“PRXY” 后面的輸入欄輸入“0.24”，單擊“OK”。再選中“Density”并雙擊，彈出一個(gè)“Density for Material Number 2”對(duì)話框，在“DENS”后面的欄中輸入隧道圍巖材料的密度“2700”，再單擊“OK”按鈕，彈出一個(gè)定義材料模型對(duì)話框。 c.復(fù)制邊坡圍巖1材料性質(zhì)：在圖4-6對(duì)話框中，用鼠標(biāo)點(diǎn)擊“Edit>copy.”，彈出一個(gè)“Copy Material Model”對(duì)話框，如圖4-10所示。在“from Material number”欄后面的下拉菜單中選取“1”，在“TO Material number”欄后面輸入“3”，單擊“Ap

33、ply”按鈕。又彈出如土4-10所示對(duì)話框，然后依次在“TO Material number”欄后面輸入“4”、“5”、“6”、“7”、“8”“9”、“10”，“11”、“12”、“13”，每輸入一個(gè)數(shù)，就單擊“Apply”按鈕一次。圖 4-10 復(fù)制本構(gòu)模型對(duì)話框 最后得到10個(gè)復(fù)制圍巖1的邊坡材料本構(gòu)模型，如圖4-11所示。圖4-11 定義強(qiáng)度折減后材料模型對(duì)話框圖4-12 定義強(qiáng)度折減系數(shù)F=1.2時(shí)圍巖材料對(duì)話框d.定義10個(gè)強(qiáng)度折減后材料本構(gòu)模型：首先定義強(qiáng)度折減系數(shù)F=1.2后邊坡圍巖材料模型，在圖4-11對(duì)話框中，在鼠標(biāo)依次雙擊“Material Model Number 3/D

34、rucper-Prager”。彈出一個(gè)“Drucker- Prager Material Number 3”，如圖4-12所示，在“Cohesion”欄添入強(qiáng)度折減系數(shù)F=1.2后邊坡圍巖材料1的內(nèi)聚力“0.75E6”，在“Fric Angle” 欄添入折減后邊坡內(nèi)摩擦角“37.7”，單擊“OK”按鈕。用相同方法定義強(qiáng)度折減系數(shù)分別為：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的邊坡圍巖材料本構(gòu)模型。u 定義強(qiáng)度折減后本構(gòu)模型目的是為了分析邊坡穩(wěn)定性。u 強(qiáng)度折減就是降低內(nèi)聚力和摩擦角，根據(jù)式4-1和式4-2進(jìn)行折減。4.3.2.

35、2 建立模型和劃分網(wǎng)格1）創(chuàng)建邊坡線模型a.輸入關(guān)鍵點(diǎn)：Main Menu> Preprocessor> Modeling> Create>Keypoints>In Active CS，彈出“Creae Keypoints in Active Cooedinate System”對(duì)話框，如圖4-13所示。在“NPT keypoint number”欄后面輸入“1”，在“X，Y，Z Location in active CS”欄后面輸入“（0，0，0）”，單擊“Apply”按鈕，這樣就創(chuàng)建了關(guān)鍵點(diǎn)1。再依次重復(fù)在“NPT keypoint number”欄后面輸入“

36、2、3、4、5、6、7、8、9”，在對(duì)應(yīng)“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在當(dāng)前坐標(biāo)系創(chuàng)建關(guān)鍵點(diǎn)對(duì)話框b.創(chuàng)建邊坡線模型：Main Menu> Preprocessor> Modeling> Create>Lines>Straight line，彈出“Creae straight lines

37、”對(duì)話框，用鼠標(biāo)依次點(diǎn)擊關(guān)鍵點(diǎn)1、2，單擊“Apply”按鈕，這樣就創(chuàng)建了直接L1，同樣分別連接關(guān)鍵點(diǎn)“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）創(chuàng)建邊坡面模型a.打開面編號(hào)顯示：Utility Menu> PlotCtrls> Numbering，彈出“Plot Numbering Controls”對(duì)話框，如圖4-15所示。選中“Aares Numbers” 選項(xiàng)，后面的文字由“off”變?yōu)椤皁n”,單擊“OK”關(guān)閉窗口

38、。圖4-15 打開面編號(hào)對(duì)話框圖4-16 邊坡面模型b.創(chuàng)建邊坡面模型：Main Menu> Preprocessor> Modeling> Create>Areas> Arbitrary> by line，彈出一個(gè)“Create Area by lines”對(duì)話框，在圖形中選取線L4、L5、L3和L11，點(diǎn)擊“Apply ”按鈕，就生成了邊坡彈性材料區(qū)域面積A1；再依次用鼠標(biāo)在圖形中選取線L1、L2、L6、L10、和L11，點(diǎn)擊“Apply ”按鈕，就生成了邊坡塑性材料區(qū)域面積A2；再依次用鼠標(biāo)在圖形中選取線L7、L8、L10和L9，點(diǎn)擊“OK ”按鈕，就

39、生成了邊坡開挖掉區(qū)域面積A3。最后得到邊坡模型的面模型，如圖4-16所示。 4）劃分邊坡圍巖2單元網(wǎng)格 a. .給邊坡圍巖2賦予材料特性：Main Menu> Preprocessor> Meshing> MeshTool，彈出“MeshTool”對(duì)話框，如圖4-17所示。在“Element Attributes”后面的下拉式選擇欄中選擇“Areas”，按“Set”按鈕，彈出一個(gè)“Areas Attributes”面拾取框，在圖形界面上拾取邊坡圍巖2區(qū)域，單擊拾取框上的“OK”按鈕，又彈出一個(gè)如圖4-18所示的“Areas Attributes”對(duì)話框，在“Material

40、number”后面的下拉式選擇欄中選取“2”，在“Element type number ” 后面的下拉式選擇欄中選取“2 PLANE82”，單擊“Apply”。圖4-17 網(wǎng)格劃分工具欄 圖4-18定義單元屬性對(duì)話框 b.設(shè)置網(wǎng)格劃分份數(shù)：在圖4-17工具欄中“Size Control”欄，用鼠標(biāo)點(diǎn)擊“l(fā)ines”后面的“Set”，彈出一個(gè)選擇對(duì)話框，用鼠標(biāo)在圖形選擇線L3和L5。彈出一個(gè)“Element Sizes on Picked Lines”，對(duì)話框，如圖4-19所示，在“No of element division”欄后面輸入“5”，單擊“Apply”按鈕，再選擇線L4和L11，又

41、彈出圖4-19對(duì)話框，在“No of element division”欄后面輸入“26”，單擊“OK”按鈕。c.劃分單元網(wǎng)格：在圖4-17網(wǎng)格劃分工具欄中單擊“Mesh”按鈕，彈出一個(gè)拾取面積對(duì)話框，拾取面積A1，單擊拾取框上的“OK”按鈕，生成邊坡圍巖2單元網(wǎng)格。圖 4-19 設(shè)置網(wǎng)格份數(shù)對(duì)話框5）劃分邊坡圍巖1單元網(wǎng)格 圖4-20 選取線對(duì)話框 圖4-21 設(shè)置網(wǎng)格分?jǐn)?shù)對(duì)話框a.設(shè)置網(wǎng)格份數(shù)：Main Menu>Preprocessor>Meshing>Size Cntrls>ManualSize>Layers>Picked Lines，彈出一個(gè)“Se

42、t Layer Controls”對(duì)話框，如圖4-20所示，用鼠標(biāo)選取線L1、L1和L6單擊“OK”按鈕。彈出一個(gè)“Area Layer Mesh Control on Picked lines”對(duì)話框，如圖4-21所示，在“No of line division”欄后面輸入“10”，單擊“OK”按鈕。相同方法設(shè)置線L8和L10分割份數(shù)為16；設(shè)置線L7和L9線的分割份數(shù)為12。 b. 給邊坡圍巖1賦予材料特性：Main Menu> Preprocessor> Meshing> MeshTool，彈出“MeshTool”對(duì)話框，如圖4-17所示。在“Element Attri

43、butes”后面的下拉式選擇欄中選擇“Areas”，按“Set”按鈕，彈出一個(gè)“Areas Attributes”面拾取框，在圖形界面上拾取面A2和A3，單擊拾取框上的“OK”按鈕，又彈出一個(gè) “Areas Attributes”對(duì)話框，在“Material number”后面的下拉式選擇欄中選取“1”，在“Element type number ” 后面的下拉式選擇欄中選取“2 PLANE82”，單擊“Apply”。 c.劃分單元網(wǎng)格：在圖3-97網(wǎng)格劃分工具欄中單擊“Mesh”按鈕，彈出一個(gè)拾取面積對(duì)話框，拾取圍巖，單擊拾取框上的“OK”按鈕，生成邊坡圍巖1單元網(wǎng)格。最后得到邊坡模型單元網(wǎng)

44、格，如圖4-22所示。圖4-22 邊坡模型單元網(wǎng)格4.3.2.3 施加約束和荷載 1）給邊坡模型施加約束a.給邊坡模型兩邊施加約束：執(zhí)行Main Menu>Solution>Define Loads>Apply>Structural>Displacement>on Nodes，彈出在節(jié)點(diǎn)上施加位移約束對(duì)話框，用鼠標(biāo)選取隧道模型兩側(cè)邊界上所有節(jié)點(diǎn)，單擊“OK”按鈕。彈出“Apply U,ROT on Nodes”對(duì)話框，如圖4-23所示，在“DOFS to be constrained”欄后面中選取“UX”，在“Apply as”欄后面的下拉菜單中選取“Con

45、stant value”，在“Displacement value”欄后面輸入“0”值，然后單擊“OK”按鈕。圖4-23 給模型兩側(cè)施加位移約束對(duì)話框 圖4-24 給模型底部施加位移約束對(duì)話框b.給模型底部施加約束：執(zhí)行Main Menu>Solution>Define Loads> Apply>Structural>Displacement>on Nodes，彈出在節(jié)點(diǎn)上施加位移約束對(duì)話框，用鼠標(biāo)選取隧道模型底部邊界上所有節(jié)點(diǎn)，單擊“OK”按鈕。彈出圖4-24所示對(duì)話框，在“DOFS to be constrained”欄后面中選取“UX、UY”，在“Ap

46、ply as”欄后面的下拉菜單中選取“Constant value”，在“Displacement value”欄后面輸入“0”值，然后單擊“OK”按鈕。u 節(jié)點(diǎn)選擇，可以先選擇節(jié)點(diǎn)上線，再選擇附在線上的節(jié)點(diǎn)。 2）施加重力加速度：Main Menu>Solution>Define Loads>Apply>Structural>Inertia>Gravity，彈出“Apply(Gravitational)Acceleration”對(duì)話框，如圖4-25所示。只需在“Global Cartesian Y-comp”欄后面輸入重力加速度值“9.8”就可以，單擊“O

47、K”按鈕，就完成了重力加速度的施加。圖4-25 施加重力加速度對(duì)話框 這時(shí)就可以得到施加約束和重力加速度后隧道有限元模型，如圖4-25所示。圖4-25 施加約束和重力荷載后邊坡模型.4 求解1）求解設(shè)置a.指定求解類型：Main Menu>Solution>Analysis Type>New Analysis，彈出一個(gè)如圖4-25所示對(duì)話框，在“Type of analysis”欄后面選中“Static”，單擊“OK”按鈕。圖4-25 指定求解類型對(duì)話框b.設(shè)置載荷步： Main Menu>Preprocessor>Loads>Analysis Type&g

48、t;Sol'n Controls，彈出一個(gè)“Solution Controls”對(duì)話框，用鼠標(biāo)單擊“Basic”選項(xiàng)，如圖4-26所示，在“Number of Substeps”欄后面輸入“5”，在“Max no. of substeps”欄后面輸入“100”，在“Min no. of substeps”欄后面輸入“1”，單擊“OK”按鈕。圖4-26設(shè)置載荷步對(duì)話框 c.設(shè)置線性搜索：Main Menu>Solution>Analysis Type>Sol'n Controls，彈出一個(gè)“Solution Controls”對(duì)話框，用鼠標(biāo)單擊“Nonlinea

49、r”選項(xiàng)，如圖4-27所示，在“Line search”欄后面下拉菜單選中“ON”，單擊“OK”。圖4-27 設(shè)置線性搜索對(duì)話框 d. 設(shè)定牛頓-拉普森選項(xiàng)：Main Menu>Solution>Analysis Type>Analysis Options，彈出一個(gè)“Static or Steady-static Analysis”對(duì)話框，如圖4-28所示，在“New-Raphson option”欄后面下拉菜單選中“Full N-R”，單擊“OK”按鈕。圖4-28 設(shè)定牛頓-拉普森選項(xiàng)e.打開大位移求解：Main Menu>Preprocessor>Loads&

50、gt;Analysis Type>Sol'n Controls>Basic，得到圖4-26所示對(duì)話框，在“Analysis Options”欄后面下拉菜單選中大位移“Large Displacement static”，單擊“OK”按鈕。f.設(shè)置收斂條件：Main Menu>Preprocessor>Loads>Load Step Opts>Nonlinear>Convergence Crit，彈出一個(gè)“Default Nonlinear Convergence Criteria”，如圖4-29所示。圖中顯示ANSYS默認(rèn)的收斂條件：分別設(shè)置了

51、力和力矩的收斂條件。圖4-29 ANSYS默認(rèn)收斂條件 為了使求解順利進(jìn)行和得到較好解，可以修改默認(rèn)收斂設(shè)置，可以分別設(shè)置力、力矩和位移收斂條件。單擊圖4-29中的“Replace”，彈出一個(gè)“Nonlinear Convergence Criteria”對(duì)話框，如圖4-30所示。在“Lab Convergence based on”欄后面第一欄中選中“Structural”，第二欄中選中“Force F”；在“TOLER Tolerance about value”欄后面輸入“0.005”；在“NORM Convergence norm”欄后面下拉菜單選中“L2 norm”；在“MNREF

52、Minimum reference value”后面欄中輸入“0.5”，單擊“OK”按鈕，這就設(shè)置好了求解時(shí)力收斂條件。圖4-30 設(shè)置力收斂條件對(duì)話框 單擊圖4-30對(duì)話框中“OK”按鈕后，彈出一個(gè)如圖4-31所示的對(duì)話框，單擊“Add”按鈕。彈出一個(gè)如圖4-31所示對(duì)話框，在“Lab Convergence based on”欄后面第一欄中選中“Structural”，第二欄中選中“Displacement U”；在“TOLER Tolerance about value”欄后面輸入“0.05”；在“NORM Convergence norm”欄后面下拉菜單選中“L2 norm”；在“MN

53、REF Minimum reference value”后面欄中輸入“1”，單擊“OK”按鈕，這就設(shè)置好了求解時(shí)力矩收斂條件。圖4-31 設(shè)置完力收斂條件對(duì)話框圖4-32 設(shè)置位移收斂條件對(duì)話框圖4-33 設(shè)置好求解收斂條件 最后得到設(shè)置好求解收斂條件，如圖4-33所示。2）邊坡在強(qiáng)度折減系數(shù)F=1時(shí)求解a.求解：Main Menu>Solution>Solve>Current LS，彈出一個(gè)求解選項(xiàng)信息和一個(gè)當(dāng)前求解載荷步對(duì)話框，檢查信息無錯(cuò)誤后，單擊“OK”，開始求解運(yùn)算，直到出現(xiàn)一個(gè)“Solution is done”的提示欄，表示求解結(jié)束。b. 保存求解結(jié)果；Util

54、ity Menu> File> Save as，彈出一個(gè)“Save Database”對(duì)話框，在“Save Database to”下面輸入欄中輸入文件名“F1.db”，單擊“OK”。3）邊坡在強(qiáng)度折減系數(shù)F=1.2時(shí)求解a.折減邊坡強(qiáng)度：首先選擇需要折減單元，再執(zhí)行Main Menu>Solution>Load Step Opts>Other>Change Mat Props>Change Mat Num，彈出一個(gè)“Change Material Number”對(duì)話框，如圖4-34所示，在“New material number”欄后面輸入新材料號(hào)“

55、3”，在“Element no. to be modified”欄后面輸入“ALL”，表示把剛才選定的單元材料該為3號(hào)材料，單擊“OK”按鈕。圖4-34 改變材料號(hào)對(duì)話框b.求解：Main Menu>Solution>Solve>Current LS，彈出一個(gè)求解選項(xiàng)信息和一個(gè)當(dāng)前求解載荷步對(duì)話框，檢查信息無錯(cuò)誤后，單擊“OK”，開始求解運(yùn)算，直到出現(xiàn)一個(gè)“Solution is done”的提示欄，表示求解結(jié)束。c. 保存求解結(jié)果；Utility Menu> File> Save as，彈出一個(gè)“Save Database”對(duì)話框，在“Save Database

56、 to”下面輸入欄中輸入文件名“”，單擊“OK”。同理，依次對(duì)強(qiáng)度折減系數(shù)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進(jìn)行求解，直到求解不收斂為止，并保存各次求解結(jié)果：、F1.8.db、。 當(dāng)強(qiáng)度折減系數(shù)F=3.0時(shí)，求解不收斂，此時(shí)求解迭代力和位移不收斂過程如圖4-36所示。邊坡穩(wěn)定性有限元分析一般采用強(qiáng)度折減方法來求得邊坡安全系數(shù)。強(qiáng)度折減根據(jù)式4-1和式4-2來進(jìn)行折減。求解不收斂是判斷邊坡不穩(wěn)定的一個(gè)準(zhǔn)則。圖4-35 F=2.8 求解收斂時(shí)迭代過程圖圖4-36 F=3.0求解迭代不收斂過程圖 4.3.2.5 后處理 伴隨強(qiáng)度

57、折減系數(shù)的增加，邊坡的塑性應(yīng)變?cè)龃螅苄詤^(qū)也隨之?dāng)U大，當(dāng)塑性區(qū)發(fā)展成一個(gè)貫通區(qū)域，邊坡就不穩(wěn)定，此時(shí)求解也不收斂。與之同時(shí)，邊坡水平位移也變大。因此，主要通過觀察后處理中邊坡塑性應(yīng)變、塑性區(qū)、位移和收斂來判斷邊坡穩(wěn)定性與否。 1）強(qiáng)度折減系數(shù)F=1時(shí)結(jié)果分析 a.讀入強(qiáng)度折減系數(shù)F=1時(shí)結(jié)果數(shù)據(jù)：Utility Menu> resume，彈出一個(gè)“Resume Database”對(duì)話框，選中剛才保存的文件“F1.db”，單擊“OK”按鈕。 b.繪制邊坡變形圖：Main Menu>General Postproc>Plot Results>Deformed Shape，彈

58、出一個(gè)“Plot Deform Shape”對(duì)話框，如圖4-37所示。選中“Def +undeformed”，單擊“OK”按鈕，得到邊坡變形圖，如圖4-39所示。圖4-37 繪制變形圖對(duì)話框4-38 節(jié)點(diǎn)解云圖繪制對(duì)話框 c.顯示邊坡X方向位移云圖：Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，彈出一個(gè)“Contour Nodal Solution Data”對(duì)話框，如圖4-38所示，用鼠標(biāo)依次點(diǎn)擊“Nodal Solution/DOF Solution/ X-Compoment of displacement”，再單擊“OK”按鈕，就得到邊坡X方向位移云圖，如圖4-40所示。此時(shí)，邊坡水平方向最大位移為58.815mm。圖4-39 F=1時(shí)邊坡變形圖圖4-40 F=1時(shí)，邊坡X方向位移云圖 d.顯示邊坡塑性應(yīng)變