CC BY
23
88 ДИАГНОСТИКА
УДК 620.1
Оценка технического состояния узлов газоперекачивающих агрегатов как модель колебательной системы
А.А. Паранук
ассистент1
rambi.paranuk@gmail.com П.С. Кунина
д.т.н., заведующая кафедрой1
А.В. Бунякин
к.т.н, доцент1
М.К. Абессоло
аспирант1
■кафедра ОНГП КубГТУ, Краснодар, Россия
В работе представлен метод моделирования колебаний основных элементов газоперекачивающего агрегата — вала, подшипника скольжения с масляным слоем и их содержащего корпуса как линейно связанной механической системы. На основе математической модели такой системы определяются ее внутренние параметры, по отклонению которых от их нормальных значений можно судить о развитии неисправностей некоторых типов. В результате анализа значений этих параметров, соответствующих реальным эксплуатационным значениям, предлагается схема лабораторного моделирования радиально-упорного подшипника скольжения, учитывающая критерии физического подобия.
Материалы и методы
Математическое моделирование.
Ключевые слова
вибрационный спектр, диагностика, линейно-колебательная система
В настоящее время надежность и эффективная эксплуатация энергоемкого оборудования для транспортировки природного газа (газоперекачивающих агрегатов — ГПА) обеспечиваются выявлением и идентификацией на раннем этапе развития признаков возможных отказов. Эти признаки в виде отклонений измеряемых параметров и величин (амплитуд и частот виброскоростей и виброускорений при колебании агрегатов) характеризуют возникновение повреждений или дефектов оборудования, которые могут привести к аварийным ситуациям.
Совершенствование методов диагностики технического состояния газоперекачивающих агрегатов по параметрам вибрации позволит предотвратить отказы оборудования с тяжелыми последствиями.
Для анализа технического состояния ГПА (системы, состоящей из центробежного нагнетателя и газотурбинного двигателя, соединенных через муфты и мультипликатор) используется модель колебательной системы трех последовательно связанных гармонических осцилляторов, один из которых жестко закреплен [1, 2]. Эта идеальная механическая конструкция принята за основу для выявления закономерности в связях между собственными частотами при вибрации ГПА. Рассмотрим модель колебательной системы ГПА (рис. 1)
Условно обозначим: нижний элемент 1 — корпус; средний элемент обозначим 2 — подшипниковый узел (опорно-упорный подшипник, вкладыш, масляный слой); верхний элемент обозначим 3 — вал. Соответственно М, М2, М3 — массы элементов; К, К2, К3 — жесткости упругих соединений.
Опишем с помощью уравнения Лагранжа эту идеальную механическую систему:
х123 — смещения элементов от положений их равновесия; точкой обозначена — (й/Л) производная по времени.
Напишем уравнение движения механической системы в форме Лагранжа:
где 1 =1, 2, 3. И получим следующую систему:
(2)
т,х, + к л-, - к, (х,-х) = О \ MjXj + к: (х, - X,) - к( (х, -xj = 0 (з)
т,х, + к,х, — к, fir, — xj = О Подстановкой в эту систему вектора столбца решения получается матрица
№
(4)
которая используется для определения характеристических (собственных) частот.
т - ¿j I А .
k. лО * ■
О -А,
О
к, -к, -mjo' + ks
(5)
(6)
Частоты находятся из характеристического (векового) уравнения, которое выражает условие равенства нулю определяется этой матрицей:
(■ т,га! + к, + к;) (■■ т.оУ + к, (-т,^ + к,) - к,!Х~т,аУ + к,) -к/( ш^Н к, + к3) = О
Это уравнение третьей степени (кубическое) относительно ю2, ввиду того, что матрица симметричная, ее собственные числа [3, 4] (частоты в квадрате) вещественны, но кроме того, они должны быть еще и положительными. Уравнение приводится к стандартной записи кубического уравнения: х3 + ах3 + Ьх + с = 0 введены следующие обозначения:
х = ю2 (7)
Рис. 1 — Модель колебательной системы ГПА
(k, + kj o = - „ ■ +
(к, 4 kj
(8)
С- -
kl k; к) т, т, т.
После замены переменной У х~ ^ уравнение приводится к неполному виду
где р = - j + Ь, q
«• 'Г-iÜÍ
о
Это проверяется непосредственной подстановкой и использованием тожества
ens' г — ^ (cos 3 + 3 cos.О
(9)
Далее, ввиду симметричности матрицы, корни должны быть вещественными. Это выражается условием:
Это так называемый «неприводимый случай», в котором формулы Кардана для корней кубического уравнения могут иметь следующее (тригонометрическое) представление:
(10)
(И)
(12)
Естественно p < 0, а условие Q < 0 при этом равносильно \cos а\ < 1, случай Q = 0, то есть cos а = 1. Это соответствует кратным корням (особым).
Таким образом: _
(13)
(14)
Здесь \ *У~Т} в интервале
0<а<к. Для общего сведения тригонометрическое представление корней имеет вид:
то есть Q<0. Случай же Q=0, то есть cos а = 1, то есть а=0 и
Л Гр 2л
V, = V, - 2.1—J— oos — Л 13 3
Это кратность корней. Проверка выкладок может быть осуществлена таким образом: для произвольных и положительных mJ23; к123 по формулам (1), (2) и (3) вычисляются частоты (в квадрате ы212). Проверяется, что все эти величины положительны. Далее они подставляются (по очереди) в характеристическое уравнение и проверяется, что оно приближенно выполнено (~ 0).
Диагностические параметры — это жесткости к,, к, к3 — заранее неизвестны. Массы же в принципе известны (корпус, движущиеся части подшипника и вал). Частоты (ю123= 2пу12) измеряются вибрографом. Из уравнений (l), (2) и (3) находятся неизвестные к,2 3 каким-либо численным методом, например, встроенным в Excel или Mathlab.
Дальше проведем эксперимент на ради-ально-упорный подшипник скольжения.
На рис. 2 представлено вращения масляного слоя радиально-упорного подшипника скольжения. Слева показаны траектория и направление вращения с оси Y к оси X; справа разложение колебаний на оси Х по вертикали. Данные получены с компрессорной станции «Донская». Интерпретируя и анализируя этот график, мы сможем написать модель этих колебаний.
Из графика мы получили следующие величины:
• условный период колебаний масляного слоя — T = 103c;
• частота вращения масляного слоя — п = 3000 об/мин;
• угловую скорость вращения;
• разность перемещений Al = 2,5 • 106м.
Для наших расчетов условно возьмем следующие величины:
• масса вала центробежного компрессора m = 50 кг;
• площадь поперечного сечения вала s = 0,5 • 10-2м2
• максимальное перемещение масляного слоя L = 8ÁL = 0,002 мм.
Выразим условную жесткость колебаний «k» между валом и масляным слоем из формулы: ю2=к/ш и получим k=m®2=4040^50=240í> н/м.
Отсюда можно найти давление, которое оказывает масляной слой на поверхность вала:
P=f/s
где, f — сила, с которой масляной слой давит на вал.
Из закона Гука(.) известно, что f=kAl. Таким образом, мы окончательно получим:
р =к - =2-10' s
210"
■=8-Ю*Па
Рис. 2 — Вращение масляного слоя радиально-упорного подшипника скольжения
0,510 -(=0,8кгс/мм2).
По справочнику такие характеристики совпадают с каучуком. Поэтому для прогнозирования состояния радиально-упор-ного подшипника скольжения можно по подобию в лабораторных условиях провести испытания на [5, 6] каучук с помощью вынужденных колебаний и следить за тем как параметр жесткости будет меняться и какие при этом возникнут дефекты на поверхности каучука. Полученные опытные результаты можно сравнить с реальными данными с целью уточнения погрешности измерений.
Итоги
Создана работающая модель, которая испы-тывается на лабораторной установке и позволяет применить новые методы к прогнозированию состояния радиально-упорного подшипника скольжения.
Выводы
Описанные авторами новым подход и математически модели позволят проводить прогнозирования состояния радиально-упорно-го подшипника скольжения.
Список использдуеомй литературы
1. Зарицкий С., Стрельченко А., Тимофеев В., Бойко А., Шайхутдинов А. Вибромониторинг и диагностика - основа достоверной информации о состоянии ГПА // Газотурбинные технологии. 2000. № 5. С. 24-28.
2. Кунина П.С., Павленко П.П.. Диагностика газоперекачивающих агрегатов с центробежными нагнетателями. Ростов-на-Дону: РГУ, 2001. 362 с.
3. Кунина П.С. Бунякин А.В. Методы анализа спектров вибрации. Ростов-на-Дону: РГУ, 2001. 98 с.
4. Кунина П.С. Бунякин А.В. Математическое моделирование спектра вибрации механического модуля турбокомпрессора методом «слабого резонанса». Ростов-на-Дону: СКНЦ, №2, 2003 .
5. Трение, изнашивание и смазка. Под ред. Крагельского И.В., Алисина В.В. М.: Машиностроение, т. 1, 1978. 400 с.
6. Карасёв В.А. Вибрационная диагностика. ГТД. М.: Машиностроение, 1978. 342 с.
то есть
ENGLISH
DIAGNOSTICS
Evaluation of the technical condition of sectors of gas compressor units as a model of oscillating system
Author:
Arambi A.Paranuk — assistant1; rambi.paranuk@gmail.com Polina S. Kunina — Sc.D., head of the department1 Alexey V. Bunyakin —Ph.D., assistant professor1 Marten K. Abessolo — graduate student1
department "Equipment of oil and gas field' of Kuban State Technological University, Krasnodar, Russian Federation
UDC 620.1
Abstract
This paper presents a method for modeling the vibrations of the main elements of pumping unit — the shaft, bearings with the oil layer and their shell as a linearly related mechanical system. On the basis of a mathematical model of this system, we can define its internal parameters, and the deviation from their normal values can be used to discuss about the development of certain types of failures. By the result of the analysis of
these parameters, we suggest a scheme of laboratory simulation of journal-thrust bearing, using the criteria of physical similarity.
Materials and methods
Mathematical modeling.
Results
It was created a working model that is being tested in a laboratory setting allows you to apply new methods to predict the state of
angular contact ball bearing slide. Conclusions
The authors described a new approach and a mathematical model predicting the state will allow to carry out angular contact ball bearing slide.
Keywords
vibration spectrum,
diagnostics,
linear oscillating system
2.
References
Zaritskiy S., Strel'chenko A., Timofeev V., Boyko A., Shaykhutdinov A. Vibromonitoring i diagnostika - osnova dostovernoy informatsii osostoyanii GPA [Vibration monitoring and diagnostics - the basis of reliable information on the gas pumping unit]. Gas turbo technology, 2000, issue 5. pp. 24-28. Kunina P.S., Pavlenko P.P. Diagnostika gazoperekachivayushchikh agregatov s tsentrobezhnymi nagnetatelyami
4.
[Diagnostics of gas pumping units with centrifugal blowers]. Rostov-on-Don: Rostov State University, 2001, 362 p. Kunina P.S. Bunyakin A.V. Metody analiza spektrov vibratsii [Methods of analysis of vibration spectrs]. Rostov-on-Don: Rostov State University, 2001, 98 p. Kunina P.S. Bunyakin A.V. Matematicheskoe modelirovanie spektra vibratsii mekhanicheskogo modulya turbokompressora metodom «slabogo
6.
rezonansa» [Mathematical modeling of the vibration spectrum of mechanical turbocharger module by "weak resonance"]. Rostov-on-Don: SKNTS, issue 2, 2003. Trenie, iznashivanie ismazka [Friction, wear and lubrication]. Ed. by Kravel'skiy I.V., Alisin V.V.. Moscow: Mashinostroenie, Vol. 1, 1978, 400 p.
Karasev V.A. Vibratsionnaya diagnostika. GTD. [Vibration diagnostics. GTD]. Moscow: Mashinostroenie, 1978, 342 p.
СИГНАЛИЗАТОРЫ УРОВНЯ ЖИДКОСТИ ПОПЛАВКОВЫЕ «СУЖ-П-И»
ПРЕДНАЗНАЧЕН для сигнализации одного, двух или трёх предельных уровней нефти и нефтепродуктов, а также других технических жидкостей, включая воду, в резервуарах и технологических аппаратах. Сигнализатор обеспечивает предотвращение перетока жидких продуктов, фиксацию предельных уровней, снижение ущерба от аварии и повышение безопасности.
ИСПОЛНЕНИЕ ВЗРЫВОЗАЩИЩЕННОЕ
Модуль преобразователя вторичного предназначения для формирования искробезопасного напряжения и коммутации исполнительных устройств, имеет маркировку «[Ех1а] 11 В» и устанавливается вне взрывоопасных зон помещений и наружных установок.
ОБЛАСТЬ ПРИМЕНЕНИЯ СИГНАЛИЗАТОРА -
нефтебазы,технологические аппараты и резервуарные парки нефтеперерабатывающих заводов и производств. Принцип действия - срабатывание контактного устройства (геркона) при достижении чувствительным элементом (поплавком) заданного (контролируемого) уровня продукта. Сигнализатор СУЖ-П-И состоит из модуля преобразователя и преобразователей первичных.
Преобразователи первичные вертикального исполнения разделяются на разборные и не разборные, и служат для сигнализации одной, двух или трёх точек контроля.
Преобразователи первичные предназначены для сигнализации предельных уровней нефти и нефтепродуктов; имеют два вида конструктивного исполнения:вертикальное и горизонтальное.
Преобразователи первичные имеют маркировку взрывозащиты «0Ех1а118Т5:» и могут устанавливаться во взрывоопасных зонах помещений и наружных установок.
О 6 6
■ щ
ОТКРЫТОЕ АКЦИОНЕРНОЕ ОБЩЕСТВО
«КРАСНОЕ ЗНАМЯ»
390043, Россия, г. Рязань, пр. Шабулина, 2 а. +7 (4912) 938-517. post@kz.ryazan.ru
