二维肋片稳态导热问题数值计算.docVIP

例题4-5 二维肋片稳态导热问题的数值计算 （1）自主编程，编程语言自定，最后提交源程序 （2）提交电子报告（word格式），包括： （a）给出空间离散示意图（网格划分） （b）节点离散方程 （c） 图示温度等值线（可以利用origin或matlab） 解：（a）空间离散示意图(由origin的graph作图) （b）节点离散方程 由上图所示得各节点的节点离散方程 结点1：T(m,n)=0.25*(T(m+1,n)+T(m-1,n)+T(m,n+1)T(m,n-1)) 结点2：T(M,n)= 1/(4+2*Bi)*(T(M,n-1)+T(M,n+1)+2*T(M-1,n)) 结点3：T(m,N)=1/(4+2*Bi)*(T(m-1,N)+T(m+1,N)+2*T(m,N-1)) 结点4：T(M,N)= 1/(2+2*Bi)*(T(M-1,1)+T(M,2)) 结点5：T(m,1)=0.25*(T(m-1,1)+T(m+1,1)+2*T(m,2)) 结点6：T(M,1)= 1/(2+2*Bi)*(T(M,N-1)+T(M,2)) （c） 温度等值线 等温线图，工况1（Bi=0.01）η= 0.9670 温度与y轴分布图（工况1） 工况2（Bi=1）η=0.1910 温度与y轴分布图（工况2） 图像分析：四幅图的显示来看，结