单元等效节点外力计算.PDF

新加坡南洋理工大学机械与宇航工程学院 许江平 单元等效节点外力计算 对于单元等效节点外力的计算，相信学习有限元的同学都比较熟悉以下的计算公式 e T e fb N bdV , b是单元所受体力 (1) Ve e T e ft N tdS , t是单元所受面力 (2) S e 这里只讨论平面问题，鉴于很多数值计算问题都会涉及面力的计算，如图(1~4)所 示。无论是自己编制有限元程序还是采用ANSYS 、ABAQUS 等软件进行平面问题 的模拟分析，了解其计算步骤都至关重要。因此，下面主要介绍受面力的单元等效 节点外力计算。 图1 Cook’s 斜梁问题 图2 Bathe’s 大应变问题 图3 受剪悬臂梁 图4 受弯悬臂梁 四边形为例：有以下简化 （参考图5 ）： 1 新加坡南洋理工大学机械与宇航工程学院 许江平 图5 平面问题 图6 笛卡尔坐标系下受面力  e e e dV tdA t J dd V 为单元体积   (3) e e dS tds S 为单元表面积，s为边界曲线C长度(l +l l l )  12 23 34 41 其中 x y    J   (4) x y      将公式(3)-(4)代入公式(2) ，可得到 e T e T T T T ft N tdS t N tds t N tds t N tds t N tds (5) S e 12 23 34 45 其中形函数矩阵为 N1 0 N 2