冰川冻土, 2022, 44(6): 1820-1832 doi: 10.7522/j.issn.1000-0240.2022.0158

寒区工程与灾害

冻融循环作用下砂土粒组变化的分形特征及其与长期强度的相关性分析

刘启,1,2, 张泽,1,2,3, 张圣嵘1,2, 恽晴飞4, 付峻松1,2

1.东北林业大学 土木工程学院/寒区科学与工程研究院/冻土研究所,黑龙江 哈尔滨 150040

2.东北林业大学 东北多年冻土区地质 环境系统教育部野外科学观测研究站/东北多年冻土区环境、道路建设与养护协同创新中心,黑龙江 哈尔滨 150040

3.中国科学院 西北生态环境资源研究院 冻土工程国家重点实验室,甘肃 兰州 730000

4.上海长凯岩土工程有限公司,上海 200093

Fractal characteristics of sand grain group changes under freeze-thaw cycles and its correlation with long-term strength

LIU Qi,1,2, ZHANG Ze,1,2,3, ZHANG Shengrong1,2, YUN Qingfei4, FU Junsong1,2

1.School of Civil Engineering / Institute of Cold Region Science and Engineering / Permafrost Institute, Northeast Forestry University, Harbin 150040, China

2.Northeast-China Observatory and Research-Station of Permafrost Geo-Environment of the Ministry of Education / Collaborative Innovation Centre for Permafrost Environment and Road Construction and Maintenance in Northeast China (CIC-PERCM), Northeast Forestry University, Harbin 150040, China

3.State Key Laboratory of Frozen Soil Engineering, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China

4.Shanghai Changkai Geotechnical Engineering Co. , Shanghai 200093, China

通讯作者: 张泽,教授,主要从事冻土工程与环境研究. E-mail: zez@nefu.edu.cn

收稿日期: 2022-05-03   修回日期: 2022-09-01  

基金资助: 国家自然科学基金委员会(NSFC)与俄罗斯基础研究基金会(RFBR)合作交流项目.  42011530083
国家自然科学基金项目.  41771078
黑龙江省交通投资集团有限公司科研项目.  JT-100000-ZC-FW-2021-0129
黑龙江省重点研发计划项目(指导类).  GZ20220095

Received: 2022-05-03   Revised: 2022-09-01  

作者简介 About authors

刘启,硕士研究生,主要从事冻土力学与冻土工程研究.E-mail:qil@nefu.edu.cn , E-mail:qil@nefu.edu.cn

摘要

在季节冻土区施工时,冷暖交替所产生的冻融循环作用是寒区工程破坏的主要原因之一,还原土体受冻融作用时的性质和状态尤为重要。制备了不同冻融循环次数的砂土试样,使用球形模板压入仪测试冻结砂土的长期强度,利用分形理论以及显著性分析的方法研究了不同冻融循环次数下两种冻结砂土的粒组与长期强度的变化规律。结果表明:对于细砂(FS),分形维数DB 与长期强度变化呈极显著正相关,其中≥0.15~0.20 mm和≥0.25~0.40 mm粒组与长期强度的拟合优度较好,为FS的优势粒级;对于中砂(MS),分形维数DB 与长期强度呈微正相关,变化呈“竖向N”趋势波动,其中≥0.30~0.40 mm和≥0.40 mm粒组与长期强度的拟合优度较好,为MS的优势粒级;其他粒组含量与长期强度变化并无明显关联。冻融作用改变了土体中粗、细颗粒粒组的含量配比。同时,随着整体粒径区间的增大,优势粒级也随之增大,表现为冻结砂土的长期强度会随着某些粒级含量的增加呈现先降低后升高的趋势。研究结果可为受冻融作用地区长期强度的测定提供理论参考。

关键词: 冻融循环 ; 长期强度 ; 冻结砂土 ; 粒组 ; 分形维数

Abstract

Seasonally frozen soils are widely distributed in China in terms of area, and the freeze-thaw cycle effect generated by the alternation of cold and warmth is one of the causes of engineering damage in cold areas during construction, and it is particularly important to restore the nature and state of the soil when it is subjected to freeze-thaw action. Therefore, sandy soil specimens with different numbers of freeze-thaw cycles were prepared, and the long-term strength of frozen sandy soil was tested using a spherical template indenter. Using fractal theory and the microstructure image processing software ImageJ, the change law of grain group and long-term strength of two frozen sandy soils under different numbers of freeze-thaw cycles were studied. The results show that: for fine sand (FS), the fractal dimension DB has a highly significant positive correlation with the long-term strength variation, among which ≥0.15~0.20 mm and ≥0.25~0.40 mm have the best fit with the long-term strength, and are the dominant grain classes of FS. For medium sand (MS), the fractal dimension DB is slightly positively correlated with the long-term strength, and the variation shows a “vertical N” trend, in which the grain size content of ≥0.30~0.40 mm and ≥0.40 mm fits better with the long-term strength, and is the dominant grain class for MS. The content of other grain groups did not correlate significantly with the long-term strength change. The freeze-thaw action changed the content ratios of coarse and fine grain agglomerates in the soil. With the increase of the overall particle size interval, the dominant particle size also increases, which shows that the long-term strength of frozen sandy soil tends to decrease and then increase with the increase of the content of some particle sizes. The results of the study can provide theoretical reference for the determination of long-term strength in areas subject to freeze-thaw action.

Keywords: freeze-thaw cycles ; long-term strength ; frozen sandy soil ; grain group ; fractal dimension

