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

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

.PDF
Скачиваний:
23
Добавлен:
06.03.2016
Размер:
1.41 Mб
Скачать

v2

gz

C

1

C .

0

0

 

 

1

2

0

2

3

Вычтем это выражение из (а):

 

v2 v2

 

 

 

 

 

 

 

 

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

g(z z

 

) C

 

 

 

 

 

 

 

 

0

 

 

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1

1

 

 

1

 

 

g(z z

 

)

v

2 v02

 

.

 

 

 

 

(в)

 

 

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C2 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдѐм C , чтобы сюда вставить. Из (2.14) и (2.15) имеем C

 

 

p0

,

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

c2

 

 

 

 

 

 

 

 

 

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

0 г . Отсюда C

 

 

 

 

 

г

 

 

. Подставим в (в):

 

 

 

 

 

 

0

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

v2

 

v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

0 1

 

 

g(z z0 )

 

 

 

 

 

 

 

0

.

 

 

 

 

 

 

 

(г)

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что величина

1

очень мала (т.к. 1 близко к нулю, а ско-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рость звука cг

велика). Кроме того, скорость течения обычно мала по срав-

нению с cг , перепад высот z z0

невелик.

Поэтому в фигурных скобках

второй член мал и выражение (г) допускает следующее приближение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

v2 v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

 

 

 

g(z z0 )

 

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cг

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

частности,

для

горизонтального канала

(когда z z0 )

будем

иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v2

v2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2cг

 

 

 

 

 

 

 

31

Глава 3

ДИНАМИКА ВЯЗКОЙ ЖИДКОСТИ

3.1. Понятие о вязкости

Реальная жидкость и газ имеют вязкость.

Вязкость – это сила трения между слоями текущей жидкости.

Именно из-за вязкости текущая жидкость, предоставленная самой себе (т.е. когда убраны причины, вызвавшие движение), останавливается.

Рассмотрим два опыта, в которых проявляется внутреннее трение (вязкость).

О п ы т 1. В вертикальную трубку с краном вни-

 

зу нальѐм воду, а на неѐ сверху осторожно доль-

 

ѐм коричневое растительное масло (рис.

3.1). В

 

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

 

ризонтальной (пунктирная линия).

 

 

 

Откроем кран так, чтобы течение было очень

 

медленным. Вскоре граница раздела примет

 

форму параболоида вращения. Слой жидкости,

 

прилегающий к стенке трубки, неподвижен. Ско-

 

рость течения остальных слоѐв увеличивается по

 

мере приближения к оси трубки.

 

 

 

О п ы т 2. Поместим в жидкость две парал-

 

лельные пластины одинаковой площади S , рас-

 

стояние между которыми y, причѐм y

 

 

 

S (рис.

 

3.2).

 

 

 

Будем перемещать вправо верхнюю пластину

Рис. 3.1

со скоростью v. Убеждаемся, что для перемеще-

 

32

 

 

 

ния пластины с постоянной скоростью, к ней нужно приложить определѐнную постоянную силу F. Раз пластина движется не ускорен-

но, значит, на неѐ действует нулевая суммарная сила. Следователь-

но, внешнюю силу F, которую мы прикладываем, уравновешивает

противоположно направленная сила трения жидкости о пластину.

Обозначим еѐ Fтр.

Выполняя этот опыт при различных S , y, v (как это делал Ньютон), можно убедиться, что

- растёт при увеличении площадипластины S; Fтр : - растёт при увеличениискорости перемещения v;

- уменьшается при увеличениирасстояния y.

Рис. 3.2

Значит,

 

 

 

F S

v

.

(3.1)

 

тр

y

 

 

 

Коэффициент называется динамической вязкостью; он зави-

сит от типа жидкости и еѐ состояния (например, от температуры). Опыт показывает, что при нагревании вязкость жидкости уменьшается, а газов – растѐт.

Относительно верхней пластины нижняя пластина движется с той же скоростью v, но в обратную сторону, влево. Значит, на ниж-

нюю пластину действует сила Fтр .

Формула (3.1) определяет не только силу трения, действующую на пластину, но и силу трения между соприкасающимися слоями жидкости.

Для слоѐв, находящихся на расстоянии dy, скорости отличаются

на dv. Поэтому формулу (3.1) можно записать так:

 

F S

dv

.

(3.2)

 

тр

dy

 

 

 

33

3.2.Течение жидкости в круглой трубе

Всоответствии с рис. 3.2, при течении в круглой трубе скорость жидкости равна нулю около стенок и максимальна на оси трубы (рис. 3.3). Иначе говоря,

при r R должно быть v 0.

(а)

при r 0 должно быть v v0 .

(б)

Найдѐм закон изменения скорости v частиц жидкости в зависимости от их расстояния r до оси трубы радиуса R.

В трубе выделим воображаемый цилиндрический объѐм жидкости радиуса r и длины L (рис. 3.4).

