复杂性 复杂性 1099 - 0526 1076 - 2787 Hindawi 10.1155 / 2020/2706091 2706091 研究文章 损伤特征和深层岩石的本构模型下频繁影响干扰在卸货的过程中静态压力高 https://orcid.org/0000 - 0002 - 0858 - 654 x 1 2 3 4 https://orcid.org/0000 - 0002 - 4233 - 2893 张ydF4y2Ba Shuai-fei 1 https://orcid.org/0000 - 0002 - 8082 - 5955 美之 2 5 https://orcid.org/0000 - 0001 - 6626 - 7937 1 3 https://orcid.org/0000 - 0002 - 0164 - 2406 鲁平 1 3 6 https://orcid.org/0000 - 0003 - 4214 - 7599 Zu-qiang 1 Mobayen 萨利赫 1 能源科学与工程学院 河南理工大学 焦作 河南454000年 中国 hpu.edu.cn 2 广西防灾工程重点实验室安全 广西大学 南宁 广西530029年 中国 gxu.edu.cn 3 州和地方联合工程实验室深矿井的瓦斯抽放与地面控制 河南理工大学 焦作 河南454000年 中国 hpu.edu.cn 4 河南煤炭安全生产的协同创新中心 焦作 河南454000年 中国 5 土木工程和建筑学院 广西大学 南宁 广西530029年 中国 gxu.edu.cn 6 资源与安全工程学院 中南大学 长沙 湖南410083年 中国 csu.edu.cn 2020年 30. 5 2020年 2020年 18 12 2019年 14 03 2020年 27 03 2020年 30. 5 2020年 2020年 版权©2020年春王等。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

在矿产资源的开采深度,岩石在高地壳的复杂力学环境压力和爆破开挖卸荷,和改进的SHPB试验系统是用来进行频繁冲击扰动试验深岩石在卸货的过程中静态压力高。首先,动态应力-应变曲线的一般特征的整个过程趋势可分为四个阶段:稳定发展阶段的微裂纹,微裂纹的无机窝发展阶段,疲劳损伤阶段,疲劳破坏阶段。然后,岩石的损伤变量是由连续的因素,应变等效原理和统计损伤理论,这是基于整个测试期间的岩石的变形特征。和派生的损伤变量方程,岩石破坏过程中被证明是合理的,和频繁的干扰下的岩石损伤本构模型在卸货的过程中静态高压力是由组合模型建立的方法。最后,验证了本构方程的合理性测试数据,显示了动态应力-应变曲线包络良好一致性与本构方程的理论曲线。

广西防灾工程重点实验室安全 2019年zdk052 河南理工大学 672707年 SJF201802 中国国家自然科学基金 51904093 主要研发和推广的特殊项目(科技)河南省 192102310247 河南理工大学青年创新基础 NSFRF180321
1。介绍

随着经济的迅速发展,人类,各种资源消耗越来越多。而肤浅的资源耗尽,深层资源的开发是急需的,特别是深海矿产资源。然而,矿产资源的开采深度是不可避免的使深部岩体在高静态压力和频繁的复杂力学环境动态干扰( 1]。此外,静态力学理论不足以解释一些复杂的深层岩石变形特征和破坏现象,因此迫切需要探索适合解释力学特性的理论和方法。

为了研究深岩石的力学特性,岩石力学领域的学者进行了大量的实验,涉及损害的特点和本构模型在静态和动态负载下很深的岩石。对测试设备的研究和数据监测技术也重视。设备的持续改进和数据监测技术研究提供一个基本的保证深层岩石的力学特性。例如,Mobayen和Tchier 2- - - - - - 4)已经提出了新颖的鲁棒自适应二阶滑模跟踪控制技术和滑模扰动观测器控制技术,和一个非奇异的快速终端滑模稳定器。所有这些奠定了基础的动态负载测试技术的改进。

深岩石的损伤特性研究主要指损伤变量,破坏过程,损伤力学等等。从能源的角度来看,法律的能量耗散的过程中岩石变形和破坏被一些学者研究,可以有效地反映损伤过程( 5- - - - - - 7]。因此,岩石的方法来描述损伤变量的定义是根据能量耗散的律法,和损伤变量的理论公式推导出估算特定机械条件下岩石的损伤程度( 8- - - - - - 10]。基于岩石的力学参数,也定义为弹性模量损伤变量的方法,最大应变方法和剩余强度的方法,和损伤变量的相应公式推断,但没有一个能反映岩石的初始损伤程度下的频繁干扰高静态绝对压力( 11- - - - - - 14]。损伤机制,部分学者们还讨论了温度的影响,应力环境,湿度对岩石破裂特征。发现高温诱导萌发和扩张中的微裂缝岩石和岩石的强度减弱( 15, 16]。当压力环境不同,裂纹萌发或扩张的方向和形式在岩石内部是不同的,所以最终的失效模式也不同 17, 18]。此外,较高的岩石含水量更容易受到损害( 19, 20.]。

至于深岩石的本构模型,它包括影响因素的研究,建立模型,方程推导等等。因为岩石地质过程的产物,其特点影响力学性能是不连续的,各向异性不均匀等。例如,冲击载荷下岩石强度远远大于静态压力下( 21- - - - - - 23),高围压下岩石脆性和韧性变化( 24- - - - - - 27]。因此,它可以看到有很多复杂的因素影响岩石本构模型,和大量的研究已经由学者。考虑温度的交互作用,残余强度和力学条件下,形变场耦合损伤模型( 28, 29日)和热损伤模型( 30., 31日)建立,提供了理论参考深部岩体的开挖和钻井工程。由于含水量的影响和水的pH值在岩石力学性能,发现化学水造成的损害的程度在一定程度上反映了岩石内部损伤,所以损伤本构模型建立了裂隙岩体在水岩化学环境下( 32, 33]。学者们还讨论了高应力下岩石的本构模型和动态干扰影响,如非线性损伤高应力软岩的蠕变本构模型,建立了非线性蠕变破坏特征( 34, 35]。与高应力岩石从脆性过渡到韧性的特点,并建立了高应力下的损伤本构模型( 36]。考虑温度和动态耦合的双重作用,岩石力学行为分析研究岩石本构模型在这种条件下( 37, 38]。