PDF (4101KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘启, 张泽, 张圣嵘, 恽晴飞, 付峻松. 冻融循环作用下砂土粒组变化的分形特征及其与长期强度的相关性分析[J]. 冰川冻土, 2022, 44(6): 1820-1832 doi:10.7522/j.issn.1000-0240.2022.0158

LIU Qi, ZHANG Ze, ZHANG Shengrong, YUN Qingfei, FU Junsong. Fractal characteristics of sand grain group changes under freeze-thaw cycles and its correlation with long-term strength[J]. Journal of Glaciology and Geocryology, 2022, 44(6): 1820-1832 doi:10.7522/j.issn.1000-0240.2022.0158

0 引言

冻土是指温度在0 ℃及以下,并含有冰的各种岩石和土壤1。根据存在时间,地球上主要分布两类冻土,即多年冻土和季节冻土。在我国,多年冻土主要分布在青藏高原、祁连山以及东北等地区2。季节冻土在我国分布更加广泛,季节冻土区由于受到四季交替,冷暖变化引起的反复冻融作用,路基会出现冻胀、沉陷等一系列灾害。可以将冻融循环理解为一种特殊的强风化作用3,导致土体物理力学性质均在不同程度上发生改变。大量学者研究发现冻融循环作用会改变土体结构4-5、模量6、孔隙率7等与土体长期强度密切相关的因素。因此在受冻融作用影响的地区建设工程设施时,尤其有必要考虑土体的长期强度变化规律问题8-9

Mandelbrot10提出了分形的概念,分形理论的出现为研究分布复杂且无规则的事物提供了一个有力的工具。已有大量学者证明岩石和土体具有良好的分形特性11-14。以往的学者针对各类土体进行了研究,取得了一些规律性的成果。例如,Li等14通过对钙质砂表面分维以及质量分维的研究发现,对于颗粒破碎的定量描述,分形维数比形状描述符更具有统计规律性。陈鑫等15基于不同的分形模型,利用压汞法,定量分析了冻融循环对黄土孔隙的改变规律。于钱米16通过建立增长路径法(GP),推导出粗粒土颗粒破碎的极限二维分形模型和极限分形维数,并计算了粗粒土的极限分形维数约为2.585。王展等17使用人工控温冻融的方法,从冻融次数和分形维数入手,揭示了冻融循环对土微团聚体稳定性改变的影响机制。周泓等18则借助球形模板压入仪探究了冻融循环作用对黄土长期强度的劣化机制。刘建坤等19通过不固结不排水剪切试验,研究了冻融作用下细粒土的不均匀分布对粗粒土力学特性的影响。张泽等20通过对冰碛亚黏土的反复融冻试验,发现粒度成分呈双向性变化,大颗粒发生分裂,小颗粒发生团聚,粒径呈现均一化发展。卜建清等21通过室内三轴试验研究,发现在经历冻融循环后,细粒土的冻胀削弱了粗粒土的强度,并建议以6次冻融循环后的力学指标作为参考值。

综上所述,在已有的成果中很少涉及土体在经历冻融循环作用后,内部的粒组与长期强度的变化规律。因此本文以冻结砂土为研究对象,使用球形模板压入仪获得了不同冻融循环次数下冻结砂土的长期强度,利用分形理论和显著性分析的方法,探究冻结砂土的不同粒组与长期强度的变化趋势,以期得到不同冻融次数下粒度含量对冻结砂土长期强度的影响规律。

1 试验方法

1.1 试样制备

本文以两种不同粒度范围的砂土为研究对象,如图1所示的试验方案设计图,其中C代表冻融循环,例如C4表示样品经历了4次冻融循环。试样进行指定次数的冻融循环试验后进行球形模板压入试验。

图1

图1   试验方案

Fig. 1   The test scheme


根据《公路土工试验规程》(JTG E40—2007)的制备要求,参与试验的砂土均过2 mm筛,颗粒级配曲线如图2所示,并依据该规程以及叙述简便,将土样相对粒度范围较小的称为细砂(FS),反之则称为中砂(MS)。试验土样的干密度为2.65×103 kg·m-3,均为饱水状态。标准环刀和样品的底面内径为61.8 mm,高为20.0 mm。样品采用制样机制备,试样制备完成后,放入饱和缸内进行抽真空饱和,同时用保鲜膜与宽胶带进行保存密封。试验样品的冻结、融化温度分别为-20 ℃与20 ℃,在冻融循环箱中进行样品的冻融,冻融过程无外界水源补给。样品采用多向冻融,即从试样外围向中心冻结。试验进行前,在所设定的冻结和融化温度下,对单个样品进行了反复的冻结和融化试验,发现当冻结和融化时间设置为2 h时,试样均能完成充分的冻结和融化,因此将样品的冻结与融化时间均设置为2 h。为模拟季节冻土区的冻融作用对土体的影响,将冻融循环次数N设为4、6、10、50、100次,对达到要求的冻融样品进行球形模板压入试验。其中冻土强度降低最急剧的阶段是在加载初期,加载8 h得到的黏聚力仅为超极限长期黏聚力的0.10~0.35倍,加载16 h后冻土的黏聚力已不再发生变化,此时的黏聚力便可视为冻土的长期黏聚力22。所以为了实用,极限长期黏聚力可近似地根据短期的,例如8 h后的试验资料,同时引入修正系数0.75来确定(不同土壤其系数有很大差别),亦或是16 h后的试验数据,将其直接视为冻土的长期黏聚力。本文为了进一步降低试验误差以及数据记录的便捷性,设定球模仪压入时间为24 h,并采用24 h时所计算的黏聚力视为本次试验土样的长期强度。为消除土体表面因颗粒不均匀而导致测试结果的离散性问题,对每个冻融循环次数同时制取6个平行样品,并取测试结果平均值作为参考。

图2

图2   颗粒级配曲线

Fig. 2   Particle size distribution curves of the soil


1.2 长期强度测定

球形模板压入仪(球模仪)是一种能够快速测出土体(冻结、未冻的黏土和冻砂等)黏聚力的试验装置(图3),测试原理与布氏硬度仪十分相似23。球模仪将土体视作理想黏结的不密实体,根据黏聚力与固定恒载下球形压在土体中沉陷深度的反比例关系进行推算。若固定恒载没有超出长期强度的极限值,则变形速度随着时间的累积而衰减,最终变形逐渐趋于稳定。

图3

图3   球模仪结构示意图

Fig. 3   Structure sketch of spherical template indenter


冻土的抗剪强度由黏聚力和内摩擦力两部分组成,并在一定条件下满足库仑定律23

τ=cθ,ω,t+σtanφθ,ω,t

式中:cθ,ω,t为冻土的黏聚力(MPa);σ为法向应力(MPa);φθ,ω,t为冻土的内摩擦角(°);θ为温度(℃);ω为含水率(%);t为时间(s)。

其中cφ根据温度、含水率、作用时间来确定。将冻结的黏性土视为理想黏塑性体,此时tanφθ,ω,t→0,式(1)可简化为

τ=cθ,ω,t

即冻土的抗剪强度实际上等于黏聚力,从塑性理论精确解答中推导其黏聚力的大小18

ct=KPπdSt

式中:ct 为随时间变化的黏聚力(MPa);P为作用在球形压头上的固定恒载(kN);K为修正系数,表示压板与土体的硬度比,塑性材料取0.18;d为球形压头的直径(mm);St 为球形压头在土中的沉入深度(mm)。

在固定恒载的作用下,球形压头沉入土中,用沉入深度St 来计算冻土的黏聚力ct。黏聚力ct 与压头接触面的最大接触压力存在一定联系,表述为ct =γq0,本研究中球形压头直径选用为22 mm24。由于受到内摩擦力的影响,将黏聚力ct 视为等效黏聚力,其中包含了内摩擦力的影响。如图4所示,随着时间的增加,球形压头沉入土中的深度受试验土样与球形压头接触面积的制约。球形压头与土体接触面积的增加,导致传入土中单位面积上的应力减少,从而加速了变形的稳定,所以利用球模仪测试土的长期强度非常节省时间,但由于得到的冻土黏聚力为等效黏聚力,是黏聚力和内摩擦系数的综合反映,而不是常见的摩尔-库仑准则中所用的单一黏聚力,因此所测的黏聚力通常要偏大23

图4

图4   球模仪球形压头压入土中示意图

Fig. 4   Sketch of ball indenter pushing into soil of spherical template indenter


使冻融试验箱的温度控制在-20 ℃,试验开始前(图3)首先检查螺栓、螺帽的松紧度,确保球模仪的正常使用。检查完成后,放置试样,使球形压头与试样上表面刚好接触,保持球模仪固定栓的松紧度。使用砝码施加竖向固定恒载,采用YHD-50型位移传感器测量球形压头的沉入深度,保持位移计的纵向竖直度,并将位移计归零。待温度稳定后,扭开球模仪的固定栓,开始进行等效黏聚力的测试,并利用计算机跟踪记录不同时刻下的沉降量。

1.3 土的分形特征

土是一种典型的多孔介质,从而导致传统的欧氏几何很难去描述15,分形理论的创立让描述土体不规则的特性有了新的发展,其主要通过分形维数来定量表征事物的复杂程度,其中分形维数包括Hausdorff维数DH 、计盒维数DB 等,由于计盒维数DB 在数学上的便捷性以及物理上的直观性,成为了分形领域普遍采用的分形维数25-26。因此,本文利用PIP 9.1型颗粒图像处理仪所生成的颗粒显微图像,来分别测算经历不同冻融循环次数下土颗粒的计盒维数DB,从而进一步分析其变化规律。

图5所示的图像处理流程,首先利用PIP 9.1型颗粒图像处理仪将经历不同冻融循环次数下土的颗粒显微图像保存记录,使用ImageJ图像处理软件进行图像类型转换、阈值调整、二值化处理等步骤。对于任意灰度图片其阈值范围为0~255,为选取一个合适的阈值将试样与其他部分分离,经过大量的测试后本文的阈值选取255。

图5

图5   图像处理流程

Fig. 5   Image processing flow


图6所示,经过二值化处理后,可以清楚地看到颗粒的形态特征与分布情况,图中黑色区域即代表颗粒。用任意尺寸r>0的小方块(小盒子)将处理后的二值化图像进行网格划分,其中Nr)代表黑色像素所占有的网格个数。

