ГОСТ ISO 16063-11-2013. Межгосударственный стандарт. Вибрация. Методы калибровки датчиков вибрации и удара. Часть 11. Первичная вибрационная калибровка методами лазерной интерферометрии

(введен в действие Приказом Росстандарта от 05.12.2013 № 2181-ст)

ИСТОРИЯ ИЗМЕНЕНИЙ НОРМАТИВНО-ТЕХНИЧЕСКИХ АКТОВ

flowchart LR
  n0["ГОСТ ИСО 5347-1-96"]:::rst
  click n0 "https://bukva.ai/docs/rukovodstvo-rb-002-97" "ГОСТ ИСО 5347-1-96. Межгосударственный стандарт. Вибрация. Калибровка датчиков вибрации и удара. Часть 1. Первичная вибрационная калибровка методами лазерной интерферометрии (приня"
  n1["ГОСТ ISO 16063-11-2013"]:::doc
  click n1 "https://bukva.ai/docs/standart-gost-iso-16063-11-2013" "ГОСТ ISO 16063-11-2013. Межгосударственный стандарт. Вибрация. Методы калибровки датчиков вибрации и удара. Часть 11. Первичная вибрационная калибровка методами лазерной интерфером"
  n0 --> n1
classDef rst fill:#FFEBEB,stroke:#8B0000,color:#000;
classDef doc fill:#EBFFEB,stroke:#006400,color:#000;
classDef notEnt fill:#FFFFEB,stroke:#FF8C00,color:#000;

Введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 5 декабря 2013 г. № 2181-ст

МЕЖГОСУДАРСТВЕННЫЙ СТАНДАРТ

ВИБРАЦИЯ

МЕТОДЫ КАЛИБРОВКИ ДАТЧИКОВ ВИБРАЦИИ И УДАРА

ЧАСТЬ 11

ПЕРВИЧНАЯ ВИБРАЦИОННАЯ КАЛИБРОВКА МЕТОДАМИ ЛАЗЕРНОЙ ИНТЕРФЕРОМЕТРИИ

Vibration. Methods for the calibration of vibration and shock transducers. Part 11. Primary vibration calibration by laser interferometry

(ISO 16063-11:1999, Methods for the calibration of vibration and shock transducers — Part 11: Primary vibration calibration by laser interferometry, IDT)

ГОСТ ISO 16063-11-2013

МКС 17.160

Дата введения 1 ноября 2014 года

Предисловие

Цели, основные принципы и основной порядок проведения работ по межгосударственной стандартизации установлены в ГОСТ 1.0-2015 «Межгосударственная система стандартизации. Основные положения» и ГОСТ 1.2-2015 «Межгосударственная система стандартизации. Стандарты межгосударственные, правила и рекомендации по межгосударственной стандартизации. Порядок разработки, принятия, обновления и отмены»

Сведения о стандарте

1 ПОДГОТОВЛЕН Федеральным государственным унитарным предприятием Всероссийский научно-исследовательский институт метрологической службы (ФГУП «ВНИИМС») на основе собственного перевода на русский язык англоязычной версии стандарта, указанного в пункте 5, с участием Открытого акционерного общества «Научно-исследовательский центр контроля и диагностики технических систем» (АО «НИЦ КД»)

2 ВНЕСЕН Федеральным агентством по техническому регулированию и метрологии

3 ПРИНЯТ Межгосударственным советом по стандартизации, метрологии и сертификации (протокол от 3 декабря 2013 г. № 62-П)

За принятие стандарта проголосовали:

4 Приказом Федерального агентства по техническому регулированию и метрологии от 5 декабря 2013 г. № 2181-ст межгосударственный стандарт ГОСТ ISO 16063-11-2013 введен в действие в качестве национального стандарта Российской Федерации с 1 ноября 2014 г.

5 Настоящий стандарт идентичен международному стандарту ISO 16063-11:1999 «Методы калибровки датчиков вибрации и удара. Часть 11. Первичная вибрационная калибровка методами лазерной интерферометрии» («Methods for the calibration of vibration and shock transducers — Part 11: Primary vibration calibration by laser interferometry», IDT).

Международный стандарт разработан Подкомитетом ISO/TC 108/SC 3 «Применение и калибровка средств измерений вибрации и удара» Технического комитета по стандартизации ISO/TC 108 «Вибрация и удар» Международной организации по стандартизации (ISO).

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

Наименование настоящего стандарта изменено относительно наименования указанного международного стандарта для приведения в соответствие с ГОСТ 1.5 (подраздел 3.6).

Настоящий стандарт подготовлен на основе применения ГОСТ Р ИСО 16063-11-2009 <*>

<*> Приказом Федерального агентства по техническому регулированию и метрологии от 5 декабря 2013 г. № 2181-ст ГОСТ Р ИСО 16063-11-2009 отменен с 1 ноября 2014 г.

6 ВЗАМЕН ГОСТ ИСО 5347-1-96

7 ПЕРЕИЗДАНИЕ. Март 2019 г.

Информация об изменениях к настоящему стандарту публикуется в ежегодном указателе «Национальные стандарты», а текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ежемесячном информационном указателе «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет (www.gost.ru)

1 Область применения

Настоящий стандарт устанавливает три метода первичной вибрационной калибровки преобразователей прямолинейного ускорения (далее — акселерометров) совместно с усилителями или без них для определения комплексного коэффициента преобразования посредством возбуждения гармонической вибрации и измерения амплитуды колебаний методами лазерной интерферометрии.

