唐振華，尹明玉，張小濤，趙鵬強

（中水東北勘測設(shè)計研究有限責任公司，吉林 長春130021）

0 引 言

堤防設(shè)計中，滲流分析是一項非常重要的內(nèi)容，通常都采用數(shù)值模擬的方法，主要包括有限差分法、有限單元法、邊界元法以及離散單元法等。由于有限元法能夠建立各種復雜的幾何模型、具有強大的前、后處理功能等特點，其越來越多的被應用于滲流分析當中。

作為成熟有限元分析程序的ANSYS 軟件中雖然不存在專門的滲流分析模塊，但由于滲流場與溫度場的有限元計算公式相同，均屬于拉普拉斯定解問題，因此可以利用ANSYS 中的熱分析模塊進行分析，只需相應參數(shù)（熱傳導系數(shù)—滲透系數(shù)、溫度邊界條件—水頭）對應采取即可。

1 計算原理

對于各向異性、非均質(zhì)的連續(xù)介質(zhì)，服從達西定律的穩(wěn)定滲流問題，可歸結(jié)為下列定解問題：

式中：H 為水頭函數(shù)；kx，ky，kz分別為x，y，z 主方向的滲透系數(shù)，坐標軸方向與滲透主方向一致；Ω 為滲流區(qū)域；S1為已知水頭值的邊界曲面；S2為給定流量邊界曲面；S3為浸潤面；S4為逸出段；q 為邊界上的單位面積流量，這里q=0 表示為無流量交換邊界，n 為邊界的外法線方向。對于各向同性的介質(zhì)，即kx=ky=kz=k，公式（4）可簡化為

ANSYS 熱分析模塊中穩(wěn)態(tài)傳熱微分方程為：

式中：T 為溫度；kx，ky，kz分別為x，y，z 主方向的介質(zhì)熱傳導率。比較滲流微分方程（式1）與傳熱微分方程（式5）可以發(fā)現(xiàn)二者非常相似，滲流分析問題可以按照溫度場問題的特殊形式考慮。

利用ANSYS 模擬滲流分析時，采用固定網(wǎng)格法求解自由水面位置（浸潤線）?；贐athe 提出的單元滲透矩陣調(diào)整法的原理上，通過改變自由面以上的單元的滲透系數(shù)，來調(diào)整這部分單元對總體滲透矩陣的貢獻，從而通過逐步迭代使假設(shè)自由面趨近于實際自由面。

2 計算實例

2.1 計算模型及參數(shù)

這次計算以某管道穿越西里尼西河為例。

通過分析設(shè)計部門提供的西里尼西河穿越縱斷面圖和地質(zhì)勘察報告之后可以發(fā)現(xiàn)，管道基本布置在第一層的雜色卵石地層中，根據(jù)《土力學》及勘察報告提供的地層顆粒組成情況，確定各地質(zhì)層的滲透系數(shù)分別為地面以上堤防土體，k1=10 m/d；卵石層，k2=100 m/d；礫砂層，k3=60 m/d；全風化花崗巖，k4=0.1 m/d；施工回填土，根據(jù)施工情況分別取卵石層滲透系數(shù)的1，1.2，1.5，2，5 和10 倍；管道，k5=0 m/d，即不透水。

計算時以50年一遇設(shè)計洪水位523.46 m 作為堤外水位，堤內(nèi)地面無水。地表以下20 m 深處巖石視為相對不透水邊界，作為下限，水平方向上堤防前后均取20 m，沿堤防軸線方向取30 m。模擬中依次建立管道、回填土體（形狀按開挖尺寸定義）、堤防、地基單元，通過改變回填土體的滲透系數(shù)，來模擬工程實施前和實施后按不同施工質(zhì)量復堤后堤防的浸潤線、最大滲透坡降和單寬滲流量（滲流計算簡圖見圖1）。

圖1 滲流計算簡圖

2.2 結(jié)果分析

2.2.1 工程實施前滲流分析

圖2給出了現(xiàn)狀堤防在工程實施前的壓力水頭及等勢線云圖，根據(jù)計算結(jié)果，堤防的最大滲透比降為0.476，管道穿越處筑堤材料的允許滲透比降為0.50～0.80，即發(fā)生50年一遇洪水時，堤防的最大滲透比降小于材料的允許滲透比降，因此堤防不會發(fā)生滲透破壞。

圖2 現(xiàn)狀堤防滲流壓力水頭及等勢線云圖（單位：m）

2.2.2 工程實施后滲流分析

根據(jù)施工情況將管道周圍回填土體和堤防的滲透系數(shù)，按原土體滲透系數(shù)的1.2，1.5，2，5 和10 倍取值，對工程實施后的堤防滲流情況進行分析，分別繪出了自由水面以下的滲流體網(wǎng)格圖、堤防內(nèi)回填土體單元網(wǎng)格和流速場圖及堤防內(nèi)原土體單元網(wǎng)格和流速場圖（圖略，以k后為k前的5 倍為例）。

計算結(jié)果表明，當回填土體的滲透系數(shù)增大時，堤防內(nèi)部滲流浸潤線略有抬高，但增幅不大。從滲流量增大幅度、最大滲透比降與回填土體滲透系數(shù)的關(guān)系圖（見圖3、圖4）可以看出，滲流量增大幅度與施工前后回填土體滲透系數(shù)的比值成線性關(guān)系，而施工前后回填土體的滲透系數(shù)之比小于7 時，堤防背水坡的滲透坡降均小于筑堤材料的允許滲透比降，不會發(fā)生滲透變形，是安全的。但在該滲透系數(shù)比大于7 時，堤防背水坡的滲透坡降將可能大于筑堤材料的允許滲透比降，造成堤防的滲透破壞。

圖4 最大滲透比降與回填土體滲透系數(shù)關(guān)系

3 結(jié) 語

通過滲流場與溫度場的相似性比較，發(fā)現(xiàn)ANSYS 軟件中的熱分析模塊可以用來進行滲流分析。計算實例表明，該軟件可以很好的模擬滲流場隨時間的變化過程，并能夠反映實際的滲流狀態(tài)。而管道穿越工程對堤防的影響，歸根結(jié)底可將回填土體的滲透系數(shù)作為衡量工程質(zhì)量的指標，研究堤防滲透比降和滲流量與管道周圍土體的滲透系數(shù)之間的變化關(guān)系可以作為工程提供控制指標。施工中可參照本結(jié)果進行相應的防滲措施以預防堤防的滲流量加大或造成滲透破壞。

［1］ 毛昶熙.滲流數(shù)值計算與程序應用［M］.南京：河海大學出版社，1999.

［2］ 齊曉華，李曉麗，張波.基于有限元法的大壩滲流分析［J］.內(nèi)蒙古農(nóng)業(yè)大學學報，2012，33（3）：211-215.

［3］ 戚藍，黃君寶，黃曉東.ANSYS 軟件在有垂直防滲墻的堤壩及地基滲流計算中的應用［J］.長沙交通學院學報，2005，21（3）：48-51.

［4］ 許尚杰，黨發(fā)寧，田威.基于ANSYS 的土壩滲流分析及其應用［J］.西安理工大學學報，2008，24（4）：467-471.

［5］ 羅啟北，郝方，肖良錦，張艷霞. ANSYS 在滲流場動態(tài)分析中的應用［J］.貴州水力發(fā)電，2005，19（1）：68-71.