周連偉

（大連市市政設(shè)計(jì)研究院有限責(zé)任公司,遼寧 大連 116011）

樁錨復(fù)合支護(hù)深基坑變形的流固耦合分析

周連偉

（大連市市政設(shè)計(jì)研究院有限責(zé)任公司,遼寧 大連 116011）

深基坑開挖降水過程中,坑內(nèi)外壓力差可能引起嚴(yán)重的工程事故。通過有限元軟件MIDAS/GTS,采用滲流應(yīng)力耦合理論及摩爾庫(kù)倫模型,結(jié)合某樁錨復(fù)合支護(hù)深基坑工程實(shí)例建立三維有限元模型分析了該深基坑工程的變形情況,并與不考慮耦合的基坑變形情況進(jìn)行對(duì)比,主要包括錨桿軸力圖、樁剪力圖與彎矩圖、基坑地表沉降、坑底回彈、側(cè)向位移等。結(jié)果表明:在兩種不同的有限元模擬條件下基坑變形形態(tài)基本一致；總的來說,流固耦合分析引起的基坑變形小于不考慮耦合分析的基坑變形,支護(hù)結(jié)構(gòu)的受力變形也較小,更接近實(shí)際情況?？紤]地下水流固耦合分析對(duì)基坑變形的影響不容忽視,對(duì)實(shí)際工程的設(shè)計(jì)優(yōu)化有一定的參考意義。

樁錨復(fù)合支護(hù)；滲流；流固耦合；變形

0 引言

樁錨支護(hù)體系是將受拉桿件的一端固定在開挖基坑的穩(wěn)定地層中,另一端與圍護(hù)樁相聯(lián)的基坑支護(hù)體系,它是在巖石錨桿理論研究比較成熟的基礎(chǔ)上發(fā)展起來的一種擋土結(jié)構(gòu),安全經(jīng)濟(jì)的特點(diǎn)使它廣泛應(yīng)用于邊坡和深基坑支護(hù)工程中[1]。隨著地下空間的大力開發(fā),基坑工程向著大、深、廣的方向發(fā)展,工程實(shí)踐往往遇到地下水問題,很多工程事故都是對(duì)地下水考慮不周引起的。目前工程中往往采取基坑降水的方法,來控制地下水影響。單純的應(yīng)用土壓力原理進(jìn)行水土合算或水土分算來分析,忽略了土體中水的流動(dòng)性[2],往往不考慮孔隙水滲流場(chǎng)與應(yīng)力場(chǎng)相互作用的影響,為深基坑的穩(wěn)定帶來很多安全隱患。

要準(zhǔn)確考慮基坑降水引起的滲流與應(yīng)力的變化,必須進(jìn)行流固耦合分析。本文以大連某樁錨復(fù)合支護(hù)深基坑工程為背景,運(yùn)用MIDAS/GTS軟件,采用數(shù)值方法分析深基坑在滲流應(yīng)力耦合情況下的受力狀態(tài)和變形,并與不考慮耦合的深基坑變形作對(duì)比,得出預(yù)期結(jié)論。

1 流固耦合計(jì)算原理

流固耦合力學(xué)的重要特征是兩相介質(zhì)之間的相互作用,變形固體在流體載荷作用下會(huì)產(chǎn)生變形或運(yùn)動(dòng)。變形或運(yùn)動(dòng)又反過來影響流體運(yùn)動(dòng),從而改變流體載荷的分布和大小,正是這種相互作用將在不同條件下產(chǎn)生形形色色的流固耦合現(xiàn)象[8]。

本文的流固耦合現(xiàn)象為基坑降水開挖過程中應(yīng)力場(chǎng)與滲流場(chǎng)的相互作用?；油馏w開挖簡(jiǎn)化為卸荷條件下卸除側(cè)向約束,同時(shí)存在地下水滲流固結(jié)的變形過程,土體中地下水滲流與固結(jié)變形都對(duì)基坑圍護(hù)墻上土壓力產(chǎn)生影響,利用巖土工程有限元軟件建立分析模型,設(shè)立特定水力邊界條件,通過數(shù)值方法模擬“流固耦合”效應(yīng)來考慮基坑開挖的變形。

