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

2025

.pdf
Скачиваний:
0
Добавлен:
15.11.2022
Размер:
1.09 Mб
Скачать
&
U1равна

21

стрелки. Режим работы двигательный. Для этого режима работы угол нагрузки , а также электромагнитный

момент считаем положительными. Так же считаем, что пусковая обмотка ротора приведена к обмотке статора известными методами . Это позволяет использовать одну и ту же систему относительных единиц для ротор-

ных и статорных величин .

В векторной диаграмме используются два угла, которые характеризуют нагрузку синхронных дви-

гателей.

Угол определяется положением вектора напряжения, приложенного к обмотке статора ,и попе-

речной осью q. Так как ротор двигателя вращается синхронно с полем прямой последовательности, принима-

&

 

ем за угол нагрузки угол между вектором напряжения прямой последовательности U1

и осью q ротора.

Угол характеризует отклонение продольной оси полюса ротора от оси МДС обмотки статора, обу-

словленное изменением нагрузки на валу электродвигателя. Так как к двигателю модели приложено напряже-

&

 

&

ние U1

, то угол определяет смещение оси ротора d относительно МДС F1, создаваемой током I1 .

На данной векторной диаграмме в отличие от

[9], E0 в двигательном режиме принята величиной

положительной, т.е. это есть составляющая напряжения, уравновешивающая ЭДС. Из векторной диаграммы составляющая напряжения

&

E0

j I1

cos X d e

j

j I1 sin X q e

j

/2

r

&

&

Z1

&

(2.34)

U1

 

 

 

I1

E0

I1 ,

 

 

Z1=r+j·Xd·cos e-j· +j·I1·Xq·sin ej·( /2- ) .

 

 

 

 

(2.35)

 

 

 

 

где: r - активное сопротивление фазы обмотки статора; Xd, Xq - индуктивные сопротивления фазы

обмотки статора по осям d и q . Учитывая, что j=ej· /2 и обозначив

 

 

 

 

 

 

 

Xac=1/2·(Xad+Xaq) ; Xu=1/2·(Xd-Xq) ,

 

 

 

 

(2.36)

 

 

 

 

 

где Xad, Xaq - сопротивления взаимоиндукции фазы обмотки статора по осям d и q , выражение для

сопротивления Z1

можно представить в виде

 

 

 

 

 

 

 

 

 

 

Z1=r+Xu·sin2 +j·(X +Xac+Xu·cos2 ) ,

 

 

 

 

(2.37)

 

 

 

 

 

где X - сопротивление рассеяния фазы обмотки статора. Определим соотношение между углами

и . Из векторной диаграммы можно написать:

 

 

 

 

 

 

 

 

 

 

U1·sin

=Xq·Iq-r·I1·cos

; U1·cos =Xd·Id+E0+r·I1·sin .

(2.38)

 

 

Из

(2.38) можно найти, что

 

 

 

 

 

 

 

 

 

 

Id

U1

cos

E0

Xq

U1 r

sin

;

 

 

 

r 2

Xd

Xq

 

 

 

 

 

 

 

 

Iq

U1

cos

E0

r

U1 X d

sin

.

(2.39)

 

 

r 2

X d

X q

 

 

 

 

 

 

 

Учитывая, что Id=I1·cos и Iq=I1·sin

, разделим второе уравнение на первое

 

22

tg

 

U1

cos

E 0

r

U1

Xd

sin

 

.

(2.40)

 

U1

cos

E 0

Xq U1

r

sin

 

 

 

 

 

 

В относительных единицах

 

 

 

 

 

 

 

 

tg

U1

cos

 

r

U1

Xd

 

sin

.

(2.41)

 

 

 

 

 

 

 

 

U1

cos

 

Xq

U1

r

sin

 

 

 

 

 

 

Сопротивление Z2 для токов обратной последовательности не остается постоянным вследствие яв-

нополюсности ротора. В первом приближении этими изменениями можно пренебречь, а сопротивление Z2

представить в виде

Z2

r j X

j X ac

ryc / S

j X yc

 

 

 

 

,

(2.42)

 

 

 

 

 

ryc / S

j X yc

X ac

 

где

ryc=1/2·(ryd + ryq) ; X yc=1/2·(X yd + X yq) ;

 

 

(2.43)

ryd , ryq - активные сопротивления роторной обмотки (пусковой) по осям d и q , приведенные к об-

мотке статора;

X yd , X yq - индуктивные сопротивления рассеяния роторной обмотки по осям d и q ;

