Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Мат.мет.АРЭС МАРЭС / Методические указания к выполнению курсовой работы / Методические указания к выполнению курсовой работы.pdf
Скачиваний:
39
Добавлен:
28.05.2015
Размер:
758.92 Кб
Скачать

10.

1.2. Математическая модель электромагнитного демпфера

Выберем систему отсчета так, как показано на рисунке 2: Ось ОХ направим вниз, приняв за начало отсчета уровень внутреннего угла верхнего полюса (пунктирная линия). Координату сердечника х будем определять расстоянием от этого уровня до его верхнего торца. Причем если верхний торец находится выше нулевого уровня, то координата сердечника отрицательная. Так, на рисунке 2, х0 < 0.

Основные геометрические размеры демпфера приведены на рисунке 3. К ним относятся: диаметр сердечника DС, внутренний диаметр расточки индуктора D, воздушный зазор z (на одну сторону), высота hП нижнего и верхнего полюса индуктора, глубина bК и высота hК паза под катушку, число витков катушки w.

hП

hП

D

z

O

x

X

DС

bK

w hK

Рис. 3. Основные размеры демпфера

Для упрощения математической модели предположим, что магнитная проницаемость стали равна бесконечности, а электрическая проводимость рав-

© А. Кудинов, 2006, ТГУ, каф.ПЭ

11.

на нулю (т.е. нет вихревых токов). На практике это может быть обеспечено применением новых прессованных ферромагнитных материалов. Длину сердечника будем считать достаточно большой для полного перекрытия обоих полюсов индуктора. Предположим также, что магнитная цепь линейна, что справедливо в силу наличия воздушного зазора. В этом случае электромагнитную силу, развиваемую демпфером в направлении оси ОХ можно определять по известной формуле

 

F(i, x) =

1

(iw)2

 

G(x)

,

(1)

 

2

x

 

 

 

 

 

 

где F

электромагнитная сила;

 

 

 

 

 

 

i

ток катушки;

 

 

 

 

 

 

w

число витков в катушке;

 

 

 

 

 

 

G(x) – зависимость проводимости магнитной цепи от положения сердечника.

Примерный вид зависимости G(x) и F(x) при i = const приведен на рисун-

ке 4.

Основные свойства функции G(x) ясны из физических соображений. Если смещать якорь из положения x = 0 вверх, т.е. в область отрицательных координат, проводимость магнитной цепи должна возрастать и достигать максимального значения Gm при x = –hП. Дальнейшее уменьшение x практически не влияет на значение магнитной проводимости. Смещение якоря вниз, в область положительных значений координат, вызывает уменьшение магнитной проводимости. При x = hК + hП якорь выходит из индуктора полностью и магнитная проводимость минимальна (G0). Дальнейшее увеличение координаты x практически не влияет на проводимость. Наибольший наклон функции G(x) должен наблюдаться при значении x близком к нулю.

© А. Кудинов, 2006, ТГУ, каф.ПЭ

12.

 

Gm

G(x)

 

G0

x

0

 

F(x)

 

Рис. 4. Зависимость проводимости магнитной цепи от положения сердечника

Аналитическое определение зависимостей G(x) и F(x) – весьма сложная задача, к тому же погрешность имеющихся метедов расчета сложных магнитных цепей достаточно велика. Поэтому целесообразно аппроксимировать зависимость G(x) некоторой простой аналитической функцией. Предлагается для этого использовать зависимость вида

G(x) = G

 

 

γ

 

+1

(2)

 

 

ax

0

+e

 

 

 

1

 

 

 

где G0 – минимальная магнитная проводимость (при выведенном сердечнике);

γ = Gm G0 – коэффициент модуляции проводимости;

G0

Gm – максимальная магнитная проводимость (при полностью введенном сердечнике);

a – параметр, определяющий наклон графика G(x).

© А. Кудинов, 2006, ТГУ, каф.ПЭ

13.

Расчет минимальной G0 и максимальной Gm проводимости магнитной цепи особой сложности не представляет.

Для вычисления G0 рассмотрим магнитную цепь обмотки при выведенном сердечнике (рисунок 5). Выделим три характерные зоны прохождения магнитного потока: I – зона, занятая обмоткой, II – зона выпучивания с углов полюсов, III – зона выпучивания с боковых поверхностей полюсов. В соответствии с этим делением будем искать три слагаемых проводимости G0:

G0 = G1 +G2 +G3

(3)

Для расчета составляющей G1 магнитной проводимости участка, занятого обмоткой, выделим на расстоянии r от оси симметрии тонкую трубку потока толщиной dr (рис. 5).

bК

 

 

 

 

 

 

hК

I

II

III

dr

r

 

 

 

 

 

D

 

D1

 

 

Рис. 5. К расчету проводимости G0

Силовая линия, проходящая через выделенную трубку, охватывает часть витков wr обмотки, равную