2 工程概況

大連某深基坑工程,基坑深10 m,長(zhǎng)40 m,寬20 m,地下水位位于地表-3 m左右。工程場(chǎng)區(qū)土層自上而下大致分為回填土、風(fēng)化土、風(fēng)化巖,土體參數(shù)及選用的材料模型見表1。由于地下水較豐富,圍護(hù)墻采用高15 m,直徑為1 m,凈間距為1 m的C30混凝土鉆孔灌注樁和1 m寬的止水帷幕。共設(shè)2道錨桿,分別設(shè)置在樁頂以下3 m和7 m處,長(zhǎng)度自上而下依次為14.5 m、10.0 m,傾角分別為33.7°、45°?；硬扇∵吔邓_挖邊支護(hù)的分步施工方案,第一步開挖3 m,設(shè)置第一排錨桿；第二步開挖4 m,設(shè)置第二排錨桿；第三步開挖3 m。第一步至第三步每步開挖之前,均考慮將地下水降到開挖面以下1 m。因MIDAS/GTS計(jì)算單元有限,三維有限元模型水平向（X方向）取80 m,縱向（Y方向）取8 m,豎向（Z方向）取30 m,進(jìn)行數(shù)值模擬,二維幾何模型及計(jì)算模型見圖1、圖2。巖土體參數(shù)及支護(hù)材料模型見表1。

表1 巖土體參數(shù)及支護(hù)材料模型

圖1 二維幾何模型（單位:mm）

圖2 計(jì)算模型

3 計(jì)算模型

土工程通常采用的土體彈塑性本構(gòu)模型有M-C模型和D-P模型,二者都考慮了靜水壓力對(duì)屈服強(qiáng)度的影響,本文采用M-C模型模擬土體。

本工程實(shí)例有下面幾個(gè)基本假定:

（1）假設(shè)含水層是水平的,并且水流服從達(dá)西定律[10]。

（2）假定降水一定時(shí)間后基坑中心處水位達(dá)到要求深度并處于暫時(shí)穩(wěn)定狀態(tài)。

（3）忽略降水時(shí)非飽和土對(duì)滲透效果的影響,模型滲透系數(shù)均采用土層飽和狀態(tài)下的滲透系數(shù)。

開挖前打樁引起的土體原位應(yīng)力狀態(tài)的改變不予考慮,土體的初始應(yīng)力假定為靜止土壓力[11]。

4 計(jì)算結(jié)果分析

按照實(shí)際施工順序,分3步模擬基坑開挖的施工過程,MIDAS通過激活和鈍化土體單元及增加結(jié)構(gòu)支護(hù)單元來模擬開挖和支護(hù),并通過定義各施工階段,形成完整的施工工況,再對(duì)整體工況進(jìn)行分析,得出計(jì)算結(jié)果。本文中流固耦合工況通過考慮開挖基坑的降水變形實(shí)現(xiàn)了滲流場(chǎng)和應(yīng)力場(chǎng)的耦合模擬,包括1個(gè)打樁階段、3個(gè)降水階段、3個(gè)降水變形階段、3個(gè)開挖階段共10個(gè)主要施工階段；而不考慮流固耦合工況則通過降水后直接開挖實(shí)現(xiàn),包括1個(gè)打樁階段、3個(gè)降水階段、3個(gè)開挖階段共7個(gè)主要施工階段。兩種情況的計(jì)算位移結(jié)果見圖3。

圖3 位移云圖

4.1 基坑水平位移

