1、abaqus地應(yīng)力平衡先說(shuō)為什么要施加地應(yīng)力：1、我們所建立的幾何模型一般和工程實(shí)際情況或尺寸相對(duì)應(yīng)、相一致，比方邊坡幾何模型和實(shí)際邊坡尺寸一致，但我們可以夸張一點(diǎn)想像，實(shí)際邊坡應(yīng)是由一個(gè)更大一點(diǎn)或更高一點(diǎn)的不受重力的初始邊坡在n年前突然受重力和類(lèi)似目前的邊界條件作用下逐漸形成了今天的尺寸大小，n年前受重力和類(lèi)似目前的邊界條件作用之前邊坡的尺寸大小，我們不得而知，如果能準(zhǔn)確知曉，我們就可以建立一個(gè)那時(shí)的幾何模型，再施加重力和邊界條件進(jìn)行計(jì)算，變形后形狀和現(xiàn)狀邊坡形狀一致，其內(nèi)力也就是初始應(yīng)力場(chǎng)或地應(yīng)力，就不用專(zhuān)門(mén)去施加地應(yīng)力了，但問(wèn)題是我們不能知曉邊坡受力前的形狀尺寸，我們現(xiàn)在的幾何模型就是邊

2、坡現(xiàn)在的實(shí)際尺寸，受力后將會(huì)變成一個(gè)更小的或與現(xiàn)狀不一致的邊坡，這不符合我們模擬現(xiàn)狀邊坡的目的。如果我們知道現(xiàn)狀邊坡的內(nèi)力，將其提取出來(lái)作為幾何模型的內(nèi)力，再和外力重力平衡，那么我們建立的模型才能算和實(shí)際模型一致。真實(shí)地知道現(xiàn)狀邊坡的內(nèi)力是很難的，我們采取的方法是，用我們所建立的幾何模型施加和實(shí)際模型一致的重力和邊界條件進(jìn)行計(jì)算，得到變形后或變得更小或與現(xiàn)狀邊坡不完全一致的邊坡內(nèi)力近似的作為現(xiàn)狀邊坡的內(nèi)力，并重新將其施加于與現(xiàn)狀邊坡一致的幾何模型，再施加重力當(dāng)然邊界條件也應(yīng)根本一致以平衡，這樣才算建立了與現(xiàn)狀模型根本一致的模型，其下的計(jì)算才成為可能。這就是所謂“地應(yīng)力平衡的含義、目的、作用。2

3、. 地應(yīng)力平衡中的外力和內(nèi)力的問(wèn)題，地應(yīng)力平衡中，顯然，重力是外力，應(yīng)力場(chǎng)是內(nèi)力，僅有外力重力，沒(méi)有內(nèi)力是不可能的，同樣，僅有內(nèi)力專(zhuān)指初始應(yīng)力場(chǎng)而不受重力也是不可能的，否那么，整個(gè)體系的力不會(huì)平衡。這就是為什么我們將提取出的內(nèi)力施加于幾何模型后必須再施加重力的原因。為的是內(nèi)力和外力平衡。) q0 F3 q6 H1 O#3. 地應(yīng)力場(chǎng)的方向問(wèn)題，有網(wǎng)友在論壇里問(wèn)，既然重力是向下，為與重力平衡，那應(yīng)力場(chǎng)的方向是不是向上呢，這同樣是我開(kāi)始接觸abaqus的疑問(wèn)，相信很初學(xué)者也有這樣的疑問(wèn)，我的理解是內(nèi)力是沒(méi)有向上、向下或者向其它方向的概念的，內(nèi)力只有拉力或壓力或剪力之分，其方向也

4、按是拉是壓是順時(shí)針或逆時(shí)針而分，內(nèi)力往往都是成對(duì)出現(xiàn)，如地應(yīng)力場(chǎng)中的應(yīng)力以壓應(yīng)力為主，取一個(gè)微元，那么壓應(yīng)力同時(shí)出現(xiàn)在向下和向上，你能說(shuō)地應(yīng)力就是向上，與重力反向嗎？aba中初始地應(yīng)力場(chǎng)平衡一般在外表水平的情況下僅僅和密度相關(guān)，密度一樣的話(huà)平衡的結(jié)果很好，別的參數(shù)改變之后經(jīng)過(guò)計(jì)算，差異很小。外表不水平的情況那么最好通過(guò)文件導(dǎo)入初應(yīng)力的情況進(jìn)行平衡。一般來(lái)講，外表不平的時(shí)候有很多因素造成誤差很大。cdstudio 版主的總結(jié)： :/forum.simwe /viewthread.php?tid=908486&extra=page%3D1%26amp%3Bfilter%3Dtyp

