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

ISSN08б8-588б

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2001, том 11, № 3, с. 76-86 ОРИГИНАЛЬНЫЕ СТАТЬИ

УДК 538.65: 621.316.73 © Г. А. Сапожников, С. В. Богословский

ОПТИМИЗАЦИЯ УПРАВЛЕНИЯ МАГНИТНЫМ ПОДВЕСОМ В ДАТЧИКАХ

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

ВВЕДЕНИЕ

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

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

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

КОНСТРУКТИВНЫЕ СХЕМЫ МАГНИТНЫХ ПОДВЕСОВ ДЛЯ ИЗМЕРИТЕЛЬНЫХ ПРИБОРОВ

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

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

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

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

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

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

Несмотря на отсутствие механического контакта между статором и взвешиваемым телом, электромагнитный подвес может стать источником внутренних возмущающих моментов. Возникновение возмущающих моментов в магнитном подвесе связано с возможным отклонением электромагнитных сил от направления к центру ротора. Это отклонение может быть вызвано двумя причинами. Первая из них связана с погрешностью обработки поверхности ротора, приводящей к несовпадению нормалей к поверхности с направлениями к геометрическому центру ротора. Вторая причина обусловливается проявлением гистерезисных свойств ферромагнитного материала ротора. Гистерезисные явления приводят к тому, что элементарные магнитные силы могут отклоняться от нормали к поверхности ротора, из-за чего момент магнитного тяжения в подвесе может возникать даже при идеальной поверхности ротора. Вследствие этого положительный эффект, достигаемый применением магнитных подвесов в прецизионных датчиках первичной информации систем управления, во многом зависит от того, насколько успешно решена задача сведения к минимуму возмущающих моментов. Применение электромагнитных резонансных подвесов (ЭМРП) в качестве измерительного подвеса ограничено нелинейным характером зависимостей силы и индукции в рабочем зазоре подвеса. Кроме того, выходная информационная величина датчика на ЭМРП — это переменное напряжение высокой частоты, измерить которое с малой погрешностью затруднительно. Однако малый вес, небольшие габариты, малое потребление энергии ЭМРП были использованы при создании комбинированного электромагнитного подвеса, состоящего из ЭМРП и внешнего контура автоматического регулирования суммы индукции в рабочих зазорах подвеса. Такой подвес имеет линейную силовую характеристику, что позволяет успешно использовать его в качестве измерительного подвеса.

Получить динамическую устойчивость ЭМРП без использования внешнего механического демпфирования и улучшить его статические характеристики можно за счет:

— управления магнитным сопротивлением магнитопроводов подвеса;

— управления магнитным сопротивлением для потоков рассеяния;

— управления индуктивным сопротивлением резонансной схемы;

— управления емкостным сопротивлением резонансной схемы;

— управления активным сопротивлением резонансной схемы;

— управления сопротивлениями связей между резонансными схемами.

Наиболее перспективными методами коррекции характеристик ЭМРП являются:

— управление индуктивным сопротивлением резонансной схемы;

— управление емкостным сопротивлением резонансной схемы.

— комплексный подход к методам формирования и коррекции статических и динамических характеристик ЭМРП, к выбору рабочей точки на резонансной характеристике с учетом влияния корректирующих цепей.

Структурная схема инвариантного к внешней силе подвеса приведена на рис. 1. Электромагнитный подвес содержит: взвешиваемый якорь 1; диаметрально противоположные электромагниты 2 с рабочими 3 и дополнительными 4 обмотками; последовательно соединенные резонансные контуры, образованные рабочими обмотками диаметрально противоположных электромагнитов и конденсаторами 5; выпрямительные мосты 6; управляемые дроссели 7 с обмотками управления 8, подключенные к дополнительным обмоткам элек-

8-

//-

Рис. 1. Инвариантный подвес

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

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

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

— более высокая жесткость и быстродействие, так как явно выраженная обратная связь делает более легким регулирование крутизны изменения токов;

— возможность обеспечения меньшего уровня возмущающих моментов;

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

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

На рис. 2 показана функциональная схема одноосного электромагнитного подвеса с дифференциальным преобразователем силы. Подвес состоит: из взвешиваемого тела (ВТ); установленных с противоположных сторон от него электромагнитов ЭМ1 и ЭМ2; бесконтактного датчика положения (ДП); усилительно-преобразовательных блоков УПБ1 и УПБ2 и усилителей мощности УМ1 и УМ2.

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

При разработке измерительных электромагнит-

