Flac3D隧道 计算实例

;===============================================
;建模型
;===============================================
new
plot create model
plot set title text '模型图
plot set back white
plot add sketch
plot show

def x_rad
xrad=sqrt(6.19*6.19-2.19*2.19)
fxrad=-xrad
end
x_rad

gen zon cylinder p0 0 0 0 p1 xrad 0 -2.19 p2 add 0 1 0 p3 0 0 6.19 p4 xrad 1 -2.19 p5 0 1 6.19 size 7 1 10 rat 1 1 1

gen zon wedge p0 0 0 -2.19 p1 xrad 0 -2.19 p2 add 0 1 0 p3 0 0 0 p4 xrad 1 -2.19 p5 0 1 0 p6 0 0 0 p7 0 1 0 &
size 5 1 7 rat 1 1 1

gen zon reflect dip 90 dd 90 ori 0 0 0
group nsoil

gen zon radcyl p0 0 0 0 p1 36.19 0 -13.311021248234213556702313453022 p2 0 1 0 p3 0 0 21.19 p4 36.19 1 &
-13.311021248234213556702313453022 p5 0 1 21.19 &
p6 36.19 0 21.19 p7 36.19 1 21.19 &
p8 xrad 0 -2.19 p9 0 0 6.19 p10 xrad 1 -2.19 p11 0 1 6.19 dim 6.19 6.19 6.19 6.19 &
size 5 1 10 10 rat 1 1 1 1.2

gen zon brick p0 0 0 21.19 p1 36.19 0 21.19 p2 add 0 1 0 p3 0 0 46.19 size 5 1 5 rat 1 1 1

gen zon brick p0 36.19 0 21.19 p1 add 25 0 0 p2 add 0 1 0 p3 add 0 0 25 size 5 1 5 rat 1 1 1

gen zon brick p0 36.19 0 -13.311021248234213556702313453022 p1 add 25 0 0 p2 add 0 1 0 &
p3 36.19 0 21.19 size 5 1 5 rat 1 1 1

gen zon brick p0 0 0 -2.19 p1 0 0 -13.311021248234213556702313453022 p2 add 0 1 0 p3 add xrad 0 0 p4 0 1 -13.311021248234213556702313453022 &
p5 add xrad 1 0 p6 36.19 0 -13.311021248234213556702313453022 p7 36.19 1 -13.311021248234213556702313453022 &
size 10 1 5 rat 1.2 1 1

gen zon brick p0 0 0 -27.19 p1 add 36.19 0 0 p2 add 0 1 0 p3 add 0 0 13.87897875176578644329768654698 size 5 1 2 rat 1 1 1

gen zon brick p0 36.19 0 -27.19 p1 add 25 0 0 p2 add 0 1 0 p3 add 0 0 13.87897875176578644329768654698 size 5 1 2 rat 1 1 1
group wsoil range group nsoil not
gen zon reflect dip 90 dd 90 ori 0 0 0 range group wsoil
;pause
;===============================================
;初始应力计算程序
;===============================================
def derive
s_mod1=E_mod1/(2.0*(1.0+p_ratio1))
b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1))
end
set E_mod1=1.2e8 p_ratio1=0.43
derive

mo elas
prop bulk b_mod1 shear s_mod1 density 2000
ini szz -92380 grad 0 0 2000
ini syy -46190 grad 0 0 1000
ini sxx -46190 grad 0 0 1000

fix x range x 61.18 61.20
fix x range x -61.20 -61.18
fix y range y -0.1 0.1
fix y range y 0.9 1.1
fix x y z range z -27.20 -27.18
set grav 0 0 -10
;===============================================
;设置监测关键点
;===============================================
hist id=1 unbal

hist id=2 gp zdis 0 0 6.19
hist id=3 gp zdis 0 0 2.19
hist id=4 gp zdis xrad 0 2.19
hist id=5 gp zdis fxrad 0 2.19
hist id=6 gp zdis 6.19 0 0
hist id=7 gp zdis -6.19 0 0
hist id=8 gp zdis 4.37699 0 4.37699
hist id=9 gp zdis -4.37699 0 4.37699

hist id=12 zone smin 0 0 6.19
hist id=13 zone smin 0 0

2.19
hist id=14 zone smin xrad 0 2.19
hist id=15 zone smin fxrad 0 2.19
hist id=16 zone smin 6.19 0 0
hist id=17 zone smin -6.19 0 0
hist id=18 zone smin 4.37699 0 4.37699
hist id=19 zone smin -4.37699 0 4.37699

hist id=22 zone ssi 0 0 6.19
hist id=23 zone ssi 0 0 2.19
hist id=24 zone ssi xrad 0 2.19
hist id=25 zone ssi fxrad 0 2.19
hist id=26 zone ssi 6.19 0 0
hist id=27 zone ssi -6.19 0 0
hist id=28 zone ssi 4.37699 0 4.37699
hist id=29 zone ssi -4.37699 0 4.37699

;hist id=32 sel recover sres mx cid=19
;hist id=33 sel recover sres mx cid=20
;hist id=34 sel recover sres mx cid=21
;hist id=35 sel recover sres mx cid=22
;hist id=36 sel recover sres mx cid=3
;hist id=37 sel recover sres mx cid=4
;hist id=38 sel recover sres mx cid=6
;hist id=39 sel recover sres mx cid=7

