WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

Pages:     || 2 | 3 |

«АДАПТИВНЫЕ АЛГОРИТМЫ ОЦЕНИВАНИЯ КООРДИНАТ БЕЗДАТЧИКОВЫХ ЭЛЕКТРОПРИВОДОВ ПЕРЕМЕННОГО ТОКА С РАСШИРЕННЫМ ДИАПАЗОНОМ РЕГУЛИРОВАНИЯ ...»

-- [ Страница 1 ] --

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ

ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО

ОБРАЗОВАНИЯ

«НОВОСИБИРСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

На правах рукописи

ВДОВИН ВЛАДИМИР ВЛАДИМИРОВИЧ

АДАПТИВНЫЕ АЛГОРИТМЫ ОЦЕНИВАНИЯ КООРДИНАТ

БЕЗДАТЧИКОВЫХ ЭЛЕКТРОПРИВОДОВ ПЕРЕМЕННОГО ТОКА

С РАСШИРЕННЫМ ДИАПАЗОНОМ РЕГУЛИРОВАНИЯ

Специальность: 05.09.03 – Электротехнические комплексы и системы Диссертация на соискание ученой степени кандидата технических наук

Научный руководитель – д. т. н., профессор В.В. Панкратов Новосибирск –

СПИСОК СОКРАЩЕНИЙ И УСЛОВНЫХ ОБОЗНАЧЕНИЙ

АД – асинхронный двигатель АДКЗР – асинхронный двигатель с короткозамкнутым ротором АДФР – асинхронный двигатель с фазным ротором АНПП – адаптивный НПП АСЗМ – адаптивная система с задающей (эталонной) моделью АСМ – асинхронизированная синхронная машина АЦП – аналого-цифровое преобразование МДП – машина двойного питания МЦС – модель цепи статора НПП – наблюдатель полного порядка ПК – преобразователь координат ПЧ – преобразователь частоты РП – регулятор потокосцепления РС – регулятор скорости РТ – регулятор тока СВУ – система векторного управления СДПМ – синхронный двигатель с постоянными магнитами ШИМ – широтно-импульсная модуляция ЭП – электропривод EKF – Extended Kalman Filter (расширенный фильтр Калмана) FOC – Field-Oriented Control (векторное управление) FOO – Full Order Observer (наблюдатель полного порядка) AFOO – Adaptive Full Order Observer (адаптивный наблюдатель полного порядка) IGBT – Insulated-gate bipolar transistor (биполярный транзистор с изолированным затвором) MRAS – Model Reference Adaptive System (адаптивная система с задающей моделью)

ОГЛАВЛЕНИЕ

Стр.

СПИСОК СОКРАЩЕНИЙ И УСЛОВНЫХ ОБОЗНАЧЕНИЙ

ВВЕДЕНИЕ

ГЛАВА 1. ВЕКТОРНОЕ УПРАВЛЕНИЕ ЭЛЕКТРИЧЕСКИМИ МАШИНАМИ

ПЕРЕМЕННОГО ТОКА В УСЛОВИЯХ НЕПОЛНЫХ ИЗМЕРЕНИЙ.............. 1.1 Математическое моделирование машин переменного тока

1.1.1 Преобразования координат

1.1.2 Баланс мощностей и определение электромагнитного момента.... 1.1.3 Математическая модель АДКЗР как объекта управления............... 1.1.4 Математическая модель АСМ

1.1.5 Математическая модель СДПМ

1.2 Векторное управление машинами переменного тока на основе идеализированного ПЧ

1.2.1 Принцип векторного управления АДКЗР

1.2.2 Принцип векторного управления СДПМ

1.2.3 Принцип векторного управления АСМ

1.2.4 Способы автоматической ориентации вращающейся системы координат и вектора управляющих воздействий по магнитному полю двигателя

1.3 Бездатчиковое векторное управление

1.3.1 Методы оперативного оценивания координат состояния............... 1.3.2 Адаптивное управление в бездатчиковых СВУ с оцениванием координат по основным рабочим гармоникам электрических величин. 1.4 Выводы

ГЛАВА 2. СИНТЕЗ И ИССЛЕДОВАНИЕ АЛГОРИТМА ВЫЧИСЛЕНИЯ

ЧАСТОТЫ ВРАЩЕНИЯ И ОПОРНОГО ВЕКТОРА ПОТОКОСЦЕПЛЕНИЙ

АДКЗР, МАЛОЧУВСТВИТЕЛЬНОГО К ДРЕЙФУ АКТИВНОГО

СОПРОТИВЛЕНИЯ СТАТОРА

2.1 Обзор адаптивных алгоритмов вычисления координат АДКЗР.............. 2.1.1 Алгоритмы типа АСЗМ

2.1.2 Алгоритм на основе НПП (Кубота)

2.1.3 Алгоритм на основе АНПП (Хинкканен)

2.1.4 Выводы

2.2 Предлагаемый метод структурного синтеза алгоритмов вычисления.... оценок координат АДКЗР на основе адаптивного наблюдателя состояния. 2.2.1 Синтез наблюдателя полного порядка

вращения ротора АД

вращения ротора АД и активного сопротивления статорной цепи......... вычисления частоты вращения и сопротивления статора

2.3 Сравнительный анализ предлагаемого алгоритма с известными............ идентификации АДКЗР к отклонениям параметров

2.5 Предварительная идентификация параметров АД

2.6 Выводы

ГЛАВА 3. СИНТЕЗ И ИССЛЕДОВАНИЕ ПОДСИСТЕМЫ ОЦЕНИВАНИЯ

ЧАСТОТЫ ВРАЩЕНИЯ И УГЛОВОГО ПОЛОЖЕНИЯ СДПМ,

МАЛОЧУВСТВИТЕЛЬНОЙ К ДРЕЙФУ АКТИВНОГО СОПРОТИВЛЕНИЯ

СТАТОРА

3.1 Обзор алгоритмов оценивания неизмеряемых координат СДПМ......... 3.1.1 Бездатчиковый алгоритм управления СДПМ, использующий «расширенную ЭДС»

3.2 Структурный синтез алгоритма оценивания координат неявнополюсного СДПМ на основе адаптивного наблюдателя состояния

3.2.1 Cинтез наблюдателя на основе прямого метода Ляпунова........... 3.2.2 Синтез адаптора по частоте вращения

3.2.2 Синтез адаптора активного сопротивления статора

3.3 Исследование чувствительности адаптивного алгоритма оценивания координат СДПМ к отклонениям параметров

3.4 Предварительная идентификация параметров СДПМ. Определение начального положения ротора. Раскрутка постоянным током



3.5 Выводы

ГЛАВА 4. СИНТЕЗ И ИССЛЕДОВАНИЕ ПОДСИСТЕМЫ ОЦЕНИВАНИЯ

ЧАСТОТЫ ВРАЩЕНИЯ И ОПОРНОГО ВЕКТОРА ПОТОКОСЦЕПЛЕНИЙ

СТАТОРА АСМ. ОБОБЩЕНИЕ МЕТОДИКИ СИНТЕЗА

4.1 Обзор алгоритмов оценивания координат АСМ

4.1.1 Ориентация вращающейся системы координат при прямом измерении углового положения ротора АСМ

4.1.2 Вычисление углового положения вектора потокосцеплений статора и частоты вращения ротора методом АСЗМ

4.1.3 Структура АСЗМ без измерения напряжений статора.................. 4.1.4 Структура АСЗМ с адаптацией по ЭДС статора

4.1.5 Вычисление углового положения вектора потокосцеплений статора и частоты вращения ротора АСМ при измерениях только переменных на выходе ПЧ

4.2 Синтез и исследование адаптивного наблюдателя координат АСМ..... 4.2.1 Синтез адаптивного наблюдателя координат АСМ

4.2.2 Исследование чувствительности алгоритма оценивания координат АСМ к отклонениям активных сопротивлений обмоток

4.2.3 Выводы

4.3 Общая методика синтеза алгоритмов оценивания состояния в условиях неполных измерений и неопределенности параметров объекта ................. 4.4 Выводы

ГЛАВА 5. ЭКСПЕРИМЕНТАЛЬНОЕ ИССЛЕДОВАНИЕ АДАПТИВНЫХ

АЛГОРИТМОВ БЕЗДАТЧИКОВОГО ВЕКТОРНОГО УПРАВЛЕНИЯ.......... 5.1 Иследование алгоритма АСЗМ

5.1.1 Эспериментальная установка для исследования АСЗМ................ 5.1.2 Результаты экспериментального исследования

5.1.3 Исследование чувствительности алгоритма АСЗМ к отклонениям параметров

5.2 Исследование алгоритма АНПП

5.3 Выводы

ЗАКЛЮЧЕНИЕ

СПИСОК ЛИТЕРАТУРЫ

ПРИЛОЖЕНИЯ

ПРИЛОЖЕНИЕ А

ПРИЛОЖЕНИЕ Б

ПРИЛОЖЕНИЕ В

ПРИЛОЖЕНИЕ Г

ПРИЛОЖЕНИЕ Д

ВВЕДЕНИЕ

Известно, что более 60% всей вырабатываемой в мире электроэнергии потребляется электроприводами (ЭП), причём большинство из них – это электроприводы общепромышленных механизмов (насосов, вентиляторов, компрессоров, транспортеров и т.п.), а также электроприводы различных технологических установок, не требующие глубокого регулирования и высокого быстродействия и изначально оснащенные электродвигателями переменного тока.

Применение в них двигателей постоянного тока (ДПТ), обладающих благоприятными регулировочными характеристиками и допускающих относительно простое построение системы управления, как правило, ограничено непропорционально высокими эксплуатационными затратами, низкой надежностью и невозможностью функционирования во взрывоопасных, загрязненных и агрессивных средах.

Возможность регулирования поступающей механической мощности требуется большинством технологических процессов. При нерегулируемом электроприводе оно обычно осуществляется путем управляемого изменения в самом технологическом процессе. В частности, в насосах, компрессорах, вентиляторах – путем открытия и закрытия заслонок, в металло – и деревообрабатывающих станках – применением механических коробок передач. Однако такие технические решения нельзя считать экономически эффективными. Практика показала, что при переходе к регулируемому электроприводу можно сэкономить до (30…60)% потребленной электроэнергии по сравнению с нерегулируемым. Кроме того, уменьшаются пусковые токи и броски момента, что продлевает срок службы механических узлов.

По этим причинам мировое сообщество активно переходит на использование регулируемого по скорости электропривода переменного тока. Это связано, прежде всего, с развитием соответствующей элементной базы: силовой полностью управляемой полупроводниковой техники, микроэлектроники и микропроцессорной техники, а также становлением теоретических основ – теории векторного управления электроприводами переменного тока, теории автоматического управления многосвязными нелинейными объектами. Все это позволяет осуществить построение высокоэффективных систем управления электроприводами переменного тока.

Основными машинами переменного тока средней и большой мощности, применяемыми в промышленности, являются асинхронные двигатели с короткозамкнутым ротором (АДКЗР), асинхронные двигатели с фазным ротором (АДФР) и синхронные двигатели с постоянными магнитами (СДПМ).

Широкое использование в электроприводах АДКЗР обуславливается их известными преимуществами – простотой конструкции ротора, не содержащего скользящих контактов и постоянных магнитов, вследствие этого – высокой технологичностью, низкой ценой и высокой надежностью, а также минимальной требовательностью к обслуживанию. Кроме того, АДКЗР не требуют организации специальных пусковых мероприятий при прямом питании от сети. Поэтому применяются, в основном, в нерегулируемом ЭП.

Исследования в области действующего электрооборудования отечественных промышленных предприятий показывают широкое применение в мощных (более 300 кВт) механизмах с тяжелыми условиями пуска (ЭП подъемнотранспортных механизмов, ЭП главного движения вращающихся цементных печей и пр.) асинхронных двигателей с фазным ротором. Применение в подобных механизмах ЭП на базе АДФР с резисторными релейно-контакторными системами управления объясняется высокой перегрузочной способностью АДФР (коэффициент перегрузки – 2,5, как у ДПТ), простотой регулирования путем введения в цепь ротора добавочных активных сопротивлений. С теоретической точки зрения для электрического регулирования частоты вращения АДФР возможно закоротить обмотку ротора и подключить статорную обмотку к преобразователю частоты (ПЧ), получив тем самым хорошо проработанный электропривод на базе асинхронного двигателя с короткозамкнутым ротором.

Однако большинство мощных АДФР имеет питание от промышленной сети 6 – 10 кВ, что потребует применения ПЧ с соответствующим выходным линейным напряжением. Текущий уровень развития силовой полупроводниковой техники требует построения такого преобразователя на базе многоуровневого преобразователя частоты, обладающего высокой стоимостью. Более дешевым вариантом является установка ПЧ в роторную цепь, линейное напряжение которой, как правило, составляет 500 – 1000 В. Статорная обмотка при этом питается от промышленной сети 6 – 10 кВ. Такое включение АДФР носит название «асинхронизированная синхронная машина» (АСМ).

СДПМ все чаще находят применение в ЭП, требующих высокого качества регулирования, вытесняя ДПТ. СДПМ обладают высокими удельными энергетическими показателями – высоким КПД, за счет отсутствия потерь в роторе, и высокими массогабаритными показателями. Кроме того, за счет использования постоянных магнитов всегда имеется номинальный магнитный поток, что позволяет исключить из пусковой диаграммы ЭП участок предварительного намагничивания, необходимый для АДКЗР. За счет отсутствия скользящих контактов СДПМ характеризуются малыми эксплуатационными затратами.

Современный ЭП переменного тока, как правило, содержит двухзвенный преобразователь частоты, выпрямитель которого нагружен на автономный инвертор напряжения (АИН), работающий в режиме широтно-импульсной модуляции (ШИМ). Указанная структура преобразователя частоты позволяет независимо от режимов работы ЭП обеспечить высокий коэффициент мощности силовой цепи. В случае необходимости длительных рекуперативных режимов работы вместо нерегулируемого выпрямителя целесообразно использовать активный выпрямитель напряжения (АВН) работающий в режиме стабилизации звена постоянного тока, который обеспечивает требуемый уровень напряжения звена постоянного тока независимо от просадок питающей сети, обеспечивая тем самым требуемый диапазон регулирования частоты вращения. АВН, а также позволяет формировать синусоидальный ток, потребляемый из сети, с требуемым cos. Для улучшения формы потребляемых токов на вход АВН устанавливается трехфазный силовой синусный фильтр [1].

В 1971 г. Ф. Блашке предложил новый принцип управления машинами переменного тока [2]. В соответствии с этим принципом электромагнитные процессы рассматриваются не в стационарной (кларковской) системе координат, а во вращающейся ориентированной системе координат Парка-Горева.

Ориентация системы координат, как правило, осуществляется по опорному вектору потокосцеплений. Для АДКЗР в качестве опорного вектора, как правило, используется вектор потокосцеплений ротора, для АДФР – вектор потокосцеплений статора, а для СДПМ – вектор потокосцеплений, созданных постоянными магнитами.

Согласно [3] для реализации систем векторного управления необходимо иметь информацию обо всем векторе состояния машины. Соответственно возникает проблема получения информации об опорном векторе потокосцеплений при отсутствии полных измерений. В разное время эту проблему решали путем измерения индукции магнитного поля в расточке статора двигателя – использовали датчики Холла – или же на основании косвенной информации о потоках, получаемой по ЭДС, наводимой в специально уложенной в статор измерительной обмотке. Для измерения угловой скорости ротора использовались тахогенераторы или импульсные датчики. Однако все это требовало замены или доработки огромного парка машин, которые не предназначались для регулирования и не были снабжены этими датчиками.

Нецелесообразность применения датчиков Холла и измерительных обмоток подтверждается тем фактом, что вследствие несинусоидальности распределения магнитного поля по воздушному зазору машины, зубцовых пульсаций магнитного потока и погрешностей установки самих датчиков Холла на геометрических осях фазных обмоток получаемые с их помощью оценки лишь приближенно пропорциональны по мгновенным значениям компонентам интегрального вектора главных потокосцеплений двигателя. Использование же измерительных обмоток основано на интегрировании наведенных в них ЭДС, что связано с известными трудностями и налагает весьма жесткие ограничения на диапазон рабочих частот по статорным переменным. В целом, все описанные выше способы получения информации о координатах состояния машин переменного тока значительно снижают надежность электропривода.

К началу 80-х годов были достигнуты значительные успехи в области современной теории автоматического управления, что позволило косвенно оценивать неизвестные составляющие вектора состояния машины. Наиболее перспективным способом получения информации о недоступных прямому измерению координатах машин переменного тока является построение быстродействующих алгоритмов или устройств оценивания – наблюдателей.

Наблюдатель – это специальная динамическая подсистема, выходные сигналы которой служат в качестве оценок неизмеряемых координат состояния объекта управления. Наблюдатель должен быть устойчивым по начальным условиям и внешним воздействиям, использовать сигналы датчиков только заранее определенных и доступных непосредственным измерениям физических величин, быть малочувствительным к ошибкам в априорной информации о параметрах объекта управления.

Разработка методов построения систем управления ЭП переменного тока с различными наблюдателями и алгебраическими вычислителями явилась предпосылкой для появления к середине 80-х годов первых ЭП, не оснащенных датчиками магнитного состояния и координат механического движения (скорости, положения ротора). Такие электроприводы получили название «бездатчиковые» (sensorless) и широко распространены в промышленности, однако до сих пор они имеют существенные ограничения по диапазону и режимам регулирования.

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

Все они совмещают вычисление оценки частоты вращения с вычислением модуля и углового положения опорного вектора потокосцеплений, а отличаются друг от друга точностью вычисления скорости, чувствительностью к дрейфу параметров, входящих в математическую модель наблюдателя, способностью функционировать в характерных областях на плоскости механических характеристик. Диапазон регулирования скорости в двигательном режиме в практических разработках бездатчиковых электроприводов не превышает 50…100:1, а в режимах генераторного торможения значительно уже.

На сегодняшний день известно большое количество разнообразных алгоритмов оперативного (текущего) оценивания координат машин переменного тока [4]. Их можно разделить на пассивные и активные алгоритмы. Активные алгоритмы предполагают введение в основной спектр напряжения или тока специальных тестовых воздействий для дальнейшего анализа реакции на них электрической машины. Эти воздействия носят высокочастотный характер, например, с частотой 6-ой и более высоких гармоник. Однако, инжектируя тестовые сигналы даже с довольно малой амплитудой, приходится мириться с дополнительными потерями в электрической машине и силовом преобразователе, что, конечно, ухудшает энергоэффективность электромеханической системы в целом и увеличивает установленную мощность силовых элементов.

Системы пассивного оценивания делятся на неадаптивные и адаптивные.

Неадаптивные системы используют либо статорную модель электромагнитных процессов либо роторную. Адаптивные же системы используют две модели – эталонную и настраиваемую, что расширяет их функциональные возможности.

Для построения таких алгоритмов используется метод функций Ляпунова.

Целью диссертационной работы является построение на единой методической основе и исследование адаптивных алгоритмов вычисления неизмеряемых координат систем векторного управления электроприводами на базе АДКЗР, СДПМ и АДФР, пригодных для их применения четырехквадрантном общепромышленном ЭП переменного тока с расширенным диапазоном регулирования.

Для достижения поставленной цели в диссертации должны быть решены следующие задачи.

1. Проанализировать известные математические модели управляемых АДКЗР, СДПМ и АСМ, используемые при построении систем регулируемого ЭП, определить рациональные формы их представления при оценивании неизмеряемых координат электропривода.

2. Разработать и исследовать адаптивные алгоритмы вычисления опорного вектора потокосцеплений и частоты вращения ротора АДКЗР, СДПМ и АСМ по основным (рабочим) составляющим электрических величин, не требующих инжекции в двигатель специальных тестовых воздействий и формально работоспособные на всей плоскости механических характеристик электропривода.

3. На основе полученных результатов сформулировать обобщенную методику синтеза алгоритмов оценивания координат регулируемых электроприводов переменного тока при измерениях электрических переменных на выходе преобразователя частоты.

Научная новизна основных результатов диссертации заключается в следующем.

1. Сформулирована новая обобщенная методика структурнопараметрического синтеза адаптивных алгоритмов текущего оценивания координат и параметров электрических машин переменного тока в условиях неполных измерений, использующая в структуре вычислителя наблюдатель электромагнитных процессов полного порядка и, в отличие от известных, обеспечивающая устойчивость процессов оценивания во всех режимах работы электропривода благодаря целенаправленному заданию соотношений между элементами матрицы «стабилизирующей добавки» и матрицы весовых коэффициентов функции Ляпунова.

2. Разработаны алгоритмы оценивания опорного вектора потокосцеплений и частоты вращения ротора АДКЗР и неявнополюсного СДПМ по основным рабочим гармоникам электрических величин с возможностью адаптации к изменениям активного сопротивления статора и вычисления его текущего значения. Алгоритмы отличаются от известных работоспособностью во всех четырех квадрантах плоскости механических характеристик ЭП без инжекции в двигатель дополнительных тестовых воздействий. Сформулирована методика расчета параметров предложенных законов адаптации наблюдателей по частоте вращения и активному сопротивлению статора, учитывающая положение рабочей точки ЭП и обеспечивающая желаемое качество процессов оценивания.

Предложен новый высокоэффективный алгоритм активной предварительной идентификации параметров схемы замещения АДКЗР, совмещенный с процессом намагничивания двигателя.

3. На основе структуры наблюдателя электромагнитных процессов полного порядка разработан новый пассивный алгоритм оценивания опорного вектора потокосцеплений статора и частоты вращения ротора АСМ, работоспособный во всех четырех квадрантах плоскости механических характеристик ЭП без формальных ограничений по частоте скольжения.

Теоретическая и практическая значимость работы.

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

Предложенная методика синтеза алгоритмов вычисления координат состояния и параметров двигателей переменного тока может быть использована для построения систем управления другими, схожими по структуре динамическими объектами.

Методы исследования.

Для решения поставленных задач используются методы современной теории автоматического управления, положения теории электропривода, аналитические методы расчета, основанные на применении аппарата дифференциальных уравнений и передаточных функций. Проверка работоспособности разработанных алгоритмов осуществляется методами цифрового моделирования в пакете программ Matlab 6.5 – Simulink 5.0 и натурного эксперимента.

На защиту выносятся следующие основные результаты.

1. Алгоритм вычисления вектора потокосцеплений ротора, частоты вращения ротора и активного сопротивления статора АДКЗР. Методика расчета коэффициентов законов адаптации. Алгоритм предварительной идентификации параметров двигателя.

2. Алгоритм вычисления направления вектора потокосцеплений от постоянных магнитов, частоты вращения ротора и активного сопротивления статора неявнополюсного СДПМ. Методика расчета коэффициентов законов адаптации.

3. Алгоритм вычисления вектора потокосцеплений статора и частоты вращения ротора АСМ.

4. Обобщенная методика синтеза алгоритмов текущего оценивания координат и параметров электроприводов переменного тока на основе адаптивной системы с наблюдателем полного порядка.

Реализация результатов работы. Результаты, полученные в диссертационной работе, приняты к внедрению в системах управления асинхронными электроприводами подъемно-транспортных механизмов и синхронных электроприводов специального назначения производства ЗАО «ЭРАСИБ» (г. Новосибирск), а также используются в учебном процессе Новосибирского государственного технического университета (НГТУ).

