含雙承壓水地層中地連墻槽壁穩(wěn)定性有限元分析

周忠群，張孟喜，王 鑫，張 強(qiáng)，薛青松

(1.上海大學(xué)力學(xué)與工程科學(xué)學(xué)院，上海 200444；2.中鐵二十局集團(tuán)第一工程有限公司，江蘇 蘇州 215151)

隨著城市地下空間不斷開發(fā)，地下連續(xù)墻技術(shù)也得到廣泛運(yùn)用，地下連續(xù)墻成槽過程中槽壁穩(wěn)定性控制一直是施工質(zhì)量控制的一項(xiàng)重點(diǎn)工作，而地連墻槽壁穩(wěn)定性分析一直是巖土工程中的一個(gè)難點(diǎn)問題，正確地分析槽壁穩(wěn)定性，對(duì)于確保工程質(zhì)量安全和工程成本控制具有十分重要的意義。影響槽壁穩(wěn)定性的因素很多，而對(duì)于穿湖地層的地連墻成槽施工，地下水的滲流作用無疑是影響其穩(wěn)定性的主要因素之一。

當(dāng)前，國(guó)內(nèi)外學(xué)者對(duì)滲流作用下的地連墻槽壁穩(wěn)定性問題作了大量的研究工作。黃鋒[1]考慮成槽開挖時(shí)的滲流作用，基于極限分析上限定理，構(gòu)建運(yùn)動(dòng)速度機(jī)動(dòng)場(chǎng)對(duì)地連墻局部和整體失穩(wěn)模式進(jìn)行了分析；劉楊等[2]針對(duì)富水軟弱地層地連墻進(jìn)行研究，考慮了承壓水作用，基于極限平衡方法分析了槽壁局部失穩(wěn)機(jī)制；Liu等[3]考慮了滲流作用，采用滲流解析模型方法對(duì)槽壁進(jìn)行局部穩(wěn)定性分析；李建高等[4]考慮了水的滲流作用對(duì)槽壁穩(wěn)定性進(jìn)行分析，采用解析方法和二維數(shù)值計(jì)算得出槽壁的安全系數(shù)；孟德平[5]考慮了地下水位的影響，采用有限元強(qiáng)度折減法對(duì)槽壁穩(wěn)定性進(jìn)行分析。以上學(xué)者雖然考慮了地下水的滲流作用，但均未對(duì)雙承壓含水層槽壁失穩(wěn)模式進(jìn)行研究，且缺乏對(duì)水的滲流場(chǎng)分析。

考慮滲流作用下地連墻槽壁穩(wěn)定性分析是一個(gè)比較復(fù)雜的問題，目前，槽壁穩(wěn)定性分析方法主要有極限平衡法、有限元法和有限差分法。其中，有限元強(qiáng)度折減法不需要假定破壞面的形狀和位置，也能定量反映出地連墻槽壁的漸進(jìn)失穩(wěn)演化過程。因此，本文依托蘇州春申湖路快速化改造工程為背景，采用ABAQUS軟件基于強(qiáng)度折減法對(duì)穿湖地層雙承壓含水層地連墻成槽施工進(jìn)行槽壁失穩(wěn)破壞模式分析，有限元強(qiáng)度折減法通過不斷調(diào)整折減系數(shù)對(duì)槽壁穩(wěn)定性進(jìn)行分析，從而確定滑裂面位置和滑動(dòng)體形狀，揭示雙承壓含水層槽壁失穩(wěn)破壞模式及破壞過程，并結(jié)合現(xiàn)場(chǎng)聲波實(shí)測(cè)進(jìn)行對(duì)比分析。

1 槽壁破壞特征及強(qiáng)度折減分析

1.1 滲流作用下地連墻失穩(wěn)破壞過程

