- 积分
- 0
- 注册时间
- 2009-2-14
- 仿真币
-
- 最后登录
- 1970-1-1
|
gen zon bri size 3 3 3
model mohr
prop bu 3e6 sh 1e6 coh 10e3 fric 15
fix z ran z -.1 .1
fix x ran x -.1 .1
fix x ran x 2.9 3.1
fix y ran y -.1 .1
fix y ran y 2.9 3.1
ini dens 2000
hist unbal
set grav 10
solve elastic
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
def n_load
pnt=gp_near(1,1,3)
loop n(1,10)
load=-1e3*n
pict='zdisp_'+string(n)
command
app nstress load ran z 2.9 3.1 x 1 2 y 1 2
solve rat 1e-6
endcommand
xtable(1,n)=-1e3*n
ytable(1,n)=gp_zdisp(pnt)
command
set plot jpg
pl create pict
pl set back white
pl set title text
'Z-Displacement VS load'+string(n)
pl add table 1 mark line xlabel 'load' ylabel 'z-displacement'
set out pict.jpg
pl hardcopy pict
endcommand
endloop
end
table 1 name 'zdisp'
n_load
实际做实验的时候,总是遇到循环加载的过程。为了得到每一次加载后模型顶部的位移图,本人考虑写了上面的一小段程序,逻辑上我看没有问题的啊。程序运行的时候,会先出现bad integer conversion.请高手帮忙指出这其中的错误!谢谢
|
|