Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Козобков, А. А. Электрическое моделирование вибраций трубопроводов

.pdf
Скачиваний:
7
Добавлен:
19.10.2023
Размер:
6.3 Mб
Скачать

§ 1. ЭЛЕКТРИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЫНУЖДЕННЫХ КОЛЕБАНИЙ ТРУБОПРОВОДА ОТ ВНЕШНЕГО ВОЗДЕЙСТВИЯ

На трубопровод в месте его крепления (рис. 90) могут дей­ ствовать силовые факторы (Q, М, Мкр). При этом трубопровод может совершать изгибные и крутильные колебания, характе­

ризуемые в месте крепления следующими величинами: у, 0, ср. Для расчета вынужденных колебании трубопровода, возникаю­ щих под воздействием внешней силы, к его модели необходимо

Рис.

90. Расчетная

Рис.

91. Схема

мо­

схема трубопрово­

дели

внешнего

воз­

да,

закрепленного

 

действия

 

 

хомутом

 

 

 

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

Электрическая модель колебаний трубопроводной системы выполняется в этом случае по рассмотренной выше методике. Она будет иметь вид, показанный на рис. 91, где / — электриче­ ская модель изгибных колебаний, II — электрическая модель крутильных колебании трубопровода. Если известны силовые факторы, действующие на трубопровод, то в соответствии с при­ нятой системой аналогий внешнее воздействие моделируется источниками тока i, а если известны геометрические факторы, то источниками напряжения и (см. рис. 91).

§2. ЭЛЕКТРИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЫНУЖДЕННЫХ КОЛЕБАНИЙ ТРУБОПРОВОДА ОТ ПУЛЬСИРУЮЩЕГО

ПОТОКА

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

Возмущающее воздействие в трубопроводных системах по­ является в местах установки вентилей, задвижек, кранов, в ме­ стах перехода с одного диаметра на другой, в поворотах трубо-

130

Рнс. 92. Схема модели внешнего воздействия с преобразованием

провода и т. д., т. е. там, где имеется акустическая неоднород­ ность или включение.

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

венной координаты.

Если давление и скорость потока задать

в виде напряжений

и токов, то, преобразовав их по определен­

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

Полная схема модели для этого случая, изображенная на рис. 92, содержит эле­ ктрическую модель изгибных колебаний участка трубопровода (/) и преобразо­ ватель (II), источники напряжения (и) и гока (/), соответствующие давлению и скорости потока.

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

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

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

Линеаризованные по И. Парному уравнения, описывающие с общепринятыми допущениями поток в трубопроводе, имеют вид

др_

М 0а

а0аО;

др

__ 1

dG

(303)

дха

