有限元法筆記(共11頁)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上文 獻(xiàn) 筆 記1、 地鐵隧道盾構(gòu)法施工過程中地層變位的三維有限元模擬1.解決的問題：預(yù)測(cè)盾構(gòu)穿越所引起的地層位移，確保已有建(構(gòu))筑物正常使用和盾構(gòu)的順利掘進(jìn)。2.基本思路：在前人工作的基礎(chǔ)上，提出一種盾構(gòu)施工三維有限元模擬方法，綜合考慮盾構(gòu)推進(jìn)、開挖面前方及隧道壁面土體的剪切擾動(dòng)、盾構(gòu)刀盤超挖空隙和盾尾空隙的閉合、盾尾注漿等因素，分析了盾構(gòu)推進(jìn)過程中隧道周圍及地表處土體的位移和變形以及橫斷面不同深度上的沉降分布規(guī)律。3.盾構(gòu)施工引起土層移動(dòng)的主要因素：a.開挖面土體開挖應(yīng)力釋放 b.刀盤超挖引起的間隙 c.盾構(gòu)與周圍土體的錯(cuò)動(dòng) d.盾構(gòu)前行產(chǎn)生盾尾空隙 e.盾尾注漿

2、 f.襯砌變形4.三維有限元模擬盾構(gòu)施工的步驟：三維有限元模擬盾構(gòu)施工的計(jì)算流程圖如圖2所示。具體的模擬步驟如下： 第 1 步：計(jì)算地層初始應(yīng)力。第 2 步：盾構(gòu)機(jī)每向前推進(jìn)1 環(huán)管片作如下變化：(1) 單元材料變化將開挖面向前推進(jìn)1 環(huán)管片的距離，同時(shí)變換單元材料，單元材料變化情況如圖3 所示，圖中虛線表示被挖掉的土單元。 (2) 荷載的模擬 位于新開挖面前方的一層土單元應(yīng)力由靜止土壓力狀態(tài)k0z，k0z，zT 變?yōu)楦飨虻葔籂顟B(tài)s，s，sT，其中k0 為靜止土壓力系數(shù)，s 為土壓艙壓力。根據(jù)式DFe = òòòBTDdxdydz計(jì)算由應(yīng)力變化產(chǎn)生的結(jié)點(diǎn)荷載，并將

3、其作用在單元結(jié)點(diǎn)上。 使開挖面上隧道外圍的一層土單元結(jié)點(diǎn)上發(fā)生盾構(gòu)刀盤超挖位移，盾尾襯砌外圍土單元結(jié)點(diǎn)發(fā)生盾尾建筑空隙位移。當(dāng)土體接觸到盾構(gòu)殼或者襯砌外壁時(shí)，再將盾構(gòu)殼和襯砌單元激活，使它們與土體共同作用。 在漿液?jiǎn)卧鈬Y(jié)點(diǎn)上加遠(yuǎn)離盾構(gòu)中心方向的注漿壓力。 在盾構(gòu)殼后部單元結(jié)點(diǎn)上加指向盾構(gòu)推進(jìn)方向的結(jié)點(diǎn)力，同時(shí)在其后部的管片單元結(jié)點(diǎn)上加反向的結(jié)點(diǎn)力，前后結(jié)點(diǎn)力的總和相等，且等于千斤頂推力的總和。 在盾尾后 1 環(huán)管片結(jié)點(diǎn)上加豎直向上的結(jié)點(diǎn)力，總和等于1 環(huán)管片內(nèi)部挖去土體的重量，模擬開挖土體引起的豎向卸荷。重復(fù)第 2 步直到計(jì)算結(jié)束。5.有限元分析模型及其參數(shù)：某盾構(gòu)隧道采用土壓平衡式盾構(gòu)。