图6

图6   二值化图像示意图

Fig. 6   Schematic diagram of binarized images: raw image of PIP 9.1 particle analyzer (a) and image after binarization (b)


若存在数D,使得当r→0时

Nrr-D

且仅存在唯一正值b,使得

Nr=br-D

对上式两边取对数,可得

lgNr=lgb+lg1/rD

进一步化简为

lgNr=-Dlgr+lgb

式中:D为分形维数,即所求的计盒维数DB

在实际问题处理过程中,首先计算出覆盖土颗粒的网格数(盒子数)Nr),由式(7)可知,以lg r为横坐标,lg Nr)为纵坐标,在坐标系上使用最小二乘法对散点图进行拟合,利用拟合直线斜率来计算DB。需要指出的是,本文选取计盒方法所得的分形维数是仅具有统计意义上的分形维数,不同于颗粒形状与材质的严格意义上的分形维数。

2 结果与分析

2.1 长期强度的测试结果

使用球形模板压入仪进行冻土的长期强度测试非常节省时间,但由于土颗粒不均匀,球形模板压入仪的球形压头在压入时的测试结果出现很大的离散性23,因此进行多次重复试验是必不可少的。在经历4、6、10、50、100次冻融循环后砂土的长期强度变化特征如图7所示,可以看到FS与MS的等效黏聚力均随着时间的增加而缓慢减小。

