автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математическое моделирование колебаний манометрических трубчатых пружин в вязкой среде

кандидата технических наук
Черенцов, Дмитрий Андреевич
город
Тюмень
год
2015
специальность ВАК РФ
05.13.18
Автореферат по информатике, вычислительной технике и управлению на тему «Математическое моделирование колебаний манометрических трубчатых пружин в вязкой среде»

Автореферат диссертации по теме "Математическое моделирование колебаний манометрических трубчатых пружин в вязкой среде"

ЧЕРЕНЦОВ ДМИТРИЙ АНДРЕЕВИЧ

На правах рукописи

Ъ

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ КОЛЕБАНИИ МАНОМЕТРИЧЕСКИХ ТРУБЧАТЫХ ПРУЖИН В ВЯЗКОЙ СРЕДЕ

Специальность 05.13.18. - Математическое моделирование, численные методы и комплексы программ.

11 НОЯ ¿015

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

005564^

Тюмень - 2015

005564333

Работа выполнена на кафедре «Транспорт углеводородных ресурсов» ФГБОУ ВО «Тюменский государственный нефтегазовый университет»

Научный руководитель: доктор технических наук, профессор

Пирогов Сергей Петрович

Официальные оппоненты: Соколов Владимир Григорьевич

доктор технических наук, профессор, ФГБОУ ВПО «Тюменский государственный архитектурно - строительный университет», заведующий кафедрой «Строительная механика»

Смолин Николай Иванович

кандидат технических наук, профессор, ФГБОУ ВПО «Государственный аграрный университет Северного Зауралья», заведующий кафедрой «Лесного хозяйства, деревообработки и прикладной механики»

Ведущая организация: ФГАОУ ВПО «УрФУ имени первого

Президента России Б.Н. Ельцина», г. Екатеринбург

Защита состоится «15» декабря 2015 г. в 14-00 часов на заседании диссертационного совета Д 212.274.14 при ФГБОУ ВО «Тюменский государственный университет» по адресу: 625003, г. Тюмень, ул. Перекопская, 15А, ауд. 410.

С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «Тюменский государственный университет» и на сайте http://d21227414.utmn.ru/alltransdef.

Автореферат разослан « 2? » ОУ^&ЬрЯ' 2015 г.

Ученый секретарь диссертационного совета

Е.А.Оленников

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность исследования. В Федеральном законе «О промышленной безопасности опасных производственных объектов» от 21.07.97 №116-ФЗ (ред. от 04.03.2013), «Концепции федеральной системы мониторинга критически важных объектов и (или) потенциально опасных объектов инфраструктуры Российской Федерации и опасных грузов» одобренной распоряжением правительства РФ от 27.08.2005 №1314-р обозначена необходимость совершенствования работ в области своевременного выявления и предупреждения угроз техногенного характера.

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

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

Цель исследования

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

Основные задачи исследования

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

2. Модификация численных методов для решения математических моделей колебательного процесса манометрических трубчатых пружин.

3. Разработка алгоритма и комплекса программ для автоматического определения параметров затухания колебаний.

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

5. Экспериментальное подтверждение достоверности полученных результатов.

Объектом исследования является чувствительный элемент приборов для измерения давления - манометрическая трубчатая пружина (МТП) в вязкой среде.

Предметом исследования являются параметры затухания колебаний МТП (коэффициент затухания колебаний и частота затухающих колебаний).

Научная новизна полученных результатов по трем областям специальности 05.13.18 - Математическое моделирование, численные методы и комплексы программ — отражена в следующем: В области математического моделирования:

• Впервые учтено влияние вязкости демпфирующей жидкости в известных математических моделях, основанных на уравнениях Лагранжа второго рода для определения параметров затухания колебаний МТП.

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

В области численных методов:

• Для впервые поставленной задачи модифицирован численный метод решения задачи колебаний манометрических пружин в вязкой среде на основе метода Бубнова - Галеркина.

В области создания и реализации комплексов программ: • Создан программный комплекс, для впервые поставленной задачи, ориентированный на численное решение задачи колебаний манометрических трубчатых пружин в вязкой среде. Программный комплекс прошел государственную регистрацию.

