★地應(yīng)力平衡圖解

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

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

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

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

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

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

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

單元號(hào)

S11

S22

S33

S12

S13

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

1

,

.

,

.

,

.

,

.

,

.

,

.

2

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

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

上面那個(gè)“單元號(hào)”這一行，而且單元號(hào)前面也沒(méi)有什么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一個(gè)part分層土地應(yīng)力平衡SHAPEPart-1是個(gè)部件，要在assembly組裝成實(shí)體instance，instance的名字一般是在部件的名字后面加“-1”，例如部件part-1形成實(shí)體part-1-1。兩個(gè)part分層土地應(yīng)力平衡SHAPE另外一種方法1：在命令行中輸入mdb.models['模型名字'].setValues(noPartsInputFile=ON)

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

地應(yīng)力平衡方法：

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

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

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

(請(qǐng)嚴(yán)格按照

這個(gè)格式，注意大小寫(xiě)的字母)

第四步：提交J0B，完成后

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

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

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

單元號(hào)

S11

S22

S33

S12

S13

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

1

,

.

,

.

,

.

,

.

,

.

,

.

2

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

.

,

.

,

.

,

.

,

.

,

.

,

.

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

上面那個(gè)“單元號(hào)”這一行，而且單元號(hào)前面也沒(méi)有什么PART名字什么的，就是1，2.

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

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

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

完成

第八步：重新提交JOB，OK

第九步：如果你還沒(méi)有成功的話，那我只能說(shuō)----------------我無(wú)語(yǔ)了。：）

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

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

XYDATA

通過(guò)Tools-XYdata-edit,thencopytoexcelormatlab,etc.2．

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

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

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

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

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

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

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

a)

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