?clc?
?clear
?
%%%%?SimpleBraytonCycleLayout
?Q=input（‘熱源功率（MW）:‘）;?????????%%26?440?38?8?87?85?136?2.6
?T4=input（‘熱源溫度（C）:‘）+273.15;
?T1=input（‘冷凝器出口溫度（C）:‘）+273.15;
%??P4=input（‘熱源壓力（MPa）:‘）*1000;
?P1=input（‘運行壓力（MPa）:‘）*1000;
?E1=input（‘透平效率（%）:‘）/100;
?E2=input（‘壓氣機效率（%）:‘）/100;
?M=input（‘工質流量（kg/s）:‘）;
?Et=input（‘壓比:‘）;%壓力比
?P4=P1*Et;

H4=refpropm（‘h‘‘T‘T4‘P‘P4‘CO2‘）;
H1=refpropm（‘h‘‘T‘T1‘P‘P1‘CO2‘）;

%%汽輪機出口5
P5=P1;
S5=refpropm（‘S‘‘T‘T4‘P‘P4‘CO2‘）;
H5=refpropm（‘h‘‘P‘P5‘S‘S5‘CO2‘）;
H45=（H5-H4）*E1;
T5=refpropm（‘t‘‘P‘P5‘H‘H4+H45‘CO2‘）;

%%換熱器熱端
P3=P4;
H34=Q*1000/M;
H3=H4-H34;
%?H3=H5-H6+H2;
T3=refpropm（‘t‘‘P‘P3‘H‘H3‘CO2‘）;

%%壓氣機出口2
P2=P4;
S2=refpropm（‘S‘‘T‘T1‘P‘P1‘CO2‘）;
H2=refpropm（‘h‘‘P‘P2‘S‘S2‘CO2‘）;
H12=（H2-H1）/E2;
T2=refpropm（‘t‘‘P‘P2‘H‘H1+H12‘CO2‘）;

%%換熱器冷端
H6=H5-H3+H2;
P6=P1;
S6=refpropm（‘S‘‘P‘P6‘H‘H6‘CO2‘）;
T6=refpropm（‘t‘‘P‘P6‘H‘H6‘CO2‘）;

%%透平輸出功
Wt=M*（-H45）/1000000;
%%壓氣機耗功率
Wc=M*（H12）/1000000;
%%系統效率
n=（Wt-Wc）/Q*100;

disp（[‘透平輸出功為‘num2str（Wt）‘MW‘]）;
disp（[‘壓氣機耗功為‘num2str（Wc）‘MW‘]）;
disp（[‘熱源功率為‘num2str（Q）‘MW‘]）;
disp（[‘系統效率為‘num2str（n）‘%‘]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1240??2017-11-23?14:27??超臨界二氧化碳布雷頓簡單循環MATLAB計算程序\SimpleBraytonCycleLayout_1.m
?????文件????????1303??2017-11-23?14:27??超臨界二氧化碳布雷頓簡單循環MATLAB計算程序\SimpleBraytonCycleLayout_2.m
?????目錄???????????0??2017-12-05?16:31??超臨界二氧化碳布雷頓簡單循環MATLAB計算程序\