目前，學(xué)者們普遍認(rèn)為邊坡失穩(wěn)破壞具有漸進(jìn)破壞特征[6-8]，根據(jù)以往大量工程實(shí)際案例以及學(xué)者們的研究發(fā)現(xiàn)，地連墻槽壁失穩(wěn)破壞和邊坡失穩(wěn)破壞存在類似特征，槽壁失穩(wěn)往往表現(xiàn)為局部塌孔向整體滑動(dòng)演化的漸進(jìn)破壞過程，并非一次瞬間完成。由于地連墻槽壁土層本身存在局部軟弱土層，且在高水位區(qū)進(jìn)行成槽開挖，滲透系數(shù)較大的粉砂層對(duì)地下水的滲流力作用影響較為明顯，往往容易最先出現(xiàn)局部應(yīng)力集中現(xiàn)象，從而使這部分土體單元最先出現(xiàn)局部塑性破壞區(qū)。隨著局部塑性破壞區(qū)不斷擴(kuò)大，局部塑性區(qū)的應(yīng)力集中會(huì)出現(xiàn)轉(zhuǎn)移和調(diào)整，新的應(yīng)力集中帶會(huì)形成，之后應(yīng)力集中帶會(huì)不斷發(fā)展擴(kuò)大，最終會(huì)在槽壁土體形成貫通的塑性破壞滑裂帶。

邊坡的滑裂帶具有軟化特性[9]，地連墻槽壁也具有類似性質(zhì)，在達(dá)到抗剪強(qiáng)度峰值會(huì)降到殘余值，含承壓水的雙粉砂層槽壁漸進(jìn)破壞過程見圖1，塑性滑裂帶起于含承壓水的薄弱粉砂層底部，即點(diǎn)a處，在點(diǎn)a和b，土體剪應(yīng)力超出峰值點(diǎn)并降到殘余值，點(diǎn)c的剪應(yīng)力剛達(dá)到峰值，而點(diǎn)d的剪應(yīng)力遠(yuǎn)低于峰值抗剪強(qiáng)度。隨著滑裂帶的漸進(jìn)擴(kuò)展，雙粉砂層之間的點(diǎn)d也將出現(xiàn)應(yīng)力集中并達(dá)到峰值點(diǎn)，滑裂帶由潛在滑裂面漸進(jìn)擴(kuò)展并形成最終貫通滑動(dòng)面?？梢姡诨罨褞У男纬蛇^程中，滑帶的強(qiáng)度從峰值降低到殘余值，這是地連墻槽壁整體滑動(dòng)漸進(jìn)形成過程中滑裂帶強(qiáng)度的演化特征。

圖1 含承壓水的雙粉砂層槽壁漸進(jìn)破壞示意

1.2 考慮滲流的強(qiáng)度折減法的基本原理

穿湖地層地連墻成槽施工，圍堰外側(cè)水位較高，產(chǎn)生的滲流力對(duì)槽壁作用較為明顯，因此，有限元模擬中必須要考慮滲流場(chǎng)作用。而實(shí)際中，土體滲流場(chǎng)和應(yīng)力場(chǎng)之間是相互聯(lián)系、相互影響的，其實(shí)質(zhì)為應(yīng)力場(chǎng)作用于土體上，使土體變形，從而改變土體的滲流特性，滲流場(chǎng)因此發(fā)生變化，滲流場(chǎng)的變化又會(huì)導(dǎo)致水頭分布的改變，使土體的有效應(yīng)力和作用在土體的滲流力發(fā)生改變，即應(yīng)力場(chǎng)發(fā)生改變[10]。因此，為了分析穩(wěn)定滲流下槽壁穩(wěn)定性，必須考慮滲流和應(yīng)力的耦合。ABAQUS軟件可以對(duì)土中水的滲流和土體應(yīng)力進(jìn)行耦合分析，孔隙水的模擬主要由滲透系數(shù)和孔隙比來設(shè)置，而靜水壓力作用通過設(shè)置浸潤(rùn)面和底面的孔隙壓力從而自動(dòng)生成線性分布孔隙水壓力，兩層承壓水層應(yīng)進(jìn)行分段設(shè)置孔隙水壓力。

