DOI： 10.11817/j.issn.1672-7207.2020.08.013

考虑干密度影响的红黏土蠕变特性及蠕变模型研究

李娜^{1, 2},陈昌富^{1, 2},朱世民^{1, 2},毛凤山^{1, 2}

（1. 湖南大学 建筑安全与节能教育部重点实验室，湖南 长沙，410082；

2. 湖南大学 土木工程学院，湖南 长沙，410082）

摘要:为深入研究密实度对红黏土蠕变特性的影响，制备不同干密度的红黏土试样，经等压固结后，利用改制的三轴蠕变试验装置开展排水三轴蠕变试验，获得不同干密度的红黏土分级加载全过程曲线，并采用“陈氏叠加法”将分级加载全过程曲线转化为分别加载曲线。根据Burgers蠕变模型参数的物理含义，结合部分蠕变试验曲线对模型参数进行辨识，确定不同干密度的红黏土试样Burgers蠕变模型参数；然后，通过确定Burgers蠕变模型参数与干密度及偏应力之间的经验关系，建立考虑干密度影响的Burgers蠕变模型，并利用此蠕变模型对未参与模型参数辨识的蠕变试验曲线进行预测。研究结果表明：本文建立的蠕变模型对不同干密度的红黏土蠕变曲线具有较好的预测效果。

关键词:红黏土；三轴蠕变；干密度；Burgers模型

中图分类号:TU446            文献标志码:A            开放科学(资源服务)标识码(OSID)

文章编号:1672-7207（2020）08-2174-09

Research on creep characteristics and creep model of red clay considering effect of dry density

LI　Na^{1, 2}, CHEN　Changfu^{1, 2}, ZHU　Shimin^{1, 2}, MAO　Fengshan^{1, 2}

(1. Key Laboratory of Building Safety and Energy Efficiency of Ministry of Education, Hunan University, Changsha 410082, China；

2. College of Civil Engineering, Hunan University, Changsha 410082, China)

Abstract: In order to study the influence of compactness on the creep properties of red clay, the red clay samples with different dry densities were prepared, and the triaxial creep test was carried out by using the modified triaxial creep test device after the isostatic consolidation. The whole process curves of red clay with different dry densities were obtained, and the whole process curves of the staged loading were transformed into the corresponding loading curves by using “Chen's loading method”. According to the physical meaning of the parameters of Burgers creep model, the Burgers creep model parameters of red clay samples with different dry densities were determined according to the partial creep test curves. Then, by determining the empirical relationship between the parameters of Burgers creep model and dry density and deviatoric stress, a modified Burgers creep model considering the influence of dry density was established.Finally,the established creep model was used to predict the creep test curves which were not involved in the model parameter identification. The results show that the improved Burgers creep model has good performance in the prediction of red clay creep curves with different dry densities.

Key words: red clay; triaxial creep; dry density; Burgers model

