Методика расчёта колебаний трубопровода с демпфирующими опорами из материала МР Текст научной статьи по специальности «Технологии материалов»

УДК 629.7.06

МЕТОДИКА РАСЧЁТА КОЛЕБАНИЙ ТРУБОПРОВОДА С ДЕМПФИРУЮЩИМИ ОПОРАМИ ИЗ МАТЕРИАЛА МР

© 2014 С. А. Безбородов1, А. М. Уланов2

1ОАО «Кузнецов», г. Самара 2Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет)

Разработана методика расчёта колебаний трубопровода произвольной формы с демпфирующими опорами из материала МР при помощи конечно-элементного пакета программ АК8У8. Методика учитывает нелинейность характеристик материала МР, позволяет определять амплитуды колебаний трубопровода и напряжения в нём. Результаты расчётов подтверждены экспериментом с использованием измерительной системы АЯАМК. Методика может быть применена для работ по доводке трубопроводов.

Трубопровод, колебания, демпфирование, материал МР, ЛЫБУБ, напряжения, метод конечных элементов.

Вибрация трубопроводов, установленных на двигателе, является фактором, определяющим долговечность их работы. Для уменьшения вибраций трубопроводов широко применяются системы конструкционного демпфирования (СКД), в частности система сухого трения в виде упру-годемпфирующей опоры (УДО) с материалом МР.

Материал МР является нелинейным материалом, то есть его жёсткость и коэффициент рассеивания энергии зависят от амплитуды сжатия и предварительной статической деформации. Он изготавливается методом холодного прессования проволочной спирали, обладает высоким демпфированием и высокой стойкостью к воздействию окружающей среды. Его свойства характеризуются следующей диаграммой (рис.1), где Р - нагрузка, х -амплитуда перемещения при статическом нагружении.

Важными количественными характеристиками, определяющими жёсткост-ные и демпфирующие свойства МР, являются:

- среднециклическая жёсткость Сср:

Сср=Ра/А,

где Ра - наибольшее значение статической нагрузки, А - соответствующее ей наибольшее значение перемещения;

У =

- коэффициент рассеивания энергии

К

П '

где К - площадь петли гистерезиса, равная количеству рассеиваемой энергии за один цикл деформации, П - условная потенциальная энергия, определяемая площадью треугольника ОАВ на рис. 1 [1].

Рис. 1. К определению параметров упруго-гистерезисной петли СКД

Существующие методики расчёта колебаний трубопроводов не учитывают сложность формы трубопровода, нелинейность свойств МР, демпфирование в опорах. В них рассматривается только относительное распределение перемещений

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

В данной работе статические и динамические характеристики МР предполагаются одинаковыми, что подтверждается сходством с динамическим экспериментом динамических характеристик, рассчитанных на основе статического эксперимента [2, 3, 4].

Для отработки методики расчёта в качестве модели трубопровода была выбрана трубка из титанового сплава ПТ-7М

с внешним диаметром 6,2 мм, внутренним диаметром 4,0 мм. Трубка имеет достаточно сложную форму, чтобы полученные для неё результаты могли быть применимы для трубопроводов произвольной формы. Оба конца трубки имеют штуцерные опоры, моделирующие связь трубопровода с агрегатами двигателя. В средней части трубки установлена опора шириной 11 мм, в которой могут устанавливаться сменные прокладки из материала МР с с различной предварительной статической деформацией Опоры жёстко прикреплены к швеллеру, который крепится на вибростенд (рис. 2).

Средняя опора с МР

Хомут для установки вибродатчика

Рис. 2. Исследуемая модель трубопровода

Построение расчётной модели для исследования колебаний трубопровода проводилось в пакете А№У8. Для задания геометрии трубки использовался упругий прямой двухузловой элемент Р1ре16. Опоры моделировались конечным элементом СотЫп14, точечная масса задавалась типом элементов Ма8821.

Параметры материала МР (жёсткость С и коэффициент рассеивания энергии у) рассчитывались на основе

