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

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

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

3、一定條件下，坡角越大，路基所占面積就越小，反之越大。在山區(qū)， 坡角越大，則路堤所需填方量越少。因此，很有必要對(duì)邊坡穩(wěn)定性進(jìn)行分析，157v1.0可編輯可修改邊坡變形破壞基本原理應(yīng)力分布狀態(tài)邊坡從其形成開(kāi)始，就處于各種應(yīng)力作用(自重應(yīng)力、構(gòu)造應(yīng)力、熱應(yīng)力等)之下。在邊坡的發(fā)展變化過(guò)程中，由于邊坡形態(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)力方向則與坡面近于垂直，并開(kāi)始出現(xiàn)水平方向的剪應(yīng)力，其總趨勢(shì)是由內(nèi)向 外增多，愈近坡腳愈高，向

4、坡內(nèi)逐漸恢復(fù)到原始應(yīng)力狀態(tài)。(2)在坡腳逐漸形成明顯的應(yīng)力集中帶。邊坡愈陡，應(yīng)力集中愈嚴(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)變形的原因。邊坡巖體變形破壞基本形式

5、邊坡在復(fù)雜的內(nèi)外地質(zhì)營(yíng)力作用下形成，又在各種因素作用下變化發(fā)展。所有邊坡都在不斷變形過(guò)程中，通過(guò)變形逐步發(fā)展至破壞。 其基本變形破壞形式主要有：松弛張裂、滑動(dòng)、崩塌、158v1.0可編輯可修改傾倒、蠕動(dòng)和流動(dòng)。影響邊坡穩(wěn)定性的因案影響邊坡穩(wěn)定性的主要因素有：(1)邊坡材料力學(xué)特性參數(shù)：包括彈性模量、泊松比、摩擦角、粘結(jié)力、容重、抗剪強(qiáng)度等參數(shù)。(2)邊坡的幾何尺寸參數(shù)：包括邊坡高度、坡面角和邊坡邊界尺寸以及坡面后方坡體的幾何形狀，即坡體的不連續(xù)面與開(kāi)挖面的坡度及方向之間的幾何關(guān)系，它將確定坡體的各個(gè)部分是否滑動(dòng)或塌落。(3)邊坡外部荷載：包括地震力、重力場(chǎng)、滲流場(chǎng)、地質(zhì)構(gòu)造地應(yīng)力等。邊坡穩(wěn)定性

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

7、動(dòng)方向的代表性剖面， 以表征滑體的基本型態(tài)；將均布力簡(jiǎn)化為集中力， 有時(shí)還將力的作用點(diǎn)簡(jiǎn)化為通過(guò)滑體 重心。極限平衡方法包括以下幾種方法：(1)瑞典圓弧滑動(dòng)法(2)簡(jiǎn)化逼肖普法(3)簡(jiǎn)布普通條分法(4)摩根斯坦-普賴(lài)斯法(5)不平衡推力傳遞法以上各種方法都是假定土體是理想塑性材料，把土條作為一個(gè)剛體，按照極限平衡的原則進(jìn)行力的分析，最大的不同之處在于對(duì)相鄰上條之間的內(nèi)力作何種假定，也就是如何增加已知條件使超靜定問(wèn)題變成靜定問(wèn)題。這些假定的物理意義不一樣，所能滿足的平衡條件也不相同，計(jì)算步驟有繁有簡(jiǎn)，使用時(shí)必須注意他們的適用場(chǎng)合。極限平衡方法關(guān)鍵是對(duì)滑體的休型和滑面的形態(tài)進(jìn)行分析、正確選用滑面的

8、計(jì)算參數(shù)以及正確引用滑體的荷載條件等。因?yàn)闃O限平衡方法完全不考慮土體本身的應(yīng)力-應(yīng)變關(guān)系，不能真 實(shí)地反映邊坡失穩(wěn)時(shí)的應(yīng)力場(chǎng)和位移場(chǎng)，因此而受到質(zhì)疑。數(shù)值分析方法數(shù)值數(shù)值分析方法考慮土體應(yīng)力應(yīng)變關(guān)系，克服了極限平衡方法完全不考慮土體本身的應(yīng)力-應(yīng)變關(guān)系缺點(diǎn)，為邊坡穩(wěn)定分析提供了較為正確和深入的概念。160v1.0可編輯可修改邊坡穩(wěn)定性數(shù)值分析方法主要包含以下幾種方法:(1)有限元法有限單元法是數(shù)值模擬方法在邊坡穩(wěn)定評(píng)價(jià)中應(yīng)用最早的方法，也是目前最廣泛使用的一種數(shù)值方法，可以用來(lái)求解彈性、彈塑性、粘彈塑性、粘塑性等問(wèn)題。目前用有限元法求解邊坡穩(wěn)定主要有兩種方法。a.有限元滑面搜索法：將邊坡體離散