Работа выполнена при поддержке Министерства образования и науки Российской Федерации, проект № 7.559.2011, гос. рег. номер НИР 01201255056.

Степень достоверности и апробация работы.

Достоверность изложенных в диссертации результатов и выводов подтверждается цифровым моделированием в пакете программ Matlab – Simulink и результатами натурного эксперимента.

Основные положения диссертационной работы докладывались и обсуждались на конференциях по итогам научной работы за 2009-2010 гг. и 2010гг. «Дни науки НГТУ – 2010» и «Дни науки НГТУ –2011», Новосибирск, НГТУ, в 2010 г. и 2011 г.; на Всероссийской научной конференции молодых ученых «Наука. Технологии. Инновации», Новосибирск, 2009 г.; Всероссийской научной конференции молодых ученых «Наука. Технологии. Инновации», Новосибирск, 2010 г. (первое место в конкурсе); XVII Международной научнопрактической конференции студентов и молодых учёных «Современные техника и технологии», Томск, 2011 г.; V Юбилейной международной научнотехнической конференции «Электромеханическое преобразование энергии», Томск, 2011 г.; XV научно-технической Международной конференции «Электроприводы переменного тока», Екатеринбург, 2012 г.; VII Международной (XVIII Всероссийской) конференции по автоматизированному электроприводу АЭП-2012, Иваново, 2012 г.; V Всероссийской научно-практической конференции «Автоматизированный электропривод и промышленная электроника», Новокузнецк, 2012 г.; работа была представлена на конкурсе молодежных научноисследовательских работ, проводимом Санкт-Петербургским государственным политехническим университетом в 2013 году; ХIV международной конференции молодых специалистов по микро/нанотехнологиям и электронным устройствам EDM-2013, Алтай, 2013 г.

Публикации. По теме диссертационной работы опубликованы 18 печатных работ, 4 из которых – в центральных журналах, рекомендованных списком ВАК, 2 – в сборниках научных трудов, 11 – в материалах и трудах научных конференций.

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

Структура и объем работы. Диссертационная работа состоит из введения, 5 глав, заключения, списка сокращений и условных обозначений, списка литературы и 5 приложений. Она содержит 211 стр. основного текста, 101 рисунок, 12 таблиц и библиографический список из 64 наименований.

Первая глава диссертационной работы посвящена математическому описанию электромагнитных процессов в машинах переменного тока. Приводится математическая модель АДКЗР, использующаяся при синтезе законов векторного управления в координатах состояния «токи статора – потокосцепления ротора». Рассматриваются математическая модель СДПМ, использующаяся при синтезе законов векторного управления в координатах состояния «токи статора – потокосцепления от постоянных магнитов». Приводится математическая модель управляемой по ротору АСМ, использующаяся при синтезе законов векторного управления в координатах состояния «токи ротора – потокосцепления статора». Сформулированы принципы векторного управления АДКЗР, СДПМ и АДФР, способы ориентирования вектора управляющих воздействий по направлению магнитного поля двигателя. Сделан обзор методов текущего оценивания неизмеряемых координат состояния на примере ЭП с АДКЗР. С учетом значительной степени разделения темпов электромеханических переходных процессов в бездатчиковых ЭП и желаемых алгоритмов оценивания неизмеряемых координат для их построения предложено использовать адаптивные модели динамики электромагнитных переменных двигателя.

Вторая глава посвящена методике структурного синтеза алгоритма оценивания опорного вектора потокосцеплений, частоты вращения ротора и активного сопротивления АДКЗР и анализу известных решений. Предлагаемый алгоритм представляет собой адаптивную систему, где в качестве настраиваемой модели используется наблюдатель полного порядка. Данная структура обладает рядом преимуществ – не содержит «открытых интеграторов», не имеет статической ошибки вычисления скорости, предлагаемый алгоритм работоспособен во всех точках плоскости механических характеристик электропривода. Предложены методики расчета коэффициентов адапторов предлагаемого алгоритма.

Приводятся результаты моделирования бездатчиковой системы ЭП.

Третья глава посвящена синтезу алгоритма вычисления положения, частоты вращения ротора неявнополюсного СДПМ и активного сопротивления статора. Алгоритм использует в своей структуре наблюдатель полного порядка и получен с помощью метода функций Ляпунова. Предложена методика расчета коэффициентов адапторов и определения начального положения ротора двигателя. Представлены результаты цифрового моделирования.

Четвертая глава посвящена методике структурного синтеза алгоритма оценивания опорного вектора потокосцеплений статора и частоты вращения ротора АСМ. Предложены несколько вариантов структурного построения вычислителя. Указаны преимущества и недостатки каждого варианта, даны рекомендации к их применению.

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

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

Приложения состоят из пяти разделов. Первый посвящен расчету параметров схемы замещения и номинальных данных АДКЗР, необходимых для цифрового моделирования. Здесь же приводится численный расчет регуляторов системы векторного управления АДКЗР. Во втором разделе представлен синтез системы векторного управления исследуемого СДПМ. Третий раздел посвящен расчету параметров схемы замещения и номинальных данных АДФР, необходимых для цифрового моделирования. Четвертый раздел содержит структурные математические модели в пакете Matlab-Simulink, используемые для цифрового исследования работоспособности разрабатываемых алгоритмов. В пятом разделе приведены акты об использовании результатов диссертационного исследования в производстве и учебном процессе.

ГЛАВА 1. ВЕКТОРНОЕ УПРАВЛЕНИЕ ЭЛЕКТРИЧЕСКИМИ

МАШИНАМИ ПЕРЕМЕННОГО ТОКА В УСЛОВИЯХ НЕПОЛНЫХ

ИЗМЕРЕНИЙ

1.1 Математическое моделирование машин переменного тока В данном разделе рассматриваются математические модели трехфазных АД, МДП и СДПМ, которые будут использоваться при решении задач синтеза и цифрового моделирования законов управления электроприводами.

При записи в соответствии со вторым законом Кирхгофа уравнений электрического равновесия обмоток m-фазной электрической машины двойного питания ( m 2 ), а также при последующих фазных преобразованиях, используем следующие общепринятые допущения [3]:

активные сопротивления фазных обмоток статора R s одинаковы;

геометрические оси фаз синусоидально распределенной обмотки стагде p n – число пар полюсов; исключение составтора разнесены на угол ляет случай m 2, когда данный угол равен ;

беличья клетка машины с короткозамкнутым ротором эквивалентирована m-фазной синусоидально распределенной обмоткой;

все m-фазных обмоток ротора машины имеют одинаковые активные сопротивления R r и число пар полюсов p n, а геометрический угол между их осями равен углу между осями фазных обмоток статора;

все параметры ротора приведены к обмотке статора;

эффект вытеснения токов в короткозамкнутом роторе АД пренебрежимо мал, поскольку частота токов ротора в современных частотнорегулируемых электроприводах ограниченна рабочим участком механической характеристики двигателя и незначительно сказывается на величине Rr ;

магнитная система машины линейна.

В соответствии с принятыми допущениями уравнения электрического равновесия обмоток двигателя (в общем случае с фазным ротором) имеют вид где usk, isk, sk, urk, irk, rk – мгновенные значения напряжений, токов и полных потокосцеплений k-х фазных обмоток статора и ротора машины.

В векторно-матиричной форме эти уравнения имеют вид торы-столбцы мгновенных значений фазных напряжений x u, токов x i, и потокосцеплений x обмоток статора и ротора машины соответственно.

Значительное упрощение математического описания процессов электромеханического преобразования энергии достигается путем линейных преобразований исходной системы уравнений, при этом осуществляется замена действительных переменных новыми при условии сохранения адекватности математического описания физическому объекту. В теории электрических машин доказано [5], что любая многофазная электрическая машина при условии равенства полных сопротивлений статора (ротора) путем невырожденного линейного преобразования переменных x sk, x rk может быть представлена эквивалентной двухфазной моделью. Это означает, что электромагнитный момент двигателя будет определяться только двумя компонентами векторов электромагнитных величин с индексами k 1,2.

Существует три основных способа преобразования системы координат векторов [5, 6, 7]:

1) Результирующий вектор является суммой векторов мгновенных значений фазных переменных, и направлены они по направлению соответствующих фаз или в обратную сторону (в зависимости от знака). При этом согласующий коэффициент такого преобразования будет иметь значение kc 1;

2) Проекции изображающего вектора на каждую из осей трехфазной системы координат равны истинным значениям соответствующих фазных переменных.

В этом случае энергетические характеристики двухфазной модели отличаются от энергетических характеристик реальной трехфазной машины в раз, отсюда kc ;

3) Преобразование, основанное на инвариантности мощностей, потребляемых из сети трехфазным двигателем и его двухфазной моделью. Тогда согласующий коэффициент будет равен k c.

В данной главе для преобразования реальных переменных трехфазной машины к реальным переменным эквивалентной двухфазной машины воспользуемся третьим способом преобразования – преобразованием из условия инвариантности потребляемой электромагнитной и электромеханической мощности. Следовательно, для того чтобы непосредственно из уравнения баланса мощностей преобразованной математической модели машины можно было вывести формулу момента, данные преобразования следует выполнять так, чтобы скалярные произведения любых двух векторных переменных не изменялись.

Поставленную задачу можно решить с помощью следующих преобразований:

С их помощью при математическом описании переходят от трехфазной машины к эквивалентной двухфазной [8]. Результирующий вектор какой-либо переменной может быть представлен как в трехфазной, так и в двухфазной (декартовой) системе координат, рисунок 1.1.

В трехфазной системе координат В декартовой системе координат Векторно-матричная форма:

Проекции вектора в декартовой системе координат связаны с проекциями в трехфазной системе согласно выражениям В матричном виде записи где kc – коэффициент согласования мощностей, обеспечивающий условие инвариантности мощности;

P – матрица фазного преобразования:

Формула обратного перехода В преобразовании (1.3) сохранены все исходные значения мощностных характеристик, а именно:

мгновенная мощность, потребляемая от источника питания цепи статора мгновенная мощность, потребляемая от источника питания цепи ротора мощность потерь в меди статора и ротора двигателя где. – евклидова норма вектора (корень квадратный из суммы квадратов его элементов);

и полная электромагнитная мощность, одна часть которой «расходуется» на создание магнитного поля, а вторая преобразуется в мощность механического движения (вывод аналитических зависимостей выполнен в подразделе 1.1.2).

Кроме того, при вращении двигателя оси одноименных обмоток фаз статора и ротора смещены в пространстве на электрический угол e p n, где – геометрический угол поворота вала машины относительно его положения, в котором оси одноименных обмоток статора и ротора совпадают. Поэтому для исключения периодичности магнитных связей между обмотками статора и ротора необходимо перейти от физических переменных, определяемых системой координат, жестко связанной со статором,, к расчетным переменным, соответствующим системе координат (1, 2), вращающейся в пространстве с произdk вольной угловой скоростью k рисунок 1.2.

Формулы перехода от неподвижной системы координат к вращающейся для статорных переменных имеют вид Для роторных переменных:

Модель машины в новом пространстве переменных получается в результате почленного домножения (1.1) и (1.2) слева на матрицы A k P и A k e P соответственно. Преобразованные уравнения равновесия напряжений принимают вид где e – электрическая частота вращения ротора, равная произведению геометрической угловой скорости на число пар полюсов p n, то есть D – матрица поворота на 90 электрических градусов 1.1.2 Баланс мощностей и определение электромагнитного момента Выведем формулу электромагнитного момента машины двойного питания, используя методику, описанную в [3, 10]. Дополнительно примем следующее допущение:

