有限元分析与应用大作业

1、有限元分析及应用大作业课程名称: 有限元分析及应用班 级: 姓 名: 试题2：图示薄板左边固定，右边受均布压力P=100Kn/m作用，板厚度为0.3cm；试采用如下方案，对其进行有限元分析，并对结果进行比较。1) 三节点常应变单元；（2个和200个单元） 2) 四节点矩形单元；（1个和50个单元） 3) 八节点等参单元。（1个和20个单元） 图2-1 薄板结构及受力图一、建模由图2-1可知，此薄板长和宽分别为2m和1.5m，厚度仅为0.3cm，本题所研究问题为平面应力问题。经计算，平板右边受均匀载荷P=33.33MPa，而左边被固定，所以要完全约束个方向的自由度，如图2-2所示。取弹性模量E=

2、2.111Pa，泊松比=0.3。P=33.33MPa图2-2 数学模型二、第一问三节点常应变单元（2个和200个单元）三节点单元类型为PLANE42，设置好单元类型后，实常数设置板厚为0.3M。采用2个单元的网格划分后的结果如图2-3,200个单元的网格划分图如图2-6所示。约束的施加方式和载荷分布如图2-2中所示。约束右边线上节点全部自由度。计算得到的位移云图分别如图2-4、7所示，应力云图如图2-5、8所示。图2-3 2个三角形单元的网格划分图图2-4 2个三角形单元的位移云图图2-5 2个三角形单元的应力云图图2-6 200个三角形单元的网格划分图图2-7 200个三角形单元的位移云图图

3、2-8 200个三角形单元的应力云图三、第二问四节点矩形单元的计算四节点单元类型为PLANE42，设置好单元类型后，实常数设置板厚为0.3M。采用1个单元的网格划分后的结果如图2-9,50个单元的网格划分图如图2-12所示。约束的施加方式和载荷分布如图2-2中所示。约束右边线上节点全部自由度。计算得到的位移云图分别如图2-10、11所示，应力云图如图2-13、14所示。图2-9 1个四边形单元的网格划分图图2-10 1个四边形单元的位移云图图2-11 1个四边形单元的应力云图图2-12 50个四边形单元的网格划分图图2-13 50个四边形单元的位移云图图2-14 50个四边形单元的应力云图四、

4、第三问八节点等参单元的计算四节点单元类型为PLANE82，设置好单元类型后，实常数设置板厚为0.3M。采用1个单元的网格划分后的结果如图2-15,50个单元的网格划分图如图2-18所示。约束的施加方式和载荷分布如图2-2中所示。约束右边线上节点全部自由度。计算得到的位移云图分别如图2-16、17所示，应力云图如图2-19、20所示。图2-15 1个八节点等参单元的网格划分图图2-16 1个八节点等参单元的位移云图图2-17 1个八节点等参单元的应力云图图2-18 20个八节点等参单元的网格划分图图2-19 20个八节点等参单元的位移云图图2-20 20个八节点等参单元的应力云图四、计算结果对比

5、单元数最大位移(m)最大应力(MPa)最小应力(MPa)三节点常应变单元20.247E-333.628.72000.246E-336.828.8四节点矩形单元10.243E-336.129.6500.244E-337.628.2八节点等参单元10.242E-335.132.7200.244E-343.727.4表2-1 计算结果对比对于三节点常应变单元,两种网格划分情况的最大位移,最大应力和最小应力在数值上都差别不是很大,但是应力分布却存在较大的差别。2单元的最大位移位于薄板的右边角附近。而200个单元的最大位移分布在整个右边上。最大应力两单元的位于左边的上角点处，而200单元的位于薄板左边的

6、上下角点处。最小应力2单元的位于左下角处，而200单元的位于左边的中部。2单元网格划分的网格结构也不具有对称性，模型存在较大的误差。对于四节点矩形单元,两种网格划分情况的最大位移,最大应力和最小应力在数值上差别不是很大。两种划分方式下的最大位移都位于薄板的右边。最大应力1单元的位于右边上，而50单元的位于薄板左边的上下角点处。但是最小应力1单元的位于左边上，而50单元的位于左边的中部。另外，应变和应力对于1个单元的从右至左呈均匀分布，而50个单元的却不是均匀分布的。对于八节点矩形单元,两种网格划分情况的最大位移,最小应力在数值上差别不是很大。最大应力却有着些许的不同，这大概是对于高阶的单元更能

7、准确的模拟上下左边角处的应力集中现象。两种划分方式下的最大位移都位于薄板的右边。最大应力1单元的位于右边上，而20单元的位于薄板左边的上下角点处。但是最小应力1单元的位于左边上，而20单元的位于左边的中部。另外，应变和应力对于1个单元的从右至左呈均匀分布，而20个单元的却不是均匀分布的。五、计算命令采用八节点矩形单元的划分20个网格命令流如下：/FILNAME,shiti9/TITLE,ANALYSIS OF PLATE STRESS!前处理/PREP7ET,l,plane82MP,EX,1,2.1E11MP,PRXY,1,0.3RECTNG,0,1.5,0,2/PNUM,AREA,1 APL

8、OT/TITLE,GEOMETRIC MODELAPLOT/PNUM,KP,1/PNUM,line,1LPLOTlesize,1,4lesize,2,5lesize,3,4lesize,4,5!mshape,1amap,1,1,2,3,4/TITLE,ELEMENTS IN MODELEPLOTFINISH!求解/SOLUANTYPE,STATIC/PNUM,LINE,1LPLOTLSEL,S,LINE, ,4NSLL,S,1D,ALL,allLSEL,S,LINE,2NSLL,S,1SF,ALL,PRES,-33300000ALLSELOUTPR,BASIC,ALLSOLVEFINISH!后

9、处理/POST1PLDISP,2PLNSOL,U,SUMPLNSOL,S,EQVFINISH!/EXIT试题 3：图示图示为一带圆孔或方孔的单位厚度（0.1M）的正方形平板，在x方向作用均布压力0.25Mpa，试用三节点常应变单元和六节点三角形单 元对平板进行有限元分析，并分别就圆孔和方孔结构对以下几种计算方案的计算结果进行比较：1）分别采用相同单元数目的三节点常应变单元和六节 点三角形单元计算；2）分别采用不同数量的三节点常应变单元计算；3）比较圆孔和方孔孔边应力水平；在Y轴上，圆孔边的应力的精确解为：x=-0.75MPa在X轴上，圆孔边的应力的精确解为：y=-0.75MPa图3-1 平板结

10、构及受力图一、建模由图3-1可知，本题所研究问题为平面应力问题，又此平板结构关于图示中X、Y轴对称，可以利用此对称性，取截面的四分之一进行分析计算。此时要约束掉下水平边的Y方向自由度和竖直左边的X方向自由度。载荷为均布压力，平均分布在右边上，大小为0.25MPa。弹性模量E=2.111Pa，泊松比=0.3.数学模型简图如图3-2所示。 图3-2 数学模型二、单元数相同类型不同的圆孔板计算结果分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算：本文采用的三节点单元类型为PLANE42，该单元是四节点类型的，但是可以退化为三节点单元。六节点单元类型选择PLANE82，该单元为八节点四边形

11、单元，但是可以退化为六节点三角形单元。设置好单元类型后，实常数设置板厚为1M。网格设置：对各个line进行网格控制，设置各线1、2、3、4、5的分割份数（Ndiv）分别为20、20、15、40、40。网格划分方式为Tri+free。划分后的网格结果如图3-3。约束的施加方式和载荷分布如图3-2中所示。分别约束线号为4和5上节点的Y和X方向自由度。在线号为1的节点上施加均布压力0.25MPa。经计算，常应变三节点三角形单元的位移云图如图3-4所示，总体应力云图（von mises stress）如图3-5所示，局部应力云图如图3-6所示。图3-3 三角形单元的网格划分图3-图3-4 位移云图图3

12、-5 总体应力云图图3-6 孔边周围局部应力云图六节点三角形单元的位移云图如图3-7所示，总体应力云图（von mises stress）如图3-8所示，局部应力云图如图3-9所示。图3-7 位移云图图3-8 总体应力云图图3-9 孔边周围局部应力云图经分析可知，不同的单元类型计算得来的应力和位移云图分布情况是基本一致的。最大位移在点1处附近，最小位移在点5处附近。最大的应力发生在点5处，最小应力发生在孔边上某段圆弧位置，该位置具体见各图中所示。DMXSMNSMXxy三节点0.309e-40.396e-50.309e-4-0.6960.259六节点0.309e-40.391e-50.309e-

13、4-0.7730.268理论值-0.750.25表3-1，计算结果对比表由上表3-1可以看出，在单元数目相同的情况，六节点三角形单元的分析精度要高于三节点常应变三角形单元。由此可得出结论，同样的单元形状和大小，高阶单元的计算精度要高于低阶单元。三、不同数量的三节点常应变单元的计算单元类型选择plane42，仅改变单元的分割段数，其他不变。即第一次加密设置各线1、2、3、4、5的分割份数（Ndiv）分别为30、30、20、60、60，第二次加密设置各线1、2、3、4、5的分割份数（Ndiv）分别为40、40、30、80、80。第一次加密后划分后网格为图3-10所示。计算后的位移云图如图3-11所

14、示，应力云图如图3-12所示。第二次加密后划分后网格为图3-13所示。计算后的位移云图如图3-14所示，应力云图如图3-15所示。图3-10 第一次加密后网格划分情况图3-11 第一次网格加密后位移云图图3-12 第一次网格加密后应力云图图3-13 加密后网格划分情况图3-14 网格加密后位移云图图3-15 网格加密后应力云图DMXSMNSMXxy不加密0.309e-40.396e-50.309e-4-0.6960.259第一次加密0.309e-40.391e-50.309e-4-0.7250.269第二次加密0.309e-40.395e-50.309e-4-0.7790.266理论值-0.7

15、50.25表2 不同密度的网格计算结果对比表由上表可知，虽然常应变三角形单元的计算结果没有高阶单元的计算结果精确，但是随着单元数目的增多，计算结果逐渐的接近高阶单元的计算结果。但是随着单元数量的增加，其计算的时间也会相应的增加。三、方孔和圆孔的应力水平比较为了使结果具有对比性，采用相同的单元类型plane82，并采用全局网格控制，大小为0.2。圆孔的计算结果如图3-16,17所示。方孔的计算结果如图3-18,19所示。圆孔的最大应力分布在圆孔与Y轴的交点附近，其大小为0.776MPa，而方孔的最大应力分布在四个角处，其大小为0.956 MPa，方孔的最大应力要比圆孔的大。并且圆孔边周围的应力变

16、化比较缓慢，而方孔变化剧烈。这些因数直接造成方孔板的安全系数没有圆孔板的大，比圆孔板易遭到破坏。图3-16 圆孔板的应力云图图3-17 圆孔边的应力云图图3-17 方孔的应力云图图3-17 方孔边的应力云图三、总结由以上的分析可以知道，在有限元分析当中，采用不同的单元（包括类型和数量），会得到不同的计算结果。相同的数量下，高阶单元的计算结果要比低阶单元的精确度高。相同的单元类型，单元的数量越多，得到的计算结果会越逼近精确值，但是计算的时间也会相应的增加。要做好有限元分析，得到需要精度的计算结果，必须得准确选择与实际力学模型相对应的单元类型，并要结合时间和硬件条件划分适合的网格数量。通过对圆孔板

17、和方孔板的计算结果的比较可知，在相同的边界条件和相同尺度的孔下，圆孔板的最大应力要小于方孔板的最大应力，方孔板在四角处的应力集中情况相当严重。在工程设计中，要避免这样的设计。四、附录圆孔板的单元为PLANE42，各线1、2、3、4、5的分割份数（Ndiv）分别为20、20、15、40、40的计算命令流：/FILNAME,shiti3/TITLE,ANALYSIS OF PLATE STRESS WITH SMALL CIRCLE/PREP7ET,l,PLANE42MP,EX,1,2.1E11MP,PRXY,1,0.3RECTNG,0,24,0,24PCIRC,3,0,0,90/PNUM,ARE

18、A,1 APLOT asba,1,2 NUMCMP,ALL /TITLE,GEOMETRIC MODELAPLOT /PNUM,KP,1LPLOTlesize,1,20lesize,2,20lesize,5,40lesize,4,40lesize,3,15mshape,1amesh,1/TITLE,ELEMENTS IN MODELEPLOTFINISH/SOLUANTYPE,STATIC/PNUM,LINE,1LPLOTLSEL,S,LINE, ,4NSLL,S,1D,ALL,UYLSEL,S,LINE,5NSLL,S,1D,ALL,UXLSEL,S,LINE,1NSLL,S,1SF,ALL

19、,PRES,-250000ALLSELOUTPR,BASIC,ALLSOLVEFINISH/POST1PLDISP,2PLNSOL,U,SUMPLNSOL,S,EQV/EXPAND,4,POLAR,HALF,90PLNSOL,S,EQVFINISH试题10:图示为钢涵洞，确定最大应力、最大位移及位置。E=210Gpa，=0.3图10-1 涵洞几何模型假如涵洞宽为1M，按空间问题进行计算，并和上述结果进行比较。同时，考虑若桥墩高由2M增加到3.5M，涵洞半径增加为无穷（即圆弧为直线）。计算最大应力，指出合理的桥洞形状曲线。一、按照平面问题计算按图10-1所示模型进行建模。并用PLANE82单元进

20、行划分网格，网格大小采用全局网格控制，划分方式采用自由方式，划分后网格的模型如图10-2所示。边界条件的施加为：约束左右两个下边的全部自由度，对涵洞上边施加均布压力70N/m。计算得到的位移和应力分布分别如图10-3、10-4所示。图10-2 涵洞的网格划分图10-3 涵洞的总体位移图图10-4 涵洞的总体应力图二、按照三维问题计算按图10-1所示平面模型，并拉伸宽度为1m进行建模。设置单元类型为solid186，控制截面上的面的大小，划分方式采用扫掠方式，划分后网格的模型如图10-5所示。边界条件的施加为：约束左右两个下表面的全部自由度，对涵洞上表面施加向下的均布压力70MPa。计算得到的位移和应力分布分别如图10-6、10-7所示。图10-5 三维涵洞的网格划分图10-6 三维涵洞的位移分布图10-7 三维涵洞的应力分布三、直角三维涵洞的计算将图10-1所示的桥墩高由2m增加到3.5m，涵洞半径增加为无穷（即圆弧为直线），宽度同样为1m。设