C TP1KA,TO1KA,TN1KA,TP11KA,TO11KA,TN11KA - то же в (KA)
C T1 - суммарный Ток Однофазного КЗ для текущего узла КЗ (отн.ед.)
C T11 - суммарный Ток Двухфазного КЗ для текущего узла КЗ (отн.ед.)
C T1KA,T11KA - то же в (КА)
XRESN=1./X1(1)
TP1=ERES/(XRESP+XRESO+XRESN)
TO1=TP1
TN1=TP1
TP11=ERES/(XRESP+XRESO*XRESN/(XRESO+XRESN))
TO11=TP11*XRESN/(XRESO+XRESN)
TN11=TP11*XRESO/(XRESO+XRESN)
T1=3.*TP1
T1KA=T1*TBAZ
T11=TP11*SQRT(3.*(1.-XRESO*XRESN/((XRESO+XRESN)**2)))
T11KA=T11*TBAZ
TP1KA=TP1*TBAZ
TO1KA=TO1*TBAZ
TN1KA=TN1*TBAZ
TP11KA=TP11*TBAZ
TO11KA=TO11*TBAZ
TN11KA=TN11*TBAZ
C Блок 29: Надо ли выполнять расчет коэффициентов распределения
C токов КЗ по ветвям схемы?
IF(NRKRTKZH.EQ.0) GO TO 89
C Блок 30: Расчет напряжений собственного аварийного режима в СНП
C для текущего узла КЗ
DO 72 I=1,KWSNP
72 UUAV(I)=0.
UUAV(KWSNP)=TN1*XRESN
M=KWSNP
NK=KYIJ
74 M=M-1
M1=M+1
NN=NK-KWZUU(M)
N=NN+1
DO 75 K=M1,KWSNP
DO 73 L=N,NK
IF(NWSNPU(K).EQ.NUSUU(L))UUAV(M)=UUAV(M)+
*UUAV(K)*YIJU(L)/YSUMU(M)
73 CONTINUE
75 CONTINUE
NK=NN
IF(M.GT.1)GO TO 74
C Блок 31: Расчет коэффициентов распределения токов нулевой
C последовательности по ветвям схемы
DO 79 I=1,KHSNP
NL=N1KSNP1(I)
NP=N2KSNP1(I)
IF(NL.EQ.0)GO TO 76
IF(NP.EQ.0)GO TO 166
GO TO 77
76 ULN=0.
77 DO 78 K=1,KWSNP
IF(NP.NE.NWSNPU(K))GO TO 78
UPN=UUAV(K)
78 CONTINUE
IF(NL.NE.0)GO TO 167
IF(N1KSNP(I).LT.NP) GO TO 169
UPN=-UPN
GO TO 167
166 UPN=0.
167 DO 168 K=1,KWSNP
IF(NL.NE.NWSNPU(K))GO TO 168
ULN=UUAV(K)
168 CONTINUE
IF(NP.NE.0)GO TO 169
IF(N2KSNP(I).LT.NL) GO TO 169
ULN=-ULN
169 IF(NL.GT.NP) GO TO 170
AKN(I)=(UPN-ULN)*XSNP(I)/TP1
GO TO 79
170 AKN(I)=(ULN-UPN)*XSNP(I)/TP1
79 CONTINUE
C Блок 32: Расчет полных токов КЗ в ветвях схемы для одноC фазного и двухфазного на землю КЗ для текущего узла КЗ
DO 81 I=1,KH
TW1(I)=0.
TW11(I)=0.
81 AKNP(I)=0.
DO 88 I=1,KH
DO 87 K=1,KHSNP
IF(N1K(I).EQ.N1KSNP(K)) GO TO 84
IF(N1K(I).NE.N2KSNP(K)) GO TO 87
IF(N2K(I).EQ.N1KSNP(K)) GO TO 85
GO TO 87
84 IF(N2K(I).NE.N2KSNP(K)) GO TO 87
85 AKNP(I)=AKN(K)
W1=ABS(AKP(I))
W2=ABS(AKO(I))
W3=ABS(AKNP(I))
IF(W1.EQ.0.)GO TO 86
ZNA=AKP(I)/W1
GO TO 90
86 ZNA=1.
90 TW1(I)=ZNA*(W1+W2+W3)/3.
TW11(I)=ZNA*SQRT(((W1*TP11+W2*TO11)*0.866)**2+
*(ABS(W1*TP11-W2*TO11)*0.5+W3*TN11)**2)/T11
87 CONTINUE
88 CONTINUE
IF(NRKRTKZH.NE.1) GO TO 89
M=0
DO 178 I=1,KH
IF(N1K(I).EQ.NWK) GO TO 179
IF(N2K(I).NE.NWK) GO TO 178
179 M=M+1
AKNP(M)=AKNP(I)
TW1(M)=TW1(I)
TW11(M)=TW11(I)
178 CONTINUE
C Блок 33: Записывание результатов расчетов для текущего узла КЗ в выC ходной файл AAAA.REZ, где AAAA - имя файла входных данных
89 WRITE(3,182)NWK
182 FORMAT(/15X,'НОМЕР УЗЛА КЗ:',I5)
WRITE(3,*)' Базисное напряжение узла КЗ, кв:',UBAZ
WRITE(3,*)' Базисный ток в узле КЗ, ка:',TBAZ
WRITE(3,*)' Результативное сопрот. схемы относительно узла КЗ:'
WRITE(3,*)' -для токов прямой (обр.) последов. (отн.ед.) '
*,XRESP
IF(NSNP.EQ.0.OR.NSOWP.EQ.0) GO TO 190
WRITE(3,*)' - для токов нулевой послед. (отн.ед.) ',XRESN
190 WRITE(3,*)' Периодическая складыв. сверхпереходн. тока КЗ (КА):'
WRITE(3,*)' - трехфазное',T3KA
WRITE(3,*)' - двухфазное',T2KA
IF(NSNP.EQ.0.OR.NSOWP.EQ.0) GO TO 191
WRITE(3,*)' - однофазное',T1KA
WRITE(3,*)' - двухфазное на землю',T11KA
191 WRITE(3,192)T3KA,TO2KA,TO2KA
192 FORMAT(' Симметрические составляющие тока КЗ (КА):'/29X,
*'Прямая Обратная Нулевая'/11X,'- трехфазное',
*F12.4/11X,'- двухфазное',F12.4,F12.4)
IF(NSNP.EQ.0.OR.NSOWP.EQ.0) GO TO 194
WRITE(3,193)TP1KA,TO1KA,TN1KA,TP11KA,TO11KA,TN11KA
193 FORMAT(11X,'- однофазное',F12.4,F12.4,F12.4/3X,
*'- двухфазн. на землю',F12.4,F12.4,F12.4)
194 IF(NRKRTKZH.EQ.0) GO TO 109
WRITE(3,201)
201 FORMAT(' Коэффициенты распределения симметрических суммар',
*'ных токов КЗ по'/' ветвям схемы (положительным направлениям в',
*' ветвях считается направление'/' от конца ветви с большим',
*' номером к концу с меньшим номером)'/
*' Ветвь схемы Прямая Обратная Нулевая')
DO 107 I=1,KHSOKR
IF(NSNP.NE.0.AND.NSOWP.NE.0) GO TO 105
WRITE(3,104)N1K1(I),N2K1(I),AKP(I),AKO(I)
104 FORMAT(I12,I5,(2F12.4))
GO TO 107
105 WRITE(3,106)N1K1(I),N2K1(I),AKP(I),AKO(I),AKNP(I)
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.