图7

图7   冻结砂土的长期强度变化趋势

Fig. 7   Trend of long-term strength of frozen sandy soil


冻融循环作用使土颗粒之间产生较大的应力,颗粒反复发生分裂与团聚,并最终达到一种稳定状态20。值得注意的是,尽管土在宏观上表现出膨胀变形,但内部仍然进行着复杂的收缩膨胀现象,土样的干密度也发生变化,土颗粒的性质发生改变,这种宏观变化也是微观结构变化的一种表征27。可以看到在200 min后,砂土的长期强度降低速度减缓,达到稳定状态,冻融前后的黏聚力降低接近50%左右,这与其他的学者28的研究结果相近。在此过程FS与MS表现出了冻结砂土的在蠕变过程中的塑性、黏性、弹性等性质。在竖向荷载的作用下,FS与MS在进行24 h的测试后,认为此时试样的内部达到了稳定状态,测试的结果可以作为砂土的长期强度参考值22。如图8(a)所示,FS在冻融循环10次时的长期强度达到最大值,6、100次时与10次接近,稍低于10次。50次时的长期强度处于中间水平,4次时的长期强度最低。在图8(b)中,MS同样在冻融循环10次时的长期强度达到最大值,在100次时的长期强度低于10次,处在第二级别。冻融循环4、6次时的长期强度接近,均低于100次,冻融循环50次时的长期强度最低。长期强度是土体内部综合性质的表征,从两种变化曲线上来看,在经历10次冻融循环后,砂土的长期强度均达到最大值,说明10次循环后,试样内部的结构达到了一种稳定状态,抵抗外部荷载的作用有所增强。

图8

图8   长期强度随冻融循环次数的变化

Fig. 8   Variations of long-term strength with number of freeze-thaw cycles


2.2 不同冻融循环次数的分形特征

DB 的计算结果如表1所示,确定系数R2均在0.99以上,说明试样具有良好的分形特性。FS与MS的DB 分别在1.3~1.7之间和1.4~1.8之间,平均值分别为1.478和1.594。

表1   不同冻融循环次数下砂土的分形维数

Table 1  Fractal dimension of sandy soils under different number of freeze-thaw cycles

土样冻融循环次数拟合方程确定系数R2DB
FS4y=-1.337x+11.6570.991.337
6y=-1.441x+12.0530.991.441
10y=-1.635x+12.9700.991.635
50y=-1.408x+12.9640.991.408
100y=-1.571x+12.5620.991.571
MS4y=-1.491x+12.3650.991.491
6y=-1.595x+12.6500.991.595
10y=-1.707x+12.0680.991.707
50y=-1.638x+12.9140.991.638
100y=-1.540x+12.5410.991.540

新窗口打开| 下载CSV


图9所示的DB 随冻融次数的变化规律图,可以看到FS与MS的DB 在10次冻融循环以内随着冻融次数增加呈增大趋势(图9中a的位置),并在10次冻融后DB 达到最大值。随着冻融次数的继续增加,MS的DB 开始持续降低,并在10次至100次冻融循环中呈现线性趋势。FS的DB 在10次至50次冻融循环过程中与MS趋势相同,但降低幅度(图9中b的位置)远高于MS,在50次至100次冻融时,DB 呈现增大趋势,100次冻融循环后,数值高于MS(图9中c的位置)。随着颗粒破碎程度的增加,不同粒径的颗粒含量有所增加,试样所占有的像素比例增大,基于二值化图像所计算的DB 也会增大。在反复的冻融后,颗粒有破碎集聚现象,不同粒级的混合程度降低,DB 也呈现降低趋势。

图9

图9   分形维数随冻融循环次数的变化

Fig. 9   Variations of fractal dimension with number of freeze-thaw cycles


3 长期强度的影响因素

3.1 DB 与长期强度的相关性分析

为进一步探究冻结砂土长期强度随DB 的变化规律,将不同冻融循环次数下的长期强度与DB 做相关性分析。如图10所示,可以看到MS的DB 与长期强度确定系数R2为0.085,呈现微正相关。FS的DB 与长期强度确定系数R2为0.917,呈现极显著正相关。DB 的大小可以集中地反映颗粒的分布分形特征,从图10(a)中可以看出,随着DB 的增大,FS的长期强度迅速增加,当DB 达到最大值时,其长期强度同样达到最大值,可以利用FS的长期强度与DB 呈现高度相关的趋势来推测其长期强度的变化规律。图10(b)中所示MS的长期强度变化在DB 的增大过程中呈波动状态,第一阶段为强度上升阶段,强度增加0.036 MPa;第二阶段为强度下降阶段,降低0.104 MPa;第三阶段与第一阶段趋势相同,但幅度明显高于一阶段,强度增加0.166 MPa。MS的长期强度随着DB 的增大并无明显的线性关系,总体上随着DB 的增大呈现“竖向N”似的波动变化。

图10

