核反应堆热工分析课程设计.doc

%另一种带入方式 %P=15.8;Nt=1.82*109;W=8916.66667; t_fin=287; L=3.66 ; m=121; n=265; n0=17; d_cs=0.0095 ; d_ci=0.00836 ; %有问题d_u=0.00819 ; %有问题s=0.0126 ;kexi=0.05;Fu=0.974; F_rn=1.35; F_zn=1.528; F_ln=1.11; F_h5e=0.95; F_qe=1.03; F_he=1.085; K_in=0.75; K_out=1.0;K_gr=1.05; f1=0.48; f2=1.02; f3=1.5; f4=1.56; f5=0.96; f6=0.48;hg=5678;%t_fout=input(请任意输入一堆芯出口温度（T=325）：t_fout=);t_fout=325;r=1;kk=0;while r=0.0001kk=kk+1;t_p=0.5*(t_fout+t_fin);if (280t_p=300) cp=19.35*t_p-350;elseif (300t_p=320)cp=34.43*t_p-4874;elseif (320t_p=340)cp=84.7*t_p-20960.4;endt1=t_fin+Fu*Nt/(W*(1-kexi)*cp);r=(t1-t_fout)/t_fout;t_fout=t1;endt_p=0.5*(t_fout+t_fin);% 计算燃料棒表面平均热流密度fz=m*n*pi*d_cs*L ; %计算堆心燃料棒的总传热面积q=Fu*Nt/fz ; %燃料棒表面平均热流密度F_qn=F_rn*F_zn*F_ln; % !有问题！q_max=q*F_qn*F_qe;ql=q*pi*d_cs ;ql_max=ql*F_qn*F_qe;%计算平均管的流速dert=0.0008;Af=m*n*(s2-pi*d_cs2/4)+4*n0*s*dert*m/2; %堆心内流通面积if(300t_p=310)v=4.25*10(-6)*t_p+1.0002*10(-4);elseif(310t_p=320)v=5.112*10(-6)*t_p-1.672*10(-4);elseif(320t_p=330)v=6.444*10(-6)*t_p-5.9344*10(-4);elseif(330t_p=0.0001kk1=kk1+1;t_fhp1=0.5*(t_fh1+t_fin);if (t_fhp1=300) cp1=19.35*t_fhp1-350;elseif (t_fhp1=320) cp1=34.43*t_fhp1-4874;elseif (t_fhp1=340) cp1=84.7*t_fhp1-20960.4;endt11=t_fin+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp1*6)*f1;r1=(t11-t_fh1)/t_fh1;t_fh1=t11;endt_fhp1=0.5*(t_fh1+t_fin);%fprintf(t_fh1=%.4fn,t_fh1);ts=346.19;%计算第一控制体出口处的包壳外壁温if(t_fhp1=300)u1=(-0.40072*t_fhp1+209.504)/106;k1=(-1.42504*t_fhp1+994.816)/1000;pr1=(-0.40072*t_fhp1+209.504)*(0.01935*t_fhp1-0.35)/(-1.42504*t_fhp1+994.816);elseif(t_fhp1=320)u1=(-1.044416*t_fhp1+402.6128)/106;k1=(-6.916*t_fhp1+2642.104)/1000;pr1=(-1.044416*t_fhp1+402.6128)*(0.03443*t_fhp1-4.874)/(-6.916*t_fhp1+2642.104);elseif(t_fhp1=340)u1=(-1.044416*t_fhp1+402.6128)/106;k1=(-6.916*t_fhp1+2642.104)/1000;pr1=(-1.044416*t_fhp1+402.6128)*(0.0847*t_fhp1-20.9604)/(-6.916*t_fhp1+2642.104);endre1=Wh*De/(Ab*u1) ;h1=0.023*re10.8*pr10.4*k1/De ;t_csh11=t_fh1+q*F_rn*f1*F_qe/h1;t_csh12=ts+25*(q*F_rn*f1*F_qe/106)0.25*exp(-P/6.2);if(t_csh11=t_csh12)t_csh1=t_csh12;end%fprintf(t_csh1=%.3fn,t_csh1); %t_csh1=t_csh1;%计算包壳内壁温度%t_cih1=input(请任意输入一包壳内壁温度（T=315）:t_cih1=);t_cih1=315;r12=1;kk11=0;while r12=0.0001kk11=kk11+1;t_cp1=0.5*(t_cih1+t_csh1);k12=0.00547*(1.8*t_cp1+32)+13.8;t12=t_csh1+ql*F_rn*F_qe*f1*log(d_cs/d_ci)/(2*pi*k12);r12=(t12-t_cih1)/t_cih1;t_cih1=t12;end%fprintf(t_cih1=%.4fn,t_cih1);%计算芯块表面温度t_uh1=t_cih1+ql*F_rn*F_qe*f1/(pi*0.5*(d_ci+d_u)*hg);%fprintf(t_uh1=%.4fn,t_uh1);%计算芯块中心温度if(t_uh1=400) %300 400ku1=0.051*t_uh1+6.02;elseif(t_uh1=500) %400 500ku1=0.0451*t_uh1+8.38;elseif(t_uh1=600) %500 600 ku1=0.0404*t_uh1+10.73elseif(t_uh1=700) %600 700ku1=0.0368*t_uh1+12.89;elseif(t_uh1=800) %700 800ku1=0.0337*t_uh1+15.06;elseif(t_uh1=900) %800 900 ku1=0.0312*t_uh1+17.06;elseif(t_uh1=1000) %900 1000ku1=0.0292*t_uh1+18.86;elseif(t_uh1=1100) %1000 1100 ku1=0.0255*t_uh1+22.56;elseif(t_uh1=1200) %1100 1200ku1=0.028*t_uh1+19.81;endkuto1=ku1+ql*F_rn*F_qe*f1/(4*pi*100) ;if(kuto1=34.93) %500 600to1=22.173*kuto1-185.8 ;elseif(kuto1=38.65) %600 700to1=26.74*kuto1-335.02;elseif(kuto1=42.02) %700 800to1=27.17*kuto1-446.88;elseif(kuto1=45.14) %800 900to1=32.05*kuto1-546.79;elseif(kuto1=48.06) %900 1000to1=34.25*kuto1-645.89;elseif(kuto1=50.61) %1000 1100to1=39.22*kuto1-884.71;elseif(kuto1=53.41)% 1100 1200to1=35.71*kuto1-707.5;elseif(kuto1=55.84) %1200 1298to1=40.32*kuto1-953.63;elseif(kuto1=58.4) %1298 1405to1=41.84*kuto1-1037;elseif(kuto1=61.95) %1405 1560to1=43.668*kuto1-1145;elseif(kuto1=66.87) %1560 1738to1=36.23*kuto1-682.25;elseif(kuto1=68.86) %1378 1876to1=69.35*kuto1-2899.2;elseif(kuto1=71.31) %1876 1990 to1=46.53*kuto1-1328.1;elseif(kuto1=74.88) %1990 2155to1=46.22*kuto1-1305.8;elseif(kuto1=0.0001kk2=kk2+1;t_fhp2=0.5*(t_fh2+t_fh1);if (t_fhp2=300) cp2=19.35*t_fhp2-350;elseif (t_fhp2=320) cp2=34.43*t_fhp2-4874;elseif (t_fhp2=340) cp2=84.7*t_fhp2-20960.4;endt22=t_fh1+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp2*6)*f2;r2=(t22-t_fh2)/t_fh2;t_fh2=t22;endt_fhp2=0.5*(t_fh2+t_fh1);%fprintf(t_fh2=%.4fn, t_fh2);%计算第二控制体出口处的包壳外壁温if(t_fhp2=300)u2=(-0.40072*t_fhp2+209.504)/106;k2=(-1.42504*t_fhp2+994.816)/1000;pr2=(-0.40072*t_fhp2+209.504)*(0.01935*t_fhp2-0.35)/(-1.42504*t_fhp2+994.816);elseif(t_fhp2=320)u2=(-1.044416*t_fhp2+402.6128)/106;k2=(-6.916*t_fhp2+2642.104)/1000;pr2=(-1.044416*t_fhp2+402.6128)*(0.03443*t_fhp2-4.874)/(-6.916*t_fhp2+2642.104);elseif(t_fhp2=340)u2=(-1.044416*t_fhp2+402.6128)/106;k2=(-6.916*t_fhp2+2642.104)/1000;pr2=(-1.044416*t_fhp2+402.6128)*(0.0847*t_fhp2-20.9604)/(-6.916*t_fhp2+2642.104);endre2=Wh*De/(Ab*u2);h2=0.023*re20.8*pr20.4*k2/De;t_csh21=t_fh2+q*F_rn*f2*F_qe/h2;t_csh22=ts+25*(q*F_rn*f2*F_qe/106)0.25*exp(-P/6.2);if(t_csh21=t_csh22)t_csh2=t_csh22;end%fprintf(t_csh2=%.4fn,t_csh2);%计算第二控制体包壳内壁温度%t_cih2=input(请任意输入第二控制体包壳内壁温度（例如T=350）:t_cih1=);t_cih2=350;r22=1;kk22=0;while r22=0.0001kk22=kk22+1;t_cp2=0.5*(t_cih2+t_csh2);k22=0.00547*(1.8*t_cp2+32)+13.8;t22=t_csh2+ql*F_rn*F_qe*f2*log(d_cs/d_ci)/(2*pi*k22);r22=(t22-t_cih2)/t_cih2;t_cih2=t22;end%fprintf(t_cih2=%.4fn,t_cih2);%计算芯块表面温度t_uh2=t_cih2+ql*F_rn*F_qe*f2/(pi*0.5*(d_ci+d_u)*hg) ;%fprintf(t_uh2=%.4fn,t_uh2);%计算芯块中心温度if(t_uh2=400) %300 400ku2=0.051*t_uh2+6.02;elseif(t_uh2=500) %400 500ku2=0.0451*t_uh2+8.38;elseif(t_uh2=600) %500 600 ku2=0.0404*t_uh2+10.73elseif(t_uh2=700) %600 700ku2=0.0368*t_uh2+12.89;elseif(t_uh2=800) %700 800ku2=0.0337*t_uh2+15.06;elseif(t_uh2=900) %800 900 ku2=0.0312*t_uh2+17.06;elseif(t_uh2=1000) %900 1000ku2=0.0292*t_uh2+18.86;elseif(t_uh2=1100) %1000 1100 ku2=0.0255*t_uh2+22.56;elseif(t_uh2=1200) %1100 1200ku2=0.028*t_uh2+19.81;endkuto2=ku2+ql*F_rn*F_qe*f2/(4*pi*100) ;if(kuto2=38.65) %600 700to2=26.74*kuto2-335.02;elseif(kuto2=42.02) %700 800to2=27.17*kuto2-446.88;elseif(kuto2=45.14) %800 900to2=32.05*kuto2-546.79;elseif(kuto2=48.06) %900 1000to2=34.25*kuto2-645.89;elseif(kuto2=50.61) %1000 1100to2=39.22*kuto2-884.71;elseif(kuto2=53.41)% 1100 1200to2=35.71*kuto2-707.5;elseif(kuto2=55.84) %1200 1298to2=40.32*kuto2-953.63;elseif(kuto2=58.4) %1298 1405to2=41.84*kuto2-1037;elseif(kuto2=61.95) %1405 1560to2=43.668*kuto2-1145;elseif(kuto2=66.87) %1560 1738to2=36.23*kuto2-682.25;elseif(kuto2=68.86) %1378 1876to2=69.35*kuto2-2899.2;elseif(kuto2=71.31) %1876 1990 to2=46.53*kuto2-1328.1;elseif(kuto2=74.88) %1990 2155to2=46.22*kuto2-1305.8;elseif(kuto2=0.0001kk3=kk3+1;t_fhp3=0.5*(t_fh3+t_fh2);if (t_fhp3=300) cp3=19.35*t_fhp3-350;elseif (t_fhp3=320) cp3=34.43*t_fhp3-4874;elseif (t_fhp3=340) cp3=84.7*t_fhp3-20960.4;endt33=t_fh2+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp3*6)*f3;r3=(t33-t_fh3)/t_fh3;t_fh3=t33;endt_fhp3=0.5*(t_fh3+t_fh2);%fprintf(t_fh3=%.3fn, t_fh3);%计算第三控制体出口处的包壳外壁温if(t_fhp3=300)u3=(-0.40072*t_fhp3+209.504)/106;k3=(-1.42504*t_fhp3+994.816)/1000;pr3=(-0.40072*t_fhp3+209.504)*(0.01935*t_fhp3-0.35)/(-1.42504*t_fhp3+994.816);elseif(t_fhp3=320)u3=(-1.044416*t_fhp3+402.6128)/106;k3=(-6.916*t_fhp3+2642.104)/1000;pr3=(-1.044416*t_fhp3+402.6128)*(0.03443*t_fhp3-4.874)/(-6.916*t_fhp3+2642.104);elseif(t_fhp3=340)u3=(-1.044416*t_fhp3+402.6128)/106;k3=(-6.916*t_fhp3+2642.104)/1000;pr3=(-1.044416*t_fhp3+402.6128)*(0.0847*t_fhp3-20.9604)/(-6.916*t_fhp3+2642.104);endre3=Wh*De/(Ab*u3);h3=0.023*re30.8*pr30.4*k3/De;t_csh31=t_fh3+q*F_rn*f3*F_qe/h3;t_csh32=ts+25*(q*F_rn*f3*F_qe/106)0.25*exp(-P/6.2);if(t_csh31=t_csh32)t_csh3=t_csh32;end%fprintf(t_csh3=%.3fn,t_csh3);%计算第三控制体包壳内壁温度%t_cih3=input(请任意输入第三控制体包壳内壁温度（例如T=385）:t_cih3=);t_cih3=385;r33=1;kk33=0;while r33=0.0001kk33=kk33+1;t_cp3=0.5*(t_cih3+t_csh3);k33=0.00547*(1.8*t_cp3+32)+13.8;t33=t_csh3+ql*F_rn*F_qe*f3*log(d_cs/d_ci)/(2*pi*k33);r33=(t33-t_cih3)/t_cih3;t_cih3=t33;end%fprintf(t_cih3=%.3fn,t_cih3);%计算芯块表面温度t_uh3=t_cih3+ql*F_rn*F_qe*f3/(pi*0.5*(d_ci+d_u)*hg);%fprintf(t_uh3=%.3fn,t_uh3);%计算芯块中心温度if(t_uh3=400) %300 400ku3=0.051*t_uh3+6.02;elseif(t_uh3=500) %400 500ku3=0.0451*t_uh3+8.38;elseif(t_uh3=600) %500 600 ku3=0.0404*t_uh3+10.73 ;elseif(t_uh3=700) %600 700ku3=0.0368*t_uh3+12.89;elseif(t_uh3=800) %700 800ku3=0.0337*t_uh3+15.06;elseif(t_uh3=900) %800 900 ku3=0.0312*t_uh3+17.06;elseif(t_uh3=1000) %900 1000ku3=0.0292*t_uh3+18.86;elseif(t_uh3=1100) %1000 1100 ku3=0.0255*t_uh3+22.56;elseif(t_uh3=1200) %1100 1200ku3=0.028*t_uh3+19.81;endkuto3=ku3+ql*F_rn*F_qe*f3/(4*pi*100) ;if kuto3=38.65 %600 700to3=26.74*kuto3-335.02;elseif kuto3=42.02 %700 800to3=27.17*kuto3-446.88;elseif(kuto3=45.14) %800 900to3=32.05*kuto3-546.79;elseif kuto3=48.06 %900 1000to3=34.25*kuto3-645.89;elseif kuto3=50.61 %1000 1100to3=39.22*kuto3-884.71;elseif kuto3=53.41 % 1100 1200to3=35.71*kuto3-707.5;elseif kuto3=55.84 %1200 1298to3=40.32*kuto3-953.63;elseif kuto3=58.4 %1298 1405to3=41.84*kuto3-1037;elseif kuto3=61.95 %1405 1560to3=43.668*kuto3-1145;elseif kuto3=66.87 %1560 1738to3=36.23*kuto3-682.25;elseif kuto3=68.86 %1378 1876to3=69.35*kuto3-2899.2;elseif kuto3=71.31 %1876 1990 to3=46.53*kuto3-1328.1;elseif kuto3=74.88 %1990 2155to3=46.22*kuto3-1305.8;elseif kuto3=0.0001kk4=kk4+1;t_fhp4=0.5*(t_fh4+t_fh3);if t_fhp4=300 cp4=19.35*t_fhp4-350;elseif t_fhp4=320 cp4=34.43*t_fhp4-4874;elseif t_fhp4=340 cp4=84.7*t_fhp4-20960.4;endt44=t_fh3+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp4*6)*f4;r4=(t44-t_fh4)/t_fh4;t_fh4=t44;endt_fhp4=0.5*(t_fh4+t_fh3);%fprintf(t_fh4=%.3fn, t_fh4);ts=346.19;%计算第四控制体出口处的包壳外壁温if(t_fhp4=300)u4=(-0.40072*t_fhp4+209.504)/106;k4=(-1.42504*t_fhp4+994.816)/1000;pr4=(-0.40072*t_fhp4+209.504)*(0.01935*t_fhp4-0.35)/(-1.42504*t_fhp4+994.816);elseif(t_fhp4=320)u4=(-1.044416*t_fhp4+402.6128)/106;k4=(-6.916*t_fhp4+2642.104)/1000;pr4=(-1.044416*t_fhp4+402.6128)*(0.03443*t_fhp4-4.874)/(-6.916*t_fhp4+2642.104);elseif(t_fhp4=340)u4=(-1.044416*t_fhp4+402.6128)/106;k4=(-6.916*t_fhp4+2642.104)/1000;pr4=(-1.044416*t_fhp4+402.6128)*(0.0847*t_fhp4-20.9604)/(-6.916*t_fhp4+2642.104);endre4=Wh*De/(Ab*u4);h4=0.023*re40.8*pr40.4*k4/De;t_csh41=t_fh4+q*F_rn*f4*F_qe/h4;t_csh42=ts+25*(q*F_rn*f4*F_qe/106)0.25*exp(-P/6.2);if(t_csh41=t_csh42)t_csh4=t_csh42;end%fprintf(t_csh4=%.3fn,t_csh4);%计算第四控制体包壳内壁温度%t_cih4=input(请任意输入第五控制体包壳内壁温度（T=386）:t_cih4=);t_cih4=386;r44=1;kk44=0;while r44=0.0001kk44=kk44+1;t_cp4=0.5*(t_cih4+t_csh4);k44=0.00547*(1.8*t_cp4+32)+13.8;t44=t_csh4+ql*F_rn*F_qe*f4*log(d_cs/d_ci)/(2*pi*k44);r44=(t44-t_cih4)/t_cih4;t_cih4=t44;end%fprintf(t_cih4=%.4fn,t_cih4);%计算芯块表面温度t_uh4=t_cih4+ql*F_rn*F_qe*f4/(pi*0.5*(d_ci+d_u)*hg);%fprintf(t_uh4=%.3fn,t_uh4);%计算芯块中心温度if t_uh4=400 %300 400ku4=0.051*t_uh4+6.02;elseif t_uh4=500 %400 500ku4=0.0451*t_uh4+8.38;elseif t_uh4=600 %500 600 ku4=0.0404*t_uh4+10.73elseif t_uh4=700 %600 700ku4=0.0368*t_uh4+12.89;elseif t_uh4=800 %700 800ku4=0.0337*t_uh4+15.06;elseif t_uh4=900 %800 900 ku4=0.0312*t_uh4+17.06;elseif t_uh4=1000 %900 1000ku4=0.0292*t_uh4+18.86;elseif t_uh4=1100 %1000 1100 ku4=0.0255*t_uh4+22.56;elseif t_uh4=1200 %1100 1200ku4=0.028*t_uh4+19.81;endkuto4=ku4+ql*F_rn*F_qe*f4/(4*pi*100);if kuto4=38.65 %600 700to4=26.74*kuto4-335.02;elseif kuto4=42.02 %700 800to4=27.17*kuto4-446.88;elseif kuto4=45.14 %800 900to4=32.05*kuto4-546.79;elseif kuto4=48.06 %900 1000to4=34.25*kuto4-645.89;elseif kuto4=50.61 %1000 1100to4=39.22*kuto4-884.71;elseif kuto4=53.41 % 1100 1200to4=35.71*kuto4-707.5;elseif kuto4=55.84 %1200 1298to4=40.32*kuto4-953.63;elseif kuto4=58.4 %1298 1405to4=41.84*kuto4-1037;elseif kuto4=61.95 %1405 1560to4=43.668*kuto4-1145;elseif kuto4=66.87 %1560 1738to4=36.23*kuto4-682.25;elseif kuto4=68.86 %1738 1876to4=69.35*kuto4-2899.2;elseif kuto4=71.31 %1876 1990 to4=46.53*kuto4-1328.1;elseif kuto4=74.88 %1990 2155to4=46.22*kuto4-1305.8;elseif kuto4=0.0001kk5=kk5+1;t_fhp5=0.5*(t_fh5+t_fh4);if t_fhp5=300 cp5=19.35*t_fhp5-350;elseif t_fhp5=320 cp5=34.43*t_fhp5-4874;elseif t_fhp5=340 cp5=84.7*t_fhp5-20960.4;endt55=t_fh4+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp5*6)*f5;r5=(t55-t_fh5)/t_fh5;t_fh5=t55;endt_fhp5=0.5*(t_fh5+t_fh5);%fprintf(t_fh5=%.4fn, t_fh5);%计算第五控制体出口处的包壳外壁温if(t_fhp5=300)u5=(-0.40072*t_fhp5+209.504)/106;k5=(-1.42504*t_fhp5+994.816)/1000;pr5=(-0.40072*t_fhp5+209.504)*(0.01935*t_fhp5-0.35)/(-1.42504*t_fhp5+994.816);elseif(t_fhp5=320)u5=(-1.044416*t_fhp5+402.6128)/106;k5=(-6.916*t_fhp5+2642.104)/1000;pr5=(-1.044416*t_fhp5+402.6128)*(0.03443*t_fhp5-4.874)/(-6.916*t_fhp5+2642.104);elseif(t_fhp5=340)u5=(-1.044416*t_fhp5+402.6128)/106;k5=(-6.916*t_fhp5+2642.104)/1000;pr5=(-1.044416*t_fhp5+402.6128)*(0.0847*t_fhp5-20.9604)/(-6.916*t_fhp5+2642.104);endre5=Wh*De/(Ab*u5);h5=0.023*re50.8*pr50.4*k5/De;t_csh51=t_fh5+q*F_rn*f5*F_qe/h5;t_csh52=ts+25*(q*F_rn*f5*F_qe/106)0.25*exp(-P/6.2);if(t_csh51=t_csh52)t_csh5=t_csh52;end%fprintf(t_csh5=%.4fn,t_csh5); %计算第五控制体包壳内壁温度%t_cih5=input(请任意输入第五控制体包壳内壁温度（T=371）:t_cih5=);t_cih5=371;r55=1;kk55=0;while r55=0.00001kk55=kk55+1;t_cp5=0.5*(t_cih5+t_csh5);k55=0.00547*(1.8*t_cp5+32)+13.8;t55=t_csh5+ql*F_rn*F_qe*f5*log(d_cs/d_ci)/(2*pi*k55);r55=(t55-t_cih5)/t_cih5;t_cih5=t55;end%fprintf(t_cih5=%.4fn,t_cih5);%计算芯块表面温度t_uh5=t_cih5+ql*F_rn*F_qe*f5/(pi*0.5*(d_ci+d_u)*hg);%fprintf(t_uh5=%.4fn,t_uh5);%计算芯块中心温度if t_uh5=400 %300 400ku5=0.051*t_uh5+6.02;elseif t_uh5=500 %400 500ku5=0.0451*t_uh5+8.38;elseif t_uh5=600 %500 600 ku5=0.0404*t_uh5+10.73;elseif t_uh5=700 %600 700ku5=0.0368*t_uh5+12.89;elseif t_uh5=800 %700 800ku5=0.0337*t_uh5+15.06;elseif t_uh5=900 %800 900 ku5=0.0312*t_uh5+17.06;elseif t_uh5=1000 %900 1000ku5=0.0292*t_uh5+18.86;elseif t_uh5=1100 %1000 1100 ku5=0.0255*t_uh5+22.56;elseif t_uh5=1200 %1100 1200ku5=0.028*t_uh5+19.81;endkuto5=ku5+ql*F_rn*F_qe*f5/(4*pi*100);if kuto5=38.65 %600 700to5=26.74*kuto5-335.02;elseif kuto5=42.02 %700 800to5=27.17*kuto5-446.88;elseif kuto5=45.14 %800 900to5=32.05*kuto5-546.79;elseif kuto5=48.06 %900 1000to5=34.25*kuto5-645.89;elseif kuto5=50.61 %1000 1100to5=39.22*kuto5-884.71;elseif kuto5=53.41 % 1100 1200to5=35.71*kuto5-707.5;elseif kuto5=55.84 %1200 1298to5=40.32*kuto5-953.63;elseif kuto5=58.4 %1298 1405to5=41.84*kuto5-1037;elseif kuto5=61.95 %1405 1560to5=43.668*kuto5-1145;elseif kuto5=66.87 %1560 1738to5=36.23*kuto5-682.25;elseif kuto5=68.86 %1738 1876to5=69.35*kuto5-2899.2;elseif kuto5=71.31 %1876 1990 to5=46.53*kuto5-1328.1;elseif kuto5=74.88 %1990 2155to5=46.22*kuto5-1305.8;elseif kuto5=0.0001kk6=kk6+1;t_fhp6=0.5*(t_fh6+t_fh5);if t_fhp6=300 cp6=19.35*t_fhp6-350;elseif t_fhp6=320 cp6=34.43*t_fhp6-4874;elseif t_fhp6=340 cp6=84.7*t_fhp6-20960.4;endt66=t_fh5+q*F_rn*F_he*F_h5e*pi*L*d_cs/(Wh*cp6*6)*f6;r6=(t66-t_fh6)/t_fh6;t_fh6=t66;endt_fhp6=0.5*(t_fh6+t_fh6);%fprintf(t_fh6=%.4fn, t_fh6);%计算第六控制体出口处的包壳外壁温if(t_fhp6=300)u6=(-0.40072*t_fhp6+209.504)/106;k6=(-1.42504*t_fhp6+994.816)/1000;pr6=(-0.40072*t_fhp6+209.504)*(0.01935*t_fhp6-0.35)/(-1.42504*t_fhp6+994.816);elseif(t_fhp6=320)u6=(-1.044416*t_fhp6+402.6128)/106;k6=(-6.916*t_fhp6+2642.104)/1000;pr6=(-1.044416*t_fhp6+402.6128)*(0.03443*t_fhp6-4.874)/(-6.916*t_fhp6+2642.104);elseif(t_fhp6=340)u6=(-1.044416*t_fhp6+402.6128)/106;k6=(-6.916*t_fhp6+2642.104)/1000;pr6=(-1.044416*t_fhp6+402.6128)*(0.0847*t_fhp6-20.9604)/(-6.916*t_fhp6+2642.104);endre6=Wh*De/(Ab*u6);h6=0.023*re60.8*pr60.4*k6/De;t_csh61=t_fh6+q*F_rn*f6*F_qe/h6;t_csh62=ts+25*(q*F_rn*f6*F_qe/106)0.25*exp(-P/6.2);if(t_csh61=t_csh62)t_csh6=t_csh62;end