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

книги из ГПНТБ / Данилевич, Я. Б. Добавочные потери в турбо- и гидрогенераторах

.pdf
Скачиваний:
7
Добавлен:
20.10.2023
Размер:
8.39 Mб
Скачать

линеаризируется введением понятия эквивалентной магнитной проницаемости.

При рассмотрении в настоящей работе добавочных потерь в крайних пакетах, вызванных торцовым магнитным потоком, при­ нимаются следующие допущения.

1.Пренебрегается кривизной сердечника статора, так как процесс затухания происходит в направлении плоской поверх­ ности торца. Кроме того, как показали расчеты, учет цилиндричности сердечника статора не приводит к сколько-нибудь сущест­ венному уточнению величин электромагнитного поля в зоне край­ него пакета.

2.Относительная магнитная проницаемость крайних пакетов имеет различную величину для разных направлений и является тензором, матрица которого имеет вид

/Н-л*

0

0 \

V 0

0

у.як/

где к — номер области. Относительная магнитная проницаемость по оси у в пределах длины одного пакета принимается постоянной и равной магнитной проницаемости по оси х. При этом выбор ее величины производится с учетом распределения основного потока.

3. Влияние ярма сердечника статора в пределах полюсного деления учитывается условным увеличением высоты зубца, кото­ рая, как показано в [44], должна быть принята равной l—3h.

Для определения электромагнитного поля в крайнем пакете, считая изменение всех величин поля во времени синусоидальным (множитель е ^ ' в формулах в дальнейшем опускается, как общий для всех выражений), необходимо решить уравнения Максвелла

rotH,.

"

Р '

 

 

(2. 22}

rot È,. :

—]шВк,

для векторов Ек и Нк электрического и магнитного полей в трех свя­ занных областях: в области крайнего пакета (1), в области воздуш­ ного промежутка между крайним и вторым пакетами (2) и в области сердечника статора (3).

Для расчета поля в крайнем пакете могут быть приняты следующие граничные условия.

1. На поверхности крайнего пакета задана индукция аксиаль­ ного поля

B2 = Re в АО, »)» ~

60

которая может быть рассчитана в соответствии с указаниями раз­

дела 2 . 2 . Так как расчетная высота зубца принята равной

l=3h,

начало координат располагается выше дна паза (рис. 2 - 12) .

Поэ­

тому, учитывая, что в зоне ярма аксиальная составляющая индук­

ции экранируется нажимной плитой, Вг(0, у)

можно

представить

в интервале —1/2 < у <

1/2 в виде ряда Фурье с периодом Т=21

с

четным продолжением

и равным

нулю

средним

значением:

 

 

 

я=1,2з, .. .

 

711UJ

 

 

 

 

 

5 „ c o s — ,

 

(2.23)

 

 

J/2

 

 

 

 

 

где Вп = у

\ В. (О, у) cos —j~ dy — коэффициенты

разложения

 

 

о

 

 

 

 

 

•ЯДО, у) в ряд Фурье.

 

 

 

 

 

 

Расчеты

показывают,

что вторая

 

половина кривой индукции

в

интервале

от 0 до 112 (рис. 2 - 12) не оказывает существенного

влияния

(с точностью до 5 % ) на величину потерь в первой поло-

вине условной высоты

зуоца

в интервале от U до — у .

2 .

При z —> со аксиальное

поле В. > 0.

3 .

Так как размеры

зубца

крайнего пакета ограничены, то при

X — +

у

составляющая

плотности тока Sx ~ 0, а при у = + у

составляющая S = 0.

Расчеты растекания токов в тонких прямоугольных пластинках показывают, что при существующих отношениях ширины зубца к его высоте величина тока в области (например, ОАВС, рис. 2 - 12) , прилегающей к краю пластинки, изменяется не более чем на 6% при удалении одной стороны зубца на бесконечность. Следова­ тельно, не внося существенной погрешности в физическую и ко­ личественную стороны явления, можно упростить задачу, выпол­ нив граничное условие Sx = 0 только на стороне х — у , считая дру­ гую сторону удаленной на бесконечность. При'этом расчет следует

