非饱和土渗流-变形耦合的数值分析

基于一维非饱和土的渗流-变形控制方程,采用Flex PDE(Partial differential equation)软件对该耦合方程组进行求解分析。该方法突破了解析法对非饱和土导水系数函数的特殊限定,适用于任意的土-水特征曲线表达式;还可考虑到饱和时的渗透系数以及孔隙率是变量。与解析解相比,该数值解表现较高的精度,具有解决非饱和土耦合问题的可行性。计算分析表明,非饱和土渗流-变形耦合作用对暂态孔隙水压力分布产生重要的影响,在降雨入渗过程中需考虑土体渗流-变形耦合的影响。降雨初期,位移随着时间明显增大,地表出现下沉,考虑耦合效应的孔隙水压力慢于非耦合情况,原因是H值为正的。随着降雨持续时间的增大,地表下沉的速度减缓,到最后变形开始稳定。位移的变化快慢与孔隙水压力变化规律相同。地表沉降量还与初始孔隙水压力分布以及H值密切相关。饱和时的渗透系数以及孔隙率对非饱和土降雨入渗以及稳态流的分布产生影响,但对其地表变形产生的影响微弱。