S - скольжение.

В СДПМ обмотка возбуждения на роторе отсутствует, а постоянный магнит в модели данного дви-

гателя заменяется обмоткой возбуждения подключенной к источнику тока, вследствие чего он не влияет на сопротивление Z2.

2.5. Дифференциальные уравнения СДПМ

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

шать колебания около своей синхронной скорости. Двигатель, работающий в асинхронном тормозном режи-

ме, оказывает демпфирующее влияние на развитие колебаний.

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

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

Запишем дифференциальные уравнения равновесия напряжений для двигателя эквивалентной моде-

ли, работающего в синхронном режиме

Ua = d /dt + r·ia ; Ub = d /dt + r·ib , (2.44)

где Ua , Ub , ia , ib - мгновенные значения напряжений и токов прямой последовательности в фазах ;

a , b- полные потокосцепления фазовых обмоток . Уравнения равновесия напряжений для пуско-

вой обмотки ротора с учетом замены ее на две эквивалентные, сосредоточенные по продольной и поперечной

 

 

 

23

осям, будут иметь вид:

 

 

0 = d

yd/dt +

ryd·iyd ; 0 = d yq/dt + ryq·iyq ,

(2.45)

где

yd , yq

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

iyd , iyq - мгновенные значения токов в роторных обмотках. Потокосцепления обмоток в (2.44),(2.45) опреде-

ляются согласно , следующими выражениями:

a=La·ia+Mab·ib+Mayd·iyd+Mayq·iyq ;

b=Mab·ia+Lb·ib+Mbyd·iyd+Mbyq·iyq ; (2.46)

yd=Mayd·ia+Mbyd·ib+Lyd·iyd +Mad·Im0;

yq=Mayq·ia+Mbyq·ib+Lyq·iyq .

Как уже отмечалось выше постоянные магниты в модели конденсаторного СДПМ заменяются экви-

валентной фиктивной обмоткой возбуждения без потерь, подключенной к источнику тока [9], поэтому в

уравнение для потокосцепления yd вводится произведение тока возбуждения Im0 и взаимоиндукции Mad

фиктивной обмотки.

Уравнения (2.44), (2.45) являются дифференциальными уравнениями с периодическими коэффици-

ентами, так как вследствие явнополюсности ротора индуктивности La, Lb и взаимные индуктивности Mab, Mayd, Mayq,Mbyd ,Mbyq - периодически изменяются в зависимости от угла определяющего положение ротора относи-

тельно статора. Чтобы избавиться от периодических коэффициентов, проведем преобразование координат.

Заменим неподвижную двухфазную обмотку статора на эквивалентную, вращающуюся и жестко связанную с осями ротора d и q . В этом случае скорость вращения преобразованных обмоток статора относительно ротора равна нулю, а индуктивности и взаимоиндуктивности будут постоянными величинами. Таким образом, в пре-

образованной системе координат уравнения равновесия напряжений обмоток двухфазного СДПМ, питаемого системой напряжений прямой последовательности в относительных единицах, будут иметь вид:

U1d=d

1d/d

-

1q· +r·i1d

;

 

U1q=d

1q/d

-

1d· +r·i1q

;

(2.47)

0=d

1yd/d +ryd·i1yd ;

 

 

0=d

1yq/d +ryq·i1yq ,

 

 

где U1d , U1q

- напряжения прямой последовательности в системе координат d,q ;

- угловая частота вращения ротора;

 

i1d ,

i1q , 1d

,

1q - соответственно тока и потокосцепления обмоток статора в новой системе коор-

динат.

Индекс "1" определяет переменные и параметры, обсловленные напряжением прямой последова-

тельности U1 .

В системе относительных единиц, индуктивности и взаимоиндуктивности численно равны индук-

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

L=L/L =L· /(L · )=X/Z =X . (2.48)

Произведение тока возбуждения Im0 на сопротивление взаимоиндкции Xad в системе относительных единиц можно заменить на степень возбужденности постоянного магнита [9].

Запишем выражения для потокосцеплений с учетом выше изложенного

1d=(Xad+X )·i1d+Xad·i1yd + ;

 

24

1q=(Xaq+X )·i1q+Xaq·i1yq ;

 

1yd=Xad·i1d+(Xad+X yd)·i1yd + ;

(2.49)

1yq=Xaq·i1q+(Xaq+X yq)·i1yq .

 

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

U1d= -U1·sin

; U1q=U1·cos

