传热学上机编程题3.doc

1、传热学上机艾科热木江艾尔肯建环 1003 u201015592传热学上机4-2#includeusing 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;coutKEPSTITTBTRBTLBTBB;double T97;for(i=0;i=0;j-)B=p;for(IT=0;IT=1;j-) TT=Tij;Tij=0.25*(Ti-1j+Tij-1+Ti+1j+Tij+1);A=Tij-TT;if(eps=0;j-)传热学上机 for(i=0;

2、iusing namespace std;void main()int N,t,NP,TM,IT,M;double Bi,Fo,d,tf,a,h,y,dx,T0;double T11,TT11;coutd;coutN;传热学上机coutNP;couttf;coutt;coutTM;couth;couty;coutT0;couta;dx=d/(2.0*(double)N);Bi=h*dx/y;Fo=a*t/(dx*dx);cout1.0/(2.0*Bi+2.0)cout=1;i-)TTi=Fo*(Ti-1+Ti+1)+(1.0-2.0*Fo)*Ti;TTN=2.0*Fo*(TN-1+tf*Bi)+(1.0-2.0*Fo-2.0*Fo*Bi)*TN;TT0=2.0*Fo*T1+(1.0-2.0*Fo)*T0;for(i=0;i11;i+)Ti=TTi;if(IT%NP=0) coutM“ second：“;for(int i=0;i=10;i+)coutTi“ “;传热学上机 coutendl; Fo=a*TM/(d/2.0)*(d/2.0);coutBi Foendl;