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

книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах

.pdf
Скачиваний:
7
Добавлен:
25.10.2023
Размер:
12.62 Mб
Скачать

может быть представлен как произведение скаляра X на единичный тензор

 

0

0 1

' 1

0

0

0

X

0 • = \ -

0

1

0

0

0

1

.0

0

1.

т. е. он эквивалентен

скаляру.

 

 

 

Заметим, во избежание недоразумений, что тензор не определя­ ется однозначно тремя своими инвариантами. Чтобы вполне опре­ делить, например, симметричный тензор, нужно, кроме трех инва­ риантов, задать еще три главных направления, а они также зада­ ются тремя числами: двумя направляющими косинусами одного из направлений, например е( 1 ) , и одним направляющим косинусом од­ ного из двух других главных направлений, т. е. мы приходим к тому, что симметричный тензор определяется шестеркой чисел, как и в случае задания его в форме (3.7), а меняется только смысл этих чисел.

4. Криволинейные координаты

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

 

Яі =

ЯЛх,

У, z),

qi=q2{x,

у, z),

q3=q3(x,

у, z)

(4.1)

таких, что обратные

функции

 

 

 

 

 

x=x(qu

q2,

<7з),

y=y(qu

q2, q3),

z=z(qu

q2, q3)

(4.2)

также

однозначны

 

и дифференцируемы, то между

тройками

чисел

х, у, z

и qi, qz, qz

устанавливается

взаимно

однозначное соответст­

вие. Поэтому положение точки в пространстве можно определить заданием трех чисел qi, q%, qz, которые называются ее криволиней­ ными координатами. Очевидно, выбор функций qi, q2, qz не может быть вполне произвольным: требование взаимной однозначности со­ ответствия между криволинейными и прямолинейными прямоуголь­ ными координатами означает, что якобиан

dqi

<>Я\

dqt

 

дх

ду

dz

 

dq2

 

dqo

(4.3)

дх

ду

dz

 

ддз

<Э<7з

дЯз

 

дх

ду

dz

 

должен отличаться от нуля.

 

 

 

Три поверхности

 

 

 

qx (х, у, £) = const, q2(x,

у, z)=const, q3(x, у, z)=const

 

называются координатными поверхностями. Линия пересечения двух координатных поверхностей есть координатная линия. Коор­ динатная линия <7i образуется пересечением координатных поверх­ ностей 92 = const и <7з = const; вдоль этой линии меняется только координата qu а координаты q2 и q3 сохраняют постоянные зна­ чения.

Введем единичные векторы Єї, е 2 , ез, направленные по касатель­ ным к координатным линиям в точке М (рис. 4). Вектор производ­ ной dr/dqi в точке М имеет, очевидно, направление Єї, ибо qi и q3 при дифференцировании сохраняют постоянные значения. По­ этому

где Li — длина вектора dr/dqi, и вообще

-*L.=Lfil(i=h

2, 3).

(4.4)

 

А так

как

 

 

 

 

дг

дх

ду

dz

 

 

 

 

 

dqt

 

 

TO

дх

ду

 

 

 

 

 

 

 

 

 

dqt

 

 

Система криволинейных ко­

 

 

 

 

(4.5)

ординат.

 

 

 

 

 

Величины

L l t Li, L 3

называются

коэффициентами Ламэ.

В данном

параграфе

мы

будем

рассматривать только ортого­

нальные криволинейные координаты, т. е. такие, что в каждой точке

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

e i , е 2 , е 3

взаимно

перпендикулярны или

иначе

 

 

 

(е„

еА ,)=0,

іфк.

(4.6)

Необходимое и достаточное условие ортогональности криволи­ нейных координат есть, очевидно, условие ортогональности векто­ ров dr/dqi и dr/dqk при 1фк, т. е.

дх

дх

, ду

ду , dz

dz

Q

(4.7)

дЯі

dqu ~

dqt

dqk

dqb

~~

 

Однако в решении конкретных задач обычно не приходится проверять ортогональность так или иначе выбранной системы коор­ динат, а наоборот, выбирать систему координат, адекватную задаче из условия ортогональности координатных линий. В таких случаях условия (4.7) выполняются автоматически и проверять их не при­ ходится.

Вектор F в системе ортогональных криволинейных координат можно представить разложением

