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

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

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

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

101

Итак, из всех гармоник п е N в вязком пределе наиболее опас­ ны осесимметричные, что подтверждается классическими результата­ ми [112, 242, 251, 293].

§ 13. Идеальножёсткопластическое течение Куэтта

Как показано в п. 7.5, необходимым условием справедливости общих оценок устойчивости, выведенных там, является неравенство Т# > 0. Это отсекает, например, идеальножёсткопластические течения (течения Сен-Венана). Однако для последних тоже можно сформули­ ровать линеаризованную краевую задачу устойчивости в возмущениях.

13.1. Постановка обобщённой задачи Орра—Зоммерфельда. Вос­ пользуемся уравнением (8.13) и выпишем его для идеальнопластиче­ ской среды. Для этого надо формально устремить число Рейнольдса

в бесконечность. В пределе Re оо будем иметь уравнение второго порядка

где безразмерный предел текучести т3 —x/Re отражает уже не отно­ шение вязких и пластических свойств среды, а характеризует лишь пластические свойства.

При т3 = 0 уравнение (13.1) сводится к хорошо изученному уравнению Рэлея. К соответствующему течению применимы теоремы из теории устойчивости идеальной жидкости — условия Рэлея и Фьортьофта—Хойланда, а также теорема Ховарда о полукруге [104] — связанные со знаком кривизны профиля скорости внутри области течения. Современный обзор вопросов о разложении по собственным функциям в задаче Рэлея, а также новые результаты, относящиеся к устойчивости плоскопараллельных течений идеальной жидкости, приведены в [197].

Известна возможная неединственность решения уравнений дви­ жения идеальнопластических тел [146]. Одному и тому же сдвиговому напряжению т3 соответствует, вообще говоря, бесконечно много зна­ чений v°l. Следовательно, для плоского течения Куэтта в качестве основного профиля можно выбрать любую функцию г;°(ж3), прини­ мающую заданные значения в точках 0 и 1. Граничные условия в возмущениях имеют при этом вид

ж3 = 0 : 0 = 0; ж3 = 1 : 0 = 0 .

( 1 3 .2 )

102

ГЛАВА 3

Будем полагать далее, что ос­ новное течение реализуется в обла­ сти 0 < Ж3 < 1 и характеризует­ ся произвольной монотонно воз­ растающей функцией г7°(ж3) с не­ прерывной производной такой, что

И< 9 и

 

 

 

/

dxз

(X), V х 3 € [0; 1].

 

 

 

<

 

 

 

К (*з)|

 

 

 

 

 

 

 

 

(13.3)

0

0,5

л3

 

Возможно

и

любое распре­

деление на отрезке 0 < ж3 < 1

Рис.

12. Чередование жёстких зон и про­

жёстких зон (рис. 12), а также вы­

слоек идеальнопластического течения

 

рожденный случай

отсутствия зон

 

 

 

течения, когда скорость деформации представляет собой сумму не­ скольких ^-функций25*. Эти случаи исключаются из рассмотрения.

13.2. Интегральные оценки устойчивости. Пусть ф — элемент комплекснозначного гильбертова пространства И2(П) с нормой (9.2), имеющий две непрерывные производные во всей области.

Умножим равенство (13.1) на комплексно-сопряжённую функцию ф и проинтегрируем по жз от 0 до 1. С учётом граничных условий (13.2) получим

4Tssll + — (1\ + S2IQ) + %Ji [(v0,/ + S2V°)|0 | 2 +

0

+ г/>|0/|2] dx^+i J1у0,ф!фйх\ = 0, (13.4)

о

где величины 1 т и 1 0 такие же, как и в задаче о вязкопластическом течении Куэтта (см. (9.3)).

'■^Известный пример, который обычно приводят, говоря о таком течении, это сдвиг пачки листов бумаги. Задавая скорости верхнего и нижнего листов, заведомо нельзя ничего сказать о профиле скорости внутри пачки, все они равновероятны. Сдвиг может проходить всего по одному листу, в этом сл\чае скорость деформации — сумма двух Ь-функции.

В силу своей неединственности основное течение в данной задаче уже в опре­ делённом смысле неустойчиво.

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

103

Приравняем действительную часть стоящего в левой части (13.4) выражения к нулю

4Tssl] + у (if + s2J02) - J у'ф'ф dx$ 0.

(13.5)

о

 

и воспользуемся неравенством Шварца (9.6) в пространстве И^(П) с нормой (9.2). Из (13.5), (9.6) следует

Теорема 13Л. Пусть a (s,rs) — произвольное собственное число краевой задачи (13.1), (13.2). Тогда

s(glpl\ -

4rssll)

а* ^

 

(13.6)

1

} + s2I 2

Достаточное условие устойчивости плоского идеальнопластиче­ ского течения Куэтта — отрицательность правой части неравен­ ства (13.6).

