刚开始学，按书上的例子改写的，但是书上没有加褥垫层，大家看看写的哪里出错了，谢谢！！！
;带承台组合桩模拟,单桩直径600mm，1、4号桩为长桩20m，2、3号桩为短桩5m
n
;承台下土体模型的建立
;1号桩
gen zone radcyl p0 (0,0,0) p1 (0.9,0,0) p2 (0,0,-20) p3 (0,0.9,0) &
p4 (0.9,0,-20) p5 (0,0.9,-20) p6 (0.9,0.9,0) p7 (0.9,0.9,-20) &
p8 (.3,0,0) p9 (0,.3,0) p10 (.3,0,-20) p11 (0,.3,-20) &
size 3 20 6 3 ratio 1 1 1 1.4
gen zone radcyl p0 (0,0,-20) p1 (0.9,0,-20) p2 (0,0,-25) p3 (0,0.9,-20) &
p4 (0.9,0,-25) p5 (0,0.9,-25) p6 (0.9,0.9,-20) p7 (0.9,0.9,-25) &
p8 (.3,0,-20) p9 (0,.3,-20) p10 (.3,0,-25) p11 (0,.3,-25) &
size 3 5 6 3 ratio 1 1 1 1.4 fill
;2号桩
gen zone radcyl p0 (1.8,0,0) p1 (2.7,0,0) p2 (1.8,0,-5) p3 (1.8,0.9,0) &
p4 (2.7,0,-5) p5 (1.8,0.9,-5) p6 (2.7,0.9,0) p7 (2.7,0.9,-5) &
p8 (2.1,0,0) p9 (1.8,0.3,0) p10 (2.1,0,-5) p11 (1.8,0.3,-5) &
size 3 5 6 3 ratio 1 1 1 1.4
gen zone radcyl p0 (1.8,0,-5) p1 (2.7,0,-5) p2 (1.8,0,-25) p3 (1.8,0.9,-5) &
p4 (2.7,0,-25) p5 (1.8,0.9,-25) p6 (2.7,0.9,-5) p7 (2.7,0.9,-25) &
p8 (2.1,0,-5) p9 (1.8,0.3,-5) p10 (2.1,0,-25) p11 (1.8,0.3,-25) &
size 3 20 6 3 ratio 1 1 1 1.4 fill
gen zone reflect normal 1 0 0 ori 1.8 0 0 range x 1.8 2.7 y 0 0.9
;3号桩
gen zone copy -1.8 1.8 0 ran x 1.8 2.7 y 0 0.9
gen zone reflect normal 0 1 0 ori 0 1.8 0 ran x 0 0.9 y 1.8 2.7
;4号桩
gen zone copy 1.8 1.8 0 ran x 0 0.9 y 0 0.9
gen zone reflect normal 0 1 0 ori 0 1.8 0 ran x 1.8 2.7 y 1.8 2.7
gen zone reflect normal 1 0 0 ori 1.8 0 0 ran x 1.8 2.7 y 0.9 2.7
;承台外土体模型的建立，6m
gen zone radtunnel p0 (0,0,0) p1 (6,0,0) p2 (0,0,-25) p3 (0,6,0) &
p4 (6,0,-25) p5 (0,6,-25) p6 (6,6,0) p7 (6,6,-25) &
p8 (2.7,0,0) p9 (0,2.7,0) p10 (2.7,0,-25) p11 (0,2.7,-25) p12 (2.7,2.7,0) p13 (2.7,2.7,-25) &
size 9 25 9 5 ratio 1 1 1 1.3
group clay ran z -20 0
group sand ran z -20 -25
;1号桩接触面
interface 1 face range cylinder end1 (0,0,0.1) end2 (0,0,-20.1) radius .31 &
cylinder end1 (0,0,0.1) end2 (0,0,-20.1) radius .29 not
interface 2 face range cylinder end1 (0,0,-19.9) end2 (0,0,-20.1) radius .31
;2号桩接触面
interface 3 face range cylinder end1 (1.8,0,0.1) end2 (1.8,0,-5.1) radius .31 &
cylinder end1 (1.8,0,0.1) end2 (1.8,0,-5.1) radius .29 not
interface 4 face range cylinder end1 (1.8,0,-4.9) end2 (1.8,0,-5.1) radius .31
;3号桩接触面
interface 5 face range cylinder end1 (0,1.8,0.1) end2 (0,1.8,-5.1) radius .31 &
cylinder end1 (0,1.8,0.1) end2 (0,1.8,-5.1) radius .29 not
interface 6 face range cylinder end1 (0,1.8,-4.9) end2 (0,1.8,-5.1) radius .31
;4号桩接触面
interface 7 face range cylinder end1 (1.8,1.8,0.1) end2 (1.8,1.8,-20.1) radius .31 &
cylinder end1 (1.8,1.8,0.1) end2 (1.8,1.8,-20.1) radius .29 not
interface 8 face range cylinder end1 (1.8,1.8,-19.9) end2 (1.8,1.8,-20.1) radius .31
;建立桩体模型
;1号桩
gen zone cyl p0 (0,0,21) p1 (.3,0,21) p2 (0,0,1) p3 (0,.3,21) &
p4 (.3,0,1) p5 (0,.3,1) &
size 1 20 6
;2号桩
gen zone cyl p0 (1.8,0,27) p1 (2.1,0,27) p2 (1.8,0,22) p3 (1.8,.3,27) &
p4 (2.1,0,22) p5 (1.8,.3,22) &
size 1 5 6
gen zone reflect normal 1 0 0 ori 1.8 0 0 ran x 1.8 2.1 y 0 0.3 z 22 27
;3号桩
gen zone copy -1.8 1.8 0 ran x 1.8 2.1 y 0 0.3 z 22 27
gen zone reflect normal 0 1 0 ori 0 1.8 0 ran x 0 0.3 y 1.8 2.1 z 22 27
;4号桩
gen zone copy 1.8 1.8 0 ran x 0 0.3 y 0 0.3 z 1 21
gen zone reflect normal 0 1 0 ori 0 1.8 0 ran x 1.8 2.1 y 1.8 2.1 z 1 21
gen zone reflect normal 1 0 0 ori 1.8 0 0 ran x 1.8 2.1 y 1.5 2.1 z 1 21
group zhuang14 range z 1 21
ini z add -21.0 range group zhuang14
group zhuang23 range z 22 27
ini z add -27.0 range group zhuang23
;建立垫层模型
;1号桩
gen zone cyl p0 (0,0,6) p1 (.3,0,6) p2 (0,0,5.7) p3 (0,.3,6) &
p4 (.3,0,5.7) p5 (0,.3,5.7) &
size 1 1 6
gen zone radcyl p0 (0,0,6) p1 (0.9,0,6) p2 (0,0,5.7) p3 (0,0.9,6) &
p4 (0.9,0,5.7) p5 (0,0.9,5.7) p6 (0.9,0.9,6) p7 (0.9,0.9,5.7) &
p8 (.3,0,6) p9 (0,.3,6) p10 (.3,0,5.7) p11 (0,.3,5.7) &
size 3 1 6 3 ratio 1 1 1 1.4
;2号桩
gen zone copy 1.8 0 0 ran z 5.7 6
gen zone reflect normal 1 0 0 ori 1.8 0 0 ran x 1.8 2.7 z 5.7 6
;3号桩
gen zone copy 0 1.8 0 ran x 0 0.9 z 5.7 6
gen zone reflect normal 0 1 0 ori 0 1.8 0 ran y 1.8 2.7 z 5.7 6
;4号桩
gen zone copy 0 1.8 0 ran x 0.9 2.7 z 5.7 6
gen zone reflect normal 0 1 0 ori 0 1.8 0 ran x 0.9 2.7 y 1.8 2.7 z 5.7 6
group DC ran z 5.7 6
interface 9 face range z 5.69 5.71
ini z add -5.7 range group DC
;建立承台模型
;1号桩
gen zone cyl p0 (0,0,6.6) p1 (.3,0,6.6) p2 (0,0,6) p3 (0,.3,6.6) &
p4 (.3,0,6) p5 (0,.3,6) &
size 1 1 6
gen zone radcyl p0 (0,0,6.6) p1 (0.9,0,6.6) p2 (0,0,6) p3 (0,0.9,6.6) &
p4 (0.9,0,6) p5 (0,0.9,6) p6 (0.9,0.9,6.6) p7 (0.9,0.9,6) &
p8 (.3,0,6.6) p9 (0,.3,6.6) p10 (.3,0,6) p11 (0,.3,6) &
size 3 1 6 3 ratio 1 1 1 1.4
;2号桩
gen zone copy 1.8 0 0 ran z 6 6.6
gen zone reflect normal 1 0 0 ori 1.8 0 0 ran x 1.8 2.7 z 6 6.6
;3号桩
gen zone copy 0 1.8 0 ran x 0 0.9 z 6 6.6
gen zone reflect normal 0 1 0 ori 0 1.8 0 ran y 1.8 2.7 z 6 6.6
;4号桩
gen zone copy 0 1.8 0 ran x 0.9 2.7 z 6 6.6
gen zone reflect normal 0 1 0 ori 0 1.8 0 ran x 0.9 2.7 y 1.8 2.7 z 6 6.6
group cap ran z 6 6.6
ini z add -5.7 range group cap
pau
save pile_geom.sav
;
model mohr range group clay
prop bulk 2.9033e6 shear 1.34e6 coh 24.7e3 fric 28.2 range group clay
model mohr range group sand
prop bulk 11.433e6 shear 2.45e6 coh 9.76e3 fric 24.3 range group sand
model elas range group DC
prop bulk 2.9033e6 shear 1.34e6 range group DC
model elas range group cap
prop bulk 2.9033e6 shear 1.34e6 range group cap
model elas range group zhuang14
prop bulk 2.9033e6 shear 1.34e6 range group zhuang14
model elas range group zhuang23
prop bulk 2.9033e6 shear 1.34e6 range group zhuang23
;0.8c、fai
interface 1 prop kn 1e8 ks 1e8 coh 19.76e3 fric 22.56
interface 2 prop kn 1e9 ks 1e9 coh 300000 fric 45
interface 3 prop kn 1e8 ks 1e8 coh 19.76e3 fric 22.56
interface 4 prop kn 1e8 ks 1e8 coh 19.76e3 fric 22.56
interface 5 prop kn 1e8 ks 1e8 coh 19.76e3 fric 22.56
interface 6 prop kn 1e8 ks 1e8 coh 19.76e3 fric 22.56
interface 9 prop kn 1e8 ks 1e8 coh 0 fric 0
interface 7 prop kn 1e8 ks 1e8 coh 19.76e3 fric 22.56
interface 8 prop kn 1e9 ks 1e9 coh 300000 fric 45
;
ini dens 1930 range group clay
ini dens 1550 range group sand
ini dens 1930 range group DC
ini dens 1930 range group cap
ini dens 1930 range group zhuang14
ini dens 1930 range group zhuang23
model null range z -0.1 0.95
fix z range z -25.1 -24.9
fix x range x -0.1 0.1
fix x range x 5.9 6.1
fix y range y -.1 .1
fix y range y 5.9 6.1
set grav 0 0 -10
hist unbal
solve
save pile0.sav