圖4分別給出了兩種情況下各施工階段基坑壁水平位移有限元計(jì)算結(jié)果。從圖中可以看出,兩種情況變形曲線的形狀基本一致,在各步開挖的施工階段中,基坑側(cè)壁的最大位移均在基坑的開挖面附近。在各開挖階段,考慮流固耦合作用的基坑壁變形要比不考慮流固耦合作用小,在開挖最終結(jié)束時(shí),最大位移減小了5 mm。從圖4（a）可知,完全流固耦合的滲流方式考慮了基坑內(nèi)外側(cè)水壓,模擬了降水變形施工階段,表現(xiàn)出孔隙水壓力逐漸消散,基坑壁水平位移回彈這一過程,符合實(shí)際情況。

圖4 基坑壁的水平位移

圖5分別給出了兩種情況下各施工階段基坑地表沉降有限元計(jì)算結(jié)果。對(duì)比圖5（a）、（b）可以看出,考慮流固耦合時(shí)基坑周圍的地表沉降規(guī)律跟不考慮流固耦合時(shí)基本一致,都在距離坑壁一定距離處出現(xiàn)最大沉降位移,對(duì)于本文的基坑實(shí)例,考慮流固耦合時(shí)最大沉降位移為12 mm,略小于不耦合開挖方式的13 mm。由圖5（a）知,每次降水后,基坑的地表沉降位移在坑壁一定距離處會(huì)有大幅度的增加,在距坑壁較遠(yuǎn)處會(huì)有一定的回彈,而受開挖影響的地表豎向位移反而較小,考慮流固耦合情況的基坑開挖詳細(xì)的表現(xiàn)了基坑地表沉降變形的主要原因,值得深究。

4.3 基坑坑底隆起

對(duì)比圖6（a）、（b）兩種情況下各施工階段基坑坑底豎向位移有限元計(jì)算結(jié)果可以明顯的看出,隨著開挖深度的增加,兩種情況的坑底隆起量增加,考慮完全流固耦合時(shí)的回彈隆起位移要小于不考慮完全耦合作用的隆起位移,至開挖結(jié)束時(shí),最大隆起位移相差3 mm左右。由圖知,兩種情況下坑底土體隆起量在距坑壁較近的范圍內(nèi)呈近似曲線增大,且曲率較大；在距坑壁較遠(yuǎn)處,土體隆起量雖繼續(xù)增加,但增速較緩慢。隨著基坑開挖深度的增加,這一現(xiàn)象也更明顯,陡增的坑底范圍在縮小,緩增的范圍則進(jìn)一步加大。由圖6（a）看出,流固耦合情況下,每一次降水后坑底都會(huì)大幅下沉,但下一步的開挖則會(huì)引起坑底的重新隆起,坑底豎向位移的變化過程能更清晰的表現(xiàn)出來,開挖過后的每一次降水都會(huì)引起坑底回彈量一定程度的回落。

圖5 基坑的地表沉降位移

圖6 基坑坑底隆起位移

4.4 錨桿軸力圖

由數(shù)值模擬結(jié)果知,錨桿在自由段軸力均勻分布,為錨固段軸力最大值,符合實(shí)際情況。圖7是基坑開挖過程中各排錨桿錨固段軸力分布圖。四圖中錨桿錨固段軸力在剛進(jìn)入錨固段時(shí)錨桿軸力減幅較大,在錨固段的近端則是緩慢減小,在錨固段的末端則減小為零。隨著基坑開挖深度的增加,各排錨桿的軸力增加,在最后一步開挖過程中,錨桿軸力達(dá)到最大值,第一排錨桿的軸力值明顯大于第二排錨桿的軸力值；通過兩種情況下各排錨桿軸力對(duì)比可知,在各段開挖過程中,考慮完全流固耦合情況下各排錨桿軸力與不完全流固耦合情況下差別很大,直至最后一步開挖時(shí),第一排錨桿錨固段軸力相差100 kN,第二排錨桿錨固段軸力相差120 kN,且每次降水后的一定時(shí)間內(nèi),各排錨桿的軸力都有小幅的下降過程。

圖7 各排錨桿錨固段軸力

4.5 樁的剪力