红黏土广泛分布于我国湖广地区、云贵高原以及长江中下游地区，在分布区进行高速公路和铁路建设施工时，需要对红黏土进行夯实或碾压使其达到设计密实度。增大密实度能提高土体的强度^[1-3]，而密实度对红黏土蠕变特性的影响直接关系到红土边坡的长期稳定性以及红土路基的长期沉降量，因此，研究干密度对红黏土蠕变特性的影响具有重要的工程意义。目前,有关干密度对土体蠕变特性影响的研究成果较少。王松鹤等^[4]对不同干密度黄土的进行三轴蠕变试验，发现通过改善排水情况、增大土体密实程度能避免产生蠕变破坏。此外，史祥等^[5]对不同干密度的分散性土进行三轴固结排水试验，发现分散性土干密度愈大，土体变形愈小，蠕变现象愈不明显。干密度对红黏土蠕变特性有很大影响，但用于定量描述该影响的研究较少。用于定量描述土体蠕变特性的蠕变模型主要包括元件模型和经验模型^[6-8]。张先伟等^[9]建立的经验型蠕变模型与Singh-Mitchell模型和Mesri模型相比，其对饱和软土蠕变特性的预测效果更好。陈昌富等^[10]通过锚-土界面剪切蠕变试验建立了能同时描述衰减蠕变和加速蠕变的经验蠕变模型。罗庆姿等^[11]建立了含分数阶导数的Merchant模型，与原Merchant模型相比，其拟合精度更高，更适用于长期变形预测。经验模型虽然简单实用，但拟合得到的模型参数和方程缺乏物理含义，而元件模型中的参数都具有明确的物理意义，且借助不同应力水平下蠕变试验曲线就能辨识元件模型的参数^[12-14]。如金德海等^[15]推导出Burgers模型的三轴蠕变理论解，并通过模型参数辨识建立了一种考虑围压对粗粒土蠕变特性影响的修正Burgers模型。本文作者利用改制的三轴蠕变试验装置，通过红黏土排水三轴蠕变试验，获得不同干密度下红黏土三轴蠕变试验曲线；然后，以试验曲线为依据，辨识Burgers蠕变模型中的参数，得到模型参数随干密度的变化规律，从而建立一种用于描述不同干密度的红黏土蠕变试验曲线的Burgers蠕变模型；最后，利用所提出的模型对试验结果进行预测分析，以验证本文模型的可靠性。

1  红黏土三轴蠕变试验

1.1　试验用料及装置

试验所用红黏土取自湖南省长沙市宁乡县某高速路段，通过XRD测得该红黏土主要矿物成分为石英、高岭石、云母和赤铁矿。由标准击实试验测得红黏土的最大干密度为1.78 g/cm³，最优含水率为18%。土样的其他物理力学参数如表1所示，颗粒级配曲线如图1所示。

表1　红黏土的基本物理力学参数

Table 1　Basic physical mechanical properties of red clay

图1　土的颗粒级配曲线

Fig. 1　Particle distribution of red clay

试验所用装置为自行设计改造的TSZ应变控制式三轴仪，如图2所示。该装置由恒载系统和数据采集系统组成，具体原理见文献[16]。

图2　三轴蠕变试验装置

Fig. 2　Tri-axial creep testing device

1.2　试样制备及试验方案

试验所用试样为圆柱形土样，其直径φ=39.1 mm，高度h=80 mm。试样制备的具体步骤为：首先，配置含水率为18%的土料；然后，计算不同干密度的土样所对应的土料质量，称取设定质量的土料进行分层击实，击实后土样的干密度分别为1.49，1.60，1.65 和1.68 g/cm³；最后，对试样进行抽气饱和，并测得试样饱和后的含水率分别为29.8%，25.1%，23.5%和22.2%，由此制备不同干密度的饱和红黏土蠕变测试样。

试样制备完成后，先让其在200 kPa围压下固结完成。然后，采用分级加载的方式进行蠕变加载，加载等级为8~10级。本文选用的蠕变稳定标准为在10 000 s内试样的轴向变形小于0.01 mm^[17]，由于红黏土强度高，压缩性低，试样在加载后基本能在24 h内达到该稳定标准，故每级荷载的加载时间统一设定为24 h。

1.3　试验结果

根据上述试验方法开展不同干密度的红黏土三轴蠕变试验，得到分级加载全过程曲线。研究表明土体在外力作用下产生的变形包括固结变形和蠕变变形^[16]，由于本文仅研究红黏土的蠕变变形，因而需消除试样在偏应力作用下的固结变形。为此，首先，采用“陈氏叠加法”^[18]将红黏土试样的三轴蠕变分级加载全过程曲线转换为分别加载曲线；然后，利用文献[16]中方法消除固结变形，即可获得试样在分别加载条件下的蠕应变与时间的关系曲线(蠕变曲线)，如图3所示。

图3　红黏土三轴蠕变分别加载曲线

Fig. 3　Triaxial creep loading curves of red clay

2  考虑干密度的红黏土蠕变模型

