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

传热学上机
传热学上机

高牥杰
2012/11/3

4-2．
#include<iostream>
usingnamespacestd;
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;
cin>>K>>EPS>>TI>>TTB>>TRB>>TLB>>TBB;
doubleT[9][7];
for(i=0;i<N;i++)//顶部底部赋初始值
{
T[0][i]=TLB;
T[M-1][i]=TRB;
}
cout<<endl;
for(j=1;j<M;j++)//两侧赋初始值
{
T[j][0]=TBB;
T[j][N-1]=TTB;
}
for(i=1;i<M-1;i++)//内部点赋初值TI
for(j=1;j<N-1;j++)
T[i][j]=TI;
for(j=N-1;j>=0;j--)
B=p;
for(IT=0;IT<K;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;i<M-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(eps<0)eps=-eps;
if(A<0)A=-A;
if(eps<=A)eps=A;
}
}
if(eps<=EPS)
{for(j=N-1;j>=0;j--)
{for(i=0;i<M;i++)
{cout<<T[i][j]<<'\t';//打印内部点初值
if(p%9==0)
cout<<endl;
p++;
}
}
cout<<"IT="<<IT;
break;
}
}
if(B==p)cout<<"不收敛"<<IT<<endl;
}

4-4．
#include<iostream>
usingnamespacestd;
voidmain()
{
intN,t,NP,TM,IT,M;
doubleBi,Fo,d,tf,a,h,y,dx,T0;
doubleT[11],TT[11];
cout<<"请输入壁厚d:";
cin>>d;
cout<<"请输入半壁厚等分数N：";
cin>>N;
cout<<"请输入控制打印节点节点温度的时间间隔NP:";
cin>>NP;
cout<<"请输入流体温度tf:";
cin>>tf;
cout<<"请输入时间间隔t:";
cin>>t;
cout<<"请输入终止计算时间TM:";
cin>>TM;
cout<<"请输入表面传热系数h:";
cin>>h;
cout<<"请输入平壁导热系数y:";
cin>>y;
cout<<"请输入初始温度T0:";
cin>>T0;
cout<<"请输入热扩散率a：";
cin>>a;
dx=d/(2.0*(double)N);
Bi=h*dx/y;
Fo=a*t/(dx*dx);
cout<<Bi<<''<<Fo<<endl;
if(Fo>1.0/(2.0*Bi+2.0))
cout<<"不稳定"<<endl;
else{
for(inti=0;i<11;i++)
T[i]=T0;
cout<<"0second:";
for(i=0;i<11;i++)
cout<<T[i]<<"";
cout<<endl;
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;i<11;i++)
T[i]=TT[i];
if(IT%NP==0)
{cout<<M<<"second：";
for(inti=0;i<=10;i++)
{
cout<<T