论坛 › G0201：FLAC3D/FLAC2D › 如何计算塑性区深度

cdf5041353 发表于 2014-7-15 15:15:08

如何计算塑性区深度

请教各位，怎么计算围岩的塑性区深度啊，比如顶拱中心点处塑性区深度如何Fish编程实现？

陌上公子 发表于 2014-7-16 11:26:27

同问啊，还有各个区域的塑性区面积怎么统计呢（顶、底、两帮）

cdf5041353 发表于 2014-7-18 11:10:42

陌上公子 发表于 2014-7-16 11:26
同问啊，还有各个区域的塑性区面积怎么统计呢（顶、底、两帮）

我知道一个计算塑性区体积的命令流，至于计算面积还没有查。另外，无论是面还是体积，我觉得都可以在顶、底、两帮进行不同分组，然后分别调用命令文件输出应该可以实现。
;;;;;;;;;;;;;;;监测塑性区体积
def suxingqu
    shearnow=1
    tensionnow=2
    shearpast=4
    tensionpast=8
    v_shear_now=0
    v_tension_now=0
    v_shear_past=0
    v_tension_past=0

    p_z=zone_head
    loop while p_z # null
   if and(z_state(p_z,0),shearnow)=shearnow then
      v_shear_now=v_tension_now+z_volume(p_z)
   endif
   if and(z_state(p_z,0)),shearpast)=shearpast then
      v_shear_past=v_shear_past+z_volume(p_z)
   endif
   if and(z_state(p_z,0),tensionpast)=tensionpast then
      v_tension_past=v_tension_past=z_volume(p_z)
      endif
   p_z=z_next(p_z)
      endloop
      ii=out('剪切破坏：’+string(v_shear_now))
      ii=out('拉伸破坏：’+string(v_tension_now))
      ii=out('过去曾剪切破坏：’+string(v_shear_past))
      ii=out('过去曾拉伸破坏：’+string(v_tension_past))
end

琅琊逝水寒 发表于 2015-1-8 19:38:59

我编过用fish求洞室塑性区深度的命令，但程序有很大的限制，只能是一个洞室，跟网格大小有关。

龙biscuit 发表于 2016-3-14 05:05:33

cdf5041353 发表于 2014-7-18 11:10
我知道一个计算塑性区体积的命令流，至于计算面积还没有查。另外，无论是面还是体积，我觉得都可以在顶、 ...

学习了，谢谢分享:lol

牛掰的网络世界 发表于 2016-6-1 15:39:11

好深奥。。。

牛掰的网络世界 发表于 2016-6-13 10:14:36

谢谢分享

查看完整版本: 如何计算塑性区深度