多裂隙类岩石试样力学性质的几何损伤流变分析方法技术

本发明专利技术属于土木工程领域，具体地，涉及到多裂隙类岩石试样力学性质的几何损伤流变分析方法。步骤如下：(1)、构建多裂隙类岩石试样代表性单元的几何模型；(2)、定义多裂隙类岩石试样代表性单元的损伤张量；(3)、计算多裂隙类岩石试样中单组裂隙的损伤张量；(4)、计算多裂隙类岩石试样中多组裂隙的损伤张量；(5)、推导多裂隙类岩石试样有效应力公式；(6)、建立多裂隙类岩石试样几何损伤流变本构模型；(7)、多裂隙类岩石试样几何损伤流变本构模型程序化；(8)、统计工程岩体的几何裂隙参数，进行工程岩体的损伤流变特性数值模拟。

Geometric damage rheological analysis method for mechanical properties of rock samples with multiple fractures

【技术实现步骤摘要】
多裂隙类岩石试样力学性质的几何损伤流变分析方法
本专利技术属于土木工程岩土工程领域，具体地，涉及到多裂隙类岩石试样力学性质的几何损伤流变分析方法。
技术介绍
裂隙在岩体中广泛分布，对岩体结构的稳定性有着重要的影响，而如何评估多裂隙岩体的损伤影响，尤其是岩体的长期行为是一个重要问题。了解多裂隙岩体的力学性质在岩石工程设计中都非常重要，包括地基、斜坡或岩石中的地下开挖。以往对含裂隙试样力学行为的研究主要集中在裂纹的扩展和聚合上。一些研究利用实验室测试研究了裂缝之间的裂隙聚合，此外，还提出了许多模拟裂纹产生和聚合的数值方法。然而，这些研究主要集中于岩石试件的裂纹扩展和聚合，而没有对工程项目的稳定性进行分析，因为这些结果很难用来预测原生含节理裂隙岩体的力学行为。大多数岩体都含有多组裂隙。如果一个岩体包括许多与结构相比尺寸相对较小的裂缝，那么可以将该岩体理想化为连续体。然而，当出现中等大小的不连续结构面时，这种影响对于岩体的力学行为变得很重要，因为这些不连续面十分复杂，并且不能以简单的方式估计其力学效应。这时，损伤力学是解决这一问题的有效途径。为此，本专利技术提供了一种多裂隙岩体试样力学性质的几何损伤流变分析方法。该方法针对含多组裂隙类岩石试样，综合考虑损伤和流变的影响，建立了节理岩体的损伤流变模型，损伤流变模型能够预测与初始蠕变特性相对应的粘弹性应变、与稳态蠕变特性相对应的粘塑性应变和反映裂隙对岩体损伤效应的损伤应变。并且推导出基于西原体模型的粘弹塑性损伤本构模型，并用有限差分软件FLAC3D进行了编程。目前国内相关多裂隙类岩石试样力学性质的研究现状如下：1、《节理岩体等效损伤流变模型初步研究》一文假设岩石为各向同性体各向同性损伤，节理面的法向和切向损伤不同，分别建立了岩石和节理面流变损伤演变函数。并假设材料无损的情况想，计算应力采用有效应力，推到了节理岩体等效损伤流变模型有限元计算公式，编制了相应的有限元程序(参见《岩土力学》2011年第12期，作者：黄耀英等)；2、《锦屏一级水电站地下厂房围岩变形破裂的三维损伤流变分析》一文采用损伤流变耦合的三维模型和分析方法，结合运用裂隙张开产生的附加变形分析法，对其工程地下厂房洞室群的稳定性进行分析，并对洞室的长期时效变形进行了预测分析(参见《岩石力学与工程学报》2012年第5期，作者：朱维申，齐银萍等)。
技术实现思路
要研究多裂隙类岩石试样的力学行为，其中一个比较困难的地方是将与时间有关的变形和岩体中裂隙分布引起的损伤影响联系起来。本专利技术采用损伤和流变耦合的方法，提供一种多裂隙类岩石试样力学性质的几何损伤流变分析方法。多裂隙类岩石试样力学性质的几何损伤流变分析方法，步骤如下：(1)、构建多裂隙类岩石试样代表性单元的几何模型；(2)、根据试样中分布裂隙的假设定义多裂隙类岩石试样代表性单元的损伤张量；(3)、计算多裂隙类岩石试样中单组裂隙的损伤张量；(4)、在确定单元的法向量和损伤张量后，通过对单元进行求和来计算多裂隙类岩石试样中多组裂隙的损伤张量；(5)、推导多裂隙类岩石试样有效应力公式；(6)、综合考虑损伤和流变的影响，建立多裂隙类岩石试样几何损伤流变本构模型；(7)、将多裂隙类岩石试样几何损伤流变本构模型程序化；(8)、统计工程岩体的几何裂隙参数，进行工程岩体的损伤流变特性数值模拟。相对于现有技术，本专利技术具有如下有益效果：1、综合考虑了流变现象和岩体中裂隙存在的影响，能够更好的分析多裂隙岩体的力学性质；2、此损伤流变模型能够预测与初始蠕变特性相对应的粘弹性应变、与稳态蠕变特性相对应的粘塑性应变和反映裂隙对岩体损伤效应的损伤应变，并可利用该模型预测地下洞室周围岩体的位移和损伤带；3、该专利技术方法可应用于地下洞室的开挖，分析其稳定性，具有较强的可行性，可用于指导实际的工程施工。附图说明图1为损伤流变模型的程序流程图；图2为岩体有效面积示意图；图3为岩体表面裂隙的观测示意图；图4为任意不连续面的单位法向量示意图；图5为裂隙数量统计示意图；图6为多裂隙岩体中地下硐室示意图。具体实施方式多裂隙类岩石试样力学性质的几何损伤流变分析方法，具体步骤如下：1、构建多裂隙类岩石试样代表性单元的几何模型，具体方法如下：本专利技术中提出了一个二阶损伤张量来表征具有若干平面裂隙岩体的损伤状态:其中Ω代表一组裂隙的面积密度，n是与节理垂直的单位向量，表示张量积。2、定义多裂隙岩石类试样代表性单元的损伤张量，具体方法如下：对岩体中分布的裂隙作如下假设：1)裂隙都是平面的；2)岩体由完整岩石的单元组成，这些单元的边界面上存在裂隙。裂隙将沿着表面扩展，不能穿透到单元中，而完整的岩石材料发生屈服；3)岩石块体尺寸由平均裂隙间距决定；4)只考虑初始损伤的影响，岩体截面有效面积的减小是导致岩体参数劣化的主要因素。如果υ和V分别是单元体和岩体的体积，则通过将单元替换为具有同体积的等效立方体，将总有效面积定义为:其中，l＝υ1/3。假设在岩体中存在N个裂隙，其中第k个裂隙具有面积ak和单位法线nk。对于此裂隙，裂隙的面积密度为：然后，该裂隙的损伤张量可以用以下形式定义：所有的N个裂隙通过方程(4)计算后进行求和，得到岩体V中的损伤张量为：3、计算多裂隙类岩石试样中单组裂隙的损伤张量，具体方法如下：如图4，若一组裂隙分别在坐标面X1，X2，X3上具有θ1，θ2，θ3大小的角度，则该组的单位法向量n＝(n1n2n3)：n＝(n1n2n3)t＝(λcosθ1cosθ2λsinθ1sinθ2-λcosθ1sinθ2)tλ＝(sin2θ1sin2θ2+cos2θ1sin2θ2+cos2θicos2θ2)-1/2(6)也可以表示为：n＝(n1n2n3)t＝(-λcosθ2sinθ3λcosθ2cosθ3λsinθ2sinθ3)tλ＝(sin2θ2sin2θ3+cos2θ2sin2θ3+cos2θ2cos2θ3)-1/2(7)或者n＝(n1n2n3)t＝(λsinθ3sinθ1-λcosθ3sinθ1λcosθ3cosθ1)tλ＝(sin2θ3sin2θ1+cos2θ3sin2θ1+cos2θ3cos2θi)-1/2(8)对于每个角度都有以下关系：其中(i，j，k)＝(1，2，3)，(2，3，1)，(3，1，2)。若立方体单元的X1、X2、X3表面分别包含N1、N2和N3个裂隙。假设其中一个表面包含最小数量的裂隙，如X3，然后我们旋转立方体单元，使得该组裂隙的单位法向量n与新轴X′3重合，如图5所示，也就是表面X′3包含最少数量的裂隙。沿X′1轴的裂隙数目估计为：<本文档来自技高网...