Следствие 1. Если

T s

7Г2 + S 2

q2 >

(13.7)

87Г252

то а* < 0.

 

Следствие 2. Если

 

q2

(13.8)

4irs

то а* < 0.

Доказательство обоих следствий вытекает из очевидного неравен­

ства S I QI I ^

(I 2 + s2ll)jl

и неравенств Фридрихса для квадратичных

функционалов 1 } ^

2ll,

i l >

l}/q.

 

 

 

Кривые 1 и 2, соответствующие следствиям I и 2, изображены

на рис. 13.

Так как

следствия

I и

2 независимы и достаточны, то

справедливо

 

 

 

 

 

 

 

 

 

Следствие 3.

Если при фиксированном s0

 

 

Ts ^

f

+ 4 .

|

\

I

(I3.9)

 

—г > min

8тг2^

4irs0

J

Vs» > 0

 

Г

 

4тгз« ’

 

mo a * (s n ,r4) < 0.

104 ГЛАВА 3

Это

следствие означает,

что,

если при некотором s0

значе­

ние тя/ц2 лежит выше кривой 2

на рис.

13,

то исходное

тече­

ние устойчиво относительно

возмущения с

волновым числом s0.

Ts / q 2

 

 

 

 

 

 

Поскольку в реальном возмуще-

 

 

 

 

 

 

нии присутствуют все гармоники

0,5

 

 

 

 

 

 

 

 

 

 

 

 

s > 0, то общей достаточной и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

независимой от s оценки устой­

 

 

 

 

 

 

 

чивости

представить не удаётся.

0,25

 

 

 

 

 

 

Из рис. 13 видно, что наибо­

V

N

 

 

 

лее неустойчивыми будут длин­

 

2

^

 

 

 

 

новолновые

вариации, увеличе­

 

 

 

 

 

ние же Ts/q 2 оказывает стабили­

0

 

 

1

2

 

S

зирующее влияние.

 

 

 

 

Дадим оценку минимальной

Рис. 13. Критические кривые 1 и 2, со­

неустойчивой длины волны Л* в

ответствующие

следствиям

1 и 2

теоремы

слое с характерными для геотек­

13.1, в плоскости (s:rs/q2)

 

 

 

тонических процессов параметра-

 

 

 

 

V ~

 

 

ми: Р ^

3000

кг/м ,

10 см/год (высокоскоростные процессы),

q ~ Ю,

, ~

107 Па (при температурах ниже

200°С) [176]. Согласно

(13.9) Л*

=

87T2rs/g2,

т. е. Л*

~ 3 • Ю20. Соответствующие размерные

длины волн превышают всякие линейные размеры Земли. Следователь­ но, ограничение (13.9) в геофизических и геотектонических задачах не оказывает существенного влияния на устойчивость идеальножесткопла­ стического течения в плоском слое.

13.3. Течение Куэтта в узком смысле и длинноволновое прибли­ жение. Рассмотрим подробнее краевую задачу (13.1), (13.2) в случае линейного профиля скоростей по толщине: v°(x3) = Ж3, 0 < Х3 < 1,

что представляет собой течение Куэтта в узком смысле. Уравнение (13.1) перепишется в виде

(а + isx3 + 4Tss2)<j>" - s2(a + isx})((> = 0,

(13.10)

