论文部分内容阅读
广义Biot理论常用于描述连续多孔介质在荷载作用下的力学行为，但其数值求解较为困难，不同研究者通常采用u～U和u～p两种有限元列式求解广义Biot理论的控制方程，但对两种解法的求解稳定性和求解效率一直缺乏系统的比较。基于Compaq Visual Fortran6．6环境，开发了二维广义Biot理论u～U与u～p形式有限元数值求解程序。大量试算表明，u～U形式刚度阵病态性较弱，其计算稳定性优于u～p形式。通过实例比较不同网格剖分规模与不同时段步数条件下u～U与u～p形式算法所花费的CFU时间，观察到u～p