;hist id=42 sel recover sres qx 0 0 6.19
;hist id=43 sel recover sres qx 0 0 2.19
;hist id=44 sel recover sres qx xrad 0 2.19
;hist id=45 sel recover sres qx fxrad 0 2.19
;hist id=46 sel recover sres qx 6.19 0 0
;hist id=47 sel recover sres qx -6.19 0 0
;hist id=48 sel recover sres qx 4.37699 0 4.37699
;hist id=49 sel recover sres qx -4.37699 0 4.37699

;===============================================
;绘结果图
;===============================================

plot create stress
plot set title text '应力云图
plot set back white
plot add con sxx outline on
plot add con szz outline on
plot add con sxz outline on

plot create strain
plot set title text '应变云图
plot set back white
plot add con ssr outline on average
plot add con ssi outline on average

plot create dis-step
plot set title text '监测点位移图
plot set back white
plot add his 2 3 4 5 6 7 8 9 vs step

plot create p-stress-step
plot set title text '监测点最大主应力
plot set back white
plot add his 12 13 14 15 16 17 18 19 vs step

plot create shear-strain-step
plot set title text '监测点剪应变
plot set back white
plot add his 22 23 24 25 26 27 28 29 vs step

;plot create moment-step
;plot set title text '监测点弯矩
;plot set back white
;plot add his 32 33 34 35 36 37 38 39 vs step

;plot create shear-step
;plot set title text '监测点剪力
;plot set back white
;plot add his 42 43 44 45 46 47 48 49 vs step

;===============================================
;求解
;===============================================
solve
save sd51.sav

restore sd51.sav
plot block state now
;--------------------------------------------------------
;定义初期支护锚秆
;--------------------------------------------------------
def fen_jiao_chu
h=9.0
;9=180/20
x11=6.19*cos((-h)*degrad)
z11=6.19*sin((-h)*degrad)
x12=9.19*cos((-h)*degrad)
z12=9.19*sin((-h)*degrad)
x21=6.19*cos((180+h)*degrad)
z21=6.19*sin((180+h)*degrad)
x22=9.19*cos((180+h)*degrad)
z22=9.19*sin((180+h)*degrad)
x31=xrad
z31=-2.19
x32=x31+3/(1.4142)
z32=z31-3/(1.4142)
x41=-xrad
z41=-2.19
x42=x41-3/(1.4142)
z42=z41

-3/(1.4142)
end
fen_jiao_chu

def cab_parm_chu
array x2_b(1,54),z2_b(1,54),x2_2b(1,54),z2_2b(1,54)
loop n(1,21)
x2_b(1,n)=6.19*cos(h*(n-1)*degrad)
z2_b(1,n)=6.19*sin(h*(n-1)*degrad)
x2_2b(1,n)=9.19*cos(h*(n-1)*degrad)
z2_2b(1,n)=9.19*sin(h*(n-1)*degrad)
end_loop
end
cab_parm_chu

def inip_chu
x_d=x2_b(1,n)
z_d=z2_b(1,n)
x_d_2=x2_2b(1,n)
z_d_2=z2_2b(1,n)
end

def ini_chu_mao
loop m(3,23)
n=m-2
inip_chu
command
sel cable id m begin x_d,0.5,z_d end x_d_2,0.5, z_d_2 nseg 10
sel cable pro emod 2.1e10 xcarea 3.8e-4 &
yten 2.5e5 gr_k 1.75e7 gr_c 2e5 range id m
end_command
end_loop
end


;--------------------------------------------------------
;开挖------安装初期支护锚杆
;--------------------------------------------------------
;开挖
ini xdis=0 ydis=0 zdis=0
def derive1
s_mod3=E_mod3/(2.0*(1.0+p_ratio3))
b_mod3=E_mod3/(3.0*(1.0-2.0*p_ratio3))
end
set E_mod3=1.2e8 p_ratio3=0.43
derive1

m m
pro bulk b_mod3 shear s_mod3 density 2000 fric 22 coh 2.5e4 dil 30
m n range group nsoil
step 200
;pause
;--------------------------------------------------------
;安装初期支护锚杆和喷混凝土
ini_chu_mao
sel cable id 2 begin x11,0.5,z11 end x12,0.5, z12 nseg 10
sel cable pro emod 2.1e10 xcarea 3.79e-4 gr_per 1. &
yten 2.2e5 gr_k 17.5e6 gr_c 20e4 range id 2
sel cable id 1 begin x31,0.5,z31 end x32,0.5, z32 nseg 10
sel cable pro emod 2.1e10 xcarea 3.79e-4 gr_per 1. &
yten 2.2e5 gr_k 17.5e6 gr_c 20e4 range id 1
sel cable id 24 begin x21,0.5,z21 end x22,0.5, z22 nseg 10
sel cable pro emod 2.1e10 xcarea 3.79e-4 gr_per 1. &
yten 2.2e5 gr_k 17.5e6 gr_c 20e4 range id 24
sel cable id 25 begin x41,0.5,z41 end x42,0.5, z42 nseg 10
sel cable pro emod 2.1e10 xcarea 3.79e-4 gr_per 1. &
yten 2.2e5 gr_k 17.5e6 gr_c 20e4 range id 25

sel shell id 55 range cyl end1 0 0 0 end2 0 1 0 rad 6.19 z -2.20 6.19
sel shell prop isotropic 3.0e10,0.2 thickness 0.35 density 2400 range id 55
;pause

solve
plot add sel recover dis
save sdcable1.sav