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

книги / Релаксационные явления в полимерах

..pdf
Скачиваний:
3
Добавлен:
13.11.2023
Размер:
23.13 Mб
Скачать

функций распределения, которое лучше всего отвечают экспери­ ментальным данным. Эти функции могут иметь по крайней мере четыре характерных вида (рис. 11.33) [38].

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

Вслучае смеси включение — включение в обоих фазах независимо друг от друга обра­ зуется плоская ортогональная структура си­

ловых линий, которая в значительной сте­

 

Рис. 11.33. Характерные

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

 

виды функций распреде­

вательной системы.

 

 

 

 

 

ления для смесей различ­

Уравнения (11.101)

и

(II. 103), отвечаю­

 

ного типа:

 

для функции h (ф)

щие последовательной

и параллельной об­

 

1 — среда— среда; 2— включе­

общенным моделям, могут быть решены ме­

 

ние-включение; а—среда —

 

включение;

4— включение —

тодом численного

 

интегрирования. Если

 

 

среда

разбить фл на п мелких участков Афд, t-, и,

 

для функции f (Ф)

 

/ — включение — включение;

придавая

Дфл.г

произвольные значения,

 

2—среда —среда; 3—среда —

взять Е*(Т), Ê*(T,q>Aii)

при п— 2 выбран­

 

включение;

4—включение —

 

 

среда

ных температурах

 

Т\

 

Тп- 2 то ЦфА,0

 

 

 

и Я(фа, i)

будут выражаться

следующими матричными соотноше

ниями:

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

f l '

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

Г,)

 

h

 

 

Г

(Tt)

=

X

h

(II. 104)

 

. Г

(Тп-2)

.

(фА V Тп -2 )

 

. к

 

 

 

 

1

 

 

1

 

V

 

 

 

 

V

 

 

Ч>Л,1

 

Я

 

 

 

г

(ТО

 

£ (ФЛ,!> Г1)

X

Яз

(II. 105)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

Е* (Тп- 2)

Е (ФД, V Tn~i)

 

.V

 

 

 

 

 

 

 

 

 

 

При этом %i = А.(фА,0ДфА,* Уравнение (11.105) представляет собой систему линейных алгебраических уравнений с п неизвест­ ными %i. Решение системы этих уравнений осуществляется до­ вольно просто, но при возрастании п этот процесс становится весьма трудоемким. На основе температурных зависимостей дей­ ствительной части Ei комплексного модуля упругости смеси Е* путем использования уравнения (II.104) могут быть определены численные значения функции распределения МфА.О для смеси

ПС с СКС-30. При этом можно взять /г — 6. Ниже приведены результаты расчета <рл и \(фд, <) для этой смеси с различной величиной объемной доли каучука (VA) :

<P.4,i

Мфл, i) •

МФл. i) •

ЭДфЛ, Ù

Мфл, 1)

0,001

0,01

0,1

 

0.3

 

0,6-

2,2-Ю -3

 

КЛ-0,8

сл о

о to

2,0 -10-2

2,0-Ю -2

2,0 - 1 0 -3

 

 

 

 

 

1

 

 

1,0-10-'

 

V„-0.6

4,9-10-2

 

 

1,2. ю - 1

1,2- ю - '

 

CN 1 О

 

 

V A -* *

 

1

 

1

2,1.10-' 2,0.10“ '

5,0.10—2

4* О

О

О to

О to

 

 

Кл-0.4

3.2 - 10 -3

1,5* Ю -2

4,2.10-'

1,1.10-'

5,0-10-3

1,0

7,6-10-'

со ю

7 о

4,5-10-'

ю

7 о

Значения Цфл, «) достигают максимума при крайних значе­

ниях Фа, г, что указывает на

распределение типа среда — среда.

Аналогичный расчет был

сделан для смеси

ПБМА с

ПММА

(в соотношении 1:1)

при VA =

0,5 и п — 6.

 

 

 

<р„ I .

0,1

0,3

0,5

0,6

0,9

1,0

МФл, f) *

1,6 •ю -2

5,0•10-*

6,0 •ю -2 8,0 •10~2

1,0•ю -' 2,2 •10-'

Полученные данные говорят о том, что эта смесь относится к типу В в А (включение — среда), что дополнительно подтверж­ дается также результатами микроскопических исследований.

