FLAC3D模拟实例 循环开挖与支护
合集下载
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
step 100
end_command
n=excav_p+cut_Liang
excav_p=excav_p+cut_Liang
;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
ytens_y=25e4 ;抗拉强度
gr_per_gp=0.11 ;水泥浆外圈长度
gr_k_gk=6e6 ;单位长度上水泥浆刚度
gr_fric_gf=30.0 ;水泥浆摩擦角
gr_coh_gc=1.0e5 ;单位长度上水泥浆的粘结力
;maogan5
xp5=82.1314
yp5=15.4709
xp5_2=84.4996
yp5_2=16.2719
;maogan6
xp6=82.25
yp6=14.0
xp6_2=84.75
yp6_2=14.0
ini xdis 0
ini zdis 0
ini yvel 0
ini xvel 0
ini zvel 0
set large
set mech ratio=5e-4
his gp ydis 90 17 40 ;40m处巷道巷道顶部中点
his gp xdis 77.75 14.75 40 ;40m处巷道巷道左帮中点
; gr_coh__gc= ;单位长度上水泥浆的粘结力
; ycomp__y= ;抗压强度
end
maogan_canshu
maosuo_canshu
def maogan(suo)_p
;锚杆位置
;maogan1
; ———————————————————————-————————————
def maogan_canshu
cable_seg=10 ;构件数
emod_e=180e9 ;弹性模量
xcarea_x=3.46185e-4 ;锚杆（索）横截面积
;pl add sel geom cable blue
;sel set damp combined
;sel liner prop slide on
;添加接触面
;gen separate 2
;inter 1 wrap 1 2
;inter 2 wrap 3 2
;inter 1 face range plane norm 0 1 0 ori 0 11.0815 0 dist 0.1
;每次开挖量
cut_0=excav_p
cut_1=excav_p+cut_liang
;开挖
command
m null range grou gzm z cut_0 cut_1
n
res ini.sav
set geometry=0.001
ini ydis 0
ini xdis 0
ini zdis 0
ini yvel 0
ini xvel 0
ini zvel 0
m m
prop bulk 4.0e9 shear 2.5e9 fri 32 coh 2.0e6 &
range grou diban-shayan
;
prop bulk 1.8e9 shear 1.2e9 fri 25 coh 1.0e6 &
range grou diban-niyan any grou hangdao any
;
prop bulk 1.2e9 shear 0.8e9 fri 22 coh 0.8e6 &
;inter 2 face range plane norm 0 1 0 ori 0 10.7193 0 dist 0.1
;interface 1 prop kn 1e10 ks 1e10
;interface 2 prop kn 1e10 ks 1e10
set large
ini ydis 0
end_command
end_if
end_loop
end
set cut_liang=8 excav_p=0 excav_p_z=80
excav_hangdao
solve
save kwhdno2.sav
call kwhdzh.txt
xp1=77.75
yp1=14.0
xp1_2=72.25
yp1_2=14.0
;maogan2
xp2=77.8686
yp2=15.4709
xp2_2=75.5004
yp2_2=16.2719
;maogan3
interface 1 prop kn 20e9 ks 10e9 tens 1e9
interface 2 prop kn 20e9 ks 10e9 tens 1e9
set mech ratio=5e-4
def excav_mc
loop n (excav_p,excav_p_z+cut_liang)
;
prop bulk 3.5e9 shear 2.5e9 fri 34 coh 1.4e6 &
range grou dingban-shayan
;添加接触面
gen separate gzm
interface 1 wrap mc gzm
interface 2 wrap dingban-fenshayan gzm
his gp xdis 82.25 14.75 40 ;40m处巷道巷道右帮中点
def excav_hangdao
loop n (excav_p,excav_p_z+cut_liang)
;每次开挖量
cut_0=excav_p
cut_1=excav_p+cut_liang
his gp xdis 82.25 14.75 40 ;40m处巷道巷道右帮Байду номын сангаас点
;监控锚杆索命令
;his sel cable force x y z
;his sel cable stress x y z
;his sel cable grout stress x y z ;监控锚杆索构件端部水泥浆应力(端部)
;开挖
command
m null range grou hangdao z cut_0 cut_1
step 50
end_command
n=excav_p+cut_Liang
;
;锚索位置
;maosuo1
xps1=77.75
yps1=13.3536
xps1_2=71.25
yps1_2=13.3536
; ycomp_y= ;抗压强度
end
def maosuo_canshu
s_cable_seg=12 ;构件数
emod__e=200e9 ;弹性模量
xcarea__x=3.46185e-4 ;锚杆（索）横截面积
;—————————————————————————————————————————————
; 锚杆（索）、初撑参数及剖面位置参数
;
;
;—————————————————————————————————————————————
if cut_1=40 then
command
save kwmc40m.sav
end_command
end_if
end_loop
end
;每次需要修改
set cut_liang=5 excav_p=0 excav_p_z=80
ytens__y=35e4 ;抗拉强度
gr_per__gp=0.11 ;水泥浆外圈长度
gr_k__gk=6e6 ;单位长度上水泥浆刚度
gr_fric__gf=30 ;水泥浆摩擦角
ini yvel 0
ini xvel 0
ini zvel 0
set large
set mech ratio=5e-4
pl cont ydis shad on outline on
his gp ydis 90 17 40 ;40m处巷道巷道顶部中点
his gp xdis 77.75 14.75 40 ;40m处巷道巷道左帮中点
;
prop bulk 3.0e9 shear 2.5e9 fri 30 coh 1.8e6 &
range grou dingban-fenshayan
;
prop bulk 1.5e9 shear 1.2e9 fri 25 coh 1.1e6 &
range grou dingban-niyan
xp3=79.0593
yp3=16.7939
xp3_2=78.0140
yp3_2=19.0649
;maogan4
xp4=80.9407
yp4=16.7939
xp4_2=81.9680
yp4_2=19.0649
;n
;res kwmc2.sav
;pl set rot 90 180 0 ;从后视图看，零点面对着
;pl set rot 90 0 0 ;从前视图看，终点面对着
pl set rot 90 75 15
pl cont syy range grou hangdao
pl add sel geom node off fill on black blue ;显示支护单元
excav_p=excav_p+cut_Liang
;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
if cut_1=40 then
command
save kwhdno40m.sav
range grou diban-gentuyan
;
prop bulk 1.9e9 shear 1.3e9 fri 24 coh 1.0e6 &
range grou diban-tniyan
;
prop bulk 0.7e9 shear 0.8e9 fri 21 coh 0.7e6 &
range grou mc any grou gzm any
;his sel cable grout slip x y z ;滑动状态（端部）（0未屈服1正在屈服2曾屈服）
;his sel cable grout disp x y z ;监控锚杆索构件端部在水泥浆中的位移
; ———————————————————————————————————
excav_mc
solve
pl block grou
pl add axes red
pl set rotation 90 0 0
save kwmc2.sav
call kwhdno.txt
;n
;res kwmc.sav
ini ydis 0
ini xdis 0
ini zdis 0
end_command
n=excav_p+cut_Liang
excav_p=excav_p+cut_Liang
;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
ytens_y=25e4 ;抗拉强度
gr_per_gp=0.11 ;水泥浆外圈长度
gr_k_gk=6e6 ;单位长度上水泥浆刚度
gr_fric_gf=30.0 ;水泥浆摩擦角
gr_coh_gc=1.0e5 ;单位长度上水泥浆的粘结力
;maogan5
xp5=82.1314
yp5=15.4709
xp5_2=84.4996
yp5_2=16.2719
;maogan6
xp6=82.25
yp6=14.0
xp6_2=84.75
yp6_2=14.0
ini xdis 0
ini zdis 0
ini yvel 0
ini xvel 0
ini zvel 0
set large
set mech ratio=5e-4
his gp ydis 90 17 40 ;40m处巷道巷道顶部中点
his gp xdis 77.75 14.75 40 ;40m处巷道巷道左帮中点
; gr_coh__gc= ;单位长度上水泥浆的粘结力
; ycomp__y= ;抗压强度
end
maogan_canshu
maosuo_canshu
def maogan(suo)_p
;锚杆位置
;maogan1
; ———————————————————————-————————————
def maogan_canshu
cable_seg=10 ;构件数
emod_e=180e9 ;弹性模量
xcarea_x=3.46185e-4 ;锚杆（索）横截面积
;pl add sel geom cable blue
;sel set damp combined
;sel liner prop slide on
;添加接触面
;gen separate 2
;inter 1 wrap 1 2
;inter 2 wrap 3 2
;inter 1 face range plane norm 0 1 0 ori 0 11.0815 0 dist 0.1
;每次开挖量
cut_0=excav_p
cut_1=excav_p+cut_liang
;开挖
command
m null range grou gzm z cut_0 cut_1
n
res ini.sav
set geometry=0.001
ini ydis 0
ini xdis 0
ini zdis 0
ini yvel 0
ini xvel 0
ini zvel 0
m m
prop bulk 4.0e9 shear 2.5e9 fri 32 coh 2.0e6 &
range grou diban-shayan
;
prop bulk 1.8e9 shear 1.2e9 fri 25 coh 1.0e6 &
range grou diban-niyan any grou hangdao any
;
prop bulk 1.2e9 shear 0.8e9 fri 22 coh 0.8e6 &
;inter 2 face range plane norm 0 1 0 ori 0 10.7193 0 dist 0.1
;interface 1 prop kn 1e10 ks 1e10
;interface 2 prop kn 1e10 ks 1e10
set large
ini ydis 0
end_command
end_if
end_loop
end
set cut_liang=8 excav_p=0 excav_p_z=80
excav_hangdao
solve
save kwhdno2.sav
call kwhdzh.txt
xp1=77.75
yp1=14.0
xp1_2=72.25
yp1_2=14.0
;maogan2
xp2=77.8686
yp2=15.4709
xp2_2=75.5004
yp2_2=16.2719
;maogan3
interface 1 prop kn 20e9 ks 10e9 tens 1e9
interface 2 prop kn 20e9 ks 10e9 tens 1e9
set mech ratio=5e-4
def excav_mc
loop n (excav_p,excav_p_z+cut_liang)
;
prop bulk 3.5e9 shear 2.5e9 fri 34 coh 1.4e6 &
range grou dingban-shayan
;添加接触面
gen separate gzm
interface 1 wrap mc gzm
interface 2 wrap dingban-fenshayan gzm
his gp xdis 82.25 14.75 40 ;40m处巷道巷道右帮中点
def excav_hangdao
loop n (excav_p,excav_p_z+cut_liang)
;每次开挖量
cut_0=excav_p
cut_1=excav_p+cut_liang
his gp xdis 82.25 14.75 40 ;40m处巷道巷道右帮Байду номын сангаас点
;监控锚杆索命令
;his sel cable force x y z
;his sel cable stress x y z
;his sel cable grout stress x y z ;监控锚杆索构件端部水泥浆应力(端部)
;开挖
command
m null range grou hangdao z cut_0 cut_1
step 50
end_command
n=excav_p+cut_Liang
;
;锚索位置
;maosuo1
xps1=77.75
yps1=13.3536
xps1_2=71.25
yps1_2=13.3536
; ycomp_y= ;抗压强度
end
def maosuo_canshu
s_cable_seg=12 ;构件数
emod__e=200e9 ;弹性模量
xcarea__x=3.46185e-4 ;锚杆（索）横截面积
;—————————————————————————————————————————————
; 锚杆（索）、初撑参数及剖面位置参数
;
;
;—————————————————————————————————————————————
if cut_1=40 then
command
save kwmc40m.sav
end_command
end_if
end_loop
end
;每次需要修改
set cut_liang=5 excav_p=0 excav_p_z=80
ytens__y=35e4 ;抗拉强度
gr_per__gp=0.11 ;水泥浆外圈长度
gr_k__gk=6e6 ;单位长度上水泥浆刚度
gr_fric__gf=30 ;水泥浆摩擦角
ini yvel 0
ini xvel 0
ini zvel 0
set large
set mech ratio=5e-4
pl cont ydis shad on outline on
his gp ydis 90 17 40 ;40m处巷道巷道顶部中点
his gp xdis 77.75 14.75 40 ;40m处巷道巷道左帮中点
;
prop bulk 3.0e9 shear 2.5e9 fri 30 coh 1.8e6 &
range grou dingban-fenshayan
;
prop bulk 1.5e9 shear 1.2e9 fri 25 coh 1.1e6 &
range grou dingban-niyan
xp3=79.0593
yp3=16.7939
xp3_2=78.0140
yp3_2=19.0649
;maogan4
xp4=80.9407
yp4=16.7939
xp4_2=81.9680
yp4_2=19.0649
;n
;res kwmc2.sav
;pl set rot 90 180 0 ;从后视图看，零点面对着
;pl set rot 90 0 0 ;从前视图看，终点面对着
pl set rot 90 75 15
pl cont syy range grou hangdao
pl add sel geom node off fill on black blue ;显示支护单元
excav_p=excav_p+cut_Liang
;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
if cut_1=40 then
command
save kwhdno40m.sav
range grou diban-gentuyan
;
prop bulk 1.9e9 shear 1.3e9 fri 24 coh 1.0e6 &
range grou diban-tniyan
;
prop bulk 0.7e9 shear 0.8e9 fri 21 coh 0.7e6 &
range grou mc any grou gzm any
;his sel cable grout slip x y z ;滑动状态（端部）（0未屈服1正在屈服2曾屈服）
;his sel cable grout disp x y z ;监控锚杆索构件端部在水泥浆中的位移
; ———————————————————————————————————
excav_mc
solve
pl block grou
pl add axes red
pl set rotation 90 0 0
save kwmc2.sav
call kwhdno.txt
;n
;res kwmc.sav
ini ydis 0
ini xdis 0
ini zdis 0