Практическая значимость исследования. С помощью разработанного метода расчета и пакета прикладных программ можно определять параметры затухания колебаний и другие технические характеристики МТП, тем самым совершенствовать конструкции вибростойких манометров. Результаты работы внедрены и используются в ОАО «Манотомь» (г. Томск).

Достоверность результатов, определяется корректностью постановки задач исследования, применением классических математических методов для построения решений, результатами численных экспериментов, а так же экспериментальными исследованиями на натурных образцах.

На защиту выносятся следующие результаты, соответствующие пунктам паспорта специальности 05.13.18 — Математическое моделирование, численные методы и комплексы программ по техническим наукам:

Пункт 2. Развитие качественных и приближенных аналитических методов исследования математических моделей.

1. Впервые учтено влияние вязкости демпфирующей жидкости в известных математических моделях, основанных на уравнениях Лагранжа второго рода для определения параметров затухания колебаний МТП.

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

Пункт 4. Реализация эффективных численных методов и алгориттюв в виде комплексов проблемно-ориентированных программ для проведения вычислительного эксперимента.

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

среде. Программный комплекс внедрен и используется в ОАО «Манотомь» (г. Томск).

Пункт 5: Комплексные исследования научных и технических проблем с применением современной технологии математического моделирования и вычислительного эксперимента.

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

Таким образом, в соответствии с формулой специальности 05.13.18 в диссертации представлены оригинальные результаты одновременно из трех областей: математического моделирования, численных методов и комплекса программ.

Апробация научных положений и результатов исследования

Основные положения и результаты исследования докладывались на следующих конференциях и семинарах: Проблемы функционирования систем транспорта: всероссийская научно-практическая конференция, посвященная 50-летию со дня основания Тюменского индустриального института (г. Тюмень, 2013); XII Научно-техническая конференция молодежи. 1-й тур. «АК Транснефть. ОАО «Сибнефтепровод». (г. Тюмень,

2013); ХХУП-ХХ1Х Международная заочная научно-практическая конференция — Научная дискуссия: вопросы технических наук. (г. Москва,

2014); На расширенном заседании кафедры «Прикладная механика» ТюмГНГУ (г. Тюмень, 2015); На научном семинаре кафедры «Транспорт углеводородных ресурсов» ТюмГНГУ (г. Тюмень, 2015); На научном семинаре кафедры «Основы теории механики и автоматического управления» ОмГТУ (г. Омск, 2015).

Публикации. По теме диссертации опубликовано 12 печатных работ, в том числе 6 в журналах из перечня ведущих периодических изданий, рекомендованных ВАК РФ, получено 2 свидетельства о регистрации программ для ЭВМ.

Структура и объем работы. Диссертационная работа состоит из введения, четырех глав, заключения, содержит 125 страниц машинописного текста, 16 таблиц, 60 рисунков, 3 приложения, библиографического списка использованной литературы из 131 наименования.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В первой главе рассмотрено устройство, принцип действия, приведены классификация и обзор наиболее распространённых конструкций МТП. Проведен обзор методов повышения виброзащиты, рассмотрены известные решения задач в области колебаний стержней и оболочек, дан обзор исследованиям в области определения коэффициента сопротивления. Показано, что методы расчета параметров затухания колебаний МТП, работающих в жидкости, отсутствуют. Кроме того отсутствуют результаты экспериментальных исследований колебаний в демпфирующей жидкости.

Все это позволило сформулировать цель и задачи исследования.

Во второй главе разработаны математические модели движения МТП в вязкой среде.

В первом параграфе приведена аналитическая модель, основанная на уравнениях Лагранжа второго рода.

В данной модели МТП представлена в виде механической системы с двумя степенями свободы, за обобщенные координаты приняты: относительный угол раскрытия пружины - ср и величина увеличения малой полуоси поперечного сечения трубки - IV.

Рис.1 Обобщенные координаты.

Система дифференциальных уравнений движения для системы с двумя степенями свободы, основанная на уравнениях Лагранжа второго рода, имеет известный вид:

