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

Курс механики сплошных сред

..pdf
Скачиваний:
19
Добавлен:
15.11.2022
Размер:
20.58 Mб
Скачать

Читатель

легко найдет ограничение на функцию ф (Т, В\, BJJ, 5 щ ) в

урав­

нении (26) и,

в частности,

в уравнении (27), вытекающее из предположения

о су­

ществовании

естественного

состояния.

 

VIII.2.4. Изотермические и изоэнтропийные движения упругой среды. Энергия деформаций. Как было сказано выше, процесс, происходящий в упругой среде, называется изотермическим, если температура остается постоянной, т. е. Т = ТВ. В этом случае сво­ бодная энергия ф, описываемая уравнением (16), зависит только от компонентов тензора деформаций Грина —Лагранжа; то же можнр сказать и об уравнениях, отражающих законы поведения, например (20). В этом случае обычно говорят, что функция щ(Ь), определяе­ мая уравнением

w (L) =

р0ф о, L),

(28)

где р0—плотность в исходной

конфигурации,

представляет собой

энергию деформации частицы, отнесенную к единице объема исход­ ной конфигурации.

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

очень важном случае, когда механические

и термические эффекты

не связаны и могут изучаться независимо.

Учитывая закон тепло­

проводности Фурье (21), можно заметить, что этот случай рассмат­ ривается как предельный, когда коэффициент теплопроводности k стремится к бесконечности.

Другой случай отсутствия взаимосвязи между этими явлениями — адиабатные процессы —также предельный, когда коэффициент теп­ лопроводности равен нулю. Здесь предпочтительнее основываться не на свободной энергии ф(Г, L), а на внутренней энергии e(s, L), которая также позволяет дать полное описание термодинамических

свойств среды. В общем виде

имеем

 

 

 

 

 

бе _rp ds

 

де dLccP

 

(29)

 

Ж — 1 Ж + С0,

“ “ «.«Ц

dt

Рассуждая, как и в VIII.2.1, приходим

к тем же результатам.

Закон поведения

может быть

записан, например,

в виде уравне­

ния (20):

 

 

 

 

 

 

 

 

2 = рр з г рт-

 

 

 

(30)

Заметим, что в уравнении (20) производная

^

берется при постоян-

ной температуре,

в то время как

в (30)

дв

берется

при s = constI.

Этот вывод, полностью аналогичный данному выше выводу для общего случая, дает возможность выявить упрощения, вытекающие из пред­

положения об адиабатности процесса.

В самом деле,

рсо = a/;.D;y,

так как внутренняя диссипация равна

нулю и ^ = 0.

Принимаем,

как обычно, г = О, так что уравнение (VII. 15), эквивалентное энер­

гетическому уравнению, в нашем случае -^- = 0.

Итак, энтропия частицы, прослеживаемой в ее движении, оста­ ется постоянной. Если допустить, что энтропия исходной конфигу­ рации остается постоянной, то отсюда будет следовать, что энтро­ пия всей системы также остается постоянной при движении, назы­ ваемом изоэнтропийным, и равной, например, s0. В этом случае при­ ходится также вводить энергию деформации (на единицу объема исходной конфигурации)

uyfl(L) = p0e(s„, L),

(31)

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

различными.

Итак, в адиабатном случае уравнения движения фактически не связаны с уравнением энергии. Обратим внимание на то, что здесь имеется аналогия с движениями идеальных жидкостей. Впрочем, идеальную жидкость можно рассматривать как особый случай гипер­ упругой среды (если отказаться от предположения о существовании естественного состояния). Как и в случае идеальной среды, можно показать, что значения удельной энтропии в «тылу» фронта удар­ ной волны больше ее значений перед фронтом ударной волны.

VII 1.3. ЛИНЕЙНАЯ ТЕОРИЯ УПРУГОСТИ

Как уже отмечалось в VI.3.3, можно построить линеаризован­ ную теорию упругости, справедливую лишь в асимптотическом смысле для малых возмущений, но которая в силу своей простоты будет полезной для многих приложений. В этом случае в соответствии с общими указаниями, данными в V.4.5, можно заменить LaP на ety

