黏土颗粒间孔隙等效孔径的方法
【技术领域】
[0001] 本发明属于岩土工程实际应用技术领域，涉及一种黏土颗粒间孔隙等效孔径的方 法。
【背景技术】
[0002] 非饱和土由土颗粒、水分、空气和收缩膜组成。而土内孔隙的大小、形状及分布是 土的结构组成的重要因素，对与土体体积相关的总体工程行为起着重要的作用。测定或估 算孔隙物理性能对预测土体强度、压缩性、渗透性也有着深远的影响。土体吸力的大小能够 表征土颗粒的紧密程度，含水状况以及非饱和土中收缩膜的曲率的大小，因此，对于使用 土-水特征曲线估算孔径分布已有较多的研究。但是，已有的研究多局限于评价非饱和土内 孔径、孔径分布与毛细压力之间的关系，并未给出计算孔隙大小的定量方法。
【发明内容】
[0003] 本发明的目的是提供一种黏土颗粒间孔隙等效孔径的方法，为最大毛细上升高度 的计算提供了一种有利的途径，解决了现有技术中存在的问题。
[0004] 本发明所采用的技术方案是，一种黏土颗粒间孔隙等效孔径的方法，具体按照以 下步骤进行：
[0005] 步骤 1，
[0006] 首先制备含黏粒土样进行土 -水特征实验，并通过MATLAB用（1)式拟合得到完整的 土-水特征曲线，以质量含水量表示，拟合公式如下：
[0007]
〉
[0008] 式中，ua_uw为基质吸力；
[0009] 残余含水率；
[0010] 0S饱和含水率；
[0011] a、m及η均为拟合参数，参数η与土体孔径分布有关，参数m与土-水特征曲线的整体 对称性有关；
[0012]步骤 2，
[0013]以（Ua-Uw)i = (Ua-Uw)i-1+50(i_l)(i = 2,3,4......)，（ua_uw)i = lOOkPa的方式选取基 质吸力于144000kPa~OkPa之间的数据，通过(4)式将基质吸力转化为相对湿度RH: _4]
⑷
[0015]由以上推导可得，毛细半径r以相对湿度表示形式如下：
[0016] (5a)
[0017] 或表示为基质吸力形式： ，：
[0_(5b) -Ι?
[0019] 式中：水蒸气的偏摩尔体积(m3/m〇l);
[0020] R为通用气体常数(J/(mol · K));
[0021] T为热力学温度(K);
[0022] Ts为毛细管内液相表面张力；
[0023] α为土颗粒与孔隙水之间的接触角；
[0024] 步骤 3，
[0025] 当RH或ua-uw为第i步增量时，单位质量土体被气体或水填充的孔隙体积$用式(6) 表不：
[0026] (6)
[0027] 式中，W为含水量。
[0028]假设固液相互作用只对土颗粒表层附近水膜密度产生显著影响，则毛细吸附区PW 本质上是常数，Pw=lg/cm3，
[0029] 通过（7)式计算单位质量土体内水所占的孔隙体积减小量及累积孔隙体积变化 量，则每单位质量固体内水的排出量：
[0030]
.(.7.)
[0031 ]式中，i = 2,3,4……
[0032] 步骤 4，
[0033] 假设接触角α = 0,通过式(5a)或式(5b)中的任何一式计算开尔文半径；
[0034] 步骤 5，
[0035] 计算实际孔隙半径< 和吸附水膜厚度表达如下：
[0036],…
， (8)
[0037] 其中吸附水膜厚度t1用Halsey方程表示如下：
[0038]
(分）
[0039] 式中，t1为相对湿度在第i增量步时所对应的吸附水膜厚度；
[0 04 0 ] τ为吸附水分子的有效直径D
[0041] 假设液体水分子所占区域的横截面面积约为Α=10,8Α2,水的摩尔体积1.8 X 10-5m3/mol，Avogadro常数Na= 6 · 02 X 1023/mol，则吸附水分子的有效直径τ:
[0042] ^ T A
，
[0043] 步骤 6，[0044] 用（10)式计算平均孔隙半径WU，
[0045] ,、 α〇)
[0046] 式中，i = 2,3,4……
[0047] 步骤 7，
[0048] 选择孔隙体积Ag大于0.01cm3所对应的孔径，根据下式计算黏土颗粒间孔隙等效 孔径：
[0049] (11)
[0050] 式中，d为黏土颗粒间孔隙等效孔径；
[0051] (4 为孔隙体积<在0. 〇lcm3以上时所对应的孔径;β为无量纲修正系数，与土 体内封闭气泡及体积小于0.01cm3孔隙等因素有关，根据大量试验数据统计，取值0.02;
[0052] 步骤 8，
[0053] 将计算得到的等效毛细直径代入公式(12)计算最大毛细上升高度：
_ (12)
[0055] ______________s2)。
[0056] 本发明的有益效果是以完整的土-水特征曲线为载体，建立在已有研究的基础之 上，通过大量实验并结合统计理论，推导得到不同黏粒含量土样的孔径的等效方法，为最大 毛细上升高度的计算提供了一种有利的途径。
【具体实施方式】
[0057] -种黏土颗粒间孔隙等效孔径的方法，具体按照以下步骤进行：
[0058] 步骤 1，
[0059] 首先制备含黏粒土样进行土 -水特征试验，并通过MATLAB用（1)式拟合得到完整的 土-水特征曲线，以质量含水量表示。拟合公式如下：
[0060]
⑴
[0061] 式中:ua_uw为基质吸力；
[0062] 残余含水率；
[0063] 03为饱和含水率；
[0064] a、m及η均为拟合参数。参数η与土体孔径分布有关，参数m与土-水特征曲线的整体 对称性有关。
[0065] 土体孔隙结构决定了处于平衡状态的蒸汽压大小。Kelvin(开尔文）公式，则在不 违背土体体系原有特性的前提下，把土体孔隙理想化为一个半径为r、土颗粒与孔隙水之间 的接触角为α、液相表面张力为Ts的毛细管，管内水-土交界曲面上的总压力变化Ua_Uw为：
[_
. (2)
[0067] 为描述气-水交界曲面上的压力变化与交界曲面上方蒸汽压之间的联系，开尔文 公式可表示为：
[0068]
(3)
[0069] 式中，R为通用气体常数(J/(mol ·Κ));
[0070] Τ为热力学温度(K);
[0071 ] u?为温度T条件下平衡状态的自由水饱和蒸汽压值(kPa);
[0072] uvl为平衡状态的溶液蒸气压值(kPa);
[0073] 为水蒸气的偏摩尔体积(m3/mol)。
[0074] 定义uvl/u?为相对湿度RH。
[0075] 步骤 2，
[0076] 根据⑶式可以建立相对湿度与基质吸力之间的函数关系，以(Ua-udF^-uJi-1 +50(i_l) (i = 2,3,4......)，（ua_Uw)i = lOOkPa的方式选取基质吸力于144000kPa~OkPa之间 的数据，通过(4)式将基质吸力转化为相对湿度RH。
、 '9[0078] 由以上