『壹』 abaqus怎麼計算土體的地應力
我也是新手這些資料你看看吧！都是在網上找的，看完應該能明白個大概！
初始地應力平衡應該滿足：
（1） 豎向位移為零或者很小，一般至少10的負6次方，甚至更小
（2） 豎向應力s22要等於土體深度*密度*g
關於地應力的平衡方法,綜合了版上的一些意見,結合了自己的想法,對於初始地應力的施加,得到了e-6的效果,方法比較簡單,與大家分享!
1.先施加重力荷載的作用,可以在cae中實現;
2.在inp文件中的output request中寫上
*el print
s,
這樣就會將施加重力荷載後的應力輸出到*.dat文件中了;
3.在*.dat文件中,將單元應力的序號及單元的應力拷出,例如
ELEMENT T FOOT- S11 S22 S33 S12
NOTE
1 1 -1.2598E+05 -1434. -3.1852E+04 892.7
2 1 -1.2249E+05 -6287. -3.2194E+04 1223.
3 1 -1.1795E+05 -497.7 -2.9611E+04 1664.
4 1 -1.1210E+05 -7240. -2.9834E+04 1992.
5 1 -1.0485E+05 579.0 -2.6068E+04 2600.
6 1 -9.5803E+04 -8272. -2.6019E+04 3031.
7 1 -8.4709E+04 1915. -2.0698E+04 4083.
8 1 -7.0634E+04 -9746. -2.0095E+04 4339.
9 1 -5.1088E+04 5401. -1.1422E+04 8519.
10 1 -2.4353E+04 -1.1150E+04 -8876. 1.2126E+04
11 1 -1.2847E+05 268.1 -3.2050E+04 738.1
12 1 -1.2786E+05 -9868. -3.4433E+04 629.1
13 1 -1.2938E+05 -4224. -3.3402E+04 502.5
14 1 -1.3039E+05 -3458. -3.3461E+04 165.9
單獨存為一個*.dat文件,
4.用excel打開該文件,將其中的1所在的列去掉,在每個單元號前面加上其instance. ,即單元編號變為: instance名稱.序號 ;注意不同的instance和part要都按照其所在的單元從小到大編號,而不是按照他們在整體單元編號來編號!
5.接下來就在excel把該文件另存為*.csv格式的文件(即帶有逗號分隔符的格式),
6.最後在inp文件的step之前寫上
*initial conditions,type=stress,input=文件名.csv即可!
這種方法不需要用python,比較簡單,希望能對大家有用!
先說為什麼要施加地應力：
1、我們所建立的幾何模型一般和工程實際情況或尺寸相對應、相一致，比如邊坡幾何模型和實際邊坡尺寸一致，但我們可以誇張一點想像，實際邊坡應是由一個更大一點或更高一點的不受重力的初始邊坡在n年前突然受重力和類似目前的邊界條件作用下逐漸形成了今天的尺寸大小，n年前受重力和類似目前的邊界條件作用之前邊坡的尺寸大小，我們不得而知，如果能准確知曉，我們就可以建立一個那時的幾何模型，再施加重力和邊界條件進行計算，變形後形狀和現狀邊坡形狀一致，其內力也就是初始應力場或地應力，就不用專門去施加地應力了，但問題是我們不能知曉邊坡受力前的形狀尺寸，我們現在的幾何模型就是邊坡現在的實際尺寸，受力後將會變成一個更小的或與現狀不一致的邊坡，這不符合我們模擬現狀邊坡的目的。如果我們知道現狀邊坡的內力，將其提取出來作為幾何模型的內力，再和外力（重力）平衡，則我們建立的模型才能算和實際模型一致。
真實地知道現狀邊坡的內力是很難的，我們採取的辦法是，用我們所建立的幾何模型施加和實際模型一致的重力和邊界條件進行計算，得到變形後或變得更小或與現狀邊坡不完全一致的邊坡內力近似的作為現狀邊坡的內力，並重新將其施加於與現狀邊坡一致的幾何模型，再施加重力（當然邊界條件也應基本一致）以平衡，這樣才算建立了與現狀模型基本一致的模型，其下的計算才成為可能。這就是所謂「地應力平衡」的含義、目的、作用。
2. 地應力平衡中的外力和內力的問題，地應力平衡中，顯然，重力是外力，應力場是內力，僅有外力重力，沒有內力是不可能的，同樣，僅有內力（專指初始應力場）而不受重力也是不可能的，否則，整個體系的力不會平衡。這就是為什麼我們將提取出的內力施加於幾何模型後必須再施加重力的原因。為的是內力和外力平衡。) q0 F3 q6 H1 O#
3. 地應力場的方向問題，有網友在論壇里問，既然重力是向下，為與重力平衡，那應力場的方向是不是向上呢，這同樣是我開始接觸abaqus的疑問，相信很初學者也有這樣的疑問，我的理解是內力是沒有向上、向下或者向其它方向的概念的，內力只有拉力或壓力或剪力之分，其方向也按是拉是壓是順時針或逆時針而分，內力往往都是成對出現，如地應力場中的應力以壓應力為主，取一個微元，則壓應力同時出現在向下和向上，你能說地應力就是向上，與重力反向嗎？
不怕各位笑話，以上幾點在高手看來是很簡單的問題，卻是我經歷了漫長而艱辛的摸索才得到的，今天也寫給初學者，不要再走我的老彎路了。
aba中初始地應力場平衡一般在表面水平的情況下僅僅和密度相關，密度一樣的話平衡的結果很好，別的參數改變之後經過計算，差別很小。表面不水平的情況則最好通過文件導入初應力的情況進行平衡。一般來講，表面不平的時候有很多因素造成誤差很大。
地應力平衡方法：
第一步：建立模型，材料，分析步（GEOSTATIC）
第二步：施加荷載，LOAD，選擇施加重力GRAVITY，在你想施加重力的方向輸入數值9.8
第三步：在命令行中輸入mdb.models['模型名字'].setValues(noPartsInputFile=ON) (請嚴格按照
這個格式，注意大小寫的字母
第四步：提交J0B，完成後
第五步：按以下步驟,Roport---Report Field Output---選中S11,S22,S33,S12,S13,S23---Name:XX.INP
---Write中選擇Field Output-------------ok!!!
[第六步：用軟體(推薦使用UltraEdit很好編輯的)打開XX.INP，保存格式內容
單元號 S11 S22 S33 S12 S13 S23 （請注意，在保存內容中沒有這一行的）
, . , . , . , . , . ,
, . , . , . , . , . , .
, . , . , . , . , . , .
, . , . , . , . , . , .
這個結果文件是最重要的，在所保存的文件中只有數值部分，沒有英文字母，沒有!
上面那個「單元號」這一行，而且單元號前面也沒有什麼PART名字什麼的，就是1，2.
.......這些數字。
第七步：在ABAQUS----Model---Edit keywords---Model-1(這就是你的Model名字)---在材料屬性後面
加上：*initial conditions,type=stress,input=xx.inp 完成,
第八步：重新提交JOB，OK
第九步：如果你還沒有成功的話，那我只能說----------------我無語了。
ABAQUS的這項功能確實很不錯。：）這個功能讓基坑開挖、隧道開挖等的初始應力，開挖後的殘余應力很好的顯示；也可以很好的模擬鐵路設計中的工後沉降的概念，在地應力平衡後，加上荷載所得沉降即為工後沉降；也很好的模擬了樁土復合地基的問題，如果沒有初始應力的模擬，使土對樁產生了擠壓應力，從而通過設定摩擦系數就可以模擬了樁與土之間的摩擦力；除此之外，在進行擋土牆計算時也需要ABAQUS的這項功能，反正很多都用得著。
希望我的建議能幫助初學者（其實我也是初學者，呵呵）很快掌握這個技巧。
地應力平衡的正確時間點
所謂地應力平衡是指, 當我們建任何東西或挖任何東西之前, 地表的位移都是零, 但是土體的應力卻存在, 這種無位移但有應力的時間點叫地應力平衡. 那麼正確的時間點應該是什麼呢, 就是不管土原來是什麼樣的, 高山, 河流, 丘陵, 平原都可以, 關鍵是在我們對它做任何擾動之前. 有些人問樁土作用如何平衡? 這個問題首先就是錯的, 因為地應力平衡的時候是不能有樁的. 一旦有樁就跨過了地應力平衡這個時間點, 數值模擬要忠實於實際情況, 不能想當然, 所以正確的選擇時間點是地應力平衡計算中最最關鍵的.
處理地應力平衡的三種方法
1.地表水平土層分層水平的情況下的地應力平衡
地表水平而且土層分層水平的情況下的地應力平衡是最簡單的情況, 事實上也是大多數計算用到的平衡方法, 這個時候可以用ABAQUS提供的*initial conditions, type=stress, geostatic方法來做. 具體怎麼寫這個語句, 請自行參考文檔.
2. 其它的情況下的地應力平衡 [' s4 x& ]% I8 A# x* c1 Y$ q
地表不水平或土分層不規則的情況下只能使用導入地應力數據的方式進行地應力平衡, 其方法的核心是給土體加重力, 不必理會其土體位移, 只是單純的提取每個單元的平均應力, 再將這個應力讀回ABAQUS作為初始應力, 以達到平衡土體位移的效果. 具體操作參考這個貼
http://forum.simwe.com/viewthread.php?tid=762963&highlight=%B5%D8%D3%A6%C1%A6%C6%BD%BA%E2,
3. 不預設地應力
所謂的不進行地應力平衡是指計算前不設定土體位移歸零, 而是記錄下土體在重力作用下每個結點的位移值, 在進行後繼計算後, 將位移結果減去重力作用下每個結點的位移值就能得到當前計算步的凈位移, 這種方法只適用於小變形分析, 不適用於幾何非線性的情況, 好處是收斂性好, 精度高, 因為記錄的是結點位移, 缺點是麻煩, 因為要提前記錄, 還要處理數據.
1法是對單元集進行賦值, 與結點集無關;
常見的問題
1. 有結構物的情況如何平衡?
有結構物的情況下就不是正確的平衡點, 如果一定要這樣平衡, 比如原來有一個隧道, 現在在上面建房子, 要估計這種情況下的追加位移是多少, 這種情況下可以用上面的方法2和3來完成.
2. 多層土的情況下的平衡?
多層土的情況可以根據不同的情況選擇方法1或2進行.
3. 有地下水的情況下的平衡?
有水的情況下最主要是要搞清楚土的三相, 按照不同的孔隙率計算土的干密度和濕密度, 這種平衡目前最好是只用方法1來完成.
4. 邊坡加重力無法平衡?
邊坡加重力無法平衡表示邊坡本身是不穩的, 請檢查模型是否合理, 塑性參數是否合理./
我是新手、自學者，總結了一點經驗，主要是給和我一樣的新手
*initial condition,type=stress,geostatic 就不說了，啃啃手冊就明白了，或者在論壇搜索下 geostatic
我主要說一下*initial conditions,type=stress,input=XXX.XXX的方法：
1、先建立好你的模型，按照你的課題把邊界條件、接觸什麼的都建好，然後使用命令mdb.models['你的模型名字']setValues(noPartsInputFile=ON)，建立沒有part的inp文件（abaqus不認帶part的istress數據）
2、建立*geostatic的step，只施加重力荷載，然後求解
3、完成之後用viewer打開odb，選report->field report,下拉菜單裡面選擇 centriod，然後依次把s11、s22、s33、s12、s13、s23點選上，去setup頁面選擇報告文件的名字（例如istres.rpt)把下面的和值之類的額外輸出點掉，要不還要手動刪除
4、編輯你的報告文件，把數據行之外的東西都刪除。（正確的格式：單元號，s11，s22，s33，s12，s13，s23），使用excel或者ultraedit在每個數據間加上一個逗號（ , ），數據之間的空格不要緊，但是要注意不要有空行或者逗號之間是空白，否則會出現element 0錯誤
5、重新建立你的step，寫入新的inp文件之後，編輯它，在材料屬性之後、step之前加入*initial conditions,type=stress,input=XXX.XXX
搞定，位移一般都被控制在10-6數量級，滿足要求。
在此感謝論壇的各位前輩，也建議新人多搜索一下論壇的經典舊帖
---------------------------------------
補充一下：
使用mdb.models['你的model的名字'].setValues(noPartsInputFile=ON) 時注意:
- 要寫你的model的名字（在CAE窗口頂部，Model右面，默認名稱為Model-1, Model-2等等），不是CAE或JOB的名稱。
- 你的model的名字前後都要加上單引號
- 這個命令要在CAE窗口底部輸入（ 先點左下角的 >>> ）
- 輸入此命令後，要在Job模塊中用write job來生成INP文件（loveestboy1981 是不是沒作這一步？）
- 生成的INP文件名和Job名一樣，而不是CAE文件或Model的名稱。
- 生成的INP文件在當前的工作目錄下（不一定是TEMP文件夾）
『貳』 abaqus 混凝土梁 施加的荷載小呢，就可以complete，荷載稍大一點點，就不可以了，可能是什麼原因
首先，你把材料定義為線性部分，然後計算，看是不是可以算到你期望的30Kn的力，如果可以，就說明模型沒有什麼大的問題
其次，如果線性部分可以算過去，那就要檢查你的非線性材料屬性是不是合適的，，你可以使用你初步定義的非線性部分算下，確定出來是不是已經達到屈服的應力，看結果，與材料定義的屬性比較
再次，如果沒有達到屈服階段，那就是計算的載荷給的太劇烈了，你可以把載荷步的個數調整多一點，然後一定要保證大變形是打開的，在載入力的時候，用一個賦值曲線來光滑你的載入過程，讓力一點點的往上加
也就是說可能的原因，第一是材料屬性，第二是力載入的過程不夠緩慢
『叄』 新人求助！abaqus鋼筋混凝土箱梁模型，加重力算不下去
我是剛剛學abaqus，用梁單元模擬了鋼筋，嵌固約束在混凝土中，不加重力載荷加集中荷載都有結果，一加重力載荷就跑不下去，只有加很小的重力值才能算，和您情況一樣，請問您解決了嗎～求告知，萬分感謝！救救孩子吧～