徐 建 忠

(浙江中水東方建設(shè)工程咨詢有限公司，浙江 杭州 310005)

0 引 言

近年來(lái)，水利行業(yè)隨著B(niǎo)IM技術(shù)的大力普及，水利工程的可視化仿真發(fā)展突飛猛進(jìn)，尤其是復(fù)雜地質(zhì)模型的仿真建模得到了前所未有的提高。運(yùn)用Inventor軟件建立水利工程的三維模型和建筑物的實(shí)體模型，再結(jié)合ABAQUS的實(shí)體布爾運(yùn)算功能，可以有效規(guī)避inventor實(shí)體模型導(dǎo)入ABAQUS時(shí)幾何信息錯(cuò)誤的產(chǎn)生。采用ABAQUS網(wǎng)格劃分和優(yōu)化實(shí)體模型，可以創(chuàng)建三維滲流有限元網(wǎng)格模型，進(jìn)而分析其穩(wěn)定性，同時(shí)，還可以將此項(xiàng)研究成果作為三維仿真建模與數(shù)值計(jì)算分析(CAD/CAE)集成技術(shù)的參考[1]。本文選擇了某水庫(kù)工程作為研究對(duì)象，探究三維建模在水利工程灌漿穩(wěn)定性分析中的應(yīng)用。

某水庫(kù)工程位于浙江地區(qū)，工程等級(jí)為Ⅲ級(jí)，土石壩是其主要水工建筑物壩型，其最大壩高為83.5 m，建筑等級(jí)提高Ⅰ級(jí)至Ⅱ級(jí)。大壩采用黏土心墻防滲。心墻頂寬5 m，上下游坡比1︰0.25。在心墻底部設(shè)置1.5 m深的黏結(jié)槽。將0.5 m厚C15混凝土墊層分別設(shè)置在黏結(jié)槽底以及基礎(chǔ)接觸面，設(shè)置兩層反濾層心墻于下游。大壩上下游壩坡比分別為1︰1.9和1︰1.7。

1 壩基防滲設(shè)計(jì)

河床沖洪積層和岸坡殘積層應(yīng)徹底清除；除兩個(gè)斷層的壩心墻縫槽嵌入基巖2.8～2.5 m外，其余部位的壩心墻縫槽一般嵌入基巖1.5 m。心墻使用固結(jié)灌漿為基礎(chǔ)，保持5 m的深度。保持8 m的固結(jié)灌漿深度在兩個(gè)斷層之間。大壩裝有三排帷幕。大壩心墻軸線是大壩主帷幕灌漿軸線所在。輔助帷幕位于主帷幕上游1.2 m處。主幕和副幕為懸掛式，主帷幕灌漿底界由0.8倍水頭控制。輔助帷幕布置在河床段和0.5倍水頭以下的右壩段。第三排帷幕設(shè)置在另一斷層處，位于主帷幕下游側(cè)，底界與主帷幕相同。大壩帷幕灌漿防滲標(biāo)準(zhǔn)為Q≤ 3 Lu，固結(jié)灌漿防滲標(biāo)準(zhǔn)為Q≤ 5 Lu。

2 穩(wěn)定性分析理論

根據(jù)達(dá)西定律，將壩區(qū)滲流場(chǎng)視為非均質(zhì)各向異性不可壓縮土的三維穩(wěn)定滲流場(chǎng)。三維模型穩(wěn)定性的滲流控制方程為[2]：

(1)

式中xj為坐標(biāo)，i=1，2，3；ky為一個(gè)二階對(duì)稱滲透率張量；h為總水頭；Q為滲透池中的源項(xiàng)或匯項(xiàng)。

計(jì)算中使用的滲流計(jì)算邊界見(jiàn)圖1，邊界條件理論如下：

(2)

圖1 滲流計(jì)算邊界示意圖

其中，已知函數(shù)是h1和h2；滲流邊界平面外法線的余弦是ni，i = 1， 2， 3；可知水頭的第一類滲流邊界條是Γ1= BG和CD，第二類滲流邊界條件是Γ2= GA、AF、FE和 DC ，滲流池中真實(shí)和假想滲流區(qū)域之間的滲流自由面是Γ3= GE，滲流逸出面是Γ4= ED，因?yàn)槭孪炔⒉涣私鉂B流逃逸點(diǎn)E的具體位置，以及所有逃逸面具體的尺寸，所以，這屬于邊界非線性滲流問(wèn)題，要經(jīng)過(guò)迭代求解可知，邊界法向流是QN，流出為正。

