不同含液量下颗粒间液桥力及形态的试验研究

2020年1月文章编号：0559-9350（2020）01-0081-11

利

SHUILI学

XUEBAO报

第51卷第1期不同含液量下颗粒间液桥力及形态的试验研究

蒲诚，刘奉银，张昭，程靖轩，赵伟（西安理工大学岩土工程研究所，陕西西安

710048）

摘要：对液桥力的研究有助于从根本上揭示非饱和土持水特性的内在机理。以等径球体颗粒间的液桥为研究对象，利用纳米多功能试验机及CCD数码相机分别测量并记录了三种粒径、六种含液量下液桥拉伸过程中力值及形态的变化，并将试验结果与理论计算结果进行对比，同时对不同类型液桥拉伸曲线产生的原因提出了猜想。结果表明：液桥力-位移曲线可以分为上升段、平稳下降段以及突然跌落段三部分；不同含液量下液桥的初始形态不同，钟摆状液桥主要出现在饱和度较低且重力影响可以忽略的条件下，随着饱和度的增加液桥初始形态变为圆柱状以及外凸状；最大液桥力与粒径及含液量成正比，断裂距离与含液量成正比与粒径成反比，且粒径对最大液桥力以及断裂距离的影响更为显著；不同研究者观测到的两类液桥力-位移曲线可能是试验初始点不同导致的。关键词：非饱和土；持水特性；液桥拉伸；试验研究中图分类号：TU43

文献标识码：A

doi：10.13243/j.cnki.slxb.20190736

1研究背景

非饱和土的持水特性是其区别于其他类型土最主要的性质，正确认识非饱和土的持水特性

是岩土工程性能评价的基础和关键，也是科学研究的热点问题。长期以来众多学者采用基质吸力、含水率、饱和度等参数来表征非饱和土的持水特性，取得了长足的进展。但是随着认识水平的提升，人们发现，不同初始条件下土样的持水特性反映的试验规律适用范围较窄、物理机制不清，加之不能观测到土样内部的变化情况，想要正确认识和从根本上揭示非饱和土持水特性的内在机制是远远不够的。作为一种典型的散体材料，当位于潮湿环境中时土体会从周围环境中吸收水分，在颗粒接触点处形成液桥，有学者认为非饱和土持水特性的变化过程从本质上可以归结为土颗粒间液桥力的演化，因而，从液桥的角度认识和解释持水特性日益为土力学界所关注。

对液桥的研究最早可以追溯到1960年代表面科学领域。随后，为了计算简单起见，Gillespie等

［2］

［3］

［4］

［5］

［1］

和Clark等将液桥的外表面假定为圆形；DeBisschop等、Pietsch分别研究了两光滑球体以及粗糙球体颗粒间液桥的液桥力。Willett等利用微分天平分别测量了等径与不等径人造蓝宝石颗粒间毛细力的变化规律；Damiano等利用小挠度的悬臂装置，通过测量悬臂的变形来间接计算一对等径颗粒间液桥力的大小；Diana等、WangJi-Peng等在两颗粒拉伸试验的基础上进一步测量了三颗粒拉伸时液桥毛细力的变化。刘建林等依据最小势能原理推导了轴对称液桥的形貌和液桥力；王学卫等从数值模拟的角度研究了考虑重力影响下平板间液桥的断裂距离；庄大伟等和土持水特性做出了积极的探索，Likos等

［13］

［11］

［9］

［10］

［7］

［8］

［6］

、朱朝飞等

［12］

分别从

［14］

试验和理论分析的角度研究了狭长平板间液桥的形态。基于液桥理论，岩土工程学界的学者对非饱

研究了固-液接触角对非饱和土强度的影响；Yang等

收稿日期：2019-10-21；网络首发时间：2020-01-20

网络出版地址：http：//kns.cnki.net/kcms/detail/11.1882.TV.20200120.0759.001.html基金项目：国家自然科学基金项目（51679198，41372304，41877278，41402258）

作者简介：蒲诚（1995-），博士生，主要从事非饱和土水力特性研究。E-mail：1909488416@qq.com

通讯作者：刘奉银（1968-），教授，主要从事黄土力学与非饱和土力学研究。E-mail：fyliu@pub.xaonline.com

—81—

