基于颗粒流理论,分别建立大尺度一维、二维颗粒流数值模型,研究低频动载应力波的传播规律。在一维颗粒链模型中,分别从震源激发方式、黏结模型、数值弥散、边界条件的角度系统地研究P波在数值模型中的传播规律,为模拟动载应力波在二维颗粒流模型中的传播提供了参考数据。大尺度一维颗粒链模型(6 400 m)研究结果表明:(1)雷克子波不存在高频角点,比正弦波更适合模拟波场震源;(2)当颗粒链模型采用接触黏结时,P波振幅随距离的增大而衰减,但波速基本不变,与理论值的相对误差为0.12%;(3)当颗粒链模型采用平行黏结时,P波波形产生强烈扭曲;(4)数值弥散分析表明,雷克子波峰值频率的增大和震中距的增大均会使弥散程度增大,但峰值频率的增大主要使波的振幅降低,震中距的增大使振幅和波速均降低。当频率f≤60.23 Hz,角频率?与波数k近似成线性关系,弥散现象可忽略;(5)边界条件分析表明,对于刚性边界,反射波与入射波的振幅大小相等,符号相反,边界处质点的位移等于0;对于自由边界,反射波与入射波的振幅完全相等,符号相同,边界处质点的振幅为入射波的2倍;对于吸收边界,反射波振幅比入射波相比大约减少了96.8%,能有效降低反射波的影响。根据以上结论,在大尺度二维六角形排列颗粒流模型(3 000 m?3 000 m)中,采用接触黏结模型,并施加吸收边界条件,分别探讨了P波和S波波阵面特征及波形与波速的变化,并定性地分析了远场条件下的场源效应。结果表明:(1)当在震源处施加集中力时,模型中出现由点震源激发的P波和S波。考虑到颗粒间的切向刚度为零,S波的存在说明颗粒几何排布对应力波在颗粒模型中的传播规律影响很大;(2)当震源频率一定时,震中距越大,P波和S波的振幅及波速均因弥散作用而降低,且S波弥散程度更大;(3)场源效应分析表明,在震源处施加集中力时,P波位移场在力的方向上幅值最大,在垂直于力的方向上无位移;S波位移场在垂直于力的方向上幅值最大,在力的方向上无位移。上述研究得到了应力波在大尺度一维、二维颗粒流模型中传播的规律,与理论分析结果较为吻合,对进一步模拟研究岩石的动态力学行为及其细观机制具有重要的参考价值。