9、為有限單元格，按照施加的荷載及邊界條件進(jìn)行有限元 計(jì)算可得到每個(gè)結(jié)點(diǎn)的應(yīng)力張量。然后假定一個(gè)滑動(dòng)面，用有限元數(shù)據(jù)給出滑動(dòng)面任一點(diǎn)的向 正應(yīng)力和剪應(yīng)力，根據(jù)摩爾一庫(kù)侖準(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年代開(kāi)始自適應(yīng)理論被引入有限元計(jì)算，

10、主導(dǎo)思想是減少前處理工作量和實(shí)現(xiàn) 網(wǎng)格離散的客觀控制?，F(xiàn)已基本建立了一般彈性力學(xué)、流體動(dòng)力學(xué)、滲流分析等領(lǐng)域的平面自 適應(yīng)分析系統(tǒng)，能使計(jì)算較為快速和準(zhǔn)確。(3)離散單元法離散單元法的突出功能是它在反映巖塊之間接觸面的滑移、分離與傾翻等大位移的同時(shí)，又能計(jì)算巖塊內(nèi)部的變形與應(yīng)力分布。因此，任何一種巖體材料都可引入到模型中，例如彈性、161v1.0可編輯可修改粘彈性或斷裂等均可考慮，故該法對(duì)塊狀結(jié)構(gòu)、層狀破裂或一般破裂結(jié)構(gòu)巖體邊坡比較合適。并且，它利用顯式時(shí)間差分法(動(dòng)態(tài)差分法)求解動(dòng)力平衡方程，求解非線性大位移與動(dòng)力問(wèn)題 比較容易。離散元法在模擬過(guò)程中考慮了邊坡失穩(wěn)破壞的動(dòng)態(tài)過(guò)程，允許巖土體存

11、在滑動(dòng)、平移、轉(zhuǎn) 動(dòng)和巖體的斷裂及松散等復(fù)雜過(guò)程，具有宏觀上的不連續(xù)性和單個(gè)巖塊休運(yùn)動(dòng)的隨機(jī)性，可以 較真實(shí)、動(dòng)態(tài)地模擬邊坡在形成和開(kāi)挖過(guò)程中應(yīng)力、位移和狀態(tài)的變化，預(yù)測(cè)邊坡的穩(wěn)定性， 因此在巖質(zhì)高邊坡穩(wěn)定性的研究中得到廣泛的應(yīng)用。(4)拉格朗日元法為了克服有限元等方法不能求解大變形問(wèn)題的缺陷，人們根據(jù)有限差分法的原理，提出了 FLA徵值分析方法。該方法較有限元法能更好地考慮巖土體的不連續(xù)和大變形特性，求解速度 較快。缺點(diǎn)是計(jì)算邊界、單元網(wǎng)格的劃分帶有很大的隨意性。(5)界面元法界面元法是一種基于累積單元變形于界面的界面應(yīng)力元法模型，建立適用于分析不連續(xù)、 非均勻、各向異性和各類(lèi)非線性問(wèn)題、場(chǎng)

12、問(wèn)題，以及能夠完全模擬各類(lèi)錨桿復(fù)雜空間布局和開(kāi) 挖擾動(dòng)的方法。有限元法用于邊坡穩(wěn)定性分析優(yōu)點(diǎn)有限元法考慮了介質(zhì)的變形特征，真實(shí)地反應(yīng)了邊坡的受力狀態(tài)。它可以模擬連續(xù)介質(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)：162v1.0可編輯可修改(1)它考慮了巖體的應(yīng)力-應(yīng)變關(guān)系，求出每一單元的應(yīng)力與變形，反映了巖體真實(shí)工作狀O(2)與極限平衡法相比，不需要進(jìn)行條

13、間力的簡(jiǎn)化，巖體自始至終處于平衡狀態(tài)。(3)不需要像極限平衡法一樣事先假定邊坡的滑動(dòng)面，邊坡的變形特性、塑性區(qū)形成都根據(jù)實(shí)際應(yīng)力應(yīng)變狀態(tài)“自然”形成。(4)若巖體的初始應(yīng)力己知，可以模擬有構(gòu)造應(yīng)力邊坡的受力狀態(tài)。(5)不但能像極限平衡法一樣模擬邊坡的整體破壞，還能模擬邊坡的局部破壞，把邊坡的整體破壞和局部破壞納入統(tǒng)一的體系。(6)可以模擬邊坡的開(kāi)挖過(guò)程，描述和反應(yīng)巖體中存在的節(jié)理裂隙、斷層等構(gòu)造面。鑒于有限元法具有如此多優(yōu)點(diǎn)，本章借助通用有限元軟件ANSY來(lái)實(shí)現(xiàn)對(duì)邊坡穩(wěn)定性分析，用具體的邊坡工程實(shí)例詳細(xì)介紹應(yīng)用 ANSY軟件分析邊坡穩(wěn)定性問(wèn)題。ANSYS邊坡穩(wěn)定性分析步驟ANSYSi坡穩(wěn)定性分

14、析一般分以下幾個(gè)步驟：1、創(chuàng)建物理環(huán)境2、建立模型，劃分網(wǎng)格，對(duì)模型的不同區(qū)域賦予特性3、加邊界條件和載荷4、求解5、后處理(查看計(jì)算結(jié)果)創(chuàng)建物理環(huán)境在定義邊坡穩(wěn)定性分析問(wèn)題的物理環(huán)境時(shí)，進(jìn)入ANSY痢處理器，建立這個(gè)邊坡穩(wěn)定性分 163v1.0可編輯可修改析的數(shù)學(xué)仿真模型。按照以下幾個(gè)步驟來(lái)建立物理環(huán)境:1、設(shè)置GU珠單過(guò)濾如果你希望通過(guò)GUI路徑來(lái)運(yùn)行ANSYS當(dāng)ANSY酸激活后第一件要做的事情就是選擇菜 單路徑：Main MenuPreferences ,執(zhí)行上述命令后，彈出一個(gè)如圖 4-1所示的對(duì)話框出現(xiàn)后， 選才S Structural 。這樣 ANSYS會(huì)根據(jù)你所選擇的參數(shù)來(lái)對(duì)G

15、UI圖形界面進(jìn)行過(guò)濾，選擇Structural以便在進(jìn)行邊坡穩(wěn)定性分析時(shí)過(guò)濾掉一些不必要的菜單及相應(yīng)圖形界面。2、定義分析標(biāo)題（/ TITLE）在進(jìn)行分析前，可以給你所要進(jìn)行的分析起一個(gè)能夠代表所分析內(nèi)容的標(biāo)題，比如“Slopestability Analysis ”,以便能夠從標(biāo)題上與其他相似物理幾何模型區(qū)別。用下列方法定義分析 標(biāo)題。命令：/TITLEGUI: Utility MenuFileChange Title3、說(shuō)明單元類(lèi)型及其選項(xiàng)（KEYOP選項(xiàng)）與ANSYS勺其他分析一樣，也要進(jìn)行相應(yīng)的單元選擇。ANSY歌件提供了 100種以上的單元類(lèi)型，可以用來(lái)模擬工程中的各種結(jié)構(gòu)和材料，各

16、種不同的單元組合在一起，成為具體的物 理問(wèn)題的抽象模型。例如，不同材料屬性的邊坡土體用PLANE8卻元來(lái)模擬。大多數(shù)單元類(lèi)型都有關(guān)鍵選項(xiàng)（KEYOPTS這些選項(xiàng)用以修正單元特性。例如， PLANE82有如下KEYOPTSKEYOPT（2）包含或抑制過(guò)大位移設(shè)置KEYOPT（3）平面應(yīng)力、軸對(duì)稱(chēng)、平面應(yīng)變或考慮厚度的平面應(yīng)力設(shè)置164v1.0可編輯可修改KEYOPT（9）用戶子程序初始應(yīng)力設(shè)置設(shè)置單元以及其關(guān)鍵選項(xiàng)的方式如下：命令：ETKEYOPTGUI: Main Menu Preprocessor Element Type Add/Edit/Delete圖4-1 GUI圖形界面過(guò)濾4定義單位

17、結(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）,則導(dǎo)出力的單位是 kg? m/s2 （相當(dāng)于牛頓N）,材料的彈性模量單位是 kg/m? s2 （相 當(dāng)于帕Pa）。命令：/UNITS165v1.0可編輯可修改5、定義材料屬性大多數(shù)單元類(lèi)型在進(jìn)行程序分析時(shí)都需要指定材料特性，ANSY鈉序可方便地定義各種材料的特性，如結(jié)構(gòu)材料屬性參數(shù)、熱性能參數(shù)、流體性能參數(shù)和電磁性能參數(shù)等。ANSYS序可定義的材料特性有以下三種：(1)線性或非線性。(2)各向同性、正交異性或非彈性。(3

18、)隨溫度變化或不隨溫度變化。因?yàn)榉治龅倪吰履Ｐ筒捎美硐霃椝苄阅Ｐ?D-P模型)，因此邊坡穩(wěn)定性分析中需要定義邊坡中不同土體的材料屬性：容重、彈性模量、泊松比、凝聚力以及摩擦角。命令：MPGUI: Main MenuPreprocessorMaterial PropsMaterial Models或 Main MenuSolutionLoad Step OptsOtherChange Mat PropsMaterial Models進(jìn)行邊坡穩(wěn)定性分析計(jì)算時(shí)，采用強(qiáng)度折減法來(lái)實(shí)現(xiàn)。首先選取初始折減系數(shù)F,然后對(duì)4-1 和式 4-2。CFtanF邊坡土體材料強(qiáng)度系數(shù)進(jìn)行折減，折減后凝聚力以及摩擦角分

