CC BY
23
УДК 621.311.22.019.3.001.18
МОДЕЛЬ УПРАВЛЯЕМОГО ВИБРОГАСИТЕЛЯ КОЛЕБАНИЙ ТРУБОПРОВОДОВ
КРАСНОВСКИЙ С.Я., ЛЕБЕДЕВ В.Л., инженеры, ШАПИН В.И., канд. техн. наук
Рассматриваются аналитическая модель и принцип управляемой виброзащиты стержневых конструкций на базе электромеханических исполнительных элементов, а также алгоритм построения автоматизированной системы виброзащиты. Приводится структурная схема и результаты физического эксперимента.
Ключевые слова: виброгаситель, частота колебаний, амплитудно-частотные характеристики, функции управления.
MODEL OF CONTROLLABLE PIPELINE VIBRATIONS DAMPER
S.Ya. KRASNOVSKIY, V.L. LEBEDEV, engineers, V.I. SHAPIN, Candidate of Engineering
The article studies an analytical model and principle of controllable vibroprotection in framing based on elecrtomechanical operating elements. The authors consider the algorithm of developing an automated vibroprotection system. The structural scheme and physical experimental results are given.
Key words: vibration damper, vibration frequency, amplitude-frequency characteristics, control function.
Известно [1, 2], что существующие в настоящее время методы оценки надежности трубопроводов построены без достаточного учета их колебательных процессов. Однако вибрационные нагрузки могут оказаться весьма существенным фактором исчерпания живучести трубопроводов, эксплуатируемых с невыявленными трещинами [3]. Особенно опасными могут быть непроектные ситуации, когда возникают резонансные или близкие к резонансным пространственные изгибно-крутильные колебания трубопровода. Основным источником вибрации является спектр колебаний теплоносителя при движении по трассам упругих трубопроводов с природоимпульсной причиной вибровоздействия параметрического и автоколебательного характера [4]. В качестве сопутствующего фактора выступают элементы опорно-подвесной системы, как способствующие при определенных условиях переходных процессов потере виброустойчивости конструкции в целом. Таким образом, снижение виброактивности трубопроводов является актуальной технической проблемой.
Известные [5-9] существующие пассивные методы виброзащиты недостаточно эффективны при изменении спектральных характеристик вибрации.
Ниже рассматривается практическая реализация автоматизированной системы виброгашения на примере поперечных колебаний стержневых конструкций и на основе магнитожидкостного демпфера.
Возможность управления эффективностью гашения вибрации виброгасителем вытекает из следующих теоретических предпосылок. Приведенная на рис. 1 схема условно изображает виброгаситель поперечных колебаний объекта, представляющего собой стержневую конструкцию в виде системы с распределенными параметрами, обладающей конечной длиной I, изгибной жесткостью ЕЛу, плотностью р и поперечным сечением Г. В сечении с координатой Хо приложена возбуждаю-
щая вибрацию сила Р(т), а в сечении с координатой Х1 установлен магнитожидкостный демпфер в качестве исполнительного элемента (ИЭ), предназначенного для гашения поперечных колебаний объекта.
Дифференциальное уравнение поперечных колебаний объекта имеет вид
д
дХ 2
EJ
dz
~дт'
дХ 2
kz
-pF
д2 z ~&ґГ
-h
дz
(1)
8(Х - Х1 ) = P ()8(Х - Х0).
Здесь, в дополнение к отмеченным механическим параметрам объекта, введены следующие обозначения: X - продольная координата объекта; т - время; г(Х,т) -функция вертикального перемещения объекта; Л - коэффициент, характеризующий затухание в автономном состоянии объекта; Лі - коэффициент, характеризующий затухание в контуре управления.
Силы сопротивления (диссипативные силы) объекта, при наличии которых происходит потеря энергии, определяются силами внешнего затухания, зависящими от сопротивления среды, и силами затухания, зависящими от вязкости материала (внутреннее затухание). Для последнего случая в первом приближении принимается гипотеза, что сила затухания пропорциональна первой степени скорости перемещения объекта в колебательном процессе. Соответственно управляемый параметр Лі определяется вязкостью магнитожидкостной среды (МЖС) электромеханического ИЭ. Здесь также принято, что сила затухания пропорциональна скорости перемещения объекта в колебательном процессе.
ИЭ
Рис. 1. Расчетная схема управляемого магнитожидкостного виброгасителя
В расчетной схеме и в уравнении (1) принято также: к - управляемый параметр, определяющий же-сткостные характеристики присоединенного магнитожидкостного ИЭ; 5(Х - X), I = 0,1 - функция Дирака, определяющая местоположение соответственно возмущающей силы и ИЭ.
Ниже выполняются преобразования с уравнением (1), определяется решение для функции перемещения и(Х, т) с доведением до амплитудно-частотной характеристики и осуществляется ее анализ в зависимости от параметров настройки ИЭ.
Для придания универсальности результатам решение и анализ удобно представить в безразмерном виде по геометрии и времени.
Введем безразмерные параметры:
Z (*. t )-
l
х=—■ tEJ
l ’ l2 V P—
(2)
После подстановки обозначений (2) и дополнительных преобразований исходное дифференциальное уравнение приобретает вид
д4г(х, I) д2г (х, t)
дх 4
hj2
hl2
д t2
dZ (х, t)
1 dZ (x, t)
yjEJpF d t
■sjEJpF d t
l3
k—Z (x,t) EJ
(3)
x8(x-x,)=—P(t)8(x-xo) .
В исходном уравнении (1) рассматривается гармоническая возмущающая сила Р(т) = Р0 cos ют, где Р0 - амплитуда возмущающей силы; ю - частота вынужденных колебаний.
В преобразованном уравнении (3) принято:
P (t) = P0cos ю l2 t - P0cos Qt,
v EJ
(4)
2 /
где 0 = ю/ - безразмерная частота возбуждения.
Решение уравнения (3) определяется в известном [5] виде произведения двух функций - перемещения и времени:
г (х о=! и (х)г- (t), (5)
I
где функция перемещений задается также в известном виде:
Uj (x) = Cjj sin XjX + Cj2 COSXjX +
+Cj 3 shX jX+Cj 4 ch X jX.
Здесь Xy - коэффициент, соответствующий j-й форме поперечных колебаний объекта (j = 1,2,... ю).
В результате прямых подстановок, а также с использованием свойства ортогональности балочных функций и свойств функции Дирака уравнение в частных производных (3) преобразуется к виду
Z ¡ +X4Z¡ +h
l2
Z,-
л/EJp F ¡
l4
\jEJp F ¡ EJ
U¡ (xi)-
(7)
-P(t)——U¡ (x0) . v !K¡EJ л
После введения дополнительных безразмерных обозначений
2^-h
v-k
l2
VEJP— l4
■ 2^i -h.
1Л/Цр— ’
l3
(8)
—■ f-P0-----------
EJ 0 KjEJ
получено обыкновенное дифференциальное уравнение вынужденных колебаний в виде
^1 + 2 [^ + (Х1)] ^1 +
/ 2\2 / ч] / N (9)
() +У^і (х,) ^ = Ю, (Хо).
Здесь параметры Е, и ^1 характеризуют диссипативные характеристики в комплексной системе объект-ИЭ виброгашения; V - жесткостные характеристики магнитожидкостного ИЭ; ґ - амплитуда возмущения; К) - известная интегральная структура, зависящая от формы колебаний объекта [5].
После группировки обозначенных слагаемых и введения обозначений
[Ч+^Ю ( хі )]=^і,
-п2
(x2) +vU¡ (xi)
равнение (9) окі ому виду
Z¡ + 2^¡Z¡ +n^Z¡ - fU¡ (x0)■
(10)
уравнение (9) окончательно преобразуется к стандартному виду
,2V -fíl ( ) (11)
+
Очевидно, что предельных своих значений коэффициенты, которые характеризуют диссипативные и же-сткостные свойства комплексной системы объект-ИЭ, достигают при размещении ИЭ в пучности колебаний соответствующей формы движения объекта, т.е. при достижении максимума передаточной функции, Щх-\) ~ 1. Действительно, в этом случае регулировочные параметры (10) достигают максимальных значений.
При многоформовой вибрации объекта решение об оптимальном размещении ИЭ решается после проведения предварительного спектрального анализа, определения энергоемких форм движения и принятия компромиссного решения по расположению устройства для регулируемого гашения колебаний.
Таким образом, используя возможности вариабельности упругожесткостных и диссипативных свойств магнитожидкостного преобразователя, представляется возможным варьировать соответствующие параметры комплексной системы - параметры ц, и щ обыкновенного дифференциального уравнения колебательного типа (11).
При этом имеются два безразмерных отстроенных регулируемых параметра:
, =Г
/4
ГрЕЗ
(12)
= Л-
2^57^
Очевидно, что эффект максимальных возможностей виброзащитной системы определяется неравенствами вида
Лі
г4ез
и(Хі)>>1 тіиі(Хі)>>т-
(13)
Действительно, выражение для амплитуды колебаний системы (амплитудно-частотная характеристика) приобретает вид 1
Л = Ю, (Хо)-Кд,
П
где Кд - коэффициент динамичности. Иначе, в развернутом виде,
ю (хо )Л
А =-----------П------
(14)
(15)
а2 ^ 2 2 + 4ц2 а2 4
п
Х і =-
(16)
Приняв стандартные обозначения в виде □ , ц
Чу ' пу ’
выражение (15) приобретет стандартный для анализа вид
1
А =
Ю1 (Хо )-2
_________і
2 г
Хі)
(17)
-4ЬІхЗ
Для выражения фазочастотной характеристики, соответственно:
) = агсід
2^4-___Ш
1-^ '
1 пі
(18)
Таким образом, принцип работы магнитожидкостного ИЭ в функции активного виброгасителя основан на управлении частотной расстройкой от режима вынужденных резонансных колебаний объекта посредством вариации жесткости связи ИЭ к и вязкостного трения иЭ Л1 (12), (16). Последнее достигается вариацией электромагнитного поля ИЭ в контуре отрицательной обратной связи.
В развернутом размерном виде выражение для амплитудно-частотной характеристики объекта при включении ИЭ (к #0, Л1 #0) приобретает вид
Р0—— и, (х0)
0 к,ез л 01
1 + к-
езг
-иі (Х1
ю2/4 Рр
ю / —
ез
> I Г* Л,,, ,, ■—, 1 + ^ и , (х'
^ЕЗрЯ |_ Л 1
ю2/4 Рр ю / --------
Ей
1 + к
-Чти, (х1 езг4 л 1
Г° 1 + к------ри, (х1
1 езГ р 1
(19)
Соответствующая функция управления может быть заложена в блок регистрации и управления для реализация частотнозависимого управления работой виброгасителя.
На рис. 2 приведена блок-схема физического экспериментального стенда, предназначенного для исследования демпфирующих свойств виброгасителя на основе магнитожидкостного исполнительного элемента.
После предварительного усиления (8) и фильтрации сигнала помех (9) с помощью регистратора (10) фиксируется переходный процесс или амплитудночастотная характеристика объекта посредством вариации частоты вынуждающей силы генератора в окрестности резонансного состояния объекта. При этом сканирование частоты вынужденных колебаний осуществляется как в сторону увеличения диапазона частот, так и в сторону уменьшения для оценки линейности или нелинейности механической колебательной системы с оценкой зоны динамической неустойчивости. Контур управления механическими колебаниями (III), предназначенный для управляемой виброзащиты, включает магнитожидкостный исполнительный элемент (ИЭ) (11), блок питания (12) и приборную исследовательскую составляющую в виде веберметра (13), амперметра (14) и вольтметра (15) электроцепи контура управления.
Автоматизация эксперимента достигается посредством компьютерной связи (16), снабженной соответственно аналого-цифровым (АЦП) и цифроаналоговыми (ЦАП) преобразователями, звуковой картой в качестве задатчика частоты сканирования и известными алгоритмами распознавания резонансного состояния объекта, например, по смене знака производной амплитудно-частотной или фазочастотной характеристик с одновременным расчетом добротности системы (оценка диссипативных характеристик) по ширине ам-плитудно-чатотной характеристики, фиксированной резонансной частоте и амплитуде колебаний.
3
А
Рис. 2. Физический экспериментальный стенд: I - контур возбуждения механических колебаний объекта; II - контур регистрации механических колебаний объекта; III - контур управления механическими колебаниями; 1 - стальная балка прямоугольного сечения; 2 -монолитная опора с шарнирами; 3 - генератор звуковых колебаний; 4 - блок подмагничивания; 5 - обмотка постоянного тока; 6 - звуковая катушка; 7 - виброприемник ускорения; 8 - усилитель; 9 - фильтр сигнала от помех; 10 - регистратор; 11 - магнитожидкостный исполнительный элемент; 12 - блок питания; 13 - веберметр; 14 - амперметр; 15 - вольтметр; 16 - компьютер
Система автоматизации виброзащиты срабатывает начиная с порогового значения амплитуды колебаний объекта в следящем режиме при достижении последней нормативных значений, заложенных в базу данных по конкретному объекту. При этом сигнал с предварительного усилителя (8) подается непосредственно в компьютер через АЦП с выполнением известной цифровой операции фильтрации по методу Фурье-анализа. Управляющий сигнал в цепи контура отрицательной обратной связи определяется цепочкой: изменение тока приводит к изменению магнитного поля магнитожидкостного преобразователя в зазоре с изменением его вязкоупругих характеристик.
При отработке командного аналогового сигнала на выходе из компьютера включается дополнительный предварительный усилитель контура управления по пороговому значению тока. Изменение тока приводит к изменению магнитного потока магнитожидкостного преобразователя в зазоре, формируя тем самым включение или изменение дополняющих объект искусственных вязкоупругих характеристик ИЭ. Процесс авторегулирования параметров заканчивается при достижении амплитуды колебаний объекта нормативных значений.
При испытаниях магнитожидкостного преобразователя использовались 2 вида магнитной жидкости: МК-32 и МГМ-80. Последняя оказалась более эффективной. Ее параметры: МЖ МГМ-80; тип МК-2-40; вязкость 0,223 Пас; объемная концентрация 14,1 %; намагниченность насыщения 57 кА/м; жидкость-носитель ПЭС-4; вязкость при 00С - 0,0858 Пас, 200С - 0,0455 Пас, 300С -0,0336 Пас, 400С - 0,0259 Пас; теплоемкость
0,445 ккал/кг-°С; температура кипения 185 0С; диапазон рабочих температур ^ < t < 105 0С.
На рис. 3, 4, 5 приведены результаты экспериментов на физическом экспериментальном стенде в виде амплитудно-частотных характеристик объекта (АЧХ). Они соответствуют одной физической модели с последовательным возбуждением двух первых форм поперечных колебаний балки: первая - симметричная с узлами колебаний на опорах; вторая - кососимметричная с узловыми точками колебаний как на опорах, так и в центре пролета.
Если АЧХ, приведенные на рис. 3, определяют исходное динамическое состояние объекта, то АЧХ, отвечающие рис. 4, соответствуют режиму гашения при расположении ИЭ в пучности первой формы колебаний, что соответствует узлу второй формы движения. Естественно, что при этом колебания второй формы полностью не гасятся. Дискретность управления от максимума исходного состояния до минимума в режиме виброзащиты обозначена знаком Е.
При переводе ИЭ в пучность второй формы колебаний (рис. 5) достигается виброзащита во всем указанном спектре резонансных частот колебаний объекта.
Таким образом, предложенная модель управляемого виброгашения с использованием магнитожидкостных исполнительных элементов определяет инновационность ее практического использования для гашения поперечных колебаний стержневых конструкций, в том числе, применительно к трубопроводам ТЭС и АЭС.
А мм/с2 10-2
Ю Гц
Рис. 3. Исходная амплитудно-частотная характеристика (АЧХ) автономного объекта: 1 кососимметричная вторая форма движения
симметричная первая форма колебаний; 2
А мм/02*10-2 200
195
190
185
180
175
170
165
160
155
150
145
140
135
130
.....
^ 7 1
Е2
у
7
=£—
10 15 20 25 30 35 40 45 50 55 60 65 70 75
ю Гц
Рис. 4. Интегральная характеристика эффекта виброзащиты по первой форме движения
Рис. 5. Гашение обеих форм поперечных колебаний объекта с использованием прежней структуры ИЭ при увеличении тока управления (ИЭ размещен в пучности колебаний второй формы)
Список литературы 5.
1. ПНАЭ Г-7-002-86. Нормы расчета на прочность обору- 6.
дования и трубопроводов атомных энергетических установок: Утв. Госатомэнергонадзор СССР. - М.: Энерго- 7.
атомиздат, 1989.
2. РД Э0-0185-00. Методика оценки технического состояния
и остаточного ресурса трубопроводов энергоблоков АЭС. - 8.
М., 1999.
3. Живучесть паропроводов стареющих тепловых электро- 9.
станций / Ю.Л. Израилев, А.Я. Копсов, Ф.А. Хромченко и
др. - М.: Изд-во «Торус пресс», 2002.
4. Красновский С.Я., Шапин В.И. Моделирование живучести трубопроводов при стационарном вибровозмущении. Вибрационные машины и технологии. В 2 ч. Ч.1: Сб. науч. тр. - Курск: Курск. гос. техн. ун-т, 2005. - С. 124-130.
Красновский Сергей Ярославович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», инженер кафедры теоретической и прикладной механики, телефон (4932) 26-97-12, [email protected]
Лебедев Владимир Леонидович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», инженер кафедры теоретической и прикладной механики, телефон (4932) 26-97-12, [email protected]
Шапин Вадим Иванович,
ГОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина», кандидат технических наук, профессор кафедры теоретической и прикладной механики, [email protected]
Тимошенко С.П., Янг Д.Х., Уивер У. Колебания в инженерном деле. - М.: Машиностроение, 1985.
А.С. СССР №1689693. Вязкоупругий демпфер / Д.Ю. Павлов, В.В. Костарёв, В.А. Ветошкин и др. 1991.
Патент РФ №1689693. Кудряков Ю.Б., Крылов Н.И., Шу-рыгин М.Н. и др. Магнитореологический виброгаситель. 1998.
Патент РФ №2088349. Качак В.В., Крылов Н.И., Кудряков Ю.Б. и др. Гаситель Колебаний. 2002.
Коварда В.В. Затухание колебаний системы с магнитожидкостным инертным элементом. Вибрационные машины и технологии: В 2 ч. Ч.2: Сб. науч. тр. - Курск: Курск. гос. техн. ун-т, 2005. - С. 161-162.
