颗粒柱坍塌运动与堆积特性的研究综述

来志强，江恩慧，赵连军，周 伟，田文祥，马 刚

(1. 黄河水利科学研究院 水利部黄河泥沙重点实验室， 郑州 450003；2. 中国科学院 地理科学与 资源研究所，北京 100101；3. 武汉大学 水资源与水电工程科学国家重点实验室， 武汉 430072)

颗粒是自然界中常见的物质存在形式，如泥沙、岩屑石块、矿石、种子、药丸等.由于其独特的物质形态，颗粒物质的动力学行为有时类似固体，有时类似流体，在时间和空间尺度上相互转化，具有多面性和复杂性.例如，岩质边坡坍塌运动、山体滑坡的发生，就是颗粒物质由固态转化为流态的复杂流变现象.此外，颗粒集合体坍塌运动与堆积现象在化工、农业和医药等多个领域均有涉及.近年来，国内外许多学者通过颗粒柱坍塌运动模型研究颗粒物质的坍塌运动与堆积特性.颗粒柱坍塌运动模型最早由Lajeunesse等[1]和Lube等[2]于2004年提出.由于此模型与众多工程实际问题相关，如在地质灾害或工业制造领域中，模型边界挡墙瞬间移除导致的颗粒柱坍塌运动可概化模拟岩质边坡坍塌、农业、化工、制药等工业中物料颗粒在筒仓或通道中的输运等问题，当挡墙缓慢移动时则可模拟土体的准静态失稳过程.近年来，颗粒柱坍塌运动模型引起了诸多领域学者们的广泛关注.

根据颗粒柱的初始形状和边界挡墙的移除方式，目前颗粒柱坍塌运动模型主要可以分为长方体单向运动模型[3-4]以及圆柱体对称运动模型[5-6].尽管两者的几何形状与运动维度较为不同，但是其运动机制相同，因此获得的结论具有共性.目前，颗粒柱坍塌运动模型还被应用于验证颗粒物质动力学特性描述理论和数值模拟方法，如浅层流模型[7-8]、基于动力学理论或平均深度理论的连续体模型[9-10]、双颗粒温度热力学理论[11]、颗粒体系总摩擦因数连续性理论[12]、有限差分连续体模型[13]、考虑流变特性的Navier-Stokes连续体模型[14]、离散单元法(DEM)[15-16]、粒子有限元法(PFEM)[17-18]、有限元法(FEM)[19]、光滑粒子流体动力学(SPH)法[20]、基于移动粒子半隐式法[21]和物质点法(MPM)[22]等.由此可见，颗粒柱坍塌运动模型的应用十分广泛，已成为研究颗粒物质动力学特性的经典模型.

本文从影响颗粒物质运动与堆积特性的主要因素，如颗粒柱初始空间形态特征、颗粒基本物理特性、模型边界及环境条件等方面系统地概述了近年来国内外有关颗粒柱坍塌运动与堆积特性的研究进展，讨论了不同因素对颗粒柱动力与堆积特性的影响规律及作用机理，在总结已有成果的基础上，结合近年来工程实践遇到的科学问题，提出目前研究成果存在的争议性及未来研究应集中的方向，旨在提高对颗粒物质运动堆积基础问题的认识.

1 颗粒柱初始空间形态特征影响

颗粒柱初始几何形态一般为长方体或圆柱体，目前学者们从颗粒柱初始宽高比和孔隙率两方面对颗粒柱的初始空间形态特征进行描述，研究不同条件下颗粒柱坍塌运动与堆积特性的变化规律.

1.1 初始宽高比

颗粒柱的初始宽高比a决定了颗粒柱的运动机制能够显著影响其运动与堆积特性，这是目前学者们达成的一致结论.定义a=H0/L0(长方体模型)、a=H0/R0(圆柱体模型)，模型示意图如图1所示.其中：H0为颗粒柱初始高度；L0为颗粒柱初始长度；b0为颗粒柱初始宽度；R0为颗粒柱初始半径；L∞为长方体颗粒柱的运动距离；R∞为圆柱体颗粒柱的运动距离. 对于由单一粒径颗粒组成的长方体(圆柱体)颗粒柱的运动距离L∞(R∞)无量纲参数[L]([R])与a存在如下关系：