производить для Ѵ4 области условного зубца (ОАВС).

4. На поверхностях раздела трех областей выполняется непре­ рывность тангенциальных составляющих Ё ж ff. При этом равенство нормальных составляющих индукции выполняется автоматически.

Введя понятие векторного потенциала в систему (2 . 22) и при­ нимая во внимание, что из-за расслоенности стали составляющая векторного потенциала по оси z равна нулю, получаем следующую •систему дифференциальных уравнений для составляющей вектор­ ного потенциала по оси х (см. приложение 2 ) :

61

d*-Àxl

*дАя1

. 1

д»Ал

мы 1

дх*

diß

oj

o>22

p

(52^4 ^

а 2 І т

1

d*Àx3

/l ü ftr3

• " Я

I " " Я , * " "3

P

Öx2

' <ty2

"T" a 3

<fe2

 

где a1 = H-yi-, "3- HyS

Используя div .4=0, можно определить À

в виде

л * = - \ - д ? * у -

( 2 - 2 5 >

Решив систему (2. 24), получим выражения для составляющих векторного потенциала в крайнем пакете в виде

 

Л « і =

2

Ь1г-*>"* +

Fi,/1"*)

X

I

 

 

 

и=1, 3, S, ...

 

 

 

 

 

 

 

 

X

- —

" 1 7

COS K„X

sin

ß,,;/,

 

 

 

 

 

 

COS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.

26)

 

 

я=1, 3, 5,...

 

 

 

 

 

 

 

 

; я

jr.x

Jtb

\„ sm

Л„і

 

 

 

 

X

- —

- —

cos p„j/;

 

 

-FT-

£

— £

Pn

COS

À„ù

 

 

 

 

 

 

 

 

 

 

 

здесь

8 І . = } / « І [ В +

(Т)2-

p

J '

 

£ n „ = _ ^ü. Eô]'„Aj

 

7t\2

. . °зя Н-уЛ

 

1 4- ~

ch f'ft,, 4-

^ V «in

f*0

^ 7„

ИуЗ /

'" 2 J '

in о

 

 

 

 

- ( i - l ^ c h , , ; / , ,

62

+

& 3 И JJLQ

 

 

 

 

 

 

(2. 27)

 

 

 

 

 

 

 

 

 

 

 

 

 

р

 

/^, Л2 — см. рис. 2-12, (i1 }

f».3

— магнитные проницаемости по оси у

в первом и втором пакетах.

 

 

 

 

 

Комплексная

амплитуда

аксиальной составляющей индукции

в крайнем пакете

 

 

 

 

 

 

 

 

 

 

 

yi

xl

 

 

 

 

в, = - дх

ду

 

 

 

 

 

 

 

 

2

1

1

 

2

ІЕГ

+

FLN/^)

