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

лекция 4

.pdf
Скачиваний:
7
Добавлен:
18.03.2015
Размер:
824.46 Кб
Скачать

Лекция 4

Наименование

 

 

 

 

 

 

 

(час)

 

Содержание теоретической части

 

Затраты времени

разделов

Номер

 

 

 

Самостоятель-

 

разделов дисциплины

 

лекции

Лекции

ная (внеаудитор-

дисциплины

 

 

 

ная) работа

 

 

 

 

 

 

 

студента

ОСНОВНЫЕ

Свойства напряжений поверхностных сил.

4

2

0,5

УРАВНЕНИЯ

Давление и его свойства.

 

 

 

 

 

ГИДРОГАЗО-

Нормальные и касательные напряжения. Закон

 

 

 

ДИНАМИКИ

парности касательных напряжений. Модель

 

 

 

 

идеальной жидкости. Давление и его свойства.

 

 

 

 

Уравнение движения (количества движения или

 

 

 

 

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

 

 

 

 

Обобщенная гипотеза Ньютона о связи между

 

 

 

 

напряжениями

и скоростями деформаций

 

 

 

 

(закон Стокса).

 

 

 

 

 

 

 

Уравнения движения Навье-Стокса для

 

 

 

 

вязкой сжимаемой жидкости.

 

 

 

 

 

Уравнения движения идеальной жидкости в

 

 

 

 

форме Эйлера и Громеки-Лэмба.

 

 

 

 

 

Интегралы

уравнений

движения

идеальной

 

 

 

 

жидкости: для безвихревого неустановившегося

 

 

 

 

движения (интеграл Лагранжа или Коши-

 

 

 

 

Лагранжа),

для

установившегося

движения

 

 

 

 

(интеграл Бернулли). Уравнение Бернулли.

 

 

 

 

Интегральная

форма

закона

сохранения

 

 

 

 

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

 

 

 

 

объема (первое уравнение Эйлера). Переход от

 

 

 

 

жидкого объема к контрольному.

 

 

 

 

 

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

 

 

 

 

тело, по состоянию потока на границах.

 

 

 

 

Полный импульс потока в сечении. Определение

 

 

 

 

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

 

 

 

 

канала со стороны текущей по нему жидкости;

 

 

 

 

учет сил давления на канал со стороны

 

 

 

 

окружающей среды.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С.М. Мухаметшин. Лекции по гидрогазодинамике

~ 1 ~

Свойства напряжений поверхностных сил. Давление и его свойства.

Возьмем какую-либо точку M сплошной среды, проведем из этой точки координатные оси и построим на них бесконечно малый тетраэдр ABCM с ребрами dx , dy , dz . Ориентировка грани ABC произвольна и задается вектором

нормали n . Обозначим площадь грани ABC через dS , а площади трёх граней тетраэдра, лежащих в координатных плоскостях через dFx , dFy и dFz .

Рис. 4.1

Составим векторное уравнение движения жидкого тетраэдра, выражающее второй закон Ньютона. При этом учтем, что положительными принято считать напряжения, направленные в сторону внешней нормали, а орты координатных осей, являющиеся нормалями к площадкам dFx , dFy и dFz , направлены внутрь

выделенного жидкого объема.

На тетраэдр ABCM будут действовать поверхностные силы, характеризуемые напряжениями px , py , pz и pn , приложенными на

соответствующих гранях тетраэдра, а также массовая сила

dR = F dm = ρFdV = ρF 13 hdS ,

где R напряжение массовой силы, м/с2; dm масса в объеме тетраэдра dV , h высота тетраэдра.

В соответствии со вторым законом Ньютона сумма сил, действующих на тетраэдр, равна произведению его массы на ускорение, то есть

С.М. Мухаметшин. Лекции по гидрогазодинамике

~ 2 ~

1

r

r

r

r

r

=

dm =

1

ρh

dS

 

3

ρhFdS + pndS - pxdFx

- pydFy

- pzdFz

dt

3

dt

(4.1)

 

 

 

 

 

 

 

 

 

 

Площади граней, параллельных координатным осям, можно вычислить как площади проекций грани ABC на соответствующие координатные плоскости:

dFx

= dS cos(n, x) = nxdS ;

 

dFy

= dS cos(n, y) = nydS ;

(4.2)

dFz

= dS cos(n, z) = nz dS .

 

Здесь nx = n × i = ni cos(n, x) = cos(n, x), аналогично ny , nz .

В результате, подставив (4.2) в (4.1) сокращаем все члены равенства (4.1) на dS и стягивая тетраэдр в точку (то есть полагая h → 0), получим

pn = nx px + ny py + nz pz .

(4.3)

С другой стороны, каждый из векторов px , py , pz может быть представлен

тремя составляющими, действующими вдоль координатных осей.

Обратимся к рис. 4.1. На грань dFx действует напряжение px , для которого

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

px = x + xy + xz ;

аналогично

py = yx + y + yz ;

(4.4)

pz = zx + zy + z .

 

Первый индекс указывает ориентацию площадки (нормаль к ней), а второй ось, на которую проецируется вектор. Для нормального напряжения второй индекс совпадает с первым, поэтому может опускаться.

Теперь, учитывая (4.4), спроецируем вектор

pn на координатные оси.

Получим:

 

pnx = nxσ x + nyτ yx + nzτ zx

 

pny = nxτ xy + nyσ y + nzτ zy

(4.5)

pnz = nxτxz + nyτ yz + nzσ z

 

Таким образом, напряженное состояние в точке определяется совокупностью трех векторов напряжения или их девятью компонентами, определенными на трех взаимно перепендикулярных площадках. Соотношение (4.5) является определением тензора (от англ. tension - напряжение).

Компоненты образуют тензор второго ранга, которому можно поставить в соответствие матрицу:

С.М. Мухаметшин. Лекции по гидрогазодинамике

~ 3 ~

 

æ

σ x

τ yx

τ zx

ö

 

pni

ç

 

 

σ y

 

 

÷

(4.6)

= çτ xy

τ zy ÷ .

 

ç

τ

xz

τ

yz

σ

z

÷

 

 

è

 

 

 

ø

 

Диагональные компоненты (с одинаковыми индексами) называются нормальными напряжениями, а остальные компоненты касательными напряжениями, или

напряжениями сдвига.

Тензор напряжений (4.6) зависит от координат и времени и образует тензорное поле.

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

τ xy yx ; τ yz zy ; τ zx xz ,

(4.7)

называемые законом парности касательных напряжений. Соответствующее доказательство можно найти, например, в [1, 2] из списка основной литературы. Из этого закона следует, что тензор напряжений (4.6) является симметричным. Последнее означает, что тензор напряжений содержит только шесть различных компонент. В соответствии с этим уменьшается число неизвестных в уравнениях движения.

Рассмотрим теперь модель идеальной жидкости. Идеальной жидкостью

называется изотропная сплошная среда, в которой отсутствуют касательные напряжения, то есть τij = 0 (i ¹ j ). Касательные напряжения в жидкости

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

Тогда, из (4.5) следует, что

pnx = nxσ x ,

pny

 

= nyσ y , pnz = nzσ z .

(4.8)

С другой стороны,

 

 

 

 

 

 

 

 

 

 

pnx

= pn × i =

 

pn

 

×

 

i

 

cos(n, x) = pnnx .

 

 

 

 

 

 

 

 

 

Аналогично,

= pnny , pnz = pnnz .

 

pny

(4.9)

Сравнивая (4.8) и (4.9), видим, что

 

 

 

 

 

 

 

 

 

pn = σ x = σ y

= σ z .

(4.10)

Жидкий объем в силу своего строения не может воспринимать растягивающие напряжения. Поэтому действующие на него нормальные напряжения должны быть сжимающими. Такое сжимающее нормальное

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

обозначается p . Следовательно,

С.М. Мухаметшин. Лекции по гидрогазодинамике

~ 4 ~

p = − pn = −σ x = −σ y = −σ z

(4.11)

Знак «минус» показывает, что давление всегда направлено по внутренней нормали к поверхности выделенного жидкого объема. Это первое свойство давления.

Обозначая px = −σ x , py = −σ y , pz = −σ z , можно записать p = px = py = pz .

Осюда следует второе свойство давления: в любой точке идеальной или

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

Но в общем случае p = p(x, y, z,t) , то есть давление может изменяться в

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

Уравнение движения жидкости в напряжениях.

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

Количество движения K массы жидкости в объеме V равно

K = òρWdV ,

V

где W - скорость движения центра масс жидкого объема dV . Тогда

dK

 

d

 

r

d

r

dW

r

d

(ρdV ).

 

=

 

òρWdV = ò

 

ρWdV = òρ

 

dV + òW

 

dt

dt

dt

dt

dt

 

V

V

V

V

 

Так как ρdV - масса жидкого объема dV , не изменяющаяся со временем, то есть

 

 

 

d

(ρdV ) = 0 ,

(4.12)

 

 

 

dt

 

 

 

 

 

r

d

(ρdV ) = 0 .

 

то соответственно и òW

 

 

dt

 

V

 

 

 

 

Тогда уравнение количества движения жидкого объема V в интегральной форме запишется в виде:

С.М. Мухаметшин. Лекции по гидрогазодинамике

~ 5 ~

òρ

dW

r

r

(4.13)

dt

dV = òρRdV + ò pndF , Н

V

V

F

 

Справа в (4.13) стоит сумма внешних массовых и поверхностных сил, действующих на жидкий объем V и ограничивающую его поверхность F .

Уравнение (4.13) представляет собой основное динамическое соотношение ГГД в интегральной форме. Его часто называют также уравнением импульсов, так как оно преобразуется к последнему простым умножением на dt .

Преобразуем поверхностный интеграл, входящий в уравнение (4.13), в объемный. С учетом (4.3):

é

cos(n, x) + py

cos(n, y) + pz

ù

(4.14)

ò pndF =òë px

cos(n, z)û dF .

FF

Ввекторном анализе доказывается, что для произвольной функции a (векторной или скалярной) справедливы следующие соотношения:

 

r

 

a

dV ;

òacos(n, x)dF = ò

x

F

V

 

 

 

 

r

ò

a

dV

òa cos(n, y)dF =

y

F

V

 

 

r

 

a

dV

òacos(n, z)dF = ò

z

F

V

 

 

 

Тогда равенство (4.14) с учетом (4.15) запишется в виде

r

é

r

 

py

 

r

ù

px

 

 

pz

ò pndF =òê

 

+

 

+

 

ú dV .

x

y

z

F

V ë

 

 

û

(4.15)

(4.16)

Подставляя (4.16) в (4.13), в силу произвольности объема V можем приравнять подынтегральные выражения. Разделив обе части уравнения на ρ ,

получим векторную форму дифференциального уравнения движения в напряжениях:

r

1

æ

r

 

 

 

py

 

 

 

r

 

 

ö

 

dW

 

 

p

 

 

 

 

 

p

z

 

 

 

(4.17)

R +

 

ç

 

x

+

 

 

 

+

 

 

 

÷ =

 

 

 

ρ

 

 

y

z

dt

 

 

è

x

 

 

 

 

 

ø

 

 

 

В проекциях на оси декартовой прямоугольной системы координат,

учитывая (4.5), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X +

1

æ

σ

x +

 

τ yx

+

τ

ö

=

u

;

 

 

 

ç

 

 

 

 

 

 

 

 

 

zx

÷

 

 

 

 

 

 

 

y

 

 

 

 

 

t

 

 

ρ è

x

 

 

 

 

 

z ø

 

 

~ 6 ~

С.М. Мухаметшин. Лекции по гидрогазодинамике

Y + ρ1 X + ρ1

æ

τ

xy

+

σ

y

+

τ

zy

ö

=

ç

 

 

 

 

÷

 

 

 

y

z

 

è

x

 

 

 

 

ø

 

æ

τ

xz +

 

τ yz

 

+

σ

 

ö

=

ç

 

 

 

 

 

 

 

z ÷

 

 

y

 

z

è

x

 

 

 

 

ø

 

υ

;

(4.18)

t

 

 

wt .

Обобщенная гипотеза Ньютона о связи между напряжениями и скоростями деформаций (закон Стокса).

Обобщение гипотезы (закона) Ньютона (1.7) для жидкого объема приводит к закону трения Стокса:

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

При этом для касательных напряжений, учитывая (4.7), можно показать, что

τ xy

=τ yx = 2μθz

 

æ

υ

 

+

u

 

ö

 

= μ ç

 

 

÷

 

 

 

 

 

è

 

x

 

 

 

 

y

 

ø

 

τ yz

=τ zy = 2μθx

 

æ

 

w

+

υ

ö

(4.19)

= μ ç

 

÷

 

 

 

 

è

 

y

 

 

 

 

z

ø

 

τ zx

=τ xz = 2μθy

 

æ

u

+

 

w

ö

 

= μ ç

z

 

 

x

÷

 

 

 

 

 

è

 

 

 

 

ø

 

Для нормальных напряжений, возникающих от действия сил вязкости,

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¢

= 2μεx

= 2μ

∂u

;

 

 

 

 

 

 

 

 

 

 

 

 

 

σ x

x

 

 

 

 

 

σ ¢y

= 2με y

= 2μ

 

υ

 

;

 

 

(4.20)

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¢

= 2μεz

= 2μ

 

∂w

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ z

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

виде:

σ x

= − p0

 

+ σ x ;

 

σ y

= − p0 + σ ′y ;

(4.21)

σ z

= − p0

 

+ σ z .

~ 7 ~

С.М. Мухаметшин. Лекции по гидрогазодинамике

где p0 - некоторая скалярная величина, аналогичная давлению в идеальной

жидкости.

В тензорном анализе доказывается, что сумма нормальных напряжений (σ x + σ y + σ z ) инвариантна относительно поворота осей координат, то есть не

зависит от направления, как и давление в идеальной жидкости.

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

p = -σ x +σ y + σ z . (4.22) 3

Тогда, решая систему (4.20) – (4.22), можно найти, что

σ x

æ u

-

1

r ö

;

= - p + 2μ ç

x

3

divW ÷

 

è

 

ø

 

σ y

æ

υ

-

 

r ö

(4.23)

= - p + 2μ ç

1 divW ÷

 

è

y

 

3

ø

 

σ z

æ w

-

1

r ö

 

= - p + 2μ ç

z

3

divW ÷

 

 

è

 

ø

 

В идеальной жидкости, когда μ = 0, из (4.23) следует (4.11).

Уравнения движения Навье-Стокса для вязкой сжимаемой жидкости.

Если в уравнении движения в напряжениях (4.18) подставить зависимости (4.19) и (4.23), то получим дифференциальные уравнения неустановившегося пространственного движения вязкой сжимаемой жидкости, называемые уравнениями Навье-Стокса (1845). Не приводя преобразований, запишем их в проекциях на координатные оси (при μ = const для всего поля течения):

du

 

 

 

1 p

 

1

 

 

r

dt

= X -

 

 

 

x

+νDu +

3ν

 

 

 

 

(div W

ρ

x

 

 

 

 

 

1

 

 

p

 

1

 

 

r

dt

= Y -

 

 

 

 

y

+νDυ +

3ν

 

 

 

(div W

 

ρ

 

y

dw

 

 

 

1

 

 

p

 

1

 

 

r

dt

= Z -

 

 

 

z

+νDw +

3ν

 

 

 

(div W

 

ρ

 

z

 

В векторной форме (4.24) примет вид:

dW

r

1

r

1

r

 

= R -

 

grad p +νDW +

 

ν grad(divW

dt

ρ

3

 

 

 

); ); ).

) , Н/кг.

(4.24)

(4.25)

Уравнение (4.25) представляет собой специфическую для вязкой сжимаемой

С.М. Мухаметшин. Лекции по гидрогазодинамике

~ 8 ~

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

Уравнения движения идеальной жидкости в форме Эйлера и Громеки-Лэмба.

Из уравнений (4.25) и (4.24), полагая ν = 0, получаем дифференциальные

уравнения Эйлера (1755) движения идеальной жидкости в векторной:

и координатной формах:

Согласно (2.16)

dW dt

dWdt = Rr - ρ1 grad p .

dudt = X - ρ1 px ; ddtυ = Y - ρ1 py ; dwdt = Z - ρ1 pz .

(

r r

)

r W

 

W +

 

=

W ×Ñ

t

 

 

 

 

В векторном анализе доказывается, что

æ

 

2 ö

r

r r

r

grad çW

÷

= (W

×Ñ)W + (W ´ rotW )

è

2

ø

 

 

 

Тогда можно записать, что

dW

æ

 

2 ö

r

r

W

= gradçW

÷

- (W ´ rotW )+

dt

è

2

ø

 

 

t

(4.26)

(4.27)

(4.28)

(4.29)

Подставляя (4.29) в (4.26), получим уравнение движения идеальной жидкости в форме Громеки-Лэмба (1881 г.):

r

1

æ

 

2 ö

r

r

W

 

R -

grad p = gradçW

÷

- (W ´ rotW )+

(4.30)

ρ

t

 

è

2

ø

 

 

 

С.М. Мухаметшин. Лекции по гидрогазодинамике

~ 9 ~

В координатной форме будем иметь

 

 

 

1 p

 

 

 

æ

 

2

X -

 

=

 

 

çW

 

 

 

 

 

 

 

 

ρ x

 

 

 

 

 

 

 

 

x è

2

 

 

 

1 p

 

 

 

æ

 

2

Y -

 

=

 

 

çW

 

 

 

 

 

 

 

ρ y

 

 

 

 

 

 

 

y è

2

 

 

1 p

 

 

 

æ

 

2

Z -

 

=

 

 

çW

 

 

 

 

 

 

 

 

 

 

ρ z

 

 

 

 

 

 

 

 

z è

2

 

ö

r

 

÷

- (W ´

ø

 

 

ö

r

´

÷

- (W

ø

 

 

ö

r

´

÷

- (W

ø

 

 

r

+ ut

 

 

rotW )x

;

 

r

+ υt

 

 

rotW )y

;

(4.31)

r ) + w rotW z t .

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

Интегралы уравнений движения идеальной жидкости

 

Пусть массовые силы имеют потенциал,

то есть R = -gradF .

Введем

некоторую функцию давления P(x, y, z) такую, что

 

 

dP = dp , то есть P = ò dp ,

(4.32)

ρ

ρ

 

Тогда, например, Px = ρ1 px . При этом

 

1

æ

 

 

 

ö

 

 

 

 

dp =

ç

p dx +

p dy +

p dz ÷

=

P dx +

P dy +

P dz = dP ,

 

ρ ρ è

x

y

z

ø

 

x

y

z

то есть получаем соотношение (4.32).

 

 

 

В векторном виде имеем, что

1

 

 

grad P =

grad p .

(4.33)

 

 

ρ

 

Тогда уравнение Громеки-Лэмба (4.30) можно представить в виде

æ

 

2 ö

r

r

W

= 0

gradçF + P + W

÷

- (W ´ rotW )+

è

2

ø

 

 

t

(4.34)

Рассмотрим случай неустановившегося безвихревого течения, когда

С.М. Мухаметшин. Лекции по гидрогазодинамике

~ 10 ~