Постановка задачи
Целью работы явилось исследование дисперсионных кривых упругой цилиндрической оболочки, заполненной идеальной несжимаемой жидкостью, в зависимости от соотношения плотностей жидкости и материала оболочки и апробация различных методов их нахождения.
Рассматривается тонкая бесконечная упругая цилиндрическая оболочка толщины hс радиусом срединной поверхности R. Динамические уравнения движения оболочки по моментной теории в безразмерных координатах можно, следуя [2], описать следующим матричным уравнением:
, (1)
где, , , , , , - вектор перемещения с осевой, окружной и радиальной компонентами, – вектор внешней нагрузки (для свободных колебаний без жидкости), – плотность материала, – модуль Юнга, – коэффициент Пуассона,
Движение идеальной жидкости в оболочке в безразмерных координатах описывается [1] уравнением для потенциала скоростей:
. (2)
Будем считать жидкость несжимаемой (скорость звука в жидкости ).
Условие на границе жидкость-оболочка имеет вид:
. (3)
Таким образом, динамические уравнения колебаний оболочки с жидкостью представляют собой систему уравнений (1)–(2) с граничным условием (3). Будем искать её решения в виде суммы решений для каждого отдельно взятого значения окружного волнового числа (числа окружных волн). Это позволит решать систему (1) –(2) независимо для каждого фиксированного .
, (4)
где – частота, – приведённое осевое волновое число(,где – волновое число).
Выбор перемещений в виде (4) означает, что чисто мнимые значения соответствуют распространяющимся волнам, чисто вещественные – апериодическим неволновым перемещениям, а полностью комплексные – затухающим волнам. Всегда имеется две волны, движущиеся в противоположных направлениях. Выбранный вид решений означает, что приволна движется в направлении оси0z и затухает в направлении 0zпри .
Поиск решения в таком виде для потенциала скоростей сводит уравнение (2) к уравнениям Бесселя дляи, как показано в [1] и подробно описано в [4], после удовлетворения условию (3), от системы остается лишь уравнение (1), в правой части которого появляется контактное давление жидкости на стенку оболочки в виде:
, (5)
где – плотность жидкости.
Перенесем правые части уравнения (1) налево, подставим в него перемещения в виде (4) и формулу (5), разделимобе части уравнения на и введём приведённую частоту. Тогда получим матричное уравнение
, (6)
где , с матрицей оператора в виде
, (7)
Дисперсионное уравнение для оболочки с жидкостью при каждом фиксированномm получается, если приравнять к нулю определитель матрицы (7). Оно имеет вид:
, (8)
где.
Дисперсионное уравнение для пустой оболочки получается из (7) , если положить.
Рассмотрим случай. Вычисляя определитель матрицыс помощью разложения по второй строке, получим дисперсионное уравнение в виде совокупности следующих уравнений
(9)
Решение первого из них даёт два чисто мнимых корня, одинаковых для оболочки с жидкостью и без неё. Это говорит о том, что две распространяющиеся волны при m=0, начинающиеся с , имеются всегда, независимо от наличия жидкости. Левая часть второго уравнения в отсутствии жидкости () является многочленом 6-го порядка по k.
Для анализа полученных результатов потребуется также знание осевых волновых чисел жесткого цилиндрического волновода. В этом случае правая часть граничного условия для потенциала скоростей (3) будет равна нулю. В безразмерных координатах получим дисперсионное уравнение, не зависящее от .
. (10)
Численное решение дисперсионного уравнения при отсутствии жидкости является несложной задачей, так как при этом ищутся корни многочлена. При наличии жидкости дисперсионное уравнение является трансцендентным из-за наличия жидкостной добавки в виде функций Бесселя. В общем виде его можно записать в виде , где и сама функция F, и аргумент k являются, вообще говоря, комплексными.
Методы, применяемые для решения дисперсионного уравнения
1. Метод дихотомии для нахождения чисто вещественных и чисто мнимых корней, а также его модификация для нахождения полностью комплексных корней, основанная на том, что в некоторой окрестности любой точки линии пересечения поверхности или с комплексной плоскостью найдутся две точки, в которых функция будет иметь разные знаки.
2. Метод поточечной минимизации модуля функции F.
3. Метод наискорейшего спуска Ньютона
4. Метод, описанный в [3], основанный на численном решении дифференциального уравнения .
Заметим, что методы 3 и 4 требуют для своей реализации стартовой точки, то есть значений k при , а эти значения совпадают со значениями для случая без жидкости, что очевидно из дисперсионных уравнений. В то же время это ограничивает возможность их применения для нахождения кривых, начинающихся при .
Обсуждение результатов
При рассмотрении дисперсионных кривых для оболочки с жидкостью многие исследователи рассматривали только кривые, соответствующие распространяющимся волнам. Однако есть работы, например [5], где приводятся более полные результаты.
Чисто вещественные дисперсионные кривые, представленные на рис. 1, показывают, что, как только пустая оболочка становится заполненной жидкостью, так при сколь угодно малых появляются прямые линии, соответствующие жесткому цилиндрическому волноводу, с точками ветвления, зависящие от отношения плотностей. В них появляются новые кривые, видоизменяющие свою форму по мере роста отношения плотностей жидкости и оболочки. Эта картина наблюдается для всех значений m=0,1,2,3.
Рис. 1. Чисто вещественные дисперсионные кривые при m=0, 1, 2, 3 для оболочки без жидкости (красная линия) и с жидкостью при различных отношениях плотностей. Голубая линия соответствует; синяя – ; фиолетовая –
На рис. 2 представлены чисто мнимые ветви дисперсионных кривых, характеризующие распространяющиеся волны. Они достаточно легко определяются любым из вышеперечисленных методов. Для сравнения красной линией показаны чисто мнимые дисперсионные кривые для оболочки без жидкости. Кривые показывают, что первая и вторая распространяющиеся волны структурного происхождения – наличие и плотность жидкости не влияют на частоту отсечки, с ростом отношения плотностей меняется лишь длина и осевая частота у первой распространяющейся волны. Третья распространяющаяся волна – результат взаимодействия оболочки с жидкостью. Частота отсечки снижается с ростом отношения плотностей жидкости и оболочки. Следует заметить, что при m=2 и при m=3 первые мнимые ветви начинаются не при , а при , но этого не видно в масштабах рисунка.
Рис. 2. Чисто мнимые дисперсионные кривые при m=0, 1, 2, 3 для оболочки без жидкости (красная линия) и с жидкостью при различных отношениях плотностей. Голубая линия соответствует; синяя – ; фиолетовая –
Нахождение полностью комплексных ветвей представляло наибольшие проблемы потому, например, что в отсутствии жидкости полностью комплексная ветвь распадается в точке так называемого, кольцевого резонанса () на две чисто вещественные. При этом в указанной точке фазовая скорость () обращается в бесконечность (касательная к графику вещественной или мнимой части принимает вертикальное положение). Естественно при этом результаты, полученные методами 3 и 4, не могут приниматься как достоверные без дополнительных проверок. Недостатки методов 3 и 4 появляются также и по причине накопления погрешности, так как это итеративные методы. Поэтому, собственно, и возникла потребность в таких простых и очевидных методах, как 1 и 2, но требующих значительных вычислительных ресурсов, так как при каждом значении частоты необходимо вычислять значения комплексной функции, содержащей функции Бесселя, в большом количестве точек комплексной области немалого размера. Тем не менее эти трудности удалось преодолеть. Методы 1–3 показали близкие результаты при , метод 4 дает правильные результаты при . При меньших значениях следует применять лишь методы 1 и 2. Сопоставление различных методов иллюстрирует рис. 3.
Рис. 3. Сопоставление комплексных дисперсионных кривых, полученных различными методами при m=1 и отношении плотностей 0.1. Красная линия – метод 1, голубая – метод 2, тёмно синяя – метод 3, оранжевая – метод 4
Полученные результаты представлены на рис. 4 и означают следующее. При наличии жидкости и с ростом её относительной плотности точка бифуркации постепенно исчезает, а полностью комплексная ветвь, начинающаяся при и заканчивающаяся при , продолжается дальше, оставаясь полностью комплексной. Продолжение (зарождение) комплексной ветви после точки кольцевого резонанса происходит из верхней чисто вещественной ветви, а нижняя чисто вещественная исчезает. Увеличение относительной плотности жидкости ведёт усилению интенсивности затухания (), так и изменению “осевой частоты” (), характер которого можно проследить на рис 4 . Во всех расчетах принималось .
Рис. 4. Вещественные и мнимые части полностью комплексных дисперсионных кривых при m=0, 1, 2 для оболочки без жидкости (красная сплошная – комплексные, красная пунктирная – вещественные или мнимые) и с жидкостью при различных отношений плотностей: синяя линия – 0.08; тёмно-синяя – 0.2; фиолетовая – 0.5; черная – 1
Выводы
При наличии жидкости появляется новый набор апериодических неволновых перемещений (чисто вещественные дисперсионные кривые), изменяющихся с ростом относительной плотности, включающий в себя перемещения жидкости в жестком волноводе.
Частота отсечки третьей распространяющейся волны (чисто мнимые дисперсионные кривые) уменьшается с ростом относительной плотности жидкости.
При наличии жидкости комплексная ветвь не распадается на две вещественные. С ростом относительной плотности интенсивность затухания волны увеличивается.
Рецензенты:
Картузов Е.И., д.т.н., профессор, профессор кафедры теоретической механики Санкт-Петербургского Государственного морского технического университета, г. Санкт-Петербург;
Хазанов В.Б., д.ф.-м.н., профессор, профессор кафедры прикладной математики Санкт-Петербургского Государственного морского технического университета, г. Санкт-Петербург.