wr = w D1 2r . D1 D

© А. Кудинов, 2006, ТГУ, каф.ПЭ

14.

Пренебрегая магнитным сопротивлением стали, по закону полного тока рассчитаем напряженность Hr и индукцию Br магнитного поля на воздушном участке силовой линии при некотором токе i обмотки:

Hr = i hwr , Br = µ0Hr

к

Площадь поперечного сечения трубки

dS = 2πr dr =πr d(2r)

и соответственно поток выделенной трубки

= Br dS = µ0 π i wr r d (2r) hк

Поток обусловливает составляющую потокосцепления dΨ, равную

dΨ = w

=

µ0π

i w 2r d(2r) =

µ0π

i w2

 

D1 2r

2 r d(2r) .

h

 

 

 

r

 

r

h

D

D

 

 

 

к

 

к

 

1

 

 

Интегрируя полученное выражение в пределах от D до D1, и учитывая, что маг-

нитная проводимость G =Ψ (i w2 ) , находим искомую составляющую прово-

димости

 

µ0π

D1

 

D1 2r 2

2r

 

µ0π

D12

D1 D D2

G1 =

h

 

 

 

 

 

 

d(2r) =

h

 

 

+

 

 

 

D

D

2

24

12

8

 

к

D

 

1

 

 

 

 

 

к

 

 

 

 

 

 

 

Учитывая, что D1 = D +2bк , последнее выражение можно упростить:

G =

µ0π

 

bк (2D +bк )

(4)

 

 

1

6

 

hк

 

 

 

 

Проводимости областей II и III вычислим, пользуясь рекомендациями, изложенными в [1]:

G2 = 0.815µ0D

 

G3 =

 

2µ0D

 

(5)

1+

hк

 

 

 

h

 

 

 

 

 

 

 

 

п

 

© А. Кудинов, 2006, ТГУ, каф.ПЭ

15.

Максимальную проводимость Gm найдем как сумму проводимости пазового рассеяния G1 и проводимости GС потоку ФС, проходящему через сердечник (рисунок 6).

Gm =G1 +GC

(6)

hП

 

Ф1

 

hК

 

 

ФС

hП

 

bК

DС

 

z

Рис. 6. К расчету проводимости Gm

В силу допущения о том, что магнитная проницаемость стали бесконечна, проводимость GС будет определяться только воздушными участками магнитной цепи. Общая протяженность воздушных участков на пути потока ФС равна удвоенному зазору z, поэтому

 

 

 

 

G =

µ0π (DC + z)hп

 

 

 

 

 

 

 

(7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

2z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т.о., на основании формул (3)…(7) запишем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ0π

 

bк (2D +bк )

 

 

 

 

 

 

2

 

 

 

G

=

 

+ µ

 

D

0.815

+

 

 

(8)

 

 

 

 

hк

 

0

 

6

 

hк

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

© А. Кудинов, 2006, ТГУ, каф.ПЭ

 

 

 

 

 

 

 

16.

G

= µ π

 

bк (2D +bк )

+

(DС + z)hп

 

(9)

 

 

 

m

0

 

6hк

 

2z

 

 

 

 

 

 

 

Последним неизвестным параметром в формуле (2) является наклон a функции G(x). Этот параметр определяет наибольшее значение производной магнитной проводимости, и максимальную электромагнитную силу, развиваемую демпфером:

 

G(x)

 

G γa eax

 

 

 

G γa

 

(Gm G0 )a

 

 

= −

0

 

 

= −

0

=

 

(1+eax )2

 

 

 

 

 

x max

 

 

x = 0

 

4

 

4

Следовательно, по формуле (1) максимальная сила

Fmax = 81 (iw)2 (Gm G0 )a

С другой стороны, по формуле Максвелла,

Fmax = BС22µSС

0

где SС – площадь поперечного сечения сердечника,

BС – средняя индукция в сердечнике при его координате x = 0:

B

=

iw

(G(0)

G )=

iw

 

 

Gm +G0

G

 

 

 

 

 

С

SС

 

1

SС

 

 

2

1

 

 

 

 

 

 

 

 

(10)

(11)

(12)

В формуле (12) поток сердечника ФС = BCSC вычисляется как разность общего потока обмотки (iwG(0)) и потока пазового рассеяния (iwG1 ). Прирав-

нивая (10) и (11), с учетом (12), найдем значение параметра a:

 

(G

+G 2G )2

 

a =

m

 

0

1

 

(13)

µ S

С

(G

G

)

 

 

 

0

m

0

 

 

Формулы (8, 9, 13) позволяют рассчитать параметры аппроксимирующей функции G(x) (2) по геометрическим размерам демпфера и использовать выражение (1) для вычисления мгновенного значения его электромагнитной силы.

[Содержание]

© А. Кудинов, 2006, ТГУ, каф.ПЭ