F^Ffr+F^+Fte

(4.8)

Числа Fx, F2, F3 называются криволинейными составляющими вектора или проекциями вектора F на оси криволинейных коорди­ нат. Беря в качестве вектора F вектор а'г, получаем из (4.8) и (4.4)

dr=-~

dqx + ^ d q 2 + -^-dq3=Lx

dqxQx - f Z.2 dq2e2 + L 3

dq3e3.

Таким образом, составляющие вектора dr есть

(4.9)

 

 

ds—L.dq,

{І=Л,

2, 3).

(4.10)

В конкретных задачах определение dsi по dqx обычно легко вы­ полняется по свойствам выбранной системы координат, поэтому коэффициенты Ламэ большей частью определяются из формулы

(4.10), а не из формулы (4.5).

 

 

 

 

 

 

 

 

Пусть

N — точка, бесконечно

близкая

к М,

так

что

dr=MN

(рис. 5). Проведем

через точки М и N по три координатных поверх­

 

 

 

 

ности. Эти

поверхности образуют

 

 

 

УУ бесконечно

малый

криволинейный

 

 

 

1

параллелепипед

с

ребрами

ds^

 

 

 

 

Площади doi, do*, daa

его

граней,

 

 

 

 

нормальных

к

направлениям

е ь

 

 

 

 

ег, е3 , соответственно

будут

 

 

 

 

 

 

do1=L2L3dq2dq3,

 

 

 

do2—L3Lxdq3dqx,

 

 

 

 

 

da3=LxL2dqxdq2.

 

(4.11)

Рис 5.

К

определению элемента

а

объем dV определится по формуле

объема

н

элемента

поверхности

 

dV—LxL2L3dqxdq2dq3.

 

(4.12)

в криволинейных координатах.

 

 

Выразим теперь основные дифференциальные операции в кри­ волинейных координатах.

По определению grad Ф, проекция этого вектора на ось qi есть дФ/dsi или, согласно (4.10),

 

 

дФ

1

дФ

 

 

 

 

 

dsi

Li

dqi

 

 

 

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

 

 

 

 

 

 

 

 

ь

Ц

dqx

' L 2

dq2 1 L 3

dq3

4

'

Отсюда видно,

что для

оператора

Гамильтона

мы

имеем

вы­

ражение

 

 

 

 

 

 

 

 

V

Ц

dqx

' L 2

dq2

Т I 3 dq3

'

 

 

 

 

 

 

 

 

 

 

Для выражения других операций поля в криволинейных коор­ динатах существуют три пути:

1) можно пользоваться формулами (2.16), (2.24) и т. д., вы­ ражающими эти операции через оператор Гамильтона; при этом

необходимо учитывать, что единичные векторы криволинейных ко­ ординат е ь Є2, ез сами есть функции координат; последнее делает этот путь очень громоздким;

2)можно исходить из определений операций поля, т. е. из вы­ ражений (2.16), (2.24) и т. д., рассматривая их в криволинейных координатах так же, как это было сделано в § 2 и 3 для прямоли­ нейных координат;

3)можно апеллировать к общей теории тензоров.

Любой из этих путей позволяет получить следующие выражения для основных дифференциальных операций поля в криволинейных координатах.

1. Криволинейные составляющие wt ускорения w точки, движу­ щейся со скоростью v обобщение формул (2.13) для случая F = v:

 

З

 

 

З

о

 

 

З

 

n , , - * L i . V A

 

dvi

V

vk

dLk

і V

v(Ok dLt

2. Расхождение

вектора

F — обобщение

формулы

(2.16)

d i v F

1

Г d{L2L3Fx)

 

. d{L3LxF2)

, d(LxL2F3)

)

 

LXL2L3

\

dqx

 

 

dq2

"r"

dq3

f'

3. Криволинейные составляющие вихря вектора F — обобщение формулы (2.22)

 

1

f

d(F3L3)

(rotF)1 = =

 

