- 积分
- 0
- 注册时间
- 2004-11-1
- 仿真币
-
- 最后登录
- 1970-1-1
|
new
def series
new_fric = 85.0
step_lim=1000
ipt = gp_near(1,5,1)
loop n (1,8)
command
prop fric = new_fric
ini sxx 0.0 syy 0.0 szz 0.0 sxy 0.0 sxz 0.0 syz 0.0
ini xvel 0.0 yvel 0.0 zvel 0.0
ini xdis 0.0 ydis 0.0 zdis 0.0
print new_fric
set mech step step_lim force 50
solve step step_lim force 50
endcommand
xtable(1,n) = new_fric
ytable(1,n) = log(abs(gp_ydisp(ipt)))
new_fric = new_fric-inc_fric
endloop
end
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
gen zone brick size 3 5 3
model mohr
prop bulk 1e8 shear 0.3e8 coh 1e4 ten 1e4
fix y range y -0.1 0.1
set grav 0 -10 0
ini dens 2500
set inc_fric = 5
series
plot table 1 both |
|