19、別見(jiàn)式C(4-1)tan(4-2)166v1.0可編輯可修改C和為邊坡土體的初始凝聚力和摩擦角。C和 進(jìn)行折減，輸入邊坡模型計(jì)算，若收斂，則此時(shí)邊坡是穩(wěn)定的；繼續(xù)增大折減系數(shù)F,直到程序恰好不收斂，此時(shí)的折減系數(shù)即為穩(wěn)定或安全系數(shù)。建立模型和劃分網(wǎng)格創(chuàng)建好物理環(huán)境，就可以建立模型。在進(jìn)行邊坡穩(wěn)定性分析時(shí)，需要建立模擬邊坡土體的 PLANE8即元。在建立好的模型各個(gè)區(qū)域內(nèi)指定特性（單元類(lèi)型、選項(xiàng)、實(shí)常數(shù)和材料性質(zhì)等） 以后，就可以劃分有限元網(wǎng)格了。通過(guò)GUI為模型中的各區(qū)賦予特性：1、選擇 Main Menu Preprocessor Meshing Mesh Attributes Picked

20、 Areas2、點(diǎn)擊模型中要選定的區(qū)域。3、在對(duì)話框中為所選定的區(qū)域說(shuō)明材料號(hào)、實(shí)常數(shù)號(hào)、單元類(lèi)型號(hào)和單元坐標(biāo)系號(hào)。4、重復(fù)以上三個(gè)步驟，直至處理完所有區(qū)域。通過(guò)命令為模型中的各區(qū)賦予特性：ASEL （選擇模型區(qū)域）MAT （說(shuō)明材料號(hào)）REAL （說(shuō)明實(shí)常數(shù)組號(hào)）TYPE （指定單元類(lèi)型號(hào)）ESYS（說(shuō)明單元坐標(biāo)系號(hào)）167v1.0可編輯可修改施加約束和荷載在施加邊界條件和荷載時(shí)， 既可以給實(shí)體模型（關(guān)鍵點(diǎn)、線、面）也可以給有限元模型（節(jié) 點(diǎn)和單元）施加邊界條件和荷載。在求解時(shí)，ANSYS程序會(huì)自動(dòng)將加到實(shí)體模型上的邊界條件和載荷轉(zhuǎn)遞到有限元模型上。邊坡穩(wěn)定性分析中，主要是給邊坡兩側(cè)和底部施

