简支梁有限元结构静力分析.doc

第一章简支梁有限元结构静力分析
0 前言
本文利用ANSYS软件中BEAM系列单元建立简支梁有限元模型，对其进行静力分析与
模态分析，来比较建模时不同单元类型的选择和网格划分精细程度不同所带来的不同结果，以便了解和认识ANSYS对于分析结果准确性的影响。
梁单元介绍
梁是工程结构中最为常用的结构形式之一。ANSYS 程序中提供了多种二维和三维的梁单元，分别具有不同的特性，是一类轴向拉压、弯曲、扭转单元，用以模拟各类结构中的平面以及空间的梁构件。常用的梁单元中 BEAM3、BEAM 23 和 BEAM 54 为二维梁单元， BEAM 4、BEAM 24、BEAM344、BEAM188 和 BEAM189 为三维梁单元。下文将简单介绍常用的梁单元 BEAM3、BEAM4、BEAM44、BEAM188 以及 BEAM189。
BEAM3单元：
图 Beam3 单元几何图形
BEAM3 是具有拉伸、压缩和弯曲的单轴 2-D 弹性梁单元。上图给出了单元的几何图形、节点位置及坐标系统。单元由两个节点、横截面面积、横截面惯性矩、截面高度及材料属性定义。初始应变通过Δ/L 给定，Δ为单元长度 L（由 I，J 节点坐标算得）与 0 应变单元长度之差。该单元在每个节点处有三个自由度，可以进行忽略环箍效应的轴对称分析，例如模拟螺栓和槽钢等。在轴对称分析中，单元的面积和惯性矩必须给出 360°范围内的值。剪切变形量SHERAR 是可选的，如给SHERAR 赋值为0则表示忽略剪切变形，当然剪切模量（GXY）只有在考虑剪切变形时才起作用。同时可以运用实常数中的 ADDMAS 命令为单位长度梁单元施加附加质量。
BEAM4单元：
图 Beam4 单元几何图形
BEAM4 是具有拉伸、压缩、扭转和弯曲的单轴 3-D 弹性梁单元。关于本单元的几何模型，节点座标及座标系统详见上图。该单元在每个节点处有六个自由度。单元属性包括应力刚化与大变形。单元方向由两或三个节点确定，实常数有横截面面积，两个方向的惯性矩
（IZZ 和 IYY），梁的高和宽，与单元轴 X 轴的方向角和扭转惯性矩（IXX），如果没有给出 IXX 的值或赋予 0 时，程序自动假设 IXX=IYY+IZZ，IXX 必须为正同时一般情况下小于弯曲惯性矩，因此最好能够给出 IXX 的值。BEAM4 单元也可以定义附加质量。
BEAM4 单元的 X 轴方向为 I 节点到 J 节点，对于两节点情况，当θ= 0°时，Y 轴平行于总体的 X-Y 平面。用户可以使用方向角θ或者第三个节点控制单元的 Y 轴方向。如果两者都定义了，则以第三个节点定义的方向为主。定义梁单元的方向除了能够控制单元截面形式外还能控制单元各个面的位置，从而能够正确施加梁荷载。
BEAM44单元：
图 Beam44 单元几何图形
BEAM44 是 3-D 弹性渐变非对称单轴梁单元，具有拉伸，压缩，旋转和弯曲的能力。从上图可以看出，该单元在每个节点有六个自由度，允许在每个节点处有不同的非对称几何特性，并且允许节点有相对于梁中心的偏移。如果不考虑该项特征的话，可以用 BEAM4 单元进行替换。如果需要考虑材料非线性，应运用 BEAM188 或者 BEAM189 单元。
BEAM44 同 BEAM3 以及 BEAM4 一样，可以在单元设置中考虑剪切变形，具有应力刚化与大变形特性。该单元可以使用 SECTYPE，SECDATA，SECOFFSET，SECWRITE 和
SECREAD 命令定义任意形状的截面。单元实常数包括截面面积，面积矩，截面边界到中心的距离，截面中心的偏移量，截面剪切常数。其中截面惯性矩是关于梁横向主轴 Y、Z 轴的，在节点 I 的扭矩如果没有特殊指定的话，程序自动默认为惯性矩的和（IZZ+IYY），节点 J 的惯性矩和扭矩没有指定的话默认为和 J 节点的扭矩相同。单元扭转刚度随着扭矩的减小而降低。
单元用偏移距离(DX, DY, DZ)定义截面中心与单元定位节点之间的位置关系。当截面中心相对于节点的偏移方向在单元坐标系的正向上时值为正。所有梁单元第二个节点的实常数
（除了 DX, DY, DZ）默认与第一个节点相同。
在土木工程应用中，建立框架模型时，梁柱构件均运用相同单元，由于侧向位移三次插
值，采用 BEAM4 和 BEAM44 更合适。
BEAM188/189单元：
图 Beam188 单元几何图形 图 Beam189 单元几何图形
BEAM188/189 单元的几何形状、节点位置、坐标体系如上图所示，BEAM188 单元的方向由整体坐标系的节