基于平面弹性2板弯曲模拟关系的薄板有限元列式.PDFVIP

第 18 卷第 3 期 计 算 力 学 学 报 V o l. 18 N o. 3 2001 年 8 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS A ugust 2001 文章编号:(2001) 基于平面弹性板弯曲模拟关系的薄板有限元列式 ——从膜单元到薄板单元 郑长良, 黄若煜, 姚伟岸, 钟万勰 (大连理工大学 工程力学系, 辽宁 大连 116023) 摘　要: 本文全面讨论了基于平面弹性板弯曲模拟关系的薄板有限单元的理论和方法。由于直接 对弯矩函数进行插值, c1 连续性的要求得以自然避免。薄板单元可以直接在 c0 连续的层面上加以 构造, 无需借用 的中厚板理论, 由之引发的闭锁问题也得以避免。本文系统地阐 R eissner M indlin 明了平面弹性膜单元与薄板弯曲单元的对应关系, 及由平面弹性膜单元向薄板弯曲单元转换的一 整套方法, 为薄板单元的构造提供了一条新的有效的途径。文中给出了对应于平面弹性膜单元 、 、 4、 8 的薄板单元, 我们称之为 板单元。 板元以挠度和转角为自由度, 便于 CST L ST Q Q M PS M PS 实际应用, 和其它板单元相比具有非常高的精度。 关键词: 板弯曲; 膜板模拟; 收敛性; 有限元 中图分类号: O 241. 82　　　　文献标识码:A 题中去。但令人遗憾的是人们都将其注意力集中到 1　引　　言 平衡模式板单元的构造上, 迄今为止, 仅有 F raeij s de V eub ek e 和Sander [8 ]、Sam ue lsson [9 ]、Sander [10 ] 构 c1 连续性的要求是 K irchho ff 板单元构造中的 造出为数十分有限的所谓平衡单元。这些平衡单元 一大障碍。为了克服这一困难, 人们提出了很多方 从满足平衡方程的弯矩场出发, 求出相应的位移场, 法, 诸如混合法、杂交法、拟协调法、广义协调法、精 其构造过程复杂, 其单元节点位移向量中含有挠度 化法等等, 在一定意义下恰当地降低连续性的要求。 平均值这样的待定未知量, 不便于实际应用, 其精度 直接在 c0 连续的层面构造板单元, 以往的作法是借 也只是差强人意。以后这方面的工作便停滞不前, 可 用R eissn erM in d lin 的中厚板理论, 如著名的离散 以说Sou thw ell 方程的特点和应有价值和潜力还没 K irchho ff 板单元。应用 R eissn erM in d lin 的中厚板 有被充分认识和发挥。近来钟万勰和姚伟岸[3 ] 开始 理论, 仅需