含面内转角自由度的四结点四边形单元.pdfVIP

苇 30卷 第 18期 山 西 建 筑 Vol30NO．18 2004年 9 月 SHA XI ARCHITE(vI、URE Sep 2004 ·I · 文章编号：1009—6825(2004)18—000103 含面 内转角 自由度 的四结点 四边形 单元 薛 照 钧 摘 要 ：利．羽十二结点等参数单元的形 函数，推导含转角 自由度 的四结 点四边形膜元的单元刚度矩阵，该单元能直接与 每 结点具有六个 自由度的空间梁柱单元相联结，能有效地解决工程实际问题： 关键词：四边形单元，转角 自由度，有限元 中图分类号：TU311．4 文献标识码 ：A 引言 可表示为 ： 在平面 内具有两个 自由度，平面外具有三个 自由度的五 自由 “ ∑(No +N +N 0) (1) 度壳单元，实际上是用一般的平面应力单元与弯 曲板单元进行组 ∑(N U +N V +N胡0) (”=i，j，k，1) (2) 合而形成的一种新单元。由于每个结点仅有五个 自由度 (缺少面 内转角自由度)，所以该单元很难与其他每个结点具有六个 自由 度的空间梁柱单元相联结 。在膜元中引入转角 自由度 ，就可以使 上述 问题得以解决 ，而且一方面可以改善单元的性态提高精度 ， 另一方面还可以为平板壳元提供有效的单元。这一做法使传统 平板壳元在相邻单元共面出现奇异的问题得以解决 ，这为平板壳 元的应用带来了极大的方便 ，能非常有效地解决工程实际问题。 对于所构造的任意四边形单元 四条边界上 的位移做如下假 设：平行于所在边方向的位移沿所在边呈线性变化 ，垂直于所在 图 1 四边 形单 兀 边方向的位移沿所在边呈三次多项式变化。这样单元 四条边界 由于该单元位移模式 中的形函数 ， ， 及N ， 上内结点的位移可以用单元 四个角点的位移表示出来 ，然后利用 具有相同结点数的等参数单元的形 函数描述整个单元 的位移场 。 ， N胡 来源于十二结点等参数单元的形函数，故满足收敛准 在该文 中，就是利用十二结点等参数单元 的形函数，将整个单元 则 ；又 由于该单元边界上 的位移可 由单元结点位移唯一确定 ，故 的位移场用含转角 自由度的四个角点的位移表示出来 ，再进行归 该单元是一种完备的协调单元。 纳与化简就可以得到含转角 自由度 的任意四边形膜元 的形函数， 2 含平面转角 自由度的四结点任意四边形膜元形函数 由此推导了含转角 自由度的四结点任意四边形膜元的单元刚度 的推导 矩阵 ．其与一般的弯曲板单元组合后 ，形成了一种适用范围更广 、 收敛更快 、精度更高及抗畸变能力更强的壳单元 。在相同精度的 情况下，使用该单元求解问题可以减少单元数量 ，降低结构总 自 由度数：该单元能直接与每个结点具有六个 自由度的空间梁柱 单元相联结 ，能有效地解决工程实际问题 ，对于一些斜交工程问 题尤为适用：通过算例验证 ，表明提 出的含转角 自由度的四结点 任意四边形单元是一种高精度 的且实用广泛的单元 。 所有公式均用美 国著名符号运算软件 Mathematica30推 导，可以避免一些人为 的运算错误 。用 Mathematica3．0推导出 图 2 十二结点四边形单元 的一种新的单元展示了一种简单而又高效的运算方法。 图2所示任意 四边形单元边界上的结点 1至 12与 图3所示 1 单元位移模式 十二结点等参数母单元 的结点 1至 12一一对应。以图 1所示 四 rK