|
(21), (22) Заявка: 2006123601/28, 03.07.2006
(24) Дата начала отсчета срока действия патента:
03.07.2006
(46) Опубликовано: 27.03.2008
(56) Список документов, цитированных в отчете о поиске:
RU 2272990 C1, 27.03.2006. RU 2231750 C1, 27.06.2004. RU 2138012 C1, 20.09.1999. US 2006122798 A1, 08.06.2006. US 5097711 A, 24.03.1992.
Адрес для переписки:
443020, г.Самара, ул. Садовая, 61, Институт проблем управления сложными системами РАН
|
(72) Автор(ы):
Боровик Сергей Юрьевич (RU), Райков Борис Константинович (RU), Секисов Юрий Николаевич (RU), Скобелев Олег Петрович (RU)
(73) Патентообладатель(и):
Институт проблем управления сложными системами Российской академии наук (RU)
|
(54) СПОСОБ ОБНАРУЖЕНИЯ КРУТИЛЬНЫХ И ИЗГИБНЫХ СМЕЩЕНИЙ ТОРЦОВ ЛОПАТОК РАБОЧЕГО КОЛЕСА ОСЕВОГО КОМПРЕССОРА ПРИ ИССЛЕДОВАНИЯХ СРЫВНЫХ ЯВЛЕНИЙ
(57) Реферат:
Изобретение относится к измерительной технике. Технический результат: упрощение технических средств бесконтактного измерения смещений торцов лопаток компрессора, снижение числа установочных отверстий и площади нарушенной поверхности статора, необходимой для установки преобразователей, повышение информативности измерений за счет учета большего числа координатных составляющих пространственных смещений торцов лопаток. Сущность: устанавливают один кластерный одновитковый вихретоковый преобразователь с тремя чувствительными элементами в виде отрезка проводника. Чувствительные элементы располагают таким образом, что проводники образуют равносторонний треугольник. Производят одновременный импульсный опрос измерительных цепей один раз за оборот ротора компрессора в момент, когда центр замка лопатки оказывается под точкой начала системы отсчета, соответствующей геометрическому центру треугольника, образованного чувствительными элементами. Фиксируют цифровые коды измерительных каналов C1, С2 и С3, преобразователя. Коды C1, C2 и С3 дополняют результатом моделирования радиального зазора между торцами лопаток и статором, полученного путем расчета деформаций рабочего колеса и статора под действием центробежных сил и температуры с использованием измерительной информации о частоте вращения ротора турбомашины и температурах рабочего тела (газа) и статора, поступающей с аппаратно-программных средств, входящих в состав стендового оборудования. По значениям C1, С2 и С3 и результату моделирования путем решения системы из трех уравнений, полученных на основе семейств градуировочных характеристик, найденных экспериментально, определяют фактические на момент опроса измерительных каналов значения угла разворота торца лопатки, смещение его в осевом направлении и в направлении вращения. 4 ил.
Изобретение относится к измерительной технике и может быть использовано для бесконтактного обнаружения крутильных и изгибных смещений торцов лопаток рабочего колеса осевого компрессора при исследовании в нем срывных явлений, например срывного флаттера.
Способ позволяет обнаруживать смещение торцов лопаток относительно квазистатического, регулярно повторяющегося его положения, предписанного принципом действия осевого компрессора при аналогичных параметрах режима работы, но без срывных явлений.
Целью предлагаемого изобретения является упрощение технических средств бесконтактного измерения смещений торцов лопаток компрессора, снижение числа установочных отверстий и площади нарушенной поверхности статора, необходимой для установки преобразователей, повышение информативности измерений за счет учета большего числа координатных составляющих пространственных смещений торцов лопаток.
Известен способ измерения параметров движения торцов лопаток ротора турбомашины и устройство для его реализации [Патент 2231750 РФ, МПК G01В 7/14, 7/30, G01Р 3/50 – №2001128250/28; Способ измерения параметров движения торцов лопаток ротора турбомашины и устройство для его реализации / Боровик С.Ю., Райков Б.К., Секисов Ю.Н., Заявл. 18.10.2001; опубл. 27.06.2004, бюл. №18], заключающийся в том, что с торцами лопаток ротора работающей турбомашины вводят во взаимодействие один дифференциальный вихретоковый преобразователь с двумя чувствительными элементами, выполненными в виде отрезка проводника, и по параметрам выходного сигнала преобразователя оценивают линейную скорость, период вращения ротора, осевое смещение ротора, угол раскрутки пера лопатки относительно начальной установки, зазор между торцом лопатки и торцом преобразователя. Недостатком указанного способа является то, что он не позволяет определить составляющую изгибного смещения лопатки вдоль оси вращения ротора турбомашины.
Предлагается способ обнаружения как крутильных, так и изгибных смещений торцов лопаток рабочего колеса осевого компрессора при исследовании в нем срывных явлений, например срывного флаттера, который заключается в том, что то с торцом каждой лопатки работающего компрессора вводят во взаимодействие кластерный одновитковый вихретоковый преобразователь с тремя чувствительными элементами в виде линейных отрезков проводников, расположенными в плоскости торца преобразователя таким образом, что проводники образуют равносторонний треугольник. Кластерный одновитковый вихретоковый преобразователь устанавливают на статоре осевого компрессора таким образом, что положение одного чувствительного элемента оказывается совпадающим с направлением оси вращения компрессора.
Чувствительные элементы кластерного вихретокового преобразователя включают через согласующие трансформаторы в составе преобразователя в измерительные цепи измерительных каналов. Импульсный опрос измерительных цепей производится одновременно один раз за оборот ротора компрессора в момент, когда центр замка лопатки оказывается под началом системы отсчета, соответствующей геометрическому центру треугольника, образованного чувствительными элементами кластерного вихретокового преобразователя. Опрашиваются все измерительные каналы кластерного вихретокового преобразователя одновременно и фиксируются соответствующие цифровые коды C1, C2 и С3.
Результаты преобразования выходных параметров кластерного вихретокового преобразователя (коды C1, C2 и С3) дополняют результатом моделирования радиального смещения Yл (радиального зазора). При возникновении срывных явлений изменения Yл сравнительно невелики. Их моделирование производится на основе расчета упругих и термических деформаций ротора и статора осевого компрессора по измеренным с помощью аппаратно-программных средств, входящих в состав стендового оборудования, значениям частоты вращения ротора компрессора, температуры рабочего тела (газа) в ступени компрессора и температуре статора компрессора.
По найденным значениям цифровых кодов C1, C2 и С3, и результату моделирования Yл, путем решения системы из трех уравнений, полученных на основе семейств градуировочных характеристик измерительных каналов кластерного вихретокового преобразователя, снятых предварительно перед испытаниями, определяют фактические на момент опроса измерительных каналов значения угла разворота торца каждой лопатки ( л), смещение его в осевом направлении (Хл) и в направлении вращения (Zл).
На фиг.1 представлена схематически ступень осевого компрессора система отсчета (X, Y, Z), где 1 – вал ротора осевого компрессора, 2 – ступица колеса, 3 – диск колеса ротора, 4 – лопатка, 5 – статор. В некоторый фиксированный момент времени ti материальная точка М, которая находится на торцевой поверхности и оси разворота лопатки, совмещается с осью Y координатной системы XYZ, начало отсчета которой (0) находится на поверхности неподвижного статора, а сама система жестко связана со статором. Именно в этот момент и производится фиксация кодов C1, C2 и С3.
На фиг.2 представлена схема, поясняющая размещение кластерного одновиткового вихретокового преобразователя на статоре осевого компрессора, где 3 – диск колеса осевого компрессора, 4 – лопатка, 5 – статор, 6, 7, 8 – первый, второй и третий чувствительные элементы кластерного одновиткового вихретокового датчика соответственно.
На фиг.3 представлена диаграмма, поясняющая этапы обнаружения крутильных и изгибных смещений торцов лопаток рабочего колеса осевого компрессора при исследованиях срывных явлений, где р – период вращения ротора компрессора, tгц1, tгц2, …, tгцnл – моменты времени, соответствующие прохождению центрами замков лопаток геометрического центра треугольника, образованного чувствительными элементами кластерного вихретоковоо преобразователя.
На фиг.4 представлены временные диаграммы, поясняющие опрос измерительных каналов кластерного одновиткового вихретокового преобразователя, где SYNC – импульс синхронизации, поступающий от стандартного датчика частоты вращения, входящего в состав стендового оборудования, С1, C2, С3 – коды, полученные в результате аналого-цифрового преобразования естественных выходных сигналов (индуктивностей) кластерного вихретокового преобразователя, tгц1, tгц2, …, tгцi…, tгцnл – моменты времени, соответствующие прохождению центрами замков лопаток 1, 2, …, i, …, nл под геометрическим центром треугольника, образованного чувствительными элементами кластерного вихретокового преобразователя.
Обнаружение крутильных и изгибных смещений торцов лопаток рабочего колеса осевого компрессора при исследованиях срывных явлений (например, срывного флаттера) предлагаемым способом осуществляется следующим образом. После подачи команды, инициирующей процесс измерения, по временному интервалу между двумя соседними импульсами датчика частоты вращения, входящего в состав стендового оборудования определяют период вращения ротора осевого компрессора p и рассчитывают моменты времени фиксации результатов преобразования всех параметров кластерного одновиткового вихретокового преобразователя, соответствующие моментам прохождения центрами замков лопаток tгц1, tгц2, …, tгцnл (nл – число лопаток осевого компрессора) геометрического центра треугольника, образованного чувствительными элементами кластерного вихретокового преобразователя.
Каждый импульс опроса обеспечивает подачу питания трех измерительных цепей с чувствительными элементами кластерного одновиткового вихретокового преобразователя и запускает процесс преобразования их выходных параметров в напряжения и далее – в цифровые коды C1, C2, С3 (фиг.4).
Учет неизмеряемой величины радиального зазора между торцом лопатки и статором осевого компрессора (Yл), влияющего на результат измерения углового смещения торца лопатки ( л), смещения его в осевом направлении (Хл) и в направлении вращения ротора (Zл), осуществляется путем моделирования деформаций ротора и статора осевого компрессора под действием центробежной силы и температуры рабочего тела (газа) внутри ступени компрессора. Моделирование осуществляется на основе информации о частоте вращения ротора компрессора, температуре рабочего тела (газа) внутри ступени компрессора и температуры статора, поступающей от аппаратно-программных средств, входящих в состав стендового оборудования.
Рассмотрим пример вычисления радиального зазора между торцами лопаток и статором осевого компрессора с помощью простейших моделей упругой и термических деформаций элементов конструкции силовой установки.
Пусть проводятся исследования осевого компрессора следующего типоразмера (фиг.1): радиус вала ротора – 0,07 м, радиус ступицы колеса ротора – 0,24 м, радиус диска колеса ротора – 0,35 м, длина лопатки – 0,05 м, внутренний радиус статора – 0,712 м. Предположим, что согласно измерительным данным, полученным от стендовой системы измерения частоты вращения ротора компрессора и температур в газовоздушном тракте и элементах конструкции компрессора, ротор компрессора вращается с частотой 3000 об/мин, а температуры рабочего тела (газа) внутри ступени компрессора и внутренней поверхности статора за время t=250 сек уменьшились относительно первоначальных на 50°С и 49°С соответственно.
Материал элементов конструкции статора и ротора был выбран идентичным – нержавеющая сталь, и характеризуется следующими параметрами: модуль Юнга – 18 ГПа, плотность – 7900 кг/м3, температурный коэффициент расширения – 11,5·10-61/°С.
Для рассмотренного выше режима работы осевого компрессора, изменение радиального зазора между торцом лопатки и внутренней поверхностью статора обусловлено упругими вытяжками элементов конструкции ротора (ступицы, диска, лопатки, см. фиг.1) под действием центробежной силы, термическим расширением элементов конструкции ротора и статора.
Расчет упругой деформации под действием центробежной силы осуществлялся отдельно для ступицы, диска и лопатки ротора осевого компрессора с использованием второго закона Ньютона с последующим интегрированием результатов в соответствии с известными зависимостями механики и теории прочности [Степин П.А. Сопротивление материалов: Учебник для немашиностроительных специальностей вузов. – 7-е изд. – М.: Высш. школа, 1983. – 303 с.]. Для указанного режима работы силовой установки и характеристик материала конструкции элементов ротора компрессора полная упругая деформация элементов ротора осевого компрессора под действием центробежной силы, равная сумме упругих деформаций ступицы, диска и лопатки, составила 0,087 мм.
Расчет термической деформации ступицы, диска и лопатки ротора осевого компрессора осуществлялся в предположении о равномерном распределении температуры по объему соответственно ступицы, диска и лопатки методом электротепловых аналогий [Ольсон Г. Динамические аналогии / Пер. с англ. Коробочкина Б.Л. – М.: Иностранная литература, 1947. – 224 с.]. Для указанного режима работы силовой установки и характеристик материала конструкции элементов ротора компрессора полная термическая деформация элементов ротора осевого компрессора, равная сумме термических деформаций ступицы, диска и лопатки, составила 0,237 мм.
Расчет термической деформации внутренней поверхности статора осевого компрессора осуществлялся с помощью упрощенной модели линейного расширения [Яворский Б.М., Селезнев Ю.А. Справочное руководство по физике для поступающих в вузы и для самообразования. – М.: Наука, 1989. – 576 с.]. Для указанного режима работы силовой установки и характеристик материала конструкции элементов статора компрессора полная термическая деформация внутренней поверхности статора осевого компрессора составила 0,401 мм.
Таким образом, согласно проведенным расчетам, радиальный зазор между торцом лопатки ротора и статором осевого компрессора за время, равное t, уменьшился на 0,251 мм и составил 1,749 мм (исходная величина зазора составляла 2 мм).
Фактические на момент опроса значения угла разворота торца каждой лопатки ( л), смещение его в осевом направлении (Хл) и в направлении вращения (Zл) находятся путем решения системы из трех уравнений, полученных на основе снятых экспериментально семейств градуировочных характеристик измерительных каналов, связывающих коды аналого-цифрового преобразования естественных выходных сигналов чувствительных элементов кластерного одновиткового вихретокового преобразователя для контролируемой лопатки (С1, С2, С3), найденные в один и тот же момент времени прохождения ее замка через геометрический центр треугольника, образованного чувствительными элементами кластерного вихретокового преобразователя, с углом разворота торца лопатки ( л), смещением его в осевом направлении (Хл), в радиальном направлении (Yл) и в направлении вращения (Zл):