2.1　研究思路

传统的Burgers模型^[19]是由Maxwell模型和Kelvin模型串联而成的4元件黏弹性流变模型，其具体结构及物理特性如图4(a)所示。传统的Burgers模型因能够较准确地描述各类岩土材料在衰减和稳态蠕变阶段的特性，在岩土工程中应用广泛，为此，本文建立考虑干密度影响的Burgers蠕变模型。

图4　Burgers蠕变模型

Fig. 4　Burgers creep models

Burgers蠕变模型的蠕应变是由串联的弹簧体对应的瞬时弹性应变ε₁、串联的黏壶体对应的等速应变ε₂以及串联的Kelvin模型对应的蠕应变ε₃组成：

(1)

对式(1)求导，获得蠕变速率与时间t的关系：

(2)

下面对模型的边界条件进行分析。

当t=0时，

(3)

当t→∞时，

(4)

由式(3)和式(4)联立可得

(5)

式中：ε₁为瞬时弹性应变；k₁为蠕变曲线在时间t趋于无穷大时对应的斜率；k₂为Kelvin模型对应的起始蠕变率；ε_{3, t→∞}为Kelvin模型对应的最终应变。

参数ε₁，ε_{3, t→∞}，k₁和k₂均可根据实测的蠕变曲线确定，如图4所示。故可通过式(5)确定Burgers模型的4个模型参数。

选取部分蠕变试验曲线进行辨识，以确定参数ε₁，ε_{3, t→∞}，k₁和k₂，进而得到Burgers模型的4个模型参数；对4个模型参数与试样干密度和偏应力之间的关系进行回归分析，由此确定模型参数与试样干密度和偏应力之间的经验关系式；最后，将上述确定的经验关系式代入Burgers蠕变方程，即可建立考虑干密度影响的Burgers蠕变模型。

2.2　模型参数的辨识

下面以干密度为1.68 g/cm³的红黏土试样的三轴蠕变试验结果为例，介绍Burgers模型4个模型参数的辨识方法。

2.2.1　弹性模量E₁的辨识　

由上述分析可知，ε₁为试样加载瞬时所产生的应变，对应蠕变曲线(见图3)在初始时刻的轴向应变。瞬时弹性应变ε₁与偏应力△σ关系如图5所示。从图5可知：瞬时弹性应变ε₁与偏应力△σ大致呈线性关系，采用线性函数进行回归分析，结合式(5)，不难得出该回归直线的斜率为1/E₁。

图5　瞬时弹性应变ε₁与偏应力水平△σ关系曲线

Fig. 5　Correlation between instantaneous elastic strain ε₁ and deviatoric stress △σ

2.2.2　黏壶系数η₁的辨识　

由图3可知：当偏应力△σ较小时，试样产生衰减蠕变，蠕变应变渐趋于恒定，蠕变速率渐趋于0 min^-1，即表明串联的黏壶体对应的等速应变ε₂=0，k₁=0 min^-1；当偏应力△σ超过某一水平(假定为σ₀)时，试样产生等速蠕变，蠕变应变最终将呈线性增加，蠕变速率将渐趋于某一定值，由式(4)可知该定值为k₁；随着偏应力△σ继续增大，试样将产生加速蠕变，并最终产生破坏^[20]。本文不考虑加速蠕变阶段，为描述试样的衰减和等速蠕变特性，引入开关函数，并令

(6)

(7)

式中：为试样产生等速蠕变的偏应力阈值。

由图3可知：当t>600 min时，蠕变曲线近似为直线，由式(4)可知该直线的斜率为k₁。等速蠕变速率k₁与偏应力水平△σ的关系如图6所示。从图6可知k₁与偏应力△σ大致呈线性关系，结合式(4)和式(7)发现该回归直线的斜率为1/η₁。此外，回归直线与偏应力轴的交点即为σ₀，本例中，σ₀=75.1 kPa。

