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

УДК 629.7.054

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

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

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

© 2015 И. В. Фоминов

Военно-космическая академия имени А.Ф. Можайского, г. Санкт-Петербург

doi: 10.18287/1998-6629-2015-14-1-92-100

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

Введение

таты которых приведены в работах [1, 2], показали, что существенного повышения чувствительности и расширения диапазона измеряемых ускорений маятникового компенсационного акселерометра (МА) можно добиться путём искусственного введения его чувствительного элемента (ЧЭ) в режим автоколебаний при помощи импульсной обратной связи. При этом ЧЭ совершает колебательные движения относительно центра динамического равновесия. Центр динамического равновесия смещается при появлении ускорения или под действием возмущающих сил и моментов. Знакопеременный сигнал на выходе нелинейного звена (НЗ) в таком режиме представляет собой меандр, обработка которого аналогична режиму с ши-ротно-импульсной модуляцией. Величина,

пропорциональная измеряемому кажущемуся ускорению, определяется следующим выражением:

n =

nii + n2i

где пь - число импульсов за интервал времени при положительном значении меандра; п21 - число импульсов за интервал времени при отрицательном значении меандра.

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

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

Постановка задачи

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

= ( /д (х, ФI), /пи ^з X /ар (Х, 7 X /из (Руаг X Руаг, Рйх) ;

7 : В х С х /а; х = х + 5х(Ф) ; 5х(Ф) = Ах0 +Ах0 (Ф, Раг) + Ух (Ртаг) ,

где / (х, Фх) - функция инерциального датчика; /пи (инз) - функция преобразова-

теля импульсов; / (X, Z) - функция адаптивного регулятора; /нз (Pvar) - функция нелинейного звена типа «петля гистерезиса»; P/lx - множество фиксируемых

параметров инерциального датчика; Pvar -

множество варьируемых параметров инерциального датчика - параметров нелинейного звена; x - измеряемая физическая величина на входе инерциального датчика; X - результат измерений инерциального датчика; 5х(Ф) - погрешность измерений инерциального датчика; Ф -высокочастотное возмущающее воздействие; Ах0 - систематическая погрешность инерциального датчика, определённая в нормальных условиях; Ах0 (Ф) -квазисистематическая погрешность инер-циального датчика, обусловленная влиянием возмущающего воздействия Ф; инз - выходной параметр нелинейного звена; Z - множество априорных данных, хранящихся в базе знаний; {B, C} - множество параметров нелинейного звена, {B, C} е Pvar; /а - частота автоколебаний

чувствительного элемента.

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

Pvar (Ч) = arg min 5х(Pvar, Ф) ,

P eQ

где А7а - интервал времени адаптации, необходимый для поиска оптимальных значений параметров Pvar и формирования

их в нелинейном звене; Qa - множество

допустимых значений параметров нелинейного звена.

Синтез адаптивного инерциального датчика

Из теории адаптивных систем известно [4], что синтез адаптивной системы осуществляется декомпозицией системы на основной контур (неадаптивный) и адаптивный контур (контур самонастрой-

ки). В этой связи целесообразно структуру адаптивного инерциального датчика, приведённую на рис. 1, условно разделить на основной контур и адаптивный контур, образованный преобразователем импульсов (ПИ), адаптивным регулятором (АР) и базой знаний (БЗ). В нижней части рис. 1 изображена схема преобразования сигнала на входе НЗ.

На рис. 1 приняты следующие обозначения: - сила инерции; - сила обратной связи; в - угол отклонения ЧЭ; Ир - сигнал с выхода датчика перемещения (ДП), пропорциональный углу перемещения; В - ширина «петли гистерезиса» НЗ; С - размах «петли гистерезиса» НЗ; А - актюатор (датчик силы).

Ф.

Р.

АР

х ЧЭ ПЗ в ДП НЗ и нз

——*

А

х

ПИ

и нз —

t

БЗ

х

Т

Рис. 1. Структура адаптивного инерциального измерителя

Параметрический синтез основного передаточную функцию замкнутой систе-

контура адаптивного инерциального из- мы инерциального датчика, функциониру-

мерителя как нелинейной системы прово- ющего в режиме автоколебаний, можно

дится путём применения метода гармони- представить в следующем виде: ческой линеаризации [5]. В этом случае

^ма ( Р ) = ■

кчэ 1кдп /нз (В, С)

Т2р2 + 2£Тр +1 + kчэ I ^п /нз (В, с) / т

(1)

где Т - постоянная времени ЧЭ; - относительный коэффициент демпфирования ЧЭ; кчэ - коэффициент передачи ЧЭ; I - длина ЧЭ; АН - величина перемещения ЧЭ; кдп - коэффициент передачи ДП;

кдс - коэффициент передачи датчика силы; Ф(В, С) - линеаризованная функция, описывающая НЗ.

Приведём выражение (1) к передаточной функции второго порядка посред-

ством гармонической линеаризации НЗ

[5]:

к

Щ р) =

где к1 =

Т2 р2 + 2£Т1р +1

кчэ 1кДпФ(В, С) 1 + кЛчФ(В, С)

(2)

Т1 =

Т 2

1 + клчФ(В, С)

^ V1 + *лчФ(В, С)

¿лч = ¿чэ 1Кп ¿дс / .

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

Дю) =| Щ (>)| =

