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

Механика сплошной среды

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

Назовем энтропией s тела величину, пропорциональную сред­ нему значению по ансамблю логарифма вероятности /:

s(p, 0)= k ( \ n f ( H , 0)) = k f f l n f d p d q ,

(2.14)

Г

 

где k — константа (постоянная Больцмана). Один из параметров 0<? назовем температурой тела 0i и потребуем, чтобы количество тепла 6'<3 имело интегрирующий множитель 0i_1 и с этим множи­ телем представляло полный дифференциал (по параметрам ц, 0) энтропии s :

6'Q =6s.

(2.15)

0i

 

Мы потребовали, по существу, чтобы введенные

количества 6'ф,

5 и 01 удовлетворяли второму закону феноменологической термо­ динамики обратимых процессов, так как (2.15) имеет форму вы­ ражения этого закона. Поэтому вывод может казаться случайным. Но ведь и количеством тепла выражение

J # 6/ dp dq

г

названо тоже из сопоставления формулы (2.12) с выражением первого закона термодинамики. Необходимо еще проверить, при каких условиях справедливо уравнение (2.15), так как количество тепла б'Q уже определено выражением (2.13), вариация 6s может быть найдена из (2.14), параметры 05 еще не определены. В стати­ стической трактовке, т. е. с учетом (2.13) и (2.14), уравнение (2.15) должно быть тождеством, как и уравнение (2.12), выража­ ющее первый закон. Из (2.14) находим вариацию 6s:

6s = — k J б (/ In/) d p d q = — k ^

In/6/ dp dqk ^ 6f dpdq,

г

"

г

г

 

или, учитывая (2.9),

 

 

 

 

 

6s =

—-/2 J In /6/ dp dq.

(2.16)

 

 

r

 

 

Внесем выражения

(2.13),

(2.16) в

(2.15) и перепишем

получен­

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

 

 

 

$ [ * 1и/(Я .

6) + ^ - я ] б / ф ^ = о .

(2.17)

Чтобы для любой вариации б/ этот интеграл равнялся нулю, в со­ ответствии с (2.9) необходима и достаточна независимость от

(р, q) подынтегрального выражения, стоящего в квадратных скоб­ ках. Обозначая это выражение через 02/0ь т. е.

6011п/+Я =02

(2.18)

и усредняя (2.18) по ансамблю, получим

02= Я —0,s.

(2.19)

В макротермодинамике такое выражение известно как выражение свободной энергии Гельмгольца 02^=-ф через внутреннюю энергию Я, энтропию s и абсолютную температуру 01= 7' в градусах Кель­ вина, если в (2.14) множитель k является постоянной Больцмана. Таким образом, уравнение (2.15) действительно выражает второй закон термодинамики, если функция вероятности / имеет выте­ кающее из (2.18) выражение

 

f = er(n-*v*T,

(2.20)

где Т — абсолютная температура и ф> — свободная энергия

 

ф=Я—Ts.

(2.21)

Последняя, конечно,

может быть функцией

внешних параметров

и температуры Т.

 

 

Ансамбль системы 5дг, определяемый функцией распределения

f (2.20), называется

каноническим ансамблем Гиббса. Найдем из

(2.20)

ф=Я+&7Х1п f>.

(2.22)

 

Одна из основных задач статистической термодинамики состо­ ит в нахождении уравнений состояния тела, т. е. связей между внешними силами Рг, внешними параметрами рг и температурой Г, а также в определении энтропии s. Покажем, что если свобод­

ная

энергия известна

как функция цг и Т, т. е. ф=ф(ц,

Т), где,

как

и прежде, ц —

совокупность (р.ь ц2........цг, ...), то

уравне­

ния состояния и энтропия вполне определены. Для этого перепи­ шем уравнения (2.12), (2.15), выражающие два основных закона термодинамики, в виде

6Я = £ Pr&Pr+ S'Q,

Г

(2.23)

T8s=8'Q

и используем их следствие

b H - n s = Y .P r6iir.

Г

Внося в это уравнение выражение 6#, вытекающее из (2.21), по­ лучим основное термодинамическое соотношение

6vj) + sbT= ^ Рг6|лг.

(2.24)

Отсюда получаем выражение для энтропии и уравнения состоя­ ния системы:

s = —дЩ&Г,

(r= 1, 2, 3, ...).

