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

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