¿1

^(1 - г2ю2 )2 + 4^12Г12Ю2

(3)

где ш - собственная частота замкнутой системы.

Из выражения (3) видно, что свойства фильтрации внешних воздействий проявляются на частотах ш > 1/ Т1. Постоянная времени Т1 зависит от параметров НЗ: В, С. Фильтрующие свойства инерциального датчика иллюстрирует рис. 2, на котором представлены логарифмические амплитудно-частотные характеристики для различ-

ных частот автоколебаний /а, полученные экспериментальным путём по исходным данным, принятым в работе [6].

Снижение погрешностей инерци-ального датчика, обусловленных влиянием нестационарных внешних воздействий Ф, возможно посредством формирования определённой частоты автоколебаний /а . Результаты исследования влияния частоты автоколебаний /а на среднеквадратиче-скую погрешность авых, обусловленную частотой вибрации /вибр основания маятникового акселерометра, приведены на рис. 3.

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

I, дБ

-20

-40

-60

-80

-100

-120

-140

—/ак=100О(5 Гц /ак =3550, 3 Гц

/ак=102,9 Гц

10

10

10

10

10

Ю, Гц

Рис.2. Логарифмические амплитудно-частотные характеристики инерциального датчика

о м/с2 x 1ü"3 .■••»""

вых ' л |и .....

8.V-6 4 2

0

500

/ибр , Гц

1500 3000

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

Условием проведения адаптации инерциального датчика является необходимость малости и неизменности измеряемой величины: |x(t)| < хдоп, |x(t)| < Хдоп.

Для акселерометров такие условия могут быть достигнуты в режиме пассивного орбитального полёта КА, в котором кажущиеся ускорения, как правило, не

превышают значений 10-5 g .

Алгоритм реализации метода имеет следующий вид:

1. Формируется массив выходных данных за время At :

={x(tk),...,*(tk + At)}, к = \K .

2. Производится стохастический

анализ массива x

Atk

и вычисление мате-

матического ожидания М[], средне-квадратического отклонения а[ХхАйк ] и значений частот в выходных данных /вых[х&к ] =1(ХАк ), где Цх&к) - оператор быстрого преобразования Фурье.

3. Определяется требуемая частота

автоколебаний /' , при которой будет

соблюдаться условие инвариантности:

