存在填隙流体颗粒离散元法理论研究

离散元法(DEM)是分析散体力学行为的数值方法，颗粒离散元法适用于单个几何形状可用圆球近似而不产生显著差异的情况。两圆盘问介质的挤压流动问题是分析两圆球间介质挤压流动的基础。当颗粒之间有法向或切向相对运动时，由于存在填隙流体，颗粒之间产生的法向的挤压力和切向阻力，是湿颗粒离散元法的理论基础。因此，研究存在填隙流体时颗粒之间的相互作用力具有重要的理论价值和实际意义。本文主要研究工作如下： 1．研究了存在壁面滑移时两刚性圆盘间刚塑性介质的挤压流动。首先对现有的库仑摩擦条件下的压力规律做了进一步的研究，找到了一种更合理的圆盘上的压力分布规律。同时，在圆盘边缘处滑移速度一定，介质的滑移速度随着半径线性变化的假设下，引入合理的速度场，得到了另一种圆盘上的压力分布规律，并对不同的摩擦条件及用不同方法计算得到的结果进行了对比。 2．以Reynolds润滑理论为基础，采用小参数法研究了存在滑移时两刚性圆盘间二阶流体的挤压流动，导出了流体的速度场和压力分布，进而求出圆盘间挤压力的解析表达式。 3．研究了存在壁面滑移时两圆盘间宾汉流体的挤压流动，得到了计算流体刚性区和挤压粘性力的公式。结果表明，圆盘上的挤压力随着滑移参数的增加而减小，当流体全部处于流动状态时，即为圆盘挤压牛顿流体时的情况；当流体全部处于刚性状态时，即为圆盘挤压理想刚塑性介质时的情况。 4．采用Reynolds润滑理论和小参数法等研究了存在填隙二阶流体时两球法向及切向运动时的运动，位移及压力分布规律，从而求出挤压力或切向阻力、阻力矩的解析解。结果表明，二阶流体的本构关系虽然比牛顿流体复杂得多，但法向挤压力与牛顿流体相比只有微小的偏离，而切向阻力和阻力矩与牛顿流体相同。所以，当两球间存在填隙二阶流体时，颗粒间的相互作用力可采用存在牛顿流体时的理论结果。 5．研究了存在填隙Herschel-Bulkley流体时两球法向运动时的挤压力，并且证明了法向挤压力的解分别可以退回到宾汉流体或幂律流体情形。 6．利用上述理论以及DEM TRUBAL程序，研究了存在填隙流体时颗粒规则排列单向压缩问题，得到了宏观应力应变曲线随流体性质、壁的运动速度等的变化规律；并对随机排列单向压缩问题进行了数值模拟，得到了宏观应力应变平均曲线；并对干、湿两种情况下颗粒的模拟结果进行了对比。 7．对存在填隙流体时离散元法的应用进行了初步探索，对块石土强夯过程和挡土墙问题等进行了简单的离散元模拟，可以看出采用存在填隙流体时球颗粒间的作用规律后，可以对一些问题定性分析。 本文从细观上得到了颗粒间的相互作用规律，为离散元法增加了新的理论模型，增加了新的思路，具有参考意义和潜在的广阔应用前景。