弹性应变能密度问题

alonelyboat

请问哪位仁兄用FLAC3D做过有关弹性应变能密度问题，望指教，不胜感激！（最好是能贴出命令流）
这是我根据弹性应变能密度公式编写的命令流，能运行但最终无法显示出云图，望论坛里的高手指教。
def energy1
array energys(10,10,10)
loop i (1,10)
   loop j (1,10)
     loop k (1,10)
      s1=sxx(i,,j,k)*sxx(i,,j,k)
      s2=syy(i,,j,k)*syy(i,,j,k)
      s3=szz(i,,j,k)*szz(i,,j,k)
      s4=2*0.25*(sxx(i,,j,k)*syy(i,,j,k)+syy(i,,j,k)*szz(i,,j,k)+sxx(i,,j,k)*szz(i,,j,k));假定泊松比为0.25
      energys(i,,j,k)=(s1+s2+s3-s4)/(2*1e6);二维弹性能计算公式,假定弹性模量为1e6
      energy3=energys(i,,j,k)
      printout
    end_loop
   end_loop
end_loop
end
energy3

jlinz2008

前面一样，最后面位置。
....
end
energy1
energy3
试试看

csuroylin

可参看 zextra 命令

alonelyboat

首先感谢两位的热心解答，jlinz2008的方法不行，csuroylin能否讲详细一点，我对Fish也就懂点皮毛，看了你说的zextra 命令，还是没明白是怎么一回事，另外能否说明如何讲应变能密度的云图显示出来，再次感谢！

alonelyboat

再次呼唤高手，请热心帮助解答，感谢！

huyonglei2010

定义一个 config zextra1

plot con zextra1