DOI： 10.11817/j.issn.1672-7207.2021.03.004
不同应力路径下砂岩真三轴试验及数值模拟
李江腾，刘双飞，赵远，郭群
(中南大学 资源与安全工程学院，湖南 长沙，410083)
摘要：利用TRW-3000室内真三轴试验系统开展不同应力路径下的真三轴加载、卸载试验，研究其相应的力学特性，在此基础上开展PFC3D数值模拟对比试验，探讨细观裂纹演化规律。研究结果表明：岩石最大、最小主应力差(σ1-σ3)与中间主应力σ2呈线性相关；基于Drucker-Prager准则拟合不同应力路径下岩石强度效果良好；与加载相比，卸载条件下岩石黏聚力c、内摩擦角φ均有所降低；PFC3D数值模拟试验破坏模式与室内试验破坏模式基本吻合；在不同应力路径下，数值模型剪切裂纹数与拉伸裂纹数均随ε1增大而增大，剪切裂纹比例曲线ε1随变化趋势呈“И”型，且当主应变ε1相同时，随着中间主应力σ2增大，各类裂纹数量减少；与加载相比，卸载时各类裂纹数量快速增加，剪切裂纹数占比降低，曲线由加载的“上凸”型转变为卸载的“下凹”型。
关键词：中间主应力；应力路径；数值模拟；裂纹演化
中图分类号：TU43 文献标志码：A 开放科学(资源服务)标识码(OSID)：
文章编号：1672-7207（2021）03-0693-08
True triaxial test and numerical simulation of sandstone in different stress paths
LI Jiangteng, LIU Shuangfei, ZHAO Yuan, GUO Qun
(School of Resources and Safety Engineering, Central South University, Changsha, 410083, China)
Abstract: TRW-3000 indoor true triaxial test system was used to carry out true triaxial loading and unloading tests in different stress paths to study the corresponding mechanical characteristics of the sandstone. On this basis, PFC3D numerical simulation comparison test was carried out to explore the evolution law of microscopic cracks. The results show that the difference (σ1-σ3) between the maximum and the minimum principal stress of the rock is linearly related to the intermediate principal stress σ2. The Drucker-Prager criterion has good effect in fitting rock strength in different stress paths. Compared to load path, the cohesion c and the internal friction angle φ of the rock are reduced under unloading conditions. The results of PFC3D numerical simulation experiment are consistent with those of laboratory experiment. In different stress paths, the numbers of shear cracks and tensile cracks of the numerical model increase with the increase of the maximum principal strain ε1. Shear crack ratio curve with the maximum principal strain ε1 trends show "И" type. Under the same maximum principal strain ε1, with the increase of the intermediate principal stress σ2, the numbers of all kinds of cracks decrease. Compared with loading, the number of various cracks increases rapidly and the proportion of shear cracks decreases. The curve changes from "upward convex" type under loading to "downward concave" type under unloading.
Key words: intermediate principal stress; stress path; numerical simulation; crack evolution
隧道、边坡、矿山等施工导致围岩应力状态发生变化，多表现为沿开挖工作面的应力降低，在此过程中，岩石表现出的力学性质直接影响工程的安全。传统研究方法多基于加载岩石力学理论，未充分考虑开挖卸荷对岩体工程的影响，得到的结果与工程实际结果存在差异[1-2]，故研究岩石卸荷力学特性十分重要。国内外学者通过室内卸荷试验对岩石卸荷特性进行了大量研究，获得了很多成果，如：汪斌等[3]对大理岩进行了单轴循环加卸载、峰前卸围压、峰后卸围压3种卸载试验，发现峰前与峰后卸载条件下大理岩破坏模式存在明显差异，峰前卸围压路径下大理岩强度对围压更敏感；刘建等[4]通过不同围压下砂岩峰前恒轴压卸围压试验，发现较加载而言，卸载破坏时岩石扩容表现更明显，卸荷路径下砂岩更容易被破坏，剪胀角自卸载点开始快速增大到极值；ZHOU等[5-6]基于理论分析及数值模拟发现深部圆形隧道裂隙带和非裂隙带的位置、数量及大小与深部岩体的卸荷速率和力学参数有显著关系；ZHAO等[7-8]研究了卸载速率对岩爆破坏特性的影响；ZHANG等[9]利用结晶砂岩进行不同速率下的卸载试验，发现试件应力圆随着卸荷速率的增大而右移，极限强度、峰值应变及残余强度也随之增大；DUAN等[10]基于室内试验及模拟试验研究了不同卸荷速率下花岗岩卸荷破坏模式及应变急变；李宏哲等[11]基于大理岩加载、卸载试验，建立了大理岩加载、卸载本构模型；赵国彦等[12]发现幂函数型摩尔强度准则可以很好地拟合不同应力路径下岩石卸荷力学强度；许江等[13]通过加轴压卸围压试验研究了煤岩变形和渗透特征的演化规律；沙鹏等[14]通过SEM电镜扫描发现大理岩卸载破坏表现为张性萌生—剪切交汇—贯通破坏的渐进过程。但由于受试验条件限制，试验多采用常规三轴加载方式，其受力状态与自然界中岩石三向受力状态存在差异，故常规三轴试验无法有效地模拟岩石开挖过程的真实受力情况。本文基于真三轴加载、卸载试验与PFC3D数值模拟试验，研究不同应力路径下砂岩加载、卸载力学特性、破坏模式、强度参数及裂纹演化规律等，对岩土工程的开挖具有一定的理论意义和实用价值。
1 室内试验
室内试验在中南大学TRW-3000岩石真三轴伺服试验机上完成，试件是边长为100 mm的标准立方体砂岩试样。将试件表面仔细打磨，使其不平行度和不垂直度均满足DL/T 5368—2007水利水电工程岩石试验规程要求。
室内试验分为真三轴加载试验(方案1)和真三轴卸载试验(方案2)2部分，初始应力见表1。其中，σ10，σ20，σ30分别为第一、二、三初始应力。方案1：以0.05 MPa·s-1的速率同时加载三向应力至初始应力，而后保持最小主应力σ3及中间主应力σ2不变，继续加载最大主应力σ1直至试件破坏。方案2：参照方案1进行，各向应力加载至预定值后继续加载最大主应力σ1至卸载点应力(取同组加载试验所获最大主应力σ1max的90%)保持稳定，以0.05 MPa·s-1的速率卸载最小主应力σ3至试件破坏。这2种方案的应力路径如图1所示。
表1 真三轴试验初始应力
Table 1 Initial stress of true triaxial test
图1 加载和卸载应力路径
Fig. 1 Loading and unloading stress paths
2 试验结果分析
图2 不同应力路径下砂岩应力-应变曲线
Fig. 2 Stress-strain curves of sandstone under different stress paths
2.1 应力-应变曲线
图2所示为σ30=5 MPa时不同中间主应力下砂岩真三轴加载、卸载试验应力-应变曲线。由图2可知在不同应力路径下，砂岩应力-应变曲线变化趋势基本相同：试件经过短暂的微裂隙压密后，进入弹性变形阶段，此时，(σ1-σ3)与ε1之间的关系曲线近似为直线；当应力达到弹性极限后，试件塑性变形增大，(σ1-σ3)与ε1之间的关系曲线变缓直至最终破坏；与加载相比，卸载条件下试件变形更加迅速，其中最大主应变ε1与最小主应变ε3表现最明显，同时，卸载末期岩石多发生突发性脆性破坏，伴有破碎岩块剥落，卸载破坏峰值强度较加载时也存在一定程度降低，说明卸载时岩石强度降低，更易破坏。中间主应力对应力-应变曲线的影响主要表现为：随着中间主应力σ2增大，(σ1-σ3)与ε1之间的关系曲线中的直线段斜率增大，试件弹性增强，弹性极限升高，破坏时峰值强度也相应增大；同时，中间主应力σ2抑制了中间主应变ε2发展，随着σ2增大，中间主应变ε2明显变小。此现象在真三轴卸载试验中表现得更明显。
图3 (σ1-σ3)与σ2拟合曲线
Fig. 3 Fitting curves of (σ1-σ3) and σ2
2.2 破坏强度
2.2.1 应力差与中间主应力的关系
岩石破坏时最大、最小主应力差(σ1-σ3)随中间主应力σ2变化情况如图3所示。由图3可知：不同应力路径下岩石破坏时最大、最小主应力差(σ1-σ3)均随着中间主应力σ2增大而增大，可以采用线性拟合，各拟合系数R2均在0.9以上，拟合结果良好；在σ30相同时，卸载拟合曲线位于加载拟合曲线下方，即卸载破坏时最大、最小主应力差(σ1-σ3)较加载时降低；σ30越大，加载破坏与卸载破坏时最大、最小主应力差(σ1-σ3)的差值越大。
2.2.2 Drucker-Prager系列准则
Drucker-Prager准则(简称D-P准则)[15]考虑了中间主应力和静水压力对岩石强度的影响，克服了Mohr-Coulomb准则(简称M-C准则)的主要缺点，公式简单，计算方便，广泛地应用于岩土力学与工程的分析计算。D-P系列准则可表示为
(1)
式中：
，为应力第1不变量；
，为应力偏量第2不变量；
和
为仅与岩石特性有关的试验常数，本文采用M-C屈服准则外角点外接圆的D-P准则[16-20]，参数
，
。
图4 Drucker-Prager拟合曲线
Fig. 4 Drucker-Prager fitting curves
表2 Drucker-Prager准则拟合参数
Table 2 Drucker-Prager criterion fitting parameters
将真三轴加载试验与真三轴卸载试验所得结果基于D-P准则进行数值拟合，拟合曲线如图4所示。由图4可以看出：不同应力路径下试验结果与D-P准则拟合结果吻合良好，拟合系数R2均达到了0.97。岩石加、卸载条件下的黏聚力c和内摩擦角φ如表2所示。由表2可知这2种应力路径下岩石力学参数存在差异：与加载相比，真三轴卸载试验中岩石的黏聚力c降低5.6%，内摩擦角φ降低4.1%，其原因为在真三轴卸载试验中，前期应力加载过程中试件内部微裂纹萌生、发展，产生微小塑性变形，导致试件存在一定的损伤破坏，此时，试件黏聚力c与初始试件相比存在一定的弱化，同时，在卸载过程中，随着σ3降低，试件卸载方向所受约束减小，垂直于卸荷方向的裂纹更易扩张，岩石内部颗粒间摩擦减弱，导致内摩擦角φ也发生了相应弱化，故与加载相比，真三轴卸载条件下试件黏聚力c、内摩擦角φ均有一定程度降低。
3 PFC3D数值模拟
3.1 PFC3D建模与细观参数
PFC3D程序中模型由一系列球体颗粒组成，通过对颗粒本身及颗粒间接触赋值相关属性来匹配岩石的宏观特性，程序中的平行黏结模型既可以传递法向力与切向力，又能传递弯矩，广泛用于岩土工程数值分析[20-24]。本文以室内试验所用试件实际尺寸为依据建立1个由51 316个颗粒组成的边长为100 mm的立方体模型，如图5所示。
图5 数值模型
Fig. 5 Numerical model
3.2 细观参数标定与验证
基于“试错法”获取一组匹配砂岩力学特性的细观参数(见表3)，而后将室内真三轴加载试验结果与模拟结果进行对比，进一步验证细观参数的可靠性，见表4。由表4可知：PFC3D模拟试验峰值应力、各向应变与室内试验结果基本吻合；除个别ε2外，其他参数相对误差均小于10%，即采用此组细观参数的模型与试样的宏观力学参数基本匹配，满足数值模拟试验要求。该组细观参数可以用于模拟室内试验所用砂岩的力学特性，在此基础上进行PFC3D数值对比试验，模拟试验加载、卸载步骤参照室内试验方案执行。
表3 PFC3D细观参数
Table 3 Mesoscopic parameters of PFC3D
表4 加载条件下模拟试验和室内试验峰值应力、应变对比
Table 4 Comparison of peak stress and strain between simulated and laboratory tests under loading conditions
3.3 模拟结果分析
本文模拟计算的应力路径及相关力学参数均与室内试验的匹配。
3.3.1 室内试验与模拟试验对比
图6所示为σ30=5 MPa，σ20=11 MPa，σ10=20 MPa时室内试验与模拟试验的应力-应变曲线。由图6可知：不同应力路径下模拟试验曲线与室内试验曲线基本吻合，但模拟试验中(σ1-σ3)与ε1之间的关系曲线前期近似呈线性上升，无明显的压密阶段。其原因是数值模型生成过程中内部颗粒分布较均匀，无明显缺陷，而室内试验所用试样内部原生微裂隙导致(σ1-σ3)与ε1之间的关系曲线存在明显压密阶段；同时，模拟试验采用同一个数值模型，不存在离散性，而室内试验各试件存在细微差异，导致相同应力路径下应力-应变曲线不完全重合。
图6 室内试验与模拟试验应力-应变曲线
Fig. 6 Stress-strain curves of laboratory test and simulation test
图7 室内试验与模拟试验破坏模式
Fig. 7 Failure modes of laboratory test and simulation test
图7所示为σ30=5 MPa，σ20=11 MPa，σ10=20 MPa时室内与模拟卸载试验破坏模式。由图7可知：在室内卸载试验中，试件发生拉伸-剪切复合破坏，破坏面主要为2条宏观贯通破坏面，破坏面存在大量剪切摩擦留下的岩石粉末以及少数劈裂破坏产生的片状小岩块；而模拟卸载试验中，模型由2条主裂纹带和一系列次生裂纹带共同作用引发破坏，但次生裂纹带发育不丰富，只有模型个别棱角处裂纹数量较多，破坏主要由2条主裂纹带主导，主裂纹带中剪切裂纹数量约占70%，以剪切破坏为主，这与室内卸载试验破坏模式也相互吻合。
3.3.2 裂纹演化
图8 PFC3D模拟试验裂纹演化
Fig. 8 Crack evolution in PFC3D simulation test
图8所示为σ30=5 MPa，σ10=20 MPa时不同中间主应力PFC3D模拟试验裂纹演化图。由图8(a)和图8(b)可知：在不同应力路径下，剪切裂纹数与拉伸裂纹数均随ε1增大而增大，且中间主应力σ2越大，破坏前相同应变ε1下各类裂纹数量越低。其原因是随着σ2增大，模型弹性增强，相同应变ε1下塑性变形所占比例降低。由图8(c)和图8(d)可知：在不同应力路径下，裂纹演化以剪切裂纹为主，剪切裂纹占比(剪切裂纹数占总裂纹数比例)均在70%以上，剪切裂纹占比曲线随ε1变化总体趋势呈“И”型；裂纹萌生阶段以剪切裂纹为主，随着ε1增大，试样拉伸裂纹增多，剪切裂纹占比降低，继续加载，裂纹扩展加快，新生成裂纹中剪切裂纹居多，剪切裂纹占比增大；加载后期，试样发生明显塑性变形，剪切裂纹、拉伸裂纹均随ε1快速发展，此时，剪切裂纹占比快速降低直至最终破坏，同时，卸载破坏时模型剪切裂纹占比随σ2增大而减小，拉伸裂纹占总裂纹比例上升。由图8(e)可知：加载与卸载条件下裂纹演化存在差异；与加载相比，卸载时剪切裂纹数与拉伸裂纹数增速明显加快，2种应力路径下模型破坏前裂纹数量的差值持续增大。其原因是模型卸载方向约束变弱，裂纹更容易萌生、扩展；同时，剪切裂纹占比曲线自卸载点开始快速下降，曲线由加载的“上凸”型转变为卸载“下凹”型，说明卸载时新生成裂纹中拉伸裂纹数量明显增多，所占比例增大，拉伸破坏增多。
4 结论
1) 卸载条件下岩石变形更迅速，卸载末期岩石多发生突发性脆性破坏；在不同应力路径下，砂岩破坏时最大、最小主应力差(σ1-σ3)与中间主应力σ2均呈线性相关；其破坏强度采用D-P准则拟合效果良好，与加载相比，卸载条件下岩石黏聚力c和内摩擦角φ均有所降低。
2) 不同应力路径下剪切裂纹数以及拉伸裂纹数均随ε1增大而增大，剪切裂纹比例曲线随ε1呈“И”型趋势变化；与加载相比，卸载时各类裂纹数量迅速增加，剪切裂纹占比快速降低，卸载时的裂纹中拉伸裂纹数量增多，占比增大。
3) 室内试验结果与数值模拟结果所示的力学特性、破坏模式、裂纹演化规律均较吻合；本文构建的数值模型、匹配的细观参数以及相关程序语言只需进行相应调整即可用于相关理论研究和工程应用，具有较强的实用性。
参考文献：
[1] 哈秋舲. 三峡工程永久船闸陡高边坡各向异性卸荷岩体力学研究[J]. 岩石力学与工程学报, 2001, 20(5): 603-618.
HA Qiuling. Study on the anisotropic unloading rock mass mechanics for the steep-high rock slope of the Three Gorges Project permanent shiplock[J]. Chinese Journal of Rock Mechanics and Engineering, 2001, 20(5): 603-618.
[2] 哈秋舲. 岩体工程与岩体力学仿真分析: 各向异性开挖卸荷岩体力学研究[J]. 岩土工程学报, 2001, 23(6): 664-668.
HA Qiuling. Simulation analysis for rock mass engineering and rock mass mechanics: the study on anisotropic excavation unloading rock mass mechanics[J]. Chinese Journal of Geotechnical Engineering, 2001, 23(6): 664-668.
[3] 汪斌, 朱杰兵, 邬爱清, 等. 锦屏大理岩加、 卸载应力路径下力学性质试验研究[J]. 岩石力学与工程学报, 2008, 27(10): 2138-2145.
WANG Bin, ZHU Jiebing, WU Aiqing, et al. Experimental study on mechanical properties of Jinping marble under loading and unloading stress paths[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(10): 2138-2145.
[4] 刘建, 李建朋. 砂岩高应力峰前卸围压试验研究[J]. 岩石力学与工程学报, 2011, 30(3): 473-479.
LIU Jian, LI Jianpeng. Experimental research on sandstone pre-peak unloading process under high confining pressure[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(3): 473-479.
[5] ZHOU X P, BI J. Zonal disintegration mechanism of cross-anisotropic rock mass around a deep circular tunnel under dynamic unloading[J]. Theoretical and Applied Fracture Mechanics, 2012, 60(1): 15-22.
[6] ZHOU Xiaoping, SHOU Yundong. Excavation-induced zonal disintegration of the surrounding rock around a deep circular tunnel considering unloading effect[J]. International Journal of Rock Mechanics and Mining Sciences, 2013, 64: 246-257.
[7] ZHAO X G, WANG J, CAI M, et al. Influence of unloading rate on the strainburst characteristics of Beishan granite under true-triaxial unloading conditions[J]. Rock Mechanics and Rock Engineering, 2014, 47(2): 467-483.
[8] 何满潮, 赵菲, 杜帅, 等. 不同卸载速率下岩爆破坏特征试验分析[J]. 岩土力学, 2014, 35(10): 2737-2747.
HE Manchao, ZHAO Fei, DU Shuai, et al. Rockburst characteristics based on experimental tests under different unloading rates[J]. Rock and Soil Mechanics, 2014, 35(10): 2737-2747.
[9] ZHANG Dongming, YANG Yushun, CHU Yapei, et al. Influence of loading and unloading velocity of confining pressure on strength and permeability characteristics of crystalline sandstone[J]. Results in Physics, 2018, 9: 1363-1370.
[10] DUAN Kang, JI Yinlin, WU Wei, et al. Unloading-induced failure of brittle rock and implications for excavation-induced strain burst[J]. Tunnelling and Underground Space Technology, 2019, 84: 495-506.
[11] 李宏哲, 夏才初, 肖维民. 锦屏水电站大理岩加卸荷本构模型研究[J]. 岩石力学与工程学报, 2010, 29(7): 1489-1495.
LI Hongzhe, XIA Caichu, XIAO Weimin. Research on constitutive model of marble of Jinping hydropower station under loading and unloading conditions[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(7): 1489-1495.
[12] 赵国彦, 戴兵, 董陇军, 等. 不同应力路径下岩石三轴卸荷力学特性与强度准则研究[J]. 岩土力学, 2015, 36(11): 3121-3127.
ZHAO Guoyan, DAI Bing, DONG Longjun, et al. Experimental research on mechanical characteristics and strength criterion of rock of triaxial unloading tests under different stress paths[J]. Rock and Soil Mechanics, 2015, 36(11): 3121-3127.
[13] 许江, 李波波, 周婷, 等. 加卸载条件下煤岩变形特性与渗透特征的试验研究[J]. 煤炭学报, 2012, 37(9): 1493-1498.
XU Jiang, LI Bobo, ZHOU Ting, et al. Experimental study of coal deformation and permeability characteristics under loading-unloading conditions[J]. Journal of China Coal Society, 2012, 37(9): 1493-1498.
[14] 沙鹏, 伍法权, 常金源. 大理岩真三轴卸载强度特征与破坏力学模式[J]. 岩石力学与工程学报, 2018, 37(9): 2084-2092.
SHA Peng, WU Faquan, CHANG Jinyuan. Unloading strength and failure pattern of marble under true triaxial test[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(9): 2084-2092.
[15] DRUCKER D C, PRAGER W. Soil mechanics and plastic analysis or limit design[J]. Quarterly of Applied Mathematics, 1952, 10(2): 157-165.
[16] 刘金龙, 栾茂田, 许成顺, 等. Drucker-Prager准则参数特性分析[J]. 岩石力学与工程学报, 2006, 25(S2): 4009-4015.
LIU Jinlong, LUAN Maotian, XU Chengshun, et al. Study on parametric characters of drucker-prager criterion[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(S2): 4009-4015.
[17] 邓楚键, 何国杰, 郑颖人. 基于M-C准则的D-P系列准则在岩土工程中的应用研究[J]. 岩土工程学报, 2006, 28(6): 735-739.
DENG Chujian, HE Guojie, ZHENG Yingren. Studies on Drucker-Prager yield criterions based on M-C yield criterion and application in geotechnical engineering[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(6): 735-739.
[18] 周凤玺, 李世荣. 广义Drucker-Prager强度准则[J]. 岩土力学, 2008, 29(3): 747-751.
ZHOU Fengxi, LI Shirong. Generalized Drucker-Prager strength criterion[J]. Rock and Soil Mechanics, 2008, 29(3): 747-751.
[19] 王先军, 陈明祥, 常晓林, 等. Drucker-Prager系列屈服准则在稳定分析中的应用研究[J]. 岩土力学, 2009, 30(12): 3733-3738.
WANG Xianjun, CHEN Mingxiang, CHANG Xiaolin, et al. Studies of application of Drucker-Prager yield criteria to stability analysis[J]. Rock and Soil Mechanics, 2009, 30(12): 3733-3738.
[20] 田文岭, 杨圣奇, 方刚. 煤样三轴循环加卸载力学特征颗粒流模拟[J]. 煤炭学报, 2016, 41(3): 603-610.
TIAN Wenling, YANG Shengqi, FANG Gang. Particle flow simulation on mechanical behavior of coal specimen under triaxial cyclic loading and unloading[J]. Journal of China Coal Society, 2016, 41(3): 603-610.
[21] 刘静, 李江腾. 基于颗粒流的大理岩三轴循环加卸载细观损伤特性分析[J]. 中南大学学报(自然科学版), 2018, 49(11): 2797-2803.
LIU Jing, LI Jiangteng. Analysis on meso-damage characteristics of marble under triaxial cyclic loading and unloading based on particle flow simulation[J]. Journal of Central South University(Science and Technology), 2018, 49(11): 2797-2803.
[22] HE Pengfei, KULATILAKE P H S W, YANG Xuxu, et al. Detailed comparison of nine intact rock failure criteria using polyaxial intact coal strength data obtained through PFC3D simulations[J]. Acta Geotechnica, 2018, 13(2): 419-445.
[23] DING Xiaobin, ZHANG Lianyang, ZHU Hehua, et al. Effect of model scale and particle size distribution on PFC3D simulation results[J]. Rock Mechanics and Rock Engineering, 2014, 47(6): 2139-2156.
[24] PARK B, MIN K B. Bonded-particle discrete element modeling of mechanical behavior of transversely isotropic rock[J]. International Journal of Rock Mechanics and Mining Sciences, 2015, 76: 243-255.
(编辑 陈灿华)
收稿日期： 2020 -04 -10; 修回日期： 2020 -06 -12
基金项目(Foundation item)：国家自然科学基金资助项目(51979293，51774322)；湖南省水利厅科技项目(2015131-5) (Projects(51979293，51774322) supported by the National Natural Science Foundation of China; Project(2015131-5) supported by the Science and Technology Program of Water Resources Department of Hunan Province)
通信作者：郭群，高级实验师，从事岩石力学研究；E-mail: guoqun69@163.com
引用格式： 李江腾, 刘双飞, 赵远, 等. 不同应力路径下砂岩真三轴试验及数值模拟[J]. 中南大学学报(自然科学版), 2021, 52(3): 693-700.
Citation: LI Jiangteng, LIU Shuangfei, ZHAO Yuan, et al. True triaxial test and numerical simulation of sandstone in different stress paths[J]. Journal of Central South University(Science and Technology), 2021, 52(3): 693-700.
