Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

CLARIFYING THE RESULTS OF INFLUENCE OF FLUID ON THE PROPAGATION OF WAVES IN ELASTIC CYLINDRICAL SHELL

Ter-Akopyants G.L. 1
1 SPbGMTU « State Marine Technical University of St.-Petersburg»
Are considered the free propagation of elastic waves in an infinite thin cylindrical shell loaded by an incompressible ideal fluid. Are analyzed bending (beam mode) and asymmetrical about the axis waveform. On the basis of the dynamic equations of the couple-stress theory of shells is obtained and is numerically solved dispersion equation. Are dependences of the axial wavenumber of frequency (dispersion curves) for the first few values of the number of circumential waves. Are found the frequency bands in which there are purely imaginary roots of the dispersion equation that corresponds to the appearance of a propagating wave. Are proposed clarification of the approximate calculation of added masses and analyzed its impact on the axial wave numbers in comparison with used previously approximate calculation and exact solution. Are compared of the results for the shell loaded by an fluid and for empty shell. It is established that the presence of fluid change frequency range of propagating waves and the cut-on frequency of the attenuated waves in the direction of lower frequencies. Are analyzed the impact of the proposed clarification when calculating the added masses on the cut-on frequency.
cylindrical shell
propagation of waves
wave number
cut-on frequency
1. Volmir A.S. Obolochki v potoke zhidkosti i gaza: Zadachi gidrouprugosti. M. Nauka.1979. 320 p.
2. Kalinin V.S., Postnow V.A. Osnovy teorii obolochek. Izd-vo LKI., L. 1979. 200 p.
3. Spravochnik po spetsialnym funktsiyam pod red. Abramovitsa M. i Stigan I. M.: Nauka. 1979. 832 p.
4. Sorokin S. V. Fluid-Structure Interaction and Structural Acoustics. Book of Lecture Notes. Technical University of Denmark. 1997. 188 p.
5. Vijay Prakash S., Venkata R. Sonti Asymptotic expansions for the structural wavenumbers of isotropic and orthotropic fluid-filled circular cylindrical shells in the intermediate frequency range. Journal of Sound and Vibration. Manuscript Draft. Manuscript Number: JSV-D-12-01440. 15 p.

Рассматривается бесконечная цилиндрическая оболочка толщины h с радиусом срединной поверхности R. Колебания такой оболочки под нагрузкой, плотность которой обозначим p1, p2, pn соответственно, можно описать, следуя [1, 2, 4], уравнениями:

Eqn142.wmf (1)

где

Eqn143.wmf – вектор перемещения с осевой, окружной и радиальной компонентами соответственно (знак «+» перед pn и знак «–» перед последней компонентой инерционного члена говорит о том, что положительным считается перемещение по направлению к центру кривизны); ρ – плотность материала оболочки; E – модуль Юнга; ν – коэффициент Пуассона, Eqn144.wmf.

При рассмотрении задачи по моментной теории оболочек матрица оператора L имеет вид

Eqn145.wmf (2)

Перейдём к безразмерной осевой координате Eqn146.wmf и домножим на R2 систему (1).

Матрица полученной системы примет вид

Eqn147.wmf (3)

Будем искать решения преобразованной таким образом системы (1) в виде

Eqn148.wmf (4)

где Eqn149.wmf

Компоненты Eqn150.wmf примем в виде

Eqn151.wmf (5)

где Ω – частота; k = λR – приведенное осевое волновое число; λ – осевое волновое число; m – окружное волновое число (число окружных волн), принимающее значения m = 1, 2, 3… Случай m = 0 (осесимметричные колебания) является предметом отдельного исследования и здесь не рассматривается.

Такой подход позволит искать решение задачи для каждого фиксированного значения окружного волнового числа m независимо.

Предположим, что оболочка заполнена жидкостью, движение которой описывается (см. [1]) уравнением для потенциала скоростей

Eqn152.wmf (6)

где cfl – скорость звука в жидкости.

