книги / Устойчивость процессов деформирования вязкопластических тел
..pdfУСТОЙЧИВОСТЬ ВЯЗКОПЛАСТИЧЕСКИХ ТЕЧЕНИЙ |
91 |
Поскольку s может быть сколь угодно близким к нулю, ТО мажоранта правой части (11.23) совпадает с мажорантой выражения
2(ЕГ=о711фп(01)2/иХ)Г=о7121фп(012)- Нетрудно доказать, что для лю
бых вещественных чисел сь ...,с^ верно неравенство ( $ ^ =iC„)2 ^
^ z J n = ic»» причём равенство всегда реализуется при с\ = ••• = с#. В нашем же случае N = оо, и общая оценка (11.23) сверху невозможна.
11.4. Устойчивость точного решения задачи о разрыве танген циального напряжения на границе вязкопластической полуплоскости.
Точное решение данной задачи в области вязкопластического течения О < жз < £(t) записывается следующим образом
v
° |
' |
= |
' |
0 < , < , ь |
|
|
|
|
О |
|
|
|
|
®° = *3(1 - r s) - ^ = J c-<2d( - |
2y i e- ^ /(4i), |
(11.24) |
||||
|
|
О |
|
|
|
|
|
|
v |
|
|
~\ж’ |
|
V°'=l~Ts~^Jо е |
’ |
t,O(0,i) = |
||||
причём т)г находится из уравнения |
|
|
|
|
||
|
|
Vr |
|
|
|
|
|
|
j e-<2d ( = |
^ ( l - T s). |
|
(11.25) |
|
|
|
О |
|
|
|
|
Для того, |
чтобы |
вязкопластическое течение |
имело |
место при |
t > 0, необходимо, как видно, чтобы заданное напряжение на границе превышало предел текучести при сдвиге. Граница жёской зоны при этом имеет вид ху — 2т)гуД 9т. е. £(t) = 2г}г>Д.
Подставим параметры выписанного точного решения в оценку устойчивости (11.21) и получим окончательно
у. In (J? + s2J02K |
87Г2Га52(1 - Ts) |
|
1 - г , - 2s2 |
||
at |
7Г2 + 4T)}s2t |
|
|
(11.26) |
|
+ |
7Г2 |
|
t |
||
|
f > 2l*»(*)l2
92 ГЛАВА 3
Из-за невозможности промажорировать независимо от Фп ква дратную скобку в (11.26) представить общую достаточную оценку устойчивости не удаётся.
Классическое решение данной задачи для линейной вязкой жид кости следует из найденного выше при rs = 0, т]г — оо, £(t) = оо и моделирует течение в приповерхностном слое воды под действием ве тра постоянной силы [112]. Граничные условия в задаче устойчивости такого течения будут записываться уже не в виде (11.18), а в форме
х3 = 0 : ф = ф" = |
0, х3 = оо : ф = ф' - 0 . |
(11.27) |
Из (11.21), (11.26) для ньютоновской жидкости следует |
|
|
^ In (J? + s2J02) ^ |
sup |«°'| - 2s2 = 1- 2s2, |
(11.28) |
at |
o<z, |
|
и достаточными условиями устойчивости будут ограниченность снизу функции (2s2 - \)t при t > 0 и её стремление к +оо при t —►оо. Для возмущений с достаточно большими длинами волн эти условия не выполняются, коротковолновые же гармоники будут, наоборот, устойчивы. Аналогичные выводы о неустойчивом характере развития длинных волн и об устойчивости «мелкой ряби» имеются и в других приложениях, связанных с гидрогеодинамикой и геотектоникой.
§ 12. Вязкопластическое круговое течение Куэтта—Тейлора
Рассмотрим пример анализа 0303 в криволинейной ортогональ ной системе координат, а в качестве основного движения выберем несжимаемое вязкопластическое течение между двумя соосно вращаю щимися цилиндрами.
12.1. Невозмущённое движение и условия его существования. Обо значим радиусы внутреннего и внешнего цилиндров а и 6, а угловые скорости их вращения ша и щ соответственно. Характерной скоростью в данной задаче будет шаа, поэтому базис обезразмеривания фактиче
ски состоит из величин {р,о;а,а}. Тогда безразмерные параметры Re и к равны Re = puaa2/ii9 к = r$/(fiua)9 а сам предел текучести при сдвиге т$ следует отнести к величине pwla2. Введём также отношения
R = Ъ/а, ш = (шь - ша)/ша.
Единственная ненулевая компонента скорости v$(r) = v°(r) в невозмущённом стационарном течении Куэтта—Тейлора имеет вид
R2 |
/ |
| \ |
v |
— -(а ; ± x l n j R ) ^ r ------ |
j ^pxrlnr + r . |
( 12. 1) |
УСТОЙЧИВОСТЬ ВЯЗКОПЛАСТИЧЕСКИХ ТЕЧЕНИЙ |
93 |
На границах г = 1 и г = Я приняты условия прилипания, |
т. е. |
«°(1)= 1, г>°(Д) = Д(1+ы).
Другие характеристики невозмущённого движения являются сле дующими функциями радиуса г
v°re =
о
<ТГ0 —
Верхние знаки в формулах (12.1)-(12.3) следует брать, если ш > О, и нижние, если ш < 0.
Граница г = £, разделяющая область вязкопластического течения и жёсткую зону, которая вращается как абсолютно твёрдое целое, определяется из условия T°(r) = rs. Область Qf состоит из тех точек тела, где Т°(г) > т3 или U°(r) > 0.
В |
плоскости параметров |
(Я; |
2\и>\/х) |
выделим |
две |
области А и |
|||
Б (рис. |
10), |
разделяемые |
кривой |
2\ш\/х = R2 - |
1 |
- 2 In R. Из вида |
|||
U0(г) (12.2) |
следует, что, |
если |
исходная |
система |
с её |
геометрией и |
физико-механическими характеристиками принадлежит области А, то вязкопластическое течение имеет место при 1 < г < Я, и жёсткая зона отсутствует. Если же изображающая точка лежит в области Б, то появляется жёсткая зона, примыкающая к внешнему цилиндру. В этом случае ЕГ(£) = 0, т.е.
( 1 2 .4 )
Заметим, что £ ^ |
1 при Я > I, т. е. жёсткая зона не может зани |
мать весь зазор между |
цилиндрами (кроме, разумеется, случая о/ = 0). |
94 |
ГЛАВА 3 |
Рис. 10. |
Области А и Б на |
Рис. II. Зависимость {(Я) (12.4) |
плоскости |
(Я: 2\ш\/к). |
при различных значениях 2\ш\!х. |
Такое было бы возможным, если на одной из границ задавалось бы динамическое условие, например, удерживающий момент. На рис. 11 приведена зависимость £(R), отражаемая формулой (12.4), при различ ных значениях 2\ш\/х.
12.2. Постановка обобщённой задачи Орра—Зоммерфельда. Вы бранное в качестве невозмущённого известное течение (12.1)—(12.3) является одномерным сдвигом. В теории линейных вязких жидкостей, изучая устойчивость установившихся плоскопараллельных сдвигов, до статочно рассматривать только двумерные возмущения в плоскости основного движения. Для вязкопластических материалов при этом надо требовать дополнительное условие (см. обобщённую теорему Сквай ра 7.1), которое тем не менее допускает широкий класс возмущений, выходящих из плоскости основного движения. Как указано в § 7.2 (слу чай В), в криволинейной ортогональной системе трёхмерная картина вариаций вообще не сводится к двумерной. Поэтому анализ устойчи вости течения Куэтта—Тейлора справедлив лишь при искусственном ограничении vz = 0.
С учётом сказанного выше, исследуем двумерную картину вари аций скоростей (vr,ty,0). Линеаризуя определяющие соотношения и подста&яяя в них соотношения Стокса, получим
(Тгг — р Re \U° |
/ |
дт ’ |
гв |
|
Re \U° |
/ \ г |
дв |
дг |
г ) ’ |
|
|
, |
2 ( |
у. |
\ |
(1 |
dvt |
vr\ |
|
|
( 12.5) |
°в* |
Р |
Re \U° |
) |
( г |
дв |
г ) |
' |
|
||
|
|
УСТОЙЧИВОСТЬ ВЯЗКОПЛАСТИЧЕСКИХ ТЕЧЕНИЙ |
95 |
Введем функцию тока ^(г,0,£) |
|
1дф |
( 12.6) |
vr = г дв 1 |
тем самым удовлетворив условию несжимаемости. Подставим далее соотношения (12.5) в линеаризованные уравнения движения, исключим давление и стандартным путём, используя (12.6), придём к уравнению относительно <ф
4* / ф,т- ф /г\ _ ГдАф + |
о |
|
V |
||
V |
и° ),ввг~[ dt |
|
+ |
|
Г ( Д ^ ь - |
|
|
(12.7)
1
+
гг
К(12.7) надо добавить граничные условия при г = 1 и г = R, если исходное движение принадлежит области А на рис. 10, либо при г — 1 и г = £, если исходное движение принадлежит области Б. Рассмотрим сначала подробно первый случай:
г |
= 1 , г = Д : |
i/>r = il)$ = 0 . |
|
( 1 2 . 8 ) |
Разложим функцию тока в ряд Фурье по в |
|
|
||
ос |
|
|
|
|
i>(r, <М) = X |
<МГ) е" ‘в+а"‘ > |
« • = |
G С |
(12.9) |
и исследуем устойчивость отдельной гармоники возмущения с номером п. Критерием устойчивости будет условие а,и < 0 для любого п.
Подставляя ряд (12.9) в (12.7), (12.8) и опуская для краткости индекс п, получим обыкновенное дифференциальное уравнение
( 12. 10)
с граничными условиями
г = 1, г = R : ф = ф' = о. |
(12.11) |
Линеаризованная краевая задача (12.10), (12.11) представляет со бой 0303 для стационарного вязкопластического течения между со осно вращающимися цилиндрами.
96 |
ГЛАВА 3 |
12.3.Интегральные оценки устойчивости. Пусть комплексная
амплитуда ф(г) — элемент пространства |
с соответствующей |
нормой |
|
R |
|
№(<)Ц2 = / \ф"\2йхз. |
(12.12) |
1 |
|
Умножим обе части (12.10) на г, затем на комплексно-сопряжённую функцию ф и проинтегрируем от 1 до R . Используя однородные граничные условия (12.11), будем иметь
L = ~[а(12+ n2ll) + inQ]Re, |
(12.13) |
||||
где |
|
|
|
|
|
L = J{ + (In* + 1)J,2 + n2(n2- |
4)J() + 4xn(K* - |
К» + К 2) , |
|||
R |
ti |
|
R |
|
|
1 |
= |
f |
\Ф(т)\ |
|
|
i i =j \Ф{т)\2r2m~'dr, |
|
1 |
|
||
|
|
|
|
|
я
p
K l = J
Q ^ Q*
Q*=j
1Фlm)\2- ^ г d r , |
K 2 = |
|
\ф\2 |
(u°Yd r, |
« 2.m - 3, |
|
R |
|
|
|
|
-J |
|
r2(U0)2 |
|
|
1 |
|
|
R |
|
|
|
|
/» |
|
|
|
|
Q** > Q** — J и °(Ф’ф)*Лг, |
|
|||
2л.о |
|
|
|
|
n V |
+ и а(ф'фЪ - I x l l j dr. |
|||
+ — l^l |
Выделим в уравнении (12.13) действительные и мнимые части
nQtt - L/Re |
n Q ± |
(12.14) |
|
i 2 + n2i 2 ’ |
“ ** = - / 1 Т Щ - |
||
|
В силу неравенства Шварца [178] в Н2[1;Л] с нормой (12.12) имеем
1<2«1 < JК и°\ф'\\ф\йг = JИt7°v9|0'|Mrfr qixi0t |
( 12.15) |
УСТОЙЧИВОСТЬ ВЯЗКОПЛАСТИЧЕСКИХ ТЕЧЕНИЙ |
97 |
||
где |
а = |
max U°(r) . |
|
|
|
||
Воспользуемся |
также очевидным неравенством п!Д\ ^ |
(I 2 + |
|
п 21о ) / 2 . И з (12.14), |
(12.15) |
следует довольно общая интегральная |
оценка устойчивости основного течения. Сформулируем её в виде следующей теоремы.
Теорема 12.1. Пусть а ( п , Я е , к , а ; , Д ) |
— произвольное собственное |
число 0303 для течения Куэтта— Тейлора. Тогда, если |
|
2L |
(12.16) |
Re ^ |
q ( i ? + и% 2) ’
то а* < 0, w исходное течение устойчиво относительно двумерных возмущений (vr\ve\Q).
Оценка (12.16) носит достаточный характер. Дальнейший анализ основан на применении неравенств Фридрихса
Ц г |
_2 г2 |
_2 х2 |
2 |
(тг2 + 1п2Л)7„2 |
“ |
|
|||
J* > — - , X > |
(12.17) |
|||
|
In2 Л ’ |
In2 Л |
|
In2 Л |
доказательства которых приводятся ниже.
Доказательство /. Рассмотрим уравнение г0" + 0' + А20/г = 0 с граничными условиями 0(1) = 0(0) = 0. Умножим это уравнение на ф и проинтегрируем от 1 до R , получим А2 = i f //о - Точное решение 0(r) = А \cos(Alnr) + i42sin(Alnr) после подстановки в граничные условия даёт А = 1ж/ In R (I = 1,2,...). Таким образом, Am*n = ж2/ In2 R
и / ^ т г 2/02/1п2Д. |
■ |
|
|
|
|
|
|
|||
Доказательство 2. Рассмотрим уравнение гф1У + 2ф'" + А2ф"/г - |
||||||||||
А20 '/т2 |
= |
0 |
с граничными |
условиями |
(1 2 .1 1 ) . |
После |
умножения |
|||
его на |
0 |
и |
интегрирования |
от 1 |
до |
R с учётом |
(1 2 .1 1 ) будем |
|||
иметь |
А2 = J l / J l |
Подставляя точное решение |
ф(г) |
= |
А\ + А^г2 + |
|||||
j43rcos(Alnr) + j44rsin(Alnr) |
в граничные условия и приравнивая к |
|||||||||
нулю характеристический определитель, получим |
|
|
|
|||||||
|
|
|
I-cos(A lnfl) = ^ |
=-7^51п(Л1п Л), |
|
(12.18) |
||||
|
|
|
|
|
i t |
+ 1 |
|
|
|
|
либо |
|
|
|
|
|
|
|
|
|
|
|
|
|
I |
+ cos(AIn Л) = - ^-jr—7—sin(AIn Л ). |
|
(12.19) |
||||
|
|
|
|
|
it —1 |
|
|
|
|
98 ГЛАВА 3
Можно показать, что минимальный положительный корень урав нений (12.18) и (12.19) лежит в пределах 7r/lnR < Xmin < 2n/\nR, поэтому заведомо справедлива оценка j \ > я*2J 2/ In2 R. ■
Доказательство 3. Рассмотрим уравнение ф"/г - ф!/г2 + А2ф/гъ—
О с граничными условиями ф(1) = Ф\К) = 0. После выполнения аналогичных предыдущим операций получим А2 = J \/J Q. Точное же решение имеет вид
ф(г) = ^А\ ех р ^ \/1 - A2 In rj |
+ А2 е х р ^ - \/1- А^ ln r ^ jr , |
|
|
если |А| < |
1, |
ф(г) = (А\ + i42lnr)r, |
если |А| = |
1, |
ф(г) = ^А\ c o s ( v ^ T lnr^ |
+ А2 sin \ / А2 —1 ln r ^ jr , |
|
|
если |А| > |
1. |
Характеристические уравнения в каждом из трёх случаев в итоге сводятся к одному
tg(>/A2 - 1 1пя) = -л/А 2 - |
I. |
(12.20) |
||
Минимальный положительный корень Ат|П уравнения (12.20) |
||||
удовлетворяет неравенству |
A^in > (я*2 + In2 R)/\n2R. Поэтому |
J 2 > |
||
(я-2 + in2 а д 2/ In2 R. Ш |
|
|
|
|
Выпишем теперь следующие мажорантные оценки |
|
|||
J2 |
J 2 |
J 2 |
1 = 0,1,2, |
|
Rz |
q |
2Го2 < s— , |
(12.21) |
|
где |
s = |
min U°(r). |
|
|
|
|
|
Воспользуемся квадратичными неравенствами (12.17), (12.21) применительно к (12.16) и сформулируем результат в виде следствия
теоремы 12.1. |
Если |
|
|
Следствие 1. |
|
|
|
R2qRe |
4ж2п2х |
ir2 + n2 \n2R |
|
2 |
^ д(я*2 + п2 In2 R) |
In2R |
|
|
я*2 |
ж2 + \п2 R[\ - |
(12.22) |
|
4п2(\ + к/в)] |
||
ж2 + In2 R |
ж2 + п2 In2 R |
УСТОЙЧИВОСТЬ ВЯЗКОПЛАСТИЧЕСКИХ ТЕЧЕНИЙ |
99 |
то а* < 0.
Для получения общей нижней оценки критического числа Рей нольдса необходимо найти минимальное по п G Ъ значение правой ча сти (12.22). Исследуем аналитически некоторые частные случаи задачи, которые дают характерные качественные и количественные результаты.
12.4. Осесимметричные возмущения. При п —0 (осесимметрич ные возмущения) из (12.22) следует достаточное условие устойчивости
2(тг2 + 1п2Д)(Д2 - 1)
Я е£ |
(12.23) |
Д21п2Д[2Д2М + х (2 Д 21пД - Д2 + !)] |
|
Для изучения влияния пластической составляющей на линейно |
|
вязкое течение надо от к |
перейти к другому безразмерному числу т*, |
не содержащему параметр вязкости. Подставляя к — rsRe в (12.23), после решения квадратного неравенства получим
Яе< /г(Я )М |
, |
/|(Я )т, |
(12.24) |
||
|
f 2(R)u>2 |
||||
/з (Я) г, |
|
|
|
||
где |
|
|
|
|
|
/, (R) = 2 (тг2 + In2 Д)(Д2 - |
1)(2Д2 In R - R2 + 1), |
|
|
||
/ 2(Л) = Я31пЛ, |
/з(Д) = R(2R2\ n R - R 2 + 1) 1пД. |
|
|||
Оценка (12.24) позволяет получить асимптотику при малых т8 |
|||||
либо при больших относительных |
скоростях вращения |
(а;2 > |
т$). |
||
Удерживая первые три члена разложения Тейлора, будем иметь |
|
||||
Яе< |
(тг2 + In2 Д)(Д2 - 1) |
(12.25) |
|||
|
Д4 |w| In2 Д |
||||
|
|
|
|
||
(тг2 + In2 Д)2(Д2 - 1)2(2Д2 In Д - Д2 + 1) т3 |
|
|
|||
2Д10 |w|3 In4 Д |
|
|
|
||
Как видно из (12.25), |
при та —* 0 существует равномерная схо |
||||
димость по т4 к вязкому |
пределу. |
Устремляя T S / U 2 к |
нулю, |
мы, |
по-прежнему, остаёмся в области А (рис. 10), жёсткой зоны при этом не возникает.
В случае узкого зазора между цилиндрами, т.е. R - 1< |
(#+ 1)/2, |
достаточное условие устойчивости имеет простой вид |
|
27Г2 |
|
Я е ^ И ( я -1) ' |
( 1 2 .2 6 ) |
100 |
ГЛАВА 3 |
Предыдущие исследования проводились в предположении о том, что исходная система и траектории предельных переходов лежат в области А (рис. 10). Если же параметры системы принадлежат области Б, и тем самым зона течения сужается от 1 до £, то граничное условие прилипания при г — R (12.11) (либо (12.8)) следует изменить на условие при г = £, на поверхности жёсткой зоны в возмущённом движении.
12.5. Коротковолновые возмущения и вязкий предел. При п > 1 (коротковолновые возмущения) из общей нижней оценки (12.22) сле дует асимптотическая оценка
- R 2qRe ^ п2 + С , |
(12.27) |
где С — не зависящая от п функция x,u,R . Следовательно, гармоники возмущения с большими номерами всегда устойчивы24^, и критические числа Рейнольдса, им соответствующие, стремятся в бесконечность.
Покажем, что вязкий предел в данной задаче существует не только в случае осесимметричных возмущений, но и для любого п. При к —0 достаточное условие устойчивости (12.22) имеет вид
R 2q o n |
к 2 |
2 |
71,2 |
7Г2 + (1 - |
4 n 2) In 2 # |
2 |
In 2 R |
|
тг2+1п2Я |
|
(12.28) |
|
тг2 + 7121п2Д |
||||
где |
|
|
2R2 М |
|
|
|
|
|
|
||
|
|
|
90 ~ д ^ т р |
|
|
Формальный минимум правой части (12.28) достигается в точке |
|||||
|
п2 = а |
|
5а2+ 1 |
а |
7Г |
|
|
- |
(12.29) |
||
|
|
|
а2 + 1 |
|
\nR |
Можно показать, что правая часть (12.29) при любом а > 0 мень ше единицы, а при а2 > 2 + >/5 (R < 4,6016) и меньше нуля. Таким образом, минимум в (12.28) достигается при п —0 и совпадает с (12.23). Предельные переходы при т3 —►0 и п = 0 из самого общего случая вязкопластического течения оказываются перестановочными. Для уз кого зазора между коаксиальными цилиндрами остаётся справедливым условие устойчивости (12.26).
24^Напоминаем, что везде в этом параграфе речь идёт об устойчивости относитель но двумерной картины возмущений, т. е. фактически об условной устойчивости течения Куэтга—Тейлора.