CN109359416A - 一种反映真实工程土石混合体分布的颗粒流数值模拟方法 - Google Patents

一种反映真实工程土石混合体分布的颗粒流数值模拟方法 Download PDF

Info

Publication number
CN109359416A
CN109359416A CN201811321036.4A CN201811321036A CN109359416A CN 109359416 A CN109359416 A CN 109359416A CN 201811321036 A CN201811321036 A CN 201811321036A CN 109359416 A CN109359416 A CN 109359416A
Authority
CN
China
Prior art keywords
soil
stone
rock mixture
true
numerical simulation
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
Application number
CN201811321036.4A
Other languages
English (en)
Other versions
CN109359416B (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.)
Hebei University of Technology
Original Assignee
Hebei University of Technology
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 Hebei University of Technology filed Critical Hebei University of Technology
Priority to CN201811321036.4A priority Critical patent/CN109359416B/zh
Publication of CN109359416A publication Critical patent/CN109359416A/zh
Application granted granted Critical
Publication of CN109359416B publication Critical patent/CN109359416B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Astronomy & Astrophysics (AREA)
  • Remote Sensing (AREA)
  • Multimedia (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种反映真实工程土石混合体分布的颗粒流数值模拟方法,该方法包括以下步骤:拍照获取工程现场土石混合体的真实图像;图像处理技术识别石块并提取其边界坐标;结合颗粒筛分及图像数据获得土石混合体的颗粒级配;试验测试土体、石块的强度及变形等参数;建立土石混合体颗粒流数值模型;标定颗粒流数值模型参数,开展土石混合体试件力学性质数值试验。本方法实施方便,可程序化智能建模,可克服试验尺寸限制,精度可满足工程要求,适用于土石混合体、混凝土等二元介质离散元模拟。

Description

一种反映真实工程土石混合体分布的颗粒流数值模拟方法
技术领域
本发明涉及岩土工程及地质工程领域,具体是一种反映真实工程土石混合体分布的颗粒流数值模拟方法。
背景技术
土石混合体是地表第四纪沉积物的主要类型之一,特别是我国西部山区分布十分广泛,如常见的残积物、坡积物及洪积物等。另外,西南山区的填方工程十分普遍,如公路、铁路路基及机场等填方工程的填筑料大多为土石混合体。由于岩石和土体的物理力学性质差异性非常大,而且自然界土石混合体中石块的形态及大小多种多样、空间分布也非常杂乱,难以通过数学等理论方法描述土石混合体的结构特征。
由于土石混合体中含有大量不同尺度的石块,其剪切、压缩等力学性质很难通过室内试验的方法获得。现今工程界常通过现场试验开展土石混合体的力学性质测试,然而现场试验存在两个明显的不足:1)试件扰动性强。现场试验试件准备常通过开挖及切削等方法实现,当土样中含有较多不规则的石块时,试件制备过程中容易出现坑洞,很难形成平整的加载面。2)费用高,工时长,影响工程进度,试验结果离散。由于现场试验试件及试验准备工时长且需要大量的设备及人力,使得现场测试费用高,且占用现场场地而影响工程进度。因此现场试验样本数量一般很少,且土石混合体的力学性质测试结果常离散性大。而室内试验受用料的最大粒径限制,需要对超过试验容许最大粒径的土石料进行必要的缩尺处理,受缩尺效应影响,室内实验结果与真实性质存在一定差异。
发明内容
针对现有技术的不足,本发明拟解决的技术问题是,提供一种反映真实工程土石混合体分布的颗粒流数值模拟方法。
本发明解决所述技术问题的技术方案是,提供一种反映真实工程土石混合体分布的颗粒流数值模拟方法,其特征在于该方法包括以下步骤:
步骤1、根据工程现场的土石混合体的物质组成特征,选取代表性采样窗口,通过正面拍照及尺寸标定技术获得采样窗口中的土石混合体的真实图像;采用土石分界粒径的理论公式获得真实图像中土石混合体的土石分界粒径;
步骤2、根据真实图像利用二值图像处理技术得到二值图像,识别二值图像中的石块并提取石块边界轮廓关键点像元坐标;
步骤3、在工程现场选取真实土石混合体,通过颗粒筛分,确定小于土石分界粒径的土体级配特征;对步骤2得到的二值图像进行分析,提取石块面积和石块等效尺寸,确定大于土石分界粒径的石块级配特征;再将两种级配特征组合,得到土石混合体的整个级配曲线;
步骤4、采用经过步骤3筛分后的小于土石分界粒径的土石混合体制备重塑土试样,并对重塑土进行力学试验,得到重塑土的强度及变形参数;根据步骤2得到的石块边界轮廓关键点像元坐标,在PFC软件中对石块进行重构,进而在石块中填充石块材料颗粒单元,从而建立土石混合体中石块数值模型;再在石块与石块之间的空隙填充土体材料颗粒单元,得到土石混合体的颗粒流数值模型;
步骤5、选取颗粒接触本构模型,并通过数值试验标定土体和石块的颗粒接触细观参数,开展土石混合体试件的压缩和剪切数值试验。
与现有技术相比,本发明有益效果在于:
1、适用范围广:本方法的核心贡献是建立了快速准确地提取二元复杂结构体的几何信息方法及将二元结构(如土石混合体中土体和石块)一分为二的研究思路,适用于土石混合体、混凝土等二元结构的数值建模及力学性质研究。
2、快速且低成本:本方法采用图像信息建模及颗粒流数值模拟的方法,既不影响工程施工也无现场测试等费用,图像处理及数值模拟技术均为计算机程序化,因此快速且低成本。
3、高精度大样本仿真模拟:本方法通过现场照片的图像识别技术确定石块边界;将分别试验测试的土体、石块物理力学参数应用于数值模拟体现了二元结构的物理力学性质的差异性;数值模拟可以开展大量样本的数值试验,有利于更加准确地统计其力学性质。
4、本方法所述级配获取方法避免了现场大型颗分试验,仅需现场拍照即可,避免了现场大型颗分试验的高成本以及试验对施工过程的影响。
5、采用大尺度数值模型研究大型土石混合体试件的力学性质,克服了室内试验试件尺寸的限制,避免了缩尺效应。
6、本方法实施方便,可程序化智能建模,可克服试验尺寸限制,精度可满足工程要求,适用于土石混合体、混凝土等二元介质离散元模拟。
附图说明
图1为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的采样窗口中的土石混合体的真实图像;
图2为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的采样窗口中的土石混合体的二值图像;
图3为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的二值图像石块边界示意图;
图4为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的石块数值模型图;
图5为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的土石混合体的颗粒流数值模型图;
图6为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的石块级配曲线图;
图7为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的土体级配曲线图;
图8为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的土石混合体试件在受力达到应力-应变曲线屈服点时的图像;
图9为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的土石混合体试件在受力达到应力-应变曲线峰值点时的图像;
图10为本发明反映真实工程土石混合体分布的颗粒流数值模拟方法一种实施例的土石混合体试件在受力达到应力-应变曲线峰后段时的图像;
(图中:1、石块;2、土体;3、采样窗口;4、石块边界;5、石块材料颗粒单元;6、土体材料颗粒单元;7、十字标尺)
具体实施方式
下面给出本发明的具体实施例。具体实施例仅用于进一步详细说明本发明,不限制本申请权利要求的保护范围。
本发明提供了一种反映真实工程土石混合体分布的颗粒流数值模拟方法(简称方法),其特征在于该方法包括以下步骤:
步骤1、根据工程现场的土石混合体的物质组成特征,对工程现场进行调查测绘,必要时可挖探坑,选取代表性采样窗口,通过正面拍照及尺寸标定技术获得采样窗口中的土石混合体的真实图像(如图1);采用土石分界粒径的理论公式获得真实图像中土石混合体的土石分界粒径;
步骤2、根据真实图像利用二值图像处理技术得到二值图像,识别二值图像中的石块并提取石块边界轮廓关键点像元坐标(如图2和3);
步骤3、在工程现场选取真实土石混合体,通过颗粒筛分,确定小于土石分界粒径的土体级配特征(如图7,包括了5个采样窗口的数据);对步骤2得到的二值图像进行分析,提取石块面积和石块等效尺寸,确定大于土石分界粒径的石块级配特征(如图6,包括了8个采样窗口的数据,发现石块级配大致可以分为三组);再将通过不同途径得到的两种级配特征组合,得到土石混合体的整个级配曲线;
步骤4、采用经过步骤3筛分后的小于土石分界粒径的土石混合体制备重塑土试样,并对重塑土进行力学试验,得到重塑土的强度及变形参数;根据步骤2得到的石块边界轮廓关键点像元坐标,在PFC(Particle Flow Code)软件中对石块进行重构,进而在石块中填充石块材料颗粒单元5,从而建立土石混合体中石块数值模型(如图4);再在石块与石块之间的空隙填充土体材料颗粒单元6,得到土石混合体的颗粒流数值模型(如图5);
步骤5、选取颗粒接触本构模型,并通过数值试验标定颗粒接触本构模型的土体和石块的颗粒接触细观参数,开展土石混合体试件的压缩和剪切数值试验等模拟。
步骤1中,正面拍照是指相机镜头垂直于拍照面;尺寸标定技术采用互相正交的十字标尺7平放于采样面内进行标定;采样窗口3为长方形或正方形,表面尽量平整,采样窗口特征尺寸大于最大石块尺寸的10倍;当采样窗口为长方形时,特征尺寸为短边长度;当采样窗口为正方形时,特征尺寸为边长;土石分界粒径的理论公式为d0=0.05Lc,d0为土石分界粒径,Lc为采样窗口(即模型区域)特征尺寸。
步骤2中,因为石块1和土体2具有不同的灰度值,利于二值图像中像素点的像素坐标和灰度值提取石块边界4;通过像素点间距与图像中的颗粒边界真实尺寸的比例关系式完成像素尺寸到真实尺寸的转化,统计即可得到石块在采样窗口内所占面积(即含量);
步骤3中,石块等效尺寸是指二值图像中石块的等面积圆的直径。
步骤4中,小于土石分界粒径的重塑土试样应与工程实际的密度一致,通过重塑土剪切、三轴压缩与土石混合体中所包含的各类石块的剪切及单轴压缩试验测试土体和石块的剪切强度、压缩强度、变形参数和其变形破坏演化过程。
步骤5中,选择颗粒接触本构模型采用可传递力和力矩的平行粘接模型;需要标定的细观参数包括颗粒密度、颗粒粒径、接触刚度、粘结强度及摩擦系数。如图8-图10显示了在法向应力800kPa作用下直剪试件的剪切破坏过程(图中空白处代表已破坏),并且测得其抗剪强度参数为摩擦角41.3°、粘聚力73.4kPa。
本发明未述及之处适用于现有技术。