建立了基于固液接触角预测边界减湿土-水特征曲线的理论模型；LuNing等颗粒间的液桥力，并从微观角度解释土-水特征曲线的变化规律。贺炜等对土-水特征曲线的影响；栾茂田等玺等

［18-19］

［17］

［16］

［15］

试验测量了两不等径

研究了颗粒间固液接触角

［20-23］

［24］

从粒间液桥角度出发推导了非饱和土土水特征曲线公式；周凤

从理论角度推用微分天平手

理论推导了颗粒间毛细作用力并探讨其与土体液塑限的内在联系；张昭等

导了一对不等径颗粒间毛细力的计算方法，并将其延伸至黏土颗粒范围；余连英等动测量了球形颗粒间液桥拉伸过程中液桥力的变化。

然而，国内对于液桥的研究主要集中在理论推导和数值模拟层面，本文利用一种全新的纳米拉伸装置，配合高清度CCD数码相机，从试验角度测量和揭示三种粒径、六种含液量下一对等径球体颗粒间液桥力大小与形态的变化规律。

2

2.1

试验步骤

试验仪器

采用美国Keysight公司研发的UTMT150纳米多功能试验机进行试验，在试验机的前

侧以及左侧分别放置高精度显微照相机，确保上下颗粒处于同一轴线上并实时记录液桥形状的变化。试验机由刚性外壳、减振台、测量装置、数据采集系统四部分组成，刚性外壳可以最大程度减小液桥断裂时由于应力突然释放导致的波动；减振台可以消除周围操作人员的声音和脚步振动对试验结果产生的影响；测量装置主要由下端的机械驱动转换器以及上端可移动的刚性机械臂组成；数据采集系统分别连接在电脑和仪器的中央控制器上，用以操作仪器并实现数据的实时采集和分析处理，仪器构造示意图如图1所示，主要力学参数如表1所示。

图1试验仪器示意图

表1

项目参数

最大荷载/mN

500

试验机的主要参数

最大伸长位移/mm

200

位移分辨率/nm

35

拉伸速率0.5μm/s~5mm/s

荷载分辨率/nN

50

2.2面

［25］

试验材料选取3组直径分别为2.5、4和5mm的等径玻璃球；采用有机溶剂丙三醇模拟颗粒

间的液桥，20℃时，丙三醇的性质如表2所示。采用无水乙醇和去离子水作为清洗剂清洗玻璃珠表

，以清除颗粒表面灰尘和消除颗粒表面静电力的影响。

表2

液体丙三醇

丙三醇与水的物理力学参数（20℃）

表面张力（/N/m）

0.063

黏滞系数（/Pa/s）

0.15

密度（/g/cm3）

1.26

—82—

表3

饱和度Sr0%Sr>100%

液桥形态钟摆状索带状毛细管状泥浆状

液桥形态分类

示意图

作用机理

颗粒接触点处形成液桥部分颗粒间隙充满液体颗粒间几乎全部充满液体液体压力大于空气压力

［26］

2.3等

［27］

液桥体积的确定选定0.05、0.1、0.2、0.5、1.0、1.5μL六种液桥体积。Fournier等

、Rossetti

根据饱和度的大小及液桥与颗粒间的联系方式，将液桥的形态分为如表3中所示的四类。

将颗粒间距D=0时颗粒间的孔隙视为孔隙体积VV，如图2中浅灰色所示；将液桥体积等同于孔隙中水的体积VW，如图2中深灰色所示；则有：

Sr=

VWV=2VV4R-πR2（1）

式中：V为液桥体积，μL；R为颗粒半径，mm。

当液桥为索带状、毛细管状以及泥浆状时，由于含液量较大饱和度较高，重力对液桥力及形态的影响一般不可忽略。而对于钟摆状液桥而言，Adams等人通过研究重力与颗粒间液桥体积的映射关系，采用无量纲的液桥体积V*与BO数的乘积来定性反映重力对钟摆状液桥力及形态的影响，认BO<0.015时，处于过渡阶段

［25］

为：当V*·BO<0.01时，重力的影响可以忽略；当V*·BO>0.015时，重力的影响不可忽略；当0.01。

BO=Dρgd

2

σ（2）（3）

V*=VR3