图1 颗粒柱坍塌运动模型示意图Fig.1 Diagram of granular column collapse movement model

(1)

最大堆积高度H∞的无量纲参数[H]与L0、a存在如下关系：

(2)

式中：c1、aL、c2、aR、c3、aH为与颗粒基本物理特性和模型边界条件有关的常数.

由式(1)和(2)可知，[L]、[R]、[H]与a之间存在线性和幂律的分段函数关系.一些学者[1-66]通过物理试验和数值模拟获得的c1、aL、c2、aR、c3、aH的分布值，如表1所示.由表1可知，当颗粒的基本物理特性和模型边界条件发生变化时，c1、c2受其影响不大且均为正数.这表明当a>aL(a>aR)时，[L]([R])与a之间的幂函数正相关关系具有较强的普适性.一些文献没有给全c3的分布值，但是从已有的成果可以看出，当a>aH时，[H]与a之间也存在正相关关系.

由于相同条件下，圆柱体和长方体的颗粒柱运动机制相同，所以下文以长方体颗粒柱为例进行分析.a决定了颗粒柱的运动机制，高、低颗粒柱不同坍塌运动机制示意图如图2所示，其中α为破坏面倾角.低颗粒柱(aaL)，其破坏面(倾角也为α)左下方部位为静止区域， 此处颗粒几乎不发生运动，位于破坏面右上方颗粒则产生坍塌运动，颗粒数目所占比例较大，此过程惯性力起主导作用.Utili等[3]和Lajeunesse等[23]认为α与颗粒物质内摩擦角φ满足Rankine理论[24]即α=45°+φ/2，而Fern等[25]却提出α与φ相等，最终Xu等[26]在物理试验中采用粒子图像测速(PIV)技术直接测得破坏面倾角α=45°+φ/2，证明了前者结论的正确性.

图2 高、低颗粒柱的不同坍塌运动机制Fig.2 Different collapse movement regimes for granular columns at small and large aspect ratios

表1 已有文献中c1、aL、c2、aR、c3、aH参数取值分布Tab.1 Distribution of c1,aL,c2,aR,c3 and aH parameter values in previous research

在低颗粒柱坍塌运动过程中，垂向应力与颗粒柱高度的相关性较大，越靠近颗粒柱底部位置的垂向应力，其分布值越大，而对于高颗粒柱，垂向应力呈“楔形”分布[27].Xu等[26]在物理试验中采用接触式压力传感器获得了低颗粒柱(a=1.8)运动过程中底部压力的分布规律，发现初始状态下小压力(小于压力平均值)的概率密度方程呈幂函数分布形式，大压力(大于平均值)的概率密度方程则满足Gaussian分布，随着坍塌运动的进行，颗粒柱底部压力分布整体减小.然而，关于高颗粒柱底部压力随其运动的变化规律、颗粒受力与颗粒柱运动机制之间存在的内在联系，目前的研究仍十分不足.

此外，a显著影响颗粒柱坍塌运动中能量的演化规律.Zenit[28]采用DEM模拟了a=0.3, 3.0, 12.4时，二维长方形颗粒柱坍塌的运动过程，发现a越大，则颗粒柱相对势能的变化量与相对动能峰值越大，堆积形态的变化也随之越大.Utili等[3]利用DEM获得了a=0.9，3.3，5.9，9.3时，三维长方体颗粒柱能量和动量通量的变化规律，研究结果表明颗粒柱坍塌运动的总能耗以颗粒摩擦耗能为主.当a增大时，系统总能耗随之增大，相比于平动动能，颗粒系统的转动动能很小，可以忽略.此外，还建立了表征颗粒体动量通量以及瞬时冲击能量通量的无量纲参数，为离散颗粒体系冲击特性的评估提供了理论基础.由于Utili等[3]模拟时采用理想化的圆球颗粒，通过对颗粒施加抗转动系数人为控制其转动能力，但能否忽略颗粒转动特性的影响仍有待进一步考证.

1.2 初始孔隙率