(2.25>

Теплоемкостью системы при постоянных значениях парамет­ ров Цг называется количество тепла, необходимое для нагревания единицы ее массы на один градус Кельвина:

6'Q

1

дН

Cv — M N6T

— M N

дт -

что при pr=const следует из (2.23); при этом MN обозначает мас­ су системы. На основании (2.21), (2.25) получаем

М ^ = Т — = — Т

(2.26)

* дТ

дТ2

 

Итак, основная задача сводится к определению свободной энер­ гии ф(ц, Т) по заданной функции Гамильтона системы Н (р, q, ц). Для конкретных физических сред это почти всегда трудная за­ дача. Выражение ф через Н в виде функции интеграла состояний Z(p, Г), определяемого интегралом по всей фазовой области Г формулой

Z(p,

T) = \e - ^ P ’« ^ kTdpdq,

(2.27)

 

г

 

находится просто из условия (2.3) нормировки функции f

(2.20)

§ e - u i - w / k T dp d q = e w k T z = 1 .

 

Г

 

 

Отсюда

 

 

ф(^,

Т)= — kT\nZ(\i, Т).

(2.28)

Для простых систем, т. е. состоящих из частиц, которые доста­ точно точно представимы материальными точками, выражение ин­ теграла Z несколько упрощается, так как для них

n = 3 N 2

t=1

Интеграл состояний представляется в виде произведения двух ин­ тегралов: одного — по области Гр от функции, зависящей только от 9, который обозначим

D (ц, Т) = (

 

 

dq)V3N

(2.29)

и другого —

 

 

 

 

С -

1

3N

 

 

IFF

*=i

dp,

(2.30)

Ц Т ) = \ е

 

г„

 

 

 

 

причем, конечно,

 

 

 

 

dq=dq1 dq2 ... dqn,

 

dp=dpxdp2

dpn,

ti— ZN.

 

Обычно предполагается, что область ГР изменения импульсов бес­ конечна, т. е. каждый импульс pi изменяется от —оо до +оо. Поэтому интеграл I представляет произведение п=ЗЛГ интегралов вида

/ г(Г )=

J е

р2*/2кТтdpi=V2nkTm,

 

 

'—оо

 

 

 

т. е. равен

 

 

 

 

I (T)=(BjT)3N/2, B 1 = 2nkm.

(2.31)

Таким образом,

 

 

 

 

Z

T) =

[B1rD(ii, Г)]3" '2.

