CN116226982A - 一种粘性土-岩石隧道开挖耦合数值方法 - Google Patents
一种粘性土-岩石隧道开挖耦合数值方法 Download PDFInfo
- Publication number
- CN116226982A CN116226982A CN202310098900.3A CN202310098900A CN116226982A CN 116226982 A CN116226982 A CN 116226982A CN 202310098900 A CN202310098900 A CN 202310098900A CN 116226982 A CN116226982 A CN 116226982A
- Authority
- CN
- China
- Prior art keywords
- rock
- cohesive soil
- triangular
- elastic modulus
- triangle
- 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 106
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000009412 basement excavation Methods 0.000 title claims abstract description 26
- 239000000463 material Substances 0.000 claims abstract description 33
- 239000002131 composite material Substances 0.000 claims abstract description 25
- 238000009826 distribution Methods 0.000 claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 12
- 230000008878 coupling Effects 0.000 claims abstract description 9
- 238000010168 coupling process Methods 0.000 claims abstract description 9
- 238000005859 coupling reaction Methods 0.000 claims abstract description 9
- 239000002689 soil Substances 0.000 claims description 55
- 230000006870 function Effects 0.000 claims description 18
- 230000008859 change Effects 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000005315 distribution function Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 5
- 238000009827 uniform distribution Methods 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000012669 compression test Methods 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 239000004927 clay Substances 0.000 abstract 4
- 125000004122 cyclic group Chemical group 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000011437 continuous method Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011439 discrete element method Methods 0.000 description 1
- 238000001125 extrusion Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种粘性土‑岩石复合地层隧道开挖有限元‑离散元耦合数值模拟方法，步骤如下：(1)建立均质围岩隧道开挖数值模型，获得任意三角形单元弹性模量；(2)对所有三角形单元随机性赋值，得到所有三角形单元的随机弹性模量；(3)计算三角形单元的弹性模量E和弹性模量期望值E0的比值，若E/E0≤x，则三角形单元设为粘性土材料，若E/E0>x，则三角形单元设为岩石材料；(4)设立二维数组表示三角形单元的三个相邻的三角形单元，对所有三角形进行循环判定以优化空间分布；(5)将粘性土材料和岩石材料分别赋予各自力学属性。本模拟方法能够建立岩石、粘性土在空间上具有随机分布特征的数值模型，且岩石、粘性土的形状和体积比亦具有随机性。
Description
技术领域
本发明属于岩土力学与岩土工程领域，具体涉及粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模拟方法。
背景技术
在城市浅表地铁隧道建设过程中，常常遭遇粘性土-岩石复合地层，呈现粘性土中包裹漂石的表观特性。岩石、粘性土在空间上具有随机分布特征，且岩石、粘性土的形状和体积亦具有随机性。粘性土-岩石复合地层中粘性土、岩石的力学特性和空间分布特征对于地铁隧道变形和地表沉降具有重要影响。
为研究粘性土-岩石复合地层隧道开挖后围岩变形及破坏响应特征，为隧道掘进和支护参数设计提供指导，常采用理论解析、室内试验、现场监测和数值模拟方法。隧道开挖后，围岩不仅会发生弹塑性连续变形，亦将发生破裂，尤其对于粘性土-岩石复合地层，粘性土自身及与岩石的交界面将发生破裂，甚至岩石也将发生破裂，这取决于岩石及粘性土的物理力学特征，破裂后的围岩将产生块体间的接触挤压和块体自身的滑移翻转大运动，理论解析解难以描述围岩裂隙网络扩展、碎裂块体接触和块体的宏观大运动等复杂过程。
对于室内试验及现场监测而言，试验结果和现场观测结果往往仅对本次试验或本观测洞段有效，且室内试验难以建立具有随机分布特征的粘性土-岩石复合地层物理模型；而现场观测存在不同洞段地质条件差异大、不同观测结果对比性差、且围岩内部信息难以获知等问题，此外，这两种方法均存在试验/观测周期长、成本高昂等缺点，难以适用于粘性土-岩石复合地层隧道开挖后的变形响应研究中。
因此，数值模拟方法成为一种较佳研究手段。一般地，可分为连续性方法、非连续方法及耦合方法。连续性方法如有限元法和有限差分法难以模拟隧道围岩裂隙网络扩展过程，而非连续性方法如离散元法和颗粒元法等难以模拟岩石块体变形、块体内部开裂等过程。因此，有限元-离散元耦合数值模拟方法FDEM应运而生，并在岩土力学及岩土工程研究中得到了广泛应用。然而，当前FDEM数值模拟研究大都针对均质完整性围岩或具有原生裂隙的围岩，尚未见到对粘性土-岩石复合地层的研究，对隧道开挖后粘性土-岩石复合地层变形及破裂机理认识不清，究其原因在于缺乏粘性土-复合地层的有效建模方法。
发明内容
为解决上述技术问题，本发明的目的是提出一种粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模拟方法，该方法能够得到粘性土-岩石复合地层隧道开挖数值模型，且粘性土、岩石具有空间分布随机性特征，此外，所建数值模型中的粘性土、岩石几何形状、体积亦具有随机性特征，符合真实隧道粘性土-岩石复合地层的分布特征，为揭示隧道开挖粘性土-岩石复合地层变形及破坏机理研究提供数值模拟基础。
本发明所采用的技术方案如下：
一种粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模拟方法，包括如下步骤：
步骤1.建立均质围岩隧道开挖数值模型，根据累积分布函数计算获得任意三角形单元的弹性模量；
步骤2.对所有三角形单元随机性赋值，得到数值模型中所有三角形单元的随机弹性模量；
步骤3.计算三角形单元的弹性模量E和弹性模量期望值E0的比值，若E/E0≤x，则三角形单元设为粘性土材料，若E/E0>x，则三角形单元设为岩石材料；其中x为比值分界线；
步骤4.设立二维数组分别表示三角形单元的三个相邻的三角形单元，对所有三角形进行循环判定，若粘性土三角形相邻的两个三角形单元或相邻的三个三角形单元为岩石，则该粘性土三角形单元设为岩石材料防止出现畸形岩石或孤立的粘性土，以实现粘性土、岩石空间分布及几何形状的优化；
步骤5.将粘性土材料和岩石材料分别赋予各自力学属性，粘性土-岩石界面的四边形节理单元参数设为粘性土材料的参数，得到粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模型。
进一步，所述步骤1中，建模软件为Gmsh。
进一步，所述步骤1中，根据累积分布函数得到任意三角形单元的弹性模量
E0为弹性模量期望值、m为非均质度系数，均由输入文件给定，F(E)为0～1间的均匀分布函数。
更进一步，所述步骤1中，F(E)的值由c语言rand()函数产生0～32767之间的随机数并除以32768得到。
更进一步，所述步骤1中，m值亦可任意调整以改变粘性土-岩石的空间分布特征。
优选的，所述步骤1中，E0和m均由输入文件给定，m可设为1.5、E0可设为1GPa。
进一步，所述步骤2中，任意三角形单元弹性模量赋值方法如下：①定义随机数参数suijishu，且suijishu＝rand()/32768.0；②三角形弹性模量E＝E0×pow{log[1.0/(1-suijishu)],1.0/m}，pow为c语言幂函数，log为c语言以自然指数e为底的对数函数。
进一步，所述步骤3中，比值分界线x可任意调整，以改变粘性土、岩石的体积占比。
进一步，所述粘性土和岩石分类方法为：将粘性土和岩石三角形属性定义为不同数值，利用三角形的E/E0值与x的大小关系，对三角形属性进行赋值。优选的，x取1。
进一步，所述步骤4中，设立二维数组i2next[0][ielem]、i2next[1][ielem]、i2next[2][ielem]分别表示三角形单元ielem的三个相邻三角形单元，若粘性土三角形相邻的两个三角形单元或相邻的三个三角形单元为岩石，则该粘性土三角形单元设为岩石材料，防止出现畸形岩石或孤立的粘性土。
进一步，所述步骤5中，粘性土材料和岩石材料的力学参数均分别通过各自的单轴压缩、巴西劈裂和三轴压缩试验标定得到。
本发明的有益效果：
本发明提出的粘性土-岩石复合地层隧道开挖耦合数值方法，能够建立具有空间随机分布特征的粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模型，且粘性土、岩石的几何形状、体积亦具有随机分布特征。此外，还能够任意改变粘性土及岩石的体积占比及空间分布特征。
附图说明
图1为粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模型建模流程。
图2为非均质度系数m＝1.5时的概率密度函数。
图3为粘性土-岩石空间分布优化前后对比图。
图4为粘性土-岩石交界面四边形节理单元参数的分配方式。
图5为建立的粘性土-岩石复合地层隧道开挖数值模型。
具体实施方式
下面结合附图对本发明涉及的粘性土-岩石复合地层隧道开挖耦合数值方法进行详细说明。本发明的内容完全不限于此。
实施例
一种粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模拟方法，如图1所示，步骤如下：
步骤1.采用Gmsh软件建立均质围岩隧道开挖数值模型，采用公式
获得任意三角形单元的弹性模量，E0为弹模期望值、m为非均质度系数，均由输入文件给定，F(E)为0～1间的均匀分布函数，由c语言rand()函数产生0～32767之间的随机数并除以32768得到，32767是由c语言rand()函数所限定的。
公式
是借助于非均质岩样建模方法得到的，对于非均质岩样，岩石弹性模量服从Weibull分布，其概率密度函数如下：
m＝1.5时的概率密度函数f(E)如图2所示。概率密度函数对应的累积分布函数如下：
根据累积分布函数可得到任意三角形单元的弹性模量E为
对于公式
m值可设为1.5、E0可设为1GPa。F(E)为0-1之间的均匀分布函数，即0-1之间的任意值的取值概率是相同的。
任意三角形单元弹性模量赋值方法为：①定义随机数参数suijishu，且suijishu＝rand()/32768.0；②三角形弹性模量E＝E0×pow{log[1.0/(1-suijishu)],1.0/m}；本实施例中，E0＝1GPa、m＝1.5，pow为c语言幂函数，log为c语言以自然指数e为底的对数函数。
步骤2.对所有三角形单元循环，得到数值模型中所有三角形单元的弹性模量。
对于建立的均质围岩隧道开挖数值模型，采用公式
对所有三角形单元弹性模量进行随机性赋值，得到岩石弹性模量服从Weibull分布的非均质围岩隧道开挖数值模型。
步骤3.将E/E0≤1的三角形单元设为粘性土材料，将E/E0>1的三角形单元设为岩石材料。
在对三角形单元弹性模量赋值完毕后，将其与期望值E0对比，若比值≤1，则将该三角形单元属性设为粘性土材料，并重新赋为粘性土材料的弹性模量；若比值>1，则该三角形单元属性设为岩石材料，并赋为岩石材料的弹性模量，如图2所示。
粘性土和岩石分类方法为：①定义粘性土三角形属性为1、岩石三角形属性为2；②对所有三角形单元采用如下判别语句：
if(E/E0≤1),i1elpr[ielem]＝1；
else i1elpr[ielem]＝2；
其中，if、else为c语言判别语句，i1elpr为三角形属性数组，ielem为三角形单元编号。
步骤4.设立二维数组i2next[0][ielem]、i2next[1][ielem]、i2next[2][ielem]分别表示三角形单元ielem的三个相邻三角形单元，对所有的三角形单元进行循环，若粘性土三角形相邻的两个三角形单元或相邻的三个三角形单元为岩石，则该粘性土三角形单元设为岩石材料，防止出现畸形岩石或孤立的粘性土，实现粘性土、岩石空间分布及几何形状的优化，如图3所示。
优化实现方法为：①定义三角形单元ielem的三个相邻三角形单元分别为i2next[0][ielem]、i2next[1][ielem]、i2next[2][ielem]；②采用如下判别语句：
if(i1elpr[ielem]＝＝1&&
((i1elpr[i2next[0][ielem]]＝＝2&&i1elpr[i2next[1][ielem]]＝＝2)||
(i1elpr[i2next[0][ielem]]＝＝2&&i1elpr[i2next[2][ielem]]＝＝2)||
(i1elpr[i2next[1][ielem]]＝＝2&&i1elpr[i2next[2][ielem]]＝＝2)||
(i1elpr[i2next[0][ielem]]＝＝2&&i1elpr[i2next[1][ielem]]＝＝2&&i1elpr[i2next[2][ielem]]＝＝2)))
i1elpr[ielem]＝2；
其中，“&&”为c语言“且”符号，“||”为c语言“或”符号，“＝＝”为c语言判别公式左侧、右侧是否相等的符号。
步骤5.将粘性土材料和岩石材料分别赋予各自力学属性，粘性土-岩石界面的四边形节理单元参数设为粘性土材料的参数，如图4所示；粘性土材料和岩石材料的力学参数均分别通过各自的单轴压缩、巴西劈裂和三轴压缩试验标定得到；粘性土与粘性土三角形界面处的四边形节理单元参数设为粘性土的参数、岩石与岩石三角形界面处的四边形节理单元参数设为岩石的参数。
粘性土-岩石界面处四边形节理单元属性分配方法为：①定义i1linbian[0][jelem]、i1linbian[1][jelem]分别为四边形单元jelem两侧的三角形单元；②采用如下判别语句：
if((i1elpr[i1linbian[0][jelem]]＝＝1&&i1elpr[i1linbian[1][jelem]]＝＝2)||
(i1elpr[i1linbian[0][jelem]]＝＝2&&i1elpr[i1linbian[1][jelem]]＝＝1))
i1elpr[jelem]＝粘性土；
根据上述步骤，得到了如图5所示的粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模型。根据本发明提出的方法，能够建立具有空间随机分布特征的粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模型，且粘性土、岩石的几何形状、体积亦具有随机分布特征，此外，还能够通过调整E/E0的比值分界线改变粘性土及岩石的体积占比、通过调整m值改变粘性土及岩石的空间分布特征。
本发明所涉及的粘性土-岩石复合地层隧道开挖耦合数值方法不仅仅局限于以上实施例中所描述的内容，而是以权利要求所限定的范围为准。本发明所属领域技术人员在该实施例的基础上所做的任何修改或补充或等效替换，都在本发明的权利要求所要求保护的范围内，如对三角形单元尺寸、隧道模型尺寸、二维数组名称、输入值m和E0及E/E0比值分界线的修改等。
以上所述，仅为本发明较佳的具体实施方式，但本发明保护的范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内所做的任何修改，等同替换和改进等，均应包含在发明的保护范围之内。
Claims (8)
1.一种粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模拟方法，其特征在于，包括如下步骤：
步骤1.建立均质围岩隧道开挖数值模型，根据累积分布函数计算获得任意三角形单元的弹性模量；
步骤2.对所有三角形单元随机性赋值，得到数值模型中所有三角形单元的随机弹性模量；
步骤3.计算三角形单元的弹性模量E和弹性模量期望值E0的比值，若E/E0≤x，则三角形单元设为粘性土材料，若E/E0>x，则三角形单元设为岩石材料；其中x为比值分界线；
步骤4.设立二维数组分别表示三角形单元的三个相邻的三角形单元，对所有三角形进行循环判定，若粘性土三角形相邻的两个三角形单元或相邻的三个三角形单元为岩石，则该粘性土三角形单元设为岩石材料，以实现粘性土、岩石空间分布及几何形状的优化；
步骤5.将粘性土材料和岩石材料分别赋予各自力学属性，粘性土-岩石界面的四边形节理单元参数设为粘性土材料的参数，得到粘性土-岩石复合地层隧道开挖有限元-离散元耦合数值模型。
2.根据权利要求1所述的方法，其特征在于：所述步骤1中，根据累积分布函数得到任意三角形单元的弹性模量
E0为弹模期望值、m为非均质度系数，均由输入文件给定，F(E)为0～1间的均匀分布函数。
3.根据权利要求2所述的方法，其特征在于：所述步骤1中，F(E)的值由c语言rand()函数产生0～32767之间的随机数并除以32768得到。
4.根据权利要求2所述的方法，其特征在于：所述步骤1中，m值亦可任意调整以改变粘性土-岩石的空间分布特征。
5.根据权利要求1所述的方法，其特征在于：所述步骤2中，任意三角形单元弹性模量赋值方法如下：①定义随机数参数suijishu，且suijishu＝rand()/32768.0；②三角形弹性模量E＝E0×pow{log[1.0/(1-suijishu)],1.0/m}，pow为c语言幂函数，log为c语言以自然指数e为底的对数函数。
6.根据权利要求1所述的方法，其特征在于：所述步骤3中，比值分界线x可任意调整，以改变粘性土、岩石的体积占比及空间分布。
7.根据权利要求1所述的方法，其特征在于，所述粘性土和岩石分类方法为：将粘性土和岩石三角形属性定义为不同数值，利用三角形的E/E0值与x的大小关系，对三角形属性进行赋值。
8.根据权利要求1所述的方法，其特征在于：所述步骤5中，粘性土材料和岩石材料的力学参数均分别通过各自的单轴压缩、巴西劈裂和三轴压缩试验标定得到。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310098900.3A CN116226982B (zh) | 2023-01-31 | 2023-01-31 | 一种粘性土-岩石隧道开挖耦合数值方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310098900.3A CN116226982B (zh) | 2023-01-31 | 2023-01-31 | 一种粘性土-岩石隧道开挖耦合数值方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116226982A true CN116226982A (zh) | 2023-06-06 |
| CN116226982B CN116226982B (zh) | 2024-05-28 |
Family
ID=86578028
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310098900.3A Active CN116226982B (zh) | 2023-01-31 | 2023-01-31 | 一种粘性土-岩石隧道开挖耦合数值方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116226982B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117094057A (zh) * | 2023-08-22 | 2023-11-21 | 中铁建大桥工程局集团第五工程有限公司 | 一种隧道围岩松动圈测试方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110077918A1 (en) * | 2009-09-25 | 2011-03-31 | Ovunc Mutlu | Method of Predicting Natural Fractures And Damage In A Subsurface Region |
| CN110309536A (zh) * | 2019-05-24 | 2019-10-08 | 中南大学 | 一种岩土三轴试验柔性薄膜边界的离散元模拟方法 |
| CN112329290A (zh) * | 2020-10-28 | 2021-02-05 | 河海大学 | 可用于施工过程模拟的有限元离散元多尺度耦合计算方法 |
-
2023
- 2023-01-31 CN CN202310098900.3A patent/CN116226982B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110077918A1 (en) * | 2009-09-25 | 2011-03-31 | Ovunc Mutlu | Method of Predicting Natural Fractures And Damage In A Subsurface Region |
| CN110309536A (zh) * | 2019-05-24 | 2019-10-08 | 中南大学 | 一种岩土三轴试验柔性薄膜边界的离散元模拟方法 |
| CN112329290A (zh) * | 2020-10-28 | 2021-02-05 | 河海大学 | 可用于施工过程模拟的有限元离散元多尺度耦合计算方法 |
Non-Patent Citations (2)
| Title |
|---|
| PENGHAI DENG等: "FDEM numerical study on the mechanical characteristics and failure behavior of heterogeneous rock based on the Weibull distribution of mechanical parameters", COMPUTERS AND GEOTECHNICS, vol. 154, 15 November 2022 (2022-11-15), pages 105138 * |
| 严成增等: "基于数字图像技术的岩土材料有限元-离散元分析", 岩土力学, vol. 35, no. 8, 31 August 2014 (2014-08-31), pages 2408 - 2414 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117094057A (zh) * | 2023-08-22 | 2023-11-21 | 中铁建大桥工程局集团第五工程有限公司 | 一种隧道围岩松动圈测试方法 |
| CN117094057B (zh) * | 2023-08-22 | 2024-11-19 | 中铁建大桥工程局集团第五工程有限公司 | 一种隧道围岩松动圈测试方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116226982B (zh) | 2024-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110779795B (zh) | 裂缝性储层地质力学建模网格单元大小确定方法 | |
| CN113534291B (zh) | 岩石力学层约束下的低渗透储层不同尺度裂缝定量预测方法 | |
| CN101739716B (zh) | 工程岩体三维空间结构建模与关键块识别方法 | |
| Ghazvinian et al. | 3D random Voronoi grain-based models for simulation of brittle rock damage and fabric-guided micro-fracturing | |
| Chen et al. | Numerical study on spalling failure of rock surrounding deep buried tunnel based on DEM | |
| CN109271738B (zh) | 一种用于获取巷道围岩Weibull分布参数的数值反演方法 | |
| Liu et al. | Stability assessment of the Three-Gorges Dam foundation, China, using physical and numerical modeling—Part I: physical model tests | |
| Zhang et al. | Computational homogenization for mechanical properties of sand cobble stratum based on fractal theory | |
| CN104458309A (zh) | 一种用于物理模拟实验中的相似材料配比确定方法 | |
| CN110702881B (zh) | 岩土材料参数变异性结果的预测方法及其应用 | |
| CN115203900B (zh) | 基于离散元的裂隙岩体建模方法 | |
| Zhao et al. | Laboratory test and numerical simulations for 3D printed irregular columnar jointed rock masses under biaxial compression | |
| CN115324560A (zh) | 利用地应力场模拟确定压裂诱发油气套管变形位置的方法 | |
| Wang et al. | Development of an improved three-dimensional rough discrete fracture network model: method and application | |
| Chen et al. | Spatial estimation of material parameters and refined finite-element analysis of rockfill dam based on construction digitization | |
| Gao et al. | Proposal of a chemical weathering model for mechanical analysis of cemented sandstone using DDA | |
| CN113378410A (zh) | 一种采动覆岩导水通道演化的模拟方法 | |
| Zhang et al. | Failure modes of weak interlayers with different dip angles in red mudstone strata, Northwest China | |
| CN116226982A (zh) | 一种粘性土-岩石隧道开挖耦合数值方法 | |
| CN115324556A (zh) | 一种压裂诱发油气套管变形风险级别的综合预测方法 | |
| CN115221727B (zh) | 一种基于含水率的岩体的数值仿真模型参数确定方法 | |
| He et al. | Research on the multi-scale DFN-DEM equivalent modelling method for jointed rock masses and the collapse law of block structures | |
| Ma et al. | Optimizing 3D granular modeling with integrated 3DEC and neper techniques for granite mechanics simulation | |
| CN110705168A (zh) | 构造应力场的模拟方法 | |
| CN117390917A (zh) | 一种考虑采空区覆岩垮落的工作面采动应力数值模拟方法 |
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 |