21、加自由度約束。命令：D施加荷載包括自重荷載以及邊坡開(kāi)挖荷載。求解接著就可以進(jìn)行求解，ANSYS程序根據(jù)現(xiàn)有選項(xiàng)的設(shè)置，從數(shù)據(jù)庫(kù)獲取模型和載荷信息并進(jìn)行計(jì)算求解，將結(jié)果數(shù)據(jù)寫(xiě)入到結(jié)果文件和數(shù)據(jù)庫(kù)中。命令：SOLVEGUI: Main MenuSolutionSolve Current LS后處理后處理的目的是以圖和表的形式描述計(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ā)生了破壞。通 過(guò)研究位移、應(yīng)變和塑性區(qū)域，來(lái)綜合判斷邊坡的穩(wěn)定性。命

22、令:/POST1GUI: Main Menu General Postproc168v1.0可編輯可修改補(bǔ)充說(shuō)明邊坡的失穩(wěn)破壞定義有很多種，對(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)。采用

23、彈塑性模型更能反應(yīng)巖體的實(shí)際工作狀態(tài)。(2)巖體所承受的荷載超過(guò)材料強(qiáng)度時(shí)，就會(huì)出現(xiàn)明顯的滑移破壞面。因此，彈塑性計(jì)算不 需要假定破壞面的形狀和位置，破壞面根據(jù)剪應(yīng)力強(qiáng)度理論自動(dòng)形成。當(dāng)整個(gè)邊坡破壞時(shí)，就 會(huì)出現(xiàn)明顯的塑性區(qū)。(3 )能綜合考慮邊坡的局部失穩(wěn)和整體失穩(wěn)破壞。169v1.0可編輯可修改ANSYS邊坡穩(wěn)定性實(shí)例分析實(shí)例描述米米邊坡圍巖2彈塑性材料邊坡圍巖1彈性材料圖4-2邊坡模型邊坡實(shí)例選取國(guó)內(nèi)某礦，該邊坡考慮彈性和塑性?xún)煞N材料，邊坡尺寸如圖4-2所示。分析目的是對(duì)該邊坡進(jìn)行穩(wěn)定性計(jì)算分析，以判斷其穩(wěn)定性和計(jì)算出安全系數(shù)，該邊坡圍巖材料屬性 見(jiàn)表4-1。表4-1邊坡模型圍巖參數(shù)類(lèi)別

24、彈性模量/Gpa泊松比v容重/ kN / m3內(nèi)聚力/Mpa摩擦角（。）圍巖2 （彈塑性）30250042圍巖1 （彈性）312700-對(duì)于像邊坡這樣縱向很長(zhǎng)的實(shí)體，計(jì)算模型可以簡(jiǎn)化為平面應(yīng)變問(wèn)題。假定邊坡所承受的170v1.0可編輯可修改外力不隨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所示。采用雙層模型，模型上部為理想彈塑性材料，下部為彈性材料，左右邊界水平位移為零

25、，下邊界豎向位移為零。雙層模型考慮土體的彈塑性變形，其塑性區(qū)的發(fā)展，應(yīng)力的分布更符合實(shí)際情況?？紤]雙層模型，塑性區(qū)下部的單元可以產(chǎn)生一定的垂直變形和水平變形，基本消除了由于邊界效應(yīng)在邊坡下部出現(xiàn)的塑性區(qū)，更好地模擬了邊坡的變形和塑性區(qū)的發(fā)展。GUI操作方法創(chuàng)建物理環(huán)境1）在【開(kāi)始】菜單中依次選取【所有程序】/【】/ ANSYProduct Launcher，得到Product Launcher對(duì)話框。2 ）選中 【File Management,在“Working Directory ”欄輸入工作目錄“D:ansysexample4-1 ”,在Job Name”欄輸入文件名“Slope”。3）

26、單擊“ RUN按鈕，進(jìn)入的GUI操作界面。4）過(guò)濾圖形界面：Main Menu Preferences ,彈出Preferences for GUI Filtering ”對(duì)話框，選中“ Structural ”來(lái)對(duì)后面的分析進(jìn)行菜單及相應(yīng)的圖形界面過(guò)濾。5）定義工作標(biāo)題：Utility Menu File Change Title ,在彈出的對(duì)話框中輸入Slopestability Analysis ,單擊OK,如圖 4-3。171v1.0可編輯可修改圖4-3定義工作標(biāo)題6）定義單元類(lèi)型:a.定義 PLANE8即元：Main Menu Preprocessor Element Type Ad

27、d/Edit/Delete ,彈 出一個(gè)單元類(lèi)型對(duì)話框，單擊“ Add”按鈕。弓t出如圖4-4所示對(duì)話框。在該對(duì)話框左面滾動(dòng)欄中選擇“Solid ”，在右邊的滾動(dòng)欄中選擇 “Quad8node 82，單擊“Apply，就定義了 “PLANE82單元。圖4-4定義PLANE82I元對(duì)話框b. 設(shè)定 PLANE8邪元選項(xiàng)：Main Menu Preprocessor Element Type Add/Edit/Delete ,彈出一個(gè)單元類(lèi)型對(duì)話框，選中“Type 2 PLANE82,單擊“Options ”按鈕，彈出一個(gè)PLANE82element Type options ”對(duì)話框，如圖4-5

