韩志强，沈仲驰，唐明君

(四川师范大学 物理与电子工程学院，四川 成都 610101)

颗粒物质是指直径大于1 μm的大量离散的固体颗粒相互作用而组成的复杂体系，是常见的物态形式，也是地球上存在最多且最为人们熟悉的物质类型之一.由于颗粒物质体系处于静止堆积时，具有典型的非线性、自组织和能量耗散等特征，会产生诸如成拱、崩塌、粮仓效应等许多奇特的现象，这些现象的产生与日常生活生产息息相关.比如，当粮食堆积高度约大于2倍底面直径后，粮仓底面所受的压强不随粮食的增加而增加，而是趋于饱和值，这就是粮仓效应.粮仓效应由英国科学家Isaac Roberts在研究粮仓底面的压强时发现[1]，直到1895年，德国学者Janssen提出的连续介质模型巧妙地将颗粒物质空间应力求解的问题进行了简化, 很好地解释了“粮仓系统”中底部压力随填充高度的增加而逐渐趋于饱和的现象[2].直到如今，许多国家仍然采用Janssen的连续介质模型作为粮仓设计规范[3-5]的主要计算思想.

为了探究颗粒物质的粮仓效应和验证Janssen模型的正确性，有很多研究小组采用不同的测量方法测量了粮仓底部压力和侧壁压力.比如，早期采用的敲击法[6]，但由于敲击容器的力度很难控制，实验数据可重复性低.后来，王晶等采用间歇沉降法与连续沉降法对容器底部压力进行了测量[7]，该方法弥补了敲击法的不足，具有一定优势.但王晶等没有对侧壁压力进行探究，同时，颗粒与器壁的摩擦力未达到最大，导致所测的底部压力数据偏大.本实验通过间歇沉降法测量粮仓的底部压力，在此基础上，设计出了能连续测量侧壁压力的实验装置，并且保证了测量值在误差允许范围内具有可重复性.

1 理论计算

下面计算底面边长为a的正方形的立方体粮仓中的颗粒堆积情况，在Janssen模型中，假设颗粒在同一高度处的横向应力强度相等，且正比于纵向应力强度.

如图1所示，建立坐标系，水平方向为x轴、竖直方向为z轴方向.颗粒物质受力情况如图2所示，则在面积为S、厚度为dz的体积元上，力平衡的条件为

图1 坐标轴建立示意图

图2 颗粒物质受力示意图

-ρSgdz+Sdp=-τ,

(1)

其中，ρ为堆积密度，p为深度为z处的压强，τ为粮仓壁与颗粒物质间的摩擦力.由于纵向应力和横向应力是正比关系，摩擦力为

τ=4κμfpagdz，

(2)

κ为转换系数[8-10]，在实验中使用沉降法颗粒与容器壁之间的摩擦达到最大静摩擦力，所以μf为颗粒与容器壁摩擦系数.对上式进行积分可得深度为z处的压强为

(3)

又根据Janssen模型假设，侧向应力正比于正向应力

(4)

则用侧壁半径为R的圆来测量力面上的压力为

(5)

2 实验装置设计

实验器材：亚克力板、升降平台、压力传感器及显示屏、PVC板、钢尺、3 mm和5 mm陶碳粒.

本实验的实验装置设计如图3所示，由以下实验操作系统组成：

(a)正视图

a.粮仓主体部分，具体尺寸如图4所示，正方体粮仓的底边边长均为4 cm，高为30 cm，活板的长度长于其他3个面的高度.底部设置为边长略小于粮仓边长的方形活塞，这样可以使活塞上下移动时不会触碰容器的内壁，也避免颗粒掉落；

图4 方筒粮仓四壁设计图

b.升降支架部分，调节升降台的旋钮控制粮仓的底板和侧面的活动板的移动，用以测量不同高度和深度时，粮仓的底板正压力和侧壁压力的大小变化；

c.压力传感器部分，使用精度为0.01 g，量程为0～500 g的力传感器，将2个压力传感器的一端分别与侧壁上的小圆片和底部的方形活塞相连，另一端再与升降台上的固定装置相连，这样力传感器就能测量到小圆片和方形活塞所受到的力的大小.实验装置的实物图如图5所示.

图5 实验装置实物图

1)底部压力测量

由于在理论中假设侧面所受静摩擦力为最大静摩擦力，而将颗粒注入方筒中后，侧面对颗粒的摩擦力不是最大静摩擦力，故直接测量的底面压力并不满足此模型，所以本文采用沉降法测量底部的压力.如图3所示，在固定主体部分的位置时，通过缓慢调节升降台，使底面支撑体下降，即侧面静摩擦力达到最大静摩擦力.然后使支撑体下降很小的高度，此时堆积的颗粒将有下降的趋势，故此时产生的摩擦力接近于最大静摩擦力，在整体堆积处于Janssen模型的状态后，使用压力传感器测量颗粒对粮仓产生的底部压力.

2)侧壁压力测量

在研究粮仓效应时，侧壁压力的精确测量是整个实验和验证模型的正确性的关键.已有采用在侧壁的不同高度上放置不同的压力测量装置，测量侧壁的压力[8]，这种测量方式的缺点在于，只能测量固定高度点的侧壁压强，不能连续测量.为了测量更多高度的侧壁压力，本文设计了能连续测量不同高度侧壁压力的实验装置.如图3所示，在装置的侧面，粮仓的侧壁和测力装置作为整体.具体做法是在侧壁中切割出圆形区域，在此圆形区域，制作不与侧壁接触的圆片，将圆片与测力装置相连.调整该组件的高度，可以连续测量出不同深度时粮仓的侧壁压力值.该实验装置不仅能够实现沉降法测量底部颗粒压力，还能够连续地测量侧壁压力.

3 实验数据采集

3.1 不同堆积高度下的底部压力测量

实验测量时，将压力传感器与数码显示器连接.颗粒半径分别为3 mm和5 mm的均匀球形陶粒，成分主要为高岭土、长石和黏土，密度为1.22 g/cm3.将陶粒注入方筒中，达到实验所需高度，将底部支撑升降台旋转固定小角度，颗粒达到极限应力状态时，测量出陶粒的堆积高度和粮仓底部正压力的对应关系.

如图6所示，测量出当陶粒半径分别为3 mm和5 mm时，陶粒的堆积高度和粮仓底面所受压力大小的关系，图中实心圆点表示在同一位置进行10多次实验测量的底部压力的平均值，实线表示的是式(3)使用Matlab软件计算模拟出的理论曲线.其中3 mm陶粒的κ为0.41，5 mm陶粒的κ为0.56，在弹性理论计算的范围(0.33～2.50)内[11]，且与其粒径变化的规律相同[8].从图6中可以看出，随着陶粒高度的增加，粮仓的底部压力增大，当高度增加到一定值时候，底部压力的大小趋于平稳.从图6(a)中可以看出，半径为3 mm的陶粒高度达到大约12 cm后，粮仓底部压力趋于平稳，达到定值.如图6(b)中所示，当半径为5 mm的陶粒高度达到大约10 cm后，粮仓底部压力趋于平稳.即，当堆积高度达到方仓底边边长(4 cm)约2倍以上以后，粮仓的底部压力达到饱和.从图6的测量结果可以看出，对于不同粒径的陶粒，粮仓的底部压力达到饱和的陶粒高度不同，这是由于2种半径不同的陶粒具有不同的转向系数和堆积密度，所以，造成了同样深度下的底部压力不同.从图6中还可以看出，实验测量值与理论计算模拟值存在偏差，这是由于每次测量时，难以形成完全相同的极限应力状态造成的，但是，理论曲线和实验测量值的变化趋势都相同.

(a)3 mm

3.2 不同高度的侧壁压力测量

实验测量粮仓侧壁压力时，先将测力面置于装置能达到的最低处，在保持底部正压力为Janssen模型的状态下，旋转控制侧壁升降台的旋钮，使侧壁装置整体上升，通过压力传感器测量，读出相同堆积高度的情况下，不同测量深度处的侧壁压力大小，2种不同的陶粒不同深度的实验测量结果如图7所示.

图7测量出了当陶粒半径分别为3 mm和5 mm时，陶粒在相同堆积高度16.1 cm下不同深度的侧壁压力传感器所受压力的大小.图7中的实心圆点表示在多次侧壁整体上升过程中不同深度的测量值的平均值，虽然圆片有一定的面积，但仍能代表这个深度的侧壁压力.图7中实线表示的是根据式(5)，用Matlab积分计算模拟的理论曲线.从图7中可以看出，在相同堆积高度的粮仓中，随着测量深度的增加，侧壁压力逐渐增大，随着深度的继续增加，侧壁压力的增长趋势越来越缓，最后趋于定值.理论模拟值和实验测量值符合得比较好.

(a)3 mm

4 结论与讨论

通过本文的讨论得出：随着方形粮仓中倒入颗粒的高度的增加，粮仓底部的压力会逐步增大，最后会趋于某定值；当陶粒的堆积高度一定时，随着测量深度的增加，方形粮仓的侧壁压力逐渐增加，且会趋于定值，出现粮仓效应.

实验使用自制的实验装置能明显地演示出粮仓效应，并能较精确地测量底部正压力和侧壁压力，验证了Janssen模型的正确性.同时，该实验装置能够连续测量侧壁压力，能较好地反映侧壁压力变化的趋势，验证了产生粮仓效应的原因，即，由于颗粒物质间的相互作用，重力方向的力被分解到了水平方向，粮仓的侧壁支撑了颗粒的部分重力，使得粮仓的底部压力趋于饱和.但由于在改变侧壁装置整体高度时装置的晃动很容易改变原本的堆积应力状态，所以侧壁压力的测量误差偏大，但并不影响对整个趋势变化的判断.