搜索ABAQUS與巖土工程（一）——地應(yīng)力平衡（初級）【轉(zhuǎn)載】
2016-11-14 by:CAE仿真在線 來源:互聯(lián)網(wǎng)
轉(zhuǎn)自:算盤坊小胖微信號
地應(yīng)力是個很神奇的話題。很多寶寶對這個概念不是很明了——誰沒有個年輕的時候呢?對吧。首先我們來回顧一下啥叫地應(yīng)力。小胖不確定這種有點理論性的東西能不能講得準(zhǔn)確。如果您認為我講得不對,請及時溝通。在確認的情況下,我會及時更正,以免誤人子弟。
(一)地應(yīng)力的含義
以前混論壇的時候,經(jīng)常會有非巖土地質(zhì)背景的寶寶問:“巖土真事逼,從來沒聽說過什么應(yīng)力平衡啥的。”有的初級巖土寶寶也會有疑問啊,問“啥叫地應(yīng)力平衡啊?我直接加個重力在模型上不就好了嗎?地應(yīng)力平衡是不是只有考慮塑性的時候才需要呀?ABAQUS真裝B啊,FLAC就沒這回事啊?”Blah Blah......對這些人(當(dāng)然包括年輕時的小胖),小胖想說:你辣么NB,你咋不上天呢?
我們先來看看ABAQUS它們家怎么說的
?
基本意思人家在自己的幫助文件里也交代清楚了,算是仁至義盡了。咋說的呢?就是寶寶們在算巖土問題的時候呀,如果直接加上重力等外力,模型的位移會過大,因此需要先賦予模型一個初始的應(yīng)力同時保證變形不過大。
明白了么?
其實小胖年輕時也一直不是很明白,總覺得有那么一層窗戶紙沒通透,心里不夠亮堂。現(xiàn)在我大體能想明白了。各位寶寶可以想象一下三軸試驗。理想狀態(tài)下試件不受重力,那么初始狀態(tài)自然不必定義初始應(yīng)力,也就是內(nèi)力,也就不需要所謂的地應(yīng)力平衡。如果考慮重力,那么重力帶來的影響遠遠小于后期加載對試樣的影響,因此也可以忽略不計。
但,全尺度的巖土問題就是另一回事了。
在地球母親的吸引力下,巖土介質(zhì)的內(nèi)力是很大的,尤其是埋深較深的時候。這時候如果有卸載,自然會影響到整體的變形和位移,這部分因素是需要考慮的。同時,巖土介質(zhì)又是一種應(yīng)力路徑相關(guān)的材料,不同的應(yīng)力路徑會造成材料不同的力學(xué)性質(zhì),因此初始的地應(yīng)力也是需要考慮的。那么年輕時的小胖就問了:直接加個重力不就完了嗎?
兩碼事。
當(dāng)前的巖土介質(zhì)中的應(yīng)力,是在長期重力作用下逐漸固結(jié)沉降后形成的。其初始的變形狀態(tài)并不是你現(xiàn)在所建模型的狀態(tài)。因此,地應(yīng)力平衡的真正目的,只是為了塑造一個基本符合事實的狀態(tài)。你平時腳底下踩的土,就是你現(xiàn)在模擬的初始態(tài),而并不是你虛擬的那個零變形零應(yīng)力的土體。
話有點繞口,但基本就是這樣子。
現(xiàn)在可以回答有些寶寶的問題了:為啥只有ABAQUS有地應(yīng)力平衡?
小胖的回答是
(二)地應(yīng)力平衡方法論——最基本方法(geostress1-m1.cae)
現(xiàn)在開始學(xué)習(xí)ABAQUS的小盆友,乃們其實應(yīng)該很有幸福感了。用過6.5/6.6的老寶貝兒應(yīng)該都能明白當(dāng)年的苦。現(xiàn)在ABAQUS中的主流方法應(yīng)該是用Geostatic分析步了。Geostatic分析步有兩個選擇,有些同學(xué)傻傻分不清楚。
看到?jīng)],一個是“Fixed”,一個是“Automatic” 。我們還是來看看達索他們家是怎么說的:
各位寶寶看到了?達索他們自己說,這個automatic啊,可以根據(jù)你定義的一個最大位移來進行地應(yīng)力平衡。最大位移的定義,inp文件里就是加個“utol=displacement tolerance”。在CAE里,就是在上邊那個圖的“Max.displacements changes”里面進行設(shè)置。
?那問題來了:這個fixed是啥意思呢?
那我們做個例子吧。
小胖簡單做了個實例,放在網(wǎng)盤里了,大家可以按照文章開頭的方法獲取。
一個邊長10m的立方體,前后左右施加位移邊界,底部施加豎向位移邊界,模型頂部自由。暫時先不考慮水,因此選用C3D8R單元。材料選用線彈性加Mohr-Coulomb塑性準(zhǔn)則。c=10kPa, phi=10°,E=10MPa,泊松比0.35,密度1.8。
我們采用CAE內(nèi)部已有的預(yù)置地應(yīng)力的方法,在Load-Predefined field-Mechanical-Geostatic stress,如下圖
各行的含義aba寶寶們應(yīng)該了然于胸了吧。注意你的地應(yīng)力方向的正負就好了:本模型中和Z正向相反,因此第二點地應(yīng)力為負。靜止側(cè)壓力系數(shù)按照v/(1-v)進行計算。注意,軟件只是按照你給定的側(cè)壓力系數(shù)計算初始地應(yīng)力,一旦計算開始后,這個參數(shù)也就沒用了,豎向和橫向的應(yīng)力關(guān)系靠泊松比來控制。
然后!
我們選擇了這種初始地應(yīng)力賦值模式,那么就相當(dāng)于默認地應(yīng)力是加在Z方向的,所以需要在assembly里旋轉(zhuǎn)一下模型!不信的寶寶可以看下不旋轉(zhuǎn)所賦的地應(yīng)力和平衡效果:
從初始地應(yīng)力S33的分布可以看出,軟件是默認Z方向為地應(yīng)力方向的,再看下平衡后的位移,雖然量值都很小,但分布規(guī)律很感人啊。這個錯誤的文件我命名為geostress1-error放在網(wǎng)盤里。
正確的姿勢應(yīng)該是在assembly-instance-rotate里面對模型進行旋轉(zhuǎn),以使重力方向沿Z軸方向 (當(dāng)然,寶寶們在part里建模時也可以直接考慮Z軸的重力方向)。我們來看下正確姿勢下的平衡效果
我們再來看下Fixed選項下的平衡效果
總覺得怪怪的有沒有。我們看下如果不用Geostatic而用Static, General,計算結(jié)果如何
結(jié)論是:一樣滴。因此,盡管我們現(xiàn)在看起來平衡效果也還不錯,但這只是因為我們的模型是最簡單的一種。稍微復(fù)雜一點的話,Gestatic(Fixed)也不會得到好的效果,和Static, General其實是一回事。
寶寶們可以自己試一下。
(三)地應(yīng)力平衡方法論——基本方法(geostress1-m2.cae)
Geostatic方法經(jīng)常搞不定稍微復(fù)雜些的模型。此時,寶寶們可以選擇另一種更普適的方法——不預(yù)設(shè)地應(yīng)力,直接加重力先獲取一個應(yīng)力,再重建一個模型導(dǎo)入該應(yīng)力作為初始地應(yīng)力進行計算。
首先計算無初始地應(yīng)力的情況下的模型,采用Static,General(Geostatic也是可以的)分析步,我們看計算完成的位移云圖
可以看出位移比較大。如果按照某些寶寶的意見不進行平衡,那么后期的開挖等施工就直接在這個有如此明顯的初始位移的模型上進行,那和我們實際要求的就不是一回事了。
然后我們繼續(xù)進行平衡,以這個odb文件應(yīng)力為初始應(yīng)力,我們在cae里新建一個job,同時在predefined field-stress中設(shè)置初始地應(yīng)力(小胖現(xiàn)在已經(jīng)基本不使用手動導(dǎo)出地應(yīng)力csv文件再導(dǎo)入了,太麻煩了),見下圖
Continue后,選取整個模型,見下圖
Specification里我們選擇from output database file,也就是剛才的odb文件
剛才的Job-1我們只進行了一個增量步的計算,因此step選1(一般都是選擇1),increment選1。然后OK。
然后新建一個job進行計算。我們看下計算結(jié)果:
最大位移只有1.3e-5 m,還算不錯。小胖一般平衡到0.5 mm以下就滿足了,因為這個精度在實際工程監(jiān)測時都不見得能達到。
注意這種平衡方法并不適合多孔介質(zhì)弾性模型(因此也就不適合臨界狀態(tài)塑性,因為這兩者是綁定的),后續(xù)我還會繼續(xù)介紹多孔介質(zhì)弾性模型下的地應(yīng)力平衡,地表不平整情況下的平衡以及有建構(gòu)筑物時的地應(yīng)力平衡,敬請期待。
??
鎮(zhèn)樓啦啦啦
(Johnson, J. B., et al. (2015). "Discrete element method simulations of Mars Exploration Rover wheel performance." Journal of Terramechanics 62: 31-40.)
————誰說巖土是夕陽產(chǎn)業(yè)呢?
?
相關(guān)標(biāo)簽搜索：ABAQUS與巖土工程（一）——地應(yīng)力平衡（初級）【轉(zhuǎn)載】 abaqus分析培訓(xùn) abaqus技術(shù)教程 abaqus巖土分析 鋼筋混凝土仿真 abaqus分析理論 abaqus軟件下載 abaqus umat用戶子程序編程 Abaqus代做 Abaqus基礎(chǔ)知識 Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓(xùn)
