Seep滲流計(jì)算V30使用手冊(cè)

1、seep V3.0離鶴炎醫(yī)來(lái)劍來(lái)二四年一月1滲流根本理論1.1水工滲流的危害及滲流分析計(jì)算的任務(wù)流體在多孔介質(zhì)中的運(yùn)動(dòng)稱為滲流。水是最為常見(jiàn)的流體， 水利水電工程中由于廣泛建造堤、壩、圍堰、水閘等擋水建筑物形成了水頭差，這些建筑物或其地基通常是透水的多孔 介質(zhì)，因此水工滲流現(xiàn)象十分普遍。水工滲流造成多方面的危害。滲流造成水庫(kù)、渠道水量損失；滲流使堤壩、圍堰土體飽 和，降低壩體的有效容重和抗剪強(qiáng)度，可能導(dǎo)致壩坡失穩(wěn)； 建筑物地基滲流對(duì)建筑物底部產(chǎn)生揚(yáng)壓力，也不利于建筑物的穩(wěn)定；滲流流速過(guò)大時(shí)，還可能造成壩體或建筑物地基的土體顆粒流失，發(fā)生滲透變形，從而使堤壩崩塌或建筑物滑移、傾覆；水庫(kù)滲流還可能

2、引起下游 地下水位升高，導(dǎo)致農(nóng)田冷浸漬害、 鹽堿化，使作物減產(chǎn)；攔污壩滲流造成地下水環(huán)境污染。水工滲流分析計(jì)算的任務(wù)就是要研究水在滲流區(qū)域的滲流流速、流量、水頭分布及浸潤(rùn)線等，從而為采取合理的滲流控制措施提供依據(jù)，以防止或減緩滲流危害。1.2達(dá)西定律19世紀(jì)50年代，法國(guó)工程師亨利達(dá)西 H.Darcy丨通過(guò)對(duì)裝在圓筒中的均質(zhì)砂土進(jìn) 展?jié)B透試驗(yàn)發(fā)現(xiàn)，通過(guò)兩個(gè)滲流斷面間的平均滲流流速，正比于兩斷面間的水頭差h,反比于滲徑長(zhǎng)度L,且與土粒構(gòu)造及流體性質(zhì)有關(guān)。這就是著名的達(dá)西定律， 可用公式表達(dá)為：vkhkdh kJ（）L ds2式中h 測(cè)壓管水頭，h z , z為位置高度，p為壓強(qiáng)，丫為水的容重。因

3、為2g2滲流的流速一般很小，流速水頭可忽略，故h z衛(wèi)。2gk反映土粒構(gòu)造及流體性質(zhì)的系數(shù)，即滲透系數(shù)，對(duì)于某一具體的流體比方水而 言，k值僅與土粒構(gòu)造有關(guān)。J滲透坡降，J 生。ds式中的負(fù)號(hào)-表示水總是流向水頭減小的方向。應(yīng)當(dāng)注意，達(dá)西定律中的流速是全斷面上的平均流速v，而不是土體孔隙中的流速 v，,v nv這兩種流速存在以下關(guān)系：(1.2.2)式中n為體積孔隙率，可見(jiàn)達(dá)西流速小于土體孔隙中的流速。還應(yīng)注意，達(dá)西定律只能適用于層流狀態(tài)的滲流運(yùn)動(dòng)。在水利工程中，除了堆石壩、堆石排水體等大孔隙介質(zhì)中的滲流為流之外，絕大多數(shù)滲流都屬于層流， 達(dá)西定律都可適用。對(duì)于非層流滲流，其流動(dòng)規(guī)律可用以下公式

4、形式表達(dá)：v kJm 上式中當(dāng)m=1時(shí)，為層流滲流；當(dāng) m=2時(shí)，為完全紊流滲流；當(dāng) 1 v m 2時(shí)，為層流到紊流 的過(guò)渡區(qū)。將式丨等號(hào)兩邊向x、y、z軸投影，便得到空間直角坐標(biāo)系中的達(dá)西公式:VxkxhxkxJxvykyhykyJyVzkzhzkzJz滲流運(yùn)動(dòng)連續(xù)性方程124從滲流場(chǎng)中取某一微分單元體如圖，其體積V dxdydz，流入左側(cè)面的水體質(zhì)量的速率為vxdydz，流出右側(cè)面的水體質(zhì)量的速率為面進(jìn)出流量之差為(vz vzdz)dxdy，那么左右 z一 vxdxdydzx同樣，對(duì)于前后面和上下面作流進(jìn)流出的流入量計(jì)算,最后累加各凈有流入量， 得到微元體內(nèi)總得流入量為VxVyVz V將上

