Например, Бобцов

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

ПРИБОРЫ ТОЧНОЙ МЕХАНИКИ

УДК 622.691.4

С. П. ПИРОГОВ, А. Ю. ЧУБА
РАСЧЕТ ЧАСТОТ СОБСТВЕННЫХ КОЛЕБАНИЙ МАНОМЕТРИЧЕСКИХ ТРУБЧАТЫХ ПРУЖИН
Представлен вывод уравнений движения манометрической трубчатой пружины. Пружина рассматривается как изогнутый стержень, совершающий колебания в плоскости кривизны центральной оси. Приведены зависимости частоты собственных колебаний трубчатых пружин от геометрических параметров.
Ключевые слова: собственные колебания, манометрическая трубчатая пружина.
Манометрические приборы часто работают в условиях вибраций, при этом их основные упругие чувствительные элементы — трубчатые пружины — совершают колебательное движение, что негативно сказывается на точности приборов.
Важной характеристикой виброустойчивости приборов является частота собственных колебаний трубчатых пружин, поэтому необходимо определить влияние их геометрических размеров на данный параметр.
Существуют запатентованные конструкции манометрических пружин, обладающих повышенной вибростойкостью. Все эти пружины имеют поперечное сечение, переменное по длине трубки. В настоящей статье предлагается метод расчета частот собственных колебаний, применимый и для пружин, имеющих переменное поперечное сечение.
Будем рассматривать трубку Бурдона как изогнутый стержень, совершающий колебания в плоскости кривизны центральной оси. Уравнения колебаний (в соответствии с принципом Даламбера) получены из равенств нулю сумм проекций на нормаль и касательную всех сил, приложенных к элементу пружины (с учетом силы инерции) [1] (рис. 1):

γ R

w M + ∂M d ϕ u ∂ϕ

N + ∂N d ϕ ∂ϕ Q + ∂Q d ϕ ∂ϕ
Рис. 1





N R



1 R

∂Q ∂ϕ

=

mi

(ϕ)

∂2 ∂t

w
2

;⎪⎫ ⎪

1 R

∂N ∂ϕ



Q R

=

mi

(ϕ)

∂2u ∂t 2

,

⎬ ⎪ ⎪⎭

Q M N
R

(1)

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 1

40 С. П. Пирогов, А. Ю. Чуба
где N — продольная сила; Q — поперечная сила; mi (ϕ) — масса единицы длины трубки
(масса поперечного сечения с координатой φ); w, u — соответственно радиальная и окружная составляющие перемещения центра тяжести поперечного сечения с координатой φ.
Система уравнений (1) для перемещений u и w имеет следующий вид:

∂ ∂ϕ

⎨⎧DR ⎩

⎛ ⎜ ⎝

∂u ∂ϕ

+

w ⎟⎞ ⎫⎬ ⎠⎭

+

∂ ∂ϕ

⎧⎪ ⎨

H

R

⎩⎪

⎛ ⎜⎜⎝

∂u ∂ϕ



∂2w ∂ϕ2

⎟⎠⎟⎞⎬⎭⎪⎫⎪

=

mi

(ϕ)

∂2u ∂t 2

⎫ ;⎪ ⎪⎪

−DR

⎛ ⎜ ⎝

∂u ∂ϕ

+

w

⎞ ⎟



+

∂2 ∂ϕ2

⎪⎧ ⎨

H

R

⎩⎪

⎛ ⎝⎜⎜

∂u ∂ϕ



∂2w ∂ϕ2

⎞⎟⎠⎟⎪⎬⎫⎭⎪

=

mi

(ϕ)

∂2w ∂t 2

;

⎬ ⎪ ⎪ ⎭⎪

(2)

( ) ( )DR =

ES (ϕ)
1− µ2 R2

,

HR

=

EJ (ϕ) KK (ϕ)
R4 1−µ2

,

mi (ϕ)

=

ρS

(ϕ) ,

