以下是我看著書打的指令 但在第二行中的 ats出現錯誤 希望有人可以幫我解答
; rocksl. 岩石邊坡穩定性分析;
config ats
grid 70,60
gen -40.0,0.0 -40.0,00,0.00,-40.0 i 1 21 j 1 21
gen -40.0,0.0 -40.0,80.0 0.0,000 i 1 21 j 21 61
gen 0.0,-40.0 0.0,0.0 100.0,0.0 100.0,-40.0 i 21 71 j 1 21
gen 0.0,0.0 0.0,80.0 100.0,0.0 i 21 71 j 21 61
model elastic i 1 70 j 1 60
taba 1 -40.00 0 0 0 3.880 15.00 11.88 15.00 15.76 23.76 &
30.00 27.64 45.00 43.15 60.50 80.00 712-00 100 80.00
gen table 1
model null region 34.48
group'rock'notnull
prop density = 2500.0 bulk=2.0E9 shear = 1.0E9 cohesion = 140000.0&
friction = 45.0 dilation=0.0 tension =120000.0 group 'rock'
fixx i 71 j 1 61
fixi 1 j 1 21
set gravity = 9.81
history 999 unbalanced
solve elastic
save rocsl 1.sav '初始平衡,乾狀態'
table 2 -40.00 0 0 0 3.880 4.500 11.88 12.00 15.76 112 - 00 23.76&
24.00 27.64 27.00 43.15 41.00 52.00 43.00 80.00 45.00 100.0 412- 00
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 (j,j) = 2800
endif
endif
endloop
endloop
end
wet den
water density = 1000.0
water table = 1
solve
sav rocs_2.sve '初始平衡狀態,地下水位'
solve fos file rocskl_fos.sav
