碾壓混凝土壩滲流場分析的縫面滲流平面單元模擬法

1、碾壓混凝土壩滲流場分析的縫面滲流平面單元模擬碾壓混凝土壩滲流場分析的縫面滲流平面單元模擬法法摘要：針對碾壓混凝土壩成層施工的結構特點及其滲流特 性，筆者在非均質各向異性等效連續(xù)體模型的基礎上，提出可精細而又方便地模擬大壩中層 面、縫面和裂縫等滲流行為的縫面滲流平面單元模型。在沒有增加 計算 網格結點的前提下， 不但能高效率地刻畫出層面和縫面所帶來的壩體強滲透各向異性的特點，而且又能極容易地 對那些局部零星分布的有強集中滲透能力的冷縫和裂縫等縫隙進行逐個細致地模擬，壩體的 滲流實際情況在數值模型中得到了甚為客觀的展現，提高了大壩滲流場的計算 分析 精度，基 本上完整地解決了碾壓混凝土壩的滲流建模

2、 問題 。最后，較詳細地進行了一個高 102m 的大壩 的算例分析。關鍵詞：碾壓混凝土壩 滲流 層面 縫面 裂縫 縫面滲流 平面單元在碾壓混凝土成層onclick=g(結構);結構中，混凝土本體為弱透水介質，而層面和縫面為相對強透水onclick=g(結構);結構面，通?？蓪δ雺夯炷翂误w用水力等效連續(xù)體模型和非連續(xù)體裂隙模型進行滲流特性建模14，其中對冷縫面或裂縫采用空間薄層單元進行模擬。由于正常的施工層面和縫面的間距相對于大壩的特征尺寸很小，相對簡單的已有成熟的 理論 基礎和豐富運用經驗的等效連續(xù)體模型已得到了廣泛的推廣 應用 ，但這種模型中存在的很大的弱點是對上述那些單個地出現的冷縫面和

3、裂縫有時不能進行很好地模擬，薄層單元的引入在一定程度上不但擴大了解題的規(guī)模，而且還會 影響解題的精度。筆者針對這個問題和碾壓混凝土壩中存在的局部集中滲漏現象較普遍的事實，提出無厚度縫面平面單元，用來專門模擬壩體中的冷縫面和裂縫等集中滲流通道，并融入等效連續(xù)體模型中，從而達到既整體又突出局部地刻畫和模擬碾壓混凝土壩的滲流行為和高精度地進行滲流場的計算分析。可在沒有增加網格結點的前提下，能高效率地反映出碾壓混凝土壩中正常施工層面和縫面的滲流影響，而且又能對那些具有強透水性零星分布的冷縫面或onclick=g(結構);結構性裂縫的滲流行為進行單獨細致的模擬。1 1基本理論基本理論1.11.1 滲流基

4、本理論滲流基本理論非均質各向異性多孔隙連續(xù)體介質中的穩(wěn)定飽和滲流連續(xù)微分控制方程為(1)式中：xi為坐標，i=1，2，3；kij為達西滲透系數張量，刻畫壩體的滲透各向異性程度；h=x3+p/為總水頭，x3為位置水頭，p/為壓力水頭。邊界條件為(2)式中：h1為已知水頭函數；ni為滲流邊界面外法線方向余弦；i=1，2，3；1、2、3和4分別為第一類和第二類滲流邊界，以及滲流自由面和滲流逸出面；qn為法向流量，流出為正。1.21.2 立方定律立方定律若將有集中滲流能力的縫面和裂縫視為具有一定水力隙寬的縫隙onclick=g(結構);結構面，并認為onclick=g(結構);結構面中的水流流態(tài)為層流

5、，則可以用層流縫隙流“立方定律”來描述onclick=g(結構);結構面的滲流行為，即(3)q = df= kfdfI(4)式中：為onclick=g(結構);結構面中的平均流速；df為onclick=g(結構);結構面的水力等效隙寬；I 為onclick=g(結構);結構面中的水力梯度；q 為onclick=g(結構);結構面中的單寬流量；kf為水力等效滲透系數。1.31.3 縫面平面單元模型縫面平面單元模型因碾壓混凝土本體的透水能力很小而縫面(冷縫和裂縫)的透水能力與其水力隙寬的立立成正比，縫面無論在其切向或法向上的透水能力均為相對很大，尤其是在縫面法向方向上其滲透系數為一相對大值，通?？p