強(qiáng)度折減法在ABAQUS中的實(shí)現(xiàn)主要通過對(duì)土層塑性參數(shù)屬性進(jìn)行場(chǎng)變量設(shè)置，通過場(chǎng)變量中多組折減系數(shù)的設(shè)置來進(jìn)行有限元強(qiáng)度折減，具體折減過程是將地連墻所在土層抗剪強(qiáng)度參數(shù)c，tanφ同時(shí)除以一個(gè)試算的折減系數(shù)Fr，得到一組新的cm、φm值，然后將其作為新的參數(shù)代入有限元進(jìn)行試算，反復(fù)多次，直至槽壁土體符合給定的臨界破壞狀態(tài)判別標(biāo)準(zhǔn)時(shí)，對(duì)應(yīng)的Fr即為槽壁的穩(wěn)定性安全系數(shù)。其中，參數(shù)cm、φm分別由以下兩式求得，即

cm=c/Fr

(1)

φm=arctan(tanφ/Fr)

(2)

1.3 槽壁失穩(wěn)的判斷準(zhǔn)則

強(qiáng)度折減法的判斷準(zhǔn)則最早是在判斷邊坡失穩(wěn)中誕生的，之后在邊坡和基坑中得到廣泛應(yīng)用[11-13]，就方法而言，判斷準(zhǔn)則在地連墻槽壁失穩(wěn)中也同樣適用，目前用于失穩(wěn)判據(jù)大致可以分3類：①收斂判據(jù)。折減后的土體強(qiáng)度參數(shù)使得有限元計(jì)算在規(guī)定的迭代次數(shù)內(nèi)不能收斂。②位移判據(jù)。當(dāng)折減系數(shù)增大到某一特定值時(shí)，槽壁位移會(huì)發(fā)生突變，即可認(rèn)為槽壁已經(jīng)失穩(wěn)破壞。③塑性區(qū)判據(jù)。通過對(duì)有可能發(fā)生破壞的區(qū)域內(nèi)廣義剪應(yīng)變等物理量的變化和分布來判斷，當(dāng)塑性區(qū)從槽壁深處某一位置貫通至地面。目前，位移判據(jù)和塑性區(qū)判據(jù)應(yīng)用最廣，位移判據(jù)可以清晰得出槽壁的穩(wěn)定安全系數(shù)，而塑性區(qū)判據(jù)可以清晰看出滑裂面的位置和形狀，能很好反映出漸進(jìn)破壞演化過程。為了全面準(zhǔn)確地分析槽壁破壞過程以及得出安全系數(shù)，本文依據(jù)位移突變和塑性區(qū)貫通判斷準(zhǔn)則對(duì)槽壁穩(wěn)定性進(jìn)行綜合判斷。

2 有限元模型計(jì)算

2.1 工程背景

為了研究地連墻穿越雙粉砂層在承壓水作用下的成槽破壞模式及破壞過程，本文依托蘇州春申湖路快速化改造工程進(jìn)行有限元模型計(jì)算。該工程以明挖隧道形式由西向東穿越陽澄西湖，湖中段基坑施工長(zhǎng)度為2 685 m，基坑深度10.54～18.17 m，隧道采用圍堰明挖法，圍堰采用拉森鋼板樁進(jìn)行隔水施工，抽水變陸地后在湖底施工，現(xiàn)場(chǎng)施工和地質(zhì)分布見圖2。湖中段基坑圍護(hù)結(jié)構(gòu)采用地下連續(xù)墻，由于在湖中進(jìn)行地連墻成槽施工，地下水較為豐富，潛水位位于地下2 m，且在④2粉土夾砂層和⑥3粉土夾砂層分別存在微承壓水和承壓水，微承壓水頭標(biāo)高為1 m，承壓水頭標(biāo)高為-2 m，地連墻在進(jìn)行成槽開挖時(shí)很容易發(fā)生塌槽現(xiàn)象，給工程質(zhì)量安全和成本控制帶來極大挑戰(zhàn)。

圖2 春申湖明挖隧道現(xiàn)場(chǎng)情況

2.2 數(shù)值模型