28、所示。在Element behavior K3 ”欄后面的下拉菜 單中選取“ Plane strain ”，其它欄后面的下拉菜單采用ANSY激認(rèn)設(shè)置就可以，單擊“ OK按 172v1.0可編輯可修改鈕。PLANE82 element type optionsOptions for PLANE82, Element Type Ret No. 1 TOC o 1-5 h z Ele merit behaviorK3Extr曰 e Ie ment outputK5Extrw surface outputK6OKCancel圖4-5 PLANE82單元庫(kù)類(lèi)型選項(xiàng)對(duì)話框通過(guò)設(shè)置PLANE8理元選項(xiàng) K3

29、為 Plane strain ”來(lái)設(shè)定本實(shí)例分析采取平面應(yīng) 變模型進(jìn)行分析。因?yàn)檫吰率强v向很長(zhǎng)的實(shí)體，故計(jì)算模型可以簡(jiǎn)化為平面應(yīng)變問(wèn)題。8節(jié)點(diǎn)PLANE8邪元每個(gè)節(jié)點(diǎn)有 UX和UY兩個(gè)自由度，比4節(jié)點(diǎn)PLANE4卻元具有更 高的精確性，對(duì)不規(guī)則網(wǎng)格適應(yīng)性更強(qiáng)。7）定義材料屬性a.定義邊坡圍巖1 材料屬性：Main Menu PreprocessorMaterial PropsMaterialModels,彈出Define Material Model Behavior ”對(duì)話框，如圖 4-6 所示。173v1.0可編輯可修改圖4-6定義材料本構(gòu)模型對(duì)話框在圖4-6中右邊欄中連續(xù)雙擊Struct

30、ural Linear ElasticIsotropic ”后，又彈出如圖 4-7 所示“Linear Isotropic Properties for Material Number 1”對(duì)話框，在該對(duì)話框中“EX后面的輸入欄輸入“ 3E10”，在“PRXY 后面的輸入欄輸入，單擊“ OK。再在選中“Density ”并雙擊，彈出如圖 4-8 所示“Density for Material Number!對(duì)話框，在DENS后面的欄中輸入邊坡土體材料的密度“2500”，單擊“ OK按鈕。再次在圖 4-6 中右邊的欄中連續(xù)雙擊“ StructuralNonlinear InelasticNon-

31、metalplasticitydrucker-prager ”后，又彈出一個(gè)如圖 4-9所示對(duì)話框。在Cohesion”欄添入邊坡圍巖材料1的內(nèi)聚力，在Fric Angle ”欄添入邊坡內(nèi)摩擦角“ 42”，單擊“ OK按鈕。174v1.0可編輯可修改圖4-7線彈性材料模型對(duì)話框Etensity for Material Number 1圖4-8材料密度輸入對(duì)話框圖4-9定義邊坡材料1DP模型對(duì)話框b.定義邊坡圍巖2材料屬性：在圖4-6對(duì)話框中，單擊“ Material New Model”，彈出一個(gè) DeMe Material ID ”對(duì)話框，在“ID”欄后面輸入材料編號(hào)“ 2”，單擊“ OK

32、按鈕。彈出一個(gè)定義材料模型對(duì)話框?qū)υ捒?，選中“ Material Model Number 2，和定義邊坡圍巖 1材料一樣，在右邊的欄中連續(xù)雙擊 Structural Linear日asticIsotropic 后，又彈出一個(gè) a Linear Isotropic Properties for Material Number 2”對(duì)話框，在該對(duì)話框中EX后175v1.0可編輯可修改面的輸入欄輸入，在“ PRXY后面的輸入欄輸入，單擊“ OK。再選中“ Densly ”并雙擊，彈出一個(gè) Density for Material Number 2”對(duì)話框，在“ DENS后面的欄中輸入隧道圍巖材料的

33、密度“ 2700”，再單擊“ OK按鈕，彈出一個(gè)定義材料模型對(duì)話框。c.復(fù)制邊坡圍巖1材料性質(zhì)：在圖4-6對(duì)話框中，用鼠標(biāo)點(diǎn)擊“Editcopy .”,彈出一個(gè)“Copy Material Model ”對(duì)話框，如圖 4-10 所示。在from Material number ”欄后面的下拉菜單中選取1,在TOMaterial number”欄后面輸入“3,單擊Apply”按鈕。又彈出如土 4-10所示對(duì)話框，然后依次在“ TOMaterial number”欄后面輸入“ 4”、“5”、“6”、“7”、“8” “9”、“10”，“11”、“12”、“13”，每輸入一個(gè)數(shù)，就單擊 “ Apply

34、” 按鈕一次。GlCopy Material ModelCopy Material Modelfrom Material numberOKCancelto Material number圖4-10復(fù)制本構(gòu)模型對(duì)話框最后彳#到10個(gè)復(fù)制圍巖1的邊坡材料本構(gòu)模型，如圖 4-11所示。176v1.0可編輯可修改d.定義10個(gè)強(qiáng)度折減后材料本構(gòu)模型：首先定義強(qiáng)度折減系數(shù)5=后邊坡圍巖材料模型,在圖4-11對(duì)話框中，在鼠標(biāo)依次雙擊“ Material Model Number 3/Drucper-Prager 。彈出一 個(gè)Drucker- Prager Material Number 3,如圖 4-12

