《太原理工大学数值传热4.2上机.doc》由会员分享，可在线阅读，更多相关《太原理工大学数值传热4.2上机.doc（9页珍藏版）》请在三一办公上搜索。
1、精选优质文档-倾情为你奉上数值传热学上机学院：电力学院班级：热能1101姓名：宋磊学号： 上机时间：2014.01.0642有一各向同性材料的方形物体，其导热系数为常数，已知各边界条件如图所示，试用高斯赛德尔迭代求其内部网格节点的温度。相对偏差为0.0002。解：2维稳态导热无内热源问题的控制微分方程t2/(x)2+t2/(y)2=0。500边界条件：上边界为500，其余边界温度均为100。 100100100积分后有其中k都相同，x=y=（x）e=（x）w=（y）n=（y）s最后有程序如下：#include#include#define M 4double main ()double aMM
2、=0,b1=0,b2=0,b3=0,b4=0,max=0;int i,k;for (i = 0; i M; i+) for (k = 0; k M; k+) scanf(%lf/n,&aik); /*初始数据输入*/ for (i = 0; i M; i+) for (k = 0; k M; k+) printf(%lf, aik); printf(n); /*初始数据输出*/dob1=a(M-1)/3(M-1)/3;b2=a(M-1)/32*(M-1)/3;b3=a2*(M-1)/3(M-1)/3;b4=a2*(M-1)/32*(M-1)/3;for (i = 1; i M-1; i+) f
3、or (k = 1; k M-1; k+) aik=(ai-1k+ai+1k+aik-1+aik+1)/4; /*算法迭代*/for (i = 0; i M; i+) for (k = 0; k (b2-a(M-1)/32*(M-1)/3) max=(b1-a(M-1)/3(M-1)/3)/b1;else max=(b2-a(M-1)/32*(M-1)/3)/b2;if (maxb3-a2*(M-1)/3(M-1)/3) max=(b3-a2*(M-1)/3(M-1)/3)/b3;if (max=0.0002); /*要求输出的精度*/return (0);运行结果16个节点时49个节点时将M的值改为7运行由以上运行结果可知该程序可有本题出发推广到一般。当网格划分不同时，只需要更改预处理中M的值即可。其中每个点的值分别代表相应节点的值。专心-专注-专业
