有限元上机作业.doc

1、有限元上机作业目 录1分布载荷作用下的悬臂梁应力计算11.1问题描述11.2计算模型11.2.1有限元模型及网格划分11.2.2边界条件11.2.3材料模型21.3计算结果21.4结果检验21.5小结32悬臂式连接环的应力与变形分析42.1问题描述42.2计算模型42.2.1有限元模型及网格划分42.2.2边界条件52.2.3材料模型52.3计算结果52.4结果检验62.5小结73具有中心孔的薄壁圆筒受均匀拉伸分析73.1问题描述73.2计算模型83.2.1有限元模型及网格划分83.2.2边界条件83.2.3材料模型83.3计算结果93.4结果检验103.5小结114刚架与弹簧混合建模与分析1

2、24.1问题描述124.2计算模型124.2.1有限元模型及网格划分124.2.2边界条件134.2.3材料模型134.3计算结果134.4结果检验144.5小结165过盈配合圆环受力分析165.1问题描述165.2计算模型175.2.1有限元模型及网格划分175.2.2边界条件175.2.3材料模型175.3计算结果185.4结果检验185.5小结196带中心圆孔板的热应力分析196.1问题描述196.2计算模型206.2.1有限元模型及网格划分206.2.2边界条件206.2.3材料模型216.3计算结果216.4小结227壳梁组合结构建模计算237.1问题描述237.2计算模型237.2

3、1有限元模型及网格划分237.2.2边界条件247.2.3材料模型247.3计算结果257.4小结268复杂3D实体建模及受力分析268.1问题描述268.2计算模型278.2.1有限元模型及网格划分278.2.2边界条件278.2.3材料模型278.3计算结果288.4小结289含裂缝结构体的应力强度因子299.1问题描述299.2计算模型299.2.1有限元模型及网格划分299.2.2边界条件309.2.3材料模型309.3计算结果309.4小结3110受热载荷的正方形烟囱建模与温度场求解3110.1问题描述3110.2计算模型3210.2.1有限元模型及网格划分3210.2.2边界条件

4、3210.2.3材料模型3210.3计算结果3311旋转车轮的建模与应力计算3311.1问题描述3311.2计算模型3311.2.1有限元模型及网格划分3311.2.2边界条件3411.2.3材料模型3411.3计算结果3412U型支架的模态分析3512.1问题描述3512.2计算模型3512.2.1有限元模型及网格划分3512.2.2边界条件3612.2.3材料模型3612.3计算结果3613异形截面梁的几何特性和扭转切应力分布3713.1问题描述3713.2计算模型3813.2.1有限元模型及网格划分3813.2.2边界条件3913.2.3材料模型3913.3计算结果39IV1 分布载荷作

5、用下的悬臂梁应力计算1.1 问题描述分析模型如图1-1 所示, 梁的横截面为矩形（长x宽x高 = 10x1x2 m2）， 受到沿长度方向线性分布载荷作用，q=(10000-1000x)N/m。梁的弹性模量200GPa, 泊松比0.3。图 1-1 悬臂梁的的计算分析模型图1.2 计算模型1.2.1 有限元模型及网格划分问题类型：平面应力；单元类型：Solid Quad 4node 182；网格划分结果：单元数500节点数561。图 1-2 计算模型及网格图1.2.2 边界条件梁左边所有自由度为零（ALL DOF=0）；沿长度方向（x方向）上部作用载荷q=(10000-1000x)N/m。1.2.

6、3 材料模型梁各向同性材料，弹性模量EX=200GPa，泊松比PRXY=0.3。1.3 计算结果图 1-3 总位移云图图 1-4 Von Mises应力云图1.4 结果检验悬臂梁受线性分布荷载作用，出现以中性轴为对称的X方向位移，中性轴上部受拉，X向位移为正；中性轴下部受压，X向位移为负。沿载荷作用方向，离梁左端越远Y负方向位移越大，基本满足悬臂梁受载荷时的应力变形状况。图 1-5 X方向位移云图图 1-6 Y方向位移云图1.5 小结梁受沿Y方向随X线性变化的荷载，变形为向下弯曲，离左端越远变形越大，最大位移量为0.268x10-4m。梁的最大Von Mises应力出现在梁左端上部，为2448

7、86Pa，是结构最弱部位。2 悬臂式连接环的应力与变形分析2.1 问题描述分析模型及尺寸如图2-1 所示, 钢质悬臂式连接环左端焊接在墙壁，另一端包含一个圆孔。工作时连接栓将施加向下50MPa均布压力作用于圆孔内下壁。图 2-1 悬臂式连接环受分布载荷作用计算分析模型图2.2 计算模型2.2.1 有限元模型及网格划分问题类型：平面应力；单元类型：Solid Brick 8Node 45；网格划分结果：单元数342，节点数639。图 2-2 计算模型及网格图2.2.2 边界条件悬臂连接环左面所有自由度为零（ALL DOF=0）。圆孔内下壁受50N/mm2均布压力作用。2.2.3 材料模型各向同性

8、钢制材料，弹性模量EX=200x103N/mm2，泊松比PRXY=0.3。2.3 计算结果图 2-3 总位移云图图 2-4 Von Mises应力云图2.4 结果检验通过加密网格重新计算模型，结果如图2-57。总位移及Von Mises应力均与原模型的结果分布相同，且数值相差不大，可认为此计算结果可信。图 2-5 加密网格图图 2-6 网格加密总位移云图图 2-7 网格加密Von Mises应力云图2.5 小结悬臂连接环整体受力变形类似于悬臂梁，变形为向下弯曲，离左端越远变形越大，最大位移量为0.41541mm。最大Von Mises应力出现在梁左端上部，为371.039MPa，是结构最弱部位

9、3 具有中心孔的薄壁圆筒受均匀拉伸分析3.1 问题描述分析模型如图3-1 所示,薄壁圆筒受均匀拉伸 ：100MPa，圆筒内半径：100mm, 圆筒外半径：110mm, 圆筒长度：500mm，中心孔半径：20mm。图3-1 具有中心孔的薄壁圆筒受均匀拉伸的计算分析模型3.2 计算模型3.2.1 有限元模型及网格划分问题类型：空间问题；单元类型：Solid Brick 8node 185；网格划分结果：单元数15552，节点数5165。图3-2 计算模型及网格图3.2.2 边界条件薄壁圆筒左面所有自由度为零（ALL DOF=0），右面受100N/mm2均布拉力力作用。3.2.3 材料模型各向同性

10、材料，弹性模量EX=200x103N/mm2，泊松比PRXY=0.3。3.3 计算结果图3-3 总位移云图图3-4 Von Mises应力云图图3-5 Von Mises应力剖切图3.4 结果检验薄壁圆筒整体受力变形类似拉杆，变形为轴向拉伸，若无圆孔薄壁上拉应力应相同，均为100MPa。由于在圆孔存在，在其附近出现以圆孔为对称中心形成X共轭拉应力集中区。在远离圆孔部位（X=450，50处），绘两圆形路径Pcir和Pcir1（r=105），其Z方向应力在沿圆周的路径上基本在100MPa附近震荡，可认为问题求解基本正确。图 3-6 圆形路径图图3-7 PCIR路径图图3-8 PCIR1路径图3.5

11、 小结薄壁圆筒整体拉伸变形轴线拉伸径向收缩，最大位移量为0.261077mm。在圆孔附近出现应力集中，最大Von Mises应力为圆孔边缘垂直于拉伸方向的对称点上，为294.023MPa，是结构最弱部位。4 刚架与弹簧混合建模与分析4.1 问题描述分析模型如图4-1 所示, 一个1000N载荷垂直施加在一个弹射器上。弹射器建造在一个外径40mm，壁厚10mm钢管上。弹性模量200GPa。弹簧刚度5 N/mm。图4-1 弹射器受力的计算分析模型4.2 计算模型4.2.1 有限元模型及网格划分问题类型：空间问题；单元类型：弹射器基架-PIPE Elast straight 16，弹射棒-Revol

12、ute Joint 7，弹簧-Spring-Damper14；网格划分结果：单元数17，节点数13。图 4-2 计算模型图4.2.2 边界条件弹射器基架接地点所有自由度为零（ALL DOF=0），弹射棒上受1000N集中力作用。4.2.3 材料模型弹射器基架为各向同性材料，弹性模量EX=200KN/mm2，泊松比PRXY=0.33，弹簧单元刚度5 N/mm。4.3 计算结果图 4-3 结构变形图图 4-4 Keypoint7垂直位移图 4-5 Von Mises应力云图4.4 结果检验弹射器受结构对称，作用外力在对称面内，整体变形也对称。检验Keypoint 5、6，其各项位移满足对称要求，同

13、时表现为向内向下的位移，与弹射器受力的实际情况吻合。图 4-6 X方向位移云图图 4-7 Y方向位移云图图 4-8 Z方向位移云图图 4-9 Keypoints 5、6 XYZ方向位移4.5 小结弹射器在力的作用点发生最大变形，最大位移量为320.795152mm。弹射器支架上Von Mises应力大小相同（结构设计合理的表现）， Von Mises应力值为385.288MPa。5 过盈配合圆环受力分析5.1 问题描述分析模型如图5-1 所示, 过盈配合圆环外径200mm, 内径100mm, 厚度50mm，过盈配合通过给环内沿施加径向位移0.5mm表示。在柱坐标下，显示径向应力和环向（周向）应

14、力。材料的弹性模量200GPa, 泊松比0.3。圆环由于对称性，离散（整个模型4分之一）为90度基本扇区。图5-1 过盈配合圆环应力的计算分析模型5.2 计算模型5.2.1 有限元模型及网格划分问题类型：空间问题；单元类型：Solid Brick 8node 185；网格划分结果：单元数216，节点数364。图 5-2 计算模型及网格图5.2.2 边界条件四分之一圆环在两截段面上为对称约束，两扇形面Z方向自由度为零（UZ=0），圆环内环以径向位移0.5mm表示过盈配合。5.2.3 材料模型各向同性材料，弹性模量EX=200KN/mm2，泊松比PRXY=0.3。5.3 计算结果图 5-3 总位移

15、云图图 5-4 Von Mises应力云图5.4 结果检验在内圆外部圆周上建立圆形路径MyP，绘制其X/Y方向应力值图，如图5-5所示。MYPX为X方向应力值，MYPY为Y方向应力值，可计算出其径向应力为（454.9462+837.962）=953.495MPa。根据实际圆环在内环外沿上计算径向应力为200GPa(0.5mm/100mm)=1000MPa。两者相差很小，认为ANSYS求解可靠。图 5-5 MyP上X/Y方向应力值图5.5 小结过盈配合圆环，最大位移位于内环边上，最大位移为0.5mm；沿径向总位移逐渐减小，最小位移在外圆周上，位移量0.318389mm。Von Mises应力分布

16、规律与位移分布相同，最大应力位于内环边上，为1128MPa；最小应力在外圆周上，为322.176MPa。6 带中心圆孔板的热应力分析6.1 问题描述分析模型如图6-1 所示，一个带中心圆孔板，孔直径为0.05m，厚度为0.01m。板由两边钢板及中间复合板粘接而成。三个方板尺寸均为：0.1m0.1m0.01m。热边界条件：圆孔面：100，上、下面为对流换热边界，暴露在空气中，空气的温度为10，换热系数：h=20W/(m2K)。先计算带中心圆孔板内的稳态温度场。在两端固定约束条件下，再进行热结构耦合分析求解。图6-1 过盈配合圆环应力的计算分析模型6.2 计算模型6.2.1 有限元模型及网格划分问

17、题类型：热应力问题（热传导与结构应力耦合问题）；单元类型：Solid Brick 8Node 70；网格划分结果：单元数227，节点数546。 图 6-2 计算模型及网格图6.2.2 边界条件热模型：板上下表面施加对流荷载，对流边界换热系数h=20W/(m2K)，外界温度10，左、右表面视为绝热。孔内弧面温度100。结构模型：板左右面所有自由度为零（ALL DOF=0），板内分布由热传导及对流产生的热应力。6.2.3 材料模型材料导热系数W/(m)比热J/(m)弹性模量GPa泊松比热膨胀系数1密度kg/m3钢板704482000.31.3E-67800复合材料板1.4654400.180.6E

18、616006.3 计算结果图 6-3 温度分布云图图 6-4 总位移云图图 6-5 第一主应力云图（中部复合材料板）图 6-6 Von Mises应力云图（两边钢板）6.4 小结材料热变形，为热传导与结构应力耦合问题，先求解热传导问题计算材料温度分布，再将求得温度代入结构应力问题求解。复合板上温度高，尽管热膨胀系数小，变形却较大，最大位移在孔洞附近及复合板和钢板焊接四角上。结构最危险部位为板边四角，最大Von Mises应力为25.2MPa。7 壳梁组合结构建模计算7.1 问题描述分析模型如图7-1 所示, 确定建筑物框架由于q=0.5 N/mm2的“雪荷载”引起的变形和应力。梁的杨氏弹性模

19、量E=200GPa，泊松比v=0.28，壳的杨氏弹性模量E=80GPa，泊松比v=0.32。结构的长度、宽度方向，梁为等边角钢，长分别为3 m、2 m。结构的高度方向，梁为工字钢，高2 m。壳-梁为空间铰链连接(可以绕Y轴转动)，梁与地面为固定连接。壳厚度5mm，梁横截面积尺寸见下表：表1 等边角钢梁的横截面积尺寸（单位：mm）边宽B板厚t5010表2工字钢梁的横截面积尺寸（单位：mm）宽度B高度H腹板厚t1翼缘厚t21502501515200010001000q=0.5 N/mm2图7-1 建筑物框架受“雪荷载”的计算分析模型7.2 计算模型7.2.1 有限元模型及网格划分问题类型：空间问题

20、单元类型：梁Beam 188，壳Shell 4 node 63；网格划分结果：单元数60，节点数51。图 7-2 计算模型及网格图7.2.2 边界条件建筑物四接地点所有自由度为零（ALL DOF=0），顶部受q=0.5 N/mm2的均布载荷作用。7.2.3 材料模型各向同性材料，梁弹性模量E=200GPa，泊松比v=0.28；壳弹性模量E=80GPa，泊松比v=0.32。7.3 计算结果图 7-3 结构总位移图图 7-4 Von Mises应力云图图 7-4 节点第一主应力云图7.4 小结结构在雪载荷作用下，梁基本无变形，上部板壳变形由两边向中间逐渐增大，中心最大位移为1563mm(太假了不

21、)。8 复杂3D实体建模及受力分析8.1 问题描述分析模型如图8-1 所示，轴承座材料弹性模量为30x106N/cm2（300GPa），泊松比为0，尺寸（老师，这个还是省省吧），在沉孔上受到1000Mpa推力和5000Mpa下向作用力（按单位来说觉得应该是10Mpa推力和50Mpa）。图8-1 轴承座应力的计算分析模型8.2 计算模型8.2.1 有限元模型及网格划分问题类型：空间问题；单元类型：Solid Tet 10Node 92；网格划分结果：单元数5602，节点数9892。图 8-2 计算模型及网格图8.2.2 边界条件轴承座安装孔受对称约束，基座底部竖直方向位移约束为零（UY=0），

22、在沉孔上受到1000Mpa推力和5000Mpa下向作用力。8.2.3 材料模型各向同性材料，弹性模量EX=30x106N/cm2，泊松比PRXY=0。8.3 计算结果图 8-3 总位移云图图 8-4 第一主应力云图附文件：第一主应力动画（3D Entity）8.4 小结轴承座在沉孔上受到推力和下向作用力，变形主要发生在轴承圈和腹板部位。从轴承座基座往上唯一逐渐增大，最大位移发生在轴承座顶部，为0.01032cm。腹板与轴承圈相接处为结构最脆弱部位，存在最大第一主应力，为106.84MPa。9 含裂缝结构体的应力强度因子9.1 问题描述试样为长度0.3m，宽度0.2m的薄板，中间有一长0.04m

23、的裂纹。（和题目里的不一样啊，但提供的ANSYS操作尺寸是这个吧。）材料的弹性模量为210GPa，泊松比为0.3，荷载为均布拉应力140MPa。如图所示。利用对称条件,在ANSYS中取计算模型的1/4建模。图 9-1 裂隙薄板应力的计算分析模型9.2 计算模型9.2.1 有限元模型及网格划分问题类型：平面应力；单元类型：Solid Quad 8node 82；网格划分结果：单元数2228，节点数6829。图 9-2 计算模型及网格图9.2.2 边界条件薄板左边及除裂纹外的下边为对称边界，薄板上部受140MPa均布力。9.2.3 材料模型各向同性材料，弹性模量EX=210GPa，泊松比PRXY=

