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

книги / Моделирование цилиндрических линейных вентильных двигателей для различных отраслей промышленности

..pdf
Скачиваний:
5
Добавлен:
12.11.2023
Размер:
10.95 Mб
Скачать

I (F)

 

 

 

 

 

 

 

F

 

 

 

 

 

k

ф

B

 

2 p D

3 w

 

 

 

 

 

m

 

 

2

 

п

(3.35)

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

6 k

ф

 

B

p D

w

 

 

 

 

 

 

 

 

m

2

п

 

 

Рис. 3.6. Зависимости усилия от тока

На рис. 3.6 приведены графики зависимости усилия от тока, полученные: по формуле (3.35) для расчёта методом электрических и магнитных цепей с учётом насыщения (1) и без учёта насыщения (2); на расчётной модели ANSYS Maxwell при неподвижном вторичном элементе (3) и при частоте работы двигателя 7 Гц (4). Гра-

61

фик 5 получен экспериментально. Значения максимального усилия, развиваемого двигателем при действующем токе в 30 А, приведены в табл. 3.4.

Таблица 3 . 4

Максимальное усилие прототипа ЦЛВД при действующем значении тока 30 А, полученное разными методами расчёта

Тип расчёта

Усилие, Н

п/п

 

 

1

Метод электрических и магнитных цепей с учётом на-

3274

 

сыщения магнитной цепи

 

 

2

Метод электрических и магнитных цепей без учёта на-

3450

 

сыщения магнитной цепи

 

 

3

ANSYS Maxwell magnetostatic при постоянном токе и

4384

 

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

 

 

4

ANSYS Maxwell transient при движущемся вторичном

3264

 

элементе и частоте питающего напряжения 7 Гц

 

 

5

Экспериментальное значение

3940

Как видно из табл. 3.4, наиболее точный результат дает метод электрических и магнитных цепей, учитывающий насыщение магнитной цепи и влияние токов в обмотках статора на магнитный поток, и расчёт методом конечных элементов в ANSYS Maxwell magnetostatic. График 3 идёт параллельно экспериментальному (5), потому что экспериментальная зависимость представляет разницу электромагнитного усилия двигателя и силы трения, когда расчётные методы силу трения не учитывали. Однако токовая характеристика двигателя, полученная экспериментально не точно отображает работу двигателя, так как снималась при постоянном токе. Как описано в п. 2.3, замер усилия по максимальному значению в ANSYS Maxwell magnetostatic соответствует пиковому значению усилия при наличии биений и зависит от фазы токов.

Таким образом, можно сделать вывод, что, несмотря на то, что расчёт усилия в Maxwell magnetostatic и по формуле (3.31) соответствует эксперименту, они не подходят для построения

62

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

После расчёта магнитной цепи можно рассчитать продольное и поперечное индуктивныесопротивления, используя формулыиз [72].

Коэффициенты продольной и поперечной реакции якоря следующие:

kфd

sin( kф / 2)

,

kф sin( kф

/ 2)

 

 

 

kфq

 

(1 cos( kф / 2))

,

 

kф sin( kф

/ 2)

 

 

 

ad

1

, aq

1

,

 

 

Rd

 

 

 

 

 

 

 

Rq

 

 

(3.36)

(3.37)

(3.38)

где для точного расчёта сопротивление Rd рассчитывается согласно рис. 3.5, а сопротивление Rq рассчитывается согласно рис. 3.7, для

чего создаётся отдельная функция рrog(FI , з0 , k0 , R , R ), которая

посредством итерационного расчёта определяет поперечные магнитные потоки.

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

L

 

m w2

k2

 

ad

 

2

p

 

 

об

 

 

0

 

 

 

 

d

 

 

p kфd

 

 

 

 

 

 

 

 

 

 

 

 

 

m w2

k 2

 

aq

 

 

 

L

 

 

об

 

2

p

 

 

 

 

 

0

q

 

 

p kфq

 

 

 

 

 

 

 

 

 

 

 

wп2 l п,

wп2 l п ,

где п – относительная пазовая магнитная проводимость,

п

 

h1

h2

 

2h3

3

b1

b1 b2

 

b1

 

(3.39)

(3.40)

(3.41)

63

где

h b wп ;

h

D1 D2

b

h

 

т

h

;

h b ;

h 0;

 

 

1

2

2

2

k

я

 

1

 

3

k

4

 

 

 

 

 

 

 

 

 

 

 

 

 

b

b ;

b b .

 

 

 

 

 

 

 

 

 

 

1

3

3

2

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.7. Схема для расчёта Rq : а – магнитная цепь;