4、盾構(gòu)機(jī)長(zhǎng)8.0 m，盾構(gòu)外徑6.34 m，盾構(gòu)殼厚7 cm，盾構(gòu)刀盤外徑比盾構(gòu)殼外徑大1 cm；管片寬1 m，厚35 cm。千斤頂總推力14 000 kN。地基土層自上而下依次為：雜填土(03.0 m)、粘質(zhì)粉土(3.09.0 m)、淤泥質(zhì)粘土(9.014.0 m)、粉質(zhì)粘土(14.028.0 m)、砂質(zhì)粘土(28.038.0 m)、粉砂(38.0 m以下)。隧道中心線埋深20 m(位于粉質(zhì)粘土層中)。有限元計(jì)算網(wǎng)格如圖4 所示，分析區(qū)域豎向深50 m，寬60 m，沿隧道縱向長(zhǎng)150 m。共5 040個(gè)單元，6 105個(gè)節(jié)點(diǎn)。 計(jì)算中土體采用鄧肯張E-u 模型，盾構(gòu)殼和管片采用線彈性材料，各層

5、土的基本物理力學(xué)指標(biāo)及鄧肯張E-u 模型參數(shù)見表 1，2。 6.有限元分析結(jié)果及所用到的曲線：6.1周圍土體的位移：a.橫斷面位移矢量圖 b.縱斷面位移矢量圖6.2地表位移：a.隧道橫斷面地面位移曲線 b.隧道縱斷面計(jì)算與實(shí)測(cè)地面沉降曲線6.3不同深度處的沉降分布:a.隧道橫斷面不同深度處的沉降曲線2、 盾構(gòu)掘進(jìn)過程土體變形特性數(shù)值模擬1. 解決的問題：運(yùn)用三維有限差分軟件FLAC3D 對(duì)盾構(gòu)掘進(jìn)工程進(jìn)行了數(shù)值仿真。2.基本思路：仿真過程中考慮了盾構(gòu)機(jī)、注漿壓力、土倉(cāng)壓力等因素，得到了地表沉降槽和土體縱向變化規(guī)律，并與實(shí)測(cè)結(jié)果進(jìn)行了對(duì)比。3.盾構(gòu)掘進(jìn)模擬過程:1、計(jì)算土體的自重應(yīng)力，得到原始應(yīng)

6、力場(chǎng)，把所有節(jié)點(diǎn)位移賦值為零；2、用model null 命令挖出y 方向30 m 左隧道內(nèi)的土體、襯砌和注漿部位的土體；3、在y=2430 m 內(nèi)的襯砌和注漿部位，換算為盾構(gòu)殼的密度、彈性模量和泊松比；4、在y=018 m 內(nèi)注漿體假定為凝狀態(tài)；5、在y=1824 m 內(nèi)注漿體假定為未凝固狀態(tài)；6、在y=024 m 內(nèi)的襯砌組部分設(shè)置未混凝土襯砌材料；7、在y=1824 m 內(nèi)注漿體半徑3.1 m、3.2 m的圓筒面上設(shè)置注漿壓力為0.2 MPa；8、在土體開挖面上即y=30 m 的隧道范圍內(nèi)設(shè)置壓力為0.2 MPa，模擬土倉(cāng)壓力。4.FLAC3D 建模:（1）材料參數(shù)：所有材料參數(shù)見表1。

7、 5.FLAC3D 數(shù)值模擬模型:參考有關(guān)有限元文獻(xiàn)的分析結(jié)果，本模型計(jì)算范圍取x 方向84 m，y 方向72 m，z 方向49.4 m，計(jì)算結(jié)果表明，該模型范圍滿足計(jì)算要求。計(jì)算模型共劃分7 284 個(gè)單元，8 112 個(gè)節(jié)點(diǎn)。如圖2，3所示。模型的邊界條件為：垂直于 x=0、x=84、y=0、y=72 的面上的法線方向被約束，底面z=0 上為固定約束。 6. 結(jié)果分析曲線：a.地表沉降 b.沿隧道軸線地表沉降 c.計(jì)算的地表沉降槽與實(shí)測(cè)比較 d.覆土厚度 8 m 時(shí)的地表隆起位移3、 盾構(gòu)千斤頂推力變化對(duì)地面變形的影響1.解決的問題：研究了盾構(gòu)推進(jìn)力的變化對(duì)地面變形的影響。2.基本思路：采

