RU2514061C2 - Способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток - Google Patents

Способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток Download PDF

Info

Publication number
RU2514061C2
RU2514061C2 RU2009123249/28A RU2009123249A RU2514061C2 RU 2514061 C2 RU2514061 C2 RU 2514061C2 RU 2009123249/28 A RU2009123249/28 A RU 2009123249/28A RU 2009123249 A RU2009123249 A RU 2009123249A RU 2514061 C2 RU2514061 C2 RU 2514061C2
Authority
RU
Russia
Prior art keywords
blades
timing
probe
events
resonant
Prior art date
Application number
RU2009123249/28A
Other languages
English (en)
Other versions
RU2009123249A (ru
Inventor
Питер РУШАРД
Джейсон Дэвид БЭК
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 РОЛЛС-РОЙС Пи-Эл-Си, Великобритания
Publication of RU2009123249A publication Critical patent/RU2009123249A/ru
Application granted granted Critical
Publication of RU2514061C2 publication Critical patent/RU2514061C2/ru

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H1/00Measuring characteristics of vibrations in solids by using direct conduction to the detector
    • G01H1/003Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
    • G01H1/006Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines of the rotor of turbo machines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Lubrication Of Internal Combustion Engines (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

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

Description

Настоящее изобретение относится к временному анализу вращающихся лопаток (лопастей), таких как в газотурбинном двигателе.
При разработке газотурбинных двигателей важно определить величину вибрации вращающихся лопаток. По измерениям вибрации можно определить напряжение, вызываемое в лопатках. Затем можно предпринять действие для предупреждения достаточно высоких напряжений, повреждающих лопатки.
Как описано, например, в патентной заявке №2002/0162395, известен тензодатчик, установленный на вращающиеся лопатки компрессора/турбины для передачи информации об амплитудах и частотах вибрации лопаток. Может иметь место один или несколько тензодатчиков на каждой лопатке и соединяться с системой радио телеметрии, установленной на роторе, которая передает результаты измерений от ротора. Однако из-за количества тензодатчиков, требуемых для полного определения вибраций, система телеметрии обычно сложна, дорога, велика и занимает много времени для установки на роторе.
Альтернативный метод для характеристики вращения лопаток это «временной законцовки лопатки» (ВТТ), в котором бесконтактные временные зонды, обычно установленные на кожухе двигателя, используются для измерения времени, за которое каждая лопатка проходит через каждый зонд. Это время сравнивается со временем, за которое лопатка прошла бы зонд, если бы не подвергалась вибрации. Это названо «ожидаемым временем прибытия», и его можно рассчитать из положения вращения конкретной лопатки на роторе в сочетании с «один раз за оборот» или «OPR» сигнал, который дает сведения о положении ротора. Сигнал OPR получается от времени, в которое индикатор на роторе проходит через опорный датчик, и его использование хорошо известно из уровня техники.
Различие между ожидаемым временем прибытия и реальным временем прибытия можно умножить на скорость законцовки лопатки для сообщения смещения лопатки относительно его ожидаемого положения. Таким образом, данные ВТТ от конкретного зонда эффективно измеряют смещение законцовки лопатки у зонда.
Благоприятно, что способ временного анализа законцовки не требует наличия системы телеметрии, поскольку зонды установлены на кожухе.
Однако из-за того, что скорость выборки зондов определяется частотой вращения ротора, она часто ниже минимально допустимой частоты дискретизации для данных вибраций. Таким образом каждый зонд берет пониженную выборку вибраций, что приводит к проблемам искажения информации. Еще одна проблема с данными ВТТ состоит в том, что они часто шумовые по сути из-за движения зондов, вызванных ограничениями по установке и толщиной кожуха.
Тем не менее при множестве зондов отсчета времени можно, в принципе, осуществлять полезный анализ вибрации, который можно преобразовать в напряжение в лопатке. Однако традиционно такой анализ требует вмешательства рук человека и в высшей степени зависит от мастерства оператора, проводящего анализ. Действительно, разные операторы, сталкивающиеся с одними и теми же данными, могут анализировать их по-разному с выдачей разных результатов.
Во многом из-за необходимости в умелом вмешательстве руки человека использование информации о тайминге законцовки обычно ограничивается фазой разработки газотурбинных двигателей.
Настоящее изобретение направлено на решение проблем субъективности оператора и воспроизводимости анализа посредством создания методов, которые можно использовать для автоматизации и осуществления анализа данных ВТТ. Преимущество автоматизированной системы анализа заключается в том, что ее можно использовать для контроля рабочего состояния оборудования в реальном времени, т.е. не только для фазы разработки двигателей.
Первый аспект данного изобретения обеспечивает способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток, установленных на роторе, причем ряд расположенных по периметру стационарных зондов тайминга, связанных с лопатками, определяет моменты, в которые лопатки проходят соответствующие зонды, причем способ включает следующие этапы.
- получение тайминга лопаток, обнаруженного зондами,
- определение для последовательных вращений узла соответствующих факторов корреляции для одной или более лопаток, причем каждый фактор корреляции определяет величину степени корреляции между таймированием лопаток, обнаруженным зондами для конкретной лопатки на конкретном вращении, и таймированием лопаток, обнаруженным зондами на предыдущем вращении; и
- определение события резонансной вибрации, когда один или более факторов корреляции пересекает установленный порог. Этап получения может включать обнаружение моментов, в которые лопатки проходят соответствующие зонды.
Предпочтительно, использование фактора корреляции может помочь избежать субъективности оператора, давая повышенную воспроизводимость и постоянство определение событий резонансных вибраций. В особенности в отношении событий синхронной вибрации идентификация, основанная на факторе корреляции, может привести к идентификации событий, которые можно просмотреть ручным анализом.
Факторы корреляции можно определить для ряда лопаток и, например, объединить в средний фактор корреляции. Случай резонансной вибрации можно затем идентифицировать, когда средний фактор корреляции пересекает установленный порог.
Способ может дополнительно включать этап обработки таймирования лопаток путем выполнения подэтапов:
выбор одного из зондов таймирования;
подбор усредняющей кривой к таймированию лопаток выбранного зонда вне идентифицированных событий резонансной вибрации;
интерполирование секций для усредняющей кривой у таймирования лопатки, выбранного зонда внутри определенных событий резонансных вибраций,
повышение усредняющей кривой с интерполированными секциями; и
вычитание повышенной усредняющей кривой из таймирования лопаток выбранного зонда для получения обнуленного таймирования лопаток для выбранного зонда.
Эта процедура позволяет удалить сдвиги постоянного тока в данных таймирования (вызванных, например, аэродинамической нагрузкой на лопатки) до осуществления последующих этапов анализа.
Способ может дополнительно включать этап фильтрации таймирования лопаток путем выполнения следующих подэтапов:
- выбор одного из зондов таймирования,
- преобразование таймирования лопатки выбранного зонда в частотную область, причем преобразованное таймирование выдает отслеженные разряды в случаях резонансных частот в графике частоты относительно времени,
- определение интегральных положений выборки как частот, соответствующих частоте вращения ротора и их кратных до кратного, соответствующего количеству лопаток в узле,
- для определенных событий резонансной вибрации, определение отслеженных порядковых позиций относительно интегральных положений выборки и
- фильтрация таймирования лопаток выбранного зонда определенных событий резонансной вибрации во временной области, характеристика фильтра, меняющегося в зависимости от соответствующих отслеженных относительных (порядковых) положений (позиций) разрядов в частотной области.
Благоприятно, что соответствующие характеристики фильтра могут затем быть постоянно и воспроизводимо применены к разным типам событий резонансной вибрации. Например, событие асинхронной вибрации посреди интегральных положений выборки обычно требует разных характеристик фильтра по событию асинхронной вибрации при интегральном положении выборки. Таким образом, фильтр может приспосабливаться к спектральному содержанию данных таймирования, что определяется анализируемой частотой сигнала.
Второй аспект изобретения соответствует этапу факультативной обработки 1-го аспекта и обеспечивает способ обрабатывающего таймирования лопаток, обнаруженного стационарным зондом таймирования, связанных с узлом вращающихся лопаток, установленных на роторе, причем зонд обнаруживает моменты, в которые лопатки проходят зонд, при этом способ включает следующие этапы:
- получение таймирования лопаток, обнаруженного зондом,
- идентификация одного или нескольких событий резонансной вибрации в данных таймирования;
- подбор усредняющей кривой к таймированию лопаток вне идентифицированных событий резонансной вибрации,
- интерполирование секций для усредняющей кривой у таймирования лопаток внутри идентифицированных событий резонансных вибраций;
- увеличение усредняющей кривой с интерполированными секциями, и
- вычитание увеличенной усредняющей кривой из таймирования лопаток для получения обнуленного таймирования лопаток.
Этап получения может включать обнаружение моментов, в которые лопатки проходят зонд.
Способ может дополнительно включать этап фильтрации таймирования лопаток путем выполнения следующих подэтапов:
- преобразование таймирования лопатки в частотную область, причем преобразованное таймирование выдает отслеженные разряды в событиях резонансных частот в графике частоты относительно времени,
- определение интегральных положений выборки как частот, соответствующих частоте вращения ротора и их кратных до кратного, соответствующего количеству лопаток в узле,
- для идентифицированных событий резонансной вибрации определение отслеженных положений разряда относительно интегральных положений выборки и
- фильтрация таймирования лопаток при идентифицированных событиях резонансной вибрации, причем характеристика фильтра меняется в зависимости от соответствующих отслеженных положений разряда в частотной области.
Третий аспект изобретения соответствует факультативным этапам фильтрации 1-го и 2-го аспекта и обеспечивает способ фильтрации таймирования лопаток, обнаруженного стационарным зондом таймирования, связанного с узлом вращающихся лопаток, установленных на роторе, причем зонд обнаруживает моменты, когда лопатки проходят зонд, причем способ включает следующие этапы:
- получение таймирования лопаток, обнаруженного зондом,
- идентификация одного или нескольких событий резонансной вибрации в данных таймирования;
- преобразование таймирования лопатки в частотную область, причем преобразованное таймирование выдает отслеженные разряды в событиях резонансных частот в графике частоты относительно времени,
- определение интегральных положений выборки как частот, соответствующих частоте вращения ротора и их кратных до кратного, соответствующего количеству лопаток в узле,
- для идентифицированных событий резонансной вибрации определение отслеженных положений разряда относительно интегральных положений выборки и
- фильтрация таймирования лопаток для идентифицированных событий резонансных вибраций в области времени, причем изменяющаяся характеристика фильтра зависит от соответствующих отслеженных положений разрядов в частотной области.
Этап получения может включать обнаружение моментов прохождения лопаток через зонды.
Следующие факультативные признаки относятся к 1-му, 2-му и 3-му аспекту. Это или каждое событие идентифицированной резонансной вибрации может быть синхронной резонансной вибрацией. Предпочтительно, этот или каждый зонд измеряет отклонения законцовки лопаток. Частоты событий вибрации могут быть недостаточно выбраны этим или каждым зондом.
После выполнения способа по 1-му, 2-му или 3-му аспекту можно извлечь амплитуды или/и частоты вибрации лопатки из таймирования лопаток с использованием методов, известных специалисту.
Четвертый аспект изобретения обеспечивает компьютерную систему для выполнения способа любого из 1-го, 2-го и 3-го аспектов. Пятый аспект изобретения предусматривает наличие компьютерной программы для выполнения способа любого из 1-го, 2-го и 3-го аспектов. Шестой аспект изобретения предусматривает компьютерный программный продукт, включающий программу согласно 5-му аспекту.
Варианты осуществления изобретения теперь будут описаны посредством примеров со ссылкой на сопроводительные чертежи, на которых:
Фиг.1 схематически иллюстрирует приспособление ВТТ;
Фиг.2 - это блок-схема, показывающая этапы процедуры при обработке данных таймирования, полученных зондами фиг.1;
Фиг. 3 показывает обычные данные смещения от зонда таймирования для одной лопатки;
Фиг.4 показывает еще одни данные о смещении для разных лопаток для одного зонда;
Фиг.5 показывает обычную нерезонансную картинку стека для зонда таймирования,
Фиг.6 показывает картинку стека фиг.5, на которую наложены нерезонансные картинки стека для ряда других зондов таймирования;
Фиг.7 показывает картинку стека для зонда таймирования фиг. 5 во время случая вибрации.
Фиг.8 показывает график корреляции смещения лопатки относительно нерезонансных смещений лопатки по времени для одного зонда. Фиг.9а показывает данные о смещении от набора зондов таймирования для одной лопатки, и фиг.9b-d соответственно показывают соответственно соответствующие картинки стека в 3-х случаях во время периода времени данных о смещении;
Фиг.10 схематически показывает детектор для идентификации события;
Фиг.11 показывает данные смещения для пары зонд/лопатка и 4 области, идентифицированные как события резонансной вибрации;
Фиг.12 показывает кривую 21, полученную усреднением данных фиг.11, и еще одну кривую 22, которая является результатом вычитания усредняющей кривой 21 из данных фиг.11 вне 4-х областей, которые идентифицированы как события резонансной вибрации;
Фиг.13 показывает обнуленные данные фиг. 11;
Фиг.14 - это блок-схема ядра свертки для получения усредняющей кривой.
Фиг.15 схематически показывает данные смещения 30 при случае вибрации и соответствующие данные частотной области 31.
Фиг.16 - это схематический график частоты относительно времени, который показывает отслеженные разряды как толстую непрерывную наклонную линию и интегральные положения выборки как горизонтальные пунктирные линии.
Фиг.17а схематически показывает теоретические данные смещения во временной области для вибраций, зарегистрированных для одного зонда, и
фиг.17b схематически показывает представление данных реального смещения во временной области, на которые повлияли шумы.
Фиг.18а показывает эффект контрольного закона, налагающего линейное изменение в частоте среза с положением бункера, приложенным к данным;
фиг.17b и фиг.18b показывает эффект контрольного закона, налагающего экспоненциальное изменение частоты среза с положением бункера, наложенного на данные фиг 17b; и
Фиг.19а показывает график частота-время, отслеженный разряд, соответствующий вибрации, показанный как толстой линией с последовательными секциями А, В, С, и
фиг.19b-d - это соответствующие графики смещения от времени, показывающие соответственно эффекты разных фильтров А-С на данные смещения, собранные одним зондом от всех лопаток ротора за период времени отслеженного разряда.
Фиг.1 схематически иллюстрирует приспособление ВТТ. Зонд OPR 1 отслеживает положение ротора 2, тогда как от 1 до n зондов ВТТ 3 дают таймирование для лопаток 4, установленных на роторе.
Фиг.2 - это блок-схема, иллюстрирующая этапы процедуры при обработке данных таймирования, полученных зондами. Прежде всего данные анализируются для идентификации событий резонансных вибраций. После определения идентификации этих событий эти данные обнуляются и преобразуются в частотную область. Преобразованные данные фильтруются для уменьшения количества шумов в сигналах. Данные затем готовы для дальнейшего анализа, такого как снятие характеристик (например, определение количества фаз и амплитуд) событий вибрации способом, известным для специалиста. Ниже мы обсуждаем подробнее этапы идентификации событий резонансной вибрации, обнуление данных и фильтрация данных.
Фиг.3 иллюстрирует обычное смещение данных от зонда таймирования для одной лопатки. Данные прошли предварительную обработку для отклонения всплесков, но однако содержат то же количество шумов, как начальные данные. Скорость ротора, увеличенная в течение сбора данных, и большое снижение постоянного тока в данных вызваны последующей аэродинамической нагрузкой на лопатки.
Данные являются относительно искаженными, но при установлении простого порога амплитуды, по меньшей мере, одно возможное событие вибрации можно идентифицировать как примерно 800 оборотов.
Однако установка порогового значения амплитуды может стать неверным подходом к идентификации вибрации. Во-первых, установка порога может вызвать ложные положительные сигналы из-за шумов в сигнале. Во-вторых, установка порога амплитуды может пропустить особенно события вибрации, вызванные синхронной реакцией (ответом) лопаток. В синхронном ответе лопатка вибрирует с частотой, которая является точным кратным частоты выборки для лопатки зондов ВТТ (т.е. точным кратным частоты OPR). Обычно такой ответ создает только сдвиг постоянного тока в данных таймирования, и этот сдвиг может быть относительно небольшим. В отличие от этого, асинхронные отзывы, для которых не существует корреляции между частотой ответа на вибрацию лопатки и частоты выборки зонда ВТТ обычно легче идентифицировать путем простой установки порога амплитуды, хотя шумы сигнала могут еще вызвать трудности, и сдвиг постоянного тока в положении лопатки требует, чтобы пороговый уровень менялся на протяжении всего маневра.
Фиг.4 иллюстрирует дополнительные обычные (типичные) данные о смещении для разных лопаток у одного и того же зонда. События вибрации отмечены стрелками. Длительность каждого события может быть очень малой, и каждая лопатка может реагировать (отвечать) по-разному в разное время, что делает чрезвычайно субъективным визуальную идентификацию событий.
Более надежный способ идентификации синхронных и асинхронных случаев вибрации основан на понимании того, что данные таймирования от каждого зонда дают характеристику «рисунок стека», и что этот рисунок нарушается для синхронных и асинхронных случаев.
Фиг.5 иллюстрирует типичный нерезонансный рисунок стека для одного зонда в отношении ротора с 20 лопатками, причем число лопаток откладывается по оси у, а смещение каждой лопатки - по оси х. Из-за механических различий в расстоянии между лопатками, вызванных допусками при изготовлении, смещения немного меняются для каждой лопатки, давая уникальный рисунок характеристики для ротора. Каждая точка рисунка стека получается усреднением нерезонансного таймирования лопаток для соответствующей лопатки
До высокой степени корреляции другие зонды дают тот же нерезонансный рисунок, как показано на фиг.6, которая иллюстрирует рисунок стека фиг.5, на который наложены нерезонансные рисунки для ряда других зондов. Каждый зонд наблюдает тот же рисунок стека, когда данные скоррелированы для периферического положения зонда
Однако по мере изменения условия вращения ротора и начала вибрации рисунки стека начинают расширяться (развертываться), как показано на фиг.7, вызывая падение степени корреляции. Аналогичный рисунок стека с расширением наблюдается для синхронных и асинхронных событий.
Значение смещения на каждом зонде для одной лопатки коррелируется для значений для каждого следующего оборота для получения значений между 0 и 1 с указанием коэффициента корреляции результата момента Пирсона для 2-х наборов данных. Это значение возводится в квадрат для получения фактора корреляции R, а именно
Figure 00000001
где n - это количество зондов, Di и di - это смещения выбранных лопаток у 1-го зонда для, соответственно, данного оборота (вращения) и предыдущего оборота, D и d - это средние смещения выбранной лопатки у зондов для соответственно данного оборота и предыдущего оборота, и SD и Sd - это стандартное отклонение смещений лопаток для соответственно данного оборота и предыдущего оборота.
Порог, например, менее 0,95 для фактора корреляции используется для определения наличия вибрации.
Обычно все лопатки одновременно реагируют на событие вибрации, так что в основном фактор нужно только рассчитать для одной лопатки. Однако факторы корреляции всех или нескольких лопаток можно сочетать, например, как средний фактор корреляции.
Фиг.8 иллюстрирует график корреляции смещения лопаток относительно нерезонансного смещения лопаток от времени для одного зонда. Для сравнения, соответствующие данные смещения также имеются для одной лопатки у зонда. Основное смещение происходит одновременно с большим падением корреляции, что указывает на возможную вибрацию.
Фиг.9а иллюстрирует данные смещения от зонда таймирования для одной лопатки и фиг.9b-d - соответствующие рисункам стека в трех случаях во время периода данных смещения. Левый рисунок стека на фиг.9b - это расстояние, когда не происходит вибрации. Как ожидается, величины смещения для каждой лопатки хорошо скоррелированы. Средний и правый рисунок стека на фиг.9с и d, соответственно, однако, являются копиями, когда нет синхронной вибрации. Эти случаи объясняются потерей корреляции в рисунках стека. Для среднего рисунка стека на фиг.9с соответствующая часть данных смещения показывает пик, который возможно был бы обнаружен простым пороговым определением данных смещения. С другой стороны, соответствующая часть данных о смещении для правого рисунка стека на фиг.9d не имеет достаточного пика и потому могла бы быть потеряна при использовании простого определения порога или подходов с визуальной идентификацией.
Фиг.10 схематически отображает детектор для идентификации события. Дискретные данные 10 от ряда зондов посылаются в устройство 11 задержки выравнивания, в котором данные от каждого зонда выравниваются в соответствующий массив, имеющий ряд ячеек памяти. Данные поступают не выровненными вовремя из-за отставания зондов относительно ротора. Однако в устройстве 11 задержки выравнивания данные регулируются (смещаются) в выравнивание путем подсчета количества лопаток между каждым зондом и наложения соответствующего фиксированного смещения, состоящего из дополнительных ячеек памяти для каждого массива. Данные можно также усреднить по более чем одному обороту.
Данные перемещаются через массивы выравнивания со скоростью, равной частоте зонда лопатки (частота OPR × количество лопаток). Количество полных оборотов, удержанных в массивах, можно изменить для увеличения или уменьшения времени отклика детектора событий.
Выровненные данные затем перемещаются к дальнейшим массивам в памяти 12, причем каждый дальнейший массив соответствует соответствующему зонду и имеет одну ячейку для каждой лопатки.
Наконец, на данные в дальнейших массивах в памяти 12 действует DSP (цифровой процессор сигналов) устройство корреляции 13 для получения результата корреляции 14. В зависимости от значения этого результата можно сигнализировать о событии вибрации.
Затем как предшественник наложения быстрого преобразования Фурье (FFT) данные о смещении для каждой пары зонд/лопатка обнуляются путем удаления DC направления постоянного тока (такого как большое направление постоянного тока, показанное на фиг.3, вызванное аэродинамической нагрузкой). Предпочтительная процедура обнуления включает подгонку усредняющей кривой к данным о смещении и вычитание этой кривой из начальных данных о смещении.
Фиг 11 иллюстрирует данные о смещении для пары зонд/лопатка. Четыре области 20 указаны и имеют предварительную идентификацию как события резонансной вибрации. Если усредняющая кривая подогнана ко всей последовательности данных и затем вычтена из всей последовательности, существует риск того, что интересующие сигналы и в частности синхронные отзывы будут ослаблены.
Таким образом, усредняющая кривая только вычитается из данных вне идентифицированных событий или, по крайней мере, только вычитается из данных вне идентифицированных синхронных резонансных событий.
Фиг.12 отображает кривую 21, полученную усреднением данных фиг.11, и дополнительная кривая 22, являющаяся результатом вычитания усредняющей кривой 21 из данных фиг.11 вне 4-х областей 20, которые идентифицированы как резонансные вибрации.
Путем наложения линейной интерполяции на усредняющую кривую 21 в 4-х областях 20, кривая 21 повышается для покрытия всей последовательности данных о смещении фиг. 10.
Повышенная кривая 21 затем может быть вычтена из данных фиг. 11, давая обнуленные данные, показанные на фиг. 13.
Усредняющая кривая может генерироваться путем прохождения данных о смещении через ядро свертки, которое оптимизировано для данных ВТТ. Длина ядра определяет число выборок до того, как 1-й выход появляется на своем выходе, и отсюда его время ожидания. Фиг. 14 - это блок-схема ядра. Фильтры Савицкого-Голая могут дать подходящие ядра. Длина фильтра (n на фиг.14) для создания усредняющей кривой может быть, например, в диапазоне 100-400. Полином Савицкого-Голая может быть 4-го порядка.
Обнуленные данные о смещении затем фильтруются адаптивным процессом согласно тому, является ли событие синхронным или асинхронным. Во-первых, FFT проводится на данных смещения каждого зонда для преобразования данных в частотную область, причем набор преобразованных данных для каждого зонда предпочтительно включает сигналы от всех лопаток ротора. Фиг.15 схематически показывает данные о смещении 30 при вибрации и соответствующие данные 31 о частотной области.
Количество дискретных бункеров в данных частотной области находится в пропорции к количеству выборок, используемых для заполнения входного буфера FFT. Количество образцов также определяет время запаздывания между получением данных их обработкой. Для обработки в реальном времени данные поступают во входной буфер FFT с частотой, равной частоте OPR, умноженной на число лопаток.
Данные частотной области 31 имеют два больших пика. Когда такие данные повторно получаются при последующих этапах времени, большие пики дают отслеженные разряды. Фиг.16 - это схематический график частота - время и показывает отслеженный разряд как толстую непрерывную наклонную линию.
Фиг.16 отображает также интегральные положения выборки как горизонтальные пунктирные линии. Интегральные положения - это частоты, соответствующие частоте вращения ротора их кратным до кратного, соответствующего числу лопаток в узле. Отдаление друг от друга интегральных положений выборки определяется делением числа точек в FFT на число лопаток на роторе.
Во временной области теоретические данные о смещении вибрации, зарегистрированные на одном зонде и для одной лопатки, могут принять форму, показанную на фиг.17а, хотя шумы вызывают данные о реальном смещении, более похожие на показанные на фиг.17b. По мере того как отзыв приближается к интегральному положению выборки, паразитная частота снижается до постоянного тока и затем возрастает по мере движения от интегрального положения выдержки. В основном, максимально интересующая точка - это точка постоянного тока, но это там, где шумы имеют самый определенный эффект.
Фильтрация должна быть применена для снижения шумов в данных, но недостаточно выбранная природа данных означает, что разные компоненты данных должны иметь разные коэффициенты фильтрации.
Таким образом, 2-й этап адаптивного процесса фильтрации состоит в фильтрации данных о смещении выбранного зонда при каждой идентификации событий вибрации во временной области, причем характеристика фильтра изменяется согласно закону контроля, который зависит от положений отслеженного порядка событий относительно интегральных положений выборки. Фильтр может быть, например, фильтром Баттерворта или Бесселя, причем весовые коэффициенты фильтра меняются во времени в зависимости от относительного положения соответствующих отслеженных порядков. Взвешивания эффективно контролируют частоту отсечки фильтра и отзыва. Значения взвешивания можно определить экспериментально или по закону математики.
Возвращаясь к временной области, фиг.18а иллюстрирует эффект закона контроля, применяемого к линейному изменению в частоте отсечки положения бункеров, наложенных на данные фиг.17b, и фиг.18b иллюстрирует эффект закона контроля экспоненциального изменения в частоте отсечки с положением бункеров. В обоих случаях шумы уменьшены, чем значительно улучшена неопределенность измерений. Далее фильтр адаптируется к содержанию спектра данных о смещении, которое определяется анализируемой частотой сигнала.
Фиг.19а отображает график частоты относительно времени для реальных данных. Отслеживаемый порядок, соответствующий вибрации, показан жирной линией, последующие секции отслеживаемого порядка названы А, В, С. В секции А, которая представляет собой большинство событий, отслеживаемый порядок расположен между двумя интегральными положениями выборки и фильтр типа А налагается на данные. К концу отслеживаемого порядка в секции В событие переключается на интегральное положение выборки, и другой фильтр типа В применяется к данным. Наконец, в секции С другой фильтр типа С применяется, когда отслеженный порядок переключается на другое положение бункера (хотя другое положение бункера трудно визуально идентифицировать). Переключения между типами фильтров можно достигнуть при динамическом изменении коэффициентов фильтров или при нескольких фильтрах, действующих постоянно и просто выбирающих требуемый выход согласно отслеживаемому порядку положения бункера относительно интегральному положению выборки. Фиг.19b-с - это графики смещения по времени (их временные шкалы выровнены со шкалами фиг.19а) и показывают соответственно эффекты разных фильтров А-С на данные смещения, собранные одним зондом от одной лопатки ротора за период времени отслеживаемого порядка.
После фильтрации данных они могут затем использоваться для выделения амплитуд и частоты отдельных лопаток. Методы такого извлечения известны для специалиста. Использование следящего фильтра, как показано, снижает неопределенность до 10%.
Хотя изобретение описано в сочетании с примерными вариантами воплощения, описанными выше, очевидны многие эквивалентные модификации и изменения, для специалистов раскрытые здесь.
Соответственно, примерные варианты воплощения изложенного выше изобретения рассматриваются как иллюстративные и не ограничивающие. Можно сделать различные изменения в описанные варианты изобретения, не отходя от существа и объема изобретения.

Claims (19)

1. Способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток, установленных на роторе, и ряд отстоящих друг от друга по периферии стационарных зондов таймирования, связанных с лопатками, которые обнаруживают моменты, когда лопатки проходят соответствующие зонды, причем способ включает следующие этапы:
- получение тайминга лопаток, обнаруженного зондами;
- определение для последовательных вращений (оборотов) узла соответствующих факторов корреляции для одной или более лопаток, причем каждый фактор корреляции определяет величину степени корреляции между таймированием лопаток, обнаруженным зондами для конкретной лопатки на конкретном вращении и таймированием лопаток, обнаруженным зондами на предыдущем вращении; и
- определение события резонансной вибрации, когда один или более факторов корреляции пересекает установленный порог.
2. Способ по п.1, дополнительно включающий этап обработки таймирования лопатки путем выполнения следующих подэтапов:
выбор одного из зондов таймирования;
подбор усредняющей кривой к таймированию лопаток выбранного зонда вне определенных событий резонансной вибрации;
интерполирование секций для усредняющей кривой у таймирования лопатки, выбранного зонда внутри определенных событий резонансных вибраций;
дополнение усредняющей кривой с интерполированными секциями; и
вычитание повышенной усредняющей кривой из тайминга лопаток выбранного зонда для получения обнуленного таймирования лопаток для выбранного зонда.
3. Способ по п.1 или 2, дополнительно включающий этап фильтрации таймирования лопаток путем выполнения следующих подэтапов:
- выбор одного из зондов таймирования;
- преобразование таймирования лопатки выбранного зонда в частотную область, причем преобразованное таймирование выдает отслеженные разряды в событиях резонансных вибраций в графике частоты относительно времени;
- определение интегральных положений выборки как частот, соответствующих частоте вращения ротора и их кратных до кратного, соответствующего количеству лопаток в узле;
- для определенных событий резонансной вибрации, определение отслеженных порядковых позиций относительно интегральных положений выборки и;
- фильтрация таймирования лопаток выбранного зонда определенных событий резонансной вибрации во временной области, характеристика фильтра, меняющегося в зависимости от соответствующих отслеженных относительных (порядковых) положений (позиций) разрядов в частотной области.
4. Способ по п.1 или 2, где этап получения содержит обнаружение моментов, когда лопатки проходят соответствующие зонды.
5. Способ по любому из пп.1, 2, где это или каждое идентифицированное (определенное) событие резонансной вибрации - это событие синхронной резонансной вибрации.
6. Способ по любому из пп.1, 2, где этот или каждый зонд измеряет отклонения законцовки лопаток.
7. Способ по любому из пп.1, 2, где частоты событий вибрации недостаточно выбраны этим или каждым зондом.
8. Способ обработки таймирования лопаток стационарным зондом таймирования, связанного с узлом вращающихся лопаток, установленных на роторе, причем зонд обнаруживает моменты, когда лопатки проходят зонд, включающий следующие этапы:
- получение таймирования лопаток, обнаруженного зондом;
- определение одного или нескольких событий резонансной вибрации в данных таймирования;
- подбор усредняющей кривой к таймированию лопаток вне определенных событий резонансной вибрации;
- интерполирование секций для усредняющей кривой у таймирования лопаток внутри определенных событий резонансных вибраций;
- дополнение усредняющей кривой с интерполированными секциями; и
- вычитание увеличенной усредняющей кривой из таймирования лопаток для получения обнуленного таймирования лопаток.
9. Способ по п.8, где этап получения включает:
- обнаружение моментов, когда лопатки проходят зонд.
10. Способ по п.8, где это или каждое идентифицированное (определенное) событие резонансной вибрации - это событие синхронной резонансной вибрации.
11. Способ по п.8, где этот или каждый зонд измеряет отклонения законцовки лопаток.
12. Способ по п.8, где частоты событий вибрации недостаточно выбраны этим или каждым зондом.
13. Способ по п.8, дополнительно включающий этап фильтрации таймирования лопаток путем выполнения следующих подэтапов:
- преобразование таймирования лопатки в частотную область, причем преобразованное таймирование выдает отслеженные разряды в событиях резонансных частот в графике частоты относительно времени;
- определение интегральных положений выборки как частот, соответствующих частоте вращения ротора и их кратных до кратного, соответствующего количеству лопаток в узле;
- для идентифицированных событий резонансной вибрации определение отслеженных положений разряда относительно интегральных положений выборки; и
- фильтрация таймирования лопаток при определенных событиях резонансной вибрации, причем характеристика фильтра меняется в зависимости от соответствующих отслеженных положений разряда в частотной области.
14. Способ фильтрации таймирования лопаток, обнаруженного стационарным зондом таймирования, связанного с узлом вращающихся лопаток, установленных на роторе, причем зонд обнаруживает моменты, когда лопатки проходят зонд, причем способ включает следующие этапы:
- получение таймирования лопаток, обнаруженного зондом;
- идентификация одного или нескольких событий резонансной вибрации в данных таймирования;
- преобразование таймирования лопатки в частотную область, причем преобразованное таймирование выдает отслеженные разряды в событиях резонансных частот в графике частоты относительно времени;
- определение интегральных положений выборки как частот, соответствующих частоте вращения ротора и их кратных до кратного, соответствующего количеству лопаток в узле;
- для идентифицированных событий резонансной вибрации определение отслеженных положений разряда относительно интегральных положений выборки; и
- фильтрация таймирования лопаток для идентифицированных событий резонансных вибраций в области времени, причем изменяющаяся характеристика фильтра зависит от соответствующих отслеженных положений разрядов в частотной области.
15. Способ по п.14, где этап получения включает:
- обнаружение моментов, когда лопатки проходят зонд.
16. Способ по п.14, где это или каждое идентифицированное (определенное) событие резонансной вибрации - это событие синхронной резонансной вибрации.
17. Способ по п.14, где этот или каждый зонд измеряет отклонения законцовки лопаток.
18. Способ по п.14, где частоты событий вибрации недостаточно выбраны этим или каждым зондом.
19. Компьютерная система для выполнения способа по любому из пп.1-18.
RU2009123249/28A 2008-06-18 2009-06-18 Способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток RU2514061C2 (ru)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GB0811073.6 2008-06-18
GBGB0811073.6A GB0811073D0 (en) 2008-06-18 2008-06-18 Timing analysis

Publications (2)

Publication Number Publication Date
RU2009123249A RU2009123249A (ru) 2010-12-27
RU2514061C2 true RU2514061C2 (ru) 2014-04-27

Family

ID=39672407

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2009123249/28A RU2514061C2 (ru) 2008-06-18 2009-06-18 Способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток

Country Status (7)

Country Link
US (1) US9341511B2 (ru)
EP (3) EP2136189B1 (ru)
AT (1) ATE505715T1 (ru)
DE (1) DE602009001056D1 (ru)
GB (1) GB0811073D0 (ru)
IL (1) IL198899A (ru)
RU (1) RU2514061C2 (ru)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007051027A1 (de) * 2007-10-25 2009-04-30 Mtu Aero Engines Gmbh Strömungsmaschine, Spaltmesssystem und Verfahren zum Ermitteln eines Rotorspaltes
GB201004559D0 (en) 2010-03-19 2010-05-05 Rolls Royce Plc Rotating blade analysis
GB2491632B (en) 2011-06-10 2013-10-30 Rolls Royce Plc Rotating blade analysis
GB201203181D0 (en) 2012-02-24 2012-04-11 Rolls Royce Plc Blade tip timing uncertainty
US9395270B2 (en) 2012-10-19 2016-07-19 Florida Power & Light Company Method and system for monitoring rotor blades in combustion turbine engine
GB201305199D0 (en) * 2013-03-21 2013-05-01 Rolls Royce Plc Blade Tip Timing
GB201309623D0 (en) 2013-05-30 2013-07-10 Rolls Royce Plc Blade tip timing
GB201309624D0 (en) 2013-05-30 2013-07-10 Rolls Royce Plc Blade tip timing
GB201309622D0 (en) 2013-05-30 2013-07-10 Rolls Royce Plc Blade tip timing
GB201402419D0 (en) * 2014-02-12 2014-03-26 Rolls Royce Plc Time reference derivation from time of arrival measurements
CN105372029A (zh) * 2014-09-01 2016-03-02 中国航空工业集团公司沈阳发动机设计研究所 一种基于叶尖定时技术发动机转子叶片振幅数据识别方法
CN109716077B (zh) * 2016-06-27 2021-04-27 比勒陀利亚大学 使用叶尖定时(btt)监测涡轮机转子叶片的方法和系统
EP3483571B1 (en) 2017-11-10 2020-04-29 Rolls-Royce Corporation Displacement determination using optical measurements
CN108051078B (zh) * 2017-12-12 2020-04-24 湖南工业大学 一种转速非恒定时叶片振动叶端定时在线监测方法及装置
US11091999B2 (en) 2018-06-12 2021-08-17 Probe Technology Services, Inc. Methods and apparatus for cement bond evaluation through production tubing
CN108871543A (zh) * 2018-08-30 2018-11-23 天津大学 恒速下叶片异步振动频率非调和傅里叶分析方法
CN109883379B (zh) * 2019-03-22 2020-07-10 西安交通大学 基于模态振型的叶片位移应变测量方法
CN109883389B (zh) * 2019-03-22 2020-08-25 西安交通大学 一种旋转叶片动应变场测量方法及其系统
CN109883380B (zh) * 2019-03-22 2020-03-31 西安交通大学 一种基于叶端定时的转子叶片位移场测量方法及其系统
CN110686764A (zh) * 2019-09-17 2020-01-14 天津大学 基于全相位差分原理的恒速叶片异步振动频率测量方法
CN112179480B (zh) * 2020-08-26 2021-11-19 西安交通大学 一种谐波叠加拟合叶尖定时振动参数识别的方法、系统及电子设备
CN112903271A (zh) * 2021-01-22 2021-06-04 中国航发沈阳发动机研究所 一种转子叶片非接触式异步振动参数辨识方法
CN113404555B (zh) * 2021-05-18 2022-06-21 西安交通大学 一种基于多个叶端定时传感器的叶片固有频率识别方法
CN113504310B (zh) * 2021-05-18 2022-12-02 西安交通大学 基于单个叶端定时传感器的叶片固有频率识别方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4573358A (en) * 1984-10-22 1986-03-04 Westinghouse Electric Corp. Turbine blade vibration detection apparatus
US5148711A (en) * 1990-11-01 1992-09-22 Westinghouse Electric Corp. Apparatus and method for removing common mode vibration data from digital turbine blade vibration data
US5206816A (en) * 1991-01-30 1993-04-27 Westinghouse Electric Corp. System and method for monitoring synchronous blade vibration
WO2004021298A2 (en) * 2002-08-30 2004-03-11 Seismic Warning Systems, Inc. Sensor apparatus and method for detecting earthquake generated p- waves and generating a responsive control signal
WO2007071912A2 (en) * 2005-12-21 2007-06-28 Rolls-Royce Plc Method of monitoring blades mounted on a shaft in a gas turbine engine
FR2900745A1 (fr) * 2006-05-05 2007-11-09 Eurocopter France Procede et dispositif de diagnostic d'un mecanisme

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4335600A (en) * 1980-11-13 1982-06-22 General Electric Company Detecting internal abnormalities in turbines
US4790189A (en) * 1987-03-02 1988-12-13 Westinghouse Electric Corp. Method for eliminating sensor drift in a vibration monitoring system
US4887468A (en) * 1988-06-03 1989-12-19 Westinghouse Electic Corp. Nonsynchronous turbine blade vibration monitoring system
US5152172A (en) * 1989-03-23 1992-10-06 Electric Power Research Institute Operating turbine resonant blade monitor
RU2125247C1 (ru) * 1993-07-20 1999-01-20 Сименс АГ Мониторная система для представления колебательных состояний множества лопаток на вращающемся рабочем колесе
US5471880A (en) * 1994-04-28 1995-12-05 Electric Power Research Institute Method and apparatus for isolating and identifying periodic Doppler signals in a turbine
FR2805893B1 (fr) 2000-03-02 2002-05-03 Snecma Procede de mesure d'adherence d'un revetement sur un substrat
JP3530474B2 (ja) 2000-09-22 2004-05-24 三菱重工業株式会社 翼振動計測方法及びこれを用いた翼振動監視システム
GB2374670B (en) * 2001-04-17 2004-11-10 Rolls Royce Plc Analysing vibration of rotating blades
DE10152026A1 (de) * 2001-10-23 2004-02-19 Mtu Aero Engines Gmbh Warnung vor Pumpgrenze oder Schaufelschaden bei einer Turbomaschine
FR2866953B1 (fr) * 2004-02-27 2006-06-30 Electricite De France Procede et systeme de mesure de vibration a la peripherie d'un corps tournant
RU2326362C1 (ru) 2006-09-07 2008-06-10 Государственное образовательное учреждение высшего профессионального образования Самарский государственный аэрокосмический университет им. академика С.П. Королева Способ диагностирования повреждений и дефектов рабочих лопаток турбомашин
US7861592B2 (en) * 2007-02-27 2011-01-04 Siemens Energy, Inc. Blade shroud vibration monitor
KR101454456B1 (ko) * 2007-04-30 2014-11-03 누보 피그노네 에스피에이 터빈 블레이드의 특성화 방법 및 시스템

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4573358A (en) * 1984-10-22 1986-03-04 Westinghouse Electric Corp. Turbine blade vibration detection apparatus
US5148711A (en) * 1990-11-01 1992-09-22 Westinghouse Electric Corp. Apparatus and method for removing common mode vibration data from digital turbine blade vibration data
US5206816A (en) * 1991-01-30 1993-04-27 Westinghouse Electric Corp. System and method for monitoring synchronous blade vibration
WO2004021298A2 (en) * 2002-08-30 2004-03-11 Seismic Warning Systems, Inc. Sensor apparatus and method for detecting earthquake generated p- waves and generating a responsive control signal
WO2007071912A2 (en) * 2005-12-21 2007-06-28 Rolls-Royce Plc Method of monitoring blades mounted on a shaft in a gas turbine engine
FR2900745A1 (fr) * 2006-05-05 2007-11-09 Eurocopter France Procede et dispositif de diagnostic d'un mecanisme

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
M.Zielinski, G.Ziller, Noncontact Vibration Measurements on Compressor Rotor Blades, Measurement Science Technology, Volume 11, Number 7, July 2000, стр.847-856; . *

Also Published As

Publication number Publication date
ATE505715T1 (de) 2011-04-15
EP2136189A2 (en) 2009-12-23
RU2009123249A (ru) 2010-12-27
EP2237003A1 (en) 2010-10-06
GB0811073D0 (en) 2008-07-23
US20090319206A1 (en) 2009-12-24
EP2237002A1 (en) 2010-10-06
DE602009001056D1 (de) 2011-05-26
IL198899A (en) 2013-11-28
EP2136189A3 (en) 2010-03-03
EP2136189B1 (en) 2011-04-13
EP2237003B1 (en) 2014-07-09
US9341511B2 (en) 2016-05-17
IL198899A0 (en) 2010-02-17
EP2237002B1 (en) 2014-07-09

Similar Documents

Publication Publication Date Title
RU2514061C2 (ru) Способ определения событий вибрации с резонансной частотой в узле вращающихся лопаток
RU2470280C2 (ru) Способ обнаружения и автоматической идентификации повреждения подшипников качения
CN107505135B (zh) 一种滚动轴承复合故障提取方法及系统
EP2199764B1 (en) Timing analysis
US7031873B2 (en) Virtual RPM sensor
US7133801B2 (en) System and methodology for vibration analysis and condition monitoring
US7822580B2 (en) Method and a system for monitoring the condition and operation of periodically moving objects
JP2003528292A (ja) 振動解析によるベアリングの状態ベースのモニタリング
US20070032966A1 (en) System and methodology for vibration analysis and conditon monitoring
CN110987438B (zh) 水轮发电机变转速过程周期性振动冲击信号检测的方法
EP2532839B1 (en) Rotating blade analysis
TW200936999A (en) Knock signal detection in automotive systems
Lin et al. A review and strategy for the diagnosis of speed-varying machinery
JP5436477B2 (ja) エンコーダ解析装置
CN114486252B (zh) 一种矢量模极大值包络的滚动轴承故障诊断方法
JP2019045401A (ja) 部分放電診断装置および部分放電診断方法
CN115683644A (zh) 航空发动机双源拍振特征识别方法
JP3598266B2 (ja) 機器異常診断方法および装置
CN113567123B (zh) 旋转机械冲击类故障自动诊断方法
CN217358982U (zh) 一种信号处理装置
RU2787309C1 (ru) Способ идентификации мультисинусоидальных цифровых сигналов
CN114993649A (zh) 一种信号处理装置及方法
JPH0883265A (ja) 振動信号解析装置
CN117870857A (zh) 一种谐频噪声分量提取方法及提取效果评估方法
CN115855481A (zh) 综合传动振动信号不明成分的阶次共振检测方法

Legal Events

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

Effective date: 20200619