请高手帮忙，关于桩跟地基土的接触问题

mk05230723

本人根据《基础与工程实例》书中p130页例9-1建立桩土接触模型时，没有得到书中的结果而是地基土网格多出一部分，求高手解答。附命令流及图形。
new
gen zone radcyl p0 (0,0,0) p1 (9,0,0) p2 (0,0,-7) p3 (0,9,0) &
                p4 (9,0,-7) p5 (0,9,-7) p6 (9,9,0) p7 (9,9,-7) &
                p8 (.2,0,0) p9 (0,.2,0) p10 (.2,0,-7) p11 (0,.2,-7) &
                size 3 10 6 15 ratio 1 1 1 1.15
gen zone radcyl p0 (0,0,-7) p1 (9,0,-7) p2 (0,0,-10) p3 (0,9,-7) &
                p4 (9,0,-10) p5 (0,9,-10) p6 (9,9,-7) p7 (9,9,-10) &
                p8 (.2,0,-7) p9 (0,.2,-7) p10 (.2,0,-10) p11 (0,.2,-10) &
                size 3 6 6 15  ratio 1 1 1 1.15 fill
gen zone radcyl p0 (0,0,-10) p1 (9,0,-10) p2 (0,0,-20) p3 (0,9,-10) &
                p4 (9,0,-20) p5 (0,9,-20) p6 (9,9,-10) p7 (9,9,-20) &
                p8 (.2,0,-10) p9 (0,.2,-10) p10 (.2,0,-20) p11 (0,.2,20) &
                size 3 6 6 15  ratio 1 1 1 1.15 fill
gen zone reflect dd 270 dip 90
group clay
interface 1 face range cylinder end1 (0,0,0)  end2 (0,0,-10.1) radius .21 &
                       cylinder end1 (0,0,0)  end2 (0,0,-10.1) radius .19 not
interface 2 face range cylinder end1 (0,0,-9.9)  end2 (0,0,-10.1) radius .21
pause
gen zone cyl p0 (0,0,11) p1 (.2,0,11) p2 (0,0,1) p3 (0,.2,11) &
             p4 (.2,0,1) p5 (0,.2,1) &
             size 3 10 6
gen zone cyl p0 (0,0,11.1) p1 (.2,0,11.1) p2 (0,0,11) p3 (0,.2,11.1) &
             p4 (.2,0,11) p5 (0,.2,11) &
             size 3 1 6
gen zone reflect dd 270 dip 90 range z 1 11.1
group pile range z 1 11.1
ini z add -11.0 range group pile
save pile_geom.sav