Flac3d6.0源代码(含巷道开挖、支护与充填开采三组代码)直接导入运行版,附汉语注释说明
Flac3d6.0源代码用于巷道支护煤层开采和充填代码可直接导入运行并有汉语注释说明每句代码含义适合初学者用来学习。 直接联系。 包含三组代码 1巷道开挖研究巷道开挖后地压对巷道稳定性的影响 2巷道锚杆支护内置锚杆衬砌命令研究巷道在支护条件下地压对巷道稳定性的影响 3工作面充填开采研究煤层开采后采空区充填与否对上覆岩层的影响特征包括应力场、位移场和破坏场等等有充填体可以研究充填体条件下的围岩变化也可以删掉充填体研究采空区根据自己研究内容调整。巷道数值模拟这事儿好多新手刚开始总被复杂的建模步骤劝退。今儿咱们直接上手三个Flac3d6.0实战案例保证你拿着代码就能跑出巷道开挖、锚杆支护和充填开采的完整分析流程。下面这三个脚本我都加了逐句注释连坐标点怎么设置都给你标明白了。第一组巷道裸奔开挖model new ; 新建模型 zone create brick size 10 20 10 ; 搞个长条状煤层模型尺寸按实际改 zone cmodel assign mohr-coulomb ; 材料本构选莫尔库伦最常用 zone property bulk 5e9 shear 3e9 friction 35 ; 给岩石赋参数 model solve elastic ; 先做弹性初始化 zone face apply stress -10e6 range group top ; 顶部施加10MPa压力 zone face fix velocity range group bottom ; 底板固定 zone excavate range position-x 4 6 position-y 0 20 position-z 4 6 ; 关键在模型中间挖出2x20x2m巷道 model solve ; 开始计算应力重分布这里重点看zone excavate命令通过position参数精确定位开挖区域。运行后直接调出contour云图能清楚看到顶板应力集中和两帮位移变化——这些数据足够写篇小论文的第三章了。Flac3d6.0源代码用于巷道支护煤层开采和充填代码可直接导入运行并有汉语注释说明每句代码含义适合初学者用来学习。 直接联系。 包含三组代码 1巷道开挖研究巷道开挖后地压对巷道稳定性的影响 2巷道锚杆支护内置锚杆衬砌命令研究巷道在支护条件下地压对巷道稳定性的影响 3工作面充填开采研究煤层开采后采空区充填与否对上覆岩层的影响特征包括应力场、位移场和破坏场等等有充填体可以研究充填体条件下的围岩变化也可以删掉充填体研究采空区根据自己研究内容调整。第二组锚杆衬砌双保险struct node create by-zone-face range group excavation ; 在巷道表面生成节点 struct cable create by-line ... ; 锚杆参数设置间距、长度自己调 struct cable property young 2e10 yield-tension 3e5 ; 钢绞线属性 struct liner create by-face ... ; 喷浆层建模 struct liner property thickness 0.2 ; 20cm厚混凝土 model solve ; 带支护的计算锚杆支护的难点在结构单元和岩体的耦合。注意struct cable要配合zone face apply使用建议先用model cycle 1000做分步计算观察锚杆轴力变化趋势。这里有个小技巧把衬砌的young模量调低两个量级可以模拟柔性支护效果。第三组充填开采生死局zone group goaf range position-x 0 3 ; 标记采空区范围 zone cmodel assign null range group goaf ; 采空区材料失效 zone group backfill range position-x 0 3 position-z 0 2 ; 充填体占位 zone cmodel assign elastic range group backfill ; 充填材料 zone property young 1e8 poisson 0.3 range group backfill ; 充填体软如橡皮泥 model solve ; 有/无充填对比这段代码的精髓在于通过null模型和elastic模型的切换实现充填与否的对比。把young模量从1e8充填改成1e2松散垮落立马能看到顶板下沉量翻倍。建议跑完用plot create bitmap导出应力云图用PS叠图对比更直观。实际用的时候记得调整网格尺寸——模型太大容易算到地老天荒。建议新手先用5x5x5的小模型试跑参数调通后再放大尺寸。遇到不收敛的情况先检查材料单位是否统一别把MPa和Pa混用再试试把阻尼系数从默认的0.8降到0.6。有次我模型算崩了最后发现是锚杆端点没扎进岩体里这种低级错误大家引以为戒啊
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2478576.html
如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!