各位大佬有做盾构开挖的吗
我是根据一篇论文练习的，分几种情况，图一是不考虑护盾和管片，图二是不考虑护盾考虑管片。
我计算出来的 考虑管片情况下的竖向位移比不考虑的情况 还大，各位能帮我看一下命令流哪问题出在哪里吗，感激不尽

（图一）
model new
zone create radial-cylinder p 0 0 0 0 p 1 20 0 0 p 2 0 30 0 p 3 0 0 20 size 2 20 5 18 dim 2.815 2.815 2.815 2.815 ratio 1 1 1 1.2 group 'rock1'
zone reflect normal 0 0 1 range group 'rock1'
zone create cylinder p 0 0 0 0 p 1 2.815 0 0 p 2 0 30 0 p 3 0 0 2.815 size 5 20 5 group 'tunnel1'
zone reflect normal 0 0 1 range group 'tunnel1'
zone reflect normal -1 0 0 range group 'rock1'
zone group 'rock2' range position-x 0 20 group 'rock1'
zone reflect normal -1 0 0 range position-x 0 20 group 'tunnel1'
zone group 'tunnel2' range position-x 0 20 group 'tunnel1'
zone cmodel assign strain-softening
zone initialize density 2448
zone property young 1e9 bulk 1.11e9 coh 0.284e9 friction 36 poisson 0.35 shear 0.29e9
model gravity 0 0 -9.81
zone gridpoint fix velocity-x range position-x -20.01 -19.99
zone gridpoint fix velocity-x range position-x 19.99 20.01
zone gridpoint fix velocity-y range position-y -0.01 0.01
zone gridpoint fix velocity-y range position-y 29.99 30.01
zone gridpoint fix velocity-z range position-z -20.01 -19.99
zone gridpoint fix velocity-x range position-z -20.01 -19.99
zone gridpoint fix velocity-y range position-z -20.01 -19.99
zone initialize stress-xx -14.2e6
zone initialize stress-yy -22.1e6
zone initialize stress-zz -21.9e6
zone face apply stress-normal -21.9e6 range position-z 19.99 20.01
model solve
model save 'insol.sav'
model restore 'insol.sav'
zone gridpoint initialize displacement-x 0.0
zone gridpoint initialize displacement-y 0.0
zone gridpoint initialize displacement-z 0.0
zone cmodel assign null range position-y 0 15 group 'tunnel1'
zone cmodel assign null range position-y 0 15 group 'tunnel2'
structure shell create by-face id 1 range position-x -2.815 2.815 position-z -2.815 2.815 cylinder end-1 0 0 0 end-2 0 15 0 radius 2.815
structure shell property isotropic 10e4 0.1 thickness 0.1
model solve ratio 1e-4
model save 'step1.sav'

