传热学上机编程题答案.docVIP

传热学上机 高牥杰 2012/11/3 4-2． #includeiostream using namespace std; void main() { int i=0,j=0,M=9,N=7,K,IT,B; int p=1; double EPS,TTB,TRB,TLB,TBB,TT,eps,TI,A; cout请输入K,EPS,TI,TTB,TRB,TLB,TBB:endl; cinKEPSTITTBTRBTLBTBB; double T[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++) //内部点赋初值TI for(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． #includeiostream using namespace std; void main() { int N,t,NP,TM,IT,M; double Bi,Fo,d,tf,a,h,y,dx,T0; double T[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(int i=0;i11;i++) T[i]=T0; cout0 second:; for(i=0;i11;i++) coutT[i] ; coutendl; for(IT=1;IT*t=TM;IT++) { M=IT*t; for(int i=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) { cout