子域積分法及組合網(wǎng)格算法在無(wú)壓滲流中的應(yīng)用.pdf

1、由于只需要進(jìn)行一次傳導(dǎo)矩陣的集成與分解，流量法在穩(wěn)定無(wú)壓滲流分析中獲得了廣泛的應(yīng)用。由于滲透系數(shù)在自由面上、下的間斷性質(zhì)，導(dǎo)致初流量法中，在與自由面相交的單元中結(jié)點(diǎn)初流量的被積函數(shù)是不連續(xù)的。采用傳統(tǒng)的Gauss積分方式時(shí)，這些單元內(nèi)數(shù)值積分結(jié)果可能出現(xiàn)單元部分區(qū)域的貢獻(xiàn)被忽略、被夸大或自由面在數(shù)值積分點(diǎn)在自由面上、下的輕微移動(dòng)導(dǎo)致結(jié)點(diǎn)初流量數(shù)值積分結(jié)果出現(xiàn)突變等問(wèn)題，降低數(shù)值計(jì)算的收斂性與穩(wěn)定性。
　　 初流量法在非穩(wěn)定流中應(yīng)

2、用時(shí)，也存在著滲透系數(shù)和貯水率在自由面上、下的間斷性質(zhì)引起的非連續(xù)被積函數(shù)。這2項(xiàng)的數(shù)值積分計(jì)算可能出現(xiàn)較大誤差，降低了數(shù)值計(jì)算的收斂性與穩(wěn)定性。對(duì)這些非連續(xù)項(xiàng)的積分計(jì)算進(jìn)行改進(jìn)，對(duì)于提高數(shù)值計(jì)算的穩(wěn)定性和收斂性，促進(jìn)初流量法的進(jìn)一步推廣應(yīng)用十分有意義。
　　 在無(wú)壓穩(wěn)定滲流分析過(guò)程中，常常會(huì)遇到某一方向尺寸與整體計(jì)算區(qū)域相比非常小的線形結(jié)構(gòu)。當(dāng)這些結(jié)構(gòu)中的水流運(yùn)動(dòng)不能被忽略的時(shí)候，需要對(duì)這些結(jié)構(gòu)需要進(jìn)行精細(xì)離散，但這些結(jié)構(gòu)的厚

3、度可達(dá)厘米級(jí)甚至毫米級(jí)，為保持計(jì)算精度必然導(dǎo)致其它區(qū)域也需要精細(xì)的離散，一方面造成網(wǎng)格離散的工作量和計(jì)算量大幅增加，另一方面造成有限元求解精度降低、收斂精度較差。針對(duì)這些問(wèn)題，提出一種合理可行的處理方法對(duì)于減小網(wǎng)格離散工作量和難度，獲取線形結(jié)構(gòu)中滲流場(chǎng)精細(xì)分布的同時(shí)提高計(jì)算的速度和精度有十分重要的意義。
　　 本文針對(duì)上述問(wèn)題進(jìn)行了研究，提出了子域積分理論、組合網(wǎng)格法在無(wú)壓滲流中的應(yīng)用。采用FORTRAN語(yǔ)言編制了相應(yīng)的計(jì)算程序

4、，并對(duì)程序及理論的可靠性進(jìn)行了驗(yàn)證。
　　 本文的主要研究?jī)?nèi)容和創(chuàng)新性成果包括以下三個(gè)方面:
　　 (1)基于子域積分方法，對(duì)無(wú)壓穩(wěn)定滲流的初流量法進(jìn)行了改進(jìn)。將與自由面相交的單元分為自由面上、下兩部分，再根據(jù)這兩個(gè)區(qū)域的具體形狀將其細(xì)化為標(biāo)準(zhǔn)形式的區(qū)域—子域，然后把每個(gè)子域看做一個(gè)常規(guī)的“單元”，在每個(gè)子域內(nèi)按照標(biāo)準(zhǔn)Gauss積分方式布置數(shù)值積分點(diǎn)，從而保證每個(gè)數(shù)值積分點(diǎn)控制的區(qū)域內(nèi)結(jié)點(diǎn)初流量的被積函數(shù)都是連續(xù)的。避免

5、了結(jié)點(diǎn)初流量被夸大、忽略或者自由面在積分點(diǎn)上、下的微小移動(dòng)導(dǎo)致的數(shù)值積分結(jié)果出現(xiàn)突變的情況，從而改善了數(shù)值計(jì)算的穩(wěn)定性和收斂性。文中給出了與自由面相交的平面4邊形和3維6面體單元的基本類(lèi)型及其具體子域劃分方式，同時(shí)給出了非基本類(lèi)型單元和基本類(lèi)型單元的變換關(guān)系。編制了基于子域積分的改進(jìn)初流量法穩(wěn)定無(wú)壓滲流程序IIFM2DS(平面滲流)和IIFM3DS(三維滲流)。
　　 (2)借鑒飽和非飽和滲流理論以及Desai等人的研究，將適用

6、于穩(wěn)定滲流的初流量法推廣到非穩(wěn)定滲流，根據(jù)虛位移原理推導(dǎo)了有限元離散格式，并采用子域積分對(duì)其中的非連續(xù)項(xiàng)的數(shù)值積分進(jìn)行改進(jìn)，編制了相應(yīng)的二維有限元程序IIFM2DT和三維有限元程序IIFM3DT。
　　 (3)將組合網(wǎng)格法應(yīng)用于無(wú)壓穩(wěn)定滲流分析中。采用兩套獨(dú)立的網(wǎng)格進(jìn)行模擬，對(duì)整體區(qū)域采用尺寸較大的粗網(wǎng)格進(jìn)行模擬（不考慮線性結(jié)構(gòu)的影響），對(duì)線形結(jié)構(gòu)則采用尺寸比較小的細(xì)網(wǎng)格進(jìn)行模擬。細(xì)網(wǎng)格通過(guò)流量修正對(duì)粗網(wǎng)格滲流場(chǎng)進(jìn)行調(diào)整，計(jì)算在