我现在做的这个模型，是一个小的管节推进，就是这个管节前有一个小的掘进挖土机，一边挖，然后后面的小管节一边推，小管节是钢筋砼材料。我在建模时，掘进挖土机就省略了，直接就建立了一个管节，然后通过model null来挖土。
然后，管节与土之间应该建立接触面，但是现在我有几个问题一直没弄清：
1，接触面究竟是应该建立在管节上还是建立在周围土体上？
2，接触面如何删除？我输入delet int 1直接把所有网格单元都删除了，很头疼。。。
3，为何我的管节移动的时候会把接触面牵拉着走呢？不能设置在接触面上滑动么？
附代码：
建立土体模型的代码略去
掘进开挖代码：
new
rest ini.sav
plo con zdis ou on magf 1
;plo add int
;====初始化====
ini xdisp=0.0 ydisp=0.0 zdisp=0
ini xvel=0.0 yvel=0.0 zvel=0.0
sel node ini xdisp=0.0 ydisp=0.0 zdisp=0.0
sel node ini xvel=0.0 yvel=0.0 zvel=0.0
ini state=0
;====以下便开始掘进施工，每次土体挖一小块，管节就用后面的千斤顶推进去，当然我是直接让他移动了，不是用推力让他向前
;====1-exca：通道22====
def gen_tube; 分别建立每一条管节
i=2
j=2
p0x=1.53*(j-1)
p0x_a=p0x+0.17
p0x_b=p0x+1.36
p1x=1.53*j
p0z=1.53*(i-1)
p0z_a=p0z+0.17
p0z_b=p0z+1.36
p1z=1.53*i
tube_number='tube'+string(i)+string(j)
n=10*i+j
command
gen zone brick p0 p0x -15 p0z p1 p1x -15 p0z p2 p0x -5.000 p0z p3 p0x -15 p0z_a size 9 20 1 group tube_number
gen zone brick p0 p0x -15 p0z_b p1 p1x -15 p0z_b p2 p0x -5.000 p0z_b p3 p0x -15 p1z size 9 20 1 group tube_number
gen zone brick p0 p0x -15 p0z_a p1 p0x_a -15 p0z_a p2 p0x -5.000 p0z_a p3 p0x -15 p0z_b size 1 20 7 group tube_number
gen zone brick p0 p0x_b -15 p0z_a p1 p1x -15 p0z_a p2 p0x_b -5.000 p0z_a p3 p0x_b -15 p0z_b size 1 20 7 group tube_number
gen separate tube_number
attach face range x p0x p1x y -15 -5 z 0.16 0.18
attach face range x p0x p1x y -15 -5 z 1.35 1.37
interface n face range group tube_number
interface n prop kn 2e7 ks 2e7 coh 15e4 fric 15 ten 15e4
fix z range y -15 -0.1 z p0z
endcommand
end
gen_tube; 管节建立完毕，初始约束为底面的z向约束
set grav 10
def exca_begin1 ;定义掘进开挖
command
model elas range group tube22 ;首先设置好管节的弹性参数
ini dens 7850 range group tube22
prop bu 1.743e11 sh 8.077e10 range group tube22
model null range x 1.53 3.06 y 0 0.5 z 1.53 3.06 ;第一步开挖，只有0.5m长
ini y add 5.5 range group tube22 ;把管节送进去
attach face
gen separate tube22;如果不用这两步，第一次管节送进去与土接触时，上面的土就会弹开（为什么。。），见图3
solve ratio 1.0e-4
endcommand
loop k(2,20)
exca='exca_stage'+'22'+string(k)
kkk=k*0.5
kkk_=kkk-0.5
command
model null range exca
ini y add 0.5 range group tube22
pause
attach face
gen separate tube22 ;同样用了这样的办法，但是已经报错了，土体都被粘起来了，如图4，图5，所以这个就是我一直纠结的点，求指教。。。
solve ratio 1.0e-4
endcommand
endloop
command
free z range group tube22
solve ratio 1.0e-4 ;管节全部送入之后，再将其z向约束去除，计算一下实际应力
endcommand
end
exca_begin1
save tube22_kongwa.sav
set log off
看到那个头没有，土体网格粘起来了，所以一动，网格就变形了。。
