FLAC3D模拟实例循环开挖与支护
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
FLAC3D模拟实例循环开挖与支护n
res ini.sav
set geometry=0.001
ini ydis0
ini xdis0
ini zdis0
ini yvel0
ini xvel0
ini zvel0
m m
prop bulk 4.0e9shear 2.5e9fri32coh 2.0e6& range grou diban-shayan
;
prop bulk 1.8e9shear 1.2e9fri25coh 1.0e6& range grou diban-niyan any grou hangdao any
;
prop bulk 1.2e9shear0.8e9fri22coh0.8e6& range grou diban-gentuyan
;
prop bulk 1.9e9shear 1.3e9fri24coh 1.0e6& range grou diban-tniyan
;
prop bulk0.7e9shear0.8e9fri21coh0.7e6& range grou mc any grou gzm any
;
prop bulk 3.0e9shear 2.5e9fri30coh 1.8e6& range grou dingban-fenshayan
;
prop bulk 1.5e9shear 1.2e9fri25coh 1.1e6& range grou dingban-niyan
;
prop bulk 3.5e9shear 2.5e9fri34coh 1.4e6& range grou dingban-shayan
;添加接触面
gen separate gzm
interface1wrap mc gzm
interface2wrap dingban-fenshayan gzm interface1prop kn20e9ks10e9tens1e9 interface2prop kn20e9ks10e9tens1e9
set mech ratio=5e-4
def excav_mc
loop n(excav_p,excav_p_z+cut_liang)
;每次开挖量
cut_0=excav_p
cut_1=excav_p+cut_liang
;开挖
command
m null range grou gzm z cut_0cut_1
step100
end_command
n=excav_p+cut_Liang
excav_p=excav_p+cut_Liang
;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
if cut_1=40then
command
save kwmc40m.sav
end_command
end_if
end_loop
end
;每次需要修改
set cut_liang=5excav_p=0excav_p_z=80
excav_mc
solve
pl block grou
pl add axes red
pl set rotation9000
save kwmc2.sav
call kwhdno.txt
;n
;res kwmc.sav
ini ydis0
ini xdis0
ini zdis0
ini yvel0
ini xvel0
ini zvel0
set large
set mech ratio=5e-4
pl cont ydis shad on outline on
his gp ydis901740;40m处巷道巷道顶部中点
his gp xdis77.7514.7540;40m处巷道巷道左帮中点
his gp xdis82.2514.7540;40m处巷道巷道右帮中点
def excav_hangdao
loop n(excav_p,excav_p_z+cut_liang)
;每次开挖量
cut_0=excav_p
cut_1=excav_p+cut_liang
;开挖
command
m null range grou hangdao z cut_0cut_1
step50
end_command
n=excav_p+cut_Liang
excav_p=excav_p+cut_Liang
;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。
if cut_1=40then
command
save kwhdno40m.sav
end_command
end_if
end_loop
end
set cut_liang=8excav_p=0excav_p_z=80
excav_hangdao
solve
save kwhdno2.sav
call kwhdzh.txt
;n
;res kwmc2.sav
;pl set rot901800;从后视图看，零点面对着
;pl set rot9000;从前视图看，终点面对着
pl set rot907515
pl cont syy range grou hangdao
pl add sel geom node off fill on black blue;显示支护单元
;pl add sel geom cable blue
;sel set damp combined