CN113361099B - 一种裂隙岩体模拟方法及系统 - Google Patents
一种裂隙岩体模拟方法及系统 Download PDFInfo
- Publication number
- CN113361099B CN113361099B CN202110623951.4A CN202110623951A CN113361099B CN 113361099 B CN113361099 B CN 113361099B CN 202110623951 A CN202110623951 A CN 202110623951A CN 113361099 B CN113361099 B CN 113361099B
- Authority
- CN
- China
- Prior art keywords
- rock mass
- fractured rock
- deformation
- mechanical model
- model
- 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
Links
Images
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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种裂隙岩体模拟方法及系统,通过对裂隙岩体几何、物理参数的不连续变形相似准则的研究,提供了确定裂隙岩体几何、物理参数不连续变形相似准则的确定方法,由于该相似准则全面反映了材料、几何、边界条件与动力条件等方面的相似性,能够更好的模拟工程实际中大尺度复杂裂隙网络岩体的结构形式与受载情况,为物理模型试验方法提供指导,实现对裂隙岩体更准确的模拟。
Description
技术领域
本发明涉及岩体力学技术领域,特别是涉及一种裂隙岩体模拟方法及系统。
背景技术
岩体由于经受过复杂的地质作用,导致其内部形成各种各样的结构面,如断层、节理、裂隙等。裂隙岩体变形性能是岩石力学理论研究的重要基础,也是评价边坡、地下洞室和隧道安全稳定的关键因素。对于工程实际中的大尺度复杂裂隙岩体结构,直接采用理论分析和数值计算的方法研究其在复杂荷载作用下的变形特征具有很大的困难,而进行现场试验耗费的成本太高,且重复性差。物理模型试验方法具有可控制性、可重复性等优点,且易于实现。
物理模型试验的主要依据是相似准则,包括材料相似、几何相似、边界条件相似和动力相似等。由于工程实际中的大尺度复杂裂隙网络岩体的结构形式和受载情况都较为复杂,现有的物理模型试验方法很难同时满足上述所有的相似条件,不能完整的反映出裂隙岩体的变形规律。
发明内容
本发明的目的是提供一种裂隙岩体模拟方法及系统,综合考虑了裂隙岩体几何、物理参数的不连续变形相似准则的建立,通过协调各参数之间的关系,实现了对裂隙岩体变形情况的准确模拟。
为实现上述目的,本发明提供了如下方案:
一种裂隙岩体模拟方法,所述方法包括:
构建裂隙岩体结构力学模型;
根据所述裂隙岩体结构力学模型计算反映裂隙岩体的不连续变形模量方程;
根据所述裂隙岩体结构力学模型建立不连续变形相似力学模型;
根据所述不连续变形模量方程和所述不连续变形相似力学模型,利用相似变换得到反映裂隙岩体变形特性的相似准则;
根据所述相似准则采用物理模型试验方法对裂隙岩体的变形情况进行模拟。
一种裂隙岩体模拟系统,所述系统包括:
原始模型模块,用于构建裂隙岩体结构力学模型;
不连续变形模量模块,用于根据所述裂隙岩体结构力学模型计算反映裂隙岩体的不连续变形模量方程;
相似模型模块,用于根据所述裂隙岩体结构力学模型建立不连续变形相似力学模型;
相似准则模块,用于根据所述不连续变形模量方程和所述不连续变形相似力学模型,利用相似变换得到反映裂隙岩体变形特性的相似准则;
模拟模块,用于根据所述相似准则采用物理模型试验方法对裂隙岩体的变形情况进行模拟。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种裂隙岩体模拟方法及系统,通过对裂隙岩体几何、物理参数的不连续变形相似准则的研究,提供了确定裂隙岩体几何、物理参数不连续变形相似准则的确定方法,由于该相似准则全面反映了材料、几何、边界条件与动力条件等方面的相似性,能够更好的模拟工程实际中大尺度复杂裂隙网络岩体的结构形式与受载情况,为物理模型试验方法提供指导,实现对裂隙岩体更准确的模拟。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种裂隙岩体模拟方法流程图;
图2为本发明实施例提供的一种二维贯通裂隙岩体结构力学模型示意图;
图3为本发明实施例提供的一种二维复杂裂隙岩体模型轴向加载条件下变形分析及相似机理分析过程示意图;
图4为本发明实施例提供的二维裂隙岩体不连续变形模量计算原始模型示意图;
图5(a)为10倍纯几何扩大模型与10倍几何扩大不连续变形相似力学模型示意图,(b)为不同岩石基质模量变化对应的10倍几何扩大不连续变形相似力学模型示意图;
图6为本发明实施例提供的一种裂隙岩体模拟系统框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
现有的对裂隙岩体变形模拟的研究主要集中在针对物理模型试验相似材料的研究以及离心模型试验的研究。然而单独考虑部分变量的相似性进行物理模型试验研究具有很大的局限性,不能完整的反映出裂隙岩体的变形规律。
工程实际中岩体含有大量的裂隙网络,单裂隙相似材料岩样的研究能为工程裂隙岩体的变形规律提供一定的理论和试验基础,但仍未能反映原始裂隙岩体模型的一般性力学特性。
本发明的目的是提供一种综合考虑裂隙岩体几何、物理参数的裂隙岩体模拟方法及系统,更好的模拟工程实际中大尺度复杂裂隙网络岩体的结构形式与受载情况,为物理模型试验方法提供指导,实现对裂隙岩体更准确的模拟。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
如图1所示,本实施例提供了一种裂隙岩体模拟方法,所述方法包括:
步骤101:构建裂隙岩体结构力学模型;
步骤102:根据所述裂隙岩体结构力学模型计算反映裂隙岩体的不连续变形模量方程;
步骤103:根据所述裂隙岩体结构力学模型建立不连续变形相似力学模型;
步骤104:根据所述不连续变形模量方程和所述不连续变形相似力学模型,利用相似变换得到反映裂隙岩体变形特性的相似准则;
步骤105:根据所述相似准则采用物理模型试验方法对裂隙岩体进行模拟。
为了简化计算步骤,在构建裂隙岩体力学模型时,考虑到二维平面应力问题,忽略围压对加载方向位移的影响,根据轴向加载条件得到二维贯通裂隙岩体结构力学模型如图2所示。图2中,岩块尺寸为L1×L2,岩石基质杨氏模量和泊松比分别为E0和μ0,承受轴向加载应力为σ1,作用在结构面上的平均应力σc=σ1cosα,结构面上的法向以及切向应力分量分别为σcn=σccosα和τcs=σcsinα。采用如图3所示二维复杂裂隙岩体模型轴向加载条件下变形分析及相似机理分析方法,建立相应二维复杂裂隙网络岩体模型作为裂隙岩体结构力学模型,其中第i条裂隙长度和倾角分别为li和αi,裂隙结构面法向刚度和切向刚度分别为和
在根据所述二维贯通裂隙岩体结构力学模型建立相应的二维复杂裂隙网络岩体模型时,首先考虑二维贯通裂隙岩体结构力学模型中裂隙端部约束影响,引入第i条裂隙端部对该裂隙变形的影响系数ri,然后考虑裂隙连通率影响,引入裂隙连通率对岩体裂隙变形的影响系数f,最终得到裂隙岩体结构力学模型。
得到裂隙岩体结构力学模型后,根据裂隙岩体结构力学模型计算反映裂隙岩体的不连续变形模量方程。先根据所述裂隙岩体结构力学模型建立轴向加载位移方程,然后根据所述轴向加载位移方程得到反映裂隙岩体的不连续变形模量方程。
具体的,先分析单贯通裂隙岩体沿轴向加载方向位移方程为:
这里u0=σ1L1/E0为岩石基质沿加载方向的位移,考虑裂隙尖端约束以及裂隙间相互作用对裂隙结构面变形的影响,图3中裂隙岩体结构力学模型中轴向加载方向的位移可以表示为:
其中,u0为岩石基质沿加载方向的位移,σ1为承受轴向加载应力,li和αi分别为第i条裂隙长度和倾角,和分别为裂隙结构面法向刚度和切向刚度,f为裂隙连通率对岩体裂隙变形的影响系数,ri为第i条裂隙端部对该裂隙变形的影响系数,L2为岩块尺寸。
接着分析上述裂隙岩体结构力学模型中裂隙岩体沿加载方向的变形模量表示为:
令不连续变形相似力学模型各材料参数、物理参数设置与裂隙岩体结构力学模型相同,并分别用上角标“′”区分表示,设置与裂隙岩体结构力学模型相同的边界条件,推导不连续变形相似力学模型对应的不连续变形模量方程如下:
其中,E′1为裂隙岩体沿加载方向变形模量,l′i和α′i分别为不连续变形相似力学模型的第i条裂隙长度和倾角,和分别为不连续变形相似力学模型的裂隙结构面法向刚度和切向刚度,L′1×L′2为不连续变形相似力学模型的岩块尺寸。
假设裂隙裂隙倾角保持不变,即αi=α′i,裂隙岩体物理、几何相似常数如下:
将公式(6)代入到公式(5)中,相似变换得到
为保证公式(7)与公式(4)一致性,需要同时满足以下条件:
将公式(6)代入到公式(8),通过上述相似变换推导反应裂隙岩体结构变形特性的模量相似准则如下:
其中,E1为裂隙岩体沿加载方向变形模量,E0为岩石基质杨氏模量,和分别为裂隙结构面法向刚度和切向刚度,和分别为不连续变形相似力学模型的裂隙结构面法向刚度和切向刚度,L1×L2为裂隙岩体力学模型的岩块尺寸,L′1×L′2为不连续变形相似力学模型的岩块尺寸,li和αi分别为第i条裂隙长度和倾角,l′i和α′i分别为不连续变形相似力学模型的第i条裂隙长度和倾角。
最后,根据所述相似准则采用物理模型试验方法对裂隙岩体进行模拟,能够更好的模拟工程实际中大尺度复杂裂隙网络岩体的结构形式与受载情况,为物理模型试验方法提供指导,实现对裂隙岩体更准确的模拟。
为了更清楚的对本实施例提供的方法进行说明,以具体某裂隙岩体结构力学模型为例对本实施例所提供的方法的合理性进行验证。
根据实际工程裂隙岩体中裂隙分布规律,借助UDEC离散元软件建立2.83m×2.83m二维贯通裂隙岩体结构力学模型,裂隙密度为10,裂隙倾角呈45°-90°正态分布,为了保证各项异性裂隙岩体在各个方向上原始模型(裂隙岩体结构力学模型)与相似模型(不连续变形相似力学模型)变形模量的相似性,采用如图4所示二维裂隙岩体不连续变形模量计算原始模型,中间以不同角度旋转截取尺寸为2m×2m岩块,进行轴向压缩数值模拟试验,并以10倍几何扩大相似模型进行不连续变形模量相似准则验证。为了方便验证,令 其中CL=0.1,取裂隙法向刚度、切向刚度相同,统一记作验证结果不失一般性。裂隙岩体结构力学模型(图、表中简称模型)具体物理参数见下表1。
表1裂隙岩体结构力学模型物理参数
图5所示为该二维贯通裂隙岩体结构力学模型轴向加载条件下不同方向变形模量与对应10倍几何扩大不连续变形相似力学模型的变形模量之间的比较,图5(a)中岩石基质弹性模量保持不变,即由公式(8)得到 通过公式(9)得到裂隙岩体10倍几何扩大不连续变形相似力学模型对应的物理参数见表2。
表2 10倍几何扩大不连续变形相似力学模型物理参数
图5(b)中,不连续变形相似力学模型岩石基质材料弹性模量值减弱,即及时,通过公式(8)和公式(9)得到变形模量相似模型对应的物理参数,以上相似模型参数是依据本实施例提出的不连续变形相似力学模型建立方法得到的,三种相似模型对应的相似常数见表3。
表3三种相似模型对应的相似常数
图5(a)采用10倍纯几何扩大模型作为对照模型,结果显示裂隙岩体原始模型与相似模型各个方向变形模量具有很好的相似性,与纯几何扩大模型差异性显著。图5(b)显示岩石基质模量变化时,相似模型与原始模型变形模量复合本发明不连续变形模量相似准则。图5数值计算结果证明本发明提供的方法能够获得与原始模型相似的变形模量,本实施例提供的不连续变形相似模型建立方法合理。
实施例2
如图6所示,本实施例提供了一种裂隙岩体模拟系统,所述系统包括:
原始模型模块M1,用于构建裂隙岩体结构力学模型;
不连续变形模量模块M2,用于根据所述裂隙岩体结构力学模型计算反映裂隙岩体的不连续变形模量方程;
相似模型模块M3,用于根据所述裂隙岩体结构力学模型建立不连续变形相似力学模型;
相似准则模块M4,用于根据所述不连续变形模量方程和所述不连续变形相似力学模型,利用相似变换得到反映裂隙岩体变形特性的相似准则;
模拟模块M5,用于根据所述相似准则采用物理模型试验方法对裂隙岩体的变形情况进行模拟。
本说明书中每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (9)
1.一种裂隙岩体模拟方法,其特征在于,所述方法包括:
构建裂隙岩体结构力学模型;
根据所述裂隙岩体结构力学模型计算反映裂隙岩体的不连续变形模量方程;
根据所述裂隙岩体结构力学模型建立不连续变形相似力学模型;
根据所述不连续变形模量方程和所述不连续变形相似力学模型,利用相似变换得到反映裂隙岩体变形特性的相似准则;所述相似准则包括:
其中,E1为裂隙岩体沿加载方向变形模量,E0为岩石基质杨氏模量,和分别为裂隙结构面法向刚度和切向刚度,和分别为不连续变形相似力学模型的裂隙结构面法向刚度和切向刚度,L1×L2为裂隙岩体力学模型的岩块尺寸,L′1×L′2为不连续变形相似力学模型的岩块尺寸,li为第i条裂隙长度,l′i为不连续变形相似力学模型的第i条裂隙长度,E′1为裂隙岩体沿加载方向变形模量,E′0为岩石基质杨氏模量;
根据所述相似准则采用物理模型试验方法对裂隙岩体的变形情况进行模拟。
2.根据权利要求1所述的一种裂隙岩体模拟方法,其特征在于,所述构建裂隙岩体结构力学模型包括:
根据二维平面应力规律,建立轴向加载条件下二维贯通裂隙岩体结构力学模型;
根据所述二维贯通裂隙岩体结构力学模型建立相应的裂隙岩体结构力学模型。
3.根据权利要求2所述的一种裂隙岩体模拟方法,其特征在于,所述根据所述二维贯通裂隙岩体结构力学模型建立相应的裂隙岩体结构力学模型包括:
根据所述二维贯通裂隙岩体结构力学模型中裂隙端部约束影响以及裂隙连通率影响建立相应的裂隙岩体结构力学模型。
4.根据权利要求1所述的一种裂隙岩体模拟方法,其特征在于,根据所述裂隙岩体结构力学模型计算反映裂隙岩体的不连续变形模量方程包括:
根据所述裂隙岩体结构力学模型建立轴向加载位移方程;
根据所述轴向加载位移方程得到反映裂隙岩体的不连续变形模量方程。
7.根据权利要求1所述的一种裂隙岩体模拟方法,其特征在于,所述根据所述裂隙岩体结构力学模型建立不连续变形相似力学模型包括:
根据所述裂隙岩体结构力学模型的材料参数、物理参数与边界条件推导得到所述不连续变形相似力学模型。
9.一种裂隙岩体模拟系统,其特征在于,所述系统包括:
原始模型模块,用于构建裂隙岩体结构力学模型;
不连续变形模量模块,用于根据所述裂隙岩体结构力学模型计算反映裂隙岩体的不连续变形模量方程;
相似模型模块,用于根据所述裂隙岩体结构力学模型建立不连续变形相似力学模型;
相似准则模块,用于根据所述不连续变形模量方程和所述不连续变形相似力学模型,利用相似变换得到反映裂隙岩体变形特性的相似准则;所述相似准则包括:
其中,E1为裂隙岩体沿加载方向变形模量,E0为岩石基质杨氏模量,和分别为裂隙结构面法向刚度和切向刚度,和分别为不连续变形相似力学模型的裂隙结构面法向刚度和切向刚度,L1×L2为裂隙岩体力学模型的岩块尺寸,L′1×L′2为不连续变形相似力学模型的岩块尺寸,li为第i条裂隙长度,l′i为不连续变形相似力学模型的第i条裂隙长度,E′1为裂隙岩体沿加载方向变形模量,E′0为岩石基质杨氏模量;
模拟模块,用于根据所述相似准则采用物理模型试验方法对裂隙岩体的变形情况进行模拟。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110623951.4A CN113361099B (zh) | 2021-06-04 | 2021-06-04 | 一种裂隙岩体模拟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110623951.4A CN113361099B (zh) | 2021-06-04 | 2021-06-04 | 一种裂隙岩体模拟方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113361099A CN113361099A (zh) | 2021-09-07 |
CN113361099B true CN113361099B (zh) | 2022-08-26 |
Family
ID=77532107
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110623951.4A Active CN113361099B (zh) | 2021-06-04 | 2021-06-04 | 一种裂隙岩体模拟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113361099B (zh) |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101899972A (zh) * | 2010-06-30 | 2010-12-01 | 中国石油大学(北京) | 砂岩油藏注水开发可预测物理模型建立方法 |
CN103940666A (zh) * | 2014-03-18 | 2014-07-23 | 中国矿业大学 | 一种模拟断续裂隙岩石力学特性的细观参数确定方法 |
CN105424435B (zh) * | 2015-11-20 | 2017-11-10 | 中国矿业大学 | 一种制作三维节理裂隙圆柱型岩石试样的方法 |
CN107368645A (zh) * | 2017-07-17 | 2017-11-21 | 华东交通大学 | 一种约束阻尼层结构振动计算方法 |
CN108363827B (zh) * | 2017-12-28 | 2020-04-17 | 清华大学 | 一种基于相似理论的混联机床静刚度分析方法 |
CN111159794A (zh) * | 2018-11-08 | 2020-05-15 | 中国石油大学(华东) | 多裂隙类岩石试样力学性质的几何损伤流变分析方法 |
-
2021
- 2021-06-04 CN CN202110623951.4A patent/CN113361099B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113361099A (zh) | 2021-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | Calibration of the continuous surface cap model for concrete | |
Nandwana et al. | Detection of the location and size of a crack in stepped cantilever beams based on measurements of natural frequencies | |
CN111666699B (zh) | 基于rev全区域覆盖的岩体工程跨尺度模拟计算方法 | |
CN104316388A (zh) | 一种对各向异性材料结构件进行疲劳寿命测定的方法 | |
CN111062162B (zh) | 一种岩土材料精确本构模型的数值建模与应用方法 | |
Huespe et al. | Crack models with embedded discontinuities | |
Taqieddin | Elasto-plastic and damage modeling of reinforced concrete | |
CN107688691B (zh) | 一种木结构用齿板连接节点受力性能的数值模拟方法 | |
CN113361099B (zh) | 一种裂隙岩体模拟方法及系统 | |
CN117171864A (zh) | 一种梁结构线性振动预测方法 | |
CN112100873A (zh) | 水工建筑物承载力的确定方法 | |
Pu et al. | A two-level nesting smoothed extended meshfree method for static and dynamic fracture mechanics analysis of orthotropic materials | |
Liu et al. | Numerical shear tests on the scale effect of rock joints under CNL and CND conditions | |
CN114460267B (zh) | 基于节理面粗糙度的井壁稳定模拟方法及装置 | |
CN113987816A (zh) | 基于蒙特卡罗随机抽样和主成分分析的隧道顶板监测方法 | |
Aktas et al. | Computer-Aided Mold Design Algorithm for Precast Concrete Elements. | |
CN118094998B (zh) | 一种动静组合应力下隧道围岩的稳定性评价方法 | |
Lin et al. | Dynamic crack propagation analysis using scaled boundary finite element method | |
Mousavian et al. | Limit state approach for structurally informed design of shells composed of interlocking blocks | |
Luo et al. | Influence of the Order Exchange of the Node Connection in the Force Analysis of Steel Structures | |
Misir et al. | Numerical Model Calibration and a Parametric Study Based on the Out-Of-Plane Drift Capacity of Stone Masonry Walls. Buildings 2023, 13, 437 | |
Purnama et al. | Effect of mesh element shape in determining the convergent model of modular retaining wall with 3-D numerical simulation | |
CN114943164A (zh) | 一种基于最小耗能原理的储层裂缝定量预测方法 | |
CN117350192A (zh) | 一种原生裂隙水压致裂全应力张量反演方法及系统 | |
Neogy et al. | Effect of partial depth cracks on the natural frequency of twisted blades: a 3-D finite element analysis |
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 |