config ats
grid 70,60
gen -40,-40 -40,0 0,0 0,-40 i 1 21 j 1 21
gen -40,0 -40,80 0,80 0,0 i 1 21 j 21 61
gen 0,-40 0,0 100,0 100,-40 i 21 71 j 1 21
gen 0,0 0,80 100,80 100,0 i 21 71 j 21 61
model elastic i 1 70 j 1 60
table 1 -40 0 0 0 4 15 12 15 16 30 24 30 28 45 43 60 80 72 100 80
gen table 1
model null region 34 48
group 'rock' notnull
model mohr group 'rock'
prop density=2500 bulk=2e9 shear=1e9 cohesion=140000 &
friction=45 dilation=0 tension=12000 group 'rock'
fix x y j 1
fix x 71 j 1 61
fix x 1 j 1 21
set gravity=9.81
history 999 unbalanced
solve elastic
save rocs 1.sav
***************************************
table 2 -40 0 0 0 4 5 12 12 16 12 24 24 28 27 43 41 52 43 80 45 100 42
def wet den
loop i (1,izones)
loop j (1,jzones)
if model(i,j)>1 then
**********************下面的一行运行不了！求指教！！！！？
xa=(x(i,j)+x(i+1,j)+x(i+1,j+1)+x(i,j+1))xc=0.25*xa
ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1))
yc=0.25*ya
if yc<table(2,xc) then
density(i,j)=2600
endif
endif
endloop
endloop
end
wet end
water density=1000
water table=1
solve
save rocs2.sav
solve fos file rock1fos.sav
