构造杂交应力单元的柔度矩阵h对角化方法h=∫vep.pdfVIP

第19卷第4期 　 计 算 力 学 学 报 　 V o l. 19, N o. 4 2002 年 11 月 Ch in ese Journa l of Com puta t iona l M echan ic s N ovem b er 2002 文章编号:(2002) 构造杂交应力单元的柔度矩阵H 对角化方法 张灿辉, 　冯　伟, 　黄　黔 (上海大学 上海市应用数学和力学研究所, 上海, 200072) 摘　要: 证明了杂交元柔度矩阵H 非奇异的充分必要条件是假设应力模式线性无关; 以及等价应力模式形成相 同的杂交元。在此基础上建立了假设应力模式的H ilb er t 子空间, 从而可以利用 Schm idt 方法简单地得到等价的 正交应力模式, 实现了柔度矩阵H 对角化, 使得杂交元形成过程中完全避免了繁杂的矩阵求逆运算, 提高了杂交 元分析的计算效率, 特别在柔度矩阵不容易显式求逆的材料非线性分析中更具有实际意义。 关键词: 杂交应力有限元; H ilb er t 应力子空间; 柔度矩阵H 对角化; 材料非线性分析 中图分类号: O 242 21　　　 文献标识码: A 一步独立假设应力场 1　引　言 M 为了克服板壳单元存在的连续位移场构造困 = = ( 1) P i i ∑ i= 1 难和剪切闭锁现象, 1964 年卞学璜[ 1 ] 利用独立假 其中, i 为假设应力模式, i 为对应应力参数, P 为 设应力场的方法提出杂交应力有限元。为了避免单 假设应力模式组成的应力矩阵。记: 元零能变形模式一般需要在假设应力场中增加适 = T - 1 (2) 当高阶应力模式, 但是过多应力模式使假设应力场 Ke G H G 包含大量应力参数, 需要花费很大工作量对单元柔 其中 度矩阵求逆, 降低了杂交元分析计算效率。为了提 T T 高计算效率, 可以减少假设应力场中应力分量个数 = e , 　 = e