Установленные методы применяют в диапазоне частот от 1 Гц до 10 кГц и в диапазоне амплитуд ускорения от 0,1 до 1000 м/с2 (в зависимости от частоты).

Неопределенность измерения в соответствии с данными методами указана в разделе 2. Метод синус-аппроксимации (метод 3) позволяет проводить калибровку на частотах ниже 1 Гц (например, на частоте 0,4 Гц, используемой в качестве опорной частоты в некоторых стандартах) с амплитудами ускорения менее 0,1 м/с2 (например, 0,04 м/с2 на частоте 1 Гц) при наличии соответствующего низкочастотного вибростенда (см. раздел 9).

Метод счета полос (метод 1) применяют для определения модуля коэффициента преобразования в диапазоне частот от 1 до 800 Гц и, в особых случаях, на более высоких частотах (см. раздел 7). Метод точек минимума (метод 2) применяют для определения модуля коэффициента преобразования в диапазоне частот от 800 Гц до 10 кГц (см. раздел 8). Метод синус-аппроксимации может быть применен для определения модуля и фазового сдвига коэффициента преобразования в диапазоне частот от 1 Гц до 10 кГц.

Методы 1 и 3 обеспечивают калибровку при фиксированных значениях амплитуд ускорения на разных частотах. Метод 2 обеспечивает калибровку для фиксированных значений амплитуд перемещений (амплитуда ускорения изменяется в зависимости от частоты).

2 Неопределенность измерения

Применение методов, установленных настоящим стандартом, должно обеспечивать неопределенность измерения, не превышающую следующие значения:

a) для модуля коэффициента преобразования:

• 0,5% измеряемого значения на опорных частотах и ускорениях (опорных точках),
• не более 1% измеряемого значения вне опорных точек;

b) для фазового сдвига коэффициента преобразования:

• 0,50 в опорных точках;
• не более 10 вне опорных точек.

Рекомендуются следующие опорные точки:

• частота: 160, 80, 40, 16 или 8 Гц (или круговая частота: 1000, 500, 250, 100 или 50 рад/с);
• ускорение (амплитуда или среднее квадратическое значение): 100, 50, 20, 10, 2 или 1 м/с2.

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

Примечание — Неопределенность измерения представляют в виде расширенной неопределенности в соответствии с [2].

3 Требования к испытательному оборудованию и средствам измерений

3.1 Общие положения

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

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

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

3.2 Генератор частоты с показывающим устройством

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

a) неопределенность измерения частоты не более 0,05% показываемого значения;

b) отклонение частоты в процессе измерений не более +/- 0,05% показываемого значения;

c) отклонение амплитуды сигнала в процессе измерений не более +/- 0,05% показываемого значения.

3.3 Вибростенд с усилителем мощности

Следует использовать вибростенд с усилителем мощности, имеющий следующие характеристики:

a) коэффициент гармоник по ускорению не более 2%;

b) малые поперечные и угловые колебания стола вибростенда, чтобы не оказывать существенного влияния на результаты калибровки. При больших амплитудах колебаний, преимущественно в низкочастотном диапазоне от 1 до 10 Гц, колебания в поперечном направлении должны быть не более 1% колебаний в основном направлении, в диапазоне от 10 Гц до 1 кГц — не более 10%; свыше 1 кГц — не более 20%;

c) собственный электронный шум испытательной установки должен быть не менее чем на 70 дБ ниже максимального значения выходного сигнала;

d) отклонение амплитуды ускорения в процессе калибровки — не более +/- 0,05% показываемого значения.

Поверхность стола вибростенда, на которую устанавливают акселерометр, не должна деформировать его основание.

3.4 Сейсмический блок (блоки) вибростенда и лазерного интерферометра

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

При использовании общего блока его масса должна быть по меньшей мере в 2000 раз больше массы подвижной системы вибростенда. При выполнении данного условия реактивная вибрация акселерометра и интерферометра не будут превышать 0,05%. Если масса сейсмического блока меньше, то колебания блока следует учитывать при расчете вибрации.

Для уменьшения искажений вследствие колебаний грунта при измерениях в диапазоне от 10 Гц до 10 кГц сейсмический блок (блоки) устанавливают на демпфированных пружинах, чтобы соответствующая составляющая неопределенности измерения не превышала 0,1%.

3.5 Лазер

Для калибровки используют гелий-неоновый лазер.

В лабораторных условиях (при атмосферном давлении 100 кПа, температуре 23 °C и относительной влажности 50%) номинальная длина волны лазера равна 0,63281 мкм.

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

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

3.6 Интерферометр

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

Необходимая ширина полосы частот f_max, Гц, может быть вычислена по амплитуде измеряемой скорости v_max, м/с, по формуле

f_max = v_max·3,16·106.

Для метода 1 (см. рисунок 1) и метода 2 (см. рисунок 2) используют простой интерферометр Майкельсона с одним фотодетектором. Для метода 3 (см. рисунок 3) используют модифицированный интерферометр Майкельсона с квадратурными выходными сигналами и с двумя фотодетекторами для приема сигнала интерферометра. Может быть использован интерферометр Майкельсона, модифицированный в соответствии с рисунком 4. Четвертьволновая пластина преобразует падающий линейно поляризованный свет в два измерительных луча со взаимно ортогональной поляризацией и фазовым сдвигом 90°. После интерференции с линейно поляризованным опорным лучом две компоненты со взаимно ортогональными поляризациями пространственно разделяют при помощи соответствующих оптических элементов (например, призмы Волластона или поляризационного светоделителя) и детектируют двумя фотодиодами.

