;一号桩外土体模型
gen zone radcyl p0 (0,0,0) p1 (.75,0,0) p2 (0,0,-22) p3 (0,.75,0) &
p4 (.75,0,-22) p5 (0,.75,-22) p6 (.75,.75,0) p7 (.75,.75,-22) &
p8 (.25,0,0) p9 (0,.25,0) p10 (.25,0,-22) p11 (0,.25,-22) &
size 2 22 6 3 ratio 1 1 1 1.25
gen zone radcyl p0 (0,0,-22) p1 (.75,0,-22) p2 (0,0,-45) p3 (0,.75,-22) &
p4 (.75,0,-45) p5 (0,.75,-45) p6 (.75,.75,-22) p7 (.75,.75,-45) &
p8 (.25,0,-22) p9 (0,.25,-22) p10 (.25,0,-45) p11 (0,.25,-45) &
size 2 11 6 3 ratio 1 1 1 1.25 fill
;二号桩外土体模型
gen zone copy 1.5 0 0
gen zone reflect normal 1 0 0 ori 1.5 0 0 ran x 1.5 2.25
;三号桩外土体模型
gen zone copy 0 1.5 0 ran x 0 0.75
gen zone reflect normal 0 1 0 ori 0 1.5 0 ran y 1.5 2.25
;四号桩外土体模型
gen zone copy 0 1.5 0 ran x 0.75 2.25
gen zone reflect normal 0 1 0 ori 0 1.5 0 ran x 0.75 2.25 y 1.5 2.25
;承台外土体模型
gen zone radtunnel p0 (0,0,0) p1 (6,0,0) p2 (0,0,-45) p3 (0,6,0) &
p4 (6,0,-45) p5 (0,6,-45) p6 (6,6,0) p7 (6,6,-45) &
p8 (2.25,0,0) p9 (0,2.25,0) p10 (2.25,0,-45) p11(0,2.25,-45) &
p12 (2.25,2.25,0) p13 (2.25,2.25,-45) &
size 9 45 9 5 ratio 1 1 1 1.3
group soil1 range z 0 -1.03
group soil2 range z -1.03 -2.63
group soil3 range z -2.63 -7.61
group soil4 range z -7.61 -13.29
group soil5 range z -13.29 -17.69
group soil6 range z -17.69 -33.55
group soil7 range z -33.55 -45
model mohr range group soil1
model mohr range group soil2
model mohr range group soil3
model mohr range group soil4
model mohr range group soil5
model mohr range group soil6
model mohr range group soil7
prop bulk 4.0e6 shear 1.25e6 fric 11.4 cohesion 13e3 ten 1e10 range group soil1
ini density 1870 range group soil1
prop bulk 6.67e6 shear 3.08e6 fric 10.8 cohesion 10e3 ten 1e10 range group soil2
ini density 1890 range group soil2
prop bulk 6.67e6 shear 3.33e6 fric 9.8 cohesion 16e3 ten 1e10 range group soil3
ini density 1920 range group soil3
prop bulk 13.33e6 shear 4.44e6 fric 12.5 cohesion 14e3 ten 1e10 range group soil4
ini density 1990 range group soil4
prop bulk 10e6 shear 6e6 fric 13.1 cohesion 12e3 ten 1e10 range group soil5
ini density 1970 range group soil5
prop bulk 16.67e6 shear 7.69e6 fric 16 cohesion 30e3 ten 1e10 range group soil6
ini density 2020 range group soil6
prop bulk 16.67e6 shear 7.69e6 fric 20.8 cohesion 43e3 ten 1e10 range group soil7
ini density 2040 range group soil7
;1号桩接触面
interface 1 face range cylinder end1(0,0,0.1) end2(0,0,-22.1) radius .26 &
cylinder end1(0,0,0.1) end2(0,0,-22.1) radius .24 not
interface 2 face range cylinder end1(0,0,-21.9) end2(0,0,-22.1) radius .26
;2号桩接触面
interface 1 face range cylinder end1(1.5,0,0.1) end2(1.5,0,-22.1) radius .26 &
cylinder end1(1.5,0,0.1) end2(1.5,0,-22.1) radius .24 not
interface 2 face range cylinder end1(1.5,0,-21.9) end2(1.5,0,-22.1) radius .26
;三号桩接触面
interface 1 face range cylinder end1(0,1.5,0.1) end2(0,1.5,-22.1) radius .26 &
cylinder end1(0,1.5,0.1) end2(0,1.5,-22.1) radius .24 not
interface 2 face range cylinder end1(0,1.5,-21.9) end2(0,1.5,-22.1) radius .26
;四号桩接触面
interface 1 face range cylinder end1(1.5,1.5,0.1) end2(1.5,1.5,-22.1) radius .26 &
cylinder end1(1.5,1.5,0.1) end2(1.5,1.5,-22.1) radius .24 not
interface 2 face range cylinder end1(1.5,1.5,-21.9) end2(1.5,1.5,-22.1) radius .26
interface 1 prop kn=2.7e7 ks=2.7e7 fric=8 coh=43e3
interface 2 prop kn=2.7e7 ks=2.7e7 fric=8 coh=43e3
;建立桩体模型
;一号桩
gen zone cyl p0 (0,0,23) p1 (.25,0,23) p2 (0,0,1) p3 (0,.25,23) &
p4 (0.25,0,1) p5 (0,0.25,1) &
size 1 22 6
;二号桩
gen zone copy 1.5 0 0 ran z 1 23
gen zone reflect normal 1 0 0 ori 1.5 0 0 ran x 1.5 2.25 z 1 23
;三号桩
gen zone copy 0 1.5 0 ran x 0 0.75 z 1 23
gen zone reflect normal 0 1 0 ori 0 1.5 0 ran y 1.5 2.25 z 1 23
;4号桩
gen zone copy 0 1.5 0 ran x 0.75 2.25 z 1 23
gen zone reflect normal 0 1 0 ori 0 1.5 0 ran x 0.75 2.25 y 1.5 2.25 z 1 23
group zhuang range z 1 23
ini z add -23.0 range group zhuang
model elas range group zhuang
prop bulk 15.56e6 shear 11.67e6
ini density 2500 range group zhuang
;建立承台模型
gen zone brick p0(0,0,0) p1(2.25,0,0) p2(0,0,-1) p3(0,2.25,0) size 1,1,1
group chengtai range x 0 2.25 y 0 2.25 z 0 -1
model elas range group chengtai
prop bulk 15.56e6 shear 11.67e6
ini density 2500 range group chengtai
gen merge 1e-4 ran cylinder end1 (0,0,-.01) end2 (0,0,.01) radius 0.16
gen merge 1e-4 ran cylinder end1 (1.5,0,-.01) end2 (1.5,0,.01) radius 0.16
gen merge 1e-4 ran cylinder end1 (0,1.5,-.01) end2 (0,1.5,.01) radius 0.16
gen merge 1e-4 ran cylinder end1 (1.5,1.5,-.01) end2 (1.5,1.5,.01) radius 0.16
fix z range z -45.1 -44.9
fix x range x -.1 .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
ini szz -19261 grad 0 0 19261 range z 0 -1.03
ini szz -49801 grad 0 0 30240 range z -1.03 -2.63
ini szz -145417 grad 0 0 95616 range z -2.63 -7.61
ini szz -258449 grad 0 0 113032 range z -7.61 -13.29
ini szz -345129 grad 0 0 86680 range z -13.29 -17.69
ini szz -665501 grad 0 0 320372 range z -17.69 -33.55
ini szz -899081 grad 0 0 233580 range z -33.55 -45
apply szz -48e3 range z .05 .15
solve
plot gp disp ran id 0 0 0
