книги / Основы математической теории термовязкоупругости
..pdfСледовательно, повышение температуры на 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Н> + 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’ |
8г = |
1 + (0.Й"1 ‘ |
(7.9) |
После этого выражения свободной энергии, энтропии и функции рассеи вания принимают вид
Т |
2Го |
2 2 |
(1 + « ^ -1) |
|
^ |
= ~т7 * + 01 2 |
! + ш.д-1 * |
(7.10) |
|
|
1=1 |
|
* |
|
И ' - = | / Л ( т ^ Р )*.
Для простоты в (7.8) мы предположили, что ср не зависит от |
Заметим, |
что в (7.4) и (7.10) участвуют однотипные выражения |
|
8у |
|
= 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 |
~ Х^й-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