基于非饱和流固耦合理论,在考虑初始孔隙水压力场生成、非饱和区饱和度和渗透系数设置、降雨入渗和出渗边界开发及非饱和土有效应力修正的基础上,对FLAC3D软件中饱和流固耦合模块进行二次开发,建立土体饱和-非饱和流固耦合计算模块,并以Del Monte砂柱排水试验为对象进行数值模拟。结果表明,采用文中方法模拟计算的孔隙水压力、饱和度和排水速度更接近试验实测数据,开发的饱和-非饱和流固耦合计算模块可靠、合理。