пренебрегаем потерями в стали, обусловленными протеканиями вихревых токов (токов Фуко) в магнитопроводе двигателя и его перемагничиванием (потерями на гистерезис), ввиду их малого влияния на динамические свойства рассматриваемых в данной работе ЭП.

Для вывода формулы момента уравнения (1.4), (1.5) почленно помножим слева на I s и I r соответственно, а затем сложим, тогда получим:

Равенству (1.6) поставим в соответствие модель баланса мощностей, согласно которой полная мощность, потребляемая по цепям статора и ротора машины, определяется следующим выражением:

где мощностям потерь в меди статора Ps и ротора Pr соответствуют два первых слагаемых из правой части (1.6); составляющая электромагнитной мощности Pэм, идущая на создание магнитного поля, соответствует третьему и четвертому слагаемым (1.6); а мощность механического движения ротора и жестко связанных с ним маховых масс Pмех, расходуемая на изменение их кинетической энергии и совершение полезной работы по определению момента сопротивления нагрузки, равна двум последним членам правой части (1.6).

Необходимо отметить, что возможность такого выделения механической мощности непосредственно связана с переходом к единой системе координат, в которой все компоненты векторов потокосцеплений не зависят от угла поворота ротора, и в этой связи Pэм не содержит составляющих, создающих электромагнитный момент [3].

В какой бы системе координат, для какой k не производился анализ баланса мощностей, при 0 всегда допустимо предположить, что где. – некоторая действительная функция времени, напряжений, токов и потокосцеплений обмоток двигателя [10].

При этом электромагнитный момент машины выражается в наиболее общей форме как Так как произведения I s Ds и Ir Dr инвариантны к углу поворота вращающейся системы координат, обобщенная формула (1.7) оказывается справедливой при любом значении, в том числе никак не связанном с частотой вращения. Так, например, полагая 1 (соответствует переходу к вращающейся системе координат, жестко связанной с ротором машины) и (неподвижная декартова система координат, жестко связанная со статором машины), из (1.9) получаем формулы электромагнитного момента, используемые в большинстве работ по векторному управлению АД, МДП и СДПМ [6,7,11, 12] где с – коэффициент, который зависит от способа фазного преобразования сисpn ; б) kc, c pn ; в) k c 1.1.3 Математическая модель АДКЗР как объекта управления Пренебрежем влиянием нелинейности кривой намагничивания на характер электромагнитных переходных процессов, то есть полагаем коэффициент Lm (главная индуктивность асинхронной машины) постоянным или меняющимся достаточно медленно. Тогда векторы потокосцеплений статора и ротора определяются по формулам [3, 13]:

где L s, L r – полные индуктивности обмоток статора и ротора соответственно;

Lm – взаимоиндуктивность обмоток статора и ротора.

Для представления математической модели электромагнитных процессов в двигателе в пространстве состояний, образованном компонентами векторов I s и r, исключим из уравнений (1.8) и (1.9) векторы где – эквивалентная индуктивность рассеяния АДКЗР, привеLe денная к статору.

Для короткозамкнутого асинхронного двигателя U r 0 – получим систему уравнений Рассмотрим математическую модель двигателя в неподвижной декартовой системе координат, причем ось совпадает с осью обмотки А статора, а ось ортогональна ей и образует с положительный угол. При использовании системы, угловая скорость вращающейся системы координат равна нулю ( k 0 ), и уравнения АДКЗР в векторной форме записи принимают вид Электромагнитный момент машины Уравнение движения одномассовой системы ЭП где M c – приведенный к валу двигателя момент сопротивления;

J – суммарный момент инерции движущихся масс, приведенный к двигателю.

Для формулирования принципов векторного управления математическую модель электромагнитных процессов АДКЗР приводят в плоской декартовой системе координат d, q, повернутой относительно неподвижной координатной системы d, q на такой угол, чтобы продольная ось d совпадала по направлению с вектором r. Поперечная ось q ортогональна оси d и образует с ней положительный угол. При использовании системы d, q скорость вращающейся системы координат равна мгновенной скорости вектора потокосцепления ротора k, и уравнения АДКЗР в скалярной форме записи принимают вид где r r – евклидова норма (модуль) вектора потокосцеплений ротора;

– постоянная времени цепи ротора АДКЗР.

Для получения математической модели АСМ как объекта векторного управления необходимо выбрать пространство состояний, в котором будет осуществляться процесс автоматического управления. Проводя аналогии с наиболее распространенной на сегодня математической моделью системы векторного управления (СВУ) АДКЗР [2], выбираем в качестве первой (векторной) координаты состояния вектор токов ротора I r – переменную, находящуюся «на стороне» силовой части ПЧ и легко доступную для непосредственного измерения; в качестве второй координаты выбираем вектор потокосцеплений статора s – наиболее «независимую» от управляющих воздействий векторную переменную, находящуюся «далеко» от силовой части ПЧ и, согласно современным тенденциям, «нежелательную» для прямого измерения [14].

С помощью выражений (1.8), (1.9) исключим из моделей статорной и роторной цепей МДП (1.4), (1.5) «лишние» переменные I s и r :

где Lr s r m – эквивалентная индуктивность рассеяния машины, приe веденная к цепи ротора;

k s m – отношение главной индуктивности к индуктивности статора.

Запишем полученную систему уравнений в следующей форме:

Система (1.12) должна рассматриваться совместно с формулой момента, которую получим также в координатах состояния I r и s :

и уравнением движения привода где J – суммарный момент инерции ротора двигателя и жестко связанных с ним маховых масс;

M c – приведенный к валу двигателя момент сопротивления нагрузки, образует математическую модель управляемой по ротору АСМ, использующуюся при синтезе законов векторного управления.

Основная задача при построении системы векторного управления АСМ заключается в выборе ориентирующего (опорного) вектора. Cамым рациональным способом является ориентация вращающейся системы координат (1, 2) по направлению вектора потокосцеплений статора. Именно такой вариант ориентирования системы координат должен обеспечить автономную линейную зависимость модуля вектора s и величины M e от соответствующих компонент вектора I r, чего нельзя сказать о других, имеющих право быть использованными обобщенных системах координат. Указанное выше преимущество позволит в дальнейшем строить СВУ АСМ как двухканальную систему подчиненного регулирования с автономными каналами управления.

Перейдем из неподвижной декартовой системы координат, во вращающуюся полеориентированную систему координат d, q. Для этого подставим k (s) в систему уравнений (1.12) где (s) (s)dt – угол поворота вектора потокосцеплений статора.

Окончательно уравнения АСМ в ориентированной по вектору s системе координат (d, q) и скалярной форме записи принимают вид [14] где Ts s – постоянная времени роторной цепи АСМ;

s s – евклидова норма (модуль вектора) потокосцеплений статора;

s ( s ) e – частота скольжения вращающейся синхронно с полем статора системы координат.

Существенным отличием синхронной машины с постоянными магнитами является отсутствие обмоток на роторе, вместо которых используются постоянные магниты. Конструктивно СДПМ выполняются либо неявнополюсными, в этом случае магнитная проводимость по осям одинакова, либо явнополюсными, в этом случае индуктивность по продольной и поперечной осям различны.

Потокосцепление статора в неподвижной системе координат, определяется согласно выражению [15] где Ld – индуктивность по продольной оси;

Lq – индуктивность по поперечной оси;

– потокосцепление от постоянных магнитов.

Переходя в систему координат d, q, вращающуюся по вектору потокосцепления от постоянных магнитов, получим следующую систему уравнений:

Уравнение электромагнитного момента Окончательно уравнения СДПМ в ориентированной по вектору системе координат (d, q) и скалярной форме записи принимают вид 1.2 Векторное управление машинами переменного тока на основе В настоящее время регулирование частоты вращения машин переменного тока осуществляется, как правило, посредством двухзыенных ПЧ на базе автономных инверторов напряжения с цифровой системой управления, выполненной на цифровых сигнальных процессорах. Для низковольтных машин используются ПЧ с двухуровневым инвертором, состоящие из выпрямителя, звена постоянного тока и трехфазного мостового инвертора напряжения. В высоковольтных ЭП используются многоуровневые инверторы напряжения, чаще всего состоящие из нескольких H-мостов, каждый из которых запитан от своего изолированного звена постоянного тока.

Для коммутации в таких преобразователях используются мощные IGBT – транзисторы. Их применение приводит к ряду особенностей с точки зрения управления и ограничения режимов работы инверторов. Прежде всего, это наличие «мертвого времени» – паузы между командой на закрытие находящегося в открытом состояние транзистора до команды на открытие другого транзистора в стойке. В момент «мертвого времени» напряжение стойки определяется знаком тока, если ток положительный ток выход стойки присоединяется к отрицательному полюсу звена через обратный диод нижнего транзистора. Если ток отрицательный, то выходная клемма присоединяется к положительному полюсу звена через обратный диод верхнего транзистора стойки. При этом выходное напряжение может отличаться от заданного значения. Это обуславливает нелинейность регулировочной характеристики ПЧ в области низких напряжений. Влияние «мертвого времени» пропорционально частоте коммутации и величине «мертвого времени». Частота коммутации определяется, как правило, максимальной выходной частотой ПЧ и требуемым уровнем пульсации тока. Минимальная величина «мертвого времени» определяется параметрами IGBT –транзистора и регламентируется в его технической документации. Для подавления нелинейности регулировочной характеристики необходимо повышать быстродействие контура тока.

Форма фазных токов при моделировании в пакете программ MatlabSimulink с учетом мертвого времени и без него представлена на рисунке 1.3.

В силу цифровой реализации системы управления имеется дискретизация по времени и по уровню. Современные цифровые сигнальные процессоры в своем составе имеют 16 битные и более высокоразрядные АЦП, что позволяет пренебречь дискретизацией по уровню.

Работа ШИМ-каналов осуществляется следующим образом. Обновление регистра сравнения (регистр, хранящий значения сигнала сравнения, используя понятия принципов вертикального управления – сигнала задания) происходит два раза за период – при достижении счетного регистра (регистр, в котором хранится значение, накопленное таймером-счетчиком, используя понятия принципов вертикального управления – опорный сигнал). Временная диаграмма работы таймера-счетчика представлена на рисунке 1.4 [16].

После формирования задания на компаратор таймер счетчика происходит запуск АЦП. По окончанию преобразования начинается расчет новых заданий, естественно, возможности процессора должны обеспечивать окончание расчета до формирования нового задания. Соответствующая временная диаграмма представлена на рисунке 1.5.

Таким образом, цифровая реализация системы управления обуславливает дискретизацию по времени и наличие внутренней временной задержки на один такт при формировании новых сигналов заданий на ШИМ.

Однако при синтезе законов и построении систем векторного управления машинами переменного тока часто пренебрегают «неидеальностями» ПЧ как регулируемого источника трехфазных симметричных напряжений. Данное допущение может считаться справедливым, поскольку быстродействие собственно широтно-импульсного преобразователя (ШИП) более чем на порядок превышает быстродействие самого внутреннего контура СВУ (контура регулирования тока). Тогда математическая модель системы регулирования напряжения может быть записана в следующем виде для управляемой по статору асинхронной машины и СДПМ, и для управляемой по ротору асинхронной машины.

В выражениях (1.15), (1.16) приняты следующие обозначения U s – сигналы задания компонент вектора напряжения статора; U r – сигналы задания компонент вектора напряжения ротора; kШИП – коэффициент усиления источника напряжения.

1.2.1 Принцип векторного управления АДКЗР Поскольку во вращающейся синхронно с магнитным полем машины системе координат все электромагнитные переменные в установившемся режиме постоянны, модель статики легко получить, приравнивая к нулю производные в первых трех уравнениях (1.13) и интегрируя предпоследнее:

где s e – частота скольжения.