+Х

ЭМ1

V!

УМ1

УПБ1

ВТ

Р2

ЭМ2

Рис. 2. Активный электромагнитный подвес

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

0.2 %, даже в случае использования в качестве материала магнитопровода пермаллоя 79НМ.

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

На рис. 3 приведена функциональная схема активного дифференциального подвеса с управляемой жесткостью. Подвес содержит взвешенный элемент 1, расположенный в зазоре между двумя электромагнитами 2, обмотки которых соединены с выходами усилительно-преобразовательных блоков 3; датчик положения 4 взвешенного элемента, подключенный к входу блока динамической коррекции 5; блок суммирования 6 и блок вычитания 7; источник опорного сигнала 8; блок индикации 9; датчики индукции 10, чувствительные элементы 11 которых размещены в зазорах между взвешенным элементом и соответствующим электромагнитом, два интегратора 12; две схемы сравнения 13; управляемый делитель 14 и управляемый усилитель 15.

Задача регулирования величины рабочего зазо-

Рис. 3. Активный дифференциальный подвес с управляемой жесткостью

ра в зависимости от величины периода (частоты) напряжения питания решена в процессе синтеза квазиоптимального нелинейного детер-

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

На рис. 4 приведена функциональная схема такого подвеса. Подвес содержит взвешенный элемент 1, расположенный между двумя электромагнитами 2, обмотки которых соединены с выходами усилительно-преобразовательных блоков 3; датчик 4 положения взвешенного элемента, подключенный к входу блока 5 динамической коррекции; первый блок 6 суммирования и первый блок 7 вычитания, первые входы которых подключены к источнику 8 опорного сигнала; блок 9 индикации; датчики 10 индукции, чувствительные элементы 11 которых размещены в зазорах между взвешенным элементом и соответствующим электромагнитом; два интегратора 12; две схемы 13 сравне-

Рис. 4. Магнитный подвес, инвариантный по зазору

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

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

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

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

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

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

йУ / & + К У = и,

где У — потокосцепление, К = Я / Ь. В общем случае

и (і) = и0 (і) • вій

2 • п

Ті)

