CN107122542B - 基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法 - Google Patents

基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法 Download PDF

Info

Publication number
CN107122542B
CN107122542B CN201710281664.3A CN201710281664A CN107122542B CN 107122542 B CN107122542 B CN 107122542B CN 201710281664 A CN201710281664 A CN 201710281664A CN 107122542 B CN107122542 B CN 107122542B
Authority
CN
China
Prior art keywords
ore
rock
particle
model
rupture
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.)
Active
Application number
CN201710281664.3A
Other languages
English (en)
Other versions
CN107122542A (zh
Inventor
金爱兵
孙金海
孙浩
郭瑛玮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN201710281664.3A priority Critical patent/CN107122542B/zh
Publication of CN107122542A publication Critical patent/CN107122542A/zh
Application granted granted Critical
Publication of CN107122542B publication Critical patent/CN107122542B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume, or surface-area of porous materials
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/006—Crack, flaws, fracture or rupture
    • G01N2203/0067—Fracture or rupture

Landscapes

  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，所述方法包括：根据椭圆及其内接多边形的几何关系，采用椭圆的内接多边形模拟实际不规则矿岩散体形状，利用颗粒流数值软件构建不规则颗粒模型；根据已知放矿室内试验中所用散体材料的物理力学性质，研究松动体和放出体的形态变化规律,确定放矿数值模型及颗粒细观力学参数；通过改变颗粒间的平行粘结强度，进行放矿数值试验确定平行粘结强度和矿岩二次破裂的关系；验证不规则颗粒模型以及放矿数值模型的准确性。不仅实现了散体颗粒在形状上与实际矿岩散体更为接近，而且满足了矿岩能够发生二次破裂的特性，提高了放矿数值模型研究的准确性。

Description

基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法
技术领域
本发明涉及地质研究技术领域，特别是指一种基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法。
背景技术
实际矿山中的矿岩形状不规则，且在放矿过程中会由于相互之间的挤压和剪切而出现二次破裂现象，形成新的不规则散体。而实际的放矿室内试验模型，仅能达到几何相似，不能做到力学相似，因此，试验中散体材料不会出现二次破裂现象。此外，在基于离散单元法的放矿数值研究中，绝大多数是采用球形颗粒模型或双颗粒模型(也被称为花生模型，Peanut Model)。上述两种模型存在如下两个缺点：(1)可靠性不强：不能很好地反映实际矿岩形状的不规则性，从而不能实现矿岩间的较高的内锁力(Inter locking Force)；(2)适用性不广：花生模型破裂后变成两个独立的球形颗粒，而球形颗粒本身为刚体，无法破裂。
发明内容
本发明要解决的技术问题是提供一种基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，能够准确模拟发生二次破裂的特征。
为解决上述技术问题，本发明的实施例提供一种基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，所述方法包括：
根据椭圆及其内接多边形的几何关系，采用椭圆的内接多边形模拟实际不规则矿岩散体形状，利用颗粒流数值软件构建具有不规则颗粒模型；
根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数；
通过改变颗粒间的平行粘结强度，进行放矿数值试验确定平行粘结强度和矿岩二次破裂的关系；
验证不规则颗粒模型以及矿数值模型的准确性。
优选的，所述根据椭圆及其内接多边形的几何关系，采用椭圆的内接多边形模拟实际不规则矿岩散体，包括：
将不规则颗粒模型实际面积Sclu表示为：
将不规则颗粒模型的等效半径rclu表示为：
其中，Sp表示内接多边形的面积，ρ为不规则颗粒模型即内接多边形内的空隙率，a和b分别表示椭圆的半长轴和半短轴，θk表示多边形的第k个内角，e表示椭圆的离心率，rcir表示椭圆外接圆的半径。
优选的，所述利用颗粒流数值软件构建不规则颗粒模型，包括：
根据rcir和rclu的关系，生成不同半径的圆形颗粒；
根据内接多边形的顶点坐标，编写程序生成表示内接多边形的墙体，并删除圆形颗粒；
根据既定空隙率ρ，在每个内接多边形内生成一定数量的圆形小颗粒并使其相互粘结，将所有墙体删除，圆形小颗粒间采用的粘结模型为平行粘结模型。
优选的，所述根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数，包括：
在颗粒流数值软件构建单口放矿模型，根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数，所述细观力学参数包括密度、摩擦系数、形状、粒径分布、刚度、平行粘结强度。
优选的，所述根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数，包括：
依据已知放矿室内试验中所用的砾石材料的宏观物理力学性质调试并确定数值模型中散体颗粒的细观力学参数。
优选的，所述依据已知放矿室内试验中所用的砾石材料的宏观物理力学性质调试并确定数值模型中散体颗粒的细观力学参数，包括：
利用下表设置颗粒密度及颗粒摩擦系数：
在不规则颗粒模型将模型中的最小圆形颗粒半径设为60mm，不规则颗粒模型的长宽比范围设为2.0～2.5，设置颗粒粒径分布为(150～1350mm)，相似比为1:1，将颗粒的法向及切向刚度设为1.0e8N/m，通过双轴压缩试验测量颗粒的宏观散体摩擦角，并与实际散体摩擦角对比，从而确定不规则颗粒模型中颗粒间的平行粘结强度。
优选的，所述通过改变颗粒间的平行粘结强度，进行放矿数值试验确定平行粘结强度和矿岩二次破裂的关系，包括：
生成符合粒径分布的圆形颗粒，将其替换为花生模型或不规则颗粒模型Cluster模型，并赋予指定的细观力学参数；
当模型中颗粒在重力作用下达到自然平衡状态后，删除放矿口的颗粒，使模型中颗粒能够不断从放矿口放出；
在整个模拟过程中，记录每个放出颗粒的初始位置、ID编号、半径等信息，以用于圈定放出体；
改变颗粒间的平行粘结强度，重复上述步骤，确定平行粘结强度和矿岩二次破裂的关系。
优选的，所述生成符合粒径分布的圆形颗粒，将其替换为花生模型或不规则颗粒模型，并赋予指定的细观力学参数，包括：
当圆形颗粒半径大于等于150mm时，其将被替换为不规则颗粒模型；而当圆形颗粒半径小于150mm时，其将被等效替换为花生模型。
优选的，所述当模型中颗粒在重力作用下达到自然平衡状态后，删除放矿口的颗粒，使模型中颗粒能够不断从放矿口放出，包括：
将垂直位移大于1m的颗粒所圈定的范围定义为松动体。
优选的，所述在整个模拟过程中，记录每个放出颗粒的初始位置、ID编号、半径等信息，以用于圈定放出体，包括：
将通过颗粒流数值软件中的FISH程序记录发生的每个粘结破裂事件，并通过离散裂隙网络模块生成小裂隙的方式存储粘结破裂事件发生的位置、大小和方向。
本发明的上述技术方案的有益效果如下：
上述方案中，不仅实现了散体颗粒在形状上与实际矿岩散体更为接近，而且满足了矿岩能够发生二次破裂的特性，提高了放矿数值模型研究的准确性。
附图说明
图1为本发明的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法流程图。
图2为本发明的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法椭圆内接多边形示意图。
具体实施方式
为使本发明要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。
如图1所示，本发明的实施例提供一种基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，所述方法包括：
步骤101：根据椭圆及其内接多边形的几何关系，采用椭圆的内接多边形模拟实际不规则矿岩散体，利用颗粒流数值软件构建不规则颗粒模型；
步骤102：根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数；
步骤103：通过改变颗粒间的平行粘结强度，进行放矿数值试验确定平行粘结强度和矿岩二次破裂的关系；
步骤104：验证不规则颗粒模型以及放矿数值模型的准确性。
本发明实施例的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，不仅实现了散体颗粒在形状上与实际矿岩散体更为接近，而且满足了矿岩能够发生二次破裂的特性，提高了放矿数值模型研究的准确性。
优选的，所述根据椭圆及其内接多边形的几何关系，采用椭圆的内接多边形模拟实际不规则矿岩散体，包括：
将不规则颗粒模型实际面积Sclu表示为：
将不规则颗粒模型的等效半径rclu表示为：
其中，Sp表示内接多边形的面积，ρ为不规则颗粒模型即内接多边形内的空隙率，a和b分别表示椭圆的半长轴和半短轴，θk表示多边形的第k个内角，e表示椭圆的离心率，rcir表示椭圆外接圆的半径。
具体的，O表示外接圆的圆心，θk表示多边形的第k个内角，(xk,yk)表示内接多边形的第k个顶点坐标，α表示椭圆长轴与x坐标轴的夹角，其变化范围为0～2π，用于改变椭圆的布置方向。
图中，(xO,yO)表示外接圆的圆心坐标，a和b分别表示椭圆的半长轴和半短轴。
椭圆的面积Se是：
Se＝πab
半径为rcir的外接圆的面积Scir是：
椭圆的内接多边形的面积Sp是：
优选的，所述利用颗粒流数值软件构建不规则颗粒模型，包括：
根据rcir和rclu的关系，生成不同半径的圆形颗粒；
根据内接多边形的顶点坐标，编写程序生成表示内接多边形的墙体，并删除圆形颗粒；
根据既定空隙率ρ，在每个内接多边形内生成一定数量的圆形小颗粒并使其相互粘结，将所有墙体删除，圆形小颗粒间采用的粘结模型为平行粘结模型。
优选的，所述根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数，包括：
在颗粒流数值软件构建单口放矿模型，根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及细观力学参数，所述细观力学参数包括密度、摩擦系数、形状、粒径分布、刚度、平行粘结强度。
具体的，在PFC2D中构建尺寸为40m×50m(宽×高)单口放矿模型。通过被固定住的直径为0.7m的黑色颗粒构成模型的侧墙和底墙，用来模拟实际采场中较为粗糙的边壁；通过底墙上的中央颗粒模拟实际的放矿口。
优选的，所述根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数，包括：
依据已知放矿室内试验中所用的砾石材料的宏观物理力学性质调试并确定数值模型中散体颗粒的细观力学参数。
优选的，所述依据已知放矿室内试验中所用的砾石材料的宏观物理力学性质调试并确定数值模型中散体颗粒的细观力学参数，包括：
利用下表设置颗粒密度及颗粒摩擦系数：
在不规则颗粒模型将模型中的最小圆形颗粒半径设为60mm，不规则颗粒模型的长宽比范围设为2.0～2.5，设置颗粒粒径分布为(150～1350mm)，相似比为1:1，将颗粒的法向及切向刚度设为1.0e8N/m，通过双轴压缩试验测量颗粒的宏观散体摩擦角，并与实际散体摩擦角对比，从而确定不规则颗粒模型中颗粒间的平行粘结强度。
具体的本模拟将依据学者Castro所做放矿室内试验中所用的砾石材料的宏观物理力学性质调试并确定数值模型中散体颗粒的细观力学参数。Castro中所用的试验材料是来自澳大利亚布里斯班Mount Coutha采石场的千枚岩砾石，其主要物理力学参数如上表所示。
(1)颗粒密度及摩擦系数的取值与表1中相应参数一致。
(2)颗粒使用不规则颗粒模型以满足散体形状的不规则性。每个不规则颗粒模型中的最小圆形颗粒半径设为60mm，Cluster模型的长宽比范围设为2.0～2.5。
(3)不规则颗粒模型中颗粒粒径分布(150～1350mm)基本与Castro模型中颗粒粒径分布(70.8～1350mm)保持一致，只是出于计算效率的考虑，最小颗粒粒径设为150mm。需要说明的是，不规则颗粒模型相似比为1:1，Castro模型相似比为1:30，因此不规则颗粒模型中的颗粒粒径为扩大30倍后的结果。
(4)根据之前学者的研究，颗粒刚度仅对模拟速度有影响，而对模拟结果无影响。因此，参考之前学者的取值，将颗粒的法向及切向刚度设为1.0e8N/m。
(5)通过双轴压缩试验测量颗粒的宏观散体摩擦角，并与实际散体摩擦角对比，从而确定不规则颗粒模型中颗粒间的平行粘结强度。散体摩擦角φ反映材料抵抗变形的能力。随着应变逐渐增大，散体试样也将逐渐达到峰值强度σ1并逐渐降低。可以根据摩尔-库伦理论，通过峰值强度σ1与围压σ3的关系计算散体摩擦角φ。
散体试样的初始宽高比为1:2。围压σ3设为5.0e5MPa，应变率保持0.005/s直至最终应变达到0.2。平行粘结强度取值范围为1.0MPa～1.0e4MPa，其余细观参数均保持不变。
一般而言，崩落矿岩在峰值强度时的散体摩擦角范围是35～45°。当不规则颗粒模型中颗粒间的平行粘结强度大于3MPa时，其散体摩擦角与实际矿岩摩擦角相符。因此，在本次放矿试验研究中，将选取如下三个取值：3.0MPa、5.0MPa和1.0e4MPa。需要说明的是，当平行粘结强度为3.0MPa和5.0MPa时，Cluster颗粒模型可以发生破裂；而当平行粘结强度取值足够大，如1.0e4MPa，颗粒模型则不会发生破裂，从而用于与Castro的放矿室内试验结果进行对比，以验证不规则颗粒模型以及放矿模型的可靠性。
优选的，所述通过改变颗粒间的平行粘结强度，进行放矿数值试验确定平行粘结强度和矿岩二次破裂的关系，包括：
生成符合粒径分布的圆形颗粒，将其替换为花生模型或不规则颗粒模型，并赋予指定的细观力学参数；
当模型中颗粒在重力作用下达到自然平衡状态后，删除放矿口的颗粒，使模型中颗粒能够不断从放矿口放出；
在整个模拟过程中，记录每个放出颗粒的初始位置、ID编号、半径等信息，以用于圈定放出体；
改变颗粒间的平行粘结强度，重复上述步骤，确定平行粘结强度和矿岩二次破裂的关系。
优选的，所述生成符合粒径分布的圆形颗粒，将其替换为花生模型或不规则颗粒模型，并赋予指定的细观力学参数，包括：
当圆形颗粒半径大于等于150mm时，其将被替换为不规则颗粒模型；而当圆形颗粒半径小于150mm时，其将被等效替换为花生模型。
优选的，所述当模型中颗粒在重力作用下达到自然平衡状态后，删除放矿口的颗粒，使模型中颗粒能够不断从放矿口放出，包括：
将垂直位移大于1m的颗粒所圈定的范围定义为松动体。
优选的，所述在整个模拟过程中，记录每个放出颗粒的初始位置、ID编号、半径等信息，以用于圈定放出体，包括：
将通过颗粒流数值软件中的FISH程序记录发生的每个粘结破裂事件，并通过离散裂隙网络模块生成小裂隙的方式存储粘结破裂事件发生的位置、大小和方向。
本发明实施例的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，当平行粘结强度取值为1.0e4MPa，即颗粒不发生破裂时，其松动体和放出体的高度与宽度曲线与Castro物理试验结果十分符合，均满足幂函数关系，拟合相关系数均接近于1，较好地验证了不规则颗粒模型以及PFC放矿模拟方法的准确性。
模拟结果中，所有松动体和放出体形状均近似一个倒置的水滴，具有上下不对称性，符合放矿水滴理论(Upside-down drop shape theory)。此外，松动体和放出体的宽度会随着粘结强度的减小而增大。因此，对于矿岩强度较小的崩落法矿山，可以适当增大放矿口间距。
矿岩强度越小，二次破裂后的平均粒径越小。矿岩破裂主要发生在剪切区域，即松动体区域和未松动区域之间的过渡区域。该结果也验证了之前学者的结论：崩落矿岩的流动主要受松动体顶部应力拱塌落以及松动体四周矿岩间的相互挤压、剪切的控制。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述方法包括：
根据椭圆及其内接多边形的几何关系，采用椭圆的内接多边形模拟实际不规则矿岩散体，利用颗粒流数值软件构建不规则颗粒模型；
根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数；
通过改变颗粒间的平行粘结强度，研究松动体和放出体的形态变化规律确定平行粘结强度和矿岩二次破裂的关系；
验证不规则颗粒模型以及放矿数值模型的准确性。
2.根据权利要求1所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述根据椭圆及其内接多边形的几何关系，采用椭圆的内接多边形模拟实际不规则矿岩散体形状，包括：
将不规则颗粒模型实际面积Sclu表示为：
将不规则颗粒模型的等效半径rclu表示为：
其中，Sp表示内接多边形的面积，ρ为不规则颗粒模型即内接多边形内的空隙率，a和b分别表示椭圆的半长轴和半短轴，θk表示多边形的第k个内角，e表示椭圆的离心率，rcir表示椭圆外接圆的半径。
3.根据权利要求2所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述利用颗粒流数值软件构建不规则颗粒模型，包括：
根据rcir和rclu的关系，生成不同半径的圆形颗粒；
根据内接多边形的顶点坐标，编写程序生成表示内接多边形的墙体，并删除圆形颗粒；
根据既定空隙率ρ，在每个内接多边形内生成一定数量的圆形小颗粒并使其相互粘结，将所有墙体删除，圆形小颗粒间采用的粘结模型为平行粘结模型。
4.根据权利要求1所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数，包括：
利用颗粒流数值软件构建单口放矿模型，根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数，所述细观力学参数包括密度、摩擦系数、形状、粒径分布、刚度、平行粘结强度。
5.根据权利要求4所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述根据已知放矿室内试验中所用散体材料的物理力学性质，确定放矿数值模型及颗粒细观力学参数，包括：
依据已知放矿室内试验中所用的砾石材料的宏观物理力学性质调试并确定数值模型中散体颗粒的细观力学参数。
6.根据权利要求5所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述依据已知放矿室内试验中所用的砾石材料的宏观物理力学性质调试并确定数值模型中散体颗粒的细观力学参数，包括：
利用下表设置颗粒密度及颗粒摩擦系数：
在不规则颗粒模型将模型中的最小圆形颗粒半径设为60mm，不规则颗粒模型的长宽比范围设为2.0～2.5，设置颗粒粒径分布为(150～1350mm)，相似比为1:1，将颗粒的法向及切向刚度设为1.0e8N/m，通过双轴压缩试验测量颗粒的宏观散体摩擦角，并与实际散体摩擦角对比，从而确定不规则颗粒模型中颗粒间的平行粘结强度。
7.根据权利要求1所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述通过改变颗粒间的平行粘结强度，进行放矿数值试验确定平行粘结强度和矿岩二次破裂的关系，包括：
生成符合粒径分布的圆形颗粒，将其替换为花生模型或不规则颗粒模型Cluster模型，并赋予指定的细观力学参数；
当模型中颗粒在重力作用下达到自然平衡状态后，删除放矿口的颗粒，使模型中颗粒能够不断从放矿口放出；
在整个模拟过程中，记录每个放出颗粒的初始位置、ID编号、半径等信息，以用于圈定放出体；
改变颗粒间的平行粘结强度，重复上述步骤，确定平行粘结强度和矿岩二次破裂的关系。
8.根据权利要求7所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述生成符合粒径分布的圆形颗粒，将其替换为花生模型或不规则颗粒模型，并赋予指定的细观力学参数，包括：
当圆形颗粒半径大于等于150mm时，其将被替换为不规则颗粒模型；而当圆形颗粒半径小于150mm时，其将被等效替换为花生模型。
9.根据权利要求7所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述当模型中颗粒在重力作用下达到自然平衡状态后，删除放矿口的颗粒，使模型中颗粒能够不断从放矿口放出，包括：
将垂直位移大于1m的颗粒所圈定的范围定义为松动体。
10.根据权利要求7所述的基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法，其特征在于，所述在整个模拟过程中，记录每个放出颗粒的初始位置、ID编号、半径等信息，以用于圈定放出体，包括：
将通过颗粒流数值软件中的FISH程序记录发生的每个粘结破裂事件，并通过离散裂隙网络模块生成小裂隙的方式存储粘结破裂事件发生的位置、大小和方向。
CN201710281664.3A 2017-04-26 2017-04-26 基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法 Active CN107122542B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710281664.3A CN107122542B (zh) 2017-04-26 2017-04-26 基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710281664.3A CN107122542B (zh) 2017-04-26 2017-04-26 基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法

Publications (2)

Publication Number Publication Date
CN107122542A CN107122542A (zh) 2017-09-01
CN107122542B true CN107122542B (zh) 2019-07-09

Family

ID=59726482

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710281664.3A Active CN107122542B (zh) 2017-04-26 2017-04-26 基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法

Country Status (1)

Country Link
CN (1) CN107122542B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108804792B8 (zh) * 2018-05-30 2020-09-22 中南大学 矿岩崩落过程模拟方法及系统
CN110083936B (zh) * 2019-04-26 2023-06-16 中煤科工集团重庆研究院有限公司 考虑形状参数的可破碎碎石颗粒生成方法
CN112730026B (zh) * 2020-12-21 2021-10-15 北京科技大学 一种矿岩散体破裂时空演化规律分析系统及方法
CN113553753B (zh) * 2021-08-10 2024-02-06 中国电建集团华东勘测设计研究院有限公司 破裂识别方法、装置、存储介质及设备
CN113899661B (zh) * 2021-11-08 2022-12-27 北京科技大学 基于矿岩颗粒体系结构特征的细小颗粒穿流特性研究方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1693662A (zh) * 2005-05-23 2005-11-09 东北大学 一种改进的无底柱分段崩落法
CN106246183A (zh) * 2016-07-25 2016-12-21 辽宁科技大学 一种确定无底柱分段崩落法崩落体形态的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1693662A (zh) * 2005-05-23 2005-11-09 东北大学 一种改进的无底柱分段崩落法
CN106246183A (zh) * 2016-07-25 2016-12-21 辽宁科技大学 一种确定无底柱分段崩落法崩落体形态的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
崩落法采矿中放出体流动特性的影响因素;孙浩 等;《工程科学学报》;20150930;第37卷(第9期);全文

Also Published As

Publication number Publication date
CN107122542A (zh) 2017-09-01

Similar Documents

Publication Publication Date Title
CN107122542B (zh) 基于矿岩破裂的松动体、放出体及崩落矿岩形状的模拟方法
CN107423498B (zh) 一种高致密度离散颗粒多相体系的建模方法
Laubie et al. Stress transmission and failure in disordered porous media
CN106960070B (zh) 一种基于有限元-离散元ct重构煤体的渗流模拟方法
Mollon et al. Fourier–Voronoi-based generation of realistic samples for discrete modelling of granular materials
CN104899393B (zh) 一种离散相增强复合材料细观结构的建模方法
Cil et al. DEM modeling of grain size effect in brittle granular soils
Tang et al. Effects of statistics of cell’s size and shape irregularity on mechanical properties of 2D and 3D Voronoi foams
CN109543290B (zh) 一种深水气井防砂筛管冲蚀数值模拟方法
CN108491639B (zh) 基于纳米颗粒碰撞接触模型的封堵页岩孔喉仿真模拟方法
Fang et al. DEM study on the microscale and macroscale shear behaviours of granular materials with breakable and irregularly shaped particles
CN106290082A (zh) 一种离散元组合颗粒及其离散元堆积试验模拟方法
CN110390130A (zh) 含水合物沉积物的降压开采出砂室内实验数值模拟方法
CN112199843A (zh) 基于颗粒结构的岩石破裂演化细观模型、方法及电子设备
Li et al. Investigation of sand production mechanisms using DEM with fluid flow
CN110152771A (zh) 物料破碎腔结构及多级嵌套式物料破碎腔结构设计方法
Wan et al. Vibration response analysis of the tail beam of hydraulic support impacted by coal gangue particles with different shapes
CN107084895B (zh) 一种基于块石形态的土石混合体击实试验模拟方法及系统
Li et al. Performance optimization of banana vibrating screens based on PSO-SVR under DEM simulations
CN116050137A (zh) 一种基于爆炸震源模型的露天矿山爆破块度预测方法
Zhou et al. Micromechanics damage modeling of brittle rock failure processes under compression
CN109916760A (zh) 一种确定不同最大粒径下堆石料最优细粒含量的方法
Remond et al. Simulation of the packing of granular mixtures of non-convex particles and voids characterization
CN110083936B (zh) 考虑形状参数的可破碎碎石颗粒生成方法
Mahboob et al. Three-dimensional simulation of granular materials by discrete element method (DEM) by considering the fracture effect of particles

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