Experimental study on anisotropy of saturated coral sand under complex stress conditions
-
摘要: 各向异性是珊瑚砂的固有属性。对饱和南沙珊瑚砂开展了一系列不排水单向剪切试验,探究了固结应力方向角α0和单向加载方向角αm对饱和珊瑚砂不排水反应的影响。α0和αm对饱和珊瑚砂的不排水反应均有显著影响,且α0与αm的耦合作用对其不排水反应的影响更为复杂。对所有试验的应力条件,饱和珊瑚砂的超静孔压ue均呈现出先增大后减小的趋势。α0不同时饱和珊瑚砂的相变强度SPT、有效内摩擦角ϕ′PT与αm的关系存在显著差异。发现饱和珊瑚砂的相变强度SPT与无量纲参数β存在事实上的线相关性,其中,β是以α0,αm为变量的余弦函数。随着广义剪应力qg的增大,饱和珊瑚砂呈现出明显的应变硬化现象。Abstract: Anisotropy is the inherent property of coral sand. A series of undrained monotonic shear tests are carried out on the saturated Nansha coral sand by using the GDS hollow cylinder torsional shear apparatus. The effects of consolidation stress direction angle α0 and monotonic loading direction angle αm on the undrained response of saturated coral sand are investigated. The test results show that α0 and αm have significant influences on the undrained response of saturated coral sand. The undrained response characteristics of coral sand will become more complex under the coupling effects of α0 and αm. For all the test conditions considered, the excess pore water pressure ue of saturated coral sand presents contraction first and then dilatancy trend. The change of the shear resistance (SPT) and its effective angle (ϕ′PT) mobilized at the phase transformation state along with αm are significantly different for various α0. There is a unique linear correlation between SPT and the dimensionless parameter β, which is a cosine function with α0 and αm as variables. With the increase of the generalized shear stress qg, obvious strain hardening phenomenon can be observed for the samples.
-
0. 引言
岩石类材料的强度不仅与围压、应力状态等因素有关,还受加载速率的影响。通常应变率越高,强度越大。在分析爆炸冲击、地震等动力荷载作用时,由于应力波的衰减,岩土体各个区域处于不同的应力状态,应变率也存在显著差异,因而需要考虑复杂应力状态下强度的率相关特性。
目前,材料的率效应研究以单轴试验为主,学者针对岩石、混凝土、冻土等材料开展了系列研究。根据是否存在极限峰值强度,单轴动态强度准则可分为J准则和S准则[1]。对于J准则,随着应变速率的增加,材料的动态强度将无限增大;而S准则认为,当应变率极大时,材料的动态强度也将趋于稳定。Lu等[2]认为动态强度由真实动态强度和惯性动态强度组成,真实动态强度存在极限值。然而,由于在极高应变率下,试验加载难度极大,可靠数据较少,目前还很难准确测定或预测极限峰值强度。而J准则参数相对较少、易于测定。在材料达到极限峰值强度前,S准则仍可采用J准则近似,因而,目前仍广泛采用J准则。
通常,J准则以对数函数[3]、指数函数[1]或幂函数[4]等形式为主,也有部分学者采用分段函数的形式,但这增加了参数确定和模型使用的难度,函数非连续光滑也使与弹塑性本构理论结合的难度增加。
为了合理考虑岩石类材料在复杂应力状态下的动态强度规律,一些学者采用等效应变率概念[5-6],将单轴强度拓展至复杂应力状态。此外,王怀亮等采用二次函数作为子午面破坏包络线,在偏平面上采用Willam-Warnke模型,使用特征应力点法对拉压子午面方程进行修正以考虑偏平面上罗德角对动态强度的影响[7]。杜修力等基于广义非线性强度[8]框架,结合混凝土材料的非线性单轴、双轴动态强度特性,得出强度参数的率效应公式,由此建立混凝土材料的三轴动态强度准则[9]。这两种理论均假定静力条件下强度参数与单轴压缩、拉伸强度的关系在动态条件下仍然成立,从而将静力准则外推至考虑应变率效应。由于强度参数均为应变率的函数,这导致强度包线表达式复杂,增加了准则使用的难度。
本文通过整理单轴及多轴动态强度试验数据,分析各个强度参数与应变率效应的相关性,在此基础上建立动态坐标系,提出考虑应变率效应的多轴强度准则。准则中各个参数物理意义明确,模型的建立充分考虑动态强度试验规律,模型特性与试验相符,为强度的率效应提供了新的认识,可以用于预测材料的动态强度特性。
1. 单轴动态强度准则
对于砂土、岩石、混凝土等材料,其强度与加载速率密切相关,受试验条件制约,目前超高应变率下强度的发展规律及机理还没有定论。为了简便起见,本文基于J准则开展研究,即应变率越高,材料强度越高。图1为半对数坐标系下,动态强度典型的试验结果。当应变率较小时,强度保持恒定,此时称为准静态。当应变率超过一定范围后,强度产生突变,强度突变的转折点定义为临界应变率
˙ε0 。根据实测数据,强度与应变率在半对数坐标系下为线性关系,其斜率采用静态强度σs 无量纲化后,表示为kεa,描述强度的增长率。为了简单描述动力强度随应变率变化的规律,对试验结果进行归一化。采用动力增长因子(DIF, dynamic increase factor)作为强度特征量,其定义为动态强度与静态强度的比值;针对应变率,构造无量纲量
˙ε/˙ε0+1 ,表征相对应变率。采用无量纲量整理试验数据发现,动力增长因子与相对应变率在半对数坐标系下为线性关系。当应变率较高时,˙ε/˙ε0+ 1→˙ε/˙ε0 ,其斜率趋近于kεa。由于材料的初始状态为静态,代入相应的条件,直线始终经过点(1,1),由此,可以采用统一的方程描述岩石类材料由静态到动态的强度特性,DIF=σdσs=1+kεaln(˙ε˙ε0+1), (1) 式中,
σd 为动态强度,˙ε 为应变率。图2为式(1)所得动力增长因子随应变率变化的曲线,当应变率相对较低时,动力增长因子基本为1;而当应变率大于临界应变率时,随着应变率增加,DIF与应变率在半对数坐标系下近似为线性关系。临界应变率不同,曲线上升段起点不同,而曲线基本相互平行,斜率一致。当kεa不同,动态强度增长速率不同,斜率越大,强度升高越快。在动力作用下,材料单轴拉伸或压缩强度随应变率变化规律一致,式(1)可以较好地描述岩石类材料的单轴动态强度特性。
2. 多轴动态强度准则
2.1 统一强度准则
由于缺乏必要的试验基础,目前对材料的动态强度特性还没有清晰的认识,尚不具备通过理论直接建立三维动态强度准则的条件。分析应变率对强度参数的影响规律,是建立动态强度准则是行之有效的方法[11]。Yao等[12]提出的统一强度准则适用于描述多种岩土类材料的三维强度特性。该准则仅包含4个相互独立的参数,均具有明确的物理意义,可以方便地通过三轴压缩和拉伸试验确定,便于分析率效应对强度参数的影响规律,因而基于该准则进行多轴动态强度准则的研究。为了反映摩擦效应及静水压力效应,统一强度准则在子午面上采用幂函数作为破坏包线[12],
q*=Mf(ˉppr)npr。 (2) 式中 q*为三轴压缩强度;p为平均正应力;pr为参考应力;n为静水压力效应指数,描述了子午面破坏曲线的弯曲程度;Mf为参考应力pr处对应的割线斜率,反映了材料的摩擦特性。采用幂函数可以统一描述线性与非线性的强度特性,且通过固定点,方便参数确定。统一强度准则通过坐标平移描述材料的黏聚力效应,坐标平移公式为[12]
ˉσij=σij+σ0δij。 (3) 式中
σij 为应力张量;δij 为Kronecher符号。ˉσij 为平移后的应力张量,σ0 为等向拉伸强度,对应子午面破坏曲线与p轴的交点坐标。强度包线平移后满足式(2),强度包线如图3所示。在π平面,为了与SMP及Mises准则相结合,当
ˉσ3 ≥0时,定义[12]q*=α√ˉI21−3ˉI2+2(1−α)ˉI13√(ˉI1ˉI2−ˉI3)/(ˉI1ˉI2−9ˉI3)−1, (4) 式中,
ˉI1 ,ˉI2 和ˉI3 分别为坐标平移后的主应力不变量。α 的不同引起同一π 平面上强度包线由SMP曲边三角形向Mises圆过渡,α 的取值范围为[0,1],当α =0时,π平面上强度包线为SMP曲边三角形;当α =1时,π平面上强度包线为Mises圆,如图4所示。2.2 子午面动态强度特性
文献[13]开展不同应变率条件下混凝土的三轴压缩试验,首先结合该数据分析子午面动态强度特性。根据经验公式,静态等向拉伸强度近似取为0.09fc[14],令pr=fc,根据统一强度准则,子午面强度在ln(q/pr)- ln[(p+
σ0 )/pr]坐标系下为线性关系,直线斜率为静水压力系数n,与纵轴的截距为ln(Mf)。图5给出不同应变率条件下的试验结果,各条强度包线基本平行,即应变率对静水压力系数n几乎没有影响。计算各强度包线与横轴的截距,得到对应的平均正应力,若将应变率10-5视为静态,计算平均正应力差值Δp,由此可以计算出不同应变率下的动态等向拉伸强度
σd0 ,σd0=σs0+Δp, (5) 整理等效动态等向拉伸强度与应变率的关系,其在半对数空间与单轴强度试验规律一致(图6),因而可以得到表达式为
σd0=σs0+kεaσs0ln(˙ε˙ε0+1), (6) 式中,kε为等向拉伸强度增长率,通常试验均针对单轴拉伸或单轴压缩,根据经验公式,可以将强度值近似转化为等向拉伸强度[14]。
采用所得的动态等向拉伸强度重新整理试验数据,在ln(q/pr)-ln[(p+
σ0 )/pr]坐标系下,不同应变率的强度包线均近似归一到同一直线(图7)。由此可得,应变率对静水压力系数及摩擦强度影响较小。这一结论与张建民等[15]基于砂土所得的试验规律一致,即应变率变化对摩擦角影响不大,不同应变率下强度包线互相平行,应变率增大只引起黏聚力增加。根据统一强度准则,黏聚力变化可以通过坐标平移考虑,同理,可以采用动态坐标系描述应变率对强度包线的影响,不同应变率下坐标平移公式如下:
ˉσij=σij+σd0δij, (7) 式中,
σd0 为动态等向拉伸强度,可由式(6)得出。在动态坐标系下,子午面强度可采用式(2)描述,将强度包线转换至真实应力空间,即可得到不同应变率下的子午面强度(图8)。随着应变率的增加,强度包线形状相似,但是逐渐向左平移,强度逐渐升高。而随着平均正应力的升高,强度增长值明显降低,即强度的率敏感性减弱。
在等p条件下,根据式(2)可以得到动态增长因子
DIF=qdqs=[1+kεp+σs0ln(˙ε/˙ε0+1)]n, (8) 由式(8)可以得到,抗剪强度动力增长因子主要受静水压力系数及平均正应力影响。图9给出了强度动力增长因子随应变率及平均正应力演化的规律,随着平均正应力的增加,动力增长因子显著降低,这与Li等[16]所得的试验规律类似。而随着n值减小,静水压力效应越发显著,率效率越不明显。
在式(2)中代入单轴压缩应力状态(fc,0,0),可以得到单轴压缩强度,结合动态坐标系可以得到不同应变率下的强度特性,其规律与图1(a)中单轴试验规律基本一致(图10)。当超过临界应变率时,DIF与应变率在半对数空间近似为线性,这说明动态强度准则的合理性。随着ks的变化,其可以描述与单轴准则相同的动态强度特性(图2(b))。
2.3 多轴动态强度特性
在子午面动态强度的基础上,结合公式(2),(4),并假定中主应力参数
α 不随应变率变化,可得动态三维强度准则:Mf(ˉppr)npr=α√ˉI21−3ˉI2+2(1−α)ˉI13√(ˉI1ˉI2−ˉI3)/(ˉI1ˉI2−9ˉI3)−1 , (9) 式中,坐标平移公式由式(7)定义,拉压应力比
α 可以通过下式计算得出α=MeMc+3(Me−Mc)M2e, (10) 式中,Me,Mc分别为三轴拉伸、压缩应力比,对应强度包线在p-q空间的斜率。图11给出不同应率下三轴压缩及拉伸的试验数据[17],根据实测结果,三轴压缩及拉伸破坏线的斜率几乎不随应变率变化,由式(10),可以认为
α 不受应变率影响。根据式(9)可以得到不同应变率下的π平面强度包线(图12)。随着应变率增加,π平面强度包线外扩。由于等效破坏应力比减小,理论上,同一π平面的强度包线存在由曲边三角形向圆演变的趋势,在极高应变率下,材料强度特性趋向于金属材料。
将单轴拉伸应力状态(0,0,-ft)代入式(9),可以得到不同应变率下的单轴拉伸强度(图13),其规律与图2中的试验结果一致,当应变率较低时,单轴拉伸强度动力增长因子趋于1,而当应变率超过临界应变率时,动力增长因子显著增加。为了研究中主应力参数对动力增长因子的影响,图13分别给出不同中主应力参数下的单轴拉伸预测结果,根据统一强度准则,随着中主应力参数
α 的增加,单轴拉伸强度会逐渐提高。然而,不同应变率下动力增长因子却始终保持不变,这也从侧面说明,忽略率效应对α 的影响,不会显著影响材料单轴或拉伸子午面的动态强度规律。3. 单轴动态强度准则试验验证
根据前人的研究成果,搜集了多种材料的单轴动态强度试验数据(图14),不同材料强度随应变率变化规律基本一致,通过调整模型中临界应变率、强度增长率等参数,可以较好的拟合各种材料的单轴压缩及拉伸强度特性。
4. 多轴动态强度准则试验验证
在静态、准静态条件下,统一强度准则已通过砾石、岩石、混凝土等材料的强度试验进行验证[12]。而在高应变率下开展真三轴试验尚有较大的技术难度,目前率效应试验以双轴压-压或拉压为主。本文将采用混凝土动态试验对强度准则进行验证。图15为三轴压缩试验的预测和实测结果[13],试验所用试件为100 mm×100 mm×100 mm立方体,试验过程中,侧向压力保持恒定,分别为0,0.25fc,0.5fc,0.75fc,试件竖向加载直至破坏,加载应变率分别为10-5,10-4,10-3,10-2。理论预测参数为Mf=1.45,n=0.4,pr=fc,
σ0 =0.09fc,kε=0.236,ε0 =10-5。由对比结果可以看出,本文的动态强度准则可以描述应变率对子午面强度包线的影响。图16为双轴压-压试验的预测和实测对比结果,试件采用双轴加载,加载过程中,始终保持2个方向荷载比例不变,应力比分别为1∶0,1∶0.25,1∶0.5,1∶0.75,1∶1。试验过程中,在不同加载应变率下进行试验,使试件破坏。由于应力比不同,各个工况应力罗德角不同,应力比为1∶1时对应三轴拉伸状态。由于混凝土试件强度相同,参数与上文一致,因需要考虑三轴拉伸应力状态,中主应力系数
α 取为0.75,预测结果较好地反映了混凝土双轴压-压强度特性(图16)。5. 结论
本文结合单轴和多轴动态强度试验数据,基于统一强度准则进行研究,得出以下主要结论:
(1)在半对数坐标系下,单轴强度与无量纲量
˙ε/˙ε0+1 成线性关系,由此可以建立单轴动态强度公式,式中仅包含临界应变率和强度增长率两个物理意义明确且易于确定的参数。(2)针对子午面动态强度,不同应变率下强度包线在双对数坐标系下近似平行,应变率效应可以归结为对强度包线的平移,即对等向拉伸强度或黏聚力的影响。在半对数坐标系下,等向拉伸强度与无量纲量
˙ε/˙ε0+1 成线性关系,规律与单轴压缩、拉伸试验相似。由此建立动态坐标系,并将统一强度准则扩展至考虑应变率效应。(3)根据试验规律,应变率变化对中主应力参数几乎没有影响,由此将动态强度准则扩展至3维应力状态。当多轴动态强度准则退化至单轴压缩或拉伸应力状态,可以描述与试验一致的动态单轴强度规律。
(4)采用多轴动态强度准则,可以预测不同应变率下混凝土的双轴试验结果,由于相关试验数据有限,真三轴条件下的动态强度还需要进一步研究验证。
-
表 1 试验方案
Table 1 Summary of test schemes
试验编号/No. 固结条件 加载条件 备注 p′0/kPa b0 Rc α0/(°) bm αm/(°) 1 100 — 1 — 0.5 0 均等固结 2 22.5 3 45 4 67.5 5 90 6 100 0.5 1.5 0 0.5 0 各向异性固结 7 22.5 8 45 9 67.5 10 90 11 100 0.5 1.5 45 0.5 0 12 22.5 13 45 14 67.5 15 90 16 100 0.5 1.5 90 0.5 0 17 22.5 18 45 19 67.5 20 90 -
[1] 马维嘉, 陈国兴, 李磊, 等. 循环荷载下饱和南沙珊瑚砂的液化特性试验研究[J]. 岩土工程学报, 2019, 41(5): 981–988. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract17775.shtml MA Wei-jia, CHEN Guo-xing, LI Lei, et al. Experimental study on liquefaction characteristics of saturated coral sand in Nansha Islands under cyclic loading[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(5): 981–988. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract17775.shtml
[2] 朱长歧, 陈海洋, 孟庆山, 等. 钙质砂颗粒内孔隙的结构特征分析[J]. 岩土力学, 2014, 35(7): 1831–1836. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201407003.htm ZHU Chang-qi, CHEN Hai-yang, MENG Qing-shan, et al. Microscopic characterization of intra-pore structures of calcareous sands[J]. Rock and Soil Mechanics, 2014, 35(7): 1831–1836. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201407003.htm
[3] 陈海洋, 汪稔, 李建国, 等. 钙质砂颗粒的形状分析[J]. 岩土力学, 2005, 26(9): 1389–1392. doi: 10.3969/j.issn.1000-7598.2005.09.008 CHEN Hai-yang, WANG Ren, LI Jian-guo, et al. Grain shape analysis of calcareous soil[J]. Rock and Soil Mechanics, 2005, 26(9): 1389–1392. (in Chinese) doi: 10.3969/j.issn.1000-7598.2005.09.008
[4] 刘洋. 砂土的各向异性强度准则: 应力诱发各向异性[J]. 岩土工程学报, 2013, 35(3): 460–468. http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract14994.shtml LIU Yang. Anisotropic strength criteria of sand: stress-induced anisotropy [J]. Chinese Journal of Geotechnical Engineering, 2013, 35(3): 460–468. (in Chinese) http://manu31.magtech.com.cn/Jwk_ytgcxb/CN/abstract/abstract14994.shtml
[5] KATO S, ISHIHARA K, TOWHATA I. Undrained shear characterstics of saturated sand under anisotropic consolidations[J]. Soils and Foundations, 2001, 41(1): 1–11. doi: 10.3208/sandf.41.1
[6] YOSHIMINE M, ISHIHARA K, VARGAS W. Effects of principal stress direction and intermediate principal stress on undrained shear behavior of sand[J]. Soils and Foundations, 1998, 38(3): 179–188. doi: 10.3208/sandf.38.3_179
[7] WONG R K S, ARTHUR J R F. Induced and inherent anisotropy in sand[J]. Géotechnique, 1985, 35(4): 471–481. doi: 10.1680/geot.1985.35.4.471
[8] SATO K, YOSHIDA N. Effect of Principal stress direction on undrained cyclic shear behaviour of dense sand[C]// Proeeedings of the Ninth International Offshore and Polar Engineering Conference, 1999: 542–547. Brest.
[9] CHEN G X, WU Q, ZHOU Z L, et al. Undrained anisotropy and cyclic resistance of saturated silt subjected to various patterns of principal stress rotation[J]. Géotechnique, 2020, 70(4): 317–331. doi: 10.1680/jgeot.18.P.180
[10] 于艺林, 张建民, 童朝霞, 等. 定轴排水剪切试验中各向异性砂土的力学响应[J]. 岩土力学, 2011, 32(6): 1637–1642. doi: 10.3969/j.issn.1000-7598.2011.06.007 YU Yi-ling, ZHANG Jian-min, TONG Zhao-xia, et al. Behavior of anisotropic mica sand under fixed principal stress axes drained shear test[J]. Rock and Soil Mechanics, 2011, 32(6): 1637–1642. (in Chinese) doi: 10.3969/j.issn.1000-7598.2011.06.007
[11] 罗强, 李晓磊, 王忠涛. 初始各向异性砂土不排水剪切特性试验研究[J]. 大连理工大学学报, 2019, 59(6): 629–637. https://www.cnki.com.cn/Article/CJFDTOTAL-DLLG201906012.htm LUO Qiang, LI Xiao-lei, WANG Zhong-tao. Experimental study of shear behavior of inherent anisotropy sand in undrained condition[J]. Journal of Dalian University of Technology, 2019, 59(6): 629–637. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DLLG201906012.htm
[12] 杨仲轩, 李相崧, 明海燕. 砂土各向异性和不排水剪切特性研究[J]. 深圳大学学报(理工版), 2009, 26(2): 158–163. doi: 10.3969/j.issn.1000-2618.2009.02.010 YANG Zhong-xuan, LI Xiang-song, MING Hai-yan. Fabric anisotropy and undrained shear behavior of granular soil[J]. Journal of Shenzhen University (Science and Engineering), 2009, 26(2): 158–163. (in Chinese) doi: 10.3969/j.issn.1000-2618.2009.02.010
[13] 陈伟, 张吾渝, 常立君, 等. 定向剪切应力路径下击实黄土各向异性试验研究[J]. 岩石力学与工程学报, 2015, 34(增刊2): 4320–4324. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2015S2082.htm CHEN Wei, ZHANG Wu-yu, CHANG Li-jun, et al. Experimental study of anisotropy of compacted loess under directional shear stress path[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(S2): 4320–4324. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2015S2082.htm
[14] AGHAJANI F H, SALEHZADEH H. Anisotropic behavior of the Bushehr carbonate sand in the Persian Gulf[J]. Arabian Journal of Geosciences, 2015, 8(10): 8197–8217. doi: 10.1007/s12517-014-1758-3
[15] HIGHT D W, GENS A, SYMES M J. The development of a new hollow cylinder apparatus for investigating the effects of principal stress rotation in soils[J]. Géotechnique, 1983, 33(4): 355–383. doi: 10.1680/geot.1983.33.4.355
[16] CHEN G X, ZHOU Z L, PAN H, et al. The influence of undrained cyclic loading patterns and consolidation states on the deformation features of saturated fine sand over a wide strain range[J]. Engineering Geology, 2016, 204: 77–93. doi: 10.1016/j.enggeo.2016.02.008
[17] Standard Test Methods for Minimum Index Density and Unit Weight of Soils and Calculation of Relative Density: ASTM D. 4254—14[S]. 2006.
[18] Standard Test Methods for Maximum Index Density and Unit Weight of Soils Using a Vibratory Table: ASTM D. 4253—14[S]. 2006.
[19] 土的工程分类标准: GB/T 50145—2007[S]. 2008. GB/T 50145— 2007 Standard for Engineering Classification of Soil[S]. 2008. (in Chinese)
[20] 孙宗勋. 南沙群岛珊瑚砂工程性质研究[J]. 热带海洋学报, 2000, 19(2): 1–8. doi: 10.3969/j.issn.1009-5470.2000.02.001 SUN Zong-xun. Engineering properties of coral sands in Nansha Islands[J]. Tropic Oceanology, 2000, 19(2): 1–8. (in Chinese) doi: 10.3969/j.issn.1009-5470.2000.02.001
-
期刊类型引用(6)
1. 张雄辉,黄孝福,黄诗渊,刘发贵,聂亮冰,黎子玄. 压剪条件下裂缝开闭对准脆性材料断裂行为的影响机制研究. 中国农村水利水电. 2025(03): 128-135 . 百度学术
2. 梁鹏,李壮,刘俊岭,王聚贤,王骏涛. 三点弯曲试验下花岗岩应变场及损伤演化研究. 地下空间与工程学报. 2023(02): 486-494 . 百度学术
3. 杨旭旭,吴岳,靖洪文. 基于超声波实测的巷道围岩裂缝扩展和强度演变规律研究. 采矿与安全工程学报. 2021(03): 528-537 . 百度学术
4. 张超,杨楚卿,白允. 岩石类脆性材料损伤演化分析及其模型方法研究. 岩土力学. 2021(09): 2344-2354 . 百度学术
5. 何泓易. Hoek-Brown强度准则在隧道围岩卸荷试验中的应用研究. 韶关学院学报. 2021(12): 19-23 . 百度学术
6. 王思,胡晶,张雪东,任晓丹,陈祖煜,张紫涛. 不同水深水下爆炸数值及离心试验研究. 哈尔滨工业大学学报. 2020(06): 78-84 . 百度学术
其他类型引用(2)