n
config fluid
set fluid off
;===================================
;建立网格模型
gen zone brick &
p0 0 0 0 p1 2 0 0 p2 0 0.5 0 p3 0 0 3 size 3 1 3
gen zone brick &
p0 2 0 0 p1 20 0 0 p2 2 0.5 0 p3 2 0 3 &
size 17 1 3 ratio 1.03 1 1
gen zone brick &
p0 2 0 3 p1 20 0 3 p2 2 0.5 3 p3 12 0 13 &
p4 20 0.5 3 p5 12 0.5 13 p6 20 0 13 &
p7 20 0.5 13 size 17 1 17 ratio 1.03 1 1
def parm
nptabl=4
end
parm
table 1 2 0 3 2 0.5 3 20 0.5 8 20 0 8
def createwaterface
loop k(1,waterpoints-1)
xx1=xtable(1,k)
yy1=miny
zz1=ytable(1,k)
xx2=xx1
yy2=maxy
zz2=zz1
xx3=xtable(1,k+1)
yy3=maxy
zz3=ytable(1,k+1)
xx4=xx3
yy4=miny
zz4=zz3
;ii=out(string(xx1)+' : '+string(zz1)+' : '+string(xx3)+' : '+string(zz3))
if xx1<maxx then
command
WATER TABLE face xx1,yy1,zz1 xx2,yy2,zz2 xx3,yy3,zz3
WATER TABLE face xx1,yy1,zz1 xx3,yy3,zz3 xx4,yy4,zz4
end_command
endif
end_loop
end
createwaterface
model elas
prop bulk 3e9 shear 1e9
ini pp 0 range x -0.1 0.1 y 0 0.5 z 0 3
ini pp 0 range x 0 2 y 0 0.5 z 2.9 3.1
ini pp 0 range x 2 7 y 0 0.5 z 3 8 plane dip -45 dd 90 ori 2 0 3 above
model fl_iso
prop perm 1e-10 poro 0.8
ini fmod 2e9 ftens -1e10 fdens 1000
;====================================
;设置边界条件
fix x y z range z -0.1 0.1
fix x range x 19.9 20.1
fix x range x -0.1 0.1
fix y
hist rat
plot con pp ou on
set fluid imp off
set fluid dt 10
solve
