Программа для определения диаметра простого трубопровода

Уравнение Бернулли для двух выбранных сечений, считая a₁ = a₂, можно написать ,

где – располагаемый напор (HS)*, м;

u – средняя скорость в трубопроводе (V), м/с; l – коэффициент гидравлического трения (LAMDA); Sl – общая длина (SUML), м; d – определяемый диаметр (D), м; S z – сумма коэффициентов местных сопротивлений (SUMDZE); – коэффициент; D₂ – диаметр второго сечения (D2), м; – коэффициент; D₁ – диаметр первого сечения (D1), м.

Необходимые для выполнения расчета другие величины обозначены: Q – расход (Q), м³/с; D_н – начальный диаметр трубопровода (DР), м; n – кинематическая вязкость (NI), м²/с;
D_э – эквивалентная шероховатость (DELTAE), м; – относительная шероховатость (K); Rе – число Рейнольдса (RЕ). На рис. 33 приводится полный текст программы, написанный на языке FORTRAN, для ЭВМ СМ-4.

___________

* Обозначение величины в программе

0001 REAL K, LAMBDA, NI

0002 TYPE 2

0003 2 FORMAT (‘ВВЕДИТЕ Q, HS, DP, D1, D2, SUML, SUMDZE, NI, DELTAE’)

0004 ACCEPT 3, Q, HS, DP, D1, D2, SUML, SUMDZE, NI, DELTAE

0005 3 FORMAT (9F10.7)

0006 PRINT 4

0007 4 FORMAT (10X, ‘ОПРЕДЕЛЕНИЕ ДИАМЕТРА’//)

0008 N=0

0009 D=DP

0010 7 K=DELTAE/D

0011 V=4*Q/(3.24*D*D)

0012 RE=V*D/NI

0013 IF (RE.GE.4000.) GOTO 5

0014 IF (RE.GT.2320.) GOTO 6

0015 LAMBDA=64./RE

0016 GOTO 8

0017 6 D=D+DP

0018 N=N+1

0019 GOTO 7

0020 5 LAMBDA=0.11*(((68./RE)+K)**0.25)

0021 8 IF (DL.EQ.0.) GOTO 15

0022 B=(D/DL)**4

0023 GOTO 16

0024 15 B=1

0025 16 IF (D2.EQ.0) GOTO 17

0026 A=(D/D2)

0027 GOTO 18

0028 17 A=1

0029 18 H=(V*V/19.62)*(LAMBDA*SUML/D+SUMDZE+A-B)

0030 IF (H.GT.HS) GOTO 9

0031 PRINT 10, Q, HS, D1, D2, D, V, RE, LAMBDA, H, N

0032 10 FORMAT (10X, ‘ДАНО:Q’, = F8.6, 12X, ‘HS=’, F6, 1, 3X, ‘, ‘D1=’, F5.3, 10X,

‘D2=’, F5.3, 10X, ‘ОПРЕДЕЛЕНО:D=’, 15.3, ‘V=’, F6.3, 10X, ‘RE=’, F9.1, 13X,

‘LAMBDA=’, F7.5, 2X, ‘H=’, F6.1, ‘N=’, 13)

0033 STOP

0034 9 D=D+DP

0035 N=N+1

0036 GOTO 7

0037 END

Рис. 33