[ ( ^ )

-p- +

ß„J e ~ ~ cos P J , .

(2. 28)

я=1,3,5, ...

 

 

 

 

 

 

 

 

Потери, выделяемые в крайнем пакете от проникновения акси­

альной составляющей

индукции поля лобовых частей, равны

 

 

 

4/2

О

Л,

 

 

 

 

9кр =

Р S

SS

\S1\^dxdydz.

(2. 29)

 

 

 

О -1/2 О

 

 

 

Здесь

 

 

 

 

 

 

 

 

 

| ^ | 2 = f 4 а • + | f Äyl

Суммарные потери во всех зубцовых частях крайнего пакета

на две стороны

турбогенератора равны

 

(2. 30)

где Z1 — число

пазов статора.

Рассмотрим несколько подробнее вопрос о выборе эквивалент­ ных величин магнитной проницаемости пакетов (х^, цу 3 , р л и (л_3.

Будем считать постоянные магнитные проницаемости, входя­ щие в (2. 28), эквивалентными, если затухание модуля индукции по (2. 28) совпадает с затуханием индукции по (2.16) в нелинейной среде при наличии основного потока.

С этой целью будем определять эквивалентную магнитную про­ ницаемость из равенства выражений (2. 28) и (2.16) с учетом (2.17). Такое определение в настоящей работе выполнено расчетным пу­

тем, причем в соответствии с (2. 17) было принято

f i y = f i . . По вы­

ражению (2. 28) для различных значений р л и цг 3

рассчитывались

кривые распределения аксиальной составляющей

магнитной ин-

63

дукции по глубине пакета, причем одновременно исследовалось

влияние изменения полюсного деления т и вида кривой Be—f

(у)

на характер зависимости B.=f (z). Как показали результаты

ра­

счета, последние два фактора практически мало сказываются на

глубине проникновения Вг.

В то же время магнитная

проницае­

мость

f i t l

первого пакета, особенно в пределах до ^ і = Ю ц 0 , резко

 

 

 

 

 

 

 

 

изменяет

крутизну затухания

индук­

 

 

 

 

 

 

 

 

ции.

Магнитная

проницаемость

^ 3

 

 

 

 

 

 

 

 

второго пакета в незначительной сте­

 

 

 

 

 

 

 

 

пени влияет

на затухание

Ве

в пер­

 

 

 

 

 

 

 

 

вом пакете. С другой стороны, по (2.

 

 

 

 

 

 

 

 

16) с учетом

(2.

17)

рассчитывались

 

 

 

 

 

 

 

 

кривые B,=f

(z) для стали ЭЗЗО. Рас­

 

 

 

 

 

 

 

 

четы показали, что магнитное состоя­

 

 

 

 

 

 

 

 

ние зубца крайнего пакета из-за

 

 

 

 

 

 

 

 

наличия основного потока с учетом

 

 

 

 

 

 

 

 

слоистости среды характеризуется ма­

 

 

 

 

 

 

 

 

лой величиной магнитной

проницае­

 

 

 

 

 

 

 

 

мости [ла =(2-^-4) р.0 , причем нижний

 

 

 

 

 

 

 

 

предел относится к режиму

холостого

 

 

 

 

 

 

 

 

хода,

а верхний — к режиму

корот­

 

 

 

 

 

 

 

 

кого замыкания. Последние величины

 

 

 

 

 

 

 

 

и были приняты в качестве эквивален­

 

 

 

 

 

 

 

 

тных значений р л при расчете

прони­

 

 

 

 

 

 

 

 

кновения

В.

в крайние

пакеты

сер­

 

 

ч

 

12

 

 

 

дечника.

 

На

рис. 2-12

приведены

 

 

40

80

100

э к с п е

Р и

м

е и т а л ь н

ы е

кривые

B,=f

(z)

 

 

20

50

 

 

 

для

турбогенератора

мощностью

Рпс. 2-13. Распределение

500 Мвт, измеренные

при номиналь­

квадрата

модуля

составляю­

ном токе статора и холостом ходе при

щих плотности тока по осп у

номинальном

напряжении

на высоте

(—) и z ( — ) на торцовой по-

примерно

90 мм от расточки статора.

поверхности крайнего

пакета.

Там

же

 

приведены

данные

расчета

1 — \Si/\2—f

(у)

при к = 2 0 мм; 2 —

 

В г = = І

(z) по (2.28)

при

ft,1 =

2fi0

І ^ І 2 = /

(1/)прнх=0; з — |Sy| a =/ (к)

для режима

холостого

хода

и р-л=

прп

у = 100

мм; 4 y\*=f

(ж)

при у = 6 0 мм; 5 — \&x\*=f

(ж) при

=4ft)

для

режима

короткого

за­

у = 100

мм;

6 — | Й Ж | 2 = / (ж)

прп

мыкания. Как следует

из

рис. 2-12,

у = 6 0

 

мм;

7 — \è\-=i

(у)

 

прп-

 

 

опытные и расчетные кривые Bz=f

(z)

х=0;

8 — |ая |2

=/ (у)

при .-с=4 мм.

 

 

 

 

 

 

 

 

хорошо согласуются друг с другом.

По

выражениям (2. 26), (2. 29)

и (2. 30)

на ЭВЦМ «Минск-22»

были выполнены расчеты, распределения квадрата

модуля

состав­

ляющих плотности тока | 5 у | 2

и \SX\2 на торцовой поверхности край­

него пакета (рис. 2-13), а также добавочных потерь в крайнем пакете применительно к режиму трехфазного короткого замыка­ ния турбогенератора мощностью 500 Мвт.

Для турбогенератора мощностью 500 Мвт полная величина по­ терь от торцового потока в крайнем пакете в режиме короткого

64

замыкания составляет: по (2.30) — 44 квт, по [81] — 15.1 квт, по [17] — 34.7 квт, по [161] — 36 квт.

Потери по сечению зубца распределены крайне неравномерно. Кривые распределения квадратов составляющих плотности тока (рис. 2-13) показывают, что вихревые токи, наведенные в зубце,

вытесняются

к его краям.

 

 

 

 

 

 

 

Т а б л и ц а

2-1

 

Распределение потерь по глубине пакета

 

 

г, мм

1 *.-с Imaxi

1 By |шах>

Чх, вт

 

 

ѵ = ± * = о,

* = -£-, » = о,

Чу,

BT

 

a'/MM1

а'/мм1

 

 

 

0

3.23

20.2

1.66

80.7

 

4

2.03

16.2

1.04

65.1

 

8

1.34

13.2

0.682

53.4

 

12

0.922

11.1

0.465

44.6

 

16

0.658

9.49

0.330

38.2

 

20

0.488

8.32

0.244

33.5

 

24

0.376

7.47

0.187

30.0

 

28

0.300

6.84

0.149

27.4

 

32

0.247

6.36

0.124

25.4

 

36

0.209

5.98

0.106

23.8

 

40

0.181

5.63

0.093

22.4

 

В табл. 2-1 приведены результаты расчета квадрата модуля составляющих I ^ J 2 и І ^ І 2 и соответствующие им потери по глу­ бине крайнего пакета применительно к режиму трехфазного ко­ роткого замыкания турбогенератора мощностью 1200 Мвт.

Потери в крайнем пакете на две стороны турбогенератора мощ­ ностью 1200 Мвт составят 130 квт.

Для выявления зависимости потерь в зубцах статора от элект­ ромагнитных нагрузок и геометрических размеров были выполнены специальные исследования* потерь на макете, который представлял собой шихтованный сердечник из электротехнической стали с двумя прилегающими подвижными ярмами из того же материала. По­ тери исследовались в пакете, состоящем из набора пластин из электротехнической стали различных марок. Исследуемый пакет помещался между ярмами. Торцовый поток создавался обмоткой воз­ буждения, уложенной внутри сердечника. Потери от торцового потока определялись как разность общих потерь в магнитопроводе макета при наличии исследуемого пакета и при замене пакета тарировочным образцом, имеющим продольную шихтовку пластин.

* Исследования проведены пнж . В . Д . Клпмовицкпм и И. G. Генендером (СибНИЭТИ).

5 Я. Б. Данилевич

65

Результаты исследований показали, что потери от торцового потока зависят от квадрата индукции торцового поля (рис. 2-14).

Сопоставление потерь при одинако-

Рис . 2-14. Зависимость удельных потерь от квадрата пндукцип торцового поля (а) и его частоты (б).

J — пакет из стали 912 (р=2.7 - Ю - 5 ом -см); 2 — пакет из стали Э310 (р = 4.8 -І0~1

ом-см);

 

 

 

 

ширина

пакетов одинакова.

 

 

 

 

На

рис. 2-14, б показана

ісимость удельных потерь от ча-

стоты.

Из

рисунка

следует,

потери пропорциональны

квад­

 

 

 

 

 

 

рату частоты. Последнее

также

0.50 Ь

 

 

 

 

 

согласуется с формулой (2. 24)

 

 

 

 

 

 

и свидетельствует

о том, что по­

 

 

 

 

 

 

тери имеют

вихревой характер

 

 

 

 

 

 

и что

отсутствует заметное эк­

 

 

 

 

 

 

ранирование

торцового

потока

 

 

 

 

 

 

вихревыми токами в пакете.

0.1о Y

 

 

 

 

 

На

макете

исследовалась,

 

 

 

 

 

кроме

того,

зависимость

удель­

 

 

 

 

 

 

I

 

I

I

8,мм

ных потерь

от ширины

пакета.

0

 

10

20

30

Результаты

исследований

по­

 

 

 

 

 

 

Рис . 2-15. Зависимость удельных

казали, что

удельные

потери

потерь от ширины пакета при индук­

при заданной величине

квадра­

ции торцового потока 1.4 тл.

 

та индукции примерно пропор­

1 — пакет

из

стали

Э12;

г — пакет из

циональны

ширине

пакета

 

 

стали Э310.

 

 

(рис. 2-15).

 

 

 

 

2.6. Влияние расщепления зубцов на потери

В настоящее время турбо- и гидрогенераторы большой мощности выполняются с прорезью в зубцах крайнего пакета сердечника ста­ тора. Как показали экспериментальные исследования, это дает воз­ можность уменьшить величину потерь в зубцах в целом, а также снизить местные потери и нагревы, которые в некоторых случаях

66

могут ограничить работу мощных машин, особенно в режимах недовозбуждения.

Рассмотрим несколько подробнее влияние прорези и ее разме­ ров на величину местных потерь в сегментах крайнего пакета.

Для исследования влияния прорези в зубцах на величину местных потерь достаточно рассмотреть их распределение лишь

всегменте крайнего пакета сердечника. При этом можно ограни­ читься расчетом потерь только в пределах половины зубцового деления. Заменив сегмент его разверткой

впрямоугольной системе координат (рис.

2-16), для поля

вихревых

токов

в листе

имеем

dB.

 

 

 

1

 

 

 

rob,S = — s r ,

d\\-S

= 0.

(2.31)

Р

dt

 

 

 

Здесь Д . может быть представлено в виде

бегущего поля на поверхности

листа

B* = Re ßг (0. У)~

(2. 32)

Заменяя бегущее поле двумя пульси­ рующими и учитывая, как и выше, что в пределах зубцового деления ампли­

туду первого из полей с достаточной точностью можно считать неизменной

[Ree

-1

i - o ) ,

а второго

малой

Рис. 2-16. К расчету распределения местных потерь в сегменте край ­ него листа сердечника статора.

ющей

о )

для

аксиальной

составля­

зубца,

получим

индукции,

падающей

 

на

поверхность

 

 

 

 

В, = Вг(0,

у)

cos

ut.

 

 

Введем функцию и,

удовлетворяющую уравнению

 

 

 

 

 

rot

и =S.

 

 

 

(2.33)

Будем считать, что вихревые токи в сегменте сдвинуты по фазе

относительно

индукции

на

90°. Тогда

система

(2.31)

сводится

к уравнению

 

 

 

 

 

 

 

 

 

 

 

 

7ut

+ ^ i

=

- J ß A 0 , u ) .

 

(2.34)

В соответствии с (2.27) найдем составляющие плотности вих­

ревых

токов

 

 

да.

 

 

 

diij.

 

 

 

 

 

 

 

s ,

= -

 

(2. 35)

 

 

 

 

1у'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5*

67

Для решения (2.34) могут быть приняты следующие граничные условия:

1) на всех границах листа нормальная к границе составляющая плотности тока равна нулю (и.=0);