图10   分形维数与长期强度的关系

Fig. 10   Relationship between fractal dimension and long-term strength


3.2 粒组含量的变化

DB 的计算方法可以看出,长期强度与不同粒径的土颗粒含量密切相关。正如土质学家谢尔盖耶夫所言,粒度成分是决定土质最重要的因素,已有学者21对经历不同冻融循环次数以及不同细粒含量下的土颗粒的强度进行了分析研究,但并没有涉及土体在经历冻融循环作用后,内部的粒组与长期强度的变化规律。为进一步探究不同冻融循环次数下土体的粒组随长期强度的变化规律,将FS与MS的颗粒尺寸分布分别以不同界限值划分为不同粒组,相关参数见表2~3

表2   FS的粒径分组

Table 2  Grain size grouping of FS

冻融循 环次数粒径分组/mm不均匀 系数Cu曲率系数Cc体积平均粒径/mm颗数平均粒径/mm
<0.05≥0.05~0.10≥0.10~0.15≥0.15~0.20≥0.20~0.25≥0.25~0.40≥0.40
FP1FP2FP3FP4FP5FP6FP7
41.783.318.5346.9826.7212.6701.770.950.170.034
60.464.6619.6619.0318.0738.1201.850.980.200.058
102.6416.3119.4813.2026.9821.3902.441.080.170.038
501.284.246.1037.3934.2316.3301.931.310.190.035
1000.971.842.1214.1819.6924.3036.891.761.240.290.032

注:FP1表示土样FS的第1粒组(PSF1),以此类推。

新窗口打开| 下载CSV


表3   MS的粒径分组

Table 3  Grain size grouping of MS

冻融循环 次数粒径分组/mm不均匀系数Cu曲率系数Cc体积平均 粒径/mm颗数平均 粒径/mm
<0.10≥0.10~0.20≥0.20~0.30≥0.30~0.40≥0.40
MP1MP2MP3MP4MP5
45.4210.008.9832.9042.702.661.650.290.032
63.8411.7416.5224.1643.752.081.170.300.061
100.474.2817.7020.9856.561.820.930.340.089
501.002.507.6017.5071.391.931.010.380.034
1003.2310.6921.8832.6231.601.711.180.280.029

注:MP1表示土样MS的第1粒组(PSM1),以此类推。

新窗口打开| 下载CSV


图11(a)可知,随着冻融循环次数的增加FP1的粒组含量变化呈波动趋势,整体含量由冻融4次的1.78%下降至冻融100次的0.97%。FP2与FP3两个粒组含量变化呈现先增加后降低趋势,分别在冻融10次以及6次时达到峰值16.31%和19.68%。FP4的粒组含量总体上呈现阻尼形态波动,经历4次冻融循环后达到最大含量46.98%,10次后含量达到最低值13.20%,冻融100次后含量与10次相接近。

图11

图11   粒组含量随冻融循环次数的变化

Fig. 11   Variations of coarse grain content with number of freeze-thaw cycles


这是由于在冻融作用的影响下试样内部发生了破碎集聚的现象。FP5的粒组含量在6次冻融循环后达到最低值18.70%,继续冻融含量持续增加,在50次后含量达到最大值34.23%。与之相反,FP6粒组含量在6次冻融循环后达到峰值38.12%,8次冻融循环后含量变化相对稳定,含量最大差值为7.97%,总体上与FP5的粒组变化呈共轭现象。FP7的粒组含量仅在土样进行100次冻融循环后出现。如图11(b)所示,MP1、MP2和MP3的粒组含量在100次的冻融循环过程中,含量最大差值为4.95%、9.42%以及14.28%,粒组区间整体呈现扩大趋势。MP4的粒组含量在冻融循环过程中持续降低,在50次时达到最低值17.50%,之后继续冻融,含量有所回升。MP5的粒组含量与MP4的粒组含量变化也呈现共轭现象,在50次冻融循环后含量达到最大值56.56%,随着冻融次数的增加,颗粒破碎,含量降至最低。

3.3 粒度成分对长期强度的影响规律

为进一步分析不同粒组含量对冻结砂土长期强度的影响程度,将冻结砂土的长期强度与不同粒组含量绘制散点图,并利用最小二乘法进行拟合做相关性分析,结果如图12~13所示。FS中FP1、FP2与FP3的粒组与长期强度的确定系数R2分别为0.0208、0.1787和0.1995,均呈现微正相关,结果表明[图12(a)~(c)]长期强度与这三种粒径级别的含量均没有明显的变化关联。当FP4粒组由13.20%变化至46.98%时,长期强度下降0.0402 MPa[图12(d)],两者确定系数R2为0.9322,呈极显著负相关。FP5粒组变化与长期强度没有较强关联,总体呈下降趋势。随着粒径组别的增大,当FP6粒组由12.67%变化至38.12%时,长期强度呈现明显的升高趋势,长期强度提升0.0394 MPa,确定系数R2为0.5627[图12(e)]。由此可见,FP4与FP5粒径的百分含量对FS的长期强度起到关键性作用,可视为此类砂土的优势粒级。即当FP4与FP6粒组在一定区间内增加时,长期强度分别呈现减小和增大趋势,中间的FP5粒径级别视为过渡区域,含量变化与长期强度并无明显直接关联。

图12

图12   FS不同粒组与长期强度的相关性

Fig. 12   Correlation between different grain group and long-term strength of FS


图13

图13   MS不同粒组与长期强度的相关性

Fig. 13   Correlation between different grain group and long-term strength of MS