35、 所示，在Cohesion”欄添入強(qiáng)度折減系數(shù)5=后邊坡圍巖材料1的內(nèi)聚力“：在Fric Angle ”欄添入折減后邊坡內(nèi)摩擦角， 單擊“ OK按鈕。177v1.0可編輯可修改用相同方法定義強(qiáng)度折減系數(shù)分別為：F=、F=、F=、F=、F=、F=、F=、F=、F=的邊坡圍巖材料本構(gòu)模型。定義強(qiáng)度折減后本構(gòu)模型目的是為了分析邊坡穩(wěn)定性。強(qiáng)度折減就是降低內(nèi)聚力和摩擦角，根據(jù)式 4-1和式4-2進(jìn)行折減。建立模型和劃分網(wǎng)格1)創(chuàng)建邊坡線模型a.輸入關(guān)鍵點(diǎn)：Main Menu Preprocessor Modeling CreateKeypointsIn Active CS , 彈出Creae Keyp

36、oints in Active Cooedinate System”對(duì)話框，如圖 4-13 所示。在NPTkeypoint number”欄后面輸入1,在X, Y, Z Location in active CS ”欄后面輸入“(0, 0, 0)”, 單擊“Apply”按鈕，這樣就創(chuàng)建了關(guān)鍵點(diǎn) 1。再依次重復(fù)在“NPT keypoint number ”欄后面 輸入“2、3、4、5、6、7、8、9”，在對(duì)應(yīng)“X,Y,Z Location in activeC6欄后面輸入 “(-800,0, 0)、(-800, -800 , 0)、(-800, -1200 , 0)、(1200, -1200 ,

37、 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 MenuPreprocessor Modeling CreateLinesStraight line ,彈出Creae straight lines ”對(duì)話框，用鼠標(biāo)依次點(diǎn)擊關(guān)鍵點(diǎn)1、2,單擊“Apply”按鈕， 178v1.0可編輯可修改這樣就創(chuàng)建了直接L1,同樣分別連接關(guān)鍵點(diǎn)“ 2、3”，“3、4”，“4、5”，“5、6”，“6、7”，“7、4”，“7、8”，“8、9”，

38、“9、1”，“9、2”，最后單擊“ OK按鈕，就得到邊坡線模型，如圖 4-14所示。ANSYSHAY El ZDD71IMES圖4-14 邊坡線模型3）創(chuàng)建邊坡面模型a.打開(kāi)面編號(hào)顯示：Utility Menu PlotCtrls Numbering ,彈出Plot NumberingControls ”對(duì)話框，如圖4-15所示。選中Aares Numbers” 選項(xiàng)，后面白文字由off 變 為“on” ,單擊“OK關(guān)閉窗口。179v1.0可編輯可修改QpidtContrulsyFMLM Plat Fojirtenng Ccnffr-cfeKF Kfirnhi： nurrhFjn：r offL

39、JhE Lire rumtErsr orrAREAnumbenF pnM2LIU 也 Lime nurriber orrMJLte r-Jude IlJtntwrBr offElem / Nttrfci rumaer rq|No iximbErnqTABh T able Narnsar orr5VAL Numert ccnCDir idu尋r off40-1 iJurmtiFring skcvm wflh|cdlor5 B nunteraiRuulot3L0T現(xiàn)Lt/u氏他加弁Ok I 七 pM I Canral I +b圖4-15打開(kāi)面編號(hào)對(duì)話框圖4-16邊坡面模型b.創(chuàng)建邊坡面模型：Mai

40、n Menu Preprocessor Modeling CreateAreas Arbitrary 180v1.0可編輯可修改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”按鈕，就生成了邊坡開(kāi)挖掉區(qū)域面積A3。最后得到邊坡模型的面模型，如圖4-16所示。4 ）劃分邊坡圍巖2

41、單元網(wǎng)格a.給邊坡圍巖 2 賦予材料特性： Main Menu PreprocessorMeshing MeshTool ,彈出“MeshTool”對(duì)話框，如圖4-17所示。在Element Attributes ”后面的下拉式選擇欄中選擇 Areas,按Set按鈕，彈出一個(gè)Areas Attributes ”面拾取框，在圖形界面上拾取邊坡圍巖2區(qū)域，單擊拾取框上的OK按鈕，又彈出一個(gè)如圖 4-18所示的“Areas Attributes 對(duì)話框，在Material number”后面的下拉式選擇欄中選取2,在Element type number 后面的下拉式選擇欄中選取“ 2 PLANE8

42、2,單擊“ Apply ”。181v1.0可編輯可修改圖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è) Hement Sizes on Picked Lines ”,對(duì)話框，如圖 4-19 所示，在No of element division ”欄后面輸入5”, 單擊“Apply按鈕，再選擇線L4和L11,又彈出圖4-19對(duì)話框，在Noof element division 欄后面輸入“ 26”，單擊“ OK按鈕。c.

43、劃分單元網(wǎng)格：在圖4-17網(wǎng)格劃分工具欄中單擊 “Mesh按鈕，彈出一個(gè)拾取面積對(duì)話 框，拾取面積A1,單擊拾取框上的“ OK按鈕，生成邊坡圍巖 2單元網(wǎng)格。182v1.0可編輯可修改tlemerrt Sizes on Picked LinesSIZEElnnnt edge 回時(shí)NDVNo af eterrenit divisionsJESEE Eferrent m 0sls on picked 帕。6(fIDLV b ursed orJy If SEE is 匕 1 ol 段廣口)KYTO7 SjZE.TJDlV CoTi be chmn湖p YeEPACE Spacing ratioAIN

