ansys隧道開挖模擬的案例
隧道開挖三維模擬-1
1、引言
早期的隧道開挖模擬常假設為二維平面應變問題，為了反映隧道開挖面引起的應力松弛，學者們提出了不同的處理方式，如convergence-confinement method（或Stress relief method）（圖1），progressive softening method（圖2）, gap method（圖3）等。這些方法沒有完整體現開挖的三維效應。
圖1
圖2
圖3
隧道開挖的三維模擬一般通過“step-by-step”的形式實現，即首先無支撐下開挖一段土體（單元移除），然后在相應位置設置襯砌管片，然后進行下一步開挖。近年來也有學者嘗試完整模擬真實的開挖過程，如開挖機具（盾構）推進、灌漿等https://doi.org/10.1002/nag.395)。本系列首先介紹step-by-step分步三維開挖模擬方法。
2、算例要點
（1）考慮彈性模量隨深度變化（Gibson地基）：利用材料參數的溫度相關性。
（2）利用蒙皮（skin）模擬襯砌。
（3）利用單元生死模擬隧道施工。
3、算例概況
該算例改編于Franzius J N. Behaviour of buildings due to tunnel inducedsubsidence[D]. Department of Civil and Enviromental Engineering. ImperialCollege of Science, Technology and Medicine. London, SW7 2BU, 2003.利用對稱性，模型如圖4所示，在地面以下30.5m開挖一個直徑4.75m的隧道，開挖沿x負向進行。
展開 abaqus隧道開挖數值模擬
abaqus隧道開挖數值模擬
隧道開挖三維模擬 ￥40
模型介紹：
隧道開挖三維模擬，首先進行地應力平衡，而后進行第一次、第二次土體開挖，并用襯砌進行土體支撐。
相關模型見付費內容。凡購買本案例的朋友針對收費內容部分有疑問，可以一起交流。
隧道開挖模擬-2
隧道開挖引起土體隆起不一定都是莫爾庫倫模型的錯，接觸面的激活也有影響。
鏈接：https://pan.baidu.com/s/1fJtEUgBLBWcgbX9SJp9M5w
提取碼：kydd
為了模擬襯砌和周圍土體的接觸，可在襯砌周邊和土體界面上設置接觸對。按照施工順序，通常在Geostatic分析步移除襯砌單元，然后在第一步開挖中移除土體，在第二步中激活襯砌單元和接觸對。但是由于土體在第一步開挖后向隧道內變形（相應節點跑到了襯砌內部），激活襯砌及接觸對后，ABAQUS會認為接觸對兩側有overclosure，會在激活步中進行調整，使得隧道和土體產生向上的變形，與實際情況不符。圖1是一個平面算例開挖后的豎向位移，圖2是激活襯砌及接觸對后的豎向位移。盡管激活步中沒有施加任何荷載，土體還是出現了明顯的上抬變形。
圖1
圖2
要解決這一問題，必須讓激活后的襯砌的節點位置與開挖后土體變形相協調。這里給出一個簡單的思路：在襯砌位置處設置兩套共節點的重疊單元，單元集1采用較小的模量，在分析中一直保留；單元集2為真實的襯砌單元，在開挖后激活。
展開 
隧道開挖模擬-3
鏈接：https://pan.baidu.com/s/1xkpfwRcxiDgopkTnrTDaTA
提取碼：dpg1
盾構的掘進有兩種方式，一是利用生死單元，即先將掘進方向的開挖土體單元移除，隨后在相應位置激活盾構機單元，賦予盾構機材料。另一種方法基于大變形，精細化模擬盾構機的前進。前者較為簡單實用，應用較多。
本系列首先介紹第一種做法，算例的幾何形狀和參數參考Lambrughi A, Rodríguez L M, Castellanza R. Development andvalidation of a 3D numerical model for TBM–EPB mechanised excavations[J].Computers and Geotechnics, 2012, 40: 97-113.
主要施工過程：
（1） 初始假設盾構機處于地基之中，建立初始應力場。
（2） 移除盾構前土體，每次開挖長度1.5m。單元移除后在開挖面施加壓力215kPa。
（3） 盾構向前移動1.5m，在盾構前激活盾構機單元，盾構尾移除相應長度單元。同時激活盾尾灌漿層單元和襯砌單元。
（4） 依次向前開挖
關注點：
（1） 為了模擬盾構的掘進，需要在重復的位置建立幾套共節點的單元，可利用以下語句實現，詳見inp文件：
*ELCOPY, OLD SET=r1, NEWSET=sr1,SHIFT NODES=0, ELEMENT SHIFT=20000；該語句將在r1單元處重疊一套單元，節點不變，單元號+20000。
（2） 本例中沒有考慮灌漿壓力，可以在灌漿層單元外表面施加灌漿壓力，方法類似于開挖面上施加壓力的做法。
展開 abaqus雙線盾構隧道開挖數值模擬案例 ￥168
該案例為多層土體的雙線盾構隧道開挖，考慮了掌子面推進力和注漿壓力，模型為完整模型，不存在跑不通，有ODB結果，購買后支持售后講解，包括如何實現注漿硬化階段，地應力平衡的意義等。
附件包含雙線盾構隧道開挖的數值模擬模型以及運算結果ODB，購買后可聯系博主進行答疑。
該模型為博主自己手搓，杜絕不明不白的教程，模型上傳的目的是為了幫助大家，避免收到各類低質教程的迫害從而浪費時間。
ansys隧道開挖
/clear
/units,si
/filnam,suidao
/title,jb
/prep7
et,1,beam3
mp,ex,1,21e9
r,1,0.28,1.83e-3,0.28
et,2,plane42
keyopt,2,3,2
mp,ex,2,1.3e9
mp,prxy,2,0.4
tb,dp,2
tbdata,1,2e5,27,
mp,dens,2,2200
k,1,0,0
k,2,-6.79,-2.29
k,3,-4.87,-4.30
k,4,4.87,-4.30
k,5,6.79,-2.29
k,6,7.16,0
k,7,0,7.16
k,8,-7.16,0
k,9,-71.6,-71.6
k,10,71.6,-71.6
k,11,71.6,71.6
k,12,-71.6,71.6
larc,7,8,1,7.16
larc,8,2,1,7.16
larc,2,3,1,3.27
larc,3,4,1,16.28
larc,4,5,1,3.27
larc,5,6,1,7.16
larc,6,7,1,7.16
a,8,1,6,7
a,8,2,3,4,5,6,1
a,6,7,8,12,11
a,6,11,10,4,5
a,4,10,9,3
a,3,9,12,8,2
/replot
/pnum,line,1
/pnum,area,1
lcomb,2,3
lcomb,5,6
lcomb,7,1
lcomb,8,9
numcmb,line
/replot
latt,1,1,1
lesize,1,,,20
lesize,2,,,20
lesize,3,,,20
lesize,4,,,20
lmesh,1,4,1
aatt,2,,2
lesize,6,,,20,5
lesize,7,,,20
lesize
展開 隧道開挖過程的CEL三維大變形有限元模擬
隨著隧道的開挖，刀盤的推理和扭矩逐漸增大，并最終達到一個穩定值，最終的的推力約為25 MN，扭矩2.5MN?m。計算得到的推力和扭矩，可以為實際工程的施工提供參考，有效地控制開挖速度，避免開挖造成顯著的地表沉降。
圖6：刀盤的推力隨開挖時間的變化過程
圖7：刀盤的扭矩隨開挖時間的變化過程
5.結論
（1）CEL大變形有限元方法可以模擬隧道開挖的整個過程，能有效地模擬開挖過程中土體的大變形；
（2）CEL大變形方法可以計算盾構機開挖過程中刀盤的推力和扭矩，能為實際工程提供參考，有效地控制盾構機的開挖速度，避免隧道發生坍塌。
5.硬件與計算時常分析
（1）計算機配置：整個CEL大變形分析是在DELL工作站上進行，配置為3.7 GHz 主頻Intel Xeon W-2255 CPU，128G memory。
（2）計算時間：CEL模型中有964712三維八節點歐拉單元，整個計算時間大概12h.
展開 abaqus模擬二維隧道開挖對輸電塔的影響 ￥60
abaqus模擬二維隧道開挖對輸電塔的影響，模型文件
Dyna求解的工程爆破模擬-隧道爆破開挖
10.添加炸藥材料及空氣材料模型、狀態方程、及ALE屬性
參見我的CSDN博客<Dyna求解的工程爆破模擬-鋼筋混凝土結構抗爆>“9.給添加炸藥材料及介質材料模型、狀態方程、及ALE屬性”。
11.添加多物質流固耦合
參見我的CSDN博客<Dyna求解的工程爆破模擬-鋼筋混凝土結構抗爆>“10.多物質流固耦合添加”。其中MASTER組件包括空氣及炸藥組件，Sla-ve組件包括巖石及隧道組件。
12.添加無反射邊界條件
參見我的CSDN博客<Dyna求解的工程爆破模擬-鋼筋混凝土結構抗爆>“11.添加無反射邊界條件”。
13.導出模型
參見我的CSDN博客<Dyna求解的工程爆破模擬-鋼筋混凝土結構抗爆>“12.導出模型”。
14.空氣及炸藥匹配材料模型、狀態方程、及ALE屬性
參見我的CSDN博客<Dyna求解的工程爆破模擬-鋼筋混凝土結構抗爆>“13.介質及炸藥匹配材料模型、狀態方程、及ALE屬性”。
15.提交計算
參見我的CSDN博客<Dyna求解的工程爆破模擬-鋼筋混凝土結構抗爆>“14.提交計算”。
展開 新手模擬隧道開挖的常用幾種方法
1.作為新手，在模擬隧道方面，最開始會受《ABAQUS在巖土工程中的應用》一書上的例子的影響，在模擬時會采取隧道圍巖和襯砌分別建立PART，然后分別定義屬性，在ASSENBLY中組裝，用TIE約束圍巖和襯砌變形。采用model change 控制圍巖開挖和襯砌加入。但由于TIE約束本身的缺陷，該方法劃分網格后，節點調整后，對于大的3D模型，求解和精度都有一些問題。
2、隨著認知的增加，漸漸你會發現，其實一般隧道的模擬不是這么做的，一般是在只建立一個part，partition出襯砌和開挖塊。對于單獨襯砌建立集合和小集合后。劃分網格，在job模塊，write input，然后在inp文件中用elcopy對襯砌集合生成另一個節點號相同但是單元號不同的集合。對該集合賦予襯砌性質。最后再采用生死單元model change。
3、采用field和temperature來變化材料性質。這個方法很明了，就是只能采用的本構模型比較簡單。實際模擬不好用。inp文件下次再傳。
4、采用input方法，直接導入模型節點，和各材料的單元。可以借鑒《ABAQUS在隧道與地下工程中的應用》一書。其實該書還是不錯的。雖然很多童鞋覺得不好、但是abaqus雖然作為CAE軟件，input文件的修改和書寫還是很關鍵的能力。還是要像ANSYS一樣啃命令流呀。查閱下自帶的幫助文件。
展開 
ABAQUS環形隧道開挖的施工模擬Step by Step ￥3
ABAQUS環形隧道開挖的施工模擬-01-10.pdf
abaqus模擬隧道開挖，土體采用修正劍橋本構 ￥50
使用abaqus模擬隧道分階段開挖，土體開挖一層襯砌激活一層，土體本構采用修正劍橋本構，模型部分設置和結果見圖片
基于LS-DYNA的臨近隧道爆破開挖模擬建模分析 附LS-DYNA自定義本構子程序基本流程下載
LS-DYNA臨近隧道爆破開挖模擬建模分析
作者：turtle（在讀碩士）
擅長：hypermesh/dyna
開挖隧道初襯和既有隧道襯砌震速云圖：
開挖隧道初襯和既有隧道襯砌應力云圖：
數值模型的建立
臨近隧道爆破開挖模擬分析，采用流固偶合方法。巖石襯砌為solid單元，空氣炸藥為solid_ALE單元。炸藥和空氣采用ALE算法，并實現流固耦合的動態分析。巖石和襯砌共節點連接，空氣和炸藥共節點連接，空氣炸藥與巖石襯砌做流固偶合。采用cm-g-us單位制，爆破時間為0.2S。
幾何模型及網格劃分：
2.結果分析
2.1 開挖隧道初襯某測點震速時程曲線
2.2 不同時刻襯砌應力云圖
下載地址：LS-DYNA自定義本構子程序基本流程
展開 ansys19.0隧道爆破模擬
帶佬們，我想問一下掃略劃分畫爆破炸藥網格時，出現了形狀限制錯誤，但又顯示劃分出來了，是什么問題啊，該怎么解決！不解決會咋樣呢！劃分整體網格時又沒出現這個