компоненты тензора деформаций при малых возмущениях, а ^ j r -—

на Dij.

V III.3.1. Законы связи напряжений с деформациями. Будем счи­ тать, что справедливо предположение о существовании естественного состояния. Если ограничиться изотермическими или изоэнтропийными процессами, когда существует удельная энергия деформации,

обозначаемая

через t^(eiy), то с

учетом предположения о малости

возмущений

величины

 

 

дии

р>

 

Р«(з)==Щ ° и

и Р“ = aUDU

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

внутренних сил. Итак, имеем соотношение

dw

СТ'У ~ дъ.

которое можно установить также путем непосредственной линеаризации одного из точных законов поведения, например (19).

Так как речь идет о линейной

теории, то

правая часть

равен­

ства (32) (которая, как известно,

равна нулю

при е|7 = 0)

должна

быть линейной относительно

Следовательно, в качестве

функ­

ции w{BiJ) достаточно выбрать форму, квадратичную по переменным е/у., которая в соответствии с известными результатами термоста­ тики должна быть положительно определенной.

Итак, пусть

 

W==~2 aijhkeifehk

(33)

Постоянные aijhk, называемые упругими постоянными среды, — суть компоненты тензора четвертого ранга—тензора модулей упругости среды. Они удовлетворяют соотношениям симметрии

aijhk = ajihk = aijkh 5=5ahkip

(34)

Из этих соотношений следует, что правая часть (33) представ­ ляет собой общую квадратичную форму от шести независимых пе­ ременных е^у, имеющую, как это легко проверить, 21 различный

