A STUDY ON NONLINEAR CONSTITUTIVE MODEL OF FROZEN SAND CONSIDERING PARTICLE BREAKAGE
-
摘要: 颗粒破碎是粒状材料在高应力状态下的一种基本现象。为了研究冻结砂土中颗粒破碎对应力应变关系的影响,将冻结砂土视为复合颗粒材料,忽略冰的压融,考虑内摩擦角随应力状态的变化,构建一个适用于冻结砂土的考虑颗粒破碎的非线性本构模型。构建过程分为三步,首先是基于三轴剪切前后颗粒分析对冻结砂土颗粒破碎模式和产生机理进行探讨;其次是基于考虑颗粒破碎的能量平衡方程,对冻土在三轴剪切试验过程中的颗粒破碎耗能进行分析,结果表明颗粒破碎耗能随轴向应变呈双曲线变化趋势;最后应用考虑颗粒破碎的剪胀方程修正沈珠江三参数非线性模型中的体积切线模量νt,得到一个考虑颗粒破碎的非线性本构模型,模型参数可以通过单轴压缩试验和常规三轴试验确定。将原模型和修正后模型的计算结果与控制温度为-6℃,围压为1 MPa、4 MPa、6 MPa、8 MPa和10 MPa时冻结砂土的试验结果进行对比,结果表明该模型能够较好的模拟冻结砂土从低围压到高围压的应变软化特征与剪胀特征。Abstract: Particle breakage is a basic phenomenon of granular materials under high stress. In order to study the effect of particle breakage on the stress-strain relationship of frozen sand, the frozen sand is regarded as composite particle material, ignoring the melting of ice, considering the change of internal friction angle with stress state, a nonlinear constitutive model considering particle breakage for frozen sand is proposed. Firstly, based on particle analysis before and after the triaxial shearing, the fracture mode and mechanism of frozen sand particles were discussed. Secondly, based on the energy balance equation considering particle breakage, the energy consumption of particle breakage during the triaxial shear test of frozen soil was analyzed. The results show that the energy consumption of particle breakage shows a hyperbolic trend with axial strain. Finally, the volumetric tangent modulus νt in the three-parameter nonlinear model proposed by Shen Zhujiang is modified by the dilatancy equation considering particle breakage, a nonlinear constitutive model considering particle breakage suitable for frozen sand is obtained and the model parameters can be determined by uniaxial compression test and conventional triaxial test. The calculation results of the original model and the modified model are compared with the test results of frozen sand with the test temperature controlled at -6℃ and the confining pressure setting at 1 MPa, 4 MPa, 6 MPa, 8 MPa and 10 MPa, respectively. The comparison results show that the proposed model can better simulate the strain softening characteristics and dilatancy characteristics of frozen sand from low confining pressure to high confining pressure.
-
Key words:
- frozen sand /
- particle breakage /
- energy principle /
- nonlinear constitutive model
-
0. 引言
边坡稳定性一直都是岩土地质研究中经久不衰的课题,地震及降雨引发边坡坍塌、滑坡以及泥石流等边坡失稳问题尤其如此。目前,单工况下边坡稳定性的研究层出不穷,多工况耦合作用下边坡稳定性的研究则报道较少。蒋中明等[1]通过自编FISH函数实现了边坡三维饱和渗流计算分析;刘春玲等[2]利用FLAC3D对某边坡地震稳定性进行了研究;齐信等[3]基于GIS技术对汶川地震强震区地震诱发滑坡与后期降雨诱发滑坡控制因子进行了耦合分析;龚文俊等[4]采用GEOSTUDIO软件以甘肃西和Ⅲ号滑坡为例研究了降雨和地震耦合作用对滑坡稳定性的影响。边坡稳定性与人类生产生活息息相关,因此,对边坡稳定性的研究显得尤为重要。
本文以成兰铁路岷江双线特大桥高边坡为例,结合多种数值分析方法的特点,开展天然工况、地震工况、降雨工况以及降雨-地震耦合工况下边坡稳定性研究和工程安全性评价,为单工况及耦合工况下边坡稳定性研究提供参考数据及技术依据。
1. 研究区概况
成兰铁路岷江双线特大桥高边坡位于四川省九寨沟县大录乡八郎沟玉瓦寨,边坡所在的2号岩堆大致呈扇形展布于斜坡上,岩堆体纵坡较陡,坡度约在30°-38°之间,坡面长满杂树,树林茂密。该区域处于我国地貌阶梯起伏最剧烈的深切陡峻高山峡谷。铁路线路沿岷江两岸而上,在该区域穿越了4条活动断裂带,共涉及规模不等的十几条活动断裂,平均滑动速率均大于0.5 mm/a以上[5]。边坡区域地层岩性纷杂,不良地质分布广,在天然条件下会引起坡内岩土侧移导致轻微滑坡[6]。岷江地震带属高烈度地震强震频发区,地震基本烈度为Ⅷ度,地震峰值加速度为0.2 g[7]。该区域内气候多变,降水量较多,易形成暴雨。
2. 模型建立
2.1 数值模拟理论基础
本文采用分阶段弹塑性求解法、容重增加法与强度折减法的结合以及非线性动力反应分析法对边坡天然工况、降雨工况、地震工况及降雨-地震耦合工况下的应力应变特征进行研究。
2.1.1 分阶段弹塑性求解法
设置摩尔-库伦模型为材料的本构模型,先设置粘聚力和抗拉强度为大值,计算至平衡后,再设置粘聚力和抗拉强度为分析所采用的值求解,生成初始的应力场[8]。
2.1.2 强度折减法
首先将边坡强度参数粘聚力c和内摩擦角ϕ同时除以折减系数F,得到一组新的粘聚力和内摩擦角值(c′和ϕ′),然后,将c′和ϕ′值作为一组新的材料参数输入计算;不断增加F值,直到边坡达到极限状态,发生失稳破坏,这时所对应的F值即为边坡的稳定性系数[9]。用公式表示如下:
tanϕ′=(tanϕ)F (1) c′=cF (2) 式中:F为折减系数;c为岩土边坡最初的粘聚力,kPa;ϕ为岩土边坡初始内摩擦角,(°);c′、ϕ′分别为经过折减后的粘聚力和内摩擦角。
2.1.3 容重增加法
保持岩土的粘聚力及内摩擦角不变,逐步增加重力加速度G反复进行计算(逐步增加重力加速度等同于增加岩土体的容重),直至边坡破坏,此时采用的重力加速度Glimit与实际重力加速度(G0,通常取值为9.8 m/s2)之比即为该边坡的安全系数[10],即:
Fs+Glimit/G0 (3) 2.2 模型尺寸及结构面划分
研究所采用的边坡模型(见图 1)长180 m,宽80 m,高110 m。平均坡度降约0.48。坡体中上部易滑坡区域的坡度约36°。为便于分析,沿易滑坡区域对称轴取剖面1-1,易滑坡区对称轴与正北方向呈51°角。模型共由4个结构面组成。结构面Ⅰ是砂岩、板岩夹灰岩与上层的块石土以及碎石土的分界面,倾角40°;结构面Ⅱ是碎石土与块石土的分界面,倾角16°;结构面Ⅲ是块石土与粉质黏土的分界面,倾角34°;结构面Ⅳ是粉质黏土与上层覆土的分界面,倾角38°。边坡岩土物理力学参数(见表 1)根据现场勘测所获取资料确定。
表 1 边坡岩石力学参数Table 1. Natural and rainfall condition parameters工况 岩土体类型 D/(kg·m-3) c/kPa ϕ/(°) E/GPa μ σt/MPa 安全系数 天然 砂岩、板岩夹灰岩 2400 3.30E+07 55.0 2.30E+10 0.26 5.00E+05 1.31875 块石土 2200 1.05E+04 40.0 8.20E+08 0.28 0.00E+00 碎石土 2100 1.11E+04 40.0 8.10E+08 0.23 0.00E+00 粉质黏土 1900 2.20E+04 20.0 1.00E+08 0.35 0.00E+00 25 mm/d(降雨量) 砂岩、板岩夹灰岩 2450 2.81E+07 46.7 2.30E+10 0.26 4.25E+05 1.00125 块石土 2320 8.93E+03 34.0 8.20E+08 0.28 0.00E+00 碎石土 2210 9.44E+03 34.0 8.10E+08 0.23 0.00E+00 粉质黏土 2030 1.87E+04 17.0 1.00E+08 0.35 0.00E+00 35 mm/d(降雨量) 砂岩、板岩夹灰岩 2470 2.60E+07 43.8 2.30E+10 0.26 3.93E+05 0.93245 块石土 2340 8.29E+03 31.4 8.20E+08 0.28 0.00E+00 碎石土 2230 8.62E+03 31.4 8.10E+08 0.23 0.00E+00 粉质黏土 2100 1.73E+04 15.7 1.00E+08 0.35 0.00E+00 45 mm/d(降雨量) 砂岩、板岩夹灰岩 2485 2.39E+07 39.9 2.30E+10 0.26 3.73E+05 0.91225 块石土 2357 7.58E+03 29.2 8.20E+08 0.28 0.00E+00 碎石土 2247 8.11E+03 29.2 8.10E+08 0.23 0.00E+00 粉质黏土 2148 1.62E+04 15.1 1.00E+08 0.35 0.00E+00 注:D-容重;c-粘聚力;ϕ-内摩擦角;E-弹性模量;μ-泊松比;σt-抗拉强度 2.3 控制条件设置
2.3.1 天然工况
天然工况即自重工况。赋予摩尔-库仑破坏准则,采用分阶段弹塑性求解法生成初始地应力场,利用自编强度折减法计算得出天然工况下边坡的安全系数。
2.3.2 降雨工况
降雨为该地区的主要工况。在天然工况的基础上,采用强度折减法的思想对岩土体的抗剪强度参数进行折减,其后利用自编强度折减法计算得出不同降雨条件下的边坡安全系数。
2.3.3 地震工况
研究边坡位于岷江地震带,地震频发。采用弹塑性模型以及摩尔-库仑破坏准则。为减少变形过大对结果的影响,体系最大不平衡力设定为100,时间步长设置为1~5,采用最小临界阻尼比0.005、最小中心频率12.8的瑞利阻力。在模型底部导入峰值加速度为3 m/s2的汶川地震波,其达到峰值的时间为9.97 s,持续时间为20 s。
2.3.4 降雨-地震耦合工况
考虑研究区特殊的环境条件,综合强度折减法和容重增加法,在降雨工况的基础上施加地震动力荷载,比较分析在不同降雨量和地震共同作用下边坡的应力应变特征及其稳定性。各耦合工况设置情况见表 2。
表 2 耦合工况设置Table 2. Multi load-case set工况 降雨条件/(mm·d-1) 地震条件 耦合工况一 25 汶川地震波 耦合工况二 35 汶川地震波 耦合工况三 45 汶川地震波 3. 各工况计算及结果分析
3.1 天然工况及降雨工况
模拟过程中,考虑不同降雨情况对边坡稳定性的影响,设置25 mm/d、35 mm/d、45 mm/d等3种降雨情况,对摩擦角、粘聚力、体积模量、剪切模量等参数分别进行0.85、0.78、0.72倍折减。同时,考虑不同降雨情况下雨水入渗的加载作用,结合容重增加法,对土体容重进行不同情况的扩大。各降雨工况岩土体参数与其安全系数见表 1。
由表 1可知,自然工况下边坡的安全系数为1.318,有较高的安全储备,边坡处于稳定状态。降雨量为25 mm/d时,安全系数为1.001,较天然工况降低24.23%,安全储备所剩无几,边坡处于临界滑坡状态,微小的不利影响即可启动滑坡。降雨量为35 mm/d时,安全系数为0.932,较天然降雨工况降低29.29%,安全储备为0,此时边坡的抗滑力小于下滑力,边坡发生破坏。降雨量为45 mm/d时,安全系数为0.912,较天然降雨工况降低30.80%,安全储备为0,较35 mm/d降雨工况,滑体产生更大的位移。
图 2显示了降雨工况下安全系数随降雨强度的变化,由图 2可知,当降雨量为0时,边坡处于稳定状态;随着降雨量的增加,边坡稳定性逐渐降低,且当降雨量小于25 mm/d时降低幅度较大,降雨量大于25 mm/d时降低幅度较小。该曲线的拟合方程为:
y=1.87664E−4x2−0.01753x+1.31933 (4) 图 3为各降雨工况下边坡的位移分布图。从图 3可知,不同降雨工况下的位移高值区均分布于坡腰中上部,但范围随降雨量增加逐步扩大。降雨量为25 mm/d时,最大位移值为39.3 mm(见图 3a),位移高值区分布于高程2683~2714 m范围内。降雨量为35 mm/d时,最大位移值为341.7 mm(见图 3b),位移高值区分布于高程2682~2719 m范围内。降雨量为45 mm/d时,最大位移值为528.1 mm(见图 3c),位移高值区分布于高程2681~2733 m范围内。
分析可知,降雨工况与天然工况相比,位移高值区分布大致相同,随降雨量增加,雨水对土体的软化和加载作用增强,位移高值区范围沿坡面逐步扩大,且边坡稳定性逐渐降低。当降雨量达到35 mm/d时,滑体位移由39.3 mm突变至341.7 mm,表征其在雨水入渗和加载的综合作用下发生局部滑坡,这与实际工程情况较吻合。
3.2 地震工况
数值模拟模型采用弹塑性模型及摩尔-库仑准则。模型阻尼采用瑞利阻尼,其最小临界阻尼比取0.005,最小中心频率取12.8。施加汶川地震动力作用,汶川地震波的峰值加速度为3 m/s2,持时为20 s。汶川地震波加速度时程曲线见图 4。
在模型的关键位置布置4个监测点,具体位置(见图 5)为:A点位于高程2740 m的坡肩处,B点位于高程2723 m的上坡腰处,C点位于高程2670 m的下坡腰处,D点位于高程2627 m的坡脚处,记录各监测点的加速度、位移(见图 6、图 7)。
由图 6可知,边坡坡面各监测点加速度响应差异较大。坡脚D点加速度峰值最小,几乎为0;坡肩A点其次,加速度在±0.6 m/s2范围内变动;坡腰C点的加速度在±1.23 m/s2范围内变动;坡腰B点的加速度峰值最大,加速度在-3.3 m/s2~+1.4 m/s2范围内变动。
由图 7可知,坡脚D点位移几乎为0;坡肩A点在地震波持时3.5 s时开始移动,11 s时位移发生突变,14 s时位移达到最大值7.3 cm;坡腰C点在3.5 s时开始移动,12 s时位移发生突变,17 s时位移达到最大值-11.5 cm;坡腰B点在1 s时开始移动,12 s时位移发生突变,18 s时达到最大值-57 cm。
分析认为,在地震波持时18 s时各监测点位移均已达到稳定状态,表征各监测点在汶川地震作用下产生塑性变形。坡腰C点所在土层(粉质黏土)与坡腰B点所在土层(砂岩、板岩夹灰岩)在地震动力作用下的位移有极大差异,表征发生较大相对位移,即C点处土层发生局部破坏,边坡处于不安全状态。
3.3 降雨与地震耦合工况
由各耦合工况坡腰B点加速度图(见图 8)可知,从0 s开始各工况坡腰B点就有较大的加速度,其中耦合工况一最小,耦合工况二其次,耦合工况三最大。耦合工况一坡腰B点处加速度从0 s时的-1.27 m/s2锐减到0,然后在-1.1~1.08 m/s2范围内变动;耦合工况二坡腰B点处加速度从0 s时的4.71 m/s2锐减到0,然后在-2.01~3.15 m/s2范围内变动;耦合工况三坡腰B点处加速度从0 s时的-7.23 m/s2锐减到0,然后在-4.02~3.15 m/s2范围内变动。
图 9显示了各耦合工况下坡腰B点的位移。由图 9可知,耦合工况一、二、三的坡腰B点位移在0~17 s的范围内大致呈线性增加。耦合工况一坡腰B点处的位移在17.8 s时达到最大值1.9 m,之后趋于稳定;耦合工况二坡腰B点处的位移在18.3 s时达到最大值5.2 m,之后趋于稳定;耦合工况三坡腰B点处的位移在18.5 s时达到最大值13.6 m,之后趋于稳定。
由降雨与耦合工况坡腰B点位移对比图(见图 10)可知,降雨-地震耦合工况下坡腰的位移明显比降雨工况大。降雨量为25 mm/d时,坡腰位移在降雨工况下为0.039 m,在耦合工况下为1.9 m;降雨量35 mm/d时,坡腰位移在降雨工况下为0.34 m,耦合工况下为5.2 m;降雨量45 mm/d时,坡腰位移在降雨工况下为0.53 m,在耦合工况下为13.6 m。
分析认为,降雨工况使得边坡土体的容重和含水率增加,粘聚力和内摩擦角降低,地震动力作用效应得到了放大,导致在降雨-地震耦合工况下,边坡发生破坏,其位移量达到单一降雨工况和单一地震工况时的数十倍。
4. 结论
结合工程实例,运用FLAC3D有限差分软件详细地对成兰铁路岷江双线特大桥高边坡进行建模,并模拟、分析了天然、降雨、地震以及降雨-地震耦合等4种工况下边坡的位移、加速度等变化特征,对坡体的稳定性进行安全性评价。
通过天然工况和降雨工况坡体稳定性分析,得到天然工况和不同降雨工况的强度折减安全系数表以及降雨量与安全系数关系式,表明随降雨量增加,边坡稳定性逐步降低,为量化分析降雨量与滑坡的关系提供基础数据。
通过对比不同降雨工况与不同耦合工况,得到降雨工况和耦合工况的关系曲线,表明降雨和地震的共同作用,会大大削弱边坡的稳定性,使边坡稳定性较降雨和地震单独作用时要低,极其容易造成滑坡。故须对类似特殊环境的边坡采取合适的工程措施,增加其稳定性。
较详细地对成兰铁路岷江双线特大桥高边坡进行数值模拟分析,为单工况以及耦合工况下边坡的稳定性研究提供理论依据。同时,也对大巴山、龙门山、龙泉山等其他滑坡频发山区的研究和预警提供借鉴。
-
表 1 与单轴压缩模量Mt与切线模量Et相关参数表
Table 1. Parameters related to the uniaxial compression modulus Mt and tangent modulus Et
σ3/MPa k km h i j 1.0 0.1700 0.1358 0.0017 4.0 0.1400 0.1219 0.0017 6.0 0.47 0.71 0.1200 0.1185 0.0017 8.0 0.1000 0.1121 0.0018 10.0 0.0980 0.0966 0.0025 表 2 与切线泊松比νt相关参数表
Table 2. Parameters related to the tangent poisson ratio νt
σ3/MPa a e f 其他参数 1.0 1.9006 0.0006 0.0227 4.0 1.0356 0.0005 0.0193 b=1.0200 6.0 0.7825 0.0004 0.0164 s=0.0009 8.0 0.603 0.0003 0.0144 t=0.0012 10.0 0.4638 0.0002 0.0125 -
[1] 屈智炯.土的细观结构与力学性关系的研究进展[J].水电站设计, 1993, 9(3):77~81. http://www.cqvip.com/Main/Detail.aspx?id=1172182QU Zhijiong. Research progress on the relationship between soil microstructure and mechanical properties[J]. Design of Hydroelectric Power Station, 1993, 9(3):77~81. (in Chinese) http://www.cqvip.com/Main/Detail.aspx?id=1172182 [2] Lee K L, Seed H B. Drained strength characteristics of sands[J]. Journal of the Soil Mechanics and Foundations Division, 1967, 93(6):117~141. http://cn.bing.com/academic/profile?id=b04cc251d41a643515143f9b4ab4a7b0&encoded=0&v=paper_preview&mkt=zh-cn [3] Vesic A S, Clough G W. Behavior of granular materials under high stresses[J]. Journal of Soil Mechanics & Foundations Div, 1968, 94(3):661~688. http://cn.bing.com/academic/profile?id=9cf8480c204047368e5b4313c99f1286&encoded=0&v=paper_preview&mkt=zh-cn [4] 杨礼宁, 姜振泉, 张卫强, 等.高温作用后砂岩力学性质研究[J].地震工程学报, 2016, 38(2):299~302. doi: 10.3969/j.issn.1000-0844.2016.02.0299YANG Lining, JIANG Zhenquan, ZHANG Weiqiang, et al. Mechanical properties of sandstone after high temperature[J]. China Earthquake Engineering Journal, 2016, 38(2):299~302. (in Chinese with English abstract) doi: 10.3969/j.issn.1000-0844.2016.02.0299 [5] 郭庆国.关于粗粒土抗剪强度特性的试验研究[J].水利学报, 1987, (5):59~65. doi: 10.3321/j.issn:0559-9350.1987.05.009GUO Qingguo. Experimental study on shear strength characteristics of coarse-grained soil[J]. Journal of Hydraulic Engineering, 1987, (5):59~65. (in Chinese) doi: 10.3321/j.issn:0559-9350.1987.05.009 [6] 蒙进, 屈智炯.高压下冰碛土的颗粒破碎及应力应变关系[J].成都科技大学学报, 1989, (1):17~22, 56. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=SCLH198901002&dbname=CJFD&dbcode=CJFQMENG Jin, QU Zhijiong. Stress-strain behaviour of glacial till under high confining pressure[J]. Journal of Chengdu University of Science and Technology, 1989, (1):17~22, 56. (in Chinese with English abstract) http://kns.cnki.net/KCMS/detail/detail.aspx?filename=SCLH198901002&dbname=CJFD&dbcode=CJFQ [7] 汪稔, 孙吉主.钙质砂不排水性状的损伤-滑移耦合作用分析[J].水利学报, 2002, 33(7):75~78. doi: 10.3321/j.issn:0559-9350.2002.07.013WANG Ren, SUN Jizhu. Damage-slide coupled interaction behavior of undrained calcareous sand[J]. Journal of Hydraulic Engineering, 2002, 33(7):75~78. (in Chinese with English abstract) doi: 10.3321/j.issn:0559-9350.2002.07.013 [8] Chavez C, Alonso E E. A constitutive model for crushed granular aggregates which includes suction effects[J]. Soils and Foundations. 2003, 43(4):215~227. doi: 10.3208/sandf.43.4_215 [9] Varadarajan A, Sharma K G, Abbas S M, et al. Constitutive model for rockfill materials and determination of material constants[J]. International Journal of Geomechanics, 2006, 6(4):226~237. doi: 10.1061/(ASCE)1532-3641(2006)6:4(226) [10] 米占宽, 李国英, 陈铁林.考虑颗粒破碎的堆石体本构模型[J].岩土工程学报, 2007, 29(12):1865~1869. doi: 10.3321/j.issn:1000-4548.2007.12.019MI Zhankuan, LI Guoying, CHEN Tielin. Constitutive model for rockfill material considering grain crushing[J]. Chinese Journal of Geotechnical Engineering, 2007, 29(12):1865~1869. (in Chinese with English abstract) doi: 10.3321/j.issn:1000-4548.2007.12.019 [11] 迟世春, 贾宇峰.土颗粒破碎耗能对罗维剪胀模型的修正[J].岩土工程学报, 2005, 27(11):1266~1269. doi: 10.3321/j.issn:1000-4548.2005.11.006CHI Shichun, JIA Yufeng. Rowe's stress-dilatancy model modified for energy dissipation of particle breakage[J]. Chinese Journal of Geotechnical Engineering, 2005, 27(11):1266~1269. (in Chinese with English abstract) doi: 10.3321/j.issn:1000-4548.2005.11.006 [12] Rowe P W. The stress-dilatancy relation for static equilibrium of an assembly of particles in contact[J]. Proceedings of the Royal Society A:Mathematical, Physical and Engineering Sciences, 1962, 269(1339):500~527. doi: 10.1098/rspa.1962.0193 [13] 李兆明.考虑颗粒破碎的粗粒土本构模型[D].大连: 大连理工大学, 2007.LI Zhaoming. Constitutive model for coarse granular soil incorporating particle breakage[D]. Dalian: Dalian University of Technology, 2007. (in Chinese with English abstract) [14] 马巍, 吴紫汪, 常小晓, 等.围压作用下冻结砂土微结构变化的电镜分析[J].冰川冻土, 1995, 17(2):152~158. http://www.cnki.com.cn/Article/CJFDTOTAL-BCDT502.008.htmMA Wei, WU Ziwang, CHANG Xiaoxiao, et al. Analysis of microstructural changes in frozen sandy soil under confining pressures using scaning electronic microscope[J]. Journal of Glaciology and Geocryology, 1995, 17(2):152~158. (in Chinese with English abstract) http://www.cnki.com.cn/Article/CJFDTOTAL-BCDT502.008.htm [15] 马玲, 齐吉琳, 余帆, 等.冻结砂土三轴试验中颗粒破碎研究[J].岩土工程学报, 2015, 37(3):544~550. http://d.old.wanfangdata.com.cn/Periodical/ytgcxb201503020MA Ling, QI Jilin, YU fan, et al. Particle crushing of frozen sand under triaxial compression[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(3):544~550. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/ytgcxb201503020 [16] 郑郧, 马巍, 邴慧.冻融循环对土结构性影响的机理与定量研究方法[J].冰川冻土, 2015, 37(1):132~137. http://d.old.wanfangdata.com.cn/Periodical/bcdt201501014ZHENG Yun, MA Wei, BING Hui. Impact of freezing and thawing cycles on the structures of soil and a quantitative approach[J]. Journal of Glaciology and Geocryology, 2015, 37(1):132~137. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/bcdt201501014 [17] Legay J M. Du sac de billes au tas de sable Etienne Guyon, Jean-Paul Tioadec Odile Jacob, Paris, 1994, 236 p., 140 F.[J]. Natures Sciences Sociétés, 1997, 5(2):77. [18] 张建民.砂土的可逆性和不可逆性剪胀规律[J].岩土工程学报, 2000, 22(1):12~17. doi: 10.3321/j.issn:1000-4548.2000.01.002ZHANG Jianmin. Reversible and irreversible dilatancy of sand[J]. Chinese Journal of Geotechnical Engineering, 2000, 22(1):12~17. (in Chinese with English abstract) doi: 10.3321/j.issn:1000-4548.2000.01.002 [19] Ueng T S, Chen T J. Energy aspects of particle breakage in drained shear of sands[J]. Géotechnique, 2000, 50(1):65~72. doi: 10.1680/geot.2000.50.1.65 [20] Höeg K, Prévost J H. Soil mechanics and plasticity analysis of strain softening[J]. Géotechnique, 1975, 25(2):279~297. doi: 10.1680/geot.1975.25.2.279 [21] 沈珠江.考虑剪胀性的土和石料的非线性应力应变模式[J].水利水运科学研究, 1986, (4):1~14. http://www.cnki.com.cn/Article/CJFDTOTAL-SLSY198604000.htmSHEN Zhujiang. A nonlinear dilatant stress-strain model for soil and rock materials[J]. Hydro-Science and Engineering, 1986, (4):1~14. (in Chinese) http://www.cnki.com.cn/Article/CJFDTOTAL-SLSY198604000.htm -