40
L1 * Bet4 / (1.0 - Bet4) - 0.0337 * L2 * r_(Bet, 3)
end;
if NsuzA[NNIT] = 1 then Cb := 0.99 - 0.2262 * rJBet, 4.1);
{Для сопла Вентури Cb = С, так как KRe = 1}
if NSuzA[NNIT] = 2 then Cb := 0.9858 - 0.196 * r_(Bet, 4.5);
{Для труб Вентури Cb = С = const, так как KRe = 1 и Re > 2.е5}
case NSuzA[NNIT] of
3: Cb = 0.984;
4: Cb = 0.995;
5: Cb = 0.985
end;
{Расчет коэффициента расхода сужающего устройства и расхода при числе Рейнольдса, стремящемся к бесконечности}
Alfa := Cb * Ее;
Qcb := 0.039986 Alfa * KSh * Kk * Eps * sqr(Dd) *
sqrt(1.e3 * Dp * ROX) / Roc;
Re := 4.e6 * Qcb * Roc / 3.6 / 3.141592653 / MUX / Dt;
{Расчет поправочного коэффициента на число Рейнольдса}
case NSuzA[NNIT] of
0: KRe = 1.0 + 1.426 / (1.0 + Cb * r_(Re, 0.75) / 64.28 /
r_Bet. 2.5));
1: KRe .= 1.0 - 0.86 / (1.0 + Cb * r_(Re, 1.15) / 923.9 /
sqr(Bet) / (33.0 * r_(Bet, 2.15) - 17.5));
2: KRe := 1.0;
3: KRe := 1.0;
4: KRe := 1.0;
5: KRe := 1.0;
end;
{Определение действительного значения числа Рейнольдса}
Re := Re * KRe;
{Расчет поправочного коэффициента на шероховатость внутренней поверхности
измерительного трубопровода с учетом числа Рейнольдса для всех сужающих
устройств, кроме труб Вентури}
if (NSuzA[NNit] <= 2) and (RSh <> 0.0) then begin
Qcb := Qcb / KSh;
if Re > 1 .e4 then begin
if Re < 1.e6 then ARe := 1.0 - sqr(Ln(Re) / 2.3026 - 6.0) / 4.0
else ARe := 1.0;
KSh := 1.0 + Bet4* R0 * ARe
end;
if (Re <= 1.e4) or (KSh <= 1.0005) then KSh := 1.0;
Qcb := Qcb KSh
end;
{Определение коэффициента истечения для труб Вентури в зависимости от числа
Рейнольдса}
if (Re < 2.e5) and (NSuzA[NNIT] > 2) then
case NSuzA[NNIT] of
3: begin
if Re <= 6.e4 then Cb := 0.957;
if (Re > 6.e4) and (Re <= 1.e5) then Cb := 0.966;
if (Re > 1.e5) and (Re <= 1.5e5) then Cb := 0.976;
if Re > 1.5e5 then Cb := 0.982