с 01.01.2023 по настоящее время
Москва, г. Москва и Московская область, Россия
Москва, г. Москва и Московская область, Россия
Статья посвящена разработке программы расчета параметров Т-образной схемы замещения трехфазных асинхронных электродвигателей. Актуальность данной работы обусловлена необходимостью получения достоверных результатов компьютерного имитационного моделирования систем электропривода в программных комплексах. Целью исследования является реализация в удобном программном интерфейсе аналитической методики расчета параметров схемы замещения асинхронного двигателя с минимальными отклонениями от каталожных данных. Для программного расчета параметров схемы замещения рассматриваются аналитические аппроксимационные и итерационные методики счета, так как они не требуют проведения экспериментов на промышленном электрооборудовании. Выполнен сравнительный анализ таких методик на примере трёхфазного асинхронного двигателя марки ВА280S4 мощностью 110 кВт. Результаты сравнительного анализа показывают, что наибольшая достоверность механической и электромеханической характеристики достигается в случае применения итерационной методики счета с минимизацией целевой функции суммы взвешенных квадратов отклонений значений контрольных точек от каталожных данных. Программа расчета параметров схемы замещения по выбранной методике счета реализована в среде Embarcadero Delphi. Рас-считываемые программой достоверные значения параметров электродвигателей позволят проводить максимально соответствующее реальным объектам компьютерное моделирование систем электропривода.
асинхронный электродвигатель, схема замещения, расчет параметров, программа расчета, итерационная методика, многокритериальный подбор, целевая функция, механическая характеристика, моделирование
Введение
В настоящее время наиболее доступным и дешевым методом синтеза и исследования статических и динамических режимов работы электромеханических систем является компьютерное имитационное моделирование в программных комплексах, как например: ETAP, NI Multisim, MATLAB/Simulink, SimInTech. В связи с тем, что 80 % всей потребляемой в промышленности электрической энергии приходится на асинхронный электропривод [1], возникает актуальная задача наиболее точного приближения параметров математической и компьютерной модели трехфазного асинхронного двигателя (АД) к реальному оборудованию в целях повышения достоверности результатов моделирования.
Производители АД предоставляют информацию о технических характеристиках электрических машин в каталогах выпускаемого оборудования. Однако каталожные данные не содержат соответствующих параметров эквивалентной схемы замещения АД, составление которой требуется для расчета электромеханических переходных процессов [2]. Определение параметров схемы замещения трехфазного асинхронного двигателя традиционно производится путем стандартных испытаний: измерение сопротивления обмоток постоянному току, измерение сопротивления обмоток на холостом ходу и измерение сопротивлений при заторможенном роторе [1, 3, 4]. Чтобы избежать нарушения режима работы электромеханической системы и не отключать АД от сети и нагрузки, в настоящее время, как свидетельствуют недавние исследования [4 – 8], актуально использование расчетных оптимизационных методов оценки параметров асинхронных машин. Используя предоставленные производителем данные в качестве контрольных точек, авторы работ [4 – 8] осуществляют в реальном времени непрерывное сравнение текущих значений идентифицируемых параметров схемы замещения с адаптивной эталонной моделью, с применением современных алгоритмов нелинейной оптимизации [4], алгоритма вихревого поиска [5], алгоритма роя сальп [6], расширенного наблюдателя Калмана [7] и нейронных сетей [7, 8]. Получаемая в результате оптимизации математическая модель двигателя индивидуальна для каждого объекта исследования и максимально ему соответствует. Но в условиях невозможности проведения тестов на промышленном оборудовании, например, при отсутствии доступа исследователей на производство, для оценки параметров схемы замещения АД требуются совершенно иные, аналитические методы счета с применением ЭВМ. К такому же выводу приходят авторы работ [2, 9].
Соответственно, целью данной работы является реализация в удобном программном интерфейсе аналитической методики расчета параметров схемы замещения асинхронного двигателя с минимальными отклонениями от каталожных данных. В ходе работе выполнен обзор преимуществ и недостатков основных методик расчета параметров схемы замещения АД, проведён их сравнительный анализ на примере асинхронного двигателя марки BA280S4 мощностью 110 кВт, в среде Embarcadero Delphi 11 разработана программа расчета параметров схемы замещения по выбранной итерационной методике.
Методы расчета параметров схемы замещения асинхронного двигателя
Математическая и компьютерная модель АД строится по уравнениям для Г-образной (рис. 1) или Т-образной схемы замещения (рис. 2) обобщенной двухфазной электрической машины, приведённой к ортогональным осям. Моделируемая электрическая цепь АД характеризуется следующими электрическими и механическими величинами, обозначенными на рис. 1 и рис. 2: İ1 – комплексная величина тока в фазе обмотки статора; İ’2 – комплексная величина приведенного тока в фазе обмотки ротора; İ0 – комплексная величина тока в ветви намагничивания; U̇1 – комплексная величина фазного напряжения, приложенного к обмотке статора; s – скольжение асинхронного двигателя.
Рассчитываемые параметры схемы замещения – это активные и индуктивные сопротивления обмоток, обозначенные на рис. 1 и рис. 2: R1 – номинальное активное сопротивление статора; X1 – индуктивное сопротивление рассеяния статора; R’1 – приведенное активное сопротивление статора; X’1 – приведенное индуктивное сопротивление рассеяния статора; R’2 – приведенное активное сопротивление ротора; X’2 – приведенное индуктивное сопротивление рассеяния ротора; R0 – номинальное активное сопротивление намагничивания (опционально); X0 – номинальное индуктивное сопротивление намагничивания.
Рис. 1. Г-образная схема замещения асинхронного двигателя
Fig. 1. L-shaped asynchronous motor equivalent circuit
Рис. 2. Т-образная схема замещения асинхронного двигателя
Fig. 2. T-shaped asynchronous motor equivalent circuit
В библиотеках большинства программных комплексов в структуре готовых блоков моделей АД используется именно Т-образная схема замещения. Поэтому в работе рассматриваются методики расчета значений активных и индуктивных сопротивлений обмоток статора, ротора и ветви намагничивания только на основании Т-образной схемы замещения.
Исходными данными для расчета параметров схемы замещения по любой методике обязательно являются номинальные параметры, приводимые в паспорте: номинальная мощность Pн, номинальное линейное напряжение Uн, номинальный ток статора Iн, номинальная скорость nн, номинальный коэффициент полезного действия (КПД) ηн, номинальный коэффициент мощности cos φн, а также кратность пускового тока iп, кратность пускового момента mп, кратность критического момента mk. Значения этих параметров дают достоверную информацию о трёх точках механической и электромеханической характеристики АД: номинальный режим (при номинальном скольжении sн), пусковой режим (при скольжении s = 1), режим максимального момента (при критическом скольжении sk). В условиях отсутствия информации о эталонных значениях сопротивлений схемы замещения, степень адекватности их вычисления может быть оценена только по близости значений тока, момента, КПД и cos φ к паспортным данным, для чего выполняется расчет механической и электромеханической характеристики АД по полученной схеме замещения. Задача вычисления 5…6 параметров схемы замещения АД по трём точкам механической характеристики не имеет однозначного решения, отсюда следует несколько подходов к вычислению, и все они основаны на ряде допущений.
Одна группа аналитических методов, представленных в статьях Сафаряна В.С. [9], Макеева М.С. [10], Мощинского Ю.А. [11], Осипова В.С. [12], относится к безитерационным аппроксимационным. Данные методики отличаются последовательностью, простотой и быстротой счета, в случае необходимости расчет может быть выполнен вручную с использованием только калькулятора. Недостатком перечисленных методик является то, что аппроксимация функции сопротивления цепи ротора в зависимости от скольжения требует дополнительных значений величин, отличных от паспортных данных для номинального режима. Такими данными выступают ток холостого хода [10], энергетические показатели (КПД и cos φ) при загрузке 25 %, 50 %, 75 %, 80 % или 125 % номинальной [11], механические, магнитные и добавочные потери [12]. Все вышеперечисленные значения чаще всего задаются произвольно, как, например, активное сопротивление ротора [9]. Эти допущения приводят к расхождениям при аппроксимации характеристик и к значительным отклонениям расчетных значений параметров от каталожных. Как показано в [13], наибольшей точностью аппроксимации отличается методика Ю.А. Мощинского, в которой приводятся зависимости активного и индуктивного сопротивления ротора от скольжения R’2(s), X’2(s) с учетом эффекта вытеснения тока. Номинальные значения сопротивлений ротора достигаются в момент пуска при единичном скольжении: R’2(1) = R’20, X’2(1) = X’20.
Для подбора параметров в широком диапазоне значений применяются итерационные методы вычисления, такие, как методика C.Г. Германа-Галкина [14], В.Б. Терёхина [15], А.А. Усольцева и Д.В. Лукичёва [16]. В приведённых методиках значения активных и реактивных сопротивлений статора и ротора вычисляются с использованием коэффициента c1, являющегося коэффициентом пропорциональности индуктивностей обмотки статора и ветви намагничивания:
Вычисления параметров схемы замещения производятся в диапазоне 1,02 < c1 < 1,2, и на каждой итерации значение c1 по формуле (1) сравнивается с исходным. Для сокращения числа итераций, как предложено в [15, 16], рассчитанное значение c1 становится исходным для следующей итерации. Такой расчет сходится за 3…4 итерации, но получаемые по схеме замещения значения контрольных параметров обычно сильно отклоняются от каталожных данных, особенно в области скольжений, превышающих критическое.
Подбор параметров схемы замещения АД для совпадения расчетной характеристики со всеми контрольными точками является задачей многокритериальной нелинейной оптимизации, поэтому к её решению применим подход, представленный в [17]. В данной методике, помимо активных и индуктивных сопротивлений статора, ротора и ветви намагничивания, оценке подлежат ещё три параметра АД: a – доля магнитных потерь ротора; h’ – приведенная глубина паза ротора; α – доля неактивной части обмотки ротора. Начальные приближения всех 9 параметров определяются по классическим формулам для Т-образной схемы замещения [17]. Итерационный подбор выполняется методом покоординатного спуска, при котором значение каждого из 9 параметров изменяется в (1 ± Δ) раз в большую и в меньшую сторону, где Δ – точность счёта, по умолчанию принимаемая Δ = 0,005. Для каждого изменённого значения параметра при всех неизменных остальных вычисляется вектор критериев (контрольных точек) – значения номинального и пускового момента и тока, критического момента, КПД и коэффициента мощности в относительных единицах (о.е.):
Из вектора критериев (2) на каждой итерации формируется и сохраняется в матрице общий критерий оптимизации – целевая функция. В качестве целевой функции выступает аддитивная свертка взвешенных квадратов отклонений значений контрольных точек от каталожных данных:
где vi, о.е. – вес i-го критерия; Aнi, Ai, о.е. – каталожное и текущее значение i-го контрольного параметра АД из (2).
Из полученной матрицы критериев (значений минимизируемой целевой функции) выбирается наименьшее значение, и соответствующее ему изменённое значение параметра схемы замещения сохраняется, заменяя собой исходное. На этом завершается одна итерация, и далее начинается следующая. Расчет останавливается, когда минимизируемая целевая функция (3) перестаёт уменьшаться, т.е. на текущей итерации нет улучшений.
Если в процессе подбора параметров какой-либо из них выходит за границы ограничений: 0 < R1; 0 < X1; 0,05X0 < R0 < 0,2X0; 0 < X0 < 1000X1; 0 < c1 < 1000X’20; 0 < R’2(s); 0 < X’2(s); 0,2 < a < 0,8; 0,005 < α < 0,2; 1 < h’ < 7; то соответствующее ему значение целевой функции (3) в матрице критериев приравнивается к 106, что практически соответствует бесконечности [17].
Результаты работы и их обсуждение
Итерационная методика многокритериальной оптимизации параметров схемы замещения АД [17] реализована в программной среде Embarcadero Delphi [18]. Блок-схема её алгоритма приведена на рис. 3.
Рис. 3. Блок-схема алгоритма разработанной программы
Fig. 3. Block diagram of the algorithm for the developed program
Результаты работы программы показали, что используемая методика, благодаря учёту активного сопротивления ветви намагничивания и эффекта вытеснения тока в роторе, отличается наибольшим количеством совпадений с каталожными данными. Для примера асинхронного двигателя марки BA280S4 мощностью 110 кВт в табл. 1 представлены результаты программного расчёта сопротивлений схемы замещения АД в сравнении с результатами расчета по методикам Мощинского Ю.А., Германа-Галкина С.Г., Терёхина В.Б., Усольцева А.А.
Таблица 1
Расчетные значения сопротивлений схемы замещения асинхронного двигателя BA280S4
Table 1
Calculated values of the BA280S4 asynchronous motor equivalent circuit resistances
|
Методика |
R1, Ом |
R’20, Ом |
X1, Ом |
X’20, Ом |
X0, Ом |
R0, Ом |
|
Многокритериальная |
0,057698 |
0,010811 |
0,081539 |
0,140675 |
3,283204 |
0,164457 |
|
Мощинский Ю.А. и др. [11] |
0,030319 |
0,012625 |
0,123251 |
0,170204 |
3,467416 |
- |
|
Герман-Галкин С.Г. [14] |
0,023253 |
0,009213 |
0,041623 |
0,041623 |
2,173841 |
- |
|
Терёхин В.Б. [15] |
0,011843 |
0,046936 |
0,079600 |
0,079600 |
4,562045 |
- |
|
Усольцев А.А. и др. [16] |
0,010948 |
0,010948 |
0,177112 |
0,088562 |
3,413093 |
- |
Рассчитанные значения контрольных параметров АД BA280S4 и их максимальные и минимальные относительные отклонения от каталога сведены в табл. 2.
Таблица 2
Значения контрольных параметров асинхронного двигателя BA280S4 и их относительные отклонения от каталога
Table 2
Values of the BA280S4 asynchronous motor check parameters and their deviations from the catalog
|
Методика\Параметр |
mн, о.е. |
iн, о.е. |
mп, о.е. |
iп, о.е. |
mk, о.е. |
ηн, о.е. |
cos φн, о.е. |
|
Каталожные значения |
1 |
1 |
2,1 |
6,4 |
2,2 |
0,951 |
0,87 |
|
Итерационная многокритериальная |
1,0433 |
1,0468 |
0,2598 |
4,9800 |
2,2616 |
0,9280 |
0,8886 |
|
Многокритериальная с учетом вытеснения тока |
0,9927 |
1,0098 |
2,1195 |
6,2812 |
2,2181 |
0,9149 |
0,8891 |
|
Мощинский Ю.А. и др. [11] |
1,0254 |
1,0285 |
1,3737 |
6,2080 |
2,2206 |
0,9633 |
0,8563 |
|
Герман-Галкин С.Г. [14] |
1,3315 |
1,2848 |
1,4996 |
12,6131 |
5,9239 |
0,9595 |
0,8936 |
|
Терёхин В.Б. [15] |
0,2752 |
0,3341 |
2,1198 |
6,6200 |
3,7893 |
0,9848 |
0,6920 |
|
Усольцев А.А. и др. [16] |
1,1002 |
1,0671 |
0,2002 |
4,3759 |
2,3471 |
0,9778 |
0,8724 |
|
Мин. отклонение εmin, о.е. |
0,0073 |
0,0098 |
0,0093 |
0,0186 |
0,0082 |
0,0089 |
0,0028 |
|
Макс. отклонение εmax, о.е. |
0,7248 |
0,6659 |
0,9047 |
0,9708 |
1,6927 |
0,0380 |
0,2046 |
Полученные результаты (выделенные темно-серым цветом) подтверждают, что используемая в программе итерационная многокритериальная методика с учетом вытеснением тока отличается приемлемой точностью счета: расхождение для двигателей средней мощности не превышает 2 % для механических и электромеханических характеристик и достигает максимум 4 % для энергетических характеристик (КПД), что укладывается в допустимый заводской разброс характеристик, равный 5 %. Как иллюстрируют приведённые на рис. 4 естественные механические характеристики, построенные для примера BA280S4, расчёт по программе [18] даёт полное совпадение всех трёх точек номинального, пускового и критического режима, в то время как остальные рассмотренные методики позволяют достичь совпадения только одной или двух контрольных точек.
Рис. 4. Естественные механические характеристики асинхронного двигателя BA280S4
Fig. 4. Natural mechanic characteristics of the BA280S4 asynchronous motor
Заключение
Разработанная компьютерная программа объемом 16 МБ – это новый, удобный инструмент, совместимый с любой 32-битной и 64-битной версией Windows, не имеющий особых требований к процессору и оперативной памяти. Программа осуществляет вывод исходных данных и результатов расчета (параметров схемы замещения, контрольных точек и относительных отклонений) в текстовые поля, в файл и на график механической характеристики.
Программа может применяться: для расчета параметров двигателей в условиях отсутствия доступа к промышленному оборудованию в целях верификации; для разработки продвинутых компьютерных моделей трехфазных асинхронных двигателей с динамическими постоянными времени ротора; для составления каталогов исходных данных по двигателям для расчетно-графических и контрольных работ; для расчета недостающих и точных значений сопротивлений обмоток в целях построения достоверных естественных и искусственных механических характеристик при исследовании регулирования асинхронных двигателей. Быстрота и доступность расчетов с помощью разработанной программы обеспечит наибольшую достоверность математического и компьютерного моделирования систем электропривода для решения учебно-методических и научно-исследовательских задач.
1. Khalaf M. M., Ali A. M. Parameters Estimation Tests of Induction Machine Using Matlab/Simulink. Journal of Physics: Conference Series. IOP Publishing, 2021;1973(1):012109.
2. Gör H. Feasibility of Six Metaheuristic Solutions for Estimating Induction Motor Reactance. Mathematics. 2024; 12(3):483.
3. Ершов М.С., Комков А.Н., Гамидов О.М. Определение параметров схемы замещения и основных характеристик короткозамкнутых асинхронных двигателей малой мощности на основе анализа сигналов тока и напряжения // Электротехника. – 2022. – № 2. – С. 63-68.
4. Simba K.G., Quilumba F.L., Granda N.V. Parameter Estimation of a Three-Phase Induction Motor from Direct Starting Stator Transient Measurements. 2020 IEEE ANDESCON, Quito, Ecuador. IEEE, 2020. Pp. 1-5.
5. Montano J., Garzón O.D., Herrera-Jaramillo D.A., Montoya O.D., Andrade F., Tobón A. Estimating the Parameters of a Three-Phase Induction Motor using the Vortex Search Algorithm. Iranian Journal of Science and Technology, Transactions of Electrical Engineering. 2023;48:337-347.
6. Zorig A., Belkheiri A., Bendjedia B., Kouzi K., Belkheiri M. New identification of induction machine parameters with a meta-heuristic algorithm based on least squares method. COMPEL - The international journal for computation and mathematics in electrical and electronic engineering. 2023;42(6):1852-1866.
7. Tang J., Yang Y., Blaabjerg F., Chen J., Diao L., Liu Z. Parameter Identification of Inverter-Fed Induction Motors: A Review. Energies. 2018;11(9):2194.
8. Омельченко Е.Я., Лымарь А.Б. Идентификация параметров схемы замещения асинхронных двигателей при помощи нейронных сетей // Электротехнические и информационные комплексы и системы. – 2023. – Т. 19. – № 4. – С. 31-44.
9. Сафарян В.С., Геворгян С.Г. Определение параметров схемы замещения асинхронной машины // Энергетика. Известия высших учебных заведений и энергетических объединений СНГ. – 2015. – № 6. – С. 20-34.
10. Макеев М.С., Кувшинов А.А. Алгоритм расчета параметров схемы замещения асинхронного двигателя по каталожным данным // Вектор науки Тольяттинского государственного университета. – 2013. – №1 (23). – С. 108-112.
11. Мощинский Ю.А., Беспалов В.Я., Кирякин А.А. Определение параметров схемы замещения асинхронных двигателей по каталожным данным // Электричество. – 1998. – №4. – С. 38-42.
12. Осипов В.С. Аналитический метод расчета параметров схемы замещения трехфазных асинхронных двигателей серии АИР // Вестник Самарского государственного технического университета. Серия: Технические науки. – 2017. – №2 (54). – С. 108-120.
13. Пиляев С.Н., Афоничев Д.Н. Обоснование параметров схемы замещения асинхронного электродвигателя // Вестник Воронежского государственного аграрного университета. – 2020. – Т. 13. – № 4(67). – С. 129-138.
14. Герман-Галкин С.Г. Matlab & Simulink. Проектирование мехатронных систем на ПК: учеб. пособие для вузов / С.Г. Герман-Галкин. Санкт-Петербург: КОРОНА Век, 2008. – 368 с.
15. Терёхин В.Б. Моделирование систем электропривода в Simulink (Matlab 7.0.1): учебное пособие. – Томск: Изд-во Томского политехнического университета, 2008. – 320 с.
16. Усольцев А.А., Лукичёв Д.В. Определение параметров модели асинхронного двигателя по справочным данным // Приборостроение. – 2008. – №10. – С. 35-41.
17. Ершов М.С., Егоров А.В., Трифонов А.А. Устойчивость промышленных электротехнических систем. М.: Недра, 2010. – 318 с.
18. Свидетельство о государственной регистрации программы для ЭВМ № 2024611993 Российская Федерация. Программа расчета параметров Т-образной схемы замещения трехфазного асинхронного электродвигателя : № 2024611028 : заявл. 22.01.2024 : опубл. 26.01.2024 / Д. В. Дзюин.