{

dq2

 

1

f

d(FxLx)

(rot F ) 2 =

L3LX

[

dq3

(rot F);3"

1

f

d(F2L2)

Lit*

 

dqi

д {F2L2) дЧз

д (F3L3) dq\

d{FxLx) dq3

4. Оператор

Лапласа

от

 

скаляра

Ф — обобщение

формулы

(2.26)

 

 

 

 

 

 

 

 

 

у , 2 ф г ^

1

Г Д / L2L3

 

_ дФ \ •

д ( L3LX

_ дФ \ .

 

LXL2L3

\ dqi {

L x

 

dqi j

" 1

" dq2 [ L 2

dq2

J '

 

 

і

д

I LX 2

дФ

 

 

 

 

'

dq3

\

L 3

dq3

 

 

5. Криволинейные составляющие вектора div T — расхождения тензора Т с криволинейными компонентами Fih — обобщение фор­ мулы (3.21)

 

 

3

 

 

 

 

 

 

(Л\\т Т\

1

V I

1

д

I L\LiLbLi

77 \

Fkk

dLk \

3 Заказ № 428

Глава

II

О С Н О В Н Ы Е

УР А В Н Е Н И Я Г И Д Р О М Е Х А Н И К И

5.Основные уравнения механики сплошной среды

Будем рассматривать некоторую (не обязательно жидкую) дви­ жущуюся сплошную среду, пронизывающую в своем движении про­ странство, в котором задана какая-либо неподвижная система ко­

ординат. Эту систему для простоты можно считать

прямолинейной

прямоугольной, не нарушая этим общности

построений. Выделим

в пространстве тетраэдр (рис.6)

МАВС,

ребра

MA, MB,

МС

которого

параллельны

осям

координат х, у, z соответст­ венно. Если 5 — площадь грани

ABC, а

п — внешняя

нормаль

к этой

плоскости, то

площа­

ди граней МВС, MAC и МАВ

будут равны соответственно

5 cos

X •

Рис. 6. К выводу динамического урав­ нения механики сплошной среды.

Силы, действующие на грани ABC, МВС, MAC и МАВ тетра­ эдра со стороны остальной среды могут быть выражены так:

phSr p*S cos

py S cos

p^Scos I

где pn, p*, py, pz — средние значения сил, приходящихся на единицу площади соответствующих граней. Эти средние значения есть не­ которые векторы, направления которых не совпадают в общем слу­ чае с направлениями нормалей к граням. Из условия равновесия сил, действующих на тетраэдр, имеем

pnS— p.v5cos1

py Scos

pzS

cos

) + V q = 0 ,

где q вектор

объемных сил (например;

силы

тяжести)

и сил

инерции на единицу объема тетраэдра;

V — его полный объем. Раз­

делив это равенство на S, устремим все ребра тетраэдра

к нулю.

Учитывая, что при этом отношение V/S

стремится

к нулю

(числи-

тель — малая

третьего порядка,

а знаменатель — второго), полу­

чаем

>

 

 

 

PH =

PA-COS ( i O ) + Py cos (nTy)-f-p;Cos

(tCz).

Векторы p.v, ру, p z можно представить в виде

 

 

Px =

Wxx +

iPxy+^Pxz,

 

 

Py =

iPyx + iPyy + kpv2,

 

 

Pz =

iPzx +

iPzy+kpz2-

 

Соста вляющие рхх, рху, Pxz вектора р.-с єсть

силы,

действующие

на единицу

поверхности площадки, перпендикулярной

оси х в на­

правлениях

осей х, у, z, т. е. рхх — нормальное,

а рху

и pXz — каса­

тельные напряжения по этой площадке.

 

 

Таким образом, в силу определения, данного в § 3, напряженное состояние в произвольной точке М пространства занятого сплош­

ной средой есть тензор

 

 

 

Рхх

Рху Pxz'

 

П = / V

Руу

PyZ • ,

(5.1)

Pzx

Pzy

Pzz.

 

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

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

Сумма сил, приложенных

к субстанции,

заполняющей

этот

объем в любой момент, включая силы инерции, должна быть

равна

нулю

 

 

 

J p(F-w)dV

+ <$pndS=0.

 

(5.2)

V

s

 

 

Здесь р — плотность среды;

w ускорение

ее частиц; F уско­

рение объемных сил (например, силы тяжести). Момент указан­ ных сил относительно любой точки пространства, например, отно­

сительно начала

координат

0,

также

должен

быть равен

нулю.

А так как момент силы

Р есть

векторное произведение этой

силы

на плечо (радиус-вектор)

г точки

ее приложения,

то это второе ус­

ловие записывается так:

 

 

 

 

 

 

 

 

J

р [г, F -

w] d V+

j> [г,

р„] dS=0.

 

(5.3)

V

 

 

 

S

 

 

 

 

По формуле

(3.23)

 

 

 

 

 

 

 

 

 

Ф Ря ^

=

J dlv ПrfV.

 

 

(5.4)

 

5

 

 

 

v

 

 

 

 

Поэтому и в силу произвольности

объема

V,

уравнение

(5.2)

будет иметь вид

 

 

 

 

 

 

 

 

 

 

p ( F - w ) + divn=0 .

 

 

(5.5)

3*

-35

Это есть первое основное уравнение механики сплошной среды. Условие (5.3) равенства нулю момента приложенных к объему сил показывает, что П в уравнении (5.5) есть симметричный тен­ зор. Действительно, согласно формулам (3.23) и (3.21), а также вследствие того, что векторное произведение дифференцируется по тем же правилам, что и обычное, имеем

 

 

 

У \

 

дх

*

ду

'

дг

 

 

 

i l l '

дх

1

ду

 

 

 

+ . дх ' P v

+ " дудт'

+

+ dzдг

p]}dV=U[r,

 

 

divn] +

[i, p,] + [j,

p,l +

[k,

pz]}d V.

Теперь уравнение

(5.3)

можно записать так:

 

 

J [г, P ( F - w ) + d i v n ] d \ / + j ( i )

p,] +

[j, p,] + [k,

 

p,]}dV=0.

