Triaxial creep characteristics and damage model for red sandstone subjected to freeze-thaw cycles under different water contents
-
摘要: 针对寒区自然环境对岩体长期稳定性的影响,以川藏铁路沿线的红砂岩为研究对象,对不同含水率的岩样进行了冻融后核磁共振检测及三轴蠕变试验。试验结果表明,含水率的增加促进了冻融循环过程中红砂岩内各尺寸孔隙的发育扩展,且冻融后红砂岩中孔的发育程度最高；在长期荷载作用下,冻融后红砂岩的时效变形特征呈现明显的软化趋势,随着含水率的增加,冻融后红砂岩的蠕变量及蠕变速率随之增大,而荷载等级及长期强度显著降低,破坏断面由单一剪切向张拉劈裂型逐渐演变。根据试验结果,基于损伤力学原理建立一种考虑冻融及含水劣化的西原模型,并通过弹塑性力学理论将其扩展为三维形式,采用通用全局优化算法对模型进行参数辨识,借此揭示了模型参数及劣化系数随冻融循环次数及含水率的变化规律。通过FLAC3D岩石三轴蠕变仿真模拟,验证模型能较好地反映冻融后不同含水率红砂岩全阶段蠕变特征。研究结果可为川藏线寒区岩土工程长期稳定性提供科学依据及技术参考。Abstract: In view of the influences of the natural environment in cold regions on the long-term stability of rock mass, the red sandstone samples along Sichuan-Tibet railway are taken as the research object, and the nuclear magnetic resonance and triaxial creep experiments are conducted on the samples subjected to freeze-thaw cycles under different water contents. The experimental results show that the increase of water content promotes the development and expansion of pores in red sandstone during freeze-thaw cycles, and the development degree of mesopore in red sandstone after freeze-thaw cycles is the highest. The time-dependent deformation characteristics of red sandstone after freeze-thaw cycles present obvious softening trend, and the creep strain and rate of red sandstone subjected to freeze-thaw cycles increase with the increase of water content, while the load capacity and the long-term strength decrease significantly, and the fracture surface evolves from single shear to tension splitting gradually. According to the experimental results, a modified Nishihara model considering freeze-thaw damage and water deterioration is proposed based on the principle of damage mechanics, which is further extended to a three-dimensional creep constitutive equation based on the elastoplastic mechanics. The parameters of the model are identified by using the universal global algorithm with the creep data, and the variation trend of the parameters and deterioration coefficients with the number of freeze-thaw cycles and water content is revealed. Finally, the model is written into FLAC3D to realize the numerical simulation of triaxial creep experiments of rock, and it is verified that the model can well reflect the creep characteristics at all creep stages of red sandstone subjected to freeze-thaw cycles under different water contents. The research results can provide scientific basis and technical reference for the long-term stability of geotechnical engineering in cold regions along Sichuan-Tibet highway.
-
Keywords:
- red sandstone /
- triaxial creep /
- damage mechanics /
- constitutive model
-
0. 引言
近年来，中国大力发展南海地区基础设施建设，钙质砂作为沿海地区的主要岩土体受到广泛使用。有别于陆源石英砂，钙质砂作为难溶类碳酸盐（多为碳酸钙）物质的集合体，具有形状不规则、内孔隙高、颗粒易破碎、棱角度高等性质[1-3]。国内外学者对于钙质砂的力学特性开展了大量研究：张晨阳等[4]分析了细粒含量对钙质砂力学特性的影响，认为钙质砂峰值偏应力与剪胀性随细粒含量的增加而逐渐减小；黄宏翔等[5]得出在相同条件下，钙质砂的抗剪强度和残余强度均大于石英砂；Rui等[6]发现不规则颗粒形状的钙质砂表现出较慢的累积孔隙压力速率，显示出更高的抗液化性。对于钙质砂的颗粒破碎特性，研究表明在低围压情况下就会发生颗粒破碎，颗粒形状对其力学性质的影响较大，在高围压情况下钙质砂颗粒破碎情况加剧，颗粒破碎对其力学性质影响愈加显著[7]。Coop等[8]指出钙质砂的颗粒破碎程度受轴向应力大小及初始颗粒形状影响；汪轶群等[9]在三轴试验中发现随着初始围压的增大，钙质砂颗粒破碎程度加大，土样整体剪胀趋势减小，且破碎引起的能量耗散增加。
上述研究为以钙质砂的颗粒特性为切入点，对其力学性质进行了研究。但主要以钙质砂作为纯净砂进行分析，而仅考虑某一粒径、某一种类砂土多与实际工程情况不符，存在一定的局限性。中国南海地区广泛分布源于岛礁中的碳酸盐土体和沿岸地区的非碳酸盐类土所结合的碳酸盐混合土[10]，其力学性质势必由两种土体共同决定。李小刚等[11]通过三轴试验发现，混合砂峰值强度和相对破碎率随石英砂含量的增加而减小，黏聚力及峰值内摩擦角随碳酸盐含量的增加而增加。Lade等[12]指出土体中粉土、粉砂等细颗粒的含量是海底斜坡、土坝和路堤等产生液化破坏的主要原因。含细颗粒土体的试验同样表明细颗粒含量、粒径分布和颗粒形状特征是影响砂土抗剪强度和液化反应的主要因素[13-15]，且细颗粒含量与粉土特性对砂土静态液化行为的影响存在耦合关系[16]。此外，砂土中不同颗粒特征的粗颗粒和细颗粒相互作用，对混合土体的整体力学性质也产生显著影响[17-18]。
综上所述，细颗粒对钙质混合土体力学性质影响显著，而现有研究中多数仅考虑钙质细粒含量对钙质砂力学特性的影响，在对钙质混合砂的研究中，细颗粒特征不同对混合砂的影响程度亦不同。基于此，本文采用角形石英粉和圆形玻璃珠作为细颗粒，分别与钙质砂混合，对钙质混合砂进行不同轴向压力条件下的环剪试验，探究不同细颗粒含量及形状对钙质混合砂峰值强度、残余强度和轴向变形的影响，对钙质混合砂的颗粒破碎特征进行研究分析，进一步探讨颗粒破碎特征对钙质混合砂性质的影响规律，从而为南海地区工程建设提供参考。
1. 试验材料及方案
1.1 试验材料
试验所用钙质砂取自南海区域某岛礁附近。试验前先将砂样烘干，采用标准筛进行筛分，得到粒径范围为0.1~1 mm的钙质砂颗粒。玻璃珠因具有与土体相似的矿物学特征[19]，多用于模拟砂土，其形状接近理想的球形；石英粉则更具棱角和不规则性。因此，采用粒径小于0.075 mm的石英粉及玻璃珠模拟混合砂中不同形状的石英细颗粒。3类基本砂样的物理指标参数见表 1，颗粒级配曲线见图 1。按试样总质量的0%，10%，20%，30%，40%，将石英粉、玻璃珠细颗粒分别与钙质砂均匀混合得到9种试样：含石英粉试样，由J10，J20，J30，J40表示；含玻璃珠试样，由Y10，Y20，Y30，Y40表示；0表示纯钙质砂试样。
表 1 基本砂样的物性指标Table 1. Physical properties of basic sand samples试样 圆度 宽长比 D10/
μmD50/
μmD60/
μm不均匀系数Cu 曲率系数Cc 颗粒相对质量密度 钙质砂 0.517 0.541 374.5 637.7 643.5 1.72 1.28 2.75 玻璃珠 0.924 0.893 35.1 50.8 58.6 1.67 0.92 2.64 石英粉 0.635 0.746 27.2 45.7 51.3 1.89 1.05 2.66 注：①圆度为单颗粒投影面积与等周长圆面积之比；宽长比为单颗粒最小Feret直径与最大Feret直径之比。②圆度、宽长比取采样（1000粒）均值。 ![]()
图 1 3种砂样颗粒级配曲线Figure 1. Grain-size distribution curves of three kinds of sand particles为定量化表示土体形态特征，多通过颗粒形状扫描仪器，根据颗粒轮廓图像建立二维形状参数（圆度、宽长比、粗糙度等），能够较好地表述颗粒形状与其力学性质之间的关系[20-22]。对于二元颗粒混合物，参考Yang等[23]提出的组合形状参数概念，认为二元混合物组合形状参数由各颗粒材料的形状特征值及其对应的质量占比决定，组合圆度、组合宽长比越接近于1，表明颗粒整体越接近于圆形、越规整。如钙质混合砂J30的组合圆度、组合宽长比计算分别为
Sc−J30=Sc−J×Fc+Sc−G×(1−Fc)=0.635×30%+0.517×70%=0.552, (1) Ar−J30=Ar−J×Fc+Ar−G×(1−Fc)=0.746×30%+0.541×70%=0.603 (2) 图像颗粒分析仪如图 2所示，扫描基本砂样轮廓图如图 3所示，混合砂样特性指标如表 2所示。
表 2 混合砂样特性指标Table 2. Characteristics of mixed sand samples砂样 相对密实度Dr/% 最大孔隙比emax 最小孔隙比emin 组合圆度Sc 组合宽长比Ar 0 60 1.109 0.812 0.517 0.541 J10 60 1.010 0.707 0.529 0.562 J20 60 0.893 0.553 0.541 0.582 J30 60 0.794 0.467 0.552 0.603 J40 60 0.743 0.430 0.564 0.623 Y10 60 0.977 0.724 0.558 0.541 Y20 60 0.856 0.605 0.598 0.576 Y30 60 0.744 0.512 0.639 0.611 Y40 60 0.685 0.491 0.680 0.647 1.2 试验方案
试验采用HJ-1型环剪仪，其试样尺寸为外径100 mm，内径60 mm，高20 mm。天然状态下土体的相对密实度在20%~80%，且天然密实度多为中密状态[24]，故考虑实际工程情况，控制钙质混合砂的相对密实度Dr =60%（中密），进行不同轴向压力条件下的环剪试验。试验步骤如下：①采用漏斗法将试样装入容器。首先，将环剪仪的上、下剪切盒固定在一起，然后将试样缓慢倒入漏斗之中，通过漏斗将砂样以零高度的落差缓慢倒入盒中，最后将试样抚平，使其高度一致。②采用浸水法饱和试样。试样装好后倒入蒸馏水，使其完全浸没，浸泡4 h确保试样饱和度达到要求。浸泡完成后将仪器组装完成。③逐级施加100，200，300 kPa的轴向压力，待固结完成后，以1°/min的剪切速率进行环向剪切试验。
2. 试验结果分析
2.1 剪切应力-剪切位移曲线分析
试验结束后，得到不同细颗粒含量及形状条件下的剪切应力-剪切位移曲线，如图 4所示。由图 4（a）可得，纯钙质砂的剪切应力随着剪切位移的增加而逐渐增加，呈近线性增长，随后，剪切应力的增长趋势变缓。当轴向压力为100，200，300 kPa时，剪切应力分别在剪切位移为10.5，11.9，13.6 mm时达到峰值强度，其值分别为100.4，177.9，251 kPa；剪切位移分别在28，30，31 mm处得到土样的残余强度，其值分别为88.23，159.4，213.51 kPa。由各应力-位移图线得：各组试样的抗剪强度均随轴向压力的增加而增大；当剪切位移相同时，轴向压力越大，试样的剪切应力越大；同时，随轴向压力的增加，各组试样峰值强度出现延后性，且各组试样的应力-位移曲线均呈现软化特性。
对不同细颗粒含量及形状条件下钙质混合砂的峰值强度进行统计并绘制图 5。如图 5（a）所示，随着细颗粒含量的增加，混合砂峰值强度逐渐降低。与纯钙质砂相比，细颗粒含量为10%时，混合砂的峰值强度略小于纯钙质砂，其中，含石英粉混合砂峰值强度最大降幅为12%，含玻璃珠混合砂峰值强度最大降幅为16%，此时细颗粒填充土体骨架的孔隙部分，土体骨架仍由钙质砂颗粒组成；细颗粒含量为20%时，含石英粉混合砂峰值强度最大降幅为17%，含玻璃珠混合砂峰值强度最大降幅为29%；细颗粒含量为30%时，峰值强度持续降低，含石英粉混合砂峰值强度最大降幅为22%，含玻璃珠混合砂峰值强度最大降幅扩大至36%；当细颗粒含量增加至40%，峰值强度发生大幅度降低，含石英粉混合砂峰值强度最大降幅为30%，含玻璃珠混合砂峰值强度最大降幅则高达47%，表明此时土体骨架由钙质砂颗粒承担转变为钙质砂、细颗粒共同承担。此外，如图 5（b），（c）所示，峰值强度受细颗粒形状影响，随组合圆度和组合宽长比增加而逐渐减小，含玻璃珠混合砂出现更大降幅。结合剪切应力-剪切位移曲线（图 4）中含玻璃珠混合砂的曲线波动更为明显，表明细颗粒形状的改变影响了骨架颗粒间的接触状态以及力链传递，进而降低土体骨架的稳定性，圆形细颗粒的加入使得土体骨架的稳定性更差。各组试样峰值强度及下降幅度如表 3所示。
![]()
图 5 不同细颗粒含量及形状下的峰值强度Figure 5. Peak intensities at different contents and shapes of fine particles表 3 不同轴向压力下试样的峰值强度及下降幅度Table 3. Peak strengths and decrease amplitudes of samples under different axial pressures轴压 类别 0 J10 J20 J30 J40 Y10 Y20 Y30 Y40 100 kPa 峰值强度/kPa 100.4 93.6 90.7 78.0 70.2 84.8 71.7 64.3 53.6 降幅/% — 7 10 22 30 16 29 36 47 200 kPa 峰值强度/kPa 177.9 156.5 148.7 142.3 139.4 148.7 133.6 120.9 105.3 降幅/% — 12 16 20 22 16 25 32 41 300 kPa 峰值强度/kPa 251.0 241.3 208.1 198.4 188.6 232.0 201.3 185.2 158.9 降幅/% — 4 17 21 25 8 20 26 37 2.2 不同细粒含量及形状下的强度指标特征
图 6为混合砂内摩擦角随组合形状参数的变化规律。纯钙质砂的内摩擦角正切值为0.74（φ=36.4°），混合砂内摩擦角正切值随细颗粒含量、组合圆度及组合宽长比增加均表现线性减小趋势。在细颗粒含量为40%时，组合宽长比和组合圆度达到最大值，混合砂的整体颗粒形状趋于规整，含石英粉混合砂内摩擦角正切值达到最小值0.57（φ=29.7°），其降幅为23%；含玻璃珠混合砂的内摩擦角正切值为0.53（φ= 27.8°），其降幅为28%，此时表现出更小的抗剪强度。原因在于，细颗粒的填充降低了土体骨架的稳定性，在粗颗粒间起到润滑的作用，降低了颗粒间的摩擦，因而内摩擦角逐渐减小；玻璃珠因颗粒形状圆滑，随着细颗粒含量的增加，混合砂的组合形状参数逐渐接近圆形颗粒的形状参数，颗粒之间运动产生的逐渐摩擦力较小，而石英粉的颗粒形状较不规则，在颗粒间的相对运动中，不利于滚动，需要较多的力进行滑动，故含玻璃珠混合砂的内摩擦角呈现更大的降幅。
![]()
图 6 不同细颗粒含量及形状下的内摩擦角Figure 6. Angles of internal friction at different contents and shapesof fine particles钙质砂为无黏性土，但由于具有不规则的颗粒形状，颗粒之间的嵌固咬合作用产生了不可忽视的表观黏聚力[25]。图 7为混合砂表观黏聚力随颗粒含量及组合形状参数的变化规律。试验发现，随着细颗粒含量的增加，混合砂表观黏聚力由纯钙质砂时的28 kPa开始逐渐降低，含石英粉混合砂的表观黏聚力降幅较小，最终降至17.2 kPa，而含玻璃珠的混合砂表观黏聚力呈现出较大的下降趋势，最终仅为1.2 kPa。其原因是，随着细颗粒含量的增加，钙质砂颗粒所形成的骨架逐渐被细颗粒填充，减少了钙质砂颗粒之间的接触，进而降低了钙质砂颗粒之间的嵌固咬合作用；当细颗粒含量为40%时，土体骨架由钙质砂及细颗粒共同承担，由于外形更接近于球形的玻璃珠细颗粒并不具有咬合能力，导致含玻璃珠混合砂的表观黏聚力大幅下降。
![]()
图 7 不同细颗粒含量及形状下的表观黏聚力Figure 7. Apparent cohesions at different contents and shapes of fine particles2.3 不同细粒含量及形状下的残余强度特征
为定量化表述钙质混合砂的残余强度特性，采用软化系数IB分析混合砂在剪切试验后的应变软化程度。IB越大，说明砂样的强度在峰值后衰减程度越大，软化越明显，试样的变形稳定性越低。其公式为[26]
IB=σdp−σdrσdp (3) 式中：σdp为剪切破坏时的峰值强度（kPa）；σdr为剪切破坏后残余强度（kPa）。
图 8为不同细颗粒含量和形状条件下混合砂的软化系数。随着轴向压力的增大，各组混合砂的软化系数均呈现先减小后增加的趋势，轴向压力为200 kPa时，各组试样软化系数最小，试样相对稳定。细颗粒含量为10%时软化系数达到最大值，混合砂的强度在峰值后变化最大，软化程度明显。细颗粒含量超过10%后，混合砂的软化系数随着细颗粒含量的增加而减小，且在细颗粒含量为40%时软化系数达到最小值。含石英粉混合砂的软化系数明显小于含玻璃珠混合砂，表明含石英粉混合砂的变形稳定性较高，剪切应力在峰值强度后的衰减较少。混合砂中玻璃珠含量为10%时，软化最为明显，软化系数达到最大值，其数值均在0.14以上，其原因主要是由于细颗粒含量的增加，充填在钙质砂颗粒所组成骨架的孔隙中，而玻璃珠因其颗粒形状较为圆滑，在颗粒运动中起到了润滑作用，使得含玻璃珠的钙质砂混合砂受到较小的力时就能发生颗粒运动，因此试样的变形稳定性较低，在强度达到峰值后有着较大的衰减程度。与张晨阳等[4]研究结果类似，在细粒含量为10%时稳定性最差，后续随着细颗粒含量的不断增加，钙质砂的软化系数逐渐降低，试样稳定性逐渐提高。
![]()
图 8 不同轴向压力条件下试样的脆性指标Figure 8. Brittleness indexes of samples at different axial pressures2.4 不同细粒含量及形状下的变形特征
钙质混合砂在环剪过程中的轴向位移变化可在一定程度上反映出其剪胀特征。图 9为轴向压力100 kPa时环剪试验的轴向位移变形结果。从图 9中可以看出，纯钙质砂在剪切开始时出现小幅度剪缩，随着剪切位移不断增加，出现明显的剪胀现象，最大剪胀量为-0.438 mm。由图 9（a）可知，当含石英粉的混合砂中细颗粒含量分别为10%，20%，30%，40%时，最大剪缩量为0.181，0.203，0.26，0.536 mm。表明随着细颗粒含量的增加，混合砂的变形特征逐渐由剪胀状态转变为剪缩状态。由图 9（b）可知，当含玻璃珠混合砂中细颗粒含量为10%，20%，30%，40%时，最大剪缩量为分别为0.261，0.331，0.421，0.604 mm。表明同一围压下随着细颗粒含量的增加，试样的剪胀性逐渐减小，剪缩状态逐渐明显，当细颗粒含量为10%~30%时，两种不同细颗粒的混合砂的轴向变形都较为接近，相差幅度较小，当细颗粒含量增加至40%时，其轴向变形的增长趋势明显增加，且含玻璃珠的混合砂的剪缩量比含石英粉的混合砂具有明显的增长趋势，且剪缩特征更加明显。
为研究不同轴向压力条件下混合砂的轴向变形情况，对细颗粒含量为20%和40%的混合砂的轴向变形进行汇总分析，如图 10所示。由图 10（a）可知，细颗粒含量为20%的混合砂在轴向压力为100 kPa时，轴向位移先逐渐增大，随后发生小幅度的减小，最后趋于稳定，两种混合砂的最大剪缩量分别为0.203，0.262 mm；当轴向压力为200 kPa时，混合砂均呈现出明显的剪缩特征，剪缩量随剪切位移增加而增加，两种混合砂的最大剪缩量分别为0.41，0.662 mm；轴向压力为300 kPa时，混合砂均呈现出较大的剪缩特征，随着剪切位移的增加，两种混合砂的最大剪缩量分别为0.723，0.834 mm。如图 10（b）所示，细颗粒含量为40%时混合砂均表现出明显的剪缩特征。随着剪切位移的增加，剪缩量呈现出不断上升的趋势，且随着轴向压力的增加，两种混合砂的最大剪缩量分别达到0.536，0.604，0.718，0.843，0.974，1.181 mm。表明随着轴向压力不断增加，混合砂的最大剪缩量逐渐增加，呈现出明显的剪缩特征，是因为随着轴向压力不断增加，试样在剪切过程中发生较大程度的颗粒破碎现象，大颗粒逐渐减少，致使混合砂的剪缩现象显著。
![]()
图 10 细颗粒含量为20%和40%的钙质混合砂的轴向变形−剪切位移关系曲线Figure 10. Axial deformation-shear displacement curves for calcareous mixed sand with fine particles contents of 20% and 40%3. 颗粒破碎特征
钙质混合砂特殊的剪切特征主要是由于其颗粒特征和颗粒破碎等多重因素共同作用的结果，为了对不同细颗粒含量及形状条件下混合砂的剪切过程进行细观分析，需建立细颗粒含量及形状与颗粒破碎之间的联系。对剪切后试样进行颗粒级配试验，得到剪切前后各组试样的颗粒级配曲线，采用Hardin[27]所提出的相对破碎指标Br进行颗粒破碎的定量化评价，相对颗粒破碎率的计算原理如图 11所示，其计算公式如下：
Br=BtBp0。 (4) 式中：Bp0为初始破碎势，即试样试验前级配曲线的积分；Bt为总破碎势，即试验前后的整体破碎势之差。
图 12为不同细颗粒含量及形状条件下钙质混合砂的破碎关系图。由图 12可知，各组混合砂的相对破碎率均小于3%。在轴向压力为100，200，300 kPa条件下，试样的最大破碎率分别为0.9%，1.3%，2.4%。表明随轴向压力的不断增加，更多的钙质砂颗粒发生破碎，混合砂的颗粒破碎度呈现出不断上升的趋势。在细颗粒含量0~40%的增长过程中，当轴向应力为100 kPa时，两种混合砂的相对破碎率由0.9%分别减小至0.41%和0.29%；当轴向应力为200 kPa时，两种混合砂的相对破碎率由1.3%分别减小至0.87%和0.78%；当轴向应力为300 kPa时，两种混合砂的相对破碎率由2.4%分别减小至1.88%和1.78%，该现象是因为细颗粒含量的增加降低了钙质砂之间的咬合程度，颗粒破碎率也随之减小。与含玻璃珠混合砂的相对破碎率相比，含石英粉混合砂的相对破碎率更高，二者差值在0.07%~0.16%，其主要原因是石英粉的颗粒形状与玻璃珠相比更为不规则，棱角度较高，对于钙质砂颗粒之间的咬合降低程度较小，致使含石英粉的钙质混合砂有着较高的颗粒破碎率。
![]()
图 12 不同细颗粒含量及形状条件下钙质混合砂的破碎关系图Figure 12. Crushing relationship of calcareous mixed sand under different contents and shapes of fine particles图 13为钙质混合砂的相对破碎率随组合形状参数变化趋势图。由图 13可知，混合砂的相对破碎率随着细颗粒含量的增加不断减小，与组合圆度及组合宽长比均表现出线性负相关关系，表明当混合砂整体的组合形状参数越偏离球形颗粒，颗粒轮廓呈现枝状时，其破碎程度越高。细颗粒含量的增加使得混合砂整体的组合形状参数向球形颗粒靠近，包围在钙质砂周围细颗粒增加，减少了钙质砂颗粒之间的接触，起到润滑作用，使得颗粒破碎程度不断降低。玻璃珠含量的不断增加，表现出更为明显的润滑作用，且混合砂整体的组合形状参数更加趋近球形颗粒的形状参数，含玻璃珠混合砂的颗粒破碎程度明显较小。
![]()
图 13 相对破碎率与组合形状参数关系图Figure 13. Relationship between relative crushing rate and composite shape parameters4. 结论
通过对不同细颗粒形状及含量条件下的钙质混合砂进行环剪试验，分析细颗粒形状及含量间的差异对钙质混合砂力学特性及破碎特性的影响，得到以下5点结论。
（1）钙质混合砂的应力-位移曲线均呈现软化特性，含圆形玻璃珠混合砂的峰值强度明显小于含角形石英粉混合砂的峰值强度，随着细颗粒含量不断增加，混合砂的峰值强度逐渐降低，且含玻璃珠混合砂峰值强度有着更大的降低趋势。
（2）钙质混合砂的强度指标随细颗粒含量、组合圆度、组合宽长比增加均呈现减小趋势。由于玻璃珠表面圆滑，致使含玻璃珠混合砂表现出更低的内摩擦角和表观黏聚力。
（3）钙质混合砂的软化系数随轴向压力的增加先减小后增大，在200 kPa时达到最小值，并随细颗粒含量增加表现出先增大后减小的变化趋势，当细颗粒含量为10%时，软化系数最大，当细颗粒含量超过10%后，混合砂的软化系数逐渐减小，且含石英粉试样的软化系数明显低于含玻璃珠试样。
（4）纯钙质砂表现出明显的剪胀，且随细颗粒含量和轴向压力增加，混合砂的剪缩量逐渐增大，当细颗粒含量达到40%时呈现出明显的剪缩，含玻璃珠混合砂的剪缩增量更大。
（5）在不同轴压下，钙质混合砂的相对破碎率与组合圆度及组合宽长比均呈线性负相关，且轴压的增加会造成相对破碎率的增加。随细颗粒含量的增加，混合砂组合形状参数趋于球形形状参数，颗粒破碎程度减小，含玻璃珠混合砂具有更小的相对破碎率。
-
![]()
图 1 冻融循环条件下不同含水率红砂岩T2谱分布情况
Figure 1. T2 distribution of red sandstone samples subjected to freeze-thaw cycles under different water contents
![]()
图 3 冻融后不同含水率红砂岩全过程蠕变曲线
Figure 3. Creep curves of red sandstone after freeze-thaw cycles under different water contents
![]()
图 4 冻融后不同含水率红砂岩蠕变速率及加速度曲线
Figure 4. Creep rate and acceleration curves of red sandstone after freeze-thaw cycles under different water contents
![]()
图 5 冻融后不同含水率红砂岩蠕变破坏形态
Figure 5. Creep failure modes of red sandstone after freeze-thaw cycles under different water contents
![]()
图 7 黏弹性冻融含水劣化系数回归曲线
Figure 7. Regression curves of viscoelastic freeze-thaw deterioration coefficient considering water content
![]()
图 8 理论计算及模拟监测曲线与试验数据对比
Figure 8. Comparison among theoretical calculation, numerical simulation curve and experimental data
表 1 冻融后不同含水率红砂岩长期强度
Table 1 Long-term strengths of red sandstone after freeze-thaw cycles under different water contents
冻融循环次数 不同含水率红砂岩的长期强度/MPa 干燥 1% 2% 3% 饱和 0 124.69 124.03 104.81 102.48 101.32 30 113.90 89.73 85.11 65.51 66.02 60 75.02 55.71 53.87 48.67 28.32 
下载: 导出CSV
表 2 非屈服应力条件下（40.8 MPa）模型辨识参数
Table 2 Parameters identified by universal global algorithm under non-yield stress (40.8 MPa)
冻融循环次数/N 含水率/% Gh/GPa K/GPa Gk/GPa ηk/(GPa·h-1) 0 0 43.36 37.16 13.42 2.41 1 37.91 32.50 10.17 2.13 2 36.63 31.39 9.09 1.86 3 22.46 19.25 8.34 1.70 4.63 18.98 16.27 7.28 1.30 30 0 13.67 11.72 6.63 1.21 1 11.42 9.79 5.62 1.17 2 8.31 7.13 5.11 1.04 3 6.87 5.89 4.85 1.02 4.63 6.43 5.51 4.52 0.83 60 0 8.10 6.94 5.42 0.98 1 5.61 4.80 3.90 0.80 2 4.32 3.70 3.30 0.65 3 4.08 3.50 2.91 0.60 4.63 3.99 3.42 2.74 0.43
下载: 导出CSV
表 3 屈服应力条件下模型辨识参数
Table 3 Parameters identified by universal global optimization algorithm under yield stress
Gh/GPa K/GPa Gk/GPa ηk/(GPa·h-1) ηp/(GPa·h-1) θ λ 4.73 4.65 13.31 3.46 4588.61 39.19 5.01 3.09 3.03 3.88 0.78 1429.21 108.35 9.45 注： 冻融循环30次,含水率为2%,屈服应力为102.0 MPa；冻融循环60次,含水率为3%,屈服应力为61.2 MPa。
下载: 导出CSV
-
[1] 薛翊国, 孔凡猛, 杨为民, 等. 川藏铁路沿线主要不良地质条件与工程地质问题[J]. 岩石力学与工程学报, 2020, 39(3): 445-468. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202003003.htm XUE Yi-guo, KONG Fan-meng, YANG Wei-min, et al. Main unfavorable geological conditions and engineering geological problems along Sichuan-Tibet railway[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(3): 445-468. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202003003.htm
[2] 陈卫忠, 谭贤君, 于洪丹, 等. 低温及冻融环境下岩体热、水、力特性研究进展与思考[J]. 岩石力学与工程学报, 2011, 30(7): 1318-1336. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201107005.htm CHEN Wei-zhong, TAN Xian-jun, YU Hong-dan, et al. Advance and review on thermo- hydro-mechanical characteristics of rock mass under condition of low temperature and freeze-thaw cycles[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(7): 1318-1336. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201107005.htm
[3] PARK J, HYUN C U, PARK H D. Changes in microstructure and physical properties of rocks caused by artificial freeze-thaw action[J]. Bulletin of Engineering Geology and the Environment, 2015, 74(2): 555-565. doi: 10.1007/s10064-014-0630-8
[4] KHANLARI G, ABDILOR Y. Influence of wet-dry, freeze-thaw, and heat-cool cycles on the physical and mechanical properties of upper red sandstones in central Iran[J]. Bulletin of Engineering Geology and the Environment, 2015, 74(4): 1287-1300. doi: 10.1007/s10064-014-0691-8
[5] 贾海梁, 项伟, 谭龙, 等. 砂岩冻融损伤机制的理论分析和试验验证[J]. 岩石力学与工程学报, 2016, 35(5): 879-895. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201605002.htm JIA Hai-liang, XIANG Wei, TAN Long, et al. Theoretical analysis and experimental verifications of frost damage mechanism of sandstone[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(5): 879-895. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201605002.htm
[6] 韩铁林, 陈蕴生, 师俊平, 等. 冻融循环下钙质砂岩力学特性及其损伤劣化机制的试验研究[J]. 岩土工程学报, 2016, 38(10): 1802-1812. doi: 10.11779/CJGE201610009 HAN Tie-lin, CHEN Yun-sheng, SHI Jun-ping, et al. Experimental study on mechanical properties and damage degradation mechanism of calcareous sandstone subjected to freeze-thaw cycles[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(10): 1802-1812. (in Chinese) doi: 10.11779/CJGE201610009
[7] CHEN T C, YEUNG M R, MORI N. Effect of water saturation on deterioration of welded tuff due to freeze-thaw action[J]. Cold Regions Science and Technology, 2004, 38(2/3): 127-136.
[8] 王俐. 不同初始含水率红砂岩冻融损伤的试验研究及其机理分析[D]. 武汉: 中国科学院武汉岩土力学研究所, 2006. WANG Li. Experimental Studies and Mechanism Analysis on Different Initial Water-Saturated Red Sandstones Under Condition of Frost and Thaw[D]. Wuhan: Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, 2006. (in Chinese)
[9] 刘海康, 张思渊, 张鑫鑫. 不同初始含水率下砂岩冻融劣化特性试验研究[J]. 科学技术与工程, 2017, 17(26): 322-327. doi: 10.3969/j.issn.1671-1815.2017.26.054 LIU Hai-kang, ZHANG Si-yuan, ZHANG Xin-xin. Experimental study of freeze-thaw deterioration specialty of sandstone in different initial moisture content[J]. Science Technology and Engineering, 2017, 17(26): 322-327. (in Chinese) doi: 10.3969/j.issn.1671-1815.2017.26.054
[10] 贺晶晶, 师俊平. 冻融后不同含水状态砂岩的剪切破坏特性[J]. 岩石力学与工程学报, 2018, 37(6): 1350-1358. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201806004.htm HE Jing-jing, SHI Jun-ping. Shear failure properties of sandstone with different moisture contents after cyclic freezing- thawing[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(6): 1350-1358. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201806004.htm
[11] 王鹏, 许金余, 方新宇, 等. 红砂岩吸水软化及冻融循环力学特性劣化[J]. 岩土力学, 2018, 39(6): 2065-2072. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201806019.htm WANG Peng, XU Jin-yu, FANG Xin-yu, et al. Water softening and freeze-thaw cycling induced decay of red-sandstone[J]. Rock and Soil Mechanics, 2018, 39(6): 2065-2072. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201806019.htm
[12] 陈国庆, 郭帆, 王剑超, 等. 冻融后石英砂岩三轴蠕变特性试验研究[J]. 岩土力学, 2017, 38(增刊1): 203-210. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2017S1030.htm CHEN Guo-qing, GUO Fan, WANG Jian-chao, et al. Experimental study of creep properties of quartz sandstone after freezing-thawing cycles[J]. Rock and Soil Mechanics, 2017, 38(S1): 203-210. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2017S1030.htm
[13] 童庆闯. 冻融循环作用下炭质页岩蠕变损伤特性研究[D]. 长沙: 长沙理工大学, 2017. TONG Qing-chuang. Study on Creep Damage Characteristics of Carbonaceous Shale under Freezing and Thawing[D]. Changsha: Changsha University of Science & Technology, 2017. (in Chinese)
[14] 杨秀荣, 姜谙男, 王善勇, 等. 冻融循环条件下片麻岩蠕变特性试验研究[J]. 岩土力学, 2019, 40(11): 4331-4340. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201911024.htm YANG Xiu-rong, JIANG An-nan, WANG Shan-yong, et al. Experimental study on creep characteristics of gneiss under freeze-thaw cycles[J]. Rock and Soil Mechanics, 2019, 40(11): 4331-4340. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201911024.htm
[15] 张峰瑞, 姜谙男, 杨秀荣, 等. 冻融循环下花岗岩剪切蠕变试验与模型研究[J]. 岩土力学, 2020, 41(2): 509-519. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202002022.htm ZHANG Feng-rui, JIANG An-nan, YANG Xiu-rong, et al. Experimental and model research on shear creep of granite under freeze-thaw cycles[J]. Rock and Soil Mechanics, 2020, 41(2): 509-519. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202002022.htm
[16] 杨圣奇, 徐卫亚, 杨松林. 龙滩水电站泥板岩剪切流变力学特性研究[J]. 岩土力学, 2007, 28(5): 895-902. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200705008.htm YANG Sheng-qi, XU Wei-ya, YANG Song-lin. Investigation on shear rheological mechanical properties of shale in Longtan Hydropower Project[J]. Rock and Soil Mechanics, 2007, 28(5): 895-902. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200705008.htm
[17] 王来贵, 赵娜, 何峰, 等. 岩石蠕变损伤模型及其稳定性分析[J]. 煤炭学报, 2009, 34(1): 64-68. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB200901014.htm WANG Lai-gui, ZHAO Na, HE Feng, et al. Rock creep damage model and its stability analysis[J]. Journal of China Coal Society, 2009, 34(1): 64-68. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB200901014.htm
[18] 杨圣奇, 徐鹏. 一种新的岩石非线性流变损伤模型研究[J]. 岩土工程学报, 2014, 36(10): 1846-1854. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201410016.htm YANG Sheng-qi, XU Peng. A new nonlinear rheological damage model for rock[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(10): 1846-1854. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201410016.htm
[19] 赵延林, 唐劲舟, 付成成, 等. 岩石黏弹塑性应变分离的流变试验与蠕变损伤模型[J]. 岩石力学与工程学报, 2016, 35(7): 1297-1308. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201607001.htm ZHAO Yan-lin, TANG Jin-zhou, FU Cheng-cheng, et al. Rheological test of separation between viscoelastic-plastic strains and creep damage model[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(7): 1297-1308. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201607001.htm
[20] 申艳军, 杨更社, 荣腾龙, 等. 岩石冻融循环试验建议性方案探讨[J]. 岩土工程学报, 2016, 38(10): 1775-1782. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201610006.htm SHEN Yan-jun, YANG Geng-she, RONG Teng-long, et al. Proposed scheme for freeze-thaw cycle tests on rock[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(10): 1775-1782. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201610006.htm
[21] 武东生, 孟陆波, 李天斌, 等. 灰岩三轴高温后效流变特性及长期强度研究[J]. 岩土力学, 2016, 37(增刊1): 183-191. WU Dong-sheng, MENG Lu-bo, LI Tian-bin, et al. Study of triaxial rheological property and long-term strength of limestone after high temperature[J]. Rock and Soil Mechanics, 2016, 37(S1): 183-191. (in Chinese)
[22] 许宏发. 软岩强度和弹模的时间效应[J]. 岩石力学与工程学报, 1997, 16(3): 246-251. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX703.006.htm XU Hong-fa. Time-dependent behaviours of strength and elasticity modulus of weak rock[J]. Chinese Journal of Rock Mechanical and Engineering, 1997, 16(3): 246-251. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX703.006.htm
-
期刊类型引用(7)
1. 刘志，杨阳. 黄原胶联合微生物加固改善混合砂蓄水特性试验研究. 土木与环境工程学报(中英文). 2025(03): 49-57 . 
百度学术
2. 梁越，冉裕星，许彬，张鑫强，何慧汝. 细颗粒含量影响渗流侵蚀规律的细观机理研究. 岩土工程学报. 2025(05): 1099-1106 . 本站查看
3. 高霞，梅雯艳，吴强，崔祥龙. 球度和扁平度对含瓦斯水合物煤体宏细观力学性质影响研究. 中国安全生产科学技术. 2025(05): 86-94 . 百度学术
4. 蔡国庆，刁显锋，杨芮，王北辰，高帅，刘韬. 基于CFD-DEM的流-固耦合数值建模方法研究进展. 哈尔滨工业大学学报. 2024(01): 17-32 . 百度学术
5. 程建毅，郭晓军，李泳. 泥石流物源土体标度分布参数与粘聚力的关系. 山地学报. 2024(03): 401-410 . 百度学术
6. 孙增春，刘汉龙，肖杨. 砂-粉混合料的分数阶塑性本构模型. 岩土工程学报. 2024(08): 1596-1604 . 本站查看
7. 孙建强. 颗粒形状对黄土干密度与抗剪强度的影响. 黑龙江工程学院学报. 2024(06): 8-15+23 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 291
- HTML全文浏览量: 55
- PDF下载量: 174
- 被引次数: 12
