n
plot block group
gen zone uwedge p0 30.18 0 -25 p1 34.59 0 -25 p2 30.18 1 -25 p3 30.18 0 0 &
p4 34.59 1 -25 p5 30.18 1 0 size 5 1 25 group 1
gen zone uwedge p0 34.59 0 0 p1 30.18 0 0 p2 34.59 1 0 p3 34.59 0 -25 &
p4 30.18 1 0 p5 34.59 1 -25 size 5 1 25 group 2
;接触面
gen separate 2
int 1 wrap 1 2
int 1 maxedge 0.5
gen zone brick p0 0 0 0 p1 30.18 0 0 p2 0 0 -25 p3 0 1 0 p4 30.18 0 -25 &
p5 0 1 -25 p6 30.18 1 0 p7 30.18 1 -25 size 30 25 1
gen zone brick p0 34.59 0 0 p1 50 0 0 p2 34.59 0 -25 p3 34.59 1 0 p4 50 0 -25 &
p5 34.59 1 -25 p6 50 1 0 p7 50 1 -25 size 16 25 1
group soil-1 range z 0 -6
group soil-2 range z -6 -10
group soil-3 range z -10 -25
;模型参数
model mohr range group soil-1
model mohr range group soil-2
model mohr range group soil-3
prop bulk 4.4e10 shear 1.6e10 fric 17.2 coh 16.8e10 tens 1e10 dil 10 range group soil-1
prop bulk 5.7e10 shear 2.1e10 fric 17.7 coh 16.8e10 tens 1e10 dil 10 range group soil-2
prop bulk 10e10 shear 3.6e10 fric 20.9 coh 21.6e10 tens 1e10 dil 10 range group soil-3
interface 1 prop kn 1e8 ks 1e8 fric 15.3 coh 3.6e10 range z 0 -6
interface 1 prop kn 1e8 ks 1e8 fric 16 coh 8.48e10 range z -6 -10
interface 1 prop kn 1e8 ks 1e8 fric 20 coh 4.68e10 range z -10 -25
ini dens 1730 range group soil-1
ini dens 1860 range group soil-2
ini dens 1920 range group soil-3
;边界条件
fix x range x -0.1 0.1
fix x range x 49.9 50.1
fix y range y -0.1 0.1
fix y range y 0.9 1.1
fix z range z -25.1 -24.9
set grav 0 0 -9.8
solve
prop bulk 19e6 shear 8e6 fric 14.8 coh 21e3 tens 1e3 range group soil-1
prop bulk 26e6 shear 10e6 fric 15.2 coh 21e3 tens 1e3 range group soil-2
prop bulk 45e6 shear 18e6 fric 17.9 coh 27e3 tens 1e3 range group soil-3
interface 1 prop kn 1e6 ks 1e6 fric 14 coh 4.6e3 range z 0 -6
interface 1 prop kn 1e6 ks 1e6 fric 13.7 coh 10.6e3 range z -6 -10
interface 1 prop kn 1e6 ks 1e6 fric 17 coh 5.9e3 range z -10 -25
hist unbal
solve
save 初始平衡.sav
res 初始平衡.sav
;pile单元
sel pile id=1 begin=(39.6,0.5,0) end=(39.6,0.5,-16) nseg=16
sel pile id=1 prop emod=2.6e10 nu=0.2 density=2200 xcarea=0.5027 per=2.51 &
xciy=0.001257 xciz=0.001257 xcj=0.002514 cs_scoh=5.78e4 &
cs_sfric=18.62 cs_sk=2.871e8 cs_nk=2.871e8 &
cs_ncoh=5.78e4 cs_nfric=18.62 cs_ngap=off
solve
save 加桩平衡.sav
res 加桩平衡.sav
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
sel node ini xdis 0 ydis 0 zdis 0
sel node ini xvel 0 yvel 0 zvel 0
sel node ini xrdis 0 yrdis 0 zrdis 0
sel node ini xrvel 0 yrvel 0 zrvel 0
set large
group 第一次开挖 range x 40 50 y 0 1 z 0 -3
;设置监测点
hist id=10 gp zdisp 31 0 0
hist id=11 gp zdisp 33 0 0
hist id=12 gp zdisp 35 0 0
hist id=13 gp zdisp 37 0 0
hist id=14 gp zdisp 39 0 0
hist id=15 gp zdisp 41 0 -3
hist id=16 gp zdisp 43 0 -3
hist id=17 gp z
disp 45 0 -3
hist id=18 gp zdisp 47 0 -3
hist id=19 gp zdisp 49 0 -3
hist id=20 gp xdisp 40 0 0
hist id=21 gp xdisp 40 0 -1
hist id=22 gp xdisp 40 0 -2
hist id=23 gp xdisp 40 0 -3
hist id=24 gp xdisp 40 0 -4
hist id=25 gp xdisp 40 0 -5
hist id=26 gp xdisp 40 0 -6
hist id=27 gp xdisp 40 0 -7
hist id=28 gp xdisp 40 0 -8
hist id=29 gp xdisp 40 0 -9
hist id=30 gp xdisp 40 0 -10
hist id=31 gp xdisp 40 0 -11
hist id=32 gp xdisp 40 0 -12
hist id=33 gp xdisp 40 0 -13
hist id=34 gp xdisp 40 0 -14
hist id=35 gp xdisp 40 0 -15
hist id=36 gp xdisp 40 0 -16
hist id=61 sel pile force fz 39.6 0.5 -0.5
hist id=62 sel pile force fz 39.6 0.5 -1.5
hist id=63 sel pile force fz 39.6 0.5 -2.5
hist id=64 sel pile force fz 39.6 0.5 -3.5
hist id=65 sel pile force fz 39.6 0.5 -4.5
hist id=66 sel pile force fz 39.6 0.5 -5.5
hist id=67 sel pile force fz 39.6 0.5 -6.5
hist id=68 sel pile force fz 39.6 0.5 -7.5
hist id=69 sel pile force fz 39.6 0.5 -8.5
hist id=70 sel pile force fz 39.6 0.5 -9.5
hist id=71 sel pile force fz 39.6 0.5 -10.5
hist id=72 sel pile force fz 39.6 0.5 -11.5
hist id=73 sel pile force fz 39.6 0.5 -12.5
hist id=74 sel pile force fz 39.6 0.5 -13.5
hist id=75 sel pile force fz 39.6 0.5 -14.5
hist id=76 sel pile force fz 39.6 0.5 -15.5
hist id=81 sel pile mom mz 39.6 0.5 -0.5
hist id=82 sel pile mom mz 39.6 0.5 -1.5
hist id=83 sel pile mom mz 39.6 0.5 -2.5
hist id=84 sel pile mom mz 39.6 0.5 -3.5
hist id=85 sel pile mom mz 39.6 0.5 -4.5
hist id=86 sel pile mom mz 39.6 0.5 -5.5
hist id=87 sel pile mom mz 39.6 0.5 -6.5
hist id=88 sel pile mom mz 39.6 0.5 -7.5
hist id=89 sel pile mom mz 39.6 0.5 -8.5
hist id=90 sel pile mom mz 39.6 0.5 -9.5
hist id=91 sel pile mom mz 39.6 0.5 -10.5
hist id=92 sel pile mom mz 39.6 0.5 -11.5
hist id=93 sel pile mom mz 39.6 0.5 -12.5
hist id=94 sel pile mom mz 39.6 0.5 -13.5
hist id=95 sel pile mom mz 39.6 0.5 -14.5
hist id=96 sel pile mom mz 39.6 0.5 -15.5
model null range group 第一次开挖
solve
save 第一次开挖.sav
res 第一次开挖.sav
group 第二次开挖 range x 40 50 y 0 1 z -3 -6
;设置监测点
hist id=110 gp zdisp 31 0 0
hist id=111 gp zdisp 33 0 0
hist id=112 gp zdisp 35 0 0
hist id=113 gp zdisp 37 0 0
hist id=114 gp zdisp 39 0 0
hist id=115 gp zdisp 41 0 -6
hist id=116 gp zdisp 43 0 -6
hist id=117 gp zdisp 45 0 -6
hist id=118 gp zdisp 47 0 -6
hist id=119 gp zdisp 49 0 -6
hist id=120 gp xdisp 40 0 0
hist id=121 gp xdisp 40 0 -1
hist id=122 gp xdisp 40 0 -2
hist id=123 gp xdisp 40 0 -3
hist id=124 gp xdisp 40 0 -4
hist id=125 gp xdisp 40 0 -5
hist id=126 gp xdisp 40 0 -6
hist id=127 gp xdisp 40 0 -7
hist id=128 gp xdisp 40 0 -8
hist id=129 gp xdisp 40 0 -9
hist id=130 gp xdisp 40 0 -10
hist id=131 gp xdisp 40 0 -11
hist id=132 gp xdisp 40 0 -12
hist id=133 gp xdisp 40 0 -13
hist id=134 gp xdisp 40 0 -14
hist id=135 gp xdisp 40 0 -15
hist id=136 gp xdisp 40 0 -16
hist id=161 sel pile force fz 39
.6 0.5 -0.5
hi
st id=162 sel pile force fz 39.6 0.5 -1.5
hist id=163 sel pile force fz 39.6 0.5 -2.5
hist id=164 sel pile force fz 39.6 0.5 -3.5
hist id=165 sel pile force fz 39.6 0.5 -4.5
hist id=166 sel pile force fz 39.6 0.5 -5.5
hist id=167 sel pile force fz 39.6 0.5 -6.5
hist id=168 sel pile force fz 39.6 0.5 -7.5
hist id=169 sel pile force fz 39.6 0.5 -8.5
hist id=170 sel pile force fz 39.6 0.5 -9.5
hist id=171 sel pile force fz 39.6 0.5 -10.5
hist id=172 sel pile force fz 39.6 0.5 -11.5
hist id=173 sel pile force fz 39.6 0.5 -12.5
hist id=174 sel pile force fz 39.6 0.5 -13.5
hist id=175 sel pile force fz 39.6 0.5 -14.5
hist id=176 sel pile force fz 39.6 0.5 -15.5
hist id=181 sel pile mom mz 39.6 0.5 -0.5
hist id=182 sel pile mom mz 39.6 0.5 -1.5
hist id=183 sel pile mom mz 39.6 0.5 -2.5
hist id=184 sel pile mom mz 39.6 0.5 -3.5
hist id=185 sel pile mom mz 39.6 0.5 -4.5
hist id=186 sel pile mom mz 39.6 0.5 -5.5
hist id=187 sel pile mom mz 39.6 0.5 -6.5
hist id=188 sel pile mom mz 39.6 0.5 -7.5
hist id=189 sel pile mom mz 39.6 0.5 -8.5
hist id=190 sel pile mom mz 39.6 0.5 -9.5
hist id=191 sel pile mom mz 39.6 0.5 -10.5
hist id=192 sel pile mom mz 39.6 0.5 -11.5
hist id=193 sel pile mom mz 39.6 0.5 -12.5
hist id=194 sel pile mom mz 39.6 0.5 -13.5
hist id=195 sel pile mom mz 39.6 0.5 -14.5
hist id=196 sel pile mom mz 39.6 0.5 -15.5
model null range group 第二次开挖
solve
save 第二次开挖.sav
res 第二次开挖.sav
group 第三次开挖 range x 40 50 y 0 1 z -6 -8
;设置监测点
hist id=210 gp zdisp 31 0 0
hist id=211 gp zdisp 33 0 0
hist id=212 gp zdisp 35 0 0
hist id=213 gp zdisp 37 0 0
hist id=214 gp zdisp 39 0 0
hist id=215 gp zdisp 41 0 -8
hist id=216 gp zdisp 43 0 -8
hist id=217 gp zdisp 45 0 -8
hist id=218 gp zdisp 47 0 -8
hist id=219 gp zdisp 49 0 -8
hist id=220 gp xdisp 40 0 0
hist id=221 gp xdisp 40 0 -1
hist id=222 gp xdisp 40 0 -2
hist id=223 gp xdisp 40 0 -3
hist id=224 gp xdisp 40 0 -4
hist id=225 gp xdisp 40 0 -5
hist id=226 gp xdisp 40 0 -6
hist id=227 gp xdisp 40 0 -7
hist id=228 gp xdisp 40 0 -8
hist id=229 gp xdisp 40 0 -9
hist id=230 gp xdisp 40 0 -10
hist id=231 gp xdisp 40 0 -11
hist id=232 gp xdisp 40 0 -12
hist id=233 gp xdisp 40 0 -13
hist id=234 gp xdisp 40 0 -14
hist id=235 gp xdisp 40 0 -15
hist id=236 gp xdisp 40 0 -16
hist id=261 sel pile force fz 39.6 0.5 -0.5
hist id=262 sel pile force fz 39.6 0.5 -1.5
hist id=263 sel pile force fz 39.6 0.5 -2.5
hist id=264 sel pile force fz 39.6 0.5 -3.5
hist id=265 sel pile force fz 39.6 0.5 -4.5
hist id=266 sel pile force fz 39.6 0.5 -5.5
hist id=267 sel pile force fz 39.6 0.5 -6.5
hist id=268 sel pile force fz 39.6 0.5 -7.5
hist id=269 sel pile force fz 39.6 0.5 -8.5
hist id=270 sel pile force fz 39.6 0.5 -9.5
hist id=271 sel pile force fz 39.6 0.5 -10.5
hist id=272 sel pile force fz 39.6 0.5 -11.5
hist id=273 sel pile force
fz 39.6 0.5 -12.5
hist id=274 sel pile force fz 39.6 0.5 -13.5
his
t id=275 sel pile force fz 39.6 0.5 -14.5
hist id=276 sel pile force fz 39.6 0.5 -15.5
hist id=281 sel pile mom mz 39.6 0.5 -0.5
hist id=282 sel pile mom mz 39.6 0.5 -1.5
hist id=283 sel pile mom mz 39.6 0.5 -2.5
hist id=284 sel pile mom mz 39.6 0.5 -3.5
hist id=285 sel pile mom mz 39.6 0.5 -4.5
hist id=286 sel pile mom mz 39.6 0.5 -5.5
hist id=287 sel pile mom mz 39.6 0.5 -6.5
hist id=288 sel pile mom mz 39.6 0.5 -7.5
hist id=289 sel pile mom mz 39.6 0.5 -8.5
hist id=290 sel pile mom mz 39.6 0.5 -9.5
hist id=291 sel pile mom mz 39.6 0.5 -10.5
hist id=292 sel pile mom mz 39.6 0.5 -11.5
hist id=293 sel pile mom mz 39.6 0.5 -12.5
hist id=294 sel pile mom mz 39.6 0.5 -13.5
hist id=295 sel pile mom mz 39.6 0.5 -14.5
hist id=296 sel pile mom mz 39.6 0.5 -15.5
model null range group 第三次开挖
solve
save 第三次开挖.sav
res 第三次开挖.sav
group 第四次开挖 range x 40 50 y 0 1 z -8 -10
;设置监测点
hist id=310 gp zdisp 31 0 0
hist id=311 gp zdisp 33 0 0
hist id=312 gp zdisp 35 0 0
hist id=313 gp zdisp 37 0 0
hist id=314 gp zdisp 39 0 0
hist id=315 gp zdisp 41 0 -10
hist id=316 gp zdisp 43 0 -10
hist id=317 gp zdisp 45 0 -10
hist id=318 gp zdisp 47 0 -10
hist id=319 gp zdisp 49 0 -10
hist id=320 gp xdisp 40 0 0
hist id=321 gp xdisp 40 0 -1
hist id=322 gp xdisp 40 0 -2
hist id=323 gp xdisp 40 0 -3
hist id=324 gp xdisp 40 0 -4
hist id=325 gp xdisp 40 0 -5
hist id=326 gp xdisp 40 0 -6
hist id=327 gp xdisp 40 0 -7
hist id=328 gp xdisp 40 0 -8
hist id=329 gp xdisp 40 0 -9
hist id=330 gp xdisp 40 0 -10
hist id=331 gp xdisp 40 0 -11
hist id=332 gp xdisp 40 0 -12
hist id=333 gp xdisp 40 0 -13
hist id=334 gp xdisp 40 0 -14
hist id=335 gp xdisp 40 0 -15
hist id=336 gp xdisp 40 0 -16
hist id=361 sel pile force fz 39.6 0.5 -0.5
hist id=362 sel pile force fz 39.6 0.5 -1.5
hist id=363 sel pile force fz 39.6 0.5 -2.5
hist id=364 sel pile force fz 39.6 0.5 -3.5
hist id=365 sel pile force fz 39.6 0.5 -4.5
hist id=366 sel pile force fz 39.6 0.5 -5.5
hist id=367 sel pile force fz 39.6 0.5 -6.5
hist id=368 sel pile force fz 39.6 0.5 -7.5
hist id=369 sel pile force fz 39.6 0.5 -8.5
hist id=370 sel pile force fz 39.6 0.5 -9.5
hist id=371 sel pile force fz 39.6 0.5 -10.5
hist id=372 sel pile force fz 39.6 0.5 -11.5
hist id=373 sel pile force fz 39.6 0.5 -12.5
hist id=374 sel pile force fz 39.6 0.5 -13.5
hist id=375 sel pile force fz 39.6 0.5 -14.5
hist id=376 sel pile force fz 39.6 0.5 -15.5
hist id=381 sel pile mom mz 39.6 0.5 -0.5
hist id=382 sel pile mom mz 39.6 0.5 -1.5
hist id=383 sel pile mom mz 39.6 0.5 -2.5
hist id=384 sel pile mom mz 39.6 0.5 -3.5
hist id=385 sel pile mom mz 39.6 0.5 -4.5
hist id=386 sel pile mom mz 39.6 0.5 -5.5
hist id=387 sel pile mom mz 39.6 0.5 -6.5
hist id=388 sel pile mom mz 39.6 0.5 -7.5
hist id=389 sel pile mom m
z 39.6 0.5 -8.5
hist id=390 sel pile mom mz 39.6 0.5 -9.5
hist id=391 sel pile mom mz 39.6 0.5 -10.5
hist id=392 sel pil
e mom mz 39.6 0.5 -11.5
hist id=393 sel pile mom mz 39.6 0.5 -12.5
hist id=394 sel pile mom mz 39.6 0.5 -13.5
hist id=395 sel pile mom mz 39.6 0.5 -14.5
hist id=396 sel pile mom mz 39.6 0.5 -15.5
model null range group 第四次开挖
solve
save 第四次开挖.sav