連續介質滲流區域非線性問題數值求解分析

2020-12-11 01:32:02潘廣良

黑龍江水利科技 2020年11期

關鍵詞：模型

潘廣良

(方正縣河湖運行保障中心，黑龍江方正 150800)

0 前言

土石壩的非線性滲流求解在于滲流自由面和溢出面的確定。文章基于移動網格法、單元滲透矩陣調整法和復合單元高斯點法，運用Microsoft Developer Studio可視化集成開發二維無壓滲流分析有限元程序，結合具體算例求解連續介質滲流區域的非線性問題，成果對滲流區域非線性問題研究具有重要意義。

1 緒論

中國目前的巖土體水力學研究，為實際應用和方便處理考慮，基本選用線性數學模型。線性數學模型通常將滲流運動和物體狀態簡化處理，因此助推巖土體水力學近些年快速發展[1]。但實際滲流由多種因素復合影響而成，絕大多數的滲流分析均呈非線性，故材料的非線性問題在巖土體水力學研究中為基本問題，具有十分重要的研究價值。

2 滲流基本原理和數學模型

連續介質滲流區域非線性問題以二維非穩定流的飽和滲流模型為例，4種表達方程式[2]如下：

1)Boussinesq方程：在自由面進行流量補給，并將緩變形式簡化為水平向滲流，在區域地下水活動研究中應用廣泛。

(1)

2)Laplace方程：自由面既可進行流量補給，同時為下降流速的限值條件，表達式為：

(2)

3)擴散方程：自由面不進行流量補給，該方程僅針對自由面微弱變化的均質土壩，表達式為：

(3)

4)固結方程：因考慮到土體具有壓縮性，常應用于粘土筑壩的固結，表達式為：

(4)

3 算例分析

3.1 移動網格法

連續介質滲流的二維有限元模型分為2種不同單元進行劃分，三角形單元和8節點四邊形等參單元[3]。算例為一均質土壩，壩面寬6m，上、下游水位分別為4m和3m，計算模型示意圖，見圖1。

圖1 計算模型示意圖

邊界條件：滲流邊界BC無補給，AD為封閉邊界。

模型進行三角形單元劃分時，共劃分18個單元，16個結點；八節點四邊形單元劃分時，共劃分為9個單元，40個結點，三角形單元網格剖分，見圖2；八結點四邊形等參元網格剖分，見圖3。

圖2 三角形單元網格剖分圖3 八結點四邊形等參元網格剖分

通過程序模擬計算，得到2不同單元的滲流自由面結點水頭計算結果值，區域滲流自由面結點水頭對比值，見表1。

表1 區域滲流自由面結點水頭對比值

3.2 單元滲透矩陣調整法[4]

該法的無壓滲流分析單元劃分原理同基本網格法。算例為一均質土壩，高8m，寬4m，上、下游水位6.0m和1m，該均質土壩模型，見圖4。

圖4 均質土壩計算模型

邊界條件：BC邊界無流量補給，流量為零，CD邊界流量可能溢出，AE為封閉邊界。本算例先選用四邊形單元剖分模型，共剖分24個單元，93個結點。四邊形單元網格剖分，見圖5。

圖5 四邊形單元網格剖分

為保證迭代的穩定性和結果的精確性，進行5次迭代后計算得出溢出點水頭值5.68m，滲流自由面各節點水頭統計值，見表2。

表2 滲流自由面結點水頭統計值

由結果可得，在相同節點不同數解值有較大差距數，故在四邊形單元基礎上再次進行三角形網格加密，共劃分單元48個，結點35個，三角形單元網格剖分，見圖6。

圖6 三角形單元網格剖分

本次求解進行10次迭代后得到滲流面溢出點水頭值為3.23m，模型的滲流自由面和節點水頭結果見圖7和圖8，滲流自由面各節點水頭統計值，見表3。

圖7 計算模型滲流自由面

圖8 節點水頭等值線圖

表3 滲流自由面結點水頭統計值

由上述結果可得，網格加密的滲流求解更加貼合實際，為更好證明自編程序的合理性，重新運用四邊形單元加密劃分模型，本次單元劃分單元150個，節點501個，模型單元劃分結果見圖9。迭代9次后得到自由面水頭值為3.78m，計算模型滲流自由面，見圖10；節點水頭等值線圖，見圖11；滲流自由面結點水頭值，見表4。