3 三維有限元建模

依據(jù)本項(xiàng)目工程的地形地貌，了解主要結(jié)構(gòu)面斷層的分布，結(jié)合建筑物設(shè)置和防滲措施，構(gòu)建出涵蓋心墻壩結(jié)構(gòu)(黏土心墻、反濾層和堆石區(qū))、灌漿帷幕、灌漿蓋板、固結(jié)灌漿等在內(nèi)的本項(xiàng)目工程的三維有限元模型，同時(shí)模擬了壩基的主要斷層和基巖，其中包括基巖、斷層、河床溶蝕帶等[3]。

首先，從項(xiàng)目的工程高線圖中找出等高線點(diǎn)數(shù)據(jù)，使用rhino軟件創(chuàng)建三維地形曲面，建立基礎(chǔ)實(shí)體，然后導(dǎo)入到Inventor 軟件中。Inventor可以創(chuàng)建壩基開(kāi)挖面、溶蝕帶實(shí)體、滲透性邊界面、斷層實(shí)體、固結(jié)灌漿實(shí)體、大壩分區(qū)實(shí)體、帷幕灌漿實(shí)體。由于模型的復(fù)雜性，如果在inventor中通過(guò)布爾運(yùn)算形成項(xiàng)目建筑實(shí)體，然后將其導(dǎo)入ABAQUS，則經(jīng)常會(huì)發(fā)生幾何錯(cuò)誤，以致于造成操作失敗。運(yùn)用超強(qiáng)的ABAQUS布爾運(yùn)算功能，導(dǎo)入由invnetor創(chuàng)建的各種實(shí)體和曲面，然后進(jìn)行布爾運(yùn)算，以此完成最后的工程建筑實(shí)體。ABAQUS 能夠準(zhǔn)確識(shí)別實(shí)體并將它們以平滑網(wǎng)格劃分。

該項(xiàng)目工程的三維模型共有422 340個(gè)單元，74 123個(gè)節(jié)點(diǎn)，主要組成部分為四面體單元。其生成的三維有限元網(wǎng)格模型，主要分為大壩模型、斷層和河床溶蝕帶模型、固結(jié)灌漿、灌漿蓋板、四個(gè)斷層模型以及防滲帷幕模型[4-5]。

4 材料滲透性、邊界條件及計(jì)算工況

4.1 材料滲透特性參數(shù)

依據(jù)主要建筑材料室內(nèi)試驗(yàn)/地層、主要結(jié)構(gòu)面綜合、鉆孔水壓試驗(yàn)結(jié)果，該項(xiàng)目工程材料分為多個(gè)材料區(qū)，材料滲透特性參數(shù)見(jiàn)表1。

表1 材料滲透特性參數(shù)

4.2 邊界條件

將計(jì)算區(qū)域周圍的截止邊界和底邊界條件視為水頭界面；在地表邊界壩軸線上游方向，水頭邊界要低于水庫(kù)上游1 551.92 m的位置；已知水頭邊界條件位于壩軸線下游側(cè)，也低于下游水位的位置，可能的滲漏逸水面要高于下游水位的位置。

4.3 計(jì)算條件

計(jì)算工況(表2)綜合了河床段基巖局部帷幕損壞對(duì)庫(kù)區(qū)水流造成的影響，標(biāo)明了此次工程三維穩(wěn)定性有限元計(jì)算的分析工況。

表2 計(jì)算工況

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

5.1 帷幕完整條件下的工況