2) на осях симметрии зубца {CD, рис. 2-16) и паза {AB, рис. 2-16) равна нулю радиальная составляющая плотности тока

Местные потери в листе вычислялись по амплитудному значе­ нию плотности тока и равны

« 7 = f {S%+ s»).

(2. 36)

Решение (2. 34) с учетом граничных

условий проводилось на

ЭЦВМ «Минск-22» методом сеток, причем в качестве начального значения функции иг{х, у) принималось аналитическое решение

уравнения (2. 34) для прямоугольной области

(приложение 3).

На рис. 2-17 применительно

к турбогенератору

мощностью

500 Мвт представлены

кривые

распределения

плотностей

токов

и местных потерь вдоль

ширины

зубцовой зоны

при

к =

^<^1,

 

 

 

 

 

'II

где к — отношение высоты шлица к высоте паза. Полученное распре­

деление токов в листе полностью

согласуется с физической карти­

ной растеканпя токов при поверхностном эффекте. Согласно урав­ нению (2. 34), начало координат может занимать любое положение

на линии верхушки зубца. Тогда,

учитывая направление токов,

а также то, что для коронки зубца

^> Sx (рис. 2-17), получаем

линии тока с небольшим угловым коэффициентом tga, где а — угол между касательной к линии тока и осью у.