Claims (7)

1.一种反映真实工程土石混合体分布的颗粒流数值模拟方法,其特征在于该方法包括以下步骤:
步骤1、根据工程现场的土石混合体的物质组成特征,选取代表性采样窗口,通过正面拍照及尺寸标定技术获得采样窗口中的土石混合体的真实图像;采用土石分界粒径的理论公式获得真实图像中土石混合体的土石分界粒径;
步骤2、根据真实图像利用二值图像处理技术得到二值图像,识别二值图像中的石块并提取石块边界轮廓关键点像元坐标;
步骤3、在工程现场选取真实土石混合体,通过颗粒筛分,确定小于土石分界粒径的土体级配特征;对步骤2得到的二值图像进行分析,提取石块面积和石块等效尺寸,确定大于土石分界粒径的石块级配特征;再将两种级配特征组合,得到土石混合体的整个级配曲线;
步骤4、采用经过步骤3筛分后的小于土石分界粒径的土石混合体制备重塑土试样,并对重塑土进行力学试验,得到重塑土的强度及变形参数;根据步骤2得到的石块边界轮廓关键点像元坐标,在PFC软件中对石块进行重构,进而在石块中填充石块材料颗粒单元,从而建立土石混合体中石块数值模型;再在石块与石块之间的空隙填充土体材料颗粒单元,得到土石混合体的颗粒流数值模型;
步骤5、选取颗粒接触本构模型,并通过数值试验标定土体和石块的颗粒接触细观参数,开展土石混合体试件的压缩和剪切数值试验。
2.根据权利要求1所述的反映真实工程土石混合体分布的颗粒流数值模拟方法,其特征在于步骤1中,正面拍照是指相机镜头垂直于拍照面;尺寸标定技术采用互相正交的十字标尺平放于采样面内进行标定。
3.根据权利要求1所述的反映真实工程土石混合体分布的颗粒流数值模拟方法,其特征在于步骤1中,采样窗口为长方形或正方形,表面平整,采样窗口特征尺寸大于最大石块尺寸的10倍;当采样窗口为长方形时,特征尺寸为短边长度;当采样窗口为正方形时,特征尺寸为边长。
4.根据权利要求1所述的反映真实工程土石混合体分布的颗粒流数值模拟方法,其特征在于步骤1中,土石分界粒径的理论公式为d0=0.05Lc,d0为土石分界粒径,Lc为采样窗口特征尺寸。
5.根据权利要求1所述的反映真实工程土石混合体分布的颗粒流数值模拟方法,其特征在于步骤3中,石块等效尺寸是指二值图像中石块的等面积圆的直径。
6.根据权利要求1所述的反映真实工程土石混合体分布的颗粒流数值模拟方法,其特征在于步骤5中,颗粒接触本构模型采用可传递力和力矩的平行粘接模型。
7.根据权利要求1所述的反映真实工程土石混合体分布的颗粒流数值模拟方法,其特征在于步骤5中,需要标定的细观参数包括颗粒密度、颗粒粒径、接触刚度、粘结强度及摩擦系数。
CN201811321036.4A 2018-11-07 2018-11-07 一种反映真实工程土石混合体分布的颗粒流数值模拟方法 Active CN109359416B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811321036.4A CN109359416B (zh) 2018-11-07 2018-11-07 一种反映真实工程土石混合体分布的颗粒流数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811321036.4A CN109359416B (zh) 2018-11-07 2018-11-07 一种反映真实工程土石混合体分布的颗粒流数值模拟方法

