请教例子中的零刚度节点怎么处理

chengjia02

命令如下，有心人帮忙解决一下
new
gen zon b p0 0 0 0 p1 16.13 0 0 p2 0 1 0 p3 0 0 10 &
          p4 16.3 1 0 p5 0 1 10 p6 11.78 0 10 p7 11.78 1 10 &
          size 32 2 20
gen zon b p0 0 0 10 p1 11.78 0 10 p2 0 1 10 p3 0 0 16.2 &
          p4 11.78 1 10 p5 0 1 16.2 p6 11.24 0 16.2 p7 11.24 1 16.2 &
          size 22 2 12
gen zon b p0 0 0 16.2 p1 8.5 0 16.2 p2 0 1 16.2 p3 0 0 25 &
          p4 8.5 1 16.2 p5 0 1 25 p6 8.1 0 25 p7 8.1 1 25 &
          size 17 2 18
gen zon b p0 8.5 0 16.2 p1 11.24 0 16.2 p2 8.88 1 16.2 p3 8.6 0 25 &
          p4 8.5 1 16.2 p5 8.6 1 25 p6 10.48 0 25 p7 10.48 1 25 &
          size 4 2 18
gen merge 1e-5 range z(0,25)
attach face all
m m
pro bulk 2e8 shear 1e8 coh 2.5e5 fric 44 tens 5e5
ini dens 2000
set grav 0 0 -9.8
fix x range x -.1 .1
fix x range x 16.12 16.14
fix y
fix z range z -.1 .1
his unbal
sol

chengjia02

没人有空理我....