a после замены независимого переменного ( = 2(s x $ -ia -2ic) сводится к уравнению Уиттекера

d2<b

4Сттг - (С + 4гс)0 = 0, с = 2T SS ~,

(13.11)

которое имеет фундаментальное решение

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

105

Ф(О = М_к;1/2(С )(с| 2 J f(v)dr,) ,

(13.12)

о.

1

 

f(v) =

 

/2(V)

 

где Со = < U , = o , C i = С к , = ь Af_jc:i / 2 ( C ) - Функция Уиттекера

[205J,

связанная с вырожденной гипергеометрической функцией Ф( 1+ гс, 2; С) соотношением

М-,с:1/2(С) = СеЧ/2Ф(1+*с,2;С) =

= <е-с/2 > + £

(1 +*с)„Сп

(13.13)

п! (п + 1)!

 

(•)„ — символ Похгаммера.

 

 

 

Подставим (13.12) в граничные условия

 

С - С о , '• Ф =

0 ;

С - C i '• 0 = 0 >

(13.14)

следующие из (13.2), и придём к характеристическому уравнению

 

Ci{s,a)

 

D(s,a) =

j

f(Q d ( = 0,

(13.15)

которое неявно описывает зависимость частоты а от волнового числа 5. Далее следует подставить (13.13) в (13.12) и (13.15) и получить уравнение, содержащее четырёхкратные степенные ряды по (о($,а) и £i(s,а). Аналитическое исследование в этом направлении, есте­ ственно, затруднено. Поэтому ниже остановимся на длинноволновом приближении в данной задаче.

Разложим в ряд Тейлора по s вблизи 5 = 0 функцию D(s,a) в соотношении (13.15) и удержим первые два ненулевые члена этого разложения. Так как Со(0,or) = (i(0,а), то D(0,а) = 0, и будем иметь

s d 2D , k

5 7 № . ) t - s r (M ) = «.

Вычислим частные производные от D по s при а

9D

ds (0,«) = 2f ( - 2ia ) ,

(13.16)

( 1 3 .1 7 )

d 2D

106

ГЛАВА 3

Подставляя (13.17) в (13.16) и учитывая (13.12), получим

Из определения функции Уиттекера (13.13) и величины c(s) (13ЛI) вытекает, что

В результате (13.18) сведётся к простому виду

tg а = is .

(13.19)

Корни уравнения (13.19) при \s\ < 1 записываются так

а = тг + г arctg s, п Е Ъ .

(13.20)

Наличие корней (13.20) с положительной действительной частью а* = 7Г71, п ^ 1 говорит о неустойчивом характере длинных (s < 1) волн. Об этом же свидетельствуют и достаточные интегральные оценки

вп. 13.2.

§14. Наследственно вязкопластическне сдвиговые течения

Среди многих моделей неныотоновских жидкостей особенно по­ пулярны в исследованиях релаксационные модели, в которых напря­ жения зависят от всей истории деформирования (полимерные или вязкоупругие жидкости с затухающей памятью [21, 224]). Ниже пред­ лагается скалярное определяющее соотношение, являющееся обобще­ нием линейной вязкоупругой жидкости и вязкопластического тела, построенного на основе модели Ильюшина—Бингама. Векторные со­ отношения между девиаторами напряжений и скоростей деформаций принимаются линейными. Таким материалам даётся название наслед­ ственно вязкопластические.

Одновременный учёт релаксации вязких свойств материала и на­ личие предела текучести позволяет в задачах о плоскопараллельном сдвиге описать изменение толщины жёсткой зоны. Заметим, что это изменение вызвано не переменными внешними данными (переменным

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

107

градиентом давления в течении Пуазейля, переменной скоростью гра­ ницы в течениях Куэтта и т.д.) [2, 11, 16, 25, 34, 84, 130, 131, 184-187, 240, 294], а именно явной зависимостью материальных функций от времени при постоянных внешних воздействиях. Приводится также постановка линеаризованной краевой задачи устойчивости одномер­ ного наследственно вязкопластического сдвига относительно малых возмущений.

14.1. Наследственно вязкопластическое течение Пуазейля в плос­ ком слое. Пусть максимальное касательное напряжение Т(х,t) и максимальная скорость скольжения U(x,t) связаны следующим обра­

зом

 

T = rsH(t) + (p - M )U .

(14.1)

Здесь H(t) — функция Хевисайда, М — интегральный

оператор

разностного типа

M U = Jt M (t-Q U (()d (

(14.2)

о

с регулярным ядром M(t) ^ 0, описывающий релаксацию вязкости. Рассматривается плоскопараллельное течение материала со ска­

лярными соотношениями (14.1), (14.2) в слое 0 < ж3 < ft, причём причиной сдвига является постоянный по времени перепад давления Vp по оси Oxj (обобщённое течение Пуазейля). В качестве бази­ са обезразмеривания проще всего выбрать тройку величин {Vp,p,ft}.

Роль числа Рейнольдса Re выполняет комбинация у/ph3Vp /р, без­ размерный же предел текучести равен rs/(hVp).

Единственное уравнение движения относительно продольной скорости v\(ж3,£) = v(x,t) в безразмерных переменных принимает следующий вид

Наследственно вязкопластическое течение реализуется в двух не свя­ занных между собой пристеночных слоях 0 < х < £(t) и 1-£(t) < х < 1, где £(t) — неизвестная граница, отделяющая зону течения от жёсткой прослойки.

Таким образом, на границах х = 0, х = £(t) и при х = 1/2 (из со­ ображений симметрии относительно оси х = 1/2) должны выполняться условия

ж = 0 : г7 = 0 ;

ж = £ ( £ ) : Т = т 5; ж = ^ : Г = 0 .

(14 .4 )

В качестве начального условия можно взять соответствующее стацио­

108

ГЛАВА 3

нарное решение в слое 0 < х < £() = 1/2 - т*, полученное при M(t) = 0. Одним из методов исследования краевой задачи (14.3), (14.4)

является применение преобразования Лапласа. Его использование в области с меняющейся со временем границей предполагает привле­ чение специальных методов и приёмов операционного исчисления26*. Таким приёмом может быть аналитическое продолжение функции в более широкую область, включающую исходную, и решение задачи в изображениях в ней. Физически это означает некоторое доопределе­ ние скорости внутри жёсткой прослойки с последующим выполнением граничного условия на подвижной границе.

Отметим здесь ожидаемое (следующее из визуальных наблюде­ ний) изменение решения v(x,t) со временем. На рис. 14 приведены схематичные профили продольной скорости в несколько последова­ тельных моментов времени.

Рис. 14. Последовательные во времени схематичные профили скорости v(x,t) в наследственно вязкопластическом течении Пуазейля

14.2. Наследственно вязкопластическое течение Куэтта. В случа плоскопараллельного течения Куэтта с заданной скоростью верхней границы V в безразмерный базис естественно включить величины

{р) V, h}. Дальнейшие рассуждения в постановочном плане аналогич­ ны проведённым в п. 14.1. Существенным же отличием здесь является отсутствие подвижной границы жёсткой зоны, что упрощает примене­ ние преобразования Лапласа

ОС'

 

/.(« ) = I f(t)c -atdt.

(14.5)

о

 

26*Особенности этих методов отметил И. А. Кийко, указавший также на невозмож­ ность непосредственного применения преобразования Лапласа к уравнению, описываю­ щему толщину жёсткой зоны х = ((t).

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

109

Единственное уравнение движения

d~v dv

(14.6)

дх2 dt

с граничными условиями х = 0 : г; = 0, ж =

1: г; =

1 в изображениях

выглядит следующим образом

 

 

 

 

v'l - a Re+v+ = 0,

г;^(0) = 0;

г;*( 1) =

(14.7)

и имеет решение

 

 

 

 

 

 

 

(14.8)

Длительное число Рейнольдса Яе*(а) связано с мгновенным Re

соотношением

1/Яе*(а) = 1/Я е-

М*(а). В

силу неотрицательности

функции М(t)

справедливо неравенство Яе*(а) > Re.

 

Исследуем далее вязкопластический предел в данной задаче. Если наследственные свойства у материала отсутствуют, т. е. М = 0, и, следо­ вательно, Л = \/аЯе, то из (14.8) можно получить известное решение в пространстве оригиналов задачи о вязкопластическом течении Куэтта. Устойчивость этого течения исследовалась ранее в §9.

Применение обратного преобразования Лапласа (преобразования Мелина) к изображению (14.8) позволяет перейти в пространство оригиналов и найти функцию v(x,t). Рассмотрим два характерных вида ядер М(£) — экспоненциальное и степенное.

а) Экспоненциальное ядро, релаксирующее до нуля М(£) =

e~t/n/(nRe). Здесь п — время релаксации. В данном случае

М*(а) =

[(1 +па)Я е]“1 и согласно (14.8) А(а) = у/(\/п + a)Re.

 

б) Степенное ядро M(t) = A/[(t + t{))2Re]\ A,t = const,

А Ф to.

Изображение M* и собственное число А имеют вид

 

 

(14.9)

atoRe

(14.10)

А =

to - А[ \ + ato ехр(а^о) E'x(-ato)]

 

где Е1(ж) — интегральная показательная функция.

 

п о

ГЛАВА 3

Из выражения (14.8) видно, что, если для некоторого ядра ре­ лаксации предельное при а —►О значение А равно нулю (напри­ мер, для степенного ядра), то пре­ дельный профиль, если он суще­ ствует, представляет собой линей­ ную функцию. Для экспоненци­ ального ядра, релаксирующего до нуля, как было найдено выше, А0 = у^йе/п. Предельные профи­ ли скоростей для некоторых А0

Рис. 15. Предельные профили скоро­

 

 

 

сти

в наследственно

вязкопластическом

lim v{x, t) = sh ( у

(14.11)

течении

Куэтта

для

экспоненциально­

го

ядра

при

различных До, причём

<—oo

sh Ao

 

 

 

 

Д()1 < ^02 < ^03 < ^>4

изображены на рис. 15. Видно, что, если А0 > 1 (время релаксации мало), то почти весь слой за исключе­ нием тонкой пристеночной области практически не движется. Если же Ао < 1 или Ao « 1, то определяющую роль играют мгновенные свойства материала.

Полный анализ длительного поведения материала в пространстве оригиналов может быть основан на разложении изображения (14.8) в ряд по а [117].

14.3. Постановка линеаризованной задачи устойчивости. Приведё постановку линеаризованной краевой задачи устойчивости одномерно­ го наследственно вязкопластического сдвига. Подставим в уравнения движения в возмущениях (6.11) проварьированные с учётом (14.1) векторные определяющие соотношения

2Т°

Р = 1,3,

 

Spp- — vpp,

 

•13 = «13 - 2 ft M(t -

<)»,з(0 dC,

(14.12)

О

 

 

и далее стандартным путём придём к одному уравнению относительно функции тока ^(жь ж3,t)

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