Представляет интерес расчет Я,(фд, <) для смесей аморфного и частично кристаллического полимера и смесей двух частично кристаллических полимеров, например для смеси ПЭ высокого давления с ПИБ (при Ра = 0,5 и п — 6) Я/(фл, <) и для смеси ПЭ с ПП (при тех же значениях VA и п) Х"(фа,{);

ФЛ , .

0,1

0,3

0,5

0,6

0,9

1,0

Л'(фA l )

1,5-ю - 2

1,8 •10-2

2,0 - ю - 2 2,1

10-2

2.2-10-' 2,5-10-*

Л"(Фл t )

1,5- 10-2

1,8-10-2

2,0- 10-2 2,1

•10—2

2.2.10-' 2,5-10-'

Для смеси

ПЭ с ПП функция А'(ф а , i)

имеет

минимум при

Фа , i — 0,5 и

возрастает

при крайних значениях

фд, <• Следова­

тельно, эта смесь образует систему двух взаимопроникающих кар­ касов (среда — среда).

Для качественного описания поведения взаимонерастворимых смесей полимеров может быть предложена механическая модель, представляющая собой последовательное включение двух моделей Кельвина с элементами ЕА ца и Ев TJB и пружины с моду­ лем До (рис, 11.34). Дифференциальное уравнение такой модели имеет вид:

ё +

Е .Е Я

 

■ л~л é =

 

 

+

е а е в ) а

 

£о П а^в

 

 

Решение его для случая периодического режима нагружения может быть записано следующим образом:

Емп — Е'

(1+а^тв)

+©’ (xA-tB) [ха{х+ТГ^+Хв (l+'I j) ]

Е' =

 

 

 

 

 

 

 

 

 

 

 

 

(l+-f i +‘f r - <e’tH

’+

 

 

 

 

 

 

 

 

 

 

 

+Ш’ Ь

(,+-ft )+4

1+-ë -)]’

 

 

 

 

 

 

 

 

tgô =

E " IE '

 

 

 

(И. 106)

 

 

 

 

 

 

 

Е" = СjD

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С =

Е0Ф (хА -

rB) (1 + Е0/Еа + Е0/Ев -

©V

B) +

 

 

 

 

+

» ( 1 “ ^

X A X B ) + ® 0

тхАх в ) +

 

 

 

 

 

 

+ M

I +

W

+ M 1+ W

]

 

 

 

 

c -

( '

+

t

+

t - ‘“W

+

 

 

 

 

 

 

 

 

 

+ “ ’ h ( ' + ' ï ) + , ” ( ' + 'ft )]

Рис. 11.34. Ме­

 

 

 

 

 

 

 

 

 

 

 

 

Из

выражёния (IL 106)

следует, что

характер

ханическая мо­

 

дель для описа­

зависимости

Ета и

tgô

от частоты

согласуется

ния

релакса­

с наблюдаемым экспериментально. Таким образом,

ционных явле-

можно считать, что эта модель качественно удов­

ний во взаимно­

нерастворимых

летворительно

описывает

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

релакса­

смесях

полиме­

ционные явления во взаимонерастворимых полимер­

ров.

ных системах.

 

 

 

 

 

влияние распределения

 

Вывод соотношений, характеризующих

компонентов полимерной композиции на ее механические свойства. Для композиций полимеров используются обобщенные механиче­ ские модели простейшего (рис. 11.31) и усложненного (рис. 11.32) видов. При этом в случае полимерных композиций важно оценить распределение объемных долей компонентов фл и фв, влияющее на динамически механические величины композиций и характер их температурно-частотных зависимостей [39]. Для t-ro элемента простейшей последовательной обобщенной модели (рис. II. 31, о) комплексный модуль упругости Юнга может быть представлен в виде:

 

к *

ш(1)

. ФК А

,