式中：d为特征长度，m，是液桥体积的函数；g为重力加速度，m/s2；Dρ为液体和外部气体的密度差值，kg/m3；σ为液体的表面张力，N/m。

BO数的概念最早用以表征重力对自由液体下落时液滴形状的影响，表示为BO=DρgD2σ，其中D为液滴的直径；随后，BO数被引入到液桥领域用来表征重力对球-板间液桥力与形态的影响，表示为BO=DρgR2σ，其中R为球的半径；而对于球-球间钟摆状液桥，Adams将其定义为式（2）的形式，其中d为特征长度，是一个无任何物理含义的值，Adams将其近似表示为d=VD，其中V为液桥体积，D为液桥的直径。本文通过CCD数码相机自带的图像处理软件ImageView对不同条件下钟摆状液桥的直径进行量测，计算发现对于本试验而言特征长度值在0.3mm至1.5mm之间变动，为了计算方便统一取d=0.001m。

对于本试验而言，Dρ=1260kg/m3，g=9.81m/s2，σ=0.067N/m，颗粒半径R分别为2.5、4和5mm。将其分别带入式（1）—（3）计算出不同粒径、不同液桥体积所对应的Bo数、无量纲的液桥体积V*及饱和度Sr如表4所示。由表4可知：试验条件分别对应重力影响可以忽略、重力影响可忽略过渡状态以及重力影响不可忽略的情形，液桥形态分别为钟摆状、索带状、毛细管状及泥浆状。2.4

试验步骤

将清洗后的玻璃珠用硬基质胶固定在测量装置的刚性机械臂及NMAT延长头上，采

用微型移液枪在下部颗粒球冠处注入相应体积的液体，控制上部颗粒向下移动，待形成液桥之后反复拉伸多次使其形态趋于轴对称。试验时采用位移控制，确保上下两个颗粒接触距离为0，关闭试验机两侧及前侧的玻璃门使试验主体处于封闭状态，卸下减振台上的紧固螺栓，保持下部颗粒稳定，使上部颗粒以100μm/s的速度向上移动，设置荷载记录触发值为100nN，同时利用CCD工业电子显微镜记录试验过程中液桥形态的变化。

—83—

表4

颗粒直径/mm

液桥体积/μL

0.050.10

2.5

0.200.51.00.050.10

4

0.200.51.00.050.10

5

0.200.51.01.51.51.5

不同液桥体积对应的Bo数及饱和度Sr

液桥形态钟摆状钟摆状钟摆状索带状毛细管状泥浆状钟摆状钟摆状钟摆状钟摆状钟摆状索带状钟摆状钟摆状钟摆状钟摆状钟摆状钟摆状

3.2×10-36.4×10-30.0130.0320.0640.096

5.91×10-41.18×10-32.36×10-35.91×10-30.01870.012

6.25×10-30.01250.06250.1250.025

1.16×10-32.31×10-34.62×10-3

0.0120.023

V*0.0250.0510.102

V*·BO

考虑重力

否否是是是是否否否近似过渡

是是否否否否近似过渡

是

饱和度/%3.717.4314.937.1111.41.452.915.8214.629.243.70.931.863.739.3218.627.974.3

4.7×10-39.4×10-30.018

3

3.1

试验结果分析

液桥力-位移曲线分析

张昭等

ψ

［20］

将粗粒样简化分散为球体颗粒，将水分形态视为液桥，认为液

桥力Fliq是液桥中由基质吸力F产生的毛细斥力以及由表面张力Fσ产生的毛细引力共同作用的结果，如式（4）；根据Y-L方程可将气液交界面处的基质吸力ψ表示为式（5），液桥颈部受力如图3所示。

孔隙液桥

毛细斥力毛细引力

图2饱和度计算示意图图3液桥颈部受力示意图

r1

R

θ

r2

图4液桥几何参数示意图

—84—

Fliq=F

σ

+F

ψ

æöψ=σç1-1÷r2Røèr1RF

σ

=2π(r2R)σ-π(r2R

)2

ψ（4）（5）

式中几何参数如图4所示，其中：r2为液桥颈部半径、r1为液桥外轮廓半径、R为颗粒半径。当

>F时，液桥力为正，液桥形成；当F

ψ

σ

ψ