总之,大量的研究已经进行了岩石损伤特性和本构模型在不同的力学环境,和获得的结果用于指导相应的工程实践,例如,刘和戴 39)提出了一种损伤本构模型描述裂隙岩体的变形和强度特征下循环单轴压缩。刘等人。 40)建立损伤本构模型描述循环荷载下岩石力学行为,基于能量耗散的法则。深的损伤特点及本构模型研究的作用下岩石频繁冲击负荷不足,特别是研究深岩石的损伤特点和本构模型下频繁影响干扰在卸货的过程中静态压力高。很深的岩石工程建设期间,矿体和围岩条件下的高地应力,尤其是采矿工程。深岩石的力学性能变化的影响下频繁扰动引起的爆破开挖的影响。同时,原岩应力的平衡状态是破碎爆破开挖过程中操作。因此,在新的应力平衡状态的形成,矿石和围岩在卸货的过程中。基于上述分析,环境的工程岩体可以概括如下:深层岩石遭受频繁的冲击干扰的过程中高压卸载。因此,试验研究深岩石的损伤特点和本构模型的基础上,进行了实际工程环境可以为岩石破坏提供理论参考和支持在深部岩体工程建设。

要解决这些问题,试验研究是进行一些新的想法。如高的轴向静载荷应用于岩石样本用于模拟的高地应力深岩,岩体工程的开挖卸载过程是由卸载模拟高静态负载在不同的利率,和负载的影响经常应用的岩石样本轴是用来模拟多个爆破造成的干扰影响。同时,生物人口增长模式的介绍,分析岩石的显微裂纹传播是第一次。通过这种方式,可以定量分析岩石的损伤特点,和深岩石的动态本构关系可以建立在相应的条件下。最后,指导深岩石工程建设的目的是实现。

2。高的静载荷卸载频繁冲击扰动试验 2.1。试样和设备 2.1.1。试样

测试的核心来自采矿巷道附近的围岩在冬瓜山铜矿的地下900米的深度,这是夕卡岩致密结构和良好的均匀的分析。确保non-parallelism non-perpendicularity标本小于0.02毫米,都精心打磨的两端。标本加工成两种尺寸(50毫米×100毫米,50毫米×50毫米),高径比2:1,它是用来测量单轴抗压强度。此外,高径比的标本1:1被用来研究高的频繁干扰测试静态应力卸载。

2.1.2。测试设备

单轴抗压试验和高静态应力卸载测试是在中南大学完成,和前研究了Instron1346型电液伺服试验机,而后者依靠改进的SHPB试验系统(结构图如图 1( 41为研究])。如图 1改进的SHPB试验系统轴向加载和卸载的能力。手动液压泵用于加载,加载时间可以用秒表记录控制加载速率。液压阀需要打开缓慢卸载压力时,此外,卸载值和时间记录控制卸载速度。

测试负载的结构图:(1)压力加载单元;(2)缓冲栏;(3)支持;(4)应变计A2;(5)传动杆;(6)岩石标本;(7)应变计A1;(8)事件栏;(9)异型冲击锤;(10)发射腔; (11) Gas gun.

2.2。测试计划

模拟高静态应力状态,基本原则是公元前由轴向压力设置尽可能接近单轴抗压强度,以不同的速率和轴向压力卸载模拟岩石开挖的过程。为了确保标本可以承受多重循环的影响,当轴向压力降低到50%公元前由轴向压力,高负荷的影响为0.5 MPa应用于标本模拟小扰动。在冲击荷载之前,岩石样品首先受到轴向静载荷。首先加载速率为0.5 MPa / s,然后更改为0.1 MPa / s。为了防止岩石样品损坏,因为加载速率过快,加载速率调整到0.1 MPa / s时,后来轴5 MPa的压力。同时,为了确保相同的冲击负载率测试,异型影响锤被放置在相同的位置在发射前腔的影响。当冲击压力是固定的,冲击载荷可以保持不变,这是应用于酒吧结束事件的影响。试验方案如表所示 1 2,平均单轴压应力的矽卡岩是126.63 MPa的数据表 1

测试方案和深矽卡岩单轴压缩下的结果。