Ф(/ в

 

 

 

$

Ф*. В

. ^

«

Г

Т

 

Ь\

*

г

“Г

 

 

 

 

£Л +

*3

J

или

/

Т-1

г I

 

 

 

 

 

 

IS »'Л>

s

 

*

 

 

 

 

 

 

 

 

 

Щ -

j=1

J=1____

 

 

 

 

 

 

 

 

 

 

“А

 

 

“В

J

 

 

 

 

Учитывая,

 

 

что

 

 

+

I

 

 

окончательно

получим

 

 

/=1

у=1

'

 

 

 

 

,

 

Юнга формулу (11.96). Для

для комплексного

модуля упругости

t-го элемента простейшей параллельной

 

обобщенной

модели

(рис. II. 31, б)

 

 

комплексный

модуль

упругости

Юнга может быть

записан в виде:

 

 

 

 

 

 

 

 

 

 

 

 

£;, = <р1!Ул+ ф(!)в^ +

 

+ фм£и+ < ’в£в +

 

+ фд + ч,№

Группируя

 

члены

с

ЕА и Ев с учетом

того,

что

+

+ <?\т)А=

 

 

т

 

 

 

 

 

у[т)в =

 

т

 

 

 

 

2 Ф(^’Л

И ф<'>в +

+

 

2

ф'/V

получим:

*

fe=l

 

 

*

 

 

*

k=l

*

 

 

 

 

 

 

 

 

 

т

 

 

т

 

 

 

 

 

 

 

 

 

* !- * л 2 ч > ? л + * а 2 ф * в

 

 

 

 

 

 

 

 

 

k=l

 

 

k=\

 

 

 

 

 

Учитывая,

 

что

2

ф)*^ =

Ф, А и 2

ф£*в = ф< в,

окончательно

получим для

комплексного

модуля упругости Юнга формулу:

 

 

 

 

 

 

Ф<, А&А+

Ф<, В^В

 

 

 

 

 

Для усложненной

последовательной

 

обобщенной

модели

(рис. II.-32, а)

 

обратная величина

комплексного модуля

упругости

Юнга (механическая податливость)

может быть записана в виде:

 

 

 

 

 

 

k

 

 

 

k

 

 

 

 

 

 

 

 

 

 

=

 

~Щ‘

где

S

 

 

 

 

При этом

 

 

 

 

 

i= 1

 

1=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

П

 

 

 

 

 

 

 

 

 

2

 

 

hb,A = VA>

2

 

 

ф1 А + Ъ,в = 1

 

i=\

 

 

 

1=1

 

 

 

 

 

для случая п -* оо,

Принимая во внимание выражения

(II. 98)

окончательно получим интегральное соотношение

 

 

 

 

 

 

 

 

 

КФл).

 

 

J-

 

 

(II. 107)

 

 

 

 

 

 

Фл^ + О+ФлИв

 

А

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

которое, с учетом того,

что J* /Чфд)dq>A=

1

и

J f (фЛ)q>A dyA= V А

 

 

 

 

 

 

 

о

 

 

 

 

 

о

 

 

переходит в формулу (II. 101), если обозначить уАЕ'А + (1 — фА)ЕВ

через £‘ (фл)*

Для

усложненной

параллельной

обобщенной

модели

(рис. II. 32, б) величина

комплексного

модуля упругости

Юнга

может быть представлена в виде:

 

 

 

 

 

k

k

 

 

 

£ * =

2 ^ i . где

2JXf — I

 

 

i=I

i= l

 

 

При этом

кk

2

A = ^A' 2

B= ^B’ 4P/, B"*■ B = 1

Принимая во внимание выражения (11.98) для случая п-*- оо, окончательно получим интегральное соотношение:

£* =

(II. 108)

1

1

которое, с учетом того, что J* Я(фд)^фЛ— 1

и J Я(фд)фд^фд=Уд,

о

о

переходит в формулу (II. 102), если обозначить (фд/£д+(1—1Фд)/£д) через £*(фд). Уравнения (II. 107) и (II. 108) могут быть решены методом численного интегрирования. Переписав их в более ком­ пактном виде, получим соответственно:

 

/(«V

 

 

1

 

 

 

 

<*ФА и

Е' — J Е*"(фд) Я (фд)d<pA

 

 

*4> л )

 

 

о

 

 

 

Далее, разбив фд на п мелких участков Дфд,< и придавая

Дфд, i произвольные

значения,

можно

взять

Е * ( Т ) , Е * ( Т ,

фВ. i)

при и — 2

выбранных температурах

Т\ . . . .

Т п - 2 ,

получая

для