24、0.3。9.3 计算结果图 9-3 总位移云图图 9-4 Von Mises应力云图图 9-5 Von Mises应力局部放大云图9.4 小结带裂隙薄板受拉，在裂隙尖角处出现应力集中。薄板大部分位置Von Mises应力强度都较低，但在裂隙尖端很小范围内出现1885MPa的最大应力。10 受热载荷的正方形烟囱建模与温度场求解10.1 问题描述假定烟囱壁由两层材料构成，外层材料的导热系数为k=0.1W/(mK)，内层材料为混凝土，导热系数为k=1.4W/(mK)。外部表面的截面尺寸为60cm60cm，内部表面的截面尺寸为30cm30cm，烟囱通道的尺寸为20cm20cm。假定烟囱内表面的温度为1

25、00，烟囱外表面暴露在空气中，空气的温度为30，换热系数为h=20W/(m2K)。图 10-1 双层烟囱的热分析模型10.2 计算模型10.2.1 有限元模型及网格划分问题类型：平面热力学问题；单元类型：Solid Quad 4node 55；网格划分结果：单元数804，节点数865。图 10-2 计算模型及网格图10.2.2 边界条件模型左边及下边为对称边界，内表面温度100，外表面温度30，对流换热系数20W/(m2K)。10.2.3 材料模型内外层材料均为各向同性材料，外层材料的导热系数为k=0.1W/(mK)，内层材料导热系数为k=1.4W/(mK)。10.3 计算结果图 10-3 温

26、度场云图11 旋转车轮的建模与应力计算11.1 问题描述车轮为沿轴向具有循环对称的特性，基本扇区为45度（8分之一模型）。材料密度=7.8E-6（kg/mm3），材料的弹性模量200GPa，泊松比0.3。为给车轮内沿施加径向位移0.5mm，给车轮施加转速n=100r/s的旋转载荷(离心力)。图 11-1 旋转车轮应力分析模型11.2 计算模型11.2.1 有限元模型及网格划分问题类型：空间问题；单元类型：Tel 4node 285；网格划分结果：单元数15306，节点数3555。图 11-2 计算模型及网格图11.2.2 边界条件车轮内外圆面Y方向位移为零，旋转轮体的前、后面为对称面，轮体有绕

27、Y轴环向转速100r/s。11.2.3 材料模型各向同性材料，弹性模量200GPa，泊松比0.3, 密度7800kg/m3。11.3 计算结果图 11-3 Von Mises应力云图附文件：Von Mises应力动画（Wheel）12 U型支架的模态分析12.1 问题描述使用分块兰索斯法提取U型支架的前10阶模态。U型支架形状尺寸如下图所示。图 12-1 U型支架模态分析模型12.2 计算模型12.2.1 有限元模型及网格划分问题类型：空间问题；单元类型：Solid 20node 95；网格划分结果：单元数2821，节点数5159。图 12-2 计算模型及网格图12.2.2 边界条件支架左面所

28、有自由度为零，振动频率范围：0-10000000Hz。12.2.3 材料模型各向同性材料，弹性模量200GPa，泊松比0.3，密度7800kg/m3。12.3 计算结果图 12-3 第一阶模态总位移云图图 12-4 第一阶模态Von Mises应力云图图 12-5 第十阶模态总位移云图图 12-6 第十阶模态Von Mises应力云图附文件：第一阶模态Von Mises应力动画（Usupport）13 异形截面梁的几何特性和扭转切应力分布13.1 问题描述计算图所示异形截面对形心轴的惯性矩。201501001202020单位：mm图 13-1 异形截面图13.2 计算模型13.2.1 有限元模型及网格划分定义异形截面尺寸参数：单元类型：Beam 2 node 188；网格划分结果：单元数10，节点数21。图 13-2 计算模型及网格图13.2.2 边界条件异形截面梁的左端所有自由度为零，右端加绕x轴旋转力偶100Nm。（为什么会是MZ呢？）13.2.3 材料模型各向同性材料，弹性模量200GPa，泊松比0.3。13.3 计算结果图 13-3 总位移云图图 13-4 总转角云图图 13-5 XY面内剪应力云图图 13-6 Von Mises应力云图41