由于颗粒的离散特性，颗粒在一定空间内堆积时总会存有孔隙，如图3所示.颗粒柱初始孔隙率n为孔隙体积占颗粒柱总体积(H0L0)的比例.Kermani等[5]利用DEM分析了n(0.40≤n≤0.45)对颗粒柱坍塌运动的影响，发现n对颗粒柱运动距离的影响不大，但对堆积高度的影响显著，n越小，初始颗粒柱密实度越大，则堆积高度越大.Kermani等[5]认为在剪切作用下密实颗粒柱会产生“剪胀”现象，松散颗粒柱则产生“剪缩”现象，两者之间存在一个临界孔隙率值，颗粒柱运动过程中经历长时间剪切作用后内部孔隙结构达到此临界状态，因此，水平运动距离受初始孔隙率影响微弱，但文献[5]没有解释为何n会对颗粒柱的堆积高度产生影响.Fern等[22]采用MPM分析了n(0.50≤n≤0.90)对颗粒柱运动特性的影响，发现n越大，颗粒柱发生运动部位的体积越大，势能转化量越大，水平运动距离就越远，相同n条件下a越小，则n的影响越明显.上述两种结论不同的原因可能在于文献[5]模拟n的分布范围过小，获得的结论代表性不强.

图3 颗粒柱内部孔隙示意图Fig.3 Diagrammatic sketch of internal porosity of granular columns

当颗粒柱在液体环境中作坍塌运动时，n的影响机制将发生改变.Rondon等[29]进行了不同n(0.38≤n≤0.45)条件下颗粒柱在液体环境中的坍塌运动物理试验，发现n对颗粒柱堆积形态的影响显著.密实颗粒柱(0.38≤n<0.42)运动时产生“剪胀”现象，部分液体被吸入颗粒柱内部，在黏滞曳力作用下颗粒柱局部运动速度减慢，运动距离减小，而n对堆积坡度没有影响，均为颗粒材料的休止角度；松散颗粒柱(0.42≤n≤0.45)运动时，将内部孔隙中的水体排出，形成颗粒柱局部流态化，使其运动速度加快，最终形成的堆积体较长，斜面坡度随着n的增大而增大.Yang等[30]从细观角度解释了上述现象的形成原因，认为密实颗粒柱密集树枝状力链网络阻碍了颗粒体的剪切滑动，导致颗粒柱进行蠕变式缓慢运动，而松散颗粒柱则更易形成“滑水”式快速运动的模式，大大降低了摩擦阻抗.Lee[31]的物理试验表明了颗粒间孔隙压力对密实、松散颗粒柱运动的影响机制不同.对于密实颗粒柱，颗粒间易形成负的超孔隙水压力，增大颗粒间的摩擦，阻碍颗粒柱坍塌运动；对于松散颗粒柱，颗粒间易形成较大的孔隙水压力，减小颗粒间的摩擦，促进颗粒柱坍塌运动.

2 颗粒基本物理特性影响

目前，关于颗粒基本物理特性对颗粒柱坍塌运动的影响研究主要集中在颗粒粒径、颗粒形状、颗粒摩擦、碰撞及变形特性、颗粒破碎与颗粒潮湿等方面.不同的颗粒基本物理特性对颗粒柱坍塌运动的影响规律不同，作用机制也不同.

2.1 颗粒粒径

相对于具有不同H0、L0的颗粒柱，同一颗粒粒径d的影响效应不同，因此采用无量纲参数L0/d反映颗粒粒径与颗粒柱尺寸的相对大小.现有研究指出当L0/d至少为10时，颗粒集合体形成的剪切带和力链才能够被精准描述[32-33]，因此颗粒柱中的颗粒粒径d应小于L0/10.

关于由单一粒径颗粒组成的颗粒柱，Lube等[2](12.7≤L0/d≤303)和Artoni等[34](14≤L0/d≤35)的物理试验表明d对运动距离和堆积形态没有影响.而Gabrieli等[35]物理试验和DEM数值模拟却发现d(14≤L0/d≤35)越大，则颗粒柱运动距离越远，堆积角度越低，当颗粒处于潮湿状态时，d的影响更为显著.Huang等[36]的物理试验结果表明当颗粒粒径增大时，颗粒柱运动距离减少，但试验所采用的不同粒径的颗粒密度有所不同，因此无法排除颗粒密度对上述规律的影响.针对上述相悖的结论，Cabrera等[37]采用DEM研究了不同a(0.25≤a≤16)条件下，L0/d(10≤L0/d≤200)对颗粒柱运动堆积特性的影响规律，发现在特定范围内d才会对颗粒柱的运动时间和运动距离产生影响.当低颗粒柱L0/d≥75或高颗粒柱L0/d≥50时，d对颗粒柱坍塌运动没有影响.