5、式展開(kāi)為VxVyxyVXVy VzVx y z式中后一項(xiàng)與前一項(xiàng)相比，小得可以忽略，故改寫成質(zhì)量M隨時(shí)間的變化速率:MnVV、，n 、，nV -nV -ttttt式中，n為土體的孔隙率，p為水的密度，V為微元體的體積。式丨右邊三項(xiàng)分別代上式即為水體質(zhì)量在微元體內(nèi)積累的速率,根據(jù)質(zhì)量守恒原理,它應(yīng)等于微元體內(nèi)水體表土體骨架、孔隙體積及流體密度的改變速率。Vzv125引入彈性壓縮理論，可導(dǎo)出式中得等價(jià)表達(dá)為:即壓縮模量127M 2x/ hg n V - tt式丨中a 土顆粒骨架的壓縮性3水的壓縮性滲流測(cè)壓管水頭m，根據(jù)質(zhì)量守恒原理，式丨與式丨相等，那么得到VX Vy Vz128假定水和土不可壓縮時(shí)

6、，上式變?yōu)?Vx-Vy-Vz0x y z129式丨為不可壓縮流體在剛體介質(zhì)中流動(dòng)的連續(xù)性方程。滲流微分方程將達(dá)西公式代入滲流連續(xù)性方程可得:hhhkxkykzgx x y y z z式為非穩(wěn)定滲流微分方程的一般形式,既適合于承壓含水層,也適合于無(wú)壓滲流。式中Ss g n 稱為單位儲(chǔ)存量尺度為1/L，其含義是：?jiǎn)挝惑w積得飽和土體在水頭下降1m時(shí)，由于土體壓縮g丨和水體膨脹gn丨所釋放出來(lái)得儲(chǔ)存水量。當(dāng)介質(zhì)為均質(zhì)各向同性即kxkykzk丨時(shí)，式丨變?yōu)?.2.11 2h2h 2h Ss_hx2y2z2k tkj kk0x x y y z z2hh22x y2hz21213式、 丨分別為各向異性、各向

7、同性時(shí)的三維穩(wěn)定滲流微分方程。對(duì)于堤壩滲流及建筑物地基滲流，在許多情況下可簡(jiǎn)化為垂直剖面上的二維滲流問(wèn)題，那么式、 可變成kxkz0x22hz2(1.2.15)有自由面變動(dòng)的滲流微分方程對(duì)于土壩或堤防在水庫(kù)或河道水位下降時(shí)的非穩(wěn)定滲流情況，自由面下降引起的土體壓縮或彈性釋放水量與自由面下降時(shí)所排出的水量相比很小，故可令式kzz中的Ss 0，因此，有自由面的三維非穩(wěn)定滲流微分方程為其二維非穩(wěn)定滲流微分方程為kz h z z上面二式雖然在形式上與穩(wěn)定滲流方程式、 完全一樣，但結(jié)合自由面變動(dòng)的邊界條件所得到的水頭分布是空間座標(biāo)與時(shí)間的函數(shù)，而不象穩(wěn)定滲流方程式的解答，只是空間座標(biāo)的函數(shù)。水工滲流問(wèn)題