对于MS[图13(a)~(c)],MP1、MP2和MP3的粒组和长期强度的确定系数R2分别为0.2744、0.4928和0.4343,呈现微负相关和负相关,长期强度随着这三种粒径级别的含量增加而缓慢降低。当MP4粒组由17.50%变化至32.90%时,长期强度由0.1178 MPa急剧降低至0.0504 MPa。而对于MP5,其变化呈相反趋势,长期强度增加0.1036 MPa。确定系数R2分别为0.8118和0.8703[图13(d)~(e)],呈现极显著负相关与极显著正相关。对于MS而言,在冻融循环的作用下,长期强度的变化与MP4与MP5两种粒组密切相关,视为优势粒级。

4 影响机制分析

比较两种砂土的粒组与长期强度相关性曲线,可以明显发现,MS的优势粒级范围(≥0.30~0.40 mm与≥0.40 mm)明显高于FS的优势粒级范围(≥0.15~0.20 mm与≥0.25~0.40 mm)。可见随着整体粒径范围的增大,长期强度会随着某些粒组在一定范围内的增加而呈现先降低后升高的趋势。这是由于当试验的颗粒构成发生改变时,试样的结构也会随之发生改变。在细颗粒达到界限含量前,试样以粗颗粒间接触以及和冰的接触为主要形式,极少部分的细颗粒以填充和被冰包裹的块状结构形式存在粗颗粒形成的骨架孔隙中,表现为以点接触的粒状堆积结构29。如图14所示,当细颗粒未充分地填充在粗颗粒的孔隙中,其孔隙率为21.298%;当细颗粒充分地填充在粗颗粒的孔隙中,其孔隙率仅为4.913%。

图14

图14   孔隙率随细粒填充的变化趋势

Fig. 14   Variation trend of porosity with fine particle filling


针对细颗粒(FSC)的阈值含量,郭志杰30以最大干密度为评价指标,对粗-细粒混合土进行了一系列的室内试验,发现了随细粒含量变化的阈值为40%与85%。在本次试验中,FS的FP4与FP6以及MS的MP4与MP5均在含量接近40%时强度发生改变,这与其的研究结果相近。对于融土来说,根据细粒含量从低到高,可划分为接触状态1至接触状态4四个层次31。从完全由粗颗粒决定力学特性到完全由细颗粒决定力学特性,即从类粗粒土过渡至类细粒土,接触状态2、3称之为中间性态土。与常规融土不同的是,在本次试验中试样均为饱和的冻结砂土,土颗粒之间由于冰的存在,在受到球形压头的作用后会发生流变现象。如图15所示,在细颗粒达到界限含量(40%)以前,以粗颗粒间的接触以及和冰的接触为主要存在形式,并以粗颗粒和冰作为主要骨架,属于粗粒土的范畴,同时细颗粒会以填充形式存在粗颗粒与冰形成的孔隙之间,并充当“滚珠”的作用32。若细颗粒含量继续增加,由接触状态1转变至接触状态4(图16),也就是由粗粒土范畴转变至细粒土范畴。而细粒土在空间中的填充度要优于粗粒土,此时试样以细粒土间的接触以及和冰的接触为主要形式,使得粗、细粒土间咬合比较紧密,粒间的胶结作用加强30,球形压头压入土体中所需要克服的阻力增加,长期强度表现出增大趋势。并且在一定范围内细颗粒含量的增加,会让细颗粒的“滚珠”作用进一步增强。同时随着整体粒径范围的增大,细颗粒的“相对粒径”也会增大,例如本次试验中FP4与MP4粒级,其含量在一定范围的降低导致土样长期强度的降低。随着粒径的增大,土样重新由粗颗粒和冰占据主导地位,长期强度又随着粗颗粒含量的增加而增大。冻融循环作用改变了粗、细粒土不同粒组的配比以及孔隙冰的变化,使得土颗粒接触状态发生改变,从而导致长期强度也随之发生变化。土体的力学性质受冻融次数与细粒含量的双重影响,对于不同条件下的试验,起主导作用的因素也不同,这也是土体性质受冻融循环作用呈现多样化的原因。

图15

图15   细粒填充变化示意图

Fig. 15   Schematic diagram of variation of fine particle filling


图16

图16   融土的接触状态示意图

Fig. 16   Schematic diagram of contact state of thawed soils


5 结论

在冻融循环的条件下,饱和冻结砂土内部的粒组成分会发生变化,同时影响着冻结砂土的长期强度,考虑冻融循环作用以及细粒含量对寒区工程建设非常重要,通过室内试验改变冻融循环次数来初步探究内部颗粒构成对冻结砂土的长期强度的影响,初步结论如下:

(1)FS与MS的长期强度均随着时间的增加呈现降低的趋势。在10次冻融循环后,冻结砂土的长期强度均达到最大值,FS在4次冻融循环后的长期强度最小,MS在50次冻融循环后的长期强度最小。长期强度是土体内部成分和结构的综合反映,两种砂土在10次冻融循环后长期强度并无明显的剧烈波动特征,若在寒区进行相关的工程设施建设,建议在10次冻融循环后土体稳定时进行。

(2)基于分形理论,计算出FS与MS的DB 在1.4~1.6间波动,并在冻融10次时均达到DB 最大值。基于显著性分析理论,冻融过程中对于长期强度与DB 的关系,FS的拟合优度要大于MS。随着DB 的增大,FS的长期强度呈线性增加变化,MS的长期强度呈“竖向N”似波动变化。

