flac3d一个多层土的例子

new
set logfile cz1000
set log on
;1# 桩长28.5米,桩直径1.0米
;-------------------------------------------------
; pile in rock for example
;-------------------------------------------------
gen zone radcyl p0 (0 0 0) p1 (8 0 0) p2 (0 0 -2.26) p3 (0 8 0) &
p4 (8 0 -2.26) p5 (0 8 -2.26) p6 (8 8 0) p7 (8 8 -2.26) &
p8 (0.5 0 0) p9 (0 0.5 0) p10 (0.5 0 -2.26) p11 (0 0.5 -2.26) &
size 3 2 6 15 ratio 1 1 1 1.15 ;细砂
group fines range z -2.26 0
gen zone radcyl p0 (0 0 -2.26) p1 (8 0 -2.26) p2 (0 0 -8.46) p3 (0 8 -2.26) &
p4 (8 0 -8.46) p5 (0 8 -8.46) p6 (8 8 -2.26) p7 (8 8 -8.46) &
p8 (0.5 0 -2.26) p9 (0 0.5 -2.26) p10 (0.5 0 -8.46) p11 (0 0.5 -8.46) &
size 3 7 6 15 ratio 1 1 1 1.15 ;淤泥质土
group siltl range z -8.46 -2.26
gen zone radcyl p0 (0 0 -8.46) p1 (8 0 -8.46) p2 (0 0 -9.76) p3 (0 8 -8.46) &
p4 (8 0 -9.76) p5 (0 8 -9.76) p6 (8 8 -8.46) p7 (8 8 -9.76) &
p8 (0.5 0 -8.46) p9 (0 0.5 -8.46) p10 (0.5 0 -9.76) p11 (0 0.5 -9.76) &
size 3 2 6 15 ratio 1 1 1 1.15 ;粉砂
group fengs range z -9.76 -8.46
gen zone radcyl p0 (0 0 -9.76) p1 (8 0 -9.76) p2 (0 0 -15.36) p3 (0 8 -9.76) &
p4 (8 0 -15.36) p5 (0 8 -15.36) p6 (8 8 -9.76) p7 (8 8 -15.36) &
p8 (0.5 0 -9.76) p9 (0 0.5 -9.76) p10 (0.5 0 -15.36) p11 (0 0.5 -15.36) &
size 3 7 6 15 ratio 1 1 1 1.15 ;中砂
group middles range z -15.36 -9.76
gen zone radcyl p0 (0 0 -15.36) p1 (8 0 -15.36) p2 (0 0 -19.66) p3 (0 8 -15.36) &
p4 (8 0 -19.66) p5 (0 8 -19.66) p6 (8 8 -15.36) p7 (8 8 -19.66) &
p8 (0.5 0 -15.36) p9 (0 0.5 -15.36) p10 (0.5 0 -19.66) p11 (0 0.5 -19.66) &
size 3 6 6 15 ratio 1 1 1 1.15 ;砂质粘性土
group cohs range z -19.66 -15.36
gen zone radcyl p0 (0 0 -19.66) p1 (8 0 -19.66) p2 (0 0 -23.66) p3 (0 8 -19.66) &
p4 (8 0 -23.66) p5 (0 8 -23.66) p6 (8 8 -19.66) p7 (8 8 -23.66) &
p8 (0.5 0 -19.66) p9 (0 0.5 -19.66) p10 (0.5 0 -23.66) p11 (0 0.5 -23.66) &
size 3 5 6 15 ratio 1 1 1 1.15 ;全风化花岗岩
group quanh range z -23.66 -19.66
gen zone radcyl p0 (0 0 -23.66) p1 (8 0 -23.66) p2 (0 0 -27.46) p3 (0 8 -23.66) &
p4 (8 0 -27.46) p5 (0 8 -27.46) p6 (8 8 -23.66) p7 (8 8 -27.46) &
p8 (0.5 0 -23.66) p9 (0 0.5 -23.66) p10 (0.5 0 -27.46) p11 (0 0.5 -27.46) &
size 3 5 6 15 ratio 1 1 1 1.15 ;强风化花岗岩
group qiangh range z -27.46 -23.66
gen zone radcyl p0 (0 0 -27.46) p1 (8 0 -27.46) p2 (0 0 -28.50) p3 (0 8 -27.46) &
p4 (8 0 -28.50) p5 (0 8 -28.50) p6 (8 8 -27.46) p7 (8 8 -28.50) &
p8 (0.5 0 -27.46) p9 (0 0.5 -27.46) p10 (0.5 0 -28.50) p11 (0 0.5 -28.50) &
size 3 2 6 15 ratio 1 1 1 1.15 ;中风化花岗岩,以上为有桩体部分的岩土体
gen zone radcyl p0 (0 0 -28.50) p1 (8 0 -28.50) p2 (0 0 -36.50) p3 (0 8 -28.50) &
p4 (8 0 -36.50) p5 (0 8 -36.50) p6 (8 8 -28.50) p7 (8 8 -36.50) &
p8 (0.5 0 -28.50) p9 (0 0.5 -28.50) p10 (0.5 0 -36.50) p11 (0 0.5 -36.50) &
size 3 10 6 15 ratio 1 1.15 1 1.15 fill ;中风化砾岩,此部分为无桩部分的岩体
group zhongh range z -36.50 -27.46
gen zone reflect dd 270 dip 90
;
interface 1

