搜索
传热学上机高牥杰2012/11/3传热学上机4-2．#includeiostreamusingnamespacestd;voidmain(){inti=0,j=0,M=9,N=7,K,IT,B;intp=1;doubleEPS,TTB,TRB,TLB,TBB,TT,eps,TI,A;cout请输入K,EPS,TI,TTB,TRB,TLB,TBB:endl;cinKEPSTITTBTRBTLBTBB;doubleT[9][7];for(i=0;iN;i++)//顶部底部赋初始值{T[0][i]=TLB;T[M-1][i]=TRB;}coutendl;for(j=1;jM;j++)//两侧赋初始值{T[j][0]=TBB;T[j][N-1]=TTB;}for(i=1;iM-1;i++)//内部点赋初值TIfor(j=1;jN-1;j++)T[i][j]=TI;for(j=N-1;j=0;j--)B=p;for(IT=0;ITK;IT++){eps=T[1][N-2]-0.25*(T[0][N-2]+T[1][N-3]+T[2][N-2]+T[1][N-1]);for(i=1;iM-1;i++){for(j=N-2;j=1;j--){TT=T[i][j];T[i][j]=0.25*(T[i-1][j]+T[i][j-1]+T[i+1][j]+T[i][j+1]);A=T[i][j]-TT;if(eps0)eps=-eps;if(A0)A=-A;if(eps=A)eps=A;}}if(eps=EPS){for(j=N-1;j=0;j--){for(i=0;iM;i++){coutT[i][j]'\t';//打印内部点初值传热学上机if(p%9==0)coutendl;p++;}}coutIT=IT;break;}}if(B==p)cout不收敛ITendl;}4-4．#includeiostreamusingnamespacestd;voidmain(){intN,t,NP,TM,IT,M;doubleBi,Fo,d,tf,a,h,y,dx,T0;doubleT[11],TT[11];cout请输入壁厚d:;cind;cout请输入半壁厚等分数N：;cinN;cout请输入控制打印节点节点温度的时间间隔NP:;cinNP;传热学上机cout请输入流体温度tf:;cintf;cout请输入时间间隔t:;cint;cout请输入终止计算时间TM:;cinTM;cout请输入表面传热系数h:;cinh;cout请输入平壁导热系数y:;ciny;cout请输入初始温度T0:;cinT0;cout请输入热扩散率a：;cina;dx=d/(2.0*(double)N);Bi=h*dx/y;Fo=a*t/(dx*dx);coutBi''Foendl;if(Fo1.0/(2.0*Bi+2.0))cout不稳定endl;else{for(inti=0;i11;i++)T[i]=T0;cout0second:;for(i=0;i11;i++)coutT[i];coutendl;for(IT=1;IT*t=TM;IT++){M=IT*t;for(inti=9;i=1;i--){TT[i]=Fo*(T[i-1]+T[i+1])+(1.0-2.0*Fo)*T[i];}TT[N]=2.0*Fo*(T[N-1]+tf*Bi)+(1.0-2.0*Fo-2.0*Fo*Bi)*T[N];TT[0]=2.0*Fo*T[1]+(1.0-2.0*Fo)*T[0];for(i=0;i11;i++)T[i]=TT[i];if(IT%NP==0){coutMsecond：;for(inti=0;i=10;i++){coutT[i];}coutendl;传热学上机}}Fo=a*TM/((d/2.0)*(d/2.0));coutBi''Foendl;}}建环1003班U201015589高牥杰