《核反应堆热工分析-课程设计》由会员分享，可在线阅读，更多相关《核反应堆热工分析-课程设计（35页珍藏版）》请在金锄头文库上搜索。
1、%另一种带入方式 %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;
2、 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_
3、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
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
5、*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+
6、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
7、*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=31
8、5）: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)*
9、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) %80
10、0 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 70
11、0to1=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.8
12、4) %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;el
13、seif(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
