СПОСОБ ИЗМЕРЕНИЯ ВИБРАЦИОННЫХ НАГРУЗОК НА ДВИГАТЕЛЕ ЛЕТАТЕЛЬНОГО АППАРАТА Российский патент 2016 года по МПК G01H17/00 

Описание патента на изобретение RU2580381C1

Область техники.

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

Уровень техники.

Известен способ выделения случайной компоненты вибраций машин периодического действия, заявленный в патенте SU 1280961 A1, с измерением координаты центра тяжести амплитудно-частотного спектра этой компоненты. Однако этот способ не предполагает получения характеристик гармонических и случайных компонент, требуемых для описания вибрации на двигателе ЛА, и применяется лишь для вибрационной диагностики машин периодического действия.

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

Известен способ измерения случайных вибронагрузок на двигателях силовой установки самолета, заключающийся в том, что регистрируют одновременно входные случайные вибрационные нагрузки в контрольных точках конструкций двигателей и выходную вибрационную нагрузку в контрольной точке упругой подвески силовой установки. Выделяют гармонические составляющие с использованием корреляционной функции. Эти составляющие определяются путем последовательной многократной фильтрации корреляционной функции смешанного процесса до момента выполнения краевых условий, выбранных по степени затухания корреляционной функции в пределах задаваемой при обработке временной реализации исходного вибрационного процесса (Пат. Российской Федерации 2036450 C1, опубл. 27.06.1995).

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

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

где R(m) - корреляционная функция;

m - дискретное преобразование времени;

k - дискретное преобразование частоты;

N - количество ординат обрабатываемой реализации вибропроцесса.

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

Этот способ имеет следующие недостатки:

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

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

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

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

Существенные признаки.

Для получения указанного технического результата в предлагаемом способе измерения вибрационных нагрузок на двигателе летательного аппарата, включающем регистрацию вибрационного процесса в контрольной точке двигателя, обработку вибрационных измерений, получение спектральной характеристики этого процесса, выделение гармонических составляющих, определение амплитудных уровней и частот гармонических составляющих, сопоставление их с аналогичными параметрами, установленными в технических требованиях, оценку вибрационной совместимости нового двигателя, работающего в составе силовой установки, зарегистрированный процесс обрабатывают как случайный с представлением результатов в графической форме. Вычисляют спектральную плотность мощности ускорения вибрационного процесса и вычисляют интегральную характеристику Qxx(f) как функцию частоты, представляющую накопленную по частоте площадь по кривой спектральной плотности мощности ускорения вибрационного процесса Wxx(f). Затем определяют амплитудные уровни и частоты гармонических вибрационных нагрузок путем выбора частот с резко выделяющимися на них уровнями спектральной плотности в Wxx(f) и выявления на этих же частотах fi скачкообразных приращений Δi ординат площади, которые характеризуют гармонические составляющие в анализируемом вибрационном процессе в интегральной характеристике Qxx(f). Определяют амплитуды ускорения гармонических составляющих вибраций зарегистрированного процесса по результатам приравнивания каждой Δi на частотах fi к среднему значению квадрата амплитуды гармонического колебания, равному величине дисперсии гармонического колебания на частоте fi, а все остальные частоты, где отсутствуют скачкообразные приращения, относят к случайным составляющим измеренного вибрационного процесса с уровнями полученной спектральной плотности мощности.

В реальных условиях полета на конструкции фюзеляжа ЛА и на двигателях возникают вибрационные процессы. Источники процессов затруднительно определить и измерить их характеристики. Для этого выполняют измерения параметров вибрационных процессов. Эти процессы могут иметь разные структуры, зависящие от наличия и комбинации в них случайных и гармонических (синусоидальных) составляющих. Для описания измеренного вибрационного процесса применяют характеристики отдельно для случайных составляющих и отдельно для гармонических составляющих. Зарегистрированный вибрационный процесс в виде временной реализации текущих ординат путем применения известных преобразований Фурье превращают в амплитудный спектр, если процесс гармонический, либо в спектр плотности мощности ускорения, если процесс случайный.

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

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