V

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.6)

Из

(5.5) следует, что первый интеграл

в (5.6) равен нулю. Атак

как объем V произволен, то уравнение (5.6) сводится

к

уравнению

 

 

 

[i.

P.v] + [J,

Py] + [k,

Р , ] = 0 .

 

(5.7)

Но

уравнение

(5.7)

и

есть

условие

симметричности тензора.

Чтобы в этом убедиться, умножим

(5.7)

векторно на

произвольный

вектор

а. Учитывая,

что

согласно

формулам

(1.3) и (1.4)

 

 

 

[і,

р * ] Х а = р Л * . а ) - і ( р Л

а)

 

 

и что аналогичные соотношения имеют место для других членов

(5.7), получаем

 

 

РЛі» a)+py (j, a)+p,(k,

а ) - і ( р х ,

a ) - j ( p , , a ) - k ( p z , a ) = 0 ,

откуда, согласно формулам

(3.14) и

(3.15), имеем

с , а ) - ( П ,

а ) = 0 .

А так как а — произвольный вектор, то это значит, что П

С = П,

т. е. тензор П совпадает со своим сопряженным, и значит он

сим­

метричен

Рху==РуХ> Pyz==Pzy> Pzx===Pxz-

Если рассматривается жидкая (капельная или газообразная) среда, то, положив в формуле (2.12) F = v, получим следующее вы­ ражение для ускорения:

Основное уравнение (5.5) для этого случая будет иметь вид

- g - + ( v , V v ) = F + - I - d l v I I .

(5.9)

Следует заметить, что

оператор (v, V ) часто

обозначают как

(v, grad), тогда уравнение

(5.9) будет иметь вид

 

- ^ - + (v, grad)v = F + -£- divH.

Важный частный случай жидкой среды есть идеальная жид­

кость — среда, неспособная

к

восприятию касательных напряже­

ний. В идеальной

жидкости

р х у

= pyx = pVz = Pzy = pzx = Pxz:=0, а век­

торы р п , рл;, ру,

pz нормальны

к граням тетраэдра соответственно

ABC, МВС, МАС, МАВ

(рис. 6), при этом

 

 

Vx =

Wxx, Py = j/>yy.

Pz=kPzz,

 

P„=iA™cos ( n T x ) - f j^yycos (тСу) + кргг cos (пГ*)-

(5.10)

Проектируя равенство (5.10) последовательно на направления

х, у, г, получаем

 

 

 

 

Pn=Pxx=Pyy=Pzz,

 

(5-11)

т. е. нормальное напряжение не зависит от направления площадки, по которой оно действует, иными словами, в данном случае тензор П есть шаровой тензор

\-р

0

01

п =

 

О

о

,

