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

请问哪位仁兄用FLAC3D做过有关弹性应变能密度问题，望指教，不胜感激！（最好是能贴出命令流） `-L?x2)U  
这是我根据弹性应变能密度公式编写的命令流，能运行但最终无法显示出云图，望论坛里的高手指教。 .zO2g8(VR  
def energy1 d+JK")$9C  
array energys(10,10,10) 2!/Kt O)i^  
loop i (1,10)
N6y9'LGG`  
   loop j (1,10) Lj1l]OD  
     loop k (1,10) 1VX3pkUET  
      s1=sxx(i,,j,k)*sxx(i,,j,k) xPm. TPj  
      s2=syy(i,,j,k)*syy(i,,j,k) uDDa>Ka#+  
      s3=szz(i,,j,k)*szz(i,,j,k) O^@8Drg
c  
      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 m@2E ~m  
      energys(i,,j,k)=(s1+s2+s3-s4)/(2*1e6);二维弹性能计算公式,假定弹性模量为1e6 "Lbsq\W>  
      energy3=energys(i,,j,k) qk(Eyp  
      printout c|Ivet>3  
    end_loop >H;i#!9,  
   end_loop bBjr hi  
end_loop Yr9'2.%Q  
end y7,fFUKl  
energy3

请高手热心解答，感谢！

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

太厉害了，不过看不懂啊

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

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

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

res elasticroadway.sav uP9b^LEoN  
config zextra 1 c)EYXo  
def elasenergy 49FP&NgK  
array ai(6) $WYt`U;*lj  
wp_z=zone_head a~KtH;7<  
loop while p_z # null +w0Wg.4V  
dum=z_fsi(p_z,ai) 5<Kt"5Z%7  
sxx=abs(z_sxx(p_z)) U]Q5};FK  
syy=abs(z_syy(p_z))   HNkOPz+d&8  
szz=abs(z_szz(p_z)) ruaZ(R[  
sxy=abs(z_sxy(p_z)) >MYxj}I4{z  
syz=abs(z_syz(p_z)) ~<1s[Hu  
sxz=abs(z_sxy(p_z)) wBt7S!>G  
exx=abs(ai(1)) ! fk W;|  
eyy=abs(ai(2)) 4~ x>]  
ezz=abs(ai(3)) c!@|yE,  
exy=abs(ai(4)) k}LIMkEa4a  
eyz=abs(ai(5)) /
)K')  
exz=abs(ai(6)) *zRig|k!H  
ee=0.5*(sxx*exx+syy*eyy+szz*ezz+sxy*exy+syz*eyz+sxz*exz) v&:[?<6-  
z_extra(p_z,1)=ee` p5]W2i.,  
p_z=z_next(p_z) nOC\ =<Nsg  
endloop a"s2N%{  
end Od.@G~  
elasenergy w5Fk#zJv  
save newelasenergy.sav

楼上高人~~~~