Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
495.pdf
Скачиваний:
68
Добавлен:
16.02.2016
Размер:
959.38 Кб
Скачать

1.

МАГНИТНЫЙ ШУМ И ВИБРАЦИЯ АСИНХРОННЫХ

 

 

 

ЭЛЕКТРОДВИГАТЕЛЕЙ

 

 

 

 

1.1. Общие положения

 

 

Вибрация электрических машин зависит от характера сил, вызывающих

её, и мест их приложения, поэтому динамическая модель машины может

быть представлена в виде поступательной (рис. 1) или вращательной колеба-

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

которых производится методами прикладной теории колебаний [2].

 

 

 

В механической поступательной системе движущая сила

 

Pм

Рм действует на массу m, прикреплённую к упругому элемен-

 

m

ту λм, с потерями на трение rм, пропорциональными скорости,

а в механической вращательной системе движущий крутящий

λм

момент Mθ действует на маховик с моментом инерции J, сое-

 

динённым с вращательным упругим элементом λθ, с потерями

 

на трение rθ, пропорциональными скорости.

 

 

Рис. 1

Уравнения сил и моментов для этих систем имеют вид

 

m d 2 y + r

dy

+ y = P ejωt ,

 

 

 

(7)

 

dt2

м

dt

λм

 

м

 

 

 

J d 2ϕ

+ r

dϕ

+

ϕ = M

 

ejωt ,

(8)

 

dt2

и

dt

 

 

λи

и

 

 

 

d 2 y

 

 

dy

 

 

 

 

jωt

 

где m dt2

- сила инерции; rм

dt

- сила трения; ум – сила упругости; Рме

приложенная к системе сила;

J d 2ϕ - момент сил инерции вращающегося ма-

 

 

 

 

 

dt2

 

 

 

 

ховика; r

dϕ - момент сил трения; φ/λ

 

– момент сил упругости;

М еjωt

и

dt

 

 

 

 

 

θ

 

θ

 

внешний момент, приложенный к системе.

 

 

Решения этих уравнений имеют вид

 

 

 

 

y =

 

 

P ejωt

 

 

 

 

 

 

P

 

 

 

 

м

 

 

 

 

 

 

=

 

 

м ,

(9)

 

r

+ jωm - j/(ωλ

 

 

)

z

 

 

м

 

м

 

 

 

 

м

 

 

 

 

 

 

 

 

 

 

ϕ =

 

 

M

ejωt

 

 

 

 

 

=

m

 

,

(10)

 

 

и

 

 

 

 

 

 

 

и

r

+ jωJ - j/ (ωλ

 

 

)

 

z

 

м

 

 

и

 

 

 

 

и

 

 

 

 

 

 

 

 

 

где zм и zθ – полные комплексные сопротивления.

33

Единица сопротивления для механической поступательной системы – ньютон-секунда на сантиметр, а для механической вращательной системы – ньютон-сантиметр-секунда на радиан.

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

b(θ, t) = f(θ, t)Λ(θ, t),

(11)

где θ – координата рассматриваемой точки воздушного зазора в момент времени t; f(θ, t) и Λ(θ, t) – мгновенные значения результирующей МДС обмоток статора и ротора и магнитной проводимости зазора.

Значение и распределение радиальных магнитных сил (в ньютонах на квадратный сантиметр) в воздушном зазоре определяется уравнением

pr = 40b2(θ, t),

 

(12)

а тангенциальных сил

 

 

pτ = a(θ, t)b(θ, t)·10-2,

(13)

где a(θ, t) – линейная токовая нагрузка статора или ротора, А/см.

 

1.2. МДС обмоток статора и ротора

 

 

МДС при её разложении в ряд Фурье имеет вид

 

н=∞

м=∞

 

f (θ, t) = f0(θ, t) + fн(и,t) +

fм(и,t) ,

(14)

н=1

м=1

 

где f0(θ, t) – результирующая МДС основной волны; fν(θ, t) – МДС ν-й гармоники обмотки статора; fµ(θ, t) – МДС µ-й гармоники обмотки ротора.

В уравнении (14)

f0(θ, t) = F0cos(pθ – ω1t – φ0r),

(15)

fν(θ, t) = Fνcos(νθ – ω1t – φ1),

(16)

fµ(θ, t) = Fµcos(µθ - ωµt – φ2).

(17)

34

Амплитуды МДС для трёхфазных обмоток равны

F =1,35

 

w1kw1

 

I

0r

,

(18)

 

 

0

 

p

 

 

 

 

 

 

 

 

 

 

 

 

F =1,35

w1kwн

I

,

 

(19)

 

 

 

н

 

p

 

 

1

 

 

 

 

 

 

 

 

 

 

 

Fм =1,35

 

w1kw1

 

I2.

(20)

µ

 

где р, ν, µ – числа пар полюсов основного поля, гармоник МДС обмотки статора и короткозамкнутого ротора; ω1, ωµ – угловые скорости основной волны и µ-й гармоники ротора относительно статора; w1 – число витков фазы обмотки статора; kw1, k– обмоточные коэффициенты для основной волны и ν-й гармоники обмотки статора; I0r, I1, I2, φ0r, φ1, φ2 – и их фазовые сдвиги отно-

