-------
;Flac 锚杆模拟
;------------------
new
gen zone radcyl size 6 5 12 8 p0 (0 0 0) p1=(15 0 0) p2=(0 15 0) p3=(0 0 15) p4=(15 15 0) p5=(0 15 15) p6=(15 0 15) p7=(15 15 15) p8=(3 0 0) p9=(0 0 3) p10=(3 15 0) p11=(0 15 3) ratio 1 1 1 1.12 fill
gen zone reflect dip 90 dd 90 orig =(0 0 0)
gen zone reflect dip 0 dd 0 origin =(0 0 0)
pause
;------------------------- 岩土参数 ----------------------
DEF pro_set
young_m = 20e9
pos_f = 0.2
tiji_m = young_m/(3*(1-2*pos_f))
jian_m = young_m/(2*(1+pos_f))
END
pro_set
model mohr
property bulk tiji_m shear jian_m friction 35 cohesion 2.5e6 tension 1.5e6
set mech ratio 1e-4
; ----------------------- 施加边界条件 ----------------
ini sxx -15.3e6
ini syy -10.6e6
ini szz -5.0e6
fix x range x -15.2 -14.8
fix x range x 14.8 15.2
fix y range y -0.1 0.1
fix y range y 14.9 15.1
apply szz -5.0e6 range z 14.8 15.2
fix z range z -15.2 -14.8
step 1
;===塑性区清零====
DEF state_ini
xdisp=0.0
ydisp=0.0
zdisp=0.0
xvel=0
yvel=0
zvel=0
pnt=zone_head
loop while pnt # null
z_state(pnt,0)=0
pnt=z_next(pnt)
endloop
END
state_ini
group excavation range cylinder end1 (0 0 0) end2 (0 15 0) rad=3
model null range group excavation
sel shell range cyl end1 (0 0 0) end2 (0 15 0) rad 3
sel shell prop iso=(25.3e9 0.3) thick=0.5
solve
save shell.sav
pl sur