数量 高度(毫米) 直径(毫米) 密度(克/厘米3 纵波速度(米/秒) 加载速度(毫米/秒) 峰值负载(kN) 单轴抗压强度(MPa)
XK1 100.02 50.06 3.78 5567年 0.03 239.08 121.53
XK2 99.36 49.86 4.26 5898年 0.03 384.10 196.82
XK3 99.24 49.56 4.06 5324年 0.03 255.80 132.67
XK4 98.96 49.98 3.84 5126年 0.03 176.93 90.23
XK5 100.02 50.02 3.72 5297年 0.03 252.89 128.76
XK6 99.04 49.94 3.94 5135年 0.03 220.17 112.46
XK7 100.08 50.04 3.90 5079年 0.03 204.31 103.94

测试计划的高静态应力卸载和频繁的干扰。

数量 高度(毫米) 直径(毫米) 密度(克/厘米3) 纵波速度(米/秒) Pre-axial压力(MPa) 轴向压力的影响(MPa) 卸货速率(MPa / s)的轴向压力 冲击压力(MPa)
XK1-1 50.34 53.60 3.10 4467年 65年 32.5 0.5 0.5
XK1-2 50.22 52.98 3.09 4549年 75年 37.5
XK1-3 49.85 52.88 3.33 4674年 85年 42.5
XK1-4 49.85 53.70 3.39 4496年 95年 47.5
XK2-1 49.47 53.08 3.56 4310年 65年 32.5 1。0 0.5
XK2-2 49.81 52.87 3.22 4625年 75年 37.5
XK2-3 50.05 53.72 3.53 4613年 85年 42.5
XK2-4 50.07 53.14 3.54 4428年 95年 47.5
XK3-1 49.74 53.15 3.19 4563年 65年 32.5 1。5 0.5
XK3-2 50.07 53.69 3.52 4146年 75年 37.5
XK3-3 49.41 53.20 3.30 4567年 85年 42.5
XK3-4 49.67 53.34 3.51 4366年 95年 47.5
XK4-1 50.30 53.36 2.94 4743年 65年 32.5 2.0 0.5
XK4-2 50.30 53.24 3.21 4571年 75年 37.5
XK4-3 49.71 53.26 2.85 4227年 85年 42.5
XK4-4 49.59 53.31 3.53 4252年 95年 47.5
2.3。测试结果

据帖子:轴向压力和卸荷速率的差异,分为16组试验。每个标本受到多种干扰影响,和相应的动态应力-应变曲线可以得到扰动后的影响。研究样本的总体变形特性,动态应力-应变曲线的包络线多个扰动影响下的每个标本被选中进行分析。在测试中,相应的动态应力-应变曲线可以得出基于从每个获得的测试数据的影响。可以选择一些特性数据,形成一组新的数据,如峰值前的峰值数据,所有数据在第一个获得的应力-应变曲线的影响,相应的数据后的残余强度阶段的峰值在最后的影响,和峰值相对应的动态应力-应变曲线中间数据的影响。应力-应变曲线的新数据组合是动态应力-应变曲线的包络曲线。图 2显示了一些代表性的动态应力-应变曲线,代表的数字干扰。图 3列出了信封的动态应力-应变曲线的卸货率0.5 MPa /年代和85 MPa的公元前由轴向压力,分别。在图 3(一个),数字代表公元前由轴向压力值,数字图 3 (b)显示卸货率的轴向压力。

矽卡岩的动态应力-应变曲线和包络曲线在高静态应力卸载的共同行动和影响干扰(轴向载荷95 MPa和卸载速率0.5 MPa / s)。

典型的动态应力-应变曲线的包络线。(一)轴向卸载率是0.5 MPa / s。(b)公元前由轴向压力是85 MPa。

从图可以看出 2,动态应力-应变曲线的初始阶段是近似线性,直线段与日益影响干扰时间逐渐缩短,和速度进入非线性阶段变得更快。似乎表明,试样弹性变形首先然后塑性变形在高静态应力卸载和频繁的干扰,并通过频繁的标本的弹性减弱干扰。因为样品的微裂缝关闭几乎在公元前由高轴向压力下,没有压实阶段,直接和标本进入弹性阶段。此外,走到单轴抗压强度极其公元前由轴向压力,导致损害标本,与此同时,损伤逐渐加重加载、卸载和频繁的干扰。因此,弹性变形阶段成为缩短和塑性变形阶段动态应力-应变曲线上逐步扩大。图 2还表明,包络的动态应力-应变曲线的变化趋势是一致的与动态应力-应变曲线在每个影响,这表明整个过程趋势的标本在同一轴向压力下,同样的卸货率,可以反映在相同的冲击强度包络在某些情况下的变化特征。

3(一个)显示相对应的峰值应力降低的信封公元前由轴向压力的增加,当卸货速率是恒定的。例如,当加载轴向压力65 MPa, 75 MPa, 85 MPa和95 MPa,相应的压力峰值应力的动态应力-应变曲线的包络线115.08 MPa, 110.58 MPa,分别为104.25 MPa和90.38 MPa。进一步分析后,相应的包络曲线的峰值应力降低比率是3.91%,5.50%和12.05%,分别时加载轴向压力从65 MPa提高到75 MPa,从75 MPa - 85 MPa,从85 MPa 95 MPa。因此,内部损伤加重公元前由高轴向压力,即轴向压力越高,较弱的岩石抵抗外部干扰。一句话,影响标本的干扰可以承受的数量减少。

3 (b)显示信封的变化趋势基本保持不变,当轴向压力是一样的,它穿过一条直线段,然后进入非线性发展阶段,表明岩石的动态变化趋势不受卸荷速率的影响。卸货速率为0.5 MPa / s时,1.0 MPa /秒,1.5 MPa / s和2.0 MPa / s,相应的动态应力峰值为104.25 MPa, 103.85 MPa,分别为103.38 MPa和103.15 MPa。动态应力峰值之间的范围从0.2%变化到0.5%,显示一个小递减趋势,可以得出的结论是,允许动态应力随卸率上升而减小。当轴向压力是肯定的,卸货速度越大,所需的时间越短卸相同的轴向压力。标本的时候受到的卸载过程高轴压短之前每个影响。同时,微裂缝的时间在标本发芽和扩大将短,导致一个更小的内部损伤积累。此外,当试样轴向压力和卸货的过程中,它有延迟的影响力量的外部影响。当卸货率较高,延迟的影响力量更加明显,导致有效的冲击应力的下降和内部疲劳损伤标本。卸货率越高,越有利于稳定的岩石。

信封的一般特征的动态应力-应变曲线可以分为四个阶段,如稳定微裂缝的发展,微裂缝的不稳定发展,疲劳损伤,疲劳失效,如图 4

典型的动态应力-应变曲线的包络线。

在微裂隙的稳定发展阶段(OA部分),一条直线的曲线发展和相应的变形模量是最大的,表明样品在弹性变形阶段的开始。在这个阶段,试样的微裂隙或稳定不变,这反映了标本阻力影响最大,这是更加难以产生大的变形。

在不稳定传播阶段的微裂隙(AB),曲线非线性的方式发展和减缓;与此同时,相应的变形模量逐渐减少,表明试样在塑性变形阶段。标本显示瞬时传播的内部微裂缝和贯穿在频繁的干扰,导致宏观失败。然而,每个动态扰动很小,所以一个影响未能使微裂纹在试样立即通过,但使它在一个不稳定的发展。

在疲劳损伤阶段(BC段),这部分的曲线显示了一个下降的趋势,下降趋势是相对较慢。在这个阶段,有大量的岩石压缩变形,但减少动态应力相对较小,显示了一定的延性。因为弹性能量存储在标本是相对较大的,它能够抵消干扰频繁的一部分。然而,试样的损伤程度与增加影响数量仍在加剧,以及弹性能量存储在标本逐渐减少,这显示了一个缓慢下降的趋势包络曲线。

在疲劳失效阶段(CD段),试样的脆性增强。一般来说,没有宏观冲击应力下的标本的失败,导致宏观失败,主要因素是公元前由高轴向压力。

3所示。损伤特征的分析 3.1。损伤变量

在这个测试中,样品内部的损伤加剧相互影响,因此,合理的定义损伤变量有助于分析和反映损伤特征在整个测试过程。由于复杂的内部结构的岩石,岩石微量元素被认为是连续和随机,和定义的损伤变量是连续的因素,应变等效原理,或者统计损伤理论。与此同时,他们也使用加权平均进行分析。

3.1.1。损伤变量D <斜体> < /斜体> <子> 1 < /订阅>基于连续系数和应变等效原理

有结构紧凑和试样的均匀性好,而且它假定isotropous在测试。基于应变等效原理,认为名义应力造成的变形作用于材料等于造成的损害有效应力作用于虚拟无损材料。应变等效原理的示意图岩石微量元素如图 5

应变等效原理的示意图岩石微量元素。

总之,损害体积总量的比率的岩石微量元素定义损伤变量 D1,公式如下: (1) D = 1 年代 ˜ 年代 , 0 D 1。

结合连续因素,应变等效原理和胡克定律,应力和应变之间的关系可以从图中获得 5: (2) ε = σ E ˜ = σ ˜ E

公式( 3从公式计算)( 1)和( 2): (3) D 1 = 1 年代 ˜ 年代 = 1 σ σ ˜ = 1 E ˜ E

3.1.2。损伤变量D <斜体> < /斜体> <子> 2 < /订阅>基于统计损伤理论

假设的缺陷岩石微量元素是独立的,随机的,他们符合威布尔分布的分布。基于统计损伤理论,唐 42]。计算岩石的损伤变量,公式如下: (4) D 2 = 1 ε α + 1 经验值 ε α , ε 0。