Publications (2)

Publication Number Publication Date
CN109359416A true CN109359416A (zh) 2019-02-19
CN109359416B CN109359416B (zh) 2023-10-24

Family

ID=65344383

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811321036.4A Active CN109359416B (zh) 2018-11-07 2018-11-07 一种反映真实工程土石混合体分布的颗粒流数值模拟方法

Country Status (1)

Country Link
CN (1) CN109359416B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109977357A (zh) * 2019-03-26 2019-07-05 长安大学 一种土石混合体填料路基压实质量和强度的控制方法
CN110516366A (zh) * 2019-08-28 2019-11-29 北京工业大学 一种基于超高效液相色谱分析中随机微珠微混合器的建模方法
CN110990911A (zh) * 2019-10-29 2020-04-10 武汉大学 一种基于颗粒簇替代的岩石晶体结构快速建模方法
CN111504862A (zh) * 2020-04-29 2020-08-07 中铁二院工程集团有限责任公司 一种基于图像技术的粗颗粒土级配确定方法
CN111666620A (zh) * 2020-05-29 2020-09-15 黄河水利委员会黄河水利科学研究院 一种颗粒材料侧向压力分布特征定量描述方法
CN111737860A (zh) * 2020-06-03 2020-10-02 山东大学 岩土体渗透破坏过程强度演化特征的等效模拟方法及系统
CN113515878A (zh) * 2021-07-07 2021-10-19 重庆交通大学 一种基于块石形状和破碎的堆积体离散元三维建模方法
CN113533146A (zh) * 2021-07-09 2021-10-22 清华大学 基于图像识别技术的堆石模拟分析方法及系统
CN114036867A (zh) * 2021-11-08 2022-02-11 四川大学 一种滑坡堆积碎石材料离散元细观参数标定方法及系统
CN116738519A (zh) * 2023-08-14 2023-09-12 吉林大学 一种小麦植株通用建模方法
CN117744412A (zh) * 2024-02-19 2024-03-22 西安交通大学 一种岩石非线性力学模拟方法、系统、设备及介质
CN117807814A (zh) * 2024-02-29 2024-04-02 天津理工大学 一种模拟现场碾压效果的路基碾压过程建模方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103605841A (zh) * 2013-11-07 2014-02-26 河海大学 一种水工沥青混凝土数值仿真模型构建方法
CN105654482A (zh) * 2015-12-29 2016-06-08 河海大学 基于数字图像技术的土石混合体表征单元体尺度的确定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103605841A (zh) * 2013-11-07 2014-02-26 河海大学 一种水工沥青混凝土数值仿真模型构建方法
CN105654482A (zh) * 2015-12-29 2016-06-08 河海大学 基于数字图像技术的土石混合体表征单元体尺度的确定方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
DUOFENG CEN: "Shear deformation and strength of the interphase between the soil–rock mixture and the benched bedrock slope surface", 《ACTA GEOTECHNICA》 *
丁秀丽等: "基于数字图像的土石混合体力学性质的颗粒流模拟", 《岩石力学与工程学报》 *
徐文杰等: "土石混合体细观结构及力学特性数值模拟研究", 《岩石力学与工程学报》 *
王鹏: "土石混合体大型三轴试验及其细观力学特定的PFC模拟", 《中国优秀硕士学位论文全文数据库》 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109977357B (zh) * 2019-03-26 2022-10-04 长安大学 一种土石混合体填料路基压实质量和强度的控制方法
CN109977357A (zh) * 2019-03-26 2019-07-05 长安大学 一种土石混合体填料路基压实质量和强度的控制方法
CN110516366A (zh) * 2019-08-28 2019-11-29 北京工业大学 一种基于超高效液相色谱分析中随机微珠微混合器的建模方法
CN110516366B (zh) * 2019-08-28 2023-04-07 北京工业大学 一种基于超高效液相色谱分析中随机微珠微混合器的建模方法
CN110990911A (zh) * 2019-10-29 2020-04-10 武汉大学 一种基于颗粒簇替代的岩石晶体结构快速建模方法
CN111504862A (zh) * 2020-04-29 2020-08-07 中铁二院工程集团有限责任公司 一种基于图像技术的粗颗粒土级配确定方法
CN111666620A (zh) * 2020-05-29 2020-09-15 黄河水利委员会黄河水利科学研究院 一种颗粒材料侧向压力分布特征定量描述方法
CN111666620B (zh) * 2020-05-29 2023-06-09 黄河水利委员会黄河水利科学研究院 一种颗粒材料侧向压力分布特征定量描述方法
CN111737860A (zh) * 2020-06-03 2020-10-02 山东大学 岩土体渗透破坏过程强度演化特征的等效模拟方法及系统
CN113515878A (zh) * 2021-07-07 2021-10-19 重庆交通大学 一种基于块石形状和破碎的堆积体离散元三维建模方法
CN113533146B (zh) * 2021-07-09 2022-07-08 清华大学 基于图像识别技术的堆石模拟分析方法及系统
CN113533146A (zh) * 2021-07-09 2021-10-22 清华大学 基于图像识别技术的堆石模拟分析方法及系统
CN114036867A (zh) * 2021-11-08 2022-02-11 四川大学 一种滑坡堆积碎石材料离散元细观参数标定方法及系统
CN114036867B (zh) * 2021-11-08 2024-04-12 四川大学 一种滑坡堆积碎石材料离散元细观参数标定方法及系统
CN116738519A (zh) * 2023-08-14 2023-09-12 吉林大学 一种小麦植株通用建模方法
CN116738519B (zh) * 2023-08-14 2023-12-01 吉林大学 一种小麦植株通用建模方法
CN117744412A (zh) * 2024-02-19 2024-03-22 西安交通大学 一种岩石非线性力学模拟方法、系统、设备及介质
CN117744412B (zh) * 2024-02-19 2024-05-10 西安交通大学 一种岩石非线性力学模拟方法、系统、设备及介质
CN117807814A (zh) * 2024-02-29 2024-04-02 天津理工大学 一种模拟现场碾压效果的路基碾压过程建模方法
CN117807814B (zh) * 2024-02-29 2024-05-10 天津大学 一种模拟现场碾压效果的路基碾压过程建模方法

