ГОСТ 8.563.3-97 Государственная система обеспечения единства измерений. Измерение расхода и количества жидкостей и газов методом переменного перепада давления. Процедура и модуль расчетов. Программное обеспечение

ГОСТ 8.563.3-97

Группа Т86

МЕЖГОСУДАРСТВЕННЫЙ СТАНДАРТ

ГОСУДАРСТВЕННАЯ СИСТЕМА ОБЕСПЕЧЕНИЯ ЕДИНСТВА ИЗМЕРЕНИЙ

Измерение расхода и количества жидкостей
и газов методом переменного перепада давления

ПРОЦЕДУРА И МОДУЛЬ РАСЧЕТОВ. ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ

State system for ensuring the uniformity of measurements.
Measurement of liquids and gases flow rate and quantity by differential
pressure method. Procedure and module of calculations. Software

МКС 17.020
ОКСТУ 0008

Дата введения 1999-01-01

Предисловие

1 РАЗРАБОТАН фирмой "Газприборавтоматика" РАО "Газпром", Всероссийским научно-исследовательским центром стандартизации, информации и сертификации сырья, материалов и веществ (ВНИЦСМВ) Госстандарта России

ВНЕСЕН Госстандартом России

2 ПРИНЯТ Межгосударственным Советом по стандартизации, метрологии и сертификации (протокол № 11-97 от 25 апреля 1997 г.)

За принятие проголосовали:

3 ВВЕДЕН В ДЕЙСТВИЕ непосредственно в качестве государственного стандарта Российской Федерации с 1 января 1999 г. постановлением Государственного комитета Российской Федерации по стандартизации, метрологии и сертификации от 11 декабря 1997 г. № 410 межгосударственный стандарт ГОСТ 8.563.3-97

4 ВЗАМЕН ГОСТ 23720-79, ГОСТ 26969-86, РД 50-213-80, МИ 2204-92, МИ 2346-95

1 ОБЛАСТЬ ПРИМЕНЕНИЯ

Настоящий стандарт устанавливает программное обеспечение расчета расхода и количества жидкостей и газов (далее - среда), а также расчета погрешностей определения расхода и количества сред.

Листинги программ, приведенные в настоящем стандарте, разработаны согласно требованиям, изложенным в ГОСТ 8.563.1 и ГОСТ 8.563.2.

2 НОРМАТИВНЫЕ ССЫЛКИ

В настоящем стандарте использованы ссылки на следующие стандарты:

ГОСТ 8.563.1-97 ГСИ. Измерение расхода и количества жидкостей и газов методом переменного перепада давления. Диафрагмы, сопла ИСА 1932 и трубы Вентури, установленные в заполненных трубопроводах круглого сечения. Технические условия

ГОСТ 8.563.2-97 ГСИ. Измерение расхода и количества жидкостей и газов методом переменного перепада давления. Методика выполнения измерений с помощью сужающих устройств

ГОСТ 30319.1-96 Газ природный. Методы расчета физических свойств. Определение физических свойств природного газа, его компонентов и продуктов его переработки

ГОСТ 30319.2-96 Газ природный. Методы расчета физических свойств. Определение коэффициента сжимаемости

ГОСТ 30319.3-96 Газ природный. Методы расчета физических свойств. Определение физических свойств по уравнению состояния

3 ПРОЦЕДУРА РАСЧЕТА РАСХОДА И КОЛИЧЕСТВА ЖИДКОСТЕЙ И ГАЗОВ

Процедура расчета расхода и количества жидкостей и газов (далее - процедура) написана на алгоритмическом языке ТУРБО ПАСКАЛЬ 7.0. Обращение к процедуре - QCALC.

3.1 Исходные данные

Исходные данные передаются в процедуру QCALC в виде глобальных параметров; для работы процедуры необходимо использовать модули Dos и Crt.

3.1.1 NNit - порядковый номер измерительного трубопровода.

3.1.2 NSubA[NNit] - номер среды (таблица 1).

Таблица 1

Наименование и номер среды

3.1.3 NMethKA[NNit] - номер метода расчета коэффициента сжимаемости природного газа по ГОСТ 30319.2 (таблица 2).

Таблица 2

Наименование и номер метода расчета

3.1.4 NSuzA[NNit] - номер типа сужающего устройства (таблица 3).

Таблица 3

Наименование и номер типа сужающего устройства

3.1.5 YR - массив концентраций 16 компонентов природного газа, молярные доли (таблица 4).

Таблица 4

Наименование компонента и его порядковый номер в массиве

3.1.6 Параметры среды

Roc - плотность природного газа при стандартных условиях, кг/м.

(стандартные условия: = 1,01325 бар, = 20 °С)*.

______________

* Условные обозначения параметров при стандартных условиях - по ГОСТ 8.563.1 и ГОСТ 8.563.2.

Ya - концентрация азота в природном газе, молярные доли.

Yy - концентрация диоксида углерода в природном газе, молярные доли.

Р - давление среды, бар.

Т - температура среды, °С.

3.1.7 Характеристики сужающего устройства

AlfaSU - температурный коэффициент линейного расширения материала сужающего устройства, 1/°С.

Dd20 - диаметр отверстия сужающего устройства при 20 °С, мм.

Rn - начальный радиус закругления входной кромки диафрагмы, мм.

ТаиР - период поверки диафрагмы, год.