Лых (/вибр, /а') = В результате анализа экспериментальных исследований эмпирическим путём получено выражение

./вых = | /вибр - Х/а | , , позволяющее

утверждать, что для инвариантности погрешности, обусловленной влиянием вибраций, необходимо, чтобы частота автоколебаний /а была кратна частоте внешних вибраций /вибр. При этом значение показателя кратности х задаётся исходя из близости значения начальной частоты автоколебаний /а0 (В0, С0), рассчитанной до проведения адаптации, к требуемой частоте автоколебаний /' .

4. Решается задача фильтрации внешних возмущающих факторов путём определения нового значения параметра Вк+1 из множества 2 : В х С х /а, при кото-

ром достигается минимум частоты /вых:

Лых ^

5. Формируется новое значение параметра В в соответствии с рекуррентным соотношением Вк+1 = Вк + Хк АВ0 до достижения условия /вых «0, где

Q

AB •

Xk f (Ч ), AS0 k Afa(Atk У 0

6. Проводится стохастический анализ последней реализации xA i и определяется её квазисистематическая составляющая погрешности измерения, равная в данных условиях математическому ожиданию Aie = M [ xA ].

7. Производится компенсация квазисистематической составляющей погрешности путём формирования значения параметра Ck+1 в соответствии с формулой Ck+i =Ck (1 + Ae).

Апробация метода проведена путём численного моделирования процесса функционирования маятникового акселерометра в режиме автоколебаний с учётом

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

т = 0,29-10-3 кг, кчэ = 4,117-10-3 с2 /м, Т = 4,636 -10-3 с, Ь = 16,323, I = 5,086 -10-3м, кдп = 2,5 -105 В/м,

кдс = 2,155 -10-3 Н/В .

Результаты исследования влияния значений параметров НЗ на частоту и амплитуду автоколебаний показаны на рис. 4.

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

зависимости Ак от параметров НЗ представлен на рис. 4, б.

fа , ГЦ

1500010000 5000ч

B, В

Ah, 10-6 м

6 0

C, В

б

Рис. 4. Графики зависимостей частоты / и амплитуды Ак автоколебаний ЧЭ от параметров нелинейного звена маятникового акселерометра

6

5

а

Анализ представленных на рис. 4 графиков показал, что для обеспечения требуемой частоты автоколебаний f при

ограничениях на их амплитуду Ah и заданный диапазон измерения существует такое оптимальное значение параметра B, при котором достигается минимум средне-квадратической ошибки (СКО) измерения ускорения min аа.

Проведённые исследования показывают, что уменьшение частоты автоколебаний ниже определённого критического значения приводит к снижению информативности выходного сигнала за временной интервал, равный периоду автоколебаний. Длительность адаптации составляет 1-2 периода автоколебаний, что при частоте в 1200 Гц соответствует тысячным долям секунды (рис. 5).

м/с

1.4

1.3

1.2

1.1

0.9

0.8

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0 с

Рис. 5. График реализации выходной характеристики в процессе адаптации

1

Результаты математического моделирования процесса функционирования адаптивного маятникового акселерометра подтвердили снижение СКО погрешности с оа =0.0239 м/с2 до оа =0.0022 м/с2. При этом в качестве входных воздействий задавались постоянное ускорение а=1 м/с2, вибрация с амплитудой авибр =0,05 м/с и

частотой /вибр =6° Гц. Время на проведение самонастройки составило 0,08 с.

Заключение

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

Библиографический список

1. Рыбаков В.И., Фоминов И.В. Способ измерения навигационных параметров подвижных объектов автоколебательными датчиками первичной информации. СПБ.: Военно-космическая академия им. А.Ф. Можайского, 2005. С. 240-244.

2. Рыбаков В.И., Фоминов И.В. Инерциальный блок измерения абсолютной угловой скорости КА // Известия высших учебных заведений. Приборостроение. 2005. Т. 49, № 7. С. 37-43.

3. Фоминов И.В., Малетин А.Н. Алгоритм самонастройки маятникового автоколебательного акселерометра при воз-

действии периодических возмущений // Известия высших учебных заведений. Приборостроение. 2011. Т. 54, № 9. С. 2833.

4. Фрадков А.Л. Адаптивное управление в сложных системах. М.: Наука, 1990. 296 с.

5. Бесекерский В.А., Попов Е.П. Теория систем автоматического управления. СПб: Профессия, 2004. 752 с.

6. Распопов В.Я. Микромеханические приборы: учебное пособие. М.: Машиностроение, 2007. 400 с.

Информация об авторе

Фоминов Иван Вячеславович, кан- i.v.fominov@gmail.com. Область научных

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

кафедры автономных систем управления, движением космических аппаратов. ВКА имени А.Ф. Можайского. E-mail:

METHOD OF PARAMETRIC OPTIMIZATION OF EACH COMPONENT OF ADAPTIVE COMPENSATNG INERTIAL SENSORS OPERATING IN A SELF-OSCILATION MODE

© 2015 I. V. Fominov

Military Space Academy named after A.F. Mozhaiskiy, St. Petersburg, Russian Federation

The paper deals with a crucial task of parametric optimization of an adaptive inertial sensor capable of varying the parameters in order to reduce errors which are caused by the influence of high-frequency transient disturbances. A method of parametric optimization of each component that makes it possible to vary the parameters of the nonlinear element of a compensating inertial sensor operating in the self-oscillation mode is developed. The method was tested with a pendulous accelerometer taken as an example. Numerical investigations confirmed the possibility of filtering external vibrations and internal noise of the inertial sensor while maintaining its sensitivity in the required measuring range due to the synthesis of the adaptive loop based on a non-linear element of the "hysteresis loop" type. The essence of the method is based on parametric synthesis of the adaptive loop with a nonlinear element. Numerical investigations show that the error of the pendulous accelerometer can be reduced, on the average, by an order of magnitude thanks to the adaptive parameter setting of the non-linear element, depending on the amplitude of the external vibrations. The above method of parametric optimization of each component can be applied to most compensating inertial sensors measuring motion parameters of spacecraft operating in conditions of non-stationary external and internal disturbances. The adaptation circuit in adaptive inertial sensors can be implemented through the use of modern microcontrollers.

Adaptive inertial sensor, pendulous accelerometer, parametric optimization, self-oscillation mode, error, external vibration.

References

1. Rybakov V.I., Fominov I.V. Sposob izmereniya navigatsionnykh parametrov podvizhnykh ob"ektov avtokolebatel'nymi datchikami pervichnoy informatsii [Method of measuring navigation parameters of moving objects by self-oscillating sensors]. Spb.: Military Space Academy named after A.F. Mozhaiskiy Publ., 2005. P. 240-244.

2. Rybakov V.I., Fominov I.V. Inertial Device for Measurement of Absolute Angular Rate of a Space Vehicle. Izvestia vuzov. Priborostroenie. 2006. V. 49, no.7. P. 37-43. (In Russ.)

3. Fominov I.V., Maletin A.N. Algorithm of self-adjustment of pendulous auto-oscillating accelerometer under the influence of high-

frequency periodic disturbances. Izvestia vuzov. Priborostroenie. 2011. V. 54, no. 9. P. 28-33. (In Russ.)

4. Fradkov A.L. Adaptivnoe upravlenie v slozhnykh sistemakh [Adaptive control in complex systems]. Moscow: Nauka Publ., 1990. 296 p.

5. Besekerskiy V.A., Popov E.P. Teoriya system avtomaticheskogo upravleniya [Theory of automatic control systems]. SPb.: Professiya Publ., 2004. 752 p.

6. Raspopov V.Y. Mikromekhanicheskie pribory: uchebnoe posobie [Micromechanical devices: study guide]. Moscow: Mashi-nostroenie Publ., 2007. 400 p.

About t

Fominov Ivan Vyacheslavovich, Candidate of Science (Engineering), doctoral student, Department of Autonomous Control Systems, Military Space Academy named after A.F.

author

Mozhaisky, St. Petersburg, Russian Federation. E-mail: i.v.fominov@gmail.com. Area of Research: navigation and control of spacecraft motion.