中国地质大学ANSYS实验2..docVIP

实验二 梁结构静力有限元分析 一、实验目的： 加深有限元理论关于网格划分概念、划分原则等的理解。 熟悉有限元建模、求解及结果分析步骤和方法。 能利用ANSYS软件对三维梁结构进行静力有限元分析。 二、实验内容： 梁的截面形状为梯形截面，各个截面尺寸相同。两端受弯矩沿中性面发生弯曲，如图2-1所示。试利用ANSYS软件对此梯形截面梁进行静力学分析，以获得沿梁AA截面的应力分布情况。 三、实验步骤： 1、分析问题： 由于此问题不是轴对称的，梁上各点位移呈圆弧状，有弯曲半径和弯曲中心，所以采用三维实体单元要比采用轴对称单元好一些。其几何形状可以通过柱坐标建立。 合理简化模型。 由于梁弯曲部分的应力不随θ变化，所以可以适当简化模型，取图2-2所示的切片。AB和CD边夹角为5°。 由于不知道切片两侧截面上轴向应力的分布情况，所以只能将弯矩M直接作用在简化模型上。在定义位移约束时仍认为切片两侧保持平面，切片两端只受纯弯矩载荷，即切片端面不受外力载荷。通过有限元分析可以得到受弯矩切片端面处的应力分布情况。因应力与所受弯矩呈线性关系，所以截面上的应力与切片两端面所受弯矩Mp紧密相关。当z值不变时，梁的截面上点A、B、C和D对称分布，所以，分析梁截面时只需取截面的一半。 (2)描述模型的边界条件。 任意节点处沿u（径向）、v（环方向）、w（轴向）的约束情况如表2-1所示。 表2-1 约束条件 1＃面（Face 1） 2＃面（Face 2） U=0（节点A） 无 V=0（所有节点） V=0.0001(rc -r)（所有节点） W=0（沿AB边） W=0（沿CD边） 切片上所有节点均被约束。A节点处，u=0可阻止切片沿r方向做刚体运动；1＃面上所有节点v=0可防止1＃面做圆周运动，对于ABCD由w=0保证切片模型的对称性；2＃面上BC保证2＃面绕r= rc面转动时，2＃面保持平面。比例系数0.0001，这是随意取的，没有特别含义。开始时，不知道rc的确切值，由于rc对应的是纯弯矩，所以A节点处的反作用力Ra为零。假设开始时，rc =60mm或rc =70mm，则两个rc 值对应的Ra分别为2001N和357N。根据线性推断，当Ra=0时有rc =72.2mm。所以，在分析过程中，取rc =72.2mm（为了分析过程简洁，所以在这里给出rc值，实际问题分析中，读者只能自己确定rc 值）。 2、建立有限元模型： 建立工作文件夹。 在运行ANSYS之前，在桌面上的姓名全拼+班号文件夹下，建立姓名全拼3文件夹，在随后的分析过程中所生成的所有文件都将保存在这个文件夹中。 在启动ANSYS时，系统将工作目录指向“姓名全拼3”文件夹。启动后，进入Referenc设置过滤为必要的菜单，操作如下： GUI: Main Menu Preferences Structural 选择单元。 由于采用柱坐标进行三维实体分析，所以选择的单元与前面不同，具体单元信息参考ANSYS帮助文件关于单元的说明。选择单元的操作如下： GUI: Main MenuPreprocessorElement TypeAdd/Edit/DeleteAddStructural SolidBrick 8-node 45 然后关闭Element Types 对话框和Element Type 菜单。下面定义单元的实常数。 定义实常数 定义单元实常数的操作如下： GUI: Main MenuPreprocessorReal ConstantsAdd/Edit/DeleteAdd 由于分析不需要定义实常数，因此可忽略提示，关闭Real Constants菜单。实常数定义完毕后，保存数据库。操作如下： GUI: ToolbarSAVE_DB (4) 定义材料属性 定义弹性模量和泊松比，操作如下： GUI: Main MenuPreprocessorMaterial PropsMaterial Models Structurallinear ElasticIsotropic 在弹出的对话框中输入材料参数： 杨氏模量： 200e9 泊松比： 0.3 (5)定义几何参数 根据切片模型，首先定义切片顶点的8个关键点，然后通过关键点生成切片实体模型。在柱坐标系中生产所需关键点。由于4个关键点是模型图上的A、B、C、D，另外4个是有同样的r和θ但没有显示出来的z轴方向上的与前4个关键点对应的关键点。因此，需要通过模型几何参数创建。 通过参数定义几何实体的操作如下： GUI: