Оценка характеристик аэрометрических систем Текст научной статьи по специальности «Физика»

УДК 629.7.051

ОЦЕНКА ХАРАКТЕРИСТИК АЭРОМЕТРИЧЕСКИХ СИСТЕМ

© 2017 В.В. Алмазов1,2, Н.Н. Макаров2, М.Ю. Сорокин12

1 Ульяновский государственный технический университет 2 АО «Ульяновское конструкторское бюро приборостроения»

Статья поступила в редакцию 15.03.2017

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

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

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

На данный момент наибольшее распространение измерения высотно-скоростных параметров (ВСП) получил аэродинамический метод,

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

Рис. 1. Условная схема измерения высотно-скоростных параметров

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

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

Алмазов Василий Вячеславович, аспирант. E-mail: almaz73.91@mail.ru

Макаров Николай Николаевич, доктор технических наук, советник генерального директора. E-mail: rto@ukbp.ru Сорокин Михаил Юрьевич, кандидат технических наук, начальник отдела, доцент кафедры «Измерительно-вычислительные комплексы». E-mail: rto@ukbp.ru

сигналов [1] к интегрированным комплексам ВСП [5], в которык привлекается информация от бортовой инерциальной навигационной системы и ряд давлений с фюзеляжа ЛА, и системам воздушных сигналов на базе многофункциональных приемников воздушнык давлений (например, на рис. 2 приведена предлагаемая компанией Goodrich система воздушнык сигналов на базе SmartProbe), в которык исключены пневмотрассы и отдельные вычислители за счет объединения в едином корпусе приемника воздушнык давлений, датчиков давлений и вычислителей. Направление создания подобнык мно-гофункциональнык измерителей аэродинамических параметров на данный момент активно прорабатывается ведущими разработчиками аэрометрических систем, и вопросы метрологического характера приобретают другой характер. Целесообразно рассмотреть вопросы совместной градуировки приемника воздушных давлений и датчиков давлений, в особенности для архитектур, похожих на систему воздушнык сигналов на базе SmartProbe.

Рис. 2. Сравнение двух систем воздушных сигналов традиционной архитектуры (слева) и предлагаемой компанией Goodrich на базе многофункциональных приемников SmartProbe (справа)

Рассмотрим погрешность измерения ВСП в зависимости от восприятия воздушных давлений различными зондовыми средствами (приемниками воздушных давлений) и погрешности применяемых датчиков давления в дозвуковом диапазоне скоростей и до высоты, равной 11 км. В соответствии с [24] абсолютная барометрическая (геопотенциальная) высота Яабс вычисляется по следующей зависимости от статического давления на высоте Рст:

Набс —

3, 5 32 7467 -Р.

0,1902631

0,000079691

(1)

где Рст - статическое давление, [мм рт.ст.].

Приборная скорость Упр в зависимости от динамического давления Рдин вычисляется по следующей формуле:

Упр — 2 7 3 9,3 1 3 79 9

м

