Механика сплошной среды
..pdfвать как равенство нулю компоненты некоторого многомерного вектора — оператора А над двумя функциями 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 |
получаем |
|
т. е. возвращаемся к задаче о баротропной жидкости.