Анализируя последнее уравнение системы (1.17) совместно с третьим уравнением (1.17), делаем вывод, что продольная составляющая вектора токов статора isd определяет магнитное состояние машины, характеризующееся r, а поперечный ток isq, умножаясь на текущее значение потокосцепления, создает электромагнитный момент АД. В этой связи isd называют намагничивающей составляющей вектора токов статора (намагничивающим током), а isq – активной, или моментообразующей, составляющей (активным или моментообразующим, током). На основе данного заключения базируется так называемый принцип векторного управления АДКЗР (Field-Oriented Control – FOC [2]), который в интерпретации [3] имеет следующую формулировку.

1. Законы управления частотно-регулируемым электроприводом, построенным на базе АДКЗР, могут быть синтезированы на методической основе, известной из теории ЭП постоянного тока с независимым возбуждением, если управляющие воздействия на регулируемый источник тока или напряжения формировать во вращающейся системе координат, ориентированной по вектору потокосцеплений ротора, а затем преобразовывать их в неподвижную (фазную) систему. При этом АДКЗР должен рассматриваться как двухканальный объект управления.

2. Воздействие на поперечную составляющую вектора токов статора АДКЗР должно использоваться для управления электромагнитным моментом в канале регулирования координат механического движения электропривода, подобно току якоря в ЭП постоянного тока.

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

Из первых двух уравнений (1.17) видно, что при питании двигателя от регулируемого источника напряжения поперечная составляющая вектора токов статора должна регулироваться посредством проекции вектора напряжений на ось q. Продольная составляющая тока регулируется проекцией вектора напряжений на ось d.

Согласно принципу векторного управления система регулирования АДКЗР может строиться по принципу подчиненного регулирования координат.

Структурная схема бездатчикового векторного ЭП с АДКЗР приведена на рисунке 1.6.

На рисунке 1.6 приняты следующие обозначения:

М – асинхронный двигатель, BA – датчики токов, PWM-VSI – выходное звено преобразователя частоты, выполненное на автономном инверторе напряжения с широко-импульсной модуляцией, PTd, PTq – регуляторы тока по осям d,q, РП – регулятор потокосцепления ротора, РС – регулятор скорости, ПК – преобразователи координат, осуществляющие следующие преобразования: преобразования токов статора от трехфазной системы A, B, C к неподвижной декартовой системе координат, ; преобразование токов статора от осей, к осям d, q с использованием непосредственного полеориентирования; преобразование напряжений от осей d, q к осям,.

Идентификатор предназначается для вычисления всех необходимых для реализации принципа векторного управления величин.

Расчет параметров используемой для цифрового моделирования асинхронной машины с короткозамкнутым ротором и синтез системы векторного управления представлены в Приложении А.

Модель статики СДПМ получим, полагая все производные равными нулю в (1.14):

На основе этих зависимостей формулируется следующий принцип векторного управления СДПМ.

1. С позиций теории автоматического управления СДПМ является двухканальным объектом управления, характеризующимся наличием перекрестных связей между каналами через ЭДС вращения.

2. Управляющие воздействия на регулируемый источник напряжений статора необходимо формировать во вращающейся системе координат, ориентированной по вектору потокосцеплений от постоянных магнитов, а затем преобразовывать их в неподвижные относительно статора координаты.

3. При этом канал воздействия на поперечную составляющую вектора токов статора isq следует использовать для управления электромагнитным моментом СДПМ.

4. Путем воздействия на продольную составляющую тока статора isd для явнополюсных машин возможно изменение магнитного состояния статора и электромагнитного момента двигателя.

Согласно принципу векторного управления система регулирования СДПМ может строиться по принципу подчиненного регулирования координат.

Структурная схема бездатчикового векторного ЭП с СДПМ приведена на рисунке 1.7.

Параметры используемой для цифрового моделирования СДПМ и синтез системы векторного управления представлены в Приложении Б.

Модель статики АСМ получим, полагая все производные равными нулю в (1.13):

Если пренебречь влиянием активного сопротивления статора, можно сделать вывод о том, что модуль вектора потокосцеплений статора напрямую определяется модулем вектора напряжения и частотой питания статора:

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

Следует также подчеркнуть, что при управлении АСМ главную взаимную индуктивность машины можно также с достаточной для практики степенью точности считать величиной постоянной, т.к. магнитное состояние АСМ в основном определяется параметрами питающей обмотку статора электрической сети и не изменяется в столь значительных диапазонах, как магнитный поток АДКЗР в двухзонных или энергосберегающих ЭП.

Анализируя выражения (1.18) и (1.19), можно сформулировать следующие принципы двухканального (векторного) управления АСМ с ориентацией управляющих воздействий по направлению результирующего вектора потокосцеплений статора [17].

1. С позиций теории автоматического управления АСМ является двухканальным объектом управления, характеризующимся наличием перекрестных связей между каналами через ЭДС скольжения.

2. Управляющие воздействия на регулируемый источник напряжений ротора необходимо формировать во вращающейся системе координат, ориентированной по вектору потокосцеплений статора, а затем преобразовывать их в неподвижные относительно ротора (фазные) координаты.

3. При этом канал воздействия на поперечную составляющую вектора токов ротора i rq следует использовать для управления электромагнитным моментом АСМ.

4. Магнитное состояние машины, определяемое в основном напряжением статора, может корректироваться в соответствии с требуемым режимом работы ЭП воздействием на продольную составляющую вектора токов ротора ird.

Величина последнего должна использоваться для формирования оптимальных режимов электромеханического преобразования энергии в АСМ, критерии оптимизации обсуждаются в работах [18, 19].

Из первых двух уравнений (1.18) видно, что при питании двигателя от управляемого источника напряжений поперечная составляющая вектора токов ротора должна регулироваться посредством проекции вектора напряжений на ось q. Продольная составляющая тока регулируется проекцией вектора напряжений на ось d.

На основании вышеизложенного функциональная схема бездатчиковой СВУ АСМ может быть изображена так, как это сделано на рисунке 1.8 (канал энергетической оптимизации на рисунке не показан) [14].

На рисунке 1.8: РС, РТ – регуляторы скорости и токов ротора; ПК – преобразователи координат.

В Приложении В приведен расчет параметров АДФР, используемого в данной работе.

1.2.4 Способы автоматической ориентации вращающейся системы координат и вектора управляющих воздействий по магнитному полю двигателя При построении систем векторного управления асинхронными ЭП используются два принципиально различных подхода, называемые непосредственным и косвенным ориентированием вектора управляющих воздействий по направлению магнитного поля двигателя (непосредственное и косвенное полеориентирование).

Непосредственное полеориентирование (Direct Field Oriented Control) заключается в следующем. По результатам обработки текущей информации о доступных прямым измерениям переменных (напряжениях, токах, скорости двигателя) производится оценивание компонент вектора опорного потокосцепления в неподвижной системе координат,, через которые затем определяются мгновенные значения cos и sin, используемые в преобразовании координат вида Направляющие косинус и синус ориентирующего вектора определяются, например, так:

Структура реализации непосредственным способом ориентации вращающейся системы координат по направлению изображающего вектора магнитного потока изображена на рисунке 1.9.

Косвенное ориентирование по полю (Indirect FOC, Feedforward FOC) производится без обработки информации о мгновенных токах и напряжениях АДКЗР путем вычисления оценки фазы вектора потокосцеплений ротора интегрированием суммы электрической частоты вращения и оценки частоты скольжения или сложением электрического угла поворота ротора с интегралом частоты скольжения Поскольку вектор управляющих воздействий можно формировать как в декартовой, так и в полярной вращающейся системе координат, возможно множество различных структурных решений устройства формирующего сигналы о положении полеориентированной системы координат относительно неподвижной. Зачастую, при косвенном полеориентировании, значения cos и sin определяются по.

Структура реализации косвенным способом ориентации вращающейся системы координат по направлению изображающего вектора магнитного потока изображена на рисунке 1.10.

В данном разделе дается краткое описание доступных как в отечественной, так и зарубежной литературе методов текущей пассивной идентификации неизмеряемых координат и параметров, необходимых для построения систем векторного управления машинами переменного тока. Также производится анализ преимуществ и недостатков этих алгоритмов.

1.3.1 Методы оперативного оценивания координат состояния Для алгоритма непосредственного полеориентирования необходима информация о направлении опорного вектора потокосцеплений. Использование датчиков Холла и измерительных обмоток ЭДС требует модернизации электрической машины и связано с рядом известных трудностей. Также для построения замкнутой системы регулирования частоты вращения ротора машин переменного тока необходима информация о ее текущей величине. В «датчиковом» варианте ЭП предполагается наличия датчика частоты вращения. Однако при создании новых систем ЭП перед разработчиками все чаще ставится условие отсутствия датчика частоты вращения, что диктуется соображениями повышения надежности, уменьшения стоимости и простоты обслуживания ЭП.

При модернизации существующего оборудования установка датчика частоты вращения часто попросту невозможна в силу отсутствия свободного места в действующем оборудовании. В этом случае возникает необходимость построения алгоритмов бездатчикового векторного управления. Данные алгоритмы должны функционировать, используя информацию о доступных прямому измерению координатах ЭП и управляющих воздействиях – токам и напряжениям.

На сегодняшний день имеется большое количество различных алгоритмов оценивания координат АД [4]. Их можно разделить на пассивные и активные алгоритмы. Активные алгоритмы предполагают использование каких-либо дополнительных тестовых воздействий (например, высокочастотного напряжения), которые вызывают дополнительные потери энергии в преобразователе частоты и двигателе, уменьшают перегрузочную способность привода и неизбежно влияют на ход технологического процесса. Алгоритмы пассивного оценивания предполагают использование только «рабочих» сигналов токов и напряжений на выходе инвертора. Пассивные алгоритмы в свою очередь делятся на неадаптивные и адаптивные.

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

Нейронные сети и генетические алгоритмы [20, 21] – характеризуются сложностью и неуверсальностью при переходе ЭП из одного режима в другой, требуют большое время для «обучения» системы. Отсутствуют рекомендации по выбору количества нейронов и активационной функции для применения к существенно нелинейным объектам, какими являются любые электрические машины переменного тока.

Алгоритмы, основанные на аппарате нечеткой логики [22, 23, 24 и др.] – отсутствует аналитическое обоснование применимости. При синтезе законов управления большую роли играет интуитивность и метод подбора.

Алгоритмы с инжекцией различных тестовых воздействий и работой на «паразитных» эффектах электрической машины [25, 26, 27 и др.] – для работы алгоритмов этого типа в основной спектр напряжения или тока вводятся специальные тестовые воздействия. Эти воздействия носят высокочастотный характер, например, с частотой 6-ой и более высоких гармоник. В дальнейшем алгоритм анализирует реакцию электрической машины на эти тестовые сигналы.

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

1.3.2 Адаптивное управление в бездатчиковых СВУ с оцениванием координат по основным рабочим гармоникам электрических величин Адаптивные системы используют две модели – эталонную и настраиваемую. Оценки координат состояния получаются с выхода эталонной модели, а для получения оценок параметров используется тот или иной алгоритм адаптации, например, сводящий к нулю векторное произведение оценок векторов потокосцеплений, полученных из эталонной и настраиваемой модели. Для получения таких алгоритмов используется метод функций Ляпунова. Эти алгоритмы отличаются друг от друга используемыми моделями и законами адаптации. Например, адаптивные модели с задающей моделью (АСЗМ) в качестве эталонной модели используют модель статора, а в качестве настраиваемой – модель цепи ротора. Главным недостатком алгоритмов типа АСЗМ является наличие в их структуре идеальных звеньев интегрирования или дифференцирования. Алгоритмы вычисления на основе наблюдателя полного порядка в качестве настраиваемой модели используют наблюдатель, а эталонной моделью является сама электрическая машина.

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

