Программа для определения диаметра простого трубопровода
Уравнение Бернулли для двух выбранных сечений, считая a1 = a2, можно написать ,
где – располагаемый напор (HS)*, м;
u – средняя скорость в трубопроводе (V), м/с; l – коэффициент гидравлического трения (LAMDA); Sl – общая длина (SUML), м; d – определяемый диаметр (D), м; S z – сумма коэффициентов местных сопротивлений (SUMDZE); – коэффициент; D2 – диаметр второго сечения (D2), м; – коэффициент; D1 – диаметр первого сечения (D1), м.
Необходимые для выполнения расчета другие величины обозначены: Q – расход (Q), м3/с; Dн – начальный диаметр трубопровода (DР), м; n – кинематическая вязкость (NI), м2/с;
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
Дата добавления: 2015-06-27; просмотров: 607;