★地應力平衡圖解

*initialconditions,type=stress,input=0.csvmdb.models['Model-1'].setValues(noPartsInputFile=ON)（將網(wǎng)格劃分精密一些，地應力平衡精度更高一些，2010-7-12注）注意：導入的inp、csv或其他文件里不得有空行，否則會出現(xiàn)element0。允許有空格。--------------------------------------------------------------------方法一：樁土地應力平衡不容易通過，可先直接指定應力：*initialconditions,type=stress,geostaticSet-pile-soil,0,0,-510000,30,0.6,0.6將樁土摩擦系數(shù)設為0，容易平衡，輸出含有應力s11、s22等的應力文件；導入應力文件，將摩擦系數(shù)提高為一個較小的數(shù)值（低于正常的摩擦系數(shù)），再計算，輸出應力文件；再提高摩擦系數(shù)，導入應力文件，計算，重復上述操作，直到摩擦系數(shù)達到正常值。上述如果不行，可以干脆將樁與的接觸改為tie，2010-8-18的模型按照上面的步驟操作不行，后來tie就可以了，且收斂的速度很快。樁土之間用摩擦容易出現(xiàn)“Theratioofdeformationspeedtowavespeedexceeds1.0000inatleastoneelement.”的錯誤。而用tie，對于動力學更容易收斂，一般不會出現(xiàn)上述錯誤提示（2010-8-19的錘擊樁模型將土分為內(nèi)外兩層，內(nèi)層網(wǎng)格密度與樁劃分的密）。方法二：樁土地應力平衡方法：（建議此法）1、按照實際工程要求建立樁土模型model-1；2、模型樹上的model-1上面copy，復制成另外一個模型model-2；3、把mode-2里的所有與土無關的東西全部刪除，包括part、assembly、set、interaction、load等；4、在load模塊的BC里約束原來土與樁接觸的面，約束住土不要進入樁孔內(nèi)（豎向不要約束）；5、注意：土不能受到重力方向的約束，包括土四周的邊界，讓土自由下沉。6、不用導入初始應力，提交model-2的job，產(chǎn)生應力文件-2，編輯成ABAQUS要求的格式；6、將應力文件-2導入model-1，提交計算model-1的job，產(chǎn)生應力文件-1，編輯成ABAQUS要求的格式，重新導入model-1，提交，平衡完畢；7、刪除輔助用的model-2及其產(chǎn)生的文件。該法通過率高，精度高（2010-8-19用錘擊樁模型試算，用此法平衡達到10e-4，而用其它方法平衡的精度達10e-3）。問題：2010-8-19的模型（錘擊樁內(nèi)外土tie的模型）用該法沒有平衡成功，還是用方法一平衡成功。地應力平衡方法：

第一步：建立模型，材料，分析步（GEOSTATIC）

第二步：施加荷載，LOAD，選擇施加重力GRAVITY，在你想施加重力的方向輸入數(shù)值9.8Y

第三步：在命令行中輸入mdb.models['模型名字'].setValues(noPartsInputFile=ON)

(請嚴格按照這個格式，注意大小寫的字母)，例如：mdb.models['Model-1'].setValues(noPartsInputFile=ON)第四步：提交J0B，完成后

第五步：按以下步驟,Roport---ReportFieldOutput---選中S11,S22,S33,S12,S13,S23---Name:XX.INP

---Write中選擇FieldOutput-------------ok!!!第六步：用軟件(推薦使用UltraEdit很好編輯的)打開XX.INP[注：先用UltraEdit打開，刪除中間的英文，再用excel打開]，保存格式內(nèi)容

單元號

S11

S22

S33

S12

S13

S23（請注意，在保存內(nèi)容中沒有這一行的）

1

,

.

,

.

,

.

,

.

,

.

,

.

2

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

這個結果文件是最重要的，在所保存的文件中只有數(shù)值部分，沒有英文字母，沒有

上面那個“單元號”這一行，而且單元號前面也沒有什么PART名字什么的，就是1，2.

