傳熱學(xué)上機(jī)編程題答案

傳熱學(xué)上機(jī)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”請(qǐng)輸入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++)〃兩側(cè)賦初始值{T[j][O]=TBB;T[j][N-l]=TTB;)for(i=l;i<M-l;i++)〃內(nèi)部點(diǎn)賦初值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]<<,憂;〃打印內(nèi)部點(diǎn)初值

if(p%9==0)cout?endl;P++；))cout?,,IT=H?IT;break;))if(B==p)cout""不收斂"<<IT<<endl;)rI7-H:\Debug\4-2.exe”請(qǐng)輸入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<<"請(qǐng)輸入壁厚d:";cin?d;cout<<“請(qǐng)輸入半壁厚等分?jǐn)?shù)N：cin?N;cout<<”請(qǐng)輸入控制打印節(jié)點(diǎn)節(jié)點(diǎn)溫度的時(shí)間間隔NP:“;cin?NP;cout<<”請(qǐng)輸入流體溫度tf:M;cin?tf;cout"”請(qǐng)輸入時(shí)間間隔t:";cin?t;cout<<“請(qǐng)輸入終止計(jì)算時(shí)間TM:";cin?TM;cout<<“請(qǐng)輸入表面?zhèn)鳠嵯禂?shù)h:";cin?h;cout<<“請(qǐng)輸入平壁導(dǎo)熱系數(shù)y:";cin?y;cout"”請(qǐng)輸入初始溫度TO:";cin?T0;cout?”請(qǐng)輸入熱擴(kuò)散率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<<"不穩(wěn)定"<<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:厚打溫間露導(dǎo)溫散66厚壁*間止面壁始擴(kuò)26壁半控爵終表于黑03人A入入入入A-木入入0.等臼度隨最一瞽06分節(jié)t:呼本系T0a:反率7MM????Thy02^^數(shù):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