(3)冻融作用使FS与MS的颗粒粒组发生改变。FS中PF4和FP6粒组与长期强度的拟合优度最好,确定系数R2分别为0.93与0.56。MS中MP4和MP5粒组与长期强度的拟合优度最好,确定系数R2分别为0.81与0.87。FS与MS的长期强度变化的“优势粒级”并不相同,但“优势粒级”会随着整体粒径区间的增大而增大,这对于受冻融循环作用的地区预测长期强度具有重要意义。

参考文献

Ma WeiWang Dayan.

Studies on frozen soil mechanics in China in past 50 years and their prospect

[J]. Chinese Journal of Geotechnical Engineering, 2012344): 625-640.

[本文引用: 1]

马巍王大雁.

中国冻土力学研究50 a回顾与展望

[J]. 岩土工程学报, 2012344): 625-640.

[本文引用: 1]

Liu GuiminZhang BoWang Liet al.

Permafrost region and permafrost area in the globe and China

[J/OL]. Earth Science, 20222022-05-03]. .

[本文引用: 1]

刘桂民张博王莉.

全球和我国多年冻土分布范围和实际面积研究进展

[J/OL]. 地球科学, 20222022-05-03]. .

[本文引用: 1]

Roman L TZhang Ze.

Effect of cycles of freezing and thawing on the physical and mechanical properties of moraine loam

[J]. Foundation and Soil Mechanics, 2010473): 96-101.

[本文引用: 1]

Shen JiejieWang QingChen Yatinget al.

Evolution process of the microstructure of saline soil with different compaction degrees during freeze-thaw cycles

[J]. Engineering Geology, 2022304106699.

[本文引用: 1]

Tang LiyunLi GangLuo Taoet al.

Mechanism of shear strength deterioration of soil-rock mixture after freeze-thaw cycles

[J]. Cold Regions Science and Technology, 2022200103585.

[本文引用: 1]

Zhou ZhiweiLi GuoyuShen Mingdeet al.

Dynamic responses of frozen subgrade soil exposed to freeze-thaw cycles

[J]. Soil Dynamics and Earthquake Engineering, 2022152107010.

[本文引用: 1]

Liu BoMa RenmingFan Haoming.

Evaluation of the impact of freeze-thaw cycles on pore structure characteristics of black soil using X-ray computed tomography

[J]. Soil and Tillage Research, 2021206104810.

[本文引用: 1]

Wang YubinHuang ZeQian Jiaxinet al.

Freeze-thaw cycles aggravated the negative effects of moss-biocrusts on hydraulic conductivity in sandy land

[J]. Catena, 2021207105638.

[本文引用: 1]

Li GuoyuWang FeiMa Weiet al.

Variations in strength and deformation of compacted loess exposed to wetting-drying and freeze-thaw cycles

[J]. Cold Regions Science and Technology, 2018151159-167.

[本文引用: 1]

Mandelbrot B B. Fractals: form, chance and dimension[M]. San Francisco, CA, USAW. H. Freeman & Co.1978.

[本文引用: 1]

Tao GaoliangZhang Jiru.

Two categories of fractal models of rock and soil expressing volume and size-distribution of pores and grains

[J]. Chinese Science Bulletin, 20095423): 4458-4467.

[本文引用: 1]

陶高粱张季如.

表征孔隙及颗粒体积与尺度分布的两类岩土体分形模型

[J]. 科学通报, 2009546): 838-846.

[本文引用: 1]

Du JunHou KepengLiang Weiet al.

Experimental study of compaction characteristics and fractal feature in crushing of coarse-grained soils

[J]. Rock and Soil Mechanics, 201334(): 155-161.

杜俊侯克鹏梁维.

粗粒土压实特性及颗粒破碎分形特征试验研究

[J]. 岩土力学, 201334(): 155-161.

Peng RuidongYang YancongJu Yanget al.

Computation of fractal dimension of rock pores based on gray CT images

[J]. Chinese Science Bulletin, 20115631): 3346-3357.

彭瑞东杨彦从鞠杨.

基于灰度CT图像的岩石孔隙分形维数计算

[J]. 科学通报, 20115626): 2256-2266.

Li XueLiu JiankunLi Jinze.

Fractal dimension, particle shape, and particle breakage analysis for calcareous sand

[J]. Bulletin of Engineering Geology and the Environment, 2022813): 106.

[本文引用: 2]

Chen XinZhang ZeLi Dongqing.

Study on the pore features of freezing-thawing loess based on different fractal models

[J]. Journal of Glaciology and Geocryology, 2020424): 1238-1248.

[本文引用: 2]

陈鑫张泽李东庆.

基于不同分形模型的冻融黄土孔隙特征研究

[J]. 冰川冻土, 2020424): 1238-1248.

[本文引用: 2]

Yu Qianmi.

Study on evolutionary regularity of particle breakage coarse-grained soil

[D]. BeijingBeijing Jiaotong University201827-34.

[本文引用: 1]

于钱米.

粗粒土颗粒破碎演化规律研究

[D]. 北京北京交通大学201827-34.

[本文引用: 1]

Wang ZhanZhang YulongYu Naet al.

Effects of freezing-thawing on characteristics and fractal dimension of soil microaggregates

[J]. Acta Pedologica Sinica, 2013501): 83-88.

[本文引用: 1]

王展张玉龙虞娜.

冻融作用对土壤微团聚体特征及分形维数的影响

[J]. 土壤学报, 2013501): 83-88.

[本文引用: 1]

Zhou HongZhang YuchuanZhang Zeet al.

Changing rule of long-term strength of frozen loess cohesion under impact of freeze-thaw cycle

[J]. Rock and Soil Mechanics, 2014358): 2241-2246.