Два выходных сигнала модифицированного интерферометра Майкельсона должны иметь сдвиг нуля не более +/- 5% амплитуды сигнала, относительную разность амплитуд менее +/- 5% и отклонение разности фаз от номинального угла 90° не более +/- 5°. Чтобы удовлетворить этим требованиям, должны быть предусмотрены средства регулировки нуля, амплитуды сигнала и сдвига фаз между двумя сигналами интерферометра.

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

Пример — Для амплитуды перемещения 2,5 мм (ей соответствует амплитуда ускорения 0,1 м/с2 на частоте 1 Гц) допуски могут быть расширены до +/- 10% для относительной разности амплитуд и до +/- 20° для отклонения сдвига фаз от номинального угла 90° (см. также примечание 1 к 9.2).

Примечание - Модифицированный интерферометр Майкельсона для методов 1, 2 или 3 может быть заменен другим подходящим двухлучевым интерферометром, например модифицированным интерферометром Маха — Цандера.

3.7 Аппаратура для счета интерференционных полос (метод 1)

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

a) диапазон частот от 1 Гц до требуемой максимальной частоты (обычно 20 МГц);

b) максимальная неопределенность измерения 0,01% показываемого значения.

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

3.8 Перестраиваемый полосовой фильтр или спектроанализатор (метод 2)

Используемая аппаратура должна иметь следующие характеристики:

a) диапазон частот от 800 Гц до 10 кГц;

b) ширина полосы анализа не более 1/6 октавы (12% среднегеометрической частоты);

c) спад частотной характеристики фильтра более 24 дБ на октаву;

d) отношение «полезный сигнал/шум» более 70 дБ (определяют по максимальному значению полезного сигнала);

e) динамический диапазон более 60 дБ.

3.9 Устройство детектирования нуля (метод 2)

Устройство детектирования нуля (не является необходимым при использовании спектроанализатора) должно иметь диапазон частот от 800 Гц до 10 кГц. Диапазон должен быть достаточен для детектирования выходного шума полосового фильтра.

3.10 Устройство измерения истинного среднеквадратичного значения выходного сигнала акселерометра

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

a) диапазон частот от 1 Гц до 10 кГц;

b) неопределенность измерения 0,1% показываемого значения.

Для вычисления амплитуды сигнала среднее квадратическое значение умножают на \sqrt{2} .

Вольтметры среднеквадратичного значения применяют в методах 1 и 2. В методе 3 для измерения напряжения следует использовать специальную аппаратуру (см. 3.13). Вольтметр среднеквадратичного значения в данном случае может быть использован как вспомогательное средство.

3.11 Устройство измерения коэффициента гармоник

Для измерения искажений должно использоваться устройство, позволяющее определить коэффициент гармоник в диапазоне до 5% и имеющее следующие характеристики:

a) диапазон частот от 1 Гц до 10 кГц с возможностью измерения до 5-й гармоники;

b) неопределенность измерения 10% показываемого значения в диапазоне коэффициента гармоник от 0,5% до 5%.

3.12 Осциллограф

Применение осциллографа не является обязательным, но он может быть использован для лучшей настройки интерферометра и контроля формы сигнала акселерометра. Диапазон частот осциллографа — от 1 Гц до как минимум 2 МГц.

3.13 Записывающее устройство с компьютерным интерфейсом (метод 3)

Метод 3 требует применения записывающего устройства с компьютерным интерфейсом, позволяющего выполнять аналого-цифровое преобразование и хранить в памяти два квадратурных выходных сигнала интерферометра и выходной сигнал акселерометра. Разрешающая способность по амплитуде, частота выборки и объем памяти должны быть достаточными для калибровки в требуемом диапазоне амплитуд с неопределенностью измерения, указанной в разделе 2. Обычно для выходного сигнала акселерометра достаточно разрядности аналого-цифрового преобразования 10 бит, а для квадратурных сигналов на выходе интерферометра — 8 бит. Возможно использование двух регистрирующих устройств: двухканального для выходных сигналов интерферометра и устройства с более высокой разрешающей способностью и более низкой частоты выборки для выходного сигнала акселерометра. Оба регистрирующих устройства должны начинать и заканчивать запись сигналов одновременно. Допустимое расхождение в моментах начала записи следует оценить с учетом требований к неопределенности измерения (см. раздел 2).

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

Пример - Для калибровки акселерометра на частоте вибрации 1 Гц и при амплитуде ускорения 0,1 м/с2 требуются объем памяти более 4 Мбайт и частота выборки более 512 кГц.

3.14 Компьютер с программой обработки данных (метод 3)

Компьютер с программой обработки данных используют при расчетах по 9.4.

3.15 Другие требования

Чтобы удовлетворить требованию к неопределенности измерения не превышать 0,5%, рекомендуется акселерометр калибровать совместно с усилителем.

При расчете неопределенности измерения (см. приложение A) должны быть приняты во внимание чувствительность акселерометра к деформации основания, коэффициент поперечного преобразования и стабильность характеристик системы «акселерометр — усилитель».

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

