CN113868897A - 一种热-力耦合条件下层状岩石统计损伤计算方法 - Google Patents
一种热-力耦合条件下层状岩石统计损伤计算方法 Download PDFInfo
- Publication number
- CN113868897A CN113868897A CN202111400440.2A CN202111400440A CN113868897A CN 113868897 A CN113868897 A CN 113868897A CN 202111400440 A CN202111400440 A CN 202111400440A CN 113868897 A CN113868897 A CN 113868897A
- Authority
- CN
- China
- Prior art keywords
- rock
- stress
- damage
- under
- formula
- 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.)
- Pending
Links
- 239000011435 rock Substances 0.000 title claims abstract description 173
- 230000006378 damage Effects 0.000 title claims abstract description 134
- 238000010168 coupling process Methods 0.000 title claims abstract description 57
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 57
- 230000008878 coupling Effects 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000009471 action Effects 0.000 claims abstract description 31
- 239000010454 slate Substances 0.000 claims abstract description 19
- 238000012360 testing method Methods 0.000 claims abstract description 18
- 238000005315 distribution function Methods 0.000 claims abstract description 4
- 238000012887 quadratic function Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 7
- 230000000977 initiatory effect Effects 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 230000003204 osmotic effect Effects 0.000 claims description 3
- 208000027418 Wounds and injury Diseases 0.000 claims description 2
- 238000012512 characterization method Methods 0.000 claims description 2
- 208000014674 injury Diseases 0.000 claims description 2
- 238000012669 compression test Methods 0.000 abstract description 3
- 230000005483 Hooke's law Effects 0.000 abstract 1
- 230000008569 process Effects 0.000 description 8
- 238000010438 heat treatment Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000005336 cracking Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005056 compaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000006735 deficit Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002277 temperature effect Effects 0.000 description 1
- 230000003685 thermal hair damage Effects 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/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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/08—Thermal analysis or thermal optimisation
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Computer Hardware Design (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种热‑力耦合条件下层状岩石统计损伤计算方法，首先给定岩石微元体强度服从Weibull分布函数，其次建立了岩石应力‑应变的轴向应力‑应变关系，最后根据广义胡克定律和统计损伤原理建立了考虑各向异性非线性强度破坏准则的岩石在三轴高温‑渗流耦合力试验条件下的应力‑应变关系。本发明应用于热‑力耦合作用下层状板岩的三轴压缩试验，提出的新模型较现有的模型综合性及适用性较强，对试验所得的应力应变曲线拟合程度较高。体现了层状板岩的变形破坏特征，充分反映了峰后阶段的趋势。
Description
技术领域
本发明涉及岩石统计损伤模型计算方法技术领域，尤其涉及一种热-力耦合条件下层状岩石统计损伤计算方法。
背景技术
随着水底隧道工程建设和运营项目的不断增加，其内部发生火灾的风险和概率也在不断攀升。火灾除了对人员造成巨大伤害，也会导致围岩支护系统的力学性能结构劣化。实际上，部分水底隧道具有穿越层状岩体等特点。若在这样的环境中发生火灾，往往会对隧道围岩变形产生更强烈的影响。加之，层状岩体中具有诸如层理、节理等原生或次生结构面，其变形和强度特征具有明显的各向异性，火灾高温作用下岩体的破坏机制及方式也明显不同于其他均质岩体。因此，研究层理岩石在高温-荷载耦合条件下的力学性能具有重要意义。
同时，目前关于损伤统计本构模型的研究相对较少，无法合理描述岩石损伤演化过程的缺陷，无法较好反映出高温作用岩石损伤的力学机制。
发明内容
本发明提供一种热-力耦合条件下层状岩石统计损伤计算方法，以克服现有的层状岩石统计损伤计算方法无法合理描述岩石损伤演化过程的缺陷。
为了实现上述目的，本发明的技术方案是：
一种热-力耦合条件下层状岩石统计损伤计算方法，所述岩石统计损伤计算方法包括以下步骤：
S1：定义损伤参数D(T,β)，所述损伤参数D(T,β)用于表征温度和层理角度对岩石受力性能的影响；
S2：定义用于表征载荷对岩石受力性能的影响的连续损伤变量D，建立在力的作用下岩石损伤的本构关系，以表征力对岩石受力性能的影响；
S3、引入修正系数，并得到温度和层理角度作用下岩石应力-应变的轴向应力-应变关系；
S4、建立热-力耦合条件下层状岩石统计损伤本构模型，以表征热-力耦合条件共同作用对岩石受力性能的影响，进而确定岩石在三轴高温-渗流耦合力试验条件下的应力-应变关系，所述应力-应变关系即岩石统计损伤关系。
进一步的，所述S1具体如下：
S11：定义用来表征在温度和层理角度E(T,β)作用下，岩石的损伤参数D(T,β)为：
式中：E(T,β)为温度T、层理角度β岩样的弹性模量；E(20℃,0°)为常温20℃、层理角度0°板岩的弹性模量；
S12：假设岩石材料微元体强度k服从Weibull分布，其概率密度函数f(k)可表示为：
式中：m代表不同温度下岩石的均匀性参数，K代表不同温度下岩石的强度特性参数，k代表岩石微元体强度；
由概率密度函数f(k)可得不同温度下岩石的强度特性参数K；
S13：表征温度和层理角度对岩石损伤参数D(T,β)的影响公式如下：
式中：m0为常温20℃、层理角度为0°板岩的岩石的均匀性参数，K0为常温20℃、层理角度为0°板岩的岩石的强度特性参数。
进一步的，所述S2具体如下：
S21、构建发生破裂的微元体数量模型NF，发生破裂的微元体数量模型NF为：
式中：N代表总体岩石微元体数量，f(x)代表Weibull的分布函数；
S22、定义用来表征在力的作用下，岩石的连续损伤变量D为：
S23、根据应力不变量表示的各向异性岩石的非线性强度准则f(σ)对微元体数量模型NF进行修正，得到：
式中：
分别为有效应力下的最大主应力和最小主应力；Bβ为各向异性系数；αi为强度参数；β为岩石试件的层理角度，σcβ为不同角度β板岩的单轴抗压强度；
S24、根据各向异性系数Bβ随单轴抗压强度σcβ的变化规律对微元体数量模型NF进行修正，得到：
式中：σtβ为不同角度β板岩的抗拉强度；
S25、给定各向异性岩石的经验非线性强度准则作为岩石微元体破坏参数，即：
f(σ)＝k(8)
S26、综合各向异性系数Bβ和各向异性岩石的非线性强度准则f(σ)构建连续损伤变量D，构建连续损伤变量D具体为：
进一步的，所述S2还包括：
S27、构建荷载统计损伤值D(σ,T,β)，构建荷载统计损伤值D(σ,T,β)具体为：
式中：D(σ,T,β)为考虑温度与各向异性效应板岩的荷载统计损伤值；σD为起裂应力。
进一步的，所述S3具体如下：
S31、引入损伤修正系数η，建立在力的作用下岩石损伤本构关系：
式中：σi为应力-渗流作用下的应力，
为应力-渗流作用下的有效应力；
S32、根据广义胡克定律，得出温度和层理角度作用下岩石应力-应变的轴向应力-应变关系εi，具体公式如下：
式中：
为三轴实验中径向方向上的第二方向有效应力和第三方向有效应力，其中第一方向有效应力为三轴实验中轴向方向上的有效应力，E为弹性模量，μ为泊松比；
由公式(11)和公式(12)推导得出岩石应力-应变的轴向应力-应变关系εi为：
进一步的，所述S4具体如下：
S41、将公式(9)代入公式(13)中，得到在考虑渗透压作用下轴向的应力-应变关系：
S42、利用应力转换方法得到第二方向有效应力
和第三方向有效应力
第二方向有效应力
和第三方向有效应力
具体公式为：
S43、结合损伤参数D(T,β)与第二方向有效应力
和第三方向有效应力
得到热-力耦合条件下层状岩石统计损伤本构模型σi，得到热-力耦合条件下层状岩石统计损伤本构模型σi具体为:
S44、结合岩石应力应变过坐标原点的特征，构建热-力耦合条件下层状岩石统计损伤本构模型对应二次函数，热-力耦合条件下层状岩石统计损伤本构模型对应二次函数具体为：
式中：A、B代表二次函数的拟合参数，εi代表应变值；
S45、结合公式(17)和公式(18)获得拟合后的热-力耦合条件下层状岩石统计损伤本构模型σi，拟合后的热-力耦合条件下层状岩石统计损伤本构模型σi具体为：
式中：εD为应力阀值点所对应的应变值。
进一步的，确定不同温度下岩石的均匀性参数m和不同温度下岩石的强度特性参数K的具体方法为：
根据岩石应力-应变的轴向应力-应变关系εi，结合边界条件σ＝σC,ε＝εC和σ＝σC,
可得：
当ε1>εD时，对公式(20)求偏导，并将峰值条件代入，可得：
联立式(20)和(21)即可求得不同温度下岩石的均匀性参数m和不同温度下岩石的强度特性参数K。
进一步的，确定二次函数的拟合参数A和B数值的具体方法为：
将公式(18)代入公式(20)可得：
当岩石的应力-应变曲线位于给定的损伤阀值区间前后，应力-应变曲线一阶导数连续，根据公式(22)又可表示其一阶导数为：
联立式(22)和(23)即可求得二次函数的拟合参数A和B数值。
进一步的，还包括利用连续损伤变量D获得高温-荷载耦合损伤Dc，获得高温-荷载耦合损伤D具体为：
Dc＝D(T,β)+D-D(T,β)D (24)。
有益效果：
1、考虑了温度和层理角度对岩石力学参数和本构模型参数的影响，本发明确立了本构模型的求解过程；
2、本发明所构建的热-力耦合条件下层状岩石统计损伤本构模型所计算的理论值与试验值相差不大，充分反映了峰后阶段的趋势，能够较好反映出不同温度作用和层理角度的板岩的应力-应变关系；
3、本发明进一步阐明了热损伤和荷载统计损伤可表征岩石的耦合损伤，并揭示了不同温度和层理角度的耦合损伤变量的演化规律。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明热-力耦合条件下层状岩石统计损伤计算方法流程图；
图2为高温-荷载耦合作用层理板岩三轴试验的流程图；
图3a为不同层理角度下弹性模量随温度变化图；
图3b为不同温度下弹性模量随层理角度变化图；
图3c为不同层理角度下峰值应力随温度变化图；
图3d为不同温度下峰值应力随层理角度变化图；
图4a为获得20℃、层理角度为45°条件下试验值与理论曲线对比图；
图4b为获得60℃、层理角度为45°条件下试验值与理论曲线对比图；
图4c为获得100℃、层理角度为45°条件下试验值与理论曲线对比图；
图4d为获得150℃、层理角度为45°条件下试验值与理论曲线对比图；
图5a为加热温度为100℃下层理角度为0°的试验值与理论曲线对比图；
图5b为加热温度为100℃下层理角度为30°的试验值与理论曲线对比图；
图5c为加热温度为100℃下层理角度为45°的试验值与理论曲线对比图；
图5d为加热温度为100℃下层理角度为60°的试验值与理论曲线对比图；
图5e为加热温度为100℃下层理角度为90°的试验值与理论曲线对比图；
图6a为本构模型中参数m、K随温度变化规律图；
图6b为本构模型中参数m、K随层理角度的变化规律图；
图7为耦合损伤变量DC的演化规律图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本实施例提供了一种热-力耦合条件下层状岩石统计损伤计算方法，如图1-7，包括以下步骤：
S1：定义损伤参数D(T,β)，所述损伤参数D(T,β)用于表征温度和层理角度对岩石受力性能的影响；
S2：定义用于表征载荷对岩石受力性能的影响的连续损伤变量D，建立在力的作用下岩石损伤的本构关系，以表征力对岩石受力性能的影响；
S3、引入修正系数，并得到温度和层理角度作用下岩石应力-应变的轴向应力-应变关系；
S4、建立热-力耦合条件下层状岩石统计损伤本构模型，以表征热-力耦合条件共同作用对岩石受力性能的影响，进而确定岩石在三轴高温-渗流耦合力试验条件下的应力-应变关系，所述应力-应变关系即岩石统计损伤关系。
在具体实施例中，所述S1具体如下：
S11：定义用来表征在温度和层理角度E(T,β)作用下，采用弹性模量作为损伤变量以表征温度和层理角度对岩石造成了损伤，定义岩石的损伤参数D(T,β)为：
式中：E(T,β)为温度T、层理角度β岩样的弹性模量；E(20℃,0°)为常温20℃、层理角度0°板岩的弹性模量；
在外荷载作用下，岩石的损伤是岩石微元体的不断破坏所导致的，假设在某一应力状态下发生破坏的微元体数量NF与总体微元体数量N的比值为统计损伤变量，即：
S12：由于Weibull概率分布容易积分、均值大于0等特点满足岩石受压破坏的统计特征，因此本文假设材料微元体强度k服从Weibull分布其概率密度函数f(k)可表示为：
式中：m代表不同温度下岩石的均匀性参数，K代表不同温度下岩石的强度特性参数，k代表岩石微元体强度；
由概率密度函数f(k)可得不同温度下岩石的强度特性参数K；
S13：假设m与K受温度和层理角度的损伤影响规律相似，表征温度和层理角度对岩石损伤参数D(T,β)的影响公式如下：
式中：m0为常温20℃、层理角度为0°板岩的岩石的均匀性参数，K0为常温20℃、层理角度为0°板岩的岩石的强度特性参数。
在具体实施例中，所述S2具体如下：
S21、岩石细观上是由随机分布的颗粒微元体结合而成，当微元体强度k超过某一强度时，岩石内部逐渐发生破裂最终形成宏观上的破坏，因此构建发生破裂的微元体数量模型NF，发生破裂的微元体数量模型NF为：
式中：N代表总体岩石微元体数量，f(x)代表Weibull的分布函数；
S22、定义用来表征在力的作用下，岩石的连续损伤变量D为：
S23、根据应力不变量表示的各向异性岩石的非线性强度准则f(σ)对微元体数量模型NF进行修正，得到：
式中：
分别为有效应力下的最大主应力和最小主应力；Bβ为各向异性系数；αi为强度参数；β为岩石试件的层理角度，σcβ为不同角度β板岩的单轴抗压强度；
S24、根据各向异性系数Bβ随单轴抗压强度σcβ的变化规律对微元体数量模型NF进行修正，得到：
式中：σtβ为不同角度β板岩的抗拉强度；
S25、给定各向异性岩石的经验非线性强度准则作为岩石微元体破坏参数，即：
f(σ)＝k(8)
S26、综合各向异性系数Bβ和各向异性岩石的非线性强度准则f(σ)构建连续损伤变量D，构建连续损伤变量D具体为：
在具体实施例中，所述S2还包括：
S27、岩石变形破坏过程实际上是其内部裂纹萌生、压密、拓展、汇合并贯通的过程，在应力值低于起裂应力σD的条件下，岩石内部裂隙不断扩展，岩石不发生损伤，因此在岩石变形过程中存在阀值点。并认为在岩石应力状态小于起裂应力时，无损伤演化或很小；当应力状态大于起裂应力时，其损伤可用式(9)表示，因此，在整个应力状态下，构建荷载统计损伤值D(σ,T,β)，构建荷载统计损伤值D(σ,T,β)具体为：
式中：D(σ,T,β)为考虑温度与各向异性效应板岩的荷载统计损伤值；σD为起裂应力。
在具体实施例中，所述S3具体如下：
S31、在三轴压缩试验的过程中，岩石试件受压头摩擦力及围压的影响，当其发生破坏后会传递剪应力和产生一定的残余强度，因此为了修正名义应力与有效应力的转化，本文引入了修正系数η，其中0<η<1。则三维各项同性应力的转化关系为引入损伤修正系数η，建立在力的作用下岩石损伤本构关系：
式中：σi为应力-渗流作用下的应力，
为应力-渗流作用下的有效应力；
S32、根据广义胡克定律，得出温度和层理角度作用下岩石应力-应变的轴向应力-应变关系εi，具体公式如下：
式中：
为三轴实验中径向方向上的第二方向有效应力和第三方向有效应力，其中第一方向有效应力为三轴实验中轴向方向上的有效应力，E为弹性模量，μ为泊松比；
由公式(11)和公式(12)推导得出岩石应力-应变的轴向应力-应变关系εi为：
在具体实施例中，所述S4具体如下：
S41、将公式(9)代入公式(13)中，得到在考虑渗透压作用下轴向的应力-应变关系：
S42、由于式(15)中的微元体强度包括有效应力，需要利用应力转换将其转化为名义应力，利用应力转换方法得到第二方向有效应力
和第三方向有效应力
第二方向有效应力
和第三方向有效应力
具体公式为：
S43、随着温度和角度的变化，板岩的弹性模量有显著的变化。结合损伤参数D(T,β)与第二方向有效应力
和第三方向有效应力
得到热-力耦合条件下层状岩石统计损伤本构模型σi，得到热-力耦合条件下层状岩石统计损伤本构模型σi具体为：
S44、损伤阀值点之后，热-力耦合作用下层状板岩损伤本构方程由式(18)给出。在损伤阀值点之前，结合岩石应力应变过坐标点(0,0)的特征，可采用二次函数进行拟合，则热-力耦合条件下层状岩石统计损伤本构模型对应二次函数具体为：
式中：A、B代表二次函数的拟合参数，εi代表应变值；
S45、结合公式(17)和公式(18)获得拟合后的热-力耦合条件下层状岩石统计损伤本构模型σi，拟合后的热-力耦合条件下层状岩石统计损伤本构模型σi具体为：
式中：εD为应力阀值点所对应的应变值。
根据试验结果可知岩石损伤阀值点处的应力即为岩石的起裂应力，起裂应力所对应的应变即为损伤阀值应变。为此，引入岩样的损伤阀值点处应力应变(εD,σD)与峰值应力处应力应变(εC,σC)作为宏观力学参数来确定高温三轴压缩试验模型参数。
在具体实施例中，确定不同温度下岩石的均匀性参数m和不同温度下岩石的强度特性参数K的具体方法为：
根据岩石应力-应变的轴向应力-应变关系εi，结合边界条件σ＝σC,ε＝εC和σ＝σC,
可得：
当ε1>εD时，对公式(20)求偏导，并将峰值条件代入，可得：
联立式(20)和(21)即可求得不同温度下岩石的均匀性参数m和不同温度下岩石的强度特性参数K。
在具体实施例中，岩石产生的变形破坏是一个连续变化的过程，且在应力应变曲线中的压缩阶段和损伤演化阶段的变化规律具有连续性。即在损伤阀值点前后的岩石变形具有一致性，在对本构模型中二次函数的拟合参数A和B数值的具体方法为：
将公式(18)代入公式(20)可得：
当岩石的应力-应变曲线位于给定的损伤阀值区间前后，应力-应变曲线一阶导数连续，具体的岩石的应力-应变曲线位于给定的损伤阀值区间前后各1％-5％，根据公式(22)又可表示其一阶导数为：
联立式(22)和(23)即可求得二次函数的拟合参数A和B的数值
在具体实施例中，还包括利用连续损伤变量D获得高温-荷载耦合损伤Dc，获得高温-荷载耦合损伤D具体为：
Dc＝D(T,β)+D-D(T,β)D (24)。
根据实验经验可得本构模型中参数m、K的取值范围，取值范围见表1。
表1
由图3a-d可知弹性模量、峰值强度与温度和层里角度具有相关性，而本发明建立的模型证明了模型中的参数mk均与温度和层理角度相关，表明模型中的参数mk的定义具有物理意义。由图4a-d和图5a-e可知试验值与理论曲线相差不大，表明本模型是较为合理的。图6a-b的描述是能够表现出m、K与温度和层理角度有着明显的相关性，图7表现出耦合损伤变量是随着轴向应力产生变化，是与试验值保持一致，说明本发明的模型的合理性。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (9)
1.一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，所述岩石统计损伤计算方法包括以下步骤：
S1：定义损伤参数D(T,β)，所述损伤参数D(T,β)用于表征温度和层理角度对岩石受力性能的影响；
S2：定义用于表征载荷对岩石受力性能的影响的连续损伤变量D，建立在力的作用下岩石损伤的本构关系，以表征力对岩石受力性能的影响；
S3、引入修正系数，并得到温度和层理角度作用下岩石应力-应变的轴向应力-应变关系；
S4、建立热-力耦合条件下层状岩石统计损伤本构模型，以表征热-力耦合条件共同作用对岩石受力性能的影响，进而确定岩石在三轴高温-渗流耦合力试验条件下的应力-应变关系，所述应力-应变关系即岩石统计损伤关系。
2.如权利要求1所述的一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，所述S1具体如下：
S11：定义用来表征在温度和层理角度E(T,β)作用下，岩石的损伤参数D(T,β)为：
式中：E(T,β)为温度T、层理角度β岩样的弹性模量；E(20℃,0°)为常温20℃、层理角度0°板岩的弹性模量；
S12：假设岩石材料微元体强度k服从Weibull分布，其概率密度函数f(k)可表示为：
式中：m代表不同温度下岩石的均匀性参数，K代表不同温度下岩石的强度特性参数，k代表岩石微元体强度；
由概率密度函数f(k)可得不同温度下岩石的强度特性参数K；
S13：表征温度和层理角度对岩石损伤参数D(T,β)的影响公式如下：
式中：m0为常温20℃、层理角度为0°板岩的岩石的均匀性参数，K0为常温20℃、层理角度为0°板岩的岩石的强度特性参数。
3.如权利要求2所述的一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，所述S2具体如下：
S21、构建发生破裂的微元体数量模型NF，发生破裂的微元体数量模型NF为：
式中：N代表总体岩石微元体数量，f(x)代表Weibull的分布函数；
S22、定义用来表征在力的作用下，岩石的连续损伤变量D为：
S23、根据应力不变量表示的各向异性岩石的非线性强度准则f(σ)对微元体数量模型NF进行修正，得到：
式中：
分别为有效应力下的最大主应力和最小主应力；Bβ为各向异性系数；αi为强度参数；β为岩石试件的层理角度，σcβ为不同角度β板岩的单轴抗压强度；
S24、根据各向异性系数Bβ随单轴抗压强度σcβ的变化规律对微元体数量模型NF进行修正，得到：
式中：σtβ为不同角度β板岩的抗拉强度；
S25、给定各向异性岩石的经验非线性强度准则作为岩石微元体破坏参数，即：
f(σ)＝k (8)
S26、综合各向异性系数Bβ和各向异性岩石的非线性强度准则f(σ)构建连续损伤变量D，构建连续损伤变量D具体为：
4.如权利要求3所述的一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，所述S2还包括：
S27、构建荷载统计损伤值D(σ,T,β)，构建荷载统计损伤值D(σ,T,β)具体为：
式中：D(σ,T,β)为考虑温度与各向异性效应板岩的荷载统计损伤值；σD为起裂应力。
5.如权利要求4所述的一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，所述S3具体如下：
S31、引入损伤修正系数η，建立在力的作用下岩石损伤本构关系：
式中：σi为应力-渗流作用下的应力，
为应力-渗流作用下的有效应力；
S32、根据广义胡克定律，得出温度和层理角度作用下岩石应力-应变的轴向应力-应变关系εi，具体公式如下：
式中：
为三轴实验中径向方向上的第二方向有效应力和第三方向有效应力，其中第一方向有效应力为三轴实验中轴向方向上的有效应力，E为弹性模量，μ为泊松比；
由公式(11)和公式(12)推导得出岩石应力-应变的轴向应力-应变关系εi为：
6.如权利要求5所述的一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，所述S4具体如下：
S41、将公式(9)代入公式(13)中，得到在考虑渗透压作用下轴向的应力-应变关系：
S42、利用应力转换方法得到第二方向有效应力
和第三方向有效应力
第二方向有效应力
和第三方向有效应力
具体公式为：
S43、结合损伤参数D(T,β)与第二方向有效应力
和第三方向有效应力
得到热-力耦合条件下层状岩石统计损伤本构模型σi，得到热-力耦合条件下层状岩石统计损伤本构模型σi具体为:
S44、结合岩石应力应变过坐标原点的特征，构建热-力耦合条件下层状岩石统计损伤本构模型对应二次函数，热-力耦合条件下层状岩石统计损伤本构模型对应二次函数具体为：
式中：A、B代表二次函数的拟合参数，εi代表应变值；
S45、结合公式(17)和公式(18)获得拟合后的热-力耦合条件下层状岩石统计损伤本构模型σi，拟合后的热-力耦合条件下层状岩石统计损伤本构模型σi具体为：
式中：εD为应力阀值点所对应的应变值。
7.如权利要求6所述的一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，确定不同温度下岩石的均匀性参数m和不同温度下岩石的强度特性参数K的具体方法为：
根据岩石应力-应变的轴向应力-应变关系εi，结合边界条件σ＝σC,ε＝εC和σ＝σC,
可得：
当ε1>εD时，对公式(20)求偏导，并将峰值条件代入，可得：
联立式(20)和(21)即可求得不同温度下岩石的均匀性参数m和不同温度下岩石的强度特性参数K。
8.如权利要求7所述的一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，确定二次函数的拟合参数A和B数值的具体方法为：
将公式(18)代入公式(20)可得：
当岩石的应力-应变曲线位于给定的损伤阀值区间前后，应力-应变曲线一阶导数连续，根据公式(22)又可表示其一阶导数为：
联立式(22)和(23)即可求得二次函数的拟合参数A和B数值。
9.如权利要求8所述的一种热-力耦合条件下层状岩石统计损伤计算方法，其特征在于，还包括利用连续损伤变量D获得高温-荷载耦合损伤Dc，获得高温-荷载耦合损伤D具体为：
Dc＝D(T,β)+D-D(T,β)D。 (24)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111400440.2A CN113868897A (zh) | 2021-11-19 | 2021-11-19 | 一种热-力耦合条件下层状岩石统计损伤计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111400440.2A CN113868897A (zh) | 2021-11-19 | 2021-11-19 | 一种热-力耦合条件下层状岩石统计损伤计算方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113868897A true CN113868897A (zh) | 2021-12-31 |
Family
ID=78985258
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111400440.2A Pending CN113868897A (zh) | 2021-11-19 | 2021-11-19 | 一种热-力耦合条件下层状岩石统计损伤计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113868897A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114923786A (zh) * | 2022-01-13 | 2022-08-19 | 大连海事大学 | 一种表达高温作用后岩石拉伸特性的损伤关系计算方法 |
| CN116611265A (zh) * | 2023-07-18 | 2023-08-18 | 北京建筑大学 | 一种深部各向异性岩石的应力和应变的预测方法及装置 |
| CN116738724A (zh) * | 2023-06-14 | 2023-09-12 | 中国国家铁路集团有限公司 | 一种围岩力学性质动态损伤本构模型的构建方法 |
| CN118010509A (zh) * | 2024-04-10 | 2024-05-10 | 中国矿业大学(北京) | 一种岩石脆性的评价方法和装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108952700A (zh) * | 2018-08-21 | 2018-12-07 | 西南石油大学 | 一种各向异性地层井壁破裂压力确定方法 |
| CN110135113A (zh) * | 2019-06-05 | 2019-08-16 | 中南大学 | 考虑尺寸效应的岩石结构面损伤统计本构模型的构建方法 |
| CN110687141A (zh) * | 2019-09-29 | 2020-01-14 | 重庆大学 | 一种钻井过程中储层页岩的渗透特性分析方法 |
| CN113155699A (zh) * | 2021-04-09 | 2021-07-23 | 大连海事大学 | 一种热-水-力共同作用的岩石统计损伤计算方法及其应用 |
| CN113392505A (zh) * | 2021-05-21 | 2021-09-14 | 大连海事大学 | 一种高温-水冷却循环作用后岩石统计损伤计算方法 |
| WO2021203491A1 (zh) * | 2020-04-07 | 2021-10-14 | 北京科技大学 | 一种采动围岩地压灾害智能预测方法及系统 |
-
2021
- 2021-11-19 CN CN202111400440.2A patent/CN113868897A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108952700A (zh) * | 2018-08-21 | 2018-12-07 | 西南石油大学 | 一种各向异性地层井壁破裂压力确定方法 |
| CN110135113A (zh) * | 2019-06-05 | 2019-08-16 | 中南大学 | 考虑尺寸效应的岩石结构面损伤统计本构模型的构建方法 |
| CN110687141A (zh) * | 2019-09-29 | 2020-01-14 | 重庆大学 | 一种钻井过程中储层页岩的渗透特性分析方法 |
| WO2021203491A1 (zh) * | 2020-04-07 | 2021-10-14 | 北京科技大学 | 一种采动围岩地压灾害智能预测方法及系统 |
| CN113155699A (zh) * | 2021-04-09 | 2021-07-23 | 大连海事大学 | 一种热-水-力共同作用的岩石统计损伤计算方法及其应用 |
| CN113392505A (zh) * | 2021-05-21 | 2021-09-14 | 大连海事大学 | 一种高温-水冷却循环作用后岩石统计损伤计算方法 |
Non-Patent Citations (5)
| Title |
|---|
| XUERU FAN，ETC: ""Dynamic Breakage Characteristics of Shale with Different Bedding Angles under the Different Ambient Temperatures"", 《ROCK MECHANICS AND ROCK ENGINEERING (2021)》, 13 April 2021 (2021-04-13) * |
| 刘港;马凤山;赵海军;冯雪磊;郭捷;: "考虑统计损伤模型的非均质节理岩体拉剪破坏机制研究", 岩土力学, no. 1, 10 July 2018 (2018-07-10) * |
| 廖安杰;孟陆波;李天斌;路强;成威;: "热力作用下层状砂岩各向异性三轴压缩试验研究", 岩石力学与工程学报, no. 1, 15 April 2019 (2019-04-15) * |
| 蒋浩鹏，等: ""基于 Weibull 分布的高温岩石统计损伤本构模型及其验证"", 《岩土力学》, vol. 42, no. 7, 31 July 2021 (2021-07-31) * |
| 郭永成;刘鑫宇;王克辉;胡鹏;: "温度-孔隙水压耦合作用下砂岩卸荷损伤本构模型研究", 水力发电, no. 05, 12 May 2020 (2020-05-12) * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114923786A (zh) * | 2022-01-13 | 2022-08-19 | 大连海事大学 | 一种表达高温作用后岩石拉伸特性的损伤关系计算方法 |
| CN116738724A (zh) * | 2023-06-14 | 2023-09-12 | 中国国家铁路集团有限公司 | 一种围岩力学性质动态损伤本构模型的构建方法 |
| CN116738724B (zh) * | 2023-06-14 | 2024-03-05 | 中国国家铁路集团有限公司 | 一种围岩力学性质动态损伤本构模型的构建方法 |
| CN116611265A (zh) * | 2023-07-18 | 2023-08-18 | 北京建筑大学 | 一种深部各向异性岩石的应力和应变的预测方法及装置 |
| CN116611265B (zh) * | 2023-07-18 | 2023-09-22 | 北京建筑大学 | 一种深部各向异性岩石的应力和应变的预测方法及装置 |
| CN118010509A (zh) * | 2024-04-10 | 2024-05-10 | 中国矿业大学(北京) | 一种岩石脆性的评价方法和装置 |
| CN118010509B (zh) * | 2024-04-10 | 2024-07-19 | 中国矿业大学(北京) | 一种岩石脆性的评价方法和装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113868897A (zh) | 一种热-力耦合条件下层状岩石统计损伤计算方法 | |
| Kodur et al. | Residual response of reinforced concrete columns exposed to design fires | |
| KR101497579B1 (ko) | 유한요소해에 기초한 취성재료의 압입파괴인성평가 방법 | |
| Bedon et al. | Numerical modelling of adhesive connections including cohesive damage | |
| Pineda | Collapse and upgrading mechanisms associated to the structural materials of a deteriorated masonry tower. Nonlinear assessment under different damage and loading levels | |
| Spiller et al. | Prediction of annealed glass window response to blast loading | |
| Ye et al. | Experimental study on cyclically-damaged steel-concrete composite joints subjected to fire | |
| Annecchiarico et al. | Micro and macro finite element modeling of brick masonry panels subject to lateral loadings | |
| CN117910269B (zh) | 压裂井水泥环径向开裂、裂纹扩展贯通风险评价调控方法 | |
| Agüera et al. | Structural response of unreinforced masonry walls | |
| CN117436264A (zh) | 一种考虑高温的节理岩石累计损伤本构模型建立的方法 | |
| Cvetkovska et al. | Fire scenario influence on fire resistance of reinforced concrete frame structure | |
| Al-Hamd et al. | Influence of loading ratio on flat slab connections at elevated temperature: A numerical study | |
| Zhong et al. | Seismic response of concrete gravity dam reinforced with FRP sheets on dam surface | |
| Kong et al. | Blast response of cracked reinforced concrete slabs repaired with CFRP composite patch | |
| CN116187050A (zh) | 一种高地应力围岩三轴条件下强度预测方法 | |
| Lan et al. | Progressive collapse resistance of RC frames subjected to localized fire | |
| Komurcu et al. | Numerical modelling of the in-plane loaded homogenized masonry walls | |
| Wichert et al. | Grouted segment joints for structures made of ultra-high performance concrete | |
| Gu et al. | Determination of mechanical parameters for elements in meso-mechanical models of concrete | |
| Ni et al. | Post-fire earthquake resistance of reinforced concrete structural walls | |
| Ahmed et al. | Finite element analysis for the response of URM walls supporting RC slab | |
| Mogili et al. | Cyclic vulnerability of RC knee connections under relatively high opening shear stresses | |
| Rodríguez et al. | High temperature effects in mortar and concrete specimens using a meso-mechanical model with fracture based zero-thickness interface elements | |
| Singamsetty et al. | Damage analysis of two-way reinforced concrete slab subjected to blast loads |
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 |