Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Устойчивость процессов деформирования вязкопластических тел

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
8.23 Mб
Скачать

УСТОЙЧИВОСТЬ ВЯЗКОПЛАСТИЧЕСКИХ ТЕЧЕНИЙ

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)

+

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 < s— ,

(12.21)

где

s =

min U°(r).

 

 

 

 

 

Воспользуемся квадратичными неравенствами (12.17), (12.21) применительно к (12.16) и сформулируем результат в виде следствия

теоремы 12.1.

Если

 

 

Следствие 1.

 

 

R2qRe

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п

где

 

 

2R2 М

 

 

 

 

 

 

 

 

90 ~ д ^ т р

 

Формальный минимум правой части (12.28) достигается в точке

 

п2 = а

 

5а2+ 1

а

 

 

-

(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^Напоминаем, что везде в этом параграфе речь идёт об устойчивости относитель­ но двумерной картины возмущений, т. е. фактически об условной устойчивости течения Куэтга—Тейлора.

Соседние файлы в папке книги