FLAC3D：[2]壳单元

1、newgen zon radc鲻戟缒男yl 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.2r4 的位置上出现了较为明显的梯度变化。

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加上轮廓线