.......這些數(shù)字。例如：*elprints,*initialconditions,type=stress,geostaticSet-all,0,0,-510000,30,0.6,0.6*initialconditions,type=stress,geostaticSet-pile,0,0,-250000,10,0.2,0.2Set-soil,0,0,-255000,15,0.6,0.6一個part分層土地應力平衡SHAPEPart-1是個部件，要在assembly組裝成實體instance，instance的名字一般是在部件的名字后面加“-1”，例如部件part-1形成實體part-1-1。兩個part分層土地應力平衡SHAPE另外一種方法1：在命令行中輸入mdb.models['模型名字'].setValues(noPartsInputFile=ON)

(請嚴格按照這個格式，注意大小寫的字母)2：提交J0B，完成后3：按以下步驟,Roport---ReportFieldOutput---選中S11,S22,S33,S12,S13,S23---Name:XX.INP---Write中選擇FieldOutput-------------ok!!!以上3個步驟可以代替以下命令：*elprints,手工輸入地應力平衡SHAPE*initialconditions,type=stress,geostaticSet-pile,0,30,-500000,10,0.2,0.2set-soil,0,30,-540000,0,0.6,0.6對于很厚的土層必須先進行地應力平衡，否則有些單元會出現(xiàn)過大的扭曲變形，最后運算不了，出現(xiàn)ERROR:FIXEDTIMEINCREMENTISTOOLARGE（即使在step中打開非線性開關也不行），必須先手工輸入地應力平衡初始條件，見上面。在inp中加入*elprints,將各單元的應力導出（*.csv），再*initialconditions,type=stress,input=*.csv導入初試應力（程序自動將取排在后面的初始應力為準，自動把第一個初始應力除去如下圖，*initialconditions,type=stress,geostatic和下面一行自動被除去），這樣迭代平衡效果非常好。對于采用c3d8等非縮減積分單元，會出現(xiàn)如下圖的應力文件：用下面的方法就OK：（不要選擇積分點，而是質(zhì)心）由于老板項目多，自己論文忙，好長一段時間都沒來SIMWE了，一上來就收到不少朋友的消息，問我的地應力平衡方法，現(xiàn)在我把整個過程詳細說明一下，希望對大家有所幫助。

地應力平衡方法：

第一步：建立模型，材料，分析步（GEOSTATIC）

第二步：施加荷載，LOAD，選擇施加重力GRAVITY，在你想施加重力的方向輸入數(shù)值9.8

第三步：在命令行中輸入mdb.models['模型名字'].setValues(noPartsInputFile=ON)

(請嚴格按照

這個格式，注意大小寫的字母)

第四步：提交J0B，完成后

第五步：按以下步驟,Roport---ReportFieldOutput---選中S11,S22,S33,S12,S13,S23---Name:XX.INP

---Write中選擇FieldOutput-------------ok!!!

第六步：用軟件(推薦使用UltraEdit很好編輯的)打開XX.INP，保存格式內(nèi)容

單元號

S11

S22

S33

S12

S13

S23（請注意，在保存內(nèi)容中沒有這一行的）

1

,

.

,

.

,

.

,

.

,

.

,

.

2

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

這個結果文件是最重要的，在所保存的文件中只有數(shù)值部分，沒有英文字母，沒有

上面那個“單元號”這一行，而且單元號前面也沒有什么PART名字什么的，就是1，2.

.......這些數(shù)字。

第七步：在ABAQUS----Model---Editkeywords---Model-1(這就是你的Model名字)---在材料屬性后面

加上：*initialconditions,type=stress,input=xx.inp

完成

第八步：重新提交JOB，OK

第九步：如果你還沒有成功的話，那我只能說----------------我無語了。：）