Для пояснения сущности изобретения на фиг. 1 представлен график спектральной плотности мощности амплитуд виброускорения (1) Wxx(fi), в единицах g2/Гц, и график интегральной характеристики (2) Qxx(f), представляющий изменение площади с ростом частоты fi под графиком (1), в единицах g2, а также показаны скачкообразные приращения ординат (3) интегральной характеристики на Qxx(fi) (2).

Способ осуществляется следующим образом.

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

Структура вибрационных процессов зависит от наличия гармонических (синусоидальных) и случайных составляющих и проявляется в частотной характеристике этих процессов. Как известно, гармонические вибрации характеризуются амплитудно-частотным спектром, отражающим распределение амплитуд вибрации xi на частотах fi. В качестве амплитуды вибрации могут быть амплитуды смещения, скорости, ускорения. Оценку характеристик смешанного вибрационного процесса, возникающего на двигателе самолета или вертолета в процессе полета, определяют путем обработки смешанного процесса, замеренного на двигателе, по алгоритмам характеристик случайного процесса получают спектральную плотность мощности ускорения вибрационного процесса Wxx(f) в виде графика, имеющего также смешанную структуру. Спектральная плотность мощности амплитуды гармонического колебания равна бесконечности на частотах этих колебаний и нулю при других значениях частоты, а интеграл спектральной плотности мощности, взятый в пределах диапазона частот (f÷fi) при (f-fi)→0, имеет конечное значение, равное среднему значению квадрата амплитуды ускорения на частоте fi, таким свойством обладают только гармонические вибрации. Аналитическое выражение гармонического слагаемого, действующего в составе смешанного процесса, по принятому интегральному преобразованию для спектральной плотности мощности ускорения определяется по формуле

где δ(f-fi) - дельта-функция;

- величина дисперсии гармонического колебания.

(См. литературу Г. Дженкинс, Д. Ваттс. «Спектральный анализ и его приложения». М.: «Мир», 1971).

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

где fн - нижняя частота исследуемого процесса;

fb - верхняя частота исследуемого процесса.

При этом в пределах диапазона частот (f-fi) интегральный спектр накапливает площадь в соответствии с формулой

и эта площадь равна .

Следовательно, интегральный спектр в точке f=fi получает скачкообразное приращение площади Δi, равное величине дисперсии гармонический колебаний , т.е.

где Ai - амплитуда гармонических колебаний на частоте fi.

На частотах гармонических составляющих значение функции Qxx(f) меняется скачком на величину , где Ai - амплитуда гармонической составляющей на частоте fi. Вычисленные значения спектральной плотности мощности ускорения вибрационного смешанного процесса Wxx(fi) на этих частотах fi обнуляются, так как на этих частотах присутствуют гармонические составляющие вибрации с частотой колебания fi и амплитудой Am i. На всех остальных частотах исследуемый процесс носит случайный характер.

Точность определения значения частоты по предлагаемому способу определяется дискретным шагом по частоте Δf при вычислении Wxx(f), рекомендуется значение Δf выбирать в пределах 0,5 Гц…1,0 Гц. Значение Ai не зависит от этого параметра.

Пример.

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

Зарегистрированный процесс обрабатывают как случайный с представлением результатов в графической форме, см. фиг. 1 (1, 2). Определяют в характеристике Wxx(fi) (1) на частотах f1, f2, f3 и f4 резко выделяющиеся уровни спектральной плотности мощности. В интегральной характеристике Qxx(f) (2) выявлены скачкообразные приращения Δi (3) ординат площади на тех же частотах f1, f2, f3 и f4. Указанные приращения характеризуют гармонические составляющие в анализируемом вибрационном процессе. Затем по формуле (1) определяют амплитуды ускорения гармонических составляющих вибрации зарегистрированного процесса:

на частоте f1 приращение Δ1=2,2 g2, что соответствует A1=2,1 g;

на частоте f2 приращение Δ2=5,1 g2, что соответствует A2=3,17 g;

на частоте f3 приращение Δ3=1,03 g2, что соответствует A3=1,43 g;

на частоте f4 приращение Δ4=2,7 g2, что соответствует A4=2,32 g.

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

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

Похожие патенты RU2580381C1

название год авторы номер документа
Способ выявления синусоидальной и случайной вибраций в составе смешанного вибрационного процесса, измеряемого на летательном аппарате 2021
  • Митенков Виктор Борисович
  • Парусова Марина Георгиевна
  • Баранова Марина Сергеевна
  • Кудашин Владимир Сергеевич
  • Саркисян Анаида Фрунзевна
  • Федулов Андрей Владимирович
  • Харатян Вова Гамлетович
