高庙子膨润土热膨胀特性试验研究

曹胜飞; 刘月妙; 谢敬礼; 闫安; 高玉峰; 佟强

doi:10.11779/CJGE202202020

高庙子膨润土热膨胀特性试验研究 English Version

曹胜飞^{1, 2,},
刘月妙^{1, 2},
谢敬礼^{1, 2},
闫安^{1, 2},
高玉峰^{1, 2},
佟强^{1, 2}

1.
国家原子能机构高放废物地质处理创新中心, 北京 100029
2.
核工业北京地质研究院中核高放废物地质处置评价技术重点实验室, 北京 100029

基金项目: 国家国防科技工业局项目（地下实验室条件下缓冲材料原位试验安装技术研究）

详细信息

作者简介:
曹胜飞(1983—)，男，博士，高级工程师，主要从事高放废物地质处置缓冲材料研究工作。E-mail：csf831016@163.com

中图分类号: TU443
计量
- 文章访问数: 201
- HTML全文浏览量: 35
- PDF下载量: 127
出版历程
- 收稿日期: 2021-05-25
- 网络出版日期: 2022-09-22
- 刊出日期: 2022-01-31

Experimental study on thermal expansion properties of GMZ bentonite

CAO Sheng-fei^{1, 2,},
LIU Yue-miao^{1, 2},
XIE Jing-li^{1, 2},
YAN An^{1, 2},
GAO Yu-feng^{1, 2},
TONG Qiang^{1, 2}

1.
CAEA Innovation Center on Geologicac Disposal of High Level Radioactive Waste, Beijing 100029, China
2.
CNNC Key Laboratory on Geological Disposal of High Level Radioactive Waste, Beijing Research Institute of Uranium Geology, Beijing 100029, China

摘要

摘要: 采用DIL806热膨胀仪，通过室内试验系统研究了不同影响因素对高庙子压实膨润土样品热膨胀性能的影响，获得了膨润土样品热膨胀系数随干密度、升温速率、含水率和环境气氛等条件的变化规律。试验结果表明：膨润土样品的热膨胀系数随着干密度的增大而增大；样品热膨胀系数随初始含水率的升高而降低；高含水率的膨润土样品，其热膨胀系数随升温速率的升高而增大，含水率为0的试验样品其热膨胀系数随升温速率的升高略有下降；对于同一干密度和含水率的试验样品，空气氛围环境下样品的热膨胀系数要高于氩气氛围环境下的热膨胀系数。研究成果对分析膨润土的热膨胀性能，评价高放废物处置缓冲材料的长期稳定性具有一定的参考依据。
- 高庙子膨润土 /
- 热膨胀系数 /
- 干密度 /
- 含水率 /
- 升温速率 /
- 气氛环境
Abstract: Laboratory tests are conducted on the compacted Gaomiaozi (GMZ) bentonite by using the DIL 806 thermal dilatometer. The influences of different factors on the thermal expansion properties of the bentonite are systematically analyzed. The thermal expansion coefficient of the bentonite is obtained with the dry density, the heating rate, the water content and the atmosphere environment. The experimental results show that the coefficient of thermal expansion increases with the dry density of the bentonite. The coefficient of thermal expansion decreases with the increase of the initial water content. The coefficient of thermal expansion increases with the increase of the heating rate for the samples with high water content, and the coefficient of thermal expansion decreases with the increase of the heating rate for the drying bentonite. Under the same dry density and water content of the samples, the coefficient of thermal expansion is higher in the air environment than that in the Ar gas environment. The research results have certain reference for analyzing the thermal expansion properties of the bentonite and evaluating the long-term stability of buffer materials for high-level radioactive waste disposal.
- GMZ bentonite /
- coefficient of thermal expansion /
- dry density /
- water content /
- heating rate /
- atmosphere environment

HTML全文

0. 引言

盐渍土具有溶陷、冻胀及腐蚀等特殊工程性质，盐渍土地区修建道路在工程上常采用换填土层处理。黑龙江省处于季冻区，路基土受到交通动荷载、冬季低温及冻融循环综合作用，路基土的力学性能和动力特性因此发生改变，主要体现在土体动剪切模量和阻尼比的变化^[1]。目前对于盐渍土的研究主要包括温度变化下水盐迁移运动研究^[2-3]，通过室内试验结果，分析水盐运动主要影响因素，建立温度梯度下盐渍土水–热–盐多场耦合模型，揭示水盐运动规律；土体冻胀和溶陷特性研究，将湿陷和溶陷加以区分，提出峰值溶陷系数以描述不同含盐量对压缩、湿陷和溶陷的影响。考虑冻胀和盐胀综合作用，分析冻胀影响因素^[4-6]；盐渍土改性剂种类和掺量以及改良盐渍土的静力性质研究，采用不同改性剂对盐渍土进行改良，并对改良盐渍土抗剪强度及破坏强度进行测定，根据试验结果，建立考虑多重因素的改良盐渍土破坏强度及抗剪强度预测模型^[7-8]。以上研究多针对未处理盐渍土，或主要考虑改良盐渍土静力性质。若通过石灰改良含盐量较低的盐渍土，使其满足路基力学性能的要求，可以降低成本并节约土壤资源。综上，目前对于季冻区低温及冻融循环作用下的石灰改良盐渍土动力特性研究偏少但其影响不容忽视。

