FLAC3D：[2]殼單元?

本節主要介紹如何在圓形盾構隧道中建立殼單元。示例中，先建立幾何模型，後加入殼單元。

建議讀者把第一步中的命令流自行拷貝到 txt 檔案中，用file--call命令呼叫自行試驗一下，這樣效果會好得多。

工具/原料

FLAC3D

方法/步驟

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

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

最後的結果見下圖。

先來了解一下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 的位置上出現了較為明顯的梯度變化。

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

因為是圓形隧道，只需要用cylinder命令，把周圈圍巖的中心包括在內就可以了。

plo sel geom

顯示殼單元

plo add sketch

加上輪廓線