32
if EL > 1.0 then begin I ERR := 3; goto 1 end;
EE:=0.0;
if NSUZA[NNit] = 0 then begin
EDM:=DT*sqrt(0.1*DP*(0.681-0.651*BET)/DAC[2]);
if DAC[1] < EDM then begin
A:=BET*(13.5-15.5*BET); B:=117.0-106.0*r_(BET,1.9);
EE:=0.1*DP*sqr(DT)*(A*DT/DAC[1]-B)/DAC[3]/sqr(DAC[1])
end; {endif}
end; {endif}
EX:=0.0;
EXMIN:==0.0025*DT/(0.1+2.3*sqr(BET)*sqr(BET));
EXMAX:=2.0*EXMIN;
if TAC[1] > EXMAX then begin IERR := 4; goto 1 end;
if TAC[1] > EXMIN then EX:=0.3;
EH:=0.0; HD:=TAC[3]/DT;
if HD > 0.003 then begin
HDC:=0.002*(TAC[2]/DT+0.4)/(0.1+2.3*sqr(BET)*sqr(BET));
if (HD<=HDC) and (HD<=0.05) then EH:=0.2
else begin
IERE.:=5; goto 1
end; {endif}
end; {endif}
ESH:=100.0*(KSH-1.0); EK:=100.0*(KK-1.0);
ERC:=sqrt(sqr(EC0+EL+EE+EH+EX)+sqr(ESH)+sqr(EK));
1:
end;{ С }
{---------------------------------------------------------------------------------------------}
{ Расчет методической составляющей погрешности определения коэффициента расширения }
Procedure EPSI;
begin { EPSI }
EEPS:=0.0;
if NSUBA[NNIT]<>23 then begin
case NSUZA[NNIT] of
0: if BET <= 0.75 then EEPS := 4.0-DP/P
else EEPS := 8.0-DP/P;
1: EEPS := 2.0*DP/P;
else
EEPS := DP*(4.0+100.0*r_(BET, 8.0))/P;
end;
end; {endif}
end; { EPSI }
{----------------------------------------------------------------------------------}
{Расчет методических погрешностей определения плотности, показателя
адиабаты и динамической вязкости}
Procedure EPH;
var
N:byte;
const
EPH3:array (1..2,1..3,1..3) ofreal=(((0.3,0.4,0.2),
(0.9,1 -0,0.6),(2.0,3.0,2.0)),((0.6,1.3,0.4),(0.6,1.1,0.6),
(2.0,3.0,2.0)));