此外，Phillips等[38]开展了由两种粒径颗粒组成颗粒柱的运动特性试验研究，发现二元粒径颗粒柱的运动距离比仅由单一粒径颗粒组成的颗粒柱运动距离远.当小粒径颗粒质量占颗粒体系总质量的比例ψ约为0.3时，颗粒柱运动距离达到最远；当ψ≤0.2时，小颗粒通过自身旋转促进大颗粒运动，随着ψ的增大，越来越多的大颗粒被小颗粒包裹，使其摩擦减少，整体运动距离增大；当ψ≥0.7时，由于相同体积的小颗粒比大颗粒产生的接触面积多，随着ψ的增大，颗粒摩擦接触增多、能耗增大，导致运动距离减小；当0.2<ψ<0.7时，颗粒系统处于以上两种作用机制的过渡阶段.此外，颗粒密度对上述规律没有影响.Degaetano等[39]通过物理试验也获得了类似规律，当ψ=0.5时，颗粒柱运动距离达到最远，颗粒的初始位置越高，其最终堆积位置越远，且大粒径颗粒比小粒径颗粒运动距离远.

一些学者还研究了由多种粒径颗粒组成的颗粒柱运动特性.Vallejo等[40]采用分形维数D[41]描述了颗粒柱中小粒径颗粒的含量，D越大表示小粒径颗粒含量越多.试验发现当D增大时，颗粒柱整体运动距离增大.其原因在于当小粒径颗粒含量增多时，颗粒整体旋转运动增强，颗粒之间剪切摩擦减弱.随后，Lai等[42]利用DEM进一步探讨了D的影响发现，随着D的增大，颗粒体断面水平运动速度由线性分布形式逐渐转变为幂律分布形式，颗粒与底板之间形成边界层效应，位于此处的颗粒受到强烈的剪切作用而获得较大的运动速度，颗粒粒径越小，受到的颗粒接触作用力越大，导致其水平运动速度越大.研究成果从细观角度揭示了小粒径颗粒对颗粒柱运动特性的“润滑”作用机制.

2.2 颗粒形状

颗粒形状提供了颗粒之间的咬合作用力，进而导致颗粒具有抗转动的特性.Tapia-McClung等[43]通过DEM将不同数目的圆球子颗粒组合成类似棒状的细长颗粒，采用子颗粒的数目来表征颗粒的长度，模拟了其组成颗粒柱坍塌的运动过程，结果发现[L]、[H]与a均存在如式(1)和(2)所示的函数关系，颗粒长度对颗粒柱运动过程中的能量演化、等效摩擦力、最终运动距离与堆积高度影响很小，上述结论与以往颗粒形状显著影响离散颗粒堆积形态和力学特性的认识相悖[44-45]，而文献[43]也没有从本质上解释清楚上述现象的发生机制.Owen等[46]分析了由复杂形状块体组成的颗粒柱运动的堆积特性，发现块体颗粒形状越不规则，颗粒柱运动距离越小，堆积高度越大，但没有明确指出颗粒形状的影响机制.Trepanier等[47]利用物理试验研究了由异形长颗粒棒组成的颗粒柱坍塌运动与堆积特性，发现当颗粒棒长度Lp与颗粒粒径d比值满足Lp/d>24时，颗粒柱保持稳定状态不发生坍塌；当Lp/d≤24时，存在颗粒柱初始高度值上下限临界值HU和HL，当H0>HU时，颗粒柱必定失稳坍塌；当H0