,

(2.50)

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

статора.

 

 

 

 

 

Угол нагрузки

 

и частота вращения ротора связаны соотношением

t

 

 

 

 

 

1

 

dt

0

,

(2.51)

0

 

 

 

 

 

где 1 = б - угловая скорость вращения магнитного поля;

0 - начальное значение угла нагрузки.

Тогда продифференцировав уравнение (2.51) получим

 

d /dt=

1- .

 

 

(2.52)

Чтобы выразить уравнение (2.52) в системе относительных единиц, разделим обе его части на ба-

зисную угловую скорость

б

 

 

 

d /d

бt= 1/ б- /

б .

 

(2.53)

Учитывая выражение (2.53) и f1= fН в системе относительных единиц (2.52) запишется в виде

d /d =1- .

 

 

(2.54)

При частотном управлении конденсаторным СДПМ, учитывая, что fб = fН , б=2· ·fН, =f/fН , уравне-

ние (2.54) будет иметь другой вид, чтобы его получить выразим из уравнения (2.52) угловую скорость ротора

 

= 1- d /dt .

(2.55)

Разделим (2.55) почленно на базисную угловую скорость

 

/ б=

1/ б- d /d бt=2· ·f/(2· ·fб)- d /d = - d /d .

(2.56)

Окончательно получим

 

d

/d = - .

(2.57)

Нарушение статической устойчивости в синхронных электродвигателях сопровождается колеба-

ниями скорости вращения ротора. Поэтому, чтобы учесть данное явление в конденсаторных СДПМ, необхо-

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

Исходное уравнение движения ротора двигателя, работающего от напряжения прямой последова-

тельности (уравнение моментов на валу), можно записать в виде

 

J/p· d /dt+Mс-M1=0 ,

(2.58)

где J - момент инерции вращающихся масс ротора и нагрузки;

 

p - число пар полюсов;

 

Mc - момент сопротивления на валу;

M1 - электромагнитный момент, создаваемый двигателем, работающим в синхронном режиме.

Вэквивалентной модели конденсаторного СДПМ механически соединены два двигателя. Поэтому

вуравнении (2.58) нужно учесть тормозной момент M2, создаваемый двигателем, работающим от напряжения обратной последтельности, т.е. момент сопротивления на валу будет равен M2+Mc .

Уравнение движения ротора эквивалентной модели с учетом сказанного, можно записать в виде

 

 

25

 

M1=M2+Mс+ J/p· d /dt

(2.59)

или

M1=M2+Mс+ J· d /dt ,

(2.60)

где

= /p - угловая скорость ротора.

 

Момент двигателя, работающего в синхронном режиме M1 , определяется изестным выражением

 

M1= 1d·i1q- 1q·i1d ,

(2.61)

где токи и потокосцепления определяются из уравнений(2.47)-(2.49).

Момент Mc включает в себя момент холостого хода двигателя, учитывающий трение, момент не-

уравновешенных масс и момент внешних сил нагрузки . Особое влияние на затухание колебаний скорости вращения ротора могут оказывать силы трения. Учет сил трения, способствующих затуханию колебаний,

производится введением в уравнение движения ротора коэффициента механического демпфирования . Мо-

мент сопротивления, учитывающий демпфирующее влияние сил трения при малых колебаниях скорости вра-

щения ротора, запишется в виде

 

Mс=Mс0+km1· ,

(2.62)

где Mс0 - составляющая момента сопротивления в установившемся синхронном режиме, включая потери холостого хода и полезный момент на валу двигателя;

km1 - коэффициент механического демпфирования;

= - 1 - изменения угловой скорости вращения ротора относительно синхронной.

Двигатель обратной последовательности создает тормозной момент M2. Величина данного момента в заданном синхронном режиме работы определяется из следующих рассуждений. Ротор двигателя обратной

последовательности по отношению к полю обратной последовательности вращается со скольжением

S2=(- 1- )/(- 1)=2-S1 ,

(2.63)

где S1 - скольжение ротора относительно поля прямой последовательности.

Так как в установившемся синхронном режиме ротор двигателя вращается синхронно с полем пря-

мой последовательности (S1=0), то из (2.63) получаем S2=2.

Зная напряжение обратной последовательности и параметры двигателя, можно рассчитать асин-

хронный тормозной момент в установившемся режиме М20 по методике, изложенной в [8].

Демпфирующее влияние двигателя обратной последовательности на развитие колебаний ротора оп-

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