face range cylin end1 (0 0 0) end2 (0 0 -2.26) radius 0.51
interface 2 face range cylin end1 (0 0 -2.26) end2 (0 0 -8.46) radius 0.51
interface 3 face range cylin end1 (0 0 -8.46) end2 (0 0 -9.76) radius 0.51
interface 4 face range cylin end1 (0 0 -9.76) end2 (0 0 -15.36) radius 0.51
interface 5 face range cylin end1 (0 0 -15.36) end2 (0 0 -19.66) radius 0.51
interface 6 face range cylin end1 (0 0 -19.66) end2 (0 0 -23.66) radius 0.51
interface 7 face range cylin end1 (0 0 -23.66) end2 (0 0 -27.46) radius 0.51
interface 8 face range cylin end1 (0 0 -27.46) end2 (0 0 -28.60) radius 0.51 &
cylin end1 (0 0 -27.46) end2 (0 0 -28.60) radius 0.49 not
interface 9 face range cylin end1 (0 0 -28.40) end2 (0 0 -28.60) radius 0.51
;
gen zone cyl p0 (0 0 29.50) p1 (0.5 0 29.50) p2 (0 0 27.24) p3 (0 0.5 29.50) &
p4 (0.5 0 27.24) p5 (0 0.5 27.24) &
size 3 2 6
gen zone cyl p0 (0 0 27.24) p1 (0.5 0 27.24) p2 (0 0 21.04) p3 (0 0.5 27.24) &
p4 (0.5 0 21.04) p5 (0 0.5 21.04) &
size 3 7 6
gen zone cyl p0 (0 0 21.04) p1 (0.5 0 21.04) p2 (0 0 19.74) p3 (0 0.5 21.04) &
p4 (0.5 0 19.74) p5 (0 0.5 19.74) &
size 3 2 6
gen zone cyl p0 (0 0 19.74) p1 (0.5 0 19.74) p2 (0 0 14.14) p3 (0 0.5 19.74) &
p4 (0.5 0 14.14) p5 (0 0.5 14.14) &
size 3 7 6
gen zone cyl p0 (0 0 14.14) p1 (0.5 0 14.14) p2 (0 0 9.84) p3 (0 0.5 14.14) &
p4 (0.5 0 9.84) p5 (0 0.5 9.84) &
size 3 6 6
gen zone cyl p0 (0 0 9.84) p1 (0.5 0 9.84) p2 (0 0 5.84) p3 (0 0.5 9.84) &
p4 (0.5 0 5.84) p5 (0 0.5 5.84) &
size 3 5 6
gen zone cyl p0 (0 0 5.84) p1 (0.5 0 5.84) p2 (0 0 2.04) p3 (0 0.5 5.84) &
p4 (0.5 0 2.04) p5 (0 0.5 2.04) &
size 3 5 6
gen zone cyl p0 (0 0 2.04) p1 (0.5 0 2.04) p2 (0 0 1.00) p3 (0 0.5 2.04) &
p4 (0.5 0 1.00) p5 (0 0.5 1.00) &
size 3 2 6
gen zone cyl p0 (0 0 29.60) p1 (0.5 0 29.60) p2 (0 0 29.50) p3 (0 0.5 29.60) &
p4 (0.5 0 29.50) p5 (0 0.5 29.50) &
size 3 1 6
gen zone reflect dd 270 dip 90 range z 1 29.60
group pile range z 1 29.60
plo sur
pause
ini z add -29.50 range z 1 29.60
; material property
;------------------------------------------------------------------------------------
macro finesand 'bulk 5.667e6 shear 3.400e6 coh 0 fric 30' ;细砂
macro silt 'bulk 3.100e6 shear 1.033e6 coh 5500 fric 8.0' ;淤泥质土
macro fengsand 'bulk 4.333e6 shear 2.600e6 coh 0 fric 25' ;粉砂
macro middlesand 'bulk 8.867e6 shear 5.32e6 coh 0 fric 34.7' ;中砂
macro cohsoil 'bulk 3.647e6 shear 2.188e6 coh 18500 fric 23.0' ;砂质粘性土
macro quanfhhgy 'bulk 4.427e6 shear 3.602e6 coh 22000 fric 25' ;全风化花岗岩
macro qiangfhhgy 'bulk 7.143e6 shear 6.522e6 coh 24000 fric 27.5' ;强风化花岗岩
macro zhongfhhgy 'bulk 7.046e6 kshear 3.0e4 qdil 1.0 qvol 0.7 shear 7.686e6 tension 4.25e5' ;中风化花岗岩
;
;细砂
model mo range group fines
prop finesand range group fines
;细砂
;------------------------------------------------------------------------------------
;淤泥质土
model mo range group sil