ABAQUS的這項功能確實很不錯。：）這個功能讓基坑開挖、隧道開挖等的初始應力，開挖后的殘余應力很好的顯示；也可以很好的模擬鐵路設計中的工后沉降的概念，在地應力平衡后，加上荷載所得沉降即為工后沉降；也很好的模擬了樁土復合地基的問題，如果沒有初始應力的模擬，使土對樁產(chǎn)生了擠壓應力，從而通過設定摩擦系數(shù)就可以模擬了樁與土之間的摩擦力；除此之外，在進行擋土墻計算時也需要ABAQUS的這項功能，反正很多都用得著。

希望我的建議能幫助初學者（其實我也是初學者，呵呵）很快掌握這個技巧。原文地址：\o"有關ABAQUS數(shù)據(jù)的后處理"有關ABAQUS數(shù)據(jù)的后處理作者：\o"舉舉"舉舉三種方法：1．

XYDATA

通過Tools-XYdata-edit,thencopytoexcelormatlab,etc.2．

通過Plugins-Tools-excelutilities直接在Excel中輸出數(shù)據(jù)和圖形3．

通過Report-XY，將結果輸出至abaqus.rpt，此文件后綴名可改為.txt，通過MATLAB讀取。樁土共同作用對于新手來說一直是個難題，本人結合3年做樁土分析的經(jīng)驗，將樁土分析中遇到的問題做個總結：

主要思想是在計算土體自重時將樁殺死，然后約束土體中樁孔的橫向自由度，建立step1，通過靜力分析步計算土體的自重應力，然后將土體的自重應力編輯后，通過initialconditions，type=stress，input=自重應力文件（可以是txt或是inp文件）導入到原模型中，然后添加step2分析步，在step2分析步中，激活樁，同時建立接觸，施加樁的自重及其他荷載，即可完成樁土的共同作用分析。具體操作過程如下：

1.首先建立模型1，如model1，在該模型里建立樁和土的幾何模型，并賦予材料屬性，將所有的樁設置為一個集合，并命名（如z1），建立step1，使用靜力分析步，注意，在editkeyword中在step1定義的后面添加兩行語句，第一行為*modelchange，remove，第二行為z1（即前面定義的樁的集合），然后在荷載里定義土體的重力加速度，約束土體的位移，同時約束土體中樁孔的橫向位移（如，z為重力加速度方向，則約束x和y的自由度），計算分析，得到土體的自重應力，編輯自重應力文件，命名為XXX.txt。

2.進入到模型1中，將模型1另存為模型2，此時，模型2和模型1是相同的，在模型2中，添加step2，編輯editkeyword，在邊界條件定義前，添加*initialconditions，type=stress，input=XXX.txt，然后在step2定義后面添加兩行語句，第一行為*modelchange，add第二行為z1（樁的集合名），然后在step2中建立樁和土的接觸，施加樁的自重，然后在step中殺死土體中樁孔的橫向自由度（對應step1中樁孔的橫向自由度約束），后面的荷載可以在step2中施加，也可以建立新的step施加。

3.補充：建立接觸的分析步（本例中為step2，初始步要設置的小的，本人一般設置為0.01，這樣有利于接觸的順利建立）；當模型較大時，地應力文件的編輯工作量很大，本人找朋友編了個小程序，可以將地應力文件自動生成abaqus要求的格式，該程序?qū)⒃诟郊薪o出。附件太大，無法上傳，悲劇[local]1[/local]\o"匯總"匯總,

\o"經(jīng)驗"經(jīng)驗查看3D實體截面的內(nèi)力，如截面合力、截面總彎矩，或者截面各方向的分力、分彎矩，按照aba-aba大俠編寫的《ABAQUS常見問題匯總》，可以設置cutsurface來實現(xiàn)。但這種方法較難，而且對每一個需要查看內(nèi)力的截面都要單獨設置，但要查看內(nèi)力的截面較多時，就更麻煩。通過odb文件可以較為簡單的獲取截面內(nèi)力。具體操作步驟如下：

a)

打開odb文件，在窗口中只保留要獲取內(nèi)力的部件，如某根柱、梁等；若保留整個模型，則獲取的截面內(nèi)力是整個模型在該截面的內(nèi)力；