毛细滞回内变量模型的高效数值积分方法

摘要

毛细滞回内变量模型是基于热动力学基本原理推导出来的土水特征本构模型,能够有效地描述干湿循环情况下非饱和土中的毛细滞回现象.研究了该模型的形式特点和数值积分方法,就经典的Euler法、4阶Runge-Kutta法和4阶Adams-Bashforth法的积分效果进行了对比.计算结果表明,Euler法精度很低,累计误差较大;4阶Adams-Bashforth法精度最高,且运算效率高于同阶的Runge-Kutta法,适合对该模型进行求解和参数标定.将基于Adams-Bashforth法的子程序集成到有限元程序中,两相渗流模拟结果的精度有了较大提高.