www.pudn.com > coal-excavate.zip > 巷道已支护.txt, change:2013-05-31,size:4361b


New 
 
;园柱体外环绕放射状网格 
Generate zone radcylinder size 5 15 10 5 ratio 1 1 1 1.1 p0=(0,0,0) p1=(5,0,0) p2=(0,30,0) p3=(0,0,5) dimension 2.7,2.7,2.7 group 1 
 
;平行六面体外环绕放射状网格 
Generate zone radtunnel size 5 15 5 5 ratio 1 1 1 1.1 p0=(0,0,0) p1=(0,0,-5) p2=(0,30,0) p3=(5,0,0) dimension 2.7,2.7,2.7 group 2 
 
;两侧岩体 
gen zone b p0 5 0 -5 p1 15 0 -5 p2 5 30 -5 p3 5 0 5 size 10 15 10 group 3 
 
;上部岩体 
gen zone b p0 0 0 5 p1 15 0 5 p2 0 30 5 p3 0 0 15 size 15 15 5 group 4 
 
;下部岩体 
gen zone b p0 0 0 -15 p1 15 0 -15 p2 0 30 -15 p3 0 0 -5 size 15 15 5 group 5 
 
;镜像网格 
Generate zone reflect normal (-1,0,0) origin=(0,0,0) 
 
;设定摩尔库伦材料参数 
model mohr 
 
;材料参数赋值 
pro bulk 2.29e9 she 1.65e9 fric 28 coh 3.3e6 ten 0.93e6 range z -15 -13.16 ;3煤 
pro bulk 11.5e9 she 9.84e9 fric 24.4 coh 5e6 ten 2.79e6 range z -13.16 -9.76 ;细砂岩 
pro bulk 6.25e9 she 3.5e9 fric 27 coh 5.4e6 ten 3.0e6 range z -9.76 -5 ;泥岩 
pro bulk 12.1e9 she 10.86e9 fric 32 coh 10.27e6 ten 25.7e6 range z -5 5.51 ;中砂岩(巷道所在层) 
pro bulk 6.12e9 she 3.54e9 fric 27 coh 5.12e6 ten 2.4e6 range z 5.51 6.87 ;泥岩 
pro bulk 2.29e9 she 1.65e9 fric 28 coh 3.4e6 ten 0.94e6 range z 6.87 7.32 ;2煤 
pro bulk 3.3e9 she 3.14e9 fric 27 coh 5.1e6 ten 2.2e6 range z 7.32 10.94 ;泥岩 
pro bulk 11.5e9 she 9.84e9 fric 25 coh 4.7e6 ten 2.79e6 range z 10.94 15 ;细砂岩 
 
;安装锚索 
def maosuo 
suo=-0.4 
loop while suo<30 
suo1=suo+0.7 
command 
sel cable begin 2.65 suo1 0 end 10.65 suo1 0 nseg 15 
sel cable begin 1.88 suo1 1.88 end 7.53 suo1 7.53 nesg 15 
sel cable begin -2.65 suo1 0 end -10.65 suo1 0 nseg  15 
sel cable begin -1.88 suo1 1.88 end -7.53 suo1 7.53 nesg 15 
sel cable begin 0 suo1 2.65 end 0 suo1 10.65 nseg  15 
end_command 
suo=suo+0.7 
end_loop 
end 
maosuo 
 
;安装锚杆 
def maogan 
gan=-0.4 
loop while gan<30 
gan1=gan+0.7 
;底部 
command 
sel cable begin 1.8 gan1 -2.7 end 1.8 gan1 -5.2 nseg 15 
sel cable begin 0.9 gan1 -2.7 end 0.9 gan1 -5.2 nseg 15 
sel cable begin 0 gan1 -2.7 end 0 gan1 -5.2 nseg 15 
sel cable begin -0.9 gan1 -2.7 end -0.9 gan1 -5.2 nseg 15 
sel cable begin -1.8 gan1 -2.7 end -1.8 gan1 -5.2 nseg 15 
;帮 
sel cable begin 2.65 gan1 -1.8 end 5.2 gan1 -1.8 nseg 15 
sel cable begin 2.65 gan1 -1.1 end 5.2 gan1 -1.1 nseg 15 
sel cable begin 2.65 gan1 -0.4 end 5.2 gan1 -0.4 nseg 15 
sel cable begin -2.65 gan1 -1.8 end -5.2 gan1 -1.8 nseg 15 
sel cable begin -2.65 gan1 -1.1 end -5.2 gan1 -1.1 nseg 15 
sel cable begin -2.65 gan1 -0.4 end -5.2 gan1 -0.4 nseg 15 
;顶 
sel cable begin -2.63 gan1 0.39 end -5.15 gan1 0.75 nseg 15 
sel cable begin -2.44 gan1 1.05 end -4.78 gan1 2.04 nseg 15 
sel cable begin -2.09 gan1 1.64 end -4.1 gan1 3.2 nseg 15 
sel cable begin -1.6 gan1 2.12 end -3.14 gan1 4.15 nseg 15 
sel cable begin -1.01 gan1 2.46 end -1.98 gan1 4.82 nseg 15 
sel cable begin -0.35 gan1 2.63 end -0.67 gan1 5.16 nseg 15 
sel cable begin 2.63 gan1 0.39 end 5.15 gan1 0.75 nseg 15 
sel cable begin 2.44 gan1 1.05 end 4.78 gan1 2.04 nseg 15 
sel cable begin 2.09 gan1 1.64 end 4.1 gan1 3.2 nseg 15 
sel cable begin 1.6 gan1 2.12 end 3.14 gan1 4.15 nseg 15 
sel cable begin 1.01 gan1 2.46 end 1.98 gan1 4.82 nseg 15 
sel cable begin 0.35 gan1 2.63 end 0.67 gan1 5.16 nseg 15 
end_command 
gan=gan+0.7 
end_loop 
end 
maogan 
 
;锚索锚杆参数 
sel cable pro emod 20e10 xcarea 3.8e-4 gr_per 0.088 yten 1e20 gr_k 1e10 gr_coh 1e20 
 
;给定边界条件 
fix z range z -15.1 -15 
fix x range x -15.1 -15 
fix x range x 14.9 15.1 
fix y range y -0.1 0.1 
fix y range y 29.9 30.1 
 
;监控变量 
his nstep=10 
hist unbal ;最大不平衡力 his 1 
his gp zdis 0 15 2.7 ; 顶板节点z轴方向位移 his 2 
his zone szz 0 15 2.7 ; 顶板单元体zz方向应力 his 3 
his gp zdis 0 15 -2.7 ; 底板节点z轴方向位移 his 4 
his zone szz 0 15 -2.7 ; 底板单元体zz方向应力 his 5 
his gp xdis -2.7 15 -1.3 ; 左帮节点x轴方向位移 his 6 
his zone sxx -2.7 15 -1.3 ; 左帮单元体xx方向应力 his 7 
his gp xdis 2.7 15 -1.3 ; 右帮节点x轴方向位移 his 8 
his zone sxx 2.7 15 -1.3 ; 右帮单元体xx方向应力 his 9 
 
;施加重力场 
set gravity 0 0 -9.81;设置初始重力加速度 
 
;设定密度 
ini density 2580 range z -13.16 15 ;岩体的密度 
ini density 1480 range z -15 -13.16 ;3#煤的密度 
 
;施加等效载荷 
app nstress -10e6 range x -15 15 y 0 30 z 15 
 
;设定边界位移速度 
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0  
ini state 0 
 
plot bl st  
set large 
solve 
 
save yizhihu.sav