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

请问哪位仁兄用FLAC3D做过有关弹性应变能密度问题，望指教，不胜感激！（最好是能贴出命令流） \e=Iw"yd  
这是我根据弹性应变能密度公式编写的命令流，能运行但最终无法显示出云图，望论坛里的高手指教。 Nx%]dOa  
def energy1 |qn`z-  
array energys(10,10,10) DR d|m
<Z  
loop i (1,10) Y`O"+Jr  
   loop j (1,10) i
r72fSe  
     loop k (1,10) ')bas#=uP  
      s1=sxx(i,,j,k)*sxx(i,,j,k) ^-*q  
      s2=syy(i,,j,k)*syy(i,,j,k) ]YO &_#  
      s3=szz(i,,j,k)*szz(i,,j,k) go6
XUe  
      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 $,mljJSQv  
      energys(i,,j,k)=(s1+s2+s3-s4)/(2*1e6);二维弹性能计算公式,假定弹性模量为1e6 :jJ;&t^^  
      energy3=energys(i,,j,k) xla64Qld  
      printout h~!KNF*XW  
    end_loop S\,~6]^T  
   end_loop 0TDcQ  
end_loop z]^u@]@NC  
end gB&8TE~Y  
energy3

楼上高人~~~~

res elasticroadway.sav Y`U[Y Hx  
config zextra 1 _#SCjFz  
def elasenergy PQ#zF&gL9t  
array ai(6) RRasX;
zK  
wp_z=zone_head h&k*i  
loop while p_z # null 3YFbT Z  
dum=z_fsi(p_z,ai) Im<i.a <`  
sxx=abs(z_sxx(p_z)) iE}Lw&x  
syy=abs(z_syy(p_z))   S5,y!K]C~  
szz=abs(z_szz(p_z)) \i2S'AblYq  
sxy=abs(z_sxy(p_z)) 'a^'f]"  
syz=abs(z_syz(p_z)) 5C
G ,l  
sxz=abs(z_sxy(p_z)) _wY<8 F*  
exx=abs(ai(1)) $YPU(y  
eyy=abs(ai(2)) GY3g`M   
ezz=abs(ai(3)) Jn+k$'6%#  
exy=abs(ai(4)) 3jJd)C R  
eyz=abs(ai(5)) ^wtr~D|  
exz=abs(ai(6)) ;ATn&  
ee=0.5*(sxx*exx+syy*eyy+szz*ezz+sxy*exy+syz*eyz+sxz*exz) 4)DI0b"  
z_extra(p_z,1)=ee` T\ cJn>kCn  
p_z=z_next(p_z) b}C6/zW  
endloop : U Yn  
end  g%.;ZlK  
elasenergy 6]mFw{6qn1  
save newelasenergy.sav

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

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

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

太厉害了，不过看不懂啊

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

请高手热心解答，感谢！