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

книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред

.pdf
Скачиваний:
25
Добавлен:
19.11.2023
Размер:
41.7 Mб
Скачать

§ 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) в виде дифференциальных соотношений.

§ 4.2. Главные, квадратичные и линейные модели

383

Действительно, проведя цепочку подстановок: (4.2.83)^(4.2.73)^(4.2.62), найдем выражение для тензора ядер релаксации:

4К (*) =

K a/}(t)-

+

£

K aa{t)AT c

а,/3=1

 

&Q, ар

 

 

 

 

 

а = т + 1

 

( Л

 

d R a p{t)

= 2 ^

^ 7

ехР

 

K ap{t) = ------ ^

^

_(7)'

 

 

 

 

_(-у)

 

 

 

 

7=1 Тсч(3

 

'оД

Введем группу новых тензоров второго ранга:

 

t

 

(n)

 

 

 

 

— exp

-

t — Тч Сб>(т)^т

7 = Ь

TV.

УУсЩ—

 

r (7)

r(7)

 

 

 

т<х/3

1af3

 

 

 

(4.2.86)

(4.2.87)

(4.2.88)

 

Дифференцируя

 

no t и исключая интеграл, получаем, что

удовлетворяют дифференциальным уравнениям первого порядка:

 

 

d w S .

w g _

c e(t)

 

7 = 1, ..., TV.

(4.2.89)

 

dt

^

r(7)

_(7)

 

 

 

rl,

1а(3

 

 

 

 

 

 

1а(3

 

 

 

 

Подставляя (4.2.86) в (4.2.63) и учитывая выражения (4.2.87), (4.2.88),

получаем следующее представление определяющих соотношений:

 

 

 

(п)

О

(п)

 

N

(4.2.90)

 

 

Т

= J (4M • • С в -

W (7)),

 

 

 

 

 

7=1

 

где

спектры вязких напряжений,

являющиеся тензорами

второго

ранга:

 

 

 

 

 

 

 

WW = £ Bi7>

S

а а арч +

2

4 7TW ^

(4.2.91)

 

а,/3=1

 

 

 

а = т + 1

 

Таким образом, действительно с помощью экспоненциальных ядер (4.2.83) определяющие соотношения механически детерминированной модели Ап вяз­ коупругих сред (4.2.63) можно представить в дифференциальной форме (4.2.89)-(4.2.91). В этой форме при переходе от интегральных соотношений

к дифференциальным появляются дополнительные неизвестные W

( 7 )

ДЛЯ

а(3

которых формулируются уравнения (4.2.89).

 

 

 

При численном решении задач дифференциальная форма (4.2.89)-(4.2.91),

как правило, более удобна, чем интегральная форма (4.2.63).

 

 

 

Подставляя выражения (4.2.86), (4.2.62) в (4.2.65), для функции диссипа­

ции

получаем следующую формулу:

 

 

 

* _ J

t t

 

 

 

w

s r

K a/3(2t - n - r2) d Y P ( r l) d Y P ( r 2)+

 

 

 

= —

^

 

 

 

 

2

 

 

 

 

 

a,P=\ 0 0

 

 

 

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П

" - +

 

 

 

386 Глава 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)

 

( 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 и наоборот.

 

 

388

Глава 4. Вязкоупругие среды

Из (4.2.1026) получаем аналогичные несколько более простые формулы для вычисления констант А&1 и tad (а = т + 1 , ..., п)\

Пп

N

ЛЬ')

с п

N

в (7')

 

= Е

= V

(4.2.111)

 

 

Г У ) '

I(т)

— ^

_ Т-(т) *

г ( т )

(7 ) _

/у/—1

тл-(т')1

 

ry/—1

1аа

°аа

 

скск

‘"аа

 

Подставляя (4.2.98) и (4.2.107) в (4.2.103), находим выражение спектраль­ ных функций ползучести для экспоненциальных ядер:

N

П(xfi(t) = Ua/3 +

 

(l - ехр ( - 4 у ) ) ,

(4.2.112)

7=1

 

ск/3

 

причем

 

N

 

 

 

(4.2.113)

lirn^ Пay t ) =

+

= П ^ .

t^+oo

 

7=1

 

 

 

 

Упражнения к § 4.2

Упражнение 1. Используя правило дифференцирования интеграла с переменным верхним пределом (см. формулы (4.1.156), (4.1.25)), путем непосредственного вычис­

ления производной по t от функционала (4.2.30), показать, что из ОТТ (т. 2, (3.3.15))

(п)

действительно следуют формулы (4.2.31) и (4.2.32) для Т и гг*.

Упражнение 2. Используя определение (4.1.70), показать, что для линейных моде­ лей Ап термореологически простых вязкоупругих сред соотношения (4.2.30)-(4.2.32) имеют вид

t t

ф = ф0(в) + ^

фа0(ф - т[л’ - т/) dI(as)(Ti) dI{p{T 2 )+

а ,/3=10 0

+Х Ё VwO' -т[Л' -т'ф) dJ{as)(T\,T2),

Рa=ri+l о о

(п)

T = J

Е

ra0(f -

 

+ J

Y . raa(t' -

 

а,/3=1

Q

 

 

a=ri+l п

 

г\

t t

 

 

 

0

 

 

 

 

 

w* = ~ \

Е

 

-т'\3’ - т!ф) d l ^ ( r ^ d l ^ {т2)~

 

а,/3=1 QQ

 

 

 

 

 

 

 

 

 

г2 t

t

 

 

 

-

j

Е

Г.Фаае ' - т [ л ' - Л ) а л аз){т\,т2),

 

 

 

 

a=ri+l п

dt

где t', т[ и

определяются по (4.1.71).

 

0

0

 

 

 

Учитывая, что

 

 

 

 

 

 

d

l a ( r ) = d L I a ( T ) d T

=

71а(т')с1т' = сИа(т'),

показать, что эти соотношения можно представить как функции приведенного вре­ мени, в частности

 

§ 4.2. Главные, квадратичные и линейные модели

389

(п)

п

Г

T[)dI^{T[) + J j h

 

Т = J

Е

7 ’ ra0(t'

T\) d l(ac ( r i' ) .

 

а, /3=1

О

a = r i+l Q

 

Упражнение 3. Используя представления (4.2.56)-(4.2.58) для тензора функций ре­ лаксации 4 R (£), показать, что представления (4.2.64) и (4.2.39) для ф эквивалентны.

Упражнение 4. Подставляя в выражение (4.2.65) для формулы (4.2.56)-(4.2.58) для тензора функций релаксации 4R (t), показать, что формулы (4.2.65) в точности совпадают с (4.2.40).

Упражнение 5. Показать, что из условия (4.2.70) следует монотонное невозрастание следующих функций релаксации:

dRaaaa(t)/dt < 0, dRal3a0(t)/dt^ 0 V t^O , а,0 = 1,2,3; а ф /3,

где Rlikl(t) — компоненты тензора 4R(t) в базисе с

Упражнение 6. Используя представления (4.2.56)-(4.2.58) для тензора функций релаксации 4R (t) и его спектральное представление (4.2.73), а также учитывая вид тензоров аа и 4Га (см. п. 2.5.8 и т. 1, § 4.10), показать, что функции rap(t) и Rap(t) для изотропных сред связаны следующими соотношениями:

R\\{t) = rx(t) + (2/3 )r2(t), R22(t) = 2r2(t), а для трансверсально-изотропных сред —

Ru{t) = i?3333W = r22(i), R22(t) + i?33(t) = 2Rllll(t) = r„(i) + 2r44(t).

Упражнение 7. Показать, что для механически-детерминированных линейных моде­ лей Ап термореологически простых сред определяющие соотношения, представлен­ ные в упр. 2 к § 4.2, можно записать в форме (4.2.60), (4.2.64), (4.2.65):

t t

 

РФ = °рФо{0) +

\

( п )

 

 

( п )

 

 

d C e(n) • • 4Н(2Г -

Т[ -тф)-- dC e(T2),

 

 

 

о о

(")

д л

,

(п)

( п )

4R(i'- т') • • dCo(T),

w*=

Т = J

dCe(r,) • • ^ 7 4К ( 2 Г - г -

г ) • • dC e(T2).

о

 

 

 

о о

at

 

 

 

 

 

 

 

 

Показать, что определяющее соотношение (4.2.63) для этой модели имеет вид

(п)

t

 

 

 

 

 

 

О (п)

4К(7 - г') • •

( п )

 

 

4R (t)/dt'.

Т =

J (4M - • С e(t) -

C e(T)ae(T)dT), (tr) = - d

Упражнение 8. Показать, что для линейных моделей Ап термореологически простых сред с экспоненциальными ядрами определяющие соотношения из упр. 7 к § 4.2 можно представить в форме (4.2.89)-(4.2.91), (4.2.95)

(п)

 

 

 

,

о

 

(п)

N

Т = J(4M -- C e - ^ W ) ,

 

 

N

 

 

m

 

r (7)

7 = 1

Jag \

^ f \

Л

( n )

 

 

a(3

 

w = T

E

'E

 

T-4 7 CL a

( C 0 -

 

7 = 1

a ,/3=1

a /3

 

E

dW^l

ae(t)

( n )

+ I ^ ( w 3 w - c ^ ) ) = 0 ,

dt

1~af3

 

 

( n )

 

) <E>(Cff - W Д ) • ~a13

(7)

(n)

 

% ( c e - w ™ ) . . 4r 0 ( Ce - W ™ )

a=m+l

390

Глава 4. Вязкоупругие среды

Упражнение 9. Доказать, что из спектральных представлений (4.2.72) и (4.2.73) следует спектральное представление (4.2.74) для определяющих соотношений линей­ ных моделей Ап (2.59).

Упражнение 10. Показать, что тензоры ядер и функций ползучести 4N (t) и 4ПД) обладают теми же свойствами симметрии (4.2.67), (4.2.71), что и тензоры ядер и функций релаксации (t) и 4R (£), и наоборот.

§ 4 . 3 . М о д е л и н е с ж и м а е м ы х в я з к о у п р у г и х т в е р д ы х ср е д и в я з к о у п р у г и х ж и д к о с т е й

4.3.1. Модели несжимаемых вязкоупругих сред

Для несжимаемых вязкоупругих сред, как всегда, дополнительно суще­ ствует условие несжимаемости, которое можно записать в любой из форм (т. 2, (3.9.1)—(3.9.6)); для моделей Ап основное термодинамическое тождество принимает вид (т. 2, (3.9.9)), аналогичный и для моделей Вп.

Подставляя функционал (4.1.35) в формулу (т. 2, (3.9.9)), с использовани­ ем правила (4.1.27) получаем определяющие соотношения для несжимаемых вязкоупругих сред:

( Н

7) Н . о

Н

= —р5ф.

Поскольку для несжимаемых сред число г независимых инвариантов

(s) ^

 

 

Ц д С (£)) уменьшается на единицу по сравнению со сжимаемыми средами, то

в каждом представлении (4.1.53), (4.1.61), (4.1.73), (4.2.1),

(4.2.20) функцио­

нала свободной энергии у функции

(п)

пробегает (г —1)

С), в) индекс 7

значений.

Так же согласованным образом уменьшается и число z совместных инва­ риантов JjS\ участвующих в качестве аргументов ядер т .

4.3.2. Главные модели A n несжимаемых изотропных вязкоупругих сред

Для главных моделей Ап изотропных вязкоупругих сред функционал ф (4.2.1) зависит только от пяти совместных инвариантов

rn

m

(n)

m

(п)

j W = I a{ C e{t)),

4 + а = 1 °(с е(т))’ « = Г 2 ;

j \ I] = С в(т) ■• СЩ ), (4.3.2)

а определяющие соотношения (4.3.1) принимают вид

 

 

(n)

v (п)

(п)

 

т = ----- 7 - G - 1 +<р1Е + <р2Се,

(4.3.3)

п — 111