*'Количество ветвей в схеме нулевой последовательности KHSNP=',I2/
*'Количество узлов в схеме нулевой последовательности,'/
*' которые имеют нулевой потенциал KWSNPNP=',I2/
*'Базисная мощность SB=',F6.0,' MVA'/
*'Информация о ветвях схемы прямой последовательности:'/
*' N1K(KH) N2K(KH) X(KH) SNG(KH)'/
*' отн.ед. МВА'/
*(2(5X,I2),F11.4,F12.0))
WRITE(3,321)(NWKZ(I),UB(I),I=1,KWKZ)
321 FORMAT('Информация о узлах короткого замыкания:'/
*' NWKZ(KWKZ) UB(KWKZ)'/
*' кв'/(I7,F13.1))
IF(NSNP.EQ.0) GO TO 324
WRITE(3,322)(N1KSNP(I),N2KSNP(I),XSNP(I),I=1,KHSNP)
322 FORMAT('Информация о ветвях схемы нулевой последовательности:'/
*' N1KSNP(KHSNP) N2KSNP(KHSNP) XSNP(KHSNP)'/
*' отн.ед.'/(2И12,F16.4))
WRITE(3,323)(NWSNPNP(I),I=1,KWSNPNP)
323 FORMAT('Узлы схемы нулевой последовательности, которые имеют '/
*'нулевой потенциал: NWSNPNP(KWSNPNP)'/(20(2X,I2)))
324 WRITE(3,325)
325 FORMAT(/' РЕЗУЛЬТАТЫ РАСЧЕТА:')
C Блок 3: Переход от сопротивлений к проводимостям ветвей СПП в том
C же массиве X(KH) и образование массива E(KH) токов источников тока; пеC реход к относительным значимостям мощностей (токов) в массиве SNG(KH);
DO 32 I=1,KH
X(I)=1./X(I)
SNG(I)=SNG(I)/SB
IF(N1K(I).EQ.0)GO TO 181
IF(N2K(I).NE.0)GO TO 180
181 E(I)=1.*X(I)
GO TO 32
180 E(I)=0.
32 CONTINUE
C Блок 4: Расчет Количества линейно независимых Узлов в СПП (переC менная KW); образование массива Нумеров в этих Узлах (массив NW(KW));
KW=0
DO 36 I=1,99
M=0
DO 35 K=1,KH
IF(I.EQ.N1K(K)) GO TO 34
IF(I.NE.N2K(K)) GO TO 35
34 M=1
35 CONTINUE
IF(M.EQ.0) GO TO 36
KW=KW+1
IF(KW.LE.70) GO TO 245
WRITE(5,244)
244 FORMAT(1X,'Размеры массивов NW(70),KWZUU(70),UDAV(70),',
*'UUAV(70),UAV(70) недостаточные для вашей задачи'/1X,
*'ПРОГРАММА TKZ ЗАВЕРШАЕТ РАБОТУ'/)
STOP
245 NW(KW)=I
36 CONTINUE
C Блок 5: Вводилась ли СНП в исходных данных?
IF(NSNP.EQ.0) GO TO 45
C Блок 6: Размножение информации о номерах концов ветвей СНП
C Переход от сопротивлений к проводимостям ветвей СНП
DO 38 I=1,KHSNP
N1KSNP1(I)=N1KSNP(I)
N2KSNP1(I)=N2KSNP(I)
38 XSNP(I)=1./XSNP(I)
C Блок 7: Засылка нулей на место тех номеров узлов в массивах
C N1KSNP(KHSNP) и N2KSNP(KHSNP), которые имеют нулевой потенциал
DO 41 I=1,KWSNPNP
DO 41 K=1,KHSNP
IF(NWSNPNP(I).NE.N1KSNP1(K)) GO TO 40
N1KSNP1(K)=0
GO TO 41
40 IF(NWSNPNP(I).NE.N2KSNP1(K)) GO TO 41
N2KSNP1(K)=0
41 CONTINUE
C Блок 8: Расчет Количества линейно независимых Узлов в СНП
C (KWSNP) и создание массива номеров этих узлов (NWSNP(KWSNP))
KWSNP=0
DO 44 I=1,99
M=0
DO 43 K=1,KHSNP
IF(I.EQ.N1KSNP1(K)) GO TO 42
IF(I.NE.N2KSNP1(K)) GO TO 43
42 M=1
43 CONTINUE
IF(M.EQ.0) GO TO 44
KWSNP=KWSNP+1
NWSNP(KWSNP)=I
44 CONTINUE
45 CONTINUE
C Блок 9: Обнуление KWKZB (Количество Узлов Краткого
C Замыкания Текущая)
KWKZB=0
C Блок 10: Начало циклической части программы TKZ, которая повторяется
C для каждого узла краткого замыкания; наращивание переменной
C KWKZB; выборка Номера Узла Краткого замыкания (переменная
C NWK) из массива NWKZ(KWКZ)
47 KWKZB=KWKZB+1
NWK=NWKZ(KWKZB)
C Блок 11: Переписывание информации из массивов N1K(KH),N2K(KH),
C X(KH),E(KH) в соответствующие рабочие массивы
DO 48 I=1,KH
N1K1(I)=N1K(I)
N2K1(I)=N2K(I)
X1(I)=X(I)
48 E1(I)=E(I)
C Блок 12: Образование массива Номеров Узлов для расчета узловых
C Напряжений NWU(KWU), в котором бегущий узел короткого замыкания
C NWK в различие от массива NW(KW) записан последним;
C переменная KWU=KW-1
KWU=0
DO 49 I=1,KW
IF(NW(I).EQ.NWK) GO TO 49
KWU=KWU+1
NWU(KWU)=NW(I)
49 CONTINUE
NWU(KW)=NWK
C Блок 13: Обнуление массивов для расчета узловых напряжений в СПП
DO 50 I=1,KW
UDAV(I)=0.
UUAV(I)=0.
50 UAV(I)=0.
C Блок 14: переменная KHB - Количество Ветвей Текущая
KHB=KH
C Блок 15: Подпрограмма ZHORTKA (Свертка)
CALL ZHORTKA(N1K1,N2K1,X1,E1,NWU,KWU,
*KHB,KYIJ,KWZUU,YSUMU,AJU,NUSUU,YIJU)
C Блок 16: Расчет результирующих параметров трехфазного и
C двухфазного КЗ для текущего узла короткого замыкания NWK:
C ERES - результирующий ток источника тока относительно узла КЗ (суммарный
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.