где и0{() и Тф являются функциями времени и выражаются через переменные состояния системы электромагнитного подвеса.

Уравнение движения взвешиваемого тела имеет вид

а2X ах „ + „

т •-— + Я~ = ^

йі2 йі

В 2 • 5 2 • ^0

т

масса взвешиваемого тела,

Я — коэффициент демпфирования, ^вн — величина внешней силы.

При построении математической модели электромагнитного преобразователя силы будем предполагать, что в выражении (2) амплитуду и период

і

возмущения можно представить суммой малых ступенчатых возмущений:

” m

u0(0 = ^u0l (t), T(t) = E Tj (t).

i=\ j=1

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

п гп

Y = E71i (t) +Е72j (t):

i=1

j=1

где У1г (¿) — реакция на г-е ступенчатое возмущение, обусловленное изменением амплитуды; У2  (^) — реакция на }-е ступенчатое возмущение,

обусловленное изменением величины периода.

Электромагнитная сила, действующая со стороны электромагнитов на взвешиваемое тело, определяется выражением [1]:

F3 (t) =-----------------------------

2 • Но •w • S

1

um

4 • Но

w

S ю"

• [l - cos 2(ю • t + a0 - ф)].

отношение:

Т

>> K.

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

U(t) = Um Sin[ t + 7o

и

Y(0) = 0, Т = То, Yo = 0, Um = Uo.

Тогда решение уравнения (1) можно записать в виде

Y (t) =

U о

(

Г sin

K 2 +

2 •п

то

-t-ф

U 0

• smф • в~

„2 4 п

K +-

T2

о

2-п

Или с учетом того, что-------->> K,

То

Y (t) =

= U о • Т о 2 •п

-sin

2 •п

То

-t-ф

U • Т

оо

2-п

smф • e

-Kt

= Yx(t) + Y2(t),

где ф = arctg

2 п

T~k

При наличии в системе демпфирования и вследствие того, что взвешиваемое тело инерционно, оно не реагирует на "быстроменяющиеся" периодические составляющие электромагнитной силы. Поэтому при анализе электромагнитного преобразователя силы необходимо выделить "медленноменяющиеся" составляющие уравнения (1). В выражении (1) параметр К в общем случае является переменной величиной и зависит от изменения величины индуктивности Ь. Однако изменение параметра К при соответствующем диапазоне изменения индуктивности, имеющем место в реальных системах, практически не влияет на динамику системы. Следовательно, параметр К в уравнении (1) можно считать постоянным. На всем диапазоне изменения величины Т периода напряжения питания справедливо следующее со-2  п

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

y(t)= ад + ад,

где Y (т) =

U^T

•sini

т + Y-ф I, Uи Т-

2 п у Т

текущие значения амплитуды и величины периода напряжения питания соответственно, у — начальная фаза, а У2 (т) определяется из решения урав-

нения:

Y2(t) + K • Y2 (т) = -а U -в •Т,

где

в =

Uо • sin

а = -

U = a AU 5(т), Т = а АТ 5(т), 5(т) — импульсная дельта-функция, AU = const, АТ = const,

/ 2п ,/ V 2п ,гг ,г ,гг

Y =^т tо, Y =—10, t0 и t0— моменты вре-

То

То

! 0 А 0 мени изменений амплитуды и периода напряжения питания соответственно, Т0 — начальное значение периода напряжения питания.

+

Kt

+

+

2

Мгновенная электромагнитная сила, действующая на взвешиваемое тело, определяется выражением

Рэ (і) =

2 • ^0

1

м

и

4^ ц0 •м2 ^ ю2

[1 --сов2

где ит — амплитуда переменного напряжения питания; ю — частота напряжения питания; ^ — число витков измерительной обмотки; S — площадь полюса; ц0 — магнитная проницаемость вакуума; ю— частота питания; «0 — начальная фаза.

Выражение для среднего значения электромагнитной силы с учетом выражений (3), (4), (6) будет иметь вид

РЭср = кр

и2 •Т2 8 -п2

+ Г22

= Кр

и12 Т12 и22 Т22

8 п2

8 п2

+ ^221

Разлагая первые два слагаемых в ряд Тейлора в точках и0, Т0, получим, выражение для электро-

магнитной силы в виде:

Рэср = °1 • АТ1 - Ь1 • АТ2 + а2 • Ди 1 - Ь2 • Ди2 +

+ а3 • ^21 Ь3 • ^22 :

где а1 = Ь1 = К

и2 Т

00

а 2 = Ь2 = КР

и 0 Т 4 п2

а3 = Ь3 = Кр .

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

Подставив уравнение (8) в уравнение движения взвешиваемого тела (2), получим следующее выражение:

т • X + Л • X = а1 • Т — Ь • Т2 + а 2 • и 1 —

- Ь •и + а •¥2 — Ь •¥2 + Р

2 и^и3 -'21 и3 22 ен

Вводя промежуточную переменную Х\ = Х и добавив к уравнению (9) уравнение (4) для двух плеч дифференциальной схемы, получим систему уравнений, характеризующую дифференциальный электромагнитный преобразователь силы [1, 2]:

Х = Хх;

Я а, Ъл а 2

X! = — X! + ^ Т --^ Т2 + ^ и! -т т т т

Ь2

-^2 и + У 2 _Ш У 2 | ~ вн •

^ 2 ~ 1 21 1 22 ~ ?

т т т т

У21 =_ КУ21 _« и/1 _в Т,;

У22 = _К  У22 _ «2  и2 _ в2  Т2 ,

где а и «2, в и в определяются выражениями (5). Причем отличие « от «2 и Д от в будет обусловлено величиной начальной фазы для каждого из плеч дифференциальной схемы в момент приложения соответствующего возмущения; и = сош!:; и2 = сош!:; ¥вн — величина внешней силы.

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

СИНТЕЗ ОПТИМАЛЬНОГО УПРАВЛЕНИЯ РАБОЧИМ ЗАЗОРОМ ПО ПЕРИОДУ ПИТАЮЩЕГО НАПРЯЖЕНИЯ

Задача синтеза оптимального управления на основе метода аналитического конструирования оптимальных регуляторов [3] предполагает, что объект описывается системой линейных дифференциальных уравнений вида:

Ь3

Рв

где Т = 2п/ю — величина периода напряжения питания.

Из выражения (7) следует, что первая составляющая среднего значения электромагнитной силы не зависит от величины начальной фазы у. Составляющая У2(т) зависит от начальной фазы у. Коэффициенты а и в определяются выражениями (5) и содержат составляющую 8т(у -ф), которая в зависимости от фазы у может принимать значения от -1 до 1. Для исследования вопросов устойчивости при дальнейших расчетах будем рассматривать худший случай, когда 8т(у_ф)=±1. Знак в этом выражении не имеет значения, поскольку величина ¥(т)в выражении (6) стоит в квадрате.

Для дифференциального электромагнитного преобразователя силы

¥ = ¥ _ ¥ =

эср эср1 эср2

X = А (і) • X + В(і) • и + Б(і) • М ,

где X — вектор фазовых координат, и — вектор управлений, М — вектор возмущений, А(і) — матрица динамических свойств системы, В(і) — матрица влияния управляющих воздействий, Б(і) — матрица влияния возмущающих воздействий.

Интегральный квадратичный критерий имеет следующий вид:

T

3 = J[x T(r). F(t). X(t) + •

0

+ U T(r) • Y(t) • U(t)]jt + + X t(T )• Fj(T )• X(T),

dK

dt

= K • A + A0 • K - KBY -1 BT • K + F

где К(0) и К(Т) — начальное ( = 0) и конечное (1 = Т) значения; К(Т) = 0.

1.1. Формируется расширенная матрица

A р =

A - BY-1B T

- F

- A0

C = exp(Ap H) =

E H . H2 . 2 V1

E—A р - ^2 A р

2

x

x

2

где Р(т), Р^Т) — неотрицательно определенные, а У(Т) — положительно определенная матрицы.

Оптимальное управление, минимизирующее критерий качества (11), можно найти по формуле:

и = -¥-1

где матрица К(1) есть решение дифференциального уравнения Риккати

_ К = К  А(1) + Ат (1)  К + Е(0 _

_ К  Б(1) У) Бт(Г) К,

К(Т) = Е(Т);

С — вектор решения линейного дифференциального уравнения

(5 = [К(1)  Б(0  У _ (1)  Бт(1) _ Ат(1)]  С(1) _

_ К (1 ) Б(1 ) М(1);

С(Т) = 0 .

Дифференциальное уравнение оптимальной системы имеет вид

X = [а(0 _ Б(0  У _' (0  Бт (1)  К(1)]  X +

+ Б(0  М(0 _ Б(0  У _ (1)  Бт (1)  С(1) =

= А0(1)  X + Б(1)  М(1) _ Б(1)  С(1).

Итеративный метод решения задачи состоит из ряда этапов.

1. Вычисляется матрица начального приближения Унач= К(0) путем решения дифференциального уравнения Риккати

т. Н А Н . 2

Е + — А р +------А р

2 р 12 р

V

Т

где шаг решения Н = -—, к — целое, выбирается в

2к

пределах к=5-20.

1.3. Последовательно вычисляются:

ехр(Ар 2Н) = С С = С2; ехр(Ар 22Н) = С2 С2 = С4;

exp(Ap 2 H) = exp(Ap T) = C2 .

1.4. Матрица exp(Ap T) разбивается на четыре блока

exp( A pT) = C2 =

Pll

P21

1.5. Вычисляется матрица

1.2. Вычисляется матрица

К (0) = ^^2_1 (Т ) P2l(Т ).

2. Вычисляются

Lo = У-1 Бт Унaч, А0 =А - Б Lo.

При этом L0 выбирается так, чтобы матрица

А0 = А — Б L0

имела собственные значения с отрицательной действительной частью.

3. Решается линейное алгебраическое уравнение относительно У0:

A О • V0 + У0 •A 0 = C1

гдеC1 = L'T YL0 + F — положительно полуоп-ределенная матрица.

Так как Унач = K(0), это уравнение является уравнением Ляпунова. Для его решения используем итеративный алгоритм [4], имеющий линейную сходимость. Можно показать, что У0 = lim У0к,

где

V)( к+1) =(dT ) V0k D2к + У0к, (к = 0,1,2,...);

V00 = H ^1;

Б

Е — Н •А 0 —

2

Н2 12

А02

х

ГЕ Н . Н2 . 2л

Е +----А 0 +-----А 0

2 0 12 0

V у

Определение величины У0 производится аналогично первому этапу.

4. Зная величину У0, по формулам

Ь к = У —1 В т А к = А — ВЬ к

к—1

определяем Ll, Аь Далее, используя формулы 3-го этапа, вычисляем У1. Аналогично определяем У2, У3,... , Ук . Вычисления прекращаем по достижении требуемой точности.

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

В функционале качества (11) выберем весовые матрицы Е и У стационарными и диагональными с коэффициентами

Р. =

1

ірег • X

2

і тах

1

і и2

р ег у тах

где Хітах— максимальное допустимое отклонение і-й фазовой координаты от заданного положения, и  тах— максимальное допустимое отклонение для у-го управления, ірег — среднее время регулирования.

Подставляя значения параметров исследуемого активного электромагнитного подвеса в выражения для элементов матриц А, В, Е, У и производя необходимое масштабирование переменных, получим следующие выражения для этих матриц.

Согласно (12), оптимальный регулятор должен работать по закону и0= -У-1 Вт К X, где матрица К — решение матричного алгебраического уравнения Риккати

КА + АТ К+ Ф - В У1 ВТ К= 0.

Решение этого уравнения может быть осуществлено с помощью итерационного алгоритма.

Расширенная матрица Ар, необходимая для вычисления начального приближения К(0), имеет вид

А р =

А — ВУ —В т

— Е

—Ат

Обозначив L = -У 1 Бт К, где L= (I}) — матрица

коэффициентов усиления регулятора по цепям обратных связей, получим уравнение регулятора

и = L X.

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

Расчеты, выполненные для конкретной системы управления, показали, что линейное управление для нелинейного объекта дает значительно худший результат, чем оптимальное линейное управление без учета нелинейности объекта. Из полученных результатов можно сделать вывод о том, что для достижения заданного качества регулирования и приближения качественных показателей нелинейной системы к качественным показателям линейной системы необходимо введение нелинейной коррекции закона (13). Задача синтеза в этом случае сводится к решению дифференциального уравнения в частных производных, для которого не существует общего решения. Однако можно искать пути приближенного решения этого уравнения.

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

X = А X + Би,

т.е. найти первое (линейное) приближение квази-оптимального управления

ил(1) = У_1 Бт К(0^X .

— Выполнить синтез оптимального линейного регулятора, вычисляя матрицу линейных обратных связей размером

L = -У 1 Бт К.

— Затем уточнить оптимальное управление

и(0 = що + що.

— После этого решается задача нахождения нелинейной составляющей

ин(1) = _У_1 Бт Ях

х|тоО] _^М

и(^0) = _У _Ч^ Б т(^) х

Х-Иад, 10]^ X(t 0) + 8^(10), ^0 ]}

8^)] = Я ^[X(to)] _д[XXt0)]  X(10)| ,

где 10 — текущий момент времени; Х(10) — измеренные значения вектора состояния системы; Я — матрица размером т х п.

Для получения решения ин(0 необходимо вычислить матрицу

Я = _(ат _ К  Б  У _1 Бт )_1 К

и сформировать матрицу нелинейных обратных связей размером

N = _У_1 Бт Я.

Для упрощения реализации получаемого регу-

д Г

лятора в выражении (14) положим = 0.

Нелинейная вектор-функция f(X) для исследуемого активного электромагнитного подвеса с регулированием рабочего зазора по величине периода напряжения питания имеет вид

Г T(X) = |0 284.3  X 32 0.

Матрица N для определения нелинейных обратных связей (7) имеет вид:

Nт = |2.95  105 1.84  102 0.122|.

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

ин = _ 5.23  104 Х32.

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

ЗАКЛЮЧЕНИЕ

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

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

1. Сапожников Г.А. Повышение эффективности гашения механических колебаний в электромагнитных подвесах // Датчики систем контроля и управления технологическими процессами. Ярославль: Ярославский политехнический институт,

1984. С. 64-68.

2. Богословский С.В. Теория и практика аэромагнитного моделирования (монография). СПб.: ГУАП, 1998. 140 с.

3. Летов А.М. Аналитическое конструирование регуляторов // Автоматика и телемеханика. 1969. т. 21, № 4. С. 436-441; № 5. С. 561-568.

4. Колесников А.А. Принципы построения высококачественных систем квазиоптимального управления // В сб. "Методы построения алгоритмических моделей сложных систем": Вып. 1. Таганрог: ТРТИ, 1976. С. 166-167.

Санкт-Петербургский государственный университет аэрокосмического приборостроения

Материал поступил в редакцию 24.04.2001.

OPTIMIZATION OF CONTROL OF MAGNETIC LEVITATION IN SENSORS

G. A. Sapozhnikov, S.V. Bogoslovsky

Saint-Petersburg State University of Aerospace Instrumentation

The results of development of sensors based on magnetic levitation are briefly described. The mathematical model of control of the power supply voltage period to change the operating clearance in a magnetic suspension is considered. To optimize control of a magnetically-suspended element the method of analytical design of regulators is utilized.