图5为液桥力-位移曲线图，由图5可知，曲线分为上升段、下降段以及突然跌落段三部分。在曲线的上升段，液桥力随着液桥拉伸距离的增大不断地增大，关于这一部分不同的研究者得出了不同的结论，Rossetti等等

［29］

［27］

、Bozkurt等

［28］

认为最大液桥力出现在分离距离较小但不为零处，然而Olivier

认为最大液桥力出现在分离距离为零处，虽然不同研究者对此持有不同的观点，但并未对此现

象进行系统的分析，本文将在第四部分对其产生的原因提出猜想。在曲线的下降段，毛细斥力随着液桥拉伸距离的增加而增加而毛细引力随液桥拉伸距离的增加逐渐减小，但液桥内部仍以引力为主，当液桥的颈部变为最窄时毛细引力减到最小毛细斥力最大进入突然跌落段；当位于跌落段时，液桥不能再承受任何形式的拉伸，液桥突然断裂颗粒分离。曲线最后的残余部分为残留在下球表面液滴的重量。

0.250.20桥力/mN0.150.100.050.00

0.0

0.5

1.0

1.5

2.0

2.50.05μL

0.10μL0.20μL0.50μL1.00μL1.50μL

0.40.3桥力/mN0.20.10.0

0.05μL0.10μL0.20μL0.50μL1.00μL1.50μL

拉伸位移/mm（a）2.5mm

0.60.5桥力/mN0.40.30.20.1

0.00.5

拉伸位移/mm（b）4mm

1.01.52.0

0.05μL0.10μL0.20μL0.50μL1.00μL1.50μL

0.00.0

0.5

1.0

拉伸位移/mm（c）5mm

1.52.0

图5液桥力-位移关系

定义曲线到达峰值点的位移为峰值位移。含液量对液桥力-曲线的影响主要以饱和度及Bo数的形式反映在对上升段峰值位移及下降段曲线形态的影响上。在曲线的上升段，饱和度对曲线形态的影响较小，曲线峰值位移随饱和度的增加不断增加。在曲线下降段，当饱和度较小且重力影响可以忽略时，曲线下降趋势为“凹”状，如2.5mm颗粒0.05μL含液量以及5mm颗粒0.05～1.0μL含液量；当饱和度较大且重力影响不可忽略时，曲线下降趋势多成直线状或“凸”状，如2.5mm颗粒

—85—

图6液桥形态变化图

0.2～1.5μL含液量以及4mm颗粒1.5μL含液量。相对于含液量，无量纲的液桥体积更能反映液桥体体积变化幅度更大，因而液桥力-位移曲线的下降段离散性及差异性也更明显。3.2

液桥形态的变化

积变化对曲线形态的影响，当液桥体积从0.05μL增加至1.5μL的过程中，2.5mm颗粒的无量纲液桥

利用CCD相机全程录像，记录不同粒径及含液量下液桥形态的变化过程。由

于篇幅的原因，本文仅以2.5mm粒径为对象分析不同含液量下液桥形态的变化规律。图6（a）至图6（f）分别反映了颗粒粒径为2.5mm，液桥体积为0.05μL、0.1μL、0.2μL、0.5μL、1μL、1.5μL等6种情况下液桥从形成到断裂的形态变化。用X分别代表a至f，图6（X-1）表示液桥的初始形态，图6（X-6）表示液桥最终断裂的形态，图6（X-2）至图6（X-5）分别反映了拉伸过程中不同时刻的液桥形态。

结合图6及表4可以发现液桥的形态是由重力以及饱和度共同影响的结果：当液桥体积为0.05μL以及0.1μL时处于为钟摆状且重力影响可以忽略的情况，液桥的初始形态符合圆环假设，在颗粒接触点处形成如图6（a-1）、图6（b-1）所示的凹液桥，在拉伸过程中液桥始终保持钟摆状，当液桥最终断裂时液桥等体积的分布在上下颗粒的球冠处；当液桥体积为0.20μL以及0.5μL时分别处于毛细状液桥重力不可忽略以及索带状液桥的情况，液桥初始形态可近似为图6（c-1）、图6（d-1）所示的圆柱状，随着拉伸距离的增加液桥由圆柱状迅速变为图6（c-4）、图6（d-3）所示的符合圆弧假定的钟摆状，随着拉伸距离的增加重力对液桥形状的影响越来越明显当处于临近断裂时液桥形状为图6（c-5）、图6（d-5）所示的上部外曲率较小下半部外曲率较大的钟摆状，最终断裂时下球残留液体