【技术保护点】
1.多裂隙类岩石试样力学性质的几何损伤流变分析方法，具体步骤如下：/n(1)、构建多裂隙类岩石试样代表性单元的几何模型；/n(2)、定义多裂隙类岩石试样代表性单元的损伤张量；/n(3)、计算多裂隙类岩石试样中单组裂隙的损伤张量；/n(4)、计算多裂隙类岩石试样中多组裂隙的损伤张量；/n(5)、推导多裂隙类岩石试样有效应力公式；/n(6)、建立多裂隙类岩石试样几何损伤流变本构模型；/n(7)、多裂隙类岩石试样几何损伤流变本构模型程序化；/n(8)、统计工程岩体的裂隙几何参数，进行工程岩体的损伤流变特性数值模拟。/n

【技术特征摘要】
1.多裂隙类岩石试样力学性质的几何损伤流变分析方法，具体步骤如下：
(1)、构建多裂隙类岩石试样代表性单元的几何模型；
(2)、定义多裂隙类岩石试样代表性单元的损伤张量；
(3)、计算多裂隙类岩石试样中单组裂隙的损伤张量；
(4)、计算多裂隙类岩石试样中多组裂隙的损伤张量；
(5)、推导多裂隙类岩石试样有效应力公式；
(6)、建立多裂隙类岩石试样几何损伤流变本构模型；
(7)、多裂隙类岩石试样几何损伤流变本构模型程序化；
(8)、统计工程岩体的裂隙几何参数，进行工程岩体的损伤流变特性数值模拟。


