C расчет наполнения верхней камеры

300 Z = ZWOD +(QR /(AMW * BWOD *4.43))**0.667

DQDDT=-(Z+HWD+HWR)*9.81*FD/ALD

QD=QD+DQDDT*DT

HD=Z+HWR

WBK=WBK+QR*DT

QR=QD-QC

WRITE(2,1200)T,QC,QD,QR,Z,HWD,HWR,HD,WBK

IF(QR.LE.0.)STOP

GOTO 100

CLOSE(1)

CLOSE(2)

 1000 FORMAT(' ВРЕМЯ РАСХ. С T РАСХ. ДЕР. РАСХ. РЕЗ УРОВ. РЕЗ '

*'ПОТЕРИ ДЕР ПОТЕРИ УР ПЬЕЗ.НАПОР ОБЪЕМ ВК ')

1100 FORMAT(' C   M3/C  M3/C M3/C  M '

*' M    M    M     M3 ')

 1200 FORMAT(9F10.1)

END

ПРИМЕЧАНИЯ:

1. Текст программы линейной интерполяции приведен в Приложении 1.

2. Трансляция выполняется с использованием стандартных процедур программного обеспечения Watcom FORTRAN.

                                                                                                           

                                                                                                

ПРИЛОЖЕНИЕ 4
ИНСТРУКЦИЯ ПО ВВОДУ ИСХОДНЫХ ДАННЫХ в программу расчета колебаний в деривации и уравнительном резервуаре камерного типа

При первом запуске исполняемого модуля RESER.EXE образуются два файла DAN и RESULT. Исходные данные вводятся в файл DAN текстовым редактором (например Notepad). Числа вводятся в файл построчно в десятичном формате (с десятичной точкой). Число от числа отделяется запятой.

Строка 1 вводит 3 числа: длину деривационного водовода, площадь деривационного водовода и площадь шахты уравнительного резервуара.

Строка 2 вводит 2 числа: коэффициент шероховатости облицовки и коэффициент дополнительного сопротивления уравнительного резервуара.

Строка 3 вводит 3 числа: отметку водослива верхней камеры, коэффициент расхода водослива, длину водослива.

Строка 4 вводит 3 числа: отметку верха нижней камеры, отметку низа нижней камеры, площадь горизонтального сечения нижней камеры

Строка 5 вводит 1 число: шаг расчета по времени.

Строка 6вводит 1 число – количество точек закона Q _тур (t) изменения расхода турбинного водовода; это целое число, оно вводится без десятичной точки.

Строка 7 вводит массив значений T(i) моментов времени закона Q _тур (t)

Строка 8 вводит массив значений расходов Q _тур (i) закона Q _тур (t)

ПРИМЕЧАНИЕ Все отметки задаются относительно уровня ВБ для верхового резервуара и относительно уровня НБ для низового резервуара.

 

ПРИМЕР ВВОДА ИСХОДНЫХ ДАННЫХ

2100.,19.6,40.

0.012,1.1

10.,0.41,15.7

-12.,-17.,150.

2.

3

0.,10.,50.

80.,0.,0.

 

ИНСТРУКЦИЯ ПО ВЫВОДИМЫМ РЕЗУЛЬТАТАМ РАСЧЕТА

 

Выполнение расчетов производится запуском программы RESKAM.exe. Результаты расчета выводятся в файл RESULT. Файл находится в директории, куда помещен исполняемый модуль RESKAM.exe.

Для просмотра могут использоваться любые текстовые редакторы. Результаты выводятся в табличном виде.

Параметры колебательного процесса расположены по столбцам:

1-ый столбец – время, с

2-ой столбец – расход турбинного водовода (водоводов)

3-ий столбец – расход деривационного водовода

4-ый столбец – расход в основании уравнительного резервуара

5-ый столбец – уровень в резервуаре

6-ой столбец – потери напора в деривационном водоводе

7-ой столбец – потери напора в уравнительном резервуаре

8-ой столбец – давление в деривации

9- ый столбец – объем воды, поступивший в верхнюю камеру

 

ПРИМЕР ТАБЛИЦЫ С РЕЗУЛЬТАТАМИ РАСЧЕТА

ВРЕМЯ РАСХ.СT РАСХ.ДЕР. РАСХ.РЕЗ УРОВ.РЕЗ ПОТЕРИ ДЕР ПОТЕРИ УР ПЬЕЗ.НАПОР ОБЪЕМ ВК

C   M3/C  M3/C M3/C  M     M    M    M     M3

  .0 80.0 80.0  .0 -2.4  2.4  .0 -2.4  .0

  2.0 64.0 79.9 15.9 -1.6  2.4  .0 -1.6  .0

  4.0 48.0 79.4 31.4  -.1  2.4  .0  -.1  .0

  6.0 32.0 78.6 46.6  2.3     2.4  .0  2.3  .0

  8.0 16.0 77.1 61.1  5.4  2.4  .0  5.4  .0

10.0  .0 75.0 75.0  9.3  2.3  .0  9.3  .0

12.0  .0 72.2 72.2 13.0  2.2  .0 13.0  .0

14.0  .0 69.7 69.7 11.9  2.0  .0 11.9 144.4

16.0  .0 67.2 67.2 11.8  1.9  .0 11.8 283.8

18.0  .0 64.7 64.7 11.8  1.7  .0 11.8 418.2

20.0  .0 62.3 62.3 11.7  1.6  .0 11.7 547.6

22.0  .0 59.9 59.9 11.7  1.5  .0 11.7 672.1

 

 

ГРАФИКИ ПЕРЕХОДНОГО ПРОЦЕССА

 

 

ПРИЛОЖЕНИЕ 5
Текст программы расчета колебаний в деривации и уравнительном резервуаре ПНЕВМАТИЧЕСКОГО типа на языке Фортран