Dp - перепад давления на сужающем устройстве, бар.

3.1.8 SodSuA[NNit] - номер способа отбора давления на диафрагме:

1) 0 - угловой;

2) 1 - фланцевый;

3) 2 - трехрадиусный.

3.1.9 Характеристики измерительного трубопровода

AlfaT - температурный коэффициент линейного расширения материала измерительного трубопровода, 1/°С.

Dt20 - внутренний диаметр измерительного трубопровода при 20 °С, мм.

Rsh - эквивалентная шероховатость материала измерительного трубопровода, мм.

3.1.10 TauAv - время, за которое определяют количество среды, ч.

3.1.11 VarRoA[NNit] - номер способа определения плотности природного газа при рабочих условиях:

1) 0 - плотность рассчитывают по ГОСТ 30319.1;

2) 1 - плотность измеряют. Rot - измеренное значение плотности.

3.2 Выходные данные

Ro - плотность среды при рабочих условиях, кг/ м.

KZ - коэффициент сжимаемости среды.

Kappa - показатель адиабаты среды.

Mu - динамическая вязкость среды, мкПа·с.

Qc - объемный расход среды, приведенный к стандартным условиям, м/ч.

Vc - объем среды (количество среды, выраженное в кубических метрах), приведенный к стандартным условиям, м.

Vm - масса среды (количество среды, выраженное в тоннах), т.

Hs[1] - высшая удельная теплота сгорания среды, МДж/ м.

Hs[2] - низшая удельная теплота сгорания среды, МДж/м.

3.3 Листинг процедуры расчета расхода и количества жидкостей и газов

В нижеприведенной процедуре вызываются две внешние программы:

1) TpNg.exe - расчет теплофизических свойств природного газа в соответствии с требованиями ГОСТ 30319.2 и ГОСТ 30319.3;

2) TpSubs.exe - расчет теплофизических свойств компонентов природного газа и продуктов его переработки в соответствии с требованиями ГОСТ 30319.1.

Типы используемых переменных: Fl: text; NNit: byte; Dd, Dt, Dd20, Dt20, RSh, Rn, TauP, AlfaT, AlfaSU, Roc, Ya, Yy, Dp, P, T, Ro, Rot, Mu, Kappa, KZ, Eps, KSh, Kk, Cb, KRe, Re, Vc, Vm, TauAv: real; NSubA, NSuzA, SodSuA, NMethKA, VarRoA: array[1..30] of byte; YR: array[1..16] of real; Hs: array[1..2] of real;

Procedure Qcalc;

var

I,IBeg,IFin: byte; Code: integer;

Bet,Bet4,Ec,Rd,Psi,Rk,L1,L2,Alfa,Qcb,ARe,R0,KCb,Qc,Vcv,Log: real;

HsS: string [10];

label

1,3;

const

RocSubs: array[1..31] of real = (0.6682,1.2601,1.8641,2.488,

HsSubs1: аrrау[1..31] of real = (37.12,65.43,93.85,122.8,123.6,0.0,

HsSubs2: array[1..31] of real = (33.43,59.87,86.37,113.4,114.1,0.0,

CalcTpNg = 'TpNg.exe'; CalcTpSubs = 'TpSubs.exe';

assign(Fl,'IRD'); rewrite(Fl);

if NSubA[NNit]=0 then begin

end

writeln(Fl,P:14,Bl,T:14);

close(Fl);

TextColor(7);

gotoxy(19,9);

writeln('|------------------------------------|');

gotoxy(19,10);

writeln('| |');

gotoxy(19,11);

write('| ');

TextColor(135);

write('Ж Д И Т Е');

TextColor(7);

writeln(` |');

gotoxy(19,12);

writeln('| |');

gotoxy(19,13);

writeln(`| В Ы П О Л Н Я Е Т С Я Р А С Ч Е Т |');

gotoxy(19,14);

writeln('| |');

gotoxy(19,15);

writeln('|____________________________________|');

if NSubA[NNit]=0 then begin

gotoxy(21,12);

swapvectors; exec(CalcTpNg,CalcTpNg); swapvectors;

TextColor(7);

gotoxy(19,9);

writeln('--------------------------------------');

gotoxy(19,10);

writeln('| |');

gotoxy(19,11);

write('| ' );

TextColor(135);

write('Ж Д И Т Е');

TextColor(7);

writeln(' |');

gotoxy(19,12);

writeln('| |');

gotoxy(19,13);

writeln('| В Ы П О Л Н Я Е Т С Я Р А С Ч Е Т |');

gotoxy(19,14);

writeln('| |');

gotoxy(19,15);

writeln('|_____________________________________|');

end

end;

assign(Fl,'IRD'); reset(Fl);

if(NSubA[NNit]=0) and (NMethKA[NNit] >= 2) then

if NSubA[NNit]=0 then begin

end

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);

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/sqrt(1.0-Bet4);

{ Расчет коэффициента расширения }

Eps:= 1.0;

if NSubA[NNit] <> 23 then begin

if NSuzA[NNit] <> 0 then begin

end;

KSh:= 1.0;

if (NSuzA[NNit] <= 2) and (RSh <> 0.0) then begin

end;

KSh:= 1.0+Bet4*R0*ARe

end;

Kk:= 1.0;

if (NSuzA[NNit]=0) and (Dd20<125.0) then begin

end;

if NSuzA[NNit]=0 then begin

L1:= 0.0; L2:= 0.0;

if SodSuA[NNit]=1 then begin

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)+

end;

if NSuzA[NNit]=1 then Cb:= 0.99-0.2262*r_(Bet,4.1);

if NSuzA[NNit]=2 then Cb:= 0.9858-0.196*r_(Bet,4.5);

case NSuzA[NNit] of

end;

Alfa:= Cb*Еc;

Qcb:= 0.039986*Alfa*KSh*Kk*Eps*sqr(Dd)*

case NSuzA[NNit] of

end;

Re:= Re*KRe;

if (NSuzA[NNit] <= 2) and (RSh <> 0.0) then begin

end;

if (Re < 2.е5) and (NSuzA[NNit] > 2) then

KCb:= 1.0;

if (Re < 2.е5) and (NSuzA[NNit] > 2) then

Qc:= Qcb*KRe*KCb; Vcv:= Qc*TauAv; Vc:= Vcv;

end; { QCalc }

function r_(A,R: real): real;

function r_; begin r_:= exp(R*ln(A)) end;

4 МОДУЛЬ РАСЧЕТА ПОГРЕШНОСТИ ОПРЕДЕЛЕНИЯ РАСХОДА
И КОЛИЧЕСТВА ЖИДКОСТЕЙ И ГАЗОВ

Все процедуры и функции модуля расчета погрешности определения расхода и количества жидкостей и газов (далее - модуль ERRQSNX) написаны на алгоритмическом языке ТУРБО ПАСКАЛЬ 7.0. Обращение к модулю ERRQSNX осуществляется с помощью оператора вызова основной процедуры QS.

4.1 Исходные данные

Исходные данные передаются в процедуры модуля ERRQSNX в виде глобальных параметров; для работы процедур модуля необходимо использовать модули Dos и Crt.

4.1.1 NSubA[NNit], NSuzA[NNit], NMethKA[NNit], SodSuA[NNit], VarRoA[NNit], Dd20, Dt20, AlfaT, AlfaSU, RSh, Rn, TauP, TauAv - cм. 3.1.

4.1.2 VarPA[NNit] - номер варианта измерения давления:

1) 0 - измеряют избыточное давление;

2) 1 - измеряют абсолютное давление.

