申请/专利权人：中国矿业大学(北京)

申请日：2024-03-21

公开（公告）日：2024-06-25

公开（公告）号：CN118243498A

主分类号：G01N3/08

分类号：G01N3/08;G01N3/06;G01N23/046;G01N29/14;G06F30/25;G06F119/14

优先权：

专利状态码：在审-公开

法律状态：2024.06.25#公开

摘要：本发明涉及采空区垮落带模拟领域，具体涉及考虑颗粒破碎的采空区垮落带破碎煤岩体压实模拟方法，包括获取实验室尺度的采空区垮落带破碎煤岩体；CT扫描并导入到PFC中，获取破碎煤岩体三维轮廓；建立外接立方墙体，投放ball；使用pebble代替ball；创建采空区墙体；投放颗粒簇；使用ball代替pebble；设置接触模型；对ball子颗粒进行赋参；建立加载墙体；设置加载速度；设置最大应力加载条件，求解并记录计算参数。本发明可以得到破碎煤岩体压实过程中的声发射变化情况及破碎煤岩体破碎情况；相比传统刚性颗粒簇计算方式，本发明能够实现建立考虑破碎煤岩体真实颗粒形状和破碎煤岩体压缩过程中破碎的胶结离散元模型，能够更加真实地反映破碎煤岩体的压实破碎过程。

主权项：1.考虑颗粒破碎的采空区垮落带破碎煤岩体压实模拟方法，其特征在于，包括以下步骤：a.自工程现场获取采空区垮落带大块煤岩体，在实验室对大块煤岩体加载破碎，获得实验室尺度的采空区垮落带破碎煤岩体；b.对采空区垮落带破碎煤岩体进行CT扫描，将CT扫描获取的STL文件导入到PFC中，使用loopforeach函数遍历所有破碎煤岩体的每一个点的坐标，获取破碎煤岩体三维轮廓；c.使用wallgenerate命令建立破碎煤岩体三维轮廓的外接立方墙体，向外接立方墙体内投放ball子颗粒，设置ball子颗粒的密度及刚度；待外接立方墙体结构运算平衡后，删除外接立方墙体外部的ball子颗粒，保留外接立方墙体内部的ball子颗粒；使用fish语言记录剩余的外接立方墙体内部的ball子颗粒的坐标，使用pebble关键字在每一个ball的位置生成一个pebble子颗粒，随后删除所有ball子颗粒及外接立方墙体，此时得到的为由pebble子颗粒组成的刚性颗粒簇，每个刚性颗粒簇代表一个破碎煤岩体；d.采用domainextent命令设定整个压实模型计算区域；创建三维采空区墙体，所述三维采空区墙体为顶部开口底面水平设置的长方体结构，所述三维采空区墙体内部形成空腔；e.将步骤c获取的所有刚性颗粒簇投放到采空区墙体内；使用loopforeach函数遍历所有刚性颗粒簇的pebble子颗粒坐标，然后利用ballcreat命令在每一个pebble子颗粒位置生成一个ball子颗粒，使用clumpdelete命令删除所有的pebble子颗粒，pebble子颗粒组成的刚性颗粒簇变成ball子颗粒组成柔性颗粒簇；设置ball子颗粒的密度及刚度，使用cycle命令使柔性颗粒簇运算平衡；f.使用cmat命令设置ball-ball之间以及ball-facet之间的默认接触为滑动摩擦阻力模型，只对同一柔性颗粒簇内的ball子颗粒施加接触，不同柔性颗粒簇之间的ball子颗粒不施加接触；使用contactmethod命令将ball-ball之间的接触转为平行粘结模型；对ball子颗粒进行赋参，包括抗拉强度，粘聚力，内摩擦角，摩擦系数，弹性模量，刚度比；g.设置整个压实模型的重力，使柔性颗粒簇利用自重下降并计算平衡；使用fish语言监测柔性颗粒簇平衡后的坐标，在所有柔性颗粒簇最外处使用wall命令建立一面加载墙体；h.采用第一公式计算柔性颗粒簇间的孔隙率pb，所述第一公式为 公式中，pm是测量圆测量的孔隙率，由measureporosity命令获取；Vm是测量圆的体积；S为采空区墙体的底面积；h0为柔性颗粒簇初始堆积高度；disp为加载墙体下降的位移，下降方向为Z方向向下，由walldisp命令获取；i.基于fish语言分别采用第二公式、第三公式、第四公式计算整个压实模型加载过程中Z方向的应力tzz、应变wezz、压实度Dc，同时使用history命令进行实时记录，所述第二公式为tzz＝FzzS；所述第三公式为所述第四公式为公式中，Fzz为加载墙体上的力，由wallforcecontact命令获取；j.在整个压实模型加载过程中使用fragment命令实时监测柔性颗粒簇的破碎情况，即破碎后柔性颗粒簇内子颗粒的数量；k.定义计算频率cal_frequency，记录频率record_frequency，初始裂纹数量liewen_record为0，采用fish语言进行判断，当record_frequency-wezz≤cal_frequency时，记录此时的声发射事件数量AE_num＝crack_num-liewen_record，同时令record_frequency＝wezz，liewen_record＝crack_num，循环记录每次的声发射事件数量，直至加载结束；l.使用wallattribute命令设置加载墙体的加载速度；m.设置整个压实模型的最大应力加载条件，使用solvefishhalt命令对整个压实模型求解。

全文数据：

权利要求：

百度查询： 中国矿业大学(北京) 考虑颗粒破碎的采空区垮落带破碎煤岩体压实模拟方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD