自由液面的土石坝平面渗流分析

用 Ansys 进行土石坝渗流分析，主要难点在于渗流的浸润线未知，需要进行跌代计算。 本算例的土石坝体型比较简单，采用非饱和渗流计算，即渗透系数为空隙压力的函数，渗透系数函数关系，如下，第一列为空隙压力值（水头 m)，第二列为渗透系数参数，渗透系数等于10^a(m/d)。

-10.00 -4.0E+00 
-9.00 -3.6E+00 
-8.00 -3.2E+00 
-7.00 -2.8E+00 
-6.00 -2.4E+00 
-5.00 -2.0E+00 
-4.00 -1.6E+00 
-3.00 -1.2E+00 
-2.00 -8.0E-01 
-1.00 -4.0E-01 
0.00 0.0E+00 
APDL 如下： 
FINI 
/TITLE, EARTHDAM SEEPAGE 
/PLOPTS,DATE,0 
*DIM,TPRE,TABLE,11,1,1,PRESS,KKPE ! 定义水压与渗透系数的关系 
TPRE(1)=-4.0E+00,-3.6E+00,-3.2E+00,-2.8E+00,-2.4E+00,-2.0E+00,-1.6E+00,-1.2E+00,-8. 
0