книги из ГПНТБ / Козобков, А. А. Электрическое моделирование вибраций трубопроводов
.pdf§ 1. ЭЛЕКТРИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЫНУЖДЕННЫХ КОЛЕБАНИЙ ТРУБОПРОВОДА ОТ ВНЕШНЕГО ВОЗДЕЙСТВИЯ
На трубопровод в месте его крепления (рис. 90) могут дей ствовать силовые факторы (Q, М, Мкр). При этом трубопровод может совершать изгибные и крутильные колебания, характе
ризуемые в месте крепления следующими величинами: у, 0, ср. Для расчета вынужденных колебании трубопровода, возникаю щих под воздействием внешней силы, к его модели необходимо
Рис. |
90. Расчетная |
Рис. |
91. Схема |
мо |
схема трубопрово |
дели |
внешнего |
воз |
|
да, |
закрепленного |
|
действия |
|
|
хомутом |
|
|
|
приложить электрическое воздействие, эквивалентное внешнему возмущению, действующему на трубопровод.
Электрическая модель колебаний трубопроводной системы выполняется в этом случае по рассмотренной выше методике. Она будет иметь вид, показанный на рис. 91, где / — электриче ская модель изгибных колебаний, II — электрическая модель крутильных колебании трубопровода. Если известны силовые факторы, действующие на трубопровод, то в соответствии с при нятой системой аналогий внешнее воздействие моделируется источниками тока i, а если известны геометрические факторы, то источниками напряжения и (см. рис. 91).
§2. ЭЛЕКТРИЧЕСКОЕ МОДЕЛИРОВАНИЕ ВЫНУЖДЕННЫХ КОЛЕБАНИЙ ТРУБОПРОВОДА ОТ ПУЛЬСИРУЮЩЕГО
ПОТОКА
Для расчета вынужденных колебаний трубопровода, возни кающих под действием пульсирующего потока, на электрическую модель необходимо подать возмущающее воздействие, эквива лентное воздействию, которое оказывает пульсирующий поток на трубопровод.
Возмущающее воздействие в трубопроводных системах по является в местах установки вентилей, задвижек, кранов, в ме стах перехода с одного диаметра на другой, в поворотах трубо-
130
провода и т. д., т. е. там, где имеется акустическая неоднород ность или включение.
Величина усилия, действующего на трубопровод в месте аку стической неоднородности, в общем случае является функцией давления и скорости потока, поэтому для определения усилия, действующего на трубопровод, необходимо знать давление и скорость потока, являющиеся функциями времени и пространст
венной координаты. |
Если давление и скорость потока задать |
в виде напряжений |
и токов, то, преобразовав их по определен |
ному закону, можно получить электрические величины, пропорциональные усилиям, действующим на трубопровод. Такое пре образование может быть выполнено при помощи специального преобразователя.
Полная схема модели для этого случая, изображенная на рис. 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 |