5、e%26amp%3Btypeid%3D470地應(yīng)力平衡的正確時(shí)間點(diǎn):所謂地應(yīng)力平衡是指, 當(dāng)我們建任何東西或挖任何東西之前, 地表的位移都是零, 但是土體的應(yīng)力卻存在, 這種無(wú)位移但有應(yīng)力的時(shí)間點(diǎn)叫地應(yīng)力平衡. 那么正確的時(shí)間點(diǎn)應(yīng)該是什么呢, 就是不管土原來(lái)是什么樣的, 高山, 河流, 丘陵, 平原都可以, 關(guān)鍵是在我們對(duì)它做任何擾動(dòng)之前. 有些人問(wèn)樁土作用如何平衡? 這個(gè)問(wèn)題首先就是錯(cuò)的, 因?yàn)榈貞?yīng)力平衡的時(shí)候是不能有樁的. 一旦有樁就跨過(guò)了地應(yīng)力平衡這個(gè)時(shí)間點(diǎn), 數(shù)值模擬要忠實(shí)于實(shí)際情況, 不能想當(dāng)然, 所以正確的選擇時(shí)間點(diǎn)是地應(yīng)力平衡計(jì)算中最最關(guān)鍵的.處理地應(yīng)力平衡的

6、三種方法:1.地表水平土層分層水平的情況下的地應(yīng)力平衡0 x- : C% z3 L0 R9 K3 |) _" d地表水平而且土層分層水平的情況下的地應(yīng)力平衡是最簡(jiǎn)單的情況, 事實(shí)上也是大多數(shù)計(jì)算用到的平衡方法, 這個(gè)時(shí)候可以用ABAQUS提供的*initial conditions, type=stress, geostatic方法來(lái)做. 具體怎么寫(xiě)這個(gè)語(yǔ)句, 請(qǐng)自行參考文檔.2. 其它的情況下的地應(yīng)力平衡  ' s4 x& % I8 A# x* c1 Y$ q地表不水平或土分層不規(guī)那么的情況下只能使用導(dǎo)入地應(yīng)力數(shù)據(jù)的方式進(jìn)行地應(yīng)力平衡, 其方法的

