根据工程模拟单排悬臂桩基坑开挖

caorenyanlian

该命令流还有问题，桩顶向基坑外侧便宜，希望大家讨论下，帮我更正更正

new
;单排桩
gen zone brick size 100 2 70 &
p0 0 0 0 p1 50 0 0 p2 0 1.4 0 p3 0 0 35.2

model mohr

prop bulk 35.67e6 shear 18e6 coh 1e10 ten 1e10 fric 35 range z 0 8.4
prop bulk 31.50e6 shear 19.14e6 coh 1e10 ten 1e10 fric 30 range z 8.4 17.4
prop bulk 30.06e6 shear 17.13e6 coh 1e10 ten 1e10 fric 28 range z 17.4 25.6
prop bulk 26.67e6 shear 14.22e6 coh 1e10 ten 1e10 fric 20 range z 25.6 27.5
prop bulk 25.26e6 shear 15.79e6 coh 1e10 ten 1e10 fric 16.5 range z 27.5 31.2
prop bulk 22.17e6 shear 13.89e6 coh 1e10 ten 1e10 fric 24 range z 31.2 33.4
prop bulk 17.3e6 shear 12.11e6 coh 1e10 ten 1e10 fric 18 range z 33.4 35.2

ini dens 1900 range z 0 8.4
ini dens 1850.9 range z 8.4 17.4
ini dens 1834.9 range z 17.4 25.6
ini dens 1744.9 range z 25.6 27.5
ini dens 1820.9 range z 27.5 31.2
ini dens 1763.5 range z 31.2 33.4
ini dens 1807.2 range z 33.4 35.2

ini szz -700e3 grad 0 0 20e3
ini syy -350e3 grad 0 0 10e3
ini sxx -350e3 grad 0 0 10e3


sel pile id=1 begin 30 0.7 35.2 end 30 0.7 22 nseg 10

sel pile prop emod=30e-9 nu=0.2 xcarea=0.2826 dens 2.5e3 &
xciy=6.36e-3 xciz=6.36e-3 &
per=1.884 &
cs_sk=24.75e6 cs_scoh=10e6 cs_sfric=15.0 &
cs_nk=24.75e6 cs_ncoh=10e6 cs_nfric=18.0 &
cs_ngap=off



fix x range x -0.1 0.1
fix x range x 49.9 50.1
fix y
fix x y z range z -0.1 0.1
set grav 0 0 -10
hist unbal

plot con szz





solve

sav 111初始应力.sav



ini xdisp 0.0 ydisp 0.0 zdisp 0.0
ini xvel 0.0 yvel 0.0 zvel 0.0
sel node ini xdisp=0.0 ydisp=0.0 zdisp=0.0
sel node ini xvel=0.0 yvel=0.0 zvel=0.0
ini state=0



prop coh 1 range z 0 8.4
prop coh 5e3 range z 8.4 17.4
prop coh 11e3 range z 17.4 25.6
prop coh 6e3 range z 25.6 27.5
prop coh 17e3 range z 27.5 31.2
prop coh 10e3 range z 31.2 33.4
prop coh 8e3 range z 33.4 35.2


sel pile prop emod=30e9 nu=0.2 xcarea=0.2826 dens 2.5e3 &
xciy=6.36e-3 xciz=6.36e-3 &
per=1.884 &
cs_sk=24.75e6 cs_scoh=15e3 cs_sfric=15.0 &
cs_nk=24.75e6 cs_ncoh=15e3 cs_nfric=18.0 &
cs_ngap=off



sel delete link range id 11
sel link id 200 2 target zone


sel link attach xrdir=free yrdir=free zrdir=free range id 200
sel link attach ydir=free zdir=free range id 200
sel link attach xdir=nydeform range id 200
sel link constit nydeform 1 area=1.0 k=5.4e11 ycomp=2.22e5 range id 200


plot surf
plot add sel geom id on cid on
plot hist 1
plot sel geom black black node off
plot add sel disp
plot add sel pile moment my

sav 113kaiwaqian.sav

model null range x 30 50 z 32 35.2
solve
sav 115.sav
model null range x 30 50 z 29 32
solve
sav 117end.sav

chengcai18

一个群的 不搞基坑 帮顶