8、的定解條件發(fā)生在有限空間流場(chǎng)內(nèi)的滲流運(yùn)動(dòng)，不僅受滲流微分方程支配，也受流場(chǎng)邊界條件和初始滲流狀態(tài)條件支配。下面以土壩非穩(wěn)定滲流狀態(tài)介紹水工滲流定解條件的幾種類型。b H一4、-藝* t乂二 一-a-3：一亠.F_G圖1.2.2土壩非穩(wěn)定滲流定解條件初始條件：hh0 x,z,t僅對(duì)非穩(wěn)定滲流1218邊界條件:水頭邊界 h 1 f1 (x,乙t)在圖中，水頭邊界包括： 上游水頭邊界：hABC H,t)下游水頭邊界：hEFGH21219下游坡滲出段水頭邊界：hDE z流量邊界：q kn 2 f2(x, z,t)n在圖中，自由面AD下降時(shí)，由自由面流進(jìn)的單寬流量為:hq cost式中B自由面的法線與鉛

9、直線的夾角卩一給水度，即單位體積土體在飽和含水情況下，當(dāng)自由面下降A(chǔ)1m后排出的空隙水量，又稱土體的排水空隙率或有效空隙率。又因q knh, hkncosnz由式A B，自由面流量邊界條件改寫成hhkn tzBC當(dāng)壩體土各向異性時(shí)，kn取kz值。因滲流自由面與大氣相通，壓力水頭為零，故在自由面上還應(yīng)滿足:2 SEEP V3.0的功能和特色1適用于各種工況如水位固定或上游水位快降、復(fù)雜邊界包括各種防滲、排水設(shè)施 和復(fù)雜土層條件下的： 堤防、土石壩滲流計(jì)算；各類不透水材料壩如混凝土壩和漿砌石壩 水閘、泵房的地基滲流計(jì)算。2采用有限元法計(jì)算，可得到滲流量、浸潤(rùn)線、等水頭線等計(jì)算結(jié)果。3采用可視化界面

10、建立原始數(shù)據(jù)文件，也可在記事本中建立原始數(shù)據(jù)文件。4計(jì)算成果既可打印成表格，也可自動(dòng)繪制出浸潤(rùn)線、等水頭線圖。3 SEEP V3.0解題要點(diǎn)3.1坐標(biāo)系的選取及計(jì)算范圍確實(shí)疋如圖所示，水平坐標(biāo)x取向右為正，豎直坐標(biāo)與當(dāng)?shù)馗叱腆w系一致。對(duì)于上游入滲段長(zhǎng)度L1及下游滲出段長(zhǎng)度 L2，一般取0.51倍的水頭或取實(shí)測(cè)值圖 3.2單元剖分a) 單元自動(dòng)剖分的分組、分塊與結(jié)點(diǎn)編號(hào)根據(jù)滲流場(chǎng)的幾何形狀、土層分布等，將整個(gè)滲流場(chǎng)劃分為一些四邊形的子域，如圖所示。這些子域稱為塊，為使結(jié)點(diǎn)優(yōu)化并滿足自由面調(diào)整迭代的需要，將所劃分的 塊又分成假設(shè)干組，每組內(nèi)由假設(shè)干塊組成。具體規(guī)定：每組內(nèi)各塊首尾相接，上下貫穿整

