б) ротор с распределенными параметрами.

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

Результаты расчета вибрации одномассового ротора с трещиной и небалансом подробно представлены в [10]. В расчетной схеме J. Schmied и E. Kramer варьировали двумя параметрами: глубиной трещины и фазой небаланса. Авторы отмечают наличие ультрагармонических резонансов на частотах 0,33ω0, 0,5ω0, где ω0 – собственная частота модельного ротора. Спектр вибрации ротора помимо оборотной частоты содержит кратные гармоники. В [10] рассмотрены только 1z, 2z, и 3z гармонические составляющие. Отмечено заметное влияние фазы небаланса на амплитуду первой гармоники. При совпадении фаз небаланса и трещины амплитуда первой гармоники оказывается значительно больше по сравнению со случаем, если трещина и небаланс находятся в противофазе. Имеются определенные различия амплитудно-частотных характеристик (АЧХ) для случаев а) синфазного и б) противофазного расположения небаланса и трещины. В случае а), по мере роста частоты вращения от нуля до резонанса, амплитуда первой гармоники растет монотонно, достигая максимума при , а в случае б) амплитуда между и сначала падает, а затем растет, и достигает максимума при Фазо-частотные характеристики (ФЧХ) также имеют различия для случаев а) и б). Для а) при переходе ротора через резонанс фаза изменяется на p радиан (как в случае обычного ротора без трещины). В случае б) фаза совершает скачек на 2p радиан в точке, расположенной между и , а в окрестности = фаза плавно изменяет свое значение на радиан. Амплитуда 2 й и 3 й гармоник от фазы небаланса не зависят.

В [2] авторы также отмечают, что вибрация, вызываемая трещиной, проявляется в основном на критических частотах вращения (оборотной составляющей и частотах равных половине и трети от оборотной составляющей и частотах равных половине и трети оборотной частоты). На рабочей частоте вращения, если она отстроена от критической, амплитуды колебаний будут невелики при трещинах достигающих половины сечения вала. Чем больше жесткость и выше собственная частота ротора, тем менее интенсивны будут его колебания. Авторы также отмечают, что при глубине трещины 15-20% механизм “дыхания’ проявляет себя достаточно слабо даже для гибких роторов. Авторы представляют результаты, из которых видно, что на критической частоте вращения амплитуды вибросмещения идеально отбалансированного ротора могут достигать существенных величин при глубинах трещины менее 10% от диаметра. В качестве диагностических признаков в работе предложено использовать амплитуды 1й и 2й гармоник. Наиболее вероятной причиной их устойчивого роста, а так же изменения фазы, является развивающаяся трещина. Дополнительные признаки трещины в [2] предлагают определять при выбеге (развороте) ротора на критических частотах 1-го и 2-го рода в виде значительного возрастания амплитуд и изменения фаз 1й и 2й гармоник по сравнению с первоначальным уровнем.

В [11] представлены результаты расчетов двухопорного ротора с распределенными параметрами на анизотропных упруго-демпферных опорах. Авторами был использован метод конечных элементов (МКЭ). Присутствие трещины моделировалось с помощью КЭ переменной жесткости. Дефицит жесткости КЭ с трещиной определялся некоторой величиной s, которую авторы называют “фактором трещины”. Однако, в работе [11] s никак не связана с конкретной глубиной трещины, что не дает возможность проводить сравнение полученных результатов с экспериментом и результатам, полученными другими авторами. При расчетах s меняется в пределах 0-0,3 (большая величина s соответствует большей глубине трещины). В [11] приведены расчетные графики АЧХ и АФХ ротора с трещиной во вращающейся системе координат. Из графиков видно, что в спектре вибрации присутствуют различные гармоники. (от первой до третьей), для 2 й и 3 й гармоник существуют параметрические резонансы. АФХ для первой гармоники по мере роста частоты вращения изменяется хаотичными скачками. В работе проводится сравнение с экспериментальными данными из другой статьи. Отмечено их качественное совпадение.

В [6] расчетные исследования проведены для ротора с равномерно распределенными параметрами. Н.Г. Шульженко рассматривал абсолютно уравновешенный ротор турбоагрегата мощностью 500 МВт. Расчеты показали, что соотношение между амплитудами гармоник существенно зависит от угла раскрытия трещины g. При g< 0,05p с ростом номера гармоники уменьшается их амплитуда. При значительных размерах трещин эта закономерность не сохраняется. Автор отмечает, что, как правило, амплитуды всех гармоник растут при увеличении g до 0,25p, а затем некоторые из них уменьшаются, достигая минимума при g = 0,5 p, а другие при этих g достигают относительного максимума. Относительные минимума и максимума амплитуд наблюдается также при g = 0,75p, а при дальнейшем росте g амплитуды одних гармоник растут, других уменьшаются. В целом результаты [6] плохо согласуются с результатами приведенными в других работах [5,10,11,12]. Так график АЧХ для каждой из пяти гармоник в диапазоне частот от 0 рад/с, до 300 рад/с изменяются хаотично. В указанном диапазоне частот АЧХ имеет много пиков и впадин. График зависимости между углом раскрытия трещины g и амплитудами гармоник также носит хаотичный (немонотонный) характер. При росте g на постоянной частоте вращения амплитуды всех рассматриваемых гармоник претерпевают существенные скачки. Кроме того, амплитуды гармоник достаточно велики по сравнению с амплитудой первой гармоники даже на частотах существенно отличающихся от частот параметрических резонансов. В других работах [5,10,11,12] соотношение амплитуд носит другой характер. Амплитуды высших гармоник почти всегда существенно меньше амплитуды первой оборотной составляющей, исключение составляют лишь зоны ультрагармонических резонансов.

