1、 问题分析由于此问题不是轴对称的,梁上各点位移呈圆弧状,有弯曲半径和弯曲中心,所以采用三维实体单元要比采用轴对称单元好一些。其几何形状可以通过柱坐标建立。1.2.1 合理简化模型由于梁弯曲部分的应力不随变化,所以可以适当简化模型,取图1-2所示的切片。AB和CD边夹角为5。由于不知道切片两侧截面上轴向应力的分布情况,所以只能将弯矩M直接作用在简化模型上。在定义位移约束时仍认为切片两侧保持平面,切片两端只受纯弯矩载荷,即切片端面不受外力载荷。通过有限元分析可以得到受弯矩切片端面处的应力分布情况。因应力与所受弯矩呈线性关系,所以截面上的应力与切片两端面所受弯矩Mp紧密相关。当z值不变时,梁的截面上
2、点A、B、C和D对称分布,所以,分析梁截面时只需取截面的一半。1.2.1 描述模型的边界条件任意节点处沿u(径向)、v(环方向)、w(轴向)的约束情况如表1-1所示。 表1-1 约束条件 1面(Face 1)2面(Face 2) U=0(节点A)无 V=0(所有节点)V=(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面保持平面。比例系数,这是随意取的,没有特别含义。开始时,不知道rc
3、的确切值,由于rc对应的是纯弯矩,所以A节点处的反作用力Ra为零。假设开始时,rc =60mm或rc =70mm,则两个rc 值对应的Ra分别为2001N和357N。根据线性推断,当Ra=0时有rc =。所以,在分析过程中,取rc =(为了分析过程简洁,所以在这里给出rc值,实际问题分析中,读者只能自己确定rc 值)。2有限元模型的建立选择单元和定义实常数、材料属性由于采用柱坐标进行三维实体分析,所以选择的单元为8节点6面体单元。由于分析不需要定义实常数,因此可以选择默认值。定义弹性模量和泊松比 :杨氏模量: 200e9 泊松比:定义几何参数根据切片模型,首先定义切片顶点的8个关键点,然后通过
4、关键点生成切片实体模型。在柱坐标系中生产所需关键点。由于4个关键点是模型图上的A、B、C、D,另外4个是有同样的r和但没有显示出来的z轴方向上的与前4个关键点对应的关键点。因此,需要通过模型几何参数创建。通过参数定义几何实体的操作如下:R1=44e-3 R2=R1+88e-3 Z1=65e-3 Z2=14e-3 定义关键点由于几何模型将在柱坐标中创建,所以首先要将坐标系转换到柱坐标。注意:当当前坐标系为柱坐标时,输入提示菜单中的X、Y和Z对应柱坐标的r、(单位为度)和Z。关键点坐标参数如下:1关键点 X=R1,Y=90,Z=02关键点 X=R1,Y=95,Z=03关键点 X=R1,Y=95,Z
5、=Z14关键点 X=R1,Y=90,Z=Z15关键点 X=R2,Y=90,Z=06关键点 X=R2,Y=95,Z=07关键点 X=R2,Y=95,Z=Z28关键点 X=R2,Y=90,Z=Z2生成切片模型通过已定义的8个关键点生成实体模型:首先连接底部的关键点,然后连接顶部的关键点。这些操作均需在笛卡儿坐标系中进行。通过连接关键点而成的线为直线,即切片的边为直边。此处需要这些边为直边,而柱坐标系中生成的线却是曲线。依生成关键点的顺序依次选择关键点,即可得到切片实体模型,结果如图所示。 图 划分网格由于划分实体网格时需要根据具体的边、线进行调整,所以需要显示出边和线。另外,为了便于定义约束,需要
6、显示线的序号并关闭背景。首先显示切片的边线,操作如下:由于只有一种单元和材料,所以在划分网格时,单元类型和参数需自动选定。为便于设置网格划分参数,可显示实体边线和关键点序号,如下图所示。 图 2.5.1设定网格划分参数。从图中可知,L7线从5关键点到1关键点,所以第一个分割出现在B点附近,最后一个分割出现在A点附近。因为希望沿径向网格密度递减,因此Spacing Ratio值为。依同样的方法设定L5、L9和L11的Spacing Ratio值。表2-1总结了每条线网格参数 表2-1 网格参数设定直线序号NDIVSPACEL1,L3,L6,L101L2,L4,L8,L125L7,L9,L118L
7、51/2.5.2划分网格。这里采用六面体单元划分模型网格。其结果如图所示。3 施加载荷并求解:此部分主要是定义模型约束,具体定义过程如下。 定义约束回忆1面边界条件: 节点A(1关键点):u=0 1面 :v=0 沿线AB(L7) :w=03.1.1定义A点(1关键点)约束通过节点定义约束,选择A处节点,。在被约束自由度列表中选UY,位移值为0。3.1.2定义1面约束首先选择1面上的节点。在被约束自由度为UY。这时在图形窗口中显示出1面被施加了周向约束。其结果如图所示3.1.3定义AB线约束通过节点定义AB线的约束,在被约束自由度选择UZ,结果如图所示。3.1.4定义2面上的约束首先回忆2面的边
8、界条件: 2面上所有节点:v=(rc-r) 5线:由于v方向边界条件为空间函数,因此需要通过定义函数来定义约束。首先编辑函数,最后在2面上定义函数边界: 1e-4*(X2+Y2)选择2面上所有节点。被约束自由度为UY,设置2面上的节点环向约束。结果如图所示。3.1.5定义CD线上约束选择CD线,被约束自由度取UZ,其结果如图所示。施加载荷并求解对模型进行求解,如下图所示。 4 查看分析结果:查看等效应力首先显示等效应力等值线图,结果如图所示。从右视图上得知,最大等效应力为147MPa,出现在对称线的底部。查看环向应力显示应力分析等值线图,其结果如图所示。查看中性轴中性轴的位置就是值为零的位置,
9、通过查看梁的中性轴可以简单判断分析结果是否合理。设置使应力为负和应力为正的区域以不同的色彩显示,即在0200MPa区域显示红色,在-200MPa0区域显示蓝色,这两种颜色的相交处为中性轴,如图所示。 有限元分析结果显示中性轴是弯曲的,这与材料力学中关于中性轴的假定相矛盾,考虑到经典理论与工程实际的差别,结果可以接受。查看径向应力显示径向应力分析结果,得到径向应力等值线图,从中可发现整个截面上的径向应力均为拉应力,结果如图所示.查看变形后图形显示出变形前后图形,从图中可知,最大位移DMX=。查看变形过程动画查看三维切片在静力下的变形过程动画如图: 从截面的变形动画中可以看出,前面定义的边界条件(
10、节点A处u=0,沿AB边和CD边w=0)是合理的。验证分析结果首先验证约束是否合理,是否满足约束。这部分已经在后处理部分得到验证。下面验证反作用力是否合理。首先列出反作用力,如图所示,由于模型没有直接承受外力,所以平衡方程中合力应该为零。径向力(FX)大小为,接近零。通过更加精确地选择rc可以将径向力变得更小。周向合力FY和轴向合力FZ也非常小,但不为零,具体原因是FX不为零,因此结构平衡方程是一个近似值。5 结论 通过上述分析结果可知,使用ANSYS分析梁结构受到弯矩时梁体应力状态时,结果比较符合分析要求。在应用ansys软件进行有限元分析时,能够根据不规则边界定义相应函数,能够通过动画显示参数的变化过程,掌握静力分析问题基本的验证技巧。参考文献【1】 王新荣,初旭宏.ANSYS有限元基础教程.北京:电子工业出版社,2011.【2】 邢静忠.ANSYS应用实例与分析.北京:科学出版社,2006. 【3】 刘德彷,李丽英,周临震等,基于ANSYS的建模及分析方法研究 .盐城工学院1009-0134(2011)10(上)-0140-04【4】 李晨,金彦,雷文平等,ANSYS中建立模型的有限元方法【J】.水利电力机械2005,27(2):42-44【5】 王风丽,宋继良,谭光宇等,在ANSYS中建立复杂有限元模型【J】.哈尔滨理工大学学报,2003,8(3):22-28
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1