д(

[С0а

дха

 

 

 

где р — давление;

G — объемный расход; М0а — акустическая масса; С0а — акустическая емкость;

131

л'а — пространственная координата;

аоп — коэффициент трения.

Решая уравнения системы (303), совместно получим

дЮ

дЮ

r a0jC0 а

dG

0.

( 3 0 4 )

дха

^ОаХа dt-

dt

 

 

Приближение уравнения (304) второго порядка методом пря­ мых Слободянского имеет вид

0,704дл"С0аЛ40а d'-G„

:[Gm_1- 2 G m + Gm+1]

 

6,75

X

 

dt2

 

 

 

 

X

d-Gm- 1

dn-Gm+I

a0a^0a

dG„,

 

(305)

dt"-

dt°-

dt

 

 

 

 

 

 

 

t-(m-i)

t-m

t(m+i)

 

 

 

Рис. 93. Схема модели пульсирующего потока о трубо­ проводе

Уравнение (305) моделируется электрической цепью, состоящей из звеньев R, L, С (рис. 93). Цепь, изображенная на рис. 93, опи­ сывается уравнением, аналогичным уравнению (305)

1,73СЭЬЭ d-im

[0л-X— Чл + 0л+1

2,75

X

\/[ d4m- |

dn-tm.и

■ Х А

dim

(306)

 

ДТ2

 

dx

 

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

 

р

>

 

G

t

1

 

 

 

и

m G = — i

т , = — т ш=

“э

 

 

 

 

1

X

mt

 

_

_0,637/И 0аД.т

_

0,637С0аД *

_

а 0аДл'

 

' Ч

- ---------7-----------

 

т с а =

---------р -----------

m R — —

^ --------

 

 

L,$

 

Оз

 

/\э

 

Сравнение уравнений (305) и (306) с учетом масштабных коэффициентов позволяет получить индикаторы подобия, опреде­ ляющие подобие процессов в объекте и в модели

тмато =

mRm0

тс тр

(307)

mpmt

тр

m<amt=

niQUlt

 

132

Электрическая модель пульсирующего потока II (рис. 94), собранная из элементов, величины которых удовлетворяют соот­ ношениям (307), возбуждается источником тока I, моделирую­ щим расход среды в начальном сечении [1]. Напряжение в узлах и токи в ветвях электрической модели пропорциональны соот­ ветственно давлению и расходу среды в моделируемом трубо­ проводе. Рассчитанные подобным образом давление и расход

•среды в местах акустических не-

Рпс. 94. Структурная схема

Рис. 95. Структурная схема

элек­

электрической модели пульси­

трической

модели вынужденных

рующего

потока:

колебаний

трубопроводной

си­

J—источник тока;

// —модель по­

 

стемы

 

тока

 

 

 

 

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

Напряжения и токи в точках хаь хач, .., ха „ электрической модели потока среды I преобразовываются в функциональных блоках II: Фь Ф2, .. ., Ф„ и передаются на электрическую модель колебаний трубопровода в точки (х, у, z ) ь (х, у, z )2, ... (х, у, z)„.

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

R(x, у, z, t) =Ар~(хй, t) +BG~ (ха, t),

(308)

где А и В — постоянные коэффициенты, характеризующие аку­ стическую неоднородность.

Таким образом, усилие складывается из двух составляющих. Первая обусловлена наличием в трубопроводе пульсационной составляющей давления р, вторая— наличием пульсационной составляющей скорости, или пропорциональной ей величине объемного расхода G^4-

Используя систему принятых соответствий для модели пуль­ сирующего потока и для модели колебаний трубопровода, за­ пишем выражение (308) в электрических величинах

**(•*> у, z, t)= au3{x)-\- Ыэ{х),

(309)

133

где

lR— ток, моделирующий возмущающее воздействие,

 

о, b — постоянные

коэффициенты,

которые

 

соответствуют

 

коэффициентам А и В уравнения (308).

 

 

 

Введем масштабные коэффициенты по постоянным величинам

 

 

тА —

А

 

 

В

 

 

 

 

 

 

 

 

а

•, тв = — .

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

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

электрические

величины

с учетом их масштабных коэффициентов, получим

 

 

 

 

R = Ар

'Ajnp

BG

mQ

 

 

 

 

(310)

 

 

 

 

 

т в mG

 

 

 

 

 

Сравнивая выражения

(308) и (310), видим,

что для их тож­

 

 

 

 

 

 

дественности

необходимо,

 

 

 

 

 

 

чтобы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,

 

о

 

1.

(311)

 

 

 

 

 

 

т Атр

т в та

 

 

 

 

 

 

Таким образом, при соб­

 

 

 

 

 

 

людении

индикаторов подо­

 

 

 

 

 

 

бия

(311)

выражение

(309)

 

 

 

 

 

 

связывает

ток iB,

который

 

 

 

 

 

 

моделирует

 

возмущающее

 

 

 

 

 

 

воздействие, и токи и напря­

 

 

 

 

 

 

жения,

моделирующие

рас­

 

 

 

 

 

 

ход

и

давление

потока в

Рис.

96. Схема действия

сил

на

трубо­

трубопроводе,

т. е. опреде­

ляет условие

передачи

воз­

 

провод

 

 

 

 

мущения

с

 

электрической

 

 

 

 

 

 

модели

пульсирующего

по­

тока на электрическую модель трубопровода.

 

 

 

 

 

 

Следует, однако, отметить, что

усилие R(x,

у,

z, t) является

в общем случае величиной

векторной и может

быть

представ­

лено в виде суммы трех векторов, являющихся проекциями век­ тора R на оси координат

R{x,

у, z, t) = iRx -\- jR tJ-\-kRz,

(312)

где Rx, Ry, Rz— проекции вектора иа оси х, у, z\

L, /, k — единич­

ные векторы (орты).

может быть записано и в другой форме

Выражение (312)

/ ? = | А*|(г cos a-j- j cos ^i-j- k cosy),

 

где a, (3 и у — углы

между вектором R и соответствующими

ортами, т. е. вектор усилия R раскладывается на три состав­ ляющие.

134

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

Проиллюстрируем это положение на примере.

Пусть дана трубопроводная система (рис. 96), на которую

действует усилие R(x, у, z, t) в точке С в плоскости хоу. Это уси­ лие можно разложить на составляющие

AH-V, у, z, t )= Щх j Ru= \R\(i cos а-\- J cos |U).

Рис. 97. Схема модели колебании трубо­

Рис. 98.

Схема

проводной системы,

изображенной на

модели

разло­

рис.

96

жения

вектора

 

 

силы R

на две

 

 

составляющие

Колебаниям трубопроводной системы (см. рис. 96) в плоско­ сти соответствует электрическая модель, изображенная на рис. 97. В сечении С электрической модели включены два источ­ ника тока, при этом току первого источника ix поставлена в со­

ответствие сила iRx, а току второго — сила iRy, или

= |Я | co sa — ,

гя = |7 ?| cos

.

Л

ПЩ

J

TtlQ

. v

Следовательно, ток iH, соответствующий силе R(x, у, z, t) должен быть разложен на две составляющие

 

Lr x

COS CL

 

iR

или

.

(313)

 

i#y

iR cos Э

 

Преобразование (313) может быть выполнено при помощи трансформатора с двумя вторичными обмотками (рис. 98).

135

Если отношения чисел витков

w.vl _

1

щ\ __

1

wxn

cos a ’

wyZ

cos р’

то при достаточно большом внутреннем сопротивлении источ­ ника тока iR, токи во вторичных обмотках трансформатора (рис. 98) окажутся равными

iR cos а и f'^cosp,

т. е. будет осуществляться преобразование (313).

Очевидно, что в самом общем случае число вторичных обмо­ ток трансформатора равно трем, так как сила R(x, у, z, t) пред-

i f l F

l i n

 

&'lr C0S<X MgCOSfi 0 Lr cos j

 

 

 

t

 

Рис. 99.

Схема модели

Рис. 100. Схема сумматора

разложения

вектора си­

 

лы R на три

составляю­

 

 

щие

 

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

Теперь более подробно остановимся на внутреннем содержа­ нии функциональных блоков (см. рис. 95).

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

На рис. 100 представлена принципиальная схема суммирую­ щего усилителя, который подключается в узле т электрической модели пульсирующего потока.

Одно звено электрической модели пульсирующего потока со­ стоит из емкости Сэ, индуктивностей Ьэ, L3 и активного сопро­ тивления R3 [1]. Напряжение в точке А электрической модели пульсирующего потока равно

i:}R;}*Т Ц),

где иэ — напряжение в узле т.

136

Предположим, что параметры лампы Л\ (рис. 100) таковы: внутреннее сопротивление Ru, 'Статический 'коэффициент усиле­ ния по напряжению рь

— параметры лампы Л2: внутреннее сопротивление R i2, ста­ тический коэффициент усиления по напряжению р2-

Определим выходное напряжение «Вых на выходе суммирую­ щего усилителя. Для этого сначала предположим, что напряже­ ние на сетке лампы Л2 равно нулю, т. е. цепь между сопротивле­ нием ДД 2 и разделительной емкостью Ср2 разорвана. Тогда пред-

Цашг

( к

|

 

 

[Viz

 

 

 

Рис. 101. Расчетная схе­

Рис. 102.

Рас­

ма сумматора при нвх2= 0

четная

схема

 

 

сумматора

при

ПдХ1 —0

ставим упрощенную эквивалентную схему суммирующего усилителя, как показано на рис. 101.

Анодный ток

__ г- __ Н-1 («вх! — ’й]Яп)

где «вх1 = Ыэ.

R a

+ R i i + R i i

 

 

 

 

 

 

или после преобразований

 

 

 

 

__________Р1ц вх1__________

 

 

Отсюда

R a +

/?<| + R i 2 (1 + pi)

 

 

 

 

 

 

______ ^эр!Ra______

Ra_

(314)

^ВЫХ1

+ R i \ + Ri2(L + pi)

R;2

R a

 

Теперь предположим, что напряжение на сетке лампы Л\ равно нулю. Обозначим

tt3-j- iaRa иijX2,

а напряжение на выходе суммирующего усилителя цВЫх2. Заменим лампу Л2 источником напряжения (рис. 102) и за­

пишем выражение для анодного тока ia2

_ рщ' + и’

32 ~ R n + R a + R i i

6

366

137

Раскрывая значение и', получим

_ (jM+ 1)Ц2ннх2~Г ‘ a%Ri2 (н + 1)

 

 

 

 

R i \ + R i l + R a

 

 

 

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

 

 

 

 

 

^

_ ____ (,ц| + 1)Ц2“их2___

 

 

 

3“

R i\ Ra + R ii

+

2)

 

 

Отсюда получим выходное напряжение

-И,.

 

(315)

('|ШХ2~

и их2,а2 (1+ Ц|) R a

Ra,

 

 

 

 

 

 

 

R i 1 Н" Ra + R ii (2 + щ )

 

 

Rii

 

Если коэффициенты деления делителей

(см.

рис. 100) k\ и к2,

то величина выходного напряжения суммирующего усилителя

 

 

WB-irt==if»H*L +

if2HHL.

 

(316)

 

 

 

к I

к2

 

 

С учетом выражений

(314) п (315)

выражение (316) примет вид

 

 

Ra

I

 

W

 

или

 

k\Ri2

b\R'tel fts(иэ+

 

 

 

 

 

 

 

=«э

IV koRio

k\R‘i<i ' +,-(ЙА^ Г >

(317)

 

Выражение (317) связывает выходное напряжение сумми­

рующего усилителя

с

током

и напряжением,

моделирующими

расход и давление в трубопроводе. Напряжение ивых пропорцио­

нально усилию ~R{x, у, z, t), воздействующему

на

трубопровод

в точке С(.V, у, z ) .

 

 

 

Напряжение иаых необходимо преобразовать в

ток, так как

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

трубопровода моде­

лируется не напряжением, а током.

 

 

 

Для преобразования напряжения ивых в ток гд, моделирую­

щий усилие R, на выходе суммирующего усилителя должен быть

включен усилитель напряжения (рис. 103)

с

последовательно

включенным сопротивлением R*. Сопротивление R* должно быть

во много раз больше сопротивления нагрузки

(входное сопро­

тивление модели в точке подключения источника гд). В этом слу­ чае усилитель напряжения совместно с сопротивлением R* будет являться источником тока.

Если коэффициент усиления усилителя напряжения равен к,

то величина тока на выходе составит

 

 

 

; _кивых

иэ

Rsk

Rak

Ra Ц2

kRa -

(318)

 

k<zR*Rj2

kiR*Rn + гэ .

R*

Rii

 

138

Сравнивая выражения

(318) и

(319), видим, что они будут

тождественны в том случае,

если

коэффициенты а

и b будут

удовлетворять условию

 

1

R^kR^ixо

 

а-

Rak

Ц2

1319)

R*Ri2

kn

*1.

R*Ri2^2

 

приведена

Полная схема функционального

преобразователя

на рис. 104.

 

 

 

 

 

 

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

систему в

местах

акустических

 

 

 

неоднородностей, может быть ре- и3+ь3кэ

 

 

алнзовано в- электрических моде­

0~—~

<

U.,0

у

лях при помощи функциональных

 

Z_

К

преобразователей (рис. 104). При

 

 

3

этом предполагалось, что вид

 

 

 

зависимости

(308),

связывающей

Рис.

103.

Схема, формирующая

параметры пульсирующего пото­

 

 

ТОК I,

ка в месте данной

акустической

 

 

 

неоднородности, и усилия, возникающего в ней, известен, т. е. известны значения коэффициентов А и В, входящих в эту за­ висимость, а также углы а, (3, у. Тогда для нахождения пара­ метров функционального преобразователя, т. е. значений коэффи­ циента усиления усилителя к, коэффициентов деления ku /е2 и т. д., предположим, что параметры электрических моделей колебаний трубопровода и пульсирующего потока известны [1], а значит и известны значения масштабных коэффициентов тр, tnG, mQ.

Из выражений

(311)

находим значения масштабных коэффи­

циентов по постоянным величинам а и b

 

 

 

 

 

 

тр

I1lh

 

 

 

 

 

 

 

тп

та

 

 

 

 

 

 

 

 

 

U.3+L3R3

 

 

~

i J h c g f l

 

 

0

1

-1 <

 

 

и-э0"

 

 

 

 

 

 

К

a ' f l i t

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

~

1

 

 

 

 

cosa

cosfi

cosy

 

Рис. 104. Схема функционального

преобра-

 

а затем и значения самих коэффициентов

 

 

 

 

а = А —-

Ь ^ В ^ - .

 

 

 

 

mQ

 

 

mQ

 

 

 

При фиксированных

параметрах

функционального

преобра­

зователя R*, Ri2, цг, Ra и параметра

RD электрической

модели

6*

139

Соседние файлы в папке книги из ГПНТБ