本节主要介绍如何在圆形盾构隧道中建立壳单元。示例中，先创建几何模型，后加入壳单元。建议读者把第一步中的命令流自行拷贝到 txt 文件中，用file--call命令调用自行试验一下，这样效果会好得多。工具/原料moreFLAC3D方法/步骤1/5分步阅读

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/5

先来了解一下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/5

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

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

4/5

plo sel geom

显示壳单元

5/5

plo add sketch

加上轮廓线

flac3d