胡崢嶸，池麗敏，閆園園

（南京市水利規(guī)劃設(shè)計(jì)院股份有限公司，江蘇南京210006）

土石圍堰安全分析及其除險(xiǎn)加固措施研究

胡崢嶸，池麗敏，閆園園

（南京市水利規(guī)劃設(shè)計(jì)院股份有限公司，江蘇南京210006）

圍堰是保證涉水工程干地施工的重要建筑物。在水力等外界因素作用下，其運(yùn)行機(jī)制復(fù)雜。為有效反映滲流作用對(duì)圍堰安全的影響，以某碼頭工程的圍堰為例，對(duì)其安全性進(jìn)行數(shù)值模擬。通過(guò)對(duì)圍堰滲流場(chǎng)和應(yīng)力場(chǎng)進(jìn)行耦合分析和研究，計(jì)算其應(yīng)力、位移、安全系數(shù)等相關(guān)參數(shù)。研究表明，流固耦合作用的計(jì)算結(jié)果有利于保證圍堰的設(shè)計(jì)安全。

圍堰；滲流；流固耦合；除險(xiǎn)加固；數(shù)值模擬

0 引言

圍堰是保證涉水工程干地施工所建的臨時(shí)性擋水建筑物，其安全問(wèn)題關(guān)乎著主體工程的施工進(jìn)度和施工安全。由于圍堰失穩(wěn)大多與滲流有關(guān)，因而防滲處理值得關(guān)注。目前，工程中圍堰的設(shè)計(jì)計(jì)算都將滲流場(chǎng)和應(yīng)力場(chǎng)分開(kāi)。但是，在應(yīng)力場(chǎng)、滲流場(chǎng)等多場(chǎng)并存的復(fù)雜環(huán)境中，由于各場(chǎng)間相互影響、相互作用，致使圍堰的運(yùn)行機(jī)制復(fù)雜。因此，對(duì)圍堰進(jìn)行穩(wěn)定分析時(shí)，考慮流固耦合作用是很有必要的。羅曉輝基于基坑穩(wěn)定流與非穩(wěn)定流的數(shù)值分析，對(duì)基坑開(kāi)挖過(guò)程中滲流場(chǎng)的變化規(guī)律以及對(duì)應(yīng)力場(chǎng)產(chǎn)生的影響進(jìn)行了探討［1］。徐則民等論述了滲流場(chǎng)與應(yīng)力場(chǎng)耦合分析的基本原理及其在邊坡穩(wěn)定性評(píng)價(jià)中應(yīng)用的理論基礎(chǔ)和技術(shù)路線(xiàn)［2］。程彬、盧靖等人采用軟件Geostudio中的Seep和Sigma兩個(gè)模塊對(duì)陜西府谷縣清水川電廠邊坡進(jìn)行耦合與非耦合分析，并根據(jù)應(yīng)力場(chǎng)、邊坡穩(wěn)定性的差異，得出耦合分析下更符合實(shí)際的結(jié)果［3］。本文以某碼頭工程圍堰為例，分析計(jì)算流固耦合作用下該圍堰的應(yīng)力、變形等參數(shù)，提出了除險(xiǎn)加固措施，評(píng)估加固效果。

1 流固耦合的基本原理

1.1 本構(gòu)模型

1993年，F(xiàn)redlund和Rahardjio將非飽和土體的二維本構(gòu)模型的應(yīng)力應(yīng)變關(guān)系表示成公式（1），也可表示成增量應(yīng)力應(yīng)變形式，如公式（2）所示［4］。

式中：τ為剪應(yīng)力；ε為法向應(yīng)變；σ為法向應(yīng)力；υ為泊松比；γ為剪應(yīng)變；H為與基質(zhì)吸力（uauw）有關(guān)的非飽和土模量；E為彈性模量；D為本構(gòu)矩陣；

根據(jù)虛功原理，內(nèi)虛功等于外虛功，則存在式（3）所示的關(guān)系式。

1.2 滲流方程

對(duì)于單元土體，根據(jù)達(dá)西定律，滲流方程表示成公式（4）。

式中：kx、ky分別為x和y方向的滲透系數(shù)；υw為滲透速度；γw為水的重度；θw為單位體積含水量；t為時(shí)間。

對(duì)單元土體進(jìn)行有限元分析，可得到一個(gè)包含增量孔隙水壓力的滲流方程，如公式（5）所示［5］。

對(duì)式（2）、式（3）和式（5）聯(lián)立求解，即實(shí)現(xiàn)耦合分析。在Geostudio中，流固耦合的實(shí)質(zhì)是以比奧固結(jié)理論為基礎(chǔ)，用有限元建立的以節(jié)點(diǎn)位移和孔隙水壓力為已知數(shù)，求解應(yīng)力場(chǎng)和滲流場(chǎng)的方程組［6］。

