分区混合元法分析平面裂纹问题-工程力学-清华大学.PDFVIP

第18卷第6 期 工 程 力 学 Vol.18 No.6 2001 年 12 月 ENGINEERING MECHANICS Dec. 2001 文章编号：1000-4750(2001)06-039-08 分区混合元法分析平面裂纹问题 傅向荣，龙驭球 (清华大学土木系，北京 100084) 摘 要：基于分区混合能量原理的分区混合元法是一种高精度有限元方法，可用于分析含裂 纹、孔洞、切口等缺陷的问题。本文推导了分区混合元法中应力元的刚度矩阵；分析了应力 元的多余零能模式，并证明了整体分析中对多余零能模式的消弭；文中还对分区混合元法中 应力元的应力项数与应力单元尺寸对分析结果的影响进行了系统的讨论。 关键词：分区混合能量原理；分区混合元法；多余零能模式；应力强度因子 中图分类号：O34 文献标识码：A 1 引言 作为一种强有力的工程数值分析工具，有限单元法在其近五十年的发展历程中不断完 善和改进。有限单元法的理论基础—— 变分原理也在实践中得以升华。从势能原理，余能 原理，Hellinger-Reissner 原理、胡一鹫原理到分区变分原理，对变量的类型选择不断拓宽。 [1-4] 从协调元，非协调元、广义协调元到分区混合元 ，对变量的约束更加实用化，合理化。 分区混合能量原理允许在求解体系的不同区域采用不同的变量，这种灵活性在处理含裂 纹，孔洞，切口等缺陷的结构时显得尤为重要。本文拟采用基于分区混合能量原理构造的 分区混合元分析平面裂纹问题，对这一新型单元的基本理论进行分析。直接导出分区混合 元法中应力元的刚度矩阵，并对其收敛性和零能模式进行分析。 2 分区混合元的构造原理 以图1 所示含裂纹板为例。 y 分区混合元以势能( C ) 区的结点位移和余 1 x C2 能( C2 ) 区的应力参数作为基本未知量。在势能 C1 区可采用常规的位移型单元，结点位移{w}为： 图1 含裂纹板 {w} [w1 w2 Lw2n 1 w2 n ]T (1) ——————————————— 收稿日期： 基金项目：国家自然科学基金，清华大学基础研究基金(JC1999002) 作者简介：傅向荣(1972.5)，男，博士生 龙驭球(1926.1)，男，中国工程院院士，教授 40 工 程 力 学 分区混合元法分析平面裂纹问题 41 42 工 程 力 学 分区混合元法分析平面裂纹问题 43 4 分区混合元法的收敛性分析 以图2