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

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

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

вать как равенство нулю компоненты некоторого многомерного вектора — оператора А над двумя функциями x=<p(x, t), Т (х, t) по параметрам х, t в ограниченной области их изменения

X E [G(/), 2(0], * o < * < * '< * r (12.40)

Для обратимых операторов компонентами А! = 0, А2 = 0 векто­

ра А= 0, относящимися к области параметров

(12.40), исключая

2, являются одно (векторное) уравнение

движения и одно урав­

нение теплопроводности. Компонентами А3

А4

из области (12.40),

исключая G(i), являются одно граничное условие (векторное) и

одно условие притока тепла. Компонентами А5

А6 А7 из области

изменения параметров (12.40) при

/ = /0 = const

являются началь­

ные условия для ф и Т: при t = to

 

 

 

А5 = <f—ф0 (х)=0,

 

А6 =

^Т— Ч>о (х)= 0,

А7= ТТ0= 0 .

 

01

 

 

Если обозначим

!!(%, ть

Ле» Л7) произвольный не­

ортогональный А вектор с компонентами, определенными в соот­

ветствующих областях изменения параметров х,

t, то все семь

(15 скалярных) уравнений запишутся в виде

 

Л(х, t) А (ф, Т) = 0.

(12.41)

Наиболее просты линейные задачи, для которых оператор А — линейный относительно (ф, Г), так как для таких задач наиболее развиты точные методы решения. В общем случае распространен метод линеаризации на основе уравнения (12.41) в вариациях: если некоторое решение ф°(х, t), Г°(х, t) уравнения (12.41) изве­ стно, то разыскивается близкое решение ф1 7м удовлетворяющее уравнению

лЧФо, 7°; бф, 67) = 0,

(12.42)

причем

Ь = бА=А(ф1, Г1)—А(ф°, 7°) = А (ф° + бф, Г° + бГ) —А(ф°, 7°).

Для уравнений МСС с обратимым SF оператор L будет линейным относительно бф, бТ Эту процедуру иногда удается применить многократно и построить последовательность (ф1 Г1) ,..., (фл, Тп)> сходящуюся к решению уравнения (12.41). Формальное доказа­ тельство теоремы существования для уравнения (12.41) во мно­ гих случаях является сложной задачей. Вопрос о единственности решения в определенной мере связан с решением уравнения в ва­ риациях (12.42). Если заданные функции, входящие в уравнение

г|А= 0 (ф0, фо х£ ...), не варьируются, то уравнение в вариациях лЦфо, Т°, бф, 67") = 0 по определению ф°, Т° имеет нулевое реще-

ние: 6ф= 0, 8Т = 0. Но если существует еще и решение, не равное тождественно нулю, то заключают о возможной неединственности и неустойчивости решения <р° Т° и неустойчивости самого движе­ ния среды. Это грубо обоснованное заключение в большинстве случаев оправдывается в опыте. Строгая постановка исследования устойчивости решения уравнений (12.41) принадлежит к числу современных вопросов, не получивших общего решения.

Каждая компонента вектора А вместе с соответствующей ком­ понентой т] имеет свою определенную область изменения парамет­

ров х,

ty например, А1 Л2 т|Ь т\2 имеют область X G G, t ^ ( t 0y t\)

и т. д. Если интегралом т]А по области его определения назвать

сумму

интегралов от произведений его компонент Ак на соответ­

ствующие компоненты щ у из которых каждый

взят по «объему»

соответствующей области определения Qky то

уравнение (12.41)

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

 

$

T i A d

Q =

0 .

и

 

 

Отсюда следует (12.41) не

только

при

совершенно произвольном

т|, но и TJ в виде некоторого оператора по х, t над произвольными вариациями искомых функций. Это приводит к вариационной по­ становке задач МСС.

Глава IV

КЛАССИЧЕСКИЕ ТЕОРИИ АЭРОГИДРОДИНАМИКА И ТЕОРИЯ УПРУГОСТИ

Простейшие в МСС тела — идеальные и классические (нью­ тоновские) жидкости, идеально упругие твердые тела. Эти тела (идеализированные модели реальных тел) обладают фундамен­ тальными свойствами реальных жидких и твердых тел, причем свойства, во многих случаях второстепенные, не учитываются. Опыт показывает, что поведение многих реальных жидкостей, газов и твердых тел в определенных условиях достаточно точно описывает­ ся уравнениями механики сплошной среды, построенными для ука­ занных идеальных тел. Методическое значение моделей состоит еще и в том, что из сопоставления с опытом получается возмож­ ность изучения отклонений свойств реальных тел от свойств мо­ делей и, значит, возможность уточнения теории.

