梁单元、梁板单元和实体块单元有限元数值模拟实例

梁单元、梁板单元和实体块单元有限元数值模拟

及材料力学计算公式之间的比较分析

1. 模型描述：

工字钢的截面尺寸见图1，轴向长100mm，一端固定，另一端加1000N的集中力作用。求此条件下，工字钢的最大挠度。

用三种模型来进行计算，再用材料力学的理论公式进行计算，比较它们的计算结果。第一种模型是用一个梁单元来计算，模型见图2（a）中最下面的模型；第二种模型是用梁单元和板单元组合起来进行计算，模型见图2（a）中中间的模型；第三种模型是用实体块单元来进行计算，模型见图2（a）中最上面的模型。

图 1 工字钢的截面尺寸（单位：mm）

(a) (b)

图 2 计算模型示意图

材料力学的理论公式如下：

图 3 图 4

P是集中力，l是工字钢的长度，H是工字钢的高度。

挠曲线方程（见图3）：???Px26EI(3l?x)

应力值计算公式（见图3和图4）：??2. 挠度比较

P(l?x)H2I

图 5 各模型的挠度计算云图

表 1观察点处的挠度值（单位mm）

单元类型 梁单元 梁板单元 实体单元 理论公式 距固定端20mm处 距固定端50mm处 距固定端100mm处 -0.08 -0.09 -0.09 -0.08 -0.47 -0.46 -0.48 -0.46 -1.49 -1.42 -1.52 -1.47