传热学实验墙角matlab导热问题.docxVIP

二维导热物体温度场的数值模拟 姓名 小明 学号 111111 班级 能动 学院 能动 一、 问题描述 有一墙角模型，尺寸如图 1所示，导热系数 0.53W/(m ? K)，墙角内外壁为第一类边界 条件。求解该模型的温度分布及导热量。 图1 q=0 q=0 二、计算原理 根据热平衡法列出节点方程， 各方向导入单元体的热量之和为零。 内节点和绝热边界点 (图1点划线上的点)的方程形式不同。 图2 tuH li-! j o A X J x 円节点 绝热边界 图2所示的内节点和绝热边界节点方程如下: 内节点： 5 5 ?1?虹 ti,j) X (tj, ti,j)x (tjij lj)y (S,j tGy 0 绝热边界点: ?1?(j ti,j)x W 2y (ti,j1 kj) x 0 (ti 2y — ij kJ y o 三、计算过程 用Matlab7.1语言编写计算程序，初取网格步长 x y 0.1m HTBfun亡tion I HTB for i=l;16 IU 1) = 30; end for j=l:12 Hl, j) = 30h end for in=6; 16 for n=6：12 I % n) = 0; end 1(16, 2:5)= [24, 18, 12, C]: B±HTH1=2：15T(2:5f 12)= [24, 18, 12, S] B±HTH 1=2：1