2.根据权利要求1所述的基于多裂隙岩体力学性质的几何损伤模型分析方法，其特征在于，构建多裂隙类岩石试样代表性单元的几何模型，具体方法如下：建立一个二阶损伤张量，即：



式中：Ω为一组裂隙的面积密度，n为裂隙的单位法向量，并且表示张量积。


3.根据权利要求1所述的基于多裂隙类岩石试样力学性质的几何损伤模型分析方法，其特征在于，定义多裂隙类岩石试样代表性单元的损伤张量，具体方法如下：假设v和V分别为单元体和岩体的体积，在V内有N条裂隙，则第k个裂隙面积为ak，它的单位法向量为nk，对于第k个裂隙，损伤张量可以用下面形式来定义：



式中l＝v1/3。
考虑所有的N条裂隙，可得多裂隙类岩石试样代表性单元的损伤张量为：





4.根据权利要求1所述的基于多裂隙类岩石试样力学性质的几何损伤模型分析方法，其特征在于，计算多裂隙类岩石试样中单组裂隙的损伤张量，具体方法如下：对于单组裂隙的多裂隙类岩石试样，其几何损伤张量为：



式中，Ni为Xi面上的裂隙数量，Nj为Xj面上的裂隙数量，Li为Xi面上的裂隙的平均长度，Lj为Xj面上的裂隙的平均长度。
任意裂隙面的单位法向量n＝(n1n2n3)为：
n＝(n1n2n3)t＝(λcosθ1cosθ2λsinθ1sinθ2-λcosθ1sinθ2)t
λ＝(sin2θ1sin2θ2+cos2θ1sin2θ2+cos2θ1cos2θ2)-1/2
同时，单位法向量还可以表达为下式：
n＝(n1n2n3)t＝(-λcosθ2sinθ3λcosθ2cosθ3λsinθ2sinθ3)t
λ＝(sin2θ2sin2θ3+cos2θ2sin2θ3+cos2θ2cos2θ3)...

【专利技术属性】
技术研发人员：杨文东，张建国，薄纯杰，王令，井文君，
申请(专利权)人：中国石油大学华东，山东正元建设工程有限责任公司，
类型：发明
国别省市：山东;37