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

книги / Основы математической теории термовязкоупругости

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

Аналогично

выражается

и функция рассеяния (7.6)

 

IV =

г г

I

 

 

(т^) 3&т(т2),

 

^^

^2)

 

(8.14)

где

 

 

N

т

т

 

 

 

 

 

 

 

<2 (*-

*1, * - ъ)

=

2 1 г

2

2

 

=

 

 

 

;=1 ^

Ь=1г=1

 

 

 

 

 

N

 

 

 

 

 

 

2

 

т1)

— ^2)*

(8.15)

 

 

 

У=1

 

 

 

 

Выражение энтропии получим, дифференцируя с обратным знаком выра­ жение (7.5) по температуре и используя соотношения (8.11)

 

Л “

9Т - )

 

То

ат +

 

 

 

 

 

 

 

 

 

 

 

 

<

 

 

с

I

 

 

 

 

 

+ а^Я8 (I — х) <1ет (х) =

С.

а ^Я8 (I — т) йгт(х),

(8.16)

 

 

где

 

М

 

7/*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

я 8 {I -

г) = 2

г

{ (2 ан)

2

 

с = - т 0- ^ - . (8.17)

 

 

г=1

г=1

** -

 

 

 

 

 

 

 

 

 

 

к = 1

 

 

 

 

 

 

Для явного выражения связи (8.4) между напряжением и деформацией

воспользуемся

разложением рациональной функции

М (к)

=

Г) (Я)/Р (Я)

на

элементарные дроби

 

 

 

 

 

 

 

 

 

Т > ( 4

^1___ I

^2

I

|

 

 

 

(8.18)

 

К — %1

 

Х - Х г

 

' '1~ X — Ят

 

 

 

 

 

 

 

 

 

после чего получим

 

 

 

 

 

 

 

 

 

° в 4 ( а у ет “ м

(й)ет=

 

 

 

(

 

(8.19)

 

 

С1

 

С2

 

 

.

 

 

 

 

 

 

_*х

а~4ет

(* ~ х)Лгт00>

 

1 — Я1с!"1

 

%2(1

 

 

 

 

 

 

 

т

'

 

 

 

где

ядро релаксации .й (I

— т) является суммой экспонент

 

 

 

 

т

 

 

 

 

 

 

 

 

 

я (* - -г) = 2

 

т).

 

 

 

 

 

(8.20)

 

Подставим теперь в (7.14) полученные выражения Т*, 5,

 

а. Тогда

получим условия, накладываемые на коэффициенты, в термодинамичес­ ких функциях

м

N

Ь* Ъ1*

М

 

2 -

 

ск = 2

•(2 ««)•

 

 

г=1

 

 

г=1

Ь=1

Информация, получаемая из рассмотрения конкретной модели, слиш­ ком велика. В термодинамические функции коэффициенты акЬ ЪкЬ Е ь Щ входят не обособленно, а в некоторой комбинации. Обозначим, напри­

мер, через $к1 и

укг (к,

I = 1, 2,..., щ)

следующие комбинации:

м

 

N

,

,

 

*мгиЕ:

— РйЬ

2

ЪкГИ

(8.22)

 

 

 

Ъ

 

Соотношения (8.21) в новых обозначениях примут вид

 

 

 

т

 

 

КРы = -

Гм,

=

2

Ры.

(8.23)

 

 

 

г=1

 

Ядра свободной энергии З3 (8.13) и рассеяния <2 (8.15) с учетом обозна­ чений (8.22) записываются в виде

771 т

^>(1 — хи

1— т2) =

2

2

Р*;Л

(''''‘)+Лг(<-г»),

 

 

л-1 г=1

т

 

 

 

 

771

 

Я (« _ п

, « - т 8) =

-

2

2

(8-24)

 

 

 

К=1 *=1

 

откуда следует, что и для этой общей модели выполняется соотношение

между Я и

вида

 

 

Я ( I

Тхэ Ь Т2) =

д Г ^

^ — ^2)*§

§ 9. Вязко-упругая среда (одномерный случай)

Вязко-упругую среду будем идентифицировать моделью с бесконечным числом упругих и вязких элементов, т. е. будем считать, что для нее име­

ется бесконечное число независимых параметров состояния Ф, е[ (ъ = 1, 2,

..., оо) или, что то же'0', &Тг

д = Г - Г 0, ет . = (е ; - о с д ) | / %

<* = 1, 2, . . оо>,

(9.1)

где Е г — некоторые постоянные.

Предполагается, что образец материала однородно деформируется вдоль оси х напряжением а (I), так что упругая составляющая его деформации г (2), т. е.

ег = е — аФ,

(9.2)

однородно распределена по объему вместе с температурой Ф (2). Основны­ ми порождающими внешними параметрами являются любые два (Ф, ет) или (д, а), причем параметры равны нулю при Ь = 0.

Среда предполагается линейной в том смысле, что все параметры со­ стояния (9.1) и а, являющиеся согласно общей теории [1] функционала­ ми Ф и 8г, получаются из них с помощью линейных операторов по вре­

мени

г

г

 

■» = •&,

гтг. = 5

— т) <1гт(х),

о = $ К (I х) & т (т),

о

о

где Е г, К — положительные убывающие функции. Требования

Е г (х) > О, Е (х) > 0, йЕг (х)/Ах > О, АЕ (х)/Ах > 0

(9.4)

имеют физический смысл: напряженность 1-то упругого элемента пропор­ циональна ег*; при постоянной деформации образца и температуре, т. е.

при ет(I) =

ег Ъ, (г), из

(9.3)

имеем

гТг (0

= Кг (2) &т, о

= Е (I) Вт,

и так как ег*,

о по знаку совпадающие с е®, релаксируют, то справедливы

неравенства

(9.4).

единицы объема Т* (О, еГ1, 8т2,..., ег«,...) в рас­

Свободную энергию

сматриваемом линейном случае необходимо представить однородной квад­ ратичной формой по 8тг с точностью до аддитивной функции температуры

Т* (д), которую мы также

будем считать квадратичной

 

 

 

 

 

 

со

 

 

 

Заметим, что вследствие введения в (9.1) постоянных множителей

в

(9.4) нет необходимости вводить коэффициенты для формы ет*

На ос­

новании (9.3) выражение (9.5) перепишем в виде

 

 

 

 

 

I I

 

 

 

 

 

^

+

у

^ V ~ %1' 1 ~

^

Й8Т (^а),

(9-6)

 

 

Оо

 

 

 

 

 

причем ядро связано с Д*

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

9^ (I — ^1, ^ — ^2) ~

2

 

^2)»

 

(9‘7)

 

 

 

г—1

 

 

 

 

 

т. е. оно

положительное, симметрично

относительно

своих аргументов

 

 

 

00

 

 

 

 

 

9Ь(х, у) = 9Ь(У, х) = 2

Пг(*)

(У)

 

 

(9*8)

 

 

 

г= 1

 

 

 

 

 

и монотонно убывает по (х, у),

поскольку ядра Е г

(х)

положительны для

0 и являются монотонно убывающими по х

 

 

 

9* {хч У) !> О»

д9ь1дх <С 0, 09^/ду < 0.

 

 

(9.9)

Представимость ядра 9* (#,

у) свободной энергии

(9.6) в виде

суммы

(9.8) является фундаментальным его свойством, приводящим ко многим важным следствиям.

Полный дифференциал ф по I на основании (9.6)

 

 

1

 

= “ 17 ЬТ + дет

тх> °) йет (тх) +

 

 

О

 

Н—2“ 55~~дГ ^ ^

^ ^ ^8т ^ ^

*

о о

 

 

Из основного термодинамического соотношения, с другой стороны, он равен

бф = — ЗдТ - ЦГЫ + обе.

(9.11)

Для получения выводов из сравнения (9.10) и (9.11) теперь мы исполь­ зуем свойство модели тела; допустим, что она является мгновенно упру­ гой, т. е. мгновенная деформация обратима и для такой деформации рассеяние равно нулю; при этом выражение не содержит слага­ емых

Щ бТ + т

бв.

 

(9.12)

Тогда, учитывая вытекающее из (9.2) соотношение

6ег =

бе -

аЬТ

 

(9.13)

и сравнивая (9.11) и (9.10), получим

 

 

 

 

г

 

 

8

+ сс ^Зь {I

0) Агт(тх),

 

 

 

 

6

 

 

б

О

 

° ) * т (тг),

(9.14)

 

г

г

 

 

 

 

 

 

^

= -

Т ^

Ж -9 (* -

Х" Ь~

(Т1) ЙеТ (Т2)•

 

 

Оо

 

 

Действительно, при варьировании параметров в мгновенном процессе

мы

можем

независимо

изменять

АТ

и

с?е,

полагая

И7* — 0.

Коэф­

фициенты при ЬТ и бе в (9.10) и (9.11)

должны

совпадать, а

вариа­

ция

Ы = 0;

предполагается,

что последний

интеграл

в (9.10) не зави­

сит

от ё(2).

В

результате

получаем 3

и

а (9.14), после чего для

любо­

го процесса

из

(9.11)

находим рассеяние в

виде

(9.14). Среду

с рас­

смотренными свойствами называем

максвелловской.

 

 

 

Сравнивая выражения а (9.3) и (9.14), находим

связь между В , 3 с

учетом (9.8)

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

в (х) =

& (я, 0)

=

2

-Кг (0)

(х),

 

 

 

 

(9.15)

 

 

 

 

 

 

г = 1

 

 

 

 

 

 

 

 

причем для энтропии получаем выражение

 

 

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

 

5 = -

^

+

— Х)(1вт(х) =

^

+ аа.

 

 

(9.16)

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

Из (9.8) и (9.15) следует, что ядро релаксации К , которое мы считаем известным, представимо суммой положительных монотонно убывающих функций. Такое разложение, конечно, неоднозначно, так как суще­ ствует бесконечное множество обладающих этим свойством полных си­ стем функций ф*.(х), линейный агрегат которых с любой точностью ап­ проксимирует заданную К (,х).

По свойствам Кг (х) любую систему ф* (х) можно нормировать так, чтобы

Пусть ядро релаксации В (х) представлено в виде конечной или беско­ нечной суммы

Я(х) = 2 С»Ф (ж).

 

(9.18)

где Сг — заданные для системы ф* и данного В константы.

 

Обозначая

 

 

 

Е г (х) = ]/сТф4 (х),

Я, (0) = |ГС[, <» = 1, 2 , . . . ) ,

(9Л9)

причем этот выбор единствен, если В *пропорционально ф*, получим, что

(9.18) и (9.15) совпадают.

Из (9.8) теперь получим однозначно

^ (*, У) = 2

С«Ф* (Х) Ф* (У) = 2 Пг (Х)К1(У)-

(9.20)

Поскольку ряд

(9.18)

равномерно сходящийся (значит и

абсолютно

равномерно сходящийся), то по свойству Д* (х) и ф* {х) ряд (9.20) обладает

теми же свойствами, так как ф* (х)<1 1. Поскольку ряд для производной

-К' (*) = 2 с {Ф: (х) = - 2 с, I ч>; (х) I

(9.21)

также равномерно и абсолютно сходящийся (по свойствам Д ь ф*), то (9.20)

можно дифференцировать по х и у и затем при интегрировании по

йх,

йу менять порядок.

И7*

Обозначим в соответствии с (9.14) ядро функции рассеяния

через И7

 

*

 

= 55 (* -

* -

*>) * т (т0 * т (т2),

(9.22)

о о

 

 

 

 

причем, очевидно,

 

 

 

 

т2) = —

~ Ть 1 ~

х*)+ ИМ* — т1. * — *а).

где И7! (#, у) — любая функция, для которой

(9.23)

 

г г

 

 

 

 

$ $ И^х {I - ТЬ

I -

т2) с&т (ТО & т (т2) = 0.

(9.24)

оо

Всоответствии с (9.20) имеем при х = I — т1? у = I — т2

(ж, у) =

2

(ж) Я*(у) + И^о (*, у),

 

^0 (*, У) = 4- 2

г

(*)^ (у) - ^ (*)^ (у)1>

(9-25)

причем условие (9.17) выполняется. Таким образом, для ядра имеем

- IV (х, у) = 2 Я, (У)

= 2 С«Ф» (У)

(9-26)

Понятно, что если вместо более простого выражения И7 (9.26) мы возьмем

(при

= 0)

 

 

И7 =

Т1» * ~ Т*) =

 

=

- т 2

[■«; (< - тх) /?, (* - т2) + Д« (( - ТХ) д; (* - т,)], (9.27)

функция рассеяния ТУ* (9.22) не изменится.

Положительность функции ф + с№12Т^ очевидна по построению яд­ ра З5 (х> у). Поскольку, согласно второму началу термодинамики ТУ* долж­ на быть неотрицательна,

г

*

 

И7* = - 2 с г ($ Фг (* ~

г) * т (т)) (5 ф;(* - т) * т (т))> 0,

(9.28)

о

о

 

это условие накладывает дополнительные ограничения на функции Я 1 (х) или фг*(х)у поскольку ет (0 может быть произвольной. Требование убы­ вания Кг (х) по модулю, т. е.

Кг (х) > О, янИ > 0,

(9.29)

мы будем предполагать выполненным для любого х ]> 0.

Полные системы функций ф* (#), позволяющие представить практи­ чески встречающиеся кривые релаксации К (х) в виде (9.18), существуют, причем они таковы, что условие (9.28) будет заведомо выполнено. Дей­ ствительно, если бы мы для определения функции рассеяния воспользо­ вались скоростями изменения во времени параметров еТ, т. е. записали бы

IV* в виде

00

 

*Г = 2(ёГ)2

(9.30)

г=1

 

и представили в виде

 

*

 

— Т) * т ( Т)»

(9‘31)

О

 

то получили бы

 

г *

 

и7’ = $5 И7 (I - х г, 1 - т2) Лгт(тх) йет (т2),

(9.32)

где

 

ИЧ*, 0)

(9.33)

Значит должно существовать преобразование (9.26) к виду (9.33), где У* (я) какая-то система функций. Это требование накладывает ограничения на системы ф* (#).

Заметим, что вместо (9.5), и (9.30) можно было бы взять более общие положительные квадратичные формы

но они приводимы к (9.5), (9.30), и потому важно только существование

или несуществование определяющего выражения е* в виде (9.31), которое можно было бы взять и в более общем виде, не выявленном при рассмот­

рении механических моделей; наиболее общим линейным представлением

Г*

и # является

 

8г через внешние параметры

 

г

г

 

е? = $ V* (* — т) ^8т (т)

~ т) ^ (т)-

(9.35)

о

о

 

Мы не станем заниматься по аналогии с данным выше построением тер­

модинамики

на основе (9.5), (9.30), где 8*

взято

в виде (9.35) и гТ1 имеет

аналогичное

выражение

 

 

 

 

 

г

г

 

 

 

е Т г =

^ Я» (* — Т)

(Т) :Ь ^ '*’*! (* ”

^ ^

^

(93Г))

 

о

о

 

 

 

Отметим лишь, что приводимость формы И7* к виду (9.28) является до­ статочным, но не необходимым требованием, налагаемым на вид функций

Фг (^)•

Достаточным условием РТ* ^ 0 является практически более удобное требование положительности каждого отдельного слагаемого в (9.28), т. е. для любого <г = 1,2,...)

г г

§ § (* — т1) Яг(* — т2) с?ет (тх) с?ет (т2) > 0.

(9.37)

Оо

 

Теорема. Требование (9.28) всегда выполняется для любого практичес­ ки задаваемого ядра релаксации Я (#), поскольку оно представимо в виде (9.18), т. е. существует по крайней мере одна система функций <р* (#), для которой условие (9.28) выполняется для любого Я (х).

Для доказательства покажем, что существует система ф* = фг0 такая, что каждое слагаемое суммы (9.19) равно соответствующему слагаемому суммы (9.33)

— С{ф«(у)ф| (я) = Уг (у)У1(х). Отсюда имеем единственное решение

 

Фг0 (*) =

фгО, Ч > °. <*'= 1, 2, 3,...),

где

— любые положительные числа.

(9.38)

 

 

Теперь из (9.18) имеем

 

 

оо

 

 

Н (х) = 2 Ске х? ;

(9.39)

 

г=1

 

и мы уже знаем, что в таком виде представимы все практически задавае­ мые ядра. Из (9.6), (9.20) и (9.28) получаем

I I

т = —йт+ т5$П

~ Т1

Й8т^Т1) ^

п

Оо

 

 

 

 

 

ТГ = $$ Щ (2* - Т1 -

тг) Ы

(т2) =

о о

^

 

 

= 2

К Сг[§ ф{0 {* — Т) ЛеТ (Т)]

> °.

 

о

 

 

где

И'о (X) = 2 К С ^ Х = -

,

(9.41)

и теорема доказана.

Достаточным условием положительности IV* (9.28) для любой допу­ стимой системы функций ф* вследствие положительности коэффициентов С% является

* г

 

— $5^ ~ т1) ^ (* — Тг) <&Т (т1) Лгт(тг) > 0

(9.42)

о о

 

для всех значений г, для которых Сг Ф 0. Возможно, что это условие и необходимо.

Выбирая другую допустимую систему, в которой

 

 

#1(0) ДгО*) = С 1ф1(ж) = Л (я),

 

(9.43)

и полагая

 

 

 

 

С4= 0,

г > 2 ,

С1 = В ( 0),

 

(9.44)

получим из (9.28) для рассеяния

 

 

 

г г

 

 

 

= -

т щ - 5$# (* - т1) Я' (« - *.) ^8Т (ТО

(то

(9.45)

 

ОО

 

 

 

и из (9.6), (9.20) для свободной энергии

 

 

Т = —

4- 2д'(0)

[$ ^ ^ ~~ т) ^8т (т)] •

 

(9.46)

 

 

о

 

 

На основании (9.3) имеем

 

 

г

 

 

 

 

^ В ' (I — т) йгт(т) =

с — В (0) ет ,

 

(9.47)

о

 

 

 

 

и потому окончательные выражения

 

 

■5 = Т?- + ^

* ’ = „ ( ■ ' - * » — я ^ ) ,

 

г

 

 

 

 

б = $Д (*_т)<?ет (т);

 

(9.48)

о

 

 

 

 

причем IV* 0 для всех допустимых В (*). Такой же результат получит­ ся, если вместо (9.43), (9.44) принять все В г одинаковыми, с точностью до множителя совпадающими с В [25]

В г (х) = А (В (х),

(9.49)

т. е. предположить, что все внутренние параметры состояния с точностью до множителей релаксируют одинаково.

§10. Сплошная вязко-упругая среда максвелловского типа

Линейная теория вязко-упругости сплошной начально изотропной среды приводит к соотношениям Больцмана между напряжениями — деформа­ циями — температурой (если, конечно, использовать приведенное время)

 

 

г

 

 

 

 

% =

{I х)йе^ (т),

0 = <\)Я1Ц — т)с?0т (т),

(10.1)

где

 

о

 

о

 

— девиаторы, д =

Г — Г0, 0 = (Нуи,

 

 

0 г = 0 — Зад.

 

 

(10.2)

О сновными порождающими

внешними параметрами, следовательно, яв­

ляются девиатор ец, скаляр 0 т и параметр состояния— температура Т =

=

Г0 +

Ф» причем в (10.1) д явно не входит, т. е. предполагается, что если

в

процессе,

начинающемся из естественного состояния {I = 0, е^ = 0,

0

= 0, д

=

0), изменяется только температура {I 0, д =/= 0), напряже­

ний не возникает. Эти же параметры являются порождающими и для функ­

ций состояния и рассеяния, причем в случае I 0,

= 0 = 0, д =^= 0

изменяются только функции состояния.

 

Можно было бы повторить рассуждения предыдущего параграфа о

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

т.е. двухкратными свертками вида

*I

Т

2То = Т2 \Л'\

(* —

1— т«)

(т1) йеа (тг) +

 

г г

5 о

 

 

 

 

 

 

 

 

 

 

 

 

+. у ^

(« -

« -

Т,) Й0Г (Т,) Й0Т (т2),

 

о о

 

 

 

 

 

(10.3)

 

I (

 

 

 

 

 

=

— Ту, I — т2) йец (тх)

(т2) +

 

 

 

— тх,1 — т2) <?0Т (тх) (Ют(т2).

Работа внешних сил (для единичного объема)

 

с^6е$;*— а60 +

 

 

 

 

(10.4)

Дифференциал Ч?* находится следующим образом:

6Т +

ЬТ = йе« (*)

{I -

т, 0) йек) (т) +

 

 

I

о

 

 

г г

 

 

 

 

 

+

60т (*)

{I.т) Й0т (т) +

Ы . А^ ^ ± . 3>е {I - Ту, г - тг2) X

 

 

0

I I

 

 

0

0

 

 

 

 

 

 

 

Из основного термодинамического соотношения

6Т + &Й = сзб0 + 8ц6е15 — \У*Ы,

(10.6)

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

*

0

I

с = ^ 5й» (* — т>°) С*),

О

I

‘* = -5Т + 3 а $ 5&в(*“ !С,О)Й0т(Т)’ (107> 0

* I

-

= Т о$ \о4

"

{1 ~ Т1’

1 ~ %г)йе« (Тх) (Тг) +

 

 

* г

 

(*— т!, * — т2) <ют (тх) а ет(т,).

 

 

-I- у $ 5 “5^ Л

 

 

о о

 

 

 

 

Сравнивая выражения

и а в (10.7) и (10.1), заключаем

 

^ в(* _ т ,0 ) = Д ( г - т ) ,

$>в(1 -т ,0 ) = В 1Ц -т ),

(10.8)

и потому для энтропии имеем выражение, совпадающее с (9.16),

 

8 =

I

 

 

 

 

- ^ - + З а ^ Т ? (г - т ) й 0 т (т) = - ^ - + Злс5.

(10.9)

 

о

 

 

 

 

Поскольку правая часть равенства (10.3), определяющего свободную энергию г|э, образована как сумма квадратов параметров состояния типа

$г, причем вследствие ее инвариантности (невозможности произведений типа йвц (т2) й§т (т2)) эта сумма состоит из двух сумм квадратов, функции Зье (х, у) и Фо (#, у) должны быть представимы суммами .

3 Ье ( х , У ) =

2

г к ( х ) г к(У),

 

 

К=1

(10.10)

 

 

оо

 

^«(ж, у )=

2

ги (х)г1к(У)>

 

 

к—1

 

и потому, согласно (10.8),

 

ОО

 

Д (ж) =

2

гкФ )гА х)’

 

 

1

(10.41)

 

 

 

/?х(д:)—

2

Ал(0) г1к(х)>

 

к=1