一维非稳态导热的数值计算.doc

《一维非稳态导热的数值计算.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、传热学C程序源二维稳态导热的数值计算2.1物理问题一矩形区域，其边长L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。2.2数学描述对上述问题的微分方程及其边界条件为：x=0，T=T1=0x=1，T=T1=0y=0，T=T1=0y=1，T=T2=1该问题的解析解：2.3数值离散2.3.1区域离散区域离散x方向总节点数为N，y方向总节点数为M，区域内任一节点用I,j表示。2.3.2方程的离散对于图中所有的内部节点方程可写为：用I,j节点的二阶中心差分代替上式中的二阶导数，得：上式整理成迭代形式：(i=2,3……,N-1)，(j=

2、2,3……,M-1)补充四个边界上的第一类边界条件得：(j=1,2,3……,M)(j=1,2,3……,M)(i=1,2,3……,N)(i=1,2,3……,N)传热学C程序源之一维非稳态导热的数值计算#include#include#defineN10#defineK11main(){ inti,j,l; floatcha; floata,x,y,Fo,Bi; floatt[N][K],b[N][K]; /*打印出题目*/ printf("ttt一维非稳态导热问题tt"); printf("tttttt----何鹏举"); prin

3、tf("题目：补充材料练习题三");  y=1;/*y代表Δτ*/ x=0.05/(N-1); a=34.89/(7800*712); Fo=(a*y)/(x*x); Bi=233*x/34.89; printf("显示格式条件:"); printf("1、Fo=%3.1f<0.5t",Fo); printf("t2、1-2Fo*Bi-2Fo=%4.2f>0",1-2*Fo*Bi-2*Fo); /*时刻为零时，赋予初场温度*/ for(i=0;i

4、   for(i=0;i0.001)   {        for(i=0;i

5、])+(1-2*Fo)*t[i][j];    t[N-1][j+1]=t[N-2][j]*(1-2*Fo*Bi-2*Fo)+2*Fo*t[N-1][j]+2*Fo*Bi*20;/*边界点温度用热平衡法推导出公式*/   }      cha=0;      for(i=0;i

6、j>=0;j--)  for(i=0;i999.99)    printf("%6.1f",t[i][j]);   else    printf("%6.2f",t[i][j]);   l=l+1;   if(l==N)   {    printf("");    l=0;   }  }  getchar();/*为了是生成的exe文件结果算的后不会立即退出，方便观看*/}