44、I2S1Z DNl&iQin arc ijjacreeis)(u*w &NG曰【工 uBiS if ruuirnfaBT uf dlrkluht (NDtV) *2 tIheet t ydyt? fesiigti三工壬)nie blUik cr htu)Cfeflir attached res nd wLiittss廠 pjqQK |-pply |CcellH葉圖4-19設(shè)置網(wǎng)格份數(shù)對(duì)話框5）劃分邊坡圍巖1單元網(wǎng)格BetGontroIs聲PickUjipAukS3.vxg;l廣EonPoivannrcivileLLecpC umIf ftx 1=，m=NG 1ITiniuiAn.,XLin.r

45、*5 Lit o t 工 ft. me1LMa.K f Xixc圖4-20選取線對(duì)話框圖4-21設(shè)置網(wǎng)格分?jǐn)?shù)對(duì)話框183v1.0可編輯可修改設(shè) 置 網(wǎng) 格 份 數(shù)： Main MenuPreprocessorMeshingSize CntrlsManualSizeLayersPicked Lines ,彈出一個(gè)Set 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 ”欄后面輸入

46、“10”, 單擊“OK按鈕。相同方法設(shè)置線L8和L10分割份數(shù)為16;設(shè)置線L7和L9線的分割份數(shù)為12。給邊坡圍巖 1 賦予材料特性：Main Menu Preprocessor Meshing MeshTool ,彈出“MeshTool”對(duì)話框，如圖4-17所示。在Element Attributes ”后面的下拉式選擇欄中選擇 Areas,按Set按鈕，彈出一個(gè)Areas Attributes ”面拾取框，在圖形界面上拾取面 A2和A3,單擊拾取框上的“OK按鈕，又彈出一個(gè) Areas Attributes ”對(duì)話框，在aMaterial number”后面的下拉式選擇欄中選取1,在El

47、ement type number ” 后面的下拉式選擇欄 中選取“ 2 PLANE82,單擊“ Applyc.劃分單元網(wǎng)格：在圖3-97網(wǎng)格劃分工具欄中單擊“ Mesh按鈕，彈出一個(gè)拾取面積對(duì)話 框，拾取圍巖，單擊拾取框上的“OK按鈕，生成邊坡圍巖1單元網(wǎng)格。最后得到邊坡模型單元網(wǎng)格，如圖 4-22所示。184v1.0可編輯可修改百1 EH4MT5血aLyn上事圖4-22邊坡模型單元網(wǎng)格施加約束和荷載）給邊坡模型施加約束a.給邊坡模型兩邊施加約束：執(zhí)行 Main MenuSolutionDefine LoadsApplyStructural Displacementon Nodes,彈出在節(jié)

48、點(diǎn)上施加位移約束對(duì)話框，用鼠標(biāo)選取隧道模型兩側(cè)邊界上 所有節(jié)點(diǎn)，單擊“ OK按鈕。彈出“ Apply U,ROTon Nodes”對(duì)話框，如圖4-23所示，在“ DOFS to be constrained ”欄后面中選取UX,在Apply as”欄后面的下拉菜單中選取Constant value,在Displacement value 欄后面輸入0值，然后單擊“OK按鈕。185v1.0可編輯可修改圖4-23給模型兩側(cè)施加位移約束對(duì)話框圖4-24給模型底部施加位移約束對(duì)話框b.給模型底部施加約束：執(zhí)行Main MenuSolutionDefine Loads ApplyStructuralDi

49、splacementon Nodes,彈出在節(jié)點(diǎn)上施加位移約束對(duì)話框，用鼠標(biāo)選取隧道模型底部邊界上所有節(jié)點(diǎn)，單擊“OK按鈕。彈出圖4-24所示對(duì)話框，在DOFSto be constrained ”欄后面 中選取“UX UY,在“Apply as”欄后面的下拉菜單中選取 Constant value ,在Displacementvalue ”欄后面輸入“ 0”值，然后單擊“ OK按鈕。節(jié)點(diǎn)選擇，可以先選擇節(jié)點(diǎn)上線，再選擇附在線上的節(jié)點(diǎn)。) 施 加 重 力 加 速 度： Main MenuSolutionDefineLoadsApplyStructuralInertiaGravity ,彈出Ap

50、ply(Gravitational)Acceleration ”對(duì)話框，如圖4-25所示。只需在“Global Cartesian Y-comp ”欄后面輸入重力加速度值就可以，單擊“ OK按鈕，就完成了重力加速度的施加。186v1.0可編輯可修改Apply (CSravitaticinal) AccelerationACEL apply (Gravitational) AccelerationACELX Global Cartesian X-comp口ACELY Global Carteriar Y-compghgACELZ Global Cartesian 2-cqmpoOK |匚 #i

51、. El JH.lp |圖4-25施加重力加速度對(duì)話框這時(shí)就可以得到施加約束和重力加速度后隧道有限元模型，如圖4-25所示。圖4-25施加約束和重力荷載后邊坡模型求解1）求解設(shè)置 187v1.0可編輯可修改a.指定求解類(lèi)型：Main MenuSolutionAnalysis TypeNewAnalysis ,彈出一個(gè)如圖 4-25所示對(duì)話框，在Type of analysis ”欄后面選中“Static ,單擊OK按鈕。ANTYPF Typs of anAl/si?作 Stade ；r rvticuiC Harmcni: TransientSpEctrur C Eigen BucklingSd