基坑開挖10 m,樁和止水帷幕組成的圍護(hù)樁的深度為15 m。圖8為兩種情況下基坑開挖過程中樁的剪力分布圖。由圖8知,在各施工階段中,支護(hù)樁的最大剪力均發(fā)生在支護(hù)樁的中部偏下1 m左右,兩種情況下剪力曲線分布的形狀基本一致,呈現(xiàn)了良好的規(guī)律性,接近于線性變化,符合實(shí)際情況。如下圖8所示,支護(hù)樁樁身剪力在樁頂位置較小,然后逐漸增大,在達(dá)到第一次的開挖深度3 m后出現(xiàn)第一個(gè)極值,后又逐漸減小達(dá)到零點(diǎn)以后又反向增大出現(xiàn)峰值后又開始減小直到第二次開挖深度7 m又出現(xiàn)極值,直至最后在樁中部左右出現(xiàn)最大剪力值,中部以下的變化規(guī)律仍大致呈現(xiàn)出樁上部剪力變化規(guī)律。對(duì)比兩種情況可知,基坑各個(gè)施工階段考慮流固耦合時(shí)樁身剪力都偏小,在最后開挖結(jié)束時(shí),剪力相差最大可達(dá)到60 kN。

圖8 支護(hù)樁剪力

4.6 樁的彎矩

圖9是基坑開挖過程中樁身彎矩分布圖。樁身彎矩的總體變化過程大致為弓形,隨著開挖深度的增加,樁的上部分樁身彎矩先增加到一個(gè)極值再減小為零后反向增加到另一個(gè)極值,再經(jīng)過一次先減小再增加到最大值的過程；而在基坑底的下半部分,樁身彎矩則大致呈現(xiàn)逐漸減小,甚至變?yōu)榱愕内厔?shì)。由剪應(yīng)力和彎矩存在一階導(dǎo)數(shù)關(guān)系知,當(dāng)剪力為零時(shí),對(duì)應(yīng)的彎矩達(dá)到最大值；由圖9知,樁身的上半部分范圍內(nèi),彎矩達(dá)到最大值時(shí),所對(duì)應(yīng)的剪力值最小或等于零,從理論上證明了本次數(shù)值模擬結(jié)論的正確性。同時(shí),對(duì)比兩圖可知,在基坑的各個(gè)施工階段,考慮流固耦合時(shí)樁身的彎矩都較小,在最后開挖結(jié)束時(shí),彎矩相差最大可達(dá)到200 kN·m。

4.7 樁的軸力

圖10是基坑開挖過程中樁身軸力分布圖。由圖10知,開挖過程中,樁主要承受壓力作用。兩圖的計(jì)算結(jié)果表明,樁的軸力分布曲線大致一致,均隨樁深度的增加出現(xiàn)先增大后減小的趨勢(shì)。隨著開挖深度的增加及施工的向前推進(jìn),樁的軸力逐漸增加,開挖深度越深,軸力增加的速度越快。支護(hù)樁各施工階段的最大軸力均出現(xiàn)在基坑各步開挖面上,即分別在-3.0 m、-7.0 m、-10.0 m位置附近。對(duì)比兩圖計(jì)算結(jié)果知,在各開挖階段,考慮流固耦合作用的樁身軸力要比不考慮流固耦合作用小,在開挖最終結(jié)束時(shí),軸力最大相差100 kN左右。

圖9 支護(hù)樁彎矩

圖10 支護(hù)樁軸力

5 結(jié)語(yǔ)

（1）在基坑各步開挖的施工過程中,基坑側(cè)壁的最大位移均在基坑的開挖面附近；基坑地表最大沉降位移在距離坑壁一定距離處；坑底土體隆起量在距坑壁較近的范圍內(nèi)呈近似曲線增大,且曲率較大；在距坑壁較遠(yuǎn)處,土體隆起量雖繼續(xù)增加,但增速較緩慢。由本文考慮的兩種情況可知,完全流固耦合的基坑開挖方式減小了上述三種變形位移量,但程度較小,影響較小。

