司宪志，宁 宇，石 崇*，黄青富，张一平

(1.河海大学 岩土工程研究所，江苏 南京 210098；2.中国电建集团昆明勘测设计研究院有限公司，云南 昆明 650051)

我国西南地区多为山地丘陵地貌，在这些区域进行机场建造时，为制造出合适的场地条件和净空条件，势必会进行深挖高填工作，产生大量的高填方边坡[1]，同时山地及丘陵地带复杂的地质条件，也给机场建设的安全性带来了挑战，填方地基及高填方体的稳定性问题成为山区机场建设中存在的核心问题[2]，加强对该高填方边坡稳定性的研究，寻求科学、合理、有效的评价及防治措施，具有重要的现实意义。数值模拟由于其成本低、效果明显的特点被广泛地应用于边坡稳定性分析[3]中，在连续数值仿真方法方面：陈正东[4]等基于有限差分软件FLAC3D软件，建立了三维边坡数值计算模型，研究了开挖后的边坡体的位移变形特征，并据此对开挖后的边坡进行了稳定性评价；潘网生[5]等基于FLAC3D拉格朗日快速差分算法模拟煤层开挖，揭示边坡变形过程及其滑动机理；言志信[6]等以含软弱层的锚固岩体边坡为研究对象，利用FLAC3D软件模拟分析了地震作用下含软弱层岩体边坡两锚固界面剪切作用及其演化。但连续数值方法中单元之间需要公用节点才能合理描述荷载力的传递，所能模拟的变形量有限，当变形增大到一定程度时，单元会发生畸变导致刚度矩阵无法求解，导致计算终止。在非连续数值仿真方法方面：李龙起[7]等利用颗粒流软件PFC2D对土质滑坡在持续降雨作用下的变形发展过程进行模拟研究，揭示了该滑坡的发展及破坏模式；……