本文以黑龙江省哈肇公路沿线盐渍土为研究对象，经石灰改良后进行循环三轴试验，研究其在不同动荷载频率、围压、温度以及冻融循环作用下的动剪切模量及阻尼比变化规律。并基于Hardin双曲线模型进行曲线拟合，分析石灰改良盐渍土动力参数在不同因素下的响应规律，为工程设计提供动力参数，并为进一步研究不同改性剂种类及不同盐分种类的改良盐渍土动力特性提供理论基础。

1. 试验材料与试验方案

1.1 试验材料

试验土样取自黑龙江省哈肇公路沿线，依据《公路土工试验规程：JTG E40-2007》，测得土样基本物理参数：天然含水率为6.85%，液限为21.2%，塑限为14.2%，最大干密度为1.79 g/cm³，最佳含水率为16%。土样的颗粒级配曲线如图 1所示，由试验结果可知该土样属于低液限黏土，土样中各离子种类及含量如表 1所示。由试验结果可知此土样总含盐量为0.5047%，其中CO₃^2-离子含量最高达0.2205%，夹杂有含量0.1086%的HCO₃^-离子和0.1134%的Cl^-离子，由此判定此土样为含有少量Cl^-的碳酸盐渍土且含盐量较低。

图 1 试验用土颗粒级配曲线

Figure 1. Grain-size distribution curves of test specimens

下载: 全尺寸图片幻灯片

表 1 盐渍土离子含量测定值

Table 1. Measured values of ion content of saline soil (%)

阳离子					阴离子
K⁺	Na⁺	Ca²⁺	Mg²⁺	总量	CO₃^2-	HCO₃^-	SO₄^2-	Cl^-	总量
0.21	0.2173	0.0244	0.0357	0.4870	0.2205	0.1086	0.0622	0.1134	0.5047

下载: 导出CSV

| 显示表格

1.2 试样制备

依据《公路路基设计规范：JTG D30—2015》选择压实度为96%，并查阅相关文献及结合工程实际，确定石灰掺量为6%^[9]。

试样制备采用击实法，烘干土样经0.25 mm土工标准筛筛分，并与蒸馏水充分混合后装入塑料袋中后闷料24 h。按照石灰掺量称取一定质量的石灰均匀拌入盐渍土中，称取所需质量的土样，分3层制成直径39.1 mm，高80 mm的圆柱试样，将其密封后进行恒温恒湿养护28 d。

1.3 试验仪器

试验采用英国GDS公司的动态三轴试验系统（GDSTAS），仪器可以施加最大荷载为10 kN，荷载频率范围为1~5 Hz，最大围压为200 MPa，温度调节范围为-20℃~20℃。冻融循环过程采用DP.CDR-2型快速冻融试验设备。

1.4 试验条件和方案

试验主要控制条件为温度和冻融循环次数。根据相关文献及该地区路基冬季低温温度范围^[10-11]，低温试验温度设置为20℃，0℃，-2℃，-4℃，-6℃，-8℃，-10℃，-16℃共8个温度变化梯度，养护后的试样在设定温度下冻结24 h后，进行温度控制的循环加载试验。冻融循环次数为0，1，3，5，7，9次，冻融循环试验的冻结温度为-16℃，融化温度为20℃。

研究动荷载频率及围压对动力特性影响时采用未冻融试样，试验温度为20℃。荷载参数^[12-14]参考公路交通荷载设置频率为1，3，5 Hz，围压为50，100，150 kPa。试验方案见表 2。

表 2 试验方案

Table 2. Test plan

试验编号	试验温度 /℃	频率/Hz	有效围压 /kPa	冻融循环次数
S1	20	1	100	0
S2~S3	20	3，5	100	0
S4~S6	20	1	50，100，150	0
S7~S13	0，-2，-4，-6， -8，-10，-16	1	100	0
S14~S18	20	1	100	1，3，5，7，9

下载: 导出CSV

| 显示表格

养护完成的试样进行试验前，需在动三轴仪器上进行CO₂排气，水头差饱和、多级反压饱和，当B值大于0.95时结束饱和。试样等压固结后，施加分级正弦波荷载，固结过程及动荷载加载过程见图 2。进行低温加载试验时，养护后的试样经饱和后需按照设定温度冻结24 h，并预先将试验仪器温度控制在设定温度后再进行试验。终止加载条件为轴向变形达到5%或破坏。

图 2 正弦形式轴向动荷载示意图

Figure 2. Schematic diagram of sinusoidal axial dynamic load

下载: 全尺寸图片幻灯片

2. 试验结果与分析

2.1 不同频率下动力参数的变化