коэффициент (квадратичная форма m-го порядка имеет т(т^~ ^ неза­

висимых коэффициентов). Таким образом, согласно (32) закон пове­ дения запишется в виде

° ij—‘aijhkehh-

(35)

Отметим также тождество

(36)

° '=ST (W

Заключение. В общем случае в линейной теории упругости удель­ ная энергия деформации представляет собой квадратичную положи­ тельно определенную форму компонентов B{j (33), поведение среды при этом определяется равенством (32) или (35). Эти законы, опре­ деляемые 21 упругим модулем с размерностью напряжений, образуют тензор четвертого ранга, называемый тензором модулей упругости среды.

V III.3.2. Двойственная формулировка. Уравнение (35) можно разрешить относительно 8,-у и выразить последние как линейные функции ъ(] (так как квадратичная форма (33) положительно опре­ делена), определив тем самым некоторую функцию ср(сг/у), для ко­ торой

Ф ( о и ) =

W (в/у) =

a i j B‘J ~ w (е<у) =

т

(37)

 

Отсюда получаем,

что

сг,7 de,7 —dw = е,7 dcrl7,

 

dtp =

d<x,7 e(7 +

 

и тогда

 

Р

= i!L

 

lJ даи

Функция ф (аху) является

преобразованием Лежандра функции

wfaj). Энергия деформации здесь непосредственно выражена через

напряжения.

(38)—это закон поведения в форме

зависимости де­

Уравнение

формаций от напряжений.

формой перемен­

Очевидно,

что Ф ^у) является квадратичной

ных а{/:

 

 

 

ф(ст«7)=тг A ijhk°ijV hk,

( Щ

где коэффициенты AiJhk удовлетворяют соотношениям симметрии, аналогичным соотношениям (34). Следовательно, в этом случае урав­

нение (38) запишется в такой

форме:

 

8i/

== Aijhk°hk)

(40)

причем коэффициенты Aijhk получаются при решении уравнений (35) относительно ex-y. Таким образом установлена двойственная форму­ лировка изложенной выше теории, когда еху и оуу поменялись ролями.

V III.3.3. Изотропные среды. Если теперь сделать дополнительное предположение о том, что среда изотропна, то получим квад­ ратичную форму ш (в*у), выраженную лишь через инварианты мат­ рицы е и являющуюся фактически [как и в уравнениях (8)] линейной функцией от е{ и еи, которую можно записать так:

Ю= 1/2{^(8,7еуу) +

2це178,7}.

(41)

Вводя шаровые и девиаторные составляющие по формулам

aij = ст' 6,7 + °?i<

<*& = 0,

 

е1/ = 8^ |/ + е(/.

Екк= 0,

 

получаем выражение

 

 

И» = 4 {3 (ЗЯ+ 2(Х) (е')* + 2(ief)j =

i {9* (е*)2 + 2це^е^},

(42)

аналогичное выражению (11).

 

квадра­

Необходимым и достаточным условием того, чтобы эта

тичная форма была положительно определенной, является требование ЗА,-|- 2ц = ЗК О, ц 0.

Таким образом, вновь приходим к результатам, полученным ранее в VI.3.3 и, в частности, к закону поведения (VI,27). Приведенные выше неравенства следуют из более жестких ограничений, выте­ кающих из данных эксперимента. Энергия деформации, выраженная, через напряжения а/у, после введения модуля Юнга и коэффициента Пуассона запишется в форме

Ф = “2 ° и ъи = ~2 { ~Е a iJ°U

Е

*

( ^ )

а закон поведения (38) вновь приобретает форму (VI,31).

e/у daи = 3eJ da5

Вводя в (43) шаровые и девиаторные составляющие, имеем также

Ф =

I — 2v

W -(■ 1+V ofiofij.

(44)

Е

Замечание. Выписывая тождества

а/у deij= 3os de-j-afy defy = da;t dofj = dqp

и используя выражения (42) и (44), вновь находим законы поведения (VI,28) и (VI,32), связывающие между собой шаровые части и девиаторы:

 

os =

3/(ef ,

ai/ = 2[xefy,

es

1 — 2v

os,

D _1 *4~v n

Е

 

eij---- £—<4b

которые с учетом тождеств (VI,33), связывающих /С, р и Е, v, эквивалентны.

VIII.4. ТЕРМОУПРУГОСТЬ

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

Возьмем за отправную точку следующий экспериментальный вы­ вод. Если нагреть упругое тело, свободное от напряжений, то оно удлинится. В рамках линеаризованной теории и в предположении, что среда изотропна, можно допустить, что деформация будет про­ порциональна 0 = 7 — 7 0, где 7 0—температура естественного состоя­ ния:

8 = а 0 1, е|у. = а 0б-у,

где а —не зависящий от температуры коэффициент, который, как в теории малых возмущений состояния с однородной температурой 7 0, будем предполагать постоянным. Величину а называют коэффи­ циентом линейного расширения, а За —коэффициентом объемного расширения.

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

Этот закон будет иметь

вид

 

 

 

е17 = ф

аи -

£ (ам) 6,7 + сс0б,7

(45)

или,

если выразить напряжения

через деформации и изменение тем­

пературы,

 

 

 

т. e.

Оц = К (ekk—3a9) S,7 + 2р (el7—a06iy),

 

оij —

 

2{.i8fy—3/(a06fy,

(46)

 

 

где

3K = 3X+ 2p — модуль

объемного расширения —сжатия.

Этот

* Для простоты ограничимся изотропными средами.

закон содержит слагаемое, пропорциональное разности температур 0, которая и вызывает тепловые напряжения.

Формулы (45) и (46) и закон теплопроводности (закон Фурье) определяют законы поведения в теории термоупругости. Представ­ ляют интерес вытекающие из этих законов выражения для термо­ динамических функций. Для функции свободной энергии имеем

р0dip= — p0s d r + ои de,7 = — p0s dT —3/Ca9 de, +

+ -id[>.e,7e// + 2pel7 eiy],.

или

d [р0г|з — су(e,7)] = — p0s dT —3/Ca0 de„

(47)

где через до(е,7) обозначена плотность энергии деформации при по­ стоянной температуре Г 0, а через е,—первый инвариант е. Величина р0ф—w является функцией от Г и инвариантов е„ еп, еш , и по­

этому из

(47) следует, что энтропия s = —

не зависит от

инва­

риантов

еп и е1П и

д (PQS) = д (3К аб) = д „

 

 

 

 

 

 

 

так

что

 

дв\

дТ

1

 

 

 

p0s = 3/Cae,+/i' (Г),

 

(48)

 

 

 

 

где

h (Т) — функция

одной лишь температуры,

а /i'(T’) — ее

произ­

водная.

Так как s

определена

только с точностью до аддитивной

постоянной, то без ограничения

общности

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

h(T0) = h' (Т0) = о.

 

дифференциального уравнения (47)

 

Тогда после интегрирования

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

 

 

 

Роф = ш — 3/(a —Т0) е,—h(T).

(49)

Отсюда легко выводится выражение для плотности внутренней энергии

Ро*= Ро (Ф+ *Т) = + 3/Са7>, + ТК (Т) — h (Т).

(50)

Эта формула позволяет дать физическую интерпретацию функции h(T). Величину дв (Г, е<7) называют удельной теплоемкостью при

постоянных деформациях; она равна количеству теплоты (в обрати­ мом процессе); полученной единицей массы при повышении темпе­ ратуры на один градус при постоянной деформации. Обозначим эту величину через Се. В среде, для которой внутренняя энергия за­ дается уравнением (5), Се зависит только от температуры, и поэтому, зная функцию Се(Г), можно принять

р0Ce = Th"(T),

h'(T) = p ^ T' ^ d T \ T h '(T )-h (T ) =

= р0 Sr. с е (Г') d r 4

(51)

и получить в явном виде уравнения (45), (46) и (47).

Заключение. В классической теории термоупругости законы по­ ведения определяют напряжения двумя упругими константами (X, р) или (.Е, v) и коэффициентом линейного расширения а (все эти вели­ чины предполагаются постоянными); приток теплоты —коэффициен­ том теплопроводности k в законе Фурье, который также предпола­ гается постоянным. Если, кроме того, известна функция Се(Т), задающая удельную теплоемкость при постоянной деформации, то можно получить выражение для термодинамических функций е, ф, s.

Замечания. 1°. Так как функция h (Г) не влияет на законы поведения, то зависимости (46) можно было бы получить непосредственно, исходя из выражения для ф в виде полинома второй степени переменных е,у и 0. Можно было бы,

например,

принять

h (Г) =

21о■Q2 при постоянном Се .

 

 

Однако

-7р- и 8/у имеют разную физическую природу, и именно поэтому пред­

ставляется

ке

 

 

теорию, которая лучше

выявляла бы зависи­

интересным построить

мость термодинамических

величин

от 0 и, в частности,

дала

бы более точные

приближения для функций

е, ф, s,

фигурирующих в уравнении

энергии которые

входят в решения задач классической теории термоупругости.

 

2°. Проведенные

выше

рассуждения имеют к тому же другое достоинство —

они открывают еще

один

путь использования законов термодинамики для описа­

ния поведения физической системы. Первый путь, в основном дедуктивный, пред­ полагает известным термодинамический потенциал (в частности, именно это было сделано в VII 1.1 и V III.2); зная термодинамический потенциал, можно построить законы состояния. Практически именно законы состояния и в более общем случае законы поведения легче проверяются в экспериментах. Предполагая эти законы известными например (46), используют затем законы термодинамики для опреде­ ления термодинамических потенциалов и получения дополнительной информации. Именно этот второй путь был использован выше.

Можно провести аналогию между использованием здесь начал термодинамики и двойственным применением фундаментального закона в классической динамике. Из наблюдения некоторых видов движений можно сформулировать законы о силах — пример теории гравитации в этом отношении особенно впечатляет. Если же известны

силы,

приложенные к

некоторой системе,

и распределение масс в ней, то основ­

ной

закон позволяет описать движение. Аналогия заключается в следующем.

Если

экспериментально

известны

законы

поведения, то из них могут быть най­

дены

термодинамические

функции.

Если

же

известны термодинамические функ­

ции,

то из них выводят

законы

поведения,

и затем изучается движение среды

с учетом заданных условий.

 

 

на простом примере различие между

3°. Представляется

интересным показать

изотермическими и адиабатными процессами. Если ограничиться для простоты

функцией h (Т) вида h (Т) =

(как было упомянуто выше), то

уравнение

 

2То

 

 

 

 

(48) запишется в такой форме:

 

 

 

 

 

p0s — 3Kci&kk

1о

 

 

 

 

 

 

Для адиабатного

процесса, когда

s =

0, уравнение (46) примет вид

 

° i j ~

Н—9* ос е--- ^

eftfe^/7_t“^Iie^7==^ efefe^7+2pe/y.

(52)

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

Читатель получит более детальное представление об этой проблеме, разбирая задачи 13 и 14. Из них станет видно, что модули упругости в адиабатных про­ цессах просто выражаются через модули упругости изотермических процессов и отношение Се /Са удельных теплоемкостей соответственно при постоянных дефор­ мациях и постоянных напряжениях.

 

 

VIII.6. ВЯЗКОУПРУГИЕ СРЕДЫ КЕЛЬВИНА — ФОЙХТА

что

За отправную точку

примем те же термодинамические гипотезы,

и в V III.2.1, с той

лишь разницей, что внутренняя диссипация

не

будет

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

поненты

Laр тензора деформаций Грина —Лагранжа относительно

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

нием плотности

свободной энергии ф (71, LaP). Если saP—матрица

напряжений Пиола —Кирхгофа, то

уравнение

(23) можно переписать

в таком виде:

 

 

 

 

 

 

 

 

ро (Sap

Po^afO^afb

 

 

 

где р0—плотность

в естественном

состоянии;

Ра$-

дф

- в соот-

dLa

 

 

 

 

 

р

ветствии с (17),

а

LaP—более простое обозначение

субстанциональ­

ной производной

 

d^ !L

 

 

 

 

Это выражение аналогично форме 2Х аУа, причем в качестве аргу­

ментов Уа

выступают

переменные Lap — скорости деформаций. Пред­

положим,

что

механизм диссипации является

нормальным (VII.3.2)

и что для упрощения функция диссипации

(Yа) — квадратичная

форма переменных Ya

 

 

 

 

 

^арХц^ар

(53)

коэффициенты

СаР^

которой удовлетворяют

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

метрии:

 

СаРЛ|х — СраЯ.д — СарцД— W^iap*

 

 

Тогда из общей теории можно вывести следующий закон поведения:

 

SccP

Ро

+^арХц

ju*4

»

 

(54)

который дает непосредственно выражение тензора

Пиола —Кирх­

гофа

в виде суммы двух слагаемых, первое из

которых можно наз­

вать

упругой составляющей,

а второе — вязкой

составляющей

этого

тензора. Заметим, что

коэффициенты Са^

могут

зависеть

от Т

и LaP. Наконец, чтобы замкнуть систему и определить для нее допустимые термодинамические процессы, следует написать закон Фурье, связывающий приток теплоты с градиентом температуры

(VII.4.2).

Ограничимся для простоты линеаризованным вариантом предыду­ щей теории, выбирая в качестве исходной конфигурации естествен­ ное состояние, свободное от напряжений, и считая среду изотропной [как и квадратичную форму (53)], а процессы —изотермическими.

Тогда напряжения о,у (которые можно отождествить с saP) представ­ ляют собой сумму некоторой линейной изотропной функции от е,у

и некоторой линейной изотропной функции от е;/ (которые можно

отождествить с величинами Lap). Обозначая через os и Оц шаровую составляющую и девиатор симметричной матрицы aiJy можно записать результат в виде

a s = 3 K (е* + 0се*),

 

°?i= 2р (е° + Qgtfj),

(55)

где К и р —положительные константы, имеющие размерность напря­ жений; 0С и 0^ — константы, имеющие размерность времени. Здесь нет возможности детально изучить эти законы. Заметим, тем не

менее,

что наличие в правой части производной по времени от дефор­

мации

(в то

время как левая часть —компонент напряжения) озна­

чает наличие явления вязкости. Выкладки,

которые привели к урав­

нению

(42),

позволяют сразу же написать

выражение для диссипа­

тивной

функции

 

 

 

 

 

+

(56)

откуда видно, что 0,

и 0^ —необходимо положительные величины.

И наконец, представляется интересным рассмотреть соотношения,

обратные законам (55). Для конкретности

остановимся на первом

из них. Имеем соотношение

 

 

°'«)«Ч > ( £ ) = З К { е * (0 + в > (0 |е х р ( У -

 

 

З К 9 ,-|(е -(0 в х р ( £ ) } ,

из которого

вытекает,

что

 

 

 

 

 

<57>

если предположить, что в моменты / < 0 среда пребывает в своем естественном состоянии. Чтобы получить представление о значении этого вывода, предположим, что сJS (t) сохраняет постоянное значе­ ние а0 для любого положительного /. Тогда

с‘ <о—

t ) } -

<58)

Если вязкость отсутствует (0^ = 0), величина es сразу становится равной ■ £ . Если же 0С=^О, то es стремится асимптотически к этому

же значению при больших значениях t тем быстрее, чем меньше 0С. Точно так же можно рассмотреть и второй из законов (55). Коэф­ фициенты / Сир, естественно назвать длительными модулями упру- гости среды Кельвина, а 0С и 0^ —время запаздывания при равно­ мерном сжатии и сдвиге соответственно.

Законы (55) можно было бы получить, опираясь непосредственно на выраже­ ние для плотности свободной энергии для линеаризованной упругой изотропной

р0г|)= ш (е,7).

Упругая

составляющая

напряжений

dw

 

тогда будет

 

 

 

 

tel/'

 

 

VIII.6. ДВОЙСТВЕННАЯ ФОРМУЛИРОВКА.

 

 

ВЯЗКОУПРУГИЕ СРЕДЫ МАКСВЕЛЛА

 

V III.6.1. Общие

замечания.

При изложении

общей теории

в VIII.2.

предполагалось, что термодинамические

параметры %р,

задающие совместно с температурой Т состояния системы, совпадают с деформациями; в этом случае сопряженные переменные г]р, опре­ деляемые уравнениями (VII.9), имеют смысл механических напряже­

ний и вычисляются по законам

состояния. С другой стороны, участ­

вующие в выражении для внутренней диссипации

Фг =

S

W i

 

х = \

 

величины Yi представляют скорости

деформаций, которые считаются

определенными предысторией,

а величины Х к (типа напряжений)

должны быть определены из дополнительных законов. Двойственный

характер

проблемы «напряжение —деформация», который несколько

раз

подчеркивался, подсказывает мысль

о возможности нового опре­

деления

термодинамически допустимых

процессов для данной сплош­

ной

среды на основе предложенных ранее, в котором следует поме­

нять

ролями

напряжения и деформации. Опишем кратко основные

положения этого нового

подхода.

 

За отправную точку

примем предположение о том, что величины

Ту T]j, . . . , г)я

образуют

полную систему термодинамических пере­

менных и что можно ввести термодинамический потенциал £2(7*,%,

*421 •••* Лл)»

полученный

преобразованием Лежандра из

свободной

энергии

(П II 1.7.4). Более

точно,

положим

 

 

 

 

п

 

 

 

 

 

2

ipf\P—Ф-

 

 

 

 

р=1

 

 

С учетом

(VII, 9) имеем

 

 

 

 

 

 

П

 

 

 

d Q = sd T +

2 X^dV

 

 

 

 

 

Р = 1

 

откуда

с помощью двойственной

формулировки законов

состояния

(VII, 9)

или

(VII, 17) получаем

 

 

 

 

 

dQ /гр

ч

 

 

 

%р~~

дг\р

^

 

Если

внутренняя диссипация имеет вид

 

 

 

 

Фх = 2