列管式换热器设计程序

clc
clear
t11=input（‘t11 ‘）;
t12=input（‘t12 ‘）;
t21=input（‘t21 ‘）;
t22=input（‘t22 ‘）;
P1=input（‘p1 ‘）;
P2=input（‘p2 ‘）;
M1=input（‘M1 ‘）;


tm1=（t11+t12）/2;
cp1=input（‘cp1 单位kJ/kg.K ‘）;
cp1=cp1*1000;
p1=input（‘p1密度 ‘）;%密度
u1=input（‘u1 ‘）;
lm1=input（‘lm1 ‘）;
Pr1=（u1*cp1）/lm1;
tm2=（t21+t22）/2;
cp2=input（‘cp2 单位kJ/kg.K‘）;
cp2=cp2*1000;
p2=input（‘p2 密度 ‘）;%密度



lm2=input（‘lm2 ‘）;
u2=input（‘u2 ‘）;
Pr2=（u2*cp2）/lm2;

nL=0.98;      %取用
Q=M1*cp1*（t11-t12）*nL;
M2=Q/（cp2*（t22-t21））;
deta_t1mc=abs（（t11-t12）-（t22-t21））/（log（abs（（t11-t12）/（t22-t21））））;
P=（t22-t21）/（t11-t21）;
R=（t11-t12）/（t22-t21）;


%（1-2）型，具体见p18
%w1=sqrt（R^2+1）/（R-1）/2;
%w2=log（（1-R）/（1-P*R））;
%w3=2*（1+sqrt（（1-P）*（1-P*R）））-P*（1+R-sqrt（R^2+1））;
%w4=2*（1+sqrt（（1-P）*（1-P*R）））-P*（1+R+sqrt（R^2+1））;