大家好，可以帮我看一下这个Fish语言定义的循环开挖支护么？支护实现了循环，但是开挖却没有，不知道为什么，后面把支护去掉又可以实现循环开挖了，不知道为什么会这样。刚学习，还希望大家可以帮指导一下…………红色部分没有实现，谢谢大家了
new
title
lined tunnel with slip
ge zo radc p0 0 0 0 p1 5 0 0 p2 0 2 0 p3 0 0 5 &
dim 1.0 ratio 1.0 1.0 1.0 1.2 size 1 8 4 4 fill group tunnel
ge zo radc p0 0 2 0 p1 5 2 0 p2 0 5 0 p3 0 2 5 &
dim 1.0 ratio 1.0 1.0 1.0 1.2 size 1 3 4 4 fill group tunnel
model elastic
pro bu 50e6 sh 18e6
in sxx -1e6 syy -1e6 szz -1e6
fix x range x -0.1 0.1
fix y range y -0.1 0.1
fix z range z -0.1 0.1
fix x z y range x 4.9 5.1
fix x y z range y 4.9 5.1
fix x y z range z 4.9 5.1
hist id=1 unb
;第一阶段无支护开挖
model null range y 0 1 group tunnel
plot create the_system
plot set rot 0 0 50
plot add axes scale 0.075 black
plot add surface yellow
plot show
solve
;
in xdisp 0.0 ydisp 0.0 zdisp 0.0
hist id=2 gp zdisp 0 1 1
hist id=3 gp zvel 0 1 1
save liner-tun1.sav
;第二阶段开挖
model null range y 1 2 group tunnel
save liner-tun1-temp.sav
solve
restore liner-tun1.sav
def kaiwa
loop a(1,b)
c=a+1
c2=float(a)+0.01
c3=a-1
id1=99+a*3
id2=id1+1
id3=id2+1
;开挖
command
model null range y a c group tunnel
end_command
command
sel liner id=a range cyl end1 0 c3 0 end2 0 c2 0 rad 1.0
;3个对称面约束
sel node local xdir 1 0 0 ydir 0 -1 0 range x -0.1 0.1
sel node fix lsys range x -0.1 0.1
sel node fix x yr zr range x -0.1 0.1
sel node local xdir 0 0 -1 ydir 0 -1 0 range z -0.1 0.1
sel node fix lsys range z -0.1 0.1
sel node fix z xr yr range z -0.1 0.1
sel node fix y xr zr range y -0.1 0.1
sel liner id=a prop iso=(10.5e9 0.25) thick=0.2 &
cs_nk=7.4e10 cs_sk=7.4e10 cs_scoh=1e20
hist id=id1 gp zdisp 0 1 1
hist id=id2 sel node zdisp id=25
hist id=id3 sel linersel coupling disp normal nd1 cid=32
plot add sel geom scale 0.005 shrink 0.0 black black
plot add hist id1 id2
plot show
step 1000
end_command
end_loop
end
set b=4
kaiwa
