Устойчивость решения по начальным данным задачи о колебаниях ледового покрова в канале

УДК 517.95 + 532.58

Устойчивость решения по начальным данным задачи о колебаниях ледового покрова в канале*

К.А. Шишмарев

Алтайский государственный университет

Stability of Ice Cover Oscillations in a Channel Problem Solution with Initial Data

K.A. Shishmarev

Altai State University

В рамках линейной теории гидроупругости рассматривается начально-краевая задача о колебаниях ледового покрова в бесконечном канале, вызванных движением внешней нагрузки. В основе математической модели лежит связанная система дифференциальных уравнений, описывающая вертикальное отклонение ледового покрова и движение жидкости в канале. Ледовый покров моделируется уравнением тонкой вязкоупругой пластины. Вязкоупругие свойства льда моделируются на основе реологического закона Кельвина -Фойгта. Жидкость в канале невязкая и несжимаемая. Связь двух задач заключена в линеаризованных условиях - кинематическом и динамическом. Система уравнений замыкается следующими условиями: жесткого защемления для пластины на стенках канала; непротекания для потенциала скорости течения; затухания колебаний на бесконечности. Исследования в данной работе посвящены проблемам корректности постановок задач, описываемых совместными уравнениями динамики вязкоупругой пластины и идеальной жидкости. В пункте 1 доказана теорема об устойчивости по начальным данным классического решения начально-краевой задачи вязкоупругих колебаний ледового покрова в канале. В пункте 2 доказан аналог теоремы пункта 1 для упругих колебаний ледового покрова. Теоремы доказаны с использованием методов получения энергетических оценок, а также специальных функциональных неравенств. Вопросы устойчивости описанных задач исследованы при конечных временах.

Ключевые слова: устойчивость, начально-краевая задача, идеальная жидкость, вязко-упругие колебания, ледовый покров, внешняя нагрузка.

БО! 10.14258,^а8и(2017)4-30

*Работа выполнена при финансовой поддержке гранта РФФИ №16-08-00291.

In this paper, an initial-boundary problem of ice cover oscillations in an infinite channel caused by a moving load is considered within the linear theory of hydroelasticity. Mathematical model is based on a combined system of differential equations describing vertical deflection of the ice cover and motion of fluid in the channel. The ice cover is modeled by the equation of a thin viscoelastic plate. Viscoelastic properties of ice are modeled by the Kelvin-Voigt rheological law. The liquid in the channel is inviscid and incompressible. The relation between the viscoelastic and hydrodynamic parts of the problem is established by the linearized kinematic and dynamical conditions. The system of equations is closed by conditions of clamped edges for the plate on the walls of the channel, conditions of impermeability for the flow velocity potential, and damping conditions for the oscillations far away from the load. The investigation in this paper is devoted to the correctness of problems formulations described by combined equations of a viscoelastic plate and an ideal fluid. In section 1, the stability theorem for the classical solution of the initial boundary-value problem of viscoelastic oscillations of the ice cover in the channel with respect to initial data is proved. In section 2, the similar theorem is proved for elastic oscillations of the ice cover. Theorems are proved using energy estimation methods and special functional inequalities. Stability of described problems is studied for finite times.

Key words: stability, initial-boundary value problem, ideal fluid, viscoelastic oscillations, ice cover, external load.

Введение. Рассмотрим следующую систему уравнений:

Ар

+ D ( 1 + т-^) = Р(х, у, г)-

= рр(х,у, о, г) - р1ди>(х, у, г) (х,у,г) е Пт,

решаемую в областях Пт = П х (0,Т), Пт = П х (0,Т), П = {-то < х < то, -L < у < L, -Н < z < 0}, П = {-то < х < то, —L < у < Ь} при следующих начально-краевых условиях:

р.у =00 (у = ±Ь); рг =0 ^ = -Н);

р(х,у, 0, г) = ,ш(х,у,г); IV = и>у = 0; (у = ±Ь);

)^(х,у, 0) = w1(x,y); 'ш^х,у, 0) = у)2(х,у);

V ^ 0, wx ^ 0, р ^ 0

Данная начально-краевая задача описывает колебания ледового покрова, "примороженного" к стенкам бесконечного канала [1,2]. Здесь (х, у, г) - эйлеровы координаты; р(х,у, z,t), w(x,y,t) -соответственно потенциал скорости течения идеальной жидкости и прогиб ледового покрова; Р(х,у,г) - внешняя нагрузка (заданная функция своих аргументов); постоянные D > 0, М > 0, Р1 > 0, д > 0, т > 0 - соответственно изгибная жесткость ледовой пластины, масса льда на единицу площади, плотность жидкости, ускорение силы тяжести, коэффициент запаздывания вязко-упругого материала Кельвина - Фойгта; V4 = А = 0X4 + 2 вх2ву2 + дх. Искомыми являются функции р(х,у, z,t), w(x,y,t).

