有限元上机实验报告材料

1、有限元法基础及应用上机报告2015南京理工大学年12月上机实验一1实验题目设计一个采用减缩积分线性四边形等参元的有限元模型，通过数值试验来研 究网格密度、位移约束条件与总刚度矩阵奇异性、沙漏扩展、求解精度的关系， 并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。总结出你的研究结论，撰写实验报告。2实验目的通过实验来研究减缩积分方案中网格密度和位移约束条件对总体刚度矩阵 奇异性和求解精度的影响，以此加深对有限元减缩积分的理解，和对减缩积分中 保证总体刚度矩阵非奇异性的认识。3建模概述先保持位移约束条件不变，研究网格密度对总体刚度矩阵奇异性和求解精度的影响，并验证采用减缩积分时保证总刚度矩阵非

2、奇异的必要条件。如下图1所示，建立一个简支和链杆的约束条件，然后不断增加网格密度，通过 ABAQUS 来计算位移和应力的变化规律。简支（两个约束）积分点（3个独立关系式）链杆（一个约束）°节点（两个自由度）4计算结果分析讨论与结论1) 1*1单元四边形减缩积分实验载荷布种/单元应力云图2) 2*1单元四边形减缩积分实验载荷单元/?£，的也好 | 熊Saw3n -|' wrefnftfWb-i 配 Mnr5Arg. 75&E* LildsMSO*1 3Z3PDO pUZIBpOO 1总善心i H .嘱谑M- *1.333d*CO «U?3a£

3、;jQ * *1.J； Jw*dO *1 3Z3e*CO - I r -n- - I T建韭TO 日:1卷:第应力云图3) 4*4单元四边形减缩积分实验n- u nu oo u n Q u n u n 1 -IHI- 1 -1 ,e+时B0+0+0 Hs+iie+JAB- 974129419 6317J5m Jf 冉之JQ，1 J 由 33 3s22? 2 11119应力云图结果分析项目自由度数比较独立关系式结果1)4*2-3=5>1*3=3奇异2)6*2-3=9>2*3=6奇异3)25*2-3=47<16*3=48非奇异5实验体会与小结单元刚度矩阵的特征：(1)对称性(2)

4、奇异性(3)主元恒正 e .(4)平面图形相似、弹性矩阵 D厚度t相同的单兀，K相同(5) K e的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两 歹I，其位置与结点位置对应。整体刚度矩阵的特征：(1)对称性(2)奇异性(3)主元恒正(4)稀疏性(5)非零元素呈带状分布。K的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足 力和力矩的平衡。为消除K的奇异性，需要引入边界条件，至少需给出能限制 刚体位移的约束条件。对于一个给定形式的单元，如果采用精确积分，则插值函数中所有项次在 |J|=常数的条件下能被精确积分，并能保证刚度矩阵的非奇异性。如果采用减缩 积分，因为插值函数中

5、只有完全多项式的项次能被精确积分，因此需要进行刚度矩阵非奇异必要条件的检查。上机实验二1实验题目图示一个简支梁平面应力模型。梁截面为矩形，高度h=160mm长度L=1000mm 厚度t=10mm 上边承受均布压力 q =1N/mn,材料E=206GPa = =0.29 0 X方向 正应力弹性力学理论解为：22二x = -3( -x2)y q-(4-y2 -)h3 4h h2 5分别应用3节点三角形单元、4节点线性等参元(完全积分、减缩积分、非 协调模式)、8节点二次等参元完全积分进行下列数值实验：1)用较粗单元网格 求解梁中部应力分量6x的最大值和上下边法向应力分量， 并对各单元计算精度 进行

6、比较分析；2)对粗网格下梁中部最大位移进行对比和分析；3)通过网格加 密对比试验3节点三角形单元和8节点二次等参元的收敛速度。总结出研究结论， 撰写实验报告。2实验目的(1)通过ABAQUSC件，用有限元法分析整个梁上的6_x的分布规律，讨论 b_y的有限元解与材料力学解的区别；(2)用有限元法求梁顶边和底边中点正应力(7 x的最大值；(3)比较3节点三角形单元、4节点线性等参元(完全积分、减缩积分、非 协调模式)、8节点二次等参元完全积分求得的参数与理论解的接近程度；(4)逐步加密单元网格，把有限元法求得的 x值与理论值进行对比，考察 解的收敛性；(5)针对以上力学模型，对比分析 3节点三角

7、形平面单元和8节点四边形平 面单元的求解精度和收敛性；(6)绘制x -误差一一计算次数曲线，并进行分析说明。3建模概述(1)启用 ABAQUS/CAE序；(2)创建部件(Module： Part ),选择 2DPlanar , Approximate size=2000 , 绘制长度为1000,高度为200的矩形；(3)创建材料和截面属性(Module： Property ),弹性模量为E=210000MPa 泊松比为0.3 ;(4)定义装配件(Module： Assembly),选择 Independent ;(5)设置分析步(Module： Step);(6)定义边界条件和载荷(Modul

8、e： Load),将左右两边界在中点分成两段， 在梁的上侧面施加1Mpa6勺载荷，选取左边界中点设置 U1=0和U2=0,选 取右边界中点设置U2=Q(7)划分网格(Module： Mesh),全局尺寸设为50,采用3节点线性平面应 力三角形单元CPS3aJ分网格；(8)提交分析作业(Module： Job);(9)后处理(Module： Visualization ):显示的应力云图；查询底边中点的最大值，与理论最大应力值比较考察其精度；绘制底边上各点的应力的曲线；(10)细化网格验证收敛性，重设单元尺寸为 20和10,给出和的应力云图， 考察收敛于理论解的程度；(11)高阶单元分析，将单元

9、类型改为CPS8单元尺寸分别取为100、50、20: 绘制和的云图；查询底边中点的最大值，对比(10)和(11)中的结果，体会高精 度有限元网格的精度和收敛性； 查询模型中从上至下的值。4计算结果分析讨论与结论(1)尺寸大小为50的3节点三角形单元计算云图如下：1) x方向应力云图3Fro U型口4 l£：:42i4Ei GMT4-OB: DD 3015tV.?y iju5一t+.TaBe-ai7-7-t/小二 'I r I -L.3KE-+CJ梁中部x方向应力分量最大值为：19.57514HqdtA -3 nnoe-oiI Sr-，川 -2 £21&-b&

10、#163;)ai riHir-tpiT 了脸 Q祠ITHEEO 4- 74=;&+in tog ju2) y方向应力云图上边界二'方向应力分布曲线:皿.上边界法向最大应力：-1.54372下边界二、方向应力分布曲线:一.UXLTrue (trance along 网卜E, =-1 l Kge 冶: Trua 0工 口口丁0400,oOQ.Tine dislarHje ateng kth：cm.9( 611 5口 75S)： T<UB B5t Qq PatMIDOD.下边界法向最大应力：0.5437243) y方向位移云图梁中部y方向最大位移：-0.157743(2)尺寸大

11、小为50的四节点完全积分 计算云图如下：1) x方向应力云图梁中部x方向应力分量最大值为：28.69352) y方向应力云图上边界,r方向应力分布曲线:梁中部y方向最大位移：-0.188218-曲15.8 口口9QD.1040.Troe distance alortg pathECl 飞飞“ Itg Dut -sfcng上边界法向最大应力：-3.1434下边界法向最大应力：2.1434下边界二'方向应力分布曲线:3）位移云图(3)尺寸大小为50的四节点减缩积分1) x方向应力云图4 工 mg 4 m ,丽，中g-l a-L 3MC+OI-I 45+n i梁中部x方向应力分量最大值为：2

12、3.33512) y方向应力云图上边界二、方向应力分布曲线:5rSL3 (A)i-a： :ra需仁 TTue Det. -stnc PJth-5-上边界法向最大应力：-1.04348下边界、x方向应力分布曲线:梁中部x方向应力分量最大值为：29.27661?. EC1 衿p Ttipp Drt 却cmg n?s-sc-da»n'下边界法向最大应力：0.0434821位移云图一 2P;.'J11.11.11111 h<la在心町七山与屯 9 _K c fl M H Q fl 3 E M- Nncr-Jnr* LI L L Fa FJ2 用甥器幽sffi哥 in梁中

13、部y方向最大位移：-0.240197(4)尺寸大小为50的四节点非协调模式积分1) x方向应力云图2) y方向应力云图上边界；'x方向应力分布曲线:True along pathE,l m7字%I: Tiue Di:t 131cli口,qh什lq上边界法向最大应力：-0.900612下边界二、方向应力分布曲线:边界法向最大应力：-0.09938853)位移云图七 HE .* 口 C 4。二 H F 嗯媵窟盥学就皿 '1£0-54.4-11 -n7? !a “-_I - I -L X « IL Edu I h I H i i%一 口匚：DG 1 口 2gH G

14、MF+口8.00 1口157D& 亡口三日-*：=： J.dbV 叠印 5WP-1LETvrniHrf; 1; btat Trn«Rnrwv 的工 six a -j . * . di.& 51173,1梁中部y方向最大位移：-0.196154(5)尺寸大小为50的八节点二次等参元完全积分1) x方向应力云图十二跖工3旬i 二曲肛su十LdfT二4口 1 k q-。* j +专底3二口J%W2 口I -.-JUC 4 寸T944口口 -1970+01. -： * .-5 5=h+i：l2) y方向应力云图梁中部x方向应力分量最大值为:29.5572gfwiu ¥

15、;Q/b *5产"/5. 522lA>a. 73%J置力 二桀花94。1:.山*， 41 soge+ai 4 L J I JEI - 订工皿口胸：2 43 <9 "W -SJDODe-DI 吊汨4】。出？9誓g 7.门3“厂上边界,r方向应力分布曲线:梁中部y方向最大位移：-0.198043上边界法向最大应力：-1.04859下边界；：x方向应力分布曲线:SCO.ODDHR,LDODTrue dstmnee atongi pathE 3 L lAnia： E% I： Irje Dud . ahng acp -5n-dcr*r下边界法向最大应力：0.0485932

16、3）位移云图修sJrTii rBI2(!rl 1 i i z-aJI-an-llA-aQJ-aJ-a Ms"痔f.lFT 辑心(s.l对 +OT -4CI_-y-g441-g-1(6)尺寸大小为20的3节点三角形单元& SE.1WB- CPGMS的1: E "m51,圈,Qs，v1) x万向应力云图心如*01 包；知*01 M jrn-o i j.H.Cla+OO “ S：abd iXi T OUKTU 货空F Ik -1 - 凸丁12*肛 IU -2.1： WQ41J 弓甘 口片梁中部x方向应力分量最大值为：25.68172) y方向应力云图上边界法向最大应力：-

17、1.35765上边界,x方向应力分布曲线:下边界、x方向应力分布曲线:True Ustancft along path5, Sil 8Ttijp Rst HtaH '<5533-3D-drawn"下边界法向最大应力：0.3576523)位移云图Mjjwomx三而Ze梁中部y方向最大位移：-0.187793(7)尺寸大小为10的3节点三角形单元1) x方向应力云图梁中部x方向应力分量最大值为：27.9066-看 EM ±&bA-+hai-l-*+gr *等芳二.-.? 比.T二W+JM:1*2“£1Wri irii o+i doac 注喻2)

18、y方向应力云图4 . .1Sun gq十网5-3r*&O5e+0上边界；:x方向应力分布曲线:上边界法向最大应力：-1.19271下边界二、方向应力分布曲线:下边界法向最大应力：0.1927053)位移云图IMlalssis 黑片铝 .Q62JB弓JjJ 4 h-5 9LL LLLi '.'I.,SEtEH1>DB CP55-l：14&JPlC = .n-LZ.g3b 口上£q£, =ta-Ja -d §Psfiiij- 7工.L J/DefcHMMU梁中部y方向最大位移：-0.195956(8)尺寸大小为20的8节点二次等参

19、元1) x方向应力云图梁中部x方向应力分量最大值为：29.50852) y方向应力云图E LO'J.IT，4Nr*口-5 ylu迪vq凸、)45.341»01 T 44,+lJ t <E43tn小.小弓Euk k -1，,1.1.4444-09 - n(ne-Dl 我场近m .E47n40L上边界,r方向应力分布曲线:梁中部y方向最大位移：-0.1993155 1 lAra： TiueDrt. JhE qns»-2EHjp上边界法向最大应力：-1.01358下边界；：x方向应力分布曲线:5.S l'ia： True Oidr. sbnn qzsS-ZD

20、-dc*r下边界法向最大应力：0.01358293）位移云图(9)尺寸大小为10的8节点四边形等参单元1) x方向应力云图梁中部x方向应力分量最大值为：27.59812) y方向应力云图“叶臼I L吟HU-5OD0*-0i I 11 -'' ilj IWMIK “fiM0 TlM上边界：x方向应力分布曲线:5, SllCfik-a- T.|- T«ub C»5t &口ng'rp9日-lOtcf上边界法向最大应力：-1.00365下边界,r方向应力分布曲线:下边界法向最大应力：0.003645013）位移云图梁中部y方向最大位移：-0.2003

21、94单元尺寸为50的不同类型单元有限元分析结果:单元类型单元尺寸6 Xo- x( %上边法向应 力分量下边法向应 力分量梁中部最大位移理论解29.4969CPS3CPS4505019.575133.63675505-1.543720.543724-0.15774328.69352.723676047-3.14342.1434-0.188218CPS4R5023.335120.88965281-1.043480.0434821-0.240197CPS4I5029.27660.746858144-0.900612 -0.0993885-0.196154CPS85029.55720.20442826

22、2-1.048590.0485932-0.198043相同类型单元不同网格密度下的有限元分析对比:单元单元尺寸(TX(TX( %上边法向应 力分量下边法向应 力分量梁中部最大位移理论解29.4969CPS3CPS3502019.575133.63675505-1.543720.543724-0.15774325.68172.723676047-1.357650.357652-0.187793CPS31027.906620.88965281-1.192710.192705-0.195956CPS85029.55720.204428262-1.048590.0485932-0.198043CPS8

23、2029.508512.93424055-1.013580.0135829-0.199315CPS81027.59810.039326166-1.003650.00364501-0.200394在同等网格密度下计算水平方向应力分量时，二次单元的计算精度要优于线性单元。在单元阶数相同的情况下，三角形单元计算精 度次于四节点单元。 止匕外，减缩积分的线性四节点单元因为存在沙漏模式，其应力的计算精度要比同为线性四节点单元的完全积分单元和非协调模式的单元低。在同等网格密度下计算法向应力分量时， 减缩积分的线性四节点单元的计 算精度最高，其次依次是八节点二次单元、非协调模式的线 性四节点单元、三 角形单

24、元、完全积分的线性四节点单元。用八节点二次单元和线性四节点非协调单元计算应力是比较合适的方案。在同等网格密度下计算位移，八节点二次单元的精度最高，非协调模式的线 性四节点单元的精度次之，这两种单元的误差都在 5%内，而其他单元的计算 精度都偏大。5实验体会与小结通过本次实验，我基本掌握了使用 ABAQU软件建模求解平面应力问题，学 会了查询应力以及位移分布情况的方法。了解了单元尺寸变化和单元类型不同对 结果精度的影响。单元尺寸越小，精度越高，收敛性越好。高精度单元收敛速度 快，运算精度高。但要根据实际情况选择合理的单元，提高计算效率。上机实验三1实验题目一个矩形平板，长1000mm宽100mm厚度5mm材料的E=200GPa k = °.3 , P = 7.82M10-6Kg/mm3。板的一对短边简支。进行下列计算分析并撰写实验报告。1）在相同的粗网格（厚度方向1层单元）下，分别用8节点线性六面体全 积分等参元、8节点六面体非协调元、20节点二次六面体等参元计算其前六阶自 由振动频率（Hz）,对计算结果用表格作汇总对比和分析（不要附图）。2）采用适当网格密度的8节点六面体非协调元，用隐式直接积分法对该平 板进行瞬态响应研究。板的上表面受对称三角形脉冲均布压力，最大值0.05MPa 要求：（1）通过计算获得平板下表面中心位置沿长边方向的正应力分