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

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

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

Следовательно, повышение температуры на 70° имеет показатель (5.11) /г4х= 4, т. е. «ускоряет течение времени» в 104 раз.

Использование принципа эквивалентности для проведения экспрессиспытаний и для других целей очень заманчиво и является одной из не­ многих возможностей прогнозирования поведения материалов в тече­ ние длительного времени. Однако условие его справедливости и точности, т. е. возможность преобразования серии кривых (5.20) к виду (5.21), должно тщательно проверяться на прямых опытах с достаточно большим диапазоном изменения истинного времени. Реалистической оценкой ошиб­ ки в значениях показателя преобразования времени тг12 (5.11) будет 10— 20%, т. е. при сокращении масштаба времени на /г12 порядков за счет по­ вышения температуры ошибка в истинном времени может быть на (0,1

— 0,2) щ 2порядка.

III.Некоторые термодинамические функции

изаконы сохранения

§ 6. Простейшие вязко-упругие модели

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

Рассмотрим термодинамику простейших моделей (Максвелла и Фохта) [23]. Для всех моделей изотермические физические соотношения одинаковы

а' -

Е г \

 

а" -

це'.

 

(6.1)

Структурные

кинематические

и

динамические соотношения различны.

Для модели

Максвелла

они

имеют

вид

е' +

е" =

е,

а' -

а"

-

а;

(6.2)

для модели Фохта

 

 

 

 

 

 

8/ =

8,; =

8,

О*

о" — 0.

(6*3)

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

Обозначим разность между температурой Т и некоторой фиксирован­ ной Го в начальный момент

О - Т - Т 0.

(6.4)

С учетом теплового расширения физические соотношения (6.1) примут вид

а' = Е (е' — ад),

а" = цё",

(6.5)

где а — коэффициент теплового (линейного) расширения.

Таким образом, термодинамические функции моделей Максвелла и Фохта: энергию, энтропию, свободную энергию можно рассматривать как функции температуры Ф и упругой деформации е \ Свободную энергию любой из моделей Т (Г, е') представим в виде разложения в ряд Тейлора, сохранив только квадратичные члены с Ф и е'

Т’ (Г, е') =

ф0 +

+

'Фге' + 'Фп&2 + 'фхг'&е' +

ф22е'2.

(6.6)

Свободная энергия

Т связана

с внутренней энергией

и и энтропией

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

 

 

 

 

 

Ч = и — Т8.

 

 

 

(6.7)

Первый закон термодинамики связывает внутреннюю энергию притока

тепла с работой внешних сил

 

 

 

 

ди — 8(2 +

оде,

 

 

 

(6.8)

где первый член правой части представляет собой приток тепла, а второй — работу.

Второй закон термодинамики утверждает существование энтропии

определяемой

равенством

 

 

Тд8 -

IУМ = б(?, ТУ*д1 >

0,

(6.9)

где И7* — функция рассеивания энергии

в единице объема (в единицу

времени). Для

обратимых процессов

ТУ* = 0, для необратимых — пред­

ставляет собой часть работы внутренних

сил, необратимо переходящих

в тепло.

 

 

 

Рассматриваемые нами модели вязко-упругого тела могут быть только мгновенно обратимыми, так как при очень больших скоростях деформации вязкие элементы практически не деформируются. Вообще же процессы в них необратимы, и работа вязких элементов полностью и необратимо переходит в тепло. Поэтому для рассеяния ТУ* имеем

В7* = <у"ё" =

цё"2 > О,

 

= о"дг"

= 1хг" дг" > 0.

(6.10)

Подставляя выражение притока тепла (6.9) в (6.8) и используя (6.7), получим основное термодинамическое соотношение, одинаковое для обеих

моделей

 

 

 

 

 

 

сПГ +

Д дТ = оде — ТУ* д1 =

о'дг',

(6.11)

поскольку для первой модели а" =

а'

= а, дг дг" =

йе', а для второй

дг" = дгг — дг,

о — о" = а'. Отсюда находим, учитывая (6.6),

— 5

=

дУ/дТ =

Ч>! + ^>12е' +

2^ПФ,

^ 4

о' =

дЧ1дг’ = г|)2

+

 

е \

 

Сравнивая полученное выражение а' с (6.5), находим постоянные

ф2 =

0,

2Ф22 ^

» ^12 =

сх>Е]

(6.13)

обозначим постоянную

 

 

 

 

Константа ^ аддитивно входит в выражение энтропии (6.12) и потому мо­ жет быть отброшена. В результате получаем следующие выражения 4% Я и и = У + № + Т08 + с:

2Т = - [— ■+ ^ Е ) № + Ев1

=

■/ О

+ Ев'%+ 2Н> + 2аТ0Еег,

 

■5 - - - ж

= ( ж + - » * Е) # + аЕг;-'

причем через е'т обозначена упругая составляющая полной деформацшГе' упругого элемента

ет = е — аб'.

(6.16)

Смысл постоянной с в (6.15) найдем из закона сохранения энергии Аи =

= 6(? + оАе, который, согласно (6.15),

(6.2) и (6.3), запишем для обеих

моделей в виде

 

6() = II-ат ЕаТйв' — оЧв".

(6.17)

Применяя его к ненапряженной модели в начальном состоянии при тем­ пературе Т — Г0 (Ф = 0), т. е. при а" = 0 и а' = Е (е' — аФ) = 0, и давая этой модели свободно расширяться за счет сообщаемого тепла 6()»

получим приращение деформации

е' из Ао'

= 0, т. е. Ае' = айФ. Следо­

вательно, из

(6.17)

 

 

 

 

6<? =

+ Еа2Т0) АТ.

 

 

(6.18)

Отсюда ясно, что с + Еа2Т 0 есть начальная

объемная теплоемкость при

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

с = ср -

Еа2Т0.

 

 

(6.19)

Из (6.15) получаем простые

выражения ф, /У и и через параметры со­

стояния моделей

 

ет

 

^

 

2 Т = -

^ Ъ

2 + Ее'1

3

=

Еет,

 

 

 

 

 

 

(6.20)

2 и =

2срд +

2аЕ (Т0+

■&)гт+

Евт].'

Уравнение (6.9) теперь дает закон изменения температуры модели при

заданном во времени притоке тепла

6(), причем ]№*А1