通過(guò)ABAQUS的非飽和滲流分析模塊進(jìn)行滲流計(jì)算。從該項(xiàng)目滲流特征結(jié)合區(qū)域整體滲透場(chǎng)圖來(lái)看，具備明顯的三維滲流特征[6]：蓄水之后，庫(kù)水通過(guò)壩基滲透帶向下游滲漏，同時(shí)繞過(guò)左右壩肩經(jīng)過(guò)兩岸滲透帶。因此可得出：(1)壩體防滲心墻以及壩基帷幕灌漿的效果非常明顯，大壩滲流場(chǎng)水頭勢(shì)等值線分布具備高密度，并在心墻內(nèi)大量減小總水頭差；(2)從兩岸滲流場(chǎng)分布看，在兩岸延伸帷幕作用下，左右岸地下水位顯著降低，尤其是兩岸壩肩段。由于地質(zhì)條件很差，沿河床和溶蝕帶上下游均有斷層穿過(guò)，對(duì)庫(kù)區(qū)滲流產(chǎn)生不利影響。經(jīng)過(guò)處理，壩體浸潤(rùn)線在心墻內(nèi)迅速下降，水位等電位線在心墻、混凝土墊層、固結(jié)灌漿帶、帷幕內(nèi)密集。由此可見(jiàn)，防滲體系的效果是明顯的。

壩體和壩基的滲流是該項(xiàng)目工程防滲設(shè)計(jì)中特別關(guān)注的內(nèi)容之一。在正常蓄水位情況下，該項(xiàng)目工程大壩、壩基等部位的滲流情況見(jiàn)表3。在計(jì)算范圍內(nèi)，庫(kù)區(qū)總滲流量為54.31 L/s，年滲流量約867.45 m3(按水庫(kù)滿庫(kù)6個(gè)月計(jì)算)，占水庫(kù)庫(kù)容的2.47%。壩段滲流量為39.03 L/s，其中心墻滲透性為4.38×10-6cm/s，壩體滲流量為0.87 L/s，僅占?jí)味慰倽B流量的2.36%，壩段壩基滲流量為37.12 L/s，左岸山體滲流量為5.02 L/s，右岸山體滲流量為10.19 L/s。

表3 大壩、壩基等部位的滲流情況

通過(guò)水力梯度值針對(duì)壩體、壩基的一些重要部位展開(kāi)分析，以便于更深入評(píng)估壩體和壩基的穩(wěn)定性、安全性。 經(jīng)計(jì)算，心墻最大水力梯度值為2.317。有必要在心墻下游側(cè)設(shè)置反濾層，以防止黏土心墻出現(xiàn)裂縫時(shí)心墻土壤持續(xù)侵蝕和流失。混凝土蓋板、固結(jié)灌漿以及帷幕的水力梯度值均低于允許值，符合滲透穩(wěn)定性要求。左壩肩斷層以及右壩肩斷層的最大水力梯度值分別為1.323、0.541，2.520，為河床段斷層的最大水力梯度值。所以，在斷層下游的裸露區(qū)域安裝適當(dāng)?shù)倪^(guò)濾層勢(shì)在必行。

5.2 局部帷幕破壞研究

表4體現(xiàn)了部分滲流條件下帷幕的破壞情況。從中可以看出，在A2、A3和A4工況條件下的總滲流量比A1工況條件(帷幕完好)分別增加了3.61 L/s、3.05 L/s和5.40 L/s，增加百分比分別為6.7%、5.6%和10.0%，說(shuō)明斷層和溶蝕帶外帷幕灌漿的質(zhì)量對(duì)樞紐區(qū)的滲流有很大影響。同時(shí)還可看出，在部分帷幕破壞的情況下，斷層和溶蝕帶的水力梯度值比帷幕完整條件下大，不利于滲流的穩(wěn)定性。

表4 部分滲流條件下帷幕的破壞情況

6 結(jié) 語(yǔ)

BIM技術(shù)的推廣和發(fā)展，為實(shí)現(xiàn)建立模擬水利工程的有限元模型提供了無(wú)限可能，同時(shí)，也強(qiáng)化了數(shù)值模擬分析在工程建設(shè)以及安全評(píng)估時(shí)的重要影響。通過(guò)三維穩(wěn)定性有限元計(jì)算分析得知，壩基是此次項(xiàng)目工程滲漏最主要集中位置，左右岸岸坡是次要位置。為了防止斷層中填料的滲透損失，有必要在三條斷層的暴露下游鋪設(shè)過(guò)濾材料。研究發(fā)現(xiàn)，假如斷層及河床溶蝕帶的帷幕受損，也會(huì)影響到總滲流，對(duì)促進(jìn)斷層和溶蝕帶的穩(wěn)定性和安全性有不利影響。所以，應(yīng)該在帷幕施工時(shí)，對(duì)幕墻的施工質(zhì)量進(jìn)行嚴(yán)格管控。