请问哪位仁兄用FLAC3D做过有关弹性应变能密度问题,望指教,不胜感激!(最好是能贴出命令流) ,J~1~fg89
这是我根据弹性应变能密度公式编写的命令流,能运行但最终无法显示出云图,望论坛里的高手指教。 *1%g=vb
def energy1 zCQv:.0L
array energys(10,10,10) stDn{x.
loop i (1,10) /J6CSk
loop j (1,10) S_QDYnF)`
loop k (1,10) MUo?ajbqOd
s1=sxx(i,,j,k)*sxx(i,,j,k) y~d W=zO
s2=syy(i,,j,k)*syy(i,,j,k) l9jcoVo.
s3=szz(i,,j,k)*szz(i,,j,k) O7K))w
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 YW'Y=*
energys(i,,j,k)=(s1+s2+s3-s4)/(2*1e6);二维弹性能计算公式,假定弹性模量为1e6 c49#aNR
energy3=energys(i,,j,k) /#
eBDo
printout 'Oc8[8
end_loop NMDNls&)k
end_loop !kIw835U
end_loop '|cuVxcE55
end i3~!ofTb
energy3