Рис. 3.3

Рис. 3.4

У цилиндра площадь боковой поверхности равна S (2 r)L, по-

этому со стороны внешней среды на боковую поверхность действует сила

Fтр (3.2) 2 rL dvdr .

На торцы цилиндра действует сила давления p1 r2 p2 r2 ( p1 p2 ) r2 .

Цилиндр движется равномерно (не ускоренно), поэтому действующая на цилиндр суммарная сила равна нулю:

2 rL dvdr ( p1 p2 ) r 2 0.

Решаем полученное дифференциальное уравнение:

2 L

dv

( p p

 

) r,

 

 

 

 

2

 

 

 

 

dr

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv

( p1 p2 )

r dr, v

( p1 p2 )

r dr,

 

 

 

 

2 L

2 L

 

 

 

 

 

 

 

34

 

 

 

 

 

 

 

 

 

 

 

 

 

( p

p

2

)

 

 

 

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

v

 

1

 

 

 

 

 

 

 

 

 

 

 

C.

(в)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2L

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения неизвестной константы C воспользуемся началь-

ным условием (а):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( p p

 

)

 

R2

 

 

 

 

 

 

( p p

2

)

 

 

 

R2

 

0

1

 

2

 

 

 

C. Отсюда C

 

 

 

1

 

 

 

 

 

 

 

 

 

.

 

2L

 

 

2

 

 

 

2L

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим в (в):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

( p1 p2 )

(R2

 

 

r 2 ).

(3.3)

 

 

 

 

 

 

 

 

 

 

4L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Условие (б) даѐт

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

( p1 p2 )

R2 .