4.1.3 NVarYA[NNit] - номер варианта задания концентраций компонентов природного газа:

1) 0 - задание полного компонентного состава (см. 3.1.5);

2) 1 - задание концентраций азота и диоксида углерода.

4.1.4 Характеристики гидравлических сопротивлений

AmountR - количество гидравлических сопротивлений до сужающего устройства (от 2 до 9).

NRA[J] - массив номеров, характеризующих тип гидравлических сопротивлений, которые расположены до сужающего устройства (1 J 9) (таблица 5).

Таблица 5

Наименование (тип) гидравлических сопротивлений и их номер

LRA[J] - отношение расстояния от сужающего устройства до гидравлического сопротивления (включая гидравлическое сопротивление за сужающим устройством) к диаметру измерительного трубопровода (1 J 10).

4.1.5 Характеристики сужающего устройства

DAC[1] - толщина материала диафрагмы, мм.

DAC[2] - предел текучести материала диафрагмы при 20 °С, МПа.

DAC[3] - модуль Юнга материала диафрагмы при 20 °С, МПа.

4.1.6 Характеристики измерительного трубопровода

ТАС[1] - эксцентриситет между осями измерительного трубопровода и сужающего устройства, мм.

ТАС[2] - расстояние от уступа до отверстия для отбора давления, мм.

ТАС[3] - высота уступа между двумя секциями измерительного трубопровода, мм.

4.1.7 Характеристики измерительного комплекса (измерение перепада давления)

SIZMDPD - номер варианта состава измерительного комплекса:

1) 0 - средство измерений с линейной функцией преобразования;

2) 1 - средство измерений с квадратичной функцией преобразования;

3) 2 - измерительный преобразователь и средство измерений с линейной функцией преобразования;

4) 3 - измерительный преобразователь и средство измерений с квадратичной функцией преобразования;

5) 4 - измерительный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования;

6) 5 - измерительный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с квадратичной функцией преобразования;

7) 6 - измерительный преобразователь, средство измерений с квадратичной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования;

8) 7 - измерительный преобразователь, корнеизвлекающий преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования.

EDPA[J] - массив погрешностей измерительного комплекса (таблицы 6, 7, 8).

Таблица 6

Наименование погрешности средства измерений или измерительного преобразователя и ее номер (J) в массиве EDPA[J]

Таблица 7

Наименование погрешности средства измерений или корнеизвлекающего преобразователя и ее номер (J) в массиве EDPA[J]

Таблица 8

Наименование погрешности средства измерений и ее номер (J) в массиве EDPA[J]

EDPA[33] - систематическая составляющая погрешности планиметра (или интегратора), %.

EDPA[34] - случайная составляющая погрешности планиметра (или интегратора), %.

EDPA[35] - верхний предел измерения перепада давления, бар.

4.1.8 Характеристики измерительного комплекса (измерение давления)

SIZMPD - номер варианта состава измерительного комплекса:

1) 0 - средство измерений с линейной функцией преобразования;

2) 1 - средство измерений с квадратичной функцией преобразования;

3) 2 - измерительный преобразователь и средство измерений с линейной функцией преобразования;

4) 3 - измерительный преобразователь и средство измерений с квадратичной функцией преобразования;

5) 4 - измерительный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования;

6) 5 - измерительный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с квадратичной функцией преобразования;

7) 6 - измерительный преобразователь, средство измерений с квадратичной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования.

EPA[J] - массив погрешностей измерительного комплекса (таблицы 9, 10).

Таблица 9

Наименование погрешности средства измерений или измерительного преобразователя и ее номер (J) в массиве EPA[J]

Таблица 10

Наименование погрешности средства измерений и ее номер (J) в массиве EPA[J]

ЕРА[22] - систематическая составляющая погрешности планиметра (или интегратора), %.

ЕРА[23] - случайная составляющая погрешности планиметра (или интегратора), %.

ЕРА[24] - основная погрешность барометра, %.

ЕРА[25] - верхний предел измерения атмосферного давления, бар.

ЕРА[26] - верхний предел измерения давления, бар.

4.1.9 Характеристики измерительного комплекса (измерение температуры)

SIZMTD - номер варианта состава измерительного комплекса:

1) 0 - термометр;

2) 1 - измерительный преобразователь и средство измерений с линейной функцией преобразования;

3) 2 - измерительный преобразователь и средство измерений с квадратичной функцией преобразования;

4) 3 - измерительный преобразователь, средство измерений с линеиной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования;

5) 4 - измерительный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с квадратичной функцией преобразования;

6) 5 - измерительный преобразователь, средство измерений с квадратичной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования;

7) 6 - измерительный преобразователь, вторичный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования.

ЕТА[J] - массив погрешностей измерительного комплекса (таблицы 11, 12, 13).

ЕТА[0] - абсолютная погрешность термометра, °С.

Таблица 11

Наименование погрешности измерительного преобразователя и ее номер (J) в массиве ETA[J]

Таблица 12

Наименование погрешности средства измерений или вторичного преобразователя и ее номер (J) в массиве ETA(J)

Таблица 13

Наименование погрешности средства измерений и ее номер (J) в массиве ETA[J]

ЕТА[32] - систематическая составляющая погрешности планиметра (или интегратора), %.

ЕТА[33] - случайная составляющая погрешности планиметра (или интегратора), %.

ЕТА[34] - верхний предел измерения температуры, °С.

ЕТА[35] - нижний предел измерения температуры, °С.

4.1.10 Характеристики измерительного комплекса (измерение плотности природного газа при стандартных условиях)

EROSA[J] - массив погрешностей измерительного комплекса (таблица 14).

Таблица 14

Наименование погрешности измерительного преобразователя и ее номер (J) в массиве EROSA[J]

EROSA[6] - относительная погрешность средства измерений, %.

EROSA[7] - EROSA[22] - относительные погрешности измерения концентраций компонентов природного газа (см. 3.1.5), %.

4.1.11 ErrTau - погрешность определения интервала времени TauAV, %.

4.1.12 DHD - абсолютная погрешность хода приводного механизма диаграммы, мин.

4.1.13 Параметры эксплуатации (измеряемые за определенный промежуток времени - месяц, год и т.д.)

T1LO, Т1НО - нижнее и верхнее значения температуры окружающей среды, определяемые измерительным преобразователем, °С.

T2LO, Т2НО - нижнее и верхнее значения температуры окружающей среды, определяемые средством измерений, °С.

UL, UH - нижнее и верхнее значения напряжения питания, В.

TLG, THG - нижнее и верхнее значения температуры среды, °С.

PLG, PHG - нижнее и верхнее значения давления среды, бар.

DPL, DPH - нижнее и верхнее значения перепада давления на сужающем устройстве, бар.

ROSL, ROSH - нижнее и верхнее значения плотности природного газа при стандартных условиях, кг/м.

YIM[J,1] - верхние значения концентраций компонентов природного газа (1 J 16, см. 3.1.5).

YIM[J,2] - нижние значения концентраций компонентов природного газа (1 J 16, см. 3.1.5).

RoL, RoH - нижнее и верхнее значения плотности природного газа при рабочих условиях, кг/м.

4.1.14 Характеристика измеряемых параметров:

IfConstRo = 0, если плотность природного газа при стандартных условиях - непрерывно измеряемая величина; IfConstRo = 1, если плотность природного газа при стандартных условиях - условно-постоянная величина;

IfConstY = 0, если концентрации компонентов природного газа - непрерывно измеряемые величины; IfConstY = 1, если концентрации компонентов природного газа - условно-постоянные величины;

IfConstP = 0, если давление среды - непрерывно измеряемая величина; IfConstP = 1, если давление среды - условно-постоянная величина;

IfConstT = 0, если температура среды - непрерывно измеряемая величина; IfConstT = 1, если температура среды - условно-постоянная величина;

IfConstDp = 0, если перепад давления - непрерывно измеряемая величина; IfConstDp = 1, если перепад давления - условно-постоянная величина.

4.1.15 Характеристика измерительного комплекса (измерение плотности природного газа при рабочих условиях)

ERowA[J] - массив погрешностей измерительного комплекса (таблицы 15, 16).

Таблица 15

Наименование погрешности измерительного преобразователя и ее номер (J) в массиве ERowA[J]

