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

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

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

Менее общим* является предположение, что u(t), ф(£), s(t)

функции («внутренних»,

«структурных»,

«микроскопических»...)

«параметров

состояния»

mk = \ik(t),

подчиняющихся эволюцион­

ным уравнениям |д = со(р, П (0),

 

 

 

 

——k.^ —= ak (И- (Q. П (0)

(*=0, 1, .. •),

(10.25")

 

at

 

 

 

 

где со(|д, П)

— выбираемые из различных

физических,

нередко

интуитивных и качественных, соображений функции. Импонирует гипотеза, что скорость роста объемной плотности величины p,(f, х)

(трещин — в твердом теле или воздушных

пузырьков — в воде

и т. п.) является какой-то функцией самой

p(f) и параметров

процесса 11(f) или реакции R (t) (10.5).

 

Соотношения (10.25") нередко дополняются гипотезой о функ­ ции рассеяния и работе

а»*(0=а>к(М-(0, И (0 )М 9 .

-±-b'A=Phd\ih (k=0,

1 ,...,л ),

(10.26)

Р

 

 

 

 

после чего из (10.14) и (10.20)

получается уравнение

 

( ~ r - + T - f -----u>k-Pk)d\ih= 0,

(10.27')

\ д^к

дцк

 

}

 

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

 

 

рк = - г - + т- г -—

“V

(Ю-27")

 

д[1к

d\ik

 

 

Параметры типа р могут относиться к числу параметров Р(т),

Рл в (10.5).

(10.25),

(10.25'),

(10.25") являются

пара­

Все функционалы

метрами термомеханической природы, и соответствующие им оп­ ределяющие уравнения приводимы к виду (10.5).

В случае (10.24)

уравнения

состояния

имеют вид согласно

(10.14)

 

 

 

 

 

а

 

&i =

л

дг|э

(10.27)

йГ ’

~д¥'

 

^

7 ’

 

Но без экспериментального изучения конкретной сплошной среды, ее структуры и свойств, которые и могут быть отражены параметрами типа (10.25'), (10.25"), нельзя указать прямых опы­ тов, которые позволили бы установить число и выражение пара­

метров состояния через характеристики процесса

t), T(t).

В частности, в равновесных опытах с постоянными

по объему и

меняющимися во времени параметрами состояния T(t)t р«(0> б'С?г выражение (10.2) на основании (10.22) преобразуется к виду (при ш* = 0)

т. е. б' Q T имеет интегрирующий множитель Г-1:

^2

(10.28)

Этот интеграл зависит не от того, по какому закону во времени: подается тепло 6<3г и изменяются параметры рл, а только от ко­ нечных и начальных значений параметров. Равенство (10.28) и служит для нахождения энтропии s(T, р). Аналогичное равенство* вытекающее для рассматриваемого объема из (10.12), служит для определения внутренней энергии а(Т, р). Но для нахождения числа и физического смысла параметров р в макроопытах могут быть полезны лишь функционалы (10.6), так как параметры в силу основного постулата должны преобразовывать функциона­ лы процесса (10.25) в функции этих параметров.

Для необратимых процессов основное термодинамическое тож­

дество (10.14)

может быть преобразовано к общему виду,

если

обобщить определение процесса я(т), реакции r(t)

и скалярного

потенциала V*

которые до этого принимались я = П,

r= R ,

I

=

Тождественные преобразования дают

 

 

 

 

 

dw = duTdr\ad§ = dij?,, + r\dTad? =

 

 

 

= dy— Tdr\ + S i = d%+ т]dT + d a i= 0 ,

 

(10.29')

где

 

 

 

 

 

 

 

 

% = и ~ - Т т ] , ф = u — aS,

 

 

(10.29)

 

% = u - T i \ — o § = % - - e i = q > — Tr\.

 

 

 

 

 

Отсюда, выбирая различные допустимые пары процессов я и

соответствующих реакций г, получим всегда

 

 

 

 

 

w‘(it (t)) =

V' (it)—r$) it (0 =

0

 

(10.30)

для произвольного л (т),

При

этом

процессы я,

реакции

г* и скалярный функционал

V1 определяются таблицей 5 функций

(я, г, V).

г можно рассматривать

 

 

 

 

 

Объекты л,

как

(4X4)

квадратные

матрицы или

семимерные

(7-мерные)

ортогональные

векторы

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 5

 

 

я

 

 

 

Г

 

V'

 

 

V (t)= v*

 

1

(Л. 3)

 

 

(Т,о)

 

и(

 

 

 

и

 

2

(л. о)

 

(Т,

—9)

 

¥

 

 

Ф= ц — а^Г

3

(Г.

9 )

 

 

(—л. <?)

 

<

 

 

*£- II

1

.3

4

 

(Т,

а)

 

(—Л> - 9 )

 

%*

 

 

X = tri а?'

5

(и,

1)

 

(Г-1, —Г-%)

 

V

 

 

 

Л

 

 

 

 

 

 

о

 

 

 

 

л0

о

 

 

(10.31)

 

 

 

 

Я = ( Я 0, Я0, JTJL,

>ЛБ) =

 

 

 

 

 

 

О ||я"»"||

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г =

г0

О

 

=

о

 

б)

( т ,/1=1, 2,3)

 

 

 

(r, г0, гъ

 

 

 

 

 

О

||rmn||

 

 

 

 

 

 

 

 

 

со скалярным умножением

 

 

 

 

 

 

 

 

гя =

 

0

0

 

.

 

.

0

0

 

 

(10.32)

 

г я = г 0(я0)-+ г0я0+ ... + гБяБ= г 0(я0)‘+ гтпятГ1

причем

о

о

 

 

 

элементы

(скаляры)

в

столбцах таблицы;

я0 г0— первые

Лтпу ттп — компоненты тензоров,

указанных

в столбцах я, г;

(я0 яь

 

Я5)

и (г0

гь

г5)

— построенные согласно § 9 по

тензорам я и 7 шестимерные ортогональные векторы

(деформаций

и напряжений,

и а, или наоборот — согласно таблице).

 

Например, в третьей строке таблицы:

 

 

 

 

 

 

 

 

л(х) =

(Г, |), г—г1—(— л> о);

 

 

 

 

 

 

 

п0= Т ,

г0= — г); птп= втп, гтп= а тп)

 

(10.33)

 

 

 

 

о о

 

.

ГтпПтп= ОтпЕтп.

 

 

 

 

 

 

Го Ы ' = — Л Т,

 

 

Уравнение

(10.30) имеет вид

 

 

 

 

 

 

 

 

 

 

 

w* =

'Фп (л (т) ) ~

г< (л) я (0 = 0 ;

 

 

Г‘ (я)я = — л' (Я) 7 (0 + <У?п(Я) ётп (t).

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

ад*

и

энтропию s

(10.20) и функцию ф (10.23), положив

 

 

 

 

 

 

г^ф ^ + Гт^ф -Е Ts, Ti\= Ts— ад*

 

 

 

(10.35)

сохранив процесс я = (Г, <§) и введя

реакцию rs= (—s,

а),

полу­

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

(10.23)

в виде

 

ф* (я)—г*(я)я(/) + ад* =

0; я =

(71 £),

гs =

(— s, а).

(10.36)

Все параметры,

входящие

в уравнение

(10.30)

и в таблицу

(я, г, V), могут быть измерены в J^-опытах; предполагается, что

именно по данным этих ^-опытов и построен функционал

У'(я),

вследствие чего (10.30) и называется

тождеством. Здесь

полная

аналогия с законом

Ньютона

/пг = Р (г(т)). «Образцы

сплошной

среды» нередко практически недоступны для .^-опытов, но неко­ торые закономерные проявления ее свойств в некоторых натурных условиях имеются. Тогда можно принять некоторую гипотезу о виде одного скалярного функционала У*(я), например энергии и*(я), и если из тождества (10.30) найти уравнения состояния, т. е. реакцию /*'(я), то решить соответствующую натурному явле­ нию задачу МСС; сравнение результатов с известными законо­ мерными проявлениями позволит оценить правильность гипотезы.

Этот путь научного познания древнее наших наук.

среды,

Определение реакции

г (я), т. е. уравнений состояния

из тождества

(10.30) по

одному

заданному

функционалу

0 (я ) ,

следовательно,

представляется

необходимым,

как и выбор

этого

функционала с учетом свойств довольно большого числа уже из­ вестных в МСС.

Вид тождества (10.36) может быть придан тождеству (10.30) применительно к любой строке таблицы (я, г, V); он отличается выделенным слагаемым ад*^50, имеющим много названий: рассе­ яние, диссипация, некомпенсированное тепло и т. д. Как уже ска­

зано, все параметры, входящие в

таблицу (я,

г,

V) и

(10.30),

в принципе могут быть измерены

в J t-опытах,

но

этого

нельзя

сказать о w* = T(sг\)^0 и об энтропии s: измерим функционал

энтропии Т1'(я), но не st(n) и ад* по отдельности. Эта неопреде­ ленность породила ряд гипотез о существовании вместе с (10.30) еще одного термодинамического уравнения (kd).

Целесообразно следующее доопределение функционалов рассе­ яния ад* и энтропии s, соответствующее сложившимся в МСС представлениям: если в конечном объеме среды возможен ста­

ционарный периодический

процесс

П(Г, # )= П (х ,

/)=П (х, / +

+ 2л/ш) с частотой со то

рассеяние

ад*

и энтропия

5 для произ­

вольного момента t и любого объема

V с границей 2 удовлетво­

ряют условиям доопределения:

 

 

 

*4+2л/(й

d t\p w ’dV = Q1] Q ,=

/4-{-2я/(0

(— qn) d2;

/4*f2Jl/0

dt\pTsdV=Q .

С

J dt

jj

it

V

tx

2

*4

V

 

 

 

 

 

(10.37)

Здесь Qi — количество тепла, выделяемое за один цикл массой объема V через 2 во внешнюю среду.

При заданном функционале ф*(я) не найдено общего решения уравнения (10.36), т. е. реакции iV(n). Но если ф'(я) имеет

дифференциал Фреше, например

 

 

 

 

 

t

 

 

 

 

(я + 6я)— ф (я) = 6 ^ = J ^

J 6я (т) dr,

(10.38)

 

 

 

 

 

и '

т

 

 

и если ф'(л)

имеет вид

 

 

 

 

 

 

 

 

У (л)=

 

 

) Я +

 

 

(10.39)

причем для

любого

процесса

функция фиу^О, то

(10.36)

имеет

частное решение, физически допустимое,

 

 

 

 

г '(я )= (

)

,

w = ~ ^ w > 0;

 

(10.40)

 

 

Dn

)x=t

 

 

 

 

иначе

 

 

 

 

 

 

 

 

 

\

DT i t

Oi]=

\ Den It

 

(10.41)

 

 

 

P

 

 

причем необходимо еще выполнение (10.37) для s.

(10.30)

возни­

На основе произвола л(т)

в тождестве в форме

кает и более строгий подход к определению реакции г'(я).

 

Поскольку тождество (10.30)

справедливо для любого процесса

я(т), оно верно и для я + бл,

и возникает вариационный принцип,

на основе которого можно искать решение

(10.30) при следующих

предположениях и обозначении со = гя.

1.Функции физически допустимых процессов я(т), реакций

г(т) с их производными

я(т), г(т)

и мощностей со(т) =г(т)я(т),

которые определены на

интервале

/о ^ т ^ < о о , принадлежат ли­

нейному функциональному пространству Ф, состоящему из кусоч­

но-непрерывных функций с разрывами в счетном

множестве

М0

фиксированных точек полуинтервала t0^ r < ti,

где

t\> t

на­

ибольшее значение параметра t.

 

со(т), г(т) и

Для двух любых наборов функций из Ф — я(т),

я'(т), г'(т), со'(т) — разности процессов, реакций

и мощностей

бя (т)= я' (т)—л (т), бг (т)= гх (я + бл)—гт (я), бсо (т)=б (г (т) я (т))=г (т) бя (т) 4- бг (т) я (т) + бг (я) бя (т)

принадлежат Ф, причем бл (x)=d/dx • бя(т).

2. Из тождества (10.30) и определения следует, что для любо­

го фиксированного я(т)

существует функционал V1(я)= V (t) е Ф,

а для любых я(т), бя(т)

— приращение 6V'=V* (я+ бя)—У '(я)=

= б у (t) е Ф.

3. По определению функционала г'(я), который при заданном

я(т) равен г(т)еФ , для любых заданных процессов я(т)

и для

произвольных процессов бя (т) справедливо тождество

 

6V" (0—бсо (0 = 0.

(10.42)

Вариационный термодинамический принцип. При заданном функционале I/', определенном на Ф, в фиксированном процессе я(т) и в произвольном ненулевом процессе бя(т) существуют первые вариации, или ненулевые линейные функционалы U(8n), Ь'(бя), такие, что

V* + бя)—V‘ (я)= U (бя (т)) + о N),

г' (яЦ-бя)г' (я)=Ь' (бя) + о ( Ш )

по некоторой норме 8N, определенной на Ф и имеющей вид

sup |kj6n (т) | + sup 162бл (т) |, т е (/0

t) \M 0,

где k\>0, ki>0 — конечные величины.

реакции г'(я) удов­

При заданном У'(я) и 6iV-»-0 функционал

летворяет функциональному тождеству относительно бя, -которое

соответствует

(10.42)

 

 

6V' (я (т))—г* (я (т)) бл (0—бг' (л (т))я (0 = 0.

(10.43)

Здесь 61^', бг' суть L', L', т. е. первые вариации

V‘ и г'.

 

Теорема .

Если множество функций я(т)

совпадает с прост­

ранством непрерывно дифференцируемых функций, то

 

t

t

 

 

6V' = I

W dxv (*, x) + a (0 бя (/„) = \ 8nu (x) dxv‘i (t, т)+ а‘1(t) 6п4/(*0),

to

Г.

 

t

 

бг' = J бя (x) dxR (t, x) + бл (t0) В (0,

 

to

t

бг'я (0 = 1 бя(/ям (0 dxR‘M (t, T) + бя;/ (/0) nhl (0 Bitkl (0, (10.44) to

где индексы i, j, k, I принимают значения 0, 1, 2, 3, вектор v(t, х) (известен при данных V\ л) и тензор-функция R(t, х) с ограни­ ченной вариацией по х определены с точностью до аддитивной постоянной (зависящей только от t).

Рассмотрев множество процессов 6я(т), проходящих через точ­ ку t0>получим из (10.43) функциональное уравнение для реакции

г '(я (т )):

v (t, т)—я (t)R(t, т)—гh(t, т)= с(0 ,

(10.45)

где c(t) — постоянный (по т) вектор

С(t)=v(t, t)— n(t)R(t, t)— rh[t, t)=v(t, t0)~n(t)R(t, t0)— rh(t, i0)

и rh(t, T ) — вектор скачков функционала г*(я(т)), который в простейшем случае равен

Г*(*, x)=T(t)h(t, т),

причем h(t, т) = 1 при x = t и h(t,

т) =0 при т <t.

 

При еще более общих предположениях установлены условия,

при которых 1/'(я) = \ i ‘ (n)n(t)

и тождество (10.30)

в виде усло­

вия ортогональности

 

 

(И (я)—г'(я))я(/) = 0

(10.46)

имеет единственным решением г'(я) = V /(я).

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

параметрами рг

(г= 1 2,

...), включающие и структурные типа

mk в уравнения

вносятся

поправки. Предполагается, что возни­

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

и что возникает дополнительный к б'QT приток тепла pq$tAV за счет проникающего облучения или источника диссипативного ха­ рактера. Например, в случае электромагнитного поля с электри­ ческим вектором Е и магнитным Н (§ 22)

6'Лр=№рЛ = (1/4я) (EdD+ HdB),

(10.47)

где D В — векторы электрической и магнитной индукции, т. е. набор параметров рг состоит из компонент двух векторов: р= = (D, В); объемная плотность источника тепла рq$ представляет джоулево тепло.

Первый и второй законы термодинамики получаются в обоб­ щенном виде:

 

р

01

01 + W» ~ div Ч + Р9р.

 

(10.48)

 

 

 

 

 

 

 

 

Рт ■$“ =

—div ^+ Р^р+ рш* ’

 

 

 

 

01

 

 

 

 

причем входящие

в них функции

или функционалы

(и,

s, w*,

qt,

) определяются процессом

изменения функций

времени

&{х),

Т(т), р(т).

 

 

 

 

 

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

ванной движущейся массы

с объемом Vg и поверхностью

2 g по­

лучается из (10.48) интегрированием по объему и из определения внутренней энергии Ug и энтропии Sg рассматриваемой массы:

Ut (t)= f pudV, Sg

рsdV

(10.49)

Из первого равенства (10.48) и (8.29)

получаем

 

dUg= J (Sild£tf—div qdt + Wg/dt + pq^dt) dV = ve

— b' A— dKg— J qnd2dM-6'lp,

или согласно (10.2')

d (Ke+ Ug) = 8 ’A + 8'Qr + б'Лр,

(10.50)

U A ^ ^ W ^ + pq^dVdt.

Здесь Kg — кинетическая энергия, б'A — работа массовых и по­

верхностных сил, б'QT — приток тепла, б'Ар — приток

энергии,

связанный с параметрами р.

с множителем

Т~х и

Интегрируя второе уравнение (10.48)

учитывая преобразование

 

 

§ j r d i v q d V = \ - L q n d 2 + Г

q grad TdV,

(10.51)

Vg

получим интегральную форму закона баланса энтропии

d S g = - ^ - L q n d m + ^ - L (w*— -yqgradT + Pq^dVdt. (10.52)

 

 

vg

Согласно

закону

Фурье —q gradT=A(gradT)2> 0 . Сумму

w*— q grad Т/Т ^ 0

называют иногда объемной плотностью не­

компенсированного тепла в термомеханических процессах. След­ ствие (10.52) в виде

s g + £ j r

qnd S - j jr M td V =

£ ( —qf - - + ^ ) dV> 0 (10.53).

*

vg

ve

называют иногда неравенством Клаузиуса—Дюгема. Оно неопре­ деленно, если не определена энтропия S (10.49) или w*.

Отметим еще одно постулируемое фундаментальное неравенство для замкнутых по деформациям изотермических процессов [13]: работа внутренних напряжений в таких циклах неотрицательна

^ SiJ dei/ = ф

0,

Т = const.

(10.54)-

§ 11. ЗАМКНУТЫЕ СИСТЕМЫ УРАВНЕНИЙ МСС

 

Для малой окрестности

физической

точки (частицы)

среды

установлены дифференциальные и интегральные уравнения сохра­ нения массы, импульса (уравнения движения), сохранения энер­ гии, баланса энтропии (уравнение притока тепла), а также урав­ нения, связывающие тензор напряжения и вектор теплового по­ тока с деформациями, температурой и немеханическими задан­ ными параметрами. Эти соотношения в принципе определяются,, и притом однозначно, непосредственно в «^-опытах для всех воз­ можных в частице процессов; поскольку все входящие в эту сис­ тему равенств параметры измеряются приборами и системе удов­ летворяют, группа параметров, названная реакцией (г), одно­ значно определяется группой процесса (я). Следовательно, для малой частицы решение существует: r(t) = г /(я(т)). Поэтому пе­ речисленная система уравнений в МСС называется замкнутой для всех внутренних точек области движения среды.

В частице состояние в любой

момент

однородно

(одинаково

во всех точках частицы), поэтому для нее

все

параметры внутри

и на границе частицы одинаковы

(деформации,

напряжения, тем­

пература, градиент температуры

и т. д.). В конечных

областях

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

В силу основного постулата МСС решение замкнутой системы уравнений МСС для внутренних точек существует при некоторых

начальных условиях (/ = /0)

и условиях на границе области дви­

жения или рассматриваемой области.

 

Адекватность методов

Лагранжа и Эйлера в МСС позволяет

пользоваться любым из них. Сначала будем следовать

методу

Лагранжа и систему координат х/ в начальный момент

времени

считать декартовой ортогональной. Тензоры напряжений

и дефор­

маций по-прежнему обозначим S, <£. Общие уравнения МСС для внутренних точек области движения следующие. ^

1.Уравнения, выражающие тензор деформаций и плотность

•среды р через вектор перемещения u(x, t) или х=<р=х+и (§ 4):

 

Ъ

gii= A U h Aki= - 1£-; рЛ =р0;

(11.1)

 

 

 

 

 

I= d efu ; ег/= ( £ ;/-6,/)/2;

Y g = A = \A * \,

 

Эти уравнения алгебраически выражают е»/ и р через первые про­

изводные

дх‘/дх! = A\j= 6/ + ди'/дх'.

 

 

2.

Уравнения движения, которые эквивалентны одному вектор­

ному уравнению

(§ 8):

 

 

 

 

V(S‘ + p (F - * ) = 0,

- ^ - = х ,

(11.2')

причем F= F(x, t) — заданная массовая сила. Уравнения записы­ ваются также в виде

 

po( ? - X O =

^ ( A 4 'S '" 0 .

 

 

(Н-2)

Структура этих трех скалярных уравнений, очевидно, такова: они

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

S ij

и их

первые

производные S% со вторыми производными

х по времени, а так­

же первыми и вторыми производными х по координатам.

 

3.

В тех случаях, когда за

независимые искомые функции при­

нимаются деформации гц или напряжения

Sl/,

должны

быть

удовлетворены еще условия совместности (§ 5)

 

 

 

 

Rlmpq—О,

 

 

 

где индексы Itnpq имеют шесть независимых комбинаций

(1221 и

т. д.). Тензор кривизны RimPq алгебраически

выражается

через

е//, е//, k и 6//, kmy причем через вторые производные от е;/ — линей­

но. Эти шесть уравнений кратко запишем в виде

 

 

 

R = R (ei/t

eijtkin)= 0.

 

 

(11.3)

Существенно, что если, как и сделаем, за независимую искомую функцию принять вектор х или перемещение и, то на основании (11.1) уравнения (11.3) обратятся в тождества и потому могут не рассматриваться.

Следовательно, для вектора и и тензора напряжений S полу­ чается только одно уравнение (11.2'), т. е. для девяти функций

х1, ха, х3 S11 S22 S33 S12= S 21 S23= S 32 S31 = S13

мы из (11.2) имеем только систему трех дифференциальных урав-