有限元实体建模和网格划分.docx

PAGE1

2.有限元实体建模和网格划分

2.1建立几何模型

由于的分析对象为房屋，其建造过程中由许多梁和截面组成，如果利用三维实体模型去模拟，往往会得不到比较好的分析结果，因此，需要进行简化，将柱子简化为梁单元，将屋子的标准层以及屋面层简化为面单元。对于屋子模型，其由面和线组成，因此只需建立线和面，建模所选取的单位制为m，在apdl的Modeling模块，创建100键点，通过创建的关键点建，如下图2所示。

图2屋子几何模型

2.2添加单元类型

在ElementType中添加单元类型为BEAM188单元、SHELL181单元、SURF154单元、MASS21单元，其中，BEAM188单元用来模拟房屋的柱子结构，SHELL181单元用来模拟各个楼层的。添加的单元类型如下图3所示。

图3添加单元类型

其中，BEAM188单元的截面形状都为长方形，并设置三种不同的截面尺寸，分别为500*500，500*300，500*250。设置SHELL181的厚度为130mm。

图4不同的BEAM单元截面尺寸

2.3添加材料

在Preprocessor→MaterialProps→MaterialModels→Structural→Linear→Elastic→Isotropic→杨氏模量（EX），波松比（PRXY）的输入以表1的材料进行添加。该房屋采用混凝土建成，设置该材料的弹性模量为30000MPa，泊松比为0.3，密度为2600kg/m3，材料参数设置如图5所示。

表1材料参数

2.4网格划分

几何模型建立完成后，还需要进行网格划分才能进行后续的有限元分析。网格质量会对后续的结构分析有着重要的影响，网格数量不足，则会导求解精度不够，网格数量过大，则会导致求解时间过长，因此需要合理设置网格数量。Preprocessor→Meshing→Meshtool→在弹出的工具框内：在ElementSize选项，设置全局划分网格尺寸为2，点击mesh。网格划分结果如图6所示

图6网格划分结果

2.5施加约束与添加载荷

载荷与约束的施加正确与否对有限元的分析结果有重要的影响，根据房屋的实际受到的约束，其底部与大地进行固定，因此约束房屋的底部节点的全部自由度。约束的添加结果如图7所示。屋面层楼板附加恒载以及重力加速度载荷，载荷的添加结果如图8所示。载荷以及约束采用命令流的方式进行施加。

图7添加约束

图8施加载荷

命令流：

!施加约束

nsel,s,loc,z,-3900

d,all,all,0

!=

/solu

allsel,all

!=

!恒载计算

time,1

!墙体线荷载

esel,s,type,,1

esel,u,cent,z,3900*5

sfbeam,all,1,pres,7

!修改外墙荷载

esel,s,cent,y,0

esel,a,cent,y,10680

esel,r,type,,1

esel,u,cent,z,3900*5

sfbeam,all,1,pres,5.2

!=

!标准层楼板附加恒载

esel,s,type,,3

esel,u,cent,z,3900*5

sfe,all,1,pres,1,-1.5e-3

!=

!屋面层楼板附加恒载

esel,s,type,,3

esel,r,cent,z,3900*5

sfe,all,1,pres,1,-4.0e-3

2.6求解结果分析

载荷与约束添加完成后，提交求解，命令流如下：

allsel,all

acel,,,9800

solve

提取分析结果，ANSYSmainmenu→GeneralPostproc→PlotResults→ContourPlot→NodalSolu→Stress→VonMisesstress→OK。屋子为混凝土材料，其强度应按第四强度理论计算，因此在查看计算结果时一般查看其等效应力（VonMises应力）的高低。最终变形结果如图9所示，应力结果如图10所示。从图中可以看出，房屋的变形为8.26mm，最大变形位置位于屋子顶部，其变形数值较小，满足要求。屋子的最大应力数值为4.2MPa，最大应力数值在屋子顶部，应力数值较小，屋子的强度满足要求。

图9变形结果

图10应力结果