Also Published As

Publication number Publication date
CN109359416B (zh) 2023-10-24

Similar Documents

Publication Publication Date Title
CN109359416A (zh) 一种反映真实工程土石混合体分布的颗粒流数值模拟方法
Liu et al. Mechanical properties of brittle rock governed by micro-geometric heterogeneity
Zhou et al. Novel grain-based model for simulation of brittle failure of Alxa porphyritic granite
Müller et al. Modeling deformation and damage of rock salt using the discrete element method
Xu et al. Mesoscopic damage and fracturing of heterogeneous brittle rocks based on three-dimensional polycrystalline discrete element method
Lv et al. Effects of internal pores on the mechanical properties of marine calcareous sand particles
Jiang et al. Contact behavior of idealized granules bonded in two different interparticle distances: An experimental investigation
Wang et al. Evaluation of the anisotropy and directionality of a jointed rock mass under numerical direct shear tests
Lu et al. Uniaxial strength and failure in sandstone containing a pre-existing 3-D surface flaw
CN105319603A (zh) 致密砂岩储层复杂网状裂缝的预测方法
Yu et al. Numerical modeling of jointed rock under compressive loading using X-ray computerized tomography
Kuang et al. Experimental and numerical study on the fragmentation mechanism of a single calcareous sand particle under normal compression
Yang et al. Correlation study between fracability and brittleness of shale-gas reservoir
Beck Conceptual approaches for the analysis of coupled hydraulic and geomechanical processes
Sun et al. Comparison of crack processes in single-flawed rock-like material using two bonded–particle models under compression
Zhou et al. Mechanical and fracture behaviors of brittle material with a circular inclusion: insight from infilling composition
Li et al. A parallel-bonded chemical corrosion model for discrete element modelling of chemically corroded limestone
Wang et al. A model of anisotropic property of seepage and stress for jointed rock mass
Liu et al. Mechanical characteristics of soil-rock mixtures containing macropore structure based on 3D modeling technology
Song et al. Mechanical behavior and fracture evolution mechanism of composite rock under triaxial compression: insights from three-dimensional DEM modeling
CN114626215A (zh) 一种基于离散元与合成岩体技术的岩质边坡建模方法
Tan et al. Simulation of rock-breaking process by drilling machine and dynamic classification of surrounding rocks
Zhang et al. Numerical study of the evolution of cohesion and internal friction in rock during the pre-peak deformation process
Zengin et al. Experimental investigation of pore-fracture relationship on failure behaviour of porous rock materials
Zhang Failure mechanism and stability analysis of rock slope

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