39
if NSuzA[NNIT] = 0 then
Eps := 1.0 - (0.41 + 0.35 * Bet4) * Dp / P / KAPPAX;
ifNSuzA[NNIT] <> 0 then begin
Psi := 1.0 - Dp / P;
Eps := KAPPAX * r_(Psi, 2.0 / KAPPAX) / (KAPPAX - 1.0) *
(1.0 - Bet4) / (1.0 - Bet4 r_(Psi, 2.0 / KAPPAX))*
(1.0 - r_(Psi, (KAPPAX - 1.0) / KAPPAX)) /
(1.0 - Psi);
Eps := sqrt(Eps)
end;
end;
{Расчет поправочного коэффициента на шероховатость внутренней поверхности
измерительного трубопровода без учета числа Рейнольдса}
KSh .= 1.0;
if (NSuzA[NNit] <= 2) and (RSh <> 0.0) then begin
ARe := 0.5; Rd := RSh / Dt; Log := Ln(Rd * 1.e4) / 2.3026;
if NSuzA(NNit) = 0 then begin
if Log <= (1.0 / 10.0 / Bet4 + 8.0) / 14.0 then begin
R0 := 0.0; goto 1;
end;
R0 := 0.07 * Log - 0.04
end
else begin
if Log <= (1.0 / 10.0 / Bet4 + 5.0) / 9.0 then begin
R0 := 0.0; goto 1;
end;
R0 := 0.045 Log - 0.025
end;
KSh :== 1.0 + Bet4 * R0 ARe
end;
1:;
{Расчет поправочного коэффициента на притупление входной кромки отверстия диафрагмы}
Kk := 1.0;
if (NSuzA[NNIT] = 0) and (DdO < 125.0) then begin
if TauP = 0.0 then Rk := Rn;
if TauP = 1.0 then Rk := 0.0292 + 0.85 * Rn;
if (TauP <> 0.0) and (TauP 0 1.0) then
Rk := 0.195 - (0.195 - Rn) * (1.0 - Exp(-TauP / 3.0)) *
3.0 / TauP;
Kk := 1.0547 - 0.0575 * Exp(-149.0 * Rk / Dd)
end;
{Расчет коэффициента истечения при числе Рейнольдса, стремящемся к бесконечности}
if NSuzA[NNit] = 0 then begin
L1 := 0.0; L2 := 0.0;
if SodSuA[NNit] = 1 then begin
L1 := 25.4 / Dt; L2 := L1;
if L1 >= 0.4333 then L1 := 0.039 else L1 := 0.09 * LI
end;
if SodSuA[NNit] = 2 then begin L1 := 0.039; L2 := 0.47 end;
Cb := 0.5959 + 0.0312 * r_(Bet, 2.1) - 0.184 * sqr(Bet4) +