graduate student from 01.01.2023 until now
Moscow, Moscow, Russian Federation
employee
Moscow, Moscow, Russian Federation
The article is devoted to developing a program for calculating the parameters of a T-shaped equivalent circuit of three-phase asynchronous electric motors. The relevance of this work is due to the need to obtain reliable results of computer simulation modelling of electric drive systems in software packages. The aim of the study is to implement in a convenient software interface an analytical technique for calculating the parameters of an asynchronous motor equivalent circuit with minimal deviations from the catalogue data. For the software calculation of the equivalent circuit parameters the paper considers analytical approximation and iterative counting techniques, since they do not require experiments on industrial electrical equipment. A comparative analysis of such techniques is performed using the example of a three-phase asynchronous motor brand BA280S4 with a power of 110 kW. The results of the compar-ative analysis show that the highest reliability of the mechanical and electromechanical characteristics is achieved in the case of using the iterative counting technique with minimizing the objective function of the sum of the weighted squares of the deviations of the control point values from the catalogue data. The authors implement the program for calculating the equivalent circuit parameters according to the selected calculation method in the Embarcadero Delphi environment. The reliable values of the electric motor parameters calculated by the program will allow for computer to model electric drive systems that best matches real objects.
asynchronous electric motor, equivalent circuit, calculation of parameters, calculation program, iterative method, multi-criteria selection, objective function, mechanical characteristic, modelling
Введение
В настоящее время наиболее доступным и дешевым методом синтеза и исследования статических и динамических режимов работы электромеханических систем является компьютерное имитационное моделирование в программных комплексах, как например: 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. Ershov M.S., Komkov A.N., Gamidov O.M. Determination of the Equivalent-Circuit Parameters and the Main Characteristics of Squirrel-Cage Low-Power Induction Motors Based on Analysis of Current and Voltage Signals. Electrical Engineering. 2022;2:63-68.
4. Simba KG, Quilumba FL, Granda NV. Parameter Estimation of a Three-Phase Induction Motor From Direct Starting Stator Transient Measurements. In: Proceedings of 2020 IEEE ANDESCON; Quito, Ecuador: IEEE: 2020. p. 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. Omelchenko E.Ya., Lymar A.B. Identification of the Parameters of an Induction Motor Equipment Circuit Using Neural Networks. Electrical and Data Processing Facilities and Systems. 2023;19(4):31-44.
9. Safaryan V.S., Gevorgyan S.G. Ascertainment of the Equivalent Circuit Parameters of the Asynchronous Machine. Energetika. Proceedings of CIS Higher Education Institutions and Power Engineering Associations. 2015;6:20-34.
10. Makeev M.S., Kuvshinov A.A. Algorithm for Calculating the Parameters of the Equivalent Circuit of the Asynchronous Engine at the Catalogue Data. Science Vector of Togliatti State University. 2013;1(23):108-112.
11. Moshchinsky Yu.A., Bespalov V.Ya., Kiryakin A.A. Determining the Asynchronous Motor Equivalent Circuit From the Catalogue Data. Electricity. 1998;4:38-42.
12. Osipov V.S. Analytical Calculation of Equivalent Circuit Parameters for Three Phase Asynchronous Motors of AIR-Series. Vestnik of Samara State Technical University. Series: Technical Sciences. 2017;2(54):108-120.
13. Pilyaev S.N., Afonichev D.N. Substantiation of Parameters of the Asynchronous Electric Motor Replacement Circuit. Vestnik of Voronezh State Agrarian University. 2020;13-4(67):129-138.
14. German-Galkin S.G. Matlab & Simulink. Design of Mechatronic Systems on a PC. Saint Petersburg: KORONA Vek; 2008.
15. Terekhin V.B. Simulating Electric Drive Systems in Simulink (Matlab 7.0.1). Tomsk: Publishing House of Tomsk Polytechnic University; 2008.
16. Usoltsev A.A., Lukichev D.V. Definition of Parameters of Model of the Asynchronous Engine on the Help Data. Journal of Instrument Engineering. 2008;10:35-41.
17. Ershov M.S., Egorov A.V., Trifonov A.A. Stability of the Industrial Electrical Systems. Moscow: Nedra; 2010.
18. Dzyuin D.V. Program for Calculating the Parameters of a T-shaped Equivalent Circuit of a Three-Phase Asynchronous Electric Motor. The Certificate on Official Registration of the Computer Program in Russia. No. 2024611993; 2024 Jan 26.