8、用三維彈塑性有限元, 在真實(shí)模擬施工情況的基礎(chǔ)上,分析由于盾構(gòu)曲線推進(jìn)、抬頭或叩頭推進(jìn)對(duì)地面變形的影響。3.分析方法:對(duì)于粘土來說, 由于粘土的特性導(dǎo)致沉陷是與時(shí)間有關(guān)的, 大多數(shù)情況下, 這種固結(jié)沉降可以采用二維模型模擬。在粘土中隧道的修建期間, 現(xiàn)場(chǎng)觀察也表明總沉陷的顯著部分, 發(fā)生在較短的期間內(nèi), 即基本上在不排水條件下發(fā)生的 。對(duì)于修建在粘土中的隧道, 一旦在初始條件(不排水)下保持穩(wěn)定, 隨著孔隙水壓力的消散, 土體強(qiáng)度增大,隧道會(huì)一直保持穩(wěn)定。鑒于上海土層的實(shí)際情況, 及本文的目的在于研究隧道開挖引起土層三維移動(dòng)情況, 所以本研究中僅考慮不排水情況。本文對(duì)SAP 程序進(jìn)行了修改,

9、使之增加了適用于土體計(jì)算的三維彈塑性功能。為了計(jì)算中的精確, 本文要用三維二十節(jié)點(diǎn)空間等參元分析整個(gè)模型, 本構(gòu)模型取為橫觀各向異性彈塑性體, 采用德魯克普拉格屈服準(zhǔn)則。開挖過程是用代表開挖一步(推進(jìn)工作面)所挖土體的單元的剛度條件以“活化減退”來模擬。這是通過將它的剛度降低到極小值來達(dá)到的。通過這些降低模量的單元,模擬出工作面的逐步推進(jìn)。在施工中, 挖土階段所挖土體采用“活化減退”的方法來模擬, 隨著降低模量單元的增多, 模擬出盾構(gòu)工作面的推進(jìn)。在盾構(gòu)工作面上施加壓力P 以模擬支護(hù)壓力。在盾構(gòu)機(jī)單元的周圍施加力F , 模擬盾構(gòu)機(jī)的工作狀態(tài)。隧道管片的施作是在工作面向前推進(jìn)了一段特定的距離或在

10、開挖后經(jīng)過一段特定的時(shí)間以后, 采用“重新活化”法模擬。在漿液?jiǎn)卧闹車┘恿藟毫 , 來模擬注漿壓力。4.計(jì)算模型及土層狀況：本文是對(duì)一個(gè)實(shí)際工程的研究, 其計(jì)算模型見圖1 , 模型中, 有限元網(wǎng)格的劃分, 是考慮了許多不同網(wǎng)格劃分的計(jì)算結(jié)果后確定的。網(wǎng)格側(cè)向邊界的選擇以它不能顯著影響沉陷預(yù)測(cè)結(jié)果為基準(zhǔn)。具體范圍如下:在隧道工作面后取9D 的長(zhǎng)度,D 為盾構(gòu)推進(jìn)的直徑, 這個(gè)距離遠(yuǎn)遠(yuǎn)超過了由文獻(xiàn) 得出的1.5 2D 的影響范圍。隧道軸線左右取9D 范圍, 自地面下的計(jì)算深度為5H ,H 為地表到拱頂?shù)木嚯x。隧道采用 4350 大刀盤土壓平衡式盾構(gòu)施工, 隧道外徑為 4.2m , 內(nèi)徑 3.

11、6m, 隧道頂覆土深度為5 6m 。隧道所處地層的地質(zhì)資料如下表所示。 5.計(jì)算結(jié)果及分析曲線: a.正常推進(jìn)隧道中心線縱截面上地表面處的隆沉情況 b.曲線推進(jìn)隧道中心線縱截面上地表面處的隆沉情況 c.抬頭推進(jìn)隧道中心線縱截面上地表面處隆沉情況 d.叩頭推進(jìn)時(shí)隧道中心線縱截面上地表面處隆沉情況 e.正常推進(jìn)隧道不同橫斷面上地表面處的變形情況 f.曲線推進(jìn)隧道不同橫斷面上地表面處變形情況 g.抬頭推進(jìn)隧道不同橫斷面上地表面處的變形情況 h.叩頭推進(jìn)時(shí)隧道不同橫斷面上地表面處的變形情況4、 盾構(gòu)隧道施工過程的有限元分析1.解決的問題：針對(duì)軟土地區(qū)的特點(diǎn), 利用大型通用有限元軟件ABAQUS 對(duì)盾構(gòu)

