张 鸿，张 榜，丰浩然，吴 灿

(1.南昌工程学院 江西省水利土木工程基础设施安全重点实验室，江西 南昌 330099；2.南昌工程学院 土木与建筑工程学院，江西 南昌 330099)

煤系地层边坡开挖后，风化作用加强，当雨水渗入使得煤系土快速崩解软化，力学强度急速降低，边坡极易发生坍塌滑坡现象。目前，对煤系土的研究较少，尚处于起步阶段，而且现有的研究主要集中在对煤系土不良地质边坡的加固处理方法[1-3]、边坡稳定性分析[4-6]、以及煤系土强度特性等宏观层面上[7-9]，从细观角度对煤系土边坡失稳的细观机理进行分析与解释的文献未有报道。煤系土边坡的失稳破坏运动是一个存在岩土体滑动、平移、转动的复杂过程，具有宏观上的不连续性和单个块体运动的随机性，用传统的宏观连续介质理论已经不能合理地分析散体过程，上述煤系土边坡风化之后的散、动特征都与传统的均匀、连续等假定冲突，将导致理论与实际的偏离。因此，要想揭示煤系土边坡的失稳机理，有必要从微细观角度进行研究[10]。

降雨入渗是边坡发生失稳破坏的主要诱发因素之一，目前考虑流固耦合的计算方法主要有3类，第1类是Euler-Euler方法，该方法是将流体和固体都视为连续介质；第2类是Lagrange-Lagrange方法，该方法是把流体和固体都假设为离散介质；第3类是Euler-Lagrange计算方法，该方法将流体视为连续介质，采用流体动力学方法（computational fluid dynamic，CFD）进行模拟计算，固体视为离散介质，采用离散单元方法（discrete element method，DEM）进行模拟计算。这种DEM-CFD流……