К примеру, ошибка в определении величины активных сопротивлений приводит к колебаниям координат ЭП «внизу» диапазона регулирования скорости (скольжения для АСМ); ошибка в определении эквивалентной индуктивности рассеяния также приводит к колебаниям координат ЭП, но уже «вверху»

диапазона регулирования скорости (скольжения); отклонение постоянной времени ротора (статора) проявляется в возникновении статической ошибки вычисления скорости [28].

Наблюдатели полного порядка (НПП) используют полную модель электромагнитных процессов электрической машины. Что позволяет оценивать большее число параметров машины. Кроме того в структуре НПП все интеграторы охвачены «естественными» отрицательными обратными связями и не нуждаются в коррекции. Впервые алгоритм данного типа был опубликован в совместных работах японских исследователей H. Kubota и K. Matsuse [29, 30] и в настоящее время активно ими развивается. Однако результаты исследования предложенного алгоритма методом цифрового моделирования показали его недостатки – неустойчивость алгоритма адаптации при низких частотах вращения и в генераторных режимах.

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

1. В первой главе диссертации приведено математическое описание асинхронной машины двойного питания как объекта управления в частотнорегулируемом электроприводе. На его основе составлены математические модели АДКЗР, управляемой по ротору АСМ и СДПМ.

2. На общей методической основе сформулированы принципы векторного управления АДКЗР в координатах ток статора – потокосцепление ротора, АСМ – в координатах ток ротора – потокосцепление статора, а также СДПМ. Определены управляющие воздействия и соответствующие им регулируемые координаты, описаны способы ориентирования вектора управляющих воздействий по направлению магнитного поля двигателя.

3. На примере ЭП с АДКЗР сделан обзор методов вычисления неизмеряемых координат состояния в бездатчиковых системах регулирования скорости. В качестве технически рационального для общепромышленных ЭП выбрана структура алгоритма адаптивного НПП.

ГЛАВА 2. СИНТЕЗ И ИССЛЕДОВАНИЕ АЛГОРИТМА ВЫЧИСЛЕНИЯ

ЧАСТОТЫ ВРАЩЕНИЯ И ОПОРНОГО ВЕКТОРА

ПОТОКОСЦЕПЛЕНИЙ АДКЗР, МАЛОЧУВСТВИТЕЛЬНОГО К

ДРЕЙФУ АКТИВНОГО СОПРОТИВЛЕНИЯ СТАТОРА

2.1 Обзор адаптивных алгоритмов вычисления координат АДКЗР На сегодняшний день наиболее распространенными как в теоретических, так и в инженерных приложениях системами пассивного оценивания частоты вращения ротора АД [31] являются алгоритмы типа MRAS – Model Reference Adaptive System (адаптивная система с задающей или эталонной моделью – АСЗМ), AFOO – Adaptive Full Order Observer (адаптивный наблюдатель полного порядка – АНПП), EKF – Extended Kalman Filter (расширенный фильтр Калмана). Выходными сигналами всех этих алгоритмов являются оценки вектора потокосцеплений и электрической частоты вращения ротора асинхронной машины.

Идентификаторы типа АСЗМ (MRAS) впервые предложенные Сolin Schauder содержат в своей структуре две различные модели, которые вычисляют две оценки одной и той же (в общем случае – векторной) переменной состояния, например вектора потокосцеплений ротора. Одна модель, не содержащая неопределённого параметра, является задающей, и вычисленная с её помощью переменная считается эталоном. Другая модель, зависящая от неопределённого параметра, является настраиваемой, причем её адаптация осуществляется путём изменения оценки неопределённой переменной, подлежащей вычислению.

Рассмотрим процедуру синтеза алгоритма MRAS, сводящего к нулю векторное произведение потокосцеплений ротора, полученных из уравнений модели цепи статора АД, которая принимается эталонной, и модели цепи ротора АД, которая является настраиваемой. Структурная схема алгоритма представлена на рисунке 2.1 [28].

Для построения идентификатора частоты вращения ротора на основе метода АСЗМ, используется математическая модель АДКЗР в координатах состояния в неподвижной декартовой системе координат,, жестко связанной со статором, имеющая вид Так как модель цепи статора содержит только одну неизвестную величину – вектор потокосцеплений ротора, ее можно использовать для его нахождения:

Полагая, что имеется точная информация о параметрах двигателя, модель цепи статора принимается в качестве эталонной модели:

Структурная схема модели цепи статора представлена на рисунке 2.2.

Модель цепи ротора содержит настраиваемый параметр – электрическую частоту вращения вала машины и является настраиваемой. Рассмотрим совместно математическую модель цепи ротора АДКЗР и настраиваемую модель:

Вычитая из уравнений объекта наблюдения (модели цепи ротора АДКЗР) уравнения настраиваемой модели получим уравнения в отклонениях:

ра;

e e – отклонение оценки скорости (настраиваемый параметр).

Так как адаптивные системы являются нелинейными, для синтеза алгоритмов идентификации применяется прямой метод Ляпунова [32, 33].

Зададимся следующей функцией Ляпунова:

где – некоторая положительная константа.

Производная функции Ляпунова в силу уравнений объекта наблюдения (2.1) имеет вид Согласно теореме Ляпунова для устойчивости алгоритма оценивания достаточно обеспечить положительную определенность функции Ляпунова и отрицательную определенность ее производной.

Первое требование обеспечено выбором функции Ляпунова:

Для обеспечения отрицательной определенности ее производной достаточно обнулить два ее последних слагаемых:

Следовательно, алгоритм адаптации принимает вид:

При этом производная функции Ляпунова будет отрицательно определенной функцией С целью повышения гибкости в выборе показателей качества процессов дополним (2.2) пропорциональной составляющей и получим закон оценивания электрической частоты вращения ротора вида Известно, что прямой метод Ляпунова позволяет получить структуру закона адаптации, однако не позволяет найти численные значения коэффициентов регулятора – адаптора. Для их определения предлагается следующий подход.

Рассмотрим векторное произведение векторов этал – модуль оценки вектора потокосцепления ротора, полученный по где r эталонной модели;

н – модуль оценки вектора потокосцепления ротора, полученный по наr страиваемой модели;

е s dt – разница углов между векторами потокосцеплений.

Полагая, что этал н rном и что – малый угол, получим структурную схему, представленную на рисунке 2.3.

Настраивая замкнутый контур как звено второго порядка с собственной частотой АСЗМ и коэффициентом формы ААСЗМ, получим выражения для коэффициентов регулятора:

Известно, что модель цепи статора является нейтрально устойчивой, так как имеющиеся в ее структуре два интегратора, не охваченные обратной связью, обуславливают нулевой корень характеристического уравнения кратности два. Для того чтобы сделать идентификатор устойчивым асимптотически, необходимо стабилизировать модель цепи статора (МЦС).

Стабилизацию МЦС осуществим охватом интегратора отрицательной обратной связью и введением в нее ПИ-регулятора [31], который обеспечит нечувствительность идентификатора к смещениям нулей сигналов с датчиков тока и напряжения. Соответствующая структурная схема представлена на рисунке 2.4.

Передаточная функция контура после стабилизации имеет вид Настраивая знаменатель передаточной функции замкнутого контур на характерестический полином второго порядка с собственной частотой СМЦС и коэффициентом формы АСМЦС, получим выражения для коэффициентов регулятора:

ЛАЧХ замкнутого контура стабилизации модели цепи статора представлена на рисунке 2.5.

Однако наличие в ЛАЧХ замкнутого контура участка с наклоном +20дБ/дек ограничивает диапазон регулирования скорости. Поэтому собственная частота контура стабилизации должна выбираться из условия:

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

Структурная схема синтезированного алгоритма оценивания опорного вектора потокосцеплений ротора и электрической частоты вращения ротора, основанного на АСЗМ, представлена на рисунке 2.6.

Цифровое моделирование вычислителя типа АСЗМ проводилось в пакете программ MATLAB/Simulink в составе системы векторного управления двигателем 4А225М4У3. Модели алгоритма представлены в Приложении Г.

На рисунке 2.7 представлены графики переходных процессов в системе векторного управления с вычислителем АСЗМ по механической частоте вращения ротора, заданию на токи в системе координат d, q, заданиям фазных токов и напряжений при работе в двигательном режиме. Электропривод отрабатывает следующие режимы – намагничивание двигателя, разгон без нагрузки на частоту 3ном 4, наброс номинальной нагрузки, сброс нагрузки.

Известно, что адаптивные алгоритмы (в том числе и MRAS) чувствительны к отклонениям параметров АДКЗР, которыми оперирует вычислитель от их реальных величин. Эти отклонения могут быть связаны как с ошибочной информацией о параметрах, полученной на этапе предварительной их идентификации, так и их дрейфом в процессе работы электропривода, в частности изменением активных сопротивлений статора и ротора в результате нагрева машины. Ошибочная информация о параметрах машины приводит к возникновению ошибок оценивания и может стать причиной неустойчивости СВУ.

Переходный процесс по электрической частоте вращения ротора и ее оценке внизу диапазона регулирования в двигательном режиме под номинальной нагрузкой на частоте ном 20 при ступенчатом отклонении сопротивления статора на 20% представлен на рисунке 2.8.

Аналогичные переходные процессы при работе электропривода в генераторном режиме представлены на рисунке 2.9.

«Развал» алгоритма оценивания и всей СВУ внизу диапазона на частоте ном 50 под номинальной нагрузкой при отклонении активного сопротивления статора на 30% показан на рисунке 2.10.

Влияние отклонения эквивалентной индуктивности на 20% вверху диапазона регулирования представлено на рисунке 2.11.

На рисунке 2.12 представлена реакция алгоритма АСЗМ на отклонение постоянной времени ротора на 20%.

Таким образом, алгоритм типа АСЗМ имеет область неустойчивой работы на малых частотах питания двигателя, что связано с необходимостью стабилизации модели цепи статора. Также он чувствителен к отклонениям параметров АДКЗР от их реальных величин, особенно к дрейфу активного сопротивления статора. Так как в качестве эталонной модели принята модель статора, а в роторную модель активное сопротивление статора не входит, текущая его пассивная идентификация в рамках рассматриваемого подхода невозможна. Однако следует отметить относительную простоту рассмотренного алгоритма, что обуславливает его широкое применение.

Рассмотрим процедуру синтеза и недостатки алгоритма, предложенного японским исследователем Hisao Kubota [29, 30]. Данный алгоритм считается наилучшим по ряду показателей [34] и в качестве настраиваемой модели использует полную модель электромагнитных процессов асинхронного двигателя в декартовой неподвижной системе координат,. В качестве эталонной модели используется сам асинхронный электродвигатель. Выходными координатами модели НПП являются оценки вектора тока статора и потокосцеплений ротора. Далее вычисляется вектор токовой ошибки – разность между оценкой тока статора, полученной из модели АДКЗР (или просто АД), и реальным (измеренным) током в двигателе. Идея метода заключается в сведении к нулю произведения вектора токовой ошибки на оценку вектора потокосцеплений ротора, посредством интегрального или пропорционально-интегрального регулятора, формирующего оценку электрической частоты вращения ротора. Структурная схема алгоритма представлена на рисунке 2.13.

Математическая модель объекта наблюдения – электромагнитных процессов АД – в неподвижной декартовой системе координат в обозначениях [29] имеет вид где A11, A12, A 21, A 22 – матрицы вида:

Данная модель принимается в качестве эталонной.

Модель наблюдателя полного порядка (построенного по структуре наблюдателя Люенбергера):

где ^ – знак оценки соответствующих величин;

L – матрица «стабилизирующей добавки», dim L 4 2, Вычитая из уравнения наблюдателя (2.3) уравнение объекта наблюдения (2.4), можно получить уравнения для ошибок:

где i I s I s – токовая ошибка;

r r – ошибка потокосцепления;