Если измеряют вибрацию верхней поверхности имитатора, то она должна быть отполирована, и пятно лазерного луча должно находиться вблизи геометрического центра этой поверхности. В случаях, когда имитатор нельзя считать абсолютно твердым телом, следует учитывать относительное движение верхней и нижней поверхностей имитатора. Например, в качестве имитатора типичного акселерометра массой 20 г может быть использован шестигранный стальной стержень длиной 12 мм с размером под ключ 16 мм. При частоте 5 кГц относительное движение между верхней и нижней поверхностями такого имитатора приводит к систематической погрешности около 0,26% при измерениях амплитуды и 4,2° при измерениях фазового сдвига.

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

4 Условия окружающей среды

Калибровку проводят при следующих условиях окружающей среды:

a) комнатная температура (23 +/- 3) °C;

b) относительная влажность не более 75%.

Должны быть приняты меры к тому, чтобы на результат измерений не оказывали влияние фоновые шум и вибрация.

5 Предпочтительные значения опорных ускорений и частот

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

a) ускорение (методы 1 и 3): 0,1; 0,2; 0,5; 1; 2; 5; 10; 50; 100; 200; 500; 1000 м/с2 (последнее значение — только для амплитуды);

b) частоты выбирают из стандартизованного ряда частот для третьоктавных полос (см. [1]) в диапазоне от 1 Гц до 10 кГц.

6 Требования, общие для всех методов

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

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

7 Метод счета полос (метод 1)

7.1 Общие положения

Метод используют для определения модуля коэффициента преобразования в диапазоне частот от 1 до 800 Гц.

Примечание — Метод 1 можно применять для калибровки в диапазоне частот выше установленного, если погрешность квантования устраняют специальными средствами (см. [2], [4]). Это позволяет проводить калибровку при заданной амплитуде ускорения (например, 100 м/с2) на более высоких частотах.

7.2 Испытания

После настройки интерферометра (см. 3.6) определяют коэффициент преобразования при заданных амплитудах ускорения и частотах вибрации (см. раздел 5) путем измерения частоты полос при помощи счетчика полос или отношения частоты вибрации к частоте полос при помощи измерения отношения частот (см. 3.7). Измерительная система для метода 1 показана на рисунке 1.

1 — генератор частоты (3.2); 2 — усилитель мощности (3.3); 3 — вибростенд (3.3); 4 — подвижная система вибростенда; 5 — имитатор; 6 — акселерометр; 7 — усилитель; 8 — интерферометр (3.6); 9 — лазер (3.5); 10 — фотодетектор; 11 — связь с измерителем отношения частот (3.7); 12 — счетчик или измеритель отношения частот (3.7); 13 — вольтметр (3.10); 14 — измеритель коэффициента гармоник (3.11); 15 — осциллограф (3.15)

Рисунок 1 — Измерительная система для метода счета полос (метод 1)

7.3 Представление результатов

Теоретические основы метода изложены в разделе B.1 (приложение B).

Рассчитывают амплитуду ускорения \hat{a} , м/с2, на входе акселерометра по показаниям счетчика полос по формуле

\hat{a} = ff_{f} \cdot 3,123 \cdot 10^{- 6}

и определяют модуль коэффициента преобразования S, В/м/с2, по формуле

S = \frac{\hat{u}}{ff_{f}} \cdot 0,3202 \cdot 10^{6} ,

где \hat{u} — амплитуда напряжения на выходе акселерометра, В;

f — частота вибрации, возбуждаемой вибростендом, Гц;

ff — частота следования интерференционных полос (то есть число полос, измеренное за достаточно продолжительный период времени, деленное на длительность периода времени), Гц.

Если используют измеритель отношения частоты, то амплитуду ускорения \hat{a} , м/с2, рассчитывают по формуле

\hat{a} = f^{2}R_{f} \cdot 3,123 \cdot 10^{- 6}

и определяют модуль коэффициента преобразования S, В/м/с2, по формуле

S = \frac{\hat{u}}{f^{2}R_{f}} \cdot 0,3202 \cdot 10^{6} ,

где Rf — отношение частоты следования полос ff к частоте вибрации f, измеренное на достаточном числе периодов вибрации (число периодов измерения зависит от частоты; например, при частоте 160 Гц период измерений должен составлять по меньшей мере 100 периодов вибрации).

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

8 Метод точек минимума (метод 2)

8.1 Общие положения

Метод используют для определения модуля коэффициента преобразования в диапазоне частот от 800 Гц до 10 кГц.

Примечание — Метод 2 основан на измерении перемещения с использованием значений аргументов функции Бесселя первого рода первого порядка в точках ее пересечения с осью абсцисс [см. раздел B.2 (приложение B)]. Тот же принцип может быть применен с использованием функции Бесселя первого рода нулевого порядка, однако в таком случае необходимо обеспечить модуляцию положения опорного зеркала (см. [5]).

8.2 Испытания

Сигнал с фотодетектора (см. 3.6) пропускают через полосовой фильтр (см. 3.8), среднегеометрическая частота которого совпадает с частотой возбуждения вибростенда. Отфильтрованный сигнал имеет точки минимума, соответствующие амплитудам перемещения, указанным в таблице 1.

Таблица 1

Амплитуды перемещения для точек минимума (\lambda = 0,63281\text{мкм})

