热传导方程的c语言程序,传热学上机C程序源答案之一维非稳态导热的数值计算...

《传热学上机C程序源答案之一维非稳态导热的数值计算》由会员分享，可在线阅读，更多相关《传热学上机C程序源答案之一维非稳态导热的数值计算(4页珍藏版)》请在人人文库网上搜索。

1、二维稳态导热的数值计算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,3,M-1)。

2、补充四个边界上的第一类边界条件得： (j=1,2,3,M)(j=1,2,3,M)(i=1,2,3,N)(i=1,2,3,N)传热学C程序源之一维非稳态导热的数值计算#include#include#define N 10#define K 11main()int i,j,l;float cha;float a,x,y,Fo,Bi;float tNK,bNK;/*打印出题目*/printf(ttt一维非稳态导热问题tt);printf(ntttttt-何鹏举n);printf(n题目：补充材料练习题三n);y=1;/*y代表*/x=0.05/(N-1);a=34.89/(7800*712);Fo=(a*y)/(x*x);Bi=233*x/34.89;printf(n显示格式条件:);printf(n1、Fo=%3.1f0nn,1-2*Fo*Bi-2*Fo);/*时刻为零时，赋予初场温度*/for(i=0;i0.001) for(i=0;i=0;j-)for(i=0;i999.99)printf(%6.1f ,tij);elseprintf(%6.2f ,tij);l=l+1;if(l=N)printf(n);l=0;getchar();/*为了是生成的exe文件结果算的后不会立即退出，方便观看*。