8
ifTauP = 1.0 then Rk := 0.0292 + 0.85 * Rn;
if(TauP <> 0.0) and (TauP <> 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;
ifSodSuA[NNit] = 2 then begin L1 := 0.039; L2 := 0.47 end;
Cb := 0.5959 + 0.0312 * r_(Bet, 2.1) - 0.184 sqr(Bet4) +
L1 Bet4 / (1.0 - Bet4) - 0.0337 * L2 * r_(Bet, 3)
end;
if NSuzA[NNit] = 1 then Cb := 0.99 - 0.2262 * r_(Bet, 4.1);
{ Для сопла Вентури Cb = С, так как KRe = 1 }
ifNSuzA[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* Ro) / Roc;
Re := 4.e6 * Qcb * Roc / 3.6 / 3.141592653 / Mu / 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 * г (Re, 1.15) / 923.9 / sqr(Bet) / (33.0 *
г_(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
elseARe := 1.0;
KSh:= 1.0 + Bet4* RO-ARe