Жидкостями в механике сплошной среды называются тела, со­ противление которых сдвигу при любой деформации стремится к нулю, если скорости деформации равны нулю в течение достаточ­ но большого промежутка времени (t->оо).

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

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

Напряженное состояние малой частицы любой среды в любой момент t характеризуется тензором напряжений 3, который в главных осях напряжений всегда имеет диагональную матрицу

ог О О О о2 0 . О 0 о3

Экстремальные касательные напряжения, действующие в этой ча­ стице в момент t по плоскостям, делящим пополам углы между плоскостями главных напряжений, равны

^ _ а1

~ _ а2 а3 _ а3 а1

112

Г

» 123

~

> 131

По определению рассматриваемая частица будет частицей жид­ кости, если при <§=const и t->оо Т12=Т2з='Тз1= 0, т. ё.

*->■ оо, а1= а 2==а3= —

или частицей твердого тела, если при <§=const

t 00» I I max ^ тоо > О

т. е., вообще говоря, при t-+oo а^ог^огз^аь

Классические теории относятся к термомеханическим процес­ сам.

§ 13. ИДЕАЛЬНЫЕ ЖИДКОСТИ И ГАЗЫ

Идеальная жидкость (газ) — это среда, в которой рассеяние отсутствует (йУ*=0) и сдвиговые сопротивления которой при лю­ бой деформации и скорости деформации в любой момент време­ ни равны нулю, т. е. для любого t

°1= а2 = а3= --Ру

(13.1)

следовательно, и девиатор напряжений равен нулю. Поскольку главные оси напряжений всегда взаимно ортогональны, то на ко­ сой площадке с нормалью v(/b /2 /з) вектор напряжений ov= P v равен

av= а г1г+ а212+ <т3/3== —р\ ,

(13.2)

т. е. он направлен по нормали к площадке и имеет величину —р. Следовательно, давление р по всем площадкам, в данный момент проходящим через данную точку среды, одинаково. Обычно пред­ полагается, что р > 0, т. е. это всегда действительно давление. Од­ нако реальные жидкости могут выдерживать и некоторые всесто­ ронние растягивающие напряжения (р< 0), и в понятие идеаль­

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

В

эйлеровых ортогональных

декартовых

координатах Xi

(i=

= 1

2, 3) тензор напряжений в

идеальной

жидкости имеет

про­

стейшее выражение

 

 

 

 

 

 

 

ос1=

—р8ц.

(13.3)

Уравнения движения сплошной среды в рассматриваемых ко­

ординатах уже найдены в § 8 в виде

(8.13):

 

 

 

Р — — рх; =

др_

(i= 1

2, 3),

(13.4)

 

ш

dxt

 

 

 

или любой из двух векторных форм, включающих (8.16),

Р—---- pF= —grad р,

at

 

 

 

(13.5)

-^ - + — gradv2—vxrotv + — gradp=F.

 

dt

2

p

 

Они называются динамическими уравнениями Эйлера.

К ним

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

 

 

 

— + р div v = 0 .

(13.6)

 

 

dt

 

Система (13.5), (13.6) есть система совместных дифференциаль­ ных уравнений в частных производных первого порядка относи­ тельно вектора скорости v и скаляров р (давления) и р (плотно­ сти). Это незамкнутая система, так как для пяти функций коор­ динат и времени vu v2 v3t р, р она дает только четыре уравнения указанного типа.

Динамические уравнения Эйлера можно записать также в лагранжевых координатах. Пусть, как и ранее, х= х(х, t) — закон движения частицы, причем х=х*е/ — ее начальный декартов ра­ диус-вектор, а

X = ф (х, /) = x-fu(x, t)

— ее текущий радиус-вектор. Поскольку в момент t лагранжева система координат x;=const (i= l, 2, 3) является, вообще гово­ ря, криволинейной, ковариантные и контравариантные базисы и метрические тензоры определяются законом движения (4.15), <4.23), (4.48)

 

 

(13.17)

т. е. а* э\ gn, gu выражаются

через частные производные

*1-.

За искомые функции можно

принять либо xi(xu х2 х3

dxk

/)> либо

#г(хь х2 х3 /), причем

gij и g'i алгебраически выражаются че­

рез компоненты тензора деформации.

Найдем выражения

контравариантных компонент S i}' тензора

напряжений 5 в момент t в базисе э -. По определению идеальной жидкости вектор истинного напряжения на площадке, построен­ ной на векторах э2 и э3 направлен по нормали к ней (а направ­ ление нормали совпадает с э1) и равен давлению р, т. е.