Таблица 16

Наименование погрешности корнеизвлекающего средства измерений и ее номер (J) в массиве ERowA[J]

ERowA[21] - верхний предел измерения плотности газа при рабочих условиях, кг/м.

ERowA(22] - нижний предел измерения плотности газа при рабочих условиях, кг/м.

4.2 Выходные данные

RO, KAPPA, MU, QC, VC, VM, HS[1], HS[2] - см. 3.2.

HSV[1] - высшая теплота сгорания (энергосодержание) среды, МДж.

HSV[2] - низшая теплота сгорания (энергосодержание) среды, МДж.

EQR - случайная составляющая относительной погрешности расчета расхода среды для измерительного трубопровода с порядковым номером NNit.

EQS - систематическая составляющая относительной погрешности расчета расхода среды для измерительного трубопровода с порядковым номером NNit.

EQ1 - полная относительная погрешность расчета расхода среды для измерительного трубопровода с порядковым номером NNit.

ЕVС - полная относительная погрешность расчета количества среды для измерительного трубопровода с порядковым номером NNit за время TauAV.

ЕН[1] - полная относительная погрешность расчета высшей теплоты сгорания среды для измерительного трубопровода с порядковым номером NNit.

ЕН[2] - полная относительная погрешность расчета низшей теплоты сгорания среды для измерительного трубопровода с порядковым номером Nnit.

4.3 Листинг модуля расчета погрешности определения расхода и количества жидкостей и газов

4.3.1 Типы используемых переменных:

FL: text; NNIT, IERR, SIZMDPD, SIZMPD, SIZMTD: byte; T1HO, T1LO, UH, UL, PHG, PLG, THG, TLG, DPH, DPL, T2HO, T2LO, ROSH,ROSL, EDPS, EDPR, EDP, EPS, EPR, EP, ETS, ETR, ET, EROSS, EROSR, EROS, EROWS, EROWR, EROW, EQS, EQR, EQ1, EVC,DD20, DT20, DD, DT, ALFAT, ALFASU, RSH, RN, TAUP, TAUAV, ERRTAU, DHD, RO, KAPPA, MU, VM, ROSC, QC, VCV, VC, RoH, RoL, RoM: real; RocStr: string[6]; LRA: array [1..10] of real; DAC, TAC: array [1..3] of real; EPA: err1; YIM: err2; HS,HSV,EH: err3; NMETHKA, NSUBA, NSUZA, NVARYA, VARPA, VarRoA, AMOUNTR, SODSUA: err7; NRA: err8; EDPA: err9; ETA: err10; EROSA, ERowA: err11;

type err1 = array [1..26] of real; err2 = array [1..16,1..2] of real; err3 = array [1..2] of real; err7 = array [1..30] of byte; err8 = array[1..10] of byte; err9 = array [1..35] of real; err10 = array [0..35] of real; err11 = array [1..22] of real;

4.3.2 Листинг модуля ERRQSNX

Unit ERRQSNX;

Interface

{-------------------------------------------------}

Procedure QS;

Implementation

uses DOS, CRT;

{--------------------------------------------------------------------}

var

Zc: real;

const