Устанавливают частоту калибровки и регулируют амплитуду перемещения стола вибростенда от нуля до значения, при котором сигнал с фотодетектора после достижения максимального значения возвращается к минимальному. В данной точке минимума № 1 амплитуда перемещения равна 0,1930 мкм. Амплитуды перемещений для точек минимума других порядков приведены в таблице 1. Измерительная система для метода точек минимума показана на рисунке 2.

1 — генератор частоты (3.2); 2 — усилитель мощности (3.3); 3 — вибростенд (3.3); 4 — подвижная система вибростенда; 5 — имитатор; 6 — акселерометр; 7 — усилитель; 8 — интерферометр (3.6); 9 — лазер (3.5); 10 — фотодетектор; 11 — анализатор частоты (3.8); 12 — полосовой фильтр, настроенный на частоту вибрации (3.8); 13 — вольтметр; 14 — вольтметр (3.10); 15 — измеритель коэффициента гармоник (3.11); 16 — осциллограф (3.15)

Рисунок 2 — Измерительная система для метода точек минимума (метод 2)

Примечание 1 — Коэффициент преобразования акселерометра может быть также определен с использованием функции Бесселя первого рода нулевого порядка посредством модуляции положения опорного зеркала с частотой, которая мала по сравнению с частотой вибрации, и подстройкой среднегеометрической частоты полосового фильтра или частотного анализатора до значения, равного частоте модуляции зеркала (см. [6]).

Примечание 2 — Модуляция положения опорного зеркала может быть также применена для улучшения эффективности метода с использованием функции Бесселя первого рода первого порядка (см. [7]).

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

8.3 Представление результатов

Теоретические основы метода изложены в разделе B.2 (приложение B).

Вычисляют амплитуду ускорения акселерометра \hat{a} , м/с2, по формуле

\hat{a} = 39,478\hat{s}f^{2}

и определяют модуль коэффициента преобразования S, В/м/с2, по формуле

S = 0,02533 \frac{\hat{u}}{\hat{s}f^{2}} ,

где \hat{u} — амплитуда напряжения на выходе акселерометра, В;

\hat{s} — амплитуда перемещения для разных точек минимума в соответствии с таблицей 1;

f — частота вибрации, возбуждаемой вибростендом, Гц.

9 Метод синус-аппроксимации (метод 3)

9.1 Общие положения

Метод используют для определения модуля и/или фазового сдвига коэффициента преобразования в диапазоне частот от 1 Гц до 10 кГц.

9.2 Испытания

Измерительная система для метода синус-аппроксимации показана на рисунках 3 и 4.

Лазерный интерферометр настраивают так, чтобы его выходные квадратурные сигналы u₁ и u₂ были в пределах допусков, указанных в 3.6.

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

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

Примечание 1 — При амплитудах перемещения менее 0,5 мкм погрешность измерений, обусловленная изменениями квадратурных сигналов в пределах допусков, установленных в 3.6, не будет превышать 0,3% для модуля коэффициента преобразования и 0,3° для фазового сдвига. Уменьшить погрешность можно путем более тщательной регулировки для достижения значений меньших, чем установлено в 3.6 (см. [8]), или применением процедуры коррекции (см. [9]).

Примечание 2 — Чтобы измерить модуль и фазовый сдвиг комплексного коэффициента преобразования акселерометра при амплитудах перемещения в нанометровом диапазоне, метод синус-аппроксимации может быть применен с использованием соответствующей гетеродинной техники, как указано в [10] и [11]. Это позволяет проводить калибровку при средних уровнях амплитуд ускорения (например, 100 м/с2) на высоких частотах вибрации (например, 20 кГц).

1 — генератор частоты (3.2); 2 — усилитель мощности (3.3); 3 — вибростенд (3.3); 4 — подвижная система вибростенда; 5 — имитатор; 6 — акселерометр; 7 — усилитель; 8 — интерферометр (3.6); 9 — лазер (3.5); 10 — фотодетекторы; 11 — записывающее устройство (3.13); 12 — вольтметр (3.10); 13 — измеритель коэффициента гармоник (3.11); 14 — осциллограф (3.15)

Рисунок 3 — Измерительная система для метода синус-аппроксимации (метод 3)

1 — опорное зеркало; 2 — подвижная система вибростенда; 3 — акселерометр; 4 — имитатор; 5 — светоделитель; 6 — поляризатор; 7 — лазер; 8 — призма Волластона; 9 — фотодетекторы

Рисунок 4 — Лазерный интерферометр с квадратурным выходом

Примечание 3 — Чтобы повысить эффективность метода, может быть использована процедура умножения значений перемещения или фазовой модуляции на весовую функцию окна (см. [12]), если при этом будут соблюдены требования к неопределенности измерения, установленные в разделе 2.

9.3 Сбор данных

Частоты среза низкочастотных и высокочастотных (если используются) фильтров должны быть выбраны таким образом, чтобы искажающие влияния от низкочастотной и высокочастотной фильтрации на результаты калибровки находились в пределах допусков (см. [8]). Частота выборки должна быть не ниже удвоенной максимальной частоты в спектре сигнала в соответствии с теоремой Найквиста — Котельникова.

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

