非饱和土粒间毛细作用的微观不连续变形分析

非饱和土粒间毛细水作用是其基质吸力的主要来源。本文基于二维非连续变形分析算法(DDA),建立了一种能够模拟不同饱和度下毛细水分布状态及计算土水特征曲线的扩展DDA算法。该算法首先通过反复迭代的方式来计算出不同饱和度下毛细水弯液面的半径,并利用圆心轨迹交汇法确定粒间毛细水弯液面与颗粒的搭接位置,从而确定各个颗粒表面被毛细水浸湿的区域。而后利用Young-Laplace方程计算出各个浸湿区域毛细负压的大小,并将其与表面张力添加到传统DDA算法控制方程中。为验证该算法,参照真实黄土中含量最多的粗颗粒建立一理想黄土微观土骨架模型,模拟了不同饱和度下土体中的毛细水分布状态,并获得了土水特征曲线。将模拟结果与实测的土水特征曲线进行对比分析,结果显示预测值和实测值基本一致,该方法能够从微观上揭示非饱和土粒间吸力作用的本质。