Расчет токов несимметрических коротких замыканий на ЭВМ (Лабораторная работа № 6), страница 9

106 FORMAT(I12,I5,(3F12.4))

107 CONTINUE

IF(NRKRTKZH.EQ.2) GO TO 109

WRITE(3,202)

202 FORMAT(' Коэффициенты распределения полного тока КЗ в',

*' поврежденных фазах по'/'   ветвям схемы при разных',

*' видах КЗ в текущем узле КЗ (положительным'/'   направлениям',

*' в ветвях считается направление от конца ветви'/

*'   с большим номером к концу с меньшим номером)'/

*'   Ветвь схемы  Трехфазной  Двухфазной  Однофазной',

*' Двухф.на землю')

DO 197 I=1,KHSOKR

IF(NSNP.NE.0.AND.NSOWP.NE.0) GO TO 198

WRITE(3,199)N1K1(I),N2K1(I),AKP(I),AKO(I)

199 FORMAT(I12,I5,F12.4,F12.4)

GO TO 197

198 WRITE(3,200)N1K1(I),N2K1(I),AKP(I),AKO(I),TW1(I),TW11(I)

200 FORMAT(I9,I5,F12.4,F12.4,F12.4,F12.4)

197 CONTINUE

C Блок 34: Вывод сообщения на экран дисплея

109 WRITE(5,400)NWK

400 FORMAT(1X,' Рассчитанные токи КЗ для узла ',I2)

C Блок 35: Все ли узлы КЗ?

IF(KWKZB.LT.KWKZ) GO TO 47

C Блок 36: Закрытие файла выходных данных на диске ЭВМ

CLOSE(UNIT=3)

C Блок 37: Выдача окончательного сообщения на экран

WRITE(5,390)Z2

390 FORMAT(1X,'Программа TKZ работу закончила. Результаты',

*' расчета записаны на диск в файл '1X,A8)

STOP

END

SUBROUTINE ZHORTKA(N1K,N2K,Y,AJ,NW,KW,KH,KYIJ,

*KHMPZWW,YSMPZWW,AJMPZWW,NWMPZWW,YHMPZWW)

C Подпрограмма свертки схемы к узлу, номер которого записан

C     последним в массиве NW(KW+1)

C Входные переменные и массивы подпрограммы

C   N1K(KH) - массив Номеров "Первых" Концов ветвей схемы

C   N2K(KH) - массив Номеров "Других" Концов ветвей схемы

C   Y(KH) - массив Проводимостей ветвей схемы (отн. ед.)

C   AJ(KH) - массив Токов источников тока ветвей схемы (отн. ед.)

C   NW(KW) - массив Номеров линейно независимых Узлов схемы;

C            последним в этом массиве записан номер узла КЗ

C   KW - Количество линейноо независимых Узлов схемы минус один

C   KH - Количество Ветвей схемы

C Выходные переменные и массивы подпрограммы

C   KYIJ - Количество ветвей с Проводимостями Yij, которые необходимо

C       запомнить при выполнении прямого хода методом Гаусса для того,

C       чтобы соблюсти расчет узловых напряжений при обратном

C       ходе (количество всех ненулевых элементов в верхней

C       треугольной матрице коэффициентов узловых уравнений

C       после окончания прямого хода)

C   KHMPZWW(KW) - массив Количеств Ветвей Многолучевых

C       Звезд Удаляемых Узлов, которые были подключенны к удаленному

C       узлу в момент его удаления (исп. при рассчете  напряжения.)

C   YSMPZWW(KW) - массив Сумм. Проводимостей ветвей Многолучевых

C       Звезд Удаляемого Узла, подключенных к удаляемым узлам

C       в момент их удаления (диагональный коэффициент системы узлоC       вых урвнений для уравнения, которое соответствует номеру удаляемого узла) (отн.ед.)

C   AJMPZWW(KW) - массив Tоков источников тока Многолучевых Звезд,

C       которые "нагнетали" ток в удаляемый узел в момент его

С       удаления (правая часть системы того узлового уравнения,

C       которое соответствует номеру удаляемого узла) (отн.ед.)

C   NWMPZWW(KYIJ) - массив Номеров Вершин Многолучевых

C       Звезд Удаляемых Узлов, на которые опирается звезда в

C       момент удаления этого узла (номера столбцов ненулевых

C       элементов в строчках верхней треугольной матрицы коэффициенC       тов системы узловых уравнений, которые соответствуют ноC       меру удаляемого узла)

C   YHMPZWW(KYIJ) - массив Проводимостей Ветвей Многолучевых

C       Звезд Удаляемых Узлов (ненулевые элементы тех строчек

C       верхней треугольной матрицы коэффициентов системы узловых

C       уравнений, которые соответствуют номеру удаляемого узла)

C       (отн.ед.)

C Рабочие массивы в подпрограмме ZHORTKA

C KHDNW(KW) - Количество Ветвей, Присоединенных к Неудаленным узлам,

C       начиная от текущего узла и кончая последним узлом

C       в списке NW(KW)

C NWMPZ(KHMPZ)-массив Номеров Узлов Многолучевой Звезды

C YHMPZ(KHMPZ)-массив Провод. Ветвей Многолучевой Звезды

C AJHMPZ(KHMPZ) - массив Токов источников тока в Ветвях

C        Многолучевой Звезды

DИMENSION N1K(100),N2K(100),Y(100),AJ(100),NW(70),

*KHMPZWW(70),YSMPZWW(70),AJMPZWW(70),NWMPZWW(300),