где Е — модуль упругости материала трубки; S (ϕ) — площадь поперечного сечения трубки;
µ — коэффициент Пуассона; J (ϕ) — момент инерции сечения; KK (ϕ) — коэффициент
Кармана (определяется методом, предложенным в работе [2]); величины S (ϕ) , J (ϕ) , KK (ϕ)
зависят от угловой координаты φ сечения и его геометрических размеров. Система уравнений (2) решается при следующих граничных условиях: в сечении жест-
кого закрепления пружины (при φ=0) касательное и нормальное перемещения и угол поворота γ поперечного сечения трубки равны нулю, а на противоположном конце (при φ=γ) изгибающий момент (М), перерезывающие и растягивающие усилия обращаются в нуль.
Собственные изгибные колебания подчиняются гармоническому закону с частотой k, поэтому решение системы уравнений (2) можно представить в виде

u (ϕ,t ) = u (ϕ)sin (kt + β) , w(ϕ,t ) = w(ϕ)sin (kt + β) ,

(3)

где k — круговая частота колебаний. При решении системы (2) применялся приближенный метод Бубнова — Галеркина.
В соответствии с этим методом зададим искомые функции составляющих перемещений:

nn
∑ ∑u (ϕ) = a1u1 + a2u2 + ... + anun = aiui ; w(ϕ) = b1w1 + b2w2 + ... + bn wn = bj wj , i=1 j=1

(4)

где a1, a2, …, an, b1, b2, …, bn — неопределенные коэффициенты; u1, u2, un, w1, w2, wn — базисные функции переменной φ.
Подстановка выражений (4) в преобразованную методом Галеркина систему (2) приводит
к системе из 2n уравнений, при этом в силу граничных условий часть слагаемых равны нулю:

∫ ∫a1

γ

⎛ ⎜



DR

0⎝

du1 dϕ



HR

du1 dϕ

⎞ dui

⎟ ⎠



d

ϕ

+

...

+

an

γ

⎛ ⎜



DR

0⎝

dun dϕ



HR

dun dϕ

⎞ ⎟ ⎠

dui dϕ

dϕ +

∫ ∫+b1

γ 0

⎛ ⎝⎜⎜



DR

w1

+

HR

d 2w1 d ϕ2

⎞ ⎠⎟⎟

dui dϕ

γ⎛ d ϕ + ... + bn 0 ⎝⎜⎜ −DR wn

+

HR

d 2wn d ϕ2

⎞ ⎠⎟⎟

dui dϕ



=

γγ
∫ ∫= −a1k 2 m0u1u j d ϕ − ... − ank 2 m0unu j dϕ ; 00

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 1

Расчет частот собственных колебаний манометрических трубчатых пружин

41

∫ ∫γ ⎛

a1

0

⎜ ⎜⎝



DR

du1 dϕ

wj

+ HR

du1 dϕ

d2wj d ϕ2

⎞ γ⎛

⎟d ⎟⎠

ϕ

+

...

+

an

0

⎜ ⎝⎜

− DR

dun dϕ

wi

+ HR

dun dϕ

d2wj d ϕ2

⎞ ⎟dϕ + ⎟⎠

∫ ∫γ ⎛

+b1

0

⎜ ⎝⎜



DR

w1w

j



HR

d 2w1 d ϕ2

d2wj d ϕ2

⎞ γ⎛

⎟d ⎟⎠

ϕ

+

...

+

bn

0

⎜ ⎝⎜



DR

wn

w

j

− HR

d 2 wn d ϕ2

d2wj d ϕ2



⎟d ⎠⎟

ϕ

=

γγ
∫ ∫= −b1k 2 m0w1wj dϕ − ... − bnk 2 m0wn wi dϕ . 00
Выберем базисные функции ui и wj:
ui (ϕ) = ϕi , i = 1, ..., n; wj (ϕ) = ϕ j+1, j = 1, ..., n.

Легко проверить, что выбранные функции удовлетворяют главным граничным усло-
виям. Получим однородную систему алгебраических уравнений порядка 2n относительно
неизвестных a1, …, an, b1, …, bn. Данная система уравнений имеет ненулевое решение только в том случае, если определитель матрицы этой системы равен нулю. Запишем его
в виде

A1,1 ... B1,n ... ... ... = 0 ,
A2n,1 ... B2n,n

(5)

где

γ γγ
∫ ∫ ∫ ( )A1,1 = (−DR − HR )dϕ + k 2 m0ϕ2dϕ , B1,n = −DRϕn+1 + HRn (n + 1) ϕn−1 dϕ ; 0 00