有些学者在进行颗粒柱DEM数值模拟时，向圆球施加抗转动系数以考虑颗粒不规则形状对颗粒转动的阻抗效应.Cleary等[49]和Kermani等[5]的研究均表明当颗粒抗转动特性增强时，颗粒转动减弱，颗粒柱运动距离减小，堆积高度增加.但抗转动系数无法真实反映颗粒形状之间的咬合互锁效应，施加抗转动系数的圆球颗粒柱运动距离往往比真实不规则形状的颗粒运动距离大很多.

2.3 颗粒摩擦碰撞及其变形特性

有关颗粒摩擦特性对颗粒柱运动距离与堆积形态影响的研究存在一些互相矛盾的结论.有些学者认为颗粒摩擦特性没有影响，如基于DEM数值模拟方法，Zenit[28]、Owen等[46]和Lo等[50]分别考察了颗粒摩擦因数μp(0.2≤μp≤1.0)的影响，均认为μp对颗粒柱的运动堆积特征影响不大.一些学者则认为μp显著影响颗粒柱的运动特性.Staron等[15,51]的DEM数值模拟结果表明，当μp(0.01≤μp≤2.00)增大时，颗粒体与滑槽底部逐渐形成静止层，导致颗粒柱运动距离减小.Cleary等[49]发现当μp(0.2≤μp≤0.5)增大时，颗粒柱运动距离减小，堆积高度增大.相关的PFEM[17]、FEM[19]和MPM[27]数值模拟研究均得到了类似结论.μp对颗粒柱运动过程中的摩擦能耗起了决定作用，特别是对于以摩擦力为主控作用力的低颗粒柱运动过程，μp的影响会更加显著，因此需结合颗粒柱初始宽高比a及其运动机制进一步探讨μp的影响规律.

图4 颗粒恢复系数示意图Fig.4 Diagrammatic sketch of restitution coefficient of particles

颗粒刚度表征颗粒在外荷载作用时抵抗自身弹性变形的能力，如图5所示，其中：Fn为颗粒法向荷载；Fs为颗粒切向荷载；kn为颗粒法向刚度；ks为颗粒切向刚度；un=Fn/kn为颗粒法向变形；us=Fs/ks为颗粒切向变形.当颗粒在法向荷载Fn和切向荷载Fs作用下，颗粒刚度(法向刚度kn和切向刚度ks)越大，则颗粒变形(法向变形un和切向变形us)越小，DEM数值模拟时间步长越小，计算成本越大.目前已有研究均认为颗粒刚度对颗粒柱的运动堆积特性影响不大，但颗粒刚度越小，DEM计算完成所需的时间就越短[16,46,52]，因此许多学者在DEM数值模拟时选取的颗粒刚度一般较实际值小，然而颗粒模拟刚度可比实际刚度具体小几个数量级，此方面的研究仍有待进行.

图5 颗粒刚度示意图Fig.5 Diagrammatic sketch of particle stiffness

2.4 颗粒破碎与颗粒潮湿

自然界岩质边坡内部存在很多裂隙，这些边坡在坍塌运动时不断发生破碎，导致其运动距离显著提高[53].此外，人工制造的脆性材料如玻璃、陶瓷、混凝土等，由于其材料内部存在诸多微裂纹，在冲击荷载作用下易发生破碎[54].以此为研究背景，Langlois等[4]采用DEM将颗粒柱中颗粒粘结为整体，并赋予不同的初始粘结强度以反映不同的裂隙发展程度，研究颗粒破碎对颗粒柱坍塌运动特性的影响.研究发现：当颗粒柱初始粘结强度降低(裂隙发展程度增强)时，颗粒破碎效应增强，颗粒柱滑动距离增大，堆积体表层光滑度提高，堆积形态不规则性降低；破碎后的颗粒累计体积分数比曲线约在25%和75%两处产生峰值，即呈双峰特征，同时[L]与a之间仍满足式(1)和(2)的函数关系.图6定性地给出了不同初始粘结强度下颗粒柱堆积的形态示意图.

图6 不同初始粘结强度下，颗粒柱的堆积形态示意图[4]Fig.6 Diagrammatic sketch of accumulation characteristics of granular columns at different initial bond strengths[4]

当颗粒之间存在少量液体(潮湿颗粒)时，颗粒的物理力学性质(如毛细力Fc的存在，如图7所示)将发生改变，进而影响颗粒集合体的动力学特性[55].毛细力Fc的表达式如下：