В данной работе под решением задачи (1) - (3) понимается пара функций w(x,y,t) и р(х,у^,г), определенных на Пт и Пт соответственно и удовлетворяющих рассматриваемой системе уравнений и краевым и начальным условиям, как непрерывные в Пт и Пт функции. Совместная система уравнений динамики тонкой вязко-упругой или упругой пластины и жидкости исследовалась численно и аналитически для разных начально-краевых условий. (см., например, [3-7] и др.). Обзор работ, посвященных задачам с неограниченной тонкой упругой пластиной, приведен в [1]. Рассматриваемая задача была исследована численными методами в случае примороженного ледового покрова к двум стенкам [5,8,9], к одной стенке [10], и в случае незакрепленного ледового покрова в канале [11]. Задачи с телом, погруженным в жидкость, при наличии ледового покрова рассмотрены в [12, 13].

Разрешимость начально-краевых задач для уравнений Эйлера и динамики тонкой пластины подробно исследованы по отдельности [14-16].

Единственность решения в задаче протекания с заданным вихрем для нестационарных уравнений Эйлера идеальной несжимаемой жидкости доказана в работе [14]. Корректность начально-краевых задач для уравнения колебаний балки исследована в [17,18]. Существование решения задачи изгиба упругой пластинки с жестким включением доказана в [15]. Однозначная разрешимость задачи сопряжения двух упругих однородных балок, одна из которых моделируется балкой Эйлера — Бернулли, а другая балкой Тимошенко, доказана в [16].

Основной результат данной работы состоит в теореме, сформулированной далее в пункте 1.

1. Случай вязко-упругих колебаний. Пусть существуют два отличных от нуля решения Wl,рl и W2, Р2 системы (1) - (3), удовлетворяющих разным начальным условиям:

Wl(x,y, 0) = w^(x,y) Wl,t(x,y, 0) = 1л2(х,у),

W2(х, у, 0) = w'1(x,y) W2,t(x, у, 0) = ^^х^у).

Функции W0 = W^(x,y) —W1(x,y) и w0 = w'2(x,y) -w'22(x,y) удовлетворяют следующему неравенству:

Р1Цр°

где

У(х,у,г) = {шХХх(х,у,г) + Ц)у(х,у,г) +^1у(х,у,г));

У(х,у, 0)= ^хх + W02yy + 2w02xy) ;

р0(х,у^) = ф1(х,у^, 0) - ф2(х,у^, 0), и || • ||2,П и || • || 1 ,п - соответствующие нормы функций в пространствах Ь2(П) и Ь^П).

Функции w = Wl - W2 и р = р1 - р2 удовлетворяют следующей задаче:

+ - + т^ V4w = РlРt - РШ

Ру = 0 (у = ±Ь), рг =0 ^ = -Н);

рх = wt ^ = 0);

w = 0, -Шу = 0

р(х,у^,г), >^(х,у,г), "х (х, у, г) ^ 0 (|х| ^ то);

^п^^у, 0) = wi0; wt

Теорема 1 Для решения задачи (5) - (10) справедливы следующие оценки:

где С - константа, зависящая от параметров задачи.

Заметим, что решение ) задачи (5) - (10) об- Интегрируя (14) по t, с учетом условий (10), при-ладает следующими свойствами: ходим к тождеству

J w

+ ^ m n -n t t

П DrV4w + D J V4wdr + Mwt + p< + Pi9 J wdr =

J \Уф,y, z, Щ2 Ш = j Ф,y,0, t)wt dn,

J VF • V<pd£l = J Fwtdn, (13) Умножим тождество (15) на wt(x, y, t) и проинте-

n п грируем по области Ид = {-R < x < R, —L <y <

где F = F(x,y) - произвольная непрерывная L}. Получим гладкая ограниченная в П функция.

Интегрируя уравнение (6) по области Пд = г i I t \

{—R < x < R, —L < y < L, —H<z< Dt J wtviwdnR +D J \J v4wdT I wt dnR+

0}, учитывая условия (7) и формулу Гаусса- Пд Пд \0 /