γ
∫ ( )A2n,1 = −DRϕn+1 + HR n (n + 1) ϕn−1 dϕ , 0

∫( ) ∫γ γ
B2n,n = −DRϕ2n+2 − H R n2 ( n + 1)2 ϕ2n−2 d ϕ + k 2 m0ϕ2n+2dϕ .
00
Условие равенства нулю определителя (5) можно рассматривать как уравнение для определения частот колебаний k. Поскольку порядок определителя равен 2n, то уравнение будет иметь 2n корней, которые являются частотами собственных колебаний трубчатой пружины. Те значения k, при которых определитель равен нулю, являются круговыми частотами собственных колебаний [3].
Результаты численного эксперимента показали, что с увеличением количества базисных функций ui и wj частота колебаний стремится к некоторому предельному значению. Согласно расчетам для получения удовлетворительных результатов по первой собственной частоте колебаний достаточно вычислять по пять базисных функций.
На основе рассмотренного способа определения частот собственных колебаний составлены алгоритм и программа для ЭВМ [4], с помощью которой проведены расчеты частот собственных колебаний. По результатам расчетов построены графики зависимости частот собственных колебаний от геометрических параметров трубчатых пружин (рис. 2): радиуса изгиба

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 1

42 С. П. Пирогов, А. Ю. Чуба
трубки (R), толщины стенки трубки (h), радиуса трубки-заготовки (r), угла поворота (γ) и отношения полуосей эллиптической трубки-заготовки (a/b).
ν, Гц R=0,0415 м а=0,0104 м
180 b=0,0022 м γ=220°
160

140

120

100 0 0,5 1 1,5 h, м

ν, Гц 400
300

R=0,0415 м а=0,01 м b=0,002 м
h=0,001 м

200

100

0 ν, Гц
400
300 200
100

100 200 γ, …°
R=0,037 м γ=250° h=0,001 м r=0,008 м

0

246

8 a/b

Рис. 2
Для пружин с переменным по длине поперечным сечением установлено, что уменьшение толщины стенки трубки от закрепленного конца к свободному, а также уменьшение радиуса трубки-заготовки от закрепленного конца к свободному приводит к увеличению частоты собственных колебаний. Сравнение манометрических пружин разных конструкций показало, что наибольшей частотой собственных колебаний обладают манометрические пружины, сечения которых изменяются от восьмеркообразного (в закреплении) до плоскоовального (на свободном конце). Частоты собственных колебаний трубчатых пружин переменного по длине поперечного сечения на 20 — 40 % превышают частоты собственных колебаний пружин обычных конструкций.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 1

Метод повышения надежности гидродинамических подшипников скольжения

43

СПИСОК ЛИТЕРАТУРЫ

1. Бидерман В. Л. Теория механических колебаний. М.: Высш. школа, 1980. 408 с.

2. Дорофеев С. М., Пирогов С. П., Самакалев С. С. Приближенное решение задачи об изгибе манометрической пружины переменного сечения // Мегапаскаль: Сб. науч. тр. Тюмень: ТюмГНГУ, 2006. № 2. С. 46—48.

3. Чуба А. Ю., Пирогов С. П., Дорофеев С. М. Определение собственных частот колебаний манометрических трубчатых пружин // Изв. вузов. Нефть и газ. 2007. № 2. С. 70—74.

4. Свид. об официальной регистрации программы ЭВМ, 2007612005 РФ. Программный комплекс „ПКРМТП“ для расчета манометрических трубчатых пружин / А. Ю. Чуба, С. С. Самакалев, С. П. Пирогов. 2007611194; Заявл. 2.04.2007; Опубл. 17.05.2007.

Сергей Петрович Пирогов Александр Юрьевич Чуба

Сведения об авторах — канд. техн. наук, доцент; Тюменский государственный нефтегазовый уни-
верситет, кафедра теоретической и прикладной механики; E-mail: piro-gow@yandex.ru — канд. техн. наук, доцент; Тюменская государственная сельскохозяйственная академия, кафедра общетехнических дисциплин; E-mail: aleksandr-chuba@mail.ru

Рекомендована кафедрой теоретической и прикладной механики ТюмГНГУ

Поступила в редакцию 24.12.10 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 1