52、bstructLirngzLMSCK ICanoe I IHelp I圖4-25指定求解類(lèi)型對(duì)話框b.設(shè)置載荷步： Main MenuPreprocessorLoadsAnalysis TypeSoln 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按鈕。188v1.0可編輯可修改圖4-26設(shè)置載荷步對(duì)話框c.設(shè)置線性搜索：Mai

53、n MenuSolutionAnalysis TypeSoln Controls，彈出一個(gè)SolutionControls ”對(duì)話框，用鼠標(biāo)單擊Nonlinear ”選項(xiàng)，如圖4-27所示，在Line search”欄后面下拉菜單選中“ ON,單擊“ OK。189v1.0可編輯可修改圖4-27設(shè)置線性搜索對(duì)話框d.設(shè)定牛頓-拉普森選項(xiàng)：Main MenuSolutionAnalysis TypeAnalysis Options ,彈出 一個(gè)Static or Steady-static Analysis ”對(duì)話框，如圖 4-28 所示，在New-Raphsonoption 欄后面下拉菜單選中“

54、 Full N-R ”，單擊“ OK按鈕。IA(T Sternly State Ar alvslsrjmirear Cpocre：HL里DIM _arqe deta-iti elT=ct=KRQT | Hav-tor-RittLCn optcriAdaptive dm centr油nH rf i怕蟲(chóng)二ILirr OptionsLUM=M1 USe bmpeJI加口r口感n口EQSLV Bqurtir iMiverToler src3.iWvfeH- raid for al except Frontal and 5pa-se So 到efePrcgnm chogen |Mjltpliar -Ho

55、n hf fw F( HcsnuiUari CG【Ff毛nskM s舊ptwl-r m臼日 or y fir Fcunurti匚n CGhE酬配Mwnoq Savs -n cffnlri on X Errmncttlrn OGP onv3idor*y Frortal bpfse andpcg seversr-TTpFTRFrSTE S ffr環(huán),Dr pnp tii-Ka回代RJuIli： imCEOM.Cf-J Hull ea： E&TtF.OMtoffgt TeirpEranre ctffirencE-aUQkJts iercand詢(xún)口 cfooa temp -U qkcmncml圖4-2

56、8設(shè)定牛頓-拉普森選項(xiàng)e.打開(kāi)大位移求解： Main MenuPreprocessorLoadsAnalysis TypeSolnControlsBasic ,得到圖4-26所示對(duì)話框，在Analysis Options ”欄后面下拉菜單選中大位190v1.0可編輯可修改移“Large Displacement static ,單擊OK 按鈕。f. 設(shè) 置收 斂條 件： Main MenuPreprocessorLoadsLoadStepOptsNonlinearConvergence Crit , 彈出個(gè) Default Nonlinear Convergence Criteria , 如圖

57、4-29所示。圖中顯示 ANSY欲認(rèn)的收斂條件：分別設(shè)置了力和力矩的收斂條件。圖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 ” 欄后面輸入；在 “ NORM Convergence n

58、orm”欄后面下拉菜單選中L2 norm;在MNREF Minimum reference value ” 后面欄中輸入，單擊“ OK按鈕，這就設(shè)置好了求解時(shí)力收斂條件。191v1.0可編輯可修改圖4-30設(shè)置力收斂條件對(duì)話框單擊圖4-30對(duì)話框中“OK按鈕后，彈出一個(gè)如圖4-31所示的對(duì)話框，單擊“Add”按鈕。彈出一個(gè)如圖4-31所示對(duì)話框，在Lab Convergence based on 欄后面第一欄中選中aStructural ,第二欄中選中“ Displacement U 在TOLER Tolerance about value 欄后面輸入；在NORConvergence norm

59、”欄后面下拉菜單選中L2 norm;在MNREMinimumreference value ”后面欄中輸入“ 1”，單擊“ OK按鈕，這就設(shè)置好了求解時(shí)力矩收斂條件。192v1.0可編輯可修改圖4-31設(shè)置完力收斂條件對(duì)話框足fsonlnear CorwergencoCriterEWirroL卜inn re汗口EimzEnrw 二Lg-月 |_ Cunvergpnce 曰 ba&3zl onVALUE R日展日門(mén)午 F -日pf LabTOL 國(guó) Toltrce 如 UL VALUE網(wǎng)口RM Ctrwflrprna noririNl： hPG Mim el m rert-renre value

60、&珀d Wy if VALUE憐Ulsnk,開(kāi)麗網(wǎng)t歸no恂產(chǎn)團(tuán)G曰n%8Q)Ues圖4-32設(shè)置位移收斂條件對(duì)話框圖4-33設(shè)置好求解收斂條件最后得到設(shè)置好求解收斂條件，如圖4-33所示。2）邊坡在強(qiáng)度折減系數(shù)F=1時(shí)求解a.求解：Main MenuSolutionSolveCurrentLS,彈出一個(gè)求解選項(xiàng)信息和一個(gè)當(dāng)前求解載荷步對(duì)話框，檢查信息無(wú)錯(cuò)誤后，單擊“OK,開(kāi)始求解運(yùn)算，直到出現(xiàn)一個(gè)Solution is done”193v1.0可編輯可修改的提示欄，表示求解結(jié)束。b.保存求解結(jié)果；Utility Menu File Save as ,彈出一個(gè)Save Database 對(duì)話框