本文采用有限元軟件ABAQUS對(duì)春申湖隧道地連墻成槽施工進(jìn)行數(shù)值模擬，基于對(duì)稱原理，可以將現(xiàn)場(chǎng)地連墻進(jìn)行簡(jiǎn)化，按現(xiàn)場(chǎng)實(shí)際1/4尺寸建立三維數(shù)值模型，考慮到模型的邊界尺寸效應(yīng)，三維模型尺寸取60 m×60 m×90 m，依據(jù)據(jù)現(xiàn)場(chǎng)實(shí)際地連墻長(zhǎng)度、厚度和深度，模型中分別取3、0.6、45 m。有限元網(wǎng)格劃分如圖3所示，模型采用孔隙流體/應(yīng)力耦合的C3D8P八節(jié)點(diǎn)六面體單元，土體采用Mohr-Coulomb屈服準(zhǔn)則模型，導(dǎo)墻采用彈性單元，模型側(cè)面邊界約束水平位移，底面邊界約束水平和豎向兩個(gè)方向的位移，上表面為自由面，孔壓邊界條件為側(cè)面和底面默認(rèn)為不透水面，水位浸潤(rùn)面孔壓設(shè)為零，底面孔壓根據(jù)浸潤(rùn)面以下土層深度呈現(xiàn)線性分布，但由于粉砂層中存在承壓水，需對(duì)潛水層和承壓水層進(jìn)行分段設(shè)置孔壓，從而更準(zhǔn)確地模擬現(xiàn)場(chǎng)實(shí)際的孔隙水壓力分布。

圖3 有限元網(wǎng)格劃分(單位：m)

根據(jù)實(shí)際地質(zhì)勘測(cè)，沿線場(chǎng)地地表下90.3 m深度范圍內(nèi)地基土構(gòu)成除填土外，其余為第四系濱海、第四系河泛、河床相沉積物，一般由黏性土、粉(砂)土組成。計(jì)算模型的土層物理力學(xué)參數(shù)信息見表1。

表1 土層物理力學(xué)參數(shù)

3 數(shù)值結(jié)果分析

3.1 有限元滲流場(chǎng)分析

成槽開挖前沿深度方向的孔隙水壓力分布見圖4a，由于該地層中在上粉砂層和下粉砂層中存在微承壓水和承壓水，且微承壓水頭為1 m，而承壓水頭為-2 m，和潛水面齊平，故圖4a中微承壓水層孔隙水壓力比鄰近層凸出一部分，說明該層水壓力比鄰近層的孔隙水壓大，而承壓水層無突變。地連墻成槽開挖后的孔壓等值線云圖見圖4b。由圖4b可知，此時(shí)槽壁滲流的浸潤(rùn)面即為孔壓為0的分界面，開挖后地連墻附近孔壓分布呈現(xiàn)漏斗狀。

圖4 沿深度孔隙水壓力分布

土體開挖過程中孔隙水有效速度矢量見圖5a，由圖5a中可知，孔隙水流主要集中在槽壁底端附近，說明滲流力在槽壁底端出現(xiàn)集中現(xiàn)象，在下粉砂層出現(xiàn)流線聚集，且流線較長(zhǎng)，流線由槽壁外側(cè)向槽壁內(nèi)側(cè)移動(dòng)，由此說明該土層流速過快，產(chǎn)生的滲流力較大，對(duì)槽壁的沖擊較大，而該土層滲流力較大是由于粉砂層滲透系數(shù)較大，且靠近槽段底端。下粉砂層底端處孔隙水有效滲透速度隨折減系數(shù)變化見圖5b。由圖5b可知，下粉砂層有效滲透速度隨著折減系數(shù)的增加而增加，增長(zhǎng)曲線呈現(xiàn)線性增長(zhǎng)，當(dāng)折減系數(shù)依次增大，土層的抗剪強(qiáng)度隨之降低，土層塑性破壞程度越來越大，導(dǎo)致土層中孔隙占比增大，滲透系數(shù)也隨之增大，孔隙水的有效滲透速度將相應(yīng)增大。由此說明，應(yīng)力場(chǎng)和滲流場(chǎng)是相互影響的，強(qiáng)度折減不僅會(huì)降低土層的抗剪強(qiáng)度，也在一定程度上改變土層的滲透系數(shù)，造成滲流力作用更加明顯，從而加速粉砂層局部破壞進(jìn)程。

