杂交有限元构造的正交假设应力场方法-Core.PDF

第27 卷 第3 期 应 用 力 学 学 报 Vol.27 No.3 2010 年9 月 CHINESE JOURNAL OF APPLIED MECHANICS Sep.2010 文章编号：1000- 4939(2010) 03-0526-06 杂交有限元构造的正交假设应力场方法 张灿辉 王东东 (厦门大学 361005 厦门) 摘要：分别对各向同性和正交各向异性材料的假设应力场进行正交化并形成相应的杂交元，由于 避免了柔度矩阵求逆运算，从而提高了杂交元分析效率。对各向同性材料杂交元直接利用柔度矩 阵特征向量导出了正交假设应力场，其正交性不依赖于材料，因而具有更好的适用性。此外由于 不需要借助本征变形模式进行迭代而避免了复杂的数值计算。对于正交各向异性材料提出了一种 材料矩阵分裂法对假设应力场进行正交化，研究结果表明，所得的正交应力场只与材料两个主方 向弹性模量的比值有关，因而不受横向泊松效应的影响。采用本文方法对2D-4 节点单元和3D-8 节点单元的常用应力场进行正交化，给出了十分简洁的结果。 关键词：杂交应力元；假设应力场；正交化方法；各向同性材料；正交各向异性材料 中图分类号：O242. 21 文献标识码： A 到正交的本征应力模式不仅形式简单而且不包含 1 引 言 任何材料参数，使得其正交性不依赖于材料而具有 更好的适用性，但需要借助本征变形模式进行数值 杂交元是在单元内部独立假设应力场和位移 迭代求解。 场的双变量有限元[1] ，与单变量位移元相比，由于 本文采用特征分析方法对各向同性材料杂交 单元性能得到了改善而受到广泛重视[2-3]，目前已得 元的应力模式进行正交化，不仅给出了不包含任何 到了很大发展[4-6] 。由于泊松效应假设应力场中同类 材料参数的简单结果，而且避免了本征函数法的数 型应力模式之间一般是耦合的，因此在形成单元过 值迭代求解；同时，对正交各向异性材料杂交元提 程中需要对柔度矩阵求逆，对于应力参数比较多 出一种应力模式正交化的材料矩阵分裂法，所得正 （例如三维单元）的假设应力场需要较大工作量， 交应力模式虽然与材料主刚度比值有关，但不受横 影响了杂交元分析的效率。为了提高杂交元分析效 向泊松效应的影响，因此具有更好的适用性。 率，众多学者已经提出了不同的应力模式正交化方 法[7-9] 。文献[7]采用Schmidt 方法对应力模式进行正 2 杂交元假设应力场的正交化 交化，由于所定义内积不包含材料参数，使得计算 方法 过程和结果十分复杂。文献[8-10]也采用 Schmidt 方法进行应力模式正交化，因为所定义内积包含材 2.1 杂交元等价定理 料参数所以得到了形式较为简单的正交应力模式， 在多场变量的有限元分析中，如果分别独立假 但由于结果包含横向泊松效应的影响，所以其正交 设单元应力场和位移场，可以得到杂交应力元刚度 性依赖于不同材料。文献[11]采用本征函数法所得 基金项目：国家自然科学基金 ；福建省科技项目(2007F3096) ；厦门市科技项目(3502 来稿日期：2009-03-03 修回日期：2009-12-22 第一作者简介：张灿辉，男，1967 年生，博士，厦门大学土木工程系，副教授；研究方向——工程结构力学的数值模拟和仿真分析。 E-mail ：chzha