В работе [8] рассматривался ротор с двумя дисками, которые расположены на небольшом расстоянии друг от друга, симметрично относительно центра ротора. Трещина находилась между дисками, то есть строго в центре ротора, и была выполнена электроискровым методом. Ротор опирался на шариковые подшипники. Между подшипниками и станиной находился слой резины. Авторами была составлена математическая модель данного ротора. В [8] представлены как расчетные, так и экспериментальные результаты. В представлены как расчетные, так и экспериментальные результаты. Полученные данные авторы приводят в виде годографов, т.е. на одном графике приводятся как АЧХ, так и ФЧХ. Данные представлены для 1й и 2й гармоник. Остановимся на расчетных результатах. При фиксированной глубине трещины T. Inagaki, H. Kanki и K. Shiraki варьировали фазу небаланса и сопоставляли годографы для ротора с трещиной и без нее, а также для постоянно раскрытой трещины. Рассмотрены как случай идеально отбалансированного ротора с трещиной, так и влияние небаланса. Полученные данные свидетельствуют, что фаза небаланса существенно влияет на амплитуду первой гармоники. Что касается второй гармоник, то на нее фаза оказывает незначительное влияние. Амплитуда второй гармоники в случае “дышащей“ трещины вдвое меньше, по сравнению со случаем постоянно раскрытой трещины. Изменение фазы при выбеге ротора происходит монотонно без скачков, что, по-видимому, обусловлено малым весовым прогибом модельного ротора и, вследствие этого, относительно малым влиянием трещины на вибрацию по сравнению с небалансом.

В [5] B. Grabowski приводит расчетные результаты для одномассового двухопорного ротора, двухопорного ротора с распределенными параметрами, на анизотропных опорах, системы турбина-генератор на 3-х опорах. Отмечено, что при росте трещины критическая частота для первой гармоники снижается, а амплитуда ее возрастает или снижается (в зависимости от фазы небаланса). Для одномассового ротора, на частоте вращения =, амплитуда первой гармоники растет монотонно с ростом трещины, в то время как амплитуда второй гармоники сначала медленно растет, а затем, приняв максимальное значение при глубине трещины » 35% от диаметра, начинает уменьшаться. Для двухопорного ротора с распределенными параметрами представлены графики амплитуд в вертикальном и горизонтальном направлении в зависимости от глубины трещины для правого, левого подшипников и для центра ротора (трещина находится ближе к правому подшипнику). Вторая гармоника при росте трещины сохраняет ту же тенденцию, что в случае с одномассовым ротором. В центре ротора амплитуда второй гармоники оказывается меньше, чем на подшипниках. Поведение же первой гармоники в центре ротора и на подшипниках существенно разнятся. Автор заключает, что амплитуды гармоник зависят от месторасположения трещины. Амплитуды тем выше, чем ближе находится трещина к месту максимального изгиба ротора. Так как ротор с распределенными параметрами имеет бесконечное множество собственных частот форм (на практике нас интересует лишь несколько первых), то каждой собственной форме соответствует свой закон изменения изгиба вдоль ротора. Амплитуда колебаний по некоторой собственной форме будет тем выше, чем больше величина ротора в месте трещины для этой собственной формы. Если же изгиб по данной форме трещины равен нулю, то колебания по этой форме не возбуждаются. Проведенные расчеты для ротора с распределенными параметрами позволяют проанализировать влияние на спектр вибрации не только глубины трещины и фазы небаланса, но и место расположения трещины вдоль оси ротора.

Анализ экспериментальных результатов В [13] представлены экспериментальные результаты, полученные на двухпролетном роторе, состоящем из двух валов, соединенных жесткой муфтой. Весь валопровод экспериментальной установки опирался на подшипники жидкого трения (всего 4 подшипника), суммарная длина валопровода составляла 2800 мм. На каждом пролете расположено по 3 диска, два из которых использовались для балансировки. На одном из пролетов был выполнен вырез, который инициировал развитие трещины. Далее трещину развивали, прикладывая статическую нагрузку к пролету. Диаметр ротора в месте трещины составлял 25 мм. Ротор с начальным вырезом принимался за эталон. Он назван авторами “ротор без трещины”. В [13] получены важные экспериментальные данные, показывающие влияние фазы небаланса на резонансные амплитуды гармоник для трещин разной глубины. Из приведенных результатов видно, что первая гармоника на первой критической скорости зависит от фазы небаланса. Зависимость эта тем сильнее, чем глубже трещина (где номер кривой соответствует определенной глубине трещины: 1 - вырезу глубиной h = 0,42 (h = h/R), 2 – трещине h=0,524, 3 - h = 0,622, 4 – h = 0,788, 5 – h = 0,905). Амплитуда второй гармоники от фазы небаланса не зависит для любой глубины трещины.

