Научная статья на тему 'Обеспечение статической и динамической точности микроэлектромеханического датчика угловой скорости'

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

CC BY
421
106
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МЭМС / ДАТЧИК УГЛОВОЙ СКОРОСТИ / СТАТИЧЕСКАЯ ХАРАКТЕРИСТИКА / ТОЧНОСТЬ / ЧУВСТВИТЕЛЬНОСТЬ / СРЕДНИЙ КВАДРАТ ОШИБКИ / МОДЕЛИРОВАНИЕ / MEMS / SENSOR OF ANGULAR VELOCITY / STEADY-STATE FEATURE / ACCURACY / SENSITIVITY / AVERAGE SQUARE OF MISTAKE / MODELING

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Волков Валерий Леонидович, Хрисанова Марина Васильевна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Волков Валерий Леонидович, Хрисанова Марина Васильевна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

PROVISION STEADY-STATE AND DYNAMIC ACCURACY OF THE MICRO ELECTROMECHANIC SENSOR OF ANGULAR VELOCITY

Purpose: Making methods of calculation to accuracy micro electro mechanic sensor on base of analysis of influence constructive parameter on steepness of steady-state feature. Design/methodology/approach: is Formed mathematical model of micro sensor of angular velocity (DUS) on base principle functioning the micro gyroscope LL-type. The Considered methods of calculation of main features DUS on base of programming in ambience Matlab. The Considered method of modeling DUS for reason statistical discovery factor of influence constructive parameter on steepness of steady-state feature. Findings: is Presented results of model test DUS in the manner of tables factor to sensitivity parameter and graphics to connecting function. On graph is shown achievement it is enough small time of connecting process of sensor. Will Revealed dominant parameters, defining instability of steepness of steady-state feature of sensor. They Are Determined optimum parameters, corresponding to minimum dynamic mistake of sensor. Research limitations/implications: Using the offered methods of calculation of features DUS gives the minimum fair square-law of mistake of sensor for given a range of measurement. Originaly/value: On base of considered methods of designing and studies DUS is worded requirements to electro mechanic parameter MEMS sensor on criterion steady-state and dynamic accuracy.

Текст научной работы на тему «Обеспечение статической и динамической точности микроэлектромеханического датчика угловой скорости»

УДК 629.7.05/06: 531.781.2 (075.8)

В.Л. Волков, М.В. Хрисанова

ОБЕСПЕЧЕНИЕ СТАТИЧЕСКОЙ И ДИНАМИЧЕСКОЙ ТОЧНОСТИ МИКРОЭЛЕКТРОМЕХАНИЧЕСКОГО ДАТЧИКА УГЛОВОЙ СКОРОСТИ

Арзамасский политехнический институт (филиал НГТУ)

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

Ключевые слова: МЭМС, датчик угловой скорости, статическая характеристика, точность, чувствительность, средний квадрат ошибки, моделирование.

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

В ряде работ по микросистемной технике затрагивались вопросы обеспечения статической и динамической точности датчиков [1, 2]. Так для микромеханического акселерометра, модель которого представлена в работах [1, 3], была показана возможность построения корректирующего устройства (КУ) с параметрами, обеспечивающими минимальные значения динамических ошибок [1]. Там же были выявлены предельно допустимые отклонения параметров КУ относительно расчетных.

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

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

Задача определения коэффициентов чувствительности параметров для одного из основных параметров ДУС - крутизны статической характеристики (К), как и для акселерометра, может быть решена методом статистических испытаний [2]. При статистических испытаниях производится моделирование работы датчика с номинальными значениями параметров при случайных вариациях этих параметров. В каждом испытании вычисляются отклонения варьируемых параметров от их номинальных значений - Щ(к), а также показателя качества К относительно своего номинального значения - 5К(к). Коэффициент чувствительности ^(к) у-го параметра в к-м испытании при этом определяется в виде

% (к) = ^. (1) А ) С (к) ( )

Среднее значение коэффициента влияния для каждого параметра вычисляется в виде:

% =± ^ • (2)

к=1 ^

где N - число статистических испытаний.Коэффициенты чувствительности показателя каче-

