Проектирование, сооружение и эксплуатация систем трубопроводного транспорта
УДК 622.691.4
РЕЗУЛЬТАТЫ РАСЧЕТА ПАРАМЕТРОВ ЗАТУХАНИЯ МАНОМЕТРИЧЕСКИХ ТРУБЧАТЫХ ПРУЖИН
RESULTS OF CALCULATION OF MANOMETRIC TUBULAR SPRINGS OSCILLATION DUMPING PARAMETERS
Д. А. Черенцов, С. П. Пирогов, Ю. Д. Земенков, А. Ю. Чуба
D. A. Cherentsov, S. P. Pirogov, Yu. D. Zemenkov, A. Yu. Chuba
Тюменский государственный нефтегазовый университет, г. Тюмень
Государственный аграрный университет Северного Зауралья, г. Тюмень
Ключевые слова: параметры затухания колебаний; манометрическая трубчатая пружина Key words: oscillation damping parameters; manometric tubular spring
При измерении избыточного давления в условиях пульсации измеряемой среды, больших амплитуд и внешних вибраций в манометрах происходит изнашивание различных механизмов, повышение погрешности измерения, что может привести к его поломке [1, 2]. Эффективным методом снижения вибрации является усиление процессов трения в конструкции с помощью демпфирующей жидкости. Ниже представлены динамические модели и результаты расчетов параметров затухания манометрических трубчатых пружин (МТП).
В первой модели МТП представлена как система с двумя степенями свободы: ф — относительный угол раскрытия пружины и w — увеличения малой полуоси поперечного сечения. Для получения дифференциальных уравнений было использовано уравнение Лагранжа второго рода [3, 4]:
где Т — кинетическая энергия, — обобщенная сила, соответствующая потенциальным силам; <2йг — обобщенная сила сопротивления (рис.1).
Выражения для определения кинетической энергии, обобщенных сил, соответствующих потенциальным силам и силам сопротивления, получены в [5, 6]. После диф-
Рис. 1. Обобщенные координаты
70
Нефть и газ № 1, 2015
ференцирования получается система двух однородных дифференциальных уравнений второго порядка с постоянными коэффициентами:
агф 4- Ьгф 4- сгср 4- с31Л/ = О а2м> 4- Ь2м/ 4- с3(р 4- с2м/ = О
где а¿у — коэффициенты инерции, Ьц — коэффициенты диссипации, а — коэффициенты жесткости.
Решая характеристическое уравнение
а1а2х4 + (ахЬ2 + а2Ь1)х3 + (ахс2 + ЬХЬ2 + а2сх)х2 + (Ьхс2 + сгЪ2)х + (схс2 - с3с3) = О
Получим значения параметров затухания в комплексном виде
хг = —пг + ¿/схл х{ = — - ¿/сх I
х2 = —я2 + \к2 |
х2 — Т12 ¿/с^у
где пг ы п2 — коэффициенты затухания (положительные); кг и к2 — частоты затухающих колебаний; х{ и х^ — числа, сопряженные с соответствующими комплексными числами.
Во второй модели МТП представлена в виде тонкостенного изогнутого стержня [7], совершающего колебания в плоскости кривизны центральной оси, на рис. 2 показан бесконечно малый элемент, вырезанный из этого стержня.
Рис. 2. Элемент стержня
Были получены дифференциальные уравнения в радиальном (м/) и продольном (и) перемещениях [8]:
„ , д2\л/ дуу д2 ( (ди д2уу\) ( /ди \л
№ 1, 2015
Нефть и газ
71
т(ср)
д2и
где Нк -
_ Е]{(р)Кк{(р), (1-Д2)Д4 '
дер I =
ди дер
д2м/ ^йр1
д ( (ди
ЕБ((р)
Ш0 = рБ{ср)', Е — модуль упругости мате-
К (1-д2)Д2'
риала трубки; / (ср) — момент инерции сечения, зависящий от угловой координаты сечения (р; Кк((р) — коэффициент Кармана, зависящий от угловой координаты сечения ср; — площадь поперечного сечения трубки, зависящая от угловой координаты ср этого сечения; /и — коэффициент Пуассона.
Главные граничные условия при ф = 0: и(0) = 0; О)(0) = 0; ^ (0) = 0
Естественные граничные условия при (ср = у): М(у) = 0; А1(у) = 0; у) = 0; или в перемещениях:
(:
Л(г)-ё«|=»
дер
дер3
Решение производилось методом Бубнова — Галеркина [9, 10] и реализовано с помощью программного комплекса «Манометр» на языке MATLAB [10].
Сравнение результатов расчетов произведено на примере трех рядов пружин постоянного сечения с соотношениями полуосей поперечного сечения 2, 5 и 8. В каждом ряду пружины отличаются толщиной стенки, которая варьируется от 0,1 мм до 1,6 мм. Пружины характеризуются следующими параметрами: радиус трубки-заготовки г = 8 мм, радиус кривизны центральной оси Я = 50 мм, центральный угол пружины у = 2200, динамическая вязкость демпфирующей жидкости / = 1 Па с. Графики зависимостей частот колебаний от параметра кривизны и тонкостенности __
/¿о = л] 12(1 — V2) — показаны на рис. 3 и 4. Кгь
350 ЗОО 250 200 150
\ у а/Ь=2
\ _ ___ .
5
--- - —^Г4:
а/Ь=й
О 10 20 30 40 50 цО
Стержневая теория - ™ -Энергетическая теория
Рис. 3. Сравнение частот затухающих колебаний
72
Нефть и газ № 1, 2015
Рис. 4. Сравнение коэффициентов затухания
Из графиков видно, что во всем рассмотренном диапазоне трубчатых пружин значения частот колебаний, полученные по энергетической теории (первая динамическая модель) несколько выше значений, полученных по стержневой модели. Как показали расчеты, коэффициент затухания в большей степени зависит от вязкости демпфирующей жидкости, для разных отношений полуосей разность полученных значений не превышает 1 %.
Сравнение расчетных значений для одних и тех же образцов, полученных данными методами, показало, что оба метода хорошо согласуются между собой при значениях параметра кривизны и тонкостенности < 10. При увеличении значений этого параметра отклонение между параметрами затухания, полученными обоими методами, увеличивается.
Список литературы
1. Андреева Л. Е. Упругие элементы приборов. - М.: Машгиз, 1962. - 456 с.
2. Пирогов С. П. Манометрические трубчатые пружины. - СПб: ООО «Недра», 2009. - 276 с.
3. Яблонский А. А., Никифорова В. М. Курс теоретической механики. -М.: Высшая школа, 1966. -439 с.
4. Бидерман В. Л. Теория механических колебаний. - М.: Высшая школа, 1980. - 480 с.
5. Чуба А. Ю. Расчет собственных частот колебаний манометрических трубчатых пружин: дис. ... канд. техн. наук/ А. Ю. Чуба. - Тюмень, 2007. - 137 с.
6. Черенцов Д. А. Определение обобщенных сил сопротивления при колебаниях манометрических трубчатых пружин, при наличии сил вязкого сопротивления // Естественные и технические науки. - 2014. - № 4 (72). - С. 12-16.
7. Беляев Н. М. Сопротивление материалов. - М.: Наука, 1965.
8. Черенцов Д. А., Пирогов С. П., Дорофеев С. М. Математическая модель манометрической пружины в вязкой среде // Вестник тюменского государственного университета. - 2014. - № 7. - С. 234.
9. Флетчер К. Численные методы на основе метода Галеркина. -М.: Мир, 1988. - 352 с.
10. Иглин С. П. Математические расчеты на базе Ма^аЬ.- СПб: БХУ - Санкт-Петербург, 2005. - 640 с.
Сведения об авторах
Черенцов Дмитрий Андреевич, ассистент кафедры «Транспорт углеводородных ресурсов», Тюменский государственный нефтегазовый университет, г. Тюмень, тел. 89222652095, е-таИ: cheгentsovda@bk. ги
Пирогов Сергей Петрович, д. т. н., профессор кафедры «Прикладная механика», Тюменский государственный нефтегазовый университет, г. Тюмень, тел. 8(3452)905785, е-таИ: piro-goweov@yandex.ru
Земенков Юрий Дмитриевич, д. т. н., профессор, заведующий кафедрой «Транспорт углеводородных ресурсов», Тюменский государственный нефтегазовый университет, г. Тюмень, тел. 8(3452)201931, е-таИ: zemenkov@tsogu. ги
Чуба Александр Юрьевич, к. т. н., доцент, государственный аграрный университет Северного Зауралья, г. Тюмень, е-таИ: aleksandr-chuba@mail.ru
Information about the authors
Cherentsov D. A., assistant of the chair «Transport of hydrocarbon Resources», Tyumen State Oil and Gas University, phone: 89222652095, e-mail: cherentsovda@bk.ru
Pirogov S. P., Doctor of Engineering, professor of the chair «Applied mechanics», Tyumen State Oil and Gas University, phone: 8(3452)905785, e-mail: piro-goweov@yandex. ru
Zemenkov Yu. D., Doctor of Engineering, professor, head of the chair «Transport of hydrocarbon Resources», Tyumen State Oil and Gas University, phone: 8(3452) 201931, e-mail: zemenkov@tsogu.ru
Chuba A. Yu., Candidate of Science in Engineering, associate professor, State Agrarian University of Northern Zauralye, e-mail: aleksandr-chuba@mail. ru
№ 1, 2015
Нефть и газ
73