1. 引言

堆沙堡是沙滩游客的娱乐活动之一，沙堡工作的核心其实在于沙子和水之间的混合比例，太少的水无法粘结沙子，太多的水不能固定成形。阿姆斯特丹大学的Daniel Bonn通过实验发现最好的比例是99%的沙子和1%的水，然而水沙比1:99在实际生活中很难实现。

栾茂田等 [1] 基于等径颗粒间液桥的受力分析来建立液桥桥面的弯液面方程，由于考虑了固液接触角，其弯液面方程由3个几何条件联立确定，是一个高度非线性的方程组，不能直接得到其解析解，需要通过数值分析的方法求得解析解。张昭等 [2] 等人在液桥毛细力演化规律的探索中对于液桥体积也是采用迭代的方式计算。

本文在栾茂田等人的基础上进行改进，不考虑固液接触角，得到简化后的弯液面方程，再利用定积分求得液桥体积，以曲率半径为中间变量，得到液桥体积和液桥毛细力的关系，从而可以求得液桥毛细力形成和破裂时的水沙比，并进一步在该范围内求得最佳的水沙比。

本文通过探究颗粒间的液桥毛细力和液桥体积来确定最佳的水沙混合比，为Daniel Bonnie的实验提供一定的理论依据，并且给出现实情况下容易实现的最佳水沙混合比。从搭建沙堡的材料配比方面着手，使沙堡保存更为长久，以此增加沙滩度假游客的趣味性，使他们得到堆建沙堡成功的自豪感与满足感，吸引大家前去动手参与活动，提高动手能力，除此之外，堆建沙堡也是一项团体活动，可以在活动中促进大家合作，培养人际沟通能力与合作能力，通过研究出来的水沙混合比，也能给沙雕艺术材料的配比提供一定的理论基础。

2. 数据来源与模型假设

为简化模型，本文做出以下假设：

1) 每个沙子都被认为是大小相等的球体；

2) 认为与毛细力相比，重力可以忽略不计；

3) 不考虑沙子颗粒的浮力。

3. 最佳水沙混合比的模型建立

3.1. 毛细力的形成过程

一堆干沙子之间由于没有吸引力所以是不成形的，我们常常是将沙子与水混合，才能捏出形状。这是因为当沙子和水混合之后，附着在沙子颗粒间的水分子，其颈部表面存在着一种促使液面收缩的力，称之为表面张力。同时还存在着基质吸力，这两种力共同构成了液桥毛细力。液桥如图一所示一堆干沙子之间由于没有吸引力所以是不成形的，我们常常是将沙子与水混合，才能捏出形状。这是因为当沙子和水混合之后，附着在沙子颗粒间的水分子，其颈部表面存在着一种促使液面收缩的力，称之为表面张力。同时还存在着基质吸力，这两种力共同构成了液桥毛细力。液桥如图1所示。

图中的蓝色部分表示液桥，它可看作是由两个弧面与两个固体颗粒围成的。液桥中带有阴影部分的椭圆形为液桥颈部，其表面存在着一种与液面相切而与边界线相垂直的促使液面收缩的力，称之为表面张力，还存在促使颗粒之间相互吸引的基质吸力，这两个力共同构成液桥毛细力，以颈部为例，其毛细力如图2所示。

Figure 1. Liquid bridge shape

图1. 液桥形状图

Figure 2. Liquid Capillary Force Diagram

图2. 液桥毛细力图

由图2可知，水分子对沙子颗粒有相互吸引的作用。若水过多，则沙子的吸水量过大，会使其表面张力强度下降，沙子就撑不住彼此而坍软，这也是沙雕很注重排水的主要原因。若水过少，则沙子的吸水量多小，导致表面张力很小。

因此，寻求一个合适的水沙混合比是非常必要的，通过张昭等人的研究发现，液桥毛细力与液桥体积有着很大的关系，液桥体积最大和最小时分别对应着毛细力的形成和断裂状态，并且在这个范围内，液桥毛细力随液桥体积的增大呈先增后减的趋势，所以液桥毛细力势必会有一个最大值。液桥毛细力最大时对应的水沙比就是我们要寻找的最佳水沙比

3.2. 液桥毛细力建立

由gorge法可知，液桥毛细力是由作用在两个沙子之间液桥颈部截面内的基质吸力

φ 与其截面周上的表面张力

γ 共同产生，可得到液桥毛细力

f

q 的表达式如下：

(1)

由Young-Laplace方程可推导得到基质吸力

φ 和表面张力

γ 的关系，如下：

φ

=

γ

(

1

r

1

−

1

r

2

) (2)

可进一步得到液桥毛细力的表达式:

f

q

=

π

γ

r

2

(

r

1

+

r

2

r

1

) (3)

3.3. 液桥体积模型建立

由栾茂田等 [1] 人对液桥弯液面方程的建立和求解得到启发，如图3所示，将弯液面方程表示如下：

Figure 3. Derivation of meniscus equation

图3. 弯液面方程推导图

图中S区域为一个直角三角形，根据勾股定理可构造弯液面方程如下:

