在离散元(DEM)商业软件PFC2D的基础上,通过将描述流体体变–压力非线性关系的Tait状态方程(EOS)引入计算流体动力学(CFD),建立模拟弱可压缩流体的CFD-DEM耦合计算模块。首先,推导CFD-DEM的控制方程:包括流体–颗粒相互作用力方程、流体运动方程、Tait状态方程和颗粒运动方程。接着,通过PFC2D自定义FISH语言和C++语言将离散化的CFD-DEM控制方程嵌入商业软件PFC2D中。最后,通过单颗粒水下自由沉降运动和一维单面排水固结试验的模拟验证该耦合模块的可行性。模拟结果表明单颗粒自由沉降速度和Stokes理论解接近,一维单面排水固结试验中孔压消散和固结度随着无因次时间因数Tv的变化均与太沙基固结理论符合较好。从而使得用离散元法全耦合分析土体不排水条件下的力学特性成为可能。