FLAC3D：[2]壳单元

1、newgen zon radcyl p0 0 0 0 p1 7.5 0 0 p2 0 36 0 p3 0 0 5 size 4 24 8 4 dim 3.2 3.2 3.2 3.2 ratio 1 1 1 1.2 gen zone reflect normal（-1 0 0）origin (0 0 0)gen zone reflect normal（0 0 -1) origin (0 0 0)sel shell id=12 ran cylinder end1=(0 0 0) end2=(0,36,0) radius=3.3plo sel geomplo add sketch

------------------------------------------------

最后的结果见下图。

2、先来了解一下radcylinder命令。见下图。

gen zon radcyl p0 0 0 0 p1 7.5 0 0 p2 0 36 0 p3 0 0 5 size 4 24 8 4 dim 3.2 3.2 3.2 3.2 ratio 1 1 1 1.2

将这个语句拆分，它是先生成了外面的整体长方体（通过P0，P1，P2，P3的确定）。 p0 7.5 0 0 p1 15 0 0 p2 7.5 36 0 p3 7.5 0 5

继而，定义了在每个方向上的网格数量：

size 4 24 8 4

只有大的长方体，也有网格数，行不行？

还得有具体的尺寸限制，于是出现了下面的语句

dim 3.2 3.2 3.2 3.2 

它指定的是内部的圆柱体的四个边长

如果想要继续进行限制，要实现比率的变化，进行

rat 1 1 1 1.2

r4 的位置上出现了较为明显的梯度变化。

3、sel shell id=12 ran cylinder end1=(0 0 0) end2=(0,36,0) radius=3.3

因为是圆形隧道，只需要用cylinder命令，把周圈围岩的中心包括在内就可以了。

4、plo sel geom

显示壳单元

5、plo add sketch

加上轮廓线