（图二1）
model new
zone create radial-cylinder p 0 0 0 0 p 1 20 0 0 p 2 0 30 0 p 3 0 0 20 size 2 20 5 18 dim 2.815 2.815 2.815 2.815 ratio 1 1 1 1.2 group 'rock1'
zone reflect normal 0 0 1 range group 'rock1'
zone create cylindrical-shell p 0 0 0 0 p 1 2.815 0 0 p 2 0 30 0 p 3 0 0 2.815 size 1 20 5 5 dim 2.465 2.465 2.465 2.465 group 'guanpian1'
zone reflect normal 0 0 1 range group 'guanpian1'
zone reflect normal -1 0 0 range group 'guanpian1'
zone group 'guanpian2' range position-x 0 2.815 group 'guanpian1'
zone create cylinder p 0 0 0 0 p 1 2.465 0 0 p 2 0 30 0 p 3 0 0 2.465 size 5 40 5 group 'tunnel1'
zone reflect normal 0 0 1 range group 'tunnel1'
zone reflect normal -1 0 0 range group 'rock1'
zone group 'rock2' range position-x 0 20 group 'rock1'
zone reflect normal -1 0 0 range position-x 0 20 group 'tunnel1'
zone group 'tunnel2' range position-x 0 20 group 'tunnel1'
zone cmodel assign strain-softening
zone initialize density 2448
zone property young 1e9 bulk 1.11e9 coh 0.284e9 friction 36 poisson 0.35 shear 0.29e9
model gravity 0 0 -9.81
zone gridpoint fix velocity-x range position-x -20.01 -19.99
zone gridpoint fix velocity-x range position-x 19.99 20.01
zone gridpoint fix velocity-y range position-y -0.01 0.01
zone gridpoint fix velocity-y range position-y 29.99 30.01
zone gridpoint fix velocity-z range position-z -20.01 -19.99
zone gridpoint fix velocity-x range position-z -20.01 -19.99
zone gridpoint fix velocity-y range position-z -20.01 -19.99
zone initialize stress-xx -14.2e6
zone initialize stress-yy -22.1e6
zone initialize stress-zz -21.9e6
zone face apply stress-normal -21.9e6 range position-z 19.99 20.01
model solve ratio 1e-4
model save 'insol.sav'
model restore 'insol.sav'
zone gridpoint initialize displacement-x 0.0
zone gridpoint initialize displacement-y 0.0
zone gridpoint initialize displacement-z 0.0
zone cmodel assign null range position-y 0 13.5 group 'tunnel1'
zone cmodel assign null range position-y 0 13.5 group 'guanpian1'

（图二2）
zone cmodel assign null range position-y 0 13.5 group 'tunnel2'
zone cmodel assign null range position-y 0 13.5 group 'guanpian2'
structure shell create by-face id 1 range position-x -2.815 2.815 position-z -2.815 2.815 cylinder end-1 0 0 0 end-2 0 13.5 0 radius 2.815
structure shell property isotropic 10e4 0.1 thickness 0.1
zone cmodel assign strain-softening range position-y 0 1.5 group 'guanpian1'
zone cmodel assign strain-softening range position-y 0 1.5 group 'guanpian2'
zone initialize density 2548
zone property young 25e9 poisson 0.15 range position-y 0 1.5 group 'guanpian1'
zone property young 25e9 poisson 0.15 range position-y 0 1.5 group 'guanpian2'
model solve ratio 1e-4
model save 'step1.sav'
model restore 'step1.sav'
zone cmodel assign null range position-y 13.5 15 group 'tunnel1'
zone cmodel assign null range position-y 13.5 15 group 'guanpian1'
zone cmodel assign null range position-y 13.5 15 group 'tunnel2'
zone cmodel assign null range position-y 13.5 15 group 'guanpian2'
structure shell create by-face id 2 range position-x -2.815 2.815 position-z -2.815 2.815 cylinder end-1 0 13.5 0 end-2 0 15 0 radius 2.815
structure shell property isotropic 10e4 0.1 thickness 0.1
zone cmodel assign strain-softening range position-y 1.5 3 group 'guanpian1'
zone cmodel assign strain-softening range position-y 1.5 3 group 'guanpian2'
zone initialize density 2548
zone property young 25e9 poisson 0.15 range position-y 1.5 3 group 'guanpian1'
zone property young 25e9 poisson 0.15 range position-y 1.5 3 group 'guanpian2'
model solve ratio 1e-4
model save 'step2.sav'

可帮做

我也不知道 可能 类似问题 我的是巷道支护加了钢管混凝土支柱，我是用犀牛导入的 效果也一样 还不如之前不加钢管的支护效果好位移 也是 变大了 可能是 需要建立接触面嘛？现在 正研究犀牛的了，要是你先解决了，劳烦通知一下。我准备开学 把这个问题弄明白咯，要是我先的话 我就告诉你 。对 qq群里有个大佬，你可以问问这人，他很吊哈哈哈

shell？怎么不直接用liner？

应该是接触面的问题，你可以好好研究手册里面的接触面