Аналогичные выражения получим для Р(2) и Р(3). Вектор напря­ жения S1 равен

Si = p i-^ T i= —рэь

аналогично выражаются S2 S3. Следовательно,

& = &кэ к= — р9*.

(13.8)

Умножая это равенство на э* и учитывая (13.7), получим

(13.9)

Соотношения (13.9) представляют просто преобразования со­ отношений (13.3) от базиса е, к базису э1. Уравнения движения сплошной среды в лагранжевых координатах, имеющие вид (8.8) или, после умножения его на э2

V/S/7 + р (F*w‘)= О

преобразуем, учитывая (13.7), (13.9) и очевидные соотношения

 

V/ W

!)= gaviP= ga OXj ,

 

 

 

F‘= F di

 

 

OXj

,

(13.10)

 

 

 

 

 

 

f

:

: : iirr dXh

:i

d 2Xh

dXh

 

w

= w3l= g1]Wb — — = P1]-----------—

 

 

6

h dxj

s

a/2

dxj

 

где Х*:, Wk — декартовы компоненты (в репере е,) векторов мас­ совой силы и ускорения. Имеем

Рg“'

d2xk

\

- ^ - = 0 .

 

дР

)

дх\

Отсюда, умножая на g/m получаем лагранжеву форму уравнений движения идеальной жидкости

-s r( ^ -x‘)+T^ =0

<1311>

Здесь, конечно,

д*хк _

3 4

дхк

я

-1

dt2

а/2 ’

dxt

— °ih

1

дик

дхс

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

Ро _=1/g или 9о__ dXi = А .

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

ординат х/ и времени t: х и х2 х3

(или щ, и2, и3), р, р. Эти урав­

нения имеют второй порядок по t

(относительно х) и первый — по

Хг (относительно X р).

 

Идеальная несжимаемая жидкость — это идеальная жидкость, плотность каждой малой частицы которой во времени не изменяет­

ся. Если

при t= t0 плотность ро была постоянной, то она в несжи­

маемой

жидкости

останется постоянной

и при t> t0;

такая жид­

кость называется

однородной. Если при

t= t0 ро=ро(хь

х2, х3), то

она такой останется и при t> t0i т. е. р = р о (х ь х2, х3) *

Уравнения

движения в форме Эйлера и условие несжимаемости в эйлеровом пространстве имеют вид (13.5) и

div v=0,

(13.13)

причем представляют собой замкнутую систему четырех дифферен­ циальных уравнений первого порядка в частных производных для четырех функций координат и времени: vu v2j и3 р.

В форме Лагранжа уравнения движения идеальной несжимае­ мой жидкости имеют вид (13.11) с условием несжимаемости