图6　等速蠕变速率k₁与偏应力水平△σ关系曲线

Fig. 6　Correlation between isokinetic creep rate k₁ and deviatoric stress △σ

2.2.3　弹性模量E₂的辨识　

根据式(1)，将蠕变应变中的瞬时弹性应变ε₁和等速蠕变段应变ε₂去除，即得到Kelvin模型对应的蠕变应变ε₃，ε₃随着时间的增大逐渐趋于某一定值ε_{3, t→∞}。统计不同偏应力△σ下Kelvin模型的最终蠕变应变ε_{3, t→∞}如表2所示。由式(4)可知，弹性模量E₂即为偏应力△σ与ε_{3, t→∞}的比值。弹性模量E₂与偏应力水平△σ关系曲线如图7所示。由图7可知E₂与△σ大致呈指数函数关系。采用式(8)进行回归分析，即可建立E₂与△σ的经验关系：

表2　偏应力与Kelvin最终蠕应变

Table 2　Deviatoric stress and Kelvin final creep strain

图7　弹性模量E₂与偏应力水平△σ关系曲线

Fig. 7　Correlation between elastic modulus E₂ and deviatoric stress △σ

(8)

式中：α和β为拟合参数，与试样的干密度有关。

2.2.4　黏壶系数η₂的辨识　

前面已获得Kelvin模型对应的蠕变应变ε₃，由式(4)可知，蠕变应变ε₃在初始时刻的曲线斜率即为k₂。起始蠕变率k₂与偏应力水平△σ的关系如图8所示。从图8可见：k₂与偏应力△σ基本呈线性关系。采用线性函数进行回归分析，由式(5)可知，该回归直线的斜率为1/η₂。

图8　起始蠕变率k₂与偏应力水平△σ关系曲线

Fig. 8　Correlation between initial creep rate k₂ and deviatoric stress △σ

按照上述方法对其他干密度的部分红黏土蠕变曲线进行分析，由此建立Burgers模型参数E₁，η₁，E₂和η₂与偏应力△σ之间的经验关系，并进一步获得用以表征上述经验关系的6个参数E₁，η₁，α，β，η₂和σ₀，见表3。

2.3　考虑干密度影响的蠕变模型的建立

前面已建立Burgers模型参数E₁，η₁，E₂和η₂与偏应力△σ之间的经验关系，并获得了用以表征上述经验关系的6个参数E₁，η₁，α，β，η₂和σ₀。若再建立出这6个参数与试样干密度ρ_d之间的经验关系，则可进一步建立Burgers模型参数E₁，η₁，E₂和η₂与偏应力△σ及干密度ρ_d之间的经验关系，由此可建立考虑干密度影响的红黏土蠕变模型。

模型参数与干密度关系曲线如图9所示。由图9可知：参数E₁，η₁和η₂与干密度ρ_d均大致呈指数函数关系；参数α与干密度ρ_d大致呈线性关系；参数β与干密度ρ_d大致呈抛物线关系。参数σ₀取平均值81 kPa，采用相应的函数形式对参数E₁，η₁，α，β，η₂和σ₀与偏应力△σ之间的关系进行回归分析，建立这6个参数与干密度ρ_d之间的经验关系(见式(9))，再将式(9)代入式(10)即可建立考虑干密度影响的Burgers蠕变模型。

表3　蠕变模型拟合参数值

Table 3　Creep model fitting parameter values

图9　模型参数与干密度关系曲线

Fig. 9　Correlation between model parameters and dry density

(9)

(10)

利用上述建立的考虑干密度影响的Burgers蠕变模型对用于模型参数辨识的红黏土蠕变曲线进行预测分析，结果如图10所示。由图10可知：当偏应力较小时，本文模型对蠕变曲线的预测效果较好；而随着偏应力增加，模型预测精度有所降低。但总体上，本文模型对参与模型参数辨识的红黏土蠕变曲线的预测效果较好。

图10　本文模型对参与参数辨识的蠕变曲线预测分析