определено (6.10)

То <и

йг

Р \ ль) — осТЕ

д&'

(6.21)

дЛ

Поскольку в любой модели произвольно мы можем задавать лишь внешние (порождающие) параметры, т. е. либо а (2) и б(}/А1, либо е (2) и Ь(11А1, либо их комбинации, то уравнение (6.21) должно быть дополнено структурными соотношениями (6.2) или (6.3) и физическими (6.5), после чего уравнение (6.21) может быть разрешено.

Свободная энергия, энтропия, внутренняя энергия являются всегда функциями параметров состояния (в наших моделях любой пары д, е'

или *&, ет, или Ф, а')- Но сами параметры состояния всегда являются не~

которыми операторами над порождающими параметрами (д, а, е), причем операторы эти и отражают физические свойства тела [1].

Оператор дифференцирования по приведенному согласно температур­ но-временной аналогии времени I' будем по-прежнему обозначать буквой Д

Д = (Ц йЬ', й11 = Ж/ах, р = [х0а?- (6.22)

Для модели Максвелла из структурных и физических соотношений

найдем

 

 

 

 

 

Е

 

 

 

 

 

 

 

 

Ет— 1 + ©а-11

ет - е

ад,

со = — ;

 

ЦЮ»

(6.23)

 

 

 

 

 

соД_1ег

 

е' = ет +

 

ад,

 

 

I

в

 

 

е" = Ех — Ет — 1 + соД'1

 

Здесь а-1 — обратный для а оператор, т. е. если

 

 

с; =

йу =

с?у/Л' ,

 

 

 

(6.24)

то

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У == $*7Л' =

 

 

 

 

(6.25)

 

о

 

 

 

 

 

 

 

причем для

операторов Д, Д"1, постоянных

и функций времени / ь ...

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

 

«1 (С1/1 +

с2/2) =

М Д + МД.

 

 

 

а -1

(С1д

+

С2Д)

= М " 1 Д

+ с2л-%,

а а -1

= а -*а = 1,...

 

Из (6.23) находим параметр состояния Ет(т .

е. и е')

 

 

е'т=

г

 

Лет (г')

 

 

 

 

^

 

 

 

 

(6.26)

 

О

 

 

 

 

 

 

 

и параметр

модели, не являющийся параметром состояния,

 

 

г

 

ЕТ (т') йт’.

 

 

 

(6.27)

е" =

(0 5

 

 

 

 

 

о;

 

 

 

 

 

 

 

Из (6.10) и (6.20) теперь находим выражения рассеяния, свободной энергии и энтропии через общую деформацию е {I') и температуру д (^')

 

 

 

е “(2< ^ ^

Лет (х[)^ет (г'),

 

 

о

о

 

 

 

г

 

(6.28)

8

Ф + <%Е^

Лет (т'),

^

 

Т оо

<*2> * т

К ).

 

 

 

 

Связь между напряжением а и деформацией е для этой модели имеет

вид

V

 

(6.29)

о = Еет = е \

д.Ет(?').

о

 

 

Для модели Фохта имеем непосредственно из структурных соотноше­ ний (6.3)

г

е" =

е' =

е = ^

— т') йг (т'),

(6.30)

 

 

 

о

V

 

 

 

 

 

 

 

гт=

е' — а# = ет = ^ к {V — т') йгт (*')»

 

 

 

 

о

 

 

где к (2) — функция Хевисайда. Из (6.20) и (6.10)

ч' = - Й - « * . + - г 'Ь

« =

+

 

 

 

 

 

 

(6.31)

и из а = сг' +

а" — Егт +

цё" находим

 

<5 =

р

I

^б?6

 

 

(6.32)

Е г т +

|х0 - ^ г ■

 

 

Отметим, что для модели Максвелла ядра релаксации напряжения Д,

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

(6.28)

 

$>(II' -

I, V -

Т)) =

 

 

 

(6.33)

<? (*' — 5, V — т,) =

 

 

 

 

связаны между собой соотношениями

 

 

 

 

25й ( Г - |, 0 ) ,

 

 

 

<? (*' -

*' -

г|) = -

А 53 (Г -

Е, * ' -

Л),

(6 ‘34)

причем энтропия имеет выражение через Д

 

 

 

 

 

V

 

 

 

8 = -

- I? 1 =

^ 5 -0

$ Д ( * ' -

(х ').

(6 .35)

 

 

 

о

 

 

 

§ 7. Обобщенные максвелловские модели

Рассмотрим теперь более сложную модель, составленную из п максвел­ ловских элементов, соединенных параллельно, как показано на рис. 13. Для этой модели кинематические структурные соотношения состоят из п уравнений

Бг + е* = е, = 1, 2,..., п).

(7.1)

Структурных динамических соотношений будет п -+■ 1, и они имеют вид

п

Ог = б$,

2

^ б.

 

«=1

 

Физические соотношения остаются прежними, причем мы исполь зуем ра­ нее принятые обозначения, опуская штрихи для приведенного времени и индекс «О» у р,*.

I = Ер'ти

о'- = Ме",

г’п = г \ — аФ,

(7.3)

е т = е — ай, <г = 1, 2 , ге>.

Связь между напряжениями и деформациями записывается в виде опера­ торного соотношения

гь

— 2

Е*&т

(О,- =

Е .

(7.4)

 

1 + со.а-1 ’

 

 

Выражение свободной энергии ф в соответствии с (6.20)

 

 

 

п

 

 

Т (е', Г) = Т'(О) +

~ 2

 

(7.5)

 

 

г=1|

 

 

Точно так же для функции рассеяния имеем обобщение формулы (6.10)

 

= 2

^ 8 ? .

 

 

 

(7.6)

Энтропия

1—1

 

 

 

 

 

 

 

 

 

 

 

с _

дЧГ

ЛЧГЩ

, V I ? '

,

(7.7)

 

вт

ял

а 2

 

 

 

 

г=1

 

 

а теплоемкостью

с'р будем называть

 

 

-

Го

 

= ср, -

Т ' (О) = 5

^

Ср (<>') Ю',

с — ср — а2 2

Ек.

 

 

(7.8)

 

 

г=1|

 

 

 

 

В формулы для свободной энергии (7.5), энтропии (7.7) и функции рас­ сеивания (7.6) подставим выражения упругих и вязких деформаций через полные

 

V

со ^ х&т

 

Втг = 1 + ю.а-1’

=

1 + (0.Й"1 ‘

(7.9)

После этого выражения свободной энергии, энтропии и функции рассеи­ вания принимают вид

Т

2Го

2 2

(1 + « ^ -1)

 

^

= ~т7 * + 01 2

! + ш.д-1 *

(7.10)

 

1=1

 

*

 

И ' - = | / Л ( т ^ Р )*.

Для простоты в (7.8) мы предположили, что ср не зависит от

Заметим,

что в (7.4) и (7.10) участвуют однотипные выражения

 

 

= 1 + ш.а-э- *

 

эквивалентные дифференциальным уравнениям относительно щ

 

щ + в>%щ = ет , <г = 1, 2,

 

Решения таких уравнений при нулевых начальных условиях имеют вид

г

 

 

Щ= г + 4 а-1ег - $

А т СО

(7.11)

*О

ивыясняют смысл стоящего перед етдробного оператора. Теперь для сво­ бодной энергии, энтропии и рассеяния рассматриваемой модели имеем [24]

 

+

1=1

А т ы А т (т2),

 

 

 

ОО

 

•У =

« + « 2

Щ$

Ат (т),

(7.12)

 

г= 1

о

 

 

п

^ ^

 

 

 

Т Г = 2

 

 

А т (тх) А Г (т,).

 

г=1

ОО

 

 

 

Уравнение связи между напряжением, деформацией и температурой

(7.4)

запишется в виде

 

 

 

71

I

 

 

 

о = 2

Ек$

А т (О-

(7.13)

1—1

О

 

 

 

Справедливость полученных соотношений можно проверить непосредственно

п одстановкой

выражений (7.10) и (7.4) в

основное

соотношение

(6.11)

аЧ +

5 ЛТ +

= ай8,

 

(7.14)

которое при этом обращается в тождество.

(6.34) факт, что ядро

инте­

Обратим внимание на отмеченный ранее

грального выражения

функции рассеяния

в (7.12)

получается из

ядра

интегрального слагаемого свободной энергии дифференцированием пос­ леднего по времени и умножением на —1.

Если вместо (7.12) пользоваться выражениями функций состояния че­

рез параметры

состояния Ф и порождающими параметрами

вида (7.5),

(7.6), (7.7), то

нужно знать функциональные связи е», с

основными

порождающими параметрами а, 8, входящими в уравнения механики сплош­ ной среды. Согласно (7.9), (7.11) получаем

е; — а# = ъ’тг = 1 +8^ = 5 А т (т),

го

&*= = со^гь — 8 — ад, <г = 1, 2, ..., тг>. д (7.15)

Единственное здесь ядро, как видно из (7.12), служит базисным для всех других ядер.

§ 8. Сложная линейная модель

Рассмотрим теперь общую модель, составленную из М упругих и N вяз­ ких элементов, которые связаны так, что нет свободных элементов и нет непосредственной связи между любыми дву“мя однотипными из них, так как такие пары образуют единичные элементы. Свободная энергия модели равна упругой энергии гуковских элементов с точностью до аддитивной функции температуры, т. е. имеет прежнее выражение (7.5). Рассеяние также равно мощности всех вязких элементов, т. е. имеет выражение (7.6).

Пусть модель определена 5 структурными динамическими

М N

 

,^3

 

<3;— В к$=. О,

= 1, 2, ..., 8)

(8.1)

 

г—1

; = 1

 

 

 

и Ь

структурными

кинематическими

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

(Ь — N + М —

- 5

+ 1)

 

 

 

 

 

М

N

 

 

 

 

^

-^гй8Гг+ 2

= 1, 2, ..., Ь),

(8.2)

 

г= 1

; = 1

 

 

 

причем по-прежнему

&т = е — аб' и

=

1. Последние

соотношения могут быть на основании физических законов (7.3) сведены

кследующим:

МN

2

+

2

= ^ 8г-

(А =

1 ,2 ,...Д ) .

(8.3)

1 = 1

*

3 = 1

^ 3

 

 

 

Система М + N + 1 уравнений (8.1),

(8.3)

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

N + 2

переменных (включая а и

8у) может быть сведена к одному операторному

■соотношению

 

 

 

 

 

 

<5=

8т,

 

 

 

 

(8.4)

где и Р — дифференциальные операторы. Их вид разбирался в предыду­ щей главе. Из системы (8.1), (8.3) или из рис. 15 видно, что эти соотноше­

ния могут быть разрешены относительно о\ и

Ог =

 

(г =

1 ,2 ,...,М ),

(8.5)

 

(у =

1 ,2 ,..., Я),

 

 

 

причем порядок операторов

и Р

совпадает, как видно из матрицы на

рис. 15, а порядок операторов (^ и Р

либо совпадает, если число кинемати­

ческих соотношений Ь меньше или равно числу упругих элементов М, либо

порядок оператора

на единицу больше порядка оператора Р, если чис­

ло кинематических соотношений Ь больше М .

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

то корни многочлена Р

(Я) действительны, отрицательны и все различны,

причем, если отсутствует упругий элемент, соединенный параллельно с •остальной системой, то среди этих корней нет нулевого. Если порядок опе­ ратора Р (Д) равен яг, то имеем пг корней

Разложим теперь рациональные функции М{ (X) и Ку (X)

М* (а) =

<*; (Х)/Р (X), ^ (X) =

 

= 1, 2

, М),

(а)/Р (X),

 

(8.7)

 

 

 

 

( /=

1, 2........ЛГ),

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

 

 

 

 

| м , ( М

-

” 2ъ |

|

т г

 

 

X-- А,!

 

 

 

(8.8>

 

 

|

|

 

 

х ^ ( ч

=

 

 

 

X—Яг

XХт

 

 

В разложениях

(8.8) отсутствуют члены

- ^ и - ^ - ,

ибо

многочлены

(^(Я) и (^' (Я) обязательно имеют хотя бы по одному нулевому корню. По­

этому дифференциальные

операторы М* (Д) и Ку (Д) можно представить

следующим образом:

 

 

 

 

 

! ^ м <

( < ‘> - ^ +

г а = г + -

 

ат1

1

—Я,т Д-1 9

^г <а> =

N,•(^ =

60, +

^

__

+ 1

^тз

ХхД"1

~ Х^й-1 *

В выражениях (8.5) для с'*, о$ мы, следовательно, должны суммировать слагаемые типа (7.11)

1 -1 :а -1 8 т = ^ Хг(*~т)<*ег (<)•

 

 

(8.10)

1

0

 

 

 

 

 

Следовательно, получаем

I

 

^

 

/

771

 

 

г1т= 2Г =

т~ 2

«м

дят (х) =

$П*(г — т) <?ет (г),

 

»

пг

 

г

 

г

(8.11)

 

 

 

е) =

2

Ь*у ^ ех*(*~т) $5Т (т)

= ^ Ру (* — х) йгт(т),

 

*

]

 

о

 

О

 

где обозначения ядер Пу (0,

Ру (0

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

(7.5)

запишем в следующем виде:

 

 

 

 

 

 

 

 

 

Т (7\ е) =

ф' (О) + у

(I -

хг, I -

х2) Лгт(тх) * т (т2),

(8.12)

 

 

 

ОО

 

 

 

где ядро Р {I — тх, I — т2) представляет собой сумму экспонент

М

771 т

Р <« ■- т,, < - т.) =

2

4 -

2 2

=

 

г=1

г

Ь*=1 *=1

 

 

М

 

 

 

=

2 а д ( * - » х ) п 4( * - т я).

 

г=1