- 积分
- 0
- 注册时间
- 2006-12-1
- 仿真币
-
- 最后登录
- 1970-1-1
|
我做了一个简单是算例,模型有2种材料(2个group),想达到这个目的:如果是group 1则密度为2700,如果是group 2则在水面上密度为2500,水面下为2600。
命令流如下:
gen zone brick size 4 4 4 p0 0 0 0 p1 20 0 0 p2 0 20 0 p3 0 0 20
gen zone brick size 4 4 4 p0 20 0 0 p1 40 0 0 p2 20 20 0 p3 20 0 20
group 1
gen zone brick size 4 4 4 p0 20 0 20 p1 40 0 20 p2 20 20 20 p3 30 0 40 &
p4 40 20 20 p5 30 20 40 p6 40 0 40 p7 40 20 40
group 2 range group 1 not
m m
pro bulk 2e8 shear 1e8 fric 45 coh 1e5 ten 1e5
set grav 0 0 -9.81
water dens 1000
;形成水位面
water table face 20,0,20 20,20,20 25,0,28
water table face 25,0,28 20,20,20 25,20,28
water table face 25,0,28 25,20,28 35,0,38
water table face 35,0,38 25,20,28 35,20,38
water table face 35,0,38 35,20,38 40,0,40
water table face 40,0,40 35,20,38 40,20,40
water table face 0,0,20 0,20,20 20,0,20
water table face 20,0,20 0,20,20 20,20,20
def ini_dens
pnt=zone_head
loop while pnt#null
if z_group(pnt)=1 then
z_density(pnt)=2700
else
if z_pp(pnt)#0 then
z_density(pnt)=2600
else
z_density(pnt)=2500
endif
endif
pnt=z_next(pnt)
endloop
end
ini_dens
模型请见附件,本人对FISH不是很懂,请问上面的函数哪里出了问题。 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|