传热学上机编程题答案

传热学上机2012/11/34-2.#include<iostream>usingnamespacestd;voidmain()(inti=0,j=0/M=9/N=7/K/n；B;intp=l;doubleEPS^BJRBJLBJBB^epsJLA;coukv”请输入K/EPS/TI,TTB/TRBzTLB/TBB:,I«endl;cin»K»EPS»TI»TTB»TRB»TLB»TBB;doubleT[9][7];for(i=0;i<N;i++)〃顶部底部赋初始值{T[O][i]=TLB;T[M-l][i]=TRB;}cout«endl;for(j=l;j<M;j++)〃两侧赋初始值{T[j][O]=TBB;T[j][N-l]=TTB;)for(i=l;i<M-l;i++)〃内部点赋初值Tlfor(j=l;j<N-l;j++)T[i]Ul=TI;for(j=N-l;j>=0;j-)B=p；for(IT=0;IT<K;IT++){eps=T[l][N-2]-0.25*(T[0][N-2]+T[l][N-3]+T[2][N-2]+T[l][N-l]);for(i=l;i<M-l;i++){for(j=N-2;j>=l;j-){TT=T[i][j];T[i]0]=O.25*(T[i-l]U]+T[i]U-l]+T[i+l]U]+T[i]U+l])；A=T[i]D]-TT;if(eps<0)eps=-eps;if(A<0)A=-A;if(eps<=A)eps=A;))if(eps<=EPS){for(j=N-l;j>=0;j-){for(i=0;i<M;i++){cout<<T[i][j]<<,忧;〃打印内部点初值

if(p%9==0)cout«endl;P++；))cout«,,IT=H«IT;break;))if(B==p)cout""不收敛"<<IT<<endl;)rI7-H:\Debug\4-2.exe”请输入k，eps,ti,ttb，trb，tlb，tbb：1000.0120.0050.0050.00200.0050.00144.6891R7.S22R4.4789150.769113.70489.0118120.48287.239870.9524144.6891R7.S22R4.4789150.769113.70489.0118120.48287.239870.95245050505062.079556.931453.771451.65887A.455561.8849S6.5R252.866173.394663.66357.492753.306270.462461.890656.505852.86862.087556.93853.77651.661IT=34Pressanykeytocontinue4-4.#include<iostream>usingnamespacestd;voidmain()(intNXNRTMJTM;doubleBijFo,d,tf,a,h,ydx,TO;doubleTtllLTTlll];cout<<"请输入壁厚d:";cin»d;cout<<“请输入半壁厚等分数N：cin»N;cout<<”请输入控制打印节点节点温度的时间间隔NP:“;cin»NP;cout<<”请输入流体温度tf:M;cin»tf;cout"”请输入时间间隔t:";cin»t;cout<<“请输入终止计算时间TM:";cin»TM;cout<<“请输入表面传热系数h:";cin»h;cout<<“请输入平壁导热系数y:";cin»y;cout"”请输入初始温度TO:";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<ll;i++)T[i]=T0;cout«M0second:";for(i=0;i<ll;i++)cout«T[i]«"cout«endl;for(IT=l;IT*t<=TM;IT++){M=IT*t;forfinti=9;i>=l;i-)(TT[i]=Fo*(T[i-l]+T[i+l])+(1.0-2.0*Fo)*T[i];)TT[N]=2.0*Fo*(T[N-l]+tf*Bi)+(1.0-2.0*Fo-2.0*Fo*Bi)*T[N];TT[0]=2.0*Fo*T[l]+(1.0-2.0*Fo)*T[0];for(i=0;i<ll;i++)T[i]=TT[i];if(IT%NP==O){cout«M«nsecond：M;for(inti=0;i<=10;i++)(cout«T[i]«H)cout«endl;

)}Fo=a*TM/((d/2.0)*(d/2.0));。回！汉cout«Bi«',«Fo«endl;。回！汉a"C:\Users\Administfator\Desktop\Debug\4-4.exe'474747475I.474747575U3

主皇皇皇皇皇皇皇皇皇0E3.、Tk、Tk、Tk、Tk、Tk474747475I.474747575U3

主皇皇皇皇皇皇皇皇皇0E3.、Tk、Tk、Tk、Tk、Tk、Tk、Tk、Tk、Tk、Tkd:厚打温间露导温散66厚壁*间止面壁始扩26壁半控爵终表于黑03人A入入入入A-木入入0.等臼度随最一瞽06分节t:呼本系T0a:反率7MM••••Thy02^^数:2124:212243648600172second：second：second：second：second：second：।202020202020202020202020202020202020202020202020।2020202020202020202020.001220.046121.022332.47772020.000120.001520.028120.372423.423840.368520.000120.000920.011920.127121.037226.218445.837320.000420.004520.04320.328421.973229.01349.941720.000220.001720.014620.106720.641223.095131.66453.2220.000120.000620.00520.035220.211721.059724.331434.13155.96035.960384second：58.325596second：52760.4139202020.000220.001720.011820.070720.363121.570925.629836.41392020.000120.000620.00420.023820.124720.562522.159326.954138.20201R8second,2B.0RR120.000220.001220.008120.042920.22R,8R8922.8095