摘要

采用广义多辛数值算法研究不可压饱和多孔弹性Timoshenko梁的流固耦合动力响应特性,构造梁动力响应方程的广义多辛形式,给出其Preissmann Box离散格式及各种广义多辛局部守恒律误差离散格式。数值模拟两端可渗透多孔弹性Timoshenko悬臂梁的动力响应过程,并分析其动力响应特性。发现两相耦合作用系数增大,梁各横截面的孔隙流体压力等效力偶、固相挠度和固相有效应力达到稳态值所需的时间缩短;梁长细比增大,所需时间加长,且挠度稳态值越接近相应经典单相弹性Euler-Bernoulli梁的静挠度值;随时间的推移,梁固相骨架承担所有外荷载,孔隙流体压力等效力偶最终将为零。表征耗散效应的参数取值减小,各种广义多辛数值误差的数量级也减小。