![](/user_photo/_userpic.png)
книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред
.pdf§ 4.2. Главные, квадратичные и линейные модели |
381 |
Соотношения (4.2.74) называют спектральным представлением линейных моделей Ап вязкоупругих сред.
Если ввести спектральное разложение (4.2.72) и для тензора h: h =
П,,ч
= Рек >то с учетом (4.2.73) неравенство (4.2.70) примет вид
а = \ |
т |
Л |
|
|
Е j t Rap(t)Ya(h)Yp(h) + |
Е JtR^ Y a ( h) < 0. |
(4.2.76) |
|
а,(3=\ |
а = т -\-1 |
|
Если значения спектральных линейных инвариантов Y^(h) принять равными нулю, то, в силу независимости спектральных инвариантов между собой, из (4.2.76) получим условие монотонного невозрастания спектральных функций релаксации:
а = 1, ..., п. |
(4.2.77) |
С помощью спектральных функций релаксации формулируют частные случаи линейных моделей Ап вязкоупругих сред. Так, для простейшей линейной модели Ап изотропных вязкоупругих сред полагают одно из двух спектральных ядер релаксации постоянным:
R\\(t) = R\\(0) = l\ + (2/3) /2 = const, dR \\/dt = 0. |
(4.2.78) |
С учетом ynp. 6 к § 4.2 и формул (4.2.46) это условие можно записать в виде связи ядер q\(t) и q2 (t)\
|
qi(t) = |
-(2 /3 ) g2(i), dra/dt = - q a{t). |
(4.2.79) |
|
Определяющие соотношения (4.2.74) в этом случае принимают вид |
||||
|
(п) |
|
(п) |
|
|
I l( T) = J R n (0)Il ( C e), |
|
||
< |
(п) |
t |
Q (п) |
(4 .2 .8 0 ) |
. |
dev Т = j j R 22(t - |
т) dev — С в(т) dr, |
|
|
|
о |
ат |
|
где обозначены ортопроекторы тензоров относительно полной ортогональной группы I, называемые девиапгорами:
dev |
(п) |
(п) |
1 (n) |
(n) |
(n) |
J (п) |
|
|
Се = Се - |
^ ( С ^ Е , |
dev Т |
= Т - |
- /j (T )E . |
(4 .2 .8 1 ) |
|||
Соотношения (4.2.80) и (4.2.81) также можно записать в виде |
|
|||||||
( п ) |
J |
|
( п ) |
|
|
( п ) |
|
|
|
D /4 - |
|
д С в ( т ) |
dr. |
(4.2.82) |
|||
Т |
= -Д ц(0)/1(С „)Е + J |
R22{t - |
т)dev ——— |
При гидростатическом сжатии, когда тензор напряжений Коши и градиент деформации шаровые:
(п) |
1 |
(п) |
Т = -рЕ , F = kB, С = -----— (kn~lu - 1)Е, |
dev С = 0, |
|
|
п —III |
|
382 |
Глава 4. Вязкоупругие среды |
поведение вязкоупругой среды, согласно простейшей модели (4.2.80), анало гично поведению чисто упругой среды:
(п) |
(п) |
Т= (J/3)i211(0)/1(C)E,
т.е. вязкоупругие деформации при гидростатическом сжатии в такой модели отсутствуют.
Под действием давления вплоть до высоких значений р многие твердые материалы действительно обладают такими свойствами, поэтому простейшую модель (4.2.80) достаточно широко применяют в МСС.
4.2.14. Экспоненциальные функции релаксации и дифференциальная форма определяющих соотношений
При аналитическом и численном решении задач удобно иметь аналитиче ский вид спектральных функций релаксации R ap(t). В предыдущем разделе было установлено, что следствием неравенства диссипации ^ 0 является свойство монотонного невозрастания спектральных функций R aa(t). Функ ции с таким свойством можно аппроксимировать суммой экспонент:
N |
|
Ray t ) = R™p + Y ^ В рр{ exp ( - t/r c$ ), |
(4.2.83) |
7—1
где E>(T ) и (т) — константы, называемые спектрами значении релаксации и
времен релаксации; R ^ |
— предельное значение функций релаксации: |
||
|
lim Rap{t) = I Q , |
(4.2.84) |
|
|
t —ИЗО |
' |
|
которое может быть и равным нулю: R ^ |
= 0. |
|
|
Константы R ^ и |
удовлетворяют условию нормировки при t = 0: |
||
|
N |
|
|
|
K f , + Y , B«t! = c<-f- |
<4-2-85) |
|
О |
^ = 1 |
|
|
где Сар = Rap(0) — спектральные (двухиндексные) модули упругости при мгновенном нагружении.
Существуют и другие способы аналитической аппроксимации функций ре лаксации, однако экспоненциальные функции обладают рядом преимуществ: 1) выбирая достаточно большое число N экспонент в (4.2.83), можно аппрок симировать практически любую функцию Rap(t)\ 2) определяющие соотно шения (4.2.74) и (4.2.75) с экспоненциальными ядрами допускают обращение (см. п. 4.2.15), причем ядра обратных функционалов также оказываются экс поненциальными; 3) ядра (4.2.83) позволяют представить определяющие со отношения (4.2.74), (4.2.75) или (4.2.59), (4.2.60) в виде дифференциальных соотношений.
384 |
Глава 4. Вязкоупругие среды |
|
|
|||
|
t t |
|
|
|
|
|
|
а =Ет + 1 ..Q Q K aa(2t - |
( п ) |
|
( п ) |
|
|
7 |
т\ - T2)dCe(Ti) |
с1Св(т2). (4.2.92) |
||||
Здесь YaC\ r ) |
— линейные |
спектральные |
инварианты |
тензора |
( п ) |
|
С д(т) (см. |
||||||
(2.5.60)): |
|
|
( п ) |
|
|
|
|
Г р ( т ) = - |
|
|
|
(4.2.93) |
|
|
а п |
се(т), |
а = 1, . . . , |
т . |
Если в формулу (4.2.92) подставить экспоненциальные ядра (4.2.87) и
преобразовать двойные интегралы следующим образом:
t t
exp (- t - T \ - T 2 ) d Y P ( r x)d Y P (T 2) = 1~а{3
0 0
exp { - —t ^ - ) d Y i C\ T\) |
exp |
t~ T 2> riP) i |
{ - —t ^ - ) d Y f< \ T 2) = |
||
Ta[3 |
|
7~af3 |
1
exp { - 1— ^)ЯУ^с \т\)йт^ (Y^C\ t ) - d — exp ( - t~:^ j ) d Y ^ C\T 2)dT2J ,
ТД/3
(4.2.94) то с учетом обозначений (4.2.88) формулу (4.2.92) можно представить в виде
N |
гп |
Вос(3( т ) |
dW о( т(3) |
dW а(7(3) |
|
|
|
|
||
» 'Ч Е (Е |
' |
+ |
|
|
||||||
®а®(3 |
dt |
) (а/3 • • |
dt |
|
|
|
||||
7=1 |
а,/3=1 |
|
|
|
|
|
|
|
|
|
|
|
|
+ |
У |
в™ |
d W tt |
|
4Г . . (ПУЙ |
|
(4.2.95) |
|
|
|
dt |
|
— |
J |
||||
|
|
|
' |
/ j |
^ OLOL |
|
dt |
|
a = m -\-1
4.2.15. Обращение определяющих соотношений для линейных моделей вязкоупругих сред
В теории вязкоупругости часто используют обратные к (4.2.59) или (4.2.63) определяющие соотношения. Для этого соотношение (4.2.63) рас сматривают как линейное интегральное уравнение Вольтерры второго рода
( п )
относительно процесса С о(т), 0 < г < t. Ядро этого уравнения К (t) пола гают непрерывно-дифференцируемым, удовлетворяющим условиям (4.2.66), (4.2.67). Из теории интегральных уравнений известно (см. [32]), что уравне ние (4.2.63) с таким ядром имеет решение всегда, и это решение формально записывают в такой же форме, как и исходное уравнение:
\ |
, |
(п) |
|
(п) |
|
|
(п) |
т |
4N (t - |
т) ■■ у (г) d r , |
(4.2.96) |
||
С» = 4П |
" - + |
|||||
|
|
|
![](/html/65386/197/html_PhmaX2xzBv.Kjra/htmlconvd-3PPgnP385x1.jpg)
386 Глава 4. Вязкоупругие среды
|
4П |
4К(х) = 4N(X) • • 4М + |
4N (X —и) • -4К(г^) du, |
0 ^ и ^ |
х, |
|
|
о |
|
|
|
или |
снова |
возвращаясь к стандартным обозначениям аргументов |
х —►t, |
||
и |
т: |
|
t |
|
|
|
4П • • 4K(i) = 4N(i) • • 4М + |
4N (i —т) • -4К(т) dr, |
(4.2.101) |
О
получаем интегральное соотношение между тензором ядер релаксации 4К (t) и тензором ядер ползучести 4N(£). Если известно ядро 4К (£), то соотношение (4.2.101) представляет собой линейное интегральное уравнение Вольтерры второго рода для вычисления ядра 4N(£) и наоборот.
Подставляя в (4.2.101) спектральные представления (4.2.86), |
(4.2.98) тен- |
||||
зоров 4К (£), 4N(t) и аналогичные разложения для тензоров 4П |
|
о |
|||
и 4М, ввиду |
|||||
взаимной ортогональности тензоров 4Га , з.а |
(см. т. 1, § 4.10), получаем |
||||
т |
Р |
|
|
|
|
У У п apKpe{t) - N ap(t)Cpe - |
N ap(t - т)Крг (т)<1т) = 0 , |
а , г |
= |
1, ..., m; |
|
/?=1 |
о |
|
|
|
(4.2.102а) |
|
|
|
|
|
|
По/Q/i^Q/Q/(t) Naa(t)Caa |
Naa(t - r)K aa(t)dr = 0, a |
= rn + |
1, ..., n, |
||
|
|
|
|
|
(4.2.1026) |
— систему скалярных интегральных уравнений для нахождения ядер Nap(t) через ядра К ар и наоборот.
Подобно тензору функций релаксации 4 R(£), введем тензор функций
ползучести 4П(£), удовлетворяющий уравнению |
|
Jt 4n(t) = 4N (t), 4П(0) = 4П, |
(4.2.103) |
тогда обратное определяющее соотношение (4.2.96) можно представить в форме Больцмана:
|
(п) |
(п) |
|
( п ) |
т |
|
(4.2.104) |
|
4П ( t - т ) -- |
^ ( т ) . |
Тензоры ядер и функций ползучести обладают теми же свойствами сим метрии (4.2.67), (4.2.71), что и тензоры 4К и 4R (см. упр. 10 к § 4.2):
4П (i) = 4П (1243Д ) = 4П (2134Д ) = 4n (3412)(t) Vi > 0 . |
(4.2.105) |
Для тензора функций ползучести 4П (Ц как и для 4R(i), можно ввести спектральное представление по формуле (4.2.73):
§ 4.2. Главные, квадратичные и линейные модели |
387 |
4п(*) = £ |
п |
Y , |
n aa(t)4Tc |
(4.2.106) |
а,/3=1 |
Р |
а=т+1 |
|
|
где Паа(£) и ПаД£) — спектральные функции ползучести. |
|
|||
Теорем а 4 .2 .3 . Если |
спектральные |
ядра |
релаксации K ap{t) |
являются |
экспоненциальными, т. е. имеют вид (4.2.87), то и спектральные ядра ползучести Nap{t) являются экспоненциальными:
N А( i ) |
ехР (“ ТёТу)’ |
(4.2.107) |
Nap(t) = Y , |
||
7=1 Щ |
ск/3 |
|
и наоборот. |
|
|
(7) Лч) |
|
Константы А^р, tyap называют спектрами значений ползучести и времен |
|
ползучести. Они, вообще говоря, не совпадают с |
и r ^ j , однако число N |
в(4.2.87) и (4.2.106) одно и то же.
▼Покажем, что если ядра K ap(t) имеют вид (4.2.87), то решением интегрального уравнения (4.2.102) являются ядра (4.2.107).
Подставляя (4.2.87) и (4.2.107) в (4.2.102), находим
N |
Б (7) |
|
|
л(7) |
|
|
|
|
|
Яре |
ехр ( - * |
|
^а(3 |
ехр ( - Ai)ж |
)')J - |
|
|
Е(Е(п af3 Ai) |
|
: Ai) |
|
|||||
/3=1 7=1 |
Тде |
|
|
Ъа(3 |
|
Ъа(3 |
|
|
N |
N |
/ л(т)о(7') |
j |
|
|
|
|
|
|
:е |
( Аа13^13е |
|
ехр ( - Г - f)exp |
(-- |
- ) d T |
= 0 . (4.2.108) |
|
|
V ЛУЛУ) |
Q |
Al) |
|
Xi\ |
|
||
|
[ 7'=1 |
1а0~0е |
Za(3 |
|
'Ре |
|
||
Вычисляя в этом уравнении интеграл |
|
|
|
|
ехр (- |
t —Т |
|
,11% |
|
|
Ai) |
^ dT= |
(ехр |
_ехр |
(4-2-109) |
|
|
Ъа(3 |
T/3e |
La(3 ‘ f3e |
Tf3e |
La(3 |
и приравнивая коэффициенты в (4.2.108) при одинаковых экспонентах, полу чаем систему уравнений
гп |
N |
л(т') |
|
ГП |
|
N |
Б |
(7 ') |
4 (7) - О |
( Пд/3 |
|
Ла/3 |
В £{ |
= 0. £ 1 |
Ai) |
Е |
/Зе |
||
Е |
^ |
_ , ( V ) |
(l) |
—t(7) |
Aaf3 —U |
||||
/3=1 т/3е |
7'= 1Т/Зе |
Ъа(3 |
|
/3=1 |
Ъа(3 |
7'= 1 Че |
а/3 |
(4.2.110) |
|
|
|
|
|
|
|
|
|
|
|
для вычисления констант |
и |
через константы |
|
и rjA |
(константы |
||||
всегда могут быть вычислены через Сар и считаются известными). |
|||||||||
При тех |
значениях B ^J, |
при которых |
существует решение |
системы (4.2.110), существует и экспоненциальное представление ядер пол зучести (4.2.107). А
По формулам (4.2.110) можно вычислить константы |
и |
через B ^J, |
Tp7J и наоборот. |
|
|