12.3 Характеристика програми
Описаний алгоритм
реалізований нами в оригінальній
програмі GMVR, яка дає
можливість виконати експлуатаційний
гідравлічний розрахунок розгалуженої
газової мережі середнього чи високого
тиску довільної конфігурації. Методика
дозволяє розрахувати всі можливі
варіанти підключення відгалужень та
відводів. Для цього у програмі введені
масиви коефіцієнтів роботи відгалужень
і відводів
,
які приймаються рівними нулю при
непрацюючому відводі і рівним одиниці,
якщо споживач в кінці відводу відбирає
газ.
Текст програми GMVR наведений у додатку Ш. Перелік вхідних ідентифікаторів програми GMVR наведений у таблиці 12.1.
Таблиця 12.1 - Перелік вхідних ідентифікаторів програми GMVR
Назва параметра
|
Позначення |
Одиниця вимірювання |
|
в алгоритмі |
в програмі |
||
Абсолютний тиск газу на початку газової мережі |
|
|
МПа |
Абсолютний тиск газу на вході споживачів |
|
|
МПа |
Середня температура газу в газовій мережі |
|
T |
К |
Максимальна швидкість руху газу на ділянках мережі |
|
WM |
м/c |
Коефіцієнт динамічної в‘язкості газу за умов перекачування |
|
|
Па с |
Відносна густина газу за повітрям |
|
|
кг/м3 |
Абсолютна еквівалентна шорсткість |
|
KE |
см |
Крок зміни витрати газу |
|
|
м3/год |
Кількість ділянок основної магістралі |
|
|
- |
Масив довжин ділянок основної магістралі |
|
|
м |
Масив внутрішніх діаметрів основної магістралі |
|
|
М
|
Масив кількості ділянок на відгалуженнях |
ki |
K(I) |
- |
Масив довжин ділянок основної магістралі відгалужень |
|
|
м |
Масив внутрішніх діаметрів основної магістралі відгалужень |
|
|
м
|
Масив довжин ділянок відводів від відгалужень |
|
|
м |
Масив внутрішніх діаметрів відводів від відгалужень |
|
|
м |
12.4 Приклад визначення пропускної здатності розгалуженої газової мережі середнього тиску за програмою gmvr
Використовуючи
програму GMVR,
визначимо пропускну здатність розгалуженої
газової мережі середнього тиску,
розрахункова схема якої наведена на
рисунку 9.1. Абсолютний тиск газу на
початку мережі становить
МПа, абсолютний тиск газу на вході
споживачів
МПА. Абсолютна еквівалентна шорсткість
внутрішньої поверхні труб
м. Фізичні умови газу за умов газових
мереж: динамічна в’язкість газу
Па·
с, відносна густина газу
за повітрям
.
Середня температура газу у мережі
К, максимальна швидкість газу в
газопроводах середнього тиску
м/с.
Визначаємо пропускну здатність системи для таких чотирьох схем роботи відгалужень:
всі відгалуження відключені;
працює перше відгалуження;
працюють перше і друге відгалуження;
працюють всі три відгалуження.
Роздруківка із результатами визначення пропускної здатності розгалуженої газової мережі середнього тиску за програмою GMVR наведена у додатку Щ.
За результатами розрахунків складаємо технологічну карту режимів роботи розгалуженої газової мережі середнього тиску при різних схемах підключення відгалужень (див. таблицю 12.2).
Таблиця 12.2 – Результати визначення пропускної здатності
розгалуженої газової мережі середнього тиску
Схема роботи відгалужень |
Витрата газу, м3/год |
||||
на початковій ділянці магістралі |
на першому відгалуженні |
на другому відгалуженні |
на третьому відгалуженні |
на кінцевій ділянці магістралі |
|
Всі відгалуження відключені |
8160 |
- |
- |
- |
8160 |
Працює перше відгалуження |
11065 |
2935 |
- |
- |
8130 |
Працюють перше і друге відгалуження |
14215 |
2935 |
3560 |
|
7720 |
Працюють три відгалуження |
16250 |
2935 |
3560 |
6910 |
2845 |
Рисунок 12.1 – Схема газопровідної системи зі структурою “складного дерева”
Додаток А
Текст програми GMN1 для проектного розрахунку газової мережі низького тиску кільцевої структури
10 OPEN "o",#1,"gmn1"
20 PRINT "Кiльцева мережа низького тиску"
30 PRINT "Дiаметри дiльниць невiдомi"
40 KK=1
50 DIM D(25,25), L(25,25), Q(25,25), DP(25,25), NM(25), DPP(25,25)
60 DIM PQ(25,25), SP(25), SPA(25),SPQ(25), SS(25), DE(25)
70 DIM DQK(25),DQS(25),DQ(25), J(25,25),RE(25,25), DS(25)
80 DATA 3.2,3.59,3.9,4.1,5.1,7.0,8.3,10.2,12.5,15, 20.9,26.3, 31.5
90 FOR R=1 TO 13
100 READ DS(R):PRINT DS(R);
110 NEXT R:PRINT
120 PRINT "Ввести NU,RO, EPS, NK, KE"
130 INPUT NU,RO,EPS,NK,KE
140 PRINT NU;RO;EPS;NK;KE
150 WRITE #1,NU,RO,EPS,NK,KE
160 PRINT "Ввести T, PN, EPS1, LO, WM, QN, DPS"
170 INPUT T,PN,EPS1,LO,WM,QN,DPS
180 PRINT T;PN;EPS1;LO;WM;QN;DPS
190 WRITE #1,T,PN,EPS1,LO,WM,QN,DPS
200 FOR K=1 TO NK
210 PRINT " NM(K) "
220 INPUT NM(K)
230 PRINT "NM(K)="NM(K)
240 WRITE #1,NM(K)
250 FOR I=1 TO NM(K)
260 PRINT "Ввести Q(K,I),L(K,I),J(K,I)"
270 INPUT Q(K,I),L(K,I),J(K,I)
280 PRINT Q(K,I);L(K,I);J(K,I)
290 WRITE #1,Q(K,I),L(K,I),J(K,I)
300 NEXT I
310 NEXT K
320 PK=PN-DPS
330 PM=.5*(PN+PK)
340 IG=DPS/1.1/LO
350 W=WM
360 DN=.036238*SQR(QN*T/PM/.000001/W)
370 IR=75.9*(KE/DN+1922*NU*DN/QN)^.25*QN^2*RO /DN^5
380 IF ((IR-IG-EPS1)>0) THEN W=W-.01:GOTO 360
390 FOR K=1 TO NK
400 FOR I=1 TO NM(K)
410 D(K,I)=.036238*SQR(ABS(Q(K,I))*T/PM/.000001/W)
420 FOR R=1 TO 13
430 IF ((D(K,I)-DS(R))>0) GOTO 470
440 DZ=DS(R)
450 D(K,I)=DZ
460 GOTO 480
470 NEXT R
480 NEXT I
490 NEXT K
500 FOR K=1 TO NK
510 SP(K)=0
520 SPA(K)=0
530 SPQ(K)=0
540 FOR I=1 TO NM(K)
550 RE(K,I)=.0354*ABS(Q(K,I))/D(K,I)/NU
560 IF (RE(K,I)>2000) GOTO 590
570 DP(K,I)=1245000!*ABS(Q(K,I))*NU*RO*L(K,I) /D(K,I)^4
580 GOTO 640
590 IF (RE(K,I)>4000) GOTO 620
600 DP(K,I)=.567*ABS(Q(K,I))^2.333*RO*L(K,I)/ D(K,I)^5.333/NU^.333
610 GOTO 640
620 DPP(K,I)=75.9*(KE/D(K,I)+1922*NU*D(K,I)/ ABS(Q(K,I)))^.25
630 DP(K,I)=DPP(K,I)*Q(K,I)*Q(K,I)*RO*L(K,I)/D(K,I)^5
640 IF (Q(K,I)>0) GOTO 660
650 DP(K,I)=-DP(K,I)
660 SP(K)=SP(K)+DP(K,I)
670 SPA(K)=SPA(K)+ABS(DP(K,I))
680 NEXT I
690 NEXT K
700 FOR K=1 TO NK
710 DE(K)=200*SP(K)/SPA(K)
720 IF (ABS(DE(K))>EPS) GOTO 750
730 NEXT K
740 GOTO 1090
750 FOR K=1 TO NK
760 FOR I=1 TO NM(K)
770 PQ(K,I)=DP(K,I)/Q(K,I)
780 SPQ(K)=SPQ(K)+PQ(K,I)
790 NEXT I
800 SS(K)=0
810 DQK(K)=-(SP(K))/(1.75*SPQ(K))
820 NEXT K
830 DQK(0)=0
840 FOR K=1 TO NK
850 FOR I=1 TO NM(K)
860 FOR M=1 TO NK
870 IF(M-J(K,I)=0)GOTO 890
880 NEXT M
890 SS(K)=SS(K)+PQ(K,I)*DQK(M)
900 NEXT I
910 DQS(K)=SS(K)/SPQ(K)
920 DQ(K)=DQK(K)+DQS(K)
930 NEXT K
940 FOR K=1 TO NK
950 FOR I=1 TO NM(K)
960 Q(K,I)=Q(K,I)+DQ(K)
970 NEXT I
980 NEXT K
990 DQ(0)=0
1000 FOR K=1 TO NK
1010 FOR I=1 TO NM(K)
1020 FOR M=1 TO NK
1030 IF(M-J(K,I)=0) GOTO 1050
1040 NEXT M
1050 Q(K,I)=Q(K,I)-DQ(M)
1060 NEXT I
1070 NEXT K
1080 GOTO 500
1090 WRITE #1,W
1100 FOR K=1 TO NK
1110 WRITE #1,DE(K)
1120 FOR I=1 TO NM(K)
1130 WRITE #1,D(K,I),Q(K,I),L(K,I),DP(K,I)
1140 NEXT I
1150 NEXT K:CLOSE
1160 PRINT "Для друку ввести будь-яку цифру"
1170 INPUT O
1180 LPRINT "Кiльцева мережа низького тиску"
1190 LPRINT "Дiаметри дiльниць невiдомi"
1200 IF KK=1 THEN GOTO 1240
1210 DIM D(25,25), L(25,25), Q(25,25), DP(25,25), NM(25), DPP(25,25)
1220 DIM PQ(25,25), SP(25), SPA(25),SPQ(25), SS(25), DE(25)
1230 DIM DQK(25), DQS(25), DQ(25), J(25,25), RE(25,25), DS(25)
1240 OPEN "i",#1,"gmn1n-rez"
1250 DATA 3.2,3.59,3.9,4.1,5.1,7.0,8.3,10.2,12.5,15,20.9, 26.3, 31.5
1260 FOR R=1 TO 13
1270 READ DS(R)
1280 PRINT DS(R);
1290 NEXT R:LPRINT
1300 INPUT #1,NU,RO,EPS,NK,KE
1310 LPRINT NU;RO;EPS;NK;KE
1320 INPUT #1,T,PN,EPS1,LO,WM,QN,DPS
1330 LPRINT T;PN;EPS1;LO;WM;QN;DPS
1340 FOR K=1 TO NK
1350 INPUT #1,NM(K)
1360 LPRINT "NM(K)="NM(K)
1370 LPRINT " Q(K,I) L(K,I) J(K,I)"
1380 FOR I=1 TO NM(K)
1390 INPUT #1,Q(K,I),L(K,I),J(K,I)
1400 LPRINT USING " ####.# ###.# ## ";Q(K,I);L(K,I); J(K,I)
1410 NEXT I
1420 NEXT K
1430 LPRINT "Результати розрахунку"
1440 INPUT #1,W
1450 LPRINT "W="W
1460 FOR K=1 TO NK
1470 INPUT #1,DE(K)
1480 LPRINT "Контур="K, "DE="DE(K)
1490 LPRINT "D(K,I) Q(K,I) L(K,I) DP(K,I) "
1500 FOR I=1 TO NM(K)
1510 INPUT #1,D(K,I),Q(K,I),L(K,I),DP(K,I)
1520 LPRINT USING " ##.## ####.## ###. #### "; D(K,I), Q(K,I),L(K,I),DP(K,I)
1530 NEXT I
1540 NEXT K
1550 CLOSE
1560 END
Додаток Б
Результати проектного розрахунку газової мережі низького тиску кільцевої структури за програмою GMN1
Кiльцева мережа низького тиску
Дiаметри дiльниць невiдомi
3.2 3.59 3.9 4.1 5.1 7 8.3 10.2 12.5 15 20.9 26.3 31.5
Ввести NU,RO, EPS, NK, KE
1.366E-05 .756 1 5 .01
Ввести T, PN, EPS1, LO, WM, QN, DPS
280 104325 .2 1320 7 50 1200
NM(K)= 8
Q(K,I) L(K,I) J(K,I)
120.5 230.0 0
98.0 220.0 0
76.5 210.0 0
56.0 200.0 0
-306.0 200.0 0
-79.5 210.0 3
-34.5 250.0 2
-11.0 220.0 2
NM(K)= 6
Q(K,I) L(K,I) J(K,I)
34.5 250.0 1
11.0 220.0 1
35.0 220.0 0
-11.0 220.0 3
-60.0 240.0 5
-12.5 250.0 4
NM(K)= 4
Q(K,I) L(K,I) J(K,I)
79.5 210.0 1
11.0 220.0 2
-195.5 210.0 0
-83.0 220.0 5
NM(K)= 4
Q(K,I) L(K,I) J(K,I)
12.5 250.0 2
12.0 240.0 0
-11.5 230.0 5
-11.5 230.0 0
NM(K)= 6
Q(K,I) L(K,I) J(K,I)
83.0 220.0 3
60.0 240.0 2
11.5 230.0 4
-79.5 230.0 0
-57.0 220.0 0
-34.5 230.0 0
Результати розрахунку
W= 2.599933
Контур= 1 DE= 6.584805E-02
D(K,I) Q(K,I) L(K,I) DP(K,I)
15.00 125.43 230.0 68
12.50 102.93 220.0 110
12.50 81.43 210.0 69
10.20 60.93 200.0 104
20.90 -301.07 200.0 -57
12.50 -65.16 210.0 -46
7.00 -27.19 250.0 -190
3.90 -3.69 220.0 -58
Контур= 2 DE= .8032928
D(K,I) Q(K,I) L(K,I) DP(K,I)
7.00 27.19 250.0 190
3.90 3.69 220.0 58
7.00 32.61 220.0 232
3.90 -3.97 220.0 -69
10.20 -59.26 240.0 -119
5.10 -14.58 250.0 -287
Контур= 3 DE= .1665969
D(K,I) Q(K,I) L(K,I) DP(K,I)
12.50 65.16 210.0 46
3.90 3.97 220.0 69
20.90 -204.91 210.0 -30
12.50 -89.28 220.0 -85
Контур= 4 DE= 2.754324E-02
D(K,I) Q(K,I) L(K,I) DP(K,I)
5.10 14.58 250.0 287
4.10 11.69 240.0 534
4.10 -8.68 230.0 -300
4.10 -11.81 230.0 -521
Контур= 5 DE= .1655761
D(K,I) Q(K,I) L(K,I) DP(K,I)
12.50 89.28 220.0 85
10.20 59.26 240.0 119
4.10 8.68 230.0 300
12.50 -82.63 230.0 -78
10.20 -60.13 220.0 -112
7.00 -37.63 230.0 -313
Додаток В
Текст програми GMN2 для розрахунку газової мережі
низького тиску кільцевої структури при відомих діаметрах
10 REM "РОЗРАХУНОК КIЛЬЦЕВИХ ГАЗОВИХ МЕРЕЖ НИЗЬКОГО ТИСКУ"
20 DIM D(30,30),L(30,30),Q(30,30),DP(30,30),PQ(30,30)
30 DIM SP(30),SPA(30),SPQ(30),SS(30),DE(30)
40 DIM DQK(30),DQS(30),DQ(30),J(30,30),RE(30,30), NM(30)
50 OPEN "o",#1,"rez"
60 INPUT "ВВЕСТИ HU,RO,EPS,NK,KE";HU,RO,EPS, NK,KE
70 WRITE #1,HU,RO,EPS,NK,KE
80 FOR K=1 TO NK
90 PRINT "NM(";K;")=";
100 INPUT NM(K)
110 WRITE #1,NM(K)
120 NEXT K
130 FOR K=1 TO NK
140 FOR I=1 TO NM(K)
150 PRINT "Ввести вн.дiaметр(см),витрату, геом.довжину(м), номер сус.кiльця "
160 PRINT "D(";K;",";I;"),";"q(";K;",";I;"),";"l(";K;",";I;"), ";"j(";K;",";I;")"
170 INPUT D(K,I),Q(K,I),L(K,I),J(K,I)
180 WRITE #1,D(K,I),Q(K,I),L(K,I),J(K,I)
190 NEXT I
200 NEXT K
210 FOR K=1 TO NK
220 SP(K)=0
230 SPA(K)=0
240 SPQ(K)=0
250 FOR I=1 TO NM(K)
260 RE(K,I)=.0354*ABS(Q(K,I))/D(K,I)/HU
270 IF RE(K,I)<=2000 GOTO 320
280 IF RE(K,I)<=4000 GOTO 340
290 F1=75.9*(KE/D(K,I)+1922*HU*D(K,I)/ABS(Q(K,I)))^.25
300 DP(K,I)=F1*Q(K,I)*Q(K,I)*RO*L(K,I)/D(K,I)^5
310 GOTO 350
320 DP(K,I)=1245000!*ABS(Q(K,I))*HU*RO*L(K,I) /D(K,I)^4
330 GOTO 350
340 DP(K,I)=.567*ABS(Q(K,I))^2.333*RO*L(K,I)/ D(K,I)^5.333/HU^.333
350 IF Q(K,I)<0 THEN DP(K,I)=-DP(K,I)
360 SP(K)=SP(K)+DP(K,I)
370 SPA(K)=SPA(K)+ABS(DP(K,I))
380 NEXT I
390 NEXT K
400 FOR K=1 TO NK
410 DE(K)=200*SP(K)/SPA(K)
420 IF ABS(DE(K))>=EPS GOTO 450
430 NEXT K
440 GOTO 810
450 FOR K=1 TO NK
460 FOR I=1 TO NM(K)
470 PQ(K,I)=DP(K,I)/Q(K,I)
480 SPQ(K)=SPQ(K)+PQ(K,I)
490 NEXT I
500 SS(K)=0
510 DQK(K)=-(SP(K))/(1.75*SPQ(K))
520 NEXT K
530 DQK(0)=0
540 FOR K=1 TO NK
550 FOR I=1 TO NM(K)
560 FOR M=1 TO NK
570 IF J(K,I)<>M THEN GOTO 600
580 SS(K)=SS(K)+PQ(K,I)*DQK(M)
590 GOTO 610
600 NEXT M
610 NEXT I
620 DQS(K)=SS(K)/SPQ(K)
630 DQ(K)=DQK(K)+DQS(K)
640 NEXT K
650 FOR K=1 TO NK
660 FOR I=1 TO NM(K)
670 Q(K,I)=Q(K,I)+DQ(K)
680 NEXT I
690 NEXT K
700 DQ(0)=0
710 FOR K=1 TO NK
720 FOR I=1 TO NM(K)
730 FOR M=1 TO NK
740 IF J(K,I)<>M GOTO 770
750 Q(K,I)=Q(K,I)-DQ(M)
760 GOTO 780
770 NEXT M
780 NEXT I
790 NEXT K
800 GOTO 210
810 FOR K=1 TO NK
820 WRITE #1,DE(K)
830 FOR I=1 TO NM(K)
840 WRITE #1,D(K,I),Q(K,I),L(K,I),DP(K,I)
850 NEXT I
860 NEXT K
870 CLOSE
880 PRINT "Для друку результатiв -клавiша 1 "
890 INPUT AA
900 IF AA<>1 THEN GOTO 1310
910 LPRINT " "
920 LPRINT "ВХIДНI ДАНI"
930 OPEN "i",#1,"rez"
940 INPUT #1,HU,RO,EPS,NK,KE
950 LPRINT "HU=";HU,"RO=";RO,"EPS=";EPS,"NK=";NK, "KE=";KE
960 LPRINT " "
970 LPRINT "МАСИВ NM"
980 FOR K=1 TO NK
990 INPUT #1,NM(K)
1000 LPRINT "NM(";K;")=";NM(K)
1010 NEXT K
1020 LPRINT "______________________________________"
1030 LPRINT "I Дiлянка I Дiaметр, I Витрата, I Довжина, I Номер сус.I"
1040 LPRINT "I I см Iкуб.м/год I м I кiльця I"
1050 LPRINT "I---------I----------I----------I----------I-----------I"
1060 FOR K=1 TO NK
1070 LPRINT "I I Kiльцe";K;" I I I"
1080 FOR I=1 TO NM(K)
1090 INPUT #1,D(K,I),Q(K,I),L(K,I),J(K,I)
1100 LPRINT USING "I ## I ###.# I ####.# I ###.# I ## I";I,D
(K,I),Q(K,I),L(K,I),J(K,I)
1110 NEXT I
1120 LPRINT "--------------------------------------------------------"
1130 NEXT K
1140 LPRINT "--------------------------------------------------------"
1150 LPRINT " "
1160 LPRINT "Результати розрахунку кiлець мережi низького тиску"
1170 LPRINT "______________________________________"
1180 LPRINT "I Дiлянка I Дiaметр, I Витрата, I Довжина, I Перепад I"
1190 LPRINT "I I см Iкуб.м/год I м I тиску,Па I"
1200 LPRINT "I---------I----------I----------I----------I----------I"
1210 FOR K=1 TO NK
1220 LPRINT "I I Кiльцe";K;" I I"
1230 INPUT #1,DE(K)
1240 LPRINT USING "I Dk,%= #.###";ABS(DE(K))
1250 FOR I=1 TO NM(K)
1260 INPUT #1,D(K,I),Q(K,I),L(K,I),DP(K,I)
1270 LPRINT USING "I ## I ###.## I ####.# I ###.# I #### I";I,D(K,I),Q(K,I),L(K,I),DP(K,I)
1280 NEXT I
1290 NEXT K
1300 LPRINT "-------------------------------------------------------"
1310 END
Додаток Г
Результати розрахунку газової мережі низького тиску
кільцевої структури при відомих діаметрах газопроводів
за програмою GMN2
HU= 1.366E-05 RO= .756 EPS= 1 NK= 5 KE= .01
МАСИВ NM
NM( 1 )= 8
NM( 2 )= 6
NM( 3 )= 4
NM( 4 )= 4
NM( 5 )= 6
__________________________________________________
I Дiлянка I Дiaметр, I Витрата, I Довжина, I Номер сус.I
I I см I куб.м/год I м I кiльця I
I------------I-------------I--------------I---------------I--------------I
I I Kiльцe 1 I I I
I 1 I 15.0 I 120.5 I 230.0 I 0 I
I 2 I 12.5 I 98.0 I 220.0 I 0 I
I 3 I 12.5 I 76.5 I 210.0 I 0 I
I 4 I 10.2 I 56.0 I 200.0 I 0 I
I 5 I 20.9 I -306.0 I 200.0 I 0 I
I 6 I 12.5 I -79.5 I 210.0 I 3 I
I 7 I 7.0 I -34.5 I 250.0 I 2 I
I 8 I 3.9 I -11.0 I 220.0 I 2 I
I I Kiльцe 2 I I I
I 1 I 7.0 I 34.5 I 250.0 I 1 I
I 2 I 3.9 I 11.0 I 220.0 I 1 I
I 3 I 7.0 I 35.0 I 220.0 I 0 I
I 4 I 3.9 I -11.0 I 220.0 I 3 I
I 5 I 10.2 I -60.0 I 240.0 I 5 I
I 6 I 5.1 I -12.5 I 250.0 I 4 I
I I Kiльцe 3 I I I
I 1 I 12.5 I 79.5 I 210.0 I 1 I
I 2 I 3.9 I 11.0 I 220.0 I 2 I
I 3 I 20.9 I -195.5 I 210.0 I 0 I
I 4 I 12.5 I -83.0 I 220.0 I 5 I
I I Kiльцe 4 I I I
I 1 I 5.1 I 12.5 I 250.0 I 2 I
I 2 I 4.1 I 12.0 I 240.0 I 0 I
I 3 I 4.1 I -11.5 I 230.0 I 5 I
I 4 I 4.1 I -11.5 I 230.0 I 0 I
I I Kiльцe 5 I I I
I 1 I 12.5 I 83.0 I 220.0 I 3 I
I 2 I 10.2 I 60.0 I 240.0 I 2 I
I 3 I 4.1 I 11.5 I 230.0 I 4 I
I 4 I 12.5 I -79.5 I 230.0 I 0 I
I 5 I 10.2 I -57.0 I 220.0 I 0 I
I 6 I 7.0 I -34.5 I 230.0 I 0 I
Результати розрахунку кiлець мережi низького тиску
___________________________________________________
I Дiлянка I Дiaметр, I Витрата, I Довжина, I Перепад I
I I см I куб.м/год I м I тиску,ПаI
I-------------I---------------I---------------I--------------I-------------I
I I Кiльцe 1 I I
I Dk,%= 0.066
I 1 I 15.00 I 125.4 I 230.0 I 68 I
I 2 I 12.50 I 102.9 I 220.0 I 110 I
I 3 I 12.50 I 81.4 I 210.0 I 69 I
I 4 I 10.20 I 60.9 I 200.0 I 104 I
I 5 I 20.90 I -301.1 I 200.0 I -57 I
I 6 I 12.50 I -65.2 I 210.0 I -46 I
I 7 I 7.00 I -27.2 I 250.0 I -190 I
I 8 I 3.90 I -3.7 I 220.0 I -58 I
I I Кiльцe 2 I I
I Dk,%= 0.803
I 1 I 7.00 I 27.2 I 250.0 I 190 I
I 2 I 3.90 I 3.7 I 220.0 I 58 I
I 3 I 7.00 I 32.6 I 220.0 I 232 I
I 4 I 3.90 I -4.0 I 220.0 I -69 I
I 5 I 10.20 I -59.3 I 240.0 I -119 I
I 6 I 5.10 I -14.6 I 250.0 I -287 I
I I Кiльцe 3 I I
I Dk,%= 0.167
I 1 I 12.50 I 65.2 I 210.0 I 46 I
I 2 I 3.90 I 4.0 I 220.0 I 69 I
I 3 I 20.90 I -204.9 I 210.0 I -30 I
I 4 I 12.50 I -89.3 I 220.0 I -85 I
I I Кiльцe 4 I I
I Dk,%= 0.028
I 1 I 5.10 I 14.6 I 250.0 I 287 I
I 2 I 4.10 I 11.7 I 240.0 I 534 I
I 3 I 4.10 I -8.7 I 230.0 I -300 I
I 4 I 4.10 I -11.8 I 230.0 I -521 I
I I Кiльцe 5 I I
I Dk,%= 0.166
I 1 I 12.50 I 89.3 I 220.0 I 85 I
I 2 I 10.20 I 59.3 I 240.0 I 119 I
I 3 I 4.10 I 8.7 I 230.0 I 300 I
I 4 I 12.50 I -82.6 I 230.0 I -78 I
I 5 I 10.20 I -60.1 I 220.0 I -112 I
I 6 I 7.00 I -37.6 I 230.0 I -313 I
Додаток Д
Текст програми GMN3 для проектного розрахунку газової мережі низького тиску розгалуженої структури
10 PRINT "Розрахунок розгалужених мереж низького тиску"
20 PRINT "Дiаметри дiлянок невiдомi"
30 DIM DS(25),Q(20),L(20),D(20),DP(20)
40 DATA 3.2,3.59,3.9,4.1,5.1,7.0,8.3,10.2,12.5,15.0,20.9,26.3,31.5,41.4,51.8
50 FOR R=1 TO 15
60 READ DS(R):PRINT DS(R);
70 NEXT R
80 PRINT
90 PRINT "Початковў данў"
100 PRINT "Введить NU,RO,T,EPS,KE
110 INPUT NU,RO,T,EPS,KE
120 PRINT "NU="NU;"RO="RO;"T="T;"EPS="EPS;"KE="KE
130 PRINT "Введўть PN,PK,QN,WN,NI"
140 INPUT PN,PK,QN,WN,NI
150 PRINT
160 LO=0
170 FOR I=1 TO NI
180 PRINT "Введўть номер дiлянки,L(I),Q(I)"
190 INPUT J$(I),L(I),Q(I)
200 LO=LO+L(I)
210 NEXT I
220 DPS=PN-PK
230 PM=.5*(PN+PK)+101325!
240 IG=DPS/1.1/LO
250 W=WN
260 DN=.036238*SQR(QN*T/PM/W/.000001)
270 IR=75.9*(KE/DN+1922*NU*DN/QN)^.25*QN^2*RO/DN^5
280 IF((IR-IG-EPS)>0) THEN W=W-.01:GOTO 260
290 FOR I=1 TO NI
300 D(I)=.036238*SQR(Q(I)*T/PM/W/.000001)
310 FOR R=1 TO 25
320 IF((D(I)-DS(R))>0) GOTO 350
330 D(I)=DS(R)
340 GOTO 360
350 NEXT R
360 NEXT I
370 SP=0
380 FOR I=1 TO NI
390 RE(I)=.0354*Q(I)/D(I)/NU
400 IF(RE(I)>2000) GOTO 430
410 DP(I)=1245000!*Q(I)*NU*RO*L(I)/D(I)^4
420 GOTO 470
430 IF(RE(I)>4000) GOTO 460
440 DP(I)=.567*Q(I)^2.333*RO*L(I)/D(I)^5.333/NU^.333
450 GOTO 470
460 DP(I)=75.9*(KE/D(I)+1922*NU*D(I)/Q(I))^.25*Q(I)^2*RO*L(I)/D(I)^5
470 SP=SP+DP(I)
480 PK(I)=PN-SP
490 NEXT I
500 PKD=PN-SP
510 PRINT "PK="PKD
520 PRINT "Якщо результати не оптимальнi,введўть iнше QN,iнакше 0"
530 INPUT Q
540 IF(Q=0) GOTO 570
550 QN=Q
560 GOTO 250
570 PRINT "PN="PN;"PK="PK;"QN="QN;"WN="WN;"NI="NI
580 PRINT "Результати розрахунку"
590 PRINT "W="W,"PKдiйсне="PKD
600 PRINT "Дiлянка L(I) Q(I) D(I) PK(I)"
610 FOR I=1 TO NI
620 PRINT J$(I);
630 PRINT USING " #### ###.# ##.## #### ";L(I),Q(I),D(I),PK(I)
640 NEXT I
650 PRINT "Для продовження розрахункўв введўть будь-яку цифру,для закўнчення 0"
660 INPUT M
670 IF M=0 GOTO 700
680 LPRINT
690 GOTO 130
700 END
Додаток Ж
Результати проектного розрахунку газової мережі
низького тиску розгалуженої структури
за програмою GMN3
Розрахунок розгалужених мереж низького тиску
Дiаметри дiлянок невiдомi
3.2 3.59 3.9 4.1 5.1 7 8.3 10.2 12.5 15 20.9 26.3 31.5 41.4 51.8
Початкові дані
Введіть NU,RO,T,EPS,KE
? 13.66e-6,.756,280,.1,.01
NU= 1.366E-05 RO= .756 T= 280 EPS= .1 KE= .01
Введіть PN,PK,QN,WN,NI
? 3000,1800,50,7,5
Введіть номер дiлянки,L(I),Q(I)
? ГРП-1,250,148.5
? 1-2,220,125
? 2-3,210,59.5
? 3-4,310,33.5
? 4-5,180,9
PK= 1911.745
Якщо результати не оптимальнi,введўть iнше QN,iнакше 0
? 0
PN= 3000 PK= 1800 QN= 50 WN= 7 NI= 5
Результати розрахунку
W= 2.609933 PKдiйсне= 1911.745
Дiлянка L(I) Q(I) D(I) DP(I) PK(I)
ГРП-1 250 148.5 15.00 100 2900
1-2 220 125.0 15.00 65 2835
2-3 210 59.5 10.20 105 2730
3-4 310 33.5 7.00 343 2388
4-5 180 9.0 3.59 476 1912
Для продовження розрахунків введіть будь-яку цифру, для закінчення 0
? 1
Введіть PN,PK,QN,WN,NI
? 2835,1800,50,7,2
Введіть номер дiлянки,L(I),Q(I)
? 2-6,270,30.5
? 6-7,170,8.5
PK= 1873.685
Якщо результати не оптимальнi,введіть iнше QN,iнакше 0
? 0
PN= 2835 PK= 1800 QN= 50 WN= 7 NI= 2
Результати розрахунку
W= 3.589932 PKдiйсне= 1873.685
Дiлянка L(I) Q(I) D(I) DP(I) PK(I)
2-6 270 30.5 7.00 252 2583
6-7 170 8.5 3.20 709 1874
Для продовження розрахунків введіть будь-яку цифру,для закінчення 0
? 0
Додаток К
Текст програми GMN4 для розрахунку газової мережі
низького тиску розгалуженої структури при відомих діаметрах
10 OPEN "o",#1,"rez"
20 OPEN "o",#2,"rez1"
30 PRINT "Розрахунок розгалужених мереж низького тиску"
40 PRINT "Дiаметри дiльниць заданi"
50 DIM L(30),Q(30),D(30),DP(30),RE(30),PK(30),J$(30)
60 PRINT "Ввести NU,RO,KE,PN,N"
70 PRINT "Початковi данi"
80 INPUT NU,RO,KE,PN,N
90 PRINT "NU="NU;"RO="RO;"Ke="KE;"PN="PN;"N="N
100 WRITE #2,NU,RO,KE,PN,N
110 PRINT "РЕЗУЛЬТАТИ РОЗРАХУНКУ":PRINT
120 FOR I=1 TO N
130 PRINT "ВВЕДIТЬ НОМЕР ДIЛЬНИЦI ,L(I),Q(I),D(I)"
140 INPUT J$(I),L(I),Q(I),D(I)
150 NEXT I
160 SP=0
170 FOR I=1 TO N
180 RE(I)=.0354*Q(I)/D(I)/NU
190 IF(RE(I)>2000) GOTO 220
200 DP(I)=1245000!*Q(I)*NU*RO*L(I)/D(I)^4
200 DP(I)=1245000!*Q(I)*NU*RO*L(I)/D(I)^4
210 GOTO 260
220 IF (RE(I)>4000) GOTO 250
230 DP(I)=.568*Q(I)^2.333*RO*L(I)/D(I)^5.333/NU^.333
240 GOTO 260
250 DP(I)=75.9*(KE/D(I)+1922*NU*D(I)/Q(I))^.25*Q(I)^2*RO*L(I)/D(I)^5
260 SP=SP+DP(I)
270 PK(I)=PN-SP
280 PRINT "УЧАСТОК="J$(I),"L="L(I);"Q="Q(I);"D="D(I);"DP="DP(I);"PK="PK(I)"
290 WRITE #1,J$(I),L(I),Q(I),D(I),DP(I),PK(I)
300 NEXT I
310 PRINT "Для продовження ввести 1, для друку -0":INPUT Z
320 IF Z<>1 GOTO 350
330 INPUT "PN,N";PN,N
340 GOTO 100
350 CLOSE:OPEN "i",#1,"rez":OPEN "i",#2,"rez1"
360 LPRINT "Розрахунок розгалужених мереж низького тиску":LPRINT
370 LPRINT "Дiаметри дiльниць заданi":LPRINT
380 WHILE NOT(EOF(2))
390 LPRINT "Початковi данi":LPRINT
400 INPUT #2,NU,RO,KE,PN,N
410 LPRINT "NU="NU;"RO="RO;"Ke="KE;"PN="PN;"N="N:LPRINT
420 LPRINT "РЕЗУЛЬТАТИ РОЗРАХУНКУ":LPRINT
430 FOR I=1 TO N
440 INPUT #1,J$,L,Q,D,DP,PK
450 LPRINT "УЧАСТОК="J$,"L="L;"Q="Q;"D="D;"DP="DP;"PK="PK
460 NEXT I:LPRINT
470 WEND:CLOSE
480 END
Додаток Л
Результати розрахунку газової мережі
низького тиску розгалуженої структури при відомих діаметрах ділянок за програмою GMN4
Розрахунок розгалужених мереж низького тиску
Дiаметри дiлянок заданi
Ввести NU,RO,KE,PN,N
Початковi данi
13.66e-6,.756,.01,3000,5
NU= 1.366E-05 RO= .756 Ke= .01 PN= 3000 N= 5
ВВЕДIТЬ НОМЕР ДIЛЯНКИ ,L(I),Q(I),D(I)
ГРП-1,250,148.5,15
1-2,220,125,15
2-3,210,59.5,10.2
3-4,310,33.5,7
4-5,180,9,3.59
РЕЗУЛЬТАТИ РОЗРАХУНКУ
Дiлянка L(I) Q(I) D(I) DP(I) PK(I)
ГРП-1 250 148.5 15.00 100 2900
1-2 220 125.0 15.00 65 2835
2-3 210 59.5 10.20 105 2730
3-4 310 33.5 7.00 343 2388
4-5 180 9.0 3.59 476 1912
Для продовження ввести 1
Ввести PN,N
2835,2
ВВЕДIТЬ НОМЕР ДIЛЯНКИ ,L(I),Q(I),D(I)
2-6,270,30.5,7
6-7,170,8.5,3.2
РЕЗУЛЬТАТИ РОЗРАХУНКУ
Дiлянка L(I) Q(I) D(I) DP(I) PK(I)
2-6 270 30.5 7.00 252 2583
6-7 170 8.5 3.20 709 1874
Для продовження ввести 1
? 0
Ok_
Додаток М
Текст програми GMN5 для проектного розрахунку
внутрішніх газових мереж житлового будинку
10 PRINT "Розрахунок внутрiшнiх газових мереж житлового будинку без номограм
20 DIM DS(8),Q(30),L(30),D(30),DP(30),DPG(30),J$(30),
DPS(30),RE(30)
30 DIM A(30),B(30),LR(30),H(30)
40 DATA 1.57,2.12,2.71,3.2,3.59,3.9,4.1,5.1
50 FOR R=1 TO 8
60 READ DS(R):PRINT DS(R);
70 NEXT R
80 PRINT
90 PRINT "Початковi данi"
100 PRINT "Ввести NU,RO,T,EPS,KE"
110 INPUT NU,RO,T,EPS,KE
120 PRINT "NU="NU;"RO="RO;"T="T;"EPS="EPS;"KE="KE
130 PRINT "Ввести напрям руху газу"
140 INPUT A$
150 PRINT "Ввести PN,DPD,QN,WN,NI"
160 INPUT PN,DPD,QN,WN,NI
170 PRINT
180 LRO=0
190 PRINT "Ввести номер дiлянки ,L(I),A(I),Q(I), вертикальнiсть дiлянки"
200 FOR I=1 TO NI
210 INPUT J$(I),L(I),A(I),Q(I),B(I)
220 LR(I)=L(I)*(1+A(I)/100)
230 LRO=LRO+LR(I)
240 NEXT I
250 PK=PN-DPD
260 PM=.5*(PN+PK)+101325!
270 IG=DPD/LRO
280 W=WN
290 DN=.036238*SQR(QN*T/PM/W/.000001)
300 IR=75.9*(KE/DN+1922*NU*DN/QN)^.25*QN^2*RO/DN^5
310 IF((IR-IG-EPS)>0) THEN W=W-.01:GOTO 290
320 FOR I=1 TO NI
330 D(I)=.036238*SQR(Q(I)*T/PM/W/.000001)
340 FOR R=2 TO 8
350 IF ((D(I)-DS(R))>0) GOTO 380
360 D(I)=DS(R-1)
370 GOTO 390
380 NEXT R
390 NEXT I
400 DPS=0
410 FOR I=1 TO NI
420 RE(I)=.0354*Q(I)/D(I)/NU
430 IF (RE(I)<2000) GOTO 470
440 IF (RE(I)>4000) GOTO 490
450 DP(I)=.516*Q(I)^2.333*RO*LR(I)/(D(I)^5.333*NU^.333)
460 GOTO 500
470 DP(I)=1132000!*Q(I)*NU*RO*LR(I)/D(I)^4
480 GOTO 500
490 DP(I)=69*(KE/D(I)+1922*NU*D(I)/Q(I))^.25*Q(I)^2*
RO*LR(I)/D(I)^5
500 DPG(I)=9.810001*B(I)*L(I)*(1.293-RO)
510 DPS(I)=DP(I)+DPG(I)
520 DPS=DPS+DPS(I)
530 NEXT I
540 PRINT "DPS="DPS
550 PRINT "Якщо результат не оптимальний введiть QN, інакше 0"
560 INPUT Q
570 IF (Q=0) GOTO 600
580 QN=Q
590 GOTO 280
600 PRINT "PN="PN;"DPD="DPD;"QN="QN;"WN="WN;"NI="NI
610 PRINT "Результати розрахунку по напряму "A$
620 PRINT "W="W,"Сумарний перепад DPS="DPS
630 PRINT " L(I) Q(I) D(I) DP(I) DPG(I) DPS(I)"
640 FOR I=1 TO NI
650 PRINT J$(I);
660 PRINT USING " ##.# ##.## ##.## ### ### ###"; L(I),Q(I),D(I),DP(I),
DPG(I),DPS(I)
670 NEXT I
680 PRINT "Для продовження ввести будь-яку цифру крiм 0"
690 INPUT M
700 IF M=0 GOTO 840
710 PRINT
720 PRINT "Введiть порядковий номер дiлянки"
730 PRINT "кiнцевий тиск якої вас цiкавить"
740 INPUT K
750 X=0
760 FOR I=1 TO K
770 X=X+DPS(I)
780 A$=A$(I)
790 NEXT I
800 DP=DPD-X
810 PK=PN-X
820 PRINT "Дiлянка"A$;"PK="PK;"DPостаточ="DP
830 GOTO 130
840 END
Додаток Н
Результати проектного розрахунку внутрішніх газових мереж житлового будинку за програмою GMN5
Розрахунок внутрiшнiх газових мереж житлового будинку без номограм
1.57 2.12 2.71 3.2 3.59 3.9 4.1 5.1
Початковi данi
Ввести NU,RO,T,EPS,KE
NU= 1.366E-05 RO= .756 T= 293 EPS= .1 KE= .01
Ввести напрям руху газу
основний
Ввести PN,DPD,QN,WN,NI
1500,350,5,7,11
Ввести номер дiлянки ,L(I),A(I),Q(I),вертикальнiсть ділянки
11-12,4.3,25,10.83,0
10-11,1.8,25,10.83,-1
9-10,3.8,25,10.83,0
8-9,4.8,25,5.81,0
7-8,3,20,5.81,-1
6-7,3,20,4.84,-1
5-6,3,20,4.84,-1
4-5,3,20,3.97,-1
3-4,3,20,2.32,-1
2-3,1.2,450,2.32,0
1-2,.8,450,2.3,0
DPS= 457.437
Якщо результат не оптимальний введiть QN, інакше 0
3
DPS= 340.069
Якщо результат не оптимальний введiть QN, інакше 0
PN= 1500 DPD= 350 QN= 3 WN= 7 NI= 11
Результати розрахунку по напряму основному
W= 3.019932 Сумарний перепад DPS= 340.069
Ділянка L(I) Q(I) D(I) DP(I) DPG(I) DPS(I)
11-12 4.3 10.83 3.59 18 0 18
10-11 1.8 10.83 3.59 8 -9 -2
9-10 3.8 10.83 3.59 16 0 16
8-9 4.8 5.81 2.12 85 0 85
7-8 3.0 5.81 2.12 51 -16 35
6-7 3.0 4.84 2.12 37 -16 21
5-6 3.0 4.84 2.12 37 -16 21
4-5 3.0 3.97 2.12 26 -16 10
3-4 3.0 2.32 1.57 38 -16 22
2-3 1.2 2.32 1.57 69 0 69
1-2 0.8 2.30 1.57 45 0 45
Для продовження ввести будь-яку цифру крiм 0
Введiть порядковий номер дiлянки
кiнцевий тиск якої вас цiкавить
3
ДiлянкаPK= 1467.914 DPостаточ= 317.9136
Ввести напрям руху газу
Ввести PN,DPD,QN,WN,NI
1468,318,3,7,9
Ввести номер дiлянки ,L(I),A(I),Q(I),вертикальнiсть дiлянки
13-9,3.2,25,7,0
14-13,3,20,7,-1
15-14,3,20,6.02,-1
16-15,3,20,5.04,-1
17-16,3,20,3.92,-1
18-17,3,20,2.45,-1
19-18,1.2,450,2.45,0
20-19,.8,450,2.3,0
DPS= 345.4353
Якщо результат не оптимальний введiть QN, інакше 0
2
DPS= 223.1782
Якщо результат не оптимальний введiть QN, інакше 0
2.5
DPS= 261.377
Якщо результат не оптимальний введiть QN, інакше 0
Результати розрахунку по відгалуженню
W= 3.209932 Сумарний перепад DPS= 261.377
Ділянка L(I) Q(I) D(I) DP(I) DPG(I) DPS(I)
13-9 3.2 7.00 2.71 24 0 24
14-13 3.0 7.00 2.71 22 -16 6
15-14 3.0 6.02 2.12 55 -16 39
16-15 3.0 5.04 2.12 39 -16 24
17-16 3.0 3.92 2.12 25 -16 9
18-17 3.0 2.45 1.57 46 -16 30
19-18 1.2 2.45 1.57 85 0 85
20-19 0.8 2.30 1.57 45 0 45
Для продовження ввести будь-яку цифру крiм 0
? 0
Додаток П
Текст програми GMN6 для розрахунку газових мереж
житлового будинку при відомих діаметрах
10 OPEN "o",#1,"rez"
20 OPEN "o",#2,"rez1":OPEN "o",#3,"rez2"
30 PRINT "Розрахунок внутрiшнiх газових мереж житлового будинку":PRINT
40 PRINT "Дiаметри дiльниць вiдомi":PRINT
50 DIM DS(8),Q(30),L(30),D(30),DP(30),DPG(30),J$(30),DPS(30),RE(30)
60 DIM A(30),B(30),LR(30),H(30)
70 PRINT
80 PRINT "Початковi данi":PRINT
90 PRINT "Введiть NU,RO,KE"
100 INPUT NU,RO,KE
110 PRINT "NU="NU;"RO="RO;"KE="KE:PRINT
120 WRITE #3,NU,RO,KE
130 PRINT "Введiть напрям руху газу"
140 INPUT A$
150 PRINT "Введiть PN,DPдоп,NI"
160 INPUT PN,DPD,NI
170 PRINT
180 FOR I=1 TO NI
190 PRINT "Введiть по ходу газу номер дiльницi,L(I),A(I),Q(I),D(I),B(I)"
200 INPUT J$(I),L(I),A(I),Q(I),D(I),B(I)
210 LR(I)=L(I)*(1+A(I)/100)
220 NEXT I
230 PK=PN-DPD
240 DPS=0
250 FOR I=1 TO NI
260 RE(I)=.0354*Q(I)/D(I)/NU
270 IF (RE(I)<2000) GOTO 310
280 IF (RE(I)>4000) GOTO 330
290 DP(I)=.516*Q(I)^2.333*RO*LR(I)/(D(I)^5.333*NU^.333)
300 GOTO 340
310 DP(I)=1132000!*Q(I)*NU*RO*LR(I)/D(I)^4
320 GOTO 340
330 DP(I)=69*(KE/D(I)+1922*NU*D(I)/Q(I))^.25*Q(I)^2*RO*LR(I)/D(I)^5
340 DPG(I)=9.810001*B(I)*L(I)*(1.293-RO)
350 DPS(I)=DP(I)+DPG(I)
360 DPS=DPS+DPS(I)
370 NEXT I
380 PRINT "PN="PN;"DPдоп="DPD;"NI="NI:PRINT
390 PRINT "Результати розрахунку по напрямку "A$:PRINT
400 PRINT "Сумарний перепад DPS="DPS:PRINT
410 WRITE #2,A$,DPS,PN,DPD,NI
420 FOR I=1 TO NI
430 PRINT "Дiльниця "J$(I);" ";"L="L(I);"a="A(I);"Q="Q(I);"D="D(I)
440 PRINT " DP="DP(I);"DPG="DPG(I);"DPS="DPS(I)
450 WRITE #1,J$(I),L(I),A(I),Q(I),D(I),DP(I),DPG(I),DPS(I)
460 NEXT I
470 PRINT "Для розрахунку вiдгалужень ввести будь-яку цифру окрiм 0"
480 INPUT M
490 IF M=0 GOTO 630
500 PRINT
510 PRINT "Введiть порядковий номер дiльницi,"
520 PRINT "кiнцевий тиск якоi вас цiкавить"
530 INPUT K
540 X=0
550 FOR I=1 TO K
560 X=X+DPS(I)
570 A1$=J$(I)
580 NEXT I
590 DP=DPD-X
600 PK=PN-X
610 PRINT "Дiльниця"A1$;"PK="PK;"DPостаточ="DP
620 GOTO 130
630 INPUT "Для друку результатiв -1";M1
640 IF M1<>1 THEN GOTO 820
650 CLOSE
660 OPEN "i",#1,"rez":OPEN "i",#2,"rez1":OPEN "i",#3,"rez2"
670 INPUT #3,NU,RO,KE
680 LPRINT "Розрахунок внутрiшнiх газових мереж житлового будинку":LPRINT
690 LPRINT "Дiаметри дiльниць вiдомi":LPRINT
700 LPRINT "Початковi данi":LPRINT
710 LPRINT "NU="NU;"RO="RO;"KE="KE:LPRINT
720 WHILE NOT(EOF(2))
730 INPUT #2,A$,DPS,PN,DPD,N
740 LPRINT "PN="PN;"DPдоп="DPD;"NI="N:LPRINT
750 LPRINT "Результати розрахунку по напрямку "A$:LPRINT
760 LPRINT "Сумарний перепад DPS="DPS:LPRINT
770 FOR I=1 TO N
780 INPUT #1,J$,L,A,Q,D,DP,DPG,DPS
790 LPRINT "Дiльниця "J$;" ";"L="L;"a="A;"Q="Q;"D="D
800 LPRINT " DP="DP;"DPG="DPG;"DPS="DPS
810 NEXT I:LPRINT
820 WEND
830 CLOSE:END
Додаток Р
Результати розрахунку внутрішніх газових мереж житлового будинку при відомих діаметрах за програмою GMN6
Розрахунок внутрiшнiх газових мереж житлового будинку
Дiаметри дiлянок вiдомi
Початковi данi
Введiть NU,RO,KE
NU= 1.366E-05 RO= .756 KE= .01
Введiть напрям руху газу
Основний
Введiть PN,DPдоп,NI
1500,350,11
Введiть по ходу газу номер дiлянки,L(I),A(I),Q(I),D(I),B(I)
11-12,4.3,25,10.83,3.59,0
10-11,1.8,25,10.83,3.59,-1
9-10,3.8,25,10.83,3.59,0
8-9, 4.8,25,5.81,2.12,0
7-8,3,20,5.81,2.12,-1
6-7,3,20,4.84,2.12,-1
5-6,3,20,4.84,2.12,-1
4-5,3,20,3.97,2.12,-1
3-4,3,20,2.32,1.57,-1
2-3,1.2,450,2.32,1.57,0
1-2,.8,450,2.3,1.57,0
DPS= 340.069
PN= 1500 DPдоп= 350 NI= 11
Результати розрахунку по напрямку основному
Сумарний перепад DPS= 340.069
Дiлянка L(I) A(I) Q(I) D(I) DP(I) DPG(I) DPS()
11-12 4.3 25 10.83 3.59 18 0 18
10-11 1.8 25 10.83 3.59 8 -9 -2
9-10 3.8 25 10.83 3.59 16 0 16
8-9 4.8 25 5.81 2.12 85 0 85
7-8 3.0 20 5.81 2.12 51 -16 35
6-7 3.0 20 4.84 2.12 37 -16 21
5-6 3.0 20 4.84 2.12 37 -16 21
4-5 3.0 20 3.97 2.12 26 -16 10
3-4 3.0 20 2.32 1.57 38 -16 22
2-3 1.2 450 2.32 1.57 69 0 69
1-2 0.8 450 2.30 1.57 45 0 45
Для розрахунку вiдгалужень ввести будь-яку цифру окрiм 0
1
Введiть порядковий номер дiлянки,
кiнцевий тиск якої вас цiкавить
3
ДiлянкаPK= 1467.914 DPостаточ= 317.9136
Введiть напрям руху газу
Введiть PN,DPдоп,NI
1468,318,8
Введiть по ходу газу номер дiлянки,L(I),A(I),Q(I),D(I),B(I)
13-9,3.2,25,7,2.71,0
14-13,3,20,7,2.71,-1
15-14,3,20,6.02,2.12,-1
16-15,3,20,5.04,2.12,-1
17-16,3,20,3.92,2.12,-1
18-17,3,20,2.45,1.57,-1
19-18,1.2,450,2.45,1.57,0
20-19,.8,450,2.3,1.57,0
PN= 1468 DPдоп= 318 NI= 8
Результати розрахунку по відгалуженню
Сумарний перепад DPS= 261.377
Дiлянка L(I) A(I) Q(I) D(I) DP(I) DPG(I) DPS()
13-9 3.2 25 7.00 2.71 24 0 24
14-13 3.0 20 7.00 2.71 22 -16 6
15-14 3.0 20 6.02 2.12 55 -16 39
16-15 3.0 20 5.04 2.12 39 -16 24
17-16 3.0 20 3.92 2.12 25 -16 9
18-17 3.0 20 2.45 1.57 46 -16 30
19-18 1.2 450 2.45 1.57 85 0 85
20-19 0.8 450 2.30 1.57 45 0 45
Для розрахунку вiдгалужень ввести будь-яку цифру окрiм 0
? 0
Ok_
Додаток С
Текст програми GMV1 для проектного розрахунку
розгалуженої газової мережі високого і середнього тиску
10 PRINT "GMV1"
20 PRINT "Аналiтичний розрахунок мережi середнього тиску"
30 DIM DD(20),Q(20),L(20),DS(20),PM(20),PMR(20)
40 DIM P1(20),A(20),D(20),P2(20),KP2(20)
50 DATA 3.2,3.59,3.9,4.1,5.1,7,8.3,10.2,15,20.9,26.3,31.5,41.4,51.8
60 FOR J=1 TO 14
70 READ DS(J)
80 PRINT DS(J);
90 NEXT J:PRINT
100 PRINT "Ввести PN,PK,NU,RO,T,W,KE,EPS,N"
110 INPUT PN,PK,NU,RO,T,W,KE,EPS,N
120 PRINT PN,PK,NU,RO,T,W,KE,EPS,N
130 PRINT "Ввести J$(I),L(I),Q(I)"
140 FOR I=1 TO N
150 INPUT J$(I),L(I),Q(I)
160 PRINT J$(I),L(I),Q(I)
170 NEXT I
180 PRINT " Результати розрахунку"
190 PRINT
200 PRINT
210 PRINT " Дiлянка Дiаметр Витрата Тиск почат. Тиск кiнц."
220 PRINT
230 PRINT
240 P1(1)=PN
250 IF(N=1)GOTO 280
260 N1=N-1
270 GOTO 290
280 N1=N
290 PRINT " Швидкiсть руху газу w="W
300 PRINT "Ввести напрям руху газу"
310 INPUT A$
320 PRINT " "A$
330 PRINT " Варiант N1="N1
340 FOR I=1 TO N1
350 PM(I)=P1(I)
360 DD(I)=.036238*SQR(Q(I)*T/PM(I)/W)
370 IF(DD(I)-3.2)>0 GOTO 390
380 D(I)=3.2:GOTO 450
390 FOR J=1 TO 14
400 IF(DD(I)-DS(J))>=0 GOTO 420
410 DZ=DS(J):GOTO 440
420 NEXT J
430 D(I)=D(10):GOTO 450
440 D(I)=DZ
450 A(I)=.000014*(KE/D(I)+1922*NU*D(I)/Q(I))^(.25)*Q(I)^2
*RO/D(I)^5
460 KP2(I)=P1(I)^2-A(I)*L(I)*1.1
470 IF(KP2(I)<=0)THEN PRINT "Збiльшiть швидкiсть":GOTO 1030
480 P2(I)=SQR(KP2(I))
490 PMR(I)=2/3*(P1(I)+P2(I)^2/(P1(I)+P2(I)))
500 IF(PM(I)-PMR(I))<=EPS GOTO 520
510 PM(I)=PMR(I):GOTO 360
520 PRINT J$(I),D(I),Q(I),P1(I),P2(I)
530 P1(I+1)=P2(I)
540 NEXT I
550 IF(N1=1)GOTO 580
560 N2=N1+1
570 GOTO 610
580 N2=1
590 P1(N2)=P1(N1)
600 GOTO 620
610 P1(N2)=P2(N1)
620 FOR I=N2 TO N
630 PM(I)=P1(I)
640 DD(I)=.036238*SQR(Q(I)*T/PM(I)/W)
650 IF(DD(I)-5.1)>0 GOTO 670
660 D(I)=5.1:GOTO 730
670 FOR J=1 TO 10
680 K=10+1-J
690 IF(DD(I)-DS(K))<=0 GOTO 710
700 DZ=DS(K):GOTO 720
710 NEXT J
720 D(I)=DZ
730 A(I)=.000014*(KE/D(I)+1922*NU*D(I)/Q(I))^(.25)*Q(I)^2
*RO/D(I)^5
740 KP2(I)=P1(I)^2-A(I)*L(I)*1.1
750 IF(KP2(I))<=0 GOTO 810
760 P2(I)=SQR(KP2(I))
770 PMR(I)=2/3*(P1(I)+P2(I)^2/(P1(I)+P2(I)))
780 IF(PM(I)-PMR(I))<=EPS GOTO 800
790 PM(I)=PMR(I):GOTO 640
800 IF(P2(I)-PK)>=0 GOTO 820
810 PRINT "Округлення у меншу сторону неможливо I="I
820 P1(I+1)=P2(I)
830 PRINT J$(I),D(I),Q(I),P1(I),P2(I)
840 NEXT I
850 IF(N2-1)=0 GOTO 900
860 IF(P2(N)-PK)<=0 GOTO 940
870 IF(N1-1)>0 GOTO 890
880 N2=1:GOTO 620
890 N1=N1-1:GOTO 330
900 IF(P2(N)-PK)<=0 THEN PRINT "Варiант неможливий":
GOTO 940
910 IF(N1-1)>0 GOTO 930
920 IF(MM=0)GOTO 1030
930 PRINT
940 PRINT "Введўть PN,W,N"
950 INPUT PN,W,N
960 PRINT PN,W,N
970 PRINT "Введiть J$(I),L(I),Q(I)"
980 FOR I=1 TO N
990 INPUT J$(I),L(I),Q(I)
1000 PRINT J$(I),L(I),Q(I)
1010 NEXT I
1020 GOTO 180
1030 END
Додаток Т
Результати проектного розрахунку
розгалуженої газової мережі середнього тиску
за програмою GMV1
3.2 3.59 3.9 4.1 5.1 7 8.3 10.2 15 20.9 26.3 31.5 41.4 51.8
Ввести PN,PK,NU,RO,T,W,KE,EPS,N
.3 .2 1.366E-05 .756 280
15 .01 .2 4
Ввести J$(I),L(I), Q(I)
ГРС-1 500 12800
1-2 1200 11000
2-3 800 8800
3-РК 700 5200
Результати розрахунку
Дiлянка Дiаметр Витрата Тиск почат. Тиск кiнц.
Швидкiсть руху газу w= 15
Ввести напрям руху газу Основний
ГРС-1 41.4 12800 .3 .298238
1-2 31.5 11000 .298238 .2850843
2-3 31.5 8800 .2850843 .2792146
3-РК 20.9 5200 .2792146 .2633824
Варiант N1= 2
ГРС-1 41.4 12800 .3 .298238
1-2 31.5 11000 .298238 .2850843
2-3 20.9 8800 .2850843 .2317126
3-РК 20.9 5200 .2317126 .2123682
Варiант N1= 1
ГРС-1 41.4 12800 .3 .298238
ГРС-1 20.9 12800 .3 .2320789
Округлення у меншу сторону неможливо I= 2
1-2 20.9 11000 .2320789 .2850843
2-3 20.9 8800 .2850843 .2317126
3-РК 20.9 5200 .2317126 .2123682
Введіть PN,W,N
.298238 15 1
Введiть J$(I),L(I),Q(I)
1-ГРП 300 1800
Результати розрахунку
Дiлянка Дiаметр Витрата Тиск почат. Тиск кiнц.
Швидкiсть руху газу w= 15
Ввести напрям руху газу
? Відвід від точки 1
1-ГРП 15 1800 .298238 .2938967
1-ГРП 10.2 1800 .298238 .264872
Введіть PN,W,N
.2850843 15 1
Введiть J$(I),L(I),Q(I)
2-ХЗ 200 2200
Швидкiсть руху газу w= 15
Ввести напрям руху газу
? Відвід від точки 2
2-ХЗ 15 2200 .2850843 .2806105
2-ХЗ 10.2 2200 .2850843 .2503492
Введіть PN,W,N
.2317126 15 1
Введiть J$(I),L(I),Q(I)
3-РЗ,200,3600
3-РЗ 200 3600
Швидкiсть руху газу w= 15
Ввести напрям руху газу
Відвід від точки 3
3-РЗ 20.9 3600 .2317126 .22911
3-РЗ 15 3600 .2317126 .2169285
Додаток У
Текст програми GMV2 для розрахунку газової мережі
високого і середнього тиску при відомих діаметрах
10 PRINT "РОЗРАХУНОК ГАЗОВОI МЕРЕЖI СЕРЕДНЬОГО ТИСКУ ПIСЛЯ НОМОГРАМ "
20 DIM J$(30),A$(30),D(30),Q(30),L(30),P1(30),P2(30),DS(30), KA(30)
30 PRINT "Ввести NU,RO,PN,PK,KE,N"
40 INPUT NU,RO,PN,PK,KE,N
50 PRINT NU,RO,PN,PK,KE,N
60 DATA 51.8,41.4,31.5,26.3,20.9,15,10.2,8.3,7,5.1
70 FOR R=1 TO 10
80 READ DS(R):PRINT DS(R);
90 NEXT R:PRINT
100 PRINT "Ввести номер дiлянки,Q(I),L(I),D(I)"
110 FOR I=1 TO N
120 INPUT A$(I),Q(I),L(I),D(I)
130 PRINT A$(I),Q(I),L(I),D(I)
140 NEXT I
150 PRINT "*****************************************************
160 PRINT " "
170 PRINT " "
180 PRINT " РОЗРАХУНОК ГАЗОВОI МЕРЕЖI "
190 PRINT" "
200 PRINT
210 PRINT"I I I I I I"
220 PRINT"I ДIЛЯНКА I ВИТРАТА I ВНУТРIШНIЙ IПОЧАТКОВИЙI КIНЦЕВИЙ I"
230 PRINT"I I ГАЗУ I ДIАМЕТР I ТИСК I ТИСК I"
240 PRINT"I___________Iкуб.м /год_I_____см_______I___МПа ___I____МПа___I"
250 PRINT "I I I I I I"
260 PRINT "Ввести напрям руху газу J$"
270 INPUT J$(I)
280 PRINT" " J$(I)
290 K=1
300 P1(1)=PN
310 FOR I=K TO N
320 KA(I)=P1(I)^2-.0000154*L(I)*(KE/D(I)+1922*NU*D(I)/Q(I))^.25*Q(I)^2*RO/D(I)^5
330 IF(KA(I)<0)GOTO 470
340 P2(I)=SQR(KA(I))
350 IF(P2(I)<PK)GOTO 470
360 PRINT " "A$(I),
370 PRINT USING " ##### I ##.# I #.#### I #.#### I";Q(I),D(I),P1(
I),P2(I)
380 P1(I+1)=P2(I)
390 NEXT I
400 FOR R=1 TO 10
400 FOR R=1 TO 10
410 IF(D(N)<=DS(R))GOTO 450
420 D(N)=DS(R)
430 K=N:NN=N
440 GOTO 310
450 NEXT R
460 PRINT "I___________I___________I______________I__________I__________I
470 PRINT "Дiаметр D("A$(I);")="D(I);"неможливий"
480 PRINT "Для продовження розрахунку ввести будь-яку цифру,iнакше 0"
490 INPUT Z
500 IF(Z=0)GOTO 600
510 PRINT "Ввести PN,N"
520 INPUT PN,N
530 PRINT PN,N
540 PRINT "Ввести номер дiлянки,Q(I),L(I),D(I)"
550 FOR I=1 TO N
560 INPUT A$(I),Q(I),L(I),D(I)
570 PRINT A$(I),Q(I),L(I),D(I)
580 NEXT I
590 GOTO 260
600 END
Додаток Ф
Результати розрахунку газової мережі середнього тиску при відомих діаметрах ділянок за програмою GMV2
РОЗРАХУНОК ГАЗОВОI МЕРЕЖI СЕРЕДНЬОГО ТИСКУ ПIСЛЯ НОМОГРАМ
Ввести NU,RO,PN,PK,KE,N
1.366E-05 .756 .3 .2 .01 4
51.8 41.4 31.5 26.3 20.9 15 10.2 8.3 7 5.1
Ввести номер дiлянки,Q(I),L(I),D(I)
ГРС-1 12800 500 41.4
1-2 11000 1200 31.5
1-3 8800 800 20.9
3-РК,5200,700,20.9
РОЗРАХУНОК ГАЗОВОI МЕРЕЖI
Ввести напрям руху газу J$
Основний
ГРС-1 12800 41.4 0.3000 0.2982
1-2 11000 31.5 0.2982 0.2851
1-3 8800 20.9 0.2851 0.2317
3-РК 5200 20.9 0.2317 0.2124
Дiаметр D(3-РК)= 15 неможливий
Для продовження розрахунку ввести будь-яку цифру,iнакше 0
Ввести PN,N
.2982 1
Ввести номер дiлянки,Q(I),L(I),D(I)
1-ГРП 1800 300 10.2
Ввести напрям руху газу J$
Вiдвiд вiд точки 1
1-ГРП 1800 10.2 0.2982 0.2648
Дiаметр D(1-ГРП)= 8.3 неможливий
Для продовження розрахунку ввести будь-яку цифру,iнакше 0
Ввести PN,N
.2851 1
Ввести номер дiлянки,Q(I),L(I),D(I)
2-ХЗ 2200 200 10.2
Ввести напрям руху газу J$
Вiдвiд вiд точки 2
2-ХЗ 2200 10.2 0.2851 0.2504
Дiаметр D(2-ХЗ)= 8.3 неможливий
Для продовження розрахунку ввести будь-яку цифру,iнакше 0
Ввести PN,N
.2317 1
Ввести номер дiлянки,Q(I),L(I),D(I)
3-РЗ 3600 200 15
Ввести напрям руху газу J$
Вiдвiд вiд точки 3
3-РЗ 3600 15.0 0.2317 0.2169
Дiаметр D(3-РЗ)= 10.2 неможливий
Для продовження розрахунку ввести будь-яку цифру,iнакше 0
0
Додаток Х
Текст програми GMVK для проектного розрахунку кільцевої газової мережі високого або середнього тиску
10 PRINT "Кiльцева мережа середнього тиску"
20 PRINT "Дiаметри дiлянок невiдомi"
30 PRINT "Ваше прiзвище"
40 INPUT N$
50 PRINT N$
60 DIM D(25,25), L(25,25), Q(25,25), AP(25,25), NM(25),A(25,25)
70 DIM PQ(25,25), SP(25), SPA(25),SPQ(25), SS(25), DE(25)
80 DIM DQK(25), DQS(25), DQ(25), J(25,25), RE(25,25), DS(25)
90 DATA 5.1,7.0,8.3,10.2,15,20.9,26.3,31.5,41.4,51.8
100 FOR R=1 TO 10
110 READ DS(R):PRINT DS(R);
120 NEXT R:PRINT
130 PRINT "Ввести NU,RO, EPS, NK, KE"
140 INPUT NU,RO,EPS,NK,KE
150 PRINT NU;RO;EPS;NK;KE
160 PRINT "Ввести T, PN, EPS1, LO, WM, QN, PK"
170 INPUT T,PN,EPS1,LO,WM,QN,PK
180 PRINT T;PN;EPS1;LO;WM;QN;PK
190 FOR K=1 TO NK
200 PRINT " NM(K) "
210 INPUT NM(K)
220 PRINT "NM(K)="NM(K)
230 PRINT "Ввести Q(K,I),L(K,I),J(K,I)"
240 FOR I=1 TO NM(K)
250 INPUT Q(K,I),L(K,I),J(K,I)
260 PRINT Q(K,I);L(K,I);J(K,I)
270 NEXT I
280 NEXT K
290 PM=.5*(PN+PK)
300 AO=(PN^2-PK^2)/(1.1*LO)
310 W=WM
320 DN=.036238*SQR(QN*T/PM/W)
330 AN=.0000154*(KE/DN+1922*NU*DN/QN)^.25*QN^2*
RO/DN^5
340 IF ((AN-AO-EPS1)>0) THEN W=W-.01:GOTO 320
350 FOR K=1 TO NK
360 FOR I=1 TO NM(K)
370 D(K,I)=.036238*SQR(ABS(Q(K,I))*T/PM/W)
380 FOR R=1 TO 10
390 IF ((D(K,I)-DS(R))>0) GOTO 430
400 DZ=DS(R)
410 D(K,I)=DZ
420 GOTO 440
430 NEXT R
440 NEXT I
450 NEXT K
460 FOR K=1 TO NK
470 SP(K)=0
480 SPA(K)=0
490 SPQ(K)=0
500 FOR I=1 TO NM(K)
510 AP(K,I)=.0000154*(KE/D(K,I)+1922*NU*D(K,I)/ ABS(Q(K,I)))^.25
520 A(K,I)=AP(K,I)*Q(K,I)*Q(K,I)*RO/D(K,I)^5
530 IF (Q(K,I)>0) GOTO 550
540 A(K,I)=-A(K,I)
550 SP(K)=SP(K)+A(K,I)*L(K,I)
560 SPA(K)=SPA(K)+ABS(A(K,I))*L(K,I)
570 NEXT I
580 NEXT K
590 FOR K=1 TO NK
600 DE(K)=200*SP(K)/SPA(K)
610 IF (ABS(DE(K))>EPS) GOTO 640
620 NEXT K
630 GOTO 980
640 FOR K=1 TO NK
650 FOR I=1 TO NM(K)
660 PQ(K,I)=A(K,I)/Q(K,I)*L(K,I)
670 SPQ(K)=SPQ(K)+PQ(K,I)
680 NEXT I
690 SS(K)=0
700 DQK(K)=-(SP(K))/(2*SPQ(K))
710 NEXT K
720 DQK(0)=0
730 FOR K=1 TO NK
740 FOR I=1 TO NM(K)
750 FOR M=1 TO NK
760 IF(M-J(K,I)=0)GOTO 780
770 NEXT M
780 SS(K)=SS(K)+PQ(K,I)*DQK(M)
790 NEXT I
800 DQS(K)=SS(K)/SPQ(K)
810 DQ(K)=DQK(K)+DQS(K)
820 NEXT K
830 FOR K=1 TO NK
840 FOR I=1 TO NM(K)
850 Q(K,I)=Q(K,I)+DQ(K)
860 NEXT I
870 NEXT K
880 DQ(0)=0
890 FOR K=1 TO NK
900 FOR I=1 TO NM(K)
910 FOR M=1 TO NK
920 IF(M-J(K,I)=0) GOTO 940
930 NEXT M
940 Q(K,I)=Q(K,I)-DQ(M)
950 NEXT I
960 NEXT K
970 GOTO 460
980 PRINT "Для друку ввести будь-яку цифру"
990 INPUT O
1000 PRINT "Результати розрахунку"
1000 PRINT "Результати розрахунку"
1010 PRINT "W="W
1020 FOR K=1 TO NK
1030 PRINT "Контур="K, "DE="DE(K)
1040 PRINT "D(K,I) Q(K,I) L(K,I) A(k,i)*L(k,i)"
1050 FOR I=1 TO NM(K)
1060 PRINT USING " ##.## ####.## ###. #.#######";D(K,I),Q(K,I),L(K,I),A(
K,I)*L(K,I)
1070 NEXT I:INPUT AZ
1080 NEXT K
1090 END
Додаток Ц
Результати розрахунку кільцевої газової мережі середнього тиску за програмою GMVK
3.2 3.59 3.9 4.1 5.1 7 8.3 10.2 15 20.9 26.3 31.5 41.4 51.8
Ввести NU,RO, EPS, NK, KE
1.366E-05 .756 1 5 .01
Ввести T, PN, EPS1, LO, WM, QN, PK
280 .3 .2 1320 15 50 .2
NM(K)= 8
Ввести Q(K,I),L(K,I),J(K,I)
120.5 230 0
98 220 0
76.5 210 0
56 200 0
-306 200 0
-79.5 210 3
-34.5 250 2
-11 220 2
NM(K)= 6
Ввести Q(K,I),L(K,I),J(K,I)
34.5 250 1
11 220 1
35 220 0
-11 220 3
-60 240 5
-12.5 250 4
NM(K)= 4
Ввести Q(K,I),L(K,I),J(K,I)
79.5 210 1
11 220 2
-195.5 210 0
-83 220 5
NM(K)= 4
Ввести Q(K,I),L(K,I),J(K,I)
12.5 250 2
12 240 0
-11.5 230 5
-11.5 230 0
NM(K)= 6
Ввести Q(K,I),L(K,I),J(K,I)
83 220 3
60 240 2
11.5 230 4
-79.5 230 0
-57 220 0
-34.5 230 0
Для друку ввести будь-яку цифру
1
Результати розрахунку
W= 15
Контур= 1 DE= .1333371
D(K,I) Q(K,I) L(K,I) A(k,i)*L(k,i)
3.59 98.52 230. 0.0107812
3.20 76.02 220. 0.0112515
3.20 54.52 210. 0.0056614
3.20 34.02 200. 0.0021970
7.00 -327.98 200. -.0031472
3.20 -101.40 210. -.0187871
3.20 -52.09 250. -.0061756
3.20 -28.59 220. -.0017414
Контур= 2 DE= .1782474
D(K,I) Q(K,I) L(K,I) A(k,i)*L(k,i)
3.20 52.09 250. 0.0061756
3.20 28.59 220. 0.0017414
3.20 30.62 220. 0.0019806
3.20 -15.31 220. -.0005451
3.20 -63.77 240. -.0087432
3.20 -14.93 250. -.0005917
Контур= 3 DE= 8.136181E-02
D(K,I) Q(K,I) L(K,I) A(k,i)*L(k,i)
3.20 101.40 210. 0.0187871
3.20 15.31 220. 0.0005451
5.10 -195.57 210. -.0061468
3.20 -82.46 220. -.0131696
Контур= 4 DE= .6614364
D(K,I) Q(K,I) L(K,I) A(k,i)*L(k,i)
3.20 14.93 250. 0.0005917
3.20 10.05 240. 0.0002754
3.20 -12.84 230. -.0004124
3.20 -13.45 230. -.0004491
Контур= 5 DE= .1535313
D(K,I) Q(K,I) L(K,I) A(k,i)*L(k,i)
3.20 82.46 220. 0.0131696
3.20 63.77 240. 0.0087432
3.20 12.84 230. 0.0004124
3.20 -80.11 230. -.0130177
3.20 -57.61 220. -.0065920
3.20 -35.11 230. -.0026813
Додаток Ш
Текст програми GMVR для визначення пропускної здатності розгалуженої газової мережі високого або середнього тиску
10 PRINT "Розрахунок пропускної здатностi розгалуженої газової мережi"
20 PRINT "високого i середнього тиску за основною формулою"
30 DATA 11e-6,.562,280,1e-4,.4,.2,7
40 PRINT " ET DE T KE PPS PKS N"
50 READ ET,DE,T,KE,PPS,PKS,N
60 PRINT ET;DE;T;KE,PPS,PKS,N
70 PRINT "WMA DQ"
80 DATA 15,10
90 READ WMA,DQ
100 PRINT WMA,DQ:PRINT
110 PRINT "Параметри основной магiстралi":PRINT
120 DATA.315,200,.315,100,0.315,300,0.263,370,0.263,280,
0.209, 590,0.209,710
130 PRINT " i D,m l,m"
140 FOR I=1 TO N
150 READ D(I),L(I)
160 PRINT USING " ## #.### #### ";I,D(I),L(I)
170 NEXT I:PRINT
180 PRINT "Кiлькicть дiлянок на магiстральних вiдводах"
190 DATA 1,3,1,1,2,1
200 FOR I=1 TO N-1
210 READ K(I)
220 PRINT " ";K(I);
230 NEXT I:PRINT:INPUT AZ
240 PRINT "Параметри вiдводiв":PRINT
250 PRINT " i j DV,m lV,m"
260 DATA.083,350,0.209,270,0.15,430,0.102,310,.07,210,.209,
400,.15,100,.15,360,.102,320
270 FOR I=1 TO N-1
280 FOR J=1 TO K(I)
290 READ DV(I,J),LV(I,J)
300 PRINT USING " ## ## #.### ####";I,J,DV(I,J),LV(I,J)
310 NEXT J
320 NEXT I:PRINT:INPUT AZ
330 PRINT "Параметри вiдводiв вiд вiдводiв"
340 PRINT " i j DW,m lW,m"
350 DATA .209,290,.07,250,.15,340
360 FOR I=1 TO N-1
370 IF K(I)=1 GOTO 420
380 FOR J=1 TO K(I)-1
390 READ DW(I,J),LW(I,J)
400 PRINT USING " ## ## #.### ####";I,J,DW(I,J),LW(I,J)
410 NEXT J
420 NEXT I:PRINT
430 PRINT "Схема включення вiдвoдiв":PRINT
440 PRINT
450 DATA 1,0,0,0,0,0
460 FOR I=1 TO N-1
470 READ S(I)
480 PRINT " ";S(I);
490 NEXT I
500 PRINT
510 PRINT "Схема включення вiдвoдiв вiд вiдвoдiв":PRINT
520 PRINT
530 DATA 0,0,0
540 FOR I=1 TO N-1
550 IF K(I)=1 GOTO 610
560 IF S(I)=0 THEN GOTO 610
570 FOR J=1 TO K(I)-1
580 READ SS(I,J)
590 PRINT " ";SS(I,J);
600 NEXT J
610 NEXT I:PRINT:INPUT AZ
620 PI=3.14159:DQ=DQ/3600
630 AL=5500000!*DE^1.3/T^3.3
640 FOR I=1 TO N
650 KA(I)=1.1*DE*T*7.781E-10*L(I)/D(I)^5
660 QMA(I)=PI*D(I)^2/4*WMA
670 NEXT I
680 FOR I=1 TO N-1
690 FOR J=1 TO K(I)
700 KAV(I,J)=1.1*DE*T*7.781E-10*LV(I,J)/DV(I,J)^5
710 QMAV(I,J)=PI*DV(I,J)^2/4*WMA
720 IF J=K(I) THEN GOTO 760
730 KAW(I,J)=1.1*DE*T*7.781E-10*LW(I,J)/DW(I,J)^5
740 QMAW(I,J)=PI*DW(I,J)^2/4*WMA
750 NEXT J
760 NEXT I
770 Q(1)=DQ
780 PP(1)=PPS
790 FOR I=1 TO N
800 Z(I)=1-AL*PP(I)
810 RE(I)=1.6456*DE*Q(I)/(D(I)*ET)
820 LA(I)=.11*(KE/D(I)+68/RE(I))^.25
830 PK(I)=SQR(PP(I)^2-KA(I)*LA(I)*Z(I)*Q(I)^2)
840 IF I=N GOTO 1110
850 IF S(I)=0 THEN QV(I,J)=0:GOTO 1080
860 IF QV(I,1)=0 THEN QV(I,1)=DQ
870 PPV(I,1)=PK(I)
880 FOR J=1 TO K(I)
890 IF QV(I,1)=>Q(I) THEN GOTO 1110
900 ZV(I,J)=1-AL*PPV(I,J)
910 REV(I,J)=1.6456*DE*QV(I,J)/(DV(I,J)*ET)
920 LAV(I,J)=.11*(KE/DV(I,J)+68/REV(I,J))^.25
930 PKV(I,J)=SQR(PPV(I,J)^2-KAV(I,J)*LAV(I,J)*ZV(I,J)*QV(I,J)^2)
940 IF(K(I)=J) GOTO 1070
950 IF (SS(I,J)=0) THEN QW(I,J)=0:GOTO 1040
960 IF QW(I,J)=0 THEN QW(I,J)=DQ
970 PPW(I,J)=PKV(I,J)
980 IF QW(I,J)=>QV(I,J) THEN GOTO 1070
990 ZW(I,J)=1-AL*PPW(I,J)
1000 REW(I,J)=1.6456*DE*QW(I,J)/(DW(I,J)*ET)
1010 LAW(I,J)=.11*(KE/DW(I,J)+68/REW(I,J))^.25
1020 PKW(I,J)=SQR(PPW(I,J)^2-KAW(I,J)*LAW(I,J)*ZW(I,J)*QW(I,J)^2)
1030 IF PKW(I,J)>PKS THEN QW(I,J)=QW(I,J)+DQ:GOTO 980
1040 QV(I,J+1)=QV(I,J)-QW(I,J):PPV(I,J+1)=PKV(I,J)
1050 IF QV(I,J+1)<=0 THEN GOTO 1110
1060 NEXT J
1070 IF PKV(I,J)>PKS THEN QV(I,1)=QV(I,1)+DQ:GOTO 880
1080 Q(I+1)=Q(I)-QV(I,1):PP(I+1)=PK(I)
1090 IF Q(I+1)<=0 THEN GOTO 1110
1100 NEXT I
1110 IF PK(I)>PKS THEN Q(1)=Q(1)+DQ:GOTO 790
1120 PRINT "Параметри роботи основноi магiстралi":PRINT
1130 PRINT " i D,m l,m Q,m^3/год Qmax,m^3/год"
1140 FOR I=1 TO N
1150 PRINT USING "## #.### #### ###### ######"; I,D(I),L(I),Q(I)*3600,QMA(I)*3600
1160 NEXT I:PRINT:INPUT AZ
1170 PRINT "Параметри роботи вiдводiв":PRINT
1180 PRINT " i j DV,m lV,m QV,m^3/год QVmax,m^3/год"
1190 FOR I=1 TO N-1
1200 FOR J=1 TO K(I)
1200 FOR J=1 TO K(I)
1210 PRINT USING"## ## #.### #### ###### ######"; I,J,DV(I,J),LV(I,J),QV(I,J)*3600,QMAV(I,J)*3600
1220 NEXT J
1230 NEXT I: PRINT:INPUT AZ
1240 PRINT "Параметри роботи вiдводiв вiд вiдводiв": PRINT
1250 PRINT " i j DW,m lW,m QW,m^3/год QWmax,m^3/год"
1260 FOR I=1 TO N-1
1270 IF K(I)=1 THEN GOTO 1310
1280 FOR J=1 TO K(I)-1
1290 PRINT USING "## ## #.### #### ###### ######";I,J,DW(
I,J),LW(I,J),QW(I,J)*3600,QMAW(I,J)*3600
1300 NEXT J
1310 NEXT I
1320 END
Додаток Щ
Результати розрахунку за програмою GMVR розгалуженої газової мережі середнього тиску
Розрахунок пропускної здатностi розгалуженої газової мережi
високого i середнього тиску за основною формулою
ET DE T KE PPS PKS N
.0000102 .585 280 .0001 .3 .2
4
WMA DQ
15 5
Параметри основної магiстралi
i D,m l,m
1 0.414 500
2 0.315 1200
3 0.209 800
4 0.209 700
Кiлькicть дiлянок на магiстральних вiдводах
1 1 1
Параметри вiдводiв
i j DV,m lV,m
1 1 0.102 300
2 1 0.102 200
3 1 0.150 200
Параметри вiдводiв вiд вiдводiв
i j DW,m lW,m
Схема включення вiдвoдiв
0 0 0
Параметри роботи основної магiстралi
i D,m l,m Q,m^3/год Qmax,m^3/год
1 0.414 500 8160 7269
2 0.315 1200 8160 4208
3 0.209 800 8160 1853
4 0.209 700 8160 1853
Параметри роботи вiдводiв
i j DV,m lV,m QV,m^3/год QVmax,m^3/год
1 1 0.102 300 0 441
2 1 0.102 200 0 441
3 1 0.150 200 0 954
Схема включення вiдвoдiв
1 0 0
Параметри роботи основної магiстралi
i D,m l,m Q,m^3/год Qmax,m^3/год
1 0.414 500 11065 7269
2 0.315 1200 8130 4208
3 0.209 800 8130 1853
4 0.209 700 8130 1853
Параметри роботи вiдводiв
i j DV,m lV,m QV,m^3/год QVmax,m^3/год
1 1 0.102 300 2935 441
2 1 0.102 200 0 441
3 1 0.150 200 0 954
Схема включення вiдвoдiв
1 1 0
Параметри роботи основної магiстралi
i D,m l,m Q,m^3/год Qmax,m^3/год
1 0.414 500 14215 7269
2 0.315 1200 11280 4208
3 0.209 800 7720 1853
4 0.209 700 7720 1853
Параметри роботи вiдводiв
i j DV,m lV,m QV,m^3/год QVmax,m^3/год
1 1 0.102 300 2935 441
2 1 0.102 200 3560 441
3 1 0.150 200 0 954
Схема включення вiдвoдiв
1 1 1
Параметри роботи основної магiстралi
i D,m l,m Q,m^3/год Qmax,m^3/год
1 0.414 500 16250 7269
2 0.315 1200 13315 4208
3 0.209 800 9755 1853
4 0.209 700 2845 1853
Параметри роботи вiдводiв
i j DV,m lV,m QV,m^3/год QVmax,m^3/год
1 1 0.102 300 2935 441
2 1 0.102 200 3560 441
3 1 0.150 200 6910 954