(1 + Щ V:

V 760/

3,5

- 1

(2)

А Иабс — JT^A Рст>

и г cm

(3)

Частная производная абсолютной высоты:

д набс д Рст

Рн'9

R-Th Рст' 9

(4)

где д =9,80665 - ускорение свободного падения, [м/с2]; Тн - температура воздуха на высоте Набс, [К]; Я = 287,05287- газовая постоянная, [Дж/кгК].

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

AVnp— dnV АРП2+АРС

(5)

Частная производная приборной скорости описывается следующим выражением:

dVm

д V„r

10,57959

дРн дРп

( 1

\ 760 )

(6)

На рис. 3 и в табл. 1 приведены погрешности измерения абсолютной барометрической высоты Навс в зависимости от погрешности измерения статического давления, которая изменяется в диапазоне от 0.1 мм рт. ст. до 0,5 мм рт. ст.

Сложнее ситуация с погрешностью измерения приборной скорости 1^р, которая зависит от погрешности измерения как полного давления Рп, так и статического давления Рст. На рис. 4 и 5 приведены погрешности измерения приборной скорости в зависимости от высоты Набс и скорости полета, выраженной через число Маха.

где Рдин=Рп - Рст , [мм рт.ст.].

Переходя к оценке абсолютной погрешности измерения абсолютной барометрической высоты /1Набс в зависимости от погрешности измерения статического давления АРст:

Рис. 3. Погрешность измерения абсолютной барометрической высоты Набс

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

пр

Таблица 1. Погрешность измерения абсолютной барометрической высоты Hабс

Набс, м Рст, мм рт. ст. ДНабс, м

АРст=0,1 мм рт. ст. ДРст=0,2 мм рт. ст. ДРст=0,3 мм рт. ст. ДРст=0,4 мм.рт. ст. ДРст=0,5 мм.рт. ст.

-500 806,15 1,06 2,12 3,17 4,23 5,29

-0 760,00 1,11 2,22 3,33 4,44 5,55

600 707,47 1,18 2,35 3,53 4,70 5,88

3000 525,86 1,50 2,99 4,49 5,98 7,48

7000 307,98 2,31 4,61 6,92 9,22 11,53

10000 198,29 3,29 6,59 9,88 13,18 16,47

через который происходит удаление посторонних частиц (вода, песок и т.п.) из камеры торможения. Однако наличие дренажа приводит к уменьшению величины полного давления [8], что приводит к появлению дополнительной погрешности измерения. На рис. 6 приведена погрешность измерения приборной скорости в зависимости от диаметра дренажного отверстия при различных скоростях набегающего потока для одного из макетов приемника полного давления. Видно, что зависимость носит нелинейный характер и описывается известными выражениями движения воздуха с учетом числа Рейнольдса. Таким образом, незначительное изменение диаметра дренажного отверстия в пределах от 0,7 до 1 мм приводит к существенным изменениям погрешности измерения полного давления (приборной скорости в данном случае). И это изменение только одного из параметров, характеризующих метрологические характеристики приемника полного давления.

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

Скорость набегающего потока, км/ч

Рис. 6. Погрешность измерения приборной скорости в зависимости от диаметра дренажного отверстия

Рис. 4. Погрешность измерения приборной скорости на Набс = 0 м

Рис. 5. Погрешность измерения приборной скорости на Набс = 10000 м

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

Итак, показано, что погрешность измерения ВСП может превышать допустимую уже за счет искажения восприятия воздушный давлений самим средством восприятия (приемником полного давления, многофункциональным измерителем аэродинамических параметров). Однако, следующим элементов в цепи преобразования, см. рис. 1, является датчик давления. Современные датчики давления для измерения ВСП имеют погрешность измерения до 0,01% от рабочего диапазона давления [9]. Выходной сигнал датчика давления зависит от многих факторов, таких как температура, плотность среды, механические внешние воздействующие факторы, напряжение питания и другие. Наибольшим вкладом в погрешность измерения вносит погрешность, вызванная изменением температуры. Для датчиков тензометрического типа, погрешность измерения, вызванная изменением температуры, составляет до 30% от рабочего диапазона давления. Менее чувствительными к изменению температуры являются вибрационно-частотные датчики на основе вибрирующего цилиндра. Их погрешность, вызванная изменением температуры, составляет до 7% от рабочего диапазона давления. Для уменьшения погрешности, вызванной изменением температуры, применяются различные методы.

Наиболее очевидным является использование датчика давления при одном значении температуры. Фиксация температуры производится путем помещения чувствительного к давлению элемента в термостат. Изоляционный материал позволяет задать внутри него значение температуры, отличное от внешней среды. Для работы на температурах от минус 60°С до 75°С требуется поддерживать чувствительный элемент как минимум при температуре на 5°С выше верхнего значения. Схема управления обеспечивает поддержание высокой температуры при помощи датчика температуры и нагревательного элемента. Недостатком данного метода является то, что термостат имеет сравнительно большие габариты. К тому же при отрицательных температурах для создания температуры внутри термостата более чем 80°С требуется большая энергия, что увеличивает потребляемую мощность датчика давления. Для беспилотныгх ЛА габаритные размеры и потребляемая мощность имеют важную роль, поэтому необходим способ, исключающий применение термостата.

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

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

Рис. 7. Схема установки для измерения градуировочных данных датчика давления

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

Для расчета давления по полученным данным X, У (значение давления и значение температуры) необходимо определить, в какой квадрат попадают данные значения. При этом фиксируются значения У1, У2 между которыми попало значение температуры и значения XI, Х2 между которыми попало значение давления. Точки пересечения столбцов и строк являются вершинами квадрата со значениями 011, 012, 021, 022. Далее по полученным значениям и значениям из таблицы рассчитывается значение давления используя линейную интерполяцию по двум осям.

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

(Х2-ХШ2-уП) (Х2-ХШ2-П)

Рис. 8. Схематичное изображение вычисления действующего давления

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

Р = ! ^¡Сц-Т 1-и) (7)