б– эквивалентная электрическая цепь

Спомощью продольной и поперечной индуктивностей можно выразить зависимость индуктивных сопротивлений от частоты:

Xd ( f ) 2 f Ld ,

(3.42)

Xq ( f ) 2 f Lq .

(3.43)

При заданных техническим заданием, мгновенной скорости υ и тяговом усилии F можно рассчитать необходимые для подачи на двигатель величины напряжения и его частоты.

Система уравнений (согласно векторной диаграмме рис. 3.8):

64

2 f ;

 

(3.44)

 

Pэм

 

m E0 Is cos( )

 

 

F

 

;

(3.45)

 

 

 

 

 

2 f

 

 

 

 

 

 

 

 

(3.46)

U sin( ) Iq Xq Id r;

 

 

 

 

 

 

 

(3.47)

U cos( ) Id Xd Iq r E0 ;

 

 

sin( );

 

(3.48)

Id Is

 

Iq Is

cos( );

 

(3.49)

 

 

 

 

 

 

 

 

 

2

2

 

(3.50)

Is

Id

Iq ;

 

 

 

 

 

 

 

 

где Is – токстатора, который связан сдействующимтокомфазыI как:

Is

3

I .

(3.51)

2

 

 

 

Техническим заданием нам заданы значения тягового усилия и мгновенной скорости для каждого момента времени работы двигателя, а мгновенное значение температуры получено от датчика температуры, но для расчёта зададимся температурой 100 °С.

Выразим из (3.45) и (3.49) по-

перечную составляющую тока Iq :

Iq

2 f F .

(3.52)

Рис. 3.8. Векторная

диаграмма ЦЛВД

 

m E0

 

 

Далее по векторной диаграмме найдём продольную составляющую тока Id , предварительно рассчитав I по формуле (3.34):

Id

Is2 Iq2 .

(3.53)

Согласно (3.46) и (3.47) напряжение, подаваемое на статор,

65

U (U sin( ))2 (U cos( ))2 ,

U (I

q

X

q

I

d

r)2 (I

d

X

d

I

q

r E )2 .

(3.54)

 

 

 

 

 

 

0

 

Согласно формуле (3.54) рассчитана зависимость фазного напряжения питания ЦЛВД от частоты при номинальном значении действующего тока 30 А. Эта зависимость показана на рис. 3.9.

Рис. 3.9. Зависимость напряжения питания от частоты при заданном токе обмотки статора 30 А

Для расчёта зависимости активного сопротивления от температуры необходимо высчитать линейные размеры проводника. Зная число витков в одном пазу и их расположение, можно рассчитать суммарную длину проводников на одну фазу статора:

 

 

 

 

I ps p w l ,

(3.55)

где средняя длина витка обмотки l dср ,

где средний диаметр

витка d

ср

D 2

(b ) a

wп

.

 

 

 

 

2

k

4

 

 

 

 

 

 

 

 

Для точности прибавим к этой величине длину статора двигателя, как длину проводников, соединяющих катушки одной фазы между собой:

66

lp lps p .

(3.56)

Далее найдём активное сопротивление обмоток одной фазы двигателя при температуре 20 ºС согласно закону Ома:

R

lp

.

(3.57)

20 С Sr

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

r(t°) R20 °C (1 ) .

(3.58)

3.3. Расчёт угловых и рабочих характеристик ЦЛВД по предлагаемой методике

При построении угловых характеристик задаёмся постоянными значениями величины и частоты подаваемого напряжения, а угол нагрузки изменяем в диапазоне от –90 до 270 градусов. Продольная и поперечная составляющие токов статора вычисляются по форму-

лам [72]:

Id

 

Xq (U cos( ) E0 ) rs U sin( )

,

(3.59)

X

d

X

q

r2

 

 

 

 

 

 

s

 

 

 

I

q

 

rs (U cos( ) E0 ) Xd U sin( )

.

(3.60)

 

 

 

 

X

d

X

q

r2

 

 

 

 

 

 

 

 

 

s

 

 

 

Для выполнения расчёта с учётом насыщения магнитной цепи

токи Id и Iq

 

 

рассчитываются

 

соответствующей

функцией

сurrent(U , , f ) .

Усилие, развиваемое двигателем, можно упрощённо рассчитать по формуле (3.45). Но эта формула может быть использована лишь для двигателей, у которых продольное и поперечное индуктивные сопротивления примерно равны, поэтому для расчёта угловых характеристик двигателя воспользуемся формулами [72].

67

