rt 我想模拟隧道施工，在土层与注浆体之间施加径向节点力，但是用下面这个东西运行缺什么都没有，请大神来解答一下罢，究竟哪里有问题
new
gen zone radcyl p0 0,0,0 p1 add 18,0,0 p2 add 0,8,0 &
p3 add 0,0,18 dim 1.75,1.75,1.75,1.75 &
ratio 1.0,0.83,1.0,1.2 size 6,6,6,10 group tunnel
gen zone cyl p0 0,0,0 p1 1.75,0,0 p2 0,8,0 p3 0,0,1.75 &
size 6,6,6 ratio 1.0,0.83,1.0 group soil
mo mo
prop bulk 33e9 she 25e9 fric 45 coh 10e6 ten 5e6
fix x range x -.1,.1
fix x range x 17.9,18.1
fix y range y -.1,.1
fix y range y 7.9,8.1
fix z range z -.1,.1
def aaaa
pnt=gp_head
loop while pnt # null
_zx=gp_xpos(pnt)
_zy=gp_ypos(pnt)
_zr=sqrt(_zx*_zx+_zy*_zy)
_zid=gp_id(pnt)
if _zr=1.75 then
f_x=f_1*_zx/_zr
f_y=f_1*_zx/_zr
command
apply xforce f_x range id=_zid
apply yforce f_y range id=_zid
endcommand
endif
pnt=gp_next(pnt)
endloop
end
f_1=22e3
aaaa
step 1
plot add ax
plot sk
plot fob range cyl end1 0,0,0 end2 0,8,0 r 1.75
plot fap
plot show