Линии тока оказываются симметричными относительно се­ редины расстояния от стенки паза до стенки шлица. В области коронки зубца образуется сильное сгущение линий тока, которое усиливается из-за большой величины индукции, падающей на ко­ ронку. У конца рассечки зубца картина растекания вихревых то­ ков несколько изменяется, линии тока становятся несимметрич­ ными (рис. 2-17). Угловой коэффициент линий тока непрерывно

возрастает от стенки паза

к середине зубца.

В результате потери

в области между концом

рассечки и дном

паза возрастают, что

приводит к увеличению нагрева крайних листов пакета статора

у дна паза.

Последнее также связано с ухудшением охлаждения

листа

вблизи

обмотки.

 

На

рис. (2-18)-f-(2-20) показано влияние

высоты расщепле­

ния зубца на распределение местных потерь

в различных сече­

ниях по высоте сегмента крайнего пакета сердечника статора турбо­ генератора мощностью 500 Мвт.

68

Потери на стенке паза убывают вдоль высоты зубца (рис. 2-18), причем характер убывания величины местных потерь практически

не

зависит

от величины к : <?кор >

Отса > Оу ,„

где QKW

— потери

по

коронке

зубца, Qmcc — потери

по рассечке

зубца,

Оу „ — по­

тери около

угла лаза.

 

 

 

.<

ö

 

 

а

 

fr, -

«а К

 

СО

со

 

300 -3000-300

200 -2000-200 -

100-woo -100

-1000 -2000

-3000

вп/2

45

 

 

 

б

200

2000

 

 

 

 

 

 

V

/

100

1000-

Ч/2

 

-1000

. впіг

 

-2000 -

/

 

Рис . 2-17. Распределение плот­ ностей тока Sx (І) Ii S g (2) и местных потерь (3) вдоль ши­ рины зубцовой зоны сегмента в различных сечеппях по высо­

те зубца при

/с=0.73.

а — вдоль коронки

зубца; б —

вдоль 3/4 высоты зубца; в — вдоль линии около дна паза.

Вш/2

в

60

1000

to

20 500 1 ѵч 3

2 /

-500

Infi

вш

Расчеты показывают (рис. 2-18), что увеличение отношения с /«=0.75 до к=1 позволяет снизить потери в зоне прежнего конца рассечки приблизительно в 1.4 раза. При этом нагрев должен снизиться в еще большей степени из-за улучшения вентиляции по щели. Дальнейшее увеличение высоты рассечки приводит к сни­ жению потерь только в области дна паза, не оказывая заметного

69

Соседние файлы в папке книги из ГПНТБ