基于四边形网格的下限原理有限元法.pdfVIP

第 31 卷 第 3 期 岩石力学与工程学报 Vol.31 No.3 2012 年3 月 Chinese Journal of Rock Mechanics and Engineering March ，2012 基于四边形网格的下限原理有限元法 李春光，朱宇飞，刘 丰，邓 琴，郑 宏 ( 中国科学院武汉岩土力学研究所 岩土力学与工程国家重点实验室，湖北 武汉 430071) 摘要：为了克服四边形网格下限法不能充分利用线性规划算法的不足，通过在四边形网格积分意义上建立平衡方 程的弱形式，从而得到线性的平衡方程，克服了常规基于插值应力场得到的方程为非线性的缺点。另外，四边形 网格比三角形网格具有更少的自变量和约束方程，使求解效率得到较大提升。算例表明：四边形网格的精度和效 率都优于三角形网格，这与普通有限元是一致的。 关键词：数值分析；下限分析；有限元法；四边形网格；线性规划 中图分类号：O 241 文献标识码：A 文章编号：1000 –6915(2012)03 –0461 –08 LOWER BOUND LIMIT ANALYSIS USING QUADRILATERAL FINITE ELEMENTS LI Chunguang ，ZHU Yufei ，LIU Feng ，DENG Qin ，ZHENG Hong (State Key Laboratory of Geomechanics and Geotechnical Engineering，Institute of Rock and Soil Mechanics ， Chinese Academy of Sciences，Wuhan，Hubei 430071 ，China) Abstract ：Triangular meshes are always used in previous lower bound analysis by finite element method. For quadrilateral meshes，the equilibrium equations would be nonlinear on past form. To overcome this shortcoming ， the weak form of equilibrium equations are established by integral over the whole elements and the equations are linear，so that，the linear programming can be used. On the other hand ，there are less variables and constraint equations for quadrilateral element than triangle element，and so using quadrilateral element is more efficient than using triangle element. Examples show that the results based on quadrilateral elements are more accurate than that of triangle elements，which are consistent with general finite elements. Key words ：numerical analysis ；lower bound analysis ；finite elements ；quadrilateral meshes ；linear programming 丰富而被工程单位广泛采纳，但极限平衡法并不总 1 引 言