CN117473797A - 充填节理岩石损伤度离散元数值模拟方法及系统 - Google Patents
充填节理岩石损伤度离散元数值模拟方法及系统 Download PDFInfo
- Publication number
- CN117473797A CN117473797A CN202311804203.1A CN202311804203A CN117473797A CN 117473797 A CN117473797 A CN 117473797A CN 202311804203 A CN202311804203 A CN 202311804203A CN 117473797 A CN117473797 A CN 117473797A
- Authority
- CN
- China
- Prior art keywords
- rock
- filling joint
- layer
- numerical model
- microcracks
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 260
- 230000006378 damage Effects 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004088 simulation Methods 0.000 title claims abstract description 29
- 230000009471 action Effects 0.000 claims abstract description 53
- 230000035515 penetration Effects 0.000 claims abstract description 49
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 20
- 239000011707 mineral Substances 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 11
- 238000007906 compression Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000006835 compression Effects 0.000 claims description 5
- 230000036962 time dependent Effects 0.000 claims description 5
- 238000011161 development Methods 0.000 abstract description 6
- 239000002245 particle Substances 0.000 description 13
- 239000000463 material Substances 0.000 description 8
- 238000004891 communication Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 208000027418 Wounds and injury Diseases 0.000 description 2
- 230000003750 conditioning effect Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 208000014674 injury Diseases 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011236 particulate material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Abstract
本发明公开了一种充填节理岩石损伤度离散元数值模拟方法及系统，方法包括S100：建立充填节理岩石数值模型，通过与实际物理试验结果对比进行细观参数标定，赋予充填节理岩石数值模型通过标定后的细观参数，确定充填节理层和岩石层的损伤权重因子；S200：建立落锤数值模型，并赋予落锤数值模型的位置坐标、密度、下落初始速度和冲击接触时间参数；S300：根据微裂纹数目以及微裂纹出现的位置数据计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度。本发明以微裂纹发育情况评估冲击作用下充填节理岩石的损伤度，同时引入损伤权重因子，可以更准确地评估充填节理岩石整体的损伤程度，更贴合实际情况，精度更高。
Description
技术领域
本发明属于充填节理岩石损伤测定技术领域，更具体地，涉及一种充填节理岩石损伤度离散元数值模拟方法及系统。
背景技术
工程岩体中随机分布有大量的节理，而节理往往含有一定的充填物。充填节理岩石不同于完整岩石和非充填节理岩石，充填节理具有低强度和大变形特性，其力学强度和变形特性更差，充填节理对整个岩石的强度和稳定性有着很大影响。实际中充填节理岩石不可避免地会遭受到爆破等冲击作用，会产生一定的损伤，而损伤程度关乎到充填节理岩石的安全稳定性，所以衡量冲击作用下充填节理岩石的损伤程度至关重要。
现有技术中往往通过波速值、抗压强度以及体积变化率等特征的变化反映损伤度，公开号为CN111208198A的中国发明专利公开了一种岩体实时波速测定及质量评价的方法，通过在不同超声波发送位置、不同时间不断按动超声波发送装置的按钮，并且不断计算每次超声波传播到每个传感器的速度，根据波速、衰减系数、波形、频率、频谱及振幅等参数变化值可反演得到岩土体的物理力学指标及细观结构特征，如节理的变化情况，进而分析矿山岩体稳定情况。公开号为CN113109182B的中国发明专利公开了一种基于体积变化率确定岩石节理动态损伤的实验方法，分别基于单次冲击和循环冲击的一维分离式霍普金森压杆SHPB实验和三维节理形貌扫描技术，可以确定岩石节理在单次冲击或循环冲击下的动态损伤状态，其与传统的基于静态或准静态力学实验的超声检测相比，更适用于岩石节理在爆破或地震作用等动荷载作用下损伤程度的测量，直接采用岩石节理本身的形貌变化反映节理的损伤程度。
但考虑到充填节理岩石可能具有复杂的内部结构，包括填充物和岩石的层状分布，可能会干扰波速值和抗压强度的变化，难以准确反映实际损伤程度；此外，充填节理岩石可能在不同部分存在非均质性，不同区域的波速值和抗压强度可能不同，难以用统一的标准评估损伤。现有研究较少从冲击作用下充填节理岩石的裂纹发育情况角度评判损伤程度，而裂纹发育情况往往决定岩体是否失稳破坏，裂纹发育情况对评价岩体稳定性具有重要参考意义。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供一种充填节理岩石损伤度离散元数值模拟方法及系统，通过构建充填节理岩石数值模型与落锤数值模型，模拟落锤数值模型冲击充填节理岩石数值模型，记录并输出冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据，根据充填节理层与岩石层细观参数的比例关系和充填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子，计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，以模型中模拟形成微裂纹的发育情况这一特征评估冲击作用下充填节理岩石的损伤度，更贴合实际情况，具有精度高的优点。
为了实现上述目的，按照本发明的第一方面，提供一种充填节理岩石损伤度离散元数值模拟方法，包括：
S100：建立充填节理岩石数值模型，通过与实际物理试验结果对比进行细观参数标定，赋予充填节理岩石数值模型通过标定后的细观参数，根据充填节理层与岩石层细观参数的比例关系和充填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子a1和a2：
；
；
其中，emod1和emod2分别表示充填节理层和岩石层细观参数弹性模量的数值，h1和h2分别表示充填节理层和岩石层的厚度值；
S200：建立落锤数值模型，并赋予落锤数值模型的位置坐标、密度、下落初始速度和冲击接触时间参数，并模拟落锤数值模型冲击充填节理岩石数值模型，记录冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据；
S300：根据微裂纹数目以及微裂纹出现的位置数据计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，当最大贯通度大于阈值时，冲击作用下充填节理岩石损伤度为最大贯通度数值；当充填节理岩石数值模型的最大贯通度小于阈值时充填节理岩石的损伤度为：
；
式中，D表示冲击作用下充填节理岩石损伤度，f1和f2分别表示落锤数值模型冲击充填节理岩石数值模型下充填节理层和岩石层的微裂纹数目，c1和c2分别表示充填节理层和岩石层的接触数目。
进一步地，步骤S300中，所述阈值的确定包括：
S301：对建立考虑矿物组分、含量及矿物排列位置的充填节理岩石数值模型进行单轴压缩模拟；
S302：输出单轴压缩过程中充填节理岩石数值模型微裂纹数目随时间的变化曲线；
S303：提取变化曲线中微裂纹数目增长速度变化量最大的数据点，该数据点时刻充填节理岩石数值模型所对应的贯通度即为阈值。
进一步地，所述阈值与充填节理岩石的矿物组分、含量以及矿物排列位置相关。
进一步地，步骤S300中，所述最大贯通度为任一平面微裂纹所形成的连线与充填节理岩石数值模型在连线方向的总长度的比值。
进一步地，步骤S100中，充填节理岩石数值模型包括充填节理层和岩石层的接触数目、细观参数及厚度之类的参数。
进一步地，步骤S100中，所述充填节理岩石数值模型接触为平行黏结模型。
按照本发明的第二方面，提供一种充填节理岩石损伤度离散元数值模拟系统，包括：
损伤权重确定模块，用于建立充填节理岩石数值模型，通过与实际物理试验结果对比进行细观参数标定，赋予充填节理岩石数值模型通过标定后的细观参数，根据充填节理层与岩石层细观参数的比例关系和充填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子a1和a2：
；
；
其中，emod1和emod2分别表示充填节理层和岩石层细观参数弹性模量的数值，h1和h2分别表示充填节理层和岩石层的厚度值；
落锤数值模拟模块，用于建立落锤数值模型，并赋予落锤数值模型的位置坐标、密度、下落初始速度和冲击接触时间参数，并模拟落锤数值模型冲击充填节理岩石数值模型，记录冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据；
损伤度计算模块，用于根据微裂纹数目以及微裂纹出现的位置数据计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，当最大贯通度大于阈值时，冲击作用下充填节理岩石损伤度为最大贯通度数值；当充填节理岩石数值模型的最大贯通度小于阈值时充填节理岩石的损伤度为：
；
式中，D表示冲击作用下充填节理岩石损伤度，f1和f2分别表示落锤数值模型冲击充填节理岩石数值模型下充填节理层和岩石层的微裂纹数目，c1和c2分别表示充填节理层和岩石层的接触数目。
按照本发明的第三方面，提供一种终端装置，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现充填节理岩石损伤度离散元数值模拟方法的步骤。
按照本发明的第四方面，提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现充填节理岩石损伤度离散元数值模拟方法的步骤。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
1.本发明的方法，通过构建充填节理岩石数值模型与落锤数值模型，模拟落锤数值模型冲击充填节理岩石数值模型，记录并输出冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据，根据充填节理层与岩石层细观参数的比例关系和充填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子，计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，以模型中模拟形成微裂纹的发育情况这一特征评估冲击作用下充填节理岩石的损伤度，更贴合实际情况，具有精度高的优点。
2.本发明通过引入损伤权重因子这一系数，基于充填节理层与岩石层细观参数弹性模量的比例关系以及厚度比，对充填节理的软弱特性进行定量分析，可以更准确地评估充填节理岩石整体的损伤程度。
3.本发明在确定冲击作用下充填节理岩石的损伤度时综合考虑微裂纹数目以及微裂纹分布特性的影响，采用对比最大贯通度与阈值的关系确定损伤度，相较于基于单一因素评估充填节理岩石损伤度的方法，所确定的损伤度具有更高的可靠性。
4.本发明通过使用颗粒流数值模拟软件，模拟落锤数值模型冲击充填节理岩石数值模型，记录并输出冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据，可以直观呈现冲击作用下充填节理岩石的微裂纹数目和分布情况，具有良好的可视化效果。
附图说明
图1为本发明实施例一种充填节理岩石损伤度离散元数值模拟方法的流程图；
图2为本发明实施例一种充填节理岩石损伤度离散元数值模拟方法的原理图；
图3为本发明一个实施例中的落锤数值模型冲击充填节理岩石数值模型示意图；
图4为本发明一个实施例中的受冲击后充填节理岩石数值模型裂纹图；
图5为本发明一个实施例中的充填节理岩石数值模型微裂纹数目随时间的变化曲线。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1-图5所示，本发明实施例提供一种充填节理岩石损伤度离散元数值模拟方法，该方法包括以下步骤：
S1：建立充填节理岩石数值模型，所述充填节理岩石数值模型的接触为平行黏结模型；记录充填节理层和岩石层的接触数目、细观参数及厚度，所述细观参数通过与实际物理试验结果比对进行标定；计算充填节理层和岩石层的损伤权重因子a1和a2：
；
；
其中，emod1和emod2分别表示充填节理层和岩石层细观参数弹性模量的数值，h1和h2分别表示充填节理层和岩石层的厚度值。本发明通过引入损伤权重因子这一系数，基于充填节理层与岩石层细观参数弹性模量的比例关系以及厚度比，对充填节理的软弱特性进行定量分析，可以更准确地评估充填节理岩石整体的损伤程度。
S2：建立落锤数值模型，建立所述落锤数值模型包括位置坐标、密度、下落初始速度和冲击接触时间参数的赋予，所述参数由CAD数据导入。
S3：模拟落锤数值模型冲击充填节理岩石数值模型，记录并输出冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据，计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度；所述最大贯通度为任一平面微裂纹所形成的连线与充填节理岩石数值模型范围内该平面总长度的比值。
S4：当充填节理岩石数值模型的最大贯通度小于阈值时，计算冲击作用下充填节理岩石的损伤度：
；
其中，D表示冲击作用下充填节理岩石损伤度，f1和f2分别表示落锤数值模型冲击充填节理岩石数值模型下充填节理层和岩石层的微裂纹数目，c1和c2分别表示充填节理层和岩石层的接触数目；
当最大贯通度大于阈值时，冲击作用下充填节理岩石损伤度为最大贯通度数值。
其中，所述阈值与充填节理岩石的矿物组分、含量以及矿物排列位置相关，确定所述阈值包括以下步骤：
步骤一，对建立考虑矿物组分、含量及矿物排列位置的充填节理岩石数值模型进行单轴压缩模拟；
步骤二，输出单轴压缩过程中充填节理岩石数值模型微裂纹数目随时间的变化曲线；
步骤三，提取变化曲线中微裂纹数目增长速度变化量最大的数据点，该数据点时刻充填节理岩石数值模型所对应的贯通度即为阈值。本发明在确定冲击作用下充填节理岩石的损伤度时综合考虑微裂纹数目以及微裂纹分布特性的影响，采用对比最大贯通度与阈值的关系确定损伤度，相较于基于单一因素评估充填节理岩石损伤度的方法，所确定的损伤度具有更高的可靠性。
本发明的方法，通过构建充填节理岩石数值模型与落锤数值模型，模拟落锤数值模型冲击充填节理岩石数值模型，记录并输出冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据，根据充填节理层与岩石层细观参数的比例关系和充填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子，计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，以模型中模拟形成微裂纹的发育情况这一特征评估冲击作用下充填节理岩石的损伤度，更贴合实际情况，具有精度高的优点。
此外，本发明实施例还提供一种充填节理岩石损伤度离散元数值模拟系统，包括：
损伤权重确定模块，用于建立充填节理岩石数值模型，通过与实际物理试验结果对比进行细观参数标定，赋予充填节理岩石数值模型通过标定后的细观参数，根据充填节理层与岩石层细观参数的比例关系和充填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子a1和a2：
；
；
其中，emod1和emod2分别表示充填节理层和岩石层细观参数弹性模量的数值，h1和h2分别表示充填节理层和岩石层的厚度值；
落锤数值模拟模块，用于建立落锤数值模型，并赋予落锤数值模型的位置坐标、密度、下落初始速度和冲击接触时间参数，并模拟落锤数值模型冲击充填节理岩石数值模型，记录冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据；
损伤度计算模块，用于根据微裂纹数目以及微裂纹出现的位置数据计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，当最大贯通度大于阈值时，冲击作用下充填节理岩石损伤度为最大贯通度数值；当充填节理岩石数值模型的最大贯通度小于阈值时充填节理岩石的损伤度为：
；
式中，D表示冲击作用下充填节理岩石损伤度，f1和f2分别表示落锤数值模型冲击充填节理岩石数值模型下充填节理层和岩石层的微裂纹数目，c1和c2分别表示充填节理层和岩石层的接触数目。
实施例1
在颗粒流数值模拟软件PFC中建立充填节理岩石数值模型：充填节理岩石数值模型的尺寸为50mm*35mm；赋予充填节理岩石数值模型接触为平行黏结模型并记录充填节理层和岩石层的接触数目c1个和c2个，分别为1856个和5972个；通过与实际物理实验结果对比进行细观参数标定，赋予充填节理岩石数值模型通过标定后的细观参数，细观参数表如下：
表1 细观参数表
| 组 | pb_emod/GPa | pb_krat | pb_ten/MPa | pb_coh/MPa | pb_fa/° | emod/GPa | krat |
| 充填节理 | 0.125 | 1.65 | 2.76 | 2.76 | 25 | 0.125 | 1.65 |
| 岩石 | 2.0 | 1.8 | 46 | 46 | 27.4 | 2.0 | 1.8 |
表中：pb_emod为平行粘结弹性模量、pb_krat为平行粘结刚度比、pb_ten为平行粘结抗拉强度、pb_coh为平行粘结强度、pb_fa平行粘结摩擦角、emod为弹性模量、krat为刚度比。
岩土是一种地质材料，它经受长期的地质构造作用，在一定的地质环境中形成一定的结构，显现出宽广和多变的材料响应范围。与一般工程材料相比，表现出：（1）结构上的不连续性、不均匀性和各向异性；（2）物理力学性质上的非线性；（3）随机性、模糊性和不确定性；（4）地质条件和工程条件的复杂性；（5）整个岩土工程系统的非线性。
基于岩石结构与力学性质上的不连续性和非线性，采用颗粒流数值模拟（particle follow code，PFC）进行建模，模拟岩石这种非连续介质材料的力学行为。颗粒流数值模拟是一种离散单元法，通过圆形颗粒介质的运动及其相互作用来模拟颗粒材料的力学特性。在这种颗粒单元研究的基础上，通过非连续的数值方法来解决包含复杂变形模式的实际问题。在具有颗粒结构特性的充填节理岩石中的应用，就是从其细观力学特征出发，将材料的力学响应问题从物理域映射到数学域内进行数值求解。与此相应，物理域内实物颗粒被抽象为数学域内的颗粒单元，并通过颗粒单元来构建和设计任意几何性状的试样，颗粒间的相互作用通过接触本构关系来实现，数值边界条件的确定和试样的若干应力平衡状态通过迭代分析进行，直到使数值试样的宏观力学特性逼近材料的真实力学行为或者工程特性。
PFC程序中内置接触模型有多种：线性模型、接触黏结模型、平行黏结模型、平直节理模型和光滑节理模型等，其中平行黏结模型可以很好地模拟岩石材料的力学行为。本发明通过使用颗粒流数值模拟软件，模拟落锤数值模型冲击充填节理岩石数值模型，记录并输出冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据，可以直观呈现冲击作用下充填节理岩石的微裂纹数目和分布情况，具有良好的可视化效果。
本实施例中，根据充填节理层与岩石层细观参数弹性模量（emod）的比例关系和充
填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子a1和a2；充填节理层
和岩石层的损伤权重因子a1和a2的计算公式为：；
；
其中：emod1和emod2表示充填节理层和岩石层细观参数弹性模量（emod）的数值，分别为0.125GPa和2.0GPa；h1和h2充填节理层和岩石层的厚度值，分别为5mm和30mm；经计算，a1和a2的数值分别为0.672和0.328。
本实施例中，构建落锤数值模型，落锤数值模型特征由CAD数据导入，落锤数值模型如图3所示，赋予落锤数值模型的位置坐标为（0，0.2m）、密度为2.7×103kg/m³、下落初始速度为0m/s和冲击接触时间参数为0.01s；
本实施例中，模拟落锤数值模型冲击充填节理岩石数值模型，编制命令流分别记录并输出落锤数值模型冲击充填节理岩石数值模型下充填节理层和岩石层的微裂纹数目f1个和f2个，分别为1072个和679个，以及微裂纹出现位置数据，微裂纹位置如图4所示；
本实施例中，计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，最大贯通度为任一平面微裂纹所形成的连线与充填节理岩石数值模型范围内该平面总长度的比值，其值为0.387；
本实施例中，计算冲击作用下充填节理岩石数值模型的阈值；
所述阈值与充填节理岩石的矿物组分、含量及矿物排列位置相关，所述阈值确定的步骤为：
步骤一、对建立考虑矿物组分、含量及矿物排列位置的充填节理岩石数值模型进行单轴压缩模拟；
步骤二、输出单轴压缩过程中充填节理岩石数值模型微裂纹数目随时间的变化曲线，如图5所示；
步骤三、提取变化曲线中微裂纹数目增长速度变化量最大的数据点，所选数据点如图5所示，该数据点时刻充填节理岩石数值模型所对应的贯通度即为阈值，该阈值为0.7。
本实施例中，所述最大贯通度小于阈值，冲击作用下充填节理岩石损伤度计算公
式为：，经计算其数值为0.425。
故本实施中例冲击作用下充填节理岩石的损伤度为0.425。
实施例2
在颗粒流数值模拟软件PFC中建立充填节理岩石数值模型：充填节理岩石数值模型的尺寸为60mm*40mm；赋予充填节理岩石数值模型接触为平行黏结模型并记录充填节理层和岩石层的接触数目c1个和c2个，分别为2562个和6137个；通过与实际物理实验结果对比进行细观参数标定，赋予充填节理岩石数值模型通过标定后的细观参数。
本实施例中，根据充填节理层与岩石层细观参数弹性模量（emod）的比例关系和充
填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子a1和a2；充填节理层
和岩石层的损伤权重因子a1和a2的计算公式为：；
；
其中：emod1和emod2表示充填节理层和岩石层细观参数弹性模量（emod）的数值，分别为0.175GPa和2.5GPa；h1和h2充填节理层和岩石层的厚度值，分别为10mm和30mm；经计算，a1和a2的数值分别为0.467和0.533。
本实施例中，构建落锤数值模型，落锤数值模型特征由CAD数据导入，赋予落锤数值模型的位置坐标为（0，0.25m）、密度为3.2×103kg/m³、下落初始速度为0m/s和冲击接触时间参数为0.01s；
本实施例中，模拟落锤数值模型冲击充填节理岩石数值模型，编制命令流分别记录并输出落锤数值模型冲击充填节理岩石数值模型下充填节理层和岩石层的微裂纹数目f1个和f2个，分别为1929个和1195个，以及微裂纹出现位置数据；
本实施例中，计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，最大贯通度为任一平面微裂纹所形成的连线与充填节理岩石数值模型范围内该平面总长度的比值，其值为0.633；
本实施例中，计算冲击作用下充填节理岩石数值模型的阈值；
所述阈值与充填节理岩石的矿物组分、含量及矿物排列位置相关，所述阈值确定的步骤为：
步骤一、对建立考虑矿物组分、含量及矿物排列位置的充填节理岩石数值模型进行单轴压缩模拟；
步骤二、输出单轴压缩过程中充填节理岩石数值模型微裂纹数目随时间的变化曲线；
步骤三、提取变化曲线中微裂纹数目增长速度变化量最大的数据点，该数据点时刻充填节理岩石数值模型所对应的贯通度即为阈值，该阈值为0.6。
本实施例中，所述最大贯通度大于阈值，冲击作用下充填节理岩石损伤度为最大贯通度，其数值为0.633。
故本实施中例冲击作用下充填节理岩石的损伤度为0.633。
实施例3
本申请还提供了一种终端装置，该终端装置可包括：
存储器，用于存储计算机程序；
处理器，用于执行计算机程序时可实现如上述任意一种确定冲击作用下充填节理岩石损伤度方法的步骤。
终端装置可以包括：处理器、存储器、通信接口和通信总线。处理器、存储器、通信接口均通过通信总线完成相互间的通信。
在本申请实施例中，处理器可以为中央处理器（Central Processing Unit，CPU）、特定应用集成电路、数字信号处理器、现场可编程门阵列或者其他可编程逻辑器件等。
处理器可以调用存储器中存储的程序，具体的，处理器可以执行异常IP识别方法的实施例中的操作。
存储器中用于存放一个或者一个以上程序，程序可以包括程序代码，程序代码包括计算机操作指令，在本申请实施例中，存储器中至少存储有用于实现确定冲击作用下充填节理岩石损伤度方法功能的程序。
实施例4
本发明的一个实施例还公开了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如本发明第一方面所述的充填节理岩石损伤度离散元数值模拟方法的步骤。
本发明实现上述实施例方法中的全部或部分步骤，可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储的介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器（ROM，Read-OnlyMemory）、随机存取存储器（RAM，Random Access Memory）、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括是电载波信号和电信信号。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (9)
1.一种充填节理岩石损伤度离散元数值模拟方法，其特征在于，包括：
S100：建立充填节理岩石数值模型，通过与实际物理试验结果对比进行细观参数标定，赋予充填节理岩石数值模型通过标定后的细观参数，根据充填节理层与岩石层细观参数的比例关系和充填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子a1和a2：
；
；
其中，emod1和emod2分别表示充填节理层和岩石层细观参数弹性模量的数值，h1和h2分别表示充填节理层和岩石层的厚度值；
S200：建立落锤数值模型，并赋予落锤数值模型的位置坐标、密度、下落初始速度和冲击接触时间参数，并模拟落锤数值模型冲击充填节理岩石数值模型，记录冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据；
S300：根据微裂纹数目以及微裂纹出现的位置数据计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，当最大贯通度大于阈值时，冲击作用下充填节理岩石损伤度为最大贯通度数值；当充填节理岩石数值模型的最大贯通度小于阈值时充填节理岩石的损伤度为：
；
式中，D表示冲击作用下充填节理岩石损伤度，f1和f2分别表示落锤数值模型冲击充填节理岩石数值模型下充填节理层和岩石层的微裂纹数目，c1和c2分别表示充填节理层和岩石层的接触数目。
2.根据权利要求1所述的一种充填节理岩石损伤度离散元数值模拟方法，其特征在于，步骤S300中，所述阈值的确定包括：
S301：对建立考虑矿物组分、含量及矿物排列位置的充填节理岩石数值模型进行单轴压缩模拟；
S302：输出单轴压缩过程中充填节理岩石数值模型微裂纹数目随时间的变化曲线；
S303：提取变化曲线中微裂纹数目增长速度变化量最大的数据点，该数据点时刻充填节理岩石数值模型所对应的贯通度即为阈值。
3.根据权利要求2所述的一种充填节理岩石损伤度离散元数值模拟方法，其特征在于，所述阈值与充填节理岩石的矿物组分、含量以及矿物排列位置相关。
4.根据权利要求3所述的一种充填节理岩石损伤度离散元数值模拟方法，其特征在于，步骤S300中，所述最大贯通度为任一平面微裂纹所形成的连线与充填节理岩石数值模型在连线方向的总长度的比值。
5.根据权利要求1-4中任一项所述的一种充填节理岩石损伤度离散元数值模拟方法，其特征在于，步骤S100中，充填节理岩石数值模型包括充填节理层和岩石层的接触数目、细观参数及厚度之类的参数。
6.根据权利要求5所述的一种充填节理岩石损伤度离散元数值模拟方法，其特征在于，步骤S100中，所述充填节理岩石数值模型接触为平行黏结模型。
7.一种充填节理岩石损伤度离散元数值模拟系统，其特征在于，包括：
损伤权重确定模块，用于建立充填节理岩石数值模型，通过与实际物理试验结果对比进行细观参数标定，赋予充填节理岩石数值模型通过标定后的细观参数，根据充填节理层与岩石层细观参数的比例关系和充填节理层与岩石层的厚度比确定充填节理层和岩石层的损伤权重因子a1和a2：
；
；
其中，emod1和emod2分别表示充填节理层和岩石层细观参数弹性模量的数值，h1和h2分别表示充填节理层和岩石层的厚度值；
落锤数值模拟模块，用于建立落锤数值模型，并赋予落锤数值模型的位置坐标、密度、下落初始速度和冲击接触时间参数，并模拟落锤数值模型冲击充填节理岩石数值模型，记录冲击作用下充填节理层和岩石层的微裂纹数目以及微裂纹出现的位置数据；
损伤度计算模块，用于根据微裂纹数目以及微裂纹出现的位置数据计算冲击作用下产生的微裂纹在充填节理岩石数值模型任一平面的最大贯通度，当最大贯通度大于阈值时，冲击作用下充填节理岩石损伤度为最大贯通度数值；当充填节理岩石数值模型的最大贯通度小于阈值时充填节理岩石的损伤度为：
；
式中，D表示冲击作用下充填节理岩石损伤度，f1和f2分别表示落锤数值模型冲击充填节理岩石数值模型下充填节理层和岩石层的微裂纹数目，c1和c2分别表示充填节理层和岩石层的接触数目。
8.一种终端装置，其特征在于，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如权利要求1至6任一项所述的充填节理岩石损伤度离散元数值模拟方法的步骤。
9.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述的充填节理岩石损伤度离散元数值模拟方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311804203.1A CN117473797B (zh) | 2023-12-26 | 2023-12-26 | 充填节理岩石损伤度离散元数值模拟方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311804203.1A CN117473797B (zh) | 2023-12-26 | 2023-12-26 | 充填节理岩石损伤度离散元数值模拟方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117473797A true CN117473797A (zh) | 2024-01-30 |
| CN117473797B CN117473797B (zh) | 2024-03-08 |
Family
ID=89625931
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311804203.1A Active CN117473797B (zh) | 2023-12-26 | 2023-12-26 | 充填节理岩石损伤度离散元数值模拟方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117473797B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101433787B1 (ko) * | 2013-04-18 | 2014-08-27 | 한국원자력연구원 | 음향방출기법을 이용한 현장암반의 손상도 추정 시스템 및 그 추정 방법 |
| CN108629126A (zh) * | 2018-05-09 | 2018-10-09 | 中国地质大学（北京） | 一种考虑宏细观缺陷耦合的岩体力学数值建模方法 |
| CN111323562A (zh) * | 2020-04-16 | 2020-06-23 | 河海大学 | 一种建立充填裂隙岩石渗流损伤软化模型的方法 |
| CN116129035A (zh) * | 2022-11-28 | 2023-05-16 | 西南交通大学 | 考虑岩块空间分布的碎裂岩体三维离散元模型构建方法 |
-
2023
- 2023-12-26 CN CN202311804203.1A patent/CN117473797B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101433787B1 (ko) * | 2013-04-18 | 2014-08-27 | 한국원자력연구원 | 음향방출기법을 이용한 현장암반의 손상도 추정 시스템 및 그 추정 방법 |
| CN108629126A (zh) * | 2018-05-09 | 2018-10-09 | 中国地质大学（北京） | 一种考虑宏细观缺陷耦合的岩体力学数值建模方法 |
| CN111323562A (zh) * | 2020-04-16 | 2020-06-23 | 河海大学 | 一种建立充填裂隙岩石渗流损伤软化模型的方法 |
| CN116129035A (zh) * | 2022-11-28 | 2023-05-16 | 西南交通大学 | 考虑岩块空间分布的碎裂岩体三维离散元模型构建方法 |
Non-Patent Citations (2)
| Title |
|---|
| 刘云贺 等: "考虑损伤的平行黏结接触模型开发及其参数影响分析", 岩土力学, vol. 43, no. 3, 31 March 2022 (2022-03-31), pages 615 - 624 * |
| 宋浪 等: "充填节理特性对岩石静态压缩力学行为影响的颗粒流模拟", 中山大学学报（自然科学版）, vol. 62, no. 5, 19 June 2023 (2023-06-19), pages 145 - 156 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117473797B (zh) | 2024-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Celik | Determination of bruise susceptibility of pears (Ankara variety) to impact load by means of FEM-based explicit dynamics simulation | |
| Somfai et al. | Elastic wave propagation in confined granular systems | |
| Zhang et al. | Quasi-static fracturing in double-flawed specimens under uniaxial loading: the role of strain rate | |
| Paluszny et al. | Fracture and impulse based finite-discrete element modeling of fragmentation | |
| Resende et al. | Micromechanical modelling of stress waves in rock and rock fractures | |
| Hu et al. | 3D continuum-discrete coupled modelling of triaxial Hopkinson bar tests on rock under multiaxial static-dynamic loads | |
| Gu et al. | Discrete element modeling of shear wave propagation using bender elements in confined granular materials of different grain sizes | |
| Bu et al. | Evaluation of the characterization of acoustic emission of brittle rocks from the experiment to numerical simulation | |
| CN108827455B (zh) | 一种节理岩体爆破振动衰减参数预测方法及装置 | |
| Ning et al. | Shear wave propagation in granular assemblies | |
| CN112883617A (zh) | 隧道衬砌监测范围计算方法、装置、设备及可读存储介质 | |
| Zhang et al. | Three-dimensional reconstruction of block shape irregularity and its effects on block impacts using an energy-based approach | |
| Wang et al. | Impact source localization and vibration intensity prediction on construction sites | |
| CN117473797B (zh) | 充填节理岩石损伤度离散元数值模拟方法及系统 | |
| Pakzad et al. | Three‐dimensional finite element simulation of fracture propagation in rock specimens with pre‐existing fissure (s) under compression and their strength analysis | |
| Matsuyama et al. | Power law statistics of force and acoustic emission from a slowly penetrated granular bed | |
| Gresil et al. | Benchmark problems for predictive fem simulation of 1-D and 2-D guided waves for structural health monitoring with piezoelectric wafer active sensors | |
| Yuan et al. | Theoretical analysis of longitudinal wave attenuation in a stressed rock with variable cross-section | |
| Pasqual et al. | Low strain integrity tests in piles–1-D and 3-D numerical modeling and comparisons with results obtained in the field | |
| Wang et al. | A Novel $ A~ Priori $ Simulation Algorithm for Absorbing Receivers in Diffusion-Based Molecular Communication Systems | |
| Watson et al. | COMSOL modelling of the acoustoelastic effect | |
| Cui et al. | Deterministic failure prediction of toughened glass when impacted by ice | |
| Cusatis et al. | 3D Lattice model for dynamic simulations of creep, fracturing and rate effect in concrete | |
| CN117875096A (zh) | 岩石应力波载荷传播特性模拟方法及相关装置 | |
| Potyondy | A bonded-disk model for rock: Relating microproperties and macroproperties |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |