有限元报告——温度场..docVIP

有限元上机报告——温度场的有限元计算 一．问题 如图一平面结构在无热源情况下，给定热边界条件，用有限元分析温度分布。 二．解决步骤 对问题的分析 采用简单的三角形单元，单元内温度假定为线性分布，即 与平面结构一样，可用单元3个顶点的温度插值单元内部温度场，有 其中 为单元的节点温度列阵，而形状函数矩阵为 简单三角形单元内假定的温度场是线性分布的，其形状函数应为 对任一个单元，如面积域为，则单元泛函数为 而 所以，泛函数 单元刚度矩阵 所以 所以 2.数据准备 如图所示，划分单元格 每节点有一个自由度，边界约束为1,2,3,4,5,6,7,12,13,18,19,24,25,30,31,33,34,35,36，温度相当于载荷分布，所以只有边界处有载荷。和之前分析步骤相同，可得数据文件INP.DAT。 3.程序运行结果 改变边界条件 单元格不变，边界条件改变如下， 则程序的运行结果为 思考与讨论 1.分析采用该种单元分析平面温度场时是否可以收敛于真实解。 对同样的三角形单元，用节点温度插值单元内部温度的形状函数，与用节点位移插值单元内部位移的形状函数是完全一样的。与位移单元的分析一样，这种单元在单元交界处温度也是连续的，满足本问题的相容性要求。这种单元内部温度T(x,y)为完全的一次多项式，可以实现任意的常温度导数的温度状态，满足插值函数的完备性要求。因而，采用这种单元分析平面温度场时，有限元分析是可以收敛于真实解的。 附录： input.TXT 1,3,36,8,20,50,2 2.1,0.3,10., 0.,0.,20.,0.,40.,0.,60.,0.,80.,0.,100.,0., 0.,20.,20.,20.,40.,20.,60.,20.,80.,20.,100.,20., 0.,40.,20.,40.,40.,40.,60.,40.,80.,40.,100.,40., 0.,60.,20.,60.,40.,60.,60.,60.,80.,60.,100.,60., 0.,80.,20.,80.,40.,80.,60.,80.,80.,80.,100.,80., 0.,100.,20.,100.,40.,100.,60.,100.,80.,100.,100.,100., 0.,20.,40.,60.,80.,100., 20.,0.,0.,0.,0.,80., 40.,0.,0.,0.,0.,60., 60.,0.,0.,0.,0.,40., 80.,0.,0.,0.,0.,20., 100.,80.,60.,40.,20.,0., 1,2,8,2,3,9,3,4,10,4,5,11,5,6,12,1,8,7,2,9,8,3,10,9,4,11,10,5,12,11, 7,8,14,8,9,15,9,10,16,10,11,17,11,12,18,7,14,13,8,15,14,9,16,15,10,17,16,11,18,17, 13,14,20,14,15,21,15,16,22,16,17,23,17,18,24,13,20,19,14,21,20,15,22,21,16,23,22,17,24,23, 19,20,26,20,21,27,21,22,28,22,23,29,23,24,30,19,26,25,20,27,26,21,28,27,22,29,28,23,30,29, 25,26,32,26,27,33,27,28,34,28,29,35,29,30,36,25,32,31,26,33,32,27,34,33,28,35,34,29,36,35, 1,2,3,4,5,6,7,12,13,18,19,24,25,30,31,32,33,34,35,36, PLANE.FOR PROGRAM MAIN DIMENSION SK(300,30),EK(12,12),Q(300),MC(55),XY(2,100),XYE(2,4),QE *(12),NX(4,100) OPEN (7,FILE=input.TXT) REWIND 7 READ (7,*) NF,NE,NN,MB,ND,LE,LS READ (7,*) E,UM,T 10 FORMAT (7I5) 12 FORMAT (3F15.2) WRITE(*,600) NF,NE,NN,MB,ND,LE,LS,E,UM,T ME=NE*NF MS=NN*NF CALL INPUT (XY,Q,NX,MC,LS,NN,MS,NE,LE,ND) WRITE(*,102)((XY(I,J),I=1,LS),J=1,NN) 102 FORMAT (10X,XY/,(2X,6F12.3))