CC BY
44
УДК 621.313 Ю.В. Шуруб
МОДЕЛИРОВАНИЕ АСИНХРОННЫХ ЭЛЕКТРОПРИВОДОВ ПРИ СЛУЧАЙНОЙ НАГРУЗКЕ
Наведено методику моделювання та вибору параметрів асинхронних електроприводів, що працюють у несиметричному режимі, при роботі із навантаженнями, які мають випадковий характер.
Приведена методика моделирования и выбора параметров асинхронных электроприводов, работающих в несимметричном режиме, при работе с нагрузкой, которая имеет случайный характер.
ВВЕДЕНИЕ Асинхронный электропривод значительной части сельскохозяйственных, строительно-отделочных машин и механизмов, имеющих случайную нагрузку, должен рассчитываться с учетом требуемых динамических и энергетических показателей [1, 2]. Во многих случаях питание трехфазного двигателя малой мощности осуществляется от однофазной сети -трехфазно-однофазный электропривод, а в качестве фазосдвигающего элемента используется конденсатор, реже полупроводниковые элементы. При этом необходимо учитывать следующие обстоятельства.
В настоящее время проектируются и выпускаются высокоиспользуемые двигатели с меньшим запасом по нагреву и перегрузочной способности. В то же время качество напряжения питающей сети часто не соответствует норме. В общем случае на входе двигателя имеется несимметричная система напряжений. При случайном характере нагрузки двигателя и напряжения сети выбор и проверка параметров двигателя и фазосдвигающих элементов представляют сложную задачу.
ПОСТАНОВКА ЗАДАЧИ Известно, что симметрия фазных токов асинхронного двигателя с постоянными параметрами фазосдвигающих элементов возможна при одной какой-либо нагрузке. Математически такое событие - попадание в точку трехмерного пространства - имеет вероятность, равную нулю [3]. В практических расчетах режим считают симметричным, если фазные токи отличаются на определенную величину.
Вероятность появления такого режима, когда значения фазных токов заключены в пределах
\IA1 IA2\ \IB1 IB2\ интеграла
Vc Ic
C1 C2
получим в виде тройного
f (IA, Ib , Ic ) = -
1
ОI ОI О t
1A 1B 1C
Vsn^D
1
xexp[-2D ^ Dv
yv=1
(2)
0y0v
где тр шу, - математическое ожидание и дис-
персия случайной величины 1^, /у соответственно; Б -определитель третьего порядка, составленный из значений нормированных коэффициентов корреляции при т = 0:
R|Xv Rvy
Kyy (0) 0y0v
; D =
1 rab rac
’BA 1 RBC
CA RCB 1
Для анализа трехфазно-однофазных электроприводов могут быть использованы математические модели случайных процессов изменения нагрузки двигателя при симметричном номинальном напряжении.
В режиме нагрузки линейный ток ДО, потребляемая мощность Р\(Ґ) и момент М(Ґ) являются функциями времени. Если, например, статистические характеристики функции І(ґ) известны, т.е. известна плотность вероятности ДІ), математическое ожидание тІ и корреляционная функция ЯІ(ї-і,ї2), то могут быть определены соответствующее среднеквадратическое отклонение и среднеквадратическое значение тока нагрузки [3].
В случае, когда приложенная к двигателю нагрузка подчинена нормальному закону распределения и стационарна, т.е.
1
(і-т )2 2о]
(3)
■А 1в2 1с2
Р = 111/ (/А,,!с ^а^в^с . (1)
1л1 -Ц ^С\
Нормальное совместное распределение трех случайных величин 1А, 1В, 1С характеризуется плотностью вероятности
f(I) = ■ ^
О jV 2п
mj (t ) = mj = const;
о2 (t) = Rj , 12) = Rt (0) = const, математическое ожидание и дисперсия квадрата тока:
2 2
mj2 = mj +Оj ;
= 2о2(2да| +О2) .
(4)
(5)
x
МОДЕЛИРОВАНИЕ ПРОЦЕСА НАГРУЖЕНИЯ ДВИГАТЕЛЯ В НЕСИММЕТРИЧНОМ РЕЖИМЕ РАБОТЫ
Наличие сведений о статистических характеристиках тока, мощности или момента позволяет распространить метод эквивалентных величин на симметричные составляющие тока и напряжения прямой, обратной и нулевой последовательностей. Эти составляющие действуют независимо и вызываемые ими эффекты могут быть суммированы по методу суперпозиции. Если изменение нагрузки характеризуется как случайный стационарный процесс, то процесс изменения величин симметричных составляющих описывается стационарными случайными функциями. Таким образом, задача сводится к определению числовых характеристик случайных величин, имеющих функциональную зависимость от другой случайной величины.
При нестационарности процесса нагружения можно обычно выделить достаточно продолжительные зоны стационарности, т.е. построить кусочностационарную вероятностную модель процесса и на этой основе пользоваться для расчета методами теории стационарных случайных процессов.
Величина математического ожидания квадрата тока в такой же мере является объективной оценкой условий нагрева, как и среднеквадратическое значение тока для электродвигателей с детерминированным режимом работы. Правомерность эквиваленти-рования нагрузки на этой основе базируется на слабой коррелированности процесса изменения нагрузки, т. е. малой величине времени т = /2 — /ь при которой автокорреляционная связь является существенной по сравнению с постоянной нагрева двигателя [4].
Обычно при исследовании процессов изменения нагрузки располагают зависимостью потребляемого тока Щ) или мощности Р\((). При неизменном напряжении питающей сети и постоянных параметрах фазосдвигающих элементов величины токов и напряжений прямой и обратной последовательностей определяются однозначно через параметры двигателя при симметричном питании. Зная закон распределения случайной величины I, можно определить числовые характеристики случайных величин 12 = ф2(1), и = фз(1). Математическое ожидание и дисперсия, например, тока прямой последовательности, определяется выражениями
щх = 1 Фх(1 У/ (I М; (6)
—
012 = 1 [(1) — шк ]/(I) • Ш, (7)
—
где фа(1) - зависимость составляющей тока прямой последовательности от реального тока в фазе (симметричного режима); /(I) - плотность распределения величины !(/).
Отметим, что плотность распределения нагрузки прямоточной дробилки зерна хорошо аппроксимируется законом Гаусса [2].
Если функция I = ф1(!) достаточно близка к линейной в области практически возможных значений случайных нагрузок, то приближенные значения оценок математического ожидания и дисперсии выражаются элементарными формулами.
Для линейной функции вида I = Ь + а1, подчиненной нормальному закону, математическое ожидание и среднеквадратическое отклонение определяются выражениями
= а • т+ Ь ; (8)
2
012 = а2 •О 1. (9)
Если же эти функции значительно отклоняются от линейных, то формулы (3) и (4) применимы только при малых дисперсиях случайных нагрузок.
После определения допустимой величины тока прямой последовательности производим сопоставления его со значением эквивалентной величины тока прямой последовательности. При этом должно быть соблюдено условие:
11Э =^т112 +012 - 1 доп , (10)
т.е. эквивалентный ток прямой последовательности
должен быть меньше или равен допустимому. Если окажется, что условие (5) не выполняется, то необходимо изменить значение расчетной нагрузки.
С учетом сказанного, расчетная величина нагрузки для выбора параметров фазосдвигающих элементов может быть представлена как
1Р = т1 + Р0 I. (11)
Значение коэффициента в зависит от "растянутости" закона распределения тока нагрузки по оси величин и установленной мощности двигателя. Если плотность вероятности /(I) известна, то вероятность превышения ординатой функции/(I) уровня mI + р^
Е = 1 /(I)еИ . (12)
mI+во
Из (12) может быть с необходимой точностью определено значение коэффициента р. Для нормально распределенной случайной величины тока все рассеивание (с точностью до долей процента) укладывается на участке mI + 3 •
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ
По предложенному методу проведен выбор параметров фазосдвигающих элементов и двигателя прямоточной дробилки зерна, проверена работоспособность электропривода при трехфазном и однофазном питании. Определены оценки процесса нагружения как в статике (математическое ожидание, дисперсия), так и в динамике (корреляционная функция, спектральная плотность). Для установления закона распределения случайных величин тока и мощности построены соответствующие гистограммы. Проверка при помощи критерия Пирсона на уровне значимости
0,05 показала, что плотность распределения нагрузки прямоточной дробилки зерна хорошо аппроксимиру-
ется нормальным законом. Коэффициент вариации нагрузки (отношение среднеквадратического отклонения к математическому ожиданию) для исследуемых диаграмм находится в пределах 0,17 до 0,36.
Корреляционная функция и спектральная плотность рассматриваемых диаграмм хорошо аппроксимируются выражениями вида
п/ ч 2 -ат о
R(t) = а e 11 cosрт.
( ) ° 5(Ю) =-----
п
а а
(13)
(14)
а 2 + (ю+р)2 а 2 + (со—в)2 где а - коэффициент, характеризующий интенсивность затухания корреляционной функции; в - коэффициент, определяющий угловую частоту колебаний.
В этом смысле коэффициенты а и в характеризуют технологический процесс, выполняемый машиной, а форма и вид корреляционной функции - стационарность и эргодичность процесса нагружения машины. Спектральный состав процесса нагружения прямоточной дробилки зерна почти равномерен.
Результаты обработки нагрузочных диаграмм подтверждают гипотезу, что при постоянной производительности нагрузку исследуемых машин можно считать стационарной случайной функцией времени, обладающей эргодическими свойствами.
Математическую модель процесса нагружения однотипных машин с переменной производительностью можно представить как стационарную неэргоди-ческую функцию, содержащую случайную составляющую, обусловленную производительностью машины, и случайную стационарную эргодическую функцию, обусловленную технологией процесса:
г (о = у+*(о, (15)
где у - случайная величина с дисперсией су; х(ґ) - эр-годическая стационарная случайная функция с характеристиками тх и Кх.
Приняв, что x(t) и у не коррелированны, можно найти математическое ожидание и корреляционную функцию процесса нагружения машины. В дальнейшем по степени разброса нагрузки и по характеру затухания корреляционной функции производится корректировка параметров фазосдвигающих элементов в зависимости от изменения производительности машины.
СПИСОК ЛИТЕРАТУРЫ
1. Корчемний М.О., Філоненко А.Ф., Юсупов Н.А. Експлуатаційна надійність електродвигунів у сільськогосподарському виробництві // Механізація та електрифікація сільського господарства. - 1991. - Вип. 73. - С. 78-80.
2. Беспалов В.Я. Динамические показатели трехфазных асинхронных двигателей в однофазной сети // Электротехника. - 2000. - № 7. - С. 13-15.
3. Вентцель Е.С., Овчаров Л.А. Прикладные задачи теории вероятностей. - М.: Радио и связь, 1983.
4. Заславская Т.Б., Железко Ю.С., Нейман В.В. О допустимых значениях токов и напряжений обратной последовательности при случайном характере их изменения // Электричество. - 1981. - № 9. - C. 58-60.
Поступила 01.10.2009
Шуруб Юрий Викторович, к.т.н.
Институт электродинамики НАН Украины Украина, 03680, Киев, пр. Победы, 56 тел. (044) 454-26-37
Yu.V. Shurub
Simulation of induction electric drives under random load.
A simulation and parameter choice technique for phase shifting elements of induction electric drives in an asymmetrical operating mode under random loads is given.
Key words - induction motor, electric drive, asymmetrical mode, random load.
