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)在的幾何模型就是邊坡現(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)力平衡”的含義、目的、作用。地應(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)力和外力平衡。 地應(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)力只有拉力或壓力或剪力之分，其方向也按是拉是壓是順時(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) ，密度一樣的話平衡的結(jié)果很好， 別的參數(shù)改變之后經(jīng)過(guò)計(jì)算，差別很小。表面不水平的情況則最好通過(guò)文件導(dǎo)入初應(yīng)力的情況進(jìn)行平衡。一般來(lái)講，表面不平的時(shí)候有很多因素造成誤差很大。cdstudio版主的總結(jié)：/viewthread.php?tid=908486&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D470TOC\o"1-5"\h\z地應(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)力平衡的三種方法 :地表水平土層分層水平的情況下的地應(yīng)力平衡 地表水平而且土層分層水平的情況下的地應(yīng)力平衡是最簡(jiǎn)單的情況 ,事實(shí)上也是大多數(shù)計(jì)算用到的平衡方法 ,這個(gè)時(shí)候可以用 ABAQUS提供的*initialconditions,type=stress,geostatic方法來(lái)做 .具體怎么寫(xiě)這個(gè)語(yǔ)句 ,請(qǐng)自行參考文檔 .其它的情況下的地應(yīng)力平衡 地表不水平或土分層不規(guī)則的情況下只能使用導(dǎo)入地應(yīng)力數(shù)據(jù)的方式進(jìn)行地應(yīng)力平衡 ,其方法的核心是給土體加重力 ,不必理會(huì)其土體位移,只是單純的提取每個(gè)單元的平均應(yīng)力 ,再將這個(gè)應(yīng)力讀回 ABAQUS作為初始應(yīng)力 ,以達(dá)到平衡土體位移的效果 .具體操作參考這個(gè)貼/viewthread.php?tid=762963&highlight=%B5%D8%D3%A6%C1%A6%C6%BD%BA%E2 ,G（i6]（c:}!@5p-q0W不預(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ì)算步的凈位移 ,這種方法只適用于小變形分析 ,不適用于幾何非線性的情況 ,好處是收斂性好 ,精度高 ,因?yàn)橛涗浀氖墙Y(jié)點(diǎn)位移 ,缺點(diǎn)是麻煩 ,因?yàn)橐崆坝涗?,還要處理數(shù)據(jù) .1法是對(duì)單元集進(jìn)行賦值 ,與結(jié)點(diǎn)集無(wú)關(guān)常見(jiàn)的問(wèn)題 有結(jié)構(gòu)物的情況如何平衡 ?有結(jié)構(gòu)物的情況下就不是正確的平衡點(diǎn) ,如果一定要這樣平衡 ,比如原來(lái)有一個(gè)隧道 ,現(xiàn)在在上面建房子 ,要估計(jì)這種情況下的追加位移是多少 ,這種情況下可以用上面的方法 2和3來(lái)完成 .多層土的情況下的平衡 多層土的情況可以根據(jù)不同的情況選擇方法 1或2進(jìn)行 .有地下水的情況下的平衡 ?有水的情況下最主要是要搞清楚土的三相 ,按照不同的孔隙率計(jì)算土的干密度和濕密度 ,這種平衡目前最好是只用方法 1來(lái)完成 .邊坡加重力無(wú)法平衡 ?邊坡加重力無(wú)法平衡表示邊坡本身是不穩(wěn)的 ,請(qǐng)檢查模型是否合理 ,塑性參數(shù)是否合理 ./thread-753043-1-12.htmlaba_aba大哥解答：你的inp文件中有 part和instance,而geo.dat中的節(jié)點(diǎn)編號(hào)前面沒(méi)有 instance名稱(chēng)。有兩種方法：。建立沒(méi)有part的inp文件，geo.dat中的單元編號(hào)前面 沒(méi)有instance名稱(chēng)。見(jiàn)/forum/viewthread.php?tid=708942。 建立有part的inp文件，geo.dat中的單元編號(hào)前面 加上instance名稱(chēng)。見(jiàn)/forum/thread-724821-1-1.html"\.R!y%M4g:O!z!y48按照第 2種方法，在你的 geo.dat文件中單元編號(hào)前面加上 Part-1-1.，例如第一行變?yōu)?Part-1-1.1,-1.50E+04,-3.50E+04,-1.50E+04,-1.25E-12,-1.92E-13,-4.17E-12/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)力平衡方法熊志勇陳功奇第一部分 地應(yīng)力平衡方法簡(jiǎn)介地應(yīng)力平衡有三種方法：*initialconditions,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)上的應(yīng)力超出屈服面之后 ,應(yīng)力轉(zhuǎn)移要通過(guò)大量的迭代才能完成 ,而且有可能出現(xiàn)解不收斂的情況。在僅考慮自重情況下只能考慮受泊松比的影響帶來(lái)的側(cè)壓力系數(shù)效應(yīng) ,因此平衡后的效果不一定很理想 ,但無(wú)疑其適用性很強(qiáng)。*initialconditions,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)就是不能夠很好平衡具有起伏表面的幾何模型 ,需知道平整后模型的上覆巖體自重。*initialconditions,type=stress,geostatic,user該方法采用用戶子程序 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ì)大家有用。模型描述：2080kg/m3E11070.3560kpa12m 12m(0,20) (30,20)(0,0)地基(30,0)30m1.2.啟動(dòng) ABAQUS單擊, CreateModelDatebase創(chuàng)建部件（ Part）在Part模塊，單擊創(chuàng)建部件按鈕，彈出如右圖的對(duì)話框，按圖輸入部件名： Part-soil；采用二維模型選擇 2DPlanar；Type選擇可變型（ Deformable）；基本特征選擇殼體（ Shell）；Approximatesize輸入 70，這個(gè)數(shù)值的大小，應(yīng)根據(jù)模型的最大尺寸來(lái)確定：稍大于最大尺寸的 2倍。比如本模型最大尺寸是 30那么我輸入了 70，但也不是絕對(duì)， 你當(dāng)然也可以輸入 65，或75等等。最后單擊 Continue,繼續(xù)下一步。按照模型尺寸 （如圖），建立模型部件，雙擊鼠標(biāo)中鍵，完成部件的建立。3.建立材料屬性（ Property）3.在Module中切換到 Property模塊，單擊 ，輸入材料名稱(chēng)（ name）:Material-soil,單擊Density,在彈出對(duì)話框中輸入：密度 2080；然后單擊Mechaniacal—Elasticity—Elastic,在彈出的對(duì)話框輸入圖，單擊 OK完成材料的定義。Name:Material-soilDescription:MaterialBehaviorsName:Material-soilDescription:GeneralMechanicalThermalOtherMaterialBehaviorsDensityDepyarRegularizationUserMaterialUserDefinedFieldUserOutputVariablesGeneralMechanicalThermalOther1Density1UseterNumbercDataElasticity ?Elastic 「Plasticity ?DamageforDuctileMetals ?DamageforTractionSeparation Laws ?DamageforFiber-Reinforced Composites?HyperelasticHyperfoamHypoelasticPorousElasticIDeformationPlasticityViscoelasticD■DampingExpansionBrittleCracking|Delete▼Suboptions▼SuboptionsType:IsotropicUsetemperature-dependentdataNumberoffieldvariables: 07|Modulitimescale(forviscoelasticity):Long-termNocompressionNotensionDataYoung's PoissonasModulus Ratio1le70.3S單擊按鈕 ，輸入名字： Section-soil,選擇Soild,Homogeneous,單擊Continue，OK,完成截面的創(chuàng)立。單擊按鈕 ，選取部件(單擊或框選，選擇后成粉紅色，表示選中)，單擊 Done單擊按鈕 ，選取部件(單擊或框選，選擇后成粉紅色，表示選中)，單擊 Done或單擊鼠標(biāo)中鍵來(lái)確定。在彈出的對(duì)話框中選中 Section-Soil，單擊OK.裝配部件( Assembly)OK.在Module選擇Assembly模塊,單擊OK.注意：1.本模型只有一個(gè)部件所以自動(dòng)選中，如有多個(gè)可按 Shift鍵全選中；nstanceType本例都無(wú)所謂，只有一個(gè)部件，但是如有多個(gè)部件，我比較喜歡選擇 Independent,因?yàn)樵谶@種情況下，所有的部件會(huì)出現(xiàn)在一個(gè)窗口，劃分網(wǎng)格更方便，并且直觀的看到不同部件連接處網(wǎng)格劃分的是否協(xié)調(diào)。創(chuàng)建分析步（ Step）在Module選擇Step模塊,單擊 ，彈出下面的對(duì)話框， 創(chuàng)建分析步 Step-1,選擇Geostatic,單擊Continue，在分析步編輯框 Basic中選擇On（如果會(huì)發(fā)生大變形的情況下要選擇 On），其他默認(rèn)，單擊 OK