—86—

0.70.6最大液桥力/mN0.50.40.30.20.10.00.0

D=2.5mmD=4.0mmD=5.0mm

断裂距离/mm1.61.41.21.00.80.60.0

D=2.5mmD=4.0mmD=5.0mm

0.2图7

0.40.6

液桥体积/μL

0.81.01.21.41.60.20.4图8

0.60.8

液桥体积/μL

1.01.21.41.6

最大液桥力-液桥体积关系断裂距离-液桥体积关系

的体积稍大于上球残留液体的体积；当液桥体积为1μL以及1.5μL时分别处于毛细管状以及泥浆状的情况，此时重力的影响非常显著，液桥初始形态为外凸形并且发现了如图6（e-1）所示的非轴对称液桥形态，随着拉伸试验的进行液桥首先转化为圆柱状，后迅速变为图6（e-5）、图6（f-5）所示的上部外曲率较小下半部外曲率较大的轴对称钟摆状，最终断裂时下球残留的液体体积明显大于上球残留液体体积。3.3

最大液桥力及断裂距离分析

将三种粒径，六种含液量下的最大液桥力和断裂距离分别绘制在

图7及图8中。由图7可以看出，颗粒间的液桥力经历了由快速增长到缓慢增长的变化过程：当液桥体积从0.05μL增加到0.2μL的过程中，随液桥体积增加液桥力增加速率较快，即使很小的液桥体积的改变也会导致液桥力的迅速增加；当液桥体积从0.2μL增加到1.5μL过程中液桥力随液体体积增加变化幅度较小。可以推论液桥力的增加必然存在一个阈值，即液桥力不会随着液体含量的增加而无限制的增加下去。相对于液桥体积，粒径对液桥力的影响更为显著，最大液桥力与粒径大小呈正比，含液量相同时粒径越大最大液桥力越大。

定义液桥突然断裂时对应的的拉伸距离为断裂距离。将断裂距离绘制在图8中，由图可以看出，随着液桥体积的增加断裂距离也经历了由快速增长到缓慢增长的过程，且断裂距离随液桥体积的变化趋势与液桥力相同：当液桥体积从0.05μL增加到0.2μL的过程中，断裂距离随液桥体积增加的速率较快；而当液桥体积从0.2μL增加到1.5μL的过程中，断裂距离随液桥体积增加的速率放缓。粒径对断裂距离的影响更为显著，断裂距离与粒径大小呈反比，含液量相同时粒径越大断裂距离越小。3.4

相关理论

基于圆环理论，假定拉伸过程中液桥外轮廓半径以及固液接触角保持恒定，通过

［16，19-20］

几何关系迭代并结合Y-L方程，国内研究者给出了钟摆状液桥液桥力的计算公式。Olivier提

出了一种认可度较高且较为简便的液桥力计算方法，如前文图5所示，认为当R>>r1>>r2且液桥高度<<2r1cosθ时，液桥力可以表示为如式（6）所示的毛细力Fcap以及黏滞力Fvis的和：

0.5最大液桥力/mN0.40.30.20.10.0

0.0

0.2图9

0.4

0.6

0.8

1.0

1.2

1.4

D=2.5mmD=4.0mmD=5.0mm

断裂距离/mm2.22.01.81.61.41.21.00.80.6

液桥体积/μL

最大液桥力-液桥体积计算对比

0.0

0.2

0.4图10

0.6

D=2.5mm实测D=4.0mm实测D=5.0mm实测D=2.5mm计算D=4.0mm计算D=5.0mm计算

0.81.01.2液桥体积/μL

1.4

1.6

断裂距离-液桥体积计算对比

—87—

ìF=F+F=2πRσcosθX+3πηR21dDX2

capvisvïliq2Ddtv

ï

-1í

2Vù2ïX=1-é1+êúïv2πRDëûî

（6）

式中：σ为液体表面张力，N/m；θ为固-液接触角；η为液体的黏滞系数，Pa·s；D为分离距离，mm；dD/dt为分离速度。

