Введение 3
1 Вспомогательные сведения 5
1.1 Уравнение неразрывности 5
1.2 Уравнение движения 6
1.3 Уравнение Лапласа 7
1.4 Интеграл Коши-Лагранжа 8
2 Колебания битого льда под действием внешней нагрузки в
движущейся системе координат 11
2.1 Колебания невязкого битого льда 11
2.2 Периодические волны распространяющиеся вдоль канала . . 21
2.3 Численные результаты 22
2.4 Модель вязкого битого льда 24
3 Аналитическое исследование классического решения задачи о нестационарных колебаниях битого льда 32
3.1 Постановка задачи 32
3.2 Единственность классического решения 33
3.3 Устойчивость решения по начальным данным 36
4 Поведение прогибов битого льда при больших временах 38
4.1 Постановка задачи 38
4.2 Метод решения 39
4.3 Анализ результатов 48
Заключение 51
Список литературы
Во многих прикладных задачах, учитывающих взаимодействие ледового покрова и жидкости, жидкости и конструкций и т.п., возникает необходимость изучения вынужденных волновых колебаний жидкости. Вынужденные колебания являются следствием приложения внешних сил к поверхности жидкости. Поверхность жидкости может рассматриваться как классическая свободная поверхность [1], флотирующая жидкость (ненулевая масса поверхности) [2] или упругая поверхность (поверхность с ненулевой массой и обладающая изгибными силами) [3]. Исследования задач волнового движения жидкости, покрытой тонкой упругой пластиной имеет практическое значение для природных (ледовый покров) и искусственных (плавающие гибкие волнорезы, плавающие аэропорты и т.п.) структур [3-5]. Классическим подходом к изучению динамики ледового покрова является линейная теория гидроупругости [6-8]. В этом подходе лед моделируется тонкой упругой пластиной в рамках теории Кирхгофа-Лява (или, в двумерном случае, балкой Эйлера Бернулли), а жидкость считается идеальной и несжимаемой. Обычно под решением таких задач понимаются функция отклонения поверхности от положения покоя и потенциал скорости течения. Форма вынужденных колебаний определяется как комбинация свободных волн с разными волновым числом и частотой волны.
Искусственное разрушение ледового покрова судном на воздушной подушке и/или разрушение льда волнами на открытой воде приводит к образованию зоны битого льда. Характеристики битого льда отличаются от характеристик свободной поверхности и от сплошного ледового покрова, что приводит к изменению поведения распространяемых волн. При движении судна по битому льду основной прикладной задачей является определение параметров, гарантирующих безопасное движение. При отсутствии упругих сил в битом льду могут возникать волны большой амплитуды, как в области движения судна, так и в отдалении от него, в частности на стенках канала и возле речных сооружений. Рассматриваемая задача является близкой по постановке к задачам движения кораблей в морских льдах по полынье, образованной впереди идущим ледоколом, или движение по приграничной зоне ледового покрова. Поэтому, при решении задач на вынужденные колебания, возникает задача определения возможных свободных колебаний, или, в случае неограниченной области - прогрессивных (распространяющихся) волн, для заданного набора физических параметров и начально-краевых условий задачи [9, 10]. Основными искомыми характеристиками прогрессивных волн являются зависимость частоты волны от волнового числа и скорости распространения волн в пространстве. Предполагается, что лед был предварительно разрушен резонансным методом судном на воздушной подушке, которое движется вдоль замороженной части реки на определенной скорости и создает напряженно-деформированное состояние ледового покрова. Вызванные напряжения могут быть достаточными для ломки льда [11,12].
Случай когда верхняя поверхность канала описывается ледовым покровом хорошо изучена в прошлом. Прогрессивные волны в замороженном канале исследованы в [10, 13, 14]. Колебания неограниченной тонкой ледовой пластины исследованы в [15,16], колебания полубесконечного ледового покрова в [17]. Рассмотренные задачи решались в рамках линейной теории гидроупругости. При исследовании вынужденных колебаний внешняя нагрузка моделировалась гладким локализованным пятном давления. Показано, что стенки канала имеют важную роль в формировании прогибов ледового покрова. Исследование колебаний ледового покрова проведено в [18,19]. Оценка влияния периодической нагрузки на ледовый покров получена в [20,21]. Вопросы корректности начально-краевых задач динамики пороупругого льда рассмотрены в [22, 23].
Основной целью данной работы является постановка и численное исследование задачи о колебаниях битого льда в канале, вызванных движением внешней нагрузки. Битый лед моделируется флотирующей жидкостью. В рассматриваемой модели учитывается условие затухания волн в отдалении от нагрузки. Битый лед имеет постоянную толщину. Жидкость в канале предполагается идеальной и несжимаемой [24]. Движение жидкости в канале является потенциальным и вызвано отклонением битого льда от состояния покоя. Краевые условия на битый лед на стенках канала отсутствуют.
В работе исследована задача о движении внешней нагрузки с постоянной скоростью вдоль центральной линии канала. Канал заполнен идеальной жидкостью и покрыт битым льдом. Жидкость в канале невязкая и несжимаемая. Битый лед моделируется флотирующей жидкостью. Течение, вызванное прогибами битого льда, потенциальное. Нагрузка моделируется гладким локализованным пятном давления.
Исследованы установившиеся прогибы битого льда в системе координат, движущейся совместно с нагрузкой. Рассмотрены случаи вязкого и невязкого битого льда. Исходная задача с помощью преобразования Фурье сведена к задаче относительно профиля прогибов льда поперек канала, которая решена разложением на систему функцию, описывающую колебания жидкости поперек канала. Показано, что при движении нагрузки с постоянной скоростью форма прогибов льда имеет волновую структуру. Наблюдаемые волны распространяются сзади нагрузки. Вязкость льда увеличивает затухание волн в отдалении от нагрузки. Наиболее сильный эффект вязкости льда наблюдается для коротких волн. При увеличении скорости движения нагрузки уменьшается амплитуда прогибов битого льда и увеличивается длина вынужденной волны. При увеличении ширины канала уменьшаются амплитуды прогибов битого льда.
Исследована нестационарная задача о колебаниях битого льда в канале под действием движущейся нагрузки. Доказаны теоремы единственности и устойчивости по начальным данным классического решения задачи. Асимптотическими методами решения интегралов с большим параметром получена аналитическая формула для прогибов битого льда при больших временах. Полученная формула дает простую интерпретацию формирования прогибов битого льда в виде суммы локальных прогибов возле нагрузки и бесконечной системы волн, распространяющейся от нагрузки с ее же скоростью. Для этой системы волн исследована сходимость соответствующих амплитуд. Амплитуды волн уменьшаются при увеличении частоты волны. Для аккуратного определения прогибов битого льда необходимо учитывать от 20 до 30 волн с наибольшей частотой. Максимальные амплитуды для рассматриваемого случая наблюдались для небольшого значения скорости движения нагрузки.
1. Кочин Н. Е., Кибель И. А., Розе Н. В. Теоретическая гидромеханика. Часть 1. — 6-е изд., перераб и дополн. - М.: Физматлит, 1963. - 584 с.
2. С. А. Габов, А. Г. Свешников Математические задачи динамики флотирующей жидкости. — Итоги науки и техн. Сер. Мат. анал., 28, ВИНИТИ, М., 1990, 3-86; J. Soviet Math., 54:4 (1991), 979-1041.
3. V. Squire, R. Hosking, A. Kerr, P. Langhorne. Moving loads on ice. — Kluwer Academic Publishers, 1996.
4. M. Kashiwagi. Transient responses of a VLFS during landing and take-off of an airplane. — Journal of Marine Science and Technology 2004;9(1):14-23.
5. I.K. Chatjigeorgiou, A.A. Korobkin, M.J. Cooker. Three-dimensional steep wave impact on a vertical cylinder. — Journal of Hydrodynamics, Ser. B 2016, 28(4), 523-533.
6. DYe Kheysin. Moving load on an elastic plate which floats on the surface of an ideal fluid. — (in Russian) Izv. Akad. Nauk SSSR, Otd. Tekh. Nauk, Mekh. i Mashinostroenie 1963 (1); 178-180.
7. DYe Kheysin. Dynamics of floating ice covers. — (in Russian) Gidrometeorologichoskoe Izdatel’stvo, Leningrad, Technical Translation FSTC-HT-23-485-69. US Army Foreign Science and Technology Center 1967.
8. P.F. Sabodash, I.G. Filippov. A dynamic problem for a thin elastic plate. — International Applied Mechanics 1967; 3(6): 28-31.
9. A.A. Korobkin, T.I. Khabakhpasheva, A.A. Papin. Waves propagating along a channel with ice cover. — Eur J Mech B/Fluids 2014;47:166-175.
10. E.A. Batyaev, T.I. Khabakhpasheva. Hydroelastic waves in channel with free ice cover. Fluid Dynamics 2015;6:84-101.
11. K.A. Shishmarev, T.I. Khabakhpasheva, A.A. Korobkin. The response of ice cover to a load moving along a frozen channel. — Applied Ocean Research. 2016. Т. 59. С. 313-326.
12. V.M. Kozin. Resonance Method of Breaking of Ice Cover. Inventions and Experiments. — Akad. Estestvoznaniya, Moscow, 2007.
13. А.А. Коробкин, А.А. Папин, К.А. Шишмарев. Аналитическое и численное исследование квазиизотермической задачи взаимодействия ледового покрова канала и поверхностных волн. — Известия АлтГУ, Барнаул, 2012. Вып. 1/2 (73). с.23-27.
14. А.А. Коробкин, А.А. Папин, К.А. Шишмарев. Поведение ледового покрова канала под действием поверхностных волн. — Известия АлтГУ, Барнаул, 2012. Вып. 1/1 (73). с.55-59.
15. V. D. Zhestkaya. Numerical solution of the problem of an ice sheet under a moving load. — Journal of Applid Mechanics and Technical Physics, 1999, V 40(4): 770-775.
16. В.Д. Жесткая, В.М. Козин. Численное решение задачи о воздействии ударного импульса на ледяной покров. — ПМТФ, 2008, Т.49, №2.
17. P. Brocklehurst. Hydroelastic waves and their interaction with fixed structures. — PhD thesis, University of East Anglia, UK, 2012.
18. I.V. Sturova, L.A. Tkacheva. Wave motion in a fluid under and inhomogeneous ice cover — Journal of Physics: Conference Series. 2017. Т. 894. № 1. С. 012092.
19. И.В. Стурова, Л.А. Ткачева. Колебания ограниченного ледяного покрова при локальном динамическом воздействии — Полярная механика. 2016. № 3. С. 997-1007.
20. Л.А. Ткачева. Колебания ледяного покрова с трещиной при воздействии периодической по времени нагрузки — Известия Российской академии наук. Механика жидкости и газа. 2017. № 2. С. 54-64.
21. L.A. Tkacheva. Vibrations of an ice sheet with crack under a time-periodic load — Fluid Dynamics. 2017. Т. 52. № 2. С. 219-229.
22. М.А. Токарева. Конечное время стабилизации решения уравнений фильтрации жидкости в пороупругой среде. — Известия Алтайского государственного университета. 2015. Т. 2. № 1. С. 153-157.
23. M.A. Tokareva. Solvability of initial boundary value problem for the equations of filtration in poroelastic media. — Journal of Physics: Conference Series. 2016. Т. 722. № 1. С. 012037.
24. V.M. Kozin, V.D. Zhestkaya, A.V. Pogorelova, S.D. Chizhumov, M.P. Dzhabailov, V.S. Morozov, A.N. Kustov. Applied problems of the dynamics of ice cover.. — Academy of Natural Sciences, Moscow; 2008.