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

*E(100),NW(70),N1KSNP1(50),N2KSNP1(50),NWSNP(50),

*N1K1(100),N2K1(100),X1(100),E1(100),NWU(70),UDAV(70),UUAV(70),

*UAV(70),KWZUU(70),NUSUU(300),YSUMU(70),AJU(70),YИJU(300),

*AKP(100),AKP1(100),AKO(100),AKN(100),AKNP(100),TW1(100),TW11(100),

*NWSNPU(50),N1KSNP2(50),N2KSNP2(50)

C          Массивы, которые задаются во входных данных:

C N1K(100) - номера "первых" концов ветвей схемы прямой последов.(СПП)

C N2K(100) - номера "вторых" концов ветвей СПП

C X(100) – сопротивления ветвей СПП (относительные единицы)

C SNG(100) - мощности генераторов в ветвях СПП (MBA)

C NWKZ(20) - номера узлов короткого замыкания  (КЗ)

C UB(20) - напряжения базисные в узлах КЗ (кв)

C N1KSNP(50) - номера "первых" концов ветвей схемы нулев. последов.(СНП)

C N2KSNP(50) - номера "вторых" концов ветвей СНП

C XSNP(50) – сопротивления ветвей СНП (относительные единицы)

C NWSNPNP(30) - номера узлов СНП, которые имеют нулевой потенциал

C

C     Массивы, которые используются в главной программе TKZ как рабочие

C E(100) - массив токов источников тока в ветвях СПП (отн.ед.)

C NW(100) - номера узлов СПП

C N1KSNP1(50) - рабочий массив для работы с информацией массива N1KSNP

C N2KSNP1(50) - рабочий массив для работы из информацией массива N2KSNP

C N1K1(100),N2K1(100),X1(100),E1(100) - рабочие массивы

C UDAV(70) - узловые напряжения "Доаварийного" режима в СПП

C UUAV(70) - узловые напряжения "Собственно Аварийного" режима в СПП

C UAV(70) - узловые напряжения "Аварийного" режима в СПП

C Блок №1: Ввод входных данных

WRITE(5,350)

350 FORMAT(1X,'ВВЕДИТЕ ИМЯ (С РАСШИРЕНИЕМ) ВАШЕГО ФАЙЛА',

*' ИСХОДНЫХ ДАННЫХ В ВИДЕ'/

*1X,'AAAA.DAT (ГДЕ СИМВОЛЫ AAAA - ЭТО ИМЯ ФАЙЛА)')

READ(5,3)Z1

3 FORMAT(A8)

OPEN(UNIT=1,FILE=Z1,STATUS='OLD')

READ(1,*)NRKRTKZH,NSNP,KH,KWKZ,KHSNP,KWSNPNP,SB,

*(N1K(I),N2K(I),X(I),SNG(I),I=1,KH),

*(NWKZ(I),UB(I),I=1,KWKZ)

C           Простые переменные, которые вводятся в исходных данных:

C NRKRTKZH-Признак расчета коэффициентов распределения ТКЗ по ветвям

C     схемы;(NRKRTKZH=0 - программа выполняет расчет и вывод в выходной

C     файл значимостей симметрических составляющих и токов КЗ в узлах

C     КЗ; NRKRTKZH=1 -выполняется то, что и при NRKRTKZH=0, плюс

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

C     которые подключены к узлу КЗ; NRKRTKZH=2 - выполняется то, что

C     при NRKRTKZH=0, плюс расчет и выдача коэффициентов распределения

C     симметрических составляющих токов КЗ для всех ветвей схемы;

C     NRKRTKZH=3 - выполняется то, что при NRKRTKZH=2,

C     плюс расчет и выдача коэффициентов распределения

C     полных токов поврежденных фаз для всех ветвей схемы)

C NSNP-Признак Схемы Нулевой Последовательности (NSNP=0-СНП в исходных

c     данных не описывается и не вводится, расчет ТКЗ для однофазных

c     и двухфазных КЗ на землю не выполняется; NSNP=1 - СНП описывается

c     в исходных данных и выполняется расчет ТКЗ для однофазных

C     и двухфазных КЗ на землю)

C KH - Количество Ветвей СПП;

C KWKZ - Количество Узлов Короткого Замыкания;

C KHSNP - Количество Ветвей СНП (при NSNP=0 принимается равным нулю);

C KWSNPNP - Количество Узлов СНП, которые имеют Нулевой Потенциал

C      (при NSNP=0 принимается равным нолю);

C SB - базисная мощность (МВА);

IF(NSNP.EQ.0) GO TO 313

READ(1,*)(N1KSNP(I),N2KSNP(I),XSNP(I),I=1,KHSNP)

READ(1,*)(NWSNPNP(I),I=1,KWSNPNP)

CLOSE(UNIT=1)

WRITE(5,*)'Ввод исходных данных закончен'

C Блок 2: Образование имени файла выходных данных и открытие его на диске

313 Z2(1:4)=Z1(1:4)

Z2(5:8)='.REZ'

OPEN(UNIT=3,FILE=Z2,STATUS='NEW')

WRITE(3,320)NRKRTKZH,NSNP,KH,KWKZ,KHSNP,KWSNPNP,SB,

*(N1K(I),N2K(I),X(I),SNG(I),I=1,KH)

320 FORMAT(10X,'В Ы Х О Д Н О Й   Ф А Й Л   П Р О Г Р А М М Ы   TKZ'/

*/' ИСХОДНЫЕ ДАННЫЕ:'//

*'Признак расчета коэффициентов распределения токов КЗ'/

*'   по ветвях схемы NRKRTKZH=',I1/

*'Признак схемы нулевой последовательности NSNP=',I1/

*'Количество ветвей в схеме прямой последовательности KH=',I2/

*'Количество узлов краткого замыкания KWKZ=',I2/