Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Диссертация Акимжанов

.pdf
Скачиваний:
69
Добавлен:
13.03.2016
Размер:
2.41 Mб
Скачать

141

Начало

Au ( ) u ( ) x

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление собственных

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

значений матрицыAu ( ) :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 ( ), 2 ( ), 3 ( ),

 

, n ( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формирование матрицы для вычисления

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определителя Вандермонда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

( )

2

( )

 

 

 

n 1

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

( )

2

( )

 

 

 

n 1

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F ( )

1

( )

2 ( )

 

 

 

n 1 ( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

( )

n

2

( )

 

 

 

n 1

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление определителя Вандермонда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( ) det F ( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формирование матриц для вычисления определителей i

( ):

 

 

 

 

 

e 1 ( )

( )

 

2

( )

 

n 1

( )

 

 

 

 

1

 

e 1 ( )

 

2

( )

n 1

( )

 

 

e 2 ( )

1

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

( )

 

2

( )

 

n 1

( )

 

 

 

 

1

 

 

e 2 ( )

 

2

( )

n 1

( )

 

 

 

 

 

2

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

F ( )

e

3 ( )

( )

 

2

( )

 

 

n 1

( )

;

F ( )

 

 

1

 

 

e

3 ( )

 

2

( )

 

n 1

( )

 

;

0

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

3

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

( )

 

 

 

 

 

 

e

n ( )

2

 

n 1

( )

 

 

 

 

1

 

 

 

( )

 

 

 

 

 

 

n

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

n

 

 

 

 

1

( )

e 1 ( )

n 1

( )

 

1

( )

 

2

( )

e 1 ( )

 

 

 

1

 

 

 

1

( )

 

 

 

1

1

 

 

 

 

1

( )

e 2 ( )

n 1

 

1

( )

 

2

( )

e 2 ( )

 

 

 

2

 

 

 

2

 

 

 

 

2

2

 

 

 

F ( )

1

( )

e 3 ( )

n 1

( ) ;

F

( )

1

( )

 

2

( )

e 3 ( )

2

 

 

3

 

 

 

3

 

(n 1)

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

e

 

 

 

 

 

 

 

( )

 

 

( )

e

 

 

1

 

( )

 

1

2

n ( )

 

 

 

 

n (

 

)

n 1

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

n

n

 

 

 

.

 

 

 

 

 

Вычисление определителей i ( ) :

 

 

 

 

 

 

 

 

 

 

 

 

0 ( ) det F0 ( ) ;

 

1 ( ) det F1 ( ) ;

 

 

 

 

 

 

 

 

 

 

 

2 ( ) det F2 ( ) ;

n 1

( ) det Fn 1 ( ) ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление функции от матрицы:

 

 

 

 

 

 

 

 

 

e u ( ) x

 

( ) A

0 ( ) ( ) A 1 ( )

2

( ) A 2 ( )

 

n 1

( ) A n 1

( )

 

 

0

u

1

u

 

 

u

 

 

 

u

 

 

 

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Конец

Рисунок П.1.6. Алгоритм расчета функций от матрицы погонных параметров

многопроводной ВЛ e u ( n )l , e u ( n )l

142

Начало

Ai ( ) i ( ) x

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление собственных

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

значений матрицы Ai ( ) :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 ( ), 2 ( ), 3 ( ),

 

, n ( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формирование матрицы для вычисления

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определителя Вандермонда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

( )

2

( )

 

 

 

n 1

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

( )

2

( )

 

 

 

n 1

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F ( )

1

( )

2 ( )

 

 

 

n 1 ( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

( )

n

2

( )

 

 

 

n 1

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление определителя Вандермонда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( ) det F ( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формирование матриц для вычисления определителей i

( ):

 

 

 

 

 

e 1 ( )

( )

 

2

( )

 

n 1

( )

 

 

 

 

1

 

e 1 ( )

 

2

( )

n 1

( )

 

 

e 2 ( )

1

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

( )

 

2

( )

 

n 1

( )

 

 

 

 

1

 

 

e 2 ( )

 

2

( )

n 1

( )

 

 

 

 

 

2

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

F ( )

e

3 ( )

( )

 

2

( )

 

 

n 1

( )

;

F ( )

 

 

1

 

 

e

3 ( )

 

2

( )

 

n 1

( )

 

;

0

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

3

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

( )

 

 

 

 

 

 

e

n ( )

2

 

n 1

( )

 

 

 

 

1

 

 

 

( )

 

 

 

 

 

 

n

n

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

n

 

 

 

 

1

( )

e 1 ( )

n 1

( )

 

1

( )

2

 

 

 

1

 

1

 

 

 

 

 

1

1

 

1

( )

e 2 ( )

n 1

( )

 

1

( )

2

 

 

 

2

 

2

 

 

 

 

 

2

2

F ( )

1

( )

e 3 ( )

n 1

( ) ;

F

( )

1

( )

2

2

 

 

3

 

3

 

 

(n 1)

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

( )

e n ( )

n 1

( )

 

 

 

1

( )

2

 

 

 

n

 

n

 

 

 

 

 

n

n

( ) ( ) ( )

( )

e 1 ( ) e 2 ( ) e 3 ( )

e n ( )

.

 

 

 

 

 

 

Вычисление определителей i ( ) :

 

 

 

 

 

 

 

 

 

 

 

 

0 ( ) det F0 ( ) ;

 

1 ( ) det F1 ( ) ;

 

 

 

 

 

 

 

 

 

 

 

 

2 ( ) det F2 ( ) ;

n 1 ( ) det Fn 1 ( ) ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисление функции от матрицы:

 

 

 

 

 

 

 

 

 

e i ( ) x

 

( ) A

0 ( ) ( ) A1 ( )

( ) A 2 ( )

 

n 1

( ) A n 1

( )

 

 

0

i

1

i

2

i

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Конец

Рисунок П.1.7. Алгоритм расчета функций от матрицы погонных параметров

многопроводной ВЛ e i ( n )l , e i ( n )l

 

 

143

 

 

Начало

 

 

Ввод

 

 

Uн ( ),

Iн ( ),

 

 

L,

NL

 

 

x L / NL

 

 

i=1,NL

 

 

i=1

 

 

 

Ui ( ) Uн ( ),

 

 

 

Ii ( ) Iн ( ).

 

 

Ui ( ) U(i 1) LNL ( ),

 

 

Ii ( ) I(i 1) LNL ( ).

 

 

Расчет коэффициентов

 

 

A( )i , B( )i , C( )i , D( )i

 

 

x=1,LNL

 

 

Расчет

 

 

e u ( ) x , e u ( ) x , e i ( ) x , e i ( ) x

Расчет напряжений и токов вдоль

 

 

трассы ВЛ с шагом x

U

x,i

( ) e u ( n ) x A ( ) e u ( n ) x B ( );

 

i

i

I

x,i

( ) e i ( n ) xC ( ) e i ( n ) x D ( ).

 

i

i

Расчет потерь активной мощности

 

 

на одном из режимов

 

 

m

l

 

 

P( k ) ri ( k ) Ii2,( k ) (x)dx,

 

 

i 1

0

 

 

Конец

Рисунок П.1.8. Базовый алгоритм расчета потерь мощности в многопроводной ВЛ

с учетом распределенности параметров

144

Начало

Расчет Z ( n ), Y ( n )

Ввод L,

U А ( n ) 1.0 j0.0; U В ( n ) 0,866 j0, 5; UC ( n ) 0,866 j0, 5.

Расчет коэффициентов

An , Bn , Cn , Dn

Расчет напряжений и токов вдоль трассы ВЛ с шагом x

U x ( n ) e u ( n ) x An e u ( n ) x Bn ;

Ix ( n ) e i ( n ) xCn e i ( n ) x Dn .

Расчет проводимостей фаз и тросов

Yx ( n ) Ix ( n ) U x ( n )

Конец

Рисунок П.1.9. Алгоритм расчета частотных характеристик многопроводной ВЛ с

учетом распределенности параметров

145

Приложение 2. Частотные характеристики ВЛ 110 кВ различной

длины, выполненных проводом АС-120 на опорах марки П110-5В в холостом

режиме

 

 

 

 

См

0.2

 

 

20 км

 

 

 

 

 

 

фаз,

0.15

 

 

 

 

 

 

 

 

 

проводимости

0.1

 

 

 

 

0.05

 

 

 

 

 

 

 

 

 

Модуль

0

 

32

47

 

 

 

 

2 7

12 17 22 27

Номер37 42 гармоники52 57 62 67 72 77 82 87 92 97

 

 

Фаза А

 

Фаза В

Фаза С

См

0.12

 

 

40 км

 

0.10

 

 

 

 

фаз,

 

 

 

 

0.08

 

 

 

 

проводимости

 

 

 

 

0.06

 

 

 

 

0.04

 

 

 

 

0.02

 

 

 

 

 

 

 

 

 

Модуль

0.00

 

32

47

 

 

 

 

2 7

12 17 22 27

Номер37 42 гармоники52 57 62 67 72 77 82 87 92 97

 

 

Фаза А

 

Фаза В

Фаза С

См

0.09

 

 

60 км

 

0.08

 

 

 

 

фаз,

0.07

 

 

 

 

0.06

 

 

 

 

проводимости

 

 

 

 

0.05

 

 

 

 

0.04

 

 

 

 

0.03

 

 

 

 

0.02

 

 

 

 

0.01

 

 

 

 

Модуль

0

 

32

47

 

 

 

 

2 7

12 17 22 27

Номер37 42 гармоники52 57 62 67 72 77 82 87 92 97

 

 

Фаза А

 

Фаза В

Фаза С

См

0.07

 

 

80 км

 

0.06

 

 

 

 

фаз,

 

 

 

 

0.05

 

 

 

 

проводимости

0.04

 

 

 

 

0.03

 

 

 

 

0.02

 

 

 

 

0.01

 

 

 

 

 

 

 

 

 

Модуль

0

 

32

47

 

 

 

 

2 7

12 17 22 27

Номер37 42 гармоники52 57 62 67 72 77 82 87 92 97

 

 

Фаза А

 

Фаза В

Фаза С

См

0.16

 

 

30 км

 

0.14

 

 

 

 

фаз,

0.12

 

 

 

 

 

 

 

 

 

проводимости

0.1

 

 

 

 

0.08

 

 

 

 

0.06

 

 

 

 

0.04

 

 

 

 

0.02

 

 

 

 

Модуль

0

12 17 22 27

32

47

 

 

 

2 7

Номер37 42 гармоники52 57 62 67 72 77 82 87 92 97

 

 

Фаза А

 

Фаза В

Фаза С

См

0.1

 

 

50 км

 

 

 

 

 

 

фаз,

0.08

 

 

 

 

 

 

 

 

 

проводимости

0.06

 

 

 

 

0.04

 

 

 

 

0.02

 

 

 

 

 

 

 

 

 

Модуль

0

12 17 22 27

32

47

 

 

 

2 7

Номер37 42 гармоники52 57 62 67 72 77 82 87 92 97

 

 

Фаза А

 

Фаза В

Фаза С

См

0.08

 

 

70 км

 

0.07

 

 

 

 

фаз,

0.06

 

 

 

 

 

 

 

 

 

проводимости

0.05

 

 

 

 

0.04

 

 

 

 

0.03

 

 

 

 

0.02

 

 

 

 

0.01

 

 

 

 

Модуль

0

12 17 22 27

32

47

 

 

 

2 7

Номер37 42 гармоники52 57 62 67 72 77 82 87 92 97

 

 

Фаза А

 

Фаза В

Фаза С

См

0.06

 

 

90 км

 

0.05

 

 

 

 

фаз,

 

 

 

 

0.04

 

 

 

 

проводимости

 

 

 

 

0.03

 

 

 

 

0.02

 

 

 

 

0.01

 

 

 

 

 

 

 

 

 

Модуль

0

12 17 22 27

32

47

 

 

 

2 7

Номер37 42 гармоники52 57 62 67 72 77 82 87 92 97

 

 

Фаза А

 

Фаза В

Фаза С

146

Приложение 3. Текст программы расчета основных и добавочных потерь в одно- и двухцепных ВЛ 110 кВ

USE MSIMSL

DIMENSION AAH(15),BBH(13)

COMPLEX UK1(8)*16,AIK1(8)*16,UK10*16,UK11*16,UK12*16,AIK10*16, /AIK11*16,AIK12*16,AL*16

DOUBLE PRECISION AIM(700,50,3),FIM(700,50,3),UM(700,50,3), /FUM(700,50,3),AIM1(700,50,3),AIM2(700,50,3),UM1(700,50,3), /UM2(700,50,3),RMSU(700,3),RMSI(700,3),FUNU(700,3),FUNI(700,3), /FU(700,3),FI(700,3),PPR1(700),PPR2(700),PD(3,700),PPP(1000,50), /PPP1(1000,50),PPP2(1000,50),PPP3(1000,50),PPP4(1000,50), /PPP5(1000,50),PPP6(1000,50),PPP7(1000,50),PPP8(1000,50),WD(2,50) INTEGER II(10)

DOUBLE PRECISION PI,FF,SS2,SS0,PP1,PP2,RPR,PRP,WD0 REAL KNSU(700,3),KNSI(700,3)

COMMON MM,M,M1,MT,M10,M20,PR,K1,K2,K3,N1,N2,N3,MPR,MTR,MMT OPEN(UNIT=1,FILE='1.ДЛИНА И ПЕРИОД')

C************************************************************** C ВЫБОР КОЛИЧЕСТВА ПРОВОДОВ

COPEN(UNIT=2,FILE='2.3.ВВОДИМЫЕ ОБЩИЕ ДАННЫЕ') OPEN(UNIT=2,FILE='2.4.ВВОДИМЫЕ ОБЩИЕ ДАННЫЕ')

C OPEN(UNIT=2,FILE='2.6.ВВОДИМЫЕ ОБЩИЕ ДАННЫЕ')

C OPEN(UNIT=2,FILE='2.7.ВВОДИМЫЕ ОБЩИЕ ДАННЫЕ')

C OPEN(UNIT=2,FILE='2.8.ВВОДИМЫЕ ОБЩИЕ ДАННЫЕ')

C************************************************************** OPEN(UNIT=4,FILE='1.ПОДСТАНЦИЯ И ЕЕ ПРИСОЕДИНЕНИЯ') OPEN(UNIT=5,FILE='ВВЕДЕННЫЕ ОБЩИЕ ДАННЫЕ') OPEN(UNIT=6,FILE='ВВЕДЕННЫЕ ПАРАМЕТРЫ РЕЖИМА')

OPEN(UNIT=8,FILE='РЕЗУЛЬТАТЫ РАСЧЕТА') OPEN(UNIT=10,FILE='EPURAU') OPEN(UNIT=11,FILE='EPURAI') OPEN(UNIT=12,FILE='R11') OPEN(UNIT=30,FILE='PPP') OPEN(UNIT=31,FILE='PPP1') OPEN(UNIT=32,FILE='PPP2') OPEN(UNIT=33,FILE='PPP3') OPEN(UNIT=34,FILE='PPP4') OPEN(UNIT=35,FILE='PPP5') OPEN(UNIT=36,FILE='PPP6') OPEN(UNIT=37,FILE='PPP7') OPEN(UNIT=38,FILE='PPP8')

OPEN(unit=101,file='3.1.ФАЗА А')

OPEN(unit=111,file='3.1.ФАЗА А (продолжение)') OPEN(unit=102,file='3.2.ФАЗА В')

OPEN(unit=122,file='3.2.ФАЗА В (продолжение)') OPEN(unit=103,file='3.3.ФАЗА С')

OPEN(unit=133,file='3.3.ФАЗА С (продолжение)') READ(4,333) NN1

READ(4,335) READ(4,334) AAH READ(4,335)

147

DO 470 JJ=1,NN1

READ(4,335) BBH READ(4,333) N1 READ(4,335) II(JJ)=N1

N20=0

333FORMAT(5I3)

334FORMAT(15A1)

335FORMAT(13A1)

336FORMAT(204E9.2)

337FORMAT(204(E8.2,1X))

N2=1

N3=N1

K1=50

K2=1

K3=50

pi=3.14159265358979

AL=CMPLX(-0.5, 0.866025) FF=50.

1 FORMAT(130(E9.3,1X))

31FORMAT(130(E10.3))

WRITE(6,102)

102FORMAT('ГАРМОНИКИ НАПРЯЖЕНИЯ И ИХ ФАЗЫ ПО ФАЗЕ 1')

DO 111 I=1,N1

READ(101,337) (UM(I,K,1),FUM(I,K,1),K=2,50)

WRITE(6,337)(UM(I,K,1),FUM(I,K,1),K=2,50) 111 CONTINUE

WRITE(6,103)

103 FORMAT('ГАРМОНИКИ ТОКА И ИХ ФАЗЫ ПО ФАЗЕ 1,KНС НАПРЯЖЕНИЯ,ТОКА,ДЕЙ

/СТВУЮЩИЕ ЗНАЧЕНИЯ НАПРЯЖЕНИЯ, ТОКА,ОСНОВНЫЕ ГАРМОНИКИ НАПРЯЖЕНИЯ,

/ТОКА ПО ФАЗЕ 1')

DO 222 I=1,N1

READ(111,337) (AIM(I,K,1),FIM(I,K,1),K=2,50), /knsu(I,1),knsi(I,1),rmsu(I,1),rmsi(I,1),funu(I,1),funi(I,1), /fu(I,1),fI(I,1) WRITE(6,337)(AIM(I,K,1),FIM(I,K,1),K=2,50), /knsu(I,1),knsi(I,1),rmsu(I,1),rmsi(I,1),funu(I,1),funi(I,1), /fu(I,1),fI(I,1)

222CONTINUE WRITE(6,105)

105FORMAT('ГАРМОНИКИ НАПРЯЖЕНИЯ И ИХ ФАЗЫ ПО ФАЗЕ 2')

DO 4 I=1,N1

READ(102,337) (UM(I,K,2),FUM(I,K,2),K=2,50)

WRITE(6,337)(UM(I,K,2),FUM(I,K,2),K=2,50) 4 CONTINUE

WRITE(6,106)

106 FORMAT('ГАРМОНИКИ ТОКА И ИХ ФАЗЫ ПО ФАЗЕ 2,КНС НАПРЯЖЕНИЯ,ТОКА,ДЕЙ

/СТВУЮЩИЕ ЗНАЧЕНИЯ НАПРЯЖЕНИЯ, ТОКА,ОСНОВНЫЕ ГАРМОНИКИ НАПРЯЖЕНИЯ,

148

/ТОКА ПО ФАЗЕ 2')

DO 5 I=1,N1

READ(122,337) (AIM(I,K,2),FIM(I,K,2),K=2,50), /knsu(I,2),knsi(I,2),rmsu(I,2),rmsi(I,2),funu(I,2),funi(I,2), /fu(I,2),fI(I,2) WRITE(6,337)(AIM(I,K,2),FIM(I,K,2),K=2,50), /knsu(I,2),knsi(I,2),rmsu(I,2),rmsi(I,2),funu(I,2),funi(I,2), /fu(I,2),fI(I,2)

5CONTINUE WRITE(6,108)

108FORMAT('ГАРМОНИКИ НАПРЯЖЕНИЯ И ИХ ФАЗЫ ПО ФАЗЕ 3')

DO 777 I=1,N1

READ(103,337) (UM(I,K,3),FUM(I,K,3),K=2,50)

WRITE(6,337)(UM(I,K,3),FUM(I,K,3),K=2,50) 777 CONTINUE

WRITE(6,107)

107 FORMAT('ГАРМОНИКИ ТОКА И ИХ ФАЗЫ ПО ФАЗЕ 3,КНС НАПРЯЖЕНИЯ,ТОКА,ДЕЙ

/СТВУЮЩИЕ ЗНАЧЕНИЯ НАПРЯЖЕНИЯ, ТОКА, ОСНОВНЫЕ ГАРМОНИКИ НАПРЯЖЕНИЯ,

/ТОКА ПО ФАЗЕ 3')

DO 888 I=1,N1

READ(133,337) (AIM(I,K,3),FIM(I,K,3),K=2,50), /knsu(I,3),knsi(I,3),rmsu(I,3),rmsi(I,3),funu(I,3),funi(I,3), /fu(I,3),fI(I,3) WRITE(6,337)(AIM(I,K,3),FIM(I,K,3),K=2,50), /knsu(I,3),knsi(I,3),rmsu(I,3),rmsi(I,3),funu(I,3),funi(I,3), /fu(I,3),fI(I,3)

888 CONTINUE WRITE(6,11)

11FORMAT(80x) DO 13 I=1,N1 DO 13 L=1,3

UM(I,1,L)=FUNU(I,L)*SQRT(2.)

FUM(I,1,L)=FU(I,L)*PI/180.

UM1(I,1,L)=UM(I,1,L)*COS(FUM(I,1,L))

UM2(I,1,L)=UM(I,1,L)*SIN(FUM(I,1,L))

13CONTINUE DO 14 I=1,N1 DO 14 L=1,3 DO 14 K=2,50

UM(I,K,L)=UM(I,K,L)*UM(I,1,L)/100.

FUM(I,K,L)=FUM(I,K,L)*PI/180.

UM1(I,K,L)=UM(I,K,L)*DCOS(FUM(I,K,L))

UM2(I,K,L)=UM(I,K,L)*DSIN(FUM(I,K,L))

14CONTINUE DO 16 I=1,N1 DO 16 L=1,3

AIM(I,1,L)=FUNI(I,L)*SQRT(2.)

FIM(I,1,L)=FI(I,L)*PI/180.

AIM1(I,1,L)=AIM(I,1,L)*COS(FIM(I,1,L))

AIM2(I,1,L)=AIM(I,1,L)*SIN(FIM(I,1,L))

149

16CONTINUE DO 117 I=1,N1 DO 117 L=1,3 DO 117 K=2,50

AIM(I,K,L)=AIM(I,K,L)*AIM(I,1,L)/100.

FIM(I,K,L)=FIM(I,K,L)*PI/180.

AIM1(I,K,L)=AIM(I,K,L)*COS(FIM(I,K,L))

AIM2(I,K,L)=AIM(I,K,L)*SIN(FIM(I,K,L))

117CONTINUE

READ(1,70) MM,MPR,MTR,DT,MT M=MPR+MTR

MMT=MM/MT

WRITE(5,70) MM,MPR,MTR,DT,MT,M,MMT IF(M.LE.6)M1=M

IF(M.GT.6)M1=6

M10=2*M

M20=4*M

DO 1500 N=N2,N3 DO 1500 K=K2,K3

PR=0.

1700 CONTINUE PR=PR+1

UK1(1)=CMPLX(UM1(N,K,1),UM2(N,K,1))

UK1(2)=CMPLX(UM1(N,K,2),UM2(N,K,2))

UK1(3)=CMPLX(UM1(N,K,3),UM2(N,K,3))

UK1(4)=CMPLX(0.0,0.0)

AIK1(1)=CMPLX(AIM1(N,K,1),AIM2(N,K,1))

AIK1(2)=CMPLX(AIM1(N,K,2),AIM2(N,K,2))

AIK1(3)=CMPLX(AIM1(N,K,3),AIM2(N,K,3))

AIK1(4)=CMPLX(0.0,0.0) IF(K.GT.1) GOTO 1111

IF(K.EQ.1.AND.PR.EQ.2) GOTO 1111 UK10=(UK1(1)+UK1(2)+UK1(3))/3. UK11=(UK1(1)+UK1(2)*AL+UK1(3)*AL**2)/3. UK12=(UK1(1)+UK1(2)*AL**2+UK1(3)*AL)/3. SKU2=SQRT(REAL(UK12)**2+AIMAG(UK12)**2)/ /SQRT(REAL(UK11)**2+AIMAG(UK11)**2)*100 SKU0=SQRT(REAL(UK10)**2+AIMAG(UK10)**2)/ /SQRT(REAL(UK11)**2+AIMAG(UK11)**2)*100 UK1(1)=UK11

UK1(2)=UK11*AL**2

UK1(3)=UK11*AL

AIK10=(AIK1(1)+AIK1(2)+AIK1(3))/3.

AIK11=(AIK1(1)+AIK1(2)*AL+AIK1(3)*AL**2)/3.

AIK12=(AIK1(1)+AIK1(2)*AL**2+AIK1(3)*AL)/3.

SKI2=SQRT(REAL(AIK12)**2+AIMAG(AIK12)**2)/

/SQRT(REAL(AIK11)**2+AIMAG(AIK11)**2)*100

SKI0=SQRT(REAL(AIK10)**2+AIMAG(AIK10)**2)/

/SQRT(REAL(AIK11)**2+AIMAG(AIK11)**2)*100

AIK1(1)=AIK11

AIK1(2)=AIK11*AL**2

AIK1(3)=AIK11*AL

 

150

1111

CONTINUE

 

CALL RASCHET(UK1,AIK1,K,N,PPP,PP1,PP2,

 

/ppp1,ppp2,ppp3,ppp4,ppp5,ppp6,ppp7,ppp8)

 

IF(K.EQ.1.AND.PR.EQ.1)PPR1(N)=PP1

 

IF(K.EQ.1.AND.PR.EQ.2)PPR2(N)=PP2

 

IF(K.EQ.1.AND.PR.EQ.1.) GOTO 1700

 

IF(PR.EQ.2) GOTO 1500

1500

CONTINUE

 

DO 1501 N=N2,N3

 

PRP=0.

 

DO 1060 K=K2,K3

 

PRP=PRP+PPP(N,K)

1060

CONTINUE

 

RPR=0.

 

DO 1061 K=2,K3

 

RPR=RPR+PPP(N,K)

1061

CONTINUE

 

SS1=(PPP(N,1)/PPR1(N)-1.)*100

 

SS2=(RPR/PPR1(N))*100

 

SS0=PPP(N,1)-PPR1(N)

 

WRITE(30,196) (PPP(N,K),K=K2,K3),PRP,PPR1(N),SS0,RPR,SS1,SS2

 

WRITE(31,196) (PPP1(N,K),K=K2,K3)

 

WRITE(32,196) (PPP2(N,K),K=K2,K3)

 

WRITE(33,196) (PPP3(N,K),K=K2,K3)

 

WRITE(34,196) (PPP4(N,K),K=K2,K3)

 

WRITE(35,196) (PPP5(N,K),K=K2,K3)

 

WRITE(36,196) (PPP6(N,K),K=K2,K3)

 

WRITE(37,196) (PPP7(N,K),K=K2,K3)

 

WRITE(38,196) (PPP8(N,K),K=K2,K3)

1501

CONTINUE

 

WD0=0.

 

DO 1052 K=K2,K3

 

WD(1,K)=0.

 

DO 1052 I=N2,N3

 

WD0=WD0+PPP(I,K)*DT/60000.

 

WD(1,K)=WD(1,K)+PPP(I,K)*DT/60000.

1052

CONTINUE

 

WD1=WD0-WD(1,1)

 

DO 1053 K=K2,K3

 

WD(2,K)=WD(1,K)/WD0*100.

1053

CONTINUE

 

WD4=0.

 

DO 1056 K=14,K3

 

WD4=WD4+WD(2,K)

1056

CONTINUE

 

WD10=0.

 

DO 1057 I=N2,N3

 

WD10=WD10+PPR1(I)*DT/60000.

1057

CONTINUE

 

DO 1054 I=N2,N3

 

PD(1,I)=0.

 

PD(2,I)=0.