7
Hs[1] := HsSubs1[NSubA[NNit]]; Hs[2] := HsSubs2[NSubA[NNIT]]
end;
readln(Fl, Ro, Kappa, Mu);
close(Fl); erase(Fl);
if (NSubA[NNit] = 0) and (VarRoA[NNit] = 1) then Ro := Rot;
KZ := P * Roc * 293.15 / Ro / (Т + 273.15) / 1.01325;
if NSubA[NNit] = 0 then str(Roc:6:4, RocStr);
{Расчет: 1) диаметров сужающего устройства и измерительного трубопровода при рабочей температуре; 2) относительного диаметра; 3) коэффициента скорости входа}
Dd := (1.0 + AlfaSU * (Т - 20.0)) * Dd20;
Dt := (1.0 + AlfaT * (Т - 20.0)) * Dt20;
Bet := Dd / Dt; Bet4 := sqr(Bet) * sqr(Bet);
Ec := 1.0/sqrt1.0 - Bet4);
{Расчет коэффициента расширения}
Eps := 1.0;
if NSubA[NNIT] <> 23 then begin
if NSuzA[NNit] = 0 then
Eps := 1.0 - (0.41 + 0.35 * Bet4) * Dp / P / Kappa;
if NSuzA[NNit] <> 0 then begin
Psi := 1.0 - Dp/ P;
Eps := Kappa * r_(Psi, 2.0 / Kappa) / (Kappa - 1.0) *
(1.0 - Bet4) / (1.0 - Bet4 * r_(Psi, 2.0 / Kappa))*
(1.0 - r_(Psi, (Kappa - 1.0) / Kappa)) /
(1.0 - Psi);
Eps := sqn(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;
ifNSuzA[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 (Dd20 < 125.0) then begin
ifTauP = 0.0 then Rk := Rn;