FLAC3D求助！！！！！

我伙惊呆

我想建立一个管线埋设在土体里的模型，土体是实体单元，管线是结构shell单元，命令流在下面，我不知道接触面和管线的radius值取的对不对，所以不知道管线是否跟土体接触了，请大神帮我看一下！！！

我伙惊呆

n
;----------------------------建立管线周围土体-------------------------
;
ge zo radcyl p0 -20.8 18.8 14.5 p1 -20.8 18.8 15.2 p2 -0.8 18.8 14.5 p3 -20.8 20.8 14.5 p4 -0.8 18.8 15.2 p5 -0.8 20.8 14.5 p6 -20.8 20.8 15.2 p7 -0.8 20.8 15.2 &
p8 -20.8 18.8 15 p9 -20.8 19.3 14.5 p10 -0.8 18.8 15 p11 -0.8 19.3 14.5 size 3 5 6 3 group pipesoil1
;管线左上部土体建立（x从-20.8到-0.8）
;
ge zo radcyl p0 -20.8 18.8 14.5 p1 -20.8 20.8 14.5 p2 -0.8 18.8 14.5 p3 -20.8 18.8 13.3 p4 -0.8 20.8 14.5 p5 -0.8 18.8 13.3 p6 -20.8 20.8 13.3 p7 -0.8 20.8 13.3 &
p8 -20.8 19.3 14.5 p9 -20.8 18.8 14 p10 -0.8 19.3 14.5 p11 -0.8 18.8 14 size 3 5 6 3 group pipesoil2
;管线左下部土体建立（x从-20.8到-0.8）
;
ge zo radcyl p0 -0.8 18.8 14.5 p1 -0.8 18.8 15.2 p2 0 18.8 14.5 p3 -0.8 20.8 14.5 p4 0 18.8 15.2 p5 0 20.8 14.5 p6 -0.8 20.8 15.2 p7 0 20.8 15.2 &
p8 -0.8 18.8 15 p9 -0.8 19.3 14.5 p10 0 18.8 15 p11 0 19.3 14.5 size 3 1 6 3 group pipesoil3
;管线左上部土体建立（x从-0.8到0）
;
ge zo radcyl p0 -0.8 18.8 14.5 p1 -0.8 20.8 14.5 p2 0 18.8 14.5 p3 -0.8 18.8 13.3 p4 0 20.8 14.5 p5 0 18.8 13.3 p6 -0.8 20.8 13.3 p7 0 20.8 13.3 &
p8 -0.8 19.3 14.5 p9 -0.8 18.8 14 p10 0 19.3 14.5 p11 0 18.8 14 size 3 1 6 3 group pipesoil4
;管线左下部土体建立（x从-0.8到0）
;此处错误了，需要修改！
;
ge zo radcyl p0 0 18.8 14.5 p1 0 18.8 15.2 p2 40 18.8 14.5 p3 0 20.8 14.5 p4 40 18.8 15.2 p5 40 20.8 14.5 p6 0 20.8 15.2 p7 40 20.8 15.2 &
p8 0 18.8 15 p9 0 19.3 14.5 p10 40 18.8 15 p11 40 19.3 14.5 size 3 10 6 3 group pipesoil5
;管线左上部土体建立(x从0到40)
;
ge zo radcyl p0 0 18.8 14.5 p1 0 20.8 14.5 p2 40 18.8 14.5 p3 0 18.8 13.3 p4 40 20.8 14.5 p5 40 18.8 13.3 p6 0 20.8 13.3 p7 40 20.8 13.3 &
p8 0 19.3 14.5 p9 0 18.8 14 p10 40 19.3 14.5 p11 40 18.8 14 size 3 10 6 3 group pipesoil6
;管线左下部土体建立(x从0到40)
;
ge zo ref nor 0 -1 0 ori -20.8 18.8 14.5 range group pipesoil1 any group pipesoil2 any group pipesoil3 any group pipesoil4 any group pipesoil5 any group pipesoil6 any
;建立对称面
group pipesoil range group pipesoil1 any group pipesoil2 group pipesoil3 any group pipesoil4 any group pipesoil5 any group pipesoil6 any
pl bl gr
;
;----------------------------建立接触面------------------------------
;
interface 1 face ran cyl end1(-20.8,18.8,14.5) end2(40,18.8,14.5) rad 0.5 &
cyl end1(-20.8,18.8,14.5) end2(40,18.8,14.5) rad 0.49 not
;
;-----------------------------建立管线（使用shell结构单元）-------------------------
;
sel shell id=3 ran cyl end1=(-20.8,18.8,14.5) end2=(40,18.8,14.5) rad=0.5
;意思在圆柱半径为0.5的圆柱体表面建立shell结构单元
pl ad sel geom
;显示建立的管线结构单元
;

我伙惊呆

有大神能帮我看下这个问题吗？谢谢了！！！！！