[软件应用交流]弹性应变能密度问题 [复制链接]

请问哪位仁兄用FLAC3D做过有关弹性应变能密度问题，望指教，不胜感激！（最好是能贴出命令流） >=O5=\`  
这是我根据弹性应变能密度公式编写的命令流，能运行但最终无法显示出云图，望论坛里的高手指教。 Smg z}  
def energy1 1{15#W  
array energys(10,10,10) &uP,w#  
loop i (1,10) 8GlRO4yd  
   loop j (1,10) q[}W&t,  
     loop k (1,10) r|av|7R  
      s1=sxx(i,,j,k)*sxx(i,,j,k) y8%QS*  
      s2=syy(i,,j,k)*syy(i,,j,k) a1#",%{I  
      s3=szz(i,,j,k)*szz(i,,j,k) yVfF *nG  
      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 Xnc?o
T+  
      energys(i,,j,k)=(s1+s2+s3-s4)/(2*1e6);二维弹性能计算公式,假定弹性模量为1e6 FIu|eW+<l  
      energy3=energys(i,,j,k) <*_DC)&79  
      printout L+K,Y:D!W  
    end_loop  \uG^w(*)  
   end_loop N.'-9hv  
end_loop Ky[s&>02  
end (! a;}V<7  
energy3

请高手热心解答，感谢！

有谁做个这方面的研究或者有所了解的，请帮忙解答一下，急啊，万分感谢！

太厉害了，不过看不懂啊

就是不厉害啊，要是厉害早就自己搞定了

楼主现在有没有搞定啊，可不可以共享一下成果啊？

你得定义一个config gpextra或者config zextra 才能显示图形

res elasticroadway.sav Fv3:J~Yf  
config zextra 1 "v+%F  
def elasenergy +IM6 GeH  
array ai(6) eR!K8W  
wp_z=zone_head 7jPn6uz>w  
loop while p_z # null ~2}^ -,  
dum=z_fsi(p_z,ai) jO5Wemqf  
sxx=abs(z_sxx(p_z)) e ls&_BPE  
syy=abs(z_syy(p_z))   VMW<?V 2Z  
szz=abs(z_szz(p_z)) u4~(0  
sxy=abs(z_sxy(p_z)) t_iZ\_8
  
syz=abs(z_syz(p_z)) W C3b_ia  
sxz=abs(z_sxy(p_z)) An=Q`Uxt/  
exx=abs(ai(1)) s/OXZ<C|  
eyy=abs(ai(2)) GI/4<J\  
ezz=abs(ai(3)) [KsVI.gn  
exy=abs(ai(4)) EK}QjY[i  
eyz=abs(ai(5)) p/?TU  
exz=abs(ai(6)) .GuZV'  
ee=0.5*(sxx*exx+syy*eyy+szz*ezz+sxy*exy+syz*eyz+sxz*exz) UpqDGd7M  
z_extra(p_z,1)=ee` k6??+b:rE  
p_z=z_next(p_z) (^= Hq'D  
endloop Ts3(,
Y  
end ^jg{MTa  
elasenergy 7.(vog"I)  
save newelasenergy.sav

楼上高人~~~~