12、隧道施工過程進(jìn)行有限元模擬。2.基本思路：用有限元法對(duì)施工過程進(jìn)行模擬, 考慮土與襯砌結(jié)構(gòu)的共同作用, 同時(shí)還考慮了土體材料的非線性及施工參數(shù)等復(fù)雜因素的影響, 有限元分析時(shí)對(duì)各種影響因素進(jìn)行合理的簡(jiǎn)化。3.施工過程的模擬：首先就開挖前的土體施加重力和初始有效應(yīng)力場(chǎng), 同時(shí)給洞室周邊土體單元結(jié)點(diǎn)施加自重作用下的平衡反力(由相同模型施加重力和初始有效應(yīng)力場(chǎng)單獨(dú)計(jì)算獲得), 以達(dá)到地壓平衡, 作為計(jì)算的起始狀態(tài);然后再把平衡反力逐級(jí)釋放, 以模擬土體的開挖過程;在平衡反力釋放的同時(shí), 利用單元生死技術(shù), 生成襯砌單元與接觸面單元, 形成土體與襯砌共同作用。各級(jí)荷載釋放系數(shù)是通過取不同的釋放系數(shù)計(jì)算

13、比較而得出的, 假定隧道開挖初期, 洞室的荷載釋放系數(shù)為0.7,即管片施作后隧道周圍土體的荷載釋放70%, 考慮襯砌結(jié)構(gòu)接頭的影響, 管片剛度折減系數(shù)為0.75。4.工程概況與參數(shù)：a.工程概況:某地鐵由東行線和西行線組成, 東行線隧道長(zhǎng)970.15m , 西行線隧道長(zhǎng)974.481m , 隧道埋深16 30m , 鋼筋混凝土襯砌管片寬1.2m 、厚280mm , 外徑5.3m。隧道施工采用兩臺(tái)土壓平衡盾構(gòu)機(jī)推進(jìn), 同時(shí)采用同步壁后注漿技術(shù)。 b.地層參數(shù):計(jì)算范圍內(nèi)的土體由砂性土與粘性土交替組成, 土層共有7 層, 在量測(cè)斷面上, 盾構(gòu)隧道位于As2 砂性土中, 各層土體參數(shù)如表1。 c.襯

14、砌結(jié)構(gòu)參數(shù):對(duì)襯砌結(jié)構(gòu), 把襯砌作為一個(gè)整體來考慮, 并對(duì)襯砌的剛度進(jìn)行一定程度的折減以考慮其接頭的影響。設(shè)計(jì)時(shí)提供的參數(shù)為:E =3.5 ×107kPa , A =0 .28m2。 5.有限元模型:土體與襯砌的相互作用采用接觸面單元模擬, 土體采用平面應(yīng)變非協(xié)調(diào)單元, 襯砌用連續(xù)體平面應(yīng)變單元模擬。計(jì)算范圍的選取是以不顯著影響計(jì)算結(jié)果為基準(zhǔn), 具體范圍如下:盾構(gòu)隧道左右各取30m , 地表往下取60m , 對(duì)左右邊界設(shè)為水平約束, 下部則認(rèn)為只有豎向約束;為考慮土體的非線性變形, 土體材料采用Drucker -Prager 屈服準(zhǔn)則。6. 結(jié)果分析：a.實(shí)測(cè)土壓力與計(jì)算結(jié)構(gòu)對(duì)比 b

15、.實(shí)測(cè)襯砌彎矩與計(jì)算結(jié)果對(duì)比 c.實(shí)測(cè)襯砌軸力與計(jì)算結(jié)果對(duì)比5、 盾構(gòu)近距離穿越已建隧道的施工影響分析1.解決的問題：施工中對(duì)既有建（構(gòu)）筑物和自身近距離施工時(shí)的影響問題。2.基本思路：針對(duì)過黃浦江行人觀光隧道從上部穿越剛剛建成的上海地鐵2 號(hào)線越江區(qū)間隧道，建立了三維有限元計(jì)算模型，研究了由于盾構(gòu)推進(jìn)而引起的地層擾動(dòng)變形的規(guī)律性，對(duì)已建隧道產(chǎn)生的施工影響進(jìn)行了分析。3.工程概述：上海地鐵 2 號(hào)線橫貫長(zhǎng)寧、靜安、黃浦及浦東新區(qū)，除浦東東方路以南大都為農(nóng)田外，其余各段所穿越地面交通繁忙、建筑物密集以及地下管線錯(cuò)綜復(fù)雜的市政環(huán)境。其中圓形區(qū)間隧道部分西起中山公園站，東至龍東路站，雙線（上、下行）