Выборку квадратурных сигналов интерферометра осуществляют на периоде измерений Tmeas. Интервал выборки \Delta t = t_{i} - t_{i - 1} поддерживают постоянным. В результате процедуры выборки получают последовательности значений квадратурных сигналов интерферометра {u1(ti)} и {u2(ti)}, а также выходного сигнала акселерометра {u(ti)}, которые заносят в память компьютера.

9.4 Обработка данных

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

9.4.2 По выборочным значениям квадратурных сигналов интерферометра {u₁(ti)} и {u₂(ti)} вычисляют последовательность значений сигнала фазовой модуляции {_{\phi_{mod}(t_{i})}} по формуле

\phi_{mod}(t_{i}) = arctg \frac{u_{2}(t_{i})}{u_{1}(t_{i})} + n\pi ,

где n = 0, 1, 2...

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

Примечание — Метод определения n описан в [13].

9.4.3 Выборочные значения {_{\phi_{mod}(t_{i})}} подгонкой параметров A, B и C методом наименьших квадратов аппроксимируют зависимостью (синус-аппроксимация)

\phi_{mod}(t_{i}) = Acos\omega t_{i} - Bsin\omega t_{i} + C ,

где i = 0, 1, 2, ..., N;

A = \phi_{M}cos\phi_{s} ;

B = \phi_{M}sin\phi_{s} ;

\omega — угловая частота вибрации, \omega = 2\pi f ;

\phi_{s} — начальный фазовый угол сигнала перемещения;

N + 1 — число синхронизированных выборочных значений на периоде измерений по 9.3.

Значение параметра C в дальнейших преобразованиях не используют.

9.4.4 По значениям параметров A и B, полученным в результате синус-аппроксимации, рассчитывают амплитуду \phi_{M} фазовой модуляции и начальный фазовый угол \phi_{s} по формулам

\phi_{M} = \sqrt{A^{2} + B^{2}} ,

\phi_{s} = arctg \frac{B}{A} .

9.4.5 Рассчитывают амплитуду ускорения \hat{a} и начальный фазовый угол ускорения \phi_{a} по формулам

\hat{a} = \pi\lambda f^{2}\phi_{M} ,

\phi_{a} = \phi_{s} + \pi .

9.4.6 Аналогично методом синус-аппроксимации выборочные значения выходного сигнала акселерометра {u(ti)} аппроксимируют зависимостью

u(t_{i}) = A_{u}cos\omega t_{i} - B_{u}sin\omega t_{i} + C_{u} ,

где A_{u} = \hat{u}cos\phi_{u} ;

B_{u} = \hat{u}sin\phi_{u} ;

\hat{u} — амплитуда выходного сигнала акселерометра;

\phi_{u} — начальный фазовый угол выходного сигнала акселерометра.

9.4.7 Рассчитывают амплитуду \hat{u} выходного сигнала акселерометра и начальный фазовый угол \phi_{u} по формулам

\hat{u} = \sqrt{A_{u}^{2} + B_{u}^{2}} ,

\phi_{u} = arctg \frac{B_{u}}{A_{u}} .

9.4.8 На основе результатов 9.4.5 и 9.4.7 рассчитывают модуль \hat{S}_{a} и фазовый сдвиг \Delta\phi комплексного коэффициента преобразования акселерометра по формулам

\hat{S}_{a} = \frac{\hat{u}}{\hat{a}} ,

\Delta\phi = (\phi_{u} - \phi_{a}) .

В протоколе калибровки результаты расчетов указывают с расширенной неопределенностью измерения, полученной в соответствии с приложением A (где для обозначения модуля коэффициента преобразования акселерометра для простоты вместо символа \hat{S}_{a} использовано обозначение S).

10 Протокол калибровки

В протоколе калибровки кроме указания использованного метода калибровки приводят следующую информацию:

a) условия окружающей среды:

• температуру акселерометра,
• температуру окружающего воздуха;

b) сведения об условиях крепления акселерометра:

• материал поверхности, на которую устанавливают акселерометр,
• момент затяжки акселерометра (если его устанавливают на шпильку),
• материал смазки, наносимой на монтажную поверхность (если смазка использовалась),
• способ закрепления кабеля,
• ориентацию акселерометра;

c) сведения об имитаторе (если использовался):

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

d) средства отражения луча лазера:

• используемый рефлектор (полированная поверхность или зеркало),
• положение лазерного пятна на отражающей поверхности;

e) настройки усилителя (если использовались), например:

• положение регулятора коэффициента усиления,
• значения частот среза фильтров;

f) результаты калибровки:

• амплитуды и частоты ускорения, на которых проводилась калибровка,
• значение коэффициента преобразования (модуля и, если измерялся, фазового сдвига),
• расширенная неопределенность измерения и используемое значение коэффициента охвата k (если было использовано значение k, отличное от двух).

Приложение A (обязательное)

СОСТАВЛЯЮЩИЕ НЕОПРЕДЕЛЕННОСТИ ИЗМЕРЕНИЯ ПРИ ПЕРВИЧНОЙ КАЛИБРОВКЕ АКСЕЛЕРОМЕТРОВ МЕТОДАМИ ЛАЗЕРНОЙ ИНТЕРФЕРОМЕТРИИ

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

A.1.1 Вычисление Urel(S) для метода 1

Относительную расширенную неопределенность модуля коэффициента преобразования Urel(S) для частот, амплитуд и настроек усилителя, при которых проводилась калибровка, вычисляют в соответствии с [2] по формулам

Urel(S) = kuc,rel(S),

