gen z radtun p0 0 0 0 p1 @bb 0 0 p2 0 0 @modeldi p3 0 @bb 0 p4 @bb 0 @modeldi p5 0 @bb @modeldi ...
p6 @bb @bb 0 p7 @bb @bb @modeldi p8 @aa 0 0 p9 0 @aa 0 p10 @aa 0 @modeldi p11 0 @aa @modeldi ...
p12 @aa @aa 0 p13 @aa @aa @modeldi size 5 20 5 3 r 1 1 1 1.4
gro clay ran z 0 -15
gro sand ran z -15 -20
gen z refl dip 90 dd 90
gen z refl nor 0 -1 0 ori 0 0 0
;?¨???£¤??
def int_pra
pos_s=-5.9
pos_x=-6.1
rad_=0.265
radn_=0.065
end
@int_pra
int 1 face ran cylin end1 0 0 0 end2 0 0 @pos_x radius @rad_ ...
cylin end1 0 0 0 end2 0 0 @pos_x radius @radn_ not
int 2 face ran cylin end1 0 0 @pos_s end2 0 0 @pos_x radius @rad_
;??§¨¢?é?á?
gen z cylin p0 0 0 7.15 p1 @rad 0 7.15 p2 0 0 1 p3 0 @rad 7.15 p4 @rad 0 1 p5 0 @rad 1 ...
size 5 15 6 r 1 1 1 gro AA
gen z refl dip 90 dd 90 ran gro AA
gen z refl nor 0 -1 0 ori 0 0 0 ran gro AA
gro pile1 ran z 1 7.15
ini z add -7 ran gro pile1
attach face
;??§¨¢?é3D¨??§
gen z radcyl p0 0 0 0.15 p1 1 0 0.15 p2 0 0 0 p3 0 1 0.15 p4 1 0 0 ...
p5 0 1 0 p6 1 1 0.15 p7 1 1 0 p8 @rad 0 0.15 p9 0 @rad 0.15 ...
p10 @rad 0 0 p11 0 @rad 0 size 5 2 6 3 r 1 1 1 1.14 gro cap1
gen z reflect dip 90 dd 0 ran gro cap1
gen z reflect normal -1 0 0 ori 0 0 0.15 ran gro cap1
int 1 face ran cylin end1 0 0 0.15 end2 0 0 0 radius @rad_ ...
cylin end1 0 0 0.15 end2 0 0 0 radius @radn_ not
int 2 face ran cylin end1 0 0 0.1 end2 0 0 -0.1 radius @rad_
int 3 face ran x -1 1 y -1 1 z 0.15
gen z b p0 0 0 0.3 p1 1 0 0.3 p2 0 0 0.15 p3 0 1 0.3 p4 1 0 0.15 ...
p5 0 1 0.15 p6 1 1 0.3 p7 1 1 0.15 size 5 2 5 r 1 1 1 gro cap2
gen z refl dip 90 dd 90 ran gro cap2
gen z refl nor 0 -1 0 ori 0 0 0.3 ran gro cap2
gro cap ran gro cap1
gro cap ran gro cap2
;material property
model mohr ran gro clay
model mohr ran gro sand
model elas ran gro pile1
model elas ran gro cap
prop bulk 4.1667e7 shear 3e6 coh 3000 fric 35 range group clay
prop bulk 4.1667e7 shear 3e6 coh 3000 fric 45 range group sand
prop bulk 3e8 shear 3e8 range group pile1
prop bulk 3e8 shear 3e8 range group cap
ini dens 2200 ran gro clay
ini dens 1800 ran gro sand
ini dens 2800 ran gro pile1
ini dens 2800 ran gro cap
;k0=u/(1-u)
ini szz 0 grad 0 0 17000 ran z 0 -20
ini sxx 0 grad 0 0 7293 ran z 0 -20
ini syy 0 grad 0 0 7293 ran z 0 -20
;
fix x ran x -8
fix x ran x 8
fix y ran y -8
fix y ran y 8
fix z ran z -20
water density 1000
set grav 0 0 -9.8
water table ori 0 0 0 nor 0 0 -1
set mech ratio 1e-5
hist id 1 unbal
solve
;
ini state 0
ini xdis 0 ydis 0 zdis 0 xv 0 yv 0 zv 0
model fl_iso
set flu on
pro perm 1e-13 poros 0.5 biot_c 1.0 ran gro clay
pro perm 1e-17 poros 0.5 biot_c 1.0 ran gro sand
ini pp 1.0e4 grad 0 0 -1.0e4
ini fmod 2e9 fdens 1000
set flu pcut on
fix pp 0 range z -0.1 0.1
fix pp 0 ran x -8
fix pp 0 ran x 8
fix pp 0 ran y -8
fix pp 0 ran y 8
;
ini state 0
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
set dyn on flu off
app nquiet squiet ran z 0.3
table 1 read dynload.txt
app szz -2e5 hist table 1 ran z 0.3
table 2 read dynload1.txt
app syz -2e5 hist table 2 ran z 0.3
ini pp 30e3 grad 0 0 -10e3 ran z -20
app nquiet dquiet squiet ran z -20
set dyn damp rayleigh 1 22.8
;
ini state 0
ini xdis 0 ydis 0 zdis 0 xv 0 yv 0 zv 0
config thermal
set thermal on
model th_iso
set thermal off
pro conduct 1.8 spec_heat 1000 ran gro clay
pro conduct 1.21 spec_heat 1406 ran gro sand
pro conduct 1.74 spec_heat 920 ran gro pile1
pro conduct 1.74 spec_heat 920 ran gro cap
def th_add
loop n(15,30)
ini_th=5+n
command
ini t @ini_th ran gro pile1
endcommand
endloop
end
@th_add
ini t 5 ran gro clay
ini t 5 ran gro sand
app vsource 30 ran gro pile1
app vsource 30 interior ran gro pile1