- 积分
- 9
- 注册时间
- 2005-10-30
- 仿真币
-
- 最后登录
- 1970-1-1
|
我编的一个求MC模型单元安全度的FISH程序,欢迎大家调试,并提出意见和建议!!
gen zone brick &
p0 0 0 0 p1 2 0 0 p2 0 1 0 p3 0 0 2 &
size 2 1 2
config zextra 1
def saf
p_z=zone_head
loop while p_z # null
s1=z_sig1(p_z)
s3=z_sig3(p_z)
ALFA1=z_prop(p_z,'fri')
c1=z_prop(p_z,'coh')
shang=s1-s3
xia=2* c1*cos(alfa1*pi/180)-(s1+s3)*sin(alfa1*pi/180)
f1=shang/xia
z_extra(p_z,1) = f1
pnt=z_next(p_z)
endloop
end
odel mohr
prop density = 2000.0 bulk = 1.0E8 shear = 3.0E7 coh 12380 tens 0 &
friction = 20 dilation = 20
fix x y z range z -0.1 0.1
fix y range y -0.1 0.1
fix y range y 0.9 1.1
set gravity 0,0, -10.0
saf
pl con zextra 1
solve
[ 本帖最后由 ssw1216 于 2007-10-7 00:34 编辑 ] |
|