C  Программа расчёта переходного процесса в деривации и уравнительном

C  резервуаре с дополнительным сопротивлением

DIMENSION UT (20), UQST (20)

OPEN (1, FILE='dan')

OPEN (2, FILE='result')

READ (1,*) ALD,FD,FR

READ (1,*) AN,DZR

READ (1,*) DT,TRAS

READ (1,*) IK

READ (1,*) (UT(I),I=1,IK)

READ (1,*) (UQST(I),I=1,IK)

C Вспомогательные расчёты

DD=SQRT(FD*4./3.14)

R=DD/4.

IF (AN.EQ.0.) GO TO 15

C=1./AN*R**0.1667

AKWD=ALD/(C**2*R*FD**2)

GOTO 20

15 AKWD=0.

20 AKWR=DZR/(19.6*FD**2)

C Присваивание начальных условий

T=0.

CALL INT11 (T,QST,A,IK,UT,UQST)

QD=QST

QR=0

HWD=AKWD*QD*ABS(QD)

HWR=QD**2/(19.6*FD**2)

Z=-HWD-HWR

HD=Z+HWR

WRITE (2,1000)

WRITE (2,1001)

WRITE (2,1002)T,QST,QD,QR,Z,HWD,HWR,HD

C Решение дифференциальных уравнений

100 T=T+DT

IF (T.GT.TRAS) STOP

CALL INT11 (T,QST,A,IK,UT,UQST)

DZDT=(QD-QST)/FR

Z=Z+DZDT*DT

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

QD=QD+DQDDT*DT

HWD=AKWD*QD*ABS(QD)

QR=QD-QST

HWR=AKWR*QR*ABS(QR)+QD**2/(19.6*FD**2)

HD=Z+HWR

C Печать строки результатов

WRITE (2,1002) T,QST,QD,QR,Z,HWD,HWR,HD

GOTO 100

 1000 FORMAT(' ВРЕМЯ РАСХ. СТ. РАСХ. ДЕР. РАСХ. УР. УРОВЕНЬ ПОТ.ДЕР. '

1' ПОТ.УР ПЬЕЗ.ДЕР')

 1001 FORMAT(' СЕК М3/С М3/С М3/С М   М '

1' М   М ')

 1002 FORMAT(4F9.0,4F9.2)

CLOSE (1)

CLOSE (2)

END

ТЕКСТ программы линейной интерполяции

 

SUBROUTINE INT11(D,V,DYDX,N,X,Y)

C

C * * * * * * * * * * * * * * * * * * * * * * * * * * * * * S

C                                                      *

C    PROGRAM FOR LINEAR INTERPOLATION FUNCTIONS   *

C                                                      *

C    Y=F(X) AND COMPUTATION OF DERIVATIONS DY/DX  *

C                                                      *

C * * * * * * * * * * * * * * * * * * * * * * * * * * * * * F

C

DIMENSION X(N),Y(N)

DO 1 I=2,N

IF(D-X(I))2,2,1

1 CONTINUE

I=N

2 I1=I-1

DX=X(I)-X(I1)

DYDX=(Y(I)-Y(I1))/DX

V=(Y(I)*(D-X(I1))-Y(I1)*(D-X(I)))/DX

RETURN

END

 

ПРИМЕЧАНИЕ Трансляция выполняется с использованием стандартных процедур программного обеспечения Watcom FORTRAN

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

 

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

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

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

Строка 3 вводит 2 числа: шаг расчета по времени и время расчета

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

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

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

 

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

2100.,20.,190.

0.012,16.

1.,200.

3

0.,12.,200.

80.,0.,0.

 

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

 

Результаты расчета выводятся в файл RESULT. Содержимое этого файла может быть просмотрено любым текстовым редактором. Результаты выводятся в табличном виде. Параметры колебательного процесса расположены по столбцам:

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

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

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

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

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

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

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

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

 

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

 

ВРЕМЯ РАСХ.СТ. РАСХ.ДЕР. РАСХ.УР. УРОВЕНЬ ПОТ.ДЕР. ' ПОТ.УР ПЬЕЗ.ДЕР

СЕК М3/С М3/С М3/С М   М ' М  М

  0. 80. 80.  0. -4.36 3.55.82 -3.55

  1. 73. 80.  7. -4.33 3.55.91 -3.42

  2. 67. 80. 13. -4.26 3.55 1.18 -3.08

  3. 60. 80. 20. -4.15 3.54 1.63 -2.53

  4. 53. 80. 26. -4.01 3.53 2.24 -1.77

  5. 47. 80. 33. -3.84 3.52 3.03 -.81

  6. 40. 79. 39. -3.63 3.49 3.97.33

  7. 33. 79. 46. -3.39 3.46 5.05 1.66

  8. 27. 78. 52. -3.11 3.41 6.26 3.15

  9. 20. 78. 58. -2.81 3.36 7.60 4.79

10. 13. 77. 64.    -2.47 3.29 9.04 6.57

11.  7. 76. 69. -2.10 3.21 10.57 8.48

12.  0. 75. 75. -1.70 3.12 12.18 10.49

13.  0. 74. 74. -1.30 3.01 11.76 10.46

14.  0. 72. 72. -.91 2.90 11.35 10.44

 

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

 

 

ПРИЛОЖЕНИЕ 3
Текст программы расчета колебаний в деривации и уравнительном резервуаре камерного типа на языке Фортран

C пpогpамма pасчета пеpеходного пpоцесса в деpивации и камеpном

C pезеpвуаpе с дополнительным сопpотивлением

DIMENSION UT (20), UQC (20), UW (15), UZ (15)

OPEN(1,FILE='DAN')

OPEN(2,FILE='RESULT')

C ВВОД ИСХОДНЫХ ДАННЫХ

READ(1,*)ALD,FD,FR

READ(1,*)AN,DZR

READ(1,*)ZWOD,AMW,BWOD

READ(1,*)ZBHK,ZHHK,FHK

READ(1,*)DT

READ(1,*)IK

READ(1,*)(UT(I),I=1,IK)

READ(1,*)(UQC(I),I=1,IK)

C ОПРЕДЕЛЕНИЕ НАЧАЛЬНЫХ УСЛОВИЙ

      DD=SQRT(FD*4./3.14)

R=DD/4.

IF (AN.EQ.0.) GO TO 15

C=1./AN*R**0.1667

AKWD=ALD/(C**2*R*FD**2)

GOTO 20

15 AKWD=0.

20 AKWR=DZR/(19.6*FD**2)

T=0.

WBK=0.

CALL INT11(T,QC,A,IK,UT,UQC)

QD=QC

QR=0.

HWR=0.

HWD=AKWD*QD*ABS(QD)+QD**2/(19.62*FD**2)

Z =- HWD

HD = Z + HWR