Файл: Козобков, А. А. Электрическое моделирование вибраций трубопроводов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 24.10.2024

Просмотров: 106

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Из выражения (172)

я = - ^ ( е 25' -

1) ~ ± (2 8 + 2 W ) -

( 1 7 3 )

Считая, что времени / = 1/6,

за которое амплитуда свободных

колебании упадет в е раз, достаточно для определения величины о, перепишем выражение (173) в виде

а = -

165

(174>

 

А*

 

Теперь определим погрешность, с которой функция Ле_6( аппрок-

симирует функцию -----

2

с учетом выражения (174) для

любого времени t. Для этого определим значение погрешности

------ л *

- 1 = е~»1/'1Ч-4«— 1.

(175)

2/ /

 

 

Значение экстремума погрешности (175) найдем из условия

 

 

^ - = 0,

 

(176)

 

 

dt

 

 

откуда

t;

1

 

 

45

 

 

 

 

 

 

Подставляя выражение (176) в выражение (175), получим

 

min max [v]= / 2

0, 11.

(177)

 

 

АП

 

 

Учитывая, что

— < 0, перепишем выражение (177)

в виде

 

dt2

 

 

 

Vmax — 0,1 1.

Таким образом, мы подобрали такое значение б, что с точностью до 0 , 1 1 можно записать решение уравнения (158) в виде

аА*

(178)

= А е 16

У

Sill £*(%()*+ <?<>)•

Можно показать, что полученному решению удовлетворяет, в ча­ стности, уравнение

д5у

dbj

 

(179).

E J ^ + 2Ъ dx^di

h> dr-

0 ,

которое представляет, как уже выше указывалось,

силу трения

в форме, наиболее удобной для моделирования.

 

70


Найдем решение уравнения (179), воспользовавшись, как п ранее, формой, предложенной Фурье

со

У = 2 у*'

 

S = 1

 

где

ys — Xs{x)Ts(t).

(180)

и в качестве граничных условий, как н в случае (159), выберем шарнирные.

Подставим решение (180) в уравнение (179) и опуская индекс s, получим

EJ X х4+ 2ь х ' VT + IV*: f = 0.

(181)

Как известно, для шарнирных граничных условий

 

 

 

Х = A sin kx,

 

(182)

где Я — постоянное число.

 

 

 

 

Подставляя выражение (182)

в уравнение (181), получим

 

но

 

но

 

=

(183)

 

 

 

 

 

Тогда из уравнения (183)

может быть получена

временная

•функция

 

 

 

 

 

 

 

Т = e~si sin (шм0 >‘-|-ср0).

(184)

Подставляя временную

функцию

в

уравнение (183), получим

уравнение для определения шмо и б:

 

 

 

Е Ш

2

Л—- ( 8 )— cujj 1 sin wM(/-

 

 

 

Ho

 

ио

 

 

J

 

26 ----ш„о “Ь(—2шм0о)

cos со 0t = 0,

 

 

Но

 

 

 

 

 

откуда

 

 

 

 

 

 

0

 

 

2ЪЬ — ^ л / ———

 

jm

 

' но

 

но

Но

 

S

 

 

 

 

 

26

Но

 

*-

 

 

 

 

, 6

 

 

 

 

2

 

Но

 

Тогда, подставляя в выражение (184) сомо и б, получим


или, учитывая выражения (180) и (182) можем написать

/И£/

-<Po

 

У— А е 1 40 sin kx sin

W)

( 1 8 5 )

 

 

 

Теперь, сравнивая выражения (185) и (178), можно заметить, что если начальные амплитуды свободных колебаний равны

то

А = А,

 

 

А2 _

кЫ

А*

Й 4 _ _

Н-о

а 16 ~

2 ц 0

16 ’

откуда

АЫ

 

 

 

 

 

8

Таким образом, мы связали коэффициент трения а, сосредо­ точенного в опоре, с коэффициентом распределенного трения Ь.

§ 7. ЭЛЕКТРИЧЕСКОЕ МОДЕЛИРОВАНИЕ ИЗГИБНЫХ КОЛЕБАНИИ ВКЛЮЧЕНИИ И НЕОДНОРОДНОСТЕЙ

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

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

смотрим электрическое моделирование

колебаний включений

и неоднородностей.

правила, позволяющие

Определим основные соотношения и

моделировать колебания трубопровода совместно с включением или неоднородностью, которые могут быть классифицированы следующим образом:

1 ) участок трубопровода с переменной массой и жесткостью (переход с одного диаметра на другой);

2 ) участок трубопровода с включенной сосредоточенной мас­ сой (задвижка, фланец, вентиль);

3) участок трубопровода с сосредоточенной жесткостью (про­ межуточная опора).

Рассмотрим эти участки.

1.Участок трубопровода

спеременной массой и жесткостью

Здесь возможны два случая: резкий (рис. 33) и плавный (рис. 34) переходы с одного диаметра на другой.

72


Для случая резкого перехода с одного диаметра на другой запишем граничные условия для сечения А

«/(+)=*/(-); ”6(+)—®(-); Q(+)=Q(->; а4(+)= м (_ ),

( i8 6 j

где величины с индексом ( + ) относятся к границе сечения А слева, с индексами (—) относятся к границе сечения А справа.

Учитывая принятую систему соответствий (91), перепишем соотношения (186)

«( + ) — “ ( - ) ! '?(Э+ ) = ''Р(Э_ ); г’( + >о,5 = *(—)о,5; /( + ) = /< - ) ,

где (<+)о,5 , /(—)о, 5

— токи, протекающие в верхних обмотках транс­

форматоров электрических моделей изгибных колебаний.

!/(+)

,,(+)

М(-) В (-)Р -ог EJ2

 

В(+)

Мм

—— 1

 

 

 

 

Ji-otEJj

 

 

Рис.

33.

Параметры,

Рис. 34. Параметры, определяю­

определяющие изгибные

щие изгибные колебания трубо­

колебания

трубопровода

провода с плавным изменением

с дискретным изменением

диаметра

 

диаметра

 

На рис. 35 показана электрическая модель колебаний трубо­ провода, изображенного на рис. 33. Здесь /—электрическая мо­ дель колебаний трубы с параметрами щн и EJь I I — электриче­ ская модель колебаний трубы с параметрами р0 2 и EJ2. Действи­ тельно, токи и напряжения в сечении А модели удовлетворяют соотношениям (186), следовательно, электрическая модель, изо­ браженная на рис. 35, моделирует изгибные колебания трубопро­ вода, изображенного на рис. 33. Параметры и число звеньев модели колебаний однородных участков трубопровода выби­ раются по рассмотренной выше методике.

Участок трубопровода с плавно меняющимися параметрами (конический переход) (см. рис. 34) может быть сведен к участку со ступенчато изменяющимися параметрами (рис. 36) так, что параметры трубы на одном участке, образующем ступень, остаются постоянными. Тогда колебания участка трубопровода с плавно изменяющимися параметрами будут моделироваться как колебания трубопровода с несколькими ступенчатыми пере­ ходами с диаметра на диаметр.

Рассмотрим замену трубопровода с плавно изменяющимися параметрами трубопроводом со ступенчато изменяющимися параметрами. Для определения числа и длины отрезков разбие­ ния найдем зависимость массы единицы длины трубопровода

73


и жесткости в функции длины, приняв за начало отсчета левоекрайнее сечение конуса х0. Масса элемента трубопровода, огра­ ниченного сечениями а п Ь, будет равна

 

 

 

д/И =

дх [D2 (х) -

сЕ (х)],

(187)

где g — плотность материала, D(x)

и

d(x) — внешний и

внут­

ренний диаметры на середине элемента Ах.

 

 

Из выражения

(187)

следует

 

 

 

 

 

 

 

_~ ~ = 0 ~~ [D1(х)

сС~(х)]

 

 

 

 

 

Ах

 

4

 

 

 

 

 

 

И

 

 

 

 

 

 

 

H+)o,s | H~W

 

 

 

 

 

 

I

T

| f

 

 

 

 

 

 

 

 

 

|Фн

 

 

 

 

 

 

 

<•(+)

I сн

 

 

 

 

 

 

Рис. 35. Схема мо­

 

Рис. 36. Замена плавного

 

дели

пзгнбных ко­

 

изменения

диаметра тру­

 

лебаний

трубопро­

 

бопровода

ступенчатым

 

вода

с дискретным

 

 

 

 

 

изменением

диа­

 

 

 

 

 

 

метра

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

/

\

1*

А.'VI

(7-VI

ЯО

г 1,-1 : \

|■>/ ■|

( 18 8 ?

 

M x ) = I i m — = — = - f [ 7 > H - r / - ( x i ] .

 

 

 

Д.1 - - .0

А х

d x

4

 

 

 

Жесткость в сечении с координатой х

 

 

 

 

 

 

EJ{x) = EnD4x) jl

d ( x )

 

(189)

 

 

 

v

'

64

l

LD(x)

 

Зависимость внутреннего диаметра конуса от координаты при

угле конуса а примет вид

 

 

 

 

 

 

 

 

d {x )= d 1-\-{x —x0)sin -|- .

(190)

Зависимость

внешнего диаметра конуса от координаты

 

 

 

D (х) =

Dx-j- (х — х0) sin

 

(191)

Тогда, подставив выражения (190)

и (191) в соотношения

(188)

74