u_{c,rel}(S) = \frac{u_{c}(S)}{S} = \frac{1}{S} \sqrt{\sum_{i} u_{i}^{2}(S)}

с коэффициентом охвата k = 2 (см. таблицу A.1).

Таблица A.1

A.1.2 Вычисление Urel(S) для метода 2

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

Таблица A.2

A.1.3 Вычисление Urel(S) и U(\Delta\phi) для метода 3

A.1.3.1 Вычисление Urel(S)

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

Таблица A.3

A.1.3.2 Вычисление U(\Delta\phi)

Расширенную неопределенность фазового сдвига коэффициента преобразования U(\Delta\phi) для частот, амплитуд и настроек усилителя, при которых проводилась калибровка, вычисляют в соответствии с [2] по формулам

U(\Delta\phi) = ku_{c}(\Delta\phi) ,

u_{c}(\Delta\phi) \sqrt[=]{\sum_{i} u_{i}^{2}(\Delta\phi)}

с коэффициентом охвата k = 2 (см. таблицу A.4).

Таблица A.4

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

A.2.1 Вычисление Urel(St)

Относительная расширенная неопределенность модуля коэффициента преобразования Urel(St), вычисленная согласно A.1.1, A.1.2 или A.1.3.1, справедлива только для частот, амплитуд и настроек усилителя, при которых была проведена калибровка. Относительную расширенную неопределенность модуля коэффициента преобразования Urel(St) для всего диапазона частот и амплитуд в любой момент времени между последовательными калибровками вычисляют по формулам

Urel(St) = kuc,rel(St),

u_{c,rel}(S_{t}) = \frac{u_{c}(S_{t})}{S} = \frac{1}{S} \sqrt{\sum_{i} u_{i}^{2}(S_{t})}

с коэффициентом охвата k = 2 (см. таблицу A.5).

Таблица A.5

A.2.2 Вычисление U(\Delta\phi_{t})

Расширенная неопределенность сдвига фаз коэффициента преобразования U(\Delta\phi) , вычисленная согласно A.1.3.2, справедлива только для частот, амплитуд и настроек усилителя, при которых была проведена калибровка. Расширенную неопределенность сдвига фаз коэффициента преобразования U(\Delta\phi_{t}) для всего диапазона частот и амплитуд в любой момент времени между последовательными калибровками вычисляют по формулам

U(\Delta\phi_{t}) = ku_{c}(\Delta\phi_{t}) ,

u_{c}(\Delta\phi_{t}) \sqrt[=]{\sum_{i} u_{i}^{2}(\Delta\phi_{t})}

с коэффициентом охвата k = 2 (см. таблицу A.6).

Таблица A.6

Приложение B (обязательное)

РАСЧЕТНЫЕ ФОРМУЛЫ

B.1 Метод 1

Длина волны \lambda основной линии в спектре излучения неона равна 0,63281 мкм при давлении 100 кПа.

Перемещение \Delta s на входе акселерометра, соответствующее расстоянию между двумя соседними полосами (между соседними максимумами или минимумами) интерферометра, определяют по длине волны \lambda по формуле

\Delta s = \frac{\lambda}{2} .

Таким образом, число максимумов за один период вибрации, и, соответственно, амплитуду перемещения \hat{s} можно определить по формулам

\frac{4\hat{s}}{\lambda/2} = \frac{f_{f}}{f} ,

\hat{s} = \frac{\lambda}{8} \cdot \frac{f_{f}}{f} ,

где f — частота вибрации, возбужденной вибростендом;

ff — частота следования полос в интерферометре.

На основе полученного значения \hat{s} амплитуду ускорения \hat{a} рассчитывают по формулам

\hat{a} = 4\pi^{2}f^{2}\hat{s} ,

\hat{a} = \frac{\pi^{2}\lambda}{2} f_{f}f .

B.2 Метод 2

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

\hat{s} = x_{n} \frac{\lambda}{4\pi} ,

\hat{a} = 4\pi^{2}f^{2}\hat{s} ,

где xn — аргументы функции Бесселя в точках ее пересечения с осью абсцисс, указанные в таблице B.1.

Таблица B.1

Значения аргументов xn функции Бесселя первого рода первого порядка в точках пересечения с осью абсцисс

B.3 Метод 3

В соответствии с [2] комплексный коэффициент преобразования преобразователя по ускорению Sa можно определить формулой

S_{a} = \hat{S}_{a}expj(\phi_{u} - \phi_{a}) , (B.1)

где \hat{S}_{a} — модуль коэффициента преобразования,

\hat{S}_{a} = \frac{\hat{u}}{\hat{a}} ; (B.2)

\hat{u} — амплитуда выходного сигнала акселерометра (или усилителя, если акселерометр калибруют совместно с усилителем);

\hat{a} — амплитуда ускорения стола вибростенда a;

\phi_{u} — начальный фазовый угол выходного сигнала;

\phi_{a} — начальный фазовый угол ускорения;

\Delta\phi — сдвиг фаз коэффициента преобразования,

\Delta\phi = \phi_{u} - \phi_{a} . (B.3)

Амплитуда \hat{a} и фаза \phi_{a} ускорения могут быть получены по значениям амплитуды перемещения \hat{s} , частоты f (угловой частоты \omega = 2\pi f ) и начальной фазы перемещения \phi_{s} по формулам

\hat{a} = \omega^{2}\hat{s} , (B.4)

