过参课设-第二类标准节流装置设计

A = 1.432677533; % 设弦法不变量
Kappa = 1.3; % 等熵指数
ReD = 3.15e7; % 雷诺数
D = 306.88; % 管道直径
P_1 = 16.4e6; % 工作压力
Delta_P =50170; % 常用流量下差压
% 迭代变量
n = 1; % 循环次数
X = zeros(); % 设弦法过程变量
Beta = zeros(); % 直径比 
C = zeros(); % 记录流出系数
Epsilon = zeros(); % 膨胀系数
Delta = zeros(); % 差值
E = zeros(); % 记录误差
% 设置初始流出系数和膨胀系数
C(1) = 1;
Epsilon(1) = 1;
% 循环体Epsilon
while (n==1 || E(n)>1*10^(-7)) % 误差小于10^(-7)时跳出循环
% 迭代次数增加
n = n+1;
% 计算中间变量 
X(n) = A/(C(n-1)*Epsilon(n-1)); % 更新设弦法中间变量
Beta(n) = (X(n)^2/(X(n)^2+1))^(0.25); % 更新直径比
% 标准孔板膨胀系数
Epsilon(n) = 1 - (0.41+0.35*(Beta(n)^4))*(Delta_P)/(Kappa*P_1);% 更新膨胀系数
% 标准孔板流出系数
C(n) =0.5959+0.0312*(Beta(n)^2.1)-0.184*(Beta(n)^8)+0.0029*(Beta(n)^2.5)*((1e6/ReD)^0.75)+2.286*(Beta(n)^4)/(D*(1-(Beta(n)^4)))-0.8560*(Beta(n)^3)/D; % 更新流出系数
Delta(n) = A - X(n)*C(n)*Epsilon(n); % 更新差值

E(n) = abs(Delta(n)/A); % 计算误差 
end