Выберем модель несжимаемой жидкости, в этом случае будем считать cfl → ∞.

Условия «непротекания» на границе жидкость-оболочка имеют вид

Eqn153.wmf (7)

Наличие жидкости учтем в составе pn уравнения (1а), а p1 и p2 считаем равными нулю.

Решение уравнения (6) будем искать в виде

Eqn154.wmf (8)

Это позволит также искать решение независимо для каждого значения m. Подстановка (8) в (6) с учётом cfl → ∞ приводит к уравнению Бесселя

Eqn155.wmf (9)

Его решение имеет вид

Eqn156.wmf (10)

где C1m и C2m – некоторые константы; Im(rλ) и Ym(rλ) – модифицированные функции Бесселя (см. [3]).

Из-за ограниченности потенциала скоростей при r → 0 имеем C2m = 0, тогда подстановка (8) в (7) с учетом (5) позволяет найти C1m, затем ϕm(R) и получить окончательно давление жидкости на стенку оболочки в виде

Eqn157.wmf (11)

где ρfl – плотность жидкости.

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

Eqn158.wmf(12)

где Eqn159.wmf Eqn160.wmf – приведённая частота.

Перенесем правые части преобразованной системы (1) налево, подставим в полученные соотношения формулы (4) с учетом (5) и (12), уберем символ Eqn161.wmf и разделим все уравнения на Eqn162.wmf. Тогда получим матричное уравнение

Eqn163.wmf (13)

где Eqn164.wmf, с матрицей оператора L* в виде

Eqn165.wmf(14)

В отсутствие жидкости (пустая оболочка) в Eqn166.wmf надо принять ρfl = 0. Приравнивание к нулю определителя матрицы L* приводит к уравнению 8-го порядка для k.

Eqn167.wmf (15)

коэффициенты которого очевидным образом получаются из (14).

При отсутствии жидкости левая часть является многочленом 8 порядка, а при наличии жидкости – нет, так как коэффициенты a4, a2 и a0 зависят от k.

Для оболочки с жидкостью обсудим вопрос об аппроксимации при k ≥ 0 выражения

Eqn168.wmf (16)

Разложение функции (16) в ряд Тейлора по степеням k при m = 1, 2, 3…позволяет при малых k в качестве нулевого приближения взять, как это принято в [4],

f(k) = 1/m. (17)

Однако для k по модулю больших 1 недостаточно даже трёх первых ненулевых членов разложения (то есть до k4), ибо отбрасываемая часть не будет являться малой, а взять большее количество членов ряда не позволяет 8-й порядок уравнения (15).

Расчёт при точных значениях fm(k) возможен с помощью применения схемы итерационного типа, когда при нахождении k на последующем по ω шаге в fm(k) вводятся значения k с предыдущего шага, а значения k при ω = 0 по-прежнему являются корнями многочлена 8 степени, так как при ω = 0, для любого m = 1, 2, 3, … Eqn169.wmf. Однако такой расчёт требует больших вычислительных ресурсов. Поэтому предлагается достаточно простая аппроксимация (18) для fm(k), не отменяющая итерационной схемы, но более экономная по ресурсам. Она получена чисто эмпирическим путём:

Eqn170.wmf (18)

Близость получающихся при этом значений fm(k) при m = 1 и m = 2 иллюстрирует рис. 1.

pic_13.wmf 

Рис. 1. Точные значения fm(k) – красная сплошная линия, нулевое приближение – оранжевая пунктирная линия, предлагаемое приближение – синяя точечная линия. Слева при m = 1; справа при m = 2

