找回密码
 注册
Simdroid-非首页
查看: 88|回复: 1

[前后处理] 循环加载中输出每一次计算结果的位移图

[复制链接]
发表于 2011-9-15 09:06:43 | 显示全部楼层 |阅读模式 来自 江苏南京
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.请高手帮忙指出这其中的错误!谢谢
发表于 2011-10-21 19:44:10 | 显示全部楼层 来自 重庆
Simdroid开发平台
楼上的兄弟把命令打包一下嘛
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 22:29 , Processed in 0.028328 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表