n
gen zone brick size 25 25 10 p0 0 0 0 p1 2.5 0 0 p2 0 2.5 0 p3 0 0 1
group Base
gen zone brick size 10 10 8 p0 1 1 1.1 p1 add 0.5 0 0 p2 add 0 0.5 0 p3 add 0 0 0.4
group Top range group Base not
range name= Base group Base
range name= Top group Top
model mohr range Base
model elas range Top
interface 1 face range plane norm 0 0 -1 origin 1 1 -1.1 dist 0.1 x (1,1.5) y (1,1.5)
ini z add -0.1 range Top
pro shear 1.19e9 bulk 3.67e9 friction 41 cohesion 7.5e5 ten 2e6 range Base
pro shear 1.06e10 bulk 1.42e10
ini den 2730 range Base
ini den 2400 range Top
interface 1 property kn 5.67e11 ks 6e11 c 0.8e6 fric 41
fix x range x -0.001 0.001
fix x range x 2.499 2.501
fix y range y -0.001 0.001
fix y range y 2.499 2.501
fix z range z -0.001 0.001
apply szz -4e5 range z 1.39 1.41 x 1 1.5 y 1 1.5
apply sxx 2e5 range x -0.01 0.01 y 0 2.5 z 0 1.0
hist unbal
hist gp zdisp (1.25 ,1.25,1.4)
hist int 1 sdis
hist int 1 sstress
set large
set grav 0 0 -9.81
step 5000
大神们，帮我看看错在哪里啊，调了一下午都没找到啊！！都是泪啊！！！
哎