Я(фд, 0 соответствующее матричное уравнение:

 

 

 

для усложненной

последовательной обобщенной модели

 

1

1

 

1

 

 

 

 

fl

 

 

ФЛ,2

фЛ .я

 

V

Фл.1

 

 

h

/ • ( Г , )

'*(Ф л .Ь

T l)

'* (Ф Л ,2 'Tl)

Т*(Фл> п' ^ l)

h

 

Л. V T2)

 

А ,Г Т2)

 

(Фл./»’

X

h

П Ъ )

ш

 

г г)

 

 

 

 

 

 

 

 

Г ( Т п - 2)

У Ч Ф Л ..- Тп - 2) ^ ( Ф л . 2 - ^ - 2 ) П Ф Л .П -Тп - 2) .

in

 

 

 

 

 

 

 

(II. 109)

где fi / (фд, t) Дфл,

1

 

1

 

1

 

 

 

 

 

 

 

 

V

 

Фл,2

Фл,2

 

Е'(Т>)

£

(<Ра . 1 . 7’ |)

£

(фЛ>2'

^ l)

Г

(Т2)

%

л, . . 7’2)

£

(Фл, 2’

тг)

P

(Tn-t)

* ( * Л ,> Л л )

Е (Ча , 2 ' Тп-'.

 

 

1

1

V

 

 

Ф л .в

 

 

* ( Ф Л , п' Т \)

X

Яз

Е (Ф л , » ' Г г )

и

 

 

 

 

 

Х п .

Е (Ф л . п' ^ п - г )

(И. по)

где Xi = Я(фA,i)Афл, V

Уравнения (11.109) и (11.110) представляют собой системы линейных алгебраических уравнений с п неизвестными Я;. В слу­ чае небольших значений п решение их не вызывает особых труд­ ностей. При достаточно больших п решение этих уравнений может производиться с помощью электроннных вычислительных машин.

ЛИТЕРАТУРА

 

 

 

1. И. Пригожнн, Введение

в термодинамику необратимых процессов, ИЛ,

I960.

 

 

Основы моле­

2. И. Г. М и х а й л о в, В. А.-Соловьев, Ю. П. Сырников,

кулярной акустики, Изд. «Наука», 1964.

гл. 2,

Изд. «Мир»,

3. Физическая акустика, под ред. У. П. Мэзона, т. И,

1960.

свойства высокополимеров,

ИЛ,

I960.

4. Т. Алфрей, Механические

5.Т. Алфрей, Е. Ф. Гарни, Динамика вязкоупругого поведения, в кн/: «Реология», под ред. Ф. Эйриха, ЙЛ, 1962.

6. Д ж. Ферр и, Вязкоупругие свойства полимеров, ИЛ, 1963.

7.Физическая акустика, под ред. У. П. Мэзона, т. II б, гл. 1—3, Изд. «Мир», 1969.

8. Г. Л. Слонимский, ЖТФ, 9, 1791 (1939); В. А. Каргин, Г. Л. Сло­ нимский, 11,‘341 (1941).

9.В. Gross, Mathematical Structure of the Theories of Viscoelasticity, Paris, 1953.

10.А. В. Тобольский, Свойства и структура полимеров, Изд. «Химия», 1964.

11.Viscoelasticity; Phenomenological Aspects, Ed. J. T. Bergen, N. Y., 1960.

.12. Я. И. Френкель, 10. H. Образцов, ЖЭТФ,«9, 1081 (1939).

13.Г. Л. Слонимский, 13, № 2, 450 (1971).

14.Ю. В. Зеленев, Автореф. канд. дисс., МГПИ им. В. И. Ленина, 1964.

15. G. М. Bartenev, Yu. V. Zelenev, Plaste u: Kaut., 18, № 1, 16 (1971).

16.А. П. Александров, Ю. С. Лазуркин, ЖТФ, 9, 1261 (1939).

17.Г. M. Бартенев, Колл, ж., 12, 408 (1950).

18.Ю. С. Лазуркин, Автореф. докт. дисс., ИФП АН СССР, 1953.

19.Г. И. Гуревич, ЖТФ, 17, 1491 (1947).

20.Г. М. Бартенев, ВМС, 6, 2155 (1964).

21. А. П. Молотков, Ю. В. Зеленев, Г. М. Бартенев, Укр. физ. ж., 12, Кя 2, 303 (1967).

22. Ю. В. Зеленев, А. П. Молотков, Ученые записки МОПИ им.

Н. К- Крупской, сер. физ., вып. 8, 147, 151 (1964).

23.Ю. В. 3 ел е н е в, сб. «Надежность и долговечность полимерных материалов

и изделий из них», Изд. МДНТП, 1969; Жури. НТО СССР, Кя 12, 9 (1969).

24.А. П. Молотков, Ю. В. Зеленев, ВМС, 10, 5 (1968).

9 ; Д „ Ж У Р к ов, Неорганические материалы, 3, № 10, 1767 (1967); Вестник

ос fnH gc Çp, № 3, 46 (1968).

97 тп п ! еленев- В- И. Абрамова, ВМС, АН, № 4, 920 (1969).

л- tu. в. Зеленев, В. И. Абрамова, О. И. Гроссман, Л. М. Элект*

ne

R0B3

|МС, А10, 1585 (1968).

 

полимеров, 5,

№ 1, 30

 

*; *•(. Бартенев,

10. В. Зеленев, Механика

 

(1969).

 

 

 

 

 

ол

тт

А ск а дек ий,

Физико-химия

полиарилатов,

Изд. «Химия»,

1968.

30.

H. М. S w а 11 w о о d, J. Appl. Phys..

15. 758 (1944).

 

 

31.

Е. G u t h, О. G о 1d, Phys. Rev., 53, 322 (1958).

 

 

32.I. M. Dewey, J. Appl. Phys., 18, 16, 578 (1947).

33.A. Г. Фокин, Г. Д. Шер мер гор, Прикладная механика и техническая физика, 13, 37 (1968).

34.И. А. Чабан, Автореф. канд. дисс., Акуст. ин-т АН СССР, 1963.

35.К. Ni no my a, J. Coll. Sci., 14, 49 (1959).

36.

S. Uemu г а, М. Така у а па gi,

J. Appl. Polymer Sci., 10, 11, 113 (1966).

37.

H. K a v a i, M. O g a v a, Kobunsi,

12, № 10, 752 (1963).

38.A. Б. Айвазов, С. П. Гунькин, Ю. В. Зеленев, Г. M. Бартенев, Проблемы прочности, № 5, 117 (1971).

39.10. В. Зеленев, Автореф. докт. дисс., МОПИ им. Н. К. Крупской, 1971.

Ill

НЕРАВНОВЕСНАЯ ТЕРМОДИНАМИКА И РЕЛАКСАЦИОННЫЕ ЯВЛЕНИЯ В ПОЛИМЕРАХ

ОСНОВЫ СТАТИСТИЧЕСКОЙ ТЕРМОДИНАМИКИ НЕОБРАТИМЫХ ПРОЦЕССОВ

Неравновесная термодинамика образует фундамент теории релаксационных процессов. Здесь кратко изложены основные по­

нятия термодинамики необратимых

процессов. Более подробно

с затронутыми вопросами можно

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

[1-4].

Предварительно полезно очень коротко вспомнить некоторые положения равновесной статистической термодинамики. Основные понятия термодинамики обычно вводятся на примере газа под поршнем. Здесь будет более уместным рассмотреть растяжение по­ лоски несжимаемой резины. Сила растяжения f зависит от длины полоски L и от ее температуры Т. Чтобы изменить длину на малую величину dL, нужно совершить работу 6R = fdL. Внутренняя энергия тела изменяется как при растяжении полоски, так и в ре­ зультате ее нагревания или охлаждения:

d£ = ÔÆ+ ÔQ

где ÔQ— количество тепла, отданное или полученное телом. Коли­ чество теплоты не является функцией состояния, т. е. в нашем случае нельзя сказать, что при данных Т и L величина Q будет иметь какое-то определенное значение. Однако при равновесных (обра­ тимых) процессах можно ввести такую функцию состояния — энтро­ пию 5, что ÔQ = TdS. Энтропия представляет собой важнейшую термодинамическую характеристику. Значение энтропии опреде­ ляется тем, что она служит мерой неупорядоченности. Очень четко это видно на рассматриваемом нами примере. При растяжении ре­ зины происходит ориентация сегментов, т. е. увеличивается упоря­ доченность. Многочисленные эксперименты показывают, что при растяжении резины происходит уменьшение энтропии.

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

теме полностью определяет

ее микросостояние. Если,

например,

наша система — это шарик,

подвешенный на пружине

(одномер­

ный осциллятор), то полное описание состояния этой системы

оз­

начает задание смещения х

и скорости v шарика. Состояние

oç-

циллятора можно изобразить

некоторой точкой

(см. рисунок).

В процессе движения точка на

плоскости х Ь

описывает неко­

торую кривую. Для гармонического осциллятора это будет эллипс. Ангармонический осциллятор будет давать более сложную кри­ вую. Чаще вместо координат х — о пользуются координатами х и р, где р — mv — количество движения или импульс. Состояние осциллятора в-координатах х р изображено точкой Г (см. рису­ нок). Движение более сложных систем будет описываться кривы­ ми в пространствах большего числа измерений. Пространство

положений и импульсов всех частиц

 

системы называется фазовым

прос­

 

транством

этой системы.

 

 

 

Точка в этом пространстве, ко­

 

торая изображает некоторое состоя­

 

ние системы, т. е. показывает, какие

 

импульсы

имеют

частицы

системы

 

и где они находятся, называется

 

фазовой точкой. Последовательность

 

точек,

которую

проходит

система

 

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

 

зовой траекторией. Для системы из

 

N частиц

фазовое

пространство

 

имеет 6N измерений, а для одно­

 

мерного

осциллятора

фазовое

про­

 

странство

вырождается

в

плос­

 

кость.

 

 

 

распределение

Фазовая траектория гармониче­

Статистическое

для системы указывает, с какой ве­

ского осциллятора в координатах

Vх или р х.

роятностью осуществляются различ­

 

ные микроскопические состояния этой системы, т. е. оно представ­ ляет собой плотность вероятности в фазовом пространстве систе­ мы. Статистическое распределение будем обозначать как р(Г),

где Г — точка в фазовом пространстве.

Зная статистическое

рас­

пределение, мы имеем принципиальную

возможность расчета лю­

бых термодинамических величин. Так, внутренняя энергия

 

£ = | р (Г )# (Г) d r

(III. I)

где функция Гамильтона системы, т. е. ее мгновенная энергия, выраженная как функция фазовых координат

й <г> - 2 ( 4 - + “| ) + Е “«

i 4

1

7

i < k

Здесь pi — импульс t-й

частицы, /я,

и щ— соответственно ее

масса и энергия во внешнем поле и

— энергия взаимодействия

i-й и Л-й частиц, dr =

dpiXdpivdpit... dzN— элемент фазового

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

f = J* р (Г) f (Г) rfr

(III. 2)

где микроскопическая сила

duik KSk dz,

Формулы типа (HI. 1) и (III. 2) будем записывать в краткой форме:

Е = {Н )

f - Ф

где угловые скобки означают статистическое усреднение. Энтропию по аналогии с любой термодинамической величиной

можно записать как результат статистического усреднения:

S = — *<1пр> (III. 3)

где k— постоянная Больцмана.

Возвращаясь снова к истолкованию энтропии как меры не­ упорядоченности, заметим, что упорядоченные состояния менее вероятны, чем неупорядоченные. В этой связи можно привести два шуточных примера: каждый сержант знает, как трудно по­ строить толпу новобранцев; в книге Джерома Джерома «Трое в

одной лодке» очень ярко описано,

как

трудно смотать

веревку

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

что

р ^

1, т. е. что

1п р О,

видим, что для упорядоченных состояний S меньше, чем для не­

упорядоченных.

найти вид

равновесной

функции

На основании (III. 3) можно

распределения, если учесть, что

равновесное

распределение — это

самое вероятное распределение, при котором энтропия макси­ мальна. Другими словами, нужно найти такую функцию, при ко­ торой энтропия принимает максимальное значение. При решении этой задачи нужно учитывать условие нормировки:

| р ( Г М Г = 1 (III. 4)

Кроме того, следует еще учитывать закон сохранения энергии, т. е. мы должны рассматривать только такие распределения, при которых внутренняя энергия будет постоянна. Это. значит, что нужно наложить еще условие на р в виде соотношения (III. 1). Задача отыскания функции, при которой некоторый интеграл от этой функции будет максимален при данных дополнительных условиях, — это задача вариационного исчисления. Для решения