7、核心是給土體加重力, 不必理會(huì)其土體位移, 只是單純的提取每個(gè)單元的平均應(yīng)力, 再將這個(gè)應(yīng)力讀回ABAQUS作為初始應(yīng)力, 以到達(dá)平衡土體位移的效果. 具體操作參考這個(gè)貼 :/forum.simwe /viewthread.php?tid=762963&highlight=%B5%D8%D3%A6%C1%A6%C6%BD%BA%E2, G( i6 ( c: ! 5 p- q0 W3. 不預(yù)設(shè)地應(yīng)力所謂的不進(jìn)行地應(yīng)力平衡是指計(jì)算前不設(shè)定土體位移歸零, 而是記錄下土體在重力作用下每個(gè)結(jié)點(diǎn)的位移值, 在進(jìn)行后繼計(jì)算后, 將位移結(jié)果減去重力作用下每個(gè)結(jié)點(diǎn)的位移值就能得到當(dāng)前計(jì)算步的凈位移, 這

8、種方法只適用于小變形分析, 不適用于幾何非線(xiàn)性的情況, 好處是收斂性好, 精度高, 因?yàn)橛涗浀氖墙Y(jié)點(diǎn)位移, 缺點(diǎn)是麻煩, 因?yàn)橐崆坝涗? 還要處理數(shù)據(jù).( b8 m% l, w+ l7   g" P4 ' s   1法是對(duì)單元集進(jìn)行賦值, 與結(jié)點(diǎn)集無(wú)關(guān);常見(jiàn)的問(wèn)題( m  f/ L' 8 l4 d& _  m1. 有結(jié)構(gòu)物的情況如何平衡?有結(jié)構(gòu)物的情況下就不是正確的平衡點(diǎn), 如果一定要這樣平衡, 比方原來(lái)有一個(gè)隧道, 現(xiàn)在在上面建房子, 要估計(jì)這種情況下的追加位移是多少

9、, 這種情況下可以用上面的方法2和3來(lái)完成.3 H% g" d  f+ |( M5 N2. 多層土的情況下的平衡?7 V" Z1 f5 Q: & / v9 C多層土的情況可以根據(jù)不同的情況選擇方法1或2進(jìn)行.0 C9 n, 0 X) s: Z: s3. 有地下水的情況下的平衡?' D9 N* S' Z7 y+ ( ?有水的情況下最主要是要搞清楚土的三相, 按照不同的孔隙率計(jì)算土的干密度和濕密度, 這種平衡目前最好是只用方法1來(lái)完成.2 t- R# B# + w7 . O; K: D0 q7 b8 f4. 邊坡加重力無(wú)法平衡?, r

10、. E, j0 n: F* L4   |; 7 邊坡加重力無(wú)法平衡表示邊坡本身是不穩(wěn)的, 請(qǐng)檢查模型是否合理, 塑性參數(shù)是否合理./ :/bbs.simwe /thread-753043-1-12.html  aba_aba 大哥解答：你的inp文件中有part和instance, 而geo.dat中的節(jié)點(diǎn)編號(hào)前面沒(méi)有instance名稱(chēng)。有兩種方法：1。建立沒(méi)有part的inp文件，geo.dat中的單元編號(hào)前面沒(méi)有instance名稱(chēng)。見(jiàn) :/ simwe /forum/viewthread.php?tid=7089422。 建立有part

11、的inp文件，geo.dat中的單元編號(hào)前面加上instance名稱(chēng)。見(jiàn) :/ simwe /forum/thread-724821-1-1.html" . R! y% M4 g: O! z! y4 k8 F8 G8 Q! h+ i- |7 按照第2種方法，在你的geo.dat文件中單元編號(hào)前面加上 Part-1-1. ， 例如第一行變?yōu)? S* v) N  S# d; _" f2 & P( U, e7 j5 k6 QPart-1-1.1,-1.50E+04,-3.50E+04,-1.50E+04,-1.25E-12,-1.92E-13,-4.1

12、7E-12: :/forum.simwe /archiver/tid-708942.html    :先建立好你的模型，按照你的課題把邊界條件、接觸什么的都建好，然后使用命令mdb.models'你的模型名字'.setValues(noPartsInputFile=ON)，建立沒(méi)有part的inp文件abaqus不認(rèn)帶part的istress數(shù)據(jù)zhuyi:要寫(xiě)你的model的名字在CAE窗口頂部，Model右面，默認(rèn)名稱(chēng)為Model-1, Model-2等等，不是CAE或JOB的名稱(chēng)。 你的model的名字前后都要加上單引號(hào)地應(yīng)力平衡方

13、法熊志勇 陳功奇第一局部 地應(yīng)力平衡方法簡(jiǎn)介地應(yīng)力平衡有三種方法：1*initial conditions,type=stress,input=FileName.csv(或inp)該方法中的文件FILENAME.INP獲取方法為:首先將邊界條件施加到模型上進(jìn)行正演計(jì)算,然后一般是將計(jì)算得到的每個(gè)單元的應(yīng)力外插到形心點(diǎn)處并導(dǎo)出6個(gè)應(yīng)力分量(也可以導(dǎo)出積分點(diǎn)處的應(yīng)力分量,視要求平衡的精確程度而定)。其所采用的幾何模型可以考慮地表起伏不平的情況以及巖土材料極其不均勻的情況,適用范圍廣。但由于外插的應(yīng)力有一定誤差,因此采用彈塑性本構(gòu)模型時(shí),可能會(huì)導(dǎo)致某些點(diǎn)的高斯點(diǎn)應(yīng)力位于屈服面以外,當(dāng)大面積的高斯點(diǎn)上

14、的應(yīng)力超出屈服面之后,應(yīng)力轉(zhuǎn)移要通過(guò)大量的迭代才能完成,而且有可能出現(xiàn)解不收斂的情況。在僅考慮自重情況下只能考慮受泊松比的影響帶來(lái)的側(cè)壓力系數(shù)效應(yīng),因此平衡后的效果不一定很理想,但無(wú)疑其適用性很強(qiáng)。2*initial conditions,type=stress,geostatic該方法需給出不同材料區(qū)域的最高點(diǎn)和最低點(diǎn)的自重應(yīng)力及其相應(yīng)坐標(biāo)。所采用的幾何模型一般較規(guī)那么,外表大致水平,地應(yīng)力平衡的好壞一般只受巖體密度的影響,無(wú)論采用彈性或彈塑性本構(gòu)模型都能很好的到達(dá)平衡,可以不必局限于僅受泊松比的影響,能夠通過(guò)考慮水平兩個(gè)方向的側(cè)壓力系數(shù)值來(lái)施加初始應(yīng)力場(chǎng)。計(jì)算速度快,收斂性好。缺點(diǎn)就是不能

15、夠很好平衡具有起伏外表的幾何模型,需知道平整后模型的上覆巖體自重。3*initial conditions,type=stress,geostatic,user該方法采用用戶(hù)子程序SIGINI來(lái)定義初始應(yīng)力場(chǎng),可以定義其為應(yīng)力分量為坐標(biāo)、單元號(hào)、積分點(diǎn)號(hào)等變量的函數(shù),要到達(dá)精確平衡需具體邊界條件,在實(shí)際中應(yīng)用較少。第二局部 地應(yīng)力平衡方法實(shí)例詳解地應(yīng)力平衡是巖土工程數(shù)值模擬分析的重要的內(nèi)容，為了讓師弟師妹們快點(diǎn)上手，我利用第一種方法做一個(gè)較簡(jiǎn)單的模型，希望對(duì)大家有用。一、 模型描述：二、地應(yīng)力平衡過(guò)程1. 啟動(dòng)ABAQUS,單擊Create Model Datebase2. 創(chuàng)立部件Part在

16、Part模塊，單擊創(chuàng)立部件按鈕，彈出如右圖的對(duì)話(huà)框，按圖輸入部件名：Part-soil；采用二維模型選擇2D Planar；Type選擇可變型Deformable；根本特征選擇殼體Shell；Approximate size輸入70，這個(gè)數(shù)值的大小，應(yīng)根據(jù)模型的最大尺寸來(lái)確定：稍大于最大尺寸的2倍。比方本模型最大尺寸是30那么我輸入了70，但也不是絕對(duì)，你當(dāng)然也可以輸入65，或75等等。最后單擊Continue,繼續(xù)下一步。按照模型尺寸如圖，建立模型部件，雙擊鼠標(biāo)中鍵，完成部件的建立。3. 建立材料屬性P roperty在Module中切換到P roperty模塊，單擊，輸入材料名稱(chēng)name:

17、 Material-soil,單擊Density,在彈出對(duì)話(huà)框中輸入：密度2080；然后單擊MechaniacalElasticityElastic,在彈出的對(duì)話(huà)框輸入圖，單擊OK完成材料的定義。單擊按鈕，輸入名字：S ection-soil,選擇Soild , Homogeneous,單擊Continue，OK,完成截面的創(chuàng)立。單擊按鈕，選取部件單擊或框選，選擇后成粉紅色，表示選中，單擊Done或單擊鼠標(biāo)中鍵來(lái)確定。在彈出的對(duì)話(huà)框中選中Section-Soil，單擊OK.4. 裝配部件Assembly在Module選擇Assembly模塊,單擊，彈出對(duì)話(huà)框，采用默認(rèn)值，單擊OK.注意：1.本

18、模型只有一個(gè)部件所以自動(dòng)選中，如有多個(gè)可按Shift鍵全選中； 2.Instance Type 本例都無(wú)所謂，只有一個(gè)部件，但是如有多個(gè)部件，我比擬喜歡選擇Independent,因?yàn)樵谶@種情況下，所有的部件會(huì)出現(xiàn)在一個(gè)窗口，劃分網(wǎng)格更方便，并且直觀的看到不同部件連接處網(wǎng)格劃分的是否協(xié)調(diào)。5. 創(chuàng)立分析步Step 在Module選擇Step模塊,單擊，彈出下面的對(duì)話(huà)框，創(chuàng)立分析步Step-1,選擇Geostatic,單擊Continue，在分析步編輯框Basic中選擇On(如果會(huì)發(fā)生大變形的情況下要選擇On)，其他默認(rèn)，單擊OK6. 施加荷載和定義邊界條件Load在Module選擇Load模塊

19、，單擊定義邊界條件，需要對(duì)模型的左、右、下底面定義邊界條件，分別命名BC-1、BC-2 、BC-3，選擇初始步Initial,Displacement/Rotation, 單擊Continue，選擇左、右、下底面邊界，單擊Done，分別選擇U1、 U1、 U2U1是水平方向，U2是豎直方向，單擊OK完成邊界條件的定義。單擊定義重力荷載，Name:Load-grv, 選擇Step-1，Mechanical，Gravity，單擊Continue。單擊Edit Region，選擇整個(gè)模型，在Component 2中輸入-9.8重力加速度，單擊OK7. 劃分網(wǎng)格Mesh劃分網(wǎng)格是一門(mén)藝術(shù)，本例只是簡(jiǎn)單

20、的劃分在Module選擇Mesh模塊，單擊設(shè)定網(wǎng)格的種子，將全局種子大小設(shè)為1，其余默認(rèn)，單擊OK.單擊，采用默認(rèn)，單擊OK。單擊，F(xiàn)amily選擇Plane Strain，其余采用默認(rèn)，單擊OK單擊，單擊Yes,完成網(wǎng)格的劃分，如以下圖：8. 在命令行中輸入mdb.models'Model'.setValues(noPartsInputFile=ON)，按回車(chē)鍵9. 在job模塊中創(chuàng)立名為Job-NoInitialCondition的分析步，提交分析。應(yīng)力圖水平位移圖10. 將分析得到的應(yīng)力場(chǎng)保存為一個(gè)文本文件。具體做法:翻開(kāi)分析得到的ODB文件，選擇菜單ReportFiel

21、d Output,在以下圖所示的對(duì)話(huà)框中，選中積分點(diǎn)上的各個(gè)應(yīng)力分量對(duì)于二維問(wèn)題，應(yīng)力分量S11、S22、S33和S12；對(duì)于三維問(wèn)題，還應(yīng)選中S13和S23。圖表 1輸入常變量S11、S22、S33和 S12單擊此對(duì)話(huà)框中的Setup標(biāo)簽頁(yè)，在Name文本框中輸入要保存的文件名b.inp，取消對(duì)Append to file 的選擇即創(chuàng)立一個(gè)新文件，在Write后面只選中Field Output如以下圖所示。圖表 2設(shè)置輸出場(chǎng)變量注意，此處輸出的是當(dāng)前增量步結(jié)束時(shí)的應(yīng)力結(jié)果，因此上述對(duì)話(huà)框頂部的Step必須是Geostatic分析步，F(xiàn)rame必須是1。如果Frame是0，會(huì)看到輸出的應(yīng)力都是0。11. 按照ABAQUS所要求的初始應(yīng)力場(chǎng)文件格式，修改上述文件b.inp中的內(nèi)容。具體方法為：用Excel翻開(kāi)上述文件b.inp，在文本文件導(dǎo)入向?qū)У牟襟E1中選擇“分隔符號(hào)，在步驟2中選擇“Tab鍵和“空格鍵，這樣b.inp中的各列數(shù)據(jù)就成為Excel表格中的各個(gè)列。 刪除表格中開(kāi)始幾行的模型信息，再刪除積分點(diǎn)編號(hào)所在的第2列數(shù)據(jù)都為數(shù)字1，只保存單元編號(hào)和各個(gè)應(yīng)力分量列，并將各個(gè)應(yīng)力分量的科學(xué)計(jì)數(shù)法格式改為顯示小數(shù)點(diǎn)后5位數(shù)字。修改前和