книги / Устойчивость процессов деформирования вязкопластических тел
..pdfУСТОЙЧИВОСТЬ ВЯЗКОПЛАСТИЧЕСКИХ ТЕЧЕНИЙ |
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 } ^ |
7Г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)