сительно напряжения.

Число ν для обмотки с целым числом пазов на полюс и фазу q равно

ν = (6q′+1)p при q′= ±1; ± 2; ± 3.

(21)

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

q′= ±q ; ± 2q; ± 3q; … .

(22)

При q> 0 гармоники вращаются в направлении главного поля (ν = 7р, 13р, 19р, …), а при q< 0 – против вращения главного поля (ν = -5р, -11р, -17р, …). Угловая скорость ν-й гармоники относительно статора

ων = ω1/(6q´ + 1).

(23)

Число µ для короткозамкнутого ротора

µ = q˝Z2 + p; q˝ = ± 1; ± 2; ± 3; … ,

(24)

где Z2 – число пазов ротора. Угловая скорость µ-й гармоники

ωµ = ω1[1 + q˝(Z2/p)(1 – s)],

(25)

где s – скольжение.

35

1.3. Магнитная проводимость и магнитное поле воздушного зазора

Магнитная проводимость воздушного зазора в общем виде может быть представлена формулой

Λ(θ, t)= Λ0 +λm1 + λm2 + ,

(26)

m1

m2

 

где Λ0 = 1/(δkc) - постоянная составляющая магнитной проводимости воздушного зазора; λm1 и λm2 – гармоники проводимости зубчатости статора при гладком роторе и зубчатости ротора при гладком статоре; δ – величина воздушного зазора; kc – коэффициент Картера [1], учитывающий зубчатость ста-

тора kc1 и ротора kc2.

Гармоники зубчатости статора постоянны во времени

 

 

 

λm1 = Λm1cosk1Z1θ при k1 = 1, 2, 3, … ,

 

(27)

а гармоники зубчатости ротора при его вращении зависят от времени

 

λm2 = Λm2cosk2Z2[θ – (ω1/p)(1 – s)t]

при k2 = 1, 2, 3, … .

(28)

Амплитудные значения Λm1 и Λm2 определяются из уравнений

 

Л

m1

= −Л

0

(1)k1 (k

c1

1)

sink1 [(kc1 1)/kc1 ]π

= Λ ξ ,

(29)

 

 

 

 

 

 

k1 [(kc1

1)/ kc1 ]/ π

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Лm2

= −Л0 (

1)k2 (kc2 1) sink2 [(kc2

1) / kc2 ]π

= Λ0ξ2 .

(30)

 

 

 

 

 

 

 

 

 

k2 [(kc2 1) / kc2 ]/ π

 

 

 

Значения функций ξ1 и ξ2 выбирают в зависимости от коэффициентов kc1

или kc2 и порядка гармоник k1 или k2 (см. рис. 2 для k=1, k=2 и k=3).

ξ12

0.35

1

0.29

0.23

 

0.17

 

0.12

 

0.058

2

0

 

0.058

 

0.12

 

0.17

3

0.23

0.29

 

0.350.9 0.961.021.081.141.2 1.261.321.381.441.5

kc1,kc2

Рис.2

36

Подставив, все составляющие разложения для МДС и проводимостей в уравнение (11) получим результирующее магнитное поле в воздушном зазоре

b = b1 + bн + bм = B1 cos( pθ- ω1t - ϕ0r ) + Bн cos(νθ−ω1t - ϕнr ) +

н

м

н

(31)

 

 

+ Bм cos(µθ−ωмt - ϕмr ),

 

 

 

м

которое состоит из основного поля b1 и гармоник статора bν и ротора bµ , имеющие зубцовые гармоники bz1 и bz2.

Для основного поля

b1= B1cos(pθ – ω1t – φ0r),

(32)

где B1 = F0/kнΛ0 Bδ – амплитудное значение; kн – коэффициент насыщения магнитной цепи, равный отношению общей МДС к МДС воздушного зазора для основного поля; Bδ – индукция в воздушном зазоре в Тл.

Гармоники обмоток статора и короткозамкнутого ротора определяются уравнениями

bν = fνΛ0 = Bνcos(νθ – ω1t – φ1),

(33)

bµ = fµΛ0 = Bµcos(µθ - ωµt – φ2),

(34)

а их амплитудные значения по формулам

 

 

 

 

 

 

B

 

= F Λ

 

=

p

kk

 

 

 

I1

 

 

B ,

(35)

 

 

 

 

 

 

I0r

 

н

 

н

 

0

 

 

ν kw1

 

 

н

д

 

B

 

= F

Λ

 

=

p

 

1

 

I2

k

 

B .

(36)

 

 

µ kw1

 

 

 

м

м

 

 

0

 

 

 

 

I1

 

 

н

д

 

Зубцовые гармоники статора и ротора определяются уравнениями

 

bz1 = Bz1cos(νzθ – ω1t – φ0r),

(37)

bz2 = (f0/kнm2 = Bz2cos(µzθ - ωµt – φ0r),

(38)

а их амплитудные значения по формулам

 

 

 

 

 

 

Bz1 = (F0/kнm1 = (F0/kн0ξ1 = Bδξ1,

(39)

 

Bz2 = (F0/kнm2 = Bδξ2.

(40)

37