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

请问哪位仁兄用FLAC3D做过有关弹性应变能密度问题，望指教，不胜感激！（最好是能贴出命令流） Q5hb0O%a  
这是我根据弹性应变能密度公式编写的命令流，能运行但最终无法显示出云图，望论坛里的高手指教。 u>81dO]H  
def energy1 L#`7FaM?  
array energys(10,10,10) ZU)BJ!L,s  
loop i (1,10) g;(_Y1YQ  
   loop j (1,10) "$]ls9-%n  
     loop k (1,10) &3WkH W   
      s1=sxx(i,,j,k)*sxx(i,,j,k) 2ve lH;  
      s2=syy(i,,j,k)*syy(i,,j,k)
a5X`j
o  
      s3=szz(i,,j,k)*szz(i,,j,k) lfXH7jL2~  
      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 c]e`m6  
      energys(i,,j,k)=(s1+s2+s3-s4)/(2*1e6);二维弹性能计算公式,假定弹性模量为1e6 KlY,NSlQ  
      energy3=energys(i,,j,k) fE
'-.nA+  
      printout Z+r%_|kZ  
    end_loop *Yj~]E0`1  
   end_loop 1% asx'^  
end_loop Qk+=znJ  
end %)BwE  
energy3

请高手热心解答，感谢！

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

太厉害了，不过看不懂啊

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

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

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

res elasticroadway.sav 5CY@R  
config zextra 1 c:?#zX  
def elasenergy c:h.J4mv  
array ai(6) h9Tf@]W   
wp_z=zone_head CfkNy[}=  
loop while p_z # null jM3Y|}+  
dum=z_fsi(p_z,ai) F$6])F  
sxx=abs(z_sxx(p_z)) O, ``\(P  
syy=abs(z_syy(p_z))   !\&7oAs=I  
szz=abs(z_szz(p_z)) j$|Yd=  
sxy=abs(z_sxy(p_z)) <iXS0k  
syz=abs(z_syz(p_z)) ? !oVf>  
sxz=abs(z_sxy(p_z)) D<_,>{$gW  
exx=abs(ai(1)) fncwe ';?  
eyy=abs(ai(2)) PNxV
W  
ezz=abs(ai(3)) }> !"SU:d  
exy=abs(ai(4)) Kbz7  
eyz=abs(ai(5)) o/x5  
exz=abs(ai(6))
j^aQ>(t(9  
ee=0.5*(sxx*exx+syy*eyy+szz*ezz+sxy*exy+syz*eyz+sxz*exz) 5|6z1{g8  
z_extra(p_z,1)=ee` pE(<XD3Q  
p_z=z_next(p_z) '&pf  
endloop *&hXJJ[+  
end {0)WS}&  
elasenergy B(en5|  
save newelasenergy.sav

楼上高人~~~~