RU2362173C1 - Способ измерения линейного ускорения в системах управления ракет и устройство для его осуществления - Google Patents

Способ измерения линейного ускорения в системах управления ракет и устройство для его осуществления Download PDF

Info

Publication number
RU2362173C1
RU2362173C1 RU2008110986/28A RU2008110986A RU2362173C1 RU 2362173 C1 RU2362173 C1 RU 2362173C1 RU 2008110986/28 A RU2008110986/28 A RU 2008110986/28A RU 2008110986 A RU2008110986 A RU 2008110986A RU 2362173 C1 RU2362173 C1 RU 2362173C1
Authority
RU
Russia
Prior art keywords
range
acceleration
δδa
error
max
Prior art date
Application number
RU2008110986/28A
Other languages
English (en)
Inventor
Вадим Иванович Рыбаков (RU)
Вадим Иванович Рыбаков
Иван Николаевич Штанов (RU)
Иван Николаевич Штанов
Original Assignee
Государственное образовательное учреждение высшего профессионального образования "Военно-космическая академия имени А.Ф. Можайского"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Государственное образовательное учреждение высшего профессионального образования "Военно-космическая академия имени А.Ф. Можайского" filed Critical Государственное образовательное учреждение высшего профессионального образования "Военно-космическая академия имени А.Ф. Можайского"
Priority to RU2008110986/28A priority Critical patent/RU2362173C1/ru
Application granted granted Critical
Publication of RU2362173C1 publication Critical patent/RU2362173C1/ru

Links

Images

Abstract

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

Description

