這個是前幾天做的一個小項目，具體工程資料不便給出。貼出來模擬方法和思路以及處理結果，供大家參考。

該模型為二維，將盾構穿越節(jié)點處簡化為一個 平面應變問題。其中開挖邊緣具最近的一根樁為1m。整個模型分為四個part——土體、等代層、襯砌、管樁。因上部橋梁為簡支且無法簡化為平面應變問題，因此將其忽略。





其中，等代層厚度的取值大家可以隨便找來一篇文獻，或是看去年新出的《abaqus在隧道及地下工程中的應用》一書。全部的part都采用實體模擬。
在設置接觸的時候，要根據不同的問題來判斷采用何種接觸。本模型所關注的問題是 盾構在穿越時樁基及其上部結構的力學響應問題，因此，涉及盾構的接觸，完全可以用綁定約束，而管樁與土體的約束，自然就不能這樣設置了。
不考慮地下水滲流和固結。其它都設為線彈性。下面是網格示意圖。

施工過程是，右側盾構先行開挖，左側盾構再行開挖。 開挖的同時加上等代層和襯砌單元。
首先是 地應力平衡的問題。地應力平衡應該是一個比較成熟的技術了，但有些剛開始接觸aba的童鞋可能還不太善于使用論壇搜索功能，總是發(fā)帖問。其實以前有很多人已經做過詳細的闡述了，很多情況下直接搜索然后依葫蘆畫葫蘆就可以了。
對于這種樁基事先存在，且有較繁瑣的接觸的情況，大家可以采用最普通的平衡方法，即加載后導出 形心點應力場。

要注意，第一次平衡可采用initial conditions，type=stress,geostatic的方式，這時候step-1最好設置為static進行平衡，并且步長一定要小。計算完成后導出應力場，然后將原命令改為 initial conditions，type=stress, input=xx.csv進行平衡。如此反復進行，當用static無法計算下去時，可以改為geostatic進行平衡。


再次強調一下：initial conditions,type=stress,geostatic下的命令行，其默認重力方向為Z（二維：Y）軸，因此最好在cae中將方向轉過來再生成inp。在assembly中選定全部模型旋轉即可。
最后一般都可以平衡到滿意精度，如果一次不行，可在導出應力場前把精確位數擴大。如下圖所示。



具體的操作方法這里不再贅述了，因為已經是老生常談了。
下面是經過多次反復平衡后的效果：


平衡效果在10e-4m量級，對于工程來說綽綽有余。