© Волков В.Л., Хрисанова М.В., 2013.

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

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

итах / „/шах-К (3)

где /тах - верхнее предельное значение входного сигнала на заданном диапазоне измерения; итах - максимальное значение выходного сигнала в вольтах.

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

5К= А К, (4)

где А=Ъит /итах - предельно допустимая относительная ошибка датчика, 5ит - предельное отклонение выходного напряжения датчика на пределе измерения.

В предположении того, что ошибка датчика возникает только из-за у-го параметра, определим его предельно допустимое отклонение б/ от номинального значения по формуле

5/ = А-итах/Чу. (5)

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

п

зип = ). (6)

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

Математическая модель датчика угловой скорости

Рассмотрим ДУС на основе вибрационного микрогироскопа ££-типа [4]. Принципиальная схема такого ДУС представлена на рис. 1. Датчик включает в себя чувствительный элемент (ЧЭ) из двух интегральных масс, выполненных в виде внутренней 1 и наружной 2 рамок. Внутренняя рамка 1 относительно внешней рамки 2 смонтирована на упругих подвесах 3, в свою очередь рамка 2 посредством упругих подвесов 4 связана с анкерами 5.

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

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

Рис. 1. Принципиальная схема чувствительного элемента ДУС

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

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

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

Рис. 2. Структурная схема ДУС:

КЧЭ - коэффициент передачи чувствительного элемента; ЖПУ - передаточная функция подвижного узла; КПП - коэффициент передачи емкостного преобразователя перемещений; ЖКУ - передаточная функция корректирующего устройства; Кос - коэффициент передачи звена обратной связи

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

WKW

W (5) = Кта-"ПУ K ПП"КУ-. (7)

1 + W К W К

1 ^ " ПУКпп" кукос

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

Параметры электростатического возбудителя вибраций, передаточные функции чувствительного элемента и подвижного узла рассчитываются на основе конструкции аналогичной микрогироскопу ADXRS фирмы Analog Device [6].

Коэффициент передачи ЧЭ имеет вид:

Кчэ (s) = 2mvx = 2mX0 со (8)

где Xo - амплитуда колебаний в режиме движения подвижного узла; m - масса подвижного узла; ю - резонансная частота.

Максимальное перемещение подвижного узла X0 при заданных параметрах упругого подвеса и электростатического преобразователя на резонансной частоте рассчитывается по формуле:

ЕЕа Ь и2 а2 О

у _ Д — Е0ЕаМПУЬМПУи ОПаХО /0\

х 0 Атах Ш Ь С { )

где 8о, в - соответственно диэлектрическая постоянная и диэлектрическая проницаемость среды между электродами; El00 - модуль упругости кремния в плоскости 100; иОП - опорное напряжение; h - величина зазора между пластинами конденсаторов; ампу, Ъмпу - длина и ширина подвижного узла; ах, Ъх, ctx - длина, ширина и толщина торсиона внутренней рамки; Q - добротность колебательной системы.

Передаточная функция подвижного узла, как колебательного звена имеет вид:

КУ (з) = 1

ms2 + Kds + Gy (10)

где Оу = Е[100]Ь1ус^ /(16а1у) - жесткость торсионов, работающих на изгиб в направлении оси

Y; Кд - коэффициент демпфирования.

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

Кд= (11)

Коэффициент передачи емкостного преобразователя перемещений имеет вид:

Кпп =ио (12)

п

Коэффициент передачи звена электростатической обратной связи имеет вид:

Кос =2в^в1кЬкиоп / п2, (13)

где 4, bk - длина и ширина электрода обратной связи.

Методика расчета основных характеристик ДУС

С целью расчета конструктивных параметров ДУС и передаточных функций в программной среде Matlab разработана программа-функция [Z,Wpu,Wche,WЪe,W]=dus(P), входными данными которой является вектор параметров Р (подвижного узла, блока электроники, конструкции):