(3)

图7 潮湿颗粒间毛细力示意图Fig.7 Diagrammatic sketch of capillary force between wet particles

液体改变颗粒柱动力学特性的机制在于：少量液体的存在将湿润颗粒表面积，导致颗粒之间产生毛细力，减弱颗粒柱的离散性，对颗粒柱的宏观运动起阻碍作用.随着w的增大，颗粒之间的控制力由摩擦力逐渐转变为毛细力；颗粒粒径决定了其表面积和质量的大小，颗粒重力、毛细力与其呈正相关关系，而颗粒重力对颗粒柱坍塌运动起到了促进作用，因此应采用颗粒毛细力与颗粒重力的比值描述w与d的综合作用.

3 模型边界及环境条件影响

颗粒柱边界一般采用刚性挡墙对其进行约束，刚性挡墙约束情况及其摩擦特性对颗粒柱的运动堆积过程具有一定的影响.此外，不同学者根据不同的研究背景分析了颗粒柱在气体掺入导致的流态化及水体环境作用下的运动堆积过程，取得了丰硕的研究成果.

3.1 挡墙约束条件及摩擦特性

挡墙是颗粒柱运动的边界条件，其不同的移动方式决定了颗粒柱拟静态或坍塌运动的模式.Owen等[46]和Mériaux[57]在物理试验和DEM数值模拟中水平缓慢移动挡墙以模拟土坡的拟静态失稳过程，发现此过程并不满足Mohr-Coulomb破坏准则，系统能耗以摩擦能耗为主，颗粒碰撞能耗较少，挡墙移动速度v(1.8≤L0/v≤26.1)对颗粒柱的运动堆积特征没有影响.Zhu等[58]的物理试验发现二元粒径颗粒组成的颗粒柱拟静态运动过程中存在一个边界层，位于边界层以上的颗粒分层滑动，位于边界层以下的颗粒保持静止.Kermani等[59]的研究结果表明，采用基于连续介质理论的光滑粒子流体动力学法和基于离散介质理论的离散单元法均可准确模拟颗粒柱拟静态运动过程，但后者对堆积尖角与颗粒前端模拟更精确，低颗粒柱堆积形态趋向于上端被截断的圆锥体，高颗粒柱堆积形态更相似于圆锥体.

挡墙的几何特征与位置对颗粒柱的坍塌运动过程也存在一定的影响.Girolami等[60]采用DEM研究了颗粒柱在无侧限墙体约束下的运动堆积规律，发现颗粒柱运动堆积过程包括3个阶段，分别为自由落体、扩展运动和逐渐静止.[L]与a之间仍满足线性及幂律函数关系，当a≤3时，低颗粒柱以摩擦力为主控力，运动速度较慢，运动过程满足Mohr-Coulomb准则，形成堆积体的坡度与颗粒材料休止角度在相同数量级，堆积形态受侧限墙体的影响不大；当a≥3时，相比于侧限挡墙存在情况，无侧限约束下颗粒柱逐渐静止阶段耗时长，最终形成堆积体沿x方向长度较长，堆积坡度约为颗粒材料休止角度的一半，y方向的动量将会更大比例地传递至x方向.

有些学者还探讨了z方向挡墙宽度b0对颗粒柱坍塌运动与堆积过程的影响.Balmforth等[61]分别进行了极宽、极窄(L0/b0分别为0.1、2.0)条件下的颗粒柱运动堆积物理试验，发现当b0增大时，[L]、[H]与a拟合得到的幂函数指数c1、c2均增大.Lacaze等[62]也开展了此方面的物理试验研究，发现当b0越接近d时，颗粒柱局部结晶化现象越显著，而b0越接近2d时，颗粒间的阻塞现象越明显.上述两种现象均会对试验结果产生影响，因此获得的试验规律不具有共性.文献[62]还指出当b0=1.2d时，结晶化和阻塞效应最小，试验能够得到颗粒柱的典型运动堆积特征，随后利用DEM重现了该物理试验过程，发现内部静止区域上部的移动层速度满足线性分布形式.