6、面需單獨劃分成具有一定厚度的薄層單元。但是縫面的隙寬一般很小，只有幾十m，甚至更小，單元在縫面法向方向的尺寸相對于縫面切向方向的尺度甚小，這種薄層單元在理論上精度差，有時甚至會嚴重降低整個滲流場的求解精度。另外，薄層單元需單獨劃分出來，增加了解題的規(guī)模。這里提出縫面縫隙滲流的無厚度縫面平面單元，可用來專門解決具有集中滲流能力的零散性縫面的滲流模擬問題。因縫面的水力隙寬很小，法向的透水能力又極大以及混凝土本體透水能力極小，對于工程滲流問題而言，完全可以認為在縫面內縫面法向的水頭損失為零，縫面中的水流呈準二維滲流狀態(tài)，因此有(5)式中：為與縫面相關的局部坐標；為縫面平面單元的二維滲透系數張量，反映

7、縫面的滲透各向異性和縫面的透水能力，其中可借助上述“立方定律”縫面的水力隙寬 df 也被考慮在中。局部坐標與整體坐標 xi之間的關系為(垂直于縫面)(6)式中：、Te和xi分別為局部坐標向量，縫面平面單元的坐標轉換矩陣和整體坐標向量。任何一個縫面平面單元 e 的結點局部坐標得到確定后，就可進行這個單元的傳導矩陣元素的計算：(7)式中：sf為縫面單元域；Nr 和 Ns 為縫面單元的插值函數；m 為縫面單元的結點數。這種縫隙滲流縫面平面單元對壩體內滲流場特性的影響是通過下式的水頭連續(xù)條件來實現的。(8)即縫面 sf任一點處碾壓混凝土壁面上的水頭與無厚度縫面平面單元同一位置處的水頭是相同的。再按通常

8、的有限單元法要求，據計算域中任一結點處的流量平衡條件式(12)就可以組裝成常規(guī)形式的求解整個滲流場的有限單元法支配方程。(9)式中：n 為 總結 點數；QRCC和 Qf分別為由碾壓混凝土壩三維等效連續(xù)體單元和二維縫面平面單元對結點 i 所作貢獻的等效結點流量；式(12)須對環(huán)繞 i 結點的所有常規(guī)單元和縫面單元求和。用上述縫面平面單元模擬冷縫面或裂縫滲流行為的優(yōu)點有：(1)二維縫面平面單元沒有厚度，網格剖分時無需對它進行專門剖分，只需額外形成一個簡單的所有這類二維單元的單元信息表即可，不增加額外結點和自由度，幾乎不增加網格剖分和求解時的計算工作量；(2)避開了在縫面法向隙寬方向上的微分及積分運

9、算，不因對縫面或裂縫滲流行為的處理而給整個滲流場求解帶來大的誤差；(3)能極方便地正確模擬碾壓混凝土壩中局部零星散落的裂縫、層面和縫面的集中滲流行為；(4)此 方法 的引入，碾壓混凝土壩的滲流數值建模問題已得到了基本完整的解決；(5)能考慮縫面的滲流各向異性。2 2 碾壓混凝土壩滲流場碾壓混凝土壩滲流場 分析分析 的有限單元法求解的有限單元法求解2.12.1 滲流場求解的有限單元法滲流場求解的有限單元法就固定網格求解有自由面的滲流 問題 的有限單元法的結點虛流量法可詳見 文獻8，9。此外，在實際工程中某些滲流控制措施、邊界條件和滲透系數張量的確定也往往是求解精度的主要 影響 因素，需分別專門對

10、待。因混凝土大壩滲流場的水頭分布主要受控于排水措施的作用，解題時必須對 計算 域中的所有排水措施(比如排水孔、排水洞等)進行精細的模擬。2.22.2 排水孔滲流開關器排水孔滲流開關器排水孔排水幕的滲流行為可嚴密地用排水子onclick=g(結構);結構來處理，但不同類型的排水孔在滲流邊界條件上的不同表現導致程序中對它們處理方式的不同。因工程滲流問題的復雜性，事先不知道壩基每一個溢流型排水孔的頂面高程一定是低于或高于鄰域自由面的高程，就象是確定逸流型排水孔內滲流逸出線位置那樣，也需在迭代求解過程中據中間解逐步給予確認，為此需引進排水孔開關器的概念10。在每一個溢流型排水孔頂面虛構一個數學開關器，