а может вызывать рост колебаний - отрицательное демпфирование. Знак демпфирования зависит от парамет-

ров роторной обмотки, а учет влияния момента М2 на развитие колебаний производится с помощью коэффи-

циента электромагнитного демпфирования km2 . Коэффициент km2 определяется из зависимости M2=f(S2) , как тангенс угла наклона для скольжений S2=1,95 и 2,05 , что соответствует следующему выражению

k m2

M 2 S2 2.05

M 2 S2

195.

.

(2.64)

 

01.

 

 

 

 

 

 

Выражение, определяющее электромагнитный демпфирующий момент имеет вид

M2= km2· . (2.65)

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

26

Момент от постоянного магнита для двигателя обратной последовательности в установившемся ре-

жиме определяется согласно [9], следующим выражением

M 30

 

m

I32

 

r

 

,

 

 

 

(2.66)

 

 

 

 

S

 

 

 

 

 

1

 

 

 

 

 

 

 

 

где m - число фаз двигателя; I3 - ток, определяемый влиянием постоянных магнитов;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 2

1

S 2

X aq

X

2

 

 

 

I3

1 S

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

(2.67)

r 2

1 S 2

 

 

 

 

 

 

 

 

X ad

X

X aq X

 

Демпфирующее влияние на развитие колебаний скорости вращения ротора от момента M3 можно учесть с помощью коэффициента демпфирования km3. Выражение для демпфирующего момента от постоян-

ных магнитов можно записать в виде

M3 = km3· . (2.68)

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

M3=f(S2) , для скольжения S2=1,05 и 0,95. Таким образом выражение для km3 запишется следующим образом

k m3

 

M 3 S2

 

105.

M 3 S2

0.95

.

 

 

 

(2.69)

 

 

 

 

 

 

 

01.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определив сумму коэффициентов демпфирования, через результрующий коэффициент

km= km1+km2+km3 ,

 

 

 

 

 

 

 

 

(2.70)

 

уравнение движения ротора (2.60) с учетом (2.61)-(2.69) запишется в виде

 

M1=Mс0+M20+km·( - 1)+J·d /dt .

 

 

 

 

(2.71)

 

Выражение (2.71) в относительных единицах примет вид

 

 

 

 

M1=Mс0+M20+km·(

-

)+H·d

/d

,

 

 

 

 

(2.72)

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

где H

 

J

 

н

 

- механическая постоянная вращающихся масс. Результирующий коэффици-

 

Pб

p

 

ент демпфирования в относительных единицах имеет вид

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

k m

k m

 

 

б

 

k m1

k m2

k m3

 

 

б

 

.

(2.73)

 

P

 

 

 

P

 

 

 

 

 

 

p

 

 

 

 

 

p

 

 

 

 

 

б

 

 

 

 

 

 

 

б

 

 

 

 

Размерность коэффициента демпфирования в физических единицах системы СИ - Н·м·сек/рад.

Таким образом, уравнения (2.47) с учетом(2.50), а также (2.57) и (2.72) составляют дифференциаль-

ных уравнений эквивалентной модели конденсаторного СДПМ. Запишем данную систему уравнений в не-

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

пользуемых при расчетах на ЭВМ. Для этого перенесем в левые части уравнений первые производные пере-

менных, то есть представим их в виде нормальных дифференциалных уравнений. С учетом изложенного в системе относительных единиц уравнения примут вид

d

1d/d =

1q· -r·i1d-U1m·sin

;

d

1q/d =-

1d· -r·i1q+U1m·cos

;

 

 

27

d

1yd/d =-ryd·iyd ;

 

d

1yq/d =-ryq·iyq ;

(2.74)

d

/d =(1/H)·(M1-Mс0-M20-km·( - )) ;

 

d

/d = - .

 

Определив из выражения (2.49) токи iyd, iyq , i1d , i1q и подставив их в (2.74), получим систему обык-

новенных дифференциальных уравнений первого порядка. Зависимыми переменными в данной системе будут

потокосцепления

1d , 1q , 1yd , 1yq , угловая частота ротора и угол нагрузки . Время - является неза-

висимой переменной.

Уравнения системы (2.74) нелинейные, так как они содержат произведения переменных и нелиней-

ные функции sin

и cos .

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

совместно с (2.49) и (2.61), составляет математическую модель конденсаторного СДПМ при работе его от источника с регулируемой частотой и напряжением.

