gen zone radcyl p0 0 0 16 p1 6 0 16 p2 0 25 16 p3 0 0 22 p4 6 25 16 p5 0 25 22 &
p6 6 0 22 p7 6 25 22 p8 4 0 16 p9 0 0 20 p10 4 25 16 p11 0 25 20 &
size 4 14 16 4 ratio 1 1 1 1.2 fill group 3
gen zone reflect nor 0 0 15 ori 0 0 16
gen zone brick p0 6 0 22 p1 40 0 22 p2 6 25 22 p3 6 0 42 size 10 14 10 ratio 1.2 1 1 group 2
gen zone brick p0 0 0 22 p1 6 0 22 p2 0 25 22 p3 0 0 42 size 8 14 10 group 1
gen zone brick p0 6 0 10 p1 40 0 10 p2 6 25 10 p3 6 0 22 size 10 14 10 ratio 1.2 1 1
gen zone brick p0 0 0 0 p1 6 0 0 p2 0 25 0 p3 0 0 10 size 8 14 10
gen zone brick p0 6 0 0 p1 40 0 0 p2 6 25 0 p3 6 0 10 size 10 14 10 ratio 1.2 1 1
group 4 range z 20 22
model mech mohr
;====================================
def paramer
global nu_1=0.35 ;can ji tu
global E_1=1e8
global nu_2=0.25 ;wei feng hua yan
global E_2=3e8
global nu_pile=0.2 ;pile
global E_3=3e10
global nu_liner=0.2 ;liner
global E_4=5.5e9
global bulk_1=E_1/(3.*(1.-2.*nu_1))
global bulk_2=E_2/(3.*(1.-2.*nu_2))
global bulk_3=E_pile/(3.*(1.-2.*nu_pile))
global bulk_4=E_liner/(3.*(1.-2.*nu_liner))
global G_1=E_1/(2.*(1.+1.*nu_1))
global G_2=E_2/(2.*(1.+1.*nu_2))
global G_3=E_pile/(2.*(1.+1.*nu_pile))
global G_4=E_liner/(2.*(1.+1.*nu_liner))
end
@paramer
;===========================================
prop coh 1e5 fric 35 bulk @bulk_1 shear @G_1 tens 3e5 range z 22 42
prop coh 1.5e5 fric 35 bulk @bulk_2 shear @G_2 tens 4.5e3 range z 0 22
ini dens 2
fix x range x -.1 .1
fix x range x 39.9 40.1
fix y range y -.1 .1
fix y range y -24.9 25.1
fix z range z -.1 .1
ini szz 0 grad 0 0 20 range z 0 42
set grav 0 0 -9.81
set small
hist unbal
solve
save hefei_1.sav
;==============================================
ini xdisp 0 ydisp 0 zdisp 0
group pile range x 0 1.6 y 0 1.8 z 22 42
gen separate range group pile
int 1 wrap pile 1
int 2 wrap pile 4
interface 1 prop coh 1e5 kn 1e8 ks 1e8
interface 1 prop coh 1e5 kn 1e8 ks 1e8
ini dens 2.5 range group pile
;===========================================
model mech null range y 0 2 group 3
hist gp zdisp 0 0 42
hist gp zdisp 5 0 42
set large
solve
save hefei_2.sav
;========================================
model mech null range y 2 4 group 3
hist gp zdisp 0 0 42
hist gp zdisp 5 0 42
set large
solve
save hefei_3
sel shell id=1 range cyl end1 0 0 16 end2 0 2.01 16 rad 4
;sel node fix x yr zr range x -.1 .1
;sel node fix y xr zr range y -.1 .1
;sel node fix z xr yr range z -.1 .1
sel shell prop isot 5500 0.2 thick 0.2 dens 2.5
solve
save hefei_4