在公式( 4), ε, α,质量代表了岩石损伤的应变,岩石完整性和岩石性质和形状的有关参数,分别。

3.1.3。损伤变量D <斜体> < /斜体>在高静态应力卸载和频繁的干扰

根据测试结果,损伤变量 D是理性的比例 D1 D2是50%,所以呢 D被定义为它们之间的加权平均: (5) D = D 1 + D 2 2 = 1 2 1 E ˜ E + 1 ε α + 1 经验值 ε α , ε 0。

3.2。损伤演化方程 3.2.1之上。方程分析损伤变量D <斜体> < /斜体> <子> 1 < /订阅>

随着连续应用因子和等效应变损伤变量的定义 D 1,必须满足以下假设。首先,岩石材料损伤和无损的微量元素组成。其次,它们的体积等于,无损的微量元素能够不可逆转地微量元素瞬间转变成伤害。最后,破坏只发生而不是横向沿轴向方向,和无损的微量元素服从胡克定律。

结合方程( 1)和( 3),损伤变量的方程可以表示为: (6) D 1 = N 年代 0 年代 0 = N

在公式( 6), N受损的微量元素的总数,是无损的数量的总和受损微量元素和微量元素,也就是说,在样品微量元素的总数,然后呢 年代0是单个元素的体积。

假设生物种群的生长过程是用来模拟岩石的损伤演化过程,种群的生存环境被认为是摇滚外部负载。此外,在种群增长模型对应于岩石的应变在外部负载,以及种群的环境容量对应于岩石微量元素的总数,所以损坏的微量元素的数量的增长速度可以表示如下: (7) d N d ε = r N 1 N

在公式( 7), γ ε自然增长率和有关的岩石,分别。

然后,公式( 7)是解决分离变量的微分方程,表达式( 8)得到: (8) N = 1 + C e r ε , C = N 0 N 0

在公式( 8), N0微量元素是初始数量的岩石损伤。

微分公式( 9)的损伤变量 D1推导出公式( 6)和公式( 7): (9) d D 1 d ε = 1 d N d ε = r N 1 N

代入公式( 8)到公式( 7)和公式( 10)的损伤变量是通过分离变量的方法解决: (10) D 1 = 1 1 + e β r ε , β = ln N 0 1

在公式( 10), β代表岩石材料的初始损伤程度。

3.2.2。方程分析损伤变量D <斜体> < /斜体> <子> 2 < /订阅>

根据统计损伤理论,岩石中微量元素的分布缺陷是独立的,随机的,它满足泊松分布。因此,缺陷发生的概率可以表示通过泊松法律长度的区间内 l和概率函数的表达公式所示( 11)如下: (11) P k l = λ l k k ! e λ l

在公式( 11), λ微量元素的数学期望是缺陷单位长度, k微量元素缺陷的数量, P( k / l)的概率函数 k缺陷出现在区间的长度 l

胞体缺陷的概率假设Δ的长度 l P1 l元素的概率),没有身体缺陷的长度 l P( l)。公式( 12)和( 13)可以推断如下: (12) P 1 Δ l = λ Δ l e λ Δ l , (13) P l + Δ l = P l P l P Δ l

公式( 14)可以从公式推导( 12)和( 13): (14) P l + Δ l P l Δ l = λ P l e λ Δ l

l= 0, P( l)= 1,Δ的价值 l趋向于零,所以双方的极限公式( 14)可以得到: (15) P l = e λ l

因此,概率函数 ϕ( l)的缺陷不少于一个微量元素区间的长度 l可以推导出: (16) ϕ l = 1 P l

的概率密度函数 ϕ( l微量元素可以从公式推导)( 16): (17) φ l = d ϕ d l = λ e λ l

假设破坏样品中微量元素失去承载力和岩石应变是无限的受损时,损伤变量的表达 D 2表示如下: (18) D 2 = 0 ε l φ l d l 0 l φ l d l = 1 λ ε + 1 e λ ε

虽然考虑到线的缺陷,被认为是平面和身体,完整性 α = 1 / λ ( 43]。因此,公式( 18)可以写成公式( 19)。 (19) D 2 = 1 ε α + 1 经验值 ε α

在公式( 19), ε是岩石损伤的应变,然后呢表示参数相关的材料和形状的岩石。

3.2.3。定义损伤变量D <斜体> < /斜体>

基于上述方程分析,损伤变量 D岩石在高静态应力卸载和频繁的干扰可以被定义为公式( 20.): (20) D = 1 2 1 + 1 1 + e β r ε ε α + 1 经验值 ε α , β = ln N 0 1 , ε 0

3.2.4。确定损伤演化方程的参数

如果损伤变量 D的参数进行分析和计算, β, r, α首先必须确定。确定的值 β r的方程 D1需要分析的参数 α由方程的分析吗 D2

(1)确定参数β和r。首先,损伤变量的值 D 1对应于每个应变的包络计算公式( 3),然后 E(非破坏性材料的弹性模量)和 E ˜ (受损材料的弹性模量)公式( 3由图)计算 6

的原理图 E E ˜ 与确定方法。

公式( 10)是推导,公式( 21)要求: (21) Y = ln 1 D 1 1 = β r ε

计算 Y值对应于每个应变的动态应力-应变信封,然后分析拟合线性公式( Y- - - - - - ε线性公式)。最后,通过拟合线性公式( 21)与物流模型,和价值观 β r能够被确定。

(2)确定参数m和α。统计损伤理论认为 α(岩石完整性)的函数 ε( 42),而 α成反比 ε在高静态应力卸载和频繁的干扰如下: (22) α = C ε , C 是一个比例常数

假设方程 ε/ C= 1 / α1持有, α1是一个参数与岩石性质和形状,并代入公式( 22),这是发现的值 α α1是相等的。因此,公式( 19)可以改写公式( 23)如下: (23) D 2 = 1 ε α 1 + 1 经验值 ε α 1

据推测受损的微量元素失去承载力和non-damaged微量元素符合胡克定律,因此,岩石应力、应变之间的关系在高静态应力卸载和频繁的干扰所示公式( 24)如下: (24) σ = E ε 1 D 2

代入公式( 23)到公式( 24)和公式( 25)给出如下: (25) σ = E ε ε α 1 + 1 经验值 ε α 1

此外,如果 σ = σ / E α 1 ε = ε / α 1 是建立、公式( 25)可以转化为无量纲的形式,如公式( 26): (26) σ = ε ε + 1 经验值 ε

根据公式( 26),岩石的动态应力-应变曲线理论不同值画在图 7

动态应力-应变曲线的原理图与一维无量纲。