Многие авторы отмечают, что амплитуда второй, а тем более третьей гармоники вне зоны их резонанса чрезвычайно мала, и не может быть измерена штатной вибрационно-измерительной техникой. Трещины, расположенные вблизи подшипников, приводят к пропорционально меньшим изменениям вибрации, чем удаление трещины. Авторами работ [12,14] предложен “метод гистограмм”, позволяющий в процессе вращения ротора обнаружить трещину глубиной 1-2% от диаметра. Суть метода состоит в следующем. На первом этапе анализа происходит суммирование и осреднение большого числа исходных вибрационных сигналов, полученных при одинаковых условиях. При большом количестве реализаций, участвующих в процессе осреднения, суммирование позволяет устранить шумовой фон, имеющий случайную природу. В качестве одной реализации авторы предлагают рассмотреть сигнал, зарегистрированный в пределах одного оборота. Вибросигналы необходимо зацифровать. Проводится синхронное суммирование для ротора без трещины (эти данные будут храниться в качестве эталона). Операция суммирования периодически повторяется в целях контроля появления трещины. К усредненным данным, полученным таким методом, применяется преобразование Фурье. Полученный набор гармоник называют гистограммой гармоник. Данный метод, по мнению авторов, позволяет не только уменьшить шумовой фон, но и устранить те гармоники, которые присутствуют при работе ротора без трещины. Следовательно, полученный сигнал обусловлен трещиной. Результаты, полученные авторами, свидетельствуют о большой чувствительности метода. На основе данного подхода была создана система контроля трещин в роторах, которая успешно прошла испытания.

Выводы

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

Например.

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

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

Различие между результатами, полученными по предложенным методикам весьма существенны.

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

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

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

Список литературы

1. Papadopoulos С. A., Dimarogonas A. D. Vibration of Cracked (Shafts in Bending. // J. of Sound and Vibration 1983, Vol 91, N4, р.583-593.

2. Зиле А.3., Иэраилев K.Л., Руденко М. H. Особенности вибрационного проявления трещины ротора турбогенератора. // Электрические станции. – 1985. – № 4. – С.26-29.

3. Фролов К.В., Израилев Ю.Л., Махутов Н.А. и др. Расчет термо-напряжений и прочности роторов и корпусов турбин.М. Машиностроение, 1988, 240 с.

4. Changh Li, Bernschoni 0., Xenophotidis N. A General Approach of the Dynamic of Cracked Shaft // Trans of ASME J. of VIbr. – 1989. – July. – V.111. – P.527-263.

5. Grabowski B. The Vibrational Behaviour of a Turbine Rotor Containing– a Transverse Crack // Trans of ASME J. of Mech. – 1980. – V.102 Nl. – P.140-146.

6. Шульженко Г. H. Определение признака развитой трещина при изгибных колебаниях весомого ротора. // Проблемы машиностроения. 1990, том 34. ее.7-13.

7. Mayes LW.,Davies W. G. R. Analysis of the Response a Multi-Rotor Bearing System Containing a Transverse Crack in Rotor. // Trans of ASME J. of Vibration... - 1984 Vol 106, p.139-145.

8. Inagaki Т., Kanki H., Shiraki K. Transverse Vibrations of a General Cracked-Rotor Bearing System. // Trans of ASМЕ J. of Mech. 1982 Vol.104 april p.345-355.

9. Papadopoulos С. A., Dimarogonas A. D, Stability of Cracked Rotor in the Coupled Vibration Mode // Trans of ASME 1988 vol.1-10 July p.356-359.

10. Schmied J. Kramer E. Vibrational behaviour of a rotor with a crossectional crack // Int. conf. Vibr. Rot. Mach. Pap. Int. Conf-Heslidtion 11-13. Sept London, 1984. – P.183-192.

11. Nelson H. D., Natarai C. The Dynamic of Rotor System with a Cracked Shaft. // Trans of ASME J. of Vibr. 1986 april p.189-196.

12. Imam I., Azzaro 3. H., Bankert R. J., Schibel J. Development of an On-line Rotor Crack Detection and Monitoring System. // Trans of ASME J - of Vibr. july 1989 Vol.11-1 p.241-250.

13. Mayes I. W., Davies W. G. R. The Vibrational Behavior of a Multi-Shaft, Multi-Bearing System in the Presence of a Propagating Transverse Crack // ASME Paper 83–DET–82. –1983.

14. Imam I. Method for on-line Detection of Incipient Cracks in Turbine - Generators Rotors. // US Patent No 4,408,294 dated October 4, 1983.