(2.32')

и потому для простых систем

 

 

 

4i = — 2- т г [ 1пГ + 1пО(ц, D + i n s j ,

(2.32)

 

Я=4>—

дТ

(2.33)

 

 

 

 

В качестве простейшего примера физической равновесной си­ стемы рассмотрим разреженный одноатомный газ, заключенный в объеме V. Атомы на больших расстояниях не взаимодействуют, имеют пути пробега в пределах всего объема V и кратковремен­ но взаимодействуют только при «столкновениях». Поэтому такой

газ можно рассматривать как идеальный, т. е. потенциальную энергию U положить равной нулю. Из (2.29)

D3N/2= J dq= J dqx dq2 dq3,

dq3N= ( i\ dql dq2 dq3)N = VN,

r<i

r<

*

так как в произведения dq каждое последовательно одно за дру­ гим расположенное произведение трех дифференциалов декарто­ вых координат различных точек образует дифференциал объе­ ма dV

dV= dqLdq2 dqs= dq4 dqb dqQ= ..

и все координаты изменяются в пределах всего объема V Следо­ вательно,

D=V2i\

(2.34)

причем V — единственный внешний параметр системы. Из (2.33) находим свободную и внутреннюю энергии газа

г|>=----N kT J In Т

In К + constJ,

(2.35)

Н - — NkT.

2

Поскольку для рассматриваемого газа

k = \

то 3kT/2 есть средняя кинетическая энергия одного атома

~2

_3

 

 

mvk

kT

(2.36)

2

2

 

 

Внешней силой для нашего объема V газа является давление —р (р> 0), так как работа бМ, сообщаемая газу, равна

б'A= —pdV

Из (2.25) и (2.35) находим выражение энтропии и уравнение со­ стояния:

s— Nk [In (VT3/2) -|- const],

 

р = —

dV

— , pV— NkT = RT

(2.37)

H

V

 

Теплоемкость при постоянном объеме (при M^=Nm) на основа­ нии (2.26) равна cv=3kj (2т).

Неточность уравнения состояния (2.37) увеличивается с уве­ личением давления р. В табл. 2 дано значение pV/RT для тех же газов при различных р. Поправка для больших давлений и тем-

 

 

 

Т а б л и ц а 2

Вещество

р = 1бар

р = 50 бар

р = 1000 бар

Не

1,0005

1,24

1,44

Аг

0,9999

0,971

1,675

N 2

0,9998

0,996

1,99

ператур может быть найдена при учете потенциала взаимодейст­ вия молекул u(q); тогда получается уравнение состояния

РУ

= 1

, В ( Т )

С ( Т )

R T

 

V

V2 _г

причем коэффициенты В(Т), С(Т) определяются путем вычисле­ ния интеграла состояний Z.

Из (2.36) найдем характерную скорость атома

vc— V 't fk = КМ-с-1,

где М — отношение массы молекулы вещества к массе атома во­ дорода (точнее 1/16 атома кислорода).

Используем еще некоторые другие формулы для средних вели­ чин, вытекающие из кинетической теории газа. Длина I свобод­ ного пробега атома, имеющего эффективный диаметр dy выража­ ется через скорость vc, число молекул в единице объема N/V и число столкновений одной молекулы в единицу времени пст:

/ = - ^ ,

nCT= ^ p - N d \ .

Пет

V

Сведения о некоторых газах даны в табл. 3 для 7’=300 К.

 

 

 

 

 

Т а б л и ц а 3

Вещество

м

о

о

пст* С1

Vc, км с 1

(1, А

/, А

Не

4,0

2,2

1936

6,5*109

1,26-

Аг

39,9

3,6

694

5 ,7 -10е

0,40

N 2

28,0

3,7

654

7,3*10®

0,475

Геометрия монокристалла определяется относительным положе­ нием атомов, точнее, их ядер, в элементарной кристаллической решетке.

Потенциальную энергию V (q) электрически нейтральной си­ стемы частиц представляют, например, в виде суммы потенциалов парного взаимодействия частиц типа потенциала «6—12».

 

^ ( ? ) = Е и (|г*а 1).

 

ft..*.

где г*,*

— радиус-вектор, соединяющий ki-ю и А2-ю частицы

г*

—г* - Обозначая г*°, q* начальный постоянный радиус-век­

тор k-n

частицы и вектор ее смещения относительно точки г*°, так

что г*=Г/г0+Я*, получим U(q) в виде функции всех qi и начальных координат точек. Для идеального кристаллического тела с перио­ дической решеткой при малых колебаниях атомов около положе­ ния равновесия (<7;= 0), в котором U минимально

( |L ) = 0 « , 1 = 1 , 2 . n=3N), (2.38)

ограничимся для малых деформаций квадратичными членами раз­ ложения U в ряд по q:

f / = ^ 0+ 4 - £

(2-39)

i,/=l

 

Постоянные взаимодействия сц и начальная энергия U0 в со­ стоянии qi= 0 зависят от параметров этого состояния ц, например, деформации решетки. Интеграл состояний Z в этом случае имеет вид (2.32'), причем D можно вычислить, допуская, что координа­ ты qi изменяются в бесконечной области (—оо<<7/<оо). При этом условие /=0 на границах Гр и Г? выполняется, так как квадратич­ ная форма (2.39) положительна.

Вычисление интеграла D (2.29) как функции констант U0 и сц может быть выполнено приведением входящей в (2.39) квад­ ратичной формы к сумме квадратов или на основе упрощений мат­ рицы Ik//1| в случае парных потенциалов. Более сложный вопрос уточнения физического смысла параметров ц(|Дг) как характери­ стик деформации кристаллической решетки монокристалла в рав­ новесном состоянии (qi=0) в простейшем случае решается на основе понятия аффинной деформации окрестности точки сплош­ ной среды (гл. II).

Поэтому рассмотрим одномерную модель монокристалла в ви­ де линейной цепочки п атомов одинаковой массы т , расположен­ ных вдоль оси х на одинаковом расстоянии а в нулевом равновес­ ном состоянии Sn° и на расстоянии \ха в деформированном равно­

весном состоянии SnK Для определенности введем границы систе­ мы: неподвижный нулевой атом с координатой л:0=0 и подвижный (/г+1)-й атом с координатой хп+и равной x°n+i= (л + 1)а в состоя­ нии Sn°, ххп+\=(п+ 1)ра — в состоянии S nu>взаимодействие перво­ го и л-го атома системы с граничными пусть будет для простоты таким же, как с соседними атомам# системы; именно из этих гра­ ничных условий следует равномерное в Sn1 распределение удлине­

ния цепочки, т.

е. расстояния между границами

х 1п+\—*°Л+1=

= (ц—1) ( я + 1)а,

по (л + 1) энергетическим ячейкам,

следователь­

но, и положение m-го атома Хтп=тп\1 а в Snl.

Статистическое состояние Sr в момент t определяется коорди­

натами X m ( t ) = X m l + qm(t)

для

Ш = 1,

2, ...,

Л, (^ о ( 0 =

0 у Qn+\(t) =

= 0). Потенциал

парного

взаимодействия

m

и

пг+1

атомов

и{Гт,т-\-\) При Т

m+\ —Хт+\

-J-А<7т, (ТП=0,

1

...,

л)

И A (Jm=

=Цт+1—Цт приводит к представлению 0 в виде

 

 

 

 

п

й ^Г т 'т +

0 =

( " +

i ) И (|ха) +

 

 

и "

 

V — 2

 

 

 

 

 

 

 

 

(ца)

 

 

 

 

т=0

 

 

 

 

 

т=0

 

(2.40)

 

 

 

 

 

 

 

 

 

и" (|ха)

йги (г) I dr2 |г=ца

Для вычисления интеграла D введем вместо цт новые перемен­ ные, поскольку и"(ра)> 0,

1т=Ят у

А ^ = |ш +1-

| т .

 

Получим из (2.29) (при я » 1)

 

 

 

 

Dip, Т) — А —^£-— ехр [ — 2ИИ£)_\

(2.41)

и " ( ц а )

F \

kT

) ’

v

где А — число. Свободную энергию, приходящуюся на один атом, согласно (2.32), и ее производную по р, т. е. растягивающую си­ лу, находим из формул

Ч»1(Ц Т)= ty/ti= и (pa) + -L kT In иГ(pa) — kT 1n 7,

(2.42)

Р„,— - ^ - = аи' (pa) + kT

2и"

аш) .

ар

(ра)

Первое слагаемое в выражении силы определяет ее при Т=0, об-

ращаясь в ноль при р = 1. Второе дает поправку на тепловые коле­ бания атомов при Т > 0.

В области устойчивых состояний цепочки «"(ра)>0. При и"'(ра)< 0 нагревание системы уменьшает силу, необходимую для удерживания деформаций е=р—1> 0.

На примере потенциала «6—12», для которого в (2.42)

аи' (ра) — 12kQy(р~7—р."13)=726в1е + О(е2),

аи!" (ра)/а* (ра)= 14 4fl~9 ^ V~‘5.v ,

видим, что для выявления тонких термомеханических эффектов представление U квадратичной формой (2.39) недостаточно: мо­ дуль упругости цепочки зависит от Т через и'"(г), которая отбро­ шена в (2.40).

КИНЕМАТИКА И ВНУТРЕННИЕ НАПРЯЖЕНИЯ

§ 3. ОСНОВНЫЕ ПОНЯТИЯ и УРАВНЕНИЯ МСС. ЛАГРАНЖЕВО И ЭЙЛЕРОВО ПРЕДСТАВЛЕНИЯ ДВИЖЕНИЯ

Будем рассматривать настолько большие тела, что весьма ма­ лые их части объема dV содержат достаточно много частиц, и по­ тому для этих малых областей тела можно ввести понятия макро­ скопических величин плотности тела, перемещения, скорости, уско­ рения, внешних сил, внутренней энергии и других в смысле сред­ них по ансамблю (§ 1). Идеализация истинного физического те­ ла в МСС состоит в том, что все рассматриваемые величины при­ нимают в качестве истинных. Количество и математическая при­ рода вводимых средних величин таковы, что с достаточной точно­ стью можно описать внутреннее состояние тела и взаимодействие между телами. В основах МСС главным образом рассматривают­ ся механические и тепловые взаимодействия и деформации малых объемов, иногда учитывается действие на них электромагнитных полей, химических реакций и др.

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

Основные математические объекты МСС суть тензоры различ­ ных порядков: нулевого — скаляры (плотность, энергия), перво­ го — векторы (радиус-вектор, поток тепла, скорость), второго — тензоры деформаций, внутренних напряжений, третьего и четвер­ того — тензоры пьезоэлектрических констант, коэффициентов вяз­ кости и упругости и др. Все эти тензоры считаются непрерывно дифференцируемыми достаточное число раз по координатам и по времени, ограничены вместе с их производными в области тела. Все они введены в XIX веке в процессе создания теории упруго­ сти, гидромеханики и других разделов теоретической физики, и затем в алгебре и геометрии была создана их общая теория.

Определения основных величин в МСС можно рассматривать и как априорные, и тогда получающиеся уравнения, подобно за­ кону Ньютона для точки mx=Fx, станут содержательными зако­ нами только на основании экспериментально обоснованных гипо-