通过分析曲线在图的相似性 7和测量动态应力-应变曲线的包络线的价值可以被估计。后的值确定的最大价值 σ 通过公式( 26),参数 α计算公式( 27),由唐(推导出 42]。 (27) α = α 1 = P E σ 马克斯

在公式( 27), P峰值应力的动态应力-应变曲线的包络线。

3.3。损伤演化规律

为了验证损伤变量的定义 D在公式( 5)的值 β, r, α代入这个公式,所以相应的动态损伤变量 D通过结合相应的应变计算动态应力-应变曲线的包络线。通过比较样本的损伤程度和分析动态variable-strain理论曲线的特点,得出这样的结论:定义的损伤变量是合理的。图 8列出了两种典型破坏variable-strain曲线的理论关系,和图 9显示相应的损伤状态两组岩石在冲击过程中。

岩石的破坏variable-strain曲线在高静态应力卸载和频繁的干扰。(一)轴向压力是85 MPa。(b)卸货速率为0.5 MPa / s。

在冲击过程中岩石的损伤状态。(一)卸货率为0.5 MPa / s,这是第七的影响。(b)加载轴向压力是85 MPa和第十的影响。

从图可以得出结论 8

动态应变的值是0时,损伤变量的值 D大于0,所以定义损伤反映初始损伤程度的岩石在高静态压力。

存在四个阶段安装损伤variable-strain曲线,适度增长,快速增加,慢慢增加,趋于稳定。整个曲线呈现“s”,损伤变量的值 D是在0和1之间。它可以有效地对应于四个发展阶段的岩石在测试,稳定发展的微裂缝,迅速传播的微裂缝、疲劳损伤和疲劳失效。

当卸货速率是恒定的,轴向压力越大,损伤变量的值越快 D趋向于1,反映pre-high轴向压力加速了岩石破坏的过程和失败。

当预加载轴向压力是恒定的,卸货率越小,损伤变量越快 D往往是1,这表明高卸率有利于改善岩石承受频繁的干扰的能力。

如图 9(一个),当卸载速率和影响时间维持不变,公元前由轴向压力越大,岩石损伤越严重。当公元前由轴向压力和影响时间保持不变,卸货率越小的轴向压力的影响过程中,岩石损伤越严重,如图 9 (b)。上述两种现象表明,帖子:高轴向压力削弱了岩石抵抗外部冲击载荷的能力,和高卸率提高岩石承受频繁干扰的能力,证明定义损伤变量的合理性。

4所示。本构模型 4.1。基本假定

基于分析结果的信封的岩石的动态应力-应变曲线在该测试中,在冲击过程中损伤演化规律相结合建立本构模型在高静态应力卸载和频繁的干扰。然而,本构模型的建立需要满足某些假设,如下:

本构关系不受惯性效应的影响在恒定速度训练( 44]。

因为摇滚元素的特征是弹粘性和统计损伤,它的目的是作为一个麦克斯韦体(由弹性元件和粘性元素)在与伤害身体 D 一个1首先,然后串联伤害身体 D 一个2 岩石的力学模型单元组合如图 10

根据应力状态,本构关系 σ ε表示如下( 45]:

(28) σ = E ε 1 D ε 0

在公式( 28), σ ε分别是伤害身体的应力和应变。

由于无损特性,粘性元素的本构关系表示如下( 46]:

(29) σ 2 = η d ε 2 d t ,

在公式( 29日), σ2 ε2分别是粘性的应力应变分量, η粘度系数, t是相对应的时间压力。

岩石的应力-应变关系元素伤害符合线性微分方程之前,这大约是考虑应变叠加的原则仍然有效( 47]。

尽管粘弹性本构方程和应变等效的原则,破坏了岩体的影响,在高静态应力卸载和频繁的干扰,本构关系可以预测 48]。

岩石的力学模型单元组合。

4.2。建立本构模型

岩石力学模型的元素,应力和应变之间的关系的粘弹性的身体,伤害身体如下: (30) σ = σ 一个 1 + σ 2 = σ 一个 1 + σ 1 = σ 一个 2 , ε = ε 一个 1 + ε 一个 2 , ε 一个 1 = ε 1 + ε 2 , σ 1 = σ 2

在公式( 30.), σ ε分别是应力、应变的组合模型, σ 一个1 ε 一个1的应力和应变是伤害身体 D 一个1分别 σ 一个2 ε 一个2的应力和应变是伤害身体 D 一个2分别 σ1 ε1应力和应变的弹性组件,分别 σ2 ε2分别是粘性的应力应变分量。

胡克身体的本构方程,粘性的身体和受损的身体被替换成公式( 30.),所以力学模型的本构方程推导出岩石组合: (31) η E 一个 2 1 D + E 一个 1 1 D + E 1 σ ˙ + E 1 E 一个 2 1 D + E 1 E 一个 1 1 D σ = η E 1 E 一个 2 1 D + E 一个 1 E 一个 2 1 D 2 ε ˙ + E 1 E 一个 1 E 一个 2 1 D 2 ε

在公式( 31日), E1, E 一个1 E 一个2是弹性元件的弹性模量,伤害身体 D 一个1和伤害身体 D 一个2分别。

的弹性模量 E损害之前应该被有效弹性模量 E(1− D)由于应变等价原理 49),但不考虑初始损伤特点首先公式(更好的解决方案的 31日),即弹性模量 E 一个1 E 一个2采用替换 E 一个1(1− D), E 一个2(1− D在公式() 31日),分别。因此,获得组合模型的本构方程: (32) η E 一个 2 + E 一个 1 + E 1 σ ˙ + E 1 E 一个 2 + E 1 E 一个 1 σ = η E 1 E 一个 2 + E 一个 1 E 一个 2 ε ˙ + E 1 E 一个 1 E 一个 2 ε