Изобретение относится к способам и устройствам в области измерительной техники, которая занимается вопросами измерения линейных ускорений в системах инерциальной навигации самолетов, ракет, кораблей, космических аппаратов и других подвижных объектов (ПО). Изобретение может применяться в указанных отраслях техники, но в большей степени соответствует системам управления ракет, которые требуют обеспечения высокой точности измерения ускорения в конце активного участка траектории движения, в момент выключения маршевого двигателя при заданной величине кажущегося ускорения и скорости.
В литературе известны способы измерения линейных ускорений, реализованные в устройствах на основе компенсационных акселерометров, заключающиеся в том, что ускорение преобразуют в отклонение подвижной системы акселерометра с последующим преобразованием отклонения в электрический сигнал, его усиливают и преобразуют в момент отрицательной обратной связи, а параметры сигнала отрицательной обратной связи (значения: тока, числа импульсов, кодового числа) используют в качестве выходного сигнала для определения ускорения. При этом необходимым условием применения способа является проведение калибровки акселерометров, в процессе которой по результатам измерений определяют значения технических характеристик (ТХ), и используют их для контроля соответствия заданным значениям тактико-технических требований (ТТТ) системы управления ПО, в том числе заданной точности измерения ускорения [2, 3], а также для дальнейшего их учета в моделях погрешностей для систем управления движением и навигации ПО.
В последнее время уделяется внимание способам и устройствам, которые обладают свойствами диагностики своих ТХ без дополнительных измерительных средств. Такими свойствами обладают автоколебательные акселерометры, которые рассматриваются в качестве прототипа предлагаемому изобретению.
Среди ряда известных [4-7] в наибольшей степени предлагаемому соответствует способ измерения линейного ускорения с использованием автоколебательного акселерометра, заключающийся в том, что ускорение преобразуется в отклонение подвижной системы акселерометра с последующим преобразованием отклонения в электрический сигнал, который усиливается и преобразуется в последовательность временных интервалов, сформированных с помощью источника стабилизированного тока и элементов нелинейного звена в виде широтно-модулированного сигнала, который преобразуют в момент импульсной отрицательной обратной связи при помощи датчика момента, а сигналы импульсной отрицательной обратной связи преобразуют в сигнал, модулированный счетными импульсами, поступающими от генератора импульсов стабилизированной частоты fc, с последующим измерением количества счетных импульсов n1i и n2i, которые преобразуют в выходные сигналы акселерометра в виде разности Δni=n1i-n2i, суммы ni=n1i+n2i=fcTi и отношения количества импульсов Ni=Δni/ni, причем величину и знак входного воздействия aвхi+Δai, соответствующего измеряемым ускорению ai, и погрешностям Δai, определяют по окончании каждого i-того периода Тi автоколебаний с помощью выходного сигнала и по формуле Ni=Kii+Δаi), (i=1…∞), где Кi=ml/КдмI0 - коэффициент преобразования входного воздействия, ml - маятниковый момент подвижной системы, Кдм и I0 - коэффициент передачи датчика момента и поступающий от стабилизированного источника ток, величина которого ограничивает диапазон измеряемого ускорения, при этом коэффициент преобразования Кi, погрешность, например смещение нуля Δai=Δa0, их стабильность определяют при помощи формулы Nikaiki+Δа0), градуировочную характеристику (ГХ) и ее стабильность (линейность) при помощи другой формулы Nijaij(aij+Δa0j), определяют заранее и каждую из указанных ТХ раздельно по результатам калибровки в условиях гравитационного поля Земли при вертикальном и горизонтальном положениях измерительной оси акселерометра, причем измерения каждой ТХ выполняют на интервале времени сотен периодов автоколебаний, чтобы применить статистическую обработку и фильтрацию результатов измерений, которые используют в дальнейшем при измерении ускорения.
В результате анализа известных устройств измерения ускорения на основе компенсационных автоколебательных акселерометров [1, 4-7] для реализации известного способа измерения ускорения был выбран состав устройства, наиболее близкий предлагаемому. Структурная схема такого устройства показана на фиг.1, где изображен акселерометр, содержащий чувствительный элемент (ЧЭ) и нелинейное звено (НЗ), при помощи которого осуществляется режим автоколебаний. В состав ЧЭ входит: инерционный элемент (ИЭ) 1, датчик угла (ДУ) 2 и 3, датчик момента (ДМ) 7, которые механически и электрически связаны между собой и образуют замкнутый контур управления вместе с элементами, реализующими НЗ. Подвижные части ДУ, ДМ и ИЭ называют подвижной системой (ПС) ЧЭ. В состав НЗ входит: усилитель-преобразователь (УП) 4, компаратор 5 и электронный ключ (ЭК) 6. В состав устройства входят также: источник стабилизированного тока (ИСТ0) 8, генератор счетных импульсов стабилизированной частоты (ГСИ) 9, реверсивный (РСИ) 10 и суммирующий (ССИ) 11 счетчики импульсов и спецвычислитель (СВ) 12.
При наличии входного воздействия авхi+Δa0 ИЭ отклоняется на угол β и на обмотке ДУ 3 возникает сигнал в виде напряжения Uдудуβ, которое после усиления в УП в виде сигнала Uуп=KупUдудуКупβ поступает на вход компаратора, выход которого связан с управляющими входами ЭК и счетчиков РСИ и ССИ. При помощи ЭК источник ИСТ0 подключается к ДМ. Образуется электрическая цепь импульсов тока Iдм, поступающих через контакты ЭК на обмотку ДМ, в котором возникает момент импульсной отрицательной обратной связи. Выбор величины тока I0 источника ИСТ0 позволяет обеспечить процесс функционирования автоколебаний ЧЭ и измерения ускорений в заданном диапазоне. Одновременно функционирует информационный канал и ГСИ подключается к счетным входам РСИ и ССИ. На выходе счетчиков выделяется разность Δni=n1i-n2i и сумма ni=n1i+n2i=fcTi количества импульсов, пропорциональных входному воздействию и периоду автоколебаний Ti, а на выходе СВ получают информацию в виде Ni=Δni/ni. Значения сигналов Δni=n1i-n2i и
ni=n1i+n2i поступают от счетчиков на информационный вход СВ, а на его выходе получают информацию в виде числа Nii(ai+Δаi), пропорционального измеряемому ускорению и погрешностям.
При проведении калибровки ГОСТ [3] рекомендует использовать ряд технических характеристик, приведенных в таблице 1, которые применительно к автоколебательным акселерометрам представим в следующем виде:
- коэффициент преобразования акселерометра определяют путем измерения ускорения известной величины при наличии погрешности смещения нуля Δai=Δa0. При заданном диапазоне ±aвх=±amax+Δa0 удобнее фиксировать значения границ диапазона, что и является предметом измерения. На выходе устройства получают: Nmax1=Ka1
(+amax+Δa0), -Nmax2=Ka2(-amax+Δa0), а вычисление значений выполняют по формулам:
Figure 00000001
- смещение нуля измеряют при ускорении аi=0, тогда на выходе устройства, а при известном коэффициенте преобразования и на входе соответственно получают:
Figure 00000002
где Кa=1/2(Ka1a2) - среднее значение коэффициента преобразования;
- градуировочная характеристика измеряется в заданном диапазоне ускорений ±aвх=±amax+Δa0, и на выходе устройства получают:
Figure 00000003
где aij+Δа0jвхj - величина ускорения и погрешность смещения нуля, которые задают в виде дискретных значений, как правило, в нескольких точках j=(1, 2…10) заданного диапазона, Кaij - коэффициент преобразования ГХ. Линейность ГХ в большей степени зависит от качества оборудования, ошибок оператора и поэтому определение ее целесообразно выполнять и паспортизировать на заводе-изготовителе в условиях прецизионного оборудования, а на этапах предстартовой подготовки и применения в системах инерциальной навигации ПО использовать паспортные значения.
Полученные результаты измерений при калибровке подвергают математической обработке с использованием аппарата математической статистики [8] для определения стабильности и линейности измеряемых технических характеристик (ТХ), для чего требуется определенный объем полезной информации, поэтому в случае с автоколебательными акселерометрами время измерений ограничивают несколькими сотнями периодов автоколебаний Тiq≈qmaxTi≈(500-1000)Ti, где qmах≈(500-1000), - количество периодов автоколебаний, которое выбирают в интересах статистики при измерении каждой ТХ в процессе калибровки.
Основным недостатком известных способов и устройств измерения линейных ускорений является несоответствие полученных в процессе калибровки ТХ реальным условиям измерения ускорения и друг другу, так как каждая из характеристик определяется на различном оборудовании и при различном влиянии условий окружающей среды за достаточно длительный промежуток времени, в течение которого изменение внешних факторов существенно влияет на величины ТХ и их нестабильности. Как следствие - низкая точность измерения ускорения по причине использования ТХ, рассчитанных заранее и с методическими ошибками. Несмотря на то, что проводят калибровку с целью уменьшения погрешностей и определения ТХ для дальнейшего их использования в процессе измерения, указанная цель не может быть достигнута в полной мере при существующих способах калибровки, которые в свою очередь не позволяют обеспечить стабильность измеряемых ТХ.
Приведенные выражения (1)-(3) показывают, что ни одна из получаемых ТХ в современных системах не может выполнить свое назначение в соответствии с определением по ГОСТу [3], так как измерения их проводят раздельно, последовательно, на различном оборудовании, на различных временных интервалах, и, следовательно, при различных внешних возмущениях, что является причиной нестабильности определяемых ТХ, их взаимосвязи и взаимообусловленности в процессе измерения ускорения и калибровки. Использование аппарата математической статистики, якобы для повышения стабильности и линейности измеряемых ТХ, не решает указанной задачи, а вычислением количественных показателей в виде математического ожидания, дисперсии или СКО каждой ТХ, подчеркивает влияние внешних возмущений на точность измерения искомых величин.
Например, величина выходного сигнала при измерении ускорения в соответствии с известным способом определяется равенством Ni=Ki(ai+Δai) и пропорциональна входному сигналу в виде измеряемого ускорения и погрешностей, возникающих вследствие влияния возмущений. Однако коэффициент преобразования Ki измеряется при калибровке в других условиях и определяется с помощью выражений (1) при наличии погрешности смещения нуля, которая измеряется при третьих условиях, а вычисляется при помощи выражений (2) и наличии того же коэффициента преобразования.
Отсюда следует, что точность измерения ускорения зависит не только от возмущающих воздействий, влияющих на величину измеряемого ускорения непосредственно, но и от стабильности значений ТХ, которые зависят также от возмущающих воздействий, и при определении каждой из ТХ раздельно будут отличаться, что вызывает взаимосвязи и взаимообусловленность измеряемых ТХ, влияние их нестабильности на точность измерения ускорения, что видно при анализе приведенных формул.
Таким образом, выходом из данного противоречия является сближение по времени условий измерения и испытаний, когда внешние возмущения проявят себя несущественно. Первым шагом в направлении повышения стабильности измеряемых ТХ при калибровке может быть обеспечение единства их определения в условиях неизменного влияния окружающей среды за счет короткого по времени однократного измерения двух воздействий известной величины ускорения, одинаковых по модулю, но разных по знаку, при помощи которых могут быть вычислены все искомые характеристики одновременно.
Следовательно, применение известных способов измерения ускорений требует проведения калибровки акселерометров, которая привносит дополнительные погрешности, материальные и временные затраты, необходимые для выполнения измерений и вычислений. Проведение измерений каждой из искомых ТХ отличаются по длительности и сложности, для чего применяют сравнительно дорогостоящее лабораторное оборудование в виде виброустойчивых наклонно-поворотных оснований, оптических делительных головок, уровней, что усложняет процесс подготовки и является существенным недостатком известных способов измерения линейных ускорений [2, 7].
Особенности проведения калибровки акселерометров указывают на возможности уменьшения инструментальных погрешностей, которые ограничены технологией производства, а при калибровке могут быть измерены и скомпенсированы до уровня стабильности источника компенсации. Особая роль принадлежит погрешности смещения нуля, которая является составной частью измерений ТХ акселерометра, но, судя по выражению ΔN0=KaΔa0, может быть измерена при известном коэффициенте преобразования до определенного уровня, значение которого ограничено измеренным значением коэффициента преобразования, величина которого для заданного диапазона измеряемого ускорения неизменна. Следовательно, известные способы не позволяют повысить величину и стабильность коэффициента преобразования на заданном диапазоне путем увеличения чувствительности применяемых средств.
Проведенные предварительные исследования автоколебательных акселерометров открыли значительный резерв в направлении повышения стабильности основных ТХ при проведении калибровки, в том числе уменьшения погрешности смещения нуля в сторону приближения к погрешности порога чувствительности.
Задачей изобретения в способе и устройстве измерения линейного ускорения на основе автоколебательного акселерометра является повышение точности измерения ускорения в заданном диапазоне за счет уменьшения погрешности смещения нуля путем ее измерения на узком диапазоне, кратном заданному с последующей компенсацией непосредственным образом до уровня остаточной, а также за счет повышения стабильностей ТХ: коэффициента преобразования, смещения нуля и ГХ, определяемых при калибровке путем одновременного определения коэффициента преобразования, смещения нуля и ГХ за счет короткого по времени и однократного измерения известных величин ускорения, что приведет к упрощению способа измерения, к сокращению материальных и временных затрат при незначительном увеличении массогабаритных показателей.
Поставленная задача изобретения в способе и устройстве по пунктам формулы 1 и 2 достигается тем, что на основе заданного значения линейности ГХ, полученного и паспортизованного в результате калибровки акселерометра на заводе-изготовителе в условиях его прецизионного оборудования и заданного диапазона измеряемого ускорения, на начальном этапе проведения калибровки выбирают границы диапазона смещения нуля ±ai0=±K0amax, которые больше величины погрешности смещения нуля Δa0, но меньше величин границ заданного диапазона ±ai0<<±amах, и кратны им, определяют коэффициент кратности этого диапазона K0=ai0/amax<<1, запоминают его и определяют величину тока I00=K0I0 диапазона смещения нуля, повышают величину коэффициента преобразования акселерометра на этом диапазоне, для чего вместо тока I0 в дополнительной обмотке датчика момента формируют ток диапазона смещения нуля I00 и регулируют его величину с помощью регулятора так, чтобы граничным значениям ускорения ±ai0=K0amax соответствовали максимальные значения выходного сигнала ±Ni0=±Nmax0 акселерометра, выполняют калибровку акселерометра на диапазоне смещения нуля в условиях гравитационного поля Земли при горизонтальном положении измерительной оси акселерометра на интервале времени в несколько сотен периодов автоколебаний, проводят измерения ускорения известной величины, например, на границах диапазона смещения нуля с учетом погрешности смещения нуля, например, путем поворота измерительной оси прибора на угол, соответствующий воздействию долей единицы ускорения
Figure 00000004
: +aвх01=+g×sinα+Δa0=+0,0001g+Δa0 и -aвх01=-g×sinα+Δa0=-0,001g+Δa0, где α - угол между измерительной осью прибора и
Figure 00000005
, полученные значения выходного сигнала Nmах01=Ka01(+ai0+Δa0) и -Nmax02=Ka02(-ai0+Δa0) суммируют алгебраически и делят на удвоенное значение ускорения границы диапазона смещения нуля (Nmax01-(-Nmах02))/2ai0=Kа0, а затем суммируют арифметически и делят пополам
(Nmax01+(-Nmax02))/2=Kа0Δa0, определяют стабильное значение коэффициента преобразования Ka0 и значение погрешности смещения нуля на выходе с учетом знака ΔN00=Ka0Δa00, погрешность смещения нуля на входе компенсируют непосредственно до уровня остаточной δΔa00 путем формирования в дополнительной обмотке датчика момента тока компенсации IK0, и регулируют его величину при условии ai0=0, пока на выходе акселерометра в процессе измерений не появится минимальное значение δΔN00=Ka0δΔa00, где δΔa00 - значение скомпенсированной погрешности смещения нуля на входе на диапазоне смещения нуля, и оставляют включенной цепь компенсации, после чего в основной обмотке датчика момента восстанавливают ток I0 заданного диапазона, повторяют калибровку акселерометра на заданном диапазоне в условиях гравитационного поля Земли при вертикальном положении измерительной оси акселерометра на интервале времени в несколько сотен периодов автоколебаний путем короткого по времени однократного измерения двух известных величин входного ускорения, например, на границах заданного диапазона с учетом скомпенсированной погрешности смещения нуля, авх1=+amax+δΔa0, aвх2=-amax+δΔa0, полученные значения выходных сигналов
Nmах1=Ka1(+amax+δΔa0) и -Nmax2a2(-amax+δΔa0) суммируют алгебраически и делят на удвоенное значение ускорения границы заданного диапазона
(Nmax1-(-Nmax2))/2amaxa, а затем суммируют арифметически и делят пополам
(Nmax1+(-Nmах2))/2=KaδΔa0, определяют стабильное значение коэффициента преобразования для заданного диапазона Кa, остаточную погрешность смещения нуля на выходе δΔN0aδΔa0 и на входе δΔa0=δΔN0/Ka, стабильность полученных значений вычисляют при помощи отношения измеренных и вычисленных алгебраической и арифметической сумм (Nmах1-Nmах2)/(Nmах1+Nmах2)=KaδΔa0/Kаamax=δΔa0/amах, где δΔa0/amах - относительная величина скомпенсированной погрешности на заданном диапазоне, а представленная в процентах δΔa0100%/amaxδΔa0 характеризует стабильность измеренной погрешности и стабильность ГХ одновременно, причем при стабильном коэффициенте преобразования Ka=const полученные значения запоминают и используют при формировании полетного задания для измерения ускорения в заданном диапазоне с помощью выходного сигнала устройства и по формуле Nai=Ka(ai+δΔa0) при меньшем значении остаточной погрешности смещения нуля δΔa0<<Δa0.
Предлагаемый способ измерения линейного ускорения при помощи устройства на основе акселерометра, работающего в режиме автоколебаний, позволяет измерить и определить значения коэффициента преобразования, остаточной погрешности смещения нуля одновременно при коротком по времени однократном измерении двух величин входного воздействия, равных значениям границ заданного диапазона, что соответствует физической природе определения ГХ. Одновременное определение искомых ТХ по результатам однократного измерения двух величин входного воздействия, равных значениям границ заданного диапазона, причем время измерений искомых параметров ограничивают временем, соответствующим нескольким периодам автоколебаний, что уменьшает нежелательное воздействие внешних возмущений посредством взаимного влияния ТХ друг на друга и имеет место у известных способов при раздельном и последовательном измерении по времени и определении каждой из указанных ТХ при помощи различных технических средств и приспособлений. Кроме того, повышение точности измерения ускорения произойдет вследствие непосредственной компенсации измеренной погрешности смещения нуля на диапазоне смещения нуля до уровня остаточной, при этом величина измеренной остаточной погрешности смещения нуля на заданном диапазоне будет меньше исходной δΔa0<<Δa0, соответствующей известным способам и устройствам.
Заявленный способ измерения линейного ускорения, предложенный по пункту 1 формулы изобретения, может быть осуществлен устройством, предложенным по пункту 2 формулы, структурная схема представлена на фиг.2, за счет того, что введены: блок управления (БУ) 13, дополнительный регулятор PTI00 с измерителем тока диапазона смещения нуля mA00, дополнительный регулятор РТI0K с измерителем тока компенсации, в датчик момента введена дополнительная обмотка, входы регуляторов электрически связаны с источником стабилизированного тока (8), выход регулятора тока компенсации РТI0K через измеритель тока и нормально разомкнутые контакты блока управления связан с дополнительной обмоткой датчика момента, а выход регулятора с измерителем тока PTI0 связан через нормально замкнутые контакты блока управления с обмоткой датчика момента, при этом на начальном этапе калибровки по команде блока управления его контакты отключают обмотку датчика момента от регулятора с измерителем тока PTI0, а подключают к этой обмотке дополнительный регулятор с измерителем тока PTI00 диапазона смещения нуля, в режиме компенсации дополнительную обмотку датчика момента при помощи контактов БУ электрически связывают с регулятором и измерителем тока компенсации РТI0K. Кроме выходного параметра Ni, устройство и способ позволяют одновременно определять и получать на выходе основные ТХ: стабильный коэффициент преобразования Кa, остаточную погрешность смещения нуля δΔa0 и относительную скомпенсированную погрешность смещения нуля γδΔa0 - для дальнейшего их учета в моделях погрешностей для систем управления движением и навигации подвижных объектов.
Принципы повышения коэффициента преобразования в узком диапазоне измеряемого ускорения в предлагаемом способе можно пояснить при помощи схемы и графиков, изображенных на фиг.2 и 3. Как было сказано выше, на время настройки по команде блока управления нормально замкнутые контакты размыкаются, а нормально разомкнутые контакты замыкаются и обмотка датчика момента переключается на регулятор PTI00 с измерителем тока диапазона смещения нуля mA00, где формируют значение тока I00=K0I0. При настройке этого режима величину тока регулируют, пока выходной сигнал не примет максимальное значение Ni0=Nmах0, пропорциональное величине входного воздействия границы диапазона смещения нуля. Выбирают значение тока I00=K0I0 и создают возможность для повышения величины коэффициента преобразования на меньшем диапазоне.
Для сравнительного анализа на фиг.3 представлены графики градуировочных характеристик: характеристика X1 соответствует заданному диапазону ±аmах=±1g, характеристика Х2 - диапазону смещения нуля с граничными значениями ±авх0=±аmaх×10-3=±1×10-3 g. Характеристики построены по результатам эксперимента, проведенного на акселерометре АЛЕ-048, переведенного в режим автоколебаний. Для построения характеристики Х2 заданный диапазон по оси абсцисс был разделен на восемь равных интервалов с ценой деления aвx0=amax×10-3 g, равной границам диапазона смещения нуля. Сравнение полученных ГХ показывает существенное отличие значений коэффициента преобразования, измеренной погрешности на выходе ΔN0=KaΔa0 (точка 4) для заданного диапазона (градуировочная характеристика XI) и на диапазоне смещения нуля ΔN00=Kа0Δa0 (точка 2, ГХ Х2) для одного и того же значения погрешности смещения нуля на входе (точка 1), а следовательно, возникает возможность компенсации измеренной погрешности до уровня остаточной δΔN00a0δΔа00 (точка 4 для ГХ Х2') на диапазоне смещения нуля. При переходе на заданный диапазон измерения, когда в основной и дополнительной обмотках датчика момента поданы токи I0 и IK0 одновременно, т.е. реализован режим компенсации погрешности смещения нуля до уровня остаточной, ГХ Х2' трансформируется в ГХ X1', a δΔN0 будет соответствовать точка 5 на фиг.3 при погрешности смещения нуля на входе δΔa0 (точка 3). Если применить одну и ту же систему компенсации в указанных случаях, то значение остаточной погрешности в узком диапазоне будет меньше на два-три порядка, чем у известных способов.
Компенсацию измеренной погрешности ΔN00=Ka0Δa0 на диапазоне смещения нуля выполняют при условии ai0=0 и при помощи регулятора с измерителем тока, которые подключают источник к дополнительной обмотке датчика момента через контакты КЗ БУ. При этом ток компенсации регулируют, пока на выходе не появится минимальное значение ΔN00=Nmin0=Kа0δΔa00, которое измеряют и запоминают, где δΔa00 - значение скомпенсированной погрешности смещения нуля по входу; когда это значение будет скомпенсировано на диапазоне смещения нуля, контакты КЗ БУ оставляют включенными.
После этого по команде БУ вместо тока I00 в дополнительной обмотке датчика момента в основной обмотке восстанавливают ток I0 заданного диапазона, повторяют калибровку акселерометра на заданном диапазоне, для чего повторяют измерения на границах заданного диапазона Nmax1=Ka1(+amax+δΔa0), -Nmax2=Ka2(-amax+δΔa0), вычисляют арифметическую (Nmax1+(-Nmax2))/2=KaδΔa0 и алгебраическую (Nmax1-(-Nmax2))/2=Ka суммы полученных сигналов, уточняют стабильное значение коэффициента преобразования Кa, погрешности смещения нуля на выходе δΔN0aδΔа0 и на входе δΔa0=δΔN0a, стабильность полученных погрешностей вычисляют при помощи отношения измеренных и вычисленных алгебраической и арифметической сумм (Nmax1-Nmax2)/(Nmax1+Nmax2)=KaδΔa0/Kaamax=δΔa0/amax, где δΔa0/amах - относительная величина скомпенсированной погрешности на заданном диапазоне, а представленная в процентах δΔa0100%/amaxδΔa0 характеризует стабильность измеренной погрешности и стабильность градуировочной характеристики одновременно, причем при стабильном коэффициенте преобразования Кa=const, полученные значения запоминают и используют при формировании полетного задания для измерения ускорения в заданном диапазоне с помощью выходного сигнала устройства и по формуле Nai=Ka(ai+δΔa0) при меньшем значении остаточной погрешности смещения нуля δΔa0<<Δa0.
Возможность достижения положительного эффекта предлагаемого способа по сравнению с известными была проверена экспериментальными исследованиями акселерометра АЛЕ-048, переключенного в режим автоколебаний, при физическом и математическом моделировании. На макете указанного устройства (структурная схема представлена на фигуре 1), при настройке на границах заданного диапазона,
±amax=±1g, было определено значение тока I0=3,5 mА, который требуется подавать в обмотку датчика момента от источника для обеспечения режима автоколебаний. Такой макет является прототипом предлагаемому устройству.
Статические испытания макета устройства, проведенные в гироскопической лаборатории 2 группы научно-исследовательской ВКА им. А.Ф.Можайского в соответствии со штатными методиками (акт испытаний приведен в приложении), позволили определить основные технические характеристики: среднее значение коэффициента преобразования
Figure 00000006
, относительную нестабильность коэффициента преобразования σKaΣ%=4,28%, нелинейность ГХ γГХ=3,46%, погрешность смещения нуля на выходе ΔN0=-0,00132 [б/р], погрешность смещения нуля на входе Δa0=-0,0133 [м/с2], нестабильность погрешности смещения нуля σaх0=1,7*10-5 [м/с2] в виде СКО. Для построения ГХ и определения ее нелинейности получение значений выходных сигналов проводилось на диапазоне ±1 g с шагом 0,1 g. Таким образом, для набора необходимой статистики потребовалось проведение 500 измерений в 20 точках с использованием наклонно-поворотного основания и оптической делительной головки. Для определения погрешности смещения нуля проведены еще 2 серии испытаний по 500 измерений при горизонтальном положении измерительной оси акселерометра. При этом суммарное время измерений составило Тiq≈pqTi≈110 с, где p - количество точек измерения, q=500 - количество периодов автоколебаний в каждой точке. Ti≈0,01 с - длительность одного периода, а, учитывая время, необходимое для подготовки оборудования к получению результатов в каждой точке измерения, суммарное время измерений составило Ти≈3000 с.
Полученные результаты калибровки прототипа позволили построить расчетную модель предлагаемого изобретения, при прочих равных условиях, и в соответствии с формулой изобретения: выбрать границы диапазона смещения нуля
±ai0=±K0amax=1×10-3 g, коэффициент кратности этого диапазона К0=ai0/amax=1×10-3 и величину тока
I00=K0I0=3,5 мкА в обмотке датчика момента, которая обеспечивает измерение ускорения на диапазоне смещения нуля, когда измерительная ось акселерометра горизонтальна. Расчетным граничным значениям диапазона смещения нуля с учетом погрешности смещения нуля на входе, авх1=+аi0+Δа0 и авх2=-ai0+Δа0 на выходе соответствуют значения выходных сигналов Nmах01=Ka01(+ai0+Δa0)=0,8873433 [б/р], -Nmах02=Ka02(-ai0+Δa0)=-0,9039534 [б/р], при помощи которых определяют значения коэффициента преобразования Ka0=[Nmax01-(-Nmax02)]/2ai0=91,29953 [c2/м], значения смещения нуля на выходе ΔN00=(Nmax01+(-Nmax02))/2=-0,008305076 [б/р] и на входе
Δa00=ΔN00/Ka0=-0,91*10-4 [м/с2].
Погрешность смещения нуля компенсируют на диапазоне смещения нуля непосредственным образом до уровня остаточной δΔa00 при условии аi0=0 и при помощи регулятора с измерителем тока компенсации, которые подключают к дополнительной обмотке датчика момента. При этом ток компенсации регулируют, пока на выходе не появится минимальное значение δΔN00a0δΔa00≈δΔN0103=-0,83*10-5 [б/р], где δΔa00≈-1/*10-6 [м/с2] - значение скомпенсированной погрешности смещения нуля на входе на диапазоне смещения нуля, то есть на три порядка больше, чем при измерении, после чего оставляют включенной цепь компенсации. Для определения значения коэффициента преобразования и погрешностей проводят калибровку акселерометра путем измерения ускорения известной величины, например, на границах заданного диапазона с учетом скомпенсированной погрешности смещения нуля, повторяют вычисления с использованием полученных значений выходных сигналов Nmахi1=Kai1(+amax+δΔa0)=0,895648334, -Nmaxi2ai2(-amax+δΔa0)=-0,895648351, вычисляют арифметическую (Nmaxi1+(-Nmaxi2))/2=КaδΔа0 и алгебраическую (Nmaxi1-(-Nmaxi2))/2amax=Ka суммы полученных сигналов, определяют стабильное значение коэффициента преобразования Кa=0,09129953 [с2/м], погрешности смещения нуля на выходе
δΔN0aδΔa0=(Nmaxi1+(-Nmaxi2))/2=-0,83*10-8 [б/p] и на входе δΔa0=δΔN0/Ka=-0,91*10-7 [м/c2]. Стабильность полученных погрешностей вычисляют при помощи отношения измеренных и вычисленных алгебраической и арифметической сумм (Nmaxi1-Nmахi2)/(Nmахi1+Nmахi2)=KаδΔa0/Kaamaх, где δΔa0/amах - относительная величина погрешности на заданном диапазоне, а представленная в процентах δΔa0100%/аmaxδΔa0=-9,27*10-7% характеризует стабильность измеренной погрешности и стабильность ГХ одновременно, причем при стабильном коэффициенте преобразования Кa=0,09129953=const, что подчеркнуто присутствием указанного коэффициента в числителе и знаменателе приведенного отношения. Сравнительные результаты, характеризующие известный и предлагаемый способы, сведены в таблицу.
Новизна предложения не следует явным образом из известного уровня техники, обеспечивает изобретательский уровень данного изобретения, которое может быть использовано для измерения линейных ускорений в системах инерциальной навигации самолетов, ракет, кораблей, космических аппаратов и других ПО.
Таким образом, повышение точности измерения ускорения и уменьшение нестабильности ТХ: погрешности смещения нуля, коэффициента преобразования, ГХ в предлагаемом способе и устройстве, удалось осуществить путем повышения стабильности ТХ: коэффициента преобразования, смещения нуля и ГХ, определяемых при калибровке, за счет одновременного определения коэффициента преобразования, смещения нуля и ГХ. Это стало возможным за счет уменьшения погрешности смещения нуля путем ее измерения на узком диапазоне - диапазоне смещения нуля, кратном заданному, с последующей компенсацией погрешности непосредственным образом до уровня остаточной, а также с помощью короткого и однократного измерения двух известных величин ускорения, соответствующих границам заданного диапазона, что привело к упрощению способа измерения, к сокращению материальных и временных затрат.
Таблица
ТХ, полученные по существующей методике и предложенным способом
ТХ, полученные по существующей методике ТХ, полученные предложенным способом
ΔN0, [б/р] -0,00132 δΔN0, [б/р] -0,83*10-8
Δa0, [м/с2] -0,0133 δΔa0, [м/c2] -0,91*10-7
Figure 00000007
, [c2/м]
0,09871999 Ka0, [с2/м] 91,29953
Кa, [с2/м] 0,09129953
σKaΣ, [%] 4,28 γδΔa0, [%] -9,27*10-7%
γГХ, % 3,46
σaвх0, [м/с2] 1,7*10-5
суммарное время измерений Тu, с ≈3000 суммарное время измерений Тu, с 20
суммарное время вычислений Tв ≈200 суммарное время вычислений Тв, с ≈20
суммарное время TΣ=Tu+Tв, с ≈3200 суммарное время TΣ=Tu+Tв, с ≈40
Получено стабильное значение коэффициента преобразования Kа=0,09129953, при одинаковой стабильности погрешности смещения нуля на входе и ГХ на уровне относительной ошибки δΔa0100%/аmaxδΔa0=-9,27*10-7%.
Для осуществления предлагаемого способа достаточно иметь горизонтируемое виброустойчивое основание и уровень с ценой деления не более 10''. Уменьшаются суммарные временные затраты TΣ=Tu+Tв, необходимые для выполнения измерений, приблизительно в 80 раз.
Источники информации
1. Коновалов С.Ф. и др. Гироскопические системы. Ч.3. (Акселерометры, датчики угловой скорости и др.) - М: ВШ, 1980, стр.4-7, 41-46.
2. А.Е.Синельников. Низкочастотные линейные акселерометры. Методы и средства поверки и градуировки. - М: И-во стандартов. 1979, с 8, 11, 15.
3. ГОСТ 18955-73. Акселерометры низкочастотные линейные. Термины и определения.
4. Жуков В.Н., Рыбаков В.И., Хегай Д.К. Принципы построения высокочувствительных миниатюрных датчиков систем управления МКА. // Изв. вузов приборостроение. 2004, №3, с.36.
5. Скалон А.И. Акселерометр с импульсной обратной связью. А.с. №794541, 07.01.81, бюл. №1, G01Р 15/08.
6. Кутуров А.Н., Кулешов В.В. Преобразователь ускорений с относительным цифровым кодом. // Изв. вузов приборостроение. 2003, №9, с.46.
7. Скалон А.И. Обобщенный анализ характеристик прецизионных датчиков механических величин, работающих в режиме автоколебаний. // Измерительная техника. - 1990. - С.7-9.
8. Кан В.Л. К вопросу об оценке погрешностей сложных приборов (комплектов). Исследования по методике оценки погрешностей измерений. Труды институтов комитета. Вып.57 (117). Стандартгиз, М.-Л., 1962, с.7-9.

Claims (2)

1. Способ измерения линейного ускорения на основе свойств автоколебательного акселерометра, заключающийся в том, что ускорение преобразуют в отклонение подвижной системы акселерометра с последующим преобразованием отклонения в электрический сигнал, его усиливают и преобразуют в последовательность временных интервалов, сформированных с помощью источника стабилизированного тока и элементов нелинейного звена в виде широтно-модулированного сигнала, который преобразуют в момент импульсной отрицательной обратной связи при помощи датчика момента, а сигналы импульсной отрицательной обратной связи преобразуют в сигнал, модулированный счетными импульсами, поступающими от генератора импульсов стабилизированной частоты fc, с последующим измерением количества счетных импульсов n1i и n2i, которые преобразуют в выходные сигналы акселерометра в виде разности Δni=n1i-n2i, суммы ni=n1i+n2i=fсТi и отношения количества импульсов Ni=Δni/ni, причем величину и знак входного воздействия aвх=ai+Δai, соответствующего измеряемым ускорению аi и погрешностям Δai, определяют по окончании каждого i-того периода Тi автоколебаний с помощью выходного сигнала и по формуле Ni=Ki(ai+Δai), (i=1…∞), где Ki=ml/KдмI0 - коэффициент преобразования входного воздействия, ml - маятниковый момент подвижной системы, Кдм и I0 - коэффициент передачи датчика момента и поступающий от стабилизированного источника ток, величина которого ограничивает диапазон измеряемого ускорения, при этом коэффициент преобразования Кi, погрешность, например смещение нуля Δai=Δa0, их стабильность, определяют при помощи формулы Nik=Kaik(ai+Δa0), градуировочную характеристику и ее нестабильность (нелинейность) при помощи другой формулы Nij=Kaij(aij+Δa0j), определяют заранее и каждую раздельно по результатам калибровки в условиях гравитационного поля Земли при вертикальном и горизонтальном положениях измерительной оси акселерометра, причем измерения каждой технической характеристики выполняют на интервале времени сотен периодов автоколебаний, чтобы применить статистическую обработку и фильтрацию результатов измерений, которые используют при измерении ускорения, отличающийся тем, что на основе заданного значения линейности градуировочной характеристики, полученного и паспортизованного в результате калибровки акселерометра на заводе-изготовителе в условиях его прецизионного оборудования и заданного диапазона измеряемого ускорения, на начальном этапе проведения калибровки выбирают границы диапазона смещения нуля ±ai0=±K0amax, которые больше величины измеренной погрешности смещения нуля Δa0, но меньше величин границ заданного диапазона ±ai0<<±amax, и кратны им, определяют коэффициент кратности этого диапазона К0i0/amax<<1, повышают величину коэффициента преобразования акселерометра на этом диапазоне, выполняют калибровку акселерометра на диапазоне смещения нуля в условиях гравитационного поля Земли при горизонтальном положении измерительной оси акселерометра путем измерения ускорения на границах диапазона смещения нуля, при этом коэффициент преобразования, погрешность смещения нуля и их стабильность, соответствующие диапазону смещения нуля определяют одновременно при коротком по времени однократном измерении двух известных величин входного ускорения, соответствующих границам диапазона смещения нуля +aвx01=+g·sinα+Δa0=+0,001g+Δa0 и -aвx01=-g·sinα+Δa0=-0,001g+Δa0, где α - угол между измерительной осью прибора и
Figure 00000008
, полученные значения выходного сигнала
Nmax01=Ka01(+ai0+Δa0) и -Nmax02a02(-ai0+Δa0) суммируют алгебраически и делят на удвоенное значение ускорения границы диапазона смещения нуля (Nmax01-(-Nmax02))/2ai0=Ka0, а затем суммируют арифметически и делят пополам (Nmax01+(-Nmax02))/2=Ka0Δa0, определяют стабильное значение коэффициента преобразования
Кa0 и значение погрешности смещения нуля на выходе ΔN00=Ka0Δa00, погрешность смещения нуля на входе компенсируют непосредственно до уровня остаточной δΔа00 на входе и δΔN00 на выходе, повторяют калибровку акселерометра на заданном диапазоне в условиях гравитационного поля Земли при вертикальном положении измерительной оси акселерометра путем короткого по времени однократного измерения двух известных величин входного ускорения, например, на границах заданного диапазона с учетом скомпенсированной погрешности смещения нуля aвх1=+amax+δΔa0, aвx2=-amax+δΔa0, полученные значения выходных сигналов Nmax1=Ka1(+amax+δΔa0) и -Nmax2a2(-amax+δΔa0) суммируют алгебраически и делят на удвоенное значение ускорения границы заданного диапазона (Nmax1-(-Nmax2))/2amax=Ka, а затем суммируют арифметически и делят пополам (Nmax1+(-Nmax2))/2=KaδΔa0, определяют стабильное значение коэффициента преобразования для заданного диапазона Ка, остаточную погрешность смещения нуля на выходе δΔN0=KaδΔa0 и на входе
δΔa0=δΔN0/Ka, стабильность полученных значений вычисляют при помощи отношения измеренных и вычисленных алгебраической и арифметической сумм
(Nmax1-Nmax2)/(Nmax1+Nmax2)=KaδΔa0/Kaamax=δΔa0/amax, где δΔa0/amax - относительная величина скомпенсированной погрешности на заданном диапазоне, а представленная в процентах δΔa0100%/amaxδΔa0 характеризует стабильность измеренной погрешности и стабильность ГХ одновременно, причем при стабильном коэффициенте преобразования Ka=const, полученные значения запоминают и используют при формировании полетного задания для измерения ускорения в заданном диапазоне с помощью выходного сигнала устройства и по формуле Nаiаi+δΔа0) при меньшем значении остаточной погрешности смещения нуля δΔa0<<Δа0.
2. Устройство для измерения ускорения, содержащее автоколебательный акселерометр, в состав которого входят: инерционный элемент, датчики угла и момента, механически связанные друг с другом, выход датчика угла электрически связан через усилитель-преобразователь и компаратор с входом электронного переключателя, через контакты которого обмотка датчика момента связана через регулятор и измеритель с источником тока, в состав устройства также входят: генератор счетных импульсов стабилизированной частоты, выход которого через контакты электронного переключателя связан со счетными входами реверсивного и суммирующего счетчиков импульсов, а их выходы, в свою очередь, связаны с информационным входом спецвычислителя, на выходе которого получают информацию о величине и знаке измеряемого ускорения, отличающееся тем, что введены: блок управления, дополнительный регулятор с измерителем тока диапазона смещения нуля, дополнительный регулятор с измерителем тока компенсации, в датчик момента введена дополнительная обмотка, входы регуляторов электрически связаны с источником стабилизированного тока, выход регулятора тока компенсации через измеритель тока и нормально разомкнутые контакты блока управления связан с дополнительной обмоткой датчика момента, а на выходе устройства одновременно определяют основные ТХ: выходной параметр Ni, стабильный коэффициент преобразования Ка, остаточную погрешность смещения нуля δΔа0 и относительную скомпенсированную погрешность смещения нуля γδΔa0 - для дальнейшего их учета в моделях погрешностей для систем управления движением и навигации подвижных объектов.
RU2008110986/28A 2008-03-21 2008-03-21 Способ измерения линейного ускорения в системах управления ракет и устройство для его осуществления RU2362173C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2008110986/28A RU2362173C1 (ru) 2008-03-21 2008-03-21 Способ измерения линейного ускорения в системах управления ракет и устройство для его осуществления

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2008110986/28A RU2362173C1 (ru) 2008-03-21 2008-03-21 Способ измерения линейного ускорения в системах управления ракет и устройство для его осуществления

Publications (1)

Publication Number Publication Date
RU2362173C1 true RU2362173C1 (ru) 2009-07-20

Family

ID=41047278

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2008110986/28A RU2362173C1 (ru) 2008-03-21 2008-03-21 Способ измерения линейного ускорения в системах управления ракет и устройство для его осуществления

Country Status (1)

Country Link
RU (1) RU2362173C1 (ru)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9008988B2 (en) 2009-10-26 2015-04-14 Leica Geosystems Ag Method of calibrating inertial sensors
RU2654941C1 (ru) * 2017-05-31 2018-05-23 Публичное акционерное общество "Московский институт электромеханики и автоматики" (ПАО "МИЭА") Способ цифровой фильтрации шумовой составляющей в инерциальных датчиках
CN112379126A (zh) * 2020-11-02 2021-02-19 西安交通大学 一种复合量程石英谐振加速度传感器

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ГОСТ 18955-73. АКСЕЛЕРОМЕТРЫ НИЗКОЧАСТОТНЫЕ ЛИНЕЙНЫЕ. ТЕРМИНЫ И ОПРЕДЕЛЕНИЯ. *
КУТУРОВ А.Н. и др. Преобразователь ускорений с относительным цифровым кодом. Изв. вузов. Приборостроение, 2003, т.46, №9, с.34-37. *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9008988B2 (en) 2009-10-26 2015-04-14 Leica Geosystems Ag Method of calibrating inertial sensors
RU2654941C1 (ru) * 2017-05-31 2018-05-23 Публичное акционерное общество "Московский институт электромеханики и автоматики" (ПАО "МИЭА") Способ цифровой фильтрации шумовой составляющей в инерциальных датчиках
CN112379126A (zh) * 2020-11-02 2021-02-19 西安交通大学 一种复合量程石英谐振加速度传感器

Similar Documents

Publication Publication Date Title
Kuhnen et al. Adaptive inverse control of piezoelectric actuators with hysteresis operators
CN101484826A (zh) 使用渡越时间原理进行光电子非接触距离测量的方法和设备
RU2362173C1 (ru) Способ измерения линейного ускорения в системах управления ракет и устройство для его осуществления
Bluemm et al. Correcting nonlinearity and temperature influence of sensors through B-spline modeling
RU2366961C1 (ru) Способ измерения линейного ускорения
JP2645464B2 (ja) エレベータの位置制御装置の調整方法
RU2361223C1 (ru) Способ измерения линейного ускорения
CN113224630B (zh) 一种延时光路锁定装置及锁定方法
CN112284366B (zh) 一种基于tg-lstm神经网络的偏振光罗盘航向角误差校正方法
US4309103A (en) Device for contact-free thickness or interval measurement
Grebennikov et al. Method to increase stability of the scale factor of a precision pendulum accelerometer with digital feedback
Grebennikov et al. Pendulum accelerometer with digital control and new functionality
Wolf Undulator field integral measurements
Jacobs et al. Development of a flexible link setup for an advanced linear control theory course
Nikiforov et al. Identification of pendulous accelerometer mathematical model taking into account parametric uncertainty
EP4130661A1 (en) Vibration-type angular velocity sensor
US11719597B2 (en) Measuring system and method for measuring the displacement of at least one point of a bridge
RU2334198C1 (ru) Способ инерциальной навигации и устройство для его осуществления
EP4339621A1 (en) Mems accelerometer systems
KR102304429B1 (ko) 관성항법장치의 교정 시스템 및 교정 방법
Chai et al. Predictive Gradient Based Control Using Hammerstein Model for MEMS Micromirrors
CN117713759A (zh) 脉冲指令滤波方法、装置、终端设备以及存储介质
Cherkasova et al. Implementation Options of Control Systems for a Cantilever Calibration Device
SU600404A1 (ru) Способ измерени давлени
UA150147U (uk) Канал автоматичного супроводження літальних апаратів за напрямком з можливістю формування і обробки їх зображення та кібернетичним захистом інформації

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20100322