экспериментальных данных, полученных при статическом эксперименте с тонкими пластинами из материала МР с различной толщиной Н (от 1,15 до 3,45 мм), плотностью р (от 1,115 до 2,028 г/см ), относительной предварительной

статической деформацией е£ = Q / Н (от 0,087 до 0,268), диаметром проволоки dw (от 0,09 до 0,12 мм) по следующим формулам [5]:

С = 14(1 + 363еЛ'5)(1 + 727е3/5) dlw2 Н"°'33(е£ /0,262)°'°63Н£ (р /0,2)1Д+0'01Я/^

у = 1,71(1 - 7,23Н013(еЛ-0,048)) (1 - 0,1(10,5е£ -1)');

d

0,138

'Л \0,057Н

„ Нидэ (0,095/■ (0,2/р)0>37+0-015Н/^

(1)

(2)

где р

т

ык

Ж г02 Нр 8

из МР к массе сделанного из стали эле- относительная плот- мента такого же объёма.

ность, равная отношению массы элемента

В эксперименте эти зависимости были получены для диапазона амплитуд относительной деформации еЛ = Л / H от 0,043 до 0,146 (более высокие значения деформации приводят к быстрому разрушению материала МР) [4].

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

По формулам (1) и (2) для определения С и у, выведенным по результатам статических испытаний тонких пластин из МР, определялись первоначальные значения жёсткости и коэффициента рассеивания энергии С(0) и у(0) при принятой начальной амплитуде сжатия МР по двум направлениям, равной Л(0) = 0,1 мм.

Далее проводился гармонический расчёт в АйБуБ, определялись амплитуды колебания трубопровода по двум направлениям ЛдЕ)РТ, Л(Г)Р в месте установки опоры с МР. По этим значениям рассчитывались уточнённые значения С(1) и у(1). Затем снова производился расчёт в

ANSYS, определялись значения A

(2) ВЕРТ

f(2)

, C(2) и yЭтот цикл повторялся до

(2)

тех пор, пока с достаточной точностью не выполнялись условия:

Л(К ) = Л (N-!)

ЛВЕРТ пВЕРТ м

Л (N) = Л(К-1)

ЛГОР ЛГОР ■

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

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

Характеристики штуцерных опор задавались через жёсткость на перемещение, жёсткость на поворот и коэффициент рассеивания энергии. Жёсткость штуцерной опоры на перемещение определялась экспериментально путём приложения силы к трубопроводу в месте его присоединения к штуцерной опоре. Получено, что жёсткость штуцерной опоры составляет С0 = 1111 Н/мм. Жёсткость штуцерной опоры на поворот принималась бесконечно большой величиной. Коэффициент рассеивания энергии ¥ оценивался по формуле

Y » 2Рц = 6,28/46 = 0,14,

где h =46 - коэффициент передачи вибрации на резонансе, полученный экспериментально.

В конечно-элементном пакете ANSYS была разработана программа для модального и гармонического расчёта трубопровода с опорой из МР. Поскольку применяемые конечные элементы в ANSYS используют вязкое трение, для задания коэффициента рассеивания энергии Y осуществлялся переход к эквивалентному коэффициенту демпфирования по площади петли гистерезиса по формуле

C„ =

YC 2жа>

где С - жёсткость, (О - резонансная частота.

Для экспериментального исследования колебаний трубопровода использовался вибрационный стенд ВДС-1500. Измерения колебаний трубопровода осуществлялись при помощи бесконтактной системы АКАМК. Она представляет собой две высокоскоростные камеры, установленные на треножнике, которые одновременно фотографируют объект. В процессе загрузки последовательности изображений определяется 3Б-форма объекта и деформации объекта. Полученная кар-

ЮР

тина деформаций может быть пересчитана в напряжения.

Расчёт собственных частот и амплитуд колебаний трубопровода в районе установки опоры с МР производился при установке различных демпфирующих

Из табл. 2, 3 видно, что погрешность в определении перемещений трубопровода в районе опоры с МР колеблется от 4 до 9 %. Погрешность по частоте составляет 6,5%. Расхождение расчёта с экспериментом по частоте и вертикальным перемещениям менее 10%, что достаточно для практических целей, так как разброс параметров МР из-за технологии изготовления также составляет порядка 10%.

Рассчитываемые значения перемещений трубопровода в районе опоры с МР позволяют уточнять значения амплитуды деформации материала МР и, соответ-

пластин из МР. Характеристики образцов приведены в табл. 1. Сравнение результатов расчётов и экспериментов для образца 1 приведены в табл. 2-3, амплитуды колебаний трубопровода в районе установки опоры с МР показаны на рис. 3-4.

ственно, жёсткости и коэффициента рассеивания энергии опоры с МР. Таким образом, можно более точно оценивать вибрационное состояние трубопровода и определять оптимальные характеристики материала МР (плотность, предварительное статическое сжатие и т.д.), применяемого в данных опорах.

Программа позволяет рассчитывать напряжения в трубопроводе на резонансных частотах. В качестве примера на рис. 5 приведена картина распределения эквивалентных напряжений в трубопроводе на резонансной частоте 207 Гц.

Таблица 1. Характеристики образцов из материала МР

№ образца Толщина, мм Свободный внешний диаметр системы, мм Диаметр отверстия в опоре, мм Диаметр трубки, мм Относительная статическая деформация е е Относительная плотность МР р

1 2,1 10,4 9,8 6,2 0,14 0.25

2 2,1 10,4 9,8 6,2 0,14 0.15

3 3,3 12,8 11,8 6,2 0,15 0.19

Таблица 2. Собственные частоты колебаний трубопровода

Амплитуды входного ускорения стола вибростенда, "1, м/с 30 50

Собственная частота, измеренная системой АЯАМШ, Гц 195,1 193,5

Собственная частота, рассчитанная в А№У8, Гц 206,6

Погрешность, % 6,5

Амплитуда колебания опоры по вертикали, мм 0,044 0,071

Таблица 3. Перемещения трубопровода по осям (в районе средней опоры с МР)

Амплитуды входного ускорения стола вибростенда, "Ш, м/с2 30 50

Резонансная частота, Гц 206,6 206,6

Амплитуда колебаний трубы около опоры, мм по вертикали /по гори- по вертикали /по горизон-

зонтали тали

измеренная системой АЯАМШ 0,052/0,04 0,074/0,065

рассчитанная в А№У8 в 1 прибл. 0,05/0,037 0,081/0,066

рассчитанная в А№У8 во 2 прибл. 0,05/0,037 0,08/0,059

Погрешность, % 4/7,5 8/9

РОБТ26

АМРЬ1ТиБЕ

иу126 иг12о

(х10**-5) 5. 5

А, Ш

и У

и

£, Нг

1;гуЬка р1ре1б дагш 1шг,30

Рис. 3. Зависимость амплитуды от частоты возбуждающей нагрузки, образец 1, W1=30м/с2 (иу - перемещения по вертикали, их - по горизонтали)

РОЭТ26

АМРЫТШЖ

иу12 6

(х10**-5)

1 иу

и; 7 /

130 190 200 210 220 230

135 155 205 215 225

±г Нг

1гуЬка р1ре16 дагт 1тг,50

Рис. 4. Зависимость амплитуды от частоты возбуждающей нагрузки, образец 1, W1=50м/с2

копаь эоштгаы

5ТЕЕ-1

БИВ =27 БПЕ(Э=207 ЬЕЙЪ ОМЬУ

ОМХ =.489Е-04 БММ =518834 ЭМХ =.550Е+07

513834 .162Е+07 .273Е+07 .384Е+07 .494Е+07

.107Е+07 .21ВЕ+07 .328Е+07 .439Е+07 .550Е+07

^:гуЬка р1ре!6 дагш 1тг, 30

Рис. 5. Распределение эквивалентных напряжений в трубопроводе на частоте 207Гц, образец 1, W1=30м/с2

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

методиками учитывает зависимость жёсткости и коэффициента рассеивания энергии от амплитуды деформации МР и обладает достаточной простотой для её использования.

Методика может быть использована при доводке трубопроводов за счёт замены существенного объёма работ по тен-зометрированию на расчётное исследование.

Библиографический список

1. Д.Е. Чегодаев, Ю.К. Пономарев. Демпфирование. Самара: СГАУ, 1997. 334 с.

2. Конструкционное демпфирование в неподвижных соединениях Калинин Н.Г., Лебедев Ю.А., Лебедева В.И. и др. / Рига: Изд-во АН Латвийской ССР, 1960. 172 с.

3. Пановко Я.Г. Внутреннее трение при колебаниях упругих систем. М.: Физ-матгиз, 1960. 196 с.

4. Э.Н. Кузьмин, Г.Я. Егоров. Исследование динамических характеристик втулочных амортизаторов // Вибрационная прочность и надёжность двигателей и систем летательных аппаратов: межвуз. сб. Куйбышев, 1975. Вып. I (68). С.54-59.

5. Уланов А.М., Швецов А.В. Механические характеристики демпферов опор трубопроводов из материала МР // Вестн. Самар. гос. аэрокосм. ун-та, 2011. №3(27). С. 94-99.

Информация об авторах

Безбородов Сергей Александрович, инженер-конструктор, отдел прочности и теплофизики, ОАО «Кузнецов». Е-шаП: Ь е2Ьогоёоу8ег@уапёех. ги. Область научных интересов: защита трубопроводов от удара и вибрации.

Уланов Александр Михайлович, доктор технических наук, доцент, профессор кафедры конструкции и проектирова-

ния двигателей летательных аппаратов, Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет). Е-шаП: а!ехи!а[email protected]. Область научных интересов: защита от вибрации и удара при помощи систем сухого трения.

METHOD OF CALCULATING VIBRATIONS IN A PIPELINE WITH DAMPING SUPPORTS MADE OF MR MATERIAL

© 2014 S. А. Bezborodov1, А. М. Ulanov2

1Open joint-stock company "Kuznetsov", Samara, Russian Federation Samara State Aerospace University, Samara, Russian Federation

A method of calculating vibrations in a pipeline of an arbitrary shape with damping supports made of MR material with the help of the finite element method ANSYS software is developed. The method takes into account the non-linearity of MR material, it makes it possible to determine the amplitudes of pipeline vibration and the stress in the pipeline. The calculation results are confirmed by an experiment with the ARAMIS measurement system. The method can be used for placement of pipelines.

Pipeline, vibration, damping, MR material, ANSYS, stress, finite element method.

References

1. Chegodaev D.E., Ponomarev Yu.K. Dempfirovanie [Damping]. Samara: SSAU Publ., 1997. 334 p.

2. Kalinin N.G., Lebedev Yu.A., Leb-edeva V.I. etc. Konstruktsionnoe dempfirovanie v nepodvizhnyh soedineniyah [Structural Damping in Fixed Joints]. Riga: Latvian Academy of Science Publ., 1960. 172 p.

3. Panovko Ya.G. Vnutrennee trenie pri kolebaniyah uprugih sistem [Inner friction in vibration of elastic systems]. Moscow: Fiz-matgiz Publ., 1960. 196 p.

4. Kuzmin E.N. Research of dynamic characteristics of sleeve vibration isolators // Vibratsionnaya prochnost i nadezhnost dvigatelei i sistem letatelnyh apparatov. Kuibyshev, 1975. No. 1(68). P. 54-59. (In Russ.).

5. Ulanov A.M., Shvetsov A.V. Mechanical characteristics of pipeline support dampers made of MR material // Vestnik SGAU. 2011. No. 3(27). P. 94-99. (In Russ.).

About the аuthors

Bezborodov Sergey Alexandrovich,

Design Engineer, Open Joint-Stock Company "Kuznetsov", Strength and Thermal Physics Department. Е-mail: bez-

[email protected]. Area of research: protection of pipelines against vibration and shock.

Ulanov Alexander Mikhailovich,

Doctor of Science (Engineering), Associate Professor, Professor of the Aircraft Design Department, Samara State Aerospace University. Е-mail: [email protected]. Area of research: рrotection against vibration and shock by dry friction systems.