MATLAB的答案

未定义的函数或变量“valg_design_rotspeed”。

1视图(30天)
罗希特帕特尔
罗希特帕特尔 2018年3月24日
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%设计参数
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
n = valg_design_rotspeed;
param_omega = n * 2 *π/ 60;
Q = valg_design_flow;
H = valg_design_nethead;
g = valg_design_gP;
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算特定速度
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
n_sQ = n * Q / H ^ ^ 0.5 0.75;
n_sstar = n_sQ / g ^ 0.75;
n_sQstar = (1000/60) * n_sstar;
集(hand_calc_sp_speed,…
“字符串”num2str(0.1 *轮(10 * n_sQ)),…
“价值”,n_sQ);
集(hand_calc_sp_speed_dimless,…
“字符串”num2str(0.1 *轮(10 * n_sstar)),…
“价值”,n_sstar);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%计算出口价值
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算出口区域
% a₂不考虑叶片厚度
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
D_2 = valg_design_D2;
R_2 = 0.5 * D_2;
D_2_red = D_2 / R_2;%应该是2
R_2_red = 0.5 * D_2_red;%定义为11
b_2 = valg_design_b2;
b_2_red = b_2 / R_2;
集(hand_value_b2,…
“字符串”num2str (b_2),…
“价值”,b_2);
a₂= R_2 * 2 *π* b_2;
A_2_red = R_2_red * 2 *π* b_2_red;% a₂/ a₂;%定义为1
param_A2 = a₂;
param_A2_red = A_2_red;
集(hand_calc_A2,…
“字符串”num2str (a₂),…
“价值”,a₂);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算Peripherial速度,U_2
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
U_2 = R_2 * param_omega;
U_2_red = U_2 / U_2;%定义为1;
集(hand_calc_U2,…
“字符串”num2str (U_2),…
“价值”U_2);
集(hand_calc_U2_red,…
“字符串”num2str (U_2_red),…
“价值”,U_2_red);
集(hand_calc_U2slip,…
“字符串”num2str (U_2),…
“价值”U_2);
集(hand_calc_U2slip_red,…
“字符串”num2str (U_2_red),…
“价值”,U_2_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算出口子午速度分量绝对速度,C_m2
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C_m2 = Q / a₂;
C_m2_red = C_m2 / U_2;
集(hand_calc_Cm2,…
“字符串”num2str (C_m2),…
“价值”,C_m2);
集(hand_calc_Cm2_red,…
“字符串”num2str (C_m2_red),…
“价值”,C_m2_red);
集(hand_calc_Cm2slip,…
“字符串”num2str (C_m2),…
“价值”,C_m2);
集(hand_calc_Cm2slip_red,…
“字符串”num2str (C_m2_red),…
“价值”,C_m2_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算相对出口速度,W_2
%设置W_2slip值
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
beta_2 = valg_design_beta2;
W_2 = C_m2 /信德(beta_2);
W_2_red = W_2 / U_2;
集(hand_calc_W2,…
“字符串”num2str (W_2),…
“价值”,W_2);
集(hand_calc_W2_red,…
“字符串”num2str (W_2_red),…
“价值”,W_2_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算出口圆周速度分量的相对速度,W_U2
%设置W_U2slip值
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
W_U2 = C_m2 /罐内(beta_2);
W_U2_red = W_U2 / U_2;
集(hand_calc_WU2,…
“字符串”num2str (W_U2),…
“价值”,W_U2);
集(hand_calc_WU2_red,…
“字符串”num2str (W_U2_red),…
“价值”,W_U2_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算C_U2
%设置W_U2slip值
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C_U2 = U_2-W_U2;
C_U2_red = C_U2 / U_2;
集(hand_calc_CU2,…
“字符串”num2str (C_U2),…
“价值”,C_U2);
集(hand_calc_CU2_red,…
“字符串”num2str (C_U2_red),…
“价值”,C_U2_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%计算进气值
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算叶轮直径,D_1
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
D_1 = valg_design_D1;
D_1i = valg_design_D1i;
D_n = valg_design_Dn;
R_1i = 0.5 * D_1i;
D_1i_red = D_1i / R_2;
D_n_red = D_n / R_2;
R_1i_red = 0.5 * D_1i_red;
R_n_red = 0.5 * D_n_red;
R_1 = 0.5 * D_1;
D_1_red = D_1 / R_2;
R_1_red = 0.5 * D_1_red;
集(hand_calc_D1,…
“字符串”num2str (D_1),…
“价值”,D_1);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算入口区域,A_1
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
A_1 =(π)/ 4 * (D_1 ^ 2-D_n ^ 2);
A2_A1 = a₂/ A_1;
A_1_red = A_2_red / A2_A1;
集(hand_calc_A1,…
“字符串”num2str (A_1),…
“价值”,A_1);
集(hand_calc_ratA,…
“字符串”num2str (A2_A1),…
“价值”,A2_A1);
param_A1 = A_1;
param_A1_red = A_1_red;
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算进口意味着简化直径,D_1m
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
D_1m =√(D_1 ^ 2 + D_1i ^ 2) * 0.5);
D_1m_red = D_1m / R_2;
集(hand_calc_D1m,…
“字符串”num2str (D_1m),…。
“价值”,D_1m);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算进口圆周速度,U_1
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
U_1_ring = param_omega * D_1 * 0.5;
U_1_ring_red = U_1_ring / U_1_ring;%定义为1
集(hand_calc_U1_ring,…
“字符串”num2str (U_1_ring),…
“价值”,U_1_ring);
集(hand_calc_U1_ring_red,…
“字符串”num2str (U_1_ring_red),…
“价值”,U_1_ring_red);
U_1_merid = param_omega * D_1m * 0.5;
U_1_merid_red = U_1_merid / U_1_ring;
集(hand_calc_U1_merid,…
“字符串”num2str (U_1_merid),…
“价值”,U_1_merid);
集(hand_calc_U1_merid_red,…
“字符串”num2str (U_1_merid_red),…
“价值”,U_1_merid_red);
U_1_boss = param_omega * D_1i * 0.5;
U_1_boss_red = U_1_boss / U_1_ring;
集(hand_calc_U1_boss,…
“字符串”num2str (U_1_boss),…
“价值”,U_1_boss);
集(hand_calc_U1_boss_red,…
“字符串”num2str (U_1_boss_red),…
“价值”,U_1_boss_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算进口子午组件的绝对速度,C_m1
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C_m1_ring = Q / A_1;
C_m1_ring_red = C_m1_ring / U_1_ring;
集(hand_calc_Cm1_ring,…
“字符串”num2str (C_m1_ring),…
“价值”,C_m1_ring);
集(hand_calc_Cm1_ring_red,…
“字符串”num2str (C_m1_ring_red),…
“价值”,C_m1_ring_red);
C_m1_merid = Q / A_1;
C_m1_merid_red = C_m1_merid / U_1_ring;
集(hand_calc_Cm1_merid,…
“字符串”num2str (C_m1_merid),…
“价值”,C_m1_merid);
集(hand_calc_Cm1_merid_red,…
“字符串”num2str (C_m1_merid_red),…
“价值”,C_m1_merid_red);
C_m1_boss = Q / A_1;
C_m1_boss_red = C_m1_boss / U_1_ring;
集(hand_calc_Cm1_boss,…
“字符串”num2str (C_m1_boss),…
“价值”,C_m1_boss);
集(hand_calc_Cm1_boss_red,…
“字符串”num2str (C_m1_boss_red),…
“价值”,C_m1_boss_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算进气角在老板,beta_1
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
beta_1_ring = atand (C_m1_ring / U_1_ring);
集(hand_calc_beta1_ring,…
“字符串”num2str (beta_1_ring),…
“价值”,beta_1_ring);
beta_1_merid = atand (C_m1_merid / U_1_merid);
集(hand_calc_beta1_merid,…
“字符串”num2str (beta_1_merid),…
“价值”,beta_1_merid);
beta_1_boss = atand (C_m1_boss / U_1_boss);
集(hand_calc_beta1_boss,…
“字符串”num2str (beta_1_boss),…
“价值”,beta_1_boss);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算其他usesful速度
% C_U1环,merid和老板
% W_1环,merid和老板
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C_U1_ring = U_1_ring - Q / (A_1 *罐内(beta_1_ring));
C_U1_merid = U_1_merid - Q / (A_1 *罐内(beta_1_merid));
C_U1_boss = U_1_boss - Q / (A_1 *罐内(beta_1_boss));
W_1_ring = C_U1_ring /信德(beta_1_ring);
W_1_merid = C_U1_merid /信德(beta_1_merid);
W_1_boss = C_U1_boss /信德(beta_1_boss);
C_U1_ring_red = C_U1_ring / U_1_ring;
C_U1_merid_red = C_U1_merid / U_1_ring;
C_U1_boss_red = C_U1_boss / U_1_ring;
W_1_ring_red = W_1_ring / U_1_ring;
W_1_merid_red = W_1_merid / U_1_ring;
W_1_boss_red = W_1_boss / U_1_ring;
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%计算滑移值
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
Z = valg_design_Z;
eps_lim = exp (-8.16 * sin (beta_2) / Z);
D_2 = valg_design_D2;
千瓦= 1;
f1 = 0.98;
如果D_1m / D_2 > eps_lim
千瓦= (1 - ((D_1m / D_2-eps_lim) / (1-eps_lim)) ^ 3;
结束
slip_coeff = f1 * (1-sqrt(信德(beta_2)) / Z ^ 0.70) *千瓦;
集(hand_calc_slip_coeff,…
“字符串”num2str (slip_coeff),…
“价值”,slip_coeff);
%计算叶片堵塞
e_2 = valg_design_tout;
block_2 = (1-e_2 * Z /(π* D_2 *信德(beta_2))) ^ 1;
%计算CU2slip
% C_U2slip_old = C_U2 - (1-slip_coeff) * U_2;
C_U2slip = U_2 * (slip_coeff-C_m2 * block_2 / (U_2 *踮(beta_2)));
C_U2slip_red = C_U2slip / U_2;
集(hand_calc_CU2slip,…
“字符串”num2str (C_U2slip),…
“价值”,C_U2slip);
集(hand_calc_CU2slip_red,…
“字符串”num2str (C_U2slip_red),…
“价值”,C_U2slip_red);
%计算beta_2slip和slip_angle
beta_2slip = atand (C_m2_red / (1-C_U2slip_red));
slip_angle = beta_2-beta_2slip;
%计算W_2slip
W_2slip = C_m2 /信德(beta_2slip);
W_2slip_red = W_2slip / U_2;
集(hand_calc_W2slip,…
“字符串”num2str (W_2slip),…
“价值”,W_2slip);
集(hand_calc_W2slip_red,…
“字符串”num2str (W_2slip_red),…
“价值”,W_2slip_red);
%计算W_U2slip
W_U2slip = C_m2 /罐内(beta_2slip);
W_U2slip_red = W_U2slip / U_2;
集(hand_calc_WU2slip,…
“字符串”num2str (W_U2slip),…
“价值”,W_U2slip);
集(hand_calc_WU2slip_red,…
“字符串”num2str (W_U2slip_red),…
“价值”,W_U2slip_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算绝对角度在出口,alpha_2 alpha_2slip
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
alpha_2 = atand (C_m2 / C_U2);
alpha_2slip = atand (C_m2 / C_U2slip);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%更新出口角度
%设置beta2slip值
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
集(hand_value_beta2,…
“字符串”num2str (beta_2),…
“价值”,beta_2);
集(hand_calc_beta2,…
“字符串”num2str (valg_design_beta2))
集(hand_calc_beta2slip,…
“字符串”num2str (beta_2slip))
%设置(hand_value_slip,……
% '字符串',num2str (slip_angle)……
%“价值”,slip_angle);
集(hand_calc_slip,…
“字符串”num2str (slip_angle),…
“价值”,slip_angle);
集(hand_calc_alpha2,…
“字符串”num2str (alpha_2),…
“价值”,alpha_2);
集(hand_calc_alpha2slip,…
“字符串”num2str (alpha_2slip),…
“价值”,alpha_2slip);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%计算最优水力效率
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
如果问> 1
一个= 0.5;
其他的
= 1;
结束
m = 0.08 * 0.15 * (1 / Q) ^ * (45 / n_sQ) ^ 0.06;
eta_opt = 1 - 0.065 * (1 / Q) ^ m - 0.23 * 0.3 (log10 (n_sQ / 23)) ^ 2 * (1 / Q) ^ 0.05;
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%更新泵特点
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
集(hand_calc_eta_opt,…
“字符串”num2str (eta_opt));
%设置计算理论
集(hand_calc_theoretical_head,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 / valg_design_gP));
%设置计算理论在酒吧
集(hand_calc_theoretical_head_bar,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 * valg_design_rhoP * 10 ^ 5));
[hand_axes_charac, calc_net_head x_max、y_max calc_pfleiderer_slip]…
= draw_charac_curves (…
hand_axes_charac,…
valg_design_flow,…
get (hand_calc_U2,“价值”),…
get (hand_calc_U1_merid,“价值”),…
get (hand_calc_CU2,“价值”),…
valg_design_D2,…
D_1,…
get (hand_value_b2,“价值”),…
valg_design_beta2,…
beta_1_merid,…
A_1,…
get (hand_calc_slip,“价值”),…
param_omega,…
valg_charac_friction,…
valg_charac_shockloss,…
valg_design_gP,…
valg_design_Z,…
valg_design_nethead…
);
%设置计算理论的头体积损失
集(hand_calc_theoretical_head_leakage,“字符串”num2str (y_max));
%设置计算理论在酒吧与体积损失
集(hand_calc_theoretical_head_leakage_bar,“字符串”num2str (y_max * valg_design_rhoP * valg_design_gP * 10 ^ 5));
%设置计算基于design_flow Pfleiderer滑滑因素
集(hand_calc_pfleiderer,“字符串”num2str (calc_pfleiderer_slip));
%设置计算净头
集(hand_calc_net_head,“字符串”num2str (calc_net_head));
%设置计算净酒吧
集(hand_calc_net_head_bar,“字符串”num2str (calc_net_head * valg_design_rhoP * valg_design_gP * 10 ^ 5));
% %设置(hand_calc_pfleiderer“字符串”,num2str (pfleiderer));
%设置(hand_axes_charac,……
% XLim, [0 x_max],……
% YLim, [0 y_max]);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%更新入口速度的阴谋
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
[hand_plot_U1_ring, hand_text_U1_ring,…
hand_plot_C1_ring hand_text_C1_ring,…
hand_plot_W1_ring hand_text_W1_ring,…
hand_plot_U1_merid,…
hand_plot_W1_merid hand_text_W1_merid,…
hand_plot_U1_boss,…
hand_plot_W1_boss hand_text_W1_boss] =…
inlet_velocity_plot_update (…
hand_plot_U1_ring hand_text_U1_ring,…
hand_plot_C1_ring hand_text_C1_ring,…
hand_plot_W1_ring hand_text_W1_ring,…
hand_plot_U1_merid,…
hand_plot_W1_merid hand_text_W1_merid,…
hand_plot_U1_boss,…
hand_plot_W1_boss hand_text_W1_boss,…
U_1_ring_red C_m1_ring_red,…
U_1_merid_red C_m1_merid_red,…
U_1_boss_red C_m1_boss_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%更新出口速度的阴谋
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
[hand_plot_U2, hand_text_U2,…
hand_plot_C2 hand_text_C2,…
hand_plot_W2 hand_text_W2,…
hand_plot_CU2 hand_text_CU2,…
hand_plot_Cm2 hand_text_Cm2,…
hand_plot_C2slip hand_text_C2slip,…
hand_plot_W2slip hand_text_W2slip,…
hand_plot_CU2slip hand_text_CU2slip,…
hand_plot_Cm2slip] =…
outlet_velocity_plot_update (…
hand_plot_U2 hand_text_U2,…
hand_plot_C2 hand_text_C2,…
hand_plot_W2 hand_text_W2,…
hand_plot_CU2 hand_text_CU2,…
hand_plot_Cm2 hand_text_Cm2,…
hand_plot_C2slip hand_text_C2slip,…
hand_plot_W2slip hand_text_W2slip,…
hand_plot_CU2slip hand_text_CU2slip,…
hand_plot_Cm2slip,…
U_2_red、C_m2_red C_U2_red C_U2slip_red);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%设计参数
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
D_2 = valg_design_D2;
D_3 = valg_design_D3;
D_4 = valg_design_D4;
b_3 = valg_design_b3;
b_4 = valg_design_b4;
a_3 = valg_design_a3;
n = valg_npoints_cross_diffuser;
m = valg_npoints_stream_diffuser;
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%计算进气值
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算进气速度绝对的圆周速度分量,C_u3
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C_U2slip =得到(hand_calc_CU2slip,“价值”);
C_U3 = D_2 * C_U2slip / D_3;
C_U3block = C_U3;
集(hand_calc_CU3,…
“字符串”num2str (C_U3),…
“价值”,C_U3);
集(hand_calc_CU3block,…
“字符串”num2str (C_U3block),…
“价值”,C_U3block);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算进气速度绝对的圆周速度分量,C_u3,
%不滑
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
C_U2 =得到(hand_calc_CU2,“价值”);
C_U3noslip = D_2 * C_U2 / D_3;
C_U3noslipblock = C_U3noslip;
集(hand_calc_CU3noslip,…
“字符串”num2str (C_U3noslip),…
“价值”,C_U3noslip);
集(hand_calc_CU3noslipblock,…
“字符串”num2str (C_U3noslipblock),…
“价值”,C_U3noslipblock);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算子午之前绝对速度扩散器的入口
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
Q = valg_design_flow;
C_m3 = Q /(π* D_3 * b_3);
集(hand_calc_Cm3,…
“字符串”num2str (C_m3),…
“价值”,C_m3);
C3 =√C_m3 ^ 2 + C_U3 ^ 2);
集(hand_calc_C3,…
“字符串”num2str (C3),…
“价值”,C3);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算速度扩散器的喉咙
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
Q = valg_design_flow;
Z_Le = valg_design_ZLe;%扩散叶片的数量
C_q3 = Q / (Z_Le * a_3 * b_3);
集(hand_calc_CQ3,…
“字符串”num2str (C_q3),…
“价值”,C_q3);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算扩散器入口前绝对流动角
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
alpha_3 = atand (C_m3 / C_U3);
集(hand_calc_alpha3,…
“字符串”num2str (alpha_3),…
“价值”,alpha_3);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%计算流条件后扩散器入口(与叶片堵塞)
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
e_3 = valg_design_e3;%在扩压器进口叶片厚度
alpha_3B = valg_design_alpha3B;在扩压器进口%桨叶角
block_3 = (1-Z_Le * e_3 /(π* D_3 *信德(alpha_3B))) ^ 1;%叶片扩压器进口堵塞
C_m3block = Q * block_3 /(π* D_3 * b_3);%子午流后扩散器的组件
入口
alpha_3block = atand (C_m3block / C_U3);%流角后扩散器入口
C3block =√C_m3block ^ 2 + C_U3 ^ 2);
集(hand_calc_Cm3block,…
“字符串”num2str (C_m3block),…
“价值”,C_m3block);
集(hand_calc_alpha3block,…
“字符串”num2str (alpha_3block),…
“价值”,alpha_3block);
集(hand_calc_C3block,…
“字符串”num2str (C3block),…
“价值”,C3block);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%更新入口速度的阴谋
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
[hand_plot_C3, hand_text_C3,…
hand_plot_CU3 hand_text_CU3,…
hand_plot_Cm3 hand_text_Cm3,…
hand_plot_C3block hand_text_C3block,…
hand_plot_CU3block hand_text_CU3block,…
hand_plot_Cm3block hand_text_Cm3block] =…
inlet_velocity_plot_diffuser_update (…
hand_plot_C3 hand_text_C3,…
hand_plot_CU3 hand_text_CU3,…
hand_plot_Cm3 hand_text_Cm3,…
hand_plot_C3block hand_text_C3block,…
hand_plot_CU3block hand_text_CU3block,…
hand_plot_Cm3block hand_text_Cm3block,…
C_U3 C_U3block,…
C_m3 C_m3block);
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%泵的特点,传统的1
% H = (1 / (1 + Cp)) (Hthinfty) -KfQ ^ 2-Ks (Q-Qopt) ^ 2
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%设置计算理论
集(hand_calc_theoretical_head_trad_1,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 / valg_design_gP));
%设置计算理论在酒吧
集(hand_calc_theoretical_head_bar_trad_1,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 * valg_design_rhoP * 10 ^ 5));
[hand_axes_charac_trad_1, calc_net_head_trad_1 x_max_trad_1、y_max_trad_1 calc_pfleiderer_slip_trad_1]…
= draw_charac_curves_trad_1 (…
hand_axes_charac_trad_1,…
valg_design_flow,…
get (hand_calc_U2,“价值”),…
get (hand_calc_U1_merid,“价值”),…
get (hand_calc_CU2,“价值”),…
valg_design_D2,…
D_1,…
get (hand_value_b2,“价值”),…
valg_design_beta2,…
beta_1_merid,…
A_1,…
get (hand_calc_slip,“价值”),…
param_omega,…
valg_charac_friction_trad_1,…
valg_charac_shockloss_trad_1,…
valg_design_gP,…
valg_design_Z,…
valg_design_nethead…
);
%设置计算理论的头体积损失
集(hand_calc_theoretical_head_leakage_trad_1,“字符串”num2str (y_max_trad_1));
%设置计算理论在酒吧与体积损失
集(hand_calc_theoretical_head_leakage_bar_trad_1,“字符串”num2str (y_max_trad_1 * valg_design_rhoP * valg_design_gP * 10 ^ 5));
%设置计算基于design_flow Pfleiderer滑滑因素
集(hand_calc_pfleiderer_trad_1,“字符串”num2str (calc_pfleiderer_slip_trad_1));
%设置计算净头
集(hand_calc_net_head_trad_1,“字符串”num2str (calc_net_head_trad_1));
%设置计算净酒吧
集(hand_calc_net_head_bar_trad_1,“字符串”num2str (calc_net_head_trad_1 * valg_design_rhoP * valg_design_gP * 10 ^ 5));
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%泵的特点,传统的2
% H = (1 / (1 + Cp)) (Hthinfty) -KfQ ^ 2 - (u1cu1 / g)
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%设置计算理论
集(hand_calc_theoretical_head_trad_2,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 / valg_design_gP));
%设置计算理论在酒吧
集(hand_calc_theoretical_head_bar_trad_2,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 * valg_design_rhoP * 10 ^ 5));
[hand_axes_charac_trad_2, calc_net_head_trad_2 x_max_trad_2、y_max_trad_2 calc_pfleiderer_slip_trad_2]…
= draw_charac_curves_trad_2 (…
hand_axes_charac_trad_2,…
valg_design_flow,…
get (hand_calc_U2,“价值”),…
get (hand_calc_U1_merid,“价值”),…
get (hand_calc_CU2,“价值”),…
valg_design_D2,…
D_1,…
get (hand_value_b2,“价值”),…
valg_design_beta2,…
beta_1_merid,…
A_1,…
get (hand_calc_slip,“价值”),…
param_omega,…
valg_charac_friction_trad_2,…
valg_charac_shockloss_trad_2,…
valg_design_gP,…
valg_design_Z,…
valg_design_nethead…
);
%设置计算理论的头体积损失
集(hand_calc_theoretical_head_leakage_trad_2,“字符串”num2str (y_max_trad_2));
%设置计算理论在酒吧与体积损失
集(hand_calc_theoretical_head_leakage_bar_trad_2,“字符串”num2str (y_max_trad_2 * valg_design_rhoP * valg_design_gP * 10 ^ 5));
%设置计算基于design_flow Pfleiderer滑滑因素
集(hand_calc_pfleiderer_trad_2,“字符串”num2str (calc_pfleiderer_slip_trad_2));
%设置计算净头
集(hand_calc_net_head_trad_2,“字符串”num2str (calc_net_head_trad_2));
%设置计算净酒吧
集(hand_calc_net_head_bar_trad_2,“字符串”num2str (calc_net_head_trad_2 * valg_design_rhoP * valg_design_gP * 10 ^ 5));
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%泵特性,计算损失
% H = Gulich,包括所有损失
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%设置计算理论
集(hand_calc_theoretical_head_trad_3,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 / valg_design_gP));
%设置计算理论在酒吧
集(hand_calc_theoretical_head_bar_trad_3,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 * valg_design_rhoP * 10 ^ 5));
[hand_axes_charac_trad_3, calc_net_head_trad_3 x_max_trad_3 y_max_trad_3]…
= draw_charac_curves_trad_3 (…
hand_axes_charac_trad_3,…
valg_design_flow,…
get (hand_calc_U2,“价值”),…
get (hand_calc_U1_merid,“价值”),…
get (hand_calc_CU2,“价值”),…
valg_design_D2,…
D_1,…
get (hand_value_b2,“价值”),…
valg_design_beta2,…
beta_1_merid,…
A_1,…
slip_coeff,…
eta_opt,…
block_2,…
valg_charac_shockloss_trad_3,…
valg_design_gP,…
valg_design_Z,…
valg_design_nethead,…
get (hand_calc_sp_speed,“价值”)…
);
%设置计算理论的头体积损失
集(hand_calc_theoretical_head_leakage_trad_3,“字符串”num2str (y_max_trad_3));
%设置计算理论在酒吧与体积损失
集(hand_calc_theoretical_head_leakage_bar_trad_3,“字符串”num2str (y_max_trad_3 * valg_design_rhoP * valg_design_gP * 10 ^ 5));
%设置计算净头
集(hand_calc_net_head_trad_3,“字符串”num2str (calc_net_head_trad_3));
%设置计算净酒吧
集(hand_calc_net_head_bar_trad_3,“字符串”num2str (calc_net_head_trad_3 * valg_design_rhoP * valg_design_gP * 10 ^ 5));
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%
%泵特点、经验
% H =…
%
% cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
%设置计算理论
集(hand_calc_theoretical_head_trad_4,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 / valg_design_gP));
%设置计算理论在酒吧
集(hand_calc_theoretical_head_bar_trad_4,…
“字符串”num2str (get (hand_calc_U2“价值”)^ 2 * valg_design_rhoP * 10 ^ 5));
[hand_axes_charac_trad_4, calc_net_head_trad_4 x_max_trad_4 y_max_trad_4]…
= draw_charac_curves_trad_4 (…
hand_axes_charac_trad_4,…
valg_design_flow,…
get (hand_calc_U2,“价值”),…
get (hand_calc_U1_merid,“价值”),…
get (hand_calc_CU2,“价值”),…
valg_design_D2,…
D_1,…
get (hand_value_b2,“价值”),…
valg_design_beta2,…
beta_1_merid,…
A_1,…
slip_coeff,…
eta_opt,…
block_2,…
valg_charac_shockloss_trad_4,…
valg_design_gP,…
valg_design_Z,…
valg_design_nethead,…
get (hand_calc_sp_speed,“价值”)…
);
%设置计算理论的头体积损失
集(hand_calc_theoretical_head_leakage_trad_4,“字符串”num2str (y_max_trad_4));
%设置计算理论在酒吧与体积损失
集(hand_calc_theoretical_head_leakage_bar_trad_4,“字符串”num2str (y_max_trad_4 * valg_design_rhoP * valg_design_gP * 10 ^ 5));
%设置计算净头
集(hand_calc_net_head_trad_4,“字符串”num2str (calc_net_head_trad_4));
%设置计算净酒吧
集(hand_calc_net_head_bar_trad_4,“字符串”num2str (calc_net_head_trad_4 * valg_design_rhoP * valg_design_gP * 10 ^ 5));
函数[h_axes, h_design_calc x_limit、y_limit pfleiderer_slip] =…
U_2 draw_charac_curves_trad_1 (h_axes Q_design, ~, ~, D_2, D_1, b_2, beta_2, ~, ~,滑,ω,fric,休克,g, Z, H)
集(gcf,“CurrentAxes”h_axes)
%计算中使用的参数
R_2 = 0.5 * D_2;
a₂= R_2 * 2 *π* b_2;
Q = 0;
step_size = Q_design / 100;
i = 1;
eta_v_opt = 0.02;
Q_Lopt = Q_design * eta_v_opt;
Ku = Q_Lopt /√(H);
ψ= 1.1 *(1 +信德(beta_2)) * (D_1 / D_2);
pfleiderer_slip = 2 * (psi / Z) * (1 / (1 - (D_1 / D_2) ^ 2));
W_U2temptemp = (Q / a₂) /罐内(beta_2);%应0(零)
C_U2temptemp = U_2-W_U2temptemp;
Htemptemp = (U_2 / g) * C_U2temptemp;
Q_L = Ku * sqrt (Htemptemp);
W_U2 = ((Q + Q_L) / a₂) /罐内(beta_2);
C_U2 = U_2-W_U2;
h_th_inf = (U_2 / g) * C_U2;% U_2 ^ 2 / g
h_th = h_th_inf / (1 + pfleiderer_slip);
h_fric = h_th - fric * (Q + Q_L) ^ 2;
dh_shock =冲击* ((Q + Q_L) - (Q_design + Q_L)) ^ 2;
h = h_fric - dh_shock;
Q = Q + step_size;
而(Q < (2 * Q_design))
W_U2temptemp = (Q / a₂) /罐内(beta_2);%应0(零)
C_U2temptemp = U_2-W_U2temptemp;
Htemptemp = (U_2 / g) * C_U2temptemp;
Q_L = Ku * sqrt (Htemptemp);
W_U2_temp = ((Q + Q_L) / a₂) /罐内(beta_2);
C_U2_temp = U_2-W_U2_temp;
h_th_inf_temp = (U_2 / g) * (C_U2_temp);
h_th_temp = h_th_inf_temp / (1 + pfleiderer_slip);
h_fric_temp = h_th_temp-fric * (Q + Q_L) ^ 2;
dh_shock_temp =冲击* ((Q + Q_L) - (Q_design + Q_L)) ^ 2;
h_temp = h_fric_temp - dh_shock_temp;
[W_U2] = [W_U2 W_U2_temp];
[h_th_inf] = [h_th_inf h_th_inf_temp];
[h_th] = [h_th h_th_temp];
[h_fric] = [h_fric h_fric_temp];
h_temp [h] = [h];
[dh_shock] = [dh_shock dh_shock_temp];
Q = Q + step_size;
我=我+ 1;
结束
Q = linspace (0, Q,长度(h_th_inf));
% (min_dh_shock min_index] = min (dh_shock);
h_design_calc = h (101);
x_limit = Q(结束);
y_limit = h_th_inf (1);
情节(h_axes Q h_th_inf,‘g’,…
问,h_th“米”,…
问,h_fric“b”,…
Q, h,“r”,…
问,dh_shock”:k”,…
[0 9.9 20.3 30.2 39.9 50.1 59.7 70.6 80.2 90.2 101]/ 3600,(75.22 - 74.03 72.67 - 71.16
68.40 66.02 61.99 56.71 - 50.54 43.65 - 35.62)/ 3,“k”,…
[Q_design Q_design],[0马克斯(h_th_inf)],“c:”)
传奇(“H_ {th \ infty}’,“H_ {th}’,“H_ {fric}’,“H”,“\ DeltaH_{冲击}’,“H_{测量}’)
集(h_axes,“XLim”[0 x_limit],“YLim”[0 y_limit]);
包含(“问(m ^ 3 / s)”);
ylabel (“H [m]”);
函数[h_axes, h_design_calc x_limit、y_limit pfleiderer_slip] =…
U_2 draw_charac_curves_trad_2 (h_axes Q_design, U_1, ~, D_2, D_1, b_2, beta_2, beta_1, A_1, ip,ω,fric,休克,g, Z, H)
集(gcf,“CurrentAxes”h_axes)
%计算中使用的参数
R_2 = 0.5 * D_2;
a₂= R_2 * 2 *π* b_2;
Q = 0;
step_size = Q_design / 100;
i = 1;
eta_v_opt = 0.02;
Q_Lopt = Q_design * eta_v_opt;
Ku = Q_Lopt /√(H);
ψ= 1.1 *(1 +信德(beta_2)) * (D_1 / D_2);
pfleiderer_slip = 2 * (psi / Z) * (1 / (1 - (D_1 / D_2) ^ 2));
W_U2temptemp = (Q / a₂) /罐内(beta_2);%应0(零)
C_U2temptemp = U_2-W_U2temptemp;
Htemptemp = (U_2 / g) * C_U2temptemp;
Q_L = Ku * sqrt (Htemptemp);
W_U2 = ((Q + Q_L) / a₂) /罐内(beta_2);
W_U1 = ((Q + Q_L) / A_1) /罐内(beta_1);
C_U2 = U_2-W_U2;
C_U1 = U_1-W_U1;
h_th_inf = (U_2 / g) * C_U2;% U_2 ^ 2 / g
h_th = h_th_inf / (1 + pfleiderer_slip);
h_fric = h_th - fric * (Q + Q_L) ^ 2;
dh_shock = abs ((U_1 * C_U1) / g);
h = h_fric - dh_shock;
Q = Q + step_size;
而(Q < (2 * Q_design))
W_U2temptemp = (Q / a₂) /罐内(beta_2);%应0(零)
C_U2temptemp = U_2-W_U2temptemp;
Htemptemp = (U_2 / g) * C_U2temptemp;
Q_L = Ku * sqrt (Htemptemp);
W_U2_temp = ((Q + Q_L) / a₂) /罐内(beta_2);
W_U1_temp = ((Q + Q_L) / A_1) /罐内(beta_1);
C_U2_temp = U_2-W_U2_temp;
C_U1_temp = U_1-W_U1_temp;
h_th_inf_temp = (U_2 / g) * (C_U2_temp);
h_th_temp = h_th_inf_temp / (1 + pfleiderer_slip);
h_fric_temp = h_th_temp-fric * (Q + Q_L) ^ 2;
dh_shock_temp = abs ((U_1 * C_U1_temp) / g);
h_temp = h_fric_temp - dh_shock_temp;
[W_U2] = [W_U2 W_U2_temp];
[h_th_inf] = [h_th_inf h_th_inf_temp];
[h_th] = [h_th h_th_temp];
[h_fric] = [h_fric h_fric_temp];
h_temp [h] = [h];
[dh_shock] = [dh_shock dh_shock_temp];
Q = Q + step_size;
我=我+ 1;
结束
Q = linspace (0, Q,长度(h_th_inf));
% (min_dh_shock min_index] = min (dh_shock);
h_design_calc = h (101);
x_limit = Q(结束);
y_limit = h_th_inf (1);
情节(h_axes Q h_th_inf,‘g’,…
问,h_th“米”,…
问,h_fric“b”,…
Q, h,“r”,…
问,dh_shock”:k”,…
[0 9.9 20.3 30.2 39.9 50.1 59.7 70.6 80.2 90.2 101]/ 3600,(75.22 - 74.03 72.67 - 71.16
68.40 66.02 61.99 56.71 - 50.54 43.65 - 35.62)/ 3,“k”,…
[Q_design Q_design],[0马克斯(h_th_inf)],“c:”)
传奇(“H_ {th \ infty}’,“H_ {th}’,“H_ {fric}’,“H”,“\ DeltaH_{冲击}’,“H_{测量}’)
集(h_axes,“XLim”[0 x_limit],“YLim”[0 y_limit]);
包含(“问(m ^ 3 / s)”);
ylabel (“H [m]”);
函数[h_axes, h_design_calc x_limit y_limit] =…
U_2 draw_charac_curves_trad_3 (h_axes Q_design, U_1, C_U2, D_2, D_1, b_2, beta_2, beta_1, A_1, slipcoeff,ω,fric,休克,g, Z, H, nq)
集(gcf,“CurrentAxes”h_axes)
%计算中使用的参数
R_2 = 0.5 * D_2;
a₂= R_2 * 2 *π* b_2;
%的一些尺寸
a1 = 0.018;%进气喉尺寸(估计)
b1 = 0.0324;%进气喉尺寸
a2 = 0.04;%出口喉咙维度(估计)
b2 = 0.013;%出口高度
Lsch = 0.14;%叶轮叶片长度(估计)
visc = 0.553 * 10 ^ 6;%运动粘度水50度
每股收益= 0.5 * 10 ^ 3;%的绝对粗糙度中度腐蚀碳钢
% (http://neutrium.net/fluid_flow/absolute-roughness/)
a3 = 0.0127;%扩散器进气喉宽
b3 = 0.0155;%扩散进喉咙高度
a4 = 0.0276;%扩散器出口喉部宽度
b4 = 0.0184;%扩散器出口喉咙高度
cp = 0.72;%扩散cp从图
Z_Le = 10;%扩压器叶片
zetaov = 1.5;%返回通道设计系数,在0.2和1.5之间。0.2是好的设计,1.0是糟糕的设计
Q = 0;
step_size = Q_design / 100;
i = 1;
Q_Lopt = 2 * 4.1 / (nq ^ 1.6) * Q_design;
Ku = Q_Lopt /√(H);
A_R = (a4 * b4) / (a3 * b3);
Dh = (2 * (a2 * b2 + a1 * b1)) / (a1 + a2 + b1 + b2);
W_U2temptemp = (Q / a₂) /罐内(beta_2);%应0(零)
C_U2temptemp = U_2-W_U2temptemp;
Htemptemp = (U_2 / g) * C_U2temptemp;
Q_L = Ku * sqrt (Htemptemp);
W_Q1 = (Q + Q_L) / (Z * a1 * b1);
W_U1 = ((Q + Q_L) / A_1) /罐内(beta_1);
C_U1 = U_1-W_U1;
C_M1 = (Q + Q_L) / A_1;
W_M1 =√C_M1 ^ 2 + (U_1-C_U1) ^ 2);
h_Las = 0.3 / (2 * g) * (W_M1-W_Q1) ^ 2;
W_av = 2 * (Q + Q_L) / (Z * (a2 * b2 + a1 * b1));
Re = W_av * Lsch / visc;
cf = 0.136 / ((log10 (0.2 * eps / Lsch + 12.5 / Re)) ^ 2.15);
cd = (cf + 0.0015) * (1.1 + 4 * b2 / D_2);
h_Laf = 2 * * Lsch / cd / g Dh * W_av ^ 2;
h_La = h_Las + h_Laf;%总叶轮损失
W_U2 = ((Q + Q_L) / a₂) /罐内(beta_2);
C_U2 = U_2-W_U2;
C_U2slip = slipcoeff * U_2-W_U2;
C_M2 = (Q + Q_L) / a₂;
C_2slip =√C_M2 ^ 2 + C_U2slip ^ 2);
C_Q3 = (Q + 0.5 * Q_L) / (Z_Le * a3 * b3);
h_23 = U_2 ^ 2 / (2 * g) * (cf + 0.0015) * (a3 / D_2 + b3 / D_2) *π^ 3 * ((C_M2 * b2) / (U_2 * D_2)) ^ 2 / (8 * (Z_Le * (a3 * b3) / (D_2))) ^ 3 * (1 + C_2slip / C_Q3) ^ 3;
h_Le = h_23 + C_Q3 ^ 2 / (2 * g) * (0.3 * (C_2slip / C_Q3-1) ^ 2 + 1-cp - (1-zetaov) / A_R ^ 2);%的总
扩散器损失
h_th_inf = (U_2 / g) * C_U2;% U_2 ^ 2 / g
h_th = (U_2 / g) * C_U2slip;
h = h_th - h_La-h_Le;
Q = Q + step_size;
而(Q < (2 * Q_design))
W_U2temptemp = (Q / a₂) /罐内(beta_2);%应0(零)
C_U2temptemp = U_2-W_U2temptemp;
Htemptemp = (U_2 / g) * C_U2temptemp;
Q_L_temp = Ku * sqrt (Htemptemp);
W_U2_temp = ((Q + Q_L_temp) / a₂) /罐内(beta_2);
W_Q1_temp = (Q + Q_L_temp) / (Z * a1 * b1);
W_U1_temp = ((Q + Q_L_temp) / A_1) /罐内(beta_1);
C_U2_temp = U_2-W_U2_temp;
C_U2slip_temp = slipcoeff * U_2-W_U2_temp;
C_U1_temp = U_1-W_U1_temp;
C_M1_temp = (Q + Q_L_temp) / A_1;
W_M1_temp =√C_M1_temp ^ 2 + (U_1-C_U1_temp) ^ 2);
h_Las_temp = 0.3 / (2 * g) * (W_M1_temp-W_Q1_temp) ^ 2;
W_av_temp = 2 * (Q + Q_L_temp) / (Z * (a2 * b2 + a1 * b1));
Re_temp = W_av_temp * Lsch / visc;
cf_temp = 0.136 / ((log10 (0.2 * eps / Lsch + 12.5 / Re_temp)) ^ 2.15);
cd_temp = (cf_temp + 0.0015) * (1.1 + 4 * b2 / D_2);
h_Laf_temp = 2 * cd_temp / g * Lsch / Dh * W_av_temp ^ 2;
h_La_temp = h_Las_temp + h_Laf_temp;%总叶轮损失
C_M2_temp = (Q + Q_L_temp) / a₂;
C_2slip_temp =√C_M2_temp ^ 2 + C_U2slip ^ 2);
C_Q3_temp = (Q + 0.5 * Q_L_temp) / (Z_Le * a3 * b3);
h_23_temp = U_2 ^ 2 / (2 * g) * (cf_temp + 0.0015) * (a3 / D_2 + b3 / D_2) *π^ 3 * ((C_M2_temp * b2) / (U_2 * D_2)) ^ 2 / (8 * (Z_Le * (a3 * b3) / (D_2))) ^ 3 * (1 + C_2slip_temp / C_Q3_temp) ^ 3;
h_Le_temp = h_23_temp + C_Q3_temp ^ 2 / (2 * g) * (0.3 * (C_2slip_temp / C_Q3_temp-1) ^ 2 + 1-cp - (1-zetaov) / A_R ^ 2);%总扩散损失
h_th_inf_temp = (U_2 / g) * (C_U2_temp);
h_th_temp = (U_2 / g) * (C_U2slip_temp);
h_temp = h_th_temp - h_La_temp-h_Le_temp;
[W_U2] = [W_U2 W_U2_temp];
[h_th_inf] = [h_th_inf h_th_inf_temp];
[h_th] = [h_th h_th_temp];
[h_Le] = [h_Le h_Le_temp];
[h_La] = [h_La h_La_temp];
h_temp [h] = [h];
Q = Q + step_size;
我=我+ 1;
结束
Q = linspace (0, Q,长度(h_th_inf));
% (min_dh_shock min_index] = min (h_Le);
h_design_calc = h (101);
x_limit = Q(结束);
y_limit = h_th_inf (1);
情节(h_axes Q h_th_inf,‘g’,…
问,h_th“米”,…
问,h_Le“b”,…
Q, h,“r”,…
问,h_La”:k”,…
[0 9.9 20.3 30.2 39.9 50.1 59.7 70.6 80.2 90.2 101]/ 3600,(75.22 - 74.03 72.67 - 71.16
68.40 66.02 61.99 56.71 - 50.54 43.65 - 35.62)/ 3,“k”,…
[Q_design Q_design],[0马克斯(h_th_inf)],“c:”)
传奇(“H_ {th \ infty}’,“H_ {th}’,“\ DeltaH_{扩散}’,“H”,“\ DeltaH_{叶轮}’,“H_{测量}’)
集(h_axes,“XLim”[0 x_limit],“YLim”[0 y_limit]);
包含(“问(m ^ 3 / s)”);
ylabel (“H [m]”);
函数[h_axes, h_design_calc x_limit y_limit] =…
U_2 draw_charac_curves_trad_4 (h_axes Q_design, U_1, C_U2, D_2, D_1, b_2, beta_2, beta_1, A_1, slipcoeff_opt, eta_opt, block_2,休克,g, Z, H, nq)
集(gcf,“CurrentAxes”h_axes)
%计算中使用的参数
R_2 = 0.5 * D_2;
a₂= R_2 * 2 *π* b_2;
Q = 0;%
step_size = (Q_design-Q) / 100;
i = 1;
Q_Lopt = 2 * 4.1 / (nq ^ 1.6) * Q_design;
Ku = Q_Lopt /√(H);
W_U2temptemp = (Q / a₂) /罐内(beta_2);%应0(零)
C_U2temptemp = U_2-W_U2temptemp;
Htemptemp = (U_2 / g) * C_U2temptemp;
Q_L = Ku * sqrt (Htemptemp);
slipcoeff = slipcoeff_opt * (1.3746 - -0.2483 * ((Q + Q_L) / (Q_design + Q_L)) ^ 3 + 0.8639 * ((Q + Q_L) / (Q_design + Q_L)) ^ 2 - 0.9919 * ((Q + Q_L) / (Q_design + Q_L)));
W_U2 = ((Q + Q_L) / a₂) /罐内(beta_2);
C_U2 = U_2-W_U2;
C_U2slip = slipcoeff * U_2-W_U2;
h_th_inf = (U_2 / g) * C_U2;% U_2 ^ 2 / g
h_th = (U_2 / g) * C_U2slip;
η= eta_opt * (1 - 0.6 * ((Q + Q_L) / -0.9 (Q_design + Q_L)) ^ 2 - 0.25 * ((Q + Q_L) / -0.9 (Q_design + Q_L)) ^ 3);
h =η* U_2 / g * (slipcoeff * U_2-W_U2 * block_2);%如果non-rotational
入口
Q = Q + step_size;
而(Q < (2 * Q_design))
W_U2temptemp = (Q / a₂) /罐内(beta_2);%应0(零)
C_U2temptemp = U_2-W_U2temptemp;
Htemptemp = (U_2 / g) * C_U2temptemp;
Q_L_temp = Ku * sqrt (Htemptemp);
slipcoeff_temp = slipcoeff_opt * (1.3746 - -0.2483 * ((Q + Q_L_temp) / (Q_design + Q_L_temp)) ^ 3 + 0.8639 * ((Q + Q_L_temp) / (Q_design + Q_L_temp)) ^ 2 * -0.9919 ((Q + Q_L_temp) / (Q_design + Q_L_temp)));
W_U2_temp = ((Q + Q_L_temp) / a₂) /罐内(beta_2);
C_U2_temp = U_2-W_U2_temp;
C_U2slip_temp = slipcoeff_temp * U_2-W_U2_temp;
h_th_inf_temp = (U_2 / g) * (C_U2_temp);
h_th_temp = (U_2 / g) * (C_U2slip_temp);
eta_temp = eta_opt * (1 - 0.6 * ((Q + Q_L_temp) / -0.9 (Q_design + Q_L_temp)) ^ 2 - 0.25 * ((Q + Q_L_temp) / -0.9 (Q_design + Q_L_temp)) ^ 3);
h_temp = eta_temp * U_2 / g * (slipcoeff_temp * U_2-W_U2_temp * block_2);
[W_U2] = [W_U2 W_U2_temp];
[h_th_inf] = [h_th_inf h_th_inf_temp];
[h_th] = [h_th h_th_temp];
h_temp [h] = [h];
Q = Q + step_size;
我=我+ 1;
结束
Q = linspace (0, Q,长度(h_th_inf));
% (min_dh_shock min_index] = min (h);
h_design_calc = h (101);
x_limit = Q(结束);
y_limit = h_th_inf (1);
情节(h_axes Q h_th_inf,‘g’,…
问,h_th“米”,…
Q, h,“r”,…
[0 9.9 20.3 30.2 39.9 50.1 59.7 70.6 80.2 90.2 101]/ 3600,(75.22 - 74.03 72.67 - 71.16
68.40 66.02 61.99 56.71 - 50.54 43.65 - 35.62)/ 3,“k”,…
[Q_design Q_design],[0马克斯(h_th_inf)],“c:”)
传奇(“H_ {th \ infty}’,“H_ {th}’,“H”,“H_{测量}’)
集(h_axes,“XLim”[0 x_limit],“YLim”[0 y_limit]);
包含(“问(m ^ 3 / s)”);
ylabel (“H [m]”);
结束
结束
结束
结束

在回答这个问题。

答案(0)

类别

找到更多的在建立深层神经网络在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!