2維土石壩滲流計算

1、土石壩滲流分析，采用非飽和土滲流參數(shù)，迭代計算浸潤線，根據(jù)前次計算結(jié)果，不斷修改單元的滲透系數(shù)和浸潤逸出點位置，直到滿足精度要求。本算例的土石壩體型比較簡單.采用非飽和滲流計算.即滲透系數(shù)為空隙壓力的函數(shù).首先建立一個數(shù)據(jù)文件PPPP.TXT，存儲滲透系數(shù)函數(shù)關(guān)系，如下。第一列為空隙壓力值（水頭M)，第二列為滲透系數(shù)指數(shù)，滲透系數(shù)等于10^A(MD)。10.004.0E009.003.6E008.003.2E007.002.8E006.

2、002.4E005.002.0E004.001.6E003.001.2E002.008.0E011.004.0E010.000.0E00土壩頂寬4M，上下游坡比均為1:2，總高12M，底寬52M。上游水深8M，下游無水。FINITITLEEARTHDAMSEEPAGEFILNAMESEEPAGE5PLOPTSDATE0DIMTPRETABLE1111PRESSKKPE!定義水壓與滲透系數(shù)的關(guān)系數(shù)組TREADTPREPPPPTXT!讀入數(shù)

3、組DIMNCONARRAY4!定義數(shù)組，用于存貯單元四個節(jié)點號PREP7SMRTOFFANTYPESTATIC!THERMALANALYSISET1PLANE55MPKXX11!飽和狀態(tài)下的滲透系數(shù)MPKXX21E4!完全干燥下的滲透系數(shù)，假設(shè)空隙水壓力小于10M時K12412K2240K300K42812K5280K6520L13L32L12L45L56L46LESIZEALL24ENDDOALLSELALLFINISHSOLUSOL

4、VEFINISH!第一次計算完畢!!迭代計算CONUTT=20!最大循環(huán)次數(shù)DD_HEAT=0.001!前后兩次計算，總水頭最大允許計算差CHUK_ST=3!出口邊界條件重新設(shè)定的起始點CHUK_MAXY2=10E5!臨時變量，用于存儲浸潤線出口坐標DOCOM_NUM1CONUTTDD_H=0POST1SET1DOI1N_MAXIFCOM_NUMGTCHUK_ST1THENDD1=N_TEMP(I)IFABS(DD1TEMP(I))GT

5、DD_HTHENDD_H=ABS(DD1TEMP(I))ENDIFENDIFN_TEMP(I)=TEMP(I)!計算節(jié)點溫度（總水頭）N_PRE(I)=N_TEMP(I)NY(I)!計算節(jié)點壓力，總水頭Y坐標ENDDOIFCOM_NUMGTCHUK_ST1THENIFDD_HLEDD_HEATTHENEXITENDIFENDIFPREP7!重新給每個單元設(shè)定材料MATNUM=2DOI1E_MAXDOKK14GETNCON(KK)ELEM