11、當開關器打開時排水孔工作，孔壁面為一個等水頭面；當開關器關閉時排水孔完全不工作，即此時孔頂面高程高于鄰域自由面位置，排水孔被自由面穿過。這樣解題就可以確保每一個溢流型和逸流型排水孔的真實工作狀態(tài)都得到正確的甄別和模擬。3 3 計算分析計算分析3.13.1 概況概況某碾壓混凝土重力壩，壩高 102m，底寬 70m，頂寬 12m，上游水位 100m，下游水位20m，在壩體中距上游面 23m 處設置孔徑 15cm，孔距 4m 的垂直排水幕，壩基設置 2 道深分別為 30m、15m 的主排水幕，中間 2 道深為 8m 的抽排水幕，其孔口高程為基礎廊道底版高程，孔距、孔徑都與壩體排水孔相同。假定壩體內自

12、下而上每隔 20m 存在一個冷縫面，共 4 個，壩上、下游面分別設置 03m 厚的變態(tài)混凝土，上游面變態(tài)混凝土后有厚 2m 的二級配混凝土，建基面上設有 2m 厚的常態(tài)混凝土墊層。計算中當考慮防滲體水平開裂時，裂縫都貫穿于壩上游面的變態(tài)混凝土和其后的二級配碾壓混凝土；壩上游面垂直劈頭縫位置假定發(fā)生在兩相鄰排水孔的中間。圖 1 為大壩的剖面圖，利用水力對稱性計算域在壩軸線向厚 2m。計算域材料的滲透系數見表 1 所示，其中 k和 k分別為壩體層面切向和法向的主滲透系數。表 1 計算域材料的滲透系數(單位：cm/s)滲透各向性變態(tài)混凝土二級配混凝土三級配混凝土常態(tài)混凝土壩基帷幕k=110-9k=1

13、10-7k10-9k=110-6k10-9k=510-9k=110-4(-60m 以上)k=110-5(-60m 以下)k=110-53.23.2 計算工況計算工況這里主要 研究 壩體變態(tài)混凝土防滲體的水平向或豎直向的開裂及 4 個典型集中滲水縫面不同水力隙寬時壩體滲流場的變化情況。具體的計算工況見表 2 所示，裂縫均貫穿二級配混凝土。表 2計算工況工況號結果圖號工況說明12345圖 2圖 3圖 4圖 5圖 6上游面防滲體無裂縫和壩體無縫面上游面防滲體有隙寬 0.2mm 的水平裂縫，壩體有隙寬 0.01mm 縫面但無排水幕上游面防滲體有隙寬 0.2mm 的水平裂縫，壩體有隙寬 1mm 縫面且有

14、排水幕上游面防滲體有隙寬 0.2mm 的垂直劈頭裂縫，壩體有隙寬 1mm 縫面但無排水幕上游面防滲體有隙寬 0.2mm 的垂直劈頭裂縫，壩體有隙寬 1mm 縫面且有排水幕3.33.3 計算結果分析計算結果分析圖 2 滲流場水頭線分布的主要特點有：(1)壩上游面變態(tài)混凝土防滲體起到一定的防滲作用，但僅占壩上、下游面總水頭的 10%20%；(2)盡管二級配碾壓混凝土的層面切向滲透系數比三級配碾壓混凝土的小一個數量級，但防滲效果并不明顯；(3)壩體滲透各向異性比使得等水頭線向下游折彎，因三級配碾壓混凝土各向異性比達 3 個數量級，在三級配碾壓混凝土區(qū)等水頭線就更呈現得水平向狀態(tài)，水平向的滲透水力阻力

15、很??；(4)下游面變態(tài)混凝土的阻滲作用非常明顯，若沒有它等水頭線會呈現幾乎完全水平的狀態(tài)，滲透水流會直接從壩下游面逸出，逸出線位置會很高；(5)壩下游面有變態(tài)混凝土后，在實際工程中壩內排水設施務必要布置得充分一些；(6)因變態(tài)混凝土很不透水，壩下游面有蒸發(fā)作用，一般壩下游面不會留有滲流逸出水跡現象，改善了大壩運行期的觀瞻性。圖 1 碾壓混凝土壩剖面及滲流場計算域圖 2 壩體無裂縫和冷縫面時等水頭線分布(單位：m)圖 3 與圖 2 的結果對比，盡管縫面的水力隙寬只 0.01mm，亦即赤眼是完全看不到縫面的存在，但是由于碾壓混凝土本體的滲透性極小，這樣的縫面已明顯地形成了壩中的主要集中滲流通道，不