16、全長(zhǎng)約為12 km。地鐵2號(hào)線圓形區(qū)間隧道部分采用土壓平衡式盾構(gòu)掘進(jìn)施工，盾構(gòu)外徑為6.34 m，長(zhǎng)度為6.54 m；隧道襯砌環(huán)厚1 m，襯砌混凝土強(qiáng)度等級(jí)C50，抗?jié)B等級(jí)為 0.8 MPa ，外層襯砌混凝土密度r 2500 kg/m3，彈模E34.5 GPa，泊松比m 0.2；內(nèi)層襯砌混凝土密度r 2 400 kg/m3，彈模E 28.0 GPa，泊松比m 0.2。外灘觀光隧道工程?hào)|起東方明珠電視塔西側(cè)的浦東出入口豎井，西至浦西外灘觀光平臺(tái)邊上的陳毅廣場(chǎng)，全長(zhǎng)646.70 m，隧道的縱向坡度最大達(dá)48 ，為國(guó)內(nèi)盾構(gòu)隧道之最。隧道結(jié)構(gòu)采用鋼筋混凝土單層襯砌結(jié)構(gòu)，襯砌厚度為1.2 m，隧道外徑為

17、7.48 m，內(nèi)徑6.76 m，采用f 7.65 m 鉸接式土壓平衡盾構(gòu)施工。地鐵 2 號(hào)線陸家嘴河南路區(qū)間隧道在穿越黃浦江江底的過程中，在浦西防汛墻底下與同期建設(shè)的外灘行人觀光隧道成51°21斜交，上、下行線隧道頂部與其凈距分別為1.57 m 和2.18 m，形成3 條隧道疊交穿越情況，相互影響的范圍約有3040 m，且2 條隧道還要在外灘防汛墻僅為11 m寬的預(yù)留孔中穿過，如圖1，2 所示。為減小兩項(xiàng)工程施工產(chǎn)生的相互影響，行人隧道盾構(gòu)后于地鐵隧道盾構(gòu)施工，施工先后順序?yàn)榈罔F二號(hào)線上行線、地鐵二號(hào)線下行線、外灘觀光隧道。但是，地鐵二號(hào)線隧道與外灘觀光隧道施工時(shí)間間隔僅三個(gè)月左右，

18、此時(shí)隧道尚處于非穩(wěn)定狀態(tài)；由于土體的不穩(wěn)定，必將產(chǎn)生相互影響，并且由于行人隧道的施工必將引起地鐵盾構(gòu)隧道上部覆土卸荷作用，如果施工不當(dāng)，很容易引起已建隧道的交叉段發(fā)生過大的隆起或沉降，從而導(dǎo)致地鐵隧道的縱向變形及管片接頭的滲漏水。4.盾構(gòu)掘進(jìn)過程和隧道襯砌的模擬:在隧道開挖之前，土層在其自重及其它外部荷載作用下已處于受力變形狀態(tài)，在隧道開挖之后，由于洞室周邊的應(yīng)力被解除，從而，必然要引起土層的應(yīng)力重分布，此時(shí)，襯砌結(jié)構(gòu)也參與應(yīng)力重分布的過程并由此引起襯砌結(jié)構(gòu)的內(nèi)力。故此，需要適當(dāng)模擬結(jié)構(gòu)的建造過程，本文采用如下的方法模擬土壓平衡盾構(gòu)的掘進(jìn)過程和管片的安裝過程：（1）施加土體自重，但關(guān)閉管片單元

