单轴压缩下横观各向同性岩石破裂过程的数值模拟
16岩石单轴压缩试验原始记录 1
平均值（MPa）
编号
荷载（kN）
1
轴向应力（MPa）
△L1/(L= )
△L2/(L= )
平均纵向应变10-2
2
轴向应力（MPa）
△L1/(L= )
△L2/(L= )
平均纵向应变10-2
3
轴向应力（MPa）
△L1/(L= )
△L2/(L= )
平均纵向应变10-2
编号
弹性模量（MPa）
岩石单轴压缩变形试验记录表
委托编号
样品编号
委托日期
岩石种类
岩石规格
检测日期
使用部位
岩石状态
岩石产地
试验依据
SL/T264-2020《水利水电工程岩石试验规程》
仪器设备
30T电液伺服万能试验机YHZ-JC-009
岩石单轴抗压强度（MPa）
底面直径（mm）
顶面直径（mm）
试件高（mm）
平均直径（mm）
破坏荷载（N）
变形模量（MPa）
单值
均值
单值
均值
1
2
3
备注
[说明]1.计算公式： ; ;
式中 ---应力，MPa; P---荷载，N；A---试件截面面积，mm2;Ee----岩石弹性模量，MPaμe----岩石弹性泊松比；σa---应力与纵向应变关系曲线上直线段起始点的应力值，MPa；σb---应力与纵向应变关系曲线上直线段终点的应力值，MPa;εha---应力为σa时的纵向应变值；εhb---应力为σb时的纵向应变值；εda---应力为σa时的横向应变值；εdb---应力为σa时的横向应变值；E50----岩石变形模量，即割线模量，MPa;σ50---抗压强度50%时的应力值，MPa;
岩石单轴压缩、拉伸、巴西劈裂数值实验模拟
2.1 软件的基本原理
RFPA 是一个以弹性力学为应力分析工具、以弹性损伤理论及其修正后的 Coulomb 破坏准则为介质变形和破坏分析模块的真实破裂过程分析系统。 其基本 思路是： 1）材料介质模型离散化成由细观基元组成的数值模型，材料介质在细观上 是各向同性的弹-脆性或脆-塑性介质； 2）假定离散化后的细观基元的力学性质服从某种统计分布规律（如 weibull 分布），由此建立细观与宏观介质力学性能的联系； 3）按弹性力学中的基元线弹性应力、应变求解方法，分析模型的应力、应 变状态。RFPA 利用线弹性有限元方法作为应力求解器； 4）引入适当的基元破坏准则（相变准则）和损伤规律，基元的相变临界点 用修正的 Coulomb 准则； 5）基元的力学性质随演化的发展是不可逆的； 6）基元相变前后均为线弹性体； 7）材料介质的裂纹扩展是一个准静态过程，忽略因快速扩展引起的惯性力 的影响。
2.2 软件的网格划分
RFPA 选取等面积四节点的四边形单元剖分计算对象。为了使问题的解答足 够精确，RFPA 方法要求模型中的单元足够小（相对于宏观介质），以能足够精
确的地反映介质的非均匀性。但它又必需足够大（包含一定数量的矿物和胶结物 颗粒，以及微裂隙、孔洞等细小缺陷），因为作为子系统的单元实际上仍是一个 自由度很大的系统，它具有远大于微观尺度的细观尺度。这以要求正是为了保证 使剖分后的单元性质尽量接近基元性质。尽管这样会增加计算量，但是问题的处 理变得简单， 而且随着计算机技术的高速发展， 计算机瓶颈的影响将会逐渐消除。 由于模型中的基元数量足够多，宏观的力学行为，本质上是介质大量基元力学行 为的集体效应。
（a）step42-01
（b）step52-02
（c）step70-06
图 3、RFPA 模拟单轴拉伸条件下的破坏过程、最大主应力场、声发射累计分布图
岩石动态剥落破裂的数值模拟
岩石动态剥落破裂的数值模拟引言岩石动态剥落破裂是地质灾害中的一种严重类型，其产生的原因多样，如地震、爆炸、水力冲击等。
对于这种问题，数值模拟方法已被广泛应用于地质工程领域，以预测和评估岩石动态破裂过程的破坏性和具体效果，以及结构的稳定性和保护性能。
本文将介绍目前常用的岩石动态破裂数值模拟方法，包括有限元法和离散元法，并分析其优劣和应用范围。
一、有限元法有限元法是解决结构力学中的问题的常用方法，包括岩石动态破裂模拟。
其基本思想是将复杂的结构分解成若干个小元素，并对每个小元素进行简化模型假设，利用数值方法对每个小元素进行求解，最后将结果组合得到全局结构的反应。
在岩石动态破裂模拟中，将峰值强度、应力波传播、岩石内损伤等问题转化为有限元数值求解问题，可大幅简化问题的求解过程。
有限元法在岩石动态破裂模拟中的应用主要涉及到以下几个方面：1、破裂过程的数值模拟：破裂过程的分析对于预测和评估破坏的具体情况至关重要，有限元法能够对破裂过程进行数值模拟；2、弹性介质中应力波传播的数值模拟：应力波传播的速度、频率对于岩石破裂具有重要影响，有限元法可以计算弹性介质中应力波传播的特征及其影响；3、岩石内部损伤行为的数值模拟：岩石内部微观结构的变化对于破裂行为的发生有着直接的影响，有限元法可以模拟并计算微观尺度上的变化。
有限元法的优点在于：1、求解过程简便快捷；2、可对各种不同类型和形状的结构进行模拟；3、适用于各种不同工况下的模拟。
其缺点在于：1、仅适用于小小尺度下，如旋转对称或轴对称问题的处理等；2、计算机资源投入较大，对于大规模结构的处理难度较大；3、需要对于每个小元素进行较好的建模。
二、离散元法离散元法是一种分子动力学模型，其首要任务是模拟模型中各种物质颗粒在自然环境下的运动行为，其模型假设是颗粒物的弹性和摩擦不存在。
离散元法最初被应用于地质动力学的问题中，由于其适用范围广、计算速度快、能够对多种不同类型的物体进行建模等优点，迅速成为岩石动态破裂模拟中最常用的方法之一。
基于双屈服条件准则的横观各向同性本构模型研究及其数值模拟
基于双屈服条件准则的横观各向同性本构模型研究及其数值模拟QU Guangxiu;REN Peng【摘要】为研究层状岩体的力学特性,提出基于双屈服条件强度准则的本构模型.基于双屈服条件强度准则,联合横观各向同性的广义虎克定律刚度矩阵建立考虑横观各向同性的本构模型,并结合岩石单轴压缩试验数据,通过最小二乘法拟合该模型的参数;实现该模型的单轴压缩试验数值模拟,并通过室内单轴压缩试验结果对数值模拟结果进行验证,分析模型的可靠性.研究结果表明:本文提出的本构模型在层状岩体的力学分析方面具有适用性,为层状岩体力学特性研究及层状岩质边坡的稳定性分析奠定了基础.【期刊名称】《铁道科学与工程学报》【年(卷),期】2019(016)006【总页数】6页(P1448-1453)【关键词】横观各向同性;本构模型;双屈服条件强度准则;数值模拟【作者】QU Guangxiu;REN Peng【作者单位】【正文语种】中文【中图分类】TU458层状岩质边坡广泛分布于我国西南地区，其明显的横观各向同性力学特性对边坡的稳定性有着显著影响，因此如何建立适用的本构模型以探究其力学行为具有重要的工程实践意义。
关于横观各向同性岩石的本构模型研究，国内外学者进行了大量研究。
刘运思等[1]通过室内试验对横观各向同性岩体的弹性参数进行了研究。
Gonzaga等[2]通过三轴压缩试验研究了如何确定横观各向同性岩石的力学参数。
ZHANG等[3−5, 11]通过不同试验手段研究了横观各向同性岩石的破坏机理，探讨了加载速率对破坏过程的影响。
熊良宵等[6−8]采用数值模拟方法，探讨了横观各向同性岩体的力学特性。
Colak等[9−10]对横观各向同性岩体的破坏强度准则进行了研究。
上述研究成果大都基于Hoek-Brown准则，描述横观各向同性岩体的强度和变形特征，并提出不同的强度准则和弹塑性本构模型，但大多研究成果仅从强度或者变形特征这种单一因素考虑横观各向同性岩体的本构模型，如何科学地描述层状岩石的强度和变形特征仍值得商榷。
单轴压缩过程中岩石变形破坏机理
单轴压缩过程中岩石变形破坏机理摘要：岩石是地球表层的主要构成物质，其变形破坏机理是地质学、地球物理学、岩土力学等领域研究的重要内容之一。
本文主要围绕单轴压缩过程中岩石的变形和破坏机理进行探讨，从微观和宏观层面分析了岩石变形破坏的机理，对于加深对岩石力学性质的认识和岩体工程设计具有一定的指导意义。
关键词：单轴压缩；岩石变形；岩石破坏；机理一、引言岩石是地球表层的主要构成物质，其地质力学性质对于地质灾害防治、矿产资源勘探和开采、地下工程建设等具有重要意义。
岩石的变形和破坏是岩石力学研究的核心内容之一，而单轴压缩实验是岩石力学中最基本、最重要的实验之一，能够模拟岩石在地球内部受到的单向应力，揭示岩石的变形和破坏机理。
本文将围绕单轴压缩过程中岩石的变形和破坏机理进行探讨，从微观和宏观层面分析了岩石变形破坏的机理，对于加深对岩石力学性质的认识和岩体工程设计具有一定的指导意义。
二、单轴压缩实验概述单轴压缩实验是岩石力学中最基本、最重要的实验之一，其主要原理是使岩石试样在一个方向上受到单向压缩应力，从而模拟岩石在地球内部受到的单向应力。
实验中通常采用圆柱形或立方体的岩石试样，通过加压机施加压力，测量岩石试样的应变和应力变化，从而研究岩石的变形和破坏机理。
三、岩石变形机理岩石的变形是指岩石在受到外力作用下发生的形状和大小的变化，包括弹性变形和塑性变形两种形式。
岩石的破坏是指岩石在受到外力作用下失去原有的力学性质和稳定性，出现破裂和破碎的现象。
1. 岩石的弹性变形岩石的弹性变形是指岩石在受到外力作用下，能够在一定范围内恢复原有的形状和大小的变化。
在单轴压缩实验中，当岩石试样受到一定的压力时，岩石试样会发生弹性变形。
当压力撤离时，岩石试样会恢复原有的形状和大小。
2. 岩石的塑性变形岩石的塑性变形是指岩石在受到外力作用下，发生不可逆的形状和大小的变化。
在单轴压缩实验中，当岩石试样受到足够大的压力时，岩石试样会发生塑性变形。
S01-45 岩石单轴压缩变形试验
岩石单轴压缩变形试验
Ｓ０１－４５ 施工单位： 监理单位： 试验单位 工程部位 纵向变形(kN) 级数 等级 荷载 (kN) 每级 应力 (MPa) 1 平均 每级 应变 10 2 平均 每级 应变 10 3 平均 每级应 变10 1 平均 每级应 变10 2 合 同 号： 试验编号： 试验日期 试验依据 横向变形(kN) 平均 每级应 变10 3 平均 每级应 变10
六级
七级
八级
九级
十级 轴心抗压强度 （MPa） 单值 平均值 抗压弹性模量 （10MPa） 单值 平均值 单值 平均值
泊松比 监理人员意见：
施工单位自检意见： 结论 日期： 试验： 复核： 试验负责人：
签字：
日期：
双裂隙复合岩层单轴压缩力学性质及损伤机理离散元模拟
第 54 卷第 10 期2023 年 10 月中南大学学报(自然科学版)Journal of Central South University (Science and Technology)V ol.54 No.10Oct. 2023双裂隙复合岩层单轴压缩力学性质及损伤机理离散元模拟李恒1，杨圣奇1, 2，孙博文2，杨景2，董志锦2，温森3(1. 中国矿业大学 力学与土木工程学院，江苏 徐州，221116；2. 中国矿业大学 深部岩土力学与地下工程国家重点实验室，江苏 徐州，221116；3. 河南大学 土木建筑学院，河南 开封，475004)摘要：针对存在内部裂隙的软硬互层复合岩层力学特性展开研究，采用PFC 2D 建立层状复合岩层离散元模型，通过对完整复合岩层试验细观参数进行校准并验证其宏观特性，研究分析平行双裂隙复合岩层的岩桥倾角与层理倾角对力学特性与破坏模式的影响，得出双裂隙复合岩层的5类最终破坏模式，揭示裂纹演化规律及位移场变化。
研究结果表明：双裂隙复合岩层峰值强度随岩桥倾角增加呈现先减小后增加的趋势，而随层理倾角θ增大表现为“U ”型变化规律；与完整试样相比，双裂隙复合岩层峰值强度降低4.4%~42.0%；双裂隙复合岩层破坏模式主要分为5种类型，分别为轴向张拉破坏(Ⅰ型)、穿层理剪切破坏(Ⅱ型)、穿层理拉伸破坏(Ⅲ型)、沿层理面滑移破坏(Ⅳ型)、沿岩层拉伸破坏(Ⅴ型)；试样的破坏模式在θ=0°~45°时随岩桥倾角增加从Ⅰ型转变为Ⅱ型最后转变为Ⅲ型破坏模式，在θ=60°~75°和θ=90°时分别表现为Ⅳ型和Ⅴ型破坏模式；在θ=0°~45°时，翼裂纹垂直于主裂隙起裂，随岩桥倾角增大，翼裂纹沿岩桥倾角偏转并扩展成宏观破裂面，导致试样最终被破坏；在θ=60°~75°时，裂纹在层理面扩展形成宏观滑移面，导致试样破坏；在θ=90°时裂纹沿层理面发生岩桥贯通，翼裂纹在岩层内扩展形成宏观破裂面并导致试样破坏。
《2024年三维条件下的岩石破裂过程分析及其数值试验方法研究》范文
《三维条件下的岩石破裂过程分析及其数值试验方法研究》篇一一、引言岩石破裂过程是地质学、岩土工程学等领域的重要研究内容。
随着科技的发展，三维条件下的岩石破裂过程分析及其数值试验方法研究逐渐成为研究热点。
本文旨在探讨三维条件下的岩石破裂过程，并对其数值试验方法进行研究，以期为相关领域的研究提供参考。
二、岩石破裂的基本理论岩石破裂是指岩石在受到外力作用时，内部应力超过其承受极限，导致岩石结构破坏的现象。
岩石的破裂过程受到多种因素的影响，如岩石的物理性质、地质构造、应力状态等。
在三维条件下，岩石的破裂过程更为复杂，需要考虑多方向应力、温度、湿度等因素的影响。
三、三维条件下的岩石破裂过程分析1. 理论模型：在三维条件下，岩石的破裂过程可以通过建立理论模型进行分析。
常用的模型包括弹性力学模型、塑性力学模型、断裂力学模型等。
这些模型可以描述岩石在不同条件下的应力、应变及破裂过程。
2. 数值模拟：通过数值模拟方法，可以更直观地了解岩石的破裂过程。
常用的数值模拟方法包括有限元法、离散元法、边界元法等。
这些方法可以模拟岩石在不同条件下的应力分布、裂纹扩展及最终破裂形态。
3. 实验观察：通过实验观察，可以更深入地了解岩石的破裂过程。
实验方法包括光学显微镜观察、电子显微镜观察、声发射监测等。
这些方法可以观察岩石的微观结构变化、裂纹扩展及破裂模式。
四、数值试验方法研究1. 有限元法：有限元法是一种常用的数值试验方法，通过将岩石划分为有限个单元，求解每个单元的应力、应变及位移等参数，从而分析岩石的破裂过程。
该方法可以处理复杂的几何形状和边界条件，但需要较高的计算成本。
2. 离散元法：离散元法是一种基于刚性块体运动的数值试验方法，适用于模拟岩体的不连续性。
该方法通过考虑块体间的接触和相互作用，模拟岩体的破裂过程。
其优点是可以处理大变形和失稳问题，但需要较细致的模型构建。
3. 结合实际：在实际应用中，可以根据研究目的和岩石性质选择合适的数值试验方法。
实验五岩石单轴压缩实验
实验五岩石单轴压缩实验一、实验目的1、掌握岩石单轴压缩实验的基本原理和方法。
2、了解岩石单轴压缩强度的测试方法。
3、通过实验了解岩石在不同压力下的力学性质。
4、了解矿山工程中用于确定岩石层强度、稳定性和采矿方法选择的基本实验方法。
二、实验原理在实验室条件下，对岩石进行单轴压缩实验，即将岩石样品置于压力机滑动块与固定块之间，施压加荷，岩石样品在压力的作用下发生变形，最终出现破裂破坏。
这种实验方法可以测定岩石样品在单轴压缩应力下断裂时的应力水平值和断裂模式，是评估岩石力学性质和确定其强度和稳定性的重要方法。
单轴压缩强度表示岩石样品在单轴压缩下破坏时的最大承受压力或应力水平。
在实验过程中，将岩石样品沿其轴向方向施以单向的压力，直到样品发生破坏，根据压力与样品断面积之比计算出样品的单轴压缩强度。
单轴压缩实验中常用的岩石模型为标准直径为50mm、高度为100mm、直径与高度比为1:2的圆柱形样品。
通过实验获取不同压力下岩石样品的应变和应力的数据，利用数据处理方法分析出样品的单轴压缩强度和岩石在不同压力下的变形和断裂模式。
三、实验步骤1、制备标准圆柱形样品在实验之前，制备标准的圆柱形样品是非常重要的，直径为50mm、高度为100mm，或者根据实际情况选择其他规格的样品。
2、测量标准圆柱形样品尺寸使用卷尺对样品的直径和高度进行测量，并记录下来，便于计算样品的断裂强度。
3、安装压力机将样品放置在压力机的压缩板上，并确保样品在滑动块与固定块之间完全垂直。
调整滑动块的位置，使其与样品顶部接触。
将固定块和滑动块夹紧，用气动或手动方式施压。
4、开始施压施加压力，开始进行单轴压缩实验，随着施压的增加，记录下实验的每一阶段应变和应力数据。
5、记录数据根据实验数据绘制出应力-应变曲线、应力-时间曲线，计算出单轴压缩强度。
模拟分析样品破裂模式。
6、进行岩石单轴压缩实验的注意事项a、施压过程应逐步增加，避免突然增压，以免样品产生损伤。
混凝土单轴压缩破坏过程的三维细观数值模拟
ａ ｔ ｍａ ｉａｌ ｎ ｖ ｌａ ｅ ．Ｔｈ ｃ ａ ｉａ ｅ ｆｒ ｎ ｅａ ｄ ｆｉ ｒ ｒ ｃ ｓ ｆｃ ｎ ｒ ｔ ｕ ｊｃ ｕ ｏ ｔ ｌ ａ ｄ ｅ ａｕ ｔ ｄ ｃ ｙ ｅｍｅ ｈ ｎ ｃｌｐ ｒｏ ｍａ ｃ ｎ ａｌ ｅ ｐ ｏ ｅ ｓ ｏ ｏ ｃｅ ｅｓ ｂ ｅ ～ ｕ
２ Ｃｈｎ ｎ ｔ ｕ ｅｏ ａｅ ｓ ｕ ｃｓａ ｄ Ｈｙ ｒｐ ｗｅ ｓａ ｃ ， ｅｊ ｇ １ ０ ４ Ｃｈｎ ． ｉａ Ｉｓｉ ｔ ｆＷ ｔｒＲｅ ｏ ｒｅ ｎ ｄ ｏ ｏ ｒＲｅｅ ｒ ｈ Ｂ ｉｎ ０ ０ ４， ｉａ） ｔ ｉ Ａｂ ｔ ａ ｔ ｓ ｒ ｃ ：Ｂａ ｅ ｎ ｍ ｅ ｏ ｍｅ ｈａ ｃ ｓ ｄ ｏ ｓ — ｃ ｎｉｓ，ｃ ｃ ｅ ｅ ｉ ａ ｎ ａ ｈｒ ｅ ｐ ｓ ｏｍｐ ｉｅ ｏｎ ｉｔｎ ｇ— ｏｎ ｒ ｔ ｓｔ ｋｅ ｓｔ ｅ － ｈａ ｅ ｃ ｏｓｔ ｓｃ ｓ ｓ ｉ ｇ ｏｆａ ｇｅ ｔ ｒ ｇａ ｅ，ｍｏ ｔ ｒ ａ ｄ ｂ ｒ ａ ｎ ｏｎｄ ｂｅ ｗｅ ｎ ｇｇ ｅ ａ ｅ ａ ｄ ｍｏ ｔ ｒ Ｔｈｅ ｒ ｎｄ ｍ ｓｔｏ ｆ ａ ｇｒ ｇａ ｅ ｉ ｔ ｅ ａ ｒ ｇ ｔ ｎ ｒ ａ ． ａ ｏ ｐｏ ｉｉ ｎ ｏ ｇ ｅ ｔ ｎ
摘 要 ：从 细观 力 学的 角度 出发 ， 混凝 土视 为 由骨料 、 浆及 二 者之 间的界 面所 组 成 的三 相 复合 材 将 砂 料 ， 用 蒙特卡 罗 方法 ， 生 骨料 在 试 件 中的 随机 位 置 ， 立 了混凝 土 圆柱 体 试件 三 维数 值模 型 。 利 产 建 进 行单 元 划分 ， 对试件 中骨料 、 砂浆 和界 面单 元的材 料 号进行 自动 识别及 赋值 。通 过 三个计 算 方案 研 究 了混凝 土单 轴压 缩 下的 力学性 能及破 坏过 程 ， 拟 了混凝 土材料 从 裂 纹的 萌 生 、 展 、贯通 直 模 扩 到宏观 裂纹 产 生导致 破坏 的过 程 。结 果表 明 ， 该数 值模 型 及计 算 方 法 能够 较 好 地反 映 混 凝 土在 单 轴 荷载 作 用下 的力 学特性 。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
图４是７个岩石试件在弹性阶段表现出的弹性 模量和岩层倾角之间的关系。由３，图４可知，虽 然峰值强度发生在口＝６０。，但是最小的弹性模量 却出现在∥＝３０。。∥＝６０。时所表现出的弹性模量 比∥＝０。和∥＝１５。还要大。∥＝７５。和∥＝９０。的 弹性模量虽然相同，但是最后的峰值强度却相差很 大。这说明了横观各向异性岩石在弹性阶段所表现 出的弹性模量的大小并不能代表其最终峰值强度的 大小。
从细观力学的角度来说，岩石晶粒的变形是很 小的，塑性变形主要在于晶粒之间的滑移。岩石弹 塑性变形过程中，其宏观体积响应只是各种微结构 之间相互作用的平均结果，变形的主要机理是微裂 纹成核和增长及其聚集而生成细观裂纹和裂纹的扩 展与传播…Ｊ。但是，由于材料单元中微观实体数目 众多，表现形态各异，要从理论上精确地描述和分 析各个损伤实体几乎是不可能的。然而电子计算机 的出现和数值模拟方法的发展为解决这一难题提供 了新的途径。因此，如果能从细观损伤力学出发， 并充分考虑到岩石的非均匀性，建立一种简单直接 的数值模型来模拟横观各向同性材料整个破裂过 程，对于研究横观各向同性岩体开挖工程的稳定性 分析以及滑坡预防等研究无疑具有重要的理论和实 际意义。本文采用东北大学岩石破裂与失稳研究中 心开发的ＲＦＰＡ２Ｄ软件，用２种不同的岩石材料来 组成不同岩层倾角的横观各向同性的岩石试件，通 过单轴加载数值模拟试验模拟横观各向同性岩石渐 进破裂的整个过程，分析了岩层与最大主应力之间 的倾角和强度之问的关系，研究了不同岩层倾角的
１０
０ ０
０．０００５
０．００１。
Ｏ．００１５
０，００２
图２ ７种不同倾角试件的应力．应变全过程曲线 Ｆｉｇ．２ Ｃｏｍｐｌｅｔｅ盯一￡ｃｕｒｖｅｓ ｏｆ ｓｅｖｅｎ ｒｏｃｋ ｓａｍｐｌｅｓ ｗｉｔｈ
ｄｉｆｆｅｒｅｎｔ ｄｉｐ ａｎｇｌｅｓ
图２是岩层倾角∥分别为０。，１５。，３０。，４５。， ６０。，７５。和９０。的岩石试件在载荷作用下的全过 程应力一应变曲线。首先可以看出，随着岩层倾角的 增大，岩石的强度有一个逐渐减小，然后逐渐增加 的过程。值得注意的是，当口＿ｏ。和／，＝１５。的时候 强度变化很小，当达到６０。时，岩石试件的强度达
横观各向同性现象在地质材料中比较常见，各
向异性性状对岩体的应力一应变分析以及破坏力学 行为有很重要的影响，国内外学者做了很多的研究 工作。最早是Ｊａｅｇｅｒ针对节理岩体沿着节理面滑动 破裂提出了相应的破坏准则…，后来很多外国学者 也提出了很多横观各向同性岩体的破坏准则【２。７ Ｊ。也 有不少学者围绕确定横观各向同性岩体变形特性和 弹性参数做了很多实验研究。Ａｎｎａ采用一种新的三 轴传感器技术得到横观各向同性岩石三维变形特性 和弹性常数【８】。Ｔｉｅｎ等用２种不同的材料人工预制 了３组不同倾角的层状岩石，研究了横观各向同性 体倾角对岩石整体强度和弹性模量的影响【９Ｊ，并针 对横观各向同性岩石破坏的２种不同模式提出了相 应的破坏准则。李同林对煤层横观各向同性体水压
ＲＦＰＡ２Ｄ与其他数值模拟软件不同之点就是能 够模拟岩石类脆性材料在应力作用下，从细观单元 的破裂到裂纹萌生扩展以及宏观破裂的整个过程， 不但可以得到破裂过程中的应力一应变的信息，还能 模拟出相应的应力场和变形场的整个演化过程以及 破裂过程中的声发射的时空分布情况。
５０
４０
３０ 山 羔 ＼ ２０ ｂ
控制参数 参数值 控制参数
参数值
凸一
羔 ＼ 魁 想 趔 彗
钻∞”∞巧加埒
Ｏ
１５
３０
４５
６０
岩层倾角／（。）
７５
９０
图３ ７种岩石试件的峰值强度和层面倾角的关系 Ｆｉｇ．３ Ｐｌｏｔｓ ｏｆ ｐｅａｋ ｓｔｒｅｎｇｔｈ ＶＳ·ｄｉｐ ａｎｇｌｅｓ ｏｆ ｓｅｖｅｎ ｒｏｃｋ
ｓａｍｐｌｅｓ
３·５
３数值模拟结果分析
５８
岩
土
力
学
２００５笠
致裂机理进行了深入的认识和探讨［１ ０１。章根德等讨 论了横观各向同性岩土介质的弹塑性及其屈服准则 的模拟，给出了相应的力学模型及其数学表达式 ｕ¨。张玉军等分析了层状岩体强度异向性对地下洞 室稳定性的影响，并提出了相应的经验公式【ｌ ２｜。由 于在应力的作用下各向异性体所表现出来的各向异 性力学行为比较复杂，所以为了确定横观各向同性 线弹塑性体模型，通常需要确定的材料常数达到２０ 多个，因此需要大量的试验来确定这些常数，包括 不同方向的拉伸、压缩、以及剪切强度试验，这样就 很难在工程实践中应用推广。
１ 概述
在边坡和地下工程的稳定性分析中，一般将其 当作各向同性介质来处理，这种做法对某些均质的 岩体可以达到工程所需要的精度，但是对于一些力 学性质呈明显各向异性的岩体，如层状的沉积岩就 不适用。沉积过程中所形成的层状结构的粘土层， 在各向同性面内各个方向的矿物成份及物理力学性 质是大体相同的，但在垂直此面方向内的力学性质 确有很大差别。类似这样沉积形成的天然地质体， 其变形破坏分析可以采用横观各向同性弹性体模型 来描述。
摘要：采用基于细观损伤力学基础上开发的ＲＦＰＡ２Ｄ数值模拟软件，用２种不同的岩石材料来组成７个不同岩层倾角的横
观各向同性的岩石试件，通过单轴加载数值模拟试验，模拟横观各向同性岩石渐进破裂的整个过程，分析了岩层与最大主应
力之间的倾角和强度之间的关系，讨论了不同岩层倾角的横观各向同性岩体的不同破裂模式及其破坏准则。
到最低值。当肛９０。时强度达到最大值，几乎是最
低强度的２倍，这也说明层状岩体具有强烈的各向 异性。
图３是不同倾角试件强度的曲线图。
万方数据
蛊 ；３．０
曼
咖Ｉ｜
蓑２．５
霎
０
１５
３０
４５
６０
７５
９０
层面倾角／（。）
图４ ７种岩石试件的杨氏模量和层面倾角的关系 Ｆｉｇ．４ Ｐｌｏｔｓ ｏｆ Ｙｏｕｎｇ’Ｓ ｍｏｄｕｌｕｓ ＶＳ．ｄｉｐ ａｎｇｌｅｓ ｏｆ ｓｅｖｅｎ ｒｏｃｋ
万方数据
第１期
梁正召等：单轴压缩下横观各向同性岩石破裂过程的数值模拟
表１组成岩石试件的２种不同介质的力学参数 Ｔａｂｌｅ １ Ｍｅｃｈａｎｉｃａｌ ｐａｒａｍｅｔｅｒｓ ｏｆ ｔｗｏ ｄｉｆｆｅｒｅｎｔ ｒｏｃｋｓ
ｍａｔｅｒｉａｌｓ
表２ ＲＦＰＡ２Ｄ中的其它计算控制参数 Ｔａｂｌｅ ２ Ｏｔｈｅｒ ｃｏｎｔｒｏｌ ｐａｒａｍｅｔｅｒｓ ｉｎ ｓｉｍｕｌａｔｉｏｎ ｉｎ ＲＦＰＡｍ
Ａｂｓｔｒａｃｔ：Ｂａｓｅｄ ｏｎ ｍｅｓｏｓｃｏｐｉｃ ｄａｍａｇｅ ｍｅｃｈａｎｉｃｓ．ｎｕｍｅｒｉｃａｌ ｃｏｄｅ ＲＦＰＡ。“ｉｓ ｄｅｖｅｌｏｐｅｄ ｔｏ ｓｉｍｕｌａｔｅ ｔｈｅ ｆａｉｌｕｒｅ ｐｒｏｃｅｓｓ ｏｆ ｓｅｖｅｎ ｔｒａｎｓｖｅｒｓｅｌｙ ｉｓｏｔｒｏｐｉｃ ｒｏｃｋ ｓａｍｐｌｅｓ ｕｎｄｅｒ ｕｎｉａｘｉａｌ ｌｏａｄｉｎｇ．Ｔｈｅ ｒｏｃｋ ｓａｍｐｌｅｓ ａｒｅ ｃｏｍｐｏｓｅｄ ｏｆ ｔｗｏ ｄｉｆｆｅｒｅｎｔ ｒｏｃｋ ｍａｔｅｒｉａｌｓ ａｎｄ ｔｈｅｙ ａｒｅ ｆｏｒｍｅｄ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ｄｉｐ ａｎｇｌｅｓ ｂｅｔｗｅｅｎ ｔｈｅ ｒｏｃｋ ｌａｙｅｒ ｏｒｉｅｎｔａｔｉｏｎ ａｎｄ ｔｈｅ ｌｏａｄｉｎｇ ｄｉｒｅｃｔｉｏｎ．Ｃｏｍｐｌｅｔｅ ｓｔｒｅｓｓ—ｓｔｒａｉｎ ｃｕｒｖｅｓ ａｒｅ ｏｂｔａｉｎｅｄ ａｎｄ ｔｈｅ ｄｅｆｏｒｍａｔｉｏｎ ａｎｄ ｆａｉｌｕｒｅ ｂｅｈａｖｉｏｒ ａｒｅ ｄｅｓｃｒｉｂｅｄ．Ａｓ ｓｈｏｗｎ ａｓ ｔｈｅ ｎｕｍｅｒｉｃａｌ ｒｅｓｕｌｔｓ，ｒｏｃｋ ｌａｙｅｒ ｄｉｐ ａｎｇｌｅ ｏｆ
ｔｒａｎｓｖｅｒｓｅｌｙ ｉｓｏｔｒｏｐｉｃ ｒｏｃｋｓ ｈａｓ ｍｕｃｈ ｇｒｅａｔ ｉｎｆｌｕｅｎｃｅ ｏｎ ｔｈｅ ｆａｉｌｕｒｅ ｐｒｏｃｅｓｓ，ｓｕｃｈ ａｓ ｐｅａｋ ｓｔｒｅｎｇｔｈ，ｆａｉｌｕｒｅ ｍｏｄｅｓ，ａｎｄ ｄｅｆｏｒｍａｔｉｏｎ
ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｅｔ ａ１．Ｉｔ ｉｓ ｓｕｉｔａｂｌｅ ｔｏ ａｐｐｌｙ ｄｉｆｆｅｒｅｎｔ ｆａｉｌｕｒｅ ｃｒｉｔｅｒｉａ ａｃｃｏｒｄｉｎｇ ｔｏ。ｄｉｆｆｅｒｅｎｔ ｆａｉｌｕｒｅ ｍｏｄｅｓ ｃａｕｓｅｄ ｂｙ ｌａｙｅｒ ｄｉｐ ａｎｇｌｅ．Ｔｈｅ ｒｅｓｕｌｔｓ ｏｆ ａｎｏｔｈｅｒ ｓｅｒｉａｌ ｏｆ ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎｓ ｉｓ ｗｅｌｌ ａｇｒｅｅｍｅｎｔ ｗｉｔｈ ｔｈｅ ｅｘｐｅｒｉｍｅｎｔａｌ ｓｔｕｄｉｅｓ．Ｔｈｅ ｄｉｆｆｅｒｅｎｃｅ ｂｅｔｗｅｅｎ ｔｈｅ ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎｓ ａｎｄ ｅｘｐｅｒｉｍｅｎｔｓ ｉｓ ｐｒｏｂａｂｌｙ ｃａｕｓｅｄ ｂｙ ｐｌａｎｅ ｓｔｒａｉｎ ｐｒｏｂｌｅｍｓ ａｎｄ ｐａｒａｍｅｔｅｒｓ ｕｓｅｄ ｉｎ ｔｈｅ ｓｉｍｕｌａｔｉｏｎｓ． Ｋｅｙｗｏｒｄｓ：Ｔｒａｎｓｖｅｒｓｅｌｙ ｉｓｏｔｒｏｐｉｃ ｒｏｃｋ；ｆａｉｌｕｒｅ ｐｒｏｃｅｓｓ；ｕｎｉａｘｉａｌ ｃｏｍｐｒｅｓｓｉｏｎ；ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ
８－－０。
８＝３０。
踞４５。
８＝９０。
图１ ２种不同介质所组成的７个不同岩层倾角的横观各向同性岩石试件
Ｆｉｇ．１ Ｓｅｖｅｎ ｔｒａｎｓｖｅｒｓｅｌｙ ｉｓｏｔｒｏｐｉｃ ｒｏｃｋ ｓａｍｐｌｅｓ ｃｏｍｐｏｓｅｄ ｏｆ ｔｗｏ ｒｏｃｋ ｍａｔｅｒｉａｌｓ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ｄｉｐ ａｎｇｌｅｓ
通常的实验室物理岩石力学实验中，由于钻取 方向和岩体本身非均匀性的存在，取得足够的不同 方位、不同角度的横观各向同性岩芯比较困难。一 般对于横观各向同性线弹性体来说，虽然由于旋转 对称性减少了独立的弹性常数的数目，其独立的弹 性常数变为５个，但是对于横观各向同性弹塑性材 料，为了确定横观各向同性材料的屈服条件，还需 要确定很多参数，共需要２０多个材料常数才能完全 确定本模型所描述的性状。但在ＲＦＰＡ’２Ｄ中，用较 少的基元力学参数就可以达到很方便设置各种不同 岩石力学模型，图１是ＲＦＰＡ２Ｄ设置的不同岩层倾 角１３的岩石试件。岩石试件中白色条带是强度较高 的岩石介质Ａ，较暗的条带是强度较低的岩石介质 Ｂ，由这２种介质组成了不同角度层理的岩石试件。 ２种不同岩石介质细观基元的力学参数见下表１。