Bi:array[1..16] of real = (0.0436,0.0894,0.1288,0.1783,0.1703,0.2345,

Procedure DP(E:err9; DPC:real; var EDPSY,EDPRA,EDP:real);

forward;

Procedure P(E:err1; PM:real; var EPSY,EPRA,EP:real);

forward;

Procedure T(E:err10; TA:real; var ETSY,ETRA,ET:real);

forward;

Procedure ROS(E:err11; RIM:err1; var EROSS,EROSR,EROS:real);

forward;

Procedure C(BET,RE,DP,KSH,KK:real; var ERC:real; var IERR:byte);

forward;

Procedure EPSI(BET,DP,P:real; var EEPS:real);

forward;

Procedure EPH(YR:err1; T,P,ROS:real; var ERO,EKAP,EVIS:real);

forward;

Procedure EPHP(N:byte; T,P:real; var ERO,EKAP,EVIS:real);

forward;

Procedure ROW(E:err11; var EROWS,EROWR,EROW:real);

forward;

Procedure QCalc(NRQ:byte; T,P,DP,YA,YY,DD0,DT0:real; YR:err1;

var BET,RE,KSH,KK,ROC,QC,VCV:real; var VCC:real; var VMC,ROX,

KAPPAX,MUX:real; var HS:err3);

forward;

Procedure HSP(RIM:err1; YA,YY:real; var ERRH1,ERRH2:real);

forward;

function r_(A,R:real):real; forward;

Procedure QS;

const

label

var

begin { QS }

PMG:= sqr(0.5*(sqrt(PHG)+sqrt(PLG)));

DVPG:= 0.0; if IfConstP <> 0 then DVPG:= 100.0*(PHG-PLG)/(PHG+PLG);

T1:=THG+273.15; T2:=TLG+273.15;

TMG:= 4.0*T1*T2/sqr(sqrt(T1)+sqrt(T2))-273.15;

DVTG:= 0.0; if IfConstT <> 0 then DVTG:= 100.0*(T1-T2)/(T1+T2);

DPM:= sqr(0.5*(sqrt(DPH)+sqrt(DPL)));

DVDP:= 0.0; if IfConstDp <> 0 then DVDP:= 100.0*(DPH-DPL)/(DPH+DPL);

if(NSubA[NNit]=0) and (VarRoA[NNit]=1) then

NRQ:= 1;

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

C(BET,RE,DPH,KSH,KK,ERC,IERR);

if IERR=0 then begin

DP(EDPA,DPM,EDPS,EDPR,EDP);

P(EPA,PMG,EPS,EPR,EP);

T(ETA,TMG,ETS,ETR,ET);

ROS(EROSA,RIM,EROSS,EROSR,EROS);

EPSI(BET,DPM,PMG,EEPS);

EPH(RIM,TMG,PMG,ROSC,ERO,EKAP,EVIS);

if (NSubA[NNit]=0) and (VarRoA[NNit]=1) then begin

end

NRQ:= 0;

TETDP:= 0.0;

if EDP <> 0.0 then begin

DPMD:= 0.001*EDP*DPM; DPM1:= DPM+DPMD; DPM2:=DPM-DPMD;

QCalc(NRQ,TMG,PMG,DPM1,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

QCalc(NRQ,TMG,PMG,DPM2,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

TETDP:= DPM*(QC1-QC2)/(2.0*DPMD)/QC;

end;

DDD:= 0.001*EDD*DD20; DD01:= DD20+DDD; DD02:= DD20-DDD;

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD01,DT20,RIM,BET,RE,KSH,KK,ROSC,

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD02,DT20,RIM,BET,RE,KSH,KK,ROSC,

TETDD:= DD20*(QC1-QC2)/(2.0*DDD)/QC;

DDT:= 0.001*EDT*DT20; DT01:= DT20+DDT; DT02:= DT20-DDT;

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT01,RIM,BET,RE,KSH,KK,ROSC,

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT02,RIM,BET,RE,KSH,KK,ROSC,

TETDT:= DT20*(QC1-QC2)/(2.0*DDT)/QC;

DKAP:= 0.001*EKAP*KAPPA; KAP1:= KAPPA+DKAP; KAP2:= KAPPA-DKAP;

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

TETKAP:= KAPPA*(QC1-QC2)/(2.0*DKAP)/QC;

DMU:= 0.001*EVIS*MU; MU1:= MU+DMU; MU2:= MU-DMU;

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

TETMU := MU*(QC1-QC2)/(2.0*DMU)/QC;

TETC:= 1.0; TETEPS:= 1.0; TETRO:= 0.5;

if(NSubA[NNit]=0) and (VarRoA[NNit]=l) then begin

TETRO:= 0.0;

end;

NRQ:= 1;

TETYA:= 0.0; TETYY:= 0.0; for I:= 1 to 16 do TETYI[I]:= 0.0;

if (NSUBA[NNIT] <> 0) or (NVARYA[NNIT]=0) then TETROS:= 1.0

else begin

TETROS:= 0.0;

if EROS <> 0.0 then begin

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROS1,

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROS2,

TETROS:= ROSC*(QC1-QC2)/(2.0*DROS)/QC;

end;

if (YA <> 0.0) and (EROSA[17] <> 0.0)then begin

QCalc(NRQ,TMG,PMG,DPM,YA1,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

QCalc(NRQ,TMG,PMG,DPM,YA2,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

end;

if(YY <> 0.0) and (EROSA[18] <> 0.0) then begin

QCalc(NRQ,TMG,PMG,DPM,YA,YY1,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

QCalc(NRQ,TMG,PMG,DPM,YA,YY2,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

end; {endif}

TETT:= 0.0;

if ET <> 0.0 then begin

QCalc(NRQ,TMG+DTMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,

QCalc(NRQ,TMG-DTMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,

end;

TETP:= 0.0;

if ЕР <> 0.0 then begin

QCalc(NRQ,TMG,PMG+DPMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,

QCalc(NRQ,TMG,PMG-DPMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,

if(NSUBA[NNIT]=0) and (NVARYA[NNIT]=0) then

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIMX,BET,RE,KSH,KK,ROS1,

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIMX,BET,RE,KSH,KK,ROS2,

EQR:= sqrt(sqr(TETDP)*(sqr(EDPR)+sqr(DVDP))+sqr(TETP)*

EQS:= sqrt(sqr(TETDP*EDPS)+sqr(TETP*EPS)+

EQ1:= sqrt(sqr(EQS)+sqr(EQR));

if ERRTAU=0.0 then ERRTAU:= 100.0*DHD/(60*TAUAV);

EVC:= sqrt(sqr(EQS)+sqr(EQR)+sqr(ERRTAU));

HSP(RIM,YA,YY,ERRH1,ERRH2);

if ERRH1=0.0 then begin

end else begin

case IERR of

'расстояние между первым перед сужающим устройством гидравлическим

сопротивлением');

'и сужающим устройством меньше допустимой величины

(7.2.4, 7.3.1 ГОСТ 8.563.1).');

'расстояние между сужающим устройством и гидравлическим сопротивлением за

сужающим устройством ');

'меньше допустимой величины (7.2.4 ГОСТ 8.563.1).');

'дополнительная погрешность коэффициента истечения, обусловленная

сокращением');

'длин прямых участков между сужающим устройством и гидравлическими сопро-

тивлениями,');

'превышает 1 % (7.2.4 ГОСТ 8.563.1).');

'эксцентриситет между осями измерительного трубопровода и сужающего

устройства превышает допустимую');

'величину (7.5.2.3 ГОСТ 8.563.1).');

'высота уступа между двумя секциями измерительного трубопровода превышает

допусти-');

'мую величину (7.5.1.4 ГОСТ 8.563.1).');

end; { QS }

var

if(E[1]<>0) or (E[2]<>0) then

else begin

end; {endif}

E2DPS:= E[8]*(T1HO-T1LO);

E3DPS:= E[9]*(UH-UL);

E6DPS:= E[12]*(PHG-PLG);

EDDPS:=sqrt(sqr(E1DPSC)+sqr(E2DPS)+sqr(E3DPS)+sqr(E[10])+

if SIZMDPD > 1 then begin

if(E[13]<>0) or (E[14]<>0) then

else begin

end; {endif}

EDP2S:=E[20]*(T2HO-T2LO);

EVDPS:=sqrt(sqr(EDP1SC)+sqr(EDP2S)+sqr(E[21])+sqr(E[22]));

if SIZMDPD=7 then begin

if (E[23]<>0) or (E[24]<>0) then

else begin

end; {endif}

EDPRS:= E[30]*(T2HO-T2LO);

ERDPS:= sqrt(sqr(EDPRSC)+sqr(EDPRS)+sqr(E[31])+sqr(E[32])):

end;

end; {endif}

YDP[1]:= E[35]/DPC; YDP[2]:=2.0*sqrt(E[35]/DPC);

0: begin

end;

1: begin

end;

2: begin

end;

3: begin

end;

4: begin

end;

5: begin

end;

6: begin

end;

7: begin

end;

end; { endcase }

EDP:=sqrt(sqr(EDPSY)+sqr(EDPRA));

var

begin { P }

if (E[1]<>0) or(E[2]<>0) then

else begin

if E[3]<>0 then

else begin

end; {endif}

E2PS:=E[8]*(T1HO-T1LO);

E3PS:=E[9]*(UH-UL);

EDPS:=sqrt(sqr(E1PSC)+sqr(E2PS)+sqr(E3PS)+sqr(E[10])+

if SIZMPD>1 then begin

if (E[12]<>0) or(E(13]<>0) then

else begin

end; {endif}

EP2S:=E[19]*(T2HO-T2LO);

EVPS:=sqrt(sqr(EP1SC)+sqr(EP2S)+sqr(E[20])+sqr(E[21]));

end;

if VARPA[NNIT]=0 then EBP:= sqr(E[24]*E[25]/PM);

0: begin

end;

1: begin

end;

2: begin

end;

3: begin

end;

4: begin

end;

5: begin

end;

6: begin

end;

end; { endcase }

EP := sqrt(sqr(EPSY)+sqr(EPRA));

var

begin { T }

end; {endif}

E2TS:= E[8]*(T1HO-T1LO);

E3TS:= E[9]*(UH-UL);

EDTS:= sqrt(sqr(E1TSC)+sqr(E2TS)+sqr(E3TS)+sqr(E[10])+

if (E[12]<>0) or (E[13]<>0) then

else begin

end; {endif}

ET2S:= E[19]*(T2HO-T2LO);

EVTS:=sqrt(sqr(ET1SC)+sqr(ET2S)+sqr(E[20])+sqr(E[21]));

if SIZMTD=6 then begin

if(E[22]<>0) or (E[23]<>0) then

else begin

end; {endif}

ETXS:=E(29]*(T2HO-T2LO);

EXTS:=sqrt(sqr(ET2SC)+sqr(ETXS)+sqr(E[30])+sqr(E[31]));

end; {endif}

end; {endif}

DET:=(E[34]-E[35])/(TA+273.15);

YT[1]:= DET; YT[2]:= 2.0*sqrt(DET);

0: begin

end;

1: begin

end;

2: begin

end;

3: begin

end;

4: begin

end;

5: begin

end;

6: begin

end;

end { endcase }

ET:= sqrt(sqr(ETSY)+sqr(ETRA));

const

var

if NSUBA[NNIT]<>0 then begin

end else begin

if(E(1]<>0) or (E[2]<>0) then

else begin

EDRSC:=sqrt(sqr(E1RSC)+sqr(E2RSC));

EROSS:=sqrt(sqr(EDRSC)+sqr(E[6])); EROSR:=E1RRC;

EROS:=sqrt(sqr(EROSS)+sqr(EROSR))

end; {endif}

var

const

label

var

const

label

begin { LinVent1 }

1: if L1=0.0 then exit; if LRA[1] < L1 then IERR:= 1

end; { LinVent1 }

var

const

label

begin { LinVent2 }

1: if LRA[I] < L1 then DL:= 0.5

end; { LinVent2 }

begin { С }

end;

3: NR:=AMOUNTR;

if (LRA[NR+1]/L1<0.5) or (LRA[NR+1]<0.5) then begin

end;

MR:= NRA[1];

L1:= AK[MR]+BK[MR]*r_(BET,CK[MR]);

MR:= NRA[2];

L2:= 0.5*(AK[MR]+BK[MR]*r_(0.7,CK[MR]));

if NSUZA[NNit]>2 then

L21:= LRA[2]-LRA[1];

if (LRA[1] < L1) and (L21 < L2) then begin

end;

EL:=0.0;

if NSUZA[NNit] <= 2 then

5: if DL>EL then EL:= DL

7: if DL>EL then EL:= DL

end;

L1:= AK[16]+BK[16]*r_(BET,CK[16]);

if LRA[NR+1]<L1 then EL:= EL+0.5;

ЕЕ:= 0.0;

if NSUZA[NNit]=0 then begin

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

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));

EEPS:= 0.0;

if NSUBA[NNIT]<>23 then begin

case NSUZA[NNIT] of

else

EEPS:= DP*(4.0+100.0*r_(BET, 8.0))/P;

end;

end; {endif}

var

const

if NSUBA[NNIT] = 0 then begin

end; { endcase }

end else begin

EPHP(N,T,P,ERO,EKAP,EVIS);

end; {endif}

const

if((9<=N) and (N<=11)) or((14<=N) and (N<=17)) or(N=21) or

end else begin

end; {endif}

if P<=PC[N,2] then EVIS:= EV[N,1]

else EVIS:= EV[N,2];

var

if(E[1]<>0)or(E[2]<>0)then

begin

end

else begin

end; {endif}

E2ROS:= E[8]*(T1HO-T1LO);

E3ROS:= E[9]*(UH-UL);

EDROS:= sqrt(sqr(E1ROSC)+sqr(E2ROS)+sqr(E3ROS)+sqr(E[10]));

if (E(11]<>0)or(E[12]<>0) then

else begin

end; {endif}

ERO2S:= Е[18]*(Т2НО-Т2LО);

EVROS:= sqrt(sqr(ERO1SC)+sqr(ERO2S)+sqr(E[19])+sqr(E[20]));

YRO[1]:= (E[21]-E[22])/RoM; YRO[2]:=(sqrt(E[21])-sqrt(E[22]))/sqrt(RoM);

EROWS:= sqrt(sqr(YRO[1]*EDROS)+sqr(YRO[2]*EVROS));

EROWR:= sqrt(sqr(YRO[1]*ElRORC)+sqr(YRO[2]*ERO1RC));

EROW:= sqrt(sqr(EROWS)+sqr(EROWR))

var

label

const

begin { QCalc }

if NRQ<>0 then begin

assign(Fl,'IRD'); rewrite(Fl);

if NSubA[NNIT]=0 then begin

end

end

writeln(Fl,NSubA[NNIT]);

writeln(Fl,P:14,Bl,T:14);

close(Fl);

TextColor(7);

gotoxy(19,9);

writeln('|-----------------------------------------------------|');

gotoxy(19,10);

writeln('| |');

gotoxy(19,ll);

write('| ');

TextColor(135);

write('Ж Д И Т Е');

TextColor(7);

writeln(' |');

gotoxy(19,12);

writeln('| |');

gotoxy(19,13);

writeln('| В Ы П О Л Н Я Е Т С Я Р А С Ч Е Т |');

gotoxy(19,14);

writeln('| |');

gotoxy(19,15);

writeln('|____________________________________________________|'),

if NSubA[NNIT]=0 then begin

end

end;

assign(Fl,'IRD'); reset(Fl);

if (NSubA[NNIT]=0) and (NMethKA[NNIT] >= 2) then

if NSubA[NNIT]=0 then begin

end

end;

end; {endif}

if (NSubA[NNit]=0) and (VarRoA[NNit]=1) and

Dd:= (1.0+AlfaSU*(Т-20.0))*Dd0;

Dt:=( 1.0+AlfaT*(Т-20.0))* Dt0;

Bet:= Dd/Dt; Bet4:= sqr(Bet)*sqr(Bet);

Ec:= 1.0/sqrt(l.0 - Bet4);

Eps:= 1.0;

if NSubA[NNIT]<>23 then begin

end;

KSh:= 1.0;

Kk:= 1.0;

if(NSuzA[NNIT]=0) and (Dd0<125.0) then begin

end;

if NSuzA[NNit]=0 then begin

Alfa:= Cb*Еs;

Qcb:= 0.039986*Alfa*KSh*Kk*Ерs*sqr(Dd)*

Re:= 4.e6*Qcb*Roc/3.6/3.141592653/MUX/Dt;

case NSuzA[NNIT] of

end;

Re:= Re*KRe;

if (NSuzA[NNit] <= 2) and (RSh <> 0.0) then begin

if (Re<2.е5) and (NSuzA[NNIT]>2) then

KCb:= 1.0;

if (Re<2.е5) and (NSuzA[NNIT)>2) then

Qc:= Qcb*KRe*KCb; Vcv:= Qc*TauAv; Vc:= Vcv;

Vm:= Vcv*Roc/1000.0; Qc:= Qc*Roc/1000.0 end;

end; { QCalc }

const

DHS: array[1..31] of real =(0.11,0.11,0.28,0.42,0.42,0.0,0.0,0.0,

HsNg1: array[l..16] of real= (37.04,64.91,92.29,119.7,119.3,

HsNg2 array[1..16] of real= (33.37,59.39,84.94,110.5,110.1,

var

if NSUBA[NNIT]=0 then begin

if(NMETHKA[NNIT]=0) or (NMETHKA[NNIT]=1) then begin

end else begin

end else begin

end; {endif}

function r_; begin r_:= exp(R*ln(A)) end;

Текст документа сверен по:

М.: ИПК Издательство стандартов, 1998