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

工 程 力 学 ——————————————— 基金项目：国家自然科学基金面上项目；国家重点基础研究发展计划973计划项目（2011CB013602） 通讯作者：李亮(1975―)，男，山西人，博士，副教授，硕导，主要从事土动力学理论和数值计算方法研究(E-mail: liliang@bjut.edu.cn). 饱和两相介质动力固结问题时域求解的精细时程积分方法 李 亮，高 超，吴利华 (北京工业大学城市与工程安全减灾教育部重点实验室，北京 100124) 摘 要：针对形式的饱和两相介质波动方程，采用精细时程积分方法计算固相位移，采用向后差分算法求解流体压力，建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例，对该方法的计算精度进行了校核。开展了该方法相关的算法特性研究，对采用不同数值积分方法计算非齐次波动方程特解项时计算精度的差异进行了对比研究，并对采用不同积分点数目的高斯积分法计算特解项条件下计算精度的差异进行了对比研究。研究结果表明：（1）该方法具有良好的计算精度；（2）计算非齐次波动方程特解项的数值积分方法中，梯形积分法的计算精度最差，高斯积分法、辛普生积分法和科茨积分法都具有较好的计算精度；（3）增加高斯积分点数目对于提高计算精度的作用并不显著。 关键词：饱和两相介质；动力固结；波动方程；精细时程积分；向后差分 中图分类号：TU354 文献标识码：A Precise time-integration method for the time-domain solution of the dynamic consolidation problem of fluid-saturated porous media LI Liang，GAO Chao，WU Li-hua (The Key Laboratory of Urban Security and Disaster Engineering ( Beijing University of Technology ), Ministry of Education, Beijing 100124, China) Abstract: A precise time-integration method for the time-domain solution of the dynamic consolidation problem of fluid-saturated porous media is proposed on the basis of dynamic formulation. The solid-phase displacement is computed by the precise time-integration algorithm, and the pore fluid pressure is computed by the backward difference algorithm. The proposed method is verified by a standard numerical example. The characteristic of the proposed method is also investigated. The computational accuracy is compared when several numerical integral methods are employed to calculate the special solution term of dynamic formulation. The computational accuracy of Gauss integral method is also compared when different numbers of integral points are employed. It is shown that (1) the proposed precise time-integration method has fine computational accuracy. (2) Gauss integral method, Simpson integral method and Coates integral method all have fine computational accuracy, while trapezoidal integral method has poor computational accuracy when used to compu