（2）錨桿在自由段軸力均勻分布,為錨固段軸力最大值,錨固段軸力在剛進(jìn)入錨固段時(shí)錨桿軸力減幅較大,在錨固段的近端則是緩慢減小,在錨固段的末端則減小為零。對(duì)比本文兩種模擬情況可知,考慮完全流固耦合情況下各排錨桿軸力與不完全流固耦合情況下差別很大,錨桿軸力較小,在同等情況下,最大減小了120 kN,說明考慮完全流固耦合情況下,可以進(jìn)一步優(yōu)化錨桿的設(shè)計(jì)。

（3）在基坑各施工階段中,支護(hù)樁的最大剪力均發(fā)生在支護(hù)樁的中部偏下1 m左右,剪力曲線分布呈現(xiàn)了良好的規(guī)律性,接近于線性變化；樁身彎矩的總體變化過程大致為弓形,呈現(xiàn)出拋物線的變化規(guī)律,樁身的上半部分范圍內(nèi),彎矩達(dá)到最大值時(shí),所對(duì)應(yīng)的剪力值最小或等于零,從理論上證明了本次數(shù)值模擬結(jié)論的正確性；樁的軸力隨樁深度的增加出現(xiàn)先增大后減小的趨勢(shì),在各開挖面出現(xiàn)軸力最大值。對(duì)比兩種模擬情況可知,考慮完全流固耦合情況下,支護(hù)樁的剪力、彎矩、軸力都偏小,且差值很大,在同等條件下,剪力最大差值為60 kN,彎矩最大差值達(dá)到200 kN·m,軸力最大差值為100 kN左右。表明考慮完全流固耦合對(duì)支護(hù)樁的受力變形有較大的影響,值得進(jìn)一步的研究。

[1]徐勇,楊挺,王心聯(lián).樁錨支護(hù)體系在大型深基坑工程中的應(yīng)用[J].地下工間與工程學(xué)報(bào),2006,2(4):646-649.

[2]張小偉,姚笑青.基坑工程變形的滲流應(yīng)力耦合分析[J].地下工間與工程學(xué)報(bào),2012,8(2):339-344.

[3]朱彥鵬,任永忠,周勇.基于滲流場(chǎng)-應(yīng)力場(chǎng)耦合作用下的深基坑降水支護(hù)結(jié)構(gòu)的位移研究[J].工程勘察,2012(7):32-36.

[4]張飛,李鏡培.考慮流固耦合效應(yīng)的基坑水土壓力計(jì)算[J].工程勘察,2011(2):1-6.

[5]何招智.深基坑工程降水與地面沉降耦合數(shù)值模擬研究[J].中國(guó)市政工程,2012,163(6):97-99.

[6]張培印,李順群.流固耦合對(duì)基坑開挖變形性狀的影響研究[J].工程勘察,2014(5):16-20.

[7]賈彩虹,王翔,王媛.考慮滲流-應(yīng)力耦合作用的基坑變形研究[J].武漢理工大學(xué)學(xué)報(bào),2010,32(1):119-122.

[8]邢景棠,周盛,崔爾杰.流固耦合力學(xué)概述[J].力學(xué)進(jìn)展,1997,27(1): 20-39.

[9]華洪勛.考慮滲流影響與分步開挖的基坑變形性狀分析研究[D].廣東廣州:廣東工業(yè)大學(xué),2008.

[10]馮懷平,張伏光,岳祖潤(rùn).考慮流固耦合作用的深基坑有限元分析[J].地下工間與工程學(xué)報(bào),2012,8(2):286-291.

[11]楊春山,何娜,張雅寧.考慮滲流-應(yīng)力耦合基坑開挖降水?dāng)?shù)值分析[J].廣東工業(yè)大學(xué)學(xué)報(bào),2013,30(4):43-48.

TU753

B

1009-7716（2017）07-0292-05

10.16799/j.cnki.csdqyfh.2017.07.089

2017-04-06

周連偉（1977-）,男,遼寧大連人,高級(jí)工程師,從事橋梁隧道設(shè)計(jì)研究工作。