C ток трехфазного КЗ) (отн.ед.)
C XRESP – результирующее сопротивление СПП относительно узла КЗ (отн.ед.)
C XRESO - результирующее сопротивление схемы обратной последC овательности (СОП) относительно узла КЗ (отн.ед.)
C TP3 - ток прямой последовательности для трехфазного КЗ (отн.ед.)
C TP2 - ток прямой последовательности для двухфазного КЗ (отн.ед.)
C TO2 - ток обратной последовательности для двухфазного КЗ (отн.ед.)
C UBAZ - Напряжение Базисное в текущем узле КЗ (кв)
C TBAZ - Ток Базисный в текущем узле КЗ (кA)
C T3KA - суммарный ток трехфазного КЗ для текущего узла КЗ (KA)
C T2KA - суммарный ток двухфазного КЗ для текущего узла КЗ (KA)
C TO2KA - ток обратной последовательности двухфазного КЗ (KA)
ERES=E1(1)/X1(1)
XRESP=1./X1(1)
XRESO=XRESP
TP3=E1(1)
TP2=TP3/2
TO2=TP2
UBAZ=UB(KWKZB)
TBAZ=SB/(SQRT(3.)*UBAZ)
T3KA=TP3*TBAZ
T2KA=T3KA*SQRT(3.)/2
TO2KA=TO2*TBAZ
C Блок 17: Расчет узловых напряжений доаварийного режима
UDAV(KW)=ERES
M=KW
NK=KYIJ
55 M=M-1
M1=M+1
NN=NK-KWZUU(M)
N=NN+1
DO 52 K=M1,KW
DO 51 L=N,NK
IF(NWU(K).EQ.NUSUU(L))UDAV(M)=UDAV(M)+UDAV(K)*YIJU(L)
51 CONTINUE
52 CONTINUE
UDAV(M)=(UDAV(M)+AJU(M))/YSUMU(M)
NK=NN
IF(M.GT.1)GO TO 55
C Блок 18: Расчет узловых напряжений собственного аварийного режима
UUAV(KW)=-UDAV(KW)
M=KW
NK=KYИJ
57 M=M-1
M1=M+1
NN=NK-KWZUU(M)
N=NN+1
DO 56 K=M1,KW
DO 58 L=N,NK
IF(NWU(K).EQ.NUSUU(L))UUAV(M)=UUAV(M)+
*UUAV(K)*YIJU(L)/YSUMU(M)
58 CONTINUE
56 CONTINUE
NK=NN
IF(M.GT.1)GO TO 57
C Блок 19: Расчет узловых напряжений аварийного режима
DO 59 I=1,KW
59 UAV(I)=UDAV(I)+UUAV(I)
C Блок 20: Расчет коэффициентов распределения токов прямой
C и обратной последовательности по ветвям схемы
DO 63 I=1,KH
NL=N1K(I)
NP=N2K(I)
IF(NL.EQ.0)GO TO 60
IF(NP.EQ.0)GO TO 152
GO TO 61
60 ULP=0
ULO=0
61 DO 62 K=1,KW
IF(NP.NE.NWU(K)) GO TO 62
UPP=UAV(K)
UPO=UUAV(K)/2
62 CONTINUE
IF(NL.EQ.0)GO TO 155
GO TO 153
152 UPP=0
UPO=0
153 DO 154 K=1,KW
IF(NL.NE.NWU(K))GO TO 154
ULP=UAV(K)
ULO=UUAV(K)/2
154 CONTINUE
155 IF(NL.GT.NP)GO TO 156
AKP(I)=(UPP-ULP-E(I)/X(I))*X(I)/TP3
AKP1(I)=AKP(I)
AKO(I)=(ULO-UPO)*X(I)/TP2
GO TO 63
156 AKP(I)=(ULP-UPP-E(I)/X(I))*X(I)/TP3
AKP1(I)=AKP(I)
AKO(I)=(UPO-ULO)*X(I)/TP2
63 CONTINUE
C Блок 21: Расчет коэффициентов распределения токов КЗ по
C ветвям схемы, которые подключены к узлу КЗ
KHSOKR=KH
DO 176 I=1,KH
N1K1(I)=N1K(I)
176 N2K1(I)=N2K(I)
IF(NRKRTKZH.NE.1)GO TO 177
M=0
DO 174 I=1,KH
IF(N1K(I).EQ.NWK)GO TO 175
IF(N2K(I).NE.NWK)GO TO 174
175 M=M+1
N1K1(M)=N1K(I)
N2K1(M)=N2K(I)
AKP(M)=AKP(I)
AKO(M)=AKO(I)
174 CONTINUE
KHSOKR=M
177 CONTINUE
C Блок 22: Нужно ли рассчитывать токи однофазных и
C двухфазных на землю КЗ для текущего узла КЗ?
IF(NSNP.EQ.0)GO TO 89
C Блок 23: Входит ли текущий узел КЗ в состав СНП?
NSOWP=0
DO 69 I=1,KWSNP
IF(NWK.NE.NWSNP(I))GO TO 69
NSOWP=1
69 CONTINUE
IF(NSOWP.EQ.0)GO TO 89
C Блок 24: Пересылка информации о СНП в рабочие массивы
DO 70 I=1,KHSNP
N1KSNP2(I)=N1KSNP1(I)
N2KSNP2(I)=N2KSNP1(I)
X1(I)=XSNP(I)
70 E1(I)=0.
C Блок 25: Перестановка номеров узлов в массиве NWSNP так,
C чтобы узел КЗ оказался на последнем месте
KWSNPU=0
DO 71 I=1,KWSNP
IF(NWSNP(I).EQ.NWK) GO TO 71
KWSNPU=KWSNPU+1
NWSNPU(KWSNPU)=NWSNP(I)
71 CONTINUE
NWSNPU(KWSNP)=NWK
C Блок 26: Подготовка КНВ
KHB=KHSNP
C Блок 27: Подпрограмма ZHORTKA(Свертка); свертка СНП
C выполняется по той же подпрограмме, что и СПП
CALL ZHORTKA(N1KSNP2,N2KSNP2,X1,E1,NWSNPU,
*KWSNPU,KHB,KYIJ,KWZUU,YSUMU,AJU,NUSUU,YIJU)
C Блок 28: Расчет результирующих параметров однофазного КЗ
C для текущего узла краткого замыкания NWK:
C XRESN - результирующее сопрот. СНП относительно узла КЗ (отн.ед.)
C TP1 - Ток Прямой последовательности для Однофазного КЗ (отн.ед.)
C TO1 - Ток Обратной последовательности для Однофазного КЗ (отн.ед.)
C TN1 - Ток Нулевой последовательности для Однофазного КЗ (отн.ед.)
C TP11 - Ток Прямой последовательности для Двухфазного КЗ (отн.ед.)
C TO11 - Ток Обратной последовательности для Двухфазного КЗ (отн.ед.)
C TN11 - Ток Нулевой последовательности для Двухфазного КЗ (отн.ед.)
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.