Fig. 10　Predictive analysis of creep curves participating in parameter identification

3  预测效果分析

为进一步验证本文提出的考虑干密度影响的Burgers模型的预测效果，对另一部分未参与模型参数辨识的红黏土蠕变曲线进行预测，结果见图11。由图11可知：本文模型对未参与模型参数辨识的红黏土三轴蠕变试验曲线的预测效果也较好。这进一步验证了本文模型的适用性和准确性。

图11　本文模型对未参与参数辨识的蠕变曲线的预测结果

Fig. 11　Prediction results of creep curves of model that are not involved in parameter identification

4  结论

1) 为深入研究干密度对红黏土蠕变特性的影响，利用自行设计改造后的三轴蠕变试验装置开展不同干密度的红黏土三轴蠕变试验，获得不同干密度的红黏土试样的分别加载曲线。

2) 基于Burgers模型4个模型参数的物理意义，对部分蠕变试验曲线进行模型参数辨识，以此建立考虑干密度影响的Burgers蠕变模型。

3) 利用建立的考虑干密度影响的Burgers蠕变模型对蠕变试验曲线进行预测分析，结果表明本文所建模型具有良好的预测效果，适用于描述该地区的红黏土蠕变特性。

4) 重塑试样破坏了红黏土的结构特性，后续可进一步开展红黏土在原状及重塑状态下的蠕变对比试验，研究红黏土的结构性对其蠕变特性的影响。

参考文献：

[1] 宋日英, 黄志全, 陈宇. 密度及水分状态对黄土状粉土变形特征的影响[J]. 工程地质学报, 2011, 19(4): 492-497.

SONG Riying, HUANG Zhiquan, CHEN Yu. Effects of density and moisture content on the deformation characteristics of loess-like silt[J]. Journal of Engineering Geology, 2011, 19(4): 492-497.

[2] 杨果林, 滕珂, 谢兰芳. 膨胀土强度的室内直剪和原位推剪对比试验研究[J]. 中南大学学报(自然科学版), 2014, 45(6): 1952-1959.

YANG Guolin, TENG Ke, XIE Lanfang. Eexperimental study of direct shear and in situ clipper on strength of expansive soil[J]. Journal of Central South University(Science and Technology), 2014, 45(6): 1952-1959.

[3] 宋继宏, 胡明鉴, 付克俭, 等. 宜巴高速岩堆不同密实度大型直剪强度特性[J]. 工程地质学报, 2012, 20(5): 687-692.

SONG Jihong, HU Mingjian, FU Kejian, et al. Strebgth behavior of rock debris heaps with different densities in Yiba expressway[J]. Journal of Engineering Geology, 2012, 20(5): 687-692.

[4] 王松鹤, 骆亚生. 复杂应力下黄土蠕变特性试验研究[J]. 岩土力学, 2009, 30(S2): 43-47.

WANG Songhe, LUO Yasheng. Research on creep characteristics of loess under complex stress[J]. Rock and Soil Mechanics, 2009, 30(S2): 43-47.

[5] 史祥, 樊恒辉, 刘刚, 等. 分散性土蠕变特性试验研究[J]. 岩土力学, 2017, 38(4): 1015-1022.

SHI Xiang, FAN Henghui, LIU Gang, et al. An experimental study of creep properties of dispersive soil[J]. Rock and Soil Mechanics, 2017, 38(4): 1015-1022.

[6] LIINGAARD M, AUGUSTESEN A, LADE P V. Characterization of models for time-dependent behavior of soils[J]. International Journal of Geomechanics, 2004, 4(3): 157-177.

[7] SINGH A, MITCHELL J K. General stress-strain-time function for soils[J]. Geotechnical Special Publication, ASCE, 1968, 94(1): 21-46.

[8] MESRI G, FEBRES-CORDERO E, SHIELDS D R, et al. Shear stress-strain-time behaviour of clays[J]. Géotechnique, 1981, 31(4): 537-552.