\phi_{a} = \phi_{s} + \pi . (B.5)

Если вибростенд производит гармоническую вибрацию, то

s = \hat{s}cos(\omega t + \phi_{s}) , (B.6)

а выходной сигнал акселерометра может быть выражен формулой

u = \hat{u}cos(\omega t + \phi_{u}) . (B.7)

Выходной сигнал первого фотодетектора может быть представлен формулой

u_{1}(t) = u_{1}cos\phi_{mod}(t) = u_{1}cos[\phi_{0} + \phi_{M}cos(\omega t + \phi_{s})] , (B.8)

где \phi_{mod} — модулированная фаза сигнала,

\phi_{mod} = \phi_{0} + \phi_{M}cos(\omega t + \phi_{s}) , (B.9)

состоящая из начальной фазы сигнала фотодетектора \phi_{0} и модуляционной составляющей \phi_{M} , амплитуда \hat{\phi}_{M} которой пропорциональна перемещению

\phi_{M} = \hat{\phi}_{M}cos(\omega t + \phi_{s}) , (B.10)

\phi_{M} = \frac{4\pi\hat{s}}{\lambda} . (B.11)

Предполагают, что сдвиг фаз между перемещением s(t) и синусоидальной фазой \phi_{M}(t) отсутствует.

Выходной сигнал второго фотодетектора находится в квадратуре по отношению к первому и выражается формулой

u_{2}(t) = \hat{u}_{2}sin\phi_{mod}(t) = \hat{u}_{2}sin[\phi_{0} + \phi_{M}cos(\omega t + \phi_{s})] , (B.12)

причем \hat{u}_{2} = u_{1} .

Квадратурные сигналы подвергают выборке с постоянным значением периода выборки \Delta t = t_{i} - t_{i - 1} в течение периода измерения Tmeas, позволяющей получить серии выборочных значений {u1(ti)} и {u2(ti)}.

По обоим квадратурным сигналам последовательно в каждый момент времени в течение периода измерений вычисляют \phi_{mod}(t_{i}) по формуле

\phi_{mod}(t_{i}) = arctg \frac{u_{2}(t_{i})}{u_{1}(t_{i})} + n\pi , (B.13)

где n = 0, 1, 2...

Требуемые для расчета процедуры, в том числе вычисления функции арктангенса (см. [13]), являются стандартными в цифровой обработке сигналов. По значениям \phi_{mod}(t_{i}) могут быть получены параметры \hat{\phi}_{M} и \phi_{s} модуляционной составляющей фазы \phi_{M}(t) подгонкой методом наименьших квадратов (синус-аппроксимации) функции

\phi_{mod}(t_{i}) = Acos\omega t_{i} - Bsin\omega t_{i} + C,i = 0,1,2,\ldots,N , (B.14)

линейной относительно параметров

A = \hat{\phi}_{M}cos\hat{\phi}_{s},B = \phi_{M}sin\phi_{s} (B.15)

и постоянной C.

N + 1 — число синхронизированных выборочных значений в квадратурных сигналах на выходе фотодетекторов на периоде измерений.

По значениям параметров A и B рассчитывают амплитуду \hat{\phi}_{M} фазовой модуляции и начальный фазовый угол \phi_{s} по формулам

\hat{\phi}_{M} = \sqrt{A^{2} + B^{2}} , (B.16)

\phi_{s} = arctg \frac{B}{A} . (B.17)

Примечание — Для упрощения записи в формулах настоящего приложения не делается различия между оценкой величины [например, \hat{\phi}_{M} в формуле (B.16)] и ее «истинным» значением [например, \hat{\phi}_{M} в формуле (B.11)]. Там, где это необходимо для выражения неопределенности измерения, такое различие необходимо учитывать (см., например, [8]).

Для получения значения амплитуды ускорения \hat{a} формулу (B.11) следует подставить в формулу (B.4). Расчетная формула имеет вид

\hat{a} = \pi\lambda f^{2}\hat{\phi}_{M} . (B.18)

Фазу ускорения рассчитывают по формуле (B.5). Амплитуда \hat{u} и фаза \phi_{u} выходного сигнала акселерометра могут быть определены таким же путем, как соответствующие параметры модуляции фазы \hat{\phi}_{M} и \phi_{s} посредством синус-аппроксимации выборки значений выходного сигнала {u(ti)} с использованием формул, аналогичных (B.14) - (B.17). Подстановкой в формулы (B.2) и (B.3) полученных значений \hat{u} , \phi_{u} , \hat{a} и \phi_{a} можно определить модуль \hat{S}_{a} и фазовый сдвиг \Delta\phi коэффициента преобразования акселерометра.

Формулы, приведенные выше, в особенности (B.8) и (B.12), предполагают идеальные условия, которые невыполнимы на практике. В реальности сигналы фотодетекторов могут иметь несколько отличающиеся друг от друга амплитуды u_{1} и \hat{u}_{2} , а фазовый сдвиг между квадратурными составляющими может быть несколько отличен от \pi/2 . Кроме того, на практике могут наблюдаться сдвиги нулевого уровня сигналов фотодетекторов. Влияние неидеальности квадратурных сигналов в пределах допустимых значений, указанных в 3.6, кратко описано в 9.2 (примечание 1) (более подробно данный вопрос исследован в [8]).

БИБЛИОГРАФИЯ