tl
prop silt range group siltl
;淤泥质土
;------------------------------------------------------------------------------------
;粉砂
model mo range group fengs ;z -9.76 -8.46
prop fengsand range group fengs
;粉砂
;------------------------------------------------------------------------------------
;------------------------------------------------------------------------------------
;中砂
model mo range group middles ;z -15.36 -9.76
prop middlesand range group middles
;中砂
;------------------------------------------------------------------------------------
;砾质粘性土
model mo range group cohs ;z -19.66 -15.36
prop cohsoil range group cohs
;砾质粘性土
;------------------------------------------------------------------------------------
;全风化花岗岩
model mo range group quanh ;z -23.66 -19.66
prop quanfhhgy range group quanh
;全风化花岗岩
;------------------------------------------------------------------------------------
;强风化花岗岩
model mo range group qiangh ;z -27.46 -23.66
prop qiangfhhgy range group qiangh
;强风化花岗岩
;------------------------------------------------------------------------------------
;中风化花岗岩
model drucker range group zhongh ;z -36.50 -27.46
prop zhongfhhgy range group zhongh
;中风化花岗岩
;------------------------------------------------------------------------------------
;
model elas range z -2.26 0.1 group pile
prop bulk 5.667e6 shear 3.400e6 range z -2.26 0.0 group pile
model elas range z -8.46 -2.26 group pile
prop bulk 3.100e6 shear 1.033e6 range z -8.46 -2.26 group pile
model elas range z -9.76 -8.46 group pile
prop bulk 4.333e6 shear 2.600e6 range z -9.76 -8.46 group pile
model elas range z -15.36 -9.76 group pile
prop bulk 8.867e6 shear 5.32e6 range z -15.36 -9.76 group pile
model elas range z -19.66 -15.36 group pile
prop bulk 3.647e6 shear 2.188e6 range z -19.66 -15.36 group pile
model elas range z -23.66 -19.66 group pile
prop bulk 4.427e6 shear 3.602e6 range z -23.66 -19.66 group pile
model elas range z -27.46 -23.66 group pile
prop bulk 7.143e6 shear 6.522e6 range z -27.46 -23.66 group pile
model elas range z -28.50 -27.46 group pile
prop bulk 7.046e6 shear 7.686e6 range z -28.50 -27.46 group pile
;
;pause
;interface property,P5-86
interface 1 prop kn 3e6 ks 1e8 fric 20 coh 30000
interface 2 prop kn 2e5 ks 1e8 fric 20 coh 30000
interface 3 prop kn 5e6 ks 1e8 fric 20 coh 30000
interface 4 prop kn 6e6 ks 1e8 fric 20 coh 30000
interface 5 prop kn 1e7 ks 1e8 fric 20 coh 30000
interface 6 prop kn 4e7 ks 1e8 fric 20 coh 30000
interface 7 prop kn 3e7 ks 1e8 fric 20 coh 30000
interface 8 prop kn 4e7 ks 1e8 fric 20 coh 30000
interface 9 prop kn 4e7 ks 1e8 fric 20 coh 30000
;------------------------------------------------------------------------------------
;
ini dens 1980 range z -2.26 0 ;细砂
ini dens 1710 range z -8.46 -2.26 ;淤泥质土
ini dens 1750 range z -9.76