未冻融，围压100 kPa，温度20℃，不同频率石灰改良盐渍土动剪切模量变化如图 3（a）所示。3种频率下，随着剪应变的增大，动剪切模量呈衰减趋势，曲线形式为上凸形。剪应变较小时，动剪切模量衰减速度较慢。剪应变较大时，动剪切模量衰减速度加快。在动剪切模量–剪应变曲线变化过程中存在转折点，Lin等^[15]认为转折点动剪切模量可以评价不同因素对于动剪切模量的影响。动荷载频率由1 Hz增加到5 Hz，转折点动剪切模量增加168%。

图 3 不同频率下动剪切模量及阻尼比拟合曲线

Figure 3. Fitting curves of dynamic shear modulus and damping ratio at different frequencies

下载: 全尺寸图片幻灯片

不同频率下阻尼比变化情况如图 3（b）所示，3种频率下，随着剪应变的增大，阻尼比逐渐增加，较低剪应变水平时阻尼比增长缓慢，较高剪应变水平时阻尼比迅速增长，频率增大引起阻尼比降低。以剪应变最小时对应的阻尼比为最小阻尼比，当动荷载频率由1 Hz增长到5 Hz，最小阻尼比下降0.06。

对不同频率下的动剪切模量–剪应变，阻尼比–剪应变曲线按照Hardin等^[16-17]的双曲线模型进行拟合，拟合公式为

${\tau _{\text{d}}} = \frac{{{\gamma _{\text{d}}}}}{{{a} + b{\gamma _{\rm{d}}}}} \text{，}$

(1)

${{G} _{\text{d}}}{\text{ = }}\frac{{{\tau _{\text{d}}}}}{{{\gamma _{\rm{d}}}}}{\text{ = }}\frac{1}{{{a} + b{\gamma _{\rm{d}}}}} \text{，}$

(2)

$\lambda = {\lambda _{\rm{max}}}{\left( {1 - \frac{{G} }{{{{G} _{\rm{max}}}}}} \right)^{n} } = {\lambda _{\rm{max}}}{\left( {1 - \frac{1}{{1 + \frac{{{\gamma _{\text{d}}}}}{c}}}} \right)^n} 。$

(3)

式中 ${\tau _{\text{d}}}$ 为动应力； ${\lambda _{\rm{max}}}$ 为最大阻尼比；a，b，c，n为拟合参数，拟合参数如表 3所示，拟合曲线如图 3所示。

表 3 不同频率及围压下动剪切模量及阻尼比的拟合参数

Table 3. Fitting parameters of dynamic shear modulus and damping ratio under different frequencies and confining pressures

频率/Hz	动剪切模量G_d			阻尼比 $\lambda$
频率/Hz	a	b	R²	${\lambda _{\max }}$	c	n	R²
1	0.0101	4.7873	0.9865	0.46	0.0025	0.6845	0.9183
3	0.0076	3.4389	0.9886	0.39	0.0096	0.3665	0.9893
5	0.0056	3.3365	0.9887	0.36	0.0169	0.3469	0.9950

围压/kPa	动剪切模量G_d			阻尼比λ
围压/kPa	a	b	R²	λ_max	c	n	R²
50	0.0145	4.7686	0.9942	0.55	0.0036	0.4693	0.9940
100	0.0101	4.7873	0.9865	0.46	0.0025	0.6845	0.9183
150	0.0080	3.4256	0.9941	0.43	0.1306	0.2455	0.9802

下载: 导出CSV

| 显示表格

当式（2）中 ${\gamma _{\text{d}}}$ →0时

${{G} _{{\text{dmax}}}} = \frac{1}{{a} } \text{；}$

(4)

当式（2）中 ${\gamma _{\text{d}}}$ →+∞时，得到

${\tau _{{\text{dult}}}} = \frac{1}{{b} } 。$

(5)

式中，G_dmax为最大动剪切模量， ${\tau _{{\text{dult}}}}$ 为最终动应力幅值。由式（3）可以得到 ${\lambda _{\rm{max}}}$ ，由式（4）可以得出不同温度下的G_dmax。不同频率下G_dmax及 ${\lambda _{\rm{max}}}$ 如表 4所示。

表 4 不同频率及围压下最大动剪切模量及阻尼比

Table 4. Maximum dynamic shear modulus and damping ratios under different frequencies and confining pressures

参数	频率/Hz			围压/kPa
参数	1	3	5	50	100	150
G_dmax/MPa	98.81	130.89	177.31	69.11	98.81	125.31
λ_max	0.46	0.39	0.36	0.55	0.46	0.42

下载: 导出CSV

| 显示表格

曲线变化的原因与试样中土粒重新排列有关，较低剪应变水平时，石灰与土体形成稳定结构使得土体变形较小，动剪切模量和阻尼比变化缓慢。较高剪应变水平时，石灰与土体形成的稳定结构逐渐破坏，土粒重新排列，黏滞性水膜在荷载作用下发生流动，产生较大变形，动剪切模量和阻尼比迅速变化。对于不同频率而言，频率增大使得加载速率提高，相同加载周次中塑性变形积累速度变慢，达到同样塑性变形所需动应力增加。同时，试样中黏滞性水膜流动和重新排列效应减弱，土粒重组时间减少，塑性变形相应减弱，塑性变形耗能降低，表现出动剪切模量增加，阻尼比降低的特性。