圖5 孔隙水有效速度分布及變化

3.2 槽壁失穩(wěn)漸進(jìn)破壞模式分析

圖6為雙粉砂層地連墻槽壁在承壓水作用下塑性區(qū)破壞演變過程圖，從圖6可以看出，有限元強(qiáng)度折減法能夠真實(shí)地反映槽壁整體失穩(wěn)滑動(dòng)的破壞過程，通過對(duì)土層進(jìn)行整體強(qiáng)度折減，塑性破壞區(qū)一開始出現(xiàn)于槽壁下粉砂層底部和上粉砂層，隨著強(qiáng)度折減次數(shù)的不斷迭加，塑性區(qū)不斷由上下粉砂層向中間擴(kuò)展增大，直到塑性破壞區(qū)重疊連通，由此槽壁出現(xiàn)大范圍整體失穩(wěn)滑動(dòng)破壞，滑動(dòng)破壞面最終呈現(xiàn)弓字狀。塑性破壞區(qū)起于雙粉砂層底部，由于粉砂層本身不太穩(wěn)定，且滲透系數(shù)較大，在承壓水作用下很容易出現(xiàn)塑性破壞，從圖5a也可以清晰看出，當(dāng)成槽開挖至槽底時(shí)，孔隙流體會(huì)出現(xiàn)流動(dòng)現(xiàn)象，從而造成孔隙水壓力重新分布，會(huì)在槽壁底部產(chǎn)生較大的滲流力，造成下粉砂層底部最先出現(xiàn)塑性區(qū)，而上粉砂層在微承壓水作用下也會(huì)很快出現(xiàn)塑性區(qū)，由于雙粉砂層起始小范圍的塑性區(qū)會(huì)造成鄰近土層抗剪強(qiáng)度不斷被削弱，從而使塑性區(qū)不斷向鄰近土層擴(kuò)展，塑性區(qū)不斷延伸重疊，最終會(huì)連通雙粉砂層，造成槽壁整體滑動(dòng)失穩(wěn)。其實(shí)，塑性破壞區(qū)一般起于薄弱土層，這里的薄弱土層包括兩個(gè)方面，一方面土層自身抗剪強(qiáng)度參數(shù)較小，土層自身存在安全隱患；另一方面，土層所受外力過大，如滲流作用下的滲流力，會(huì)對(duì)土層產(chǎn)生很大的沖擊作用，從而造成土層的不穩(wěn)定性。由以上分析得出，地連墻槽壁穩(wěn)定狀態(tài)會(huì)隨著時(shí)間的推移不斷發(fā)生變化，槽壁失穩(wěn)破壞會(huì)經(jīng)歷從局部失穩(wěn)到整體失穩(wěn)的破壞演化過程，穩(wěn)定狀態(tài)會(huì)經(jīng)歷穩(wěn)定、欠穩(wěn)定和不穩(wěn)定三個(gè)過程。圖7為槽壁下粉砂層底端水平位移隨折減系數(shù)變化。由圖7可知，在折減系數(shù)達(dá)到1.35時(shí)，位移發(fā)生突變，此時(shí)槽壁處于不穩(wěn)定狀態(tài)，往往容易發(fā)生局部失穩(wěn)破壞現(xiàn)象，故槽壁穩(wěn)定性安全系數(shù)為1.35。

圖6 塑性破壞區(qū)發(fā)展過程

圖7 水平位移隨折減系數(shù)變化

3.3 雙粉砂層開挖相互耦合影響分析