где Р - измеренное значение давления (вычисленное), [мм рт.ст.]; Т - значение периода выходного сигнала по давлению, [мкс]; и - значение напряжения датчика температуры, [мВ]; С, - коэффициенты аппроксимирующей функции.

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

Представленные методы компенсации погрешности позволяют минимизировать значение погрешности, вызванное изменением температуры. Однако погрешность измерения 0,01% от рабочего диапазона давления достичь удается не на всех датчиках. Для датчиков давления на основе вибрирующего цилиндра присущ недостаток в виде, так называемого, «перескока частоты». Суть явления заключается в скачкообразном изменении выходного

сигнала по давлению при переходе через критические точки давления. На рис. 9 представлена идеальная амплитудно-частотная характеристика (АЧХ) преобразователя давления.

Рис. 9. АЧХ преобразователя давления

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

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

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

Рис. 10. АЧХ преобразователя давления ниже точки перескока частоты

Рис. 11. АЧХ преобразователя давления выше точки перескока частоты

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

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

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

1. Кравцов, В.Г. Анализ современных бортовых средств измерения высотно-скоростных параметров / В.Г. Кравцов, Н.В. Алексеев// Приборы и системы. Управление, контроль, диагностика. 2000. №8. С. 47-50.

2. ГОСТ 4401-81. Атмосфера стандартная. Параметры. -М.: Изд-во стандартов, 1982. 182 с.

3. ГОСТ 3295-75. Таблицы гипсометрические для геопотенциальных высот до 50000 м. Параметры. - М.: Изд-во стандартов, 1974. 75 с.

4. ГОСТ 5212-74. Таблица аэродинамическая. Динамические давления и температуры торможения воздуха для скорости полета от 10 до 4000 км/ч. Параметры. -М.: Изд-во стандартов, 1974. 239 с.

5. Олейник, А.И. Архитектура перспективного информационного комплекса высотно-скоростных параметров полёта самолёта // Перспективы науки. 2009. №1. С. 36-38.

6. Gracey, W. Wind-tunnel investigation of a number of total-pressure tubes at high angles of attack subsonic, transonic and supersonic speeds //. Technical Report 1303. 1956.

7. Петунин, А.Н. Методы и техника измерений параметров газового потока (приемники давления и скоростного напора). - М.: «Машиностроение», 1972. 332 с.

8. Пушков, С.Г. Погрешность восприятия полного давления приемником с протоком в камере торможения / С.Г. Пушков, Е.Г. Харин, В.Р. Кожурин, Л.Л Ловицкий // Авиакосмическое приборостроение. 2010. №5. С. 3-8.

9. Соколов, Л.В. Отечественные инновационные разработки датчиков давления / Л.В. Соколов, М.Ю. Сорокин // Научно-технический сборник «Оборонная техника». 2015. №11-12. С. 87-92.

ASSESSMENT OF AERO METRIC SYSTEMS CHARACTERISTICS

© 2017 V.V. Almazov1'2, N.N. Makarov2, M.Yu. Sorokin1'2

1 Ulyanovsk State Technical University 2 JSC "Ulyanovsk Design Bureau of Instrument Manufacturing

The paper describes the main sources of measurement error of altitude and speed parameters and the ways to reduce measurement uncertainty of altitude and speed parameters, mainly due to the combined calibration of the air pressure probes and pressure sensors.

Key words: pressure sensor, air pressure probe, aero metric system, measurement error, multi-function meter

Vasiliy Almazov, Post-graduate Student. E-mail: almaz73.91@mail.ru; Nikolay Makarov, Doctor of Technical Sciences, Advisor of General Director. E-mail: rto@ukbp.ru; Mikhail Sorokin, Candidate of Technical Sciences, Head of the Department, Associate Professor at the Department "Measuring and Computing Complexes". E-mail: rto@ukbp.ru