分数:有限元大作业院系:专业:年级:学号:姓名:2017年12月13日第一题21210MPa,0.2m,,,0.01m ,210GPa,0.3.(2)p p AB CD EF ==等厚度薄板受到分均布载荷其厚度为杆横截面积为所有材料的弹性模量为泊松比为试用有限元法对该系统进行离散并进行单元和节点编号.(1)写出离散后各单元的刚度矩阵;写出整体刚度矩阵和刚度方程;(3)求解各节点的节点位移;(4)求解支座反力;(5)用ANSYS 软件验证以上求解结果并写出详细的命令流程.一、操作命令/NOPR ! Suppress printing of UNDO process/PMACRO ! Echo following commands to logFINISH ! Make sure we are at BEGIN level/CLEAR,NOSTART ! Clear model since no SAVE found! WE SUGGEST YOU REMOVE THIS LINE AND THE FOLLOWING STARTUP LINES /input,menust,tmp,''/GRA,POWER/GST,ON/PLO,INFO,3/GRO,CURL,ON /CPLANE,1/REPLOT,RESIZE WPSTYLE,,,,,,,,0!*/NOPRKEYW,PR_SET,1KEYW,PR_STRUC,1KEYW,PR_THERM,0KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0KEYW,MAGEDG,0KEYW,MAGHFE,0KEYW,MAGELC,0KEYW,PR_MULTI,0/GO!*!*/PREP7!*ET,1,BEAM188!*ET,2,PLANE183!*KEYOPT,2,1,0KEYOPT,2,3,3KEYOPT,2,6,0!*!*R,1,0.2,!*!*MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,210e9MPDATA,PRXY,1,,0.3SECTYPE, 1, BEAM, RECT, , 0 SECOFFSET, CENTSECDATA,0.1,0.1,0,0,0,0,0,0,0,0,0,0 K,1,,,,K,2,3,,,K,3,3,4,,K,4,3,8,,K,5,9,8,,K,6,9,4,,K,7,12,0,,K,8,3,4,,K,9,9,4,,/DIST,1,1.08222638492,1/REP,FASTLSTR, 1, 3LSTR, 2, 3LSTR, 7, 6LSTR, 8, 9LSTR, 9, 5LSTR, 5, 4LSTR, 4, 8FLST,2,4,4FITEM,2,4FITEM,2,5FITEM,2,6FITEM,2,7AL,P51XCM,_Y,AREAASEL, , , , 1CM,_Y1,AREACMSEL,S,_Y!*CMSEL,S,_Y1AATT, 1, 1, 2, 0, 1 CMSEL,S,_YCMDELE,_YCMDELE,_Y1!*FLST,2,1,5,ORDE,1FITEM,2,1AESIZE,P51X,0.6,MSHAPE,0,2DMSHKEY,0!*CM,_Y,AREAASEL, , , , 1CM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_Y!*AMESH,_Y1!*CMDELE,_YCMDELE,_Y1CMDELE,_Y2!*FLST,5,7,4,ORDE,2 FITEM,5,1FITEM,5,-7CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,S,_Y!*!*CMSEL,S,_Y1LATT,1,1,1, , , ,1 CMSEL,S,_YCMDELE,_YCMDELE,_Y1!*FLST,5,7,4,ORDE,2 FITEM,5,1FITEM,5,-7CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,,_Y!*LESIZE,_Y1,0.6, , , , , , ,1 !*FLST,2,7,4,ORDE,2 FITEM,2,1FITEM,2,-7LMESH,P51X/PNUM,KP,0/PNUM,LINE,0/PNUM,AREA,0/PNUM,VOLU,0/PNUM,NODE,1/PNUM,TABN,0/PNUM,SVAL,0/NUMBER,0!*/PNUM,ELEM,0/REPLOT!*/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FASTNPLOT/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,0.924021086472,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1 /REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FAST/DIST,1,1.08222638492,1/REP,FASTFLST,4,2,1,ORDE,2FITEM,4,1FITEM,4,247CP,2,ALL,P51XFLST,4,2,1,ORDE,2FITEM,4,2FITEM,4,264CP,21,ALL,P51XFLST,2,3,3,ORDE,3FITEM,2,1FITEM,2,-2FITEM,2,7!*/GODK,P51X, , , ,0,ALL, , , , , ,FLST,2,2,4,ORDE,2FITEM,2,6FITEM,2,-7/GO!*SFL,P51X,PRES,1e7,FINISH/SOL/STATUS,SOLUSOLVE)/GOP ! Resume printing after UNDO process )! We suggest a save at this point二、计算结果支座反力NODE u v1 -0.14895E+07 -0.19808E+072 0 0.53105E+073 -0.65044E+07 0.86703E+07节点位移NODE UX UY UZ USUM1 0.0000 0.0000 0.0000 0.0000理论计算节点113 0.27034E-01-0.10597E-01 0.0000 0.29036E-01理论界点721 0.26855E-01-0.15443E-01 0.0000 0.30978E-01理论节点622 0.27018E-01-0.11019E-01 0.0000 0.29179E-01理论节点532 0.26850E-01-0.15049E-01 0.0000 0.30780E-01理论节点453 0.0000 0.0000 0.0000 0.0000理论计算节点363 0.0000 0.0000 0.0000 0.0000理论节点2第二题LNG储罐罐体温度场数值模拟过程1.罐底温度分析该部分研究的对象是储罐的罐底,从几何形状上来讲满足轴对称;罐底内壁温度即为底部LNG 温度,该温度对应罐底静压力下的饱和温度,运用克劳修斯-克拉贝龙方程进行计算得到底部LNG 的温度为122.37 K;外边界为对流换热边界,环境温度为20 ℃。
2.求解步骤第一步:建立工作文件名和工作标题(1)Utility Menu:File>Change Job name,输入文件名GUAN DI(2)Utility Menu:File>Change Title, 输入GUAN DI WEN DU FEN XI;(3) Main Menu:Preferences> Preferences for GUI Filtering,选择Thermal,单击OK关闭对话框。
(4)Main Menu:Plot Ctrls >Style>Colors>Reverse Video 将背景颜色改成白色;第二步:定义单元类型(1)Utility Menu:Preferences>Element Type>Add/Edit/Delete,出现Element Type对话框。
(2)单击Add按钮,出现Library of Element Type对话框。
在Library of Element Type 列表中选择Thermal Solid,Quad 4node 55,在Library of Element Type输入框中输入1。
单击OK按钮,关闭对话框。
(3)单击Element Type对话框中的Options按钮,出现PLANE55 Element Type options 对话框,在Element behavior k3下拉列表框中选择axisymmetric,其余选项均采用默认设置,如下图1所示,单击OK关闭该对话框。