P=[Ar1,Adop,Bdop,h,Aco,Bo,AЪo,Ayo,C,Dr2,Lr2,lk,Ъk,Atx,Btx,Aty,Bty,Ctx,Cty,Ky,Uop,T\,

где в соответствии с конструкцией ЧЭ микрогироскопа [6] обозначено: Ar1 - длина стороны внутренней рамки; Adop - длина дополнительной части внутренней рамки; Bdop - ширина дополнительной части внутренней рамки; h - величина зазора между пластинами конденсаторов (в программе ф0); Aco - длина центрального окна; Bo - ширина окон; Abo - длина бокового окна; Ayo - длина углового окна; C - толщина подвижного элемента; Dr2 - ширина внешней рамки; Lr2 - длина стороны внешней рамки; lk - длина пальца конденсатора чувствительности; bk - ширина пальца конденсатора чувствительности; Atx - длина торсиона внутренней рамки; Btx - ширина торсиона внутренней рамки; Ctx - толщина торсиона внутренней рамки; Aty - длина торсиона внешней рамки; Bty - ширина торсиона внешней рамки; Cty - толщина торсиона внешней рамки.

Дополнительно введены обозначения: Ky - коэффициент усилителя мощности, Uop -опорное напряжение, T - постоянная времени блока электроники.

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

Wpu - подвижного узла; Wche - чувствительного элемента; Wbe - блока электроники; W - полная передаточная функция ДУС и вектор расчетных параметров ДУС:

Z=[m, w, Gy,Koc,K],

где m - масса чувствительного элемента; w - круговая частота электростатического привода; Gy - жесткость одного торсиона внешней рамки; Koc - коэффициент обратной связи; K - коэффициент передачи датчика.

Текст программы-функции dus(P) на языке Matlab составлен в соответствии с теоретическими формулами механики, формулами для расчета коэффициента чувствительного элемента, первичного преобразователя, датчика обратной связи [3, 4], формулами для вычисления передаточных функций. Примеры такого рода программ приводились в работах [2, 5].

Использование программы-функции dus(P) осуществляется путем предварительного задания входных данных - вектора параметров P. Все исходные параметры ЧЭ, необходимые для расчета сведены в файл input.dat текстового формата ascii. Для расчета выходных параметров выполняется вызов dus(P) из управляющей программы control dus:

Po=load( 'input.dat ');

P= Po;

[Z,Wpu,Kche,Wbe,W]=dus(P) % вызов программы-функции.

При расчете параметров Ctx, Cty, определяющих диапазон измерения датчика, сначала рассчитывается коэффициент передачи датчика K, соответствующий диапазону измерения (3), и подбираются численные значения Ctx, Cty, удовлетворяющие этому диапазону.

Пример. При исходных данных ДУС для заданного диапазона измерения 150 град/с определены результирующие параметры. При использовании программы-функции dus рассчитаны все параметры ДУС (основные параметры приведены в табл. 1).

Таблица 1

Основные параметры ДУС для диапазона ±150 град/с

Параметр Kche m, кг w, рад/с Gy Koc K Kpp Uop, В Umax, В T, с

Значение 2,03e-9 4,03e-7 1,68e+3 2,03 8,85e-10 2.292 4e+10 20 6 1.6e-3

Параметр Wmax Atx Btx Ctx Aty Bty Cty lk bk

Значение 1,9 e-10 3,24e-4 1,0e-5 2e-6 1,0e-4 1,0e-5 2e-6 1e-4 1e-5

Соответствующие полученным параметрам передаточные функции представлены в табл. 2.

Таблица 2

Передаточные функции ДУС, соответствующие структурной схеме рис. 2

Передаточная функция подвижного узла Передаточная функция КУ Полная передаточная функция датчика

Т1Г, ч 4,03s2 +1,28-104 s + 2,03-107 W (s) =- s , ч 4,54 -10~7 s2 +1,44 -10~3 s + 2,292

4-107s2 +1,28-Ш"3 s + 2,03 "0(s) 2,8-1045 s3 + 2,0-10~7 s2 + 6,3-10 s +1

Далее решая задачу обоснования требований к параметрам ДУС необходимо определить коэффициенты чувствительности параметров (вектора параметров P), на основной показатель качества - крутизну статической характеристики K. Используем стандартную формулу чувствительности (2) с учетом (1).

Вариация 5f j-го параметра задается относительно номинального значения параметра fj, соответствующего заданному диапазону измерения датчика. Вариация 5f производится в пределах стандартного 5% допуска на исходные параметры по равномерному закону распределения. Проводятся многочисленные испытания, в которых фиксируются как вариации 5f, так и соответствующие изменения крутизны. Оценка коэффициента чувствительности j-го параметра производится на основе среднестатистического значения (3).

На основе метода статистического моделирования составлена управляющая программа dus_kj.m для расчета коэффициентов чувствительности параметров, текст которой представлен далее:

Po=load( input.daf У, % ввод исходных данных

[Z,Wpu,Wche,Wbe,W]= dus(Po); % расчет параметров ДУС для номинальных значений

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Po

K=dcgain(W); % крутизна статической характеристики для заданного диапазона

% Вычисление Ko(k) - крутизны при вариации s(k) j-го параметра

j=i;

while j<23

P= Po; an(j)=P(j); k=1;

while k<1000 % задание цикла в 1000 статистических испытаний P(j)=0.95*an(j)+0.1*an(j)*rand(1,1); % вариация j-го параметра [Z,Wpu,Wche,Wbe,W]= dus(P); % вычисление параметров ДУС Ko(k)=dcgain(W); s(k)=P(j); % запоминание K и j-го параметра в массив e(k)=Ko(k)-K; e1(k)=s(k)-Po(j); % расчет отклонений

kjk(k)= e(k)/e1(k); % коэффициент влияния j-го параметра в k-м испытании k=k+1; end

kj(j)=mean(kjky; % среднее значение коэффициента влияния j-го параметра

j=j+1; end

save(^result2.txi,P\ykf ,}-ascii}) % запись результатов в текстовый файл.

Эта программа вызывает подпрограмму-функцию dus(P), предназначенную для расчета параметров ДУС. При использовании dus_kj.m по формуле найдены коэффициенты чувствительности (коэффициенты влияния) параметров fi ДПИ на выходную характеристику K (3). Были получены результаты по коэффициентам влияния и допускам 5j на параметры для различных предельных погрешностей датчика. Результаты для различных предельных погрешностей датчика (0,1 % погрешность и 1 % погрешность) сведены в табл. 3.

Анализ результатов табл. 3 показал, что при 1% предельной погрешности датчика предъявляются вполне реальные требования по стабильности параметров. Особые требова-

7 8

ния на уровне 10" - 10" м предъявляются только к параметрам торсионов внутренней рамки по ширине и толщине.

Таблица 3

Коэффициенты чувствительности ^ и допуски 8; параметров

Параметр f Статистический номинал f ^rSK/bf 5i (при 0,1% точности) 5i (при 1% точности)

P(1)=Ar1 1.5e-3 5.5504687e+3 1.0809898e-6 1.0780415e-5

P(2)=Adop 1.215e-4 2.9833036e+2 2.0111933e-5 2.0111949e-4

P(3)=Bdop 8.1e-5 4.4750308e+2 1.3407729e-5 1.3407685e-4

P(4)=h 2.0e-5 1.1458764e+5 5.2361667e-8 5.2361667e-7

P(5)=Aco 1.242e-3 -7.9524027e+2 -7.5448896e-6 -7.5400269e-5

P(6)=Bo 2.16e-4 -8.5489977e+3 -7.0183666e-7 -7.0180283e-6

P(7)=Abo 8.1e-4 -7.9576011e+2 -7.5399607e-6 -7.5427061e-5

P(8)=Ayo 1.35e-4 -1.5911840e+3 -3.7707770e-6 -3.7704028e-5

P(9)=C 7,238e-5 1.5812170e+4 3.7945457e-7 3.7773373e-6

P(10)=Dr2 1.9e-4 6.1885089e+3 9.6953889e-7 9.6953889e-6

P(11)=Zr2 1.62e-3 7.2581277e+2 8.2665947e-6 8.2665947e-5

P(12)=lk 1e-4 -2.2845917e+4 -2.6262899e-7 -2.6217972e-6

P(13)=bk 1e-5 -2.3344354e+5 -2.5702146e-8 -2.6181010e-7

P(14)=Atx 3,24e-4 1.0750782e+4 5.5809893e-7 5.5643837e-6

P(15)=Btx 1e-5 -1.1013787e+5 -5.4477174e-8 -5.4691715e-7

P (16)= Aty 1.0e-4 3.8200672e+1 1.5706530e-4 1.5706530e-3

P(17)=By 1e-5 3.8200672e+2 1.5706530e-5 1.5706530e-4

P(18)=Ctx 2e-6 -1.7396985e+6 -3.4488735e-9 -3.5071179e-8

P(19)=Cty 2.0e-6 3.6232268e-9 1.6559824e+6 -2.1959583e+7

P(20)=Ku 9.0e+4 2.3614328e-20 2.5408303e+17 -6.9695992e+017

P(21)=Uop 2.0e+1 1.1458764e-1 5.2361667e-2 5.2361667e-1

Определение параметров корректирующего устройства

В структурной схеме датчика (см. рис. 2) для улучшения динамических свойств датчика предусмотрено корректирующее устройство с передаточной функцией Wky, соответствующей ПИД регулятору.

Wy (s) = Ts + T2 + T3/ s, (15)

где T1, T2, 73, - постоянные времени корректирующего устройства.

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

п

ошибки СКО при работе микроакселерометра с применением ПИД КУ (на уровне 10" , что для подобных приборов достаточно неплохой показатель).

При использовании ПИД корректирующего устройства в ДУС необходимо продолжить исследование его положительных свойств и определить достижимость наибольшей динамической точности датчика в зависимости от параметров КУ. В работе [7] утверждается, что с целью синхронизации работы узлов микросистемного датчика коэффициенты ПИД КУ должны быть одинаковыми, но доказательство этого утверждения не приводится. Очевидно, что при равенстве коэффициентов 71, 72, T3 важные качества переходного процесса датчика получатся далеко не оптимальными.

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

Запишем передаточную функцию прямой цепи в соответствии со структурной схемой рис. 2 на основании (7) с учетом передаточной функции КУ (15).

WnP (s) _ WnyKПП^КУ _

1

Js2 + K,s + G,

K ПП \TiS + T2 + T s

(16)

V + Vу "Л" ' 5 ) Используя эту передаточную функцию можно получить необходимые и достаточные условия для получения графика переходного процесса с малой колебательностью (или вообще с отсутствием таковой).

Наиболее очевидным и простым решением является подбор параметров КУ таким образом, чтобы нули КУ совпали с полюсами подвижного узла. При этом необходимо выполнить условия: Т2/Т1=Кду//, Г3/Г1= Gy// и тогда скобки (э2 + Т2/Т1 ^ + Г3/Г1) и (э2 + Кду// s + ОуЦ) сокращаются. В результате из (17) получим более простое выражение для передаточной функции прямой цепи.

(*) = Щ/О/• . (17)

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

W (s) = K4

кШУ( J • s)

к

1

(18)

1 + КосКппТ /(/ • 5) Кос /(КосЩ )1 5 +1

Обозначая ККУ=Т1//, как коэффициент передачи КУ, а То=(КппКосКку) — как постоянную времени датчика, получим передаточную функцию замкнутой системы в виде:

Г (5) =Кчэ 1 . (19)

Кос то5 +1

Из теории известно, что полученной передаточной функции датчика соответствует полоса пропускания равная юо=1 /ТО, и время переходного процесса равное ¿ПП=3ТО. Реализация ТО при ранее вычисленных параметрах КПП, КОС возможна только за счет Т1. Параметр Т1 определяет эффективность действия дифференцирующей составляющей КУ. Другие параметры Т2, Т3 оказываются тесно связанными с Т1 соотношениями: Т2 = Т1Кду/1/, Т3 = ТуОу// (в отличие от работы [7], где Т1з Т2, Т3 предполагались равными).

Моделирование ДУС в МаАаЪ 81шиИпк

С целью реального исследования работоспособности ДУС и проверки его характеристик создана адекватная модель в визуальной среде Ыа^аЪ БШиЫпк [8] при реализации рассчитанных ранее параметров.

Схема моделирования ДУС в Ыа^аЪ 8тиМпк представлена на рис. 3. Схема предусматривает вычисление ошибки отработки входного скачкообразного сигнала структурой ДУС. Причем при дальнейшей, уже программной обработке массива ошибок, полученных в блоке е, рассчитывается экспериментальное значение среднего квадрата ошибки.

T2

Scope

©HD

e

Рис. 3. Схема моделирования ДУС в Matlab Simulink

При моделировании получен график переходного процесса при максимальном входном сигнале 150 град/с (см. рис. 4), показывающий апериодичность достижения максимального выходного сигнала 6 В за время примерно 710" с.

Рис. 4. Переходный процесс датчика

На основе моделирования по схеме 3 произведена оценка СКО, результаты которой для данного примера соответствуют минимальному значению: СКО=3,58-10"6.

Выводы

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

Одним из теоретических результатов в данной работе является доказательство возможности получения динамической характеристики микродатчика с апериодическим переходным процессом с достаточно малым временем переходного процесса (¿ПП=3ТО, где постоянная времени ТО теоретически может быть выбрана сколь угодно малой).

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

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

1. Вавилов, В.Д. Оптимизация параметров микромеханического акселерометра / В.Д. Вавилов,

B.Л. Волков, А.В. Улюшкин // Труды НГТУ им. Р.Е. Алексеева. - Н. Новгород, 2010. № 3 (82).

C. 308-314.

2. Волков, В.Л. Обоснование требований к параметрам микромеханического акселерометра // Труды НГТУ им. Р.Е. Алексеева. - Н. Новгород, 2011. № 2 (87). С. 288-295.

3. Вавилов, В.Д. Теоретические основы микросистемных акселерометров и гироскопов / В.Д. Вавилов; НГТУ. - Н. Новгород, 2011. - 210 с.

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

5. Вавилов, В.Д. Компьютерное моделирование характеристик микросистемных датчиков /

B.Д. Вавилов; НГТУ. - Н. Новгород, 2007. - 80 с.

6. Шахнович, И. МЭМС гироскопы - единство выбора. Электроника НТБ. 2007. № 1.

C. 76-85.

7. Вавилов, В.Д. Оптимизация характеристик микроакселерометра со смещенной осью качания маятника // Труды НГТУ им. Р.Е. Алексеева. - Н. Новгород, 2012. № 2 (95). С. 309-318.

8. Дьяконов, В.П. Simulink 5/6/7: самоучитель / В.П. Дьяконов. - М.: ДМК Пресс, 2008. - 784 с.

Дата поступления в редакцию 12.04.2013

V.L. Volkov, M.V. Chrisanova

VI. PROVISION STEADY-STATE AND DYNAMIC ACCURACY OF THE MICRO ELECTROMECHANIC SENSOR OF ANGULAR VELOCITY

Arzamas polytechnical institute (branch) Nizhny Novgorod state technical universiy

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

n.a. R.Y. Alexeev

Purpose: Making methods of calculation to accuracy micro electro mechanic sensor on base of analysis of influence constructive parameter on steepness of steady-state feature.

Design/methodology/approach: is Formed mathematical model of micro sensor of angular velocity (DUS) on base principle functioning the micro gyroscope LL-type. The Considered methods of calculation of main features DUS on base of programming in ambience Matlab. The Considered method of modeling DUS for reason statistical discovery factor of influence constructive parameter on steepness of steady-state feature.

Findings: is Presented results of model test DUS in the manner of tables factor to sensitivity parameter and graphics to connecting function. On graph is shown achievement it is enough small time of connecting process of sensor. Will Revealed dominant parameters, defining instability of steepness of steady-state feature of sensor. They Are Determined optimum parameters, corresponding to minimum dynamic mistake of sensor.

Research limitations/implications: Using the offered methods of calculation of features DUS gives the minimum fair square-law of mistake of sensor for given a range of measurement.

Originaly/value: On base of considered methods of designing and studies DUS is worded requirements to electro mechanic parameter MEMS sensor on criterion steady-state and dynamic accuracy.

Key words: MEMS, sensor of angular velocity, steady-state feature, accuracy, sensitivity, average square of mistake, modeling.

i Надоели баннеры? Вы всегда можете отключить рекламу.