RU2782683C1
СПОСОБ ИСПЫТАНИЙ БОРТОВОЙ АППАРАТУРЫ КОСМИЧЕСКОГО АППАРАТА НА ВИБРАЦИОННЫЕ ВОЗДЕЙСТВИЯ 2011
  • Орлов Александр Сергеевич
  • Орлов Сергей Александрович
RU2476845C2
СПОСОБ ИСПЫТАНИЙ ОБОРУДОВАНИЯ НА МЕХАНИЧЕСКИЕ ВОЗДЕЙСТВИЯ 2009
  • Усанов Алексей Юрьевич
  • Орлов Сергей Александрович
  • Орлов Александр Сергеевич
RU2399032C1
Способ вибродиагностики технического состояния газотурбинных двигателей на ресурсосберегающих режимах с применением теории инвариантов 2020
  • Шигапов Ильяс Ильгизович
  • Попов Николай Николаевич
  • Казаринов Александр Николаевич
  • Сенной Николай Николаевич
  • Соколов Антон Григорьевич
  • Голубев Константин Геннадьевич
RU2754479C1
СПОСОБ ИЗМЕРЕНИЯ СЛУЧАЙНЫХ ВИБРОНАГРУЗОК НА ДВИГАТЕЛЯХ СИЛОВОЙ УСТАНОВКИ САМОЛЕТА 1992
  • Конычев В.И.
  • Митенков В.Б.
  • Важнов Ю.Л.
RU2036450C1
Способ вибродиагностики технического состояния газотурбинных двигателей на ресурсосберегающих режимах с применением теории инвариантов 2020
  • Шигапов Ильяс Ильгизович
  • Попов Николай Николаевич
  • Казаринов Александр Николаевич
  • Сенной Николай Николаевич
  • Соколов Антон Григорьевич
  • Голубев Константин Геннадьевич
RU2754476C1
СПОСОБ ДИНАМИЧЕСКИХ ИСПЫТАНИЙ КОНСТРУКЦИЙ И СИСТЕМ НА МЕХАНИЧЕСКИЕ И ЭЛЕКТРОННЫЕ ВОЗДЕЙСТВИЯ 2014
  • Зельманов Самуил Соломонович
  • Крылов Владимир Владимирович
RU2569636C2
УСТРОЙСТВО ДЛЯ ИСПЫТАНИЙ НА СЛУЧАЙНУЮ ВИБРАЦИЮ 1991
  • Баширов З.А.
  • Валеев А.А.
  • Гавриленко О.Н.
RU2008628C1
СПОСОБ ВИБРАЦИОННЫХ ИСПЫТАНИЙ ОБЪЕКТОВ 1990
  • Кузьмин Э.Н.
  • Захарова Н.Ф.
  • Синякина Л.П.
RU1773164C
Способ контроля технического состояния подшипников качения 2015
  • Лукин Виктор Федорович
  • Сенной Николай Николаевич
  • Спирькин Андрей Владимирович
  • Селезнев Алексей Андреевич
RU2623177C2

Иллюстрации к изобретению RU 2 580 381 C1

Реферат патента 2016 года СПОСОБ ИЗМЕРЕНИЯ ВИБРАЦИОННЫХ НАГРУЗОК НА ДВИГАТЕЛЕ ЛЕТАТЕЛЬНОГО АППАРАТА

Изобретение относится к метрологии, в частности к способам вибрационной диагностики. Способ измерения вибрационных нагрузок на двигателе летательного аппарата предполагает измерение вибраций с последующим определением спектральной плотности мощности ускорения вибрационного процесса Wxx(f), рассчетом интегральной характеристики Qxx(f), представляющей накопленную по частоте площадь под кривой спектральной плотности мощности. Затем определяют амплитудные уровни и частоты вибрационных нагрузок путем выбора частот с пиковыми уровнями спектральной плотности в Wxx(f) и выявления скачкообразных приращений Δi ординат площади. Определяют амплитуды ускорения по результатам приравнивания каждой Δi на частотах fi к среднему значению квадрата амплитуды гармонического колебания. Все остальные частоты, где отсутствуют скачкообразные приращения, относят к случайным составляющим измеренного вибрационного процесса с уровнями полученной спектральной плотности мощности. Сопоставляют их с аналогичными параметрами, установленными в технических требованиях, и оценивают вибрационную совместимость нового двигателя, работающего в составе силовой установки. Технический результат - повышение точности измерений, сокращение времени на диагностику. 1 ил.