г а /дТ\ ди дТ _ dt\дф) + дер д(р ®ф'

й /дЦ\ ди дТ _ (

чй + дй? ~~ дй? ~ где I - время; Т - кинетическая энергия, и - потенциальная энергия, <1<р> ~ обобщенные силы сопротивления по обобщенным координатам.

Выражение для определения потенциальной энергии было получено в [1,2] и имеет вид:

U =

Ehy

£lw2 _ 2b~<pw0 + b2A3<p2 +

(2)

где E - модуль упругости; ц - коэффициент Пуассона материала пружины;

Rh

X - характеристика манометрической пружины, Х =

S S S £

4 4 4 4

A^lj^ds; A2=¿Jycods; A3=¿Jy2ds; n = ±Jil2ds. 0 0 0 0 A1,A2lA3,m,n,ü. и ш - коэффициенты, зависящие от формы поперечного сечения МТП.

Выражение для определения кинетической энергии было получено в [2] и имеет вид:

Т = 2phR3Y (у - 4 sin у + 2y cos у + 2y j Вхф2 + 2phRy + w02, (3)

где p - плотность материала трубки; h - толщина стенки трубки; у - центральный

S S S

угол МТП; Bj = /04 ds ,В2- j*(o2ds, В2 = fluids, ws, К6 и m6 - коэффициенты, зависящие от формы поперечного сечения пружины.

Сила сопротивления движению в жидкости направлена по нормали к центральной оси МТП и пропорциональна её нормальной скорости.

1. Андреева Л. Е. Упругие элементы приборов / Л. Е. Андреева. - 2 изд., перераб. и доп. -Москва : Машиностроение, 1981. — 392 с.

2. Чуба А. Ю. Расчет собственных частот колебаний манометрических трубчатых пружин: дис.

... канд. техн. наук: 05.13.18 : защищена 13.11.07 / А. Ю. Чуба. - Тюмень, 2007. — 137 с.

Выражение для определения силы сопротивления получено с помощью метода неопределенных коэффициентов.

Р = СхБхуУр2уп (4)

где Сх — коэффициент лобового сопротивления; 5 - площадь миделевого сечения; р - плотность жидкости; и - скорость колебаний МТП; V - кинематическая вязкость; х, у, г и п - неизвестные показатели степени. Для симметричных тел коэффициент лобового сопротивления определяется формулой Стокса - Сх = 24//?е. Значения показателей степени были подобраны таким образом, чтобы колебания оставались линейными, выражение для определения силы сопротивления примет вид: 6ржЯ2(у2 \ . А .

а» \ 2 ) т

где а, - характерный линейный размер (расстояние между стенками колбы манометра и большой осью пружины), рж - динамическая вязкость жидкости; А = — величина, зависящая от формы сечения.

Выражения для определения обобщенных сил сопротивления:

-аф[- 2-+-i-¡ = -—{T + C0SYr-

(бл,"2'

<?»„« = V%-—+ ..... •/" ' I =

(6)

m

Продифференцировав потенциальную и кинетическую энергии, система (1) была получена в виде системы двух однородных дифференциальных уравнений второго порядка с постоянными коэффициентами:

( ахф + Ьгф + сг<р + c3w = О, (a2w + b2w + с3ф + c2w = 0. где а1(а2 - коэффициенты инерции, bltb2 - коэффициенты диссипации, d, с2, с3 - коэффициенты жесткости, зависящие от геометрических характеристик и физических свойств материала МТП:

at = 4phR3y (у - 4 sin у + 2у cos у + 2y"j Ви a2 = 4phRy +

6дж/?2 /у2 \ пА

Й! =-íí2— i-H-cosr , b2 = 12¡imyR —,

a \ 2 1 m

ЕНГ (А1 У2 п \ 2Я(1 -Ц2)\т2 12 т2)'

Если все 4 корня окажутся комплексными с отрицательными действительными частями, то будут наблюдаться затухающие колебания (малое вязкое сопротивление).

х-у = —щ + ¿к^

где щ и п2 - коэффициенты затухания (положительные); к± и к2 - частоты затухающих колебаний; х{ и х^ - числа, сопряженные с соответствующими комплексными числами.

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

Общие уравнение движения получено путем суммирования двух затухающих колебаний:

где £)1( 02, Ег и Е2 - амплитуды колебаний определяемые из начальных условий; а1,а2,р1и [12 - значения угла сдвига фазы, определяемые из начальных условий.

С помощью полученного метода можно аналитически определить первые два параметра затухания колебаний МТП, а так же оценить скорость затухания.

Во втором параграфе приводится модель, в которой МТП представляется как изогнутый стержень, совершающий колебания в плоскости кривизны центральной оси. Перемещение бесконечно малого элемента вырезанного из криволинейного стержня (рис.2) будет складываться из продольной и и поперечной IV составляющей.

Сила сопротивления движению в жидкости - распределенная нагрузка Я, направлена по нормали к оси стержня и, в связи с малыми скоростями движения, пропорциональна первой степени её нормальной скорости. Силы от внутреннего давления на верхнюю и нижнюю части поперечного сечения бесконечно малого элемента, представлены как распределённые нагрузки

= —щ —

х2 = —п2 + ¿к2 (' х2 = П2 1к2;

(8)

(р = фС1) + <р(2) = /^е-"1' этС/с^ + а+ 02е~ПгС + а2)

Шо = а)0Ю + Шо(2) = ^е""1' бшО^Г + + Е2е~п^ 5ш(Л2£ + /?2)

,} °6)

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

Л' -

д д<р 2 дМ й<р \

М~1крТ {

дМ й<р д<р 2

дМ Л®

Рис.2 Бесконечно малый элемент стержня. Система уравнений движения МТП в вязкой среде в перемещениях и и

IV имеет вид: 32м/

- .. / Ь\ Зи/ д2 ,

дс2

дс д<р2

ди 32иЛ) ( /ди 4 аЬ

(9)

где т(<р) - масса единичной длины трубки; Нк = £>я = ЕБУ 2;

32и 3 ( (ди Згш\1 3 г /ди \1

т 3^ - з^ГЧ^ " ^л " ^ Гя Ь+"Л=

__

Е - модуль упругости материала трубки; ]{<р) - момент инерции поперечного сечения трубки; 5(|р) - площадь поперечного сечения трубки.

Система (9) имеет следующие граничные условия: у основания манометрической пружины (в сечении, где пружина жестко закреплена) -Ф=0 продольные и поперечные перемещения, а также угол поворота поперечного сечения трубки равны нулю. На свободном конце ф=у, поперечные, растягивающие усилия, а так же изгибающий момент равны нулю.

В третьей главе разработаны численное решение системы (9), комплекс программ для ЭВМ. Здесь же исследовано влияние геометрических

характеристик МТП и вязкости демпфирующей жидкости на колебательное движение МТП в вязкой среде.

В первом параграфе приведен численный метод решения системы дифференциальных уравнений (9) основанный на методе Бубнова-Галеркина. Искомые перемещения представлены в виде:

(10)

и ~ T/;1(<p)a1(t) + i/>2(<P)a2(t) + "■tyn(fp)<hi(t) ~ ^ ^(»„(t),

é=I п

w » CiC^JbiCt) + fcfoO^CO + -í„(?0bn(t) « Y, Sn(.<p)bn{i),

i=l

где alt a2 ... an, blt b2 ... bn - неизвестные функции переменной t;

ipv ф2 ... xpv, <2 ... - заданные базисные функции переменной ф. Так как перемещения (10) не являются точным решением уравнений движения, то подставляя их в (9) вместо нуля получим невязки R, и R2. В соответствии с методом Бубнова - Галеркин требуется выполнение условий ортогональности невязок и R2 базисным функциям и t/j¿, которые в свою очередь должны удовлетворять главным граничным условиям. Заданные базисные функции переменной ср имеют вид:

= i = l,...,n.

v) = (pi+1; j = i.....п.

Функции переменой t имеют вид:

an(t) = сп • bn(t)=fnert,

где сп, /п и г — неизвестные постоянные.

Условие ортогональности сводится к системе 2п однородных дифференциальных уравнений с постоянными коэффициентами:

у у 7 2 2 у \

ггкшо Jw^ + r/i/?| СЛ^ + ЛС J ^T^dV + fiU j(i<id<p) + -

О 0 0 о /

(l' У У 2 2 У

r2fnmо J + r/„0 J CnCid-P + /„с J J^^dcp + fnH J ) + -

ООО о

У У У 2 2 У \

Аш01 МпЛ<р + г№ | + лс | ^г^г^ + л« / ) + -

О 0 0 о /

(У У У У \

г2/л | СпС^'Р + гМI + /„с I + №I <п<пЛч> I + •

О 0 0 о /

\ о о

У У

Л г[3фпд2<п. \ „

+ Сп I ^ "с" 1 з^ =

о о

о о

С»)

о о

,12 Г , . . г\дфхдф1 [дфгд^,

+ | г2схт0 J ¿<Р ~ с,С J - с1Н J —~^<1<р | + •

о

у

¿)(р з<р

. Г Г дфп 31/)! Г дфп дгрг .

+ ( г2спш0 ] ^ ^ - с„С ] - спН ] ] = 0.

о

},,, г(дф1дфп ¡дф1дфп \

0 0 о

У У У

Г , , ^ г С дфпдгрл Г д\рп дтрп

+ г

о

Условием существования нетривиального решения является равенство нулю главного определителя системы. Из этого условия находим г. Запишем определитель в виде:

^1,1 - Ci.ii

^2П,1 ■■■ С2пп

= о, (12)

где = г2т0 /; +гр ц ^ср+с +н /;

Г дфп Гдфпд

С1п = Н I —- ^.йер — С I —2 —г йср,

)д<р ) дф дер*

о о

у у у у

О / Ыг*Ч> + + с / ^Г^М* + Н1 <¿2**

ООО о

У У У

О [ Фф. с1ср — С ( ^-^йер - Н [ ^^ер.

^ } дер дер * ) дер дер ^

р2,п - Г2т.,

(■2п,п — Г2 то

Порядок определителя равен 2п. Элементы определителя содержат г2, значит, определитель эквивалентен алгебраическому уравнению степени 4п, поэтому уравнение с!е1(А)=0 должно иметь 4п корней.

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

Во втором параграфе определятся коэффициент Кармана для МТП с помощью полубезмоментной теории оболочек по выражениям предложенным Э.Л. Аксельрадом.

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

В четвертом параграфе приводится описание программного комплекса «МАНОМЕТР». На рис.3 показан интерфейс программы «Расчет затухающих колебаний».

Рис.3. Интерфейс программы «Расчет затухающих колебаний» В пятом параграфе производится сравнение результатов двух моделей. Полученные значения параметров затухания колебаний для одних и тех же образцов дают близкие значения при параметре кривизны менее 10, с увеличением параметра кривизны разница между параметрами затухания увеличивается.

В шестом параграфе оценивается влияние геометрических характеристик МТП и вязкости демпфирующей жидкости на параметры затухания колебаний. На рис.4(а-г) показано влияние геометрических характеристик МТП на коэффициент затухания.

а) Влияние толщины стенки б) Влияние отношения полуосей

в) Влияние центрального угла г) Влияние радиуса

Рис.4 Влияние геометрических характеристик МТП на коэффициент затухания

На рис.5(а-г) показано влияние геометрических характеристик МТП на частоты.

2 (Мазут) 1,5(Глицерин ) 0,1(Моторное масло)

0.1(Моторно« »«ело) г 1,5(Глицерин )

г(Маэут)

а) Влияние толщины стенки

б) Влияние отношения полуосей

ОД(Моторное масло) 1.5(Глицерин

200

Центральный угол, град.

Радиус продольной оси, мм

■ 2(Мэзут) 1,5(Глицерим) 0Д(Моторное масло)

в) Влияние центрального угла г) Влияние радиуса

Рис.5 Влияние геометрических характеристик МТП на частоты колебаний Исследование влияния геометрических характеристик МТП и свойств демпфирующей жидкости на параметры затухания колебаний показало, что увеличение толщины стенки ведет к увеличению частоты затухающих колебаний и уменьшению коэффициента затухания, а увеличение отношения полуосей, центрального угла и радиуса продольной оси ведет к уменьшению частот затухающих колебаний и увеличению коэффициента затухания. Частоты затухающих колебаний зависят от геометрических характеристик трубчатых пружин. Коэффициент затухания определяется динамической вязкостью демпфирующей жидкости и в меньшей степени зависит от геометрических характеристик трубчатых пружин.

В седьмом параграфе приводится оценка влияния вязкости демпфирующей жидкости на скорость затухания колебаний. В зависимости от динамической вязкости демпфирующей жидкости (табл.1) у 5 различных образцов МТП было получено время затухания колебаний. Результаты сведены в табл.2

Таблица 1

Значение вязкостен жидкостей

№ Жидкость Вязкость, Па с

1 Моторное масло БАЕ 40 0,319

2 Касторовое масло 0,985

3 Глицерин 1,49

4 Мазут 2,022

Таблица 2

Время затухания колебаний

Номер образца Время затухания колебаний, с

Моторное масло БАЕ 40 Касторовое масло Глицерин Мазут

1 11,1 3,2 2,3 1,7

2 21,4 6,1 4,4 3,4

3 29,8 8,5 6,1 4,7

4 53,6 15,3 11,1 8,5

5 47,9 13,6 9,8 7,5

Результаты расчетов показали, что при увеличении вязкости демпфирующей жидкости в 3-6 раза, время затухания колебаний уменьшается в 1,3-6,5 раз. Наиболее эффективным является применение глицерина, в силу его доступности и хорошей демпфирующей способности.

В восьмом параграфе произведено исследование апериодического движения. Для тех же пяти образцов было определено значение динамической вязкости, при котором наблюдается апериодическое движение. Результаты расчетов приведены в табл.3.

Таблица 3

Значение динамической вязкости при апериодическом движении, Па с

Номер образца 1 2 3 4 5

Значение

динамической 1000 1200 1200 1400 1600

вязкости

Расчеты показали, что значение динамической вязкости не соответствует реальным вязкостям жидкостей, см. табл. I, можно сделать вывод, что при демпфировании колебаний МТП жидкостью будут наблюдаться только затухающие колебания.

В четвертой главе приведены экспериментальные исследования параметров затухания колебаний МТП с целью определения достоверности созданных методов расчета.

В первом параграфе приведены описание опытных образцов МТП и приборов для определения параметров затухания колебаний МТП, описание стенда для измерения параметров затухания колебаний.

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

Рис.6. Стенд для измерения параметров затухания Сравнение значений показало, что отклонения частот от теоретических значений не превышали 10% для анализатора вибрации и 3%

для вискозиметра. Отклонения значений коэффициента затухания больше

для анализатора вибрации не более 15% и для вискозиметра не более 8%. Это

подтверждает хорошую точность теоретических результатов.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1. Учтено влияние вязкости демпфирующей жидкости в известных математических моделях, основанных на уравнениях Лагранжа второго рода для определения параметров затухания колебаний манометрических пружин в вязкой среде

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

3. Модифицирован численный метод решения задачи колебаний МТП в вязкой среде на основе метода Бубнова-Галеркина.

4. Разработан алгоритм, реализованный в комплексе программ «МАНОМЕТР», служащий для определения параметров затухания колебаний.

5. Исследование влияния геометрических характеристик трубчатых пружин и свойств демпфирующей жидкости на параметры затухания колебаний показало, что увеличение толщины стенки ведет к увеличению частоты затухающих колебаний и уменьшению коэффициента затухания, а увеличение отношения полуосей, центрального угла и радиуса продольной оси ведет к уменьшению частот затухающих колебаний и увеличению коэффициента затухания. Частоты затухающих колебаний зависят от геометрических характеристик трубчатых пружин. Коэффициент затухания определяется динамической вязкостью демпфирующей жидкости и в меньшей степени зависит от геометрических характеристик трубчатых пружин. При увеличении вязкости демпфирующей жидкости в 3-6 раза, время затухания колебаний уменьшается в 1,3-6,5 раз.

6. Экспериментальные исследования показали, что отклонение расчетных значений частот колебаний от экспериментальных значений во всем диапазоне не превышает 10% для анализатора вибрации и 3% для

вискозиметра, отклонение коэффициента затухания 15% и 8% соответственно.

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в рецензируемых журналах, рекомендованных ВАК РФ

1. Черенцов Д.А., Пирогов С.П., Дорофеев С.М. Математическая модель манометрической пружины в вязкой среде.//Вестник тюменского государственного университета №7 2014, с. 234.

2. Черенцов Д.А. Определение обобщенных сил сопротивления в уравнениях Лагранжа второго рода, описывающих колебания манометрической пружины при наличии сил вязкого трения.// Естественные и технические науки. 2014. №4 (72). С.12-15.

3. Черенцов Д.А., Пирогов С.П. Исследование затухающих колебаний манометрической трубчатой пружины в вязкой среде.//Известия высших учебных заведений. Нефть и газ №5 2014, с.66-72.

4. Черенцов Д.А., Пирогов С.П., Земенков Ю.Д., Чуба А.Ю. Решение задачи колебаний манометрической трубчатой пружины в вязкой среде.//Естественные и технические науки. 2014. №11-12. С.19-22.

5. Черенцов Д.А., Пирогов С.П., Земенков Ю.Д., Чуба А.Ю. Результаты расчета параметров затухания манометрических трубчатых пружин // Известия высших учебных заведений. Нефть и газ №1 2015, с.70-74.

6. Черенцов Д.А. Расчет параметров затухающих колебаний манометрических трубчатых пружин. / Д.А. Черенцов, С.П. Пирогов, С.М. Дорофеев, А.Ю. Чуба // Вестник Тюменского государственного университета. Серия "Физико-математическое моделирование, нефть, газ, энергетика", 2015. Том 1. N 1(1). с.136-145.

Статьи в других изданиях

7. Черенцов Д.А. Методы уменьшения влияния вибраций на точность измерения давления / Д.А. Черенцов, Е.А. Ищук, А.Ю. Чуба // Материалы всероссийской научно-практической конференции студентов, аспирантов и молодых ученых: Новые технологии - нефтегазовому региону/ Тюменский государственный нефтегазовый университет. - 2013. - .321-323.

8. Черенцов Д.А., Чуба А.Ю., Пирогов С.П. Модель колебаний манометрической трубчатой пружины в вязкой жидкости.//Материалы

всероссийской заочной научно-практической конференции: Проблемы функционирования систем транспорта/ Тюменский государственный нефтегазовый университет.-2013.-373-377.

9. Черенцов Д.А., Пирогов С.П. Определение параметров затухающих колебаний манометрических трубчатых пружин, при наличии сил вязкого сопротивления.//Научно-теоретический и практический журнал. №20 (99) 2014.

Ю.Черенцов Д.А., Пирогов С.П., Чуба А.Ю. Программный комплекс «МАНОМЕТР» для расчета параметров затухания колебаний манометрических трубчатых пружин.//Сборник статей по материалам ХХУП-ХХ1Х международной заочной научно-практической конференции: Научная дискуссия: вопросы технических наук./ «Международный центр науки и образования», 2014.-174с.

11.Черенцов Д.А. Оценка влияния геометрических характеристик и вязкости демпфирующей жидкости на параметры затухания колебаний манометрических трубчатых пружин / Черенцова С.А., Ищук Е.А., Пирогов С.П. // Материалы всероссийской научно-практической конференции «Механика и процессы управления». Тюмень 2015 г. Изд.: ТюмГНГУ.

12.Черенцов Д.А. Экспериментальное исследование параметров затухания колебаний манометрических трубчатых пружин работающих в вязкой среде с помощью анализатора вибрации АШ14. / Черенцова С.А., Ищук Е.А., Пирогов С.П. // Материалы всероссийской научно-практической конференции «Механика и процессы управления». Тюмень 2015 г. Изд.: ТюмГНГУ.

Свидетельства о регистрации программы для ЭВМ

13.Свидетельство об официальной регистрации программы ЭВМ 2014618861 РФ. Программный комплекс «ПРЗКМП» для расчета манометрических трубчатых пружин/ Черенцов Д.А. - 2014616833; Заявл. 15.07.14; Опубл. 1.09.14.

14.Свидетельство об официальной регистрации программы ЭВМ 2015615645 РФ. Программный комплекс для расчета манометрических пружин (МАНОМЕТР) / Черенцов Д.А., Пирогов С.П., Чуба А.Ю. - №2015612166; Заявл.: 25.03.15; Опубл.: 22.05.15.

Подписано в печать 13.10.2015. Формат 60x90 1/16.Усл.печ.л.1,5 Тираж 120 экз. Заказ № 322.

Библиотечно-издательский комплекс федерального государственного бюджетного образовательного учреждения высшего образования «Тюменский государственный нефтегазовый университет», 625000, Тюмень, ул. Володарского, 38.

Типография библиотечно-издательского комплекса 625039, Тюмень, ул. Киевская, 52.