Максимальные основное усилие F m и реактивное усилие Fdqm рассчитываются по формулам

 

 

 

 

 

F m

 

 

 

 

 

 

 

 

 

 

3 U E0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 f

(X

d

X

q

 

r

2 )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

(X

d

X

2

 

2 r

2

X

d

r2

 

 

X

q

)2

 

 

r

2 (2 X

2

X

d

 

X

q

r2 )2 ,

 

 

q

 

s

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

s

 

F

 

 

3 U 2

(X

 

d

X

q

)

 

 

 

 

 

(X 2

 

 

r2 )

(X

2 r2 ) .

 

4

f (X

 

 

X

 

r

2 )2

 

 

 

 

 

 

dqm

 

d

q

 

 

 

 

 

 

 

 

d

 

 

 

 

 

s

 

 

 

 

q

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тормозные основное усилие F Т и реактивное усилие FdqТ оп-

ределяются по формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

3 E2

r

 

 

 

 

 

 

 

 

 

 

 

 

(X 2

r2 ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

f

(X

 

 

 

X

 

 

r

2 )2

 

 

 

 

 

 

 

 

 

Т

 

d

 

q

 

 

 

 

 

 

 

q

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

3 U 2

 

r

 

 

 

 

 

 

 

 

 

 

 

 

(X

 

 

X

 

 

)2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

f (X

 

 

 

X

 

 

 

r

2 )2

 

 

 

 

 

 

 

 

 

 

dqТ

 

d

q

 

 

 

 

 

 

d

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Углы смещения основного

 

 

 

и реактивного dq усилий рас-

считываются по формулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

(2

X

q

2

 

X

d

 

X

q

r2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

,

 

 

 

 

 

 

 

 

X

 

X

 

 

 

 

2

r2

 

X

 

 

 

r2

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

q

2

 

 

d

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

(X

d

 

X

q

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

arctg

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

X

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dq

 

 

 

 

 

 

 

 

 

 

 

d

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полные основное усилие F и реактивное усилие Fdq определяются по формулам:

F F m sin( ) F Т ,

Fdq Fdqm sin(2 ( dq )) FdqТ .

68

В итоге получим значение синхронного усилия Fс :

 

Fс F Fdq .

(3.61)

Проанализировав графики 13 на рис. 3.10, можно сделать вывод, что целесообразнее вести расчёт усилия по формуле (3.45), а не по (3.61), потому что график по формуле (3.45) заметнее сходится с графиком, полученным методом конечных элементов, который считается наиболее надёжным, так как учитывает все процессы, происходящие в двигателе.

Рис. 3.10. Угловые характеристики ЦЛВД при частоте питающего напряжения 10 Гц, полученные расчётами: 1 – с учётом насыщения магнитной цепи; 2 – без учёта насыщения магнитный цепи;

3 – методом конечных элементов

На графиках можно заметить, что максимальное усилие в генераторном режиме становится больше, чем в двигательном. Причиной тому влияние активного сопротивления. В идеализированной синхронной машине, активное сопротивление которой крайне мало по сравнению с индуктивным, с ростом частоты пропорционально растёт и ЭДС и общее сопротивление обмотки. Но так как ЦЛВД работает на низких частотах, то его активное сопротивление иногда превышает индуктивное, потому в рабочем диапазоне частот можно наблюдать этот эффект на угловую характеристику [69]. Если под-

69

нять частоту питающего напряжения заметно больше, например до промышленной, то эффект будет незаметен, так как индуктивное сопротивление стало намного больше активного. Наглядное сравнение угловых характеристик ЦЛВД частот 1, 10 и 50 Гц представлено на рис. 3.11.

Рис. 3.11. Угловые характеристики ЦЛВД для частот питающего напряжения 1, 10, 50 Гц

На рис. 3.11 видно, что угол нагрузки θ, при котором ЦЛВД достигает максимальное усилие, имеет разные значения в зависимости от величины частоты питающего напряжения, это также обусловлено влиянием активного сопротивления. При частоте в 1 Гц двигатель развивает максимальное усилие при угле нагрузки 0 , потому как активное сопротивление много больше индуктивного, что противоречит пониманию идеализированной синхронной машины. При промышленной же частоте в 50 Гц максимальное усилие двигателя достигается при угле 90 , что соответствует идеализированной синхронной машине, в которой при расчётах активным сопротивлением обмоток пренебрегают.

Потерями в стали пренебрегаем в связи с низким диапазоном рабочих частот двигателя, также пренебрегаем механическими потерями. Таким образом, полезная мощность P2 Pэм , активная

70

Соседние файлы в папке книги