(3.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

4 L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поэтому

( p1 p2 )

 

v0

и формула (3.3) запишется так:

 

 

R2

 

 

 

 

4 L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v v

1

 

 

 

 

 

 

.

 

 

 

(3.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

R2

 

 

 

 

 

 

 

 

Найдѐм поток вектора скорости сквозь поперечное сечение трубы (т.е. найдѐм объѐмный расход жидкости), применив формулу (2.2):

V v dS v0

 

 

r 2

 

1

 

 

dS.

R

2

( S )

( S )

 

 

 

 

Так как (S ) – круг, вычисления удобно вести в полярных координатах. Тогда

dS r d dr,

(S) : 0 2 ,

0 r R.

V v0

 

 

r 2

 

 

2

R

 

r 2

 

 

R2

1

 

 

 

 

r d dr v0

d 1

 

 

 

 

r dr

v0 .

R

2

R

2

( S )

 

 

 

 

 

0

0

 

 

 

 

2

Подставив сюда значение (3.4), получим

 

V

R 4

 

p p

 

 

 

 

 

1

2

 

(3.6)

 

 

8

 

L

 

 

 

 

 

 

 

 

 

Формула Пуазейля

 

По формуле Пуазейля можно определить вязкость жидкости:

R4 p1 p2 . 8V L

35

Числовые значения в правой части определяются из условий эксперимента.

3.3. Ламинарное и турбулентное течения

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

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

или ламинарным.

Если сечение горизонтальной трубки увеличивать или увеличи-

вать скорость течения, появится не-регулярное движение частиц

жидкости: окрашенная струйка сначала начнѐт дрожать,

а потом

хаотично перемешиваться с прозрачными стру-

Рис. 3.5

ями. Такое течение называют турбулентным.

Формула Пуазейля справедлива только для ламинарного течения.

Введѐм безразмерную величину

 

Re

vR

 

(3.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Число Рейнольдса

 

 

 

 

 

 

 

 

Переход ламинарного течения в турбулентное происходит при некотором критическом числе Рейнольдса Reкр . Это значение силь-

но зависит от формы входной части трубы:

1200 Re

кр

 

20000 .

 

 

 

для закруглённого входа

При стационарном турбулентном течении скорость в данной точке

36

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

значение направлено вдоль оси трубы. Средняя скорость остаѐт-

vср

ся постоянной по всему сечению, и только в тонком пограничном слое у стенки трубы падает до нуля (рис. 3.6).

Рис. 3.6

На практике при турбулентном течении используется формула

 

 

 

 

 

 

 

 

 

 

 

p p

2

 

 

 

 

(vср )турб2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(v

ср

)2

 

 

R

 

p1

p2

,

 

 

 

(3.8)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

турб

 

 

 

 

 

k

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в которой k безразмерный гидравлический коэффициент.

Найдѐм соответствующую

 

 

формулу для ламинарного течения.

Так как V (v

 

)

лам

 

S,

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ср

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

(3.6)

 

 

 

R4

 

p p

 

 

1

 

 

 

(v

)

лам

 

 

 

 

 

 

 

 

 

1

2

 

 

,

 

 

 

ср

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

8

 

L

 

 

R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2

 

 

p p

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(v

ср

)

лам

 

 

 

 

 

 

 

 

1

 

.

 

 

 

(3.9)

 

 

 

 

 

 

 

 

8

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сравнение формул (3.8) и (3.9) показывает, что повышение скорости прокачки жидкости по трубам потребует при турбулентном течении большего перепада дав-

ления

p1 p2

, чем при ламинарном. Гра-

Рис. 3.7

L

 

 

 

фик зависимости скорости прокачки от перепада давления показан на рис. 3.7.

3.4.Тело в потоке вязкой жидкости

Вреальной, вязкой жидкости результирующая сила, действующая на тело, не равна нулю, т.е. поток всегда действует на тело с некоторой силой. Если ось симметрии тела направлена вдоль потока, то сила, с которой поток действует на тело, называется силой ло-

бового сопротивления.

Рассмотрим обтекание потоком шара радиуса r. Зависимость силы Fс лобового сопротивления от числа Рейнольдса показана на

рис. 3.8.

37

При малых скоростях течения, когда Re 102 , сила пропорциональна скорости, F ~ v. Это происходит потому, что из-за наличия

тонкого пограничного слоя на шар действуют силы вязкости. Толщина пограничного слоя оценивается по формуле

 

 

r

 

 

 

.

 

 

 

Re

Рис. 3.8

Например, в конце линейного участка, где Re 102 , толщина10r – в 10 раз меньше радиуса шара. Вне этого слоя жидкость те-

чѐт так же, как и невязкая, симметрично обтекая шар спереди и сзади.

При малых числах Рейнольдса сила лобового сопротивления, действующая на шар, рассчитывается по формуле Стокса:

Fc 6 vr.

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

m dvdt mg FA Fc .

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

случае v const,

dv / dt 0, поэтому

mg FA Fc 0.

(а)

38

Но m V

 

 

 

4

r3 , F

 

 

V

 

 

 

4

r 3 ,

F 6 vr. Подставив в

 

 

 

 

 

ш ш

 

ш

 

3

 

A

 

ж ш

 

ж

 

3

 

c

(а), получим ш 43 r3 ж 43 r3 6 vr 0, отсюда

2r 2 ( ш ж ). 9v

Определив в эксперименте равномерную скорость v, по этой формуле рассчитываем вязкость.

При больших скоростях потока, когда Re 102 , симметрия обтека-

ния относительно вертикальной оси (на рис. 3.9 – пунктирная линия) нарушается – позади шара происходит отрыв линий тока и образуется завихрѐнное пространство. Симметрии давления в точках M и M 1 уже нет.

По формуле Бернулли имеем

v2

 

 

v2

 

.

 

 

 

 

 

 

0 p

0

 

K p

K

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В точке K поток останавливается, v

 

0, поэтому

p

 

 

v2

 

.

K

K

0 p

 

 

 

 

 

 

 

 

2

0

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3.9

В точке K1

давление равно p0 .

Поэтому результирующее давле-

ние

 

на

шар

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

потока

будет пропорционально

 

 

 

 

v2

 

 

v2

S, где S r2 – площадь

p p

0

 

 

0

, а сила пропорциональна

0

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сечения шара, перпендикулярного потоку. Следовательно, сила лобового сопротивления равна

v2

Fc CS 20 ,

где C коэффициент лобового сопротивления.

39

3.5. Уравнение неразрывности

В пространстве с движущейся жидкостью поместим воображаемую неподвижную замкнутую поверхность (dS), проницаемую для

частиц жидкости (рис. 3.10). Область пространства, оказавшуюся внутри (dS), обозначим через (dV ). Итак, (dV ) неподвижная об-

ласть в пространстве, сквозь которую движутся частицы жидкости.

Пусть в момент времени t

в области

dV имеется масса жидкости

dV.

 

За время dt

эта масса увеличится на

 

d( dV )

 

dV const

 

d dV ,

(а)

 

 

где (x, y, z, t).

 

(б)

Поделив (а) на dt, получим массу dm,

Рис. 3.10

притекающую в область dV за единицу времени:

dm

d dV

 

d

dV.

 

(в)

 

 

 

 

dt

dt

 

 

С другой стороны, масса жидкости, вытекающей из области dV за единицу времени, определяется по формуле (2.6). Значит, величины (2.6) и (в) одинаковы, но противоположны по знаку:

 

 

 

 

 

 

t

dV div( v ) dV ,

отсюда

 

 

d

 

 

 

 

div( v ) 0.

(г)

 

dt

Из (б) найдѐм полную производную по времени:

 

d

 

dx

 

 

dy

 

 

dz

.

 

 

 

 

 

 

 

 

dt

x dt

 

y dt

z dt

t

Однако наша область dV

 

неподвижна, координаты x,

чек не зависят от времени:

 

dx

0,

dy

0,

dz

0. Поэтому

 

dt

dt

dt

 

 

 

 

 

 

 

 

 

 

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

40

y, z еѐ то-

d . dt t