資源簡介
該程序組共包含10個子程序,全部為計算水和水蒸汽性質(zhì)的子程序,采用的是國際公式化委員會制定的水和水蒸氣熱力性質(zhì)(IFC67)公式。各子程序分別是:
TSK.m 求某下壓力飽和溫度。
PSK.m 某溫度下飽和壓力。
HS.m 已知比焓、比熵,求其它性質(zhì)。
PX.m 已知壓力、干度,求其它性質(zhì)。
PV.m 已知壓力、比熵,求其它性質(zhì)。
PTG.m 已知壓力、溫度,求飽和汽、過熱蒸汽的性質(zhì)。
PTF.m 已知壓力、溫度,求飽和水、過冷水的性質(zhì)。
PT.m 已知壓力、溫度,求其它性質(zhì)。
PS.m 已知壓力、比熵,求其它性質(zhì)。
PH.m 已知壓力、比焓,求其它性質(zhì)
該程序是汽輪機(jī)設(shè)計和熱力系統(tǒng)設(shè)計的好幫手,減少了查表的麻煩。

代碼片段和文件信息
function?[X?PT?V]?=?HS(H?S)
%************************************************
%?已知比焓、比熵求過熱蒸汽、飽和蒸汽、濕蒸汽性質(zhì)函數(shù)
%??Parameters:?P??????壓力?MPa?
%??????????????T??????溫度??
%??????????????V??????比容??
%??????????????H??????焓??
%??????????????S??????熵??
%??????????????X??????干度(過熱度)
%?****作者:王雷?zrqwl2003@126.com*****************
%?************************************************
G=[-32230.91316763.0198-2684.70102126.0153826.08691173-0.515969833];
w=[0.08275963?-1.2121626?-73.614096.6074894-0.4065599];
HG?=?0;
for?i?=?1?:1:?6
????HG?=?HG?+?G(i)?*?S?^?(i?-?1);
end
if?(H?>?HG)?
????X?=?1;
????ZH?=?H?/?4186.8;
????ZS?=?0.238846?*?S;
????n?=?1;
????while(1)
????????ZT?=?(sqrt(0.035081?+?0.085?*?(ZH?-?0.4949))?-?0.1873)?/?0.085;
????????ZP?=?exp(13.064?-?9.072?*?ZS?+?3.3982?*?log(ZT));
????????T?=?1000?*?ZT?-?273.15;
????????P?=?98.0665?*?ZP;
????????TS?=?TSK(P);
????????if?(T?>=?TS)?
????????????ZHB?=?ZH;
????????????ZSB?=?ZS;
????????else
????????????T?=?TS;
????????????ZT?=?(T?+?273.15)?/?1000;
????????????ZHB?=?((0.085?*?ZT?+?0.1873)?^?2?-?0.035081)?/?0.085?+?0.4949;
????????????ZSB?=?(13.064?-?log(ZP)?+?3.3982?*?log(ZT))?/?9.072;
????????end
????????X?=?T?-?TS;
????????if?(X?>=?1)
????????else
????????????X?=?X?+?1;
????????end
????????[V?HB?SB]=PTG(P?T?);
????????if?((abs(H?-?HB)?0.0005)?&&?(abs(S?-?SB)?0.000001))
????????????break;
????????end
????????if?(n?>?1)
????????????ZH?=?ZHB?+?(H?-?HB)?*?(ZHB?-?ZHA)?/?(HB?-?HA);
????????????ZS?=?ZSB?+?(S?-?SB)?*?(ZSB?-?ZSA)?/?(SB?-?SA)?;
????????else
????????????ZH?=?ZHB?+?0.0002388?*?(H?-?HB);
????????????ZS?=?ZSB?+?0.2388?*?(S?-?SB);
????????end
????????n?=?n?+?1;
????????HA?=?HB;
????????SA?=?SB;
????????ZHA?=?ZHB;
????????ZSA?=?ZSB;
????????ZHB?=?ZH;
????????ZSB?=?ZS;
????end
else
????TB?=?H?/?S;
????while(1)
????????GF?=?0;
????????for?i?=?1?:1:?5
????????????GF?=?GF?+?w(i)?*?(TB?*?0.01)?^?(i?-?1);
????????end
????????T?=?(H?-?GF)?/?S?-?273.15;
????????if?(abs(T?-?TB)?0.001)
????????????break;
????????end
????????TB?=?T;
????end
????P?=?PSK(T);
????[VG?HG?SG]=PTG(PT);
????[VF?HF?SF]=PTF(PT);
????X?=?(H?-?HF)?/?(HG?-?HF);
????V?=?VF?+?X?*?(VG?-?VF);
end
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????2276??2013-03-11?18:28??IFC67\HS.m
?????文件???????1901??2012-07-13?21:31??IFC67\PH.m
?????文件???????2050??2012-07-13?21:31??IFC67\PS.m
?????文件????????564??2012-07-13?21:31??IFC67\PSK.m
?????文件????????692??2012-07-13?21:32??IFC67\PT.m
?????文件???????2181??2012-07-13?21:32??IFC67\PTF.m
?????文件???????4353??2012-07-13?21:32??IFC67\PTG.m
?????文件???????1097??2012-07-13?21:32??IFC67\PV.m
?????文件????????555??2012-07-13?21:32??IFC67\PX.m
?????文件????????456??2012-07-13?21:32??IFC67\TSK.m
?????目錄??????????0??2014-04-29?15:35??IFC67
-----------?---------??----------?-----??----
????????????????16125????????????????????11
評論
共有 條評論