[9] 张先伟, 王常明. 饱和软土的经验型蠕变模型[J]. 中南大学学报(自然科学版), 2011, 42(3): 791-796.

ZHANG Xianwei, WANG Changming. Empirical creep model for saturated soft soil[J]. Journal of Central South University(Science and Technology), 2011, 42(3): 791-796.

[10] 陈昌富, 刘俊斌, 徐优林, 等. 锚-土界面剪切蠕变试验及其经验模型研究[J]. 岩土工程学报, 2016, 38(10): 1762-1768.

CHEN Changfu, LIU Junbin, XU Youlin, et al. Tests on shearing creep of anchor-soil interface and its empirical model[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(10): 1762-1768.

[11] 罗庆姿, 陈晓平, 王盛, 等. 软黏土变形时效性的试验及经验模型研究[J]. 岩土力学, 2016, 37(1): 66-75.

LUO Qingzi, CHEN Xiaoping, WANG Sheng, et al. An experimental study of time-dependent deformation behaviour of soft soil and its empirical model[J]. Rock and Soil Mechanics, 2016, 37(1): 66-75.

[12] 夏才初, 孙钧. 蠕变试验中流变模型辨识及参数确定[J]. 同济大学学报(自然科学版), 1996, 24(5): 498-503.

XIA Caichu, SUN Jun. Distinction of rheological model and determination of parameters on creep tests[J]. Journal of Tongji University(Natural Science), 1996, 24(5): 498-503.

[13] 孙钧. 岩石流变力学及其工程应用研究的若干进展[J]. 岩石力学与工程学报, 2007, 26(6): 1081-1106.

SUN Jun. Rock rheological mechanics and its advance in engineering applications[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(6): 1081-1106.

[14] 杨振伟, 金爱兵, 周喻, 等. 伯格斯模型参数调试与岩石蠕变特性颗粒流分析[J]. 岩土力学, 2015, 36(1): 240-248.

YANG Zhenwei, JIN Aibing, ZHOU Yu, et al. Parametric analysis of Burgers model and creep properties of rock with particle flow code[J]. Rock and Soil Mechanics, 2015, 36(1): 240-248.

[15] 金德海, 徐明. 考虑围压对粗粒土蠕变特性影响的Burgers模型参数修正方法研究[J]. 工程力学, 2016, 33(12): 135-142.

JIN Dehai, XU Ming. Parameter correction method of Burgers model for coarse-grained soil considering confining pressure[J]. Engineering Mechanics, 2016, 33(12): 135-142.

[16] 陈昌富, 朱世民, 毛凤山, 等. 红黏土固结-蠕变特性及其耦合模型[J]. 工程地质学报, 2019, 27(4): 723-728.

CHEN Changfu, ZHU Shimin, MAO Fengshan, et al. Characterization and modelling of coupled consolidation-creep behavior of red clay[J]. Journal of Engineering Geology, 2019, 27(4): 723-728.

[17] 孙钧. 岩土材料流变及其工程应用[M]. 北京: 中国建筑工业出版社, 1999: 55-61.

SUN Jun. Rheological behavior of geomaterials and its engineering applications[M]. Beijing: China Architecture and Building Press, 1999: 55-61.

[18] TJONG-KIE T. Determination of the rheological parameters and the hardening coefficients of clays[C]// Rheology and Soil Mechanics Berlin, Heidelberg: Springer Berlin Heidelberg, 1966: 256-272.

[19] MANSOURI H, AJALLOEIAN R. Mechanical behavior of salt rock under uniaxial compression and creep tests[J]. International Journal of Rock Mechanics and Mining Sciences, 2018, 110: 19-27.

[20] 范广勤. 岩土工程流变力学[M]. 北京: 煤炭工业出版社, 1992: 2-4.

FAN Guangqin. Rheological mechanics of geotechnical engineering[M]. Beijing: Coal Industry Press, 1992: 2-4.

（编辑  陈灿华）