6.施加荷載和定義邊界條件（ Load6.施加荷載和定義邊界條件（ Load）在Module選擇Load模塊，單擊 定義邊界條件，需要對(duì)模型的左、右、下底面定義邊界條件，分別命名 BC-1、BC-2、BC-3，選擇初始步（Initial）,Displacement/Rotation,單擊Continue，選擇左、右、下底面邊界，單擊 Done，分別選擇 U1、U1、U2（U1是水平方向， U2是豎直方向） ，單擊 OK完成邊界條件的定義。

單擊 定義重力荷載， Name:Load-grv,選擇Step-1，Mechanical，Gravity，單擊Continue。單擊EditRegion，選擇整個(gè)模型，在 Component2中輸入 -9.8（重力加速度），單擊 OK力加速度），單擊 OK7.劃分網(wǎng)格（ Mesh）（ 劃分網(wǎng)格是一門(mén)藝術(shù) ，本例只是簡(jiǎn)單的劃分）單擊 ，采用默認(rèn)，單擊 OK。7.劃分網(wǎng)格（ Mesh）（ 劃分網(wǎng)格是一門(mén)藝術(shù) ，本例只是簡(jiǎn)單的劃分）單擊 ，采用默認(rèn)，單擊 OK。在Module選擇 Mesh模塊，單擊 設(shè)定網(wǎng)格的種子，將全局種子大小設(shè)為 1，其余默認(rèn)，單擊 OK.單擊 ，F(xiàn)amily選擇PlaneStrain單擊 ，F(xiàn)amily選擇PlaneStrain，其余采用默認(rèn)，單擊 OK單擊 ，單擊Yes,完成網(wǎng)格的劃分，如下圖：在命令行中輸入 mdb.models['Model'].setValues(noPartsInputFile=ON)，按回車(chē)鍵在job模塊中創(chuàng)建名為 Job-NoInitialCondition應(yīng)力圖水平位移圖將分析得到的應(yīng)力場(chǎng)保存為一個(gè)文本文件。具體做法 :打開(kāi)分析得到的 ODB文件，選擇菜單 Report→FieldOutput,在下圖所示的對(duì)話框中，選中積分點(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ì)話框中的 Setup標(biāo)簽頁(yè)，在 Name文本框中輸入要保存的文件名 b.inp，取消對(duì) Appendtofile的選擇（即創(chuàng)建一個(gè)新文件） ，在Write后面只選中 FieldOutput（如下圖所示） 。圖表2設(shè)置輸出場(chǎng)變量注意，此處輸出的是當(dāng)前增量步結(jié)束時(shí)的應(yīng)力結(jié)果，因此 上述對(duì)話框頂部的 Step必須是Geostatic分析步，F(xiàn)rame必須是1。如果Frame是 0，會(huì)看到輸出的應(yīng)力都是 0按照 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），只