[ dxj = А = \.

(13.14)

£ = 1 Ы = 1 или I дх,-

 

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

причем

 

 

р=р( Р),

(13.15)

 

 

 

 

 

 

 

C * = -- ^ > 0 ,

Л . > 0 .

(13.16)

 

=V—

dp

dp

 

Величина с

называется

скоростью звука.

 

V

dp

 

 

 

 

встреча­

В начале курса

на примере идеального газа мы уже

лись с уравнением состояния p=RpT, которое для изотермических

процессов совпадает с (13.15), причем р(р)

есть однородная ли­

нейная функция, а для адиабатических — приводится к

виду

(13.15), причем р(р) — степенная функция

с показателем

у> 1»

так что в обоих случаях

 

 

(£)’ у > 1 .

(13.17)

* В этом случае жидкость называется н е о д н о р о д н о й .

Уравнения движения и сохранения массы (13.5), (13.6) (в эй­ леровых координатах) или (13.11), (13.12) (в лагранжевых) за­ мыкаются для баротропной жидкости пятым соотношением (13.15).

Введением функции давления при потенциальной силе F (§ 8)

P ( p ) = J j £ - = J i l d p , F = - g ra d < rp

(13.18)

из (13.4), (13.5) получаем замкнутую систему для V, р

— F+ grad Р (р) = 0,

+ р div v= О

 

(13.19)

+ grad ^P + ^ p- f j

= vxrot v.

Аналогичная система в лагранжевых координатах получается из уравнений (13.11), (13.12), (13.15).

Для слабо сжимаемых жидкостей при небольших давлениях скорость звука с в ряде случаев может считаться постоянной. При этом Р=с2In р; уравнения сохранения массы и динамические урав­ нения будут содержать только функции v и Р, причем первое из них имеет вид

1

dP .

л

---------- с2

h div v = 0 .

dt

 

Отбрасывая в уравнениях (13.19) квадратичные члены v2 vX

ч . .

дР

получим систему

XrotvH

 

-~ - + gradP—F = 0 ,

01

div v + — -^ - = 0,

c2 dt

приводящуюся к волновому уравнению с источником

— - ^ - = ДР—divF,

(13.19')

с2 dt2

'

из которого и следует, что c=~]/dp/dp есть скорость распростране­ ния слабых возмущений (скорость звука).

Идеальный разреженный или совершенный газ — это идеаль­ ная жидкость, подчиняющаяся уравнению состояния Клапейрона

причем внутренняя энергия единицы массы прямо пропорциональ­ на температуре

и = с0Т = —

(13.21)

0R

Вначале курса (§ 2) доказано, что такими свойствами обладает одноатомный газ при давлениях, не превышающих сотен атмосфер,

причем в системе C G S постоянная c,;=3£/2m (k — постоянная Больцмана, m — масса атома), R=k/m.

Многие инертные газы с достаточным приближением подчи­ няются уравнениям (13.20), (13.21). Например, для воздуха R= =2,87-104 см2/с2*К.

Уравнение (13.20) представляет собой соотношение между на­ пряжениями, деформациями и температурой для рассматриваемо­ го тела. Выражение внутренней энергии и (13.21) через два един­

ственных

для

этой среды независимых параметра

состояния

р и

Т (w от

р не

зависит) получено в статистической

механике

и из

опыта.

Запишем закон сохранения энергии (§ 10)

р 8u=8'Q + &]8еи-.

Так как

y i A ' i = t;ii = diV V,

то

из (13.3)

 

 

 

 

Sl'i8Eif= o i/vif8t = — р div v dt.

 

Из

(13.6)

1

dp

 

 

,.

 

 

div v = ----------— .

 

 

 

p

dt

 

Учитывая все это, закон сохранения энергии запишем в виде

 

p8a=6'Q + р

или cvp8T=8'Q-\-p-^~.

(13.22)

 

Р

 

Р

 

Отсюда следует, что cv — коэффициент теплоемкости при посто­ янном объеме (p=const). Внося в (13.22) следствие (13.20)

б Inр= 8 Inр—61п7\

получим

p(cv + R)8T=8'Q + 8p,

(13.23)

т. е. cv + R есть коэффициент теплоемкости ср при постоянном дав­ лении (p=const). Таким образом, получилась формула Майера

cp- c v=R.

(13.24)

Из первого и второго законов термодинамики следует

р(6ц—T6s)= p-^~,

Р

откуда с учетом (13.21) имеем

6 s = i - 6 L _ tf J L ,

Тр

т. е. определяется энтропия

s= cvln -^- + const,

(13.25)

Pv

 

где у=cp/cv — число Пуассона (для воздуха у=1, 4).

К неизвестным функциям v, р, р, входящим в незамкнутую си­ стему уравнений Эйлера (13.5) и (13.6), добавилась еще одна функция — температура Т или энтропия s, связанные с р, р со­ отношениями (13.20) и (13.25). Но теперь закон сохранения энер­ гии дает еще одно дифференциальное уравнение

р -^ - + divq—р dln£- = 0.

 

(13.26)

dt

 

dt

 

 

 

Поток тепла q для большинства изотропных сред

связан

с по­

лем температуры Т законом Фурье

 

 

 

 

q = —Xgrad7\

 

(13.27)

где X — коэффициент теплопроводности, вообще говоря, извест­

ным образом зависящий

от Т

Внося значения и

(13.21)

и q

(13.27) в (13.26) и используя

уравнение

сохранения

массы

(13.6), получим уравнение теплопроводности

 

 

 

рcv — —

div (X grad Т) + р div v=

0.

(13.28)

Теперь система уравнений (13.5), (13.6), (13.20) и (13.28) для век­ тора v и скаляров р, р, Т стала замкнутой для рассматриваемого идеального газа. Вместо Т за искомую функцию можно принять и или s, или другую, выражающуюся через р, р, Т термодинамиче­

скую функцию (§ 10), например энтальпию (теплосодержание) i=cPT

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

батическим),

то в

уравнении баланса энтропии (§ 10),

кроме

ш*= 0, надо

положить и 6'Q=0, и значит, энтропия частицы бу­

дет постоянной во

времени (но

может быть различной у

разных

частиц). Из

(13.25)

при s=const

получаем

 

т. е. возвращаемся к задаче о баротропной жидкости.