2.2 不同围压下动力参数的变化

未冻融，20℃，1Hz加载的石灰改良盐渍土动剪切模量在不同围压下的变化情况如图 4（a）所示，不同围压曲线变化趋势与频率相同，围压增加，转折点动剪切模量增加，阻尼比下降。转折点动剪切模量增加84%，最小阻尼比降低0.05。围压增加，试样中土颗粒发生位移，孔隙得到填充，土体结构性得到增强，强度得到提升。不同围压下拟合曲线如图 4所示，拟合参数如表 3所示，不同围压下G_dmax及 ${\lambda _{\rm{max}}}$ 如表 4所示。

图 4 不同围压下动剪切模量及阻尼比拟合曲线

Figure 4. Fitting curves of dynamic shear modulus and damping ratio under different confining pressures

下载: 全尺寸图片幻灯片

由不同频率及围压下石灰改良盐渍土动剪切模量和阻尼比的拟合参数 ${\lambda _{\rm{max}}}$ ，得到最大阻尼比，并通过式（4）计算不同频率及围压下石灰改良盐渍土最大动剪切模量，如表 4所示。

2.3 温度对动力参数影响及温度修正系数

不同温度下的动剪切模量–剪应变，阻尼比–剪应变曲线同样按照式（1）~（3）进行拟合，拟合参数如表 5所示，拟合曲线见图 5。

表 5 不同温度下动剪切模量及阻尼比的拟合参数

Table 5. Fitting parameters of dynamic shear modulus and damping ratio under different temperatures

温度/℃	动剪切模量G_d			阻尼比 $\lambda$
温度/℃	a	b	R²	${\lambda _{\max }}$	c	n	R²
20	0.0101	4.7873	0.9865	0.4546	0.0025	0.6845	0.9183
0	0.0073	4.1682	0.9965	0.3928	0.0018	0.5971	0.9384
-2	0.0041	3.2572	0.9728	0.3718	0.0016	0.6851	0.9707
-4	0.0031	2.7513	0.9908	0.2818	0.0015	0.6153	0.9330
-6	0.0028	2.5054	0.9989	0.2170	0.0013	0.5750	0.9742
-8	0.0024	2.5351	0.9966	0.1915	0.0010	0.6831	0.9409
-10	0.0024	2.1925	0.9811	0.1441	0.0008	0.5857	0.9024
-16	0.0024	1.7242	0.9247	0.0700	0.0006	0.4172	0.9567

下载: 导出CSV

| 显示表格

图 5 不同温度下动剪切模量及阻尼比拟合曲线

Figure 5. Fitting curves of dynamic shear modulus and damping ratio under different temperatures

下载: 全尺寸图片幻灯片

由不同温度下石灰改良盐渍土动剪切模量和阻尼比的拟合参数 ${\lambda _{\rm{max}}}$ 得到最大阻尼比，并通过式（4）计算不同温度下石灰改良盐渍土最大动剪切模量，如表 6所示。

表 6 不同温度下最大动剪切模量及阻尼比

Table 6. Maximum dynamic shear modulus and damping ratios under different temperatures

温度/℃	G_dmax/MPa	${\lambda _{\rm{max}}}$
20	98.81	0.46
0	137.17	0.39
-2	242.13	0.37
-4	320.51	0.28
-6	357.14	0.22
-8	413.22	0.19
-10	414.94	0.14
-16	423.73	0.07

下载: 导出CSV

| 显示表格

随着温度的降低，动剪切模量增加，阻尼比降低。温度由20℃降低至-16℃，转折点动剪切模量增加263 MPa（75→338 MPa），阻尼比减小0.10。随着温度降低，土中水分不断结冰，土体结构更加致密使黏结力增大。抗动荷载能力及抗变形能力增加，动剪切模量增大，达到相同动应力幅值所需能量降低，阻尼比减小。

根据不同温度下动力参数试验结果，提出温度修正系数：

