基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】

2017-12-07  by:CAE仿真在線  來源:互聯(lián)網(wǎng)


【問題描述】:

某地鐵盾構(gòu)隧道管片襯砌內(nèi)徑為5.4m,外徑為D=6m,埋深為2D。從上至下,根據(jù)土層的物性參數(shù)不同將其分為3層,各層的材料參數(shù)和層厚為:

第1層:厚8m,E=3.94Mpa,v=0.35,ρ=18.28kN/m^3

第2層(隧道所在層):厚18m,E=20.6Mpa,v=0.3,ρ=20.62kN/m^3

第3層:厚15m,E=500Mpa,v=0.33,ρ=21.6kN/m^3

施工中掘削面頂進壓力為0.3MPa,盾尾注漿壓力為0.15MPa。

試采用ANSYS模擬此過程。

【建模要點】:

1、建模過程充分使用對稱性建模的方便,使用到的對稱性命令為 arsym

2、網(wǎng)格劃分輔助mesh200的使用,建模思路為通過建立面,采用mesh200劃分面,拉伸面成體,從而形成實體單元。

3、注意在第2步采用面拉伸成體單元后,體單元材料屬性的重新賦值。

4、自重應(yīng)力場的求解。

5、利用重啟動以及生死單元來模擬盾構(gòu)掘進的過程。

【建模過程】:

1、首先建立隧道附近的四分之一模型,注意網(wǎng)格的局部細(xì)分。

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)果圖片1

2、利用對稱性,建立二分之一隧道模型,并建立隧道上方和下方土體模型。

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)果圖片2

3、利用對稱性,建立整個隧道平面模型

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)果圖片3

4、利用面拉伸成體的思路,通過輔助單元建立實體單元,這里實體單元采用soild186進行模擬。

注意拉伸時的一個額外命令的使用:

extopt,aclear,1

該命令意思也即是在拉伸完成后刪除母體單元mesh200

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)果圖片4

5、由于在拉伸時候都是默認(rèn)的材料號為1,拉伸完成后需根據(jù)不同的位置,選擇不同的土體進行材料參數(shù)的改變。

6、約束條件的設(shè)置,本次約束取土地地面為全約束,各側(cè)邊約束為平行法向方向固定約束。頂面除四周邊界線有約束外,其余地方皆無約束。整個示意圖如下:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)果圖片5

7、自重應(yīng)力場的求解。

8、利用重啟動和單元生死功能,模擬盾構(gòu)掘進過程。


【結(jié)果查看】

1、自重應(yīng)力場求解后的相關(guān)結(jié)果

Y方向應(yīng)力云圖:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)果圖片6

位移云圖:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)果圖片7

2、地層位移分析

地層位移本例是指相對于自重固結(jié)下的位移,因而在分析地層位移時需要采用荷載工況的手段減去第一步在自重計算下的結(jié)構(gòu)位移。例如在計算開挖第一步時的計算命令流如下:

/POST1

LCDEF,1,1, ,

LCDEF,2,3, ,

LCASE,2,

LCOPER,SUB,1, , ,

開挖第一步的位移云圖如下:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys分析案例圖片8

開挖第四步的位移云圖如下:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys分析案例圖片9

開挖第七步的位移云圖如下:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys分析案例圖片10

隧道貫通后的位移云圖如下:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys分析案例圖片11

從結(jié)果可見,地層位移的最大值出現(xiàn)在隧道的拱頂和仰拱處,其中拱頂產(chǎn)生最大下沉約4.5cm,仰拱處產(chǎn)生最大上隆約4.9cm。

3、地表沉降分析

通過選擇地表表面的節(jié)點和單元,可以繪制地表的沉降分布云圖。標(biāo)準(zhǔn)的命令流格式如下:

CSYS,0

WPCSYS,-1

NSEL,S,LOC,Y,15 !選擇地表節(jié)點

ESLN,S !選擇附在節(jié)點上的單元

EPLOT

PLNSOL,U,Y,0,1

通過更改工況號,即可得到不同工況作用下的地表沉降情況,下圖分布表示了第二步地層位移分析對應(yīng)的地表沉降情況。

開挖第一步地表沉降位移云圖:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys分析案例圖片12

開挖第四步地表沉降位移云圖:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys分析案例圖片13

開挖第七步地表沉降位移云圖:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)構(gòu)分析圖片14

隧道貫通后的地表沉降位移云圖如下:

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)構(gòu)分析圖片15

通過云圖可以看出,在盾構(gòu)機的頂進作用下,地表發(fā)生了前隆后沉的現(xiàn)象。隨著盾構(gòu)機的推進,開挖面后方的沉降槽越來越寬,最大沉降量也越來越大,隆起區(qū)域也逐步向前推進。地表最大沉降量約為1.2cm,最大隆起量約為0.3cm。

4、管片結(jié)構(gòu)受力分析

進行管片受力分析時,不需要進行工況組合,直接將沒每步的計算結(jié)果讀入內(nèi)存即可。如繪制第一步開挖完畢后管片襯砌的Mises等效應(yīng)力云圖的命令流為:

VSEL,S,LOC,X,2.3,2.7

VSEL,R,LOC,Z,0,-3

ESLV,S !選取第1環(huán)管片

SET,3,LAST,1,

PLNSOL,S,EQV,0,1

基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】ansys結(jié)構(gòu)分析圖片16

通過選取不同的管片襯砌,及讀入不同的計算步,可以計算出各個開挖步結(jié)束時的管片襯砌的等效Mises應(yīng)力云圖。


開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:基于ANSYS某地鐵盾構(gòu)隧道掘進過程數(shù)值模擬分析【轉(zhuǎn)發(fā)】 Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓(xùn) 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進行在線報名,我們會及時給您回復(fù)電話,謝謝!
驗證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點擊咨詢 點擊咨詢
項目QQ咨詢:點擊咨詢
email:kf@1cae.com