(

о

0

~р\

где — р — инвариантное по

отношению

к

ориентации площадки

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

Если

П — шаровой тензор,

то, согласно

формулам

для

дивер­

генции

тензора,

d i v l l = gradp

и уравнение

(5.9)

для

идеальной

жидкости будет выглядеть так:

 

 

 

 

 

 

 

- 5 T + ( v , V ) v = F — J - g r a d j * .

 

 

(5.12)

Это уравнение

называется уравнением Эйлера

(точнее,

уравне­

ниями Эйлера называются три проекции этого векторного уравне­ ния на три координатные оси).

Если жидкость не идеальна, и в ней могут существовать ка­ сательные напряжения, то соотношение (5.11) уже не имеет места.

Ho pxx+Pyy + Pzz есть инвариант тензора

П, т. е. скаляр. Положим

в этом случае

 

Р = --Ъ-(Рхх+Руу+Ргг)

(5-13)

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

Рхх=—Р+*хх,

Руу=—Р+*уу.

Ргг=—Р

+ 1г*>

Рху ==/^'y.v= = Zxy= = ^у і-1 Pyz'

Pzy ' ^уг= =

^гуі

 

/'г.ї —Агг= = xzx

==тл"г •

 

Теперь уравнение (5.9) можно записать в виде

d v ' (г>, g r a d ) x » = F — f g r a d p + - f div7\

(5.14)

где

^xy

7=={ XyX Lyy

г-

''zx ^y ••zz

— тензор, порожденный отклонениями свойств жидкости от иде­ альной. Скаляр р в этом уравнении по-прежнему называется дав­ лением.

Второе основное уравнение сплошной среды дает закон сохране­ ния массы. Рассмотрим действие этого закона в условиях проте­ кания жидкости через тот же объем V, который имелся в виду при выводе первого основного уравнения. За время dt через поверх­ ность 5, ограничивающую этот объем, вытекает наружу масса жид­ кости

<J>(n,

pv)dSdt.

s

 

С другой стороны, за то же время масса жидкости в объеме V изменяется на

1-%-dVdt. v

Общее приращение массы за время dt есть сумма этих двух вы­ ражений. По закону сохранения вещества это приращение равно нулю, следовательно,

j-|*-rfVcW+<j>(n, pv)dSdt=0.

Сокращая это уравнение на dt и преобразуя затем поверхност­ ный интеграл в объемный по формуле (2.27), получаем

J ^ + d i v P v = 0 .

(5-15)

Это уравнение называется уравнением неразрывности. Оно мо­ жет быть записано также в форме

д?

•V • grad p-j-p div v = 0 ,

dt

 

или, согласно формуле (2.10),

1

RFP ' - d i v v = 0 .

(5.16)

p

dt

 

Обращаясь к физическому смыслу формул (2.9) и (2.10), за­ ключаем, что dp/dt характеризует изменение плотности жидкой ча­ стицы при ее перемещении в потоке. Если жидкость несжимаема, то плотность каждой ее частицы остается постоянной, т. е.

- ^ = = 4 - + v - g r a d p = 0 .

(5.17)

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

 

div v = 0 .

 

При этом, конечно, ни dp/dt, ни grad р, вообще говоря,

в нуль

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

жидких

частиц не означает неизменность плотности жидкости в данной точке пространства, если плотность этих частиц почему-либо раз­ лична.

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

А. Несжимаемые жидкости:

1)жидкость постоянной плотности;

2)общий случай несжимаемой жидкости. Б. Сжимаемые жидкости:

1)баротропная жидкость. Плотность ее зависит только от дав­ ления р = р(р). Простейший пример — упругая капельная жидкость, подчиняющаяся при сжатии закону Гука;

2)бароклинные жидкости (в частности, все газы). Плотность

еезависит не только от давления.

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

6. Уравнения Навье—Стокса

Векторное уравнение (5.14), эквивалентное трем скалярным, и скалярное уравнение (5.15) образуют систему четырех скалярных уравнений, в общем случае с одиннадцатью неизвестными: и, v, w (так мы обозначим здесь компоненты вектора скорости v по осям X, у, Z), р, р, Гхх, tyy, Tzz, txy, Тх г , Xyz. Общие ЗЭКОНЫ МеХЭНИКИ, ПрИВОдящие к этим уравнениям, никаких других независимых уравнений

Соседние файлы в папке книги из ГПНТБ