求助flac3d中怎么用FISH编写应变能代码 [复制链接]

求助：FLAC3D中怎么用fish编写应变能代码

什么叫做应变能？
可以定义一个单元上的额外变量，应变的话可以用ssi应变增量和ssr应变速率，遍历单元进行提取，参与函数计算。当然计算完了也可以输出结果

就是能量，一定意义相当于做功，没有找到合适的代码和公式进行计算、、方便加下qq便于以后交流吗，非常感谢，我的qq1316833994

eyrftjhutyjtjurujr

Wiki了一下，应变能的定义为应变能密度对体积的积分，应变能密度在各向同性弹性条件下表达式为：
                    U = 1/2*(sig_x*e_x + sig_y*e_y + sig_z*e_z + tau_xy*gama_xy + tau_yz*gama_yz + tau_zx*gama_zx)
如果用主应力表示为：
                    U = 1/2(sig1*e1 + sig2*e2 + sig3*e3)
弹性状态下，应力应变一一对应由胡克定律可以把上式改写成：
                    U = ((1-2v)/6E)*(sig1 + sig2 +sig3)^2 + ((sig1-sig2)^2 + (sig2-sig3)^2 + (sig1-sig3)^2)/12G
                    E 杨氏模量 v 泊松比 G 剪切模量
上式第一部分为体积应变能，第二部分为偏应变能。
用FISH的话可以设置一个z_extra变量，在计算过程中可以实时显示，代码基于FLAC3D 4.0 或5.0
Def Stain_Energy_elast
    E = ....
    v = .....
    G = E/(2*(1+v))
    Sum_U_elast = 0.
    p_z = zone_head
       loop while zone # null
             Uev  = ((1-2v)/6E)*(z_sig1(p_z) + z_sig2(p_z) + z_sig3(p_z))^2
             Ued = ((z_sig1(p_z) - z_sig2(p_z))^2 + (z_sig2(p_z) - z_sig3(p_z))^2 + (z_sig3(p_z) - z_sig1(p_z))^2)/12/G
             z_extra(p_z,1) = (Uev + Ued)*z_volume(p_z)
             Sum_U_elast = Sum_U_elast + z_extra(p_z,1)
          p_z = z_next(p_z)
       endloop
end
@Strain_Energy_elast
hist @Strain_Energy_elast
以上只是考虑到弹性条件下的例子，仅供参考，塑性下需要考虑能量耗散，建议自己查阅文献。
当然如果你用FLAC3D 5.0算得话就不用这么麻烦了，5.0自带了内变量，弹性： z_wevol(ptr), z_weshear(ptr),z_wetot(ptr)
塑性： z_wevol(ptr), z_weshear(ptr),z_wetot(ptr)， 直接调用ok啦

紫色毛笔:Wiki了一下，应变能的定义为应变能密度对体积的积分，应变能密度在各向同性弹性条件下表达式为：
                    U = 1/2*(sig_x*e_x + sig_y*e_y + sig_z*e_z + tau_ .. (2016-08-26 00:08) 

牛逼啊，虽然没用，看起来就不错啊

好东西，非常感谢啊，

        

能不能教教5.0版本怎么写考虑能量耗散的弹性应变能代码哇，不会写，谢谢您