DEVELOPMENT AND MODELING OF THE THREE-ELEMENT SYSTEM OF CONTROL OF ELECTROMAGNETIC BEARING
Journal: VESTNIK OF ASTRAKHAN STATE TECHNICAL UNIVERSITY. SERIES: MANAGEMENT, COMPUTER SCIENCE AND INFORMATICS
Volume 2015 № 4
, 2015
Rubrics: MANAGEMENT, MODELING, AUTOMATION
Authors:
1.
Nizhny Novgorod State Technical University named after R. E. Alexeev
2.
Nizhny Novgorod State Technical University named after R. E. Alexeev
3.
Nizhny Novgorod State Technical University named after R. E. Alexeev
Type:
Article
Pages:
from 14
to 23
Status:
Published
Received:
05.10.2015
Accepted:
25.10.2015
Published:
25.10.2015
Subject area:
VAC 05.12.2013 Системы, сети и устройства телекоммуникаций
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 621.3
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 621.3
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
Language:
Russian
Keywords:
Matlab/Simulink, electromagnetic bearing, three-loop control system, current loop, proportional-integral controller, modular optimum, Matlab/Simulink software
Abstract (English):
A three-loop control system of radial (and axial) electromagnetic bearing for powerful rotating machines, working air gap of which between the stator and rotor are tenths of millimeter if rotor weighs about one ton, is proposed. In the open state the system is unstable and has the features, which must be considered while synthesizing the regulators. Setting up the system is made by sequential optimization of control loops in accordance with the general principles of construction of the systems of the subordinate regulation. In accordance with the mathematical model of the electromagnetic bearing in the form of transfer functions and the corresponding block diagrams, the simulation models of the three-loop control system in the software package Matlab Simulink are constructed. Initially the inner loop with proportional-integral regulator was studied. The introduction of the internal circuit improved the system performance, but the simulation showed that the three-loop control system of the electromagnetic bearing has a static error. To make the system isoc, an integral controller is introduced into the outer contour of the situation. The transient analysis in the simulation showed that the electromagnetic bearings, equipped with the three-loop control system, have a high static rigidity and high speed. The failures of the rotor under the influence of external disturbing forces allow us to speak about the dependence of the dynamic stiffness of the electromagnetic bearing on the time constant of integration in the external position loop. An example of the use of research results can be a gas blower with 16 MW of the compressor station "Pomarskaya" of the pipeline "Urengoy - Pomary - Uzhgorod".
A three-loop control system of radial (and axial) electromagnetic bearing for powerful rotating machines, working air gap of which between the stator and rotor are tenths of millimeter if rotor weighs about one ton, is proposed. In the open state the system is unstable and has the features, which must be considered while synthesizing the regulators. Setting up the system is made by sequential optimization of control loops in accordance with the general principles of construction of the systems of the subordinate regulation. In accordance with the mathematical model of the electromagnetic bearing in the form of transfer functions and the corresponding block diagrams, the simulation models of the three-loop control system in the software package Matlab Simulink are constructed. Initially the inner loop with proportional-integral regulator was studied. The introduction of the internal circuit improved the system performance, but the simulation showed that the three-loop control system of the electromagnetic bearing has a static error. To make the system isoc, an integral controller is introduced into the outer contour of the situation. The transient analysis in the simulation showed that the electromagnetic bearings, equipped with the three-loop control system, have a high static rigidity and high speed. The failures of the rotor under the influence of external disturbing forces allow us to speak about the dependence of the dynamic stiffness of the electromagnetic bearing on the time constant of integration in the external position loop. An example of the use of research results can be a gas blower with 16 MW of the compressor station "Pomarskaya" of the pipeline "Urengoy - Pomary - Uzhgorod".
Keywords:
Matlab/Simulink, electromagnetic bearing, three-loop control system, current loop, proportional-integral controller, modular optimum, Matlab/Simulink software
Matlab/Simulink, electromagnetic bearing, three-loop control system, current loop, proportional-integral controller, modular optimum, Matlab/Simulink software
Состояние проблемы Основным элементом многих энергетических машин является ротор, вращающийся в подшипниковых опорах. Рост скорости вращения и необходимость снижения потерь мощных роторных машин требуют разработки принципиально нового вида подшипников, в которых для создания опорных реакций используются магнитные и электрические поля. Среди них наибольший практический интерес представляют активные магнитные подшипники (АМП). Применение электромагнитных подшипников (ЭМП) в конструкции мощных электродвигателей позволяет также увеличить срок их службы и расширить эксплуатационные возможности создаваемых агрегатов [1]. Несмотря на очевидные преимущества АМП, они не находят широкого применения в силу ряда причин, основная из которых - сложность технической реализации ЭМП и, как следствие, их высокая стоимость. Высокая сложность и стоимость связаны, как правило, с вопросами технической реализации системы автоматического управления (САУ), без которой данный вид опор в принципе функционировать не может. Для структурного и параметрического синтеза системы управления используется математическая модель одной оси ЭМП [2]: (1) где U1, U3 - напряжение питания обмоток электромагнитов, В; I1, I3 - ток питания обмоток электромагнитов, А; L1, L3 - индуктивность обмоток электромагнитов, Гн; R1, R3 - активное сопротивление обмоток электромагнитов, Ом; kE1, kE3 - коэффициенты ЭДС, В∙с/м; y - возможное перемещение ротора в поле электромагнитов, м; m - масса ротора, приходящаяся на одну ось подшипника, кг; kэм - коэффициент электромагнитной силы, Н/А; kF - коэффициент, связывающий приращение электромагнитного усилия и смещение ротора (коэффициент положительной обратной связи по перемещению), Н/м; G - вес ротора, вращающегося в поле электромагнитов, Н; Fв - равнодействующая внешних сил, действующих на ротор, Н (те же самые обозначения используются и далее, в том числе и в рисунках). Приведенной математической модели оси ЭМП (1) соответствует структурная схема на рис. 1. Рис. 1. Структурная схема одной оси электромагнитного подшипника: UР - выходное напряжение регулятора, В; kп - коэффициент передачи силового преобразователя; Тп - постоянная времени силового преобразователя, с; kE - коэффициент ЭДС, В∙с/м; L - индуктивность обмоток электромагнитов, Гн; m - масса ротора, приходящаяся на одну ось подшипника, кг; y - возможное перемещение ротора в поле электромагнитов, м Постановка задачи В общем случае при создании систем управления электромагнитным подвесом ротора преследуется цель обеспечить высокое быстродействие, точность поддержания ротора машины в центральном положении, демпфирование возможных вибраций. Сложность современных систем управления ЭМП обусловлена методами, которые применяются при их синтезе. Известно несколько принципов построения систем управления ЭМП: - одноконтурные системы управления, синтезированные методом решения обратных задач динамики; - двухконтурные и трехконтурные системы управления, построенные по принципу подчиненного регулирования координат. Нами, на примере газового нагнетателя мощностью 16 МВт перекачивающей станции «Помарская» газопровода «Уренгой - Помары - Ужгород», рассмотрена возможность построения трехконтурной системы управления радиальными ЭМП для мощных вращающихся агрегатов. Рабочие зазоры между статором и ротором в подобных установках составляют десятые доли миллиметра при весе ротора нагнетателя около одной тонны. Методы исследования Структурная схема трехконтурной системы управления представлена на рис. 2. Настройка системы производится путем последовательной оптимизации контуров регулирования. Применение такого подхода требует выделения в объекте управления ряда координат (скорость, перемещение, ток), по которым будут замкнуты соответствующие контуры регулирования. Система управления, синтезируемая по такому принципу, характеризуется простотой технической реализации (средствами как аналоговой, так и микропроцессорной техники) и простотой ограничения регулируемых координат (скорость, перемещение, ток). В качестве датчика главной обратной связи используется бесконтактный датчик, измеряющий перемещение ротора в поле электромагнитов. Коэффициент передачи датчика kд.п, В/м, выбран таким образом, чтобы при максимально возможном смещении ротора y = 0,2 мм сигнал датчика был равен 10 В. Для выделения сигнала обратной связи по скорости перемещения сигнал датчика положения ротора дифференцируется и вводится в систему с коэффициентом передачи kосс, В/с-1. Рис. 2. Структурная схема трехконтурной системы управления электромагнитным подшипником: Uз - напряжение задания на входе регулятора, В; Wр.п, Wр.с, Wр.т - передаточные функции регуляторов положения, скорости, тока; kд.т - коэффициент передачи датчика тока, В/А; kосс - коэффициент передачи датчика скорости На начальном этапе производится исследование внутреннего контура тока. В разомкнутом состоянии система является неустойчивой. Введение внутреннего контура тока позволит снизить чувствительность синтезируемой системы к внешним возмущениям, а также увеличить быстродействие системы за счёт более быстрого протекания переходных процессов. Следует отметить, что постоянная времени внутреннего контура тока достаточно велика (~ 50 мс), что приводит к необходимости большой форсировки по источнику питания [3]. В соответствии со структурной схемой на рис. 3 внутренний контур образуется регулятором тока, силовым преобразователем, объектом регулирования и обратной связью по току через датчик тока. Рис. 3. Структурная схема внутреннего контура тока Контур тока настраивается на модульный оптимум, тогда передаточная функция регулятора тока имеет вид [3]: , где - коэффициент передачи силового преобразователя; - коэффициент передачи датчика тока; L - индуктивность электромагнита, Гн; r - активное сопротивление электромагнита, Ом; Тμ - малая постоянная времени во внутреннем контуре тока, с. В качестве регулятора тока в системе подчиненного регулирования электромагнитного подвеса ротора использован пропорционально-интегральный регулятор (ПИ-регулятор). Настройки ПИ-регулятора определяются следующими параметрами: - коэффициент усиления: - постоянная времени интегрирования, с: Техническая реализация такого регулятора возможна как средствами аналоговой техники на операционном усилителе, так и программно - на микропроцессорном устройстве. В соответствии со структурной схемой на рис. 3 построена имитационная модель в программном комплексе Matlab Simulink для анализа переходных процессов во внутреннем контуре тока с ПИ-регулятором (рис. 4). Рис. 4. Имитационная модель контура тока с ПИ-регулятором: kп - коэффициент усиления регулятора тока; Ти - постоянная времени интегрирования регулятор тока, с Следует отметить, что для стабилизации выходного напряжения силового преобразователя на уровне 220 В в модели установлен блок «Ограничение напряжения». Тогда переходные процессы напряжения и тока будут иметь вид, представленный на рис. 5 и 6. Рис. 5. Выходное напряжение силового преобразователя, В Рис. 6. Переходные процессы во внутреннем контуре с ПИ-регулятором тока Анализ переходного процесса (рис. 6) показывает, что внутренний контур тока с ПИ-регулятором является астатическим как по заданию, так и по возмущению. Перерегулирование, при оптимизации на модульный оптимум, составляет 4,3 %. Подобный характер переходного процесса объясняется тем, что выходное напряжение силового преобразователя ограничивается на уровне 220 В [3]. Найдена передаточная функция регулятора скорости, обеспечивающего модульный оптимум замкнутому контуру скорости [4]: . Следовательно, для настройки контура скорости на модульный оптимум используется пропорциональный регулятор (П-регулятор) с коэффициентом передачи: . При настройке внешнего контура положения на модульный оптимум найдена передаточная функция регулятора положения, определяемая выражением . Следовательно, во внешнем контуре использован П-регулятор с коэффициентом передачи . В соответствии со структурной схемой на рис. 2 и настройками регуляторов, система управления реализована в виде соответствующей имитационной модели в программном комплексе Matlab Simulink (рис. 7) [5]. Моделирование позволит оценить качество переходных процессов при реализации системы управления с внутренним контуром тока. Рис. 7. Имитационная модель трехконтурной системы управления электромагнитным подшипником: kр.п, kр.с, kр.т - коэффициенты усиления регуляторов положения, скорости и тока соответственно Вид переходного процесса перемещения ротора в поле электромагнитов при выбранных параметрах регуляторов представлен на рис. 8. Рис. 8. Переходный процесс перемещения ротора при единичном задании В рассматриваемой системе разность между максимальным значением выходной координаты y и ее установившимся значением yуст составляет σ = 10 % . Время переходного процесса tп.п = 0,002 с. После окончания переходных процессов в системе устанавливается режим, когда выходная координата y максимально приблизилась к заданному значению yуст. Изменение режима работы системы возникает в результате прикладываемого к ней внешнего воздействия. Внешнее возмущение в виде скачкообразной нагрузки наиболее характерно для систем стабилизации, поэтому основная задача состоит в том, чтобы система была астатической по возмущению. Статическая нагрузка, определяемая массой вращающейся части (m = 1000 кг), известна и равна весу нагнетателя G = 10 000 Н. Максимальное усилие Fв должно быть больше статической нагрузки с целью обеспечения обработки динамических нагрузок и выхода ЭМП в рабочий режим. Обычно Fв принимают равным двум весам ротора. Следует оценить переходный процесс, возникающий в системе при отработке возмущающего воздействия Fв = 20 000 Н (рис. 9). Рис. 9. Переходный процесс по возмущению при Fв = 20 000 Н Анализ переходного процесса показывает, что в трехконтурной системе управления ЭМП имеет место статическая ошибка. Для достижения астатизма системы в прямую цепь САР включен интегральный регулятор (И-регулятор) с постоянной интегрирования Ти = 8∙Тμ, с. При выбранном значении постоянной времени Ти включение И-регулятора исключает статическую ошибку, возникающую в системе при внешнем воздействии, и практически не оказывает влияния на устойчивость системы. Вид преобразованной имитационной модели астатической системы управления представлен на рис. 10. Рис. 10. Преобразованная имитационная модель трехконтурной системы управления электромагнитным подшипником с И-регулятором во внешнем контуре положения Результаты исследования Результаты исследования приведены на рис. 11. Рис. 11. Переходный процесс по возмущению при Fв = 20 000 Н: 1 - статическая трехконтурная САУ; 2 - астатическая трехконтурная САУ Из рис. 11 видно, что система с интегральным звеном в прямом канале регулирования стала астатической по возмущению. При выбранной структуре системы управления переходный процесс по управляющему воздействию будет иметь вид, представленный на рис. 12. Рис. 12. График переходного процесса по управляющему воздействию Из рис. 12 видно, что в рассматриваемой системе разность между максимальным значением выходной координаты y и ее установившимся значением yуст составляет σ = 40 %. Следует отметить, что в системах стабилизации допустимо перерегулирование свыше 50 % и оно не влияет на их устойчивость. Время переходного процесса, т. е. интервал времени, по истечении которого отклонение не превышает 5 %, tп.п = 0,007 с. После окончания переходных процессов в системе устанавливается режим, когда выходная координата y максимально приблизилась к заданному значению yуст. Выводы Анализ переходных процессов показывает, что электромагнитные подшипники, оснащенные трехконтурной системой управления, обладают большой статической жесткостью и высоким быстродействием. Провалы ротора под действием внешних возмущающих сил позволяют говорить о зависимости динамической жесткости электромагнитного подшипника от постоянной времени Ти во внешнем контуре положения. Полученные результаты могут служить теоретической основой для решения актуальных задач, связанных с движением роторов в поле электромагнитных подшипников, таких как исследование устойчивости, анализ динамики переходных режимов (пуск-останов, смена режима работы), отклик системы на ударное (импульсное) возмущение.
1. Zhuravlev Yu. N. Aktivnye magnitnye podshipniki: teoriya, raschet, primenenie / Yu. N. Zhuravlev. SPb.: Politehnika, 2003. 206 s.
2. Sarychev A. P. Matematicheskaya model' rotora dlya analiza upravleniya magnitnymi podshipnikami / A. P. Sarychev, I. G. Rukovicyn // Tr. NPP VNIIEM. 2008. T. 107. S. 11-15.
3. Naumec A. E. Optimizaciya kontura regulirovaniya toka elektromagnitnogo podshipnika / A. E. Naumec, R. S. Taganov, G. M. Miryasov, V. G. Titov, A. N. Kuzmenkov // Tr. Nizhegorod. gos. tehn. un-ta im. R. E. Alekseeva. 2013. № 5 (102). S. 298-303.
4. Kuzmenkov A. N. Upravlenie magnitnymi podshipnikami moschnyh vraschayuschihsya agregatov / A. N. Kuzmenkov, G. M. Miryasov, V. G. Titov // Tr. VIII Mezhdunar. (XIX Vseros.) konf. po avtomat. elektroprivodu AEP-2014. Saransk: Izd-vo Mordov. un-ta, 2014. T. 1. S. 226-228.
5. Chernyh I. V. Modelirovanie elektrotehnicheskih ustroystv v MatLAB, Simulink i SimPowerSystems / I. V. Chernyh. M.: DMKPress; SPb.: Piter, 2008. 288 s.
