传热学上机编程题答案
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
iwv.tas »4.4St
£1.I9EK
宓泗$
53^860
ED
280
120.482
37_33?a ?0.9524
G2.0B7S
sc.ns
Ca.776
El .GG1
GO
ZOO
59
SB50
59
56
50
SB
56
J T =34Press any
key to cont Lnue
44
#in clude<iostream>
5®
50
SP
5U
120.47S
87-2327 70.9442
62.0795
56.T3H
53.7714
51.6589
SU
144吊帥
107«522 I4.47B9
70.4555
61.1949
56. sm
52.8661
血
1GB.7G*
112.794 A9.D11A
73 .2^46
的一如
E7.4927
則
144. m
T[i]=T0;
cout<< for(i=0;i<11;i++) cout<<T[i]<<;cout<<endl;
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];
潣瑵？请输入K,EPS,TI,TTB,TRB,TLB,TBB：v<endl; cin> >K>>EPS>>TI>>TTB»TRB>>TLB»TBB;
double T[9][7];
for(i=0;i<N;i++)//顶部底部赋初始值
{
T[0][i]=TLB;
T[M -1][i]=TRB;
}
cout«e ndl;
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++)
A=T[i][j] -TT;
if(eps<0) eps=-eps;
if(A<0)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]vv'\t：//
读书破万卷下笔如有神--
if(p%9==0)
cout«e ndl;
p++;
}
}
COUtv<IT=vvlT;
break;
}
}
晩？？?挠畯？尼不收敛vvlTvvendl;
}
■"H:■4 -2 .exe
倩输人
K.EPE.Tf.TTB.TRD.TLB.TBB：
iflH
0.01
20.09
re.00
50.09
bU.MkJ
59
sew
£0
传热学上机
高牥杰
2012/11/3
读书破万卷下笔如有神
4-2
#in clude<iostream>
using n amespace std;
void mai n()
{
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;
using n amespace std;
void mai n()
{
int N,t,NP,TM,IT,M;
double Bi,Fo,d,tf,a,h,y,dx,TO;
double T[11],TT[11];
潣瑵？请输入壁厚d:;
cin> >d;
潣瑵？请输入半壁厚等分数N:;
ci n»N;
潣瑵？请输入控制打印节点节点温度的时间间隔NP:;
ci源自文库>> NP;
读书破万卷下笔如有神
潣瑵？请输入流体温度tf:;
cin> >tf;
潣瑵？请输入时间间隔t:;
cin> >t;
潣瑵?请输入终止计算时间TM:; cin>>TM;
潣瑵?请输入表面传热系数h:; cin>>h;
潣瑵?请输入平壁导热系数y:; cin>>y;
潣瑵?请输入初始温度T0:;
cin>>T0;
{ 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]);
for(int i=0;i<=10;i++)
{cout<<T[i]<<;
}
cout<<endl;
}
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：;
潣瑵?请输入热扩散率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))潣瑵?不稳定<<endl;
else{
for(int i=0;i<11;i++)
£1.I9EK
宓泗$
53^860
ED
280
120.482
37_33?a ?0.9524
G2.0B7S
sc.ns
Ca.776
El .GG1
GO
ZOO
59
SB50
59
56
50
SB
56
J T =34Press any
key to cont Lnue
44
#in clude<iostream>
5®
50
SP
5U
120.47S
87-2327 70.9442
62.0795
56.T3H
53.7714
51.6589
SU
144吊帥
107«522 I4.47B9
70.4555
61.1949
56. sm
52.8661
血
1GB.7G*
112.794 A9.D11A
73 .2^46
的一如
E7.4927
則
144. m
T[i]=T0;
cout<< for(i=0;i<11;i++) cout<<T[i]<<;cout<<endl;
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];
潣瑵？请输入K,EPS,TI,TTB,TRB,TLB,TBB：v<endl; cin> >K>>EPS>>TI>>TTB»TRB>>TLB»TBB;
double T[9][7];
for(i=0;i<N;i++)//顶部底部赋初始值
{
T[0][i]=TLB;
T[M -1][i]=TRB;
}
cout«e ndl;
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++)
A=T[i][j] -TT;
if(eps<0) eps=-eps;
if(A<0)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]vv'\t：//
读书破万卷下笔如有神--
if(p%9==0)
cout«e ndl;
p++;
}
}
COUtv<IT=vvlT;
break;
}
}
晩？？?挠畯？尼不收敛vvlTvvendl;
}
■"H:■4 -2 .exe
倩输人
K.EPE.Tf.TTB.TRD.TLB.TBB：
iflH
0.01
20.09
re.00
50.09
bU.MkJ
59
sew
£0
传热学上机
高牥杰
2012/11/3
读书破万卷下笔如有神
4-2
#in clude<iostream>
using n amespace std;
void mai n()
{
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;
using n amespace std;
void mai n()
{
int N,t,NP,TM,IT,M;
double Bi,Fo,d,tf,a,h,y,dx,TO;
double T[11],TT[11];
潣瑵？请输入壁厚d:;
cin> >d;
潣瑵？请输入半壁厚等分数N:;
ci n»N;
潣瑵？请输入控制打印节点节点温度的时间间隔NP:;
ci源自文库>> NP;
读书破万卷下笔如有神
潣瑵？请输入流体温度tf:;
cin> >tf;
潣瑵？请输入时间间隔t:;
cin> >t;
潣瑵?请输入终止计算时间TM:; cin>>TM;
潣瑵?请输入表面传热系数h:; cin>>h;
潣瑵?请输入平壁导热系数y:; cin>>y;
潣瑵?请输入初始温度T0:;
cin>>T0;
{ 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]);
for(int i=0;i<=10;i++)
{cout<<T[i]<<;
}
cout<<endl;
}
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：;
潣瑵?请输入热扩散率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))潣瑵?不稳定<<endl;
else{
for(int i=0;i<11;i++)