e e e – ошибка идентификации электрической частоты вращения ротора.

Авторами [29] предлагается использовать следующую положительно– определенную функцию Ляпунова:

При этом ее производная согласно [26] будет иметь вид:

Авторами [29] предлагается использовать нулевую «стабилизирующую»

добавку Для обеспечения отрицательной определенности производной функции Ляпунова необходимо обнулить сумму второго и третьего слагаемых (2.5):

Тогда закон адаптации электрической частоты вращения имеет вид Для улучшения динамики наблюдения [29, 30] (2.6) дополняют пропорциональной частью. Следовательно, алгоритм идентификации электрической частоты вращения ротора (как параметра системы (2.3)) имеет вид Цифровое моделирование идентификатора AFOO производилось в пакете программ MATLAB/Simulink в составе системы векторного управления двигателем 4А225М4У3, модели идентификатора представлены в Приложении Г.

На рисунке 2.14 представлены графики переходных процессов в системе векторного управления, замкнутой посредством алгоритма НПП при нулевой стабилизирующей добавке L 0, по механической частоте ротора, заданию на токи в системе координат d, q, заданиям фазных токов и напряжений при работе в двигательном режиме. Электропривод отрабатывает следующие режимы – намагничивание двигателя, разгон без нагрузки на частоту 3ном 4, наброс номинальной нагрузки, сброс нагрузки.

Однако известно [33], что данный алгоритм обладает областью неустойчивой работы в генераторных режимах на малых скоростях.

«Развал» алгоритма в генераторном режиме на частоте 3ном 100 представлен на рисунке 2.15.

Для исследования области неустойчивости воспользуемся подходом, предложенным в работе [35].

Рассмотрим математическую модель асинхронного двигателя во вращающейся системе координат, направленной по вектору потокосцепления ротора, Уравнения наблюдателя в той же системе координат при L 0 с учетом алгоритма адаптации (2.6) могут быть представлены в виде Вычитая из уравнения наблюдателя уравнение объекта, получим уравнения для ошибок:

Соответствующий характеристический полином будет иметь вид где коэффициенты характеристического полинома:

Для анализа устойчивости воспользуемся алгебраическим критерием Гурвица. Матрица Гурвица для характеристического полинома пятого порядка имеет вид [36] Определители Гурвица:

Все определители Гурвица будут положительны и система устойчива в силу выражений для коэффициентов характеристического полинома при условии:

Найдем границы области, в которой это условие не выполняется, в плоскости механических характеристик ЭП. Для этого необходимо решить уравнение Переходя на плоскость механических характеристик и используя равенство получим следующее выражение для коэффициента a5 :

где 1 M e, r, 2 M e, r – границы области устойчивости системы:

В верхней полуплоскости механических характеристик область неустойчивой работы представлена на рисунке 2.16.

Соответствующие коэффиценты наклона границ области неустойчивости:

Таким образом, область неустойчивой работы алгоритма находится в генераторных режимах работы на малых скоростях и граничит с характеристикой динамического торможения двигателя. В работе [35] для обеспечения работы привода в этой области предлагается организовать контур зависимого ослабления поля, тем самым сдвинуть область неустойчивости выше требуемой рабочей точки, как показано на рисунке 2.17.

При этом наибольший коэффициент ослабления поля, для исследуемого двигателя составляет:

Переходные процессы при ослаблении поля в области неустойчивости алгоритма представлены на рисунке 2.18. Электропривод отрабатывает следующие режимы – намагничивание двигателя, разгон без нагрузки на частоту ном/100 наброс номинальной генераторной нагрузки, последующий разгон на частоты ном/50 и 3ном/100. Как видно, процессы стали устойчивыми.

Как видно, процессы стали устойчивыми. Однако такой подход не всегда приемлем. Во-первых, требуется организация дополнительного контура ослабления поля, который должен формировать задание на поток машины в зависимости от нахождения рабочей точки на плоскости механических характеристик.

Во-вторых, ослабление поля приводит к увеличению тока двигателя, а следовательо к уменьшению перегрузочной способности привода и увеличению потерь во всех его силовых элементах.

Наличие области неустойчивой работы данного алгоритма вызвано невыполнением условий отрицательной определенности производной функции Ляпунова. Во-первых, допущена ошибка, либо не учтено одно слагаемое, входящее в производную функции Ляпунова, предложенной в [26], которая реально имеет вид [37] Во-вторых, первое слагаемое производной функции Ляпунова не является строго отрицательным при L 0.

Таким образом, известный алгоритм AFOO имеет область неустойчивой работы на малых частотах питания АД, в генераторных режимах. Кроме того, в доступных публикациях нет методики выбора коэффициентов закона идентификации электрической частоты вращения ротора. Также АНПП чувствителен к отклонениям параметров АД от реальных их величин, особенно к дрейфу активного сопротивления статора. Предлагаемый в [30] алгоритм идентификации активного сопротивления статора не дает рекомендаций по выбору коэффициента регулятора и дополнительно увеличивает область неустойчивой работы в генераторных режимах.

Рассмотрим алгоритм финского исследователя Marko Hinkkanen [38, 39, 40]. Данный алгоритм отличается от алгоритма Hisao Kubota [29] наличием ненулевой стабилизирующей добавки вида где – коэффициент, зависящий от текущей частоты вращения согласно формуле где, и – некоторые положительные константы.

Для идентификации частоты вращения авторами [40] предлагается использовать следующий алгоритм:

где – угол стабилизации генераторных режимов внизу диапазона регулирования.

Для стабилизации генераторных режимов рекомендуется следующий алгоритм формирования угла:

где max – максимальное значение угла стабилизации;

– граничная частота для включения режима стабилизации.

Однако в опубликованных работах [38, 39, 40] авторами не дается методики и рекомендаций по выбору величин ',, max,, а также коэффициентов адаптора частоты kи и k п. Для используемого в этих работах двигателя мощностью 2,2 кВт авторами предлагаются следующие значения величин [40].

Параметры двигателя представлены в таблице 2.1.

Таблица 2. Приведенная индуктивность рассеяния 0,0209 Гн Цифровое моделирование алгоритма Hinkkanen проводилось в пакете программ MATLAB/Simulink в составе системы векторного управления, модели идентификатора представлены в Приложении Г.

На рисунке 2.19 представлены графики переходных процессов в системе векторного управления замкнутой по идентификатору Hinkkanen, по механической частоте ротора, заданию на токи в системе координат d, q, заданиям фазных токов и напряжений при работе в двигательном режиме. Электропривод отрабатывает следующие режимы – намагничивание двигателя, разгон без нагрузки на частоту 3ном 4, наброс номинальной нагрузки, сброс нагрузки.

Работа электропривода на частоте ном 20 представлена на рисунке 2.20.

Как видно электропривод обеспечивает приемлемое качество переходных процессов и астатизм по возмущению.

Попытки же применить приведенные выше параметры к ЭП существенно большей мощности не увенчались успехом: по-прежнему, СВУ была неустойчивой на малых скоростях в режимах генераторного торможения.

Таким образом, все известные алгоритмы вычисления неизменяемых координат состояния АДКЗР имеют области неустойчивой работы на малых частотах питания при генераторной нагрузке. Также авторами не предлагаются подходы по расчету коэффициентов законов адаптации.

Поэтому в настоящее время актуальна задача синтеза алгоритма вычисления координат состояния АД, устойчивого во всей плоскости механических характеристик.

2.2 Предлагаемый метод структурного синтеза алгоритмов вычисления оценок координат АДКЗР на основе адаптивного наблюдателя состояния Электромагнитные процессы в АД могут быть описаны следующими уравнениями [41] (совпадают с (2.3)):

C E 0 – матрица выхода;

переменных;

A – собственная матрица объекта, B Le – матрица управления;

– единичная матрица;

R r, R s, k r, Le, Tr, Lm – параметры двигателя;

e – электрическая частота вращения ротора двигателя.

Математическая модель НПП в неподвижной системе координат [42]:

где ^ – знак оценки соответствующих величин;

L T – матрица «стабилизирующей добавки»;

A – собственная матрица наблюдателя, имеет следующий вид:

Здесь все параметры двигателя, кроме активных сопротивлений статора и ротора, считаются известными достаточно точно по результатам их предварительной идентификации.

Вычитая из уравнений наблюдателя (2.8) уравнения объекта (2.7), получим уравнение динамики в отклонениях:

где x x – вектор ошибок;

s r T – вектор отклонений параметров;

e e – отклонение электрической частоты вращения двигателя;

s Rs Rs – отклонение активного сопротивления статора;

r Rr Rr – отклонение активного сопротивления ротор;

A – матрица влияния вектора отклонений параметров на вектор ошибок:

Согласно методу Красовского зададимся следующей функцией – кандидатом в функции Ляпунова [32]:

где H, H – симметричные матрицы весовых коэффициентов.

Предлагается использовать матрицы весовых коэффициентов вида где hi, h, h,, s, r – некоторые положительные константы.

Производная функции Ляпунова:

Алгоритм адаптации сформируем в виде где G – искомая матрица адаптации.

Тогда уравнения динамики НПП в отклонениях записываются как Производная функции Ляпунова при подстановке (2.9) в (2.10) примет вид квадратичной формы

H A LC A LCT H

Как известно, для анализа положительной или отрицательной определенности квадратичной формы можно использовать критерий Сильвестра [32]. Согласно нему для обеспечения отрицательной определенности производной функции Ляпунова при выполнении приведенных ниже дополнительных условий достаточно равенства при этом из производной функции Ляпунова будут исключены отклонения параметров:

V T H A LC A LCT H

Матрица этой квадратичной формы H A LC A LC H определяется как где S11, S12, S 21, S 22 – следующие матрицы:

Развернуто матрица адаптации имеет вид Соответствующие законы адаптации:

Таким образом, к весовым коэффициентам функции Ляпунова и матрице «стабилизирующей» добавки будем предъявлять следующие требования.

1) Обеспечение положительной определенности функции Ляпунова, то есть положительной определенности квадратичной формы (2.9).

2) Весовые коэффициенты и матрицы стабилизирующей добавки должны обеспечить отрицательную определенность производной функции Ляпунова (2.10) и, в частности, квадратичной формы (2.11). Следует подчеркнуть, что использование авторами [29, 43] нулевой «стабилизирующей»

матрицы это условие не обеспечивает.

3) Выбор соотношений весовых коэффициентов должен обеспечивать реализуемость алгоритмов адаптации в условиях отсутствия датчиков магнитного состояния машины.

2.2.2 Алгоритм текущей идентификации электрической частоты вращения Для исключения из закона адаптации частоты вращения ротора нереализуемых в рамках бездатчикового управления слагаемых, включающих в себя Dr, предлагается принять следующее соотношение:

Кроме того, для отрицательной определенности производной функции Ляпунова (глобальной устойчивости алгоритма адаптации) при принятом соотношении (2.12) автором предлагается принять соотношение и использовать стабилизирующую добавку При этом производная функции Ляпунова для малых является отрицательно определенной функцией.

При принятых соотношениях закон адаптации НПП по электрической частоте вращения ротора имеет вид Дополняя (2.14) пропорциональной частью, получим пропорциональноинтегральный закон формирования оценки электрической частоты вращения двигателя Известно, что прямой метод Ляпунова позволяет получить структуру закона адаптации, однако не позволяет найти численные значения коэффициентов регулятора – адаптора. Предлагается следующий подход по их определению [41].

Рассмотрим векторное произведение векторов, входящее в (2.15):

где I s – модуль вектора оценки тока статора;

r – модуль вектора оценки потокосцепления ротора;

I s – модуль вектора токов статора;

– угол между вектором оценки токов статора и вектором оценки потокосцепления ротора, см. рисунок 2.21;

i – угол между вектором тока статора и вектором его оценки.