-8.46 ;粉砂
ini dens 2000 range z -15.36 -9.76 ;中砂
ini dens 1930 range z -19.66 -15.36 ;砂质粘性土
ini dens 1940 range z -23.66 -19.66 ;全风化花岗岩
ini dens 1970 range z -27.46 -23.66 ;强风化花岗岩
ini dens 1990 range z -36.50 -27.46 ;中风化花岗岩
;
ini dens 1980 range z -2.26 0.1 group pile ;细砂
ini dens 1710 range z -8.46 -2.26 group pile ;淤泥质土
ini dens 1750 range z -9.76 -8.46 group pile ;粉砂
ini dens 2000 range z -15.36 -9.76 group pile ;中砂
ini dens 1930 range z -19.66 -15.36 group pile ;砂质粘性土
ini dens 1940 range z -23.66 -19.66 group pile ;全风化花岗岩
ini dens 1970 range z -27.46 -23.66 group pile ;强风化花岗岩
ini dens 1990 range z -28.50 -27.46 group pile ;中风化花岗岩
;
model null range z -0.1 0.15
fix z range z -8.1 -7.9
fix x range x -8.1 -7.9
fix x range x 7.9 8.1
fix y range y -0.1 0.1
fix y range y 7.9 8.1
set grav 0 0 -10
;---------------------------------------------------------
;细砂
ini szz 0 grad 0 0 19800 range z -2.26 0
ini syy 0 grad 0 0 6600 range z -2.26 0
ini sxx 0 grad 0 0 6600 range z -2.26 0
;---------------------------------------------------------
;淤泥质土
ini szz -6102 grad 0 0 17100 range z -8.46 -2.26
ini syy -3285.69 grad 0 0 9207.69 range z -8.46 -2.26
ini sxx -3285.69 grad 0 0 9207.69 range z -8.46 -2.26
;---------------------------------------------------------
;粉砂
ini szz -2718 grad 0 0 17500 range z -9.76 -8.46
ini syy -906 grad 0 0 5833.33 range z -9.76 -8.46
ini sxx -906 grad 0 0 5833.33 range z -9.76 -8.46
;---------------------------------------------------------
;中砂
ini szz 21682 grad 0 0 20000 range z -15.36 -9.76
ini syy 7227.33 grad 0 0 6666.67 range z -15.36 -9.76
ini sxx 7227.33 grad 0 0 6666.67 range z -15.36 -9.76
;---------------------------------------------------------
;砾质粘性土
ini szz 10930 grad 0 0 19300 range z -19.66 -15.36
ini syy 3643.33 grad 0 0 6433.33 range z -19.66 -15.36
ini sxx 3643.33 grad 0 0 6433.33 range z -19.66 -15.36
;---------------------------------------------------------
;全风化花岗岩
ini szz 12896 grad 0 0 19400 range z -23.66 -19.66
ini syy 2830.83 grad 0 0 4258.54 range z -23.66 -19.66
ini sxx 2830.83 grad 0 0 4258.54 range z -23.66 -19.66
;---------------------------------------------------------
;强风化花岗岩
ini szz 19994 grad 0 0 19700 range z -27.46 -23.66
ini syy 3528.38 grad 0 0 3476.47 range z -27.46 -23.66
ini sxx 3528.38 grad 0 0 3476.47 range z -27.46 -23.66
;---------------------------------------------------------
;中风化花岗岩
ini szz 22740 grad 0 0 19800 range z -36.50 -27.46
ini syy 2526.67 grad 0 0 2200 range z -36.50 -27.46
ini sxx 2526.67 grad 0 0 2200 range z -36.50 -27.46
;---------------------------------------------------------
hist unbal
;
;pause
solve rat 1.0e-6
save pile0.sav