x

2

+

(

r

1

+

r

2

−

y

)

2

=

r

1

2 (4)

将等式移项去括号等变形处理，可得到如下解析式：

f

(

x

)

=

r

1

+

r

2

−

(

r

1

2

−

x

2

) (5)

得到弯液面方程后，可将液桥体积表示成弯液面绕x轴旋转一周与接触颗粒的体积之差，如图4所示：

Figure 4. Derivation of liquid bridge volume

图4. 液桥体积推导图

图中

d 是颗粒间距的一半，

β 是液体的填充角。则得到液桥体积如下：

V

q

=

∫

−

(

d

+

R

−

R

cos

β

)

d

+

R

−

R

cos

β

π

[

f

(

x

)

]

2

d

x

−

V

1

−

V

2 (6)

由数学知识，可得球冠的体积

V

1

,

V

2 为：

{

V

1

=

π

(

R

−

R

cos

β

)

(

3

(

R

sin

β

)

2

+

(

R

−

R

cos

β

)

2

)

6

V

2

=

V

1 (7)

由数学知识，两个圆相切，切线只有一条，且连接圆心和切点的连线一定是垂直于切线的，所以图中一定有点

O，

O

2，

P 三点共线，由三角形相似，则有

R

cos

β

R

+

d

=

R

R

+

r

1，

R

sin

β

r

1

+

r

2

=

R

R

+

r

1 (8)

再根据

(

R

cos

β

)

2

+

(

R

sin

β

)

2

=

R

2，可得到颗粒间距一半

d 与两个曲率半径的关系，如下：

(9)

所以联立以上各式可得到液桥体积的最终表达式如下：

V

q

=

(

2

x

3

r

1

(

r

2

+

r

1

2

)

*

R

1

2

+

x

3

3

R

3

−

π

R

2

(

R

1

−

x

3

)

(

2

r

2

−

2

x

3

R

1

+

R

1

2

)

3

)

R

1

3 (10)

其中

R

1

=

R

+

r

1，

r

=

r

1

+

r

2，

x

3

=

R

2

+

2

r

1

r

2，

R 为颗粒半径。可以看出，液桥体积如果不考虑固液接触角，最终只与曲率半径

r

1，

r

2 有关。

3.4. 含水量模型建立

查阅资料得到含水量与液体体积的关系 [3] 表达式如下：

V

q

=

8

π

ρ

s

R

3

3

n

c

ρ

ω (11)

式中，

ρ

s 为沙子密度；

ρ 为海水密度；

n

c 为配位数，常取值8，

ω 为含水量。

根据含水量，可以得到水沙的比例为：

p

ω

=

ω

1

−

ω (12)

结合(9)和(10)式，得到液桥体积和水沙比的关系为：

p

ω

=

3

n

c

ρ

V

q

8

π

ρ

s

R

3

−

3

n

c

ρ

V

q (13)

4. 最佳水沙混合比的液桥模型求解

这里我们选取了4种不同颗粒半径的沙子。

步骤一：利用matlab求解，以曲率半径

r

1，

r

2 为自变量，采用暴力搜索的方法，得到不同液桥体积下对应的毛细力，关系如图5所示：

Figure 5. The relationship between liquid bridge capillary force and liquid bridge volume

图5. 液桥毛细力与液桥体积的关系

毛细力随液桥体积呈现先增后减的趋势，与张昭等人的研究结论相符合。也进一步验证我们的推导是正确的。

步骤二：选取沙子半径

R

=

7

×

10

−

5，绘制其毛细力与含水量的关系曲线。

由图6可以看出，有3个关键点，分别是A (0.058, 0)，P (0.072, 210)，B (0.0921, 0)。很显然，A，B两点分别对应着液桥毛细力的形成和断裂。P点是液桥毛细力最大的状态。由该图可以看到，毛细力随含水量的增加是呈先增后减的一个趋势，当颗粒吸水达到饱和后，毛细力便达到了最大值，随着含水量的慢慢增多，会破坏液桥，导致毛细力下降。最后消失为0。

Figure 6. Capillary force and water content diagram

图6. 毛细力与含水量的关系图

通过式(10)简单计算，A点对应的水沙比为16.24，B点对应的水沙比为9.85，可近似地取水沙比的范围为[1:10:1:16]。P点对应的最佳含水量为0.072即最佳水沙比例为

P

w

=

1

:

12.

8

˙。取近似为1:13。

5. 结束语

本文建立的液桥模型在张昭 [2] 等人的基础上引入了弯液面方程，简化了液桥毛细力的计算。其中弯液面方程是根据栾茂田等人关于液桥体积计算的推导过程进行改进得到，不考虑固液接触角，并且引入了相似三角形得到线段的比例，从而简便计算公式。本文通过研究，在理论上提供了现实情况下容易实现的最佳水沙混合比以及水沙比的范围，从材料配比入手，使人们更加容易地堆建沙堡并保存长久，能在一定程度上增加人们的满足感和喜悦感，培养团队的合作能力。