19、，使之不參與計(jì)算，這樣將得到洞室開挖前土層的初始應(yīng)力場(chǎng)；（2）關(guān)閉開挖隧道部分的單元，同時(shí)激活襯砌單元，以模擬隧道的開挖和管片的安裝；（3）同時(shí)在開挖面上施加面力，用以模擬土壓平衡狀態(tài)，面力大小為：P =KoH，式中 g 為土層的加權(quán)平均重度（kN/m3）；H 為隧道中心處埋深（m）； Ko 為水平壓力系數(shù)。5.工程地質(zhì)：地鐵2號(hào)線圓形區(qū)間隧道，工程沿途場(chǎng)地較為平緩，地面標(biāo)高約為+2.0 m左右，表層為褐黃色粘土，沿線主要穿越的土層有：2灰色砂質(zhì)粉土層，飽和、松散稍密、夾薄層粉砂、振動(dòng)后易液化，易發(fā)生流砂，屬中壓縮性土；灰色淤泥質(zhì)粉質(zhì)粘土層，飽和、流塑軟塑，屬高壓縮性土；灰色淤泥質(zhì)粘土層，飽和

20、、流塑、軟塑、夾少量薄層粉砂，屬高壓縮性土；1-1灰色粘土層，很濕、軟塑可塑、受擾動(dòng)后沉降大，屬高偏中壓縮性土；1-2灰色粉質(zhì)粘土層，很濕、軟塑、含水量高、孔隙比大、靈敏度高，受擾動(dòng)后沉降大、局部夾薄層粉砂，屬中壓縮性土。盾構(gòu)主要穿越灰色淤泥質(zhì)粉質(zhì)粘土層及灰色粘土層。該兩層土的強(qiáng)度低，含水量高，有較大的流變性。隧道所穿越的主要土層的地基土的力學(xué)指標(biāo)如下表1。 6. 有限元模型的建立：計(jì)算模型尺寸為 70 m×60 m×40 m，即沿地鐵2 號(hào)線方向分別取70 m，沿行人觀光隧道方向取60 m，沿深度方向取40 m。計(jì)算只考慮初襯管片，而未計(jì)及內(nèi)襯的影響。土體材料按照理想彈塑

21、性介質(zhì)來考慮，選取8 結(jié)點(diǎn)實(shí)體單元來模擬；由于管片材料的剛度較大，一般認(rèn)為在彈性范圍內(nèi)工作，故選取彈性的殼單元模擬。有限元網(wǎng)格劃分示意圖如圖3 所示。 模型的邊界條件，除了上表面為自有面外，其余4 個(gè)側(cè)面和底部均施加法向約束。盾構(gòu)每推進(jìn)步長(zhǎng)即每環(huán)襯砌步步長(zhǎng)，為了減少計(jì)算量，本文采用了簡(jiǎn)化的推進(jìn)步長(zhǎng)，即每步推進(jìn)6 m。計(jì)算共分10個(gè)荷載步，分別驗(yàn)算行人隧道盾構(gòu)在推進(jìn)到地鐵2號(hào)線之前、到達(dá)時(shí)以及到達(dá)后土層和2 號(hào)線隧道的變形和應(yīng)力變化情況，如圖4 所示。7. 結(jié)果分析曲線：a.行人隧道地表面橫向沉降曲線 b.2 號(hào)線上行線隧道底部地面隆起曲線 c.卸荷引起隧道“反彈變形”圖 d.沿行人隧道軸向地表面沉降曲線 e.不同載荷步時(shí)沿行人隧道軸向地表面沉降曲線6、 交疊隧道盾構(gòu)法施工土體變形的三維數(shù)值模擬1.解決的問題：以正在建設(shè)中的上海市軌道交通明珠線二期工程上、下行線近距離交疊區(qū)間隧道盾構(gòu)施工為研究對(duì)象, 按“ 先下后上” 的盾構(gòu)推進(jìn)過程, 采用三維彈塑性有限元ANSYS 程序軟件模擬了交疊隧道土層位移以及地表沉降曲面在盾構(gòu)推進(jìn)中的發(fā)展變化, 得出了地面最大沉降量在上行隧道推進(jìn)后將出現(xiàn)大幅度增長(zhǎng), 且在推進(jìn)前期沉降增幅最為顯著的結(jié)論.給出了隧道周圍土體塑性區(qū)的大小與分布, 論證了要將地面沉降控制在