;pause
;
model elas range group pile
prop bulk 13.9e9 shear 10.4e9 range grou

p pile
ini dens 2500 range group pile
;pause
;plot block density
;
;pause
call find_add.fis
solve rat 1.0e-6
save pile1.sav
;
ini state 0
ini xdis 0 ydis 0 zdis 0
;
; monitor vertical loading at pile cap
def zs_top
ad = top_head
zftot = 0.0
loop while ad # null
gp_pnt = mem(ad+1)
zf = gp_zfunbal(gp_pnt)
zftot = zftot + zf
ad = mem(ad)
endloop
zs_top = zftot / 0.1414
end
;pause
fix z range z 0.05 0.15 group pile
;
def ramp
while_stepping
if step < ncut then
udapp = float(step) * udmax / float(ncut)
else
udapp = udmax
endif
ad = top_head
loop while ad # null
gp_pnt = mem(ad+1)
gp_zvel(gp_pnt) = udapp
ad = mem(ad)
endloop
end
;
hist gp zdis 0 0 0
hist gp zvel 0 0 0
hist zs_top
hist zone szz 0 0 -0.1
;
set mech damp comb
set udmax=-1e-8 ncut 30000
step 225000
save pile2.sav
;
rest pile1.sav
;
apply szz -7.85e6 range z 0.05 0.15 group pile
;200t,1.57e6;300t,2.355e6;400t,3.14e6;500t,3.925e6;600t,4.71;
;700t,5.495e6;800t,6.28e6;900t,7.065e6;1000t,7.85e6(Pa)
ini state 0
ini xdis 0.0 ydis 0.0 zdis 0.0
;
hist gp zdis 0 0 0
hist gp zdis 0 0 -28.50
;
solve rat 1.e-6
;
save pilef2.sav
;-------------------------------------------------------------------------------------------------------------------------
;
;
;
;
ini xdis 0.0 ydis 0.0 zdis 0.0
;
; monitor horizontal loading along pile shaft
def xf_p1
ad = top_head
xftot = 0.0
loop while ad # null
gp_pnt = mem(ad+1)
xf = gp_xfunbal(gp_pnt)
xftot = xftot + xf
ad = mem(ad)
endloop
xf_p1 = xftot
end
call p-y.fis
call servo.fis
set vel_min 5e-8 vel_max 1e-6 unb_min 10 unb_max 50 horapp 5e-8
fix z range z 0.05 0.15 group pile
;
hist gp xdis 0,0,0
hist gp xvel 0,0,0
hist horapp
hist tot_reac
hist xf_p1
;
step 60000
save pile2h.sav
ret