饱和两相介质动力固结问题时域求解的精细时程积分方法

来源：计算力学学报杂志2017年第1期北京时间：

作者：李亮,高超,吴利华

单位：北京工业大学 城市与工程安全减灾教育部重点实验室, 北京 100124,北京工业大学 城市与工程安全减灾教育部重点实验室, 北京 100124,北京工业大学 城市与工程安全减灾教育部重点实验室, 北京 100124

摘要：针对u-p形式的饱和两相介质波动方程，采用精细时程积分方法计算固相位移u，采用向后差分算法求解流体压力p，建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例，对该方法的计算精度进行了校核。开展了该方法相关算法特性的研究，对采用不同数值积分方法计算非齐次波动方程特解项计算精度的差异进行了对比研究，并对采用不同积分点数目的高斯积分法计算特解项条件下计算精度的差异进行了对比研究。研究结果表明，（1）该方法具有良好的计算精度。（2）计算非齐次波动方程特解项的数值积分方法中，梯形积分法的计算精度最差，高斯积分法、辛普生积分法和科茨积分法都具有较好的计算精度。（3）增加高斯积分点数目对于提高计算精度的作用并不显著。

关键词：饱和两相介质,动力固结,u-p波动方程,精细时程积分,向后差分,

基金资助：国家自然科学基金面上项目（51178011）；国家重点基础研究发展计划973计划（2011CB013602）资助项目