На рис. 2–4 приведены корни уравнения (15) при m = 1, 2, 3 для точных значений fm(k), приближенных по формуле (17), и приближенных по формуле (18). В расчетах принималось Eqn171.wmf, Eqn172.wmf. При использовании точных значений fm(k) и приближенных по (18) получаются близкие дисперсионные кривые, неразличимые в масштабах графиков. Сопоставление их с кривыми, полученными с помощью нулевого приближения, показывает, что частоты, на которых зарождаются распространяющиеся волны (появляются чисто мнимые осевые волновые числа k), успешно определяются в том и в другом случае, слегка меняется лишь форма ветвей дисперсионных кривых. Зато частоты отсечки затухающих волн (исчезновение полностью комплексных осевых волновых чисел) меняются существенно. Поэтому вряд ли следует пренебрегать упомянутыми выше уточнениями в компонентах матрицы, связанных с наличием жидкости. Для сравнения на всех рис. 2–4 приводятся дисперсионные кривые для оболочки без жидкости. Для компактности приводятся только положительные значения вещественных и мнимых частей, отрицательные просто имеют противоположный знак. Отметим, что распространяющимся в направлении оси ox волнам соответствуют чисто мнимые волновые числа со знаком « + », волнам, затухающим в направлении ox, – полностью комплексные волновые числа с отрицательной вещественной частью и положительной мнимой, а затухающим в направлении ox неволновым перемещениям соответствуют чисто вещественные отрицательные волновые числа. На рис. 2–4 красной линией обозначены k для оболочки с жидкостью при fm(k) по формуле (17), синей – для оболочки с жидкостью при fm(k) по формуле (18), зелёной – для оболочки без жидкости. Вещественные и мнимые части полностью комплексных k изображены сплошной линией, а чисто вещественные и чисто мнимые корни – пунктирной.

Рис. 2–4 приводят к следующим выводам для рассмотренного диапазона частот.

При m = 1 затухающие волны зарождаются при ω = 0 и прекращаются при ω = 0,420.42; 0,74; 1,02 соответственно для каждого случая (с жидкостью по (17), с жидкостью по (18), без жидкости). Первые две распространяющиеся волны зарождаются при ω = 0 и ω = 0,61, одинаковых для всех случаев, а третья – при ω = 1,09; 1,09; 1,4 соответственно.

При m = 2 две затухающие волны зарождаются при ω = 0, первая прекращается при ω = 0,011; 0,011; 0,019, а вторая прекращается при ω = 0,59; 0,83; 1,10 соответственно для каждого случая. Две распространяющиеся волны в рассмотренном диапазоне частот зарождаются: первая при ω = 0,011; 0,011; 0,019 соответственно для каждого случая, а вторая при ω = 1,20, одинаковых для всех рассмотренных случаев.

pic_18.wmf pic_14.wmf

Рис. 2. Осевые волновые числа k при m = 1 (балочные колебания) в зависимости от частоты ω

pic_16.wmf pic_15.wmf

Рис. 3. Осевые волновые числа k при m = 2 (неосесимметричные колебания) в зависимости от частоты ω. На внутреннем рисунке те же кривые при малых ω до 0,03

pic_17.wmfpic_19.wmf

Рис. 4. Осевые волновые числа k при m = 3 (неосесимметричные колебания) в зависимости от частоты ω. На внутреннем рисунке те же кривые при малых ω до 0,06

При m = 3 две затухающие волны зарождаются при ω = 0 и прекращаются: первая при ω = 0,035; 0,035; 0,055; вторая при ω = 0,73; 0,93; 1,18 соответственно для каждого случая. Две распространяющиеся волны в рассмотренном диапазоне частот зарождаются: первая при ω = 0,034; 0,034; 0,054 соответственно для каждого случая; вторая при ω = 1,79, одинаковых для всех рассмотренных случаев.

Отметим, что в [5] аналогичная задача рассмотрена для ортотропной оболочки.

Рецензенты:

Фёдоров А.С., д.т.н., профессор кафедры сопротивления материалов Санкт-Петербургского государственного морского технического университета, г. Санкт-Петербург;

Хазанов В.Б., д.ф.-м.н., профессор кафедры прикладной математики Санкт-Петербургского государственного морского технического университета, г. Санкт-Петербург.

Работа поступила в редакцию 26.07.2013.