Решение системы (1) осуществляется известными численными методами, например методом простых итераций [Самарский А.А., Гулин А.В. Численные методы. – М: Наука, 1989. – 432 с.]
Формула изобретения
Способ обнаружения крутильных и изгибных смещений торцов лопаток рабочего колеса осевого компрессора при исследовании срывных явлений, например срывного флаттера, заключающийся в том, что с торцами лопаток работающего компрессора вводят во взаимодействие композицию из трех одновитковых вихретоковых преобразователей с чувствительными элементами в виде отрезка проводника, расположенных определенным образом на статоре компрессора и включенных в измерительные цепи измерительных каналов через согласующие трансформаторы; импульсный опрос измерительных цепей с одновитковыми вихретоковыми преобразователями производят одновременно один раз за оборот ротора компрессора в момент, когда центр замка лопатки оказывается под точкой, соответствующей началу выбранной системы отсчета; по цифровым кодам, соответствующим значениям выходных параметров одновитковых вихретоковых преобразователей, путем решения системы уравнений, полученных на основе снятых перед испытанием семейств градуировочных характеристик преобразователей, представляющих собой зависимости цифровых кодов от линейных перемещений торцов лопаток по координатным осям X, Y и Z, получают оценку перемещений торцов лопаток в радиальном, осевом направлениях и в направлении вращения ротора компрессора, отличающийся тем, что применяют кластерный одновитковый вихретоковый преобразователь с тремя чувствительными элементами в виде отрезков проводников, расположенными в плоскости торца общего корпуса преобразователя таким образом, что проводники образуют равносторонний треугольник; устанавливают кластерный одновитковый вихретоковый преобразователь на статоре осевого компрессора таким образом, что положение одного чувствительного элемента оказывается совпадающим с направлением оси вращения компрессора; принимают за начало системы отсчета геометрический центр равностороннего треугольника, образованного чувствительными элементами кластерного одновиткового вихретокового преобразователя; моделируют радиальный зазор между торцами лопаток и статором компрессора, изменение которого обусловлено упругими и термическими деформациями элементов конструкции ротора и статора, на основе результатов измерения частоты вращения ротора компрессора, температуры рабочего тела (газа) в компрессоре и температуры статора компрессора, полученных от аппаратно-программных средств, входящих в состав стендового оборудования; вводят дополнительную координату – угол разворота лопатки, в семейство градуировочных характеристик измерительных каналов кластерного одновиткового вихретокового преобразователя; определяют фактические на момент опроса значения угла разворота торца каждой лопатки, смещение торца каждой лопатки в осевом направлении и в направлении вращения для полученных цифровых кодов в измерительных каналах с кластерным одновитковым вихретоковым преобразователем и рассчитанного с помощью модели радиального зазора между торцами лопаток и статором осевого компрессора.
РИСУНКИ
|
|