?試設(shè)計(jì)某實(shí)驗(yàn)平臺(tái)所需要使用的一臺(tái)氣水換熱器。

設(shè)計(jì)校核計(jì)算通過(guò)matlab實(shí)現(xiàn)，請(qǐng)求斧正！

clear,clc %申明所有壓力單位統(tǒng)一為Pa %% 預(yù)備計(jì)算 %%預(yù)備計(jì)算 %%假設(shè)換熱器效率n=0.8，任務(wù)書(shū)中要求的是大于等于0.7，x假設(shè)為為冷測(cè)水的進(jìn)口流量 %初始值設(shè)定 t1_1=90;t2_1=25;n=0.8 x=890; t1_11 = t1_1 - n*(t1_1 - t2_1);%%計(jì)算出口氣溫 %%由查表得出的空氣的物性參數(shù) cp1 = 1003+0.02*43+4*0.00001*43^2; %%由查表得出的水的物性參數(shù) cp2 = 4184.4-0.6964*21+1.06964*21+1.036*10^(-2)*21^2; qm1=(1000+20*35)/3600;%36為學(xué)號(hào) qm2=x/3600; W1=qm1*cp1; W2=qm2*cp2; C=W1/W2;%計(jì)算C %由冷測(cè)進(jìn)口流量和C*的關(guān)系可判斷，C*在x=1600以內(nèi)是永遠(yuǎn)大于0.5的 t2_11=t2_1+W1/W2*(t1_1-t1_11); %Wmax側(cè)（冷側(cè)）取算數(shù)平均溫度 tm2=(t2_1+t2_11)/2 %兩流體的平均溫差為 deta_t1m=((t1_1-tm2)-(t1_11-tm2))/log((t1_1-tm2)/(t1_11-tm2)); %熱側(cè)平均溫度為 tm1=tm2+deta_t1m; %根據(jù)平均溫度求物性參數(shù) mu1=1.50619*10^(-6)*(tm1+273)^1.5;%空氣的動(dòng)力黏度miu lam1=2.456*10^(-4)*(tm1+273)^0.823;%空氣導(dǎo)熱系數(shù)λ1 Pr1=mu1*cp1/lam1;%計(jì)算熱側(cè)的普朗特?cái)?shù)Pr1 mu2=10^(230.298/(tm2+126.203)-4.5668);%水的動(dòng)力黏度 lam2=0.598+1.373*10^(-3)*tm2-5.333*10^(-6)*(tm2^2);%水的導(dǎo)熱系數(shù) Pr2=mu2*cp2/lam2%計(jì)算冷側(cè)的普朗特?cái)?shù)Pr2 mixsum=[x;mu1;lam1;Pr1;mu2;lam2;Pr2]'; %% 產(chǎn)品結(jié)構(gòu)規(guī)劃及其計(jì)算 %參數(shù)說(shuō)明 冷端 熱端 %間斷長(zhǎng)度 ls1 ls2 %翅片間距 Pf1 Pf2 %板間距 s1 s2 %翅片厚度 Deltaf1 Deltaf2 %翅片高度 h1 h2 %翅片層數(shù) m1 m2 %翅片厚度 Deltap %側(cè)板厚度 Deltas %封條寬度 bs1 bs2 %熱側(cè)進(jìn)口寬度 l2 %冷側(cè)出口 l1 %熱側(cè)出口 l3ls1=240 ; ls2=60 ; %間斷長(zhǎng)度 Pf1=50 ; Pf2=40 %翅片間距s1= 150 ; s2= 50 ; %板間距 Deltaf1=5.5 ; Deltaf2=3.35 ; % 翅片厚度 h1=75.5 ; h2=65 ; %翅片高度m1=370 ; m2=400 ; %翅片層數(shù) Deltap=13 ; % Deltap2= %翅片厚度 Deltas=40 ; % Deltas2= % 側(cè)板厚度 bs1=150 ; bs2= 90 ; %封條寬度 l2= 1000 % 熱側(cè)進(jìn)口寬度 l1= 1000 % 冷側(cè)出口%l3= %熱側(cè)出口 l3 = m1*s1+m2*s2+2*m2*Deltap; %計(jì)算非流動(dòng)方向長(zhǎng)度 de1=2*(Pf1-2*Deltaf1)*(s1-Deltaf1)/((Pf1-Deltaf1)+(s1-Deltaf1)+(s1-Deltaf1)*Deltaf1/ls1); de2=2*(Pf2-2*Deltaf2)*(s2-Deltaf2)/((Pf2-Deltaf2)+(s2-Deltaf2)+(s2-Deltaf2)*Deltaf2/ls2); %分別計(jì)算當(dāng)量直徑 phi1=(s1-Deltaf1)/(Pf1+2*Deltaf1+s1); phi2=(s2-Deltaf2)/(Pf2+2*Deltaf2+s2);%計(jì)算翅片面積比 Ay1=(l2-2*bs1)*(l3-2*Deltas); Ay2=(l1-1.5*bs2)*(l3-2*Deltas);%計(jì)算迎風(fēng)面積 Vp1=(l1-1.5*bs2)*(l2-2*bs1)*m1*s1; Vp2=(l1-1.5*bs2)*(l2-2*bs1)*m2*s2;%計(jì)算板間體積 %傳熱面積密度為 %以各板間體積為基準(zhǔn) A1=2*((Pf1-Deltaf1)+(s1-Deltaf1)); A2=2*((Pf2-Deltaf2)+(s2-Deltaf2)); beta1=A1/Vp1; beta2=A2/Vp2; %總傳熱面積為 A1=beta1*Vp1;%這里的單位是平方毫米 A2=beta2*Vp1; %最小流通面積為 Ac1=de1*A1/(4*l1); Ac2=de2*A2/(4*l2); %每側(cè)孔度為 sig1=Ac1/Ay1; sig2=Ac2/Ay2; Ap=(2*m1+2)*(l1-1.5*bs2)*(l2-2*bs1);%一次傳熱面積（隔板導(dǎo)熱面積）(㎡) %% 計(jì)算質(zhì)量流速、雷諾數(shù)、對(duì)流表面?zhèn)鳠嵯禂?shù)、翅片效率和表面效率 %質(zhì)量流速為 gm1=qm1/Ac1; gm2=qm2/Ac2; %計(jì)算雷諾數(shù)Re Re1=gm1*de1/mu1; Re2=gm2*de2/mu2; %對(duì)于空氣或氣體工質(zhì)，層流區(qū)Re<=1000時(shí)，適用Weiting擬合關(guān)系式（3-26） %溫流區(qū)Re>=2000時(shí)，適用Weiting擬合關(guān)系式（3-27） %對(duì)于過(guò)渡區(qū)需要根據(jù)下式確定Re^*_j, %%%%%%%%%%！！！！！！！！！！！！ %設(shè)置判斷語(yǔ)句 %判斷Re1,根據(jù)雷諾數(shù)爬滿斷面熱側(cè)空氣的傳熱系數(shù)j1 if Re1<=1000j1=0.483*(ls1/de1)^0.162 * ((Pf1-Deltaf1)/(s1-Deltaf1))^(-0.184)*Re1^(-0.536); else if Re1>=2000j1=0.242*(ls1/de1)^(-0.322)*(Deltaf1/de1)^0.089*Re1^(-0.368);elseRe_j=61.9*((ls1/de1)^0.952)*((Pf1-Deltaf1)/(s1-Deltaf1))^(-1.1)*Re1^(-0.53);if Re1>=Re_jj1=0.242*(ls1/de1)^(-0.322)*((Pf1-Deltaf1)/(s1-Deltaf1))^(0.089)*Re1^(-0.368);else j1=0.483*(ls1/de1)^0.384 * ((Pf1-Deltaf1)/(s1-Deltaf1))^(-0.184)*Re1^(-0.368);endend end %Re_j=61.9*((ls1/de1)^0.952)*((pf1-Deltaf1)/(s1-Deltaf1))^(-1.1)*Re^(-0.53); %根據(jù)題意，Re1>Re_j，因此熱側(cè)空氣傳熱系數(shù)可用Weiting擬合公式（3-26）確定j %否則用(3-27)確定 %j1=0.242*(ls1/de1)^(-0.322)*((pf1-Deltaf1)/(s1-Deltaf1))^(0.089)*Re^(-0.368); %對(duì)于冷測(cè)的傳熱系數(shù)采用下式計(jì)算 Pr=2;%Pr暫定為2，需核實(shí)，！！！！！！！！！！！！！！！ j2=0.287*Re1^(-0.322)*Pr^(0.167);%%這里的Pr不知道取冷端還是熱端。還是平均數(shù)；!!!! %對(duì)流傳熱面系數(shù) alp1=j1*gm1*cp1/Pr1^(2/3); alp2=j2*gm2*cp2/Pr2^(2/3); %翅片參數(shù) m1=sqrt(2*alp1*(1+Deltaf1/ls1)/(lam1*Deltaf1)); m2=sqrt(2*alp2*(1+Deltaf2/ls2)/(lam2*Deltaf2)); nf1=tanh(m1*h1)/(m1*h1); nf2=tanh(m2*h2)/(m2*h2); %兩側(cè)翅片表面效率為 n01=1-phi1*(1-nf1);%%書(shū)上例題公式有Af1，不知其所指？？？？？已解決，使用等式后邊的公式進(jìn)行計(jì)算 n02=1-phi2*(1-nf2); %% 計(jì)算壁面熱阻和傳熱系數(shù) lamw=169; Rew=Deltap/(lamw*Ap);%這里的λw是指什么 %對(duì)于氣-水換熱器，暫不考慮污垢熱阻，故有 KA1=1/(n01*alp1*A1)+Rew+1/(n02*alp2*A2);%KA1是指KA的倒數(shù) KA=1/KA1; %因?yàn)镵A=K1A1=K2A2,當(dāng)以熱流體側(cè)的總傳熱面積A1為基準(zhǔn)時(shí)， %對(duì)應(yīng)的換熱系數(shù)為 K1=KA/A1;%% 計(jì)算傳熱單元數(shù)、換熱器效率和出口流體溫度 W1=qm1*cp1; W2=qm2*cp2; C=W1/W2;%C*=Wmin/Wmax; NTU=KA/W1;%Wmin=W1; %兩流體各自非混流的叉流換熱器效率，此處按德雷克近似關(guān)系式計(jì)算單個(gè)芯體的效率，有 ni=1-exp(((NTU^0.22)/C)*exp((-C*NTU^0.78)-1)); %換熱器總效率為 n1 n1=(((1-C*ni)/(1-ni))^2-1/(((1-C*ni)/(1-ni))^2-C));%n1為設(shè)計(jì)換熱效率。n為假設(shè)效率，第一次計(jì)算n=0。74 %此處需要判斷是否符合當(dāng)初假定的效率！！！！！！ %換熱器的傳熱熱流量為 Phi=n1*W1*(t1_1-t2_1); %流體出口溫度為 t1_11=t1_1-Phi/W1; t2_11=t2_1-Phi/W2; %與假設(shè)相對(duì)誤差為 detan=abs((n1-n)/n); %補(bǔ)充數(shù)組換熱量與效率的關(guān)系。mixsum2=[]; %% 計(jì)算阻力 detap1=100;%預(yù)設(shè)的detap1 %比體積計(jì)算 p1_1=200*1000; p2_1=100*1000;%任務(wù)書(shū)所給的初識(shí)進(jìn)口壓力； p1_11=p1_1-detap1; p2_11=p2_1-detap1; R = 287;%是一個(gè)常數(shù)，定值 v1_1=R*t1_1/p1_1; v1_11=R*t1_11/p1_11; vm1=(v1_1+v1_11)/2;%計(jì)算熱側(cè)的參數(shù)，均為國(guó)際單位。 v2_1=R*t2_1/p2_1; v2_11=R*t2_11/p2_11; vm2=(v2_1+v2_11)/2;%計(jì)算冷測(cè)的參數(shù)； %計(jì)算熱側(cè)阻力 %由于書(shū)中圖3-12中查詢K_1和K_11較繁瑣，故將其曲線用線性關(guān)系近似代替 K1_1=(-3/7)*sig1+1.342857143; K1_11=(-13/8)*sig1+0.825; K2_1=(-3/7)*sig2+1.342857143; K2_11=(-13/8)*sig2+0.825; %設(shè)置判斷語(yǔ)句 %判斷Re1,根據(jù)雷諾數(shù)爬滿斷面熱側(cè)空氣的f if Re1<=1000f1=7.661*(ls1/de1)^0.384 * ((Pf1-Deltaf1)/(s1-Deltaf1))^(-0.092)*Re1^(-0.712); else if Re1>=2000f1=1.136*(ls1/de1)^(-0.781)*(Deltaf1/de1)^0.534*Re1^(-0.198);elseRe_f=41*((ls1/de1)^0.772)*((Pf1-Deltaf1)/(s1-Deltaf1))^(-0.179)*Re1^(-1.04);if Re1>=Re_jf1=7.661*(ls1/de1)^0.384 * ((Pf1-Deltaf1)/(s1-Deltaf1))^(-0.092)*Re1^(-0.712);else f1=1.136*(ls1/de1)^(-0.781)*(Deltaf1/de1)^0.534*Re1^(-0.198);endend end %此時(shí)已經(jīng)求得f1；計(jì)算實(shí)際的detap1，實(shí)際的命名為detap1_1; %A_1=1-sig1^2+K1_1; %B_1=2*(v1_11/v1_1-1); %C_1=(4*f1*2*l1*vm1)/(de1*v1_1); %D_1=(1-sig1^2+K1_11)*(v1_11/v1_1); %E_1=gm1^2*v1_1/2; %detap1_1=E_1*(A_1+B_1+C_1-D_1); detap1_1=((1-sig1^2+K1_1)+2*(v1_11/v1_1-1)+(4*f1*2*l1*vm1)/(de1*v1_1)-((1-sig1^2+K1_11)*v1_11/v1_1))*(gm1^2*v1_1/2); %計(jì)算冷測(cè)阻力 %冷測(cè)為雙流程，對(duì)于單流程預(yù)設(shè)detap20=80pa，采用相同的算法 %設(shè)置判斷語(yǔ)句 %判斷Re1,根據(jù)雷諾數(shù)爬滿斷面熱側(cè)空氣的f if Re2<=1000f2=7.661*(ls2/de2)^0.384 * ((Pf2-Deltaf2)/(s2-Deltaf2))^(-0.092)*Re2^(-0.712); else if Re2>=2000f2=1.136*(ls2/de2)^(-0.781)*(Deltaf2/de2)^0.534*Re2^(-0.198);elseRe_f=41*((ls2/de2)^0.772)*((Pf2-Deltaf2)/(s2-Deltaf2))^(-0.179)*Re2^(-1.04);if Re2>=Re_jf2=7.661*(ls2/de2)^0.384 * ((Pf2-Deltaf2)/(s2-Deltaf2))^(-0.092)*Re2^(-0.712);else f2=1.136*(ls2/de2)^(-0.781)*(Deltaf2/de2)^0.534*Re2^(-0.198);endend end %此時(shí)已經(jīng)求得f2；計(jì)算實(shí)際的detap2，實(shí)際的命名為detap2_1; %A_2=1-sig2^2+K2_1; %B_2=2*(v2_11/v2_1-1); %C_2=(4*f2*2*l2*vm1)/(de2*v2_1); %D_2=(1-sig2^2+K2_11)*(v2_11/v2_1); %E_2=gm1^2*v2_1/2; %detap2_1=E_2*(A_2+B_2+C_2-D_2); detap2_1=(((1-sig2^2+K2_1)+2*(v2_11/v2_1-1)+(4*f2*2*l2*vm2)/(de2*v2_1)-(((1-sig2^2+K2_11)*v2_11)/v2_1))*(gm2^2*v2_1/2)); detap2_2=detap2_1*2;%雙流程壓降為detap2_2 % %% 強(qiáng)度校核%換熱器熱側(cè)空氣的入口溫度為90°，由于鋁材料隨著溫度升高強(qiáng)度指標(biāo)有所下降，故選取100°時(shí)的強(qiáng)度進(jìn)行強(qiáng)度校核。%防銹鋁LF21-M的需應(yīng)力%sig_Al = 75*1000*1000;sig_Al=23.75*1000*1000;%計(jì)算翅片厚度，校核結(jié)果的厚度用Deltaf1_1表示Deltaf1_1=((p1_1*2.5)/(sig_Al*1))+0.01;%這里乘1是指phi，phi是一個(gè)常數(shù)，取1.Deltaf2_1=((p2_1*2.5)/(sig_Al*1))+0.01;%算出來(lái)的東西單位是毫米mm。%計(jì)算隔板厚度，校核結(jié)果用Deltaf1_2表示Deltaf1_2=2.5*sqrt((6*p1_1)/(8*sig_Al))+0.05;Deltaf2_2=2*sqrt((6*p1_1)/(8*sig_Al))+0.05;

總結(jié)

以上是生活随笔為你收集整理的板翅式换热器（锯齿形）matlab的的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。