[本文引用: 2]

周泓张豫川张泽.

冻融作用下冻结黄土黏聚力长期强度变化规律

[J]. 岩土力学, 2014358): 2241-2246.

[本文引用: 2]

Liu JiankunYu QianmiLiu Jingyuet al.

Influence of non-uniform distribution of fine soil on mechanical properties of coarse-grained soil

[J]. Chinese Journal of Geotechnical Engineering, 2017393): 562-572.

[本文引用: 1]

刘建坤于钱米刘景宇.

细粒土不均匀分布对粗粒土力学特性的影响

[J]. 岩土工程学报, 2017393): 562-572.

[本文引用: 1]

Zhang ZeMa WeiQi Jilin.

Structure evolution and mechanism of engineering properties change of soils under effect of freeze-thaw cycle

[J]. Journal of Jilin University (Earth Science Edition), 2013436): 1904-1914.

[本文引用: 2]

张泽马巍齐吉琳.

冻融循环作用下土体结构演化规律及其工程性质改变机理

[J]. 吉林大学学报(地球科学版), 2013436): 1904-1914.

[本文引用: 2]

Bu JianqingWang Tianliang.

Influences of freeze-thaw and fines content on mechanical properties of coarse-grained soil

[J]. Chinese Journal of Geotechnical Engineering, 2015374): 608-614.

[本文引用: 2]

卜建清王天亮.

冻融及细粒含量对粗粒土力学性质影响的试验研究

[J]. 岩土工程学报, 2015374): 608-614.

[本文引用: 2]

Mazurov G P. Physical and mechanical properties of permafrost[M]. Liang Huisheng, Wu Qijian, trans. BeijingChina Coal Industry Publishing House1980135-142.

[本文引用: 2]

马祖罗夫 Г П. 冻土物理力学性质[M]. 梁惠生, 伍期建, 译. 北京煤炭工业出版社1980135-142.

[本文引用: 2]

Zhang ZeMa WeiZhang Zhongqionget al.

Application of spherical template indenter to long-term strength tests for frozen soil

[J]. Rock and Soil Mechanics, 20123311): 3516-3520.

[本文引用: 4]

张泽马巍张中琼.

球形模板压入仪在冻土长期强度测试中的应用

[J]. 岩土力学, 20123311): 3516-3520.

[本文引用: 4]

Yun QingfeiXie ChunleiZhang Zeet al.

The solution of the double-sphere model and experimental research of the long-term shear strength of frozen sand based on spherical template indenter test

[J]. Journal of Glaciology and Geocryology, 2022442): 485-494.

[本文引用: 1]

恽晴飞谢春磊张泽.

球模仪测试冻结砂长期抗剪强度的双球模型解答及其试验研究

[J]. 冰川冻土, 2022442): 485-494.

[本文引用: 1]

Ivanovici M.

Fractal dimension of color fractal images with correlated color components

[J]. IEEE Transactions on Image Processing, 2020288069-8082.

[本文引用: 1]

Li Yurong.

Fractal dimension estimation for color texture images

[J]. Journal of Mathematical Imaging and Vision, 2020621): 37-53.

[本文引用: 1]

Xu XuezuWang JiachengZhang Lixin. Physics of frozen soil[M]. BeijingScience Press2001169.

[本文引用: 1]

徐学祖王家澄张立新. 冻土物理学[M]. 北京科学出版社2001169.

[本文引用: 1]

Chang DanLiu JiankunLi Xu.

Experimental study on yielding and strength properties of silty sand under freezing-thawing cycles

[J]. Chinese Journal of Rock Mechanics and Engineering, 2015348): 1721-1728.

[本文引用: 1]

常丹刘建坤李旭.

冻融循环下粉砂土屈服及强度特性的试验研究

[J]. 岩石力学与工程学报, 2015348): 1721-1728.

[本文引用: 1]

Chen YulongZhang Yuning.

Experimental study of effects of non-plastic fines on liquefaction properties of saturated sand

[J]. Rock and Soil Mechanics, 2016372): 507-516.

[本文引用: 1]

陈宇龙张宇宁.

非塑性细粒对饱和砂土液化特性影响的试验研究

[J]. 岩土力学, 2016372): 507-516.

[本文引用: 1]

Guo Zhijie.

Effect of fine soil content on physical and mechanical properties of mixed coarse- and fine-grained soil

[D]. BeijingBeijing Jiaotong University201825-27.

[本文引用: 2]

郭志杰.

细粒含量对粗-细粒混合土物理力学特性的影响

[D]. 北京北京交通大学201825-27.

[本文引用: 2]

Wu QiChen GuoxingZhou Zhenglonget al.

Experimental investigation on liquefaction resistance of fine-coarse-grained soil mixtures based on theory of intergrain contact state

[J]. Chinese Journal of Geotechnical Engineering, 2018403): 475-485.

[本文引用: 1]

吴琪陈国兴周正龙.

基于颗粒接触状态理论的粗细粒混合料液化强度试验研究

[J]. 岩土工程学报, 2018403): 475-485.

[本文引用: 1]

Bi ShengChen GuoxingZhou Zhenglonget al.

Experimental study on influences of fines content and consolidation stress on shear modulus and damping ratio of saturated sand

[J]. Chinese Journal of Geotechnical Engineering, 201739(): 48-52.

[本文引用: 1]

毕昇陈国兴周正龙.

细粒含量及固结应力对饱和砂土动剪切模量和阻尼比影响试验研究

[J]. 岩土工程学报, 201739(): 48-52.

[本文引用: 1]

/