Остроградского, получим после предельного пе- [о f

рехода при R ^ ж: + M J W dnR + Pi J <wt dnR+

Пд Пд

A<dn=l <z (x,y, 0,t) dn= /* \

n п + pig I wdr I wt dnR = V0wt dnR.

= [ wt(x,y,t) dn= dfw(x,y,t) dn = 0. Пд ' Пд

П П Последовательно интегрируя по частям слагаете мые в тождестве (16), используя условия (7), (9) Привлекая условия (10), выводим (11). „

,, /г-\ и свойство (12) и совершая в итоговом тождестве

Умножим теперь уравнение (6) на < и проин-

' предельный переход по R, получим тегрируем полученное равенство по Пд . С уче-

т°м (формулы 1^г1>нсса-Остроградского и условий M Г W2dn+Pl f \V<\*dn+ — dy(t)+Dd2Z(t) + (7) после предельного перехода при R ^ ж при- J t ' J 1 2 dt 2 dt2

ходим к равенству П n

_ d2

2

dn

— J \V<\2 dn + J <p(V< • n) dm = + pg dto j u2 dn = DY(t)+pig f u2 dn+J V% dn,

П дП

= ^У \V<\2 dtt + J <(x,y, 0,t)wt dn = 0, где

П П П

z $) = ! (и2хх + Ь2уу + 2п1у )Ж.

из которого следует (12). п

Умножим теперь уравнение (6) на F(х, у) и Используя неравенство Коши, получим проинтегрируем полученное равенство по Пд.

С учетом формулы Гаусса - Остроградского I V0wtdП < М I V02¿П.

и условий (7) после предельного перехода при ^ 2 ^ 2М ^

R ^ ж приходим к равенству

Используя последнюю оценку, усилим правую - УF ¿П + ! F (х, у)(У^  п) ¿д П = часть тождества (17) П дП М [ 2

— J V< • VF dn + J Fwt dn = 0,

n П

M J w2dn + pi J\V<\2dn+

П

Dr d D d2 plg d2 f 2

из которого следует (12). + 2 dtY(t) + 2 <302 Z(t) + 2 dt2J U dn ~

Уравнение (5) представим в виде П

, t <T DY(t) + n,n I u2 dn + — I V02

д С

— I DrV4w + D V4wdr + Mwt+

< DY

t

П

Используя представление t

2

+pi<(x,y, 0,^ + pigjwdr I =°. (14) w2(x,y,t) = 2 iw(x,y,r)wT (x,y,r)dr

00

t

и неравенство Гельдера, получим

w2dП < Т / w2¿П I ¿т.

У < 6С2е

СТ

у и <

25С

Dт

2 ес1т,

Усиливая правую часть (18) с помощью (19) и условий теоремы, и полагая

где С2 = (2МС1 + 1)/2МС1. Возвращаясь в (19) и (18), получим

У (И) = М

4 ¿П | ¿т < М еС1Т,

VI¿П | ¿т +1 DтУ(И),

w2dП <

Т6С-

М

2 ес1т,

\Ур\2 ¿п <

25С'

Р1

2 ес1т.

,'2 ТР1д\

С = тах

приходим к неравенству

^ | У е-0'1* ¿И

011 ¿и (2 рдд

и2с1П

Теорема доказана.

2. Случай упругих колебаний. Пусть существует два отличных от нуля решения w1, р1 и w2, р2 системы (1) - (3) в случае т = 0, удовле-+творяющие разным начальным условиям

wl(x,y, 0) = w'1(x,y) wl,t(x,y, 0) = w'2(x,y),

W2(x,y, 0) = чл\(х,у) ^и2,*(х, у, 0) = w'2(x,y).

(20) Функции w0 = w\(x,y)—w12(x,y) и w0 = w2(x,y) — Неравенство (20) проинтегрируем по И от 0 до 'wl(x, У) уд°влетв°ряют следующему неравенству: Учитывая, что

+С1е-

с" ¿И

6

< 2Ме-

С^

dZ

¿И

0, 4 / u2dП ¿И .)

п

*=о

выводим

у е-0"1*1 + еdL ( 2 DZ(tl) + 1 р1 / u2dП | +

+С2 / е-С1*

С1*1 (1 DZ(tl) + 2Р11 и2^ + 2DZ(t) + 2и2¿П | <

М

D,

^\\w0\\2,п + у\И\2,п + у\\у(х,у,0)\\ 1 ,п < 6,

Докажем аналог теоремы 1 для упругих колебаний, заменив условие (4) на (21). Функции V = w1 — w2 и р = р1 — р2 удовлетворяют следующему уравнению:

Mwtt + DУ V = —р1 р* — Р1 gw

и уравнениям и условиям (6) - (10). Заметим, что решение (р, V) рассматриваемой задачи при т = 0 обладает свойствами (11) - (13). Последовательно умножая уравнение (22) на V*, интегрируя результат по Пд и совершая предельный переход по R с учетом условий (7) - (8), приходим к тождеству

с с

< У(0) — 2МС (е-С" — 1) < У(0) + 2МС Dтлr, . 6 /2МС1 + 1

= тУ (0)+2мс < ч -ММС-

Первое, третье и четвертое слагаемые левой части полученного неравенства неотрицательны. Поэтому

¿1(1 ^^ р^ / u2dП) <6 (

(2МС1 + 1

2МС1

Откуда выводим

26ТС2 2 26ТС2

Z(И) < —^г-2, I u2dП <

D

± (2 г«? ¿п+|/ ^

п

+2! У

пп

Интегрируя тождество (23) по И, получим

У w22(x,y,t)dП + р I Ур2(х,у,г,ИЦП+

+У(х,у,И) ¿П + У w2(x,y,t)dП

п

р\д

*=0

= M / (Х,У, °)dU+ р f 0)dtt+

п п

+ j/Y(x,y, 0) dn+ Pfjw2 (x,y, 0)dn. (24) пп Привлекая условия теоремы, из (24) получим

У w2(xx,y,t)dn < M5, j < 25,

п п

Jy(x,y,t) dn < 2S, Jw2(x,y,t)dn < 5.

Теорема доказана.

Заключение. Для решения рассмотренной системы уравнений динамики ледовой пластины и жидкости в канале доказана устойчивость по начальным данным. Рассмотрены случаи вязко-упругих и упругих колебаний.

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

1. Squire V., Hosking R., Kerr A., Langhorne P. Moving loads on ice. — Kluwer Academic Publishers, 1996.

2. Hydroelasticity in marine technology. Edited by S. Malenica, N. Vladimir and I.Senjanovic. VIDICI d.o.o., 2015.

3. Sturova I.V. Unsteady three-dimensional sources in deep water with an elastic cover and their applications // J. Fluid Mech. — V. 730.— 2013.

4. Жесткая В.Д., Джабраилов М.Р. Численное решение задачи о движении нагрузки по ледяному покрову с трещиной // ПМТФ. — V. 49. — 2008. — № 3.

5. Шишмарев К.А. Математические вопросы моделирования взаимодействия ледового покрова и гидроупругих волн // Известия Алтайского гос. ун-та. — 2015. — № 1/1(85). DOI: 10.14258/izvasu(2015)1.1-22.

6. Шишмарев К.А. Постановка задачи о вяз-коупругих колебаниях ледовой пластины в канале в результате движения нагрузки // Известия Алтайского гос. ун-та. — 2015. — № 1/2(85). DOI: 10.14258/izvasu(2015)1.2-35.

7. Matiushina A.A., Pogorelova A.V., Kozin V.M. Effect of Impact Load on the Ice Cover During the Landing of an Airplane // International Journal of Offshore and Polar Engineering. — V. 26. — 2016.— № 1.

8. Korobkin A., Khabakhpasheva T., Papin A. Waves propagating along a channel with ice cover // European Journal of Mechanics B/Fluids. — V. 47. — 2014.

9. Shishmarev K., Khabakhpasheva T., Korobkin A. The response of ice cover to a load moving along a frozen channel // Applied Ocean Research. — V. 59. — 2016.

10. Brocklehurst P., Korobkin A.A., Parau E.I. Interaction of hydro-elastic waves with a vertical wall // Journal Enginering Mathematic. — V. 68. — 2010.

11. Batyaev E.A., Khabakhpasheva T.I. Hydroelastic waves in channel with free ice cover. Fluid Dynamics, 2015, №6.

12. Ткачева Л.А. Колебания цилиндрического тела, погруженного в жидкость, при наличии ледяного покрова // ПМТФ. — Т. 56. — 2012, № 6.

13. Шишмарев К.А. Математическая модель взаимодействия ледового покрова и гидродинамического диполя в канале // Известия Алтайского гос. ун-та. — 2017. — № 1 (93). DOI: 10.14258/izvasu(2017)1-30.

14. Vaigant V.A., Papin A.A. On the uniqueness of the solution of the flow problem with a given vortex // JMathematical notes. — 2014. — V. 96(6).

15. Хлуднев А.М. Об изгибе упругой пластины с отслоившимся тонким жестким включением // Сиб. журн. индустр. матем., — 2011. — V. 4(1).

16. Neustroeva N.V., Lazarev N.P. Junction problem for Euler-Bernoulli and Timoshenko elastic beams // ib. Elektron. Mat. Izv., V. 3. — 2016. -№ 7.

17. Lu H., Sun L., Sun J. Existence of positive solutions to a non-positive elastic beam equation with both ends fixed // Boundary Value Problems. — 2012. — № 56.

18. Basson M., de Villiers M., van Rensburg N.F.J. Solvability of a Model for the Vibration of a Beam with a Damping Tip Body // Journal of Applied Mathematics. — 2014.