因为样品受到高静压,如果 t= 0,初始条件 ε(0)= ε0 σ(0)= F建立了。因此,公式( 33通过拉普拉斯变换得到: (33) σ t + t 0 = E 一个 2 E 1 + E 一个 1 E 1 + E 一个 1 + E 一个 2 ε t + t 0 E 1 2 E 一个 2 2 η E 1 + E 一个 1 + E 一个 2 2 0 t ε τ + t 0 e E 1 E 一个 1 + E 1 E 一个 2 / η E 1 + E 一个 1 + E 一个 2 t + t 0 τ d τ

在公式( 33), t0静态压力的加载时间,什么时候 t0= 0, ε( t0= 0)= 0 σ( t0= 0)= 0是正确的。

如果 ε( t + t0)= ε0+ ε r( t0)= ε0+ ct,在那里 c是恒定应变率不变,公式( 33)能够转化为公式( 34): (34) σ t + t 0 = E 一个 2 E 1 + E 一个 1 E 1 + E 一个 1 + E 一个 2 ε 0 + ε r t ε 0 E 1 E 一个 2 2 E 一个 1 + E 一个 2 E 1 + E 一个 1 + E 一个 2 e 一个 t 0 e 一个 ε r t / c + t 0 + η c E 一个 2 2 E 一个 1 + E 一个 2 2 e 一个 t 0 e 一个 ε r t / c + t 0 E 1 E 一个 1 + E 1 E 一个 2 ε r t η c E 1 + E 一个 1 + E 一个 2 e 一个 t 0 , 地点: (35) 一个 = E 1 E 一个 1 + E 一个 2 η E 1 + E 一个 1 + E 一个 2

在公式( 34), ε r( t)下的应变扰动静应力之间的相互作用和影响, ε0静态应力产生的压力吗 F

由于应变等价原理( 43), E 一个1 E 一个2在公式( 34)所取代 E 一个1(1− D), E 一个2(1− D),本构方程计算如下: (36) σ t + t 0 = E 一个 2 1 D E 1 + E 一个 1 1 D E 1 + E 一个 1 1 D + E 一个 2 1 D ε 0 + ε r t ε 0 E 1 E 一个 2 2 1 D 2 1 D E 一个 1 + E 一个 2 E 1 + E 一个 1 1 D + E 一个 2 1 D e k t 0 e k ε r t / c + t 0 + η c E 一个 2 2 1 D 2 1 D 2 E 一个 1 + E 一个 2 2 e k t 0 e k ε r t / c + t 0 1 D E 1 E 一个 1 + E 1 E 一个 2 ε r t η c E 1 + E 一个 1 1 D + E 一个 2 1 D e k t 0 , 地点: (37) D = D t + t 0 = 1 ε 0 + ε r t α + 1 经验值 ε 0 + ε r t α , k = E 1 1 D E 一个 1 + E 一个 2 η E 1 + E 一个 1 1 D + E 一个 2 1 D

4.3。测试结果的验证

计算岩石的本构方程在高静态应力卸载和频繁的干扰,的值 E1, E 一个1, E 一个2,, α, η,β, r需要由分析实测数据和计算进行审判。本构方程的压力 ε r( t),静态加载时间的测量值 t0和恒应变速率 c测试数据,在哪里 ε r( t)是对应的应变动态应力-应变曲线的包络线,和 t0之前被公元前由静态负载的平均时间多重影响。

根据测试数据和岩石的本构方程,得到以下结果:的价值 E1类似于变形模量下的应力-应变曲线的初始阶段首先影响在实验过程中,前者可以取代后者。的价值 E 一个1能够表达的动态变形模量 E d, E d定义为割线模量的加权平均,II型割线模量和变形模量加载部分的信封,用于反映了在加载阶段的岩石压缩变形特征。此外,有一个比例关系 E 一个2 E d,这是 h= E 一个2/ E d,在那里 h是动态扩展因数。在测试,冲击载荷的作用时间短,岩石的动态和静态相结合合力增加迅速,和相应的应变降低或保持不变,表达了一种弹性现象的动态应力-应变曲线,因此,计算瞬时动态变形模量突然变化。因此,比率 h( h≥1)的瞬态动态变形模量 E 一个2和动态变形模量 E d突变后被定义为动态膨胀系数和的值 h可以通过测量数据计算。的四个参数 β, r, α相关参数的损伤演化方程,他粘滞系数吗 η根据第四个基本假设,推导出和它的值范围是500∼1000 GPa·s。

通过分析测试数据,两组数据的动态应力-应变曲线是典型的信封选择试验计算,如表所示 3

实验参数和理论拟合参数。

卸货速率(MPa·s−1) F(MPa) E1(GPa) E 一个1(GPa) E 一个2(GPa) h(GPa·s) α(10−3) c(10−3年代−1) β r t0(年代) 动态膨胀系数( h)
0.5 65年 184.66 135.11 358.04 800年 1.21 0.91 20.05 3.05 2216年 130年 2.65
0.5 75年 208.03 165.89 432.96 800年 2.08 0.83 23.39 2.87 2122年 150年 2.61
0.5 85年 148.21 124.60 505.88 800年 1.48 1.07 20.67 1.18 1515年 170年 4.06
0.5 95年 217.70 150.52 514.77 800年 0.82 0.50 19.85 1.67 1865年 190年 3.42
1 85年 144.50 125.13 472.97 800年 1.54 1.11 22.54 1.87 1835年 170年 3.78
1。5 85年 151.56 119.56 280.98 800年 1.38 1.05 21.67 2.10 1657年 170年 2.35
2 85年 225.99 152.87 502.96 800年 1.23 0.68 18.98 1.13 1221年 170年 3.29

参数表 3被替换成岩石本构( 36),以及相应的理论动态应力-应变曲线拟合和比较与实验信封进行分析,如图 11 12

理论和实验的比较图包络曲线轴向卸载下的应力-应变率为0.5 MPa / s。(一)预加载轴向压力65 MPa。(b)预加载轴向压力75 MPa。(c)预加载轴向压力85 MPa。(d)预加载轴向压力95 MPa。

理论和实验的比较图的包络曲线应力-应变的公元前由轴向压缩下85 MPa (a)轴向压力卸荷速率为0.5 MPa / s。(b)轴向压力卸荷速率为1.0 MPa / s。(c)轴向压力卸荷速率为1.5 MPa / s,和(d)轴向压力卸荷速率为2.0 MPa / s。

如数据所示 11 12,当轴向卸载或公元前由轴向压缩率保持不变,有一个很好的一致性理论动态应力-应变曲线和实验动态应力-应变包络曲线。因此,本构方程准确预测岩石的动态应力和应变之间的关系在整个生产过程中在高静态应力卸载和频繁的干扰,并间接反映了岩石的动态变形特性。此外,产量和最大应力的岩石本构方程能够有效预测,以及岩石变形对应的屈服应力,这些灾害预防提供理论依据的深层岩石开挖工程。

5。结论

深部岩体在高静态应力下降负载的机械环境和频繁的在开挖扰动,影响和损伤特点及本构模型研究的很深的岩石进行,主要结论如下:

信封建立了动态应力-应变曲线可以反映岩石的整体变化趋势下频繁扰动在卸货的过程中静态应力高,以及稳定发展阶段,无机窝扩张,岩石的疲劳损伤和疲劳失败。

基于连续因素,应变等效原理和统计损伤理论,频繁的干扰下的岩石损伤变量在卸货的过程中定义静态应力高,和损伤变量的方法来确定参数推导出。与此同时,损伤变量定义是由实验数据证明是合理的和失败的历史岩石。

根据信封的动态应力-应变曲线的变化特征,岩石本构模型下的频繁干扰在卸货的过程中建立高静态应力与某些假设,推导出了相应的本构方程。最后,实验曲线和理论曲线的包络线进行比较,并发现他们有良好的一致性。

5.1。不足和建议

的过程中研究深岩石的损伤特点和本构模型下的频繁影响干扰在卸货的过程中静态应力高,发现有一些缺点,如卸货率控制是困难的,岩性很简单,本构方程涉及到很多参数,研究结果并没有被深岩石工程验证。因此,以下几方面需要进一步改进或在未来深化。

提高卸货率控制设备和减少卸货率控制困难。同时,深岩石的损伤特点和本构模型需要进一步讨论通过扩大岩石类型和增加高温等因素的影响。

为了提高本构方程的应用,有必要进行进一步的研究,优化各种参数之间的关系,并简化本构方程的表达式。

基于特定的力学特性的矿体和围岩深部岩体工程、破坏特征和本构模型之间的定量关系深,深的岩石应进一步探讨。

数据可用性

文章中的所有相关数据和支持文件的信息。感兴趣的研究人员可以通过电子邮件来联系 wczy115728@163.com

的利益冲突

作者宣称没有利益冲突。

确认

这部分工作是支持广西重点实验室系统工程的防灾工程安全拨款2019 zdk052, 672707年格兰特博士基金会下的河南理工大学,河南理工大学重点实验室开放基金资助下SJF201802,国家自然科学基金资助下的中国51904093,主要研发和推广的特殊(科技)项目下的河南省拨款192102310247,和青春的探索创新的基础在格兰特NSFRF180321河南理工大学。论文的作者是非常感谢上面的支持。

Y。 F。 H。 实验和数值研究脆性岩石的震支座行为受到动荷载相结合 国际机械科学杂志》上 2020年 175年 105520年 10.1016 / j.ijmecsci.2020.105520 Mobayen 年代。 Tchier F。 一种新颖的鲁棒自适应二阶滑模跟踪控制技术不确定的动力系统匹配和无与伦比的干扰 国际期刊的控制、自动化和系统 2017年 15 3 1097年 1106年 10.1007 / s12555 - 015 - 0477 - 1 2 - s2.0 - 85016133324 Mofid O。 Mobayen 年代。 Khooban m . H。 基于自适应滑模扰动观测器控制一类分数阶混沌系统的同步 国际期刊的自适应控制和信号处理 2018年 33 462年 474年 Mobayen 年代。 Tchier F。 满秩时快速终端滑模稳定器对一类不确定非线性系统基于扰动观测器 Scientia Iranica 2017年 24 3 1410年 1418年 10.24200 / sci.2017.4123 2 - s2.0 - 85051051844 T。 X。 D。 R。 H。 非线性行为和循环荷载下岩石损伤断裂模型基于能量耗散原理 工程断裂力学 2019年 206年 330年 341年 10.1016 / j.engfracmech.2018.12.010 2 - s2.0 - 85058365791 唐ydF4y2Ba y L。 问:H。 j·G。 x。 z . C。 d . M。 水泥砂浆的单轴压缩行为及其基于能量理论的损伤本构模型 材料 2019年 12 1 14 10.3390 / ma12081309 2 - s2.0 - 85065646458 F。 杨ydF4y2Ba J。 X。 年代。 峰值应力应变能存储索引的岩爆倾向的岩石材料 国际岩石力学和采矿科学杂志》上 2019年 117年 76年 89年 10.1016 / j.ijrmms.2019.03.020 2 - s2.0 - 85063113280 t . B。 y . P。 问:F。 J。 岩石的本构模型基于能量耗散和转换的原则 阿拉伯地球科学杂志》 2019年 12 1 14 10.1007 / s12517 - 019 - 4678 - 4 2 - s2.0 - 85073639465 T。 H。 J。 Y。 变形和破坏过程的能量分析砂岩和损伤本构模型 KSCE土木工程杂志》上 2019年 23 2 513年 524年 10.1007 / s12205 - 018 - 0789 - 9 2 - s2.0 - 85058423241 Y。 M。 K。 冷冻盐水砂土的本构模型基于能量耗散理论 国际期刊的可塑性 2016年 78年 84年 113年 10.1016 / j.ijplas.2015.10.008 2 - s2.0 - 84954162936 Y。 X.-T。 研究损伤演化的完整和有节的大理石受到循环真三轴加载 工程断裂力学 2019年 215年 224年 234年 10.1016 / j.engfracmech.2019.05.011 2 - s2.0 - 85065845207 K。 Z。 B。 损伤特征的红色砂岩使用单轴压缩实验 RSC的进步 2018年 8 70年 40267年 40278年 10.1039 / c8ra06972g 2 - s2.0 - 85058147978 H。 c·J。 m·G。 x B。 岩石统计损伤本构模型考虑残余强度 国际地质力学杂志 2017年 17 1 9 10.1061 /(第3期)gm.1943 - 5622.0000680 2 - s2.0 - 85008871198 H。 H。 G。 B。 G。 一个三维的统计损伤本构模型 机械科学与技术杂志》上 2015年 29日 1 71年 77年 10.1007 / s12206 - 014 - 1211 - 2 2 - s2.0 - 84921027952 艾哈迈德 m F。 瓦格 U。 艾尔沙德 M。 罗杰斯 j . D。 热处理对动态特性的影响选择岩石类型来自巴基斯坦的盐范围 阿拉伯地球科学杂志》 2018年 11 1 13 10.1007 / s12517 - 018 - 4058 - 5 2 - s2.0 - 85057506560 Y。 G。 丽安 J。 H。 C。 温度和湿度对粘合强度的影响和shotcrete-surrounded岩石的破坏机制 建筑和建筑材料 2016年 124年 1109年 1119年 10.1016 / j.conbuildmat.2016.08.126 2 - s2.0 - 84986000942 x R。 M . M。 y . M。 y Q。 一个实验调查在岩石节理剪切疲劳损伤机制基本保持在商业周期见顶前循环加载条件 国际期刊的疲劳 2018年 106年 175年 184年 10.1016 / j.ijfatigue.2017.10.007 2 - s2.0 - 85031104248 l . X。 w·B。 问:B。 问:H。 g . H。 J。 岩石的损伤演化机制减少爆破引起的隧道深处 隧道与地下空间技术 2016年 58 257年 270年 10.1016 / j.tust.2016.06.004 2 - s2.0 - 84973447792 太阳 X。 H。 l M。 W。 对声发射试验研究砂岩岩爆的特点与不同的水分含量 中国科学技术科学 2016年 59 10 1549年 1558年 10.1007 / s11431 - 016 - 0181 - 8 2 - s2.0 - 84988568891 程ydF4y2Ba X X。 P。 Z。 损伤的微观结构和力量改变地下花岗岩wet-dry周期 Symmetry-Basel 2018年 10 1 13 10.3390 / sym10120716 2 - s2.0 - 85058993280 L.-p。 C。 L.-z。 动态特性研究的深粉砂岩在高静态压力和频繁的动态干扰 土木工程的发展 2019年 2019年 1 12 10.1155 / 2019/2087684 2 - s2.0 - 85070441556 Y。 F。 动态响应和脆性岩石震支座组合载荷作用下的破坏机理实验 岩石力学和岩石工程 2018年 51 3 747年 764年 10.1007 / s00603 - 017 - 1364 - 2 2 - s2.0 - 85034829953 K。 问:B。 G。 j . C。 J。 动态力学和断裂行为的砂岩多轴载荷下使用三轴霍普金森酒吧 岩石力学和岩石工程 2019年 52 7 2175年 2195年 10.1007 / s00603 - 018 - 1691 - y 2 - s2.0 - 85059540326 H.-b。 F。 Y。 Y。 H.-n。 数值调查水压下的岩石动态强度和失效行为监禁在SHPB试验 国际岩石力学和采矿科学杂志》上 2018年 108年 43 57 10.1016 / j.ijrmms.2018.05.008 2 - s2.0 - 85048481515 h . B。 F。 Y。 岩石力学响应和失效机理的静水力加压下震支座影响相结合 国际机械科学杂志》上 2020年 165年 105年 219年 10.1016 / j.ijmecsci.2019.105219 S.-Q。 Y。 F。 Gui 杨绍明。关铭 强度、可变形性和x射线ct机观察不同围压下的深埋地下的大理石 岩石力学和岩石工程 2016年 49 11 4227年 4244年 10.1007 / s00603 - 016 - 1040 - y 2 - s2.0 - 84976488518 T。 D。 Y。 Jianlin X。 Sufang 年代。 小明 l 研究变形特性和混合模式断裂韧性断裂准则的石膏夹层Yunying盐围压下洞穴 天然气的科学与工程》杂志上 2018年 58 1 14 10.1016 / j.jngse.2018.07.020 2 - s2.0 - 85050883084 x L。 Karakus M。 花岗岩的热力耦合损伤模型 国际岩石力学和采矿科学杂志》上 2018年 103年 195年 204年 10.1016 / j.ijrmms.2018.01.030 2 - s2.0 - 85044474615 X。 F。 Z。 岩石的形变场耦合损伤本构模型基于Hoek-Brown强度准则 损伤力学的国际期刊 2018年 27 8 1213年 1230年 10.1177 / 1056789517726838 2 - s2.0 - 85043337446 X.-l。 Karakus M。 F。 Z.-z。 岩石热损伤本构模型考虑损伤阈值和残余强度 中南大学学报 2018年 25 10 2523年 2536年 10.1007 / s11771 - 018 - 3933 - 2 2 - s2.0 - 85056285889 l X。 X。 程ydF4y2Ba 年代。 热损伤统计本构模型对深岩考虑初始压实阶段和残余强度 数学问题在工程 2019年 2019年 10 9035396 10.1155 / 2019/9035396 2 - s2.0 - 85070462209 苗族 美国J。 H。 m F。 首歌 y F。 j . T。 破碎岩石的损伤本构模型和变量水化学环境 阿拉伯地球科学杂志》 2018年 11 1 14 10.1007 / s12517 - 017 - 3373 - 6 2 - s2.0 - 85040346766 中方通过 t。 H。 Deleruyelle F。 一个。 本构模型的时间依赖的行为部分饱和岩石 电脑和土工技术 2016年 78年 123年 133年 10.1016 / j.compgeo.2016.05.004 2 - s2.0 - 84969234366 P。 y D。 y . X。 h·P。 b . X。 研究非线性损伤为高压软岩蠕变本构模型 环境地球科学 2016年 75年 1 8 10.1007 / s12665 - 016 - 5699 - x 2 - s2.0 - 84971508157 F。 程ydF4y2Ba J。 Q。 岩盐的非线性蠕变损伤模型 损伤力学的国际期刊 2019年 28 5 758年 771年 10.1177 / 1056789518792649 2 - s2.0 - 85052507292 杨绍明。关铭 Z.-Y。 下岩石动态损伤本构模型和高应力状态 冲击和振动 2017年 2017年 6 8329545 10.1155 / 2017/8329545 2 - s2.0 - 85029209613 z L。 H。 j·G。 下的花岗岩力学行为和损伤本构模型耦合的温度和动态加载 岩石力学和岩石工程 2018年 51 10 3045年 3059年 10.1007 / s00603 - 018 - 1523 - 0 2 - s2.0 - 85048048584 T。 P。 J。 X。 力学行为和热处理砂岩冲击载荷作用下的损伤本构模型 IEEE访问 2018年 6 72047年 72062年 10.1109 / access.2018.2881729 2 - s2.0 - 85056703415 Y。 F。 断断续续的损伤本构模型连接循环单轴压缩下岩石 国际岩石力学和采矿科学杂志》上 2018年 103年 289年 301年 10.1016 / j.ijrmms.2018.01.046 2 - s2.0 - 85044750490 x。 j·G。 唐ydF4y2Ba y L。 问:H。 损伤本构模型基于能量耗散对完整岩石循环荷载作用 国际岩石力学和采矿科学杂志》上 2016年 85年 27 32 10.1016 / j.ijrmms.2016.03.003 2 - s2.0 - 84960879047 f . Q。 x B。 x L。 特征的初步实验研究岩石受到三维耦合的静态和动态负载 中国岩石力学与工程学报 2011年 30. 1178年 1190年 c。 一场灾难的过程中岩石断裂 1993年 中国,北京 中国煤炭工业出版社 Y。 矿物岩石的损伤理论淤泥 1989年 Jorhat、印度 东北理工学院 x B。 d S。 岩石冲击动力学 1994年 长沙,中国 中南科技大学出版社 y L。 C。 s Y。 讨论粘弹性岩石和土壤的连续损伤本构模型 Hehai大学学报 1997年 25 114年 116年 g . Q。 x M。 Viscoelascity理论 1996年 合肥,中国 中国科学技术大学出版社 木下光男 年代。 佐藤 K。 Kawakita M。 在岩石脉冲载荷作用下的力学行为 北海道大学工学院的公告 1977年 8 51 62年 J。 问:K。 f . L。 岩石爆破理论模型和数值模拟 1999年 中国,北京 科学出版社 Lemaitre J。 Chaboche j·L。 固体材料的力学 1990年 英国剑桥 剑桥大学出版社