${G'_{\rm{dmax}}} = {\alpha _{\rm{G}}}{{G} _{\rm{dmax}}} \text{，}$

(6)

${\lambda '_{\rm{dmax}}} = {\alpha _\lambda }{\lambda _{\rm{max}}} 。$

(7)

式中， ${\alpha _{\rm{G}}}$ 为动剪切模量的温度修正系数， ${\alpha _\lambda }$ 为阻尼比的温度修正系数。基于孙锐等^[18]对于负温下的黏土动力参数预测公式，对 ${\alpha _G}$ 和 ${\alpha _\lambda }$ 进行拟合分析：

${\alpha _G} = {\alpha _1}\left( {\frac{{{{A} _1} - {{A} _2}}}{{1 + {{\text{e}}^{(T - {T_{\text{0}}})/{r_1}}}}} + {{A} _2}} \right) + {\alpha _2}({{y} _0} + {{A} _0}{e^{{{ - {T} } \mathord{\left/ {\vphantom {{ - {T} } {{{r} _2}}}} \right. } {{{r} _2}}}}}) \text{，}$

(8)

${\alpha _\lambda } = {y_0} + {A'_0}{{\text{e}}^{{T \mathord{\left/ {\vphantom {T {{{r'}_2}}}} \right. } {{{r'}_2}}}}} 。$

(9)

拟合参数如表 7所示，拟合曲线如图 6所示。得到温度修正系数表达式为

$\begin{array}{c} {\alpha _G} = 0.53\left( {\frac{{2.76 + 0.29}}{{1 + {{\text{e}}^{{{\left( {T + 0.76} \right)} \mathord{\left/ {\vphantom {{\left( {T + 0.76} \right)} {2.26}}} \right. } {2.26}}}}}} - 0.29} \right) + \\ 0.66(4.48-3.12{\text{e}}^{-T/-4.32})\text{ }\text{，} \end{array}$

(10)

${\alpha _\lambda } = - 0.019 + 0.91{{\text{e}}^{{T \mathord{\left/ {\vphantom {T {11.42}}} \right. } {11.42}}}} 。$

(11)

表 7 温度修正系数的拟合参数

Table 7. Fitting parameters of temperature correction coefficient

A₁	A₂	T₀	r₁	A₀	y₀	r₂
2.76	-0.29	-0.76	2.26	-4.32	4.48	-3.12

${\alpha _1}$	${\alpha _2}$	R²	${y'_0}$	${A'_0}$	${r'_2}$	R²
0.53	0.66	0.991	-0.019	0.91	11.42	0.990

下载: 导出CSV

| 显示表格

图 6 动剪切模量和阻尼比温度修正系数拟合曲线

Figure 6. Fitting curves of temperature correction coefficient for dynamic shear modulus and damping ratio

下载: 全尺寸图片幻灯片

由图 6可知，温度修正系数转折点位于-6℃附近，T < -6℃时，温度修正系数下降缓慢；T > -6℃时，温度修正系数迅速下降。

2.4 冻融循环作用对动力参数影响及冻融修正系数

在不同冻融循环次数下，动剪切模量–剪应变、阻尼比–剪应变曲线按式（1）~（3）拟合，a₂，b₂，c₂，n₂为拟合参数，拟合参数如表 8所示，拟合曲线如图 7所示。

表 8 不同冻融循环次数动剪切模量及阻尼比的拟合参数

Table 8. Fitting parameters of dynamic shear modulus and damping ratio under different freeze-thaw cycles

冻融循环次数	动剪切模量G_d			阻尼比λ
冻融循环次数	a₂	b₂	R²	${\lambda _{\rm{max}}}$	c₂	n₂	R²
0	0.0102	4.7873	0.9865	0.4649	0.0025	0.6972	0.8884
1	0.0145	6.1150	0.9960	0.5675	0.0040	0.7006	0.9551
3	0.0214	6.1155	0.9984	0.6641	0.0049	0.6692	0.9631
5	0.0227	6.9325	0.9951	0.6827	0.0061	0.7094	0.9558
7	0.0262	7.8786	0.9953	0.7022	0.0075	0.7610	0.9514
9	0.0283	7.6030	0.9983	0.7124	0.0085	1.0048	0.9741

下载: 导出CSV

| 显示表格

图 7 不同冻融循环次数下动剪切模量及阻尼比拟合曲线

Figure 7. Fitting curves of dynamic shear modulus and damping ratio under different freeze-thaw cycles

下载: 全尺寸图片幻灯片

由不同冻融循环次数下石灰改良盐渍土动剪切模量和阻尼比的拟合参数 ${\lambda _{\rm{max}}}$ ，得到最大阻尼比，并通过式（4）计算不同冻融循环次数下石灰改良盐渍土最大动剪切模量，如表 9所示。

表 9 不同冻融循环次数下最大动剪切模量及阻尼比

Table 9. Maximum dynamic shear modulus and damping ratios under different freeze-thaw cycles

冻融循环次数	G_dmax/MPa	${\lambda _{\rm{max}}}$

0	98.81	0.47
1	70.08	0.57
3	46.53	0.66
5	43.80	0.68
7	39.26	0.70
9	35.65	0.71

下载: 导出CSV

| 显示表格

由图 7（a）可知，动剪切模量在冻融循环作用后显著下降，尤其是第3次冻融后，动剪切模量下降40 MPa（75→35 MPa），降幅约为47%。当冻融循环次数超过3次后，动剪切模量衰减趋于稳定。第9次冻融循环作用后，较高剪应变水平时，试件的动剪切模量衰减至较低数值，强度几乎损失殆尽。随着冻融循环作用次数的增加，阻尼比不断增大，3次冻融循环作用后，阻尼比变化不明显。

土体冻结伴随着盐胀作用，土体中孔隙增加。融化后土体中的孔隙仍然存在，石灰胶结土体的稳定结构受到一定程度的破坏，导致土体结构较冻融前更加松散。在动荷载作用下土体塑性变形增加，动剪切模量大幅降低。达到同样动应力幅值需要的能量增加使得阻尼比提升。随着冻融循环作用次数的增加，特别是冻融循环次数超过3次后，试样土体结构中孔隙大幅增加，石灰胶结土体的稳定结构遭到破坏。同时试样土体内的楔形力不断增长，内部的应力累计超过土体自身黏结力导致土体结构破坏，但由于没有来自外界的水分补充，这种破坏效应逐渐停止，动剪切模量降低到较低状态后趋于稳定。

根据不同冻融循环次数下动力参数试验结果，提出冻融修正系数：

${G'_{\rm{dmax}}} = {\beta _G}{{G} _{\rm{dmax}}} \text{，}$

(12)

${\lambda '_{\max }} = {\beta _\lambda }{\lambda _{\max }} \text{，}$

(13)

式中， ${\beta _{\rm{G}}}$ 为动模量的冻融修正系数， ${\beta _\lambda }$ 为阻尼比的冻融修正系数。

根据严晗等^[19]、张向东等^[20]提出的冻融修正系数公式，对本试验的冻融修正系数进行拟合分析，拟合公式为

${\beta _G} = {{a} _3}{{b} _3}{{\text{e}}^{\frac{{ - N}}{{{c_3}}}}} \text{，}$

(14)

${\beta _\lambda } = {{a} _4} - {{b} _4}{e^{{c_4}N}} 。$

(15)

拟合参数如表 10所示，拟合曲线如图 8所示。

表 10 冻融循环修正系数的拟合参数

Table 10. Fitting parameters of freeze-thaw cycle correction coefficient

${a_3}$	${b_3}$	${c_3}$	${R^2}$
0.56	1.87	0.38	0.991

${a_4}$	${b_4}$	${c_4}$	${R^2}$
1.53	-0.50	0.65	0.980

下载: 导出CSV

| 显示表格

图 8 动剪切模量和阻尼比冻融修正系数拟合曲线

Figure 8. Fitting curves of freeze-thaw correction coefficient for dynamic shear modulus and damping ratio

下载: 全尺寸图片幻灯片

得到冻融修正系数：

${\beta _G} = 0.56 - 1.87{{\text{e}}^{\frac{{ - N}}{{0.38}}}} \text{，}$

(16)

${\beta _\lambda } = 1.53 - 0.50{{\text{e}}^{0.65N}} 。$

(17)

由图 8可知，冻融转折点为第3次冻融循环，第3次冻融循环之前，曲线变化急促，当冻融循环次数大于3时，冻融修正系数曲线趋于平缓。这也同样印证冻融循环作用使得土体内水分冻结，融化后对土体产生不可逆破坏，土体动剪切模量大幅下降。由于缺少外界水分的补充，在冻融循环3次之后，动剪切模量修正系数的衰减趋于停止。

3. 结论

针对季冻区低温冻融环境下，采用GDS冻土动三轴试验系统，对不同试验温度、不同冻融循环次数下的石灰改良盐渍土的动力参数进行研究，并提出两种修正系数，根据建立的动力参数的预测模型，为工程设计中提供所需的动力参数，同时对长期冻融循环劣化作用后石灰改良盐渍土路基稳定性分析提供参考。

（1）随剪应变的增加石灰改良盐渍土的动剪切模量不断衰减，存在动剪切模量转折点，经过转折点后动剪切模量衰减速度加快；阻尼比随剪应变的增加而增加。动荷载频率及围压增加，石灰改良盐渍土的动剪切模量增加，阻尼比降低。

（2）温度降低，石灰改良盐渍土动剪切模量转折点不断提高，阻尼比下降。提出温度修正系数，对不同温度下动剪切模量和阻尼比进行修正，并对温度修正系数进行曲线拟合，基于此提出温度修正系数公式，该系数在-6℃存在转折点，转折点后低温影响逐渐减弱。

（3）冻融1次后，动剪切模量转折点大幅降低，冻融3次后，动剪切模量转折点进一步降低并趋于稳定，随着冻融循环次数增加阻尼比小幅上升。提出冻融修正系数，对不同冻融循环次数下动剪切模量和阻尼比进行修正，并对冻融修正系数进行曲线拟合，拟合曲线在3次冻融循环处存在转折点，冻融循环的劣化作用存在极限，转折点后影响变小直至忽略不计。提出冻融循环修正系数，基于未冻融动力参数对不同冻融次数下动力参数进行预测。

图 1 DIL806光学热膨胀仪照片

Figure 1. Photo of thermal dilatometer

下载: 全尺寸图片幻灯片

图 2 样品测试原理示意图

Figure 2. Sketch of testing principle

下载: 全尺寸图片幻灯片

图 3 热膨胀系数随样品干密度变化曲线

Figure 3. Relationship between coefficient of thermal expansion and dry density of samples

下载: 全尺寸图片幻灯片

图 4 热膨胀系数随样品含水率变化曲线

Figure 4. Relationship between coefficient of thermal expansion and water content of samples

下载: 全尺寸图片幻灯片

图 5 热膨胀系数随升温速率变化曲线

Figure 5. Relationship between coefficient of thermal expansion and heating rate

下载: 全尺寸图片幻灯片

图 6 残余含水率随升温速率变化曲线

Figure 6. Relationship between residual water content and heating rate

下载: 全尺寸图片幻灯片

图 7 热膨胀系数随气氛环境变化曲线

Figure 7. Relationship between coefficient of thermal expansion and heating atmosphere

下载: 全尺寸图片幻灯片

图 8 残余含水率随干密度变化曲线（不同气氛环境）

Figure 8. Relationship between residual water content and dry density of samples (at different heating atmospheres)

下载: 全尺寸图片幻灯片

图 9 热膨胀率随试样温度变化曲线

Figure 9. Relationship between thermal expansion rate and temperature

下载: 全尺寸图片幻灯片

表 1 高庙子钠基膨润土X射线衍射法矿物成分

Table 1 Mineral compositions of GMZ-Na-bentonite (%)

土样编号	蒙脱石	石英	长石	方英石	斜发沸石	云母
GMZ24	47.7	26.3	16.6	0.9	7.3	1.2

下载: 导出CSV

表 2 试验样品测试方案

Table 2 Test schemes of samples

序号	试样干密度/(g·cm^-3)	含水率/%	升温速率/(℃·min^-1)	气氛环境
1	1.5~1.9	0，5.6，9.34，13.31	5	氩气
2	1.5~1.8	0，5.6，9.34，13.31，15.17	5	氩气
3	1.5~1.7	0，9.34，13.31	1，5，10	氩气
4	1.5~1.7	0，9.34，13.31	5	氩气空气

下载: 导出CSV

参考文献(15)

[1]	王驹, 陈伟明, 苏锐, 等. 高放废物地质处置及其若干关键科学问题[J]. 岩石力学与工程学报, 2006, 25(4): 801–812. doi: 10.3321/j.issn:1000-6915.2006.04.015 WANG Ju, CHEN Wei-ming, SU Rui, et al. Geological disposal of high-level radioactive waste and its key scientific issues[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(4): 801–812. (in Chinese) doi: 10.3321/j.issn:1000-6915.2006.04.015
[2]	刘月妙, 徐国庆, 刘淑芬. 高放废物地质处置库缓冲/回填材料性能测定[J]. 辐射防护, 1998, 18(4): 290–295. https://www.cnki.com.cn/Article/CJFDTOTAL-FSFH804.004.htm LIU Yue-miao, XU Guo-qing, LIU Shu-fen. A study on buffer/backfill materials for hlw geological repository[J]. Radialization Protection, 1998, 18(4): 290–295. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-FSFH804.004.htm
[3]	操龙飞, 徐光, 邓鹏, 等. 钢的热膨胀特性研究[J]. 北京科技大学学报, 2014, 36(5): 639–643. https://www.cnki.com.cn/Article/CJFDTOTAL-BJKD201405011.htm CAO Long-fei, XU Guang, DENG Peng, et al. Study on thermal expansion properties of steels[J]. Journal of University of Science and Technology Beijing, 2014, 36(5): 639–643. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BJKD201405011.htm
[4]	姚武, 郑欣. 配合比参数对混凝土热膨胀系数的影响[J]. 同济大学学报(自然科学版), 2007, 35(1): 77–81, 87. doi: 10.3321/j.issn:0253-374X.2007.01.016 YAO Wu, ZHENG Xin. Effect of mix proportion on coefficient of thermal expansion of concrete[J]. Journal of Tongji University (Natural Science), 2007, 35(1): 77–81, 87. (in Chinese) doi: 10.3321/j.issn:0253-374X.2007.01.016
[5]	陈薇, 杜红秀. 高温对C80高性能混凝土热膨胀性能及其微结构的影响[J]. 中国科技论文, 2017, 12(13): 1477–1481. doi: 10.3969/j.issn.2095-2783.2017.13.006 CHEN Wei, DU Hong-xiu. Effect of high temperature on thermal expansion and microstructure of C80 high performance concrete[J]. China Sciencepaper, 2017, 12(13): 1477–1481. (in Chinese) doi: 10.3969/j.issn.2095-2783.2017.13.006
[6]	刘海涛, 周辉, 胡大伟, 等. 含层理砂岩热膨胀系数的试验研究[J]. 岩土力学, 2017, 38(10): 2841–2846. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201710010.htm LIU Hai-tao, ZHOU Hui, HU Da-wei, et al. Experiment study of thermal expansion coefficient of sandstone with beddings[J]. Rock and Soil Mechanics, 2017, 38(10): 2841–2846. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201710010.htm
[7]	马占国, 唐芙蓉, 戚福周, 等. 高温砂岩热膨胀系数变化规律试验研究[J]. 采矿与安全工程学报, 2017, 34(1): 121–126. https://www.cnki.com.cn/Article/CJFDTOTAL-KSYL201701019.htm MA Zhan-guo, TANG Fu-rong, QI Fu-zhou, et al. Experimental study on thermal expansion coefficient changing rule of sandstone under high temperature[J]. Journal of Mining & Safety Engineering, 2017, 34(1): 121–126. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KSYL201701019.htm
[8]	赵洪宝, 谌伦建. 石灰岩热膨胀特性试验研究[J]. 岩土力学, 2011, 32(6): 1725–1730. doi: 10.3969/j.issn.1000-7598.2011.06.022 ZHAO Hong-bao, CHEN Lun-jian. Experimental study of thermal expansion property of limestone[J]. Rock and Soil Mechanics, 2011, 32(6): 1725–1730. (in Chinese) doi: 10.3969/j.issn.1000-7598.2011.06.022
[9]	AKESSON U. Extensometer Measurement of the Coefficient of Thermal Expansion of Rock[R]. Stockholm: Swedish National Testing and Research Institute, 2004.
[10]	陈皓, 吕海波, 陈正汉. 高庙子膨润土在高温高压下的强度特性研究[J]. 岩土工程学报, 2018, 40(增刊1): 28–33. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2018S1006.htm CHEN Hao, LÜ Hai-bo, CHEN Zheng-han. Strength properties of GMZ bentonite under high temperatures and pressure[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(S1): 28–33. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2018S1006.htm
[11]	刘俊新, 唐伟, 李军润, 等. 高温及碱性条件对高庙子钠基膨润土膨胀力的影响[J]. 岩土力学, 2021, 42(8): 2160–2172, 2184. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202108012.htm LIU Jun-xin, TANG Wei, LI Jun-run, et al. An experimental research on swelling pressure of GMZ Na-bentonite submitted to the strong alkali-heat environment[J]. Rock and Soil Mechanics, 2021, 42(8): 2160–2172, 2184. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202108012.htm
[12]	叶为民, 王琼, 潘虹, 等. 高压实高庙子膨润土的热传导性能[J]. 岩土工程学报, 2010, 32(6): 821–826. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201006003.htm YE Wei-min, WANG Qiong, PAN Hong, et al. Thermal conductivity of compacted GMZ01 bentonite[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(6): 821–826. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201006003.htm
[13]	谢敬礼, 马利科, 高玉峰, 等. 北山花岗岩岩屑-膨润土混合材料导热性能研究[J]. 岩土力学, 2018, 39(8): 2823–2828, 2843. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201808014.htm XIE Jing-li, MA Li-ke, GAO Yu-feng, et al. Thermal conductivity of mixtures of Beishan bentonite and crushed granite[J]. Rock and Soil Mechanics, 2018, 39(8): 2823–2828, 2843. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201808014.htm
[14]	谈云志, 李辉, 王培荣, 等. 膨润土受热作用后的水-力性能研究[J]. 岩土力学, 2019, 40(2): 489–496. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201902010.htm TAN Yun-zhi, LI Hui, WANG Pei-rong, et al. Hydro-mechanical performances of bentonite respond to heat-treated history[J]. Rock and Soil Mechanics, 2019, 40(2): 489–496. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201902010.htm
[15]	唐朝生, 崔玉军, TANG A H, 等. 土体干燥过程中的体积收缩变形特征[J]. 岩土工程学报, 2011, 33(8): 1271–1279. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201108023.htm TANG Chao-sheng, CUI Yu-jun, TANG A M, et al. Volumetric shrinkage characteristics of soil during drying[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(8): 1271–1279. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201108023.htm

施引文献(10)

期刊类型引用(8)

1.	王炜，窦锦钟，袁青云，韩磊，孙旻. 抽水蓄能电站绿色施工技术. 施工技术(中英文). 2024(18): 129-133+153 . 百度学术
2.	张昕阳，申玉生，常铭宇，王浩鱇，潘笑海，王岩岩. 基于GA-BP神经网络的隧道围岩相似材料配合比设计. 现代隧道技术. 2024(06): 82-91 . 百度学术
3.	廖光明. 胶凝砂砾石坝抗滑稳定可靠度分析方法研究. 湖南水利水电. 2023(02): 22-25 . 百度学术
4.	丁泽霖，王肖楠，张宏洋，王婧，包闯. 复杂地基重力坝与胶凝砂砾石坝变形破坏对比分析. 工程科学与技术. 2023(04): 197-206 . 百度学术
5.	杨锐婷，万露. 折线形胶凝砂砾石坝三维静力有限元分析. 吉林水利. 2023(07): 14-20 . 百度学术
6.	丁泽霖，徐良杰，王肖楠. 基于多元优化的百米级胶凝砂砾石坝设计方法. 水力发电学报. 2023(12): 119-131 . 百度学术
7.	杨玉坚，刘艳华，陈雪冬. 胶凝砂砾石坝导流设计方法研究. 水利建设与管理. 2023(S1): 35-39 . 百度学术
8.	刘欢. 基于AutoBANK的流固耦合场沙田水库坝体渗流变形特征研究. 水利科学与寒区工程. 2022(08): 26-30 . 百度学术