地應(yīng)力反演的案例
地應(yīng)力反演相互合作交流
地應(yīng)力反演依據(jù)工區(qū)測點地應(yīng)力數(shù)據(jù)反演應(yīng)力場，可以相互交流合作。
請教地應(yīng)力反演問題
在flac3d中做地應(yīng)力反演一般施加是什么樣的標(biāo)準(zhǔn)荷載進行反演？
2。如何知道反演后測點的應(yīng)力值與實測點的應(yīng)力值一直，及如何取得某一點的應(yīng)力值？
3。在flac3d中怎樣才能進行子模型分析？
這幾個問題可能有些提法不是很好，由于flac3d還不是很了解，希望大家?guī)蛶兔鉀Q。現(xiàn)急待解決實際問題。謝謝！
有意交流者發(fā)郵件britht-h@163.com
三軸壓縮，考慮了初始地應(yīng)力和沒考慮初始地應(yīng)力的區(qū)別
考慮了初始應(yīng)力，加載的是力加載
沒有考慮初始應(yīng)力，加載的是位移加載
’
ABAQUS熱應(yīng)力分析 附ABAQUS中初始地應(yīng)力的施加下載
軋輥與Cu層的熱傳導(dǎo)系數(shù)
下載地址：ABAQUS中初始地應(yīng)力的施加
APDL地應(yīng)力平衡 ￥8.9
APDL地應(yīng)力模型，將樁的單元殺死后對土壤施加重力加速度，命令流見附件
ABAQUS地應(yīng)力平衡
對于涉及開挖、填充的動態(tài)巖土工程問題，初始應(yīng)力場是正確模擬其施工過程的先決條件。初始應(yīng)力的加載必須滿足地應(yīng)力平衡，而地應(yīng)力平衡其實就是為了使數(shù)值模擬獲得一個存在初始應(yīng)力，而無初始應(yīng)變的狀態(tài)。當(dāng)考慮自重場是產(chǎn)生地應(yīng)力場的主要因素時，顯然重力是外力，初始應(yīng)力場是內(nèi)力，將提取出的內(nèi)力施加于數(shù)值模擬后再施加重力，是內(nèi)力和外力平衡，從而獲得較精確的沒有受到人為因素干擾情況下的數(shù)值模型的初始狀態(tài)。
在 ABAQUS 有限元軟件中，提供了 5 種定應(yīng)力平衡方法，這里我就比較常用的，使用比較普遍的一種進行講解。用的例子還是上節(jié)ANSYS的例子。
接下來，小編就一個簡單的算例就行地應(yīng)力平衡。在一個均質(zhì)的地基上有一根立柱。
柱子的尺寸為1m×2m×10m，密度2400kg/m3，彈性模量2.8E8Pa（這里我為了方便觀察規(guī)律，降低這個彈性模量），泊松比0.167；均質(zhì)地基的尺寸為30m×30m×15m，密度為2000kg/m3，彈性模量3E7Pa，泊松比0.3。模型見圖1.
1、將模型導(dǎo)入ABAQUS中。
2、創(chuàng)建靜力荷載步
3、殺死柱子單元
4、施加約束（四周施加法向連桿約束，底部施加全約束）
5、施加自重荷載
6、提交計算
7、查看后處理結(jié)果
豎向位移
豎向應(yīng)力
8、應(yīng)力場輸出
9、在工作目中找到stress_out.rpt文件
10、經(jīng)過格式的轉(zhuǎn)化（文件的后面不要有空行）
11、在之前的inp文件中Step-1前面加上*Initial conditions,Type=stress,Input=stress_in.rpt這句命令，記得一定要把stress_in.rpt放在工作目錄里面。
展開 ANSYS地應(yīng)力平衡方法
進行地應(yīng)力平衡的原因如下：我們建立的幾何模型一般都和工程實際情況一致，例如邊坡的幾何模型與邊坡實際尺寸相一致。但是由于邊坡的沉降和徐變作用，可以想像到，現(xiàn)在的邊坡應(yīng)該是由一個體積更大的原始邊坡在很久以前由于受到重力作用和邊界約束條件，逐漸形成了現(xiàn)今的邊坡形態(tài)。但是對于那個原始的邊坡形態(tài)，我們不得而知。假如能準(zhǔn)確知曉，我們就能夠建立原始邊坡的幾何模型，接著對邊坡施加重力和邊界條件，受力后邊坡形態(tài)應(yīng)該和現(xiàn)在的邊坡相一致，其內(nèi)力就是初始應(yīng)力場（地應(yīng)力），這樣就不用專門施加地應(yīng)力了。但現(xiàn)實情況是我們不能知曉原始邊坡的形態(tài)。現(xiàn)在的邊坡幾何模型就是其實際形態(tài)，受力之后將會變成一個與現(xiàn)狀不一致的邊坡，這不符合現(xiàn)在的實際情況。如果我們計算出現(xiàn)今邊坡的內(nèi)力，并將其作為邊坡的初始應(yīng)力場，再去和外力平衡，這樣我們建立的模型就和現(xiàn)實邊坡情況相一致了。
對于涉及開挖、回填的動態(tài)巖土工程問題，地應(yīng)力平衡是正確模擬施工過程的前提條件。初始應(yīng)力的加載必須滿足地應(yīng)力平衡，而地應(yīng)力平衡就是為了使地基僅存在初始應(yīng)力，而不存在初始應(yīng)變。當(dāng)?shù)鼗灾厥钱a(chǎn)生地應(yīng)力場的主要因素時，重力是外力，初始應(yīng)力場是內(nèi)力，將提取出的內(nèi)力施加于模型后再施加重力，此時內(nèi)力和外力平衡，該狀態(tài)就是工程建設(shè)的初始狀態(tài)。
在ANSYS中的地應(yīng)力該怎么平衡呢？
初應(yīng)力(Initial Stress)可以指定為一種“荷載”進行施加,但僅在靜態(tài)分析和全瞬態(tài)分析中使用，可以用于線性分析或非線性分析。
展開 多孔介質(zhì)的地應(yīng)力平衡
根據(jù)單向固結(jié)理論，有效應(yīng)力=總應(yīng)力-孔隙壓應(yīng)力。在非飽和問題中，ABAQUS采用的是：有效應(yīng)力=總應(yīng)力-飽和度*孔隙水壓力。又由于ABAQUS對孔壓正負(fù)號的規(guī)定和土力學(xué)中一致，而對應(yīng)力正負(fù)號的規(guī)定和土力學(xué)中剛好相反，故上面公式中的“-”號要改為“+”號。對《abaqus在土木工程中的應(yīng)用》一書中的這個例子，土體頂面的總應(yīng)力是為零的，但孔壓為-10kpa，查-10kpa對應(yīng)的飽和度得0.9789，再用上面的公式計算就得到了-9.789kpa。個人感覺不要死套公式，對每一點都分別算出總應(yīng)力和孔隙水壓力，再按上式計算有效應(yīng)力。最后一點：《abaqus在土木工程中的應(yīng)用》一書中給出的初始有效應(yīng)力的計算公式有點問題，手冊上面說的很詳細，手冊里面是把土體分為完全干燥區(qū)，半飽和區(qū)和完全飽和區(qū)三部分來考慮的。把手冊上相關(guān)的部分傳上來，希望對大家有所幫助。
多孔介質(zhì)的地應(yīng)力平衡.rar
展開 地應(yīng)力平衡方法以及注意事項
注意：只有采用彈塑性本構(gòu)模型時需要地應(yīng)力平衡，彈性本構(gòu)不需要地應(yīng)力平衡！
第一部分 地應(yīng)力平衡方法簡介
地應(yīng)力平衡主要有五種方法：
（1）自動平衡：第一步創(chuàng)建分析步geostatic ，這種方法注意只能在第一步只能有土和重力的情況下能使用，有其他部件或者接觸時計算不能收斂，效果是最好的，方便簡單！
（2）*initial conditions,type=stress,geostatic
該方法需給出不同材料區(qū)域的最高點和最低點的自重應(yīng)力及其相應(yīng)坐標(biāo)。所采用的幾何模型一般較規(guī)則,表面大致水平,地應(yīng)力平衡的好壞一般只受巖體密度的影響,無論采用彈性或彈塑性本構(gòu)模型都能很好的達到平衡,可以不必局限于僅受泊松比的影響,能夠通過考慮水平兩個方向的側(cè)壓力系數(shù)值來施加初始應(yīng)力場。計算速度快,收斂性好。缺點就是不能夠很好平衡具有起伏表面的幾何模型,需知道平整后模型的上覆巖體自重。
高版本在CAE里也能操作
用計算器算出每個分界面上的應(yīng)力和坐標(biāo)對應(yīng)填入，也比較方便 不需要修改關(guān)鍵字
（3）*initial conditions,type=stress,geostatic,user
導(dǎo)入ODB里的方法，也比較簡單，高版本可在截面上操作，不需要修改關(guān)鍵字
你放入ODB后，填入第一步不需要填名稱 就是1 ，增量步就是你第一步計算的最后一個增量步
（4）*initial conditions,type=stress,input=FileName.csv(或inp)
該方法中的文件FILENAME.INP獲取方法為:首先將已知邊界條件施加到模型上進行正演計算,然后一般是將計算得到的每個單元的應(yīng)力外插到形心點處并導(dǎo)出6個應(yīng)力分量(也可以導(dǎo)出積分點處的應(yīng)力分量,視要求平衡的精確程度而定)。
展開 有限元特輯II 初始地應(yīng)力平衡
Motivation
在 ABAQUS中，提供了 5 種定應(yīng)力平衡方法，分別是
1. ( AUTOBALANCE) 自動平衡法;
2. * INITIAL CONDITIONS，TYPE = STRESS，GEOSTATIC;
3. * INITIAL CONDITIONS，TYPE = STRESS，F(xiàn)ILE = file，INC =inc;
4. * INITIAL CONDITIONS，TYPE = STRESS，INPUT =XX． DAT;
5. * INITIAL CONDITIONS，TYPE = STRESS，GEOSTATIC，USER。
以上 5 種方法并不是每一種適用于所有的巖土模型，方法從易到難。
方法①的自動平衡法，它省去了自重應(yīng)力以及生成相應(yīng)初應(yīng)力文件和導(dǎo)入的麻煩。在( GEOSTATIC) 地應(yīng)力中選擇自動增量步就能使用自動地應(yīng)力平衡功能，還能指定允許的位移變化容限。不過自動地應(yīng)力平衡支持有限的幾種材料如彈性，塑性等，而其起單元也有一定的要求。
方法②為關(guān)鍵字定義初始地應(yīng)力法，這種方法需要給出不同材料區(qū)域的最高和最低點的自重應(yīng)力及其相應(yīng)坐標(biāo)。所采用的幾何模型一般較規(guī)則，表面水平，能夠通過考慮水平兩個方向的側(cè)壓力系數(shù)值來施加初始應(yīng)力場。關(guān)鍵字定義初始地應(yīng)力法只適合土體表面水平的土體，而初始地應(yīng)力提取法由于外插的應(yīng)力有一定的誤差，因此對于材料是彈塑性的復(fù)雜土體，應(yīng)力轉(zhuǎn)移要通過大量的迭代才能完成，而其有可能出現(xiàn)不收斂的情況，平衡效果可能不是很理想。
方法③是 ODB 導(dǎo)入法，這種方法可使用之前算過的 ODB 文件結(jié)果，也就是說提前計算一個初始應(yīng)力ODB 文件，定義初始應(yīng)力時直接指定 ODB 文件即可。
展開 abaqus地應(yīng)力平衡的幾個注意事項
首先大家來比較一下采用geostatic分析步時增量步分別為automatic與fixed計算出來的孔隙水壓力、豎向有效應(yīng)力、豎向位移結(jié)果對比，可以發(fā)現(xiàn)，除了豎向位移不同，計算出來的孔隙水壓力、豎向有效應(yīng)力完全相同，可以認(rèn)為增量步為automatic達到平衡，fixed計算出來的豎向位移未有達到平衡。
于是我再采用geostatic分析步時增量步fixed計算出來的odb文件，添加到預(yù)定于場里，可以發(fā)現(xiàn)，除了豎向位移不同，計算出來的孔隙水壓力、豎向有效應(yīng)力完全相同，此時的位移達到平衡。
注意：此模型全部只有土體一個部件，所以用geostatic分析步時增量步automatic能計算出較好的結(jié)果，但是若有其他樁體墻體等部件或者接觸時，geostatic分析步時計算基本不會收斂，希望高手能給予交流，若有其他部件時，平衡真的很難，我平時的做法是采用一般靜力分析步算出odb結(jié)果文件，再添加到預(yù)定于場中，采用一般靜力分析步，該方法大多數(shù)時候平衡效果并不是很好，負(fù)4次方級算理想結(jié)果，只有重復(fù)導(dǎo)入，重復(fù)次數(shù)越多，達到滿意結(jié)果，但是需要耐心，因為不知道會要多少次，甚至幾十次上百次了。
還有就是想說的是，如果你用geostatic分析步時增量步automatic計算時，如果有其他壓力或者其他荷載時，也是很難收斂，甚至基本是無法收斂。以上暫時也沒什么解決的方法，大家地應(yīng)力平衡盡量選擇用geostatic分析步時增量步automatic，如果有其他部件時就用odb或者csv方法，慢慢調(diào)試吧，同時希望大神能提出更好的方法！
展開 
網(wǎng)上找的別人的心得-地應(yīng)力平衡方法
：）這個功能讓基坑開挖、隧道開挖等的初始應(yīng)力，開挖后的殘余應(yīng)力很好的顯示；也可以很好的模擬鐵路設(shè)計中的工后沉降的概念，在地應(yīng)力平衡后，加上荷載所得沉降即為工后沉降；也很好的模擬了樁土復(fù)合地基的問題，如果沒有初始應(yīng)力的模擬，使土對樁產(chǎn)生了擠壓應(yīng)力，從而通過設(shè)定摩擦系數(shù)就可以模擬了樁與土之間的摩擦力；除此之外，在進行擋土墻計算時也需要ABAQUS的這項功能，反正很多都用得著。
希望我的建議能幫助初學(xué)者（其實我也是初學(xué)者，呵呵）很快掌握這個技巧。
地應(yīng)力分析步中的參數(shù)設(shè)置
*GEOSTATIC分析步幾個參數(shù)（關(guān)鍵還是箭頭所指的參數(shù)），很多人不理解，現(xiàn)在找到幫助文檔，給大家共勉，一起學(xué)習(xí)討論。
箭頭參數(shù)設(shè)置：
箭頭所指的參數(shù)：
Set this parameter equal to the tolerance for the maximum change of displacements. Abaqus/Standard will ensure that the maximum absolute value of a displacement at a node is smaller than the tolerance times the characteristic element length in the model. If this parameter is used without any value specified, the default value of 10–5 is used. If this parameter is omitted, no restrictions are imposed on the displacement values.
設(shè)置此參數(shù)等于容忍位移的最大變化。 ABAQUS/標(biāo)準(zhǔn)將確保位移的在一個節(jié)點的最大絕對值大于公差倍模型中的特征元件長度小。如果沒有指定的任何值使用此參數(shù)，則使用10-5的默認(rèn)值。如果省略該參數(shù)，則不作限制的位移值。
時間增量步設(shè)置：
1. Initial time increment. This value will be modified as required. If this entry is zero or is not specified, a default value that is equal to the total time period
展開 ABAQUS地應(yīng)力平衡（standard）-告別笨拙的導(dǎo)入法 ￥5
此帖僅講述ABAQUS/standard中的地應(yīng)力平衡，explicit大變形的地應(yīng)力平衡另開貼。
該方法適合絕大多數(shù)問題，同時也分享一些小技巧
以下通過預(yù)埋管樁軸對稱模型講解具體實現(xiàn)：
土體高10m(頂面y=0,底面y=-10)，土體浮密度0.65g/cm3，網(wǎng)格單元為CAX4R。
模型圖
結(jié)果分析：
重力平衡后，土體沉降在1×10-5左右，豎向應(yīng)力分布準(zhǔn)確，效果良好。當(dāng)然，如果花時間優(yōu)化模型，可以獲得更好的平衡效果。
關(guān)于ANSYS中初始地應(yīng)力施加方法的介紹
近日，水哥收到一粉絲對初始地應(yīng)力施加這塊的疑問，恰逢今天時間較多，便說說在ANSYS中如何施加初始地應(yīng)力。
針對巖土工程相關(guān)的分析而言，初始地應(yīng)力這個概念比較重要，所謂初始地應(yīng)力，也即是在我們對巖土進行任何外部操作之前，例如基坑開挖、邊坡開挖、隧道開挖等，其本身內(nèi)所存在的真實應(yīng)力，也可稱之為初始應(yīng)力場。
初始應(yīng)力場是平衡的，這也是經(jīng)常聽到的一個概念，初始地應(yīng)力平衡，更簡單來講，就是我們在進行分析之前的位移清零，應(yīng)力不清零。
為什么進行這樣一步操作？
答：為了使模型更加的符合實際。
ANSYS中對于初始地應(yīng)力的平衡沒有類似設(shè)計軟件那般（例如Midas NX）方便，也即是我們在做類似基坑開挖分析之前，首先要進行地應(yīng)力平衡操作，ANSYS中主要分為兩步進行：
第一步：原始應(yīng)力場計算，導(dǎo)出地應(yīng)力文件；
第二步：新建模型，導(dǎo)入地應(yīng)力文件，施加重力，平衡地應(yīng)力。
關(guān)于導(dǎo)入與導(dǎo)出的命令流，ANSYS以前老版本是采用Iswrite與isfile命令，新版本可采用Inistate命令，但是Iswrite和isfile依然可以用，只是幫助文件已經(jīng)沒了這兩個命令的解釋，其用法同inistate，可具體查詢Help。
通過上兩步操作，能達到位移清零，真實應(yīng)力不清零的效果，下面以一個小例子來進行說明操作過程。
某二維地塊，長度50m，高度20m，需進行基坑開挖操作，操作之前，需進行初始地應(yīng)力的平衡，材料為中風(fēng)化砂質(zhì)泥巖，彈性模量取1200Mpa,粘聚力取450Kpa，摩擦角取30度，采用DP材料模型。
展開 