2 工程實(shí)例分析

2.1 工程概況

某工業(yè)園區(qū)通用碼頭工程岸線(xiàn)總長(zhǎng)1 910m，在12＃泊位基坑開(kāi)挖至坑底設(shè)計(jì)標(biāo)高（＋4.6m）時(shí)，前期的地質(zhì)勘探孔處出現(xiàn)泉眼，數(shù)量較多，進(jìn)而形成管涌，圍堰發(fā)生滲透破壞，基坑被水淹沒(méi)，其計(jì)算斷面如圖1所示。

該圍堰從上到下依次分為5層，根據(jù)工程勘察報(bào)告，各土層材料參數(shù)如表1所示。

圖1 圍堰計(jì)算斷面（單位：m）Fig.1 Cofferdam calculation section（Unit:m）

表1 各土層材料參數(shù)表Tab.1 Soil layer material parameter

2.2 計(jì)算模型的邊界條件

模型計(jì)算時(shí)，左右側(cè)邊界為法向支座約束，最底部邊界采用固定約束，頂部為自由邊界。

2.3 流固耦合計(jì)算流程

流固耦合的實(shí)現(xiàn)過(guò)程為：利用Seep／w模塊進(jìn)行初始孔隙水壓力計(jì)算，并將結(jié)果導(dǎo)入到sigma模塊，計(jì)算初始應(yīng)力場(chǎng)。然后，根據(jù)初始孔隙水壓力和初始應(yīng)力，計(jì)算滲流場(chǎng)與應(yīng)力場(chǎng)耦合，并將得到的結(jié)果導(dǎo)入下一步的耦合分析中。如此重復(fù)計(jì)算，即可得到每一步的計(jì)算結(jié)果。具體計(jì)算流程如圖2所示。

2.4 滲流分析

一般認(rèn)為，基坑底面的水力坡降對(duì)坑底抗?jié)B穩(wěn)定有關(guān)鍵性影響［7］。該工程滲流計(jì)算結(jié)果如圖3和圖4所示。由圖3和圖4可知，圍堰下游坡腳處的水力梯度較大，為0.241，是基坑最易發(fā)生滲透破壞的部位，且坡腳處滲流流速較大。參考文獻(xiàn)，管涌型土的臨界水力梯度可按公式（6）計(jì)算［8］。

圖2 流固耦合計(jì)算流程圖Fig.2 Fluid-structure interaction calculation flow

式中：GS為土體比重；n為土的孔隙率；d5、d20分別為小于該粒徑的土粒含量，分別為5%和20%。

根據(jù)式（6），計(jì)算出粉沙土層的臨界水力梯度為0.139。因此，可以判斷圍堰后基坑會(huì)發(fā)生管涌，影響施工進(jìn)度及安全。

2.5 應(yīng)力分析

圖3 滲流計(jì)算結(jié)果Fig.3 Seepage calculation results

圖4 圍堰下游水力坡度Fig.4 Cofferdam downstream hydraulic slope

圖6 耦合作用下剪應(yīng)力（單位：kPa）Fig.6 Shear stress under coupling（Unit:kPa）

圖5和圖6分別為流固耦合及非耦合作用下圍堰剪應(yīng)力變化云圖。由圖5和圖6可知，耦合作用下和整體剪應(yīng)力大于非耦合作用下的整體剪應(yīng)力。這說(shuō)明，邊坡因滲流作用產(chǎn)生的最大剪應(yīng)力比靜水壓力作用產(chǎn)生的最大剪應(yīng)力要大，邊坡在滲流作用下，因剪切應(yīng)力增大，更容易產(chǎn)生剪切破壞。

2.6 變形分析

圖7和圖8分別為流固耦合及非耦合作用下圍堰剪位移變化云圖。由圖7和圖8可知，耦合作用下的整體位移大于非耦合作用下的整體位移。這說(shuō)明，結(jié)構(gòu)計(jì)算時(shí)考慮滲流作用圍堰變形比靜水壓力作用產(chǎn)生的位移要大。

圖7 無(wú)耦合下總位移（單位：m）Fig.7 The total disp lacem ent under no coup ling（Unit:m）

圖8 耦合作用下總位移（單位：m）Fig.8 The total disp lacement under coup ling（Unit:m）

2.7 穩(wěn)定分析

圖9和圖10分別為流固耦合及非耦合作用下圍堰邊坡穩(wěn)定計(jì)算結(jié)果。由上述分析可知，考慮流固耦合作用，圍堰整體應(yīng)力和變形增大，更容易產(chǎn)生破壞，因此安全系數(shù)降低。其主要原因是，當(dāng)水流入土體內(nèi)部時(shí)，滑動(dòng)力增大，使土體性質(zhì)發(fā)生改變，抗剪強(qiáng)度降低。水在這時(shí)起到了“潤(rùn)滑劑”的作用［9］，與前述分析結(jié)論一致。因此，在進(jìn)行圍堰設(shè)計(jì)計(jì)算時(shí)，若合理考慮流固耦合作用，計(jì)算結(jié)果有利于保證圍堰加固設(shè)計(jì)安全。

