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

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

1、二维稳态导热的数值计算一、实验名称：二维稳态导热的数值计算二、实验内容图2示出了一矩形区域，其边长L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。三、程序代码#includevoidmain(){inti,j;floatT[6][7];floatt[6][7]={1,1,1,1,1,1,1};loop:for(i=0;i<=5;i++)for(j=0;j<=6;j++)T[i][j]=t[i][j];for(i=1;i<=4;i++)for(j=1;j<=5;j++)t[i][j]=((25*t[i-1][

2、j]+25*t[i+1][j]+36*t[i][j-1]+36*t[i][j+1])/122);for(i=1;i<=4;i++)for(j=1;j<=5;j++){if((t[i][j]-T[i][j])/t[i][j]>0.)gotoloop;}for(i=0;i<=5;i++){for(j=0;j<=6;j++)printf("%f,",t[i][j]);printf("");}}四、运行结果