function [ ] = mjjs( )

%用于计算蒸发器的程序

Q0=6.24; % 制冷量 kw

t0=-28; %蒸发温度 ℃

mm=148.0; %循环制冷量 kg/h

t1=-21; %空气送风温度 ℃ % 送风相对湿度 95%

t2=-18; %空气回风温度 ℃ % 回风相对湿度 90%

% 含湿量及h值由hd函数计算得到

d1=0.53 %送风空气含湿量

d2=0.67 %回风空气含湿量

h1=-19.82 %送风焓值

h2=-16.46 %回风焓值

va=0.7202 % 空气的比体积，按 送、回风平均温度确定

% 风量的计算 m^3/h

qva=(Q0*va*3600)/(h2-h1)

sprintf('计算得到的风量为%f m^3/h，请选择风机后输入标准风量。',qva)

qva=input('请输入标准风量(m^3.h):');

% 送风时空气通过蒸发器的比焓差 kj/kg

drh=(Q0*3600)/(qva/va)

% 更新风量单位 为m^3/s

qva=qva/3600;

% 送风空气的比焓值

h1=h2-drh

sprintf('请根据求得h1值，查送风温度，并进行校正 原始温度为%f℃',t1)

t1=input('最新的送风温度(℃)：');

% 蒸发器的结构的初步规划

% 取Φ10×1的铜管，外套铝片。铝片间距b=10mm，厚度drta=0.2mm，管排数N=5，每排管子数目n=12

% 管子三角形排列，管间距s=0.028m

% 铝片间距 m

b=0.01;

% 铝片厚度 m

drta= 0.0002;

% 管排数

N=4

% 每排管子数目n

n=18

% 管间距 s m

s=0.028;

%蛇形管内径 m

Di=0.008;

%蛇形管外径 m

d0=0.01;

sprintf('管子及肋片采用默认');

% 最窄流通面积Fe与迎风面积Fy之比

Fe_Fy = (s-d0)*(b-drta)/(s*b)

% 最下流速取 4.5m/s

u=4.5;

% 最窄流通面积 m^2

Fe=qva/u

% 迎风面积 m^2

Fy=Fe/Fe_Fy

% 蒸发器高度

A= n*s

% 厚度

B=N*s*cos(pi/6)

% 每根管子的长度

ll=Fy/A

% 管子总长度

l=ll*n*N

% 六角形肋片单侧表面积

F= (6*s*0.5) *(s*cot(pi/3)*0.5) - (pi*d0*d0*0.25)

% 每米管子上面的肋片数

nf=1/b

% 每米肋片表面积

Ffl =2*F*nf

% 每米管子上面铜管的表面积

Frl=pi*d0*(1-nf*drta)

% 每米的总外表面积

Fofl = Ffl+Frl

% 总外表面积

Fof=Fofl*l

% 总的肋片面积

Ff=Ffl*l

% 总的铜管外边面积

Fr=Frl*l

% 总的内表面积

Fi=pi*Di*l

% 外表面积与内表面积之比

beita=Fof/Fi

% beita=6; %总的表面积与内表面积之比

% 空气侧放热表面传热系数 w/(m^2 ?℃)

hof=112.4;

sprintf('空气侧的表面传热系数与肋片参数有关，如果修改肋片参数，请重新计算');

% %质量流速kg/(m^2·s)

vm=(mm/3600)/(N*0.25*pi*Di^2)

r=225.706; %汽化热 kj/kg

x2=1; %离开蒸发器时的干度

rog=7.96; %制冷剂饱和蒸汽密度

ro1=1374.1; %制冷剂饱和液体密度

% Di=0.008; %蛇形管内径 m

miu1=297.7E-6;%液体粘度 pa·s

k1=107.5E-3; %制冷剂液体导热率 w/(m·k)

Ffl=2.2; %制冷剂性质的无量纲数

pr1=3.068; %制冷剂液体普朗特数

tao=10; %融霜周期 h

%Fof=15.855; %总的外表面积 m^2

%Ff=12.6; %肋片面积 m^2

%Fr=