- 积分
- 0
- 注册时间
- 2007-10-5
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2007-11-2 20:28:04
|
显示全部楼层
来自 福建厦门
能量等值线云图问题出在哪里?
现用一简单FISH例子,作能量等值线输出,命令流可以运算,但最后却不能输出云图,问题肯定出在新云图层定义上,对此不得其解,请达人指点.
config
grid 50,50
model mo
prop density=1500.0 bulk=2E8 shear=1E8 tens 0.5e6 dil 10 fric 30 cohesion 2e6
model null i 20 25 j 20 25
set large
def energy1
array energys(50,50)
loop i (1,50)
loop j (1,50)
s1=sxx(i,j)*sxx(i,j)
s2=syy(i,j)*syy(i,j)
s3=2*0.25*sxx(i,j)*syy(i,j);假定泊松比为0.25
s4=2*(1+0.25)*sxy(i,j)
energys(i,j)=(s1+s2-s3+s4)/(2*1e6);二维弹性能计算公式,假定弹性模量为1e6
energy3=energys(i,j)
printout
end_loop
end_loop
end
fix x i 1
fix x i 51
fix x y j 1
set gravity=10.0
apply syy -20e6 j 51
his 1 energy3
pause
solve
plot hold his 1 |
|