农业工程中，经常将固体化肥或种子等颗粒堆积于以特定频率旋转的离心盘上使其分撒于指定范围.以此为研究背景，Warnett等[63-64]利用物理试验研究了圆柱体颗粒柱在底板旋转状态下的坍塌运动过程.存在临界转动频率fr，当转动频率f

(4)

式中：C为与颗粒基本物理特性、模型边界条件有关的常数.

Nikooei等[65]初步开展了可侵蚀床面上的颗粒柱坍塌运动堆积特性模拟研究，相比于同等条件下的刚性床面，当颗粒柱在可侵蚀床面上运动时更多的势能转化为水平动能，进而促进颗粒柱的水平运动距离.尽管可侵蚀床面降低了颗粒柱的最大运动速度，但颗粒柱整体质量的增加使得颗粒柱的水平动量和冲击破坏力均有所增大.有关颗粒柱在可侵蚀床面上的运动堆积机制研究才刚刚引起学者们的关注，此方面的研究仍需加强.

有关挡墙摩擦特性对颗粒柱运动堆积特征的影响规律，目前研究存在相互矛盾的结论.Crosta等[19]的FEM数值模拟结果表明，挡墙摩擦因数μw对颗粒柱运动过程的影响很小，但当颗粒柱位于可侵蚀底板时，颗粒柱运动距离减小，此时μw越大，颗粒柱运动距离越小，堆积体前端越薄.Lube等[2,6]和Lajeunesse等[1,23]通过物理试验也获得了μw对颗粒柱运动堆积特性影响不大的结论.而Zhang等[17]的PFEM数值模拟结果则表明，当μw(0.36≤μw≤0.58)增大时，颗粒柱的运动距离减小，堆积高度增大，且a越大，μw的影响越显著.因此有必要在μw较大的分布范围内，结合其他因素如a、μp、d等影响，细化μw对颗粒柱运动堆积特征的影响规律，以获得适用性更为广泛的结论.

3.2 气体与水体作用

滑坡碎屑流运动过程中极易将其下方及前方的压缩气体裹入，从而产生流态化现象.流态化现象是滑坡碎屑流高速远程特性形成的主要原因.以此为研究背景，Roche等[66]从颗粒柱底部注入空气使其内部产生气体孔隙压力，以减少颗粒间由摩擦效应产生的流态化现象，随后分析了掺气后颗粒体运动堆积过程，发现流态化现象对大颗粒(d=330 μm)组成的颗粒柱运动影响很小.这是因为气体孔隙压力在其运动后迅速衰减，颗粒之间的摩擦效应重新得到恢复.而小颗粒(d=75 μm)组成的颗粒柱在流态化效应影响下的运动距离增大，这是由于运动过程中颗粒柱内部孔隙压力始终保持在较大值，颗粒之间的摩擦阻力较小，颗粒柱的运动距离无量纲参数[R]=λ1aλ2(λ1、λ2为常数)，其中λ1随颗粒摩擦因数的增大而增大，λ2则不受颗粒摩擦因数的影响，但随着a的增大而减小.由于流态化现象产生的振动波，导致颗粒柱最终的堆积形态呈现墨西哥帽型.

颗粒柱在有水环境中的坍塌运动机制与其在无水环境中不同，主要区别在于颗粒柱在无水环境中运动时前端松散，总是存在远离颗粒主体的散落颗粒；而在有水环境中颗粒柱开始运动时，其右上部开始运动的部位会形成涡旋，随后沿颗粒体的坍塌方向运动，颗粒体运动前端在水体作用下变得较为密实，没有远离颗粒主体的散落颗粒，水体的存在阻碍了颗粒体的坍塌运动，使其运动距离降低，如图8所示.以海底滑坡为研究背景，景路等[67]和Jing等[52,68]采用流固耦合DEM-CFD数值模拟方法探讨了低(a=1)和高(a=8)颗粒柱在水体环境中的坍塌运动过程.低颗粒柱的运动过程与上述过程相似；高颗粒柱的水中坍塌运动过程较为复杂.首先，位于上表层的颗粒主要以自由落体的形式运动，水平位移较小，涡旋运动不明显；随后，颗粒开始堆积并进行以剪切摩擦为主的运动，颗粒体前端变得厚实，同时大量水体涡旋开始发育，颗粒体在其作用下运动速度变慢，并逐渐形成平铺型堆积体；在涡旋冲蚀作用下，表层松散颗粒被带走，最终形成的堆积体表层凹凸不平，其中大部分初始势能发生转化，部分能量被水体黏性所消耗.相比于相同条件下的无水环境，颗粒体表面离散介质受水流影响较大.此外，当流体的斯托克斯数分布在0.06～6 430、颗粒与流体密度比分布在1.63～46.99时，颗粒柱在流体中的坍塌运动过程可分3个阶段，分别为自由下落阶段、流体惯性力作用阶段和黏滞力作用阶段.随着流体的斯托克斯数的减少，颗粒柱坍塌运动机制由颗粒滑移主导逐渐转变为颗粒悬浮主导.