Формула изобретения RU 2 580 381 C1

Способ измерения вибрационных нагрузок на двигателе летательного аппарата, включающий регистрацию вибрационного процесса в контрольной точке двигателя, обработку вибрационных измерений, получение спектральной характеристики этого процесса, выделение гармонических составляющих, определение амплитудных уровней и частот гармонических составляющих вибрационных нагрузок, сопоставление их с аналогичными параметрами, установленными в технических требованиях, оценку вибрационной совместимости нового двигателя, работающего в составе силовой установки, отличающийся тем, что зарегистрированный процесс обрабатывают как чисто случайный с представлением результатов в графической форме, вычисляют спектральную плотность мощности ускорения вибрационного процесса Wxx(f) и по полученным ее значениям вычисляют интегральную характеристику Qxx(f) как функцию частоты, представляющую накопленную по частоте площадь по кривой спектральной плотности мощности ускорения вибрационного процесса Wxx(f), затем определяют амплитудные уровни и частоты гармонических вибрационных нагрузок путем выбора частот с резко выделяющимися на них уровнями спектральной плотности в Wxx(f) и выявления на этих же частотах fi скачкообразных приращений Δi ординат площади, которые характеризуют гармонические составляющие в анализируемом вибрационном процессе интегральной характеристики Qxx(f), определяют амплитуды ускорения гармонических составляющих вибраций зарегистрированного процесса по результатам приравнивания каждой Δi на частотах fi к среднему значению квадрата амплитуды гармонического колебания, равному величине дисперсии гармонического колебания на частоте fi, а все остальные частоты, где отсутствуют скачкообразные приращения, относят к случайным составляющим измеренного вибрационного процесса с уровнями полученной спектральной плотности мощности.

Документы, цитированные в отчете о поиске Патент 2016 года RU2580381C1

СПОСОБ ИЗМЕРЕНИЯ СЛУЧАЙНЫХ ВИБРОНАГРУЗОК НА ДВИГАТЕЛЯХ СИЛОВОЙ УСТАНОВКИ САМОЛЕТА 1992
  • Конычев В.И.
  • Митенков В.Б.
  • Важнов Ю.Л.
RU2036450C1
СПОСОБ ИСПЫТАНИЙ КОСМИЧЕСКОГО АППАРАТА НА ВИБРОАКУСТИЧЕСКИЕ ВОЗДЕЙСТВИЯ 2007
  • Орлов Александр Сергеевич
  • Орлов Сергей Александрович
RU2354948C1
СПОСОБ ОПРЕДЕЛЕНИЯ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК ТЕНЗОМЕТРИЧЕСКОГО ПРЕОБРАЗОВАТЕЛЯ ДАВЛЕНИЯ (ВАРИАНТЫ) 2011
  • Семенов Александр Сергеевич
  • Бушуев Олег Юрьевич
RU2466368C1
СПОСОБ ВИБРОДИАГНОСТИКИ МЕХАНИЗМОВ ПО ХАРАКТЕРИСТИЧЕСКОЙ ФУНКЦИИ ВИБРАЦИИ 2012
  • Костюков Владимир Николаевич
  • Науменко Александр Петрович
  • Бойченко Сергей Николаевич
  • Кудрявцева Ирина Сергеевна
RU2517772C1
СПОСОБ ПРОИЗВОДСТВА ВАФЕЛЬ (ВАРИАНТЫ) 2010
  • Квасенков Олег Иванович
RU2427213C1
WO 2001033182 A1, 10.05.2001
JP 2004069480 A, 04.03.2004
US 6363789 B1, 02.04.2002..

RU 2 580 381 C1

Авторы

Митенков Виктор Борисович

Саркисян Анаида Фрунзевна

Баранова Марина Сергеевна

Власова Тамара Аркадьевна

Даты

2016-04-10Публикация

2015-03-11Подача