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

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

基于精细时程积分的饱和两相介质波动问题时域解法

针对u-p形式的饱和两相介质的波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了基于精细时程积分技术的饱和两相介质波动问题的时域求解方法。针对标准算例,将该方法的计算结果与可视为标准结果的Zienkiewicz隐式算法的计算结果进行比较分析,二者符合较好,表明了该方法具有良好的计算精度。同时,该方法的计算过程为交替迭代求解,避免了在每个时间分析步上求解耦联方程组,因而具有较高的计算效率。该方法具有时域显式计算方法的基本特点,是进行饱和两相介质动力问题计算与分析的一种有效方法。

全井钻柱系统耦合振动多体动力学模型的建立与算例分析

基于绝对节点坐标法建立全井钻柱系统的多体动力学模型,研究系统的耦合振动现象;将大长细比柔性钻柱离散为绝对节点坐标梁单元,讨论梁单元格式,并研究井口、钻头处边界及钻柱与井壁的接触/摩擦模型,给出包含绝对节点坐标梁单元的钻柱系统运动方程;采用向后差分法求解微分-代数方程组,开发多体动力学求解器及相应的前、后处理器。通过直井、定向井算例分析了全井钻柱系统的轴向、扭转、横向耦合振动特性,结果表明本文提出的力学建模和数值分析方法可实时捕捉到钻柱系统的耦合振动现象,能够在钻柱系统动力学研究和工程应用中发挥作用。

空间站柔性展开机构仿真方法研究

空间站展开机构是典型的大型多柔体系统,其展开收拢过程涉及非线性大变形以及接触碰撞问题。首先采用绝对节点坐标法作为柔性体的建模方法,建立了梁单元的有限元单元格式,然后基于拉格朗日方法和Hertz接触理论建立了考虑接触碰撞问题的柔性多体系统动力学方程,采用向后差分方法数值求解柔性多体系统的微分代数方程,最后以自由梁、简支梁和柔性单摆进行数学仿真验证。结果表明,该方法可用作我国空间站大型柔性太阳电池翼的设计分析参考。