图8 颗粒柱在无水、有水环境中坍塌运动过程示意图[67]Fig.8 Diagrammatic sketch of movement and accumulation process for dry and underwater granular columns[67]

4 结论

颗粒柱坍塌运动与堆积特性的研究是认识颗粒物质动力学特性的基础性工作，不同领域的学者从很多影响因素方面开展了大量研究.通过对相关国内外研究进行归纳，总结了颗粒柱初始形态特征、颗粒基本物理特性、模型边界及环境条件等因素对颗粒柱运动堆积特性的影响规律与作用机制，其中一些研究结论已经形成统一的认识.针对目前的研究不足和存有分歧的结论，以下几个方面仍需进一步探索.

(1) 有关颗粒柱初始孔隙率、颗粒摩擦因数、挡墙摩擦因数对颗粒柱坍塌运动与堆积特性的影响规律及作用机制还没有形成统一的认识.目前，不同学者得到的研究结论不同甚至相悖，其原因可能在于所探讨的影响因素分布范围小，数值模拟研究受限于自身的计算原理与本构模型，分析其影响规律时没有充分结合颗粒柱初始几何尺度和其他颗粒的物理特性，关于此方面还需深入研究.

(2) 有关颗粒柱坍塌运动过程中内部力链与底部压力分布随初始宽高比变化规律的研究较少.受限于测试技术，目前相关物理试验仅能获取颗粒柱底部的压力分布，难以追踪其内部力链的动态变化规律，而DEM数值模拟能够弥补此不足.结合物理试验和DEM数值模拟从力学角度分析高、低颗粒柱运动机制的形成差异，揭示颗粒受力特性与颗粒柱运动模式的内联机制，对于提高颗粒动力学特性的认识十分必要.

(3) 复杂颗粒形状对颗粒柱坍塌运动与堆积特性的影响机制.颗粒形状是颗粒物质最基本的物理特性之一，现有研究涉及的颗粒形状较为规则简单或者仅对圆球施加抗转动系数，无法全面地反映颗粒形状的影响效应.研究复杂形状颗粒在运动过程中受力与运动特性，建立颗粒形状表征参数与颗粒柱运动堆积参数[L]、[R]、[H]之间的函数关系，这方面研究将有助于从理论上完善对颗粒材料动力学特性的描述.

(4) 滑坡碎屑流在运动堆积过程中对运动路径存在强烈的刮铲效应.目前，绝大部分研究仅关注颗粒柱在刚性床面上的运动堆积过程，涉及颗粒柱在可侵蚀床面上运动堆积机制的研究极少.刮铲效应对碎屑流运动特性和堆积范围的影响机制仍有待深入研究.此外，火山岩碎屑流中经常包含低密度浮石和高密度岩块，而目前涉及颗粒密度对颗粒柱坍塌运动与堆积特性影响机制的研究却十分不足.侵蚀床面和颗粒密度对颗粒柱坍塌运动与堆积特性的影响机制是后续研究中值得关注的问题.

(5) 散体岸坡在动水作用下的失稳运动机制.目前的研究多集中于颗粒柱在静水环境中的运动堆积过程，很少关注水流运动条件对颗粒柱运动与堆积特性的影响机制，此方面需给与重视.研究成果将对深入了解地震等动荷载作用下海底泥石流的运动机制、我国西南多震库区泥石流、库岸滑动等灾害运动机制及其颗粒堆积后形成的堰塞坝特征等方面具有重要的意义.