Принимая допущение I s I s, и что i – малые углы, получим Для электропривода, работающего в первой зоне регулирования, можно принять что r rном. В этом случае структурная схема линеаризованного НПП имеет вид, приведенный на рисунке 2.22.

Настраивая характерестический полином замкнутого контура на нормированый полином второго порядка с собственной частотой НПП и коэффициентом формы AНПП, получим выражения для коэффициентов регулятора – адаптора:

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

Более того, данный подход выходит за рамки пассивного оценивания, так как любой сигнал инжекции будет влиять на ход технологического процесса и вызовет дополнительные потери в силовых элементах ПЧ и АД. По этим причинам одновременная идентификация сопротивления роторной цепи и частоты вращения ротора в данной работе не рассматривается.



Pages:     || 2 | 3 |


Похожие работы:

«ИЗ ФОНДОВ РОССИЙСКОЙ ГОСУДАРСТВЕННОЙ БИБЛИОТЕКИ Лунин, Николай Николаевич Мошенничество по уголовному законодательству России: уголовно­правовая характеристика и квалификация Москва Российская государственная библиотека diss.rsl.ru 2006 Лунин, Николай Николаевич.    Мошенничество по уголовному законодательству России: уголовно­правовая характеристика и квалификация  [Электронный ресурс] : Дис. . канд. юрид. наук  : 12.00.08. ­ Орел: РГБ, 2006. ­ (Из фондов Российской...»

«Еремина Мария Геннадьевна Оценка воздействия иммунозависимых дерматозов на качество жизни лиц трудоспособного возраста 14.02.05 – социология медицины Диссертация на соискание ученой степени кандидата медицинских наук Научный руководитель доктор медицинских наук,...»

«УДК 524.5-7; 52-17 Хоперсков Сергей Александрович ЭВОЛЮЦИЯ ДИСКОВЫХ ГАЛАКТИК: ИССЛЕДОВАНИЕ ИЕРАРХИИ СТРУКТУР 01.03.02 – Астрофизика и звездная астрономия ДИССЕРТАЦИЯ на соискание ученой степени кандидата физико-математических наук Научный руководитель д. ф.-м. н., проф., чл. корр. РАН Шустов Борис Михайлович Москва – Содержание Введение........................»

«Джаксумбаева Ольга Ильинична ПОДСИСТЕМА ПОДДЕРЖКИ ПРИНЯТИЯ РЕШЕНИЙ ДЛЯ ПЛАНИРОВАНИЯ СОЦИАЛЬНЫХ ВЫПЛАТ ИНФОРМАЦИОННОЙ СИСТЕМЫ ТИПА ЭЛЕКТРОННЫЙ СОЦИАЛЬНЫЙ РЕГИСТР НАСЕЛЕНИЯ Специальность 08.00.13 – Математические и инструментальные методы экономики Диссертация на соискание ученой степени кандидата экономических наук Научный руководитель : Д.э.н., доцент Халин Владимир Георгиевич 2 Санкт – Петербург - 2014 Оглавление ВВЕДЕНИЕ ГЛАВА...»

«Белякова Анастасия Александровна Холодноплазменный хирургический метод лечения хронического тонзиллита 14.01.03 — болезни уха, горла и носа Диссертация на соискание ученой степени кандидата медицинских наук Научный руководитель : член-корр. РАН, доктор медицинских наук, профессор Г.З. Пискунов Москва– СОДЕРЖАНИЕ СПИСОК СОКРАЩЕНИЙ ВВЕДЕНИЕ ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ...»

«АСАНБАЕВА АЛУА ЕРМЕКОВНА Историческая преемственность в духовной культуре казахского общества XV-XVIII вв.: теоретико-методологические проблемы 07.00.09 – Историография, источниковедение, методы исторического исследования Диссертация на соискание ученой степени кандидата исторических наук Научный руководитель доктор исторических наук, доцент Оразбаева А.И. Республика Казахстан Алматы, СОДЕРЖАНИЕ ВВЕДЕНИЕ.. 1...»

«Молчанова Ирина Сергеевна ПОНЯТИЕ РИСКА В ДОГОВОРЕ ИМУЩЕСТВЕННОГО СТРАХОВАНИЯ 12.00.03 – Гражданское право; предпринимательское право; семейное право; международное частное право Диссертация на соискание ученой степени кандидата юридических наук Научный руководитель : Лебедев Константин Константинович доцент, кандидат юридических наук, доцент кафедры коммерческого права Санкт-Петербургского государственного...»

«Вакуленко Андрей Святославович ОБЩЕСТВЕННОЕ МНЕНИЕ В СОЦИАЛЬНО–ИСТОРИЧЕСКОМ ПРОЦЕССЕ 09.00.11 – социальная философия Диссертация на соискание ученой степени кандидата философских наук Научный руководитель : доктор философских наук, профессор Зорин Александр Львович Краснодар – 2014 Содержание ВВЕДЕНИЕ.. ГЛАВА Теоретико–методологические основы изучения I. общественного мнения.. 1.1. Полисемантичность...»

«ХАХАЛИНА АНАСТАСИЯ АЛЕКСАНДРОВНА МОЛЕКУЛЯРНО-ГЕНЕТИЧЕСКИЙ АНАЛИЗ МУТАЦИЙ В ГЕНАХ gyrA и gyrB, СВЯЗАННЫХ С УСТОЙЧИВОСТЬЮ MYCOBACTERIUM TUBERCULOSIS К ФТОРХИНОЛОНАМ 03.02.03 – микробиология Диссертация на соискание ученой степени кандидата биологических наук Научные руководители: кандидат медицинских...»

«Робенкова Татьяна Викторовна ПСИХОТИПОЛОГИЧЕСКИЕ АСПЕКТЫ АДАПТАЦИИ СТУДЕНТОВ КОЛЛЕДЖА 03.00.13 – физиология Диссертация на соискание ученой степени кандидата медицинских наук Научный руководитель : доктор биологических наук, профессор В.Н. Васильев Томск - 2003 ОГЛАВЛЕНИЕ. ВВЕДЕНИЕ..7 ГЛАВА 1. ОБЗОР ЛИТЕРАТУРЫ.. 1.1.Современный подход к проблеме адаптации студентов. 1.1.1. Роль стресса в...»

«Абызов Алексей Александрович ОБЕСПЕЧЕНИЕ БЕЗОТКАЗНОСТИ ЭЛЕМЕНТОВ ХОДОВЫХ СИСТЕМ БЫСТРОХОДНЫХ ГУСЕНИЧНЫХ МАШИН ПРИ ПРОЕКТИРОВАНИИ НА ОСНОВЕ МОДЕЛИРОВАНИЯ ПРОЦЕССОВ ЭКСПЛУАТАЦИИ И ФОРМИРОВАНИЯ ОТКАЗОВ Специальности: 05.05.03 – Колесные и гусеничные машины 01.02.06 – Динамика, прочность...»

«БОГУШ Глеб Ильич Коррупция и международное сотрудничество в борьбе с ней Специальность 12. 00. 08 – уголовное право и криминология; уголовно-исполнительное право Диссертация на соискание ученой степени кандидата юридических наук Научный руководитель – Лауреат государственной премии СССР, Заслуженный деятель науки Российской Федерации, доктор юридических наук, профессор Н. Ф. Кузнецова Москва -...»

«АЗИНА Ольга Александровна МИФОДИЗАЙН КАК СРЕДСТВО ФОРМИРОВАНИЯ ИМИДЖА СОВРЕМЕННОЙ РОССИИ В ЕВРОПЕЙСКОМ ОБЩЕСТВЕННОМ МНЕНИИ Специальность 22.00.04. – Социальная структура, социальные институты и процессы ДИССЕРТАЦИЯ на соискание ученой степени кандидата социологических наук Научный руководитель : Мамедов А.К. доктор социологических наук, профессор Москва - Содержание Введение Глава 1....»

«НАСАН-ОЧИР ЭРДЭНЭ-ОЧИР ВОЕННОЕ ДЕЛО ДРЕВНИХ КОЧЕВНИКОВ МОНГОЛИИ (II тыс. до н.э. – III век до н.э.) Специальность 07.00.06 - археология Диссертация на соискание ученой степени кандидата исторических наук Научный руководитель – доктор исторических наук, профессор Ю.С. Худяков Новосибирск – ОГЛАВЛЕНИЕ Стр....»

«АБДУХАНОВА НАТАЛЬЯ ГЕННАДЬЕВНА ОРГАНИЗАЦИОННО-ЭКОНОМИЧЕСКИЙ МЕХАНИЗМ ИСПОЛЬЗОВАНИЯ ЛИЗИНГОВЫХ ОПЕРАЦИЙ В ЖИЛИЩНОКОММУНАЛЬНОМ КОМПЛЕКСЕ Специальность 08.00.05 – Экономика и управление народным хозяйством: экономика, организация и управление предприятиями, отраслями, комплексами (строительство) Диссертация на соискание ученой степени кандидата экономических...»

«А.И.Клюкина МЕТОДОЛОГИЯ И ПРАКТИКА РЕАЛИЗАЦИИ ПЕДАГОГИЧЕСКОГО ПОТЕНЦИАЛА МУЗЕЯ (На примере музеев естественнонаучного профиля) Специальность 13.00.05 – Теория, методика и организация социально-культурной деятельности Диссертация на соискание ученой степени доктора педагогических наук Санкт – Петербург 2014 год 1 СОДЕРЖАНИЕ Введение Стр. Актуальность и степень разработанности проблемы. Цель. Задачи. Объект. Предмет исследования. Методология, методика, организация, база...»

«БОГОПОЛЬСКИЙ Павел Майорович ИСТОРИЯ РЕКОНСТРУКТИВНОЙ ХИРУРГИИ ПИЩЕВОДА В РОССИИ Диссертация на соискание ученой степени доктора медицинских наук 07.00.10 – История науки и техники (медицинские науки) Научные консультанты: д.м.н. С.А. Кабанова д.м.н. проф. М.М. Абакумов Москва – 2014 г. ОГЛАВЛЕНИЕ Страницы Введение 5– Глава I. Исследования по истории развития...»

«Бардаченко Алексей Николаевич КРИМИНАЛИСТИЧЕСКОЕ ИССЛЕДОВАНИЕ СЛЕДОВ ТЕРМИЧЕСКОЙ РЕЗКИ НА ПРЕГРАДАХ Специальность 12.00.12 – криминалистика; судебно-экспертная деятельность; оперативно-розыскная деятельность Диссертация на соискание ученой степени кандидата юридических наук Научный руководитель : доктор юридических наук, профессор Ручкин Виталий Анатольевич Волгоград - СОДЕРЖАНИЕ ВВЕДЕНИЕ.....»

«УДК: 612.015.13:611.33-018.73 Животова Елена Юрьевна УЧАСТИЕ РЕГУЛЯТОРНЫХ ПЕПТИДОВ В ПОДДЕРЖАНИИ ТКАНЕВОГО ГОМЕОСТАЗА СЛИЗИСТОЙ ОБОЛОЧКИ ЖЕЛУДКА 03.03.04 – клеточная биология, цитология, гистология Диссертация на соискание ученой степени доктора медицинских наук Научный консультант...»

«Газиева Ирина Александровна ИММУНОПАТОГЕНЕТИЧЕСКИЕ МЕХАНИЗМЫ ФОРМИРОВАНИЯ ПЛАЦЕНТАРНОЙ НЕДОСТАТОЧНОСТИ И РАННИХ РЕПРОДУКТИВНЫХ ПОТЕРЬ 14.03.09 – Клиническая иммунология, аллергология Диссертация на соискание ученой степени доктора биологических наук Научный консультант : Чистякова Гузель Нуховна доктор...»






 
2014 www.av.disus.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Программы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.