广义Biot理论u～U与u～p形式算法计算效率比较

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