圖9 無(wú)耦合下圍堰邊坡穩(wěn)定分析Fig.9 The slope stability analysis of cofferdam under no coup ling

圖10 耦合作用下圍堰邊坡穩(wěn)定分析Fig.10 The slope stability analysis of cofferdam under coupling

2.8 防滲分析

滲流作用是影響圍堰安全的控制性因素，因此需對(duì)該圍堰進(jìn)行防滲處理。在水頭差的作用下，粉沙土層滲透系數(shù)較大，抗?jié)B性弱，是產(chǎn)生滲透破壞的主要因素。因此，為了阻斷粉細(xì)沙層的透水，采用水泥攪拌樁對(duì)圍堰進(jìn)行止水加固。防滲墻加固深度需穿過(guò)粉沙層，進(jìn)入黏土層。由此確定攪拌樁有效樁長(zhǎng)12m，單排布置，直徑為700mm。圖11為防滲加固后的滲流計(jì)算結(jié)果。由圖11可知，加固后浸潤(rùn)線(xiàn)降低，坡腳水力梯度明顯減小，最大為0.015，坡腳處滲流流速顯著降低。這表明防滲加固效果較好。

圖11 防滲加固后滲流計(jì)算結(jié)果Fig.11 Seepage calcu lation results after reinforcement

3 結(jié)語(yǔ)

（1）利用geostudio實(shí)現(xiàn)圍堰滲流、應(yīng)力變形、邊坡穩(wěn)定的數(shù)值模擬，可及時(shí)辨識(shí)圍堰的安全性態(tài)。

（2）考慮流固耦合作用對(duì)圍堰安全的影響十分合理，計(jì)算結(jié)果有利于保證其設(shè)計(jì)安全。

（3）對(duì)防滲加固措施進(jìn)行數(shù)值模擬，可為加固效果的評(píng)估提供科學(xué)依據(jù)。

［1］羅曉輝.深基坑開(kāi)挖滲流與應(yīng)力耦合分析［J］.工程勘察，1996（6）：37－41.

［2］徐則民，黃潤(rùn)秋.基于水一力耦合理論的超深隧道圍巖滲透性預(yù)測(cè)［J］.成都理工學(xué)院學(xué)報(bào)，2001，28（2）.

［3］程彬，盧靖.基于Geostudio的邊坡滲流場(chǎng)與應(yīng)力場(chǎng)耦合［J］.山西建筑，2010（3），146－147.

［4］萬(wàn)文.水位循環(huán)升降下庫(kù)岸邊坡滲流與應(yīng)力耦合分析及穩(wěn)定性研究［D］.南昌：南昌大學(xué)，2012.

［5］倪梅三，華躍.基于Geostudio的某土石壩穩(wěn)定性流固耦合分析［J］.礦業(yè)快報(bào)，2008（2）：23－25.

［6］John Krahn.Seepage Modeling w ith SEEP／W［S］.Canada：EO－SLOPE International，LTD，2004.

［7］陶明星，劉建民.基坑滲流的數(shù)值模擬與分析［J］.工程勘察，2006（1）：23－25.

［8］盧廷浩.土力學(xué)［M］.南京：河海大學(xué)出版社，2004.

［9］張強(qiáng)勇，李曉靜.紫坪鋪高邊坡穩(wěn)定性的流固耦合計(jì)算分析［J］.人民長(zhǎng)江，2005，36（11）：45－47.

[責(zé)任編輯 楊明慶]

Study on Earth－rock Cofferdam Safety Analysis and Reinforcement Measures

HU Zheng－rong，CHI Li－min；YAN Yuan－yuan
（Nanjing Water Planning and Designing Institute CO.Ltd，Nanjing 210006，Jiangsu，China）

Cofferdam is the important building to ensure the wading engineering performed in the dry land.Under the influence of external factors such as hydraulic power，its running mechanism is complex.In order to effectively reflect the influence of seepage effect on cofferdam safety，this paper takes the cofferdam of some dock project for example，and processes the numerical simulation on its safety.It calculates the related parameters such as stress，displacement and safety coefficient by coupling analysis and research of the cofferdam seepage field and stress field.Study shows that the calculation result of fluid－structure interaction is helpful to ensure the safety of the cofferdam design.

Cofferdam；seepage；fluid－structure interaction；reinforcement

TV551.3

A

10.13681／j.cnki.cn41－1282／tv.2016.04.004

2016-07-17

胡崢嶸(1983-)，女，湖北紅安人，工程師，碩士，主要從事水工結(jié)構(gòu)設(shè)計(jì)。