不同開挖深度下粉砂層等效塑性應(yīng)變?cè)茍D見圖8。由圖8可知，開挖15 m時(shí)，經(jīng)過不斷的強(qiáng)度折減，槽壁會(huì)在粉砂層最先出現(xiàn)局部塑性破壞區(qū)；開挖45 m時(shí)，在經(jīng)過多次重復(fù)強(qiáng)度折減下，槽壁會(huì)在雙粉砂層出現(xiàn)明顯的塑性破壞區(qū)。兩種工況下上粉砂層的等效塑性應(yīng)變對(duì)比見表2。由表2可知，兩種工況下粉砂層塑性應(yīng)變趨勢(shì)基本一致，最大的塑性應(yīng)變都是出現(xiàn)在粉砂層底部，且隨著下粉砂層的開挖會(huì)對(duì)上粉砂層的塑性應(yīng)變產(chǎn)生一定影響，粉砂層頂部塑性應(yīng)變由0.12變?yōu)?.27，變化幅度較小，而粉砂層底部由0.67變?yōu)?.65，變化幅度較大，由此說明，下粉砂層開挖會(huì)對(duì)上粉砂層產(chǎn)生一定的影響，且對(duì)粉砂層底部塑性應(yīng)變比頂部影響大。綜上所述，粉砂層相比于其他黏性土層對(duì)土體強(qiáng)度折減更為敏感，在土體強(qiáng)度不斷折減降低的條件下，更容易發(fā)生槽壁失穩(wěn)破壞。

圖8 不同開挖深度下粉砂層等效塑性應(yīng)變?cè)茍D

表2 不同開挖深度下上粉砂層等效塑性應(yīng)變

3.4 開挖深度對(duì)槽壁安全系數(shù)的影響

考慮滲流和不考慮滲流作用時(shí)槽壁安全系數(shù)隨開挖深度變化見圖9。由圖9可知，兩種曲線下，考慮滲流作用比不考慮滲流作用下的曲線整體偏下，說明滲流作用對(duì)成槽開挖中槽壁穩(wěn)定性影響較大，考慮滲流作用時(shí)，實(shí)際工程中對(duì)高水位區(qū)的成槽開挖一定要嚴(yán)格控制地下水位造成的滲流影響，尤其在對(duì)承壓含水層地連墻成槽開挖，一定要采用減壓降水措施。此外，隨著成槽開挖深度的不斷增加，槽壁穩(wěn)定性安全系數(shù)會(huì)下降很快，當(dāng)開挖深度為15 m時(shí)，穩(wěn)定性安全系數(shù)為1.61，開挖深度為45 m時(shí)，穩(wěn)定性安全系數(shù)下降至1.35，降幅達(dá)16.1%。由此可見，開挖深度對(duì)地連墻成槽開挖影響較大，實(shí)際工程中對(duì)于超深地連墻的開挖一定要嚴(yán)格控制成槽施工對(duì)槽壁土體造成的擾動(dòng)影響。

圖9 槽壁安全系數(shù)隨開挖深度變化

3.5 現(xiàn)場(chǎng)聲波實(shí)測(cè)對(duì)比

為了驗(yàn)證數(shù)值分析結(jié)果的準(zhǔn)確性，本文選取對(duì)應(yīng)地層抓槽試驗(yàn)和正式成槽開挖的聲波實(shí)測(cè)圖和數(shù)值模擬結(jié)果進(jìn)行對(duì)比分析，由于抓槽試驗(yàn)開挖了15 m，故數(shù)值結(jié)果選取開挖15 m時(shí)槽壁等效塑性等值云圖作為對(duì)比參照。數(shù)值模擬與聲波實(shí)測(cè)對(duì)比見圖10。

圖10 數(shù)值模擬與聲波實(shí)測(cè)對(duì)比

從圖10a可以看出，成槽開挖15 m數(shù)值模擬云圖中塑性應(yīng)變等值線主要集中分布于④2粉土夾砂層，顯示出塑性破壞區(qū)主要出現(xiàn)在粉土夾砂層中，聲波實(shí)測(cè)圖表現(xiàn)為在④2粉土夾砂層中發(fā)生局部塌孔現(xiàn)象，這是由于粉土夾砂層本身不穩(wěn)定，滲透系數(shù)較大，在微承壓水的作用下容易發(fā)生滲透破壞，從而造成局部塌孔。

