CC BY
5Электронный научный журнал "Математическое моделирование, компьютерный и натурный эксперимент в естественных науках" http://mathmod.esrae.ru/ URL статьи: mathmod.esrae.ru/16-54 Ссылка для цитирования этой статьи:
Тер-Акопянц Г.Л., Тер-Акопянц Л.Г. Об особенностях распространения волн при изгибных колебаниях армированной цилиндрической оболочки, заполненной жидкостью // Математическое моделирование, компьютерный и натурный эксперимент в естественных науках. 2017. №4_
УДК 539.3
ОБ ОСОБЕННОСТЯХ РАСПРОСТРАНЕНИЯ ВОЛН ПРИ ИЗГИБНЫХ КОЛЕБАНИЯХ АРМИРОВАННОЙ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ,
ЗАПОЛНЕННОЙ ЖИДКОСТЬЮ
Тер-Акопянц Г.Л.1, Тер-Акопянц Л. Г.2
^анкт-Петербургский Государственный электротехнический университет «ЛЭТИ» им. В.И. Ульянова (Ленина) (СПбГЭТУ «ЛЭТИ»), Россия, Санкт-Петербург, g.ter.a@mail.ru 2Россия, Санкт-Петербург, g.ter.a@mail.ru
ABOUT THE FEATURES OF WAVE PROPAGATION OF A REINFORCED
CYLINDRICAL FLUID-FILLED SHELL IN THE BEAM-MODE OF
VIBRATIONS
Ter-Akopyants G.L1, Ter-Akopyants L. G2
1Saint Petersburg Electrotechical University "LETI", Russia, S ankt-Petersburg, g.ter. a@mail .ru 2Russia, Sankt-Petersburg, g.ter.a@mail.ru
Аннотация. В работе рассмотрено распространение упругих волн в армированной в окружном направлении цилиндрической оболочке, заполненной жидкостью, для изгибного режима колебаний. Сопоставлены дисперсионные кривые и модальные коэффициенты для ортотропной оболочки без жидкости, с несжимаемой жидкостью и со сжимаемой жидкостью. Это позволило выявить влияние ортотропии и жидкости на частоты отсечки и на превалирующий тип перемещений.
Ключевые слова: армированная цилиндрическая оболочка заполненная жидкостью; распространение волн; дисперсионные кривые, модальные коэффициенты.
Abstract. The paper considers the propagation of elastic waves in the circumferential-reinforced cylindrical fluid-filled shell for beam-mode of vibrations. The dispersion curves and modal coefficients are compared for orthotropic shell without fluid, for orthotropic shell filled with an incompressible fluid and for orthotropic shell filled with a compressible fluid. This allowed us to identify the influence of orthotropy and of fluid on the cut-on frequency and on the prevailing type of displacements.
Keywords: reinforced cylindrical fluid-filled shell; wave propagation; dispersion curves; modal coefficients.
В различных технических и промышленных отраслях широко используются трубопроводы большой протяжённости, по которым с невысокой скоростью перекачиваются жидкости. При этом шум и вибрация, создаваемая внешним устройством в одном конце трубопровода, могут передаваться на значительные расстояния. Устранение этих нежелательных явлений требует исследования процессов распространения упругих волн в трубопроводе, заполненном жидкостью. В настоящее время широкое применение имеют элементы трубопроводов из полимерных, армированных материалов, позволяющих снизить вес и стоимость при сохранении прочностных характеристик. Но динамическое поведение такой анизотропной трубы (цилиндрической оболочки) при распространении упругих волн может отличается от изотропной, что и является предметом настоящего исследования. Целью работы является выявление особенностей характера распространения упругих волн в заполненной невязкой жидкостью длинной тонкой цилиндрической оболочке, армированной в окружном направлении. В качестве математической модели используется бесконечно длинная ортотропная оболочка, заполненная несжимаемой или сжимаемой идеальной жидкостью. Полученные результаты могут быть применены также к оболочке со спиральным армированием, когда угол намотки относительно оси близок к 900, если пренебречь эффектами "закручивания" распространяющихся волн.
Впервые исследование волновых процессов в изотропной цилиндрической оболочке с жидкостью на основе анализа дисперсионных кривых было проведено, по-видимому, в [1]. Отметим также ряд более современных работ, касающихся волновых процессов в анизотропных цилиндрических оболочках, заполненных жидкостью [2-5]. Методика нахождения модальных коэффициентов описана в [6]. В работе [7] проанализированы особенности волновых процессов и поведение модальных коэффициентов, найденных методами, сходными с [6], для изгибного режима.
Системы динамических уравнений равновесия цилиндрической оболочки имеет вид:
ôN„ 1 ôSxff , дV - +--хв= ph
2
дх R дв dt
1 dNe dSex 1 дМв 1 дН , дV ПЛ
в++ —г-- +--= ph—2е , (1)
R дв дх R2 дв R дх дt2 дм 2 д2н 1 д2Мв Ne , д2и
дх* R дхдв R2 дв2 R ±п д? где N, Ыв, Sxв, Sвx- погонные усилия, а Ых, Мв, Н- погонные изгибающие и крутящий моменты, р - плотность материала, цп - давление жидкости на стенку оболочки.
r
Для ортотропной оболочки решения системы (1), записанной в перемещениях и учитывающей наличие жидкости внутри оболочки, будем искать в виде:
( и, Л (ит со$>(тв) Л
х
_ х-гОI
V иг у
К ып(тд)
у!¥т СОЪ(тв)у
(2)
где к*- осевое волновое число, О - частота, т - фиксированное число окружных волн.
Переход к безразмерным координатам, введение приведённого осевого волнового числа к _ к* • R и приведенной частоты с, связанной с О формулой
со2 _ р 2(1 1/'2^21 )R2О2, приводит для каждого фиксированного т к следующей
Е + Е2
линейной однородной системе для амплитуд перемещений, называемых модальными коэффициентами:
(иЛ ( о Л
L •
V W
(3)
V" у V"; где Ьп _ Ь1 к2 - Ьат 2 + с
Ll2 _(Ь]^21 + Ьа )тк,
^3 _ ь1^21& ,
L2l _-(Ь2^12 + Ьа )тк,
L22 _ Ьак2 - Ь2т2 + а(- Ь2т2 + 4Ьак2 )+с2,
L23 _ -Ь2т + а{Ь2Уи + 4Ьа )тк2 - Ь2т3},
_-Ь2^12 & =
L32 _ -Ь2т + а|(Ь1У21 + 4Ьа )тк2 - Ь2т3},
L33 _ -Ь2 + а{- Ь1к4 + (Ь1у21 + Ь2у12 + 4Ьа )т2к2 - Ь2т4 }+ с2 ,2
(1 + —(&С)' р п
а _ -Ш- Ь _ 2Е1 Ь _ 2Е2 Ь _ 2а(1 -У12^21) . Е _ Е Е _ Е -
12R2' 1 Е1 + Е2' 2 Е1 + Е2' а Е1 + Е2 ' 1 ^ 2 0 продольный и окружной модули Юнга; у21 , у12 - коэффициенты Пуассона,
Е Е
а --; п
- толщина, R-срединный радиус оболочки;
Е1(1 + ^21)+ Е2 (1 + ^12 )
ие
(к ,0) =
3
т
функция
т3т
Л
к2+^ с//'
0
1
к2+^
0
3
т+1
Л
к2 + -с
0
у/
Бесселя; р и р/ - плотности оболочки и жидкости, с и с/ - скорости звука в
оболочке и в жидкости соответственно.
Изгибному режиму колебаний соответствует значение т = 1. Дисперсионное уравнение получается из равенства нулю определителя матрицы Ь.
det Ь = 0. (4)
Графические изображения его решений (зависимостей к = к (0)) представляют собой дисперсионные кривые. При выборе перемещений в виде (2) распространяющимся волнам соответствуют чисто мнимые значения к.
В качестве примера рассмотрена тонкая цилиндрическая оболочка ^/К=0.05) из полимерного материала с 10% стальным армированием в окружном направлении. При этом отношение окружного и продольного модулей Юнга принималось равным 1:25, коэффициент Пуассона основного материала равным 0.4, р/ : р = 0.621, с: с/ = 1.333. В используемой модели
ортотропного материала использовалось классическое [8] для ортотропной теории упругости соотношение Е1у21 = Е2у12 , что, вероятно, не вполне справедливо для композитного материала. Однако было проверено, что зависимость волновых чисел от коэффициентов Пуассона весьма слаба, и, даже если принять у21 = v12 = 0.4, то на качественные выводы это не повлияет.
Для решения поставленной задачи были рассмотрены и сопоставлены дисперсионные кривые и модальные коэффициенты для изотропной и ортотропной оболочек без жидкости, с несжимаемой жидкостью и со сжимаемой жидкостью.
Рисунки 1-2 позволяют проанализировать влияние окружного армирования для пустой оболочки. На них представлены дисперсионные кривые и модальные коэффициенты для изотропной и ортотропной оболочки без жидкости. Сопоставление результатов приведено в таблице 1. (Символами Т и X обозначается возрастание и убывание соответствующей величины с ростом 0).
Влияние окружного армирования проявляется главным образом в том,
что:
- снижается частота зарождения второй распространяющейся волны (оранжевая дисперсионная кривая);
2
- при этом с ростом частоты переход превалирующих преимущественно продольных перемещений в преимущественно окружные в изотропном случае, при наличии армирования не наблюдается;
- повышается частота зарождения третьей распространяющейся волны (зелёная дисперсионная кривая);
- при этом с ростом частоты происходит переход радиально-окружных перемещений в преимущественно окружные, а не в преимущественно продольные, как в изотропном случае.
Рис. 1. Дисперсионные кривые и модальные коэффициенты для изотропной оболочки без
жидкости
Рис.2. Дисперсионные кривые и модальные коэффициенты для ортотропной оболочки без
жидкости
Таблица 1.
Частоты отсечки и превалирующие перемещения для оболочки без жидкости
1-красная 2-жёлтая 3-зелёная
изотр. орт. изотр. орт. изотр. орт.
Частота 0 0 0.55 0.30 1.45 1.95
зарождения
Превалирующий
тип перемещений а) при
зарождении б) в конце V-W и и V-W V-W
интервала для с W W V и и V
Характер V i V i U i U t U t
изменения W t W t V t V i V i
модальных коэффициентов с W i W i
ростом с.
Рисунки 3-4 позволяют проанализировать влияние окружного армирования для оболочки, заполненной несжимаемой жидкостью. На них представлены дисперсионные кривые и модальные коэффициенты для изотропной и ортотропной оболочки с несжимаемой жидкостью. Сопоставление результатов приведено в таблице 2. (Символами tt ii обозначается интенсивное возрастание / убывание соответствующей величины с ростом с.) Отметим, что наличие несжимаемой жидкости внутри изотропной оболочки приводит к снижению частоты отсечки третьей распространяющейся волны (зелёная дисперсионная кривая). В остальном качественное влияние ортотропии полностью аналогично случаю пустой оболочки. Это объясняется тем, что в несжимаемой жидкости упругие волны не распространяются.
Рис.3. Дисперсионные кривые и модальные коэффициенты для изотропной оболочки с
несжимаемой жидкостью
Рис.4. Дисперсионные кривые и модальные коэффициенты для ортотропной оболочки с
несжимаемой жидкостью
Таблица 2.
Частоты отсечки и превалирующие перемещения для оболочки с несжимаемой
жидкостью
1 -красная 2-жёлтая 3-зелёная
изотр. орт. изотр. орт. изотр. орт.
Частота 0 0 0.55 0.30 1.00 1.45
зарождения
Превалирующий
тип
перемещений
а) при
зарождении V-W V-W U U V V
б) в конце
интервала для с W W V U U V
Характер V i V i U ii U i U t U t
изменения W t W t V tt V t Vi Vi
модальных
коэффициентов
с ростом с.
Рисунки 5-6 позволяют проанализировать влияние окружного армирования для оболочки, заполненной сжимаемой жидкостью. На них также представлены дисперсионные кривые и модальные коэффициенты. Результаты приведены в таблице 3. Прежде всего, отметим, что наличие сжимаемой жидкости, способной распространять упругие волны, приводит к появлению в рассматриваемом частотном диапазоне четвёртой распространяющейся волны жидкостного происхождения (голубая дисперсионная кривая). Такая же дисперсионная кривая (показано на рис. 5-6 чёрной точечной линией) описывает распространение волны в жидкости в абсолютно жёстком акустическом волноводе [9] . Влияние окружного армирования на первую и вторую распространяющиеся волны (красная и оранжевая дисперсионные кривые) такое же, как в отсутствие жидкости или с несжимаемой жидкостью. Третья распространяющаяся волна в изотропном случае появляется как волна структурного происхождения и характеризуется преимущественно окружными перемещениями, доля которых с ростом частоты уменьшается в пользу продольных, но это происходит только до частоты (о = 1.75, на которой зелёная и голубая дисперсионные кривые сближаются (veering). На этой частоте происходит интенсивное взаимодействие жидкости с оболочкой по радиальной координате, обмен энергией между оболочечной и жидкостной волнами (см [10]), в результате чего зелёная дисперсионная кривая, начиная с этой частоты, описывает преимущественно жидкостную волну, характеризующуюся нарастающими радиальными перемещениями, а голубая - преимущественно
структурную с преобладанием осевых перемещений. При этом происходят достаточно резкие изменения модальных коэффициентов. А в случае окружного армирования увеличивается частота зарождения третьей структурной волны (зелёная дисперсионная кривая) и эффекты, связанные со сближением дисперсионных кривых не наблюдаются.
Рис.5. Дисперсионные кривые и модальные коэффициенты для изотропной оболочки со сжимаемой жидкостью (точками показана дисперсионная кривая для абсолютно жёсткого
акустического волновода)
Рис.6. Дисперсионные кривые и модальные коэффициенты для ортотропной оболочки со сжимаемой жидкостью (точками показана дисперсионная кривая для абсолютно жёсткого акустического волновода, кружками - дисперсионные кривые для ортотропной оболочки с
несжимаемой жидкостью)
Таблица 3.
Частоты отсечки и превалирующие перемещения для оболочки со сжимаемой
жидкостью
1-красная 2-жёлтая 3-зелёная 4-голубая
изотр. орт. изотр. орт. изотр. орт. изотр. орт.
частота зарождения 0 0 0.55 0.30 1.00 1.45 1.40 1.40
превалирующий
тип перемещений а) при
зарождении б) в конце V-W V-W U U V V V-W V
интервала для с W W V U U V U W
Характер v ii V i U ii U i при U t при V i
изменения модальных коэффициентов с W tt Wt Vtt Vt w < 1.75 U tt Vii Vi w < 1.75 Utt Vii Wt
ростом с. W i при w > 1.75 Uii Vii W tt W ii при w > 1.75 Ut Vi W t
В заключение отметим, что во многих отраслях техники, где применяются трубопроводы большой протяженности, длинноволновые колебания, которые возникают при малых значениях Im(k), нежелательны, а окружное армирование позволяет сместить их в другой частотный диапазон. Также неблагоприятно сказываются на работе трубопровода резкие изменения преобладающего типа перемещений, чего можно избежать, если устранить зону сближения дисперсионных кривых за счёт окружного армирования. Таким образом, полученные результаты могут быть полезны в вопросах борьбы с шумом и вибрацией трубопроводов.
Литература
1. Fuller C.R., Fahy F.J. Characteristics of wave propagation and energy distributions in cylindrical elastic shells filled with fluid // Journal of Sound and Vibration (1982) 81(4), pp. 501-518.
2. Jeong Ho You, K. Inaba. Fluid-structure interaction in water-filled thin pipes of anisotropic composite materials // Journal of Fluids and Structures, 36(0), 2013. 1. pp. 162-173.
3. Kunte M.V., A. Sarkar, and V.R. Sonti. Asymptotic expansions for the coupled wavenumbers in an infinite orthotropic flexible fluid-filled cylindrical shell // Journal of the Acoustical Society of America, 131(6), 2012. pp. 4272-4282.
4. Xi Z.C., Liu G.R., Lam K.J., Shang H.M. Dispersion and characteristic surfaces of waves in laminated composite circular cylindrical shells // J. Acoust. Soc. Am. 2000. V.108 (5), pp.2179-2186.
5. Prakash Vijay S., Sonti V.R. Asymptotic expansions for the wavenumbers in orthotopic fluid-filled circular cylindrical shells for intermediate fluid loading // Journal of Acoustical Society of America, 135(6), June 2014. pp. 3198-3208.
6. Liu Yu-Cheng, Yun-Fan Hwang, Jin-Huang Huang. Dispersion relations and modal patterns of wave in a cylindrical shell // Wave process in classical and new solids, Prof. Pasquale Giovine (Ed.), 2012, ISBN: 978-953-51-0821-4, in tech, DOI: 10.5772/50477.
7. Тер-Акопянц Г.Л. Дисперсионные кривые и модальные коэффициенты при распространении волн в оболочке с жидкостью // Естественные и технические науки. 2015. № 6(84). C.77-81.
8. Лехницкий С.Г. Теория упругости анизотропного тела. М.: Наука,1977. 416 с.
9. Тер-Акопянц Г.Л., Тер-Акопянц Л.Г. Об идентификации структурных и жидкостных волн в упругой цилиндрической оболочке с жидкостью на основе анализа дисперсионных кривых // Международный журнал прикладных и фундаментальных исследований. 2016. № 10-2. С. 197-200; URL: http://www.applied-research.ru/ru/article/view?id=10316 (дата обращения: 24.10.2016).
10. Филиппенко Г.В. Энергия волн изгибного типа в бесконечной цилиндрической оболочке, заполненной сжимаемой жидкостью // Современное машиностроение: Наука и образование: материалы 5-й Международной научно-практической конференции. СПб.: Изд-во Политехн. ун-та, 2016. с. 276-286.
