ВОЛНЫ КОНЕЧНОЙ АМПЛИТУДЫ В ОДНОРОДНОЙ ЖИДКОСТИ

advertisement
ПРИКЛАДНАЯ МЕХАНИКА И ТЕХНИЧЕСКАЯ ФИЗИКА. 2009. Т. 50, N-◦ 5
67
УДК 532.59:539.3:534.12
ВОЛНЫ КОНЕЧНОЙ АМПЛИТУДЫ В ОДНОРОДНОЙ ЖИДКОСТИ
С ПЛАВАЮЩЕЙ УПРУГОЙ ПЛАСТИНОЙ
А. Е. Букатов, А. А. Букатов
Морской гидрофизический институт НАН Украины, 99011 Севастополь, Украина
E-mail: newisland@list.ru
Методом многих масштабов получены уравнения для трех нелинейных приближений
волнового возмущения в однородной идеальной несжимаемой жидкости, покрытой тонкой упругой пластиной, учитывающие нелинейность ускорения вертикальных смещений пластины при ее изгибе. На основе этих уравнений построены асимптотические
разложения до величин третьего порядка малости для потенциала скорости движения
жидких частиц и возмущения поверхности пластина — жидкость (изгиба пластины),
формируемых бегущей периодической волной конечной амплитуды. Проведен анализ зависимости волновых характеристик от модуля упругости и толщины пластины, длины
и наклона волны начальной основной гармоники.
Ключевые слова: поверхностные волны, изгибно-гравитационные волны, колебания
плавающей пластины.
Введение. В линейной постановке исследование волновых возмущений в жидкости
с плавающей упругой пластиной проведено в работах [1–7]. Волны конечной амплитуды
в однородной жидкости с упругой ледяной пластиной без учета нелинейности ускорения
ее вертикальных смещений рассмотрены в [8, 9]. В работе [10] проведен анализ зависимости характеристик прогрессивных поверхностных волн конечной амплитуды от толщины
и нелинейности ускорения вертикальных колебаний абсолютно гибкой пластины (битого
льда). В настоящей работе с использованием метода многомасштабных асимптотических
разложений исследуется влияние плавающей упругой пластины на распространение периодических волн конечной амплитуды с учетом нелинейности ускорения вертикальных
смещений пластины, обусловленных ее изгибом.
Постановка задачи. Пусть на поверхности однородной идеальной несжимаемой
жидкости, заполняющей неограниченный бассейн постоянной глубины H, плавает тонкая упругая пластина. В предположении, что движение жидкости является потенциальным, а колебания пластины безотрывными, рассмотрим влияние упругости и толщины
пластины на распространение изгибно-гравитационных
волн конечной амплитуды. В без√
размерных переменных x = kx1 , z = kz1 , t = kg t1 (k — волновое число) задача сводится
к решению уравнения Лапласа
∆ϕ = 0,
−∞ < x < ∞, −H 6 z 6 ζ
(1)
для потенциала скорости ϕ(x, z, t) с граничными условиями на поверхности пластина —
жидкость (z = ζ)
∂ 4ζ
∂ h 1 ∂ϕ 2 ∂ϕ i
∂ϕ 1 h ∂ϕ 2 ∂ϕ 2 i
D1 k 4 4 + κk
−
+ζ −
+
+
=0
(2)
∂x
∂z 2 ∂x
∂t
∂t
2 ∂x
∂z
и на дне бассейна (z = −H)
∂ϕ
= 0,
(3)
∂z
ПРИКЛАДНАЯ МЕХАНИКА И ТЕХНИЧЕСКАЯ ФИЗИКА. 2009. Т. 50, N-◦ 5
68
а также с начальными условиями
∂ζ
= 0.
∂t
ζ = f (x),
(4)
Здесь
D
D1 = ,
ρg
Eh3
D=
,
12(1 − ν 2 )
κ=h
ρ1
,
ρ
E, h, ρ1 , ν — модуль нормальной упругости, толщина, плотность, коэффициент Пуассона
пластины соответственно; ρ — плотность жидкости; g — ускорение свободного падения.
Потенциал скорости ϕ(x, z, t) и возмущение поверхности пластина — жидкость ζ(x, t) связаны кинематическим условием
∂ζ
∂ζ ∂ϕ ∂ϕ
−
+
= 0.
(5)
∂t ∂x ∂x
∂z
В динамическом условии (2) выражение, содержащее множитель κ, описывает инерцию
вертикальных смещений пластины, при этом первое слагаемое обусловлено нелинейностью
ускорения ее вертикальных смещений при изгибе.
Уравнения для нелинейных приближений. В предположении справедливости
разложений
ζ = εζ0 ,
ϕ = εϕ0 ,
ζ0 = ζ1 + εζ2 + ε2 ζ3 + O(ε3 ),
f = εf0 ,
ϕ0 = ϕ1 + εϕ2 + ε2 ϕ3 + O(ε3 ),
f0 = f1 + εf2 + ε2 f3 + O(ε3 )
методом многих масштабов [11] задача (1)–(5) сводится к определению нелинейных приближений порядка εn из уравнений
∆ϕn = 0,
−∞ < x < ∞,
−H 6 z 6 0;
(6)
∂ 4 ζn
∂ 2 ϕn
∂ϕn
−
κk
+ ζn −
= Fn∗ ,
z = 0;
(7)
4
∂x
∂z ∂T0
∂T0
∂ζn ∂ϕn
+
= Ln ,
z = 0;
(8)
∂T0
∂z
∂ϕn
= 0,
z = −H;
(9)
∂z
∂ζn
ζn = fn (x),
= Gn ,
t = 0,
(10)
∂T0
где fn (x) — приближение порядка εn для начального смещения пластина — жидкость;
D1 k 4
T0 = t,
T1 = εt,
T2 = ε2 t,
Fn∗ = Fn + Fn0 ,
F1 = F10 = L1 = G1 = 0,
n = 1, 2, 3,
∂ 2 ϕ1
∂ϕ1
∂ζ1 ∂ϕ1 1 h ∂ϕ1 2 ∂ϕ1 2 i
F2 = ζ1
+
+
−
+
+ κkN,
∂T0 ∂z ∂T1 ∂T0 ∂z
2
∂x
∂z
N=
∂ζ1 ∂ 2 ϕ1
∂ 3 ϕ1
∂ 2 ϕ1
+
ζ
+
,
1
∂T0 ∂z 2
∂T0 ∂z 2 ∂T1 ∂z
F3 = ζ1 N1 +
L2 =
∂ζ1 ∂ϕ1
∂ 2 ϕ1
∂ζ1
− ζ1
−
,
2
∂x ∂x
∂z
∂T1
G2 = −
∂ζ1
,
∂T1
1 2 ∂ 3 ϕ1
∂ 2 ϕ1
∂ϕ2 ∂ϕ1 ∂ϕ1 ∂ϕ2 ∂ϕ1 ∂ϕ2
ζ1
+
ζ
+
+
−
−
+ N2 + κkN3 ,
2
2
2 ∂T0 ∂z
∂T0 ∂z ∂T1
∂T2
∂x ∂x
∂z ∂z
69
А. Е. Букатов, А. А. Букатов
N1 =
N2 =
∂ 2 ϕ1
∂ 2 ϕ2
∂ 2 ϕ1 ∂ϕ1 ∂ 2 ϕ1 ∂ϕ1
+
−
−
,
∂T1 ∂z ∂T0 ∂z ∂z ∂x ∂x
∂z 2 ∂z
∂ζ2 ∂ϕ1
∂ζ1 ∂ 2 ϕ1
∂ζ1 ∂ϕ1 ∂ζ1 ∂ϕ1 ∂ϕ1
∂ζ1 ∂ϕ2
+
+ ζ1
+
−
,
2
∂T0 ∂z
∂T0 ∂z
∂T0 ∂z
∂T1 ∂z
∂x ∂z ∂x
∂ 3 ϕ1
∂ζ1 ∂ 2 ϕ2 ∂ζ2 ∂ 2 ϕ1
∂ 2 ϕ2
∂ζ1 ∂ 2 ϕ1
∂ 2 ϕ1
1 2 ∂ 4 ϕ1
N3 = ζ1 N4 + ζ1
+ζ2
+
+
+
+
+
,
2 ∂T0 ∂z 3
∂T0 ∂z 2 ∂T0 ∂z 2 ∂T0 ∂z 2 ∂z ∂T1 ∂T1 ∂z 2 ∂z ∂T2
N4 =
L3 = ζ1 N5 −
∂ 3 ϕ2
∂ζ1 ∂ 3 ϕ1
∂ 3 ϕ1
+
+
,
∂T0 ∂z 2 ∂T0 ∂z 3
∂z 2 ∂T1
G3 = −
∂ζ1
∂ζ2
−
,
∂T2 ∂T1
∂ 2 ϕ1 ∂ζ2 ∂ϕ1 ∂ζ1 ∂ϕ2
∂ζ2
∂ζ1
∂ϕ1 ∂ζ1 2
1 2 ∂ 3 ϕ1
ζ1
−
ζ
+
+
−
−
+
,
2
2
∂z 3
∂z 2
∂x ∂x
∂x ∂x
∂T1 ∂T2
∂z ∂x
∂ζ1 ∂ 2 ϕ1
∂ 2 ϕ2
∂ 2 ϕ1 ∂ϕ1
0
−
,
F
=
−κk
,
2
∂x ∂z ∂x
∂z 2
∂x ∂z ∂x
h
∂ζ1 ∂ 2 ϕ1 i
∂ 2 ϕ1 ∂ϕ2 ∂ζ1 ∂ϕ1 ∂ϕ1 ∂ 2 ϕ2
0
+
+
+
,
F3 = −κk ζ1 N6 +
∂x ∂z ∂x
∂x ∂z
∂x ∂x ∂z
∂x ∂z 2
∂ 2 ϕ 2 ∂ϕ ∂ 3 ϕ
1
1
1
N6 =
+
.
∂x ∂z
∂x ∂x ∂z 2
Заметим, что наличие слагаемых F20 , F30 в правых частях динамических условий (7) для
второго (n = 2) и третьего (n = 3) приближений обусловлено учетом нелинейности ускорения вертикальных смещений пластины.
Выражения для потенциала скорости и возмущения поверхности пластина — жидкость. Уравнения (6)–(10) получены для общего случая неустановившихся
возмущений конечной амплитуды. В случае бегущих периодических волн решение соотношений (6)–(10) будем искать, задавая fn (x) в соответствующем виде. Для этого первое
приближение (n = 1) возмущения поверхности пластина — жидкость ζ1 представим в форме волны
N5 =
ζ1 = cos θ,
θ = x + τ T0 + β(T1 , T2 ),
(11)
направление движения которой противоположно направлению оси x. Тогда из кинематического условия (8) находим
∂ϕ1
= τ sin θ,
z = 0.
(12)
∂z
Удовлетворяя граничному условию (9) на дне бассейна, потенциал скорости ϕ1 представим
в виде
ϕ1 = b0 ch (z + H) sin θ.
(13)
Подставляя соотношение (13) в (12), получаем b0 = τ (sh H)−1 . Следовательно,
ϕ1 = b1 sin θ,
b1 = τ (sh H)−1 ch (z + H).
(14)
С учетом (11), (14) из динамического условия (7) находим дисперсионное соотношение
τ 2 = (1 + D1 k 4 )(1 + κk th kH)−1 th kH.
Подставляя (11), (14) в правые части выражений (7), (8), при условии отсутствия основной
гармоники получаем
ζ2 = a2 cos 2θ,
ϕ2 = b2 sin 2θ,
(15)
ПРИКЛАДНАЯ МЕХАНИКА И ТЕХНИЧЕСКАЯ ФИЗИКА. 2009. Т. 50, N-◦ 5
70
где
a2 = 3τ 2 η2 /(4µ2 th H),
b2 = τ ν2 ch 2(z + H)/(4µ2 ch 2H th H),
ν2 = τ 2 (5 th H − cth H + 2κk) − 2(1 + 16D1 k 4 ),
η2 = (th H − cth H − 2κk) th 2H,
µ2 = (1 + 16D1 k 4 ) th 2H − 2τ 2 (1 + 2κk th 2H).
Подставляя решение (11), (14) для первого приближения и решение (15) для второго приближения в правые части динамического (7) и кинематического (8) условий задачи для
третьего (n = 3) приближения и исключая из них слагаемые, порождающие секулярность,
находим
ζ3 = a3 cos 3θ,
a3 = τ 2 η3 /µ3 ,
b3 = τ ν3 ch 3(z + H)/(3µ3 ch 3H),
η3 = (l1 − 3κkl2 ) th 3H − l2 ,
l1 = κkl11 + l12 ,
ϕ3 = b3 sin 3θ,
ν3 = 3τ 2 l1 − l2 (1 + 81D1 k 4 ),
l2 = (1/2)a2 (3 cth H + 6 cth 2H) − (3/2) cth H cth 2H + 5/8,
l11 = a2 (5 cth 2H − (1/2) cth H) + (1/2)(cth H − 5 cth 2H) cth H − 1/8,
l12 = a2 (11/2 − cth H cth 2H) + ((1/2) cth H cth 2H − 15/8) cth H,
µ3 = ( 1 + 81D1 k 4 ) th 3H − 3τ 2 (1 + 3κk th 3H),
σ0 = (1/2)[l3 − l4 (κk + cth H)−1 ],
β = τ σ0 T2 ,
l3 = (1/2)a2 (cth H + 2 cth 2H) − (1/2) cth H cth 2H − 3/8,
l4 = κkl41 − l42 ,
l41 = a2 (cth 2H − (5/2) cth H) + (1/2)(cth H − cth 2H) cth H − 3/8,
l42 = a2 (1/2 + cth H cth 2H) − (1/2)(cth H cth 2H − 5/4) cth H.
В результате возмущение поверхности пластина — жидкость ζ и потенциал скорости
движения жидкости ϕ до величин третьего порядка малости определяются из выражений
ζ=
3
X
n
ε an cos nθ,
ϕ=
n=1
θ = x + σt,
3
X
εn bn sin nθ,
n=1
2
σ = τ (1 + ε σ0 ),
a1 = 1.
√
В размерных переменных (ζ = ζ/k, ϕ = ϕ kg/k 2 , ε = ak (a — амплитуда начальной
гармоники)) имеем
ζ = a cos θ + a2 ka2 cos 2θ + a3 k 2 a3 cos 3θ,
p
p
p
ϕ = a g/k b1 sin θ + a2 kg b2 sin 2θ + a3 k kg b3 sin 3θ,
(16)
где
θ = kx + σ1 (1 + σ 0 )t,
σ 0 = a2 k 2 σ0 ,
σ1 = τ
p
kg,
τ 2 = (1 + D1 k 4 )(1 + κk th kH)−1 th kH.
Фазовую скорость волновых возмущений определим из формулы
p
v = τ g/k (1 + σ 0 ).
В приближении коротких волн (kH 1) решение упрощается, так как коэффициенты
в (16) имеют вид
r
3κk
1
1 + D1 k 4
a2 =
(1 + D1 k 4 ),
b2 = −
e2kz
[1 − (14 + 15κk)D1 k 4 ],
2∆2
2∆2
1 + κk
71
А. Е. Букатов, А. А. Букатов
a3 = [1 + 2κk + 24(κk)2 − AD1 k 4 − A1 (D1 k 4 )2 ]/(4∆2 ∆3 ),
b3 = e3kz
q
(1 + D1 k 4 )/(1 + κk) [13 + 7κk + 12(κk)2 − BD1 k 4 + 2B1 (D1 k 4 )2 ]/(24∆2 ∆3 ),
A = 13 − 4κk − 63(κk)2 ,
A1 = 14 − 2κk − 39(κk)2 ,
B = 499 − 529(κk) − 1014(κk)2 ,
B1 = 1899 + 4461κk + 2601(κk)2 ,
∆2 = 1 + 3κk − 2(7 + 6κk)D1 k 4 ,
∆3 = 1 + 4κk − 3(13 + 12κk)D1 k 4 ,
σ0 = −{3 − 7κk − 12(κk)2 − 2[21 + 41κk + 21(κk)2 ]D1 k 4 }/[8(1 + κk)∆2 ].
Анализ результатов. Решение (16) справедливо вне малых областей резонансных
значений волнового числа k = k1 , k = k2 , являющихся положительными действительными
корнями уравнений µ2 = 0 и µ3 = 0 соответственно. Данные уравнения сохраняют вид
и в случае пренебрежения нелинейностью ускорения вертикальных смещений пластины.
При kH 1 p
значения k1 , k2 удовлетворяют
уравнениям ∆2 = 0, ∆3 = 0. Отсюда находим
p
4
4
оценки k1 < 1/(4D1 ), k2 < 1/(9D1 ), k2 < k1 .
На рис. 1 показаны зависимости резонансных значений волнового числа от толщины
ледяной пластины при H = 100 м (рис. 1,а) и от глубины слоя жидкости при h = 0,5 м
(рис. 1,б), полученные при ρ1 /ρ = 0,87, ν = 0,34. На рис. 1 видно, что при фиксированной глубине слоя жидкости с уменьшением цилиндрической жесткости пластины значения k1 , k2 увеличиваются, стремясь к бесконечности в жидкости с открытой поверхностью (h = 0) и в жидкости, покрытой абсолютно гибкой (h 6= 0, E = 0) пластиной (битым
льдом). Влияние глубины слоя жидкости выражается в увеличении значений k1 , k2 . Однако
при уменьшении жесткости пластины глубина H, при которой это влияние существенно,
также уменьшается.
Для оценки влияния толщины и модуля упругости плавающей ледяной пластины на
амплитудно-фазовые характеристики волнового возмущения проведены численные расчеты при тех же значениях параметров E, ν, ρ1 /ρ, что и для рис. 1.
Анализ результатов расчетов показывает, что характер воздействия цилиндрической
жесткости ледяной пластины на структуру возмущений определяется глубиной слоя жидкости H, длиной λ = 2π/k и наклоном ε = ak волны начальной основной гармоники
à
k, ì -1
1,0
k, ì -1
0,12
á
2
0,08
0,5
1
0,04
2
1
0
0,5
1,0 h, ì
0
50
Рис. 1. Зависимости резонансных значений волнового числа от толщины пластины (а) и глубины бассейна (б) при ρ1 /ρ = 0,87, ν = 0,34:
сплошные линии — k2 , штриховые — k1 ; 1 — E = 3 · 109 Н/м2 , 2 — E = 5 · 108 Н/м2
100 H, ì
ПРИКЛАДНАЯ МЕХАНИКА И ТЕХНИЧЕСКАЯ ФИЗИКА. 2009. Т. 50, N-◦ 5
72
à
á
z/a
z/a
2
1
1
1
1
0
0
2
_1
_1
_2
0
100
200
300 x, ì
2
0
20
40
60 x, ì
Рис. 2. Профили возмущения поверхности пластина — жидкость вдоль оси x:
а — H = 10 м, h = 0,7 м, k = 0,04 м−1 , t = 150 c, ε = 0,08, б — H = 100 м, h = 0,15 м, k =
0,165 м−1 , t = 70 c, ε = 0,33; сплошные линии — без учета нелинейности вертикальных
смещений пластины, штриховые — с учетом нелинейности вертикальных смещений
пластины; 1 — E = 3 · 109 Н/м2 , 2 — E = 0
v, ì/ñ
300
js 0j .103
4
200
2
3
2
100
1
1
0
0,2
0,4
0,6
Рис. 3
0,8
1,0 k, ì -1
0
0,2
0,4
0,6
0,8
Рис. 4
Рис. 3. Зависимость величины сдвига фазы колебаний от волнового числа при
h = 0,5 м:
сплошная линия — без учета нелинейности вертикальных смещений пластины, штриховая — с учетом нелинейности вертикальных смещений пластины
Рис. 4. Зависимость фазовой скорости от волнового числа при различной толщине пластины:
сплошные линии — без учета нелинейности вертикальных смещений пластины, штриховые — с учетом нелинейности вертикальных смещений пластины; штрихпунктирные
линии — основная линейная гармоника; 1 — h = 0,5 м, 2 — h = 1,0 м
1,0 k, ì -1
73
А. Е. Букатов, А. А. Букатов
(рис. 2). Видно, что упругость пластины может оказывать влияние не только на амплитуду колебаний, уменьшая ее, но и на форму волнового профиля. При этом в коротковолновом
диапазоне (см. рис. 2,б) при увеличении упругости пластины влияние высших гармоник
также увеличивается. Кроме того, учет нелинейности ускорения вертикальных смещений
приводит к сдвигу фазы колебаний в направлении движения волны, в случае если пластина абсолютно гибкая (E = 0), и отставанию фазы колебаний в случае упругой пластины.
Количественную оценку изменения фазы колебаний за счет учета нелинейности ускорения
позволяют получить рис. 3, 4, где показаны распределения величины сдвига |σ 0 | фазы колебаний и фазовой скорости v соответственно (E = 3 · 109 H/м2 ). Следует отметить, что
кривые на рис. 3, 4 приведены только для участка изгибной ветви дисперсионной зависимости фазовой скорости. Влияние нелинейности ускорения в области гравитационной ветви,
где сила тяжести больше силы упругости, качественно такое же, как и в случае абсолютно
гибкой пластины (битый лед), рассмотренном в работе [10]. На изгибно-гравитационном
участке дисперсионной кривой фазовая скорость может как увеличиваться, так и уменьшаться за счет нелинейности вертикальных ускорений.
Заключение. На основе уравнений динамики нелинейных волн в однородной идеальной несжимаемой жидкости с плавающей тонкой упругой пластиной методом многомасштабных асимптотических разложений выведены дифференциальные уравнения в
частных производных для определения трех нелинейных приближений решения задачи
о пространственно-временной эволюции произвольного начального возмущения поверхности пластина — жидкость (изгиба пластины) конечной амплитуды. В случае бегущих
периодических волн для потенциала скорости движения жидких частиц и отклонения поверхности пластина — жидкость от невозмущенного состояния построены решения этих
уравнений в виде асимптотических рядов до величин третьего порядка малости. Проведен
анализ зависимости амплитудно-фазовых характеристик волновых возмущений от модуля
упругости и толщины пластины, длины и наклона волны начальной основной линейной
гармоники. Показано, что упругость пластины обусловливает не только уменьшение амплитуды волны изгиба, но и изменение формы профиля волны вдоль направления ее перемещения. В коротковолновом диапазоне влияние упругости приводит к увеличению вклада
высших гармоник. При учете нелинейности ускорения вертикальных смещений упругой
пластины наблюдается отставание фазы колебания, в то время как в случае абсолютно
гибкой пластины имеет место обратное явление.
ЛИТЕРАТУРА
1. Хейсин Д. Е. Нестационарная задача о колебаниях бесконечной пластинки, плавающей на
поверхности идеальной жидкости // Изв. АН СССР. Механика и машиностроение. 1962. № 1.
С. 163–167.
2. Хейсин Д. Е. Динамика ледяного покрова. Л.: Гидрометеоиздат, 1967.
3. Wadhams P. Attenuation of swell by sea ice // J. Geophys. Res. 1973. V. 78, N 18. P. 3552–3563.
4. Букатов А. Е. Влияние продольного сжатия на неустановившиеся колебания упругой пластинки, плавающей на поверхности жидкости // Прикл. механика. 1981. Т. 17, № 1. С. 93–98.
5. Squire V. A. A theoretical, laboratory and field study of ice coupled waves // J. Geophys. Res.
1984. V. 89, N C5. P. 1079–1089.
6. Schulkes R. M. S. M., Sneyd A. D. Time-dependent response of floating ice to a steadily
moving load // J. Fluid Mech. 1988. V. 186. P. 25–46.
7. Daffy D. G. The response of floating ice to a moving, vibrating load // Cold Red. Sci. Tech.
1991. V. 20. P. 51–64.
74
ПРИКЛАДНАЯ МЕХАНИКА И ТЕХНИЧЕСКАЯ ФИЗИКА. 2009. Т. 50, N-◦ 5
8. Гольдштейн Р. В., Марченко А. В. О длинных волнах в системе ледяной покров — жидкость при наличии ледового сжатия // Электрофизические и физико-механические свойства
льда: Сб. науч. тр. Л.: Гидрометеоиздат, 1989. С. 188–205.
9. Гладун О. М., Федосенко В. С. Нелинейные установившиеся колебания упругой пластины, плавающей на поверхности жидкости конечной глубины // Изв. АН СССР. Механика
жидкости и газа. 1989. № 3. С. 146–154.
10. Букатов А. Е., Букатов А. А. Нелинейные поверхностные волны в бассейне с плавающим
битым льдом // Мор. гидрофиз. журн. 2002. № 5. С. 34–46.
11. Найфе А. Х. Методы возмущений. М.: Мир, 1976.
Поступила в редакцию 7/VIII 2008 г.,
в окончательном варианте — 31/X 2008 г.
Download