從圖10b可以看出，不考慮微承壓水的條件下成槽開挖45 m時(shí)，數(shù)值結(jié)果并未出現(xiàn)圖8b中兩層明顯塑性破壞區(qū)，而正式成槽開挖聲波實(shí)測(cè)也只在下層粉土夾砂層出現(xiàn)大面積塌槽，而上層粉土夾砂層并未出現(xiàn)明顯塌孔，這是由于在此前進(jìn)行了抓槽試驗(yàn)，得出微承壓水是造成粉砂層塌孔的主要原因，正式施工前采用減壓降水措施，降低了微承壓水頭，從而提高了④2粉土夾砂層的成槽穩(wěn)定性。實(shí)際上，強(qiáng)度折減下的槽壁塑性破壞區(qū)范圍往往比現(xiàn)場(chǎng)聲波實(shí)測(cè)顯示的粉砂層局部塌孔的要大很多，因?yàn)楝F(xiàn)場(chǎng)聲波實(shí)測(cè)只是整體強(qiáng)度折減下出現(xiàn)整體滑動(dòng)的前期局部破壞某一階段，并不是強(qiáng)度折減下計(jì)算不收斂的最終結(jié)果，在這一局部階段，數(shù)值分析結(jié)果表現(xiàn)為槽壁局部小范圍塑性區(qū)破壞，聲波實(shí)測(cè)結(jié)果表現(xiàn)為粉砂層局部塌孔或整層失穩(wěn)塌槽。經(jīng)實(shí)測(cè)對(duì)比分析所得，數(shù)值模擬結(jié)果和聲波實(shí)測(cè)基本吻合，既驗(yàn)證了數(shù)值模擬的準(zhǔn)確性，也揭示出粉砂層在承壓水作用下較其他黏性土層更容易發(fā)生局部塑性失穩(wěn)的事實(shí)。

4 結(jié) 論

本文開展?jié)B流作用下雙承壓含水層地連墻成槽開挖的數(shù)值模型分析，通過強(qiáng)度折減法對(duì)地連墻破壞模式和槽壁穩(wěn)定性進(jìn)行研究，主要得到以下結(jié)論：

(1)考慮滲流作用的粉砂地層地連墻成槽開挖，滲透系數(shù)較大的粉砂層滲透速度較大，所受滲透力較大，容易發(fā)生滲透破壞，且粉砂層有效滲透速度會(huì)隨著土體強(qiáng)度折減系數(shù)的增加呈現(xiàn)線性增加趨勢(shì)，說明土體的應(yīng)力場(chǎng)與滲流場(chǎng)存在相互影響，粉砂層局部失穩(wěn)破壞往往是由于土體強(qiáng)度降低和滲透力增大綜合影響產(chǎn)生的。

(2)確定了滲流作用下雙承壓含水層在成槽開挖中的破壞模式，在雙粉砂層底端最先出現(xiàn)塑性區(qū)，隨后塑性區(qū)由雙粉砂層不斷向中間擴(kuò)展，最后雙粉砂層會(huì)連通形成貫通的塑性滑裂帶，滑裂面形狀最終呈現(xiàn)弓字狀，槽壁整體滑動(dòng)失穩(wěn)往往是由于粉砂層局部失穩(wěn)造成的。

(3)地連墻成槽開挖對(duì)雙粉砂層等效塑性應(yīng)變存在相互影響，隨著下粉砂層的開挖會(huì)導(dǎo)致上粉砂層等效塑性應(yīng)變?cè)龃?，且?duì)粉砂層的影響表現(xiàn)為上部影響小，底部影響最大。此外，考慮滲流作用時(shí)槽壁穩(wěn)定性系數(shù)整體會(huì)比不考慮滲流作用時(shí)大，隨著開挖深度的不斷增大，槽壁穩(wěn)定性安全系數(shù)會(huì)不斷降低。

(4)結(jié)合實(shí)際工程的聲波實(shí)測(cè)結(jié)果進(jìn)行對(duì)比分析，數(shù)值計(jì)算結(jié)果和聲波實(shí)測(cè)較為吻合，從而一定程度上驗(yàn)證了數(shù)值模擬的準(zhǔn)確性，也更好地揭示了粉砂層在承壓水作用下較其他黏性土層更容易發(fā)生局部塑性失穩(wěn)的事實(shí)。