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

请问哪位仁兄用FLAC3D做过有关弹性应变能密度问题，望指教，不胜感激！（最好是能贴出命令流）
us0{y7(p  
这是我根据弹性应变能密度公式编写的命令流，能运行但最终无法显示出云图，望论坛里的高手指教。 :=J^"c  
def energy1 a(a2xa  
array energys(10,10,10) Jj:Bi&C  
loop i (1,10) eM1=r:jgE  
   loop j (1,10) @_L:W1[  
     loop k (1,10) l)m]<EX  
      s1=sxx(i,,j,k)*sxx(i,,j,k) iem@K  
      s2=syy(i,,j,k)*syy(i,,j,k) t V:oBT*  
      s3=szz(i,,j,k)*szz(i,,j,k) xyvG+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 ^hN.FIzM  
      energys(i,,j,k)=(s1+s2+s3-s4)/(2*1e6);二维弹性能计算公式,假定弹性模量为1e6 E;>BcP
t5  
      energy3=energys(i,,j,k) NF0=t}e  
      printout dZ"B6L!^(  
    end_loop itpljh  
   end_loop M< 1rQW'  
end_loop !idQ-&  
end TA"4yri=7x  
energy3

请高手热心解答，感谢！

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

太厉害了，不过看不懂啊

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

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

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

res elasticroadway.sav Z~.3)6,z  
config zextra 1 U =()T}b>  
def elasenergy {FavF 9O  
array ai(6) W .c:Pulg  
wp_z=zone_head *+E9@r=HF  
loop while p_z # null fRTQ5V  
dum=z_fsi(p_z,ai) TT(dCHft  
sxx=abs(z_sxx(p_z)) =a_
B'^`L  
syy=abs(z_syy(p_z))   }tIIA"dZ  
szz=abs(z_szz(p_z)) 0 w"&9+kV  
sxy=abs(z_sxy(p_z)) ?1I0VA']  
syz=abs(z_syz(p_z)) (> v1)*r  
sxz=abs(z_sxy(p_z)) iG[an*#X  
exx=abs(ai(1)) hb>,\46}  
eyy=abs(ai(2)) wHQyMq^  
ezz=abs(ai(3)) GlDl0P,*r  
exy=abs(ai(4)) X3KPN  
eyz=abs(ai(5)) Dl3Df u8  
exz=abs(ai(6)) >Fzu]G4]  
ee=0.5*(sxx*exx+syy*eyy+szz*ezz+sxy*exy+syz*eyz+sxz*exz) !/w<F{cl  
z_extra(p_z,1)=ee` Q6cF<L`bW  
p_z=z_next(p_z) (wNL,<%~  
endloop &&}'  
end |WiK*  
elasenergy r+BPz%wM=O  
save newelasenergy.sav

楼上高人~~~~