11、個(gè)滲流場(chǎng)：每組內(nèi)的塊由上而下逐塊編號(hào)，且從左到右一組一組地進(jìn)展，但塊的結(jié)點(diǎn)標(biāo)號(hào)是任意的建議按照從左到右，從上到下的順序。圖中由結(jié)點(diǎn)3，7，11，12，15連成的折線為事先假定的浸潤(rùn)線。(451H13 坷2(3_2(3 (34()8 丄5 .12(533 (3(49 (3 138(5(31516 100)19 )專(3【7(39(3町(31?10214(520(322圖 b) 四條邊剖分點(diǎn)布置本軟件進(jìn)展?jié)B流有限元計(jì)算的根本微小單元為三角形單元，為此還需將上述四邊形塊單元剖分為三角形單元。如圖所示，從滲流場(chǎng)中取出任意一塊，圖中n1, n2, n3, n4分別表示相應(yīng)邊上剖分點(diǎn)數(shù)剖分點(diǎn)數(shù)決定三角形單

12、元的數(shù)目或大小，且按順時(shí)針?lè)较蚓幪?hào)，n1與n3、n2與n4分別為兩組對(duì)邊，兩對(duì)邊上的剖分點(diǎn)數(shù)可相等，也可不相等，但為 防止出錯(cuò)，最好使各對(duì)邊的點(diǎn)數(shù)相等或相差小一些。剖分點(diǎn)的疏密應(yīng)按水力條件來(lái)確定，在以下部位的剖分單元應(yīng)加密：水力坡降大的部位、壩軸線附近、截水槽及上游鋪蓋進(jìn)水處， 下游出口及排水溝附近。此外，剖分點(diǎn)數(shù)的擬定還必須滿足以下三點(diǎn)： 三角形單元總數(shù)500 ； 門3 厲門2,門4陽(yáng)門2n 1,門3 ； 浸潤(rùn)線變動(dòng)范圍內(nèi)n1 n3n1(3 n4n2(3n3(6)87圖322c) 剖分順序剖分順序依次為：劃分塊單元、分組并編組號(hào)、編塊號(hào)并核對(duì)塊單元的結(jié)點(diǎn)數(shù)、編結(jié)點(diǎn)號(hào)、在各塊邊上布置剖分點(diǎn)數(shù)。

13、同一塊內(nèi)的滲透系數(shù)必須一樣，但同一滲透系數(shù)土層可劃分為假設(shè)干塊。每個(gè)塊都必須是四邊形塊。在滲流區(qū)域邊界上有時(shí)不可防止地出現(xiàn)三角形塊，此時(shí)在該塊的邊界邊上任取一點(diǎn)作為結(jié)點(diǎn)，這樣一來(lái)，三角形塊可以看成是蛻化的四邊形塊。4 SEEP V3.0的操作步驟與文件內(nèi)容介紹4.1 DQBfc件的產(chǎn)生詳細(xì)的操作步驟如下：第一步：運(yùn)行CCED.EXE文件，彈出工程及視圖相關(guān)參數(shù)界面(如圖4.1.1):圖工程及視圖參數(shù)界面工程相關(guān)參數(shù)：(1) 填入“分塊組數(shù)和“土層種類數(shù)；(2) 等勢(shì)線間距有5%和10%兩個(gè)數(shù)可供選擇，分別表示在滲流場(chǎng)分布圖中相鄰兩條等 勢(shì)線之間的水頭勢(shì)能相差 5%和10%;(3) 滲流情況：如

14、果計(jì)算非穩(wěn)定滲流，需要填寫水位降落經(jīng)歷的時(shí)間、初始水位和降后 水位；如果計(jì)算穩(wěn)定滲流，需要選擇是否計(jì)算滲流量。視圖相關(guān)參數(shù)：(1) “最大坐標(biāo)和“最小坐標(biāo)：是您在繪圖過(guò)程中所能輸入的最大最小坐標(biāo)值，為 了方便繪圖，建議將斷面的最大X、Y坐標(biāo)值分別加上30,將斷面的最小 X、Y坐標(biāo)值分別減去30后輸入；(2) “1m=( )P 表示繪圖時(shí)一米長(zhǎng)度所顯示的象素個(gè)數(shù)，如果圖形中有比擬短的線條，可填入較大值，比方10 ;如果圖形中線條都較長(zhǎng)，建議采用默認(rèn)值（3）“點(diǎn)的選取半徑默認(rèn)為3個(gè)像素，設(shè)置大些可以更好地捕捉點(diǎn),線條的選取，建議采用默認(rèn)值；“字體大小一般默認(rèn)為10。4；但過(guò)大會(huì)影響(如圖 4.1.

15、2):-Ini x|第二步：?jiǎn)螕簟肮こ碳耙晥D參數(shù)界面中確實(shí)定按鈕，進(jìn)入圖形繪制界面文件（日編輯查看回幫助耳（斗000皿廠50伽）Line From：默認(rèn)狀態(tài)上下:進(jìn)展圖形繪制:圖圖形繪制界面圖形繪制：在右下角“Line From 的后面直接輸入節(jié)點(diǎn)坐標(biāo)并按回車鍵;如果節(jié)點(diǎn)坐標(biāo)有錯(cuò)誤，可以點(diǎn)擊【進(jìn)展修改;（3）點(diǎn)擊，用左鍵選中線條，再單擊右鍵可以刪除線條切換到狀態(tài)三進(jìn)展單元剖分:上加大巴減小“點(diǎn)的選取半徑，或者將圖形放大，然后點(diǎn)擊進(jìn)展剖分填寫，最后按塊的編號(hào)順序進(jìn)節(jié)點(diǎn)編號(hào)：按編號(hào)順序直接點(diǎn)擊該節(jié)點(diǎn)，如果不能選中該節(jié)點(diǎn)，可以點(diǎn)擊“點(diǎn)的選取半徑；剖分控制：先選中塊的四條邊如果不能選中該邊，可以點(diǎn)擊展

16、其它塊的剖分。數(shù)據(jù)填寫和文件保存：(1)輸入每層土的水平和垂直滲透系數(shù)K x、注意：每輸入一層土的滲透系數(shù)，都要點(diǎn)擊確定按鈕進(jìn)展數(shù)據(jù)保存；填寫土層滲透系數(shù)：點(diǎn)擊Ky cm/s，剖分塊信息修改：如果剖分信息有錯(cuò)誤，點(diǎn)擊 的編號(hào)，再單擊修改按鈕；曰二進(jìn)展修改先選中Reg下某一塊設(shè)置相關(guān)參數(shù)：點(diǎn)擊 匸翻開(kāi)相關(guān)參數(shù)設(shè)置對(duì)話框，檢查“分組數(shù)、 數(shù)間距和“滲流情況是否正確， 時(shí)，需填寫“滲流斷面穿過(guò)的塊號(hào)；“等勢(shì)線百分注意：當(dāng)計(jì)算工況為穩(wěn)定滲流且計(jì)算滲流量的彈出dqb文件保存對(duì)話框，選擇保存路徑后點(diǎn)擊保存（S）保存dqb.dat文件：點(diǎn)擊按鈕，注意：計(jì)算時(shí)文件名必須為 dqb.dat ；生成斷面輪廓圖：如

17、果需要斷面輪廓圖，直接在圖中點(diǎn)擊鼠標(biāo)左鍵選中所有要保存的輪廓線后，點(diǎn)擊彈出保存對(duì)話框，填入文件名并保存。4.2運(yùn)行計(jì)算點(diǎn)擊 耀按鈕，彈出計(jì)算操作界面(如圖4.2):圖4.2計(jì)算操作界面(1) 點(diǎn)擊查看 DQB按鈕，檢查產(chǎn)生的DQB文件是否正確(可以參考說(shuō)明書中 4.3 DQB DAT 文件的內(nèi)容和格式)；(2) 點(diǎn)擊剖分單元按鈕,生成滲流計(jì)算需要的原始數(shù)據(jù)文件G- DAT和剖分圖NETO.BLN;(3) 點(diǎn)擊顯示剖分圖按鈕,進(jìn)展單元剖分同時(shí)生成網(wǎng)格圖NET.BLN；(4) 點(diǎn)擊編輯數(shù)據(jù)按鈕,翻開(kāi)計(jì)算原始數(shù)據(jù)文件 G - DAT,注意：在該文件末尾另起一行 輸入0并保存；(5) 點(diǎn)擊計(jì)算按鈕，生

18、成結(jié)果文件SEEP- OUT和滲流場(chǎng)圖形文件 GS- BLN；(6) 點(diǎn)擊等勢(shì)線圖啟動(dòng) SURFER繪制滲流場(chǎng)分布圖：a、單擊 MAP菜單中的 LOAD BASE MAP出現(xiàn)IMPORT FILE視窗，雙擊 BOUND BLN, 單擊OK隨即出現(xiàn)土壩橫斷面圖；b、單擊 MAP菜單中的LOAD BASE MAPH現(xiàn)IMPORT FILE視窗，雙擊 GS- BLN 單 擊OK隨即出現(xiàn)滲流場(chǎng)圖形；c、將土壩橫斷面圖與滲流場(chǎng)圖形進(jìn)展疊加：用鼠標(biāo)將滲流場(chǎng)圖形拖至堤防橫斷面圖，用鼠標(biāo)同時(shí)框選堤防橫斷面圖和單個(gè)滲流場(chǎng)圖形，單擊MAP出現(xiàn)MAP菜單,單擊MAP菜單中的 OVERLAY MAP即完成圖形疊加；d

19、、進(jìn)展文字編輯、坐標(biāo)軸設(shè)定等，最終得到完美的滲流場(chǎng)分布圖。(7) 假設(shè)想查看計(jì)算結(jié)果 SEEP- OUT那么點(diǎn)擊計(jì)算報(bào)告按鈕。4.3 DQB. DAT文件的內(nèi)容和格式產(chǎn)生的數(shù)據(jù)文件DQB.DAT內(nèi)容格式如下：第一行：SEEP1, SEEP2, SEEP3第二行：ZU, KUAI , DIAN , BIAN , TUCENG , SLDM , SLK , SD , BFS第三行：KXX 1： KUAI , 1 : 9第四行：ZB 1 : DIAN , 1 : 2第五行：SLCS 1： TUCENG , 1 : 3第六行：BJXX 1: BIAN , 1: 3第七行：SLKH 1: SLK第八行：

20、SLDMXX 1: SLDM , 1: 3第九行：SJSW 1: 2, 1 : 2數(shù)據(jù)說(shuō)明：1SEEP 1判別是否保存滲流繪圖數(shù)據(jù)， Y 為肯定， N 為否認(rèn)。SEEP 2判別是否自動(dòng)剖分?jǐn)?shù)據(jù)， Y 為肯定， N 為否認(rèn)。SEEP 3設(shè)定為N 。2ZU分塊組數(shù)計(jì)算區(qū)域內(nèi)剖分的塊依次從左到右分為假設(shè)干組，同一組間應(yīng)為上下 關(guān)系，即第I塊的下邊n3與第1+1塊的上邊n 1應(yīng)為同一邊，每組內(nèi)由上而下順序編 號(hào)；KUAI 四邊形塊的數(shù)目；DIAN 塊角點(diǎn)數(shù)；BIAN 含各類邊界的邊數(shù)；TUCENG 不同滲透系數(shù)的土層種類數(shù)；SLDM 需要計(jì)算滲流量的滲流斷面?zhèn)€數(shù)，如果是零表示不計(jì)算滲流量；SLK 需要

21、計(jì)算滲流量的斷面所穿過(guò)的總計(jì)塊數(shù)，如果是零表示不計(jì)算滲流量；SD 計(jì)算非穩(wěn)定滲流時(shí)水庫(kù)水位降落的時(shí)段數(shù)，如果是零表示計(jì)算穩(wěn)定滲流；BFS 計(jì)算等勢(shì)線的百分?jǐn)?shù)間距值，默認(rèn)為 10。3KXX(1:ME,1:9)塊信息數(shù)組每塊9個(gè)分量，按塊號(hào)從小到大順序填寫。第1 4分量表示塊從左上角起按順時(shí)針?lè)较虻乃膫€(gè)角點(diǎn)點(diǎn)號(hào)；第58分量依次表示塊從上邊起按順時(shí)針?lè)较?條邊的剖分點(diǎn)數(shù)；第9分量為塊的滲透系數(shù)區(qū)號(hào)。4ZB 1： DIAN，1： 2塊的結(jié)點(diǎn)坐標(biāo)數(shù)組按塊結(jié)點(diǎn)序號(hào)順序表示每一結(jié)點(diǎn)的X，Y坐標(biāo)。5SLCS 1 : TUCENG，1： 3依次為每層土的水平滲透系數(shù)Kx單位：cm/s，垂直滲透系數(shù) Ky單位：c

22、m/s，給水度 ，0.11 k該式中滲透系數(shù) k的單位為：m/day。6BJXX(1:BN,1:3) 邊界信息數(shù)組每條邊界有三個(gè)分量，第一分量為邊界邊所在塊的塊號(hào)；第二分量為該邊界邊在塊中的邊序號(hào)；第三分量為邊的邊界類型號(hào)，規(guī)定為：1 滲流自由面邊界斜、心墻壩僅為上部自由面3半有壓流邊界5斜、心墻壩下游壩殼自由面邊界6均質(zhì)土壩下游壩坡或排水體滲出段邊界7上游水頭邊界8下游水頭邊界9斜、心墻土壩斜墻心墻后坡滲出段邊界7SLKH 1 : SLK需計(jì)算滲流量的過(guò)流斷面上的塊號(hào)。8SLDMXX 1 : SLDM，1 : 3滲流斷面的塊單元信息第一分量為某過(guò)流斷面所經(jīng)過(guò)的塊數(shù)；第二分量為該過(guò)流斷面與塊邊

23、界相交的信息，水工滲流計(jì)算一般為 1;第三分量為該過(guò)流斷面在塊中位置的節(jié)點(diǎn)鏈序號(hào)，一般為1或2。9SJSW 1 : 2，1： 2起降時(shí)間天、水位；降后時(shí)間天、水位。4.4 G.DAT數(shù)據(jù)內(nèi)容本軟件可將四邊形塊單元自動(dòng)剖分為三角形單元，同時(shí)由DQB.DAT自動(dòng)產(chǎn)生 G.DAT,G.DAT是本軟件進(jìn)展?jié)B流計(jì)算的原始數(shù)據(jù)文件，其數(shù)據(jù)內(nèi)容格式如下：第一行：SEEP1, SEEP2, SEEP3第二行：NE，NN，NNU，NND，NNM1，NNM2，NNF，TUCENG，AA，AA1，AA2，BB1，BB，BB2，NNP，NQ，SD，NW，NG，NNM3，LU，LD，MA第三行：Hmax, H1，H2，

24、H3，EPS，EPS1，DT，Yaa，Xaao o o o o o第十四行：SW 1 : SD第十五行：SJJG 1： SD第十六行：SLCS 1： TUCENG , 1 : 3第十七行：QTE(1 : NQ)第十八行：QEN(1 : Qe)解釋：1SEEP1, SEEP2, SEEP3同 DQB.DAT2NE 三角形單元總數(shù)o o o o o o3H max 上游最高水位H1 上游水位H2 下游水位H3 井水位EPS 自由面迭代允許誤差 (0.10.001)EPS1 井附近迭代允許誤差 (0.50.005)DT 水位變化總歷時(shí)單位：天O O O O O O O14SW 1： SD各時(shí)段末上游

25、水位。15SJJG 1 : SD各時(shí)段時(shí)間間隔天數(shù)16SLCS 1: TUCENG , 1 : 3依次填寫每層土的水平滲透系數(shù)Kx單位：m/day；垂直滲透系數(shù)Ky單位：m/day;給水度 ，0.1177、k式中k為滲透系數(shù)，單位：m/day。17QTE(1 : NQ)計(jì)算滲流量斷面的單元總數(shù)18QEN(1 : Qe)計(jì)算滲流量斷面的所有單元號(hào)。5算例例一均質(zhì)土壩穩(wěn)定滲流計(jì)算某均質(zhì)土壩橫斷面尺寸、上下游水位及滲透系數(shù)如圖5.1所示，要求計(jì)算穩(wěn)定滲流場(chǎng)解題步驟如下： 運(yùn)行CCED.EXE文件，彈出工程及視圖相關(guān)參數(shù)界面，在“分塊組數(shù)中填入5，在“土層種類數(shù)中填入2,選擇“穩(wěn)定滲流，并選擇“計(jì)算滲

26、流量；單擊視圖相關(guān)參數(shù)，“最大坐標(biāo)中X填寫170, Y填寫100, “最小坐標(biāo)中 X填 寫-30, Y填寫10,其它的按照默認(rèn)值；(3)單擊確定按鈕，進(jìn)入圖形繪制界面，參考“4.1 DQB文件的產(chǎn)生中具體的操作步驟建立均質(zhì)壩橫斷面剖分圖(如圖5.2);圖5.2均質(zhì)壩橫斷面剖分圖 保存dqb.dat文件，內(nèi)容如下：Y y N5 8 16 7 2 1 2 0 101 4 5 2 3 2 3 2 13 6 7 4 2 3 4 4 24 7 8 5 4 2 4 2 16 9 10 7 4 4 4 3 27 10 11 8 4 2 4 2 19 12 13 10 8 4 8 4 210 13 14 11

27、 8 2 8 2 113 15 16 14 4 2 4 2 10.000 45.00042.000 45.00042.000 40.00067.000 65.00067.000 45.00067.000 40.000121.000 48.000118.000 45.000118.000 40.000140.000 45.000140.000 40.0000.000002000 0.000002000 0.0471573350.000010000 0.000010000 0.0593474561 1 72 4 72 1 7 4 1 76 1 16 2 88 1 8點(diǎn)擊圖形繪制界面中的山；按鈕，彈

28、出計(jì)算操作界面:a、b、c、d、點(diǎn)擊查看DQB按鈕，檢查產(chǎn)生的 DQB文件是否正確（可以參考說(shuō)明書中 4.3 DQB.DAT文 件的內(nèi)容和格式）； 點(diǎn)擊剖分單元按鈕,生成滲流計(jì)算需要的原始數(shù)據(jù)文件 點(diǎn)擊顯示剖分圖按鈕,進(jìn)展剖分同時(shí)生成網(wǎng)格圖 點(diǎn)擊纟G- DAT和剖分圖 NET0.BLIN NET.BLN弟輯數(shù)據(jù)按鈕,翻開(kāi)計(jì)算原始數(shù)據(jù)文件 G- DAT注意：在該文件末尾另起一行輸入e、f、點(diǎn)擊等勢(shì)線圖啟動(dòng) 選擇MAF菜單中的 隨即出現(xiàn)土壩橫斷面圖；選擇MAP菜單中的 隨即出現(xiàn)滲流場(chǎng)圖形；將土壩橫斷面圖與滲流場(chǎng)圖形進(jìn)展疊加：用鼠標(biāo)將滲流場(chǎng)圖形拖至堤防橫斷面圖，同時(shí)框選堤防橫斷面圖和單個(gè)滲流場(chǎng)圖形，

29、選擇MAP菜單中的OVERLAY MAP即完成圖形疊加;0并保存； 點(diǎn)擊計(jì)算按鈕，生成結(jié)果文件SEEP- OUT和滲流場(chǎng)圖形文件 GS- BLN;SURFER觀察滲流場(chǎng)分布圖：LOADBASEMAP出現(xiàn) IMPORT FILE 視窗，雙擊 BOUND BLN,單擊 OKLOAD BASE MAP出現(xiàn) IMPORT FILE 視窗，雙擊 GS BLN,單擊 OK圖5.3均質(zhì)壩滲流場(chǎng)分布圖翻開(kāi)計(jì)算報(bào)告SEEP- OUT,得到單位壩長(zhǎng)的日滲漏流量QQ1=0.0313 卅/day/m。例二堤防非穩(wěn)定滲流計(jì)算堤防、堤基橫斷面尺寸及滲透系數(shù)如圖5.4所示，下游水位為31.0m,河水位由38.0m經(jīng)4.5天

30、降至36.0m，又經(jīng)3.5天降至34.0m。要求計(jì)算河水位降落過(guò)程中的堤防滲流場(chǎng)。圖5.4堤防橫斷面解題步驟可以參考例一，建立堤防橫斷面剖分圖如圖5.5圖5.5堤防橫斷面剖分圖 產(chǎn)生本問(wèn)題的 DQB- DAT文件內(nèi)容如下：Y y N3 5 12 7 2 0 0 2 101 5 6 2 4 5 4 5 23 7 8 4 13 3 13 3 14 8 9 5 13 3 15 3 15 9 10 6 15 5 15 5 29 11 12 10 4 5 4 5 20.000 30.0000.000 20.00030.000 38.00020.000 34.00010.000 30.00010.000 20.00053.000 33.5000.000010000 0.000010000 0.0593474560.000002000 0.000002000 0.0471573351 1 73 4 72 4 72 1 12 2 63 2 85 1 80 38.000 8.000 34.000注意：GS.BLN中包含了初始水位 38.0、下降瞬時(shí)水位 36.0、下降瞬時(shí)水位 34.0共3個(gè)水位時(shí)的滲流場(chǎng)繪圖數(shù)據(jù)，故應(yīng)將GS.BLN文件拆分成3個(gè)繪圖數(shù)據(jù)文件 GS1.