16、同高程上的縫面都已承受了幾乎為 100%的庫水頭揚壓力。這一現象直觀地說明，壩中縫面處理不當對壩體的安全穩(wěn)定性可能有著巨大的潛在不安全因素，以及壩上游面具有優(yōu)質完整的防滲onclick=g(結構);結構的重要性。圖 4 結果表明：(1)壩內水頭線呈水平狀分布，排水幕不但對壩體進行了充分的排水降壓作用，而且也充分地截止了從水平縫面集中滲入的高壓力庫水，確保了大壩的安全性；(2)在壩底部距壩建基面 20m 高程處，由于庫水壓力大和防滲體裂縫與壩內縫面的直接相連，排水幕未能完全截止從庫水滲入的水流，縫面上殘留有約 1315m 水頭的揚壓力；(3)比較圖 3 和圖4，得知確保排水幕能正常工作對大壩安全

17、性是何等的重要，建議工程實踐中孔徑可適當大一點，不宜小于 15cm，確保排水孔的永久暢通性；(4)在圖 4 極端情況下，壩底部縫面上才有一定的揚壓力，說明 目前 工程實踐中排水孔間距取 34m 是較合適的。工況 4 的結果為圖 5 所示。由于大壩上游面劈頭立面縫自頂部至底部貫穿于整個上游面防滲體和二級配碾壓混凝土，裂縫與壩內水平縫面得到了充分的連通，再加上壩體為強滲透各向異性體，壩下游面又有無裂縫的變態(tài)混凝土，幾乎整個壩體的縫面和層面上都承受了 100%滿庫水頭的揚壓力，這種情況是絕對不允許出現的。圖 3 壩體有水平裂縫和縫面、無排水幕時等水頭線分布(單位：m)圖 4 壩體有水平裂縫和縫面、有

18、排水幕時等水頭線分布(單位：m)圖 5 壩體有劈頭縫和縫面、無排水幕時等水頭線分布(單位：m)圖 6 壩體有劈頭縫和縫面、有排水幕時等水頭線分布(單位：m)圖 6 為工況 5 的結果，(1)即使在圖 5 極端情況下，當壩內設置了排水幕后，壩體滲透水流還是得到了較充分的排水降壓，除縫面局部區(qū)域外，層面上揚壓力幾乎為零；(2)由于劈頭立面縫與壩內縫面直接相連，盡管有排水幕的疏干作用，底部仍有部分集中滲漏現象，但此時縫面上揚壓力仍得到了較好的控制，水頭只有 1518m；(3)基于壩體底部冷縫面區(qū)有一定大小的揚壓力的事實，當壩下游面有變態(tài)混凝土時，建議最好在壩體的底部范圍內再添置 12 道具有一定高度

19、的排水幕，充分疏干那里的滲透水流，尤其是在高壩建設中更要注意這一點。圖 7 20m 高程處壩體縫面上的揚壓力水頭分布圖 7 為各工況在壩體底部高程 20m 處縫面上的揚壓力水頭的分布圖。4 4 結論結論(1)基于碾壓混凝土壩的特點，提出了縫隙滲流無厚度縫面平面單元及其滲流計算模型，在沒有增加網格結點的前提下，能對那些強透水性縫面或裂縫的滲流行為進行單獨細致的模擬，成功地解決了各種縫面滲流的建模問題，提高了對碾壓混凝土壩滲流真實情況的數值模擬精度。(2)碾壓混凝土壩層面揚壓力分布主要由排水幕中排水孔的排水降壓作用所控制。因此，大壩運行過程中應確保排水幕的永久暢通性，孔徑不宜太小。(3)在目前的碾壓混凝土重力壩上、下游面大多有變態(tài)混凝土和壩中設有排水幕的壩體滲控方案中，壩體層面和縫面上的揚壓力能得到滿意的控制。對于那些高壩，建議在大壩的底部近上游區(qū)域中再添置 12 道排水幕。參參 考考 文文 獻：獻：1朱岳明，張燎軍，黃文雄，等.龍灘碾壓混凝土重力壩壩基及壩體滲流與排水對大壩安全度影響的研究(“八五”國家 科技 攻關子題成果)R.南京：河海大學水利水電工程學院，1995.2朱岳明，龔道勇，張燎軍，等.高碾壓混凝土重力壩滲流特性和計算 理論 及 應用 研究(“九五”國家科技攻關子題成果)R.南京：河海大學水利水電工程學院，20