3. СИНХРОННЫЙ УСТАНОВИВШИЙСЯ РЕЖИМ РАБОТЫ

3.1. Расчет напряжений прямой и обратной последовательностей

В предыдущей главе были получены уравнения для напряжений пря-

мой и обратной последовательностей. Используя данные соотношения, разра-

ботаем методику расчета напряжений U1 и U2 . Расчеты будем проводить в относительных единицах. Зная параметры исследуемого двигателя, подводи-

мое напряжение, возбужденность постоянного магнита, напряжения прямой и

обратной последовательностей определяются в следующем порядке:

1)

задается углом

;

2)

находим угол

из равенства (2.41);

3)

рассчитываем сопротивление токам прямой последовательности Z1

из уравнения (2.37);

 

4)

рассчитываем сопротивление токам обратной последовательности

Z1 из уравнения (2.42);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

28

 

 

 

5) находим токи прямой I1

и обратной I2 последовательностей соответ-

ственно из уравнений (2.31) и (2.32);

 

 

 

 

6) определяем U1 ,U2

и угол 1

из уравнений (2.33).

 

 

Для определения напряжений U1 и U2, необходимо знать U1 и

U, где

U=

1+ , а также Z C . Сопротивление Z C

определяется из равенства (2.14).

Полное сопротивление конденсатора

Z C входящее в (2.14), рассчитывается

по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZC

j XC /

 

 

R C .

 

(3.1)

 

 

Z1 - находится из следующих рассуждений. В соответствии с [7]

можно записать

 

 

 

 

 

 

 

 

 

 

 

Z B

 

3

r

X au

sin 2

j

3

 

X

X ac

X au cos2 .

(3.2)

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда из (2.16) с учетом (2.37)

 

 

 

 

 

 

Z1

Z B

Z1

1

 

r

j

X .

(3.3)

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Величины

U и U1 неизвестны, поэтому для расчета U1 и U2 использу-

ем метод последовательных приближений. На первом этапе расчета принима-

ем U1 =U=

; U =

;

1 = 0 . Определяем из пункта 6 U1 , U2 и угол

1 .

Найденные значения

1

и U1 подставляем в расчетные формулы и снова рас-

считываем 1

и U1 . Расчет повторяется до тех пор, пока найденные 1

и U1

будут отличаться от предыдущих значений не больше, чем на 1 %. Затем сно-

ва задается значение , и расчеты повторяются.

Изложенная выше методика позволяет рассчитывать напряжения пря-

мой и обратной последовательностей также для моделей конденсаторных и однофазных синхронных реактивных двигателей (СРД) [7]. Расчет U1 и U2 для модели конденсаторного СРД проводится с подстановкой = 0 в формулах

(2.22),(2.23) и (2.25), а для однофазного СРД еще необходимо принять равным нулю емкость конденсатора С в (3.1).

29

Данная методика расчета напряжений реализована в виде программы для персональных ЭВМ, на языке ФОРТРАН. Схема алгоритма приведена на рис.3.1.Текст программы приведен в приложении 2.

Характерной особенностью СДПМ с радиально расположенными маг-

нитом и пусковым устройством является меньшее значение индуктивного со-

противления обмотки статора по продольной оси, чем по поперечной Xad

Xaq .Это объясняется тем, что в синхронных машинах данной конструкции значительная часть продольного потока якоря из-за малой магнитной прони-

цаемости магнита замыкается по путям рассеяния полюсных башмаков, маг-

нитное сопротивление которых значительно превышает магнитное сопротив-

ление воздушного зазора.

Начало

1. Ввод r, XS,Xad

Xaq, Xd, Xq, ryd, ryq, X yq, X yd, ,

S, k, kZC1,

2.= -90

=· /180

A

3. i = 1, 36, 1

4. 1 =

 

U =

Конец

30

U =

5. Расчет , Xac, Xu, Z1, ryc, X yc,

A, B, C, Z2, Zc,

B

Рис.3.1. Структурная схема расчета напряжений прямой и обратной последовательностей

B

6.Расчет A, B, I1, I2, U1, U2, A1, B1,

,U

7.Расчет A1, B1, C1, D1

 

НЕТ

8. A1 < B1

ДА

 

 

 

 

НЕТ

 

 

 

 

9. C1 < D1

 

 

 

 

 

 

10. U11 = U1

 

 

ДА

2 = 1

 

 

 

 

 

 

 

 

11. Печать U1,

 

 

 

 

U2, ,

 

 

 

 

 

 

 

 

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]