Lian

［31］

*

给出了无量纲的断裂距离Druat与无量纲液桥体积V*及固液接触角θ之间的关系如下：

*

Drupt=æ1+θöV

2øè

*13

（7）

*

式中：Drupt=DruptR；θ取弧度值。

由于式（6）的限定条件，本文仅对含液量0.05μL不同粒径颗粒的液桥力-位移曲线进行模拟（如图9），将无量纲的断裂距离转化为实际断裂距离并与本文试验结果对比绘制在图10中。

可以看到，式（6）对于最大液桥力的预测较准且无量纲液桥体积越小预测的精度越高，但其无法预测液桥力-位移曲线的上升段以及突然跌落段，对于液桥力-位移曲线的缓慢下降段预测的精度不高，这也是目前所有基于圆环理论推导得出的液桥力计算公式的普遍困局。式（7）能较好地预测断裂距离随含液量增加而增加，随粒径增大而减小的趋势，但不能反映在重力影响较大情况下断裂距离增加速率随含液量增加逐渐放缓的趋势。

4液桥力-位移曲线上升段机理猜想

对于液桥力-位移曲线的上升段，不同的研究者持有不同的观点，部分研究者认为最大液桥力出

现在分离距离为零的位置，另一部分研究者认为最大液桥力出现在分离距离较小但不为零的位置。本文试图从试验方法的角度对曲线上升段产生的原因提出猜想。

以往的研究者大多采用微分天平法以或悬臂梁法测量液桥力。以悬臂梁法为例，试验时先将左侧球体颗粒固定在小刚度悬臂梁顶端，在右侧颗粒球冠处注入相应体积的液体将其固定在可移动的刚性悬臂端，保持左侧小刚度悬臂梁挠度为零，移动右侧刚性悬臂使颗粒间形成液桥，通过拉伸过程中小刚度悬臂梁挠度的变化计算液桥力。当采用悬臂梁法时，试验操作人员对初始点的选择对液桥力-位移曲线有很大的影响：如图11（a），当采用颗粒紧密接触临界位置为试验起始点时，液桥力-位移曲线则会观察到上升段；如图11（b），当采用小刚度梁不产生压变形的临界位置为试验起始点时，则液桥力-位移曲线没有上升段。随着试验技术的不断改进，高精度的刚性拉伸试验机开始广泛运用在液桥拉伸试验中。如图11（c），将试验机的核心部分简化为一个高灵敏度的弹簧和可移动的刚性悬臂，通过拉伸过程中弹簧的变形值来计算液桥力的变化，为了避免操作失误引起的刚性悬臂在下降过程中对下部弹簧造成损坏，在试验机设计时往往会预设部分弹簧伸长量以保护试验装置，而试验者通常采用位移控制上部颗粒向下移动以颗粒紧密接触临界位置为试验起始点，在试验开始时下部弹簧储存有一部分压变形。

综上所述，当试验以不发生压变形为初始点时，颗粒与液桥处于力平衡状态，但此时颗粒间的接触距离可能并不为零；当以颗粒紧密接触为初始点时，颗粒间的接触距离为零但此时液桥中液体对另一颗粒的顶端作用有压应力，随着液桥的不断拉伸，作用在颗粒顶端的压应力逐渐减小曲线呈现上升状态，即两种液桥力-位移曲线产生差异的根源原因可能在于其试验初始点的选择不同。

—88—

刚性悬臂小刚性悬臂（a）悬臂梁法以位移零点为起始点（b）悬臂梁法以力平衡点为起始点

刚性悬臂弹簧（c）刚性试验机法图11

不同试验方法示意图

5结论与展望

使用纳米多功能试验机进行等径球体颗粒间的液桥拉伸试验，同时借助CCD相机观测液桥形态

的变化，从试验角度揭示了液桥力的变化规律，将试验结果与理论结果进行对比并初步探讨了液桥力-位移曲线上升段的产生原因，得到以下结论及展望。

液桥力-位移曲线可以分为上升段、平稳下降段以及突然跌落段三部分。当饱和度较小时且重力影响可以忽略时曲线下降段为“凹”状，当饱和度较大且重力影响不可忽略时曲线下降多成直线状或“凸”状，且相对于含液量，无量纲的液桥体积的差异更能反映液桥体积变化对曲线形态的影响。

液桥的初始形态受饱和度及重力的共同影响。钟摆状液桥主要出现在饱和度较低且重力影响可以忽略的情况下；当饱和度增加且重力影响不可忽略时，液桥的初始形态变为圆柱状及外凸状，并随着拉伸试验的进行转变为上部外曲率较小下部外曲率较大的钟摆状液桥。

随着液桥体积的增加，最大毛细力和断裂距离均经历了从快到慢增加的过程，低含液量下液桥力及断裂距离对液体体积的变化较为敏感。相对于液桥体积，粒径对液桥力及断裂距离的影响更为显著，最大液桥力与粒径大小呈正比，含液量相同时，粒径越大最大液桥力越大；断裂距离与粒径大小呈反比，含液量相同时，粒径越大断裂距离越小。

不同研究者观测到的两种类型液桥力-位移曲线，本文认为可能是由于试验初始点的不同导致的，关于这一猜想有待后续研究加以验证。

参考文献：

［1］GILLESPIET，SETTINERIWJ.Theeffectofcapillaryliquidontheforceofadhesionbetweensphericalsolid

particles［J］.JournalofColloidandInterfaceScience，1967，24（2）：199-202.

—89—

［2］CLARKWC，HAYNESJM，MASONG.Liquidbridgesbetweenasphereandaplane［J］.ChemicalEngineer⁃［3］DeBISSCHOPFRE，RIGOLEWJL.Aphysicalmodelforliquidcapillarybridgesbetweenadsorptivesolid［4］PIETSCHWB.Tensilestrengthofgranularmaterials［J］.Nature，1968，217（5130）：736-737.

muir，2000，16（24）：9396-9405.

spheres：thenodoidofplateau［J］.JournalofColloidandInterfaceScience，1982，88（1）：117-128.ingScience，1968，23：810-812.

［5］WILLETTCD，ADAMSMJ，JOHNSONSA，etal.Capillarybridgesbetweentwosphericalbodies［J］.Lang⁃［6］DAMIANOR，XAVIERP，STEFAANJ.R.Simons.Ruptureenergyandwettingbehaviorofpendularliquid

（1）：161-169.

bridgesinrelationtothesphericalagglomerationprocess［J］.JournalofColloidandInterfaceScience，2003，261

［7］DIANALievano，SACHINVelankar，JOSEPHJMcCarthy.Theruptureforceofliquidbridgesintwoandthree［8］WANGJP，ELENAG，BERTRANDF.Capillaryforceandruptureoffunicularliquidbridgesbetweenthree［9］刘建林，李广帅，聂志欣.轴对称液桥的形貌与液桥力［J］.西华大学学报（自然科学版），2010，29（3）：［10］王学卫，于洋.重力影响下板间液桥断裂距离研究［J］.实验力学，2012，27（1）：70-76.

［11］庄大伟，杨艺菲，胡海涛，等.竖直平板间液桥形状的观测与预测模型开发［J］.化工学报，2016，67（6）：［12］朱朝飞，贾建援，付红志，等.狭长平行板间液桥形态及受力研究［J］.工程力学，2016，33（6）：222-229.［13］LIKOSWJ，LUN.Hysteresisofcapillarystressinunsaturatedgranularsoil［J］.JournalofEngineeringMechan⁃［14］YANGS，LUTH.Studyofsoil-watercharacteristiccurveusingmicroscopicsphericalparticlemodel［J］.Pedo⁃［15］LUNING，JEREMYLechma，KELLYTMiller.Experimentalverificationofcapillaryforceandwaterretention［16］贺炜，赵明华，陈永贵，等.土-水特征曲线滞后现象的微观机制与计算分析［J］.岩土力学，2010，31（4）：［17］栾茂田，李顺群，杨庆.非饱和土的基质吸力和张力吸力［J］.岩土工程学报，2006，28（7）：863-868.［18］周凤玺，张海威，张家齐.基于粒间毛细作用探讨界限含水量［J］.兰州理工大学学报，2018，44（3）：［19］周凤玺，曹小林，马强.颗粒间的毛细作用以及吸应力特征曲线分析［J］.岩土力学，2017，38（7）：［20］张昭，刘奉银，齐吉琳，等.粗颗粒间液桥毛细力演化规律的动态计算方法［J］.岩土力学，2016，37（8）：［21］张昭，刘奉银，齐吉琳，等.考虑固-液接触角影响的粗颗粒间液桥毛细力计算方法［J］.水利学报，2016，［22］张昭，刘奉银，李荣建，等.黏土颗粒间液桥从形成至断裂时毛细力的演化规律［J］.岩土力学，2016，37

（10）：2845-2850，2860.

［23］张昭，刘奉银，张国平，等.不等径湿颗粒与液桥相互作用的微观水力特性［J］.水利学报，2013，44（7）：［24］余莲英，周丹洋，徐春晖，等.基于拟土颗粒DEM数值模型的静态液桥力测试方法［J］.中国农业大学学

报，2017，22（11）：68-74.

［25］ADAMSMJ，JOHNSONSA，SEVILLEJPK，etal.Mappingtheinfluenceofgravityonpendularliquidbridges［26］FOURNIERZ，GEROMICHALOSD，HEMINGHAUSS，etal.Mechanicalpropertiesofwetgranularmaterial

［J］.JournalofPhysics，2005，17：477-491.

［27］ROSSETTID，SIMONSSJR.Amicroscaleinvestigationofliquidbridgesinthesphericalagglomerationprocess

betweenrigidspheres［J］.Langmuir，2002，18（16）：6180-6184.810-817.

47（9）：1197-1207.2263-2270，2278.2036-2042.115-118.1078-1083.

betweenuneven-sizedspheres［J］.JournalofEngineeringMechanics，2008，134（5）：385-395.sphere，2012，22（1）：103-111.ics，2004，130（6）：646-655.2224-2229.1-5.

sphericalbodies［J］.PowderTechnology，2016，305（6）：89-98.particlesystems［J］.PowderTechnology，2017，313（2）：18-26.

—90—

［J］.PowderTechnology，2003，130：49-55.

［28］BOZKURTMG，FRATTAD，LIKOSWJ.Capillaryforcesbetweenequallysizedmovingglassbeads：anexper⁃［29］OLIVIERPitois，PASCALMoucheront，XAVIERChateau.Liquidbridgebetweentwomovingspheres：Anexper⁃

imentalstudyofviscosityeffects［J］.JournalofInterfaceScience，2000，231：26-31.imentalstudy［J］.CanadianGeotechnicalJournal，2017，54：1300-1309.

Experimentalstudyonliquidforceandprofilebetweentwomovingsphereparticlesunder

differentliquidcontentcondition

PUCheng，LIUFengyin，ZHANGZhao，CHENGJingxuan，ZHAOWei

（InstituteofGeotechnicalEngineering，Xi’anUniversityofTechnology，Xi’an

710048，China）

Abstract：Thestudyonliquidbridgeishelpfultorevealtheinnermechanismofwaterretentionpropertiesofunsaturatedsoil.Takingtheliquidbridgebetweenequalparticlesasobject，theforcevalueandprofilechangeofliquidbridgeunderthreeparticlediametersandsixliquidcontentsweremeasuredandrecordedbyusingNanoUTMandCCDcamera，andthetestresultswerecomparedwiththetheoreticalcalculation

results.Meanwhile，thesuppositionofdifferenttypesoftensilecurveswasproposed.Theresultsshowthatthetensilecurvecanbedividedintoascentstage，descentstageandrupturestage，respectively，andtheinitialprofileisvariousunderdifferentliquidcontentcondition，thependularoccursonlywhenthesatura⁃tionisextremelysmallandtheinfluenceofgravitycanbeneglected，withtherisingofliquidvolume，theingofparticlesizeandtheliquidcontent，whiletherupturedistanceincreasewiththeincreaseofliquidpointduringexperiment.

contentanddecreaseofparticlesize，theinfluenceofparticlesizeonthemaximumforceandtherupturedistanceismoresignificant.Finally，differentkindsofcurvesarecausedbythediscrepancyofstartingKeywords：unsaturatedsoil；hydrauliccharacteristic；liquidbridgetensile；experimentresearch

initialprofileturnintocylindricalandconvex.Inaddition，themaximumliquidforceincreasewiththeris⁃

（责任编辑：耿庆斋）

—91—