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

Книги / Теоретические основы гидравлики. Ртищева А.С. 2007 г

.pdf
Скачиваний:
265
Добавлен:
12.06.2014
Размер:
1.62 Mб
Скачать

30

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

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

Напряжение называется нормальным σ, если оно направлено по нормали к анализируемой поверхности.

Напряжение называется касательным τ, если оно направлено по касательной к анализируемой поверхности.

Нормальное напряжение вдоль оси x (действующее на площадку, перпендикулярную оси x) обозначают σx. Аналогично σy, σz.

Касательные напряжения записываются с двумя индексами, например, τxy, первый индекс указывает, какой оси перпендикулярна выделенная площадка, а второй – в направлении какой координатной оси действует рассматриваемое касательное напряжение.

Таким образом, в направлении каждой координатной оси действуют нормальная и две касательных составляющих напряжения. Так в направлении осей x, y, z действуют следующие составляющие напряжения соответственно:

 

σx

τ yx

τzx

.

(3.40)

 

τxy

σ y

τzy

 

 

 

τxz

τ yz

σ z

 

 

Выражение (3.40) представляет собой тензор напряжений.

(3.41)

τxy =τ yx ; τxz =τzx ; τzy =τ yz .

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

диагонали

 

 

 

 

 

σx

τxy

τxz

.

(3.42)

 

τxy

σ y

τ yz

 

 

 

τxz

τ yz

σz

 

 

В невязкой жидкости касательные напряжения равны нулю

(3.43)

τxy =τ yz =τxz = 0.

Также уравнение (3.43) справедливо для вязкой, но покоящейся жидкости.

Для этих случаев также верно

 

 

 

 

σ x =σ y =σz = −p ,

(3.44)

где p – давление, действующее на выделенный элемент противоположно по знаку напряжениям σ, поскольку давление создает напряжение сжатия, а

положительными принято считать растягивающие напряжения.

 

Из формулы (3.44) следует, что

(3.45)

p = −

1

(σ x +σ y +σ z ).

 

 

3

 

Таким образом, в вязкой покоящейся или невязкой (покоящейся или движущейся) жидкости силы трения не возникают, напряжения определяются

31

гидростатическим давлением. Для этих случаев тензор напряжений принимает вид

 

p

0

0

 

 

 

 

(3.46)

 

 

 

 

0

p

0

 

 

 

.

 

0

0

p

 

 

 

 

 

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

зависимостями:

= −p + ∆py ; σ z = −p + ∆pz ;

(3.47)

σ x = −p + ∆px ; σ y

где px, py, pz – дополнительные давления в направлении координатных осей x, y, z, обусловленное влиянием вязкости.

Таким образом, закон трения Ньютона (1.4) для плоского потока можно обобщить на пространственное течение.

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

 

2

r

w

x

 

 

w

 

 

 

 

wy

 

 

σx = −p

 

µdivW + 2µ

 

 

 

 

x

+

 

 

 

 

 

;

 

3

x

 

; τxy = 2µεxy = µ

y

x

 

 

 

 

 

 

 

 

 

 

 

 

 

2

r

wy

 

 

w

x

 

 

 

w

z

 

 

 

 

(3.48)

σ y = −p

 

µdivW + 2µ

 

 

 

; τxz = 2µεxz = µ

 

 

+

 

 

;

3

y

 

z

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

2

r

w

z

 

wy

 

 

 

w

 

 

 

 

 

 

σz = −p

 

µdivW + 2µ

 

 

 

 

 

 

 

 

 

 

z

 

 

 

3

z

 

; τ yz = 2µεyz = µ

z

+

y

.

 

 

 

 

 

 

 

 

 

 

 

r

wx

 

wy

 

wz

 

divW =

 

+

 

 

+

 

.

(3.49)

x

 

y

z

Уравнения (3.48) носят название закон трения Стокса.

 

Для несжимаемой жидкости

divW = 0 , поэтому выражения

(3.48)

упрощаются.

3.2.2. Дифференциальное уравнение неразрывности

Выделим в движущемся потоке элементарный объем в форме прямоугольного параллелепипеда (рис. 3.6) и определим изменение массы dm жидкости в выделенном объеме за элементарный промежуток времени dτ. Это изменение массы определяется разностью между втекающей и вытекающей массой жидкости через грани элементарного объема. Поскольку объем выделенного элемента V = dxdydz остается неизменным с течением времени, то

изменение массы жидкости dm может быть обусловлено лишь изменением ее

плотности dρ

 

d(ρV )

 

 

 

 

dm

 

dρ

dρ

 

 

dτ

=

dτ

=V dτ

= dτ

dxdydz .

(3.50)

32

Рис. 3.6. Выделение элементарного объема в движущемся потоке

Определим массу жидкости, втекающую в выделенный объем за единицу времени.

Жидкость втекает через грани ABFE, AEHD и EFGH в следующих количествах: через грань ABFE ρwx dydz ; через грань AEHD ρwy dxdz ; через

грань EFGH ρwz dxdy .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Жидкость

вытекает

через

грани DCGF,

 

BCGF и ABCD в следующих

количествах: через грань DCGH

 

ρwx

+

(ρw

x

)

 

 

 

 

 

через грань BCGF

 

 

x

 

dx dydz ;

 

 

(ρwy )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ρw

z

 

 

 

ρwy +

 

 

dy dxdz

; через грань ABCD

ρwz

+

 

 

 

dz dxdy .

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Просуммировав количества втекающей и вытекающей жидкости по всем

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

 

 

 

dρ

 

 

 

 

(ρw

x

)

 

 

(ρwy )

 

(ρw

z

)

 

 

 

 

 

 

 

 

 

dxdydz = −

 

 

 

 

 

+

 

 

 

 

 

 

 

 

+

 

 

 

 

 

dxdydz ;

 

(3.51)

 

 

dτ

x

 

 

 

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dρ

(ρw

)

 

 

 

 

(ρwy )

 

(ρw

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dτ +

 

x x

 

+

 

 

 

 

 

+

 

 

z z

 

= 0 .

 

 

 

 

(3.52)

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

Уравнение

(3.52)

 

называется

 

 

 

дифференциальным

уравнением

неразрывности или сплошности. Оно также может быть записано в виде

 

 

 

 

 

 

 

 

dρ

+ div(ρWr)= 0.

 

 

 

 

 

 

 

 

 

 

 

(3.53)

 

 

 

 

 

 

 

 

dτ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для потоков несжимаемой жидкости (стационарных и нестационарных)

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wx

+

 

 

wy

+

wz

 

= 0 ;

 

 

 

 

 

 

 

 

 

 

(3.54)

 

 

 

 

 

 

 

 

x

 

 

 

 

dy

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

divW = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.55)

33

Для стационарных потоков газа (сжимаемой жидкости) уравнение

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

 

 

(ρwy )

 

 

 

 

 

 

(ρw

)

 

 

 

(ρw

)

 

 

 

x x

 

+

 

 

+

z z

 

= 0

;

(3.56)

 

 

y

 

 

 

 

 

div(ρW )= 0 .

 

 

 

(3.57)

Необходимо отметить, что существуют и другие способы вывода уравнения неразрывности.

3.2.3.Дифференциальные уравнения переноса количества движения. Уравнения Эйлера и Навье-Стокса

Рассмотрим движение сплошной среды, предполагая скорость, плотность, напряжения и массовые силы непрерывными функциями времени и координат. В декартовой системе координат выделим элемент в виде прямоугольного параллелепипеда (рис. 3.6). К выделенному объему применяется закон сохранения количества движения, в соответствии с которым, изменение за определенный промежуток времени количества движения жидкости в элементарном объеме, равно импульсу внешних сил, действующих на этот объем.

Изменение количества движения в выделенном объеме происходит за счет изменения плотности жидкости и ее скорости, так и за счет разницы между втекающим и вытекающим количеством движения через границы этого объема.

Уравнение движения удобно рассматривать в проекциях на координатные оси x, y, z. Для получения проекций уравнения движения на каждую координатную ось применим закон сохранения количества движения к граням выделенного объема, перпендикулярным соответствующей оси.

Рассмотрим сначала грань ABFE, перпендикулярную оси x, и определим

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

нее за время dτ

в

направлении координатных осей x, y, z:

 

 

I xx (x)= ρwx wx dydzdτ =ixxdydzdτ ;

 

 

I xy (x)= ρwx wy dydzdτ = ixy dydzdτ ;

(3.58)

I xz (x)= ρwx wz dydzdτ = ixz dydzdτ .

где ixx , ixy , ixz составляющие потока импульса через грань, имеющую

площадь, равную 1 и перпендикулярную оси x, за единицу времени в направлении координатных осей x, y, z соответственно.

34

Потоки количества движения через грань DCGH можно выразить

зависимостями:

 

(ρw

w

 

)

 

 

i

 

 

 

I xx (x + dx)

 

ρwx wx +

 

 

 

 

 

 

=

x

 

x

 

 

dx dydzdτ = ixx +

 

xx

dx dydzdτ ;

 

 

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

(ρwx wy )

 

 

ixy

 

 

I xy (x + dx)= ρwx wy +

 

 

 

 

 

 

 

 

 

 

(3.59)

 

 

 

 

 

 

 

 

x

 

 

 

 

dx dydzdτ = ixy +

x

dx dydzdτ ;

 

 

 

 

 

 

 

 

 

 

 

I xz (x + dx)

 

ρwx wz +

(ρw

w

 

)

 

 

i

 

 

 

=

x

 

z

 

 

dx dydzdτ = ixz +

 

xz

dx dydzdτ .

 

 

 

 

x

 

 

 

 

 

 

x

 

 

Разница

 

между

втекающим

(положительным) и вытекающим

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

I xx = Ixx (x)Ixx (x + dx)= −

(ρwx wx )dxdydzdτ ;

 

 

x

 

I xy = I xy (x)I xy (x + dx)= −

(ρwx wy )

dxdydzdτ ;

(3.60)

 

 

x

 

I xz = Ixz (x)I xz (x + dx)= −

(ρwx wz )dxdydzdτ .

 

 

x

 

Аналогично определяются изменения составляющих количества движения через грани перпендикулярные осям y и z:

I yx = I yx (y)I yx (y + dy)= −(ρwy wx )dxdydzdτ ;

y

I yy = I yy (y)I yy (y + dy)= −(ρwy wy )dxdydzdτ ;

y

I yz = I yz (y)I yz (y + dy)= −(ρwy wz )dxdydzdτ .

y

I zx = I zx (z)I zx (z + dz)= −(ρwz wx )dxdydzdτ ;

z

I zy = I zy (z)I zy (z + dz)= −(ρwz wy )dxdydzdτ ;

z

I zz = I zz (z)I zz (z + dz)= −(ρwz wz )dxdydzdτ .

z

(3.61)

(3.62)

35

Изменение за время dτ количеств движения в выделенном объеме I x , I y , I z в направлении осей x, y, z, происходящее за счет изменения плотности жидкости и ее скорости, выразим соотношениями вида

I x = (ρτwx )dxdydzdτ ;

I y =

(ρwy )

dxdydzdτ ;

(3.63)

τ

 

 

 

 

Iz =

(ρwz )dxdydzdτ .

 

 

τ

 

 

На выделенный объем, в общем случае, действуют внешние массовые (гравитационные или инерционные) и поверхностные силы. Найдем проекции импульса массовых сил на оси x, y, z за время dτ .

Для внешней массовой силы R = i X + jY + kZ , отнесенной к единице объема, имеем

IRx = Xdxdydzdτ ;

IRy =Ydxdydzdτ ;

IRz = Zdxdydzdτ .

(3.64)

Проанализируем поверхностные силы. На грань ABFE, перпендикулярную

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

Px (x), а на противоположную

грань DCGH (и в противоположном направлении) – сила Px (x + dx).

 

r

r

 

 

r

r

 

 

pr

 

 

P

(x)= p

 

dydz ;

P

(x + dx)= p

 

+

x

dx dydz .

(3.65)

 

 

x

x

 

x

 

x

 

x

 

 

Вуравнении (3.65) prx вектор результирующего напряжения,

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

Вектор импульса результирующей поверхностной силы,rдействующей на перпендикулярные оси x грани выделенного объема, I Px определяется

выражением

px

 

 

IrPx = −

dxdydzdτ .

(3.66)

 

 

x

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

Аналогично получаются выражения для импульсов результирующей поверхностной силы, действующей на грани, перпендикулярные осям y и z:

r

 

py

 

 

I Py

= −

 

 

dxdydzdτ .

(3.67)

y

 

 

 

IrPz

= −

pz

 

dxdydzdτ .

(3.68)

z

 

 

 

36

Векторы напряжений px , py , pz , действующих на перпендикулярные осям x, y, z соответственно грани, можно выразить через составляющие тензора

напряжений

r

 

 

r

σ x + jτxy + kτxz ;

 

px = i

 

r

r

τxy + jσ y + kτ yz ;

 

py

= i

(3.69)

r

r

τxz + jτ yz + kσ z .

 

pz

= i

 

Суммарный импульс поверхностных сил

IrP = − irr σxxk τxzx

 

τxy

 

 

τ

 

 

 

r ∂τxy

+

 

 

 

+

 

xz

+

 

 

 

 

 

 

 

 

y

z

 

j

x

 

 

 

 

 

 

 

 

 

τ

yz

 

σ

 

 

 

 

+

 

 

+

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

σ

y

 

τ

yz

 

 

 

 

 

+

 

+

 

 

+

 

 

 

 

 

 

 

 

 

 

y

z

 

 

 

 

 

 

 

 

× dxdydzdτ .

(3.70)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, для различных составляющих закона сохранения

количества движения, получаем

 

 

 

 

 

 

 

 

τxy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(I xx + ∆I yx + ∆I zx )= I Rx

σ

 

x

 

 

 

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xz

 

 

 

 

 

 

 

 

 

 

I x

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

 

dxdydzdτ ;

 

 

 

x

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(I xy + ∆I yy + ∆I zy )= I Ry

 

τ

xy

 

 

 

 

 

σ

y

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yz

 

 

 

 

 

 

 

 

 

 

I y

 

 

 

 

+

 

 

 

 

+

 

 

 

 

 

 

dxdydzdτ ;

 

 

 

x

 

 

y

 

 

z

 

 

 

(3.71)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(I xz + ∆I yz + ∆I zz )= I Rz

 

τ

xz

 

 

 

 

 

τ yz

 

 

 

 

 

 

σ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

I z

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

 

dxdydzdτ .

 

 

 

 

x

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Преобразуя эти уравнения получаем

 

 

уравнения движения в напряжениях:

(ρw

 

 

)

 

 

(ρw

w

 

 

)

 

(ρwx wy )

 

(ρw

 

w

 

 

)

 

 

 

 

 

 

 

 

σ

 

 

 

 

 

τxy

 

τ

 

 

 

 

 

 

x

 

 

+

 

x

 

x

 

 

+

 

 

 

 

+

 

x

 

 

z

 

 

 

= X

+

 

 

 

 

x

+

 

 

 

+

 

xz ;

 

 

τ

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

z

 

 

 

(ρwy )

+

 

(ρwx wy )

+

(ρwy wy )

+

(ρwy wz )

 

=Y +

τxy

 

+

σ y

+

τ yz

;

(3.72)

 

τ

 

 

 

 

 

x

 

 

 

 

 

 

 

y

 

 

z

 

 

 

 

 

 

 

 

 

x

 

 

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂(ρw

 

)

 

 

(ρw

 

w

 

)

 

 

(ρwy wz )

 

(ρw

 

w

 

)

 

 

 

 

 

 

 

τ

 

 

 

 

 

 

τ yz

 

σ

 

 

 

 

 

 

 

z

 

 

+

 

x

 

 

z

 

 

+

 

 

 

 

+

 

z

 

 

z

 

 

 

= Z

+

 

 

 

xz

 

+

 

 

 

+

 

 

z .

 

 

τ

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

y

z

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

37

 

Левая часть уравнений (3.72) может быть преобразована как

 

 

(ρw

 

)

 

 

(ρw

w

 

)

 

 

(ρwx wy )

 

 

(ρw

w

)

 

 

w

 

 

ρ

 

w

 

 

 

x

 

 

+

 

 

 

 

x

 

 

x

 

+

 

 

 

 

 

 

 

 

 

+

 

 

 

x

z

 

= ρ

 

x + wx

 

 

 

+ ρwx

 

x +

 

τ

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

(ρwy )

 

z

 

 

 

 

τ

 

 

x

 

 

(ρw

)

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

w

 

 

 

(ρw

)

 

 

 

 

 

+ wx

 

 

 

 

 

x

 

+ ρwy

 

 

 

x

+ wx

 

 

 

 

 

 

 

 

+ ρwz

 

 

x + wx

 

z

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

z

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

.

(3.73)

 

w

 

 

 

 

 

w

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

+ w

 

x

+ w

 

x

 

+ w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ρ

 

 

 

 

x

 

 

y

 

 

z

 

 

 

x +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ

 

 

 

x

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

(ρw

 

)

 

 

(ρwy )

 

 

(ρw

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ wx

 

 

 

 

+

 

 

 

 

x

 

 

+

 

 

 

 

 

 

+

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично преобразовываются и два других уравнения.

В соответствии с уравнением неразрывности (3.52) последняя сумма в выражении (3.73), стоящая в квадратных скобках, равна нулю. Обозначим

Dwx

=

wx

+ w

 

wx

+ w

 

wx

+ w

 

wx

,

(3.74)

 

 

 

 

 

dτ

τ

x x

y y

z z

где Dwdτx субстанциальная производная.

Таким образом, получаем уравнения движения в форме уравнений Навье-

Стокса:

 

 

Dw

x

 

 

 

 

 

 

 

 

 

p

 

 

w

 

 

 

 

2

 

 

 

r

 

 

 

 

 

 

w

 

 

 

wy

 

 

 

 

 

 

 

 

 

= X

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

x

 

 

 

divW

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

 

 

+

 

 

 

 

µ

 

 

 

 

+

 

 

 

 

 

 

 

 

+

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ

 

 

 

 

 

 

 

 

 

 

 

 

dτ

 

 

 

 

 

w

 

 

 

x

 

 

 

 

x

 

 

 

 

x

 

 

 

3

 

 

 

 

 

 

 

y

 

 

y

 

 

 

 

dx

;

 

 

 

 

 

 

 

 

x

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

µ

 

 

 

 

+

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dwy

 

 

 

 

 

 

 

 

 

p

 

w

 

 

 

 

wy

 

 

 

 

 

 

 

wy

 

 

 

 

2

 

 

 

r

 

 

ρ

 

 

 

 

 

 

 

=Y

 

 

 

+

 

 

 

 

 

 

 

 

 

 

x

+

 

 

 

 

 

+

 

 

 

 

 

2

 

 

 

 

 

 

divW

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dτ

 

 

 

y

 

x

µ

y

 

dx

 

 

 

 

µ

 

y

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

(3.75)

 

 

 

 

 

 

wy

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

µ

 

 

 

 

 

 

+

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dw

z

 

 

 

 

 

 

 

 

p

 

w

 

 

 

w

 

 

 

 

 

wy

 

 

 

w

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

= Z

 

 

 

 

 

 

 

 

 

 

 

µ

 

 

 

x

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

+

 

 

 

+

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

dτ

 

 

 

z

 

 

 

 

 

z

 

dx

 

 

 

µ

z

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

2

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

µ

2

 

 

 

 

 

 

 

divW

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

энергию (энергия отнесена к единице массы среды). Изменение за

38

Система уравнений (3.75) описывает движение вязкой сжимаемой

жидкости (газа). Для несжимаемой жидкости divW = 0 , поэтому уравнения Навье-Стокса упрощаются

 

 

 

 

 

Dwx

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

ρ

 

 

= X p

+ µ

wx

+

wx +

 

wx

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dτ

 

 

 

 

 

 

x

 

 

 

 

x

2

 

 

 

y

2

 

 

 

 

 

 

z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dwy

 

 

 

 

 

 

p

 

 

 

2 wy

 

2 wy

 

 

 

 

 

2 wy

 

 

 

 

 

 

ρ

 

 

 

 

 

 

=Y

 

 

+ µ

 

 

 

 

 

 

+

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

dτ

 

 

 

 

 

 

y

 

 

 

 

x

 

 

 

y

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Dwz

 

 

 

 

 

 

p

 

 

 

2

wz

 

 

2

wz

 

 

 

 

 

 

2

wz

 

 

 

 

 

 

 

 

 

 

= Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

 

+ µ

 

 

 

 

2

 

+

 

 

 

 

 

2

+

 

 

 

 

 

2

 

.

 

 

 

 

 

 

dτ

z

 

x

 

 

y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для идеальной среды µ = 0 из уравнений (3.76) получим

ρ

Dw

x

 

 

 

= X

p

; ρ

Dwy

 

=Y

p

;

 

 

ρ

Dw

z

= Z

p

.

 

dτ

 

 

 

 

 

x

 

dτ

y

 

 

 

dτ

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

wx

 

 

 

 

 

 

 

wx

 

 

 

 

 

 

wx

 

 

 

 

 

 

wx

 

 

 

 

1

 

 

 

 

 

 

 

 

 

1 p

 

 

 

 

 

 

 

+ wx

 

 

+ wy

 

 

+ wz

 

 

=

 

 

X

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

y

z

 

ρ

 

 

ρ x

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wy

 

 

 

+ wx

 

wy

 

+ wy

 

wy

 

+ wz

wy

=

1

 

Y

 

1

 

p

;

 

 

 

τ

 

 

 

 

 

x

 

 

y

 

z

ρ

 

 

ρ y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wz

+ wx

wz

 

+ wy

 

wz

 

+ wz

wz

 

=

1

 

 

Z

 

 

1

p .

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

x

 

 

 

 

 

y

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

ρ z

 

 

(3.76)

(3.77)

(3.78)

Эти уравнения называются уравнениями Эйлера и описывают движение идеальной сплошной среды.

3.2.4. Дифференциальное уравнение энергии

Рассмотрим движение сплошной среды, полагая скорость, плотность, напряжения и массовые силы непрерывными функциями времени и координат. В декартовой системе координат выделим элемент в виде прямоугольного параллелепипеда (рис. 3.6). К выделенному объему применяется закон сохранения полной энергии E, включающий внутреннюю u и кинетическую

W 2

2

определенный промежуток времени количества энергии жидкости в элементарном объеме равно работе L внешних сил, действующих на этот объем, подводимому (отводимому) количеству теплоты Qτ и выделившемуся количеству теплоты при наличии внутренних источников тепла QVτ :

ρ

DE dxdydz = L + Q

+ Q

.

 

 

dτ

τ

V τ

 

(3.79)

 

 

 

 

39

 

Работу внешних массовых сил выразим формулой вида

 

Lм = (Xwx +Ywy + Zwz )dxdydzdτ .

(3.80)

Определим работу нормальных и касательных сил, действующих на грани выделенного объема. Действие на грань ABFE нормальных напряжений

обуславливает работу

(x)= −σx wx dydzdτ .

 

 

 

Lσx

 

(3.81)

На грань противоположную грань DCGH:

 

L

(x + dx)=

σ

w

x

+

(σx wx )

dx dydzdτ .

 

 

 

σx

 

x

 

 

x

 

(3.82)

 

 

 

 

 

 

 

Работа сжимающих напряжений считается отрицательной, а растягивающих – положительной.

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

перпендикулярные оси x:

 

(σx wx )

 

 

L

= L

(x)+ L

(x + dx)=

dxdydzdτ .

 

 

 

σx

σx

σx

 

x

(3.83)

 

 

 

 

Аналогично определяется результирующая работа касательных напряжений, действующих на грани, перпендикулярные оси x:

L

=

 

(τxy wy )

dxdydzdτ ;

(3.84)

 

 

τxy

 

 

x

 

 

 

 

 

 

L

=

 

(τxz wz )

dxdydzdτ .

(3.85)

 

 

τxz

 

 

x

 

 

 

 

Аналогичные рассуждения проводятся применительно к граням, перпендикулярным осям y и z.

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

 

(σ

 

w

 

+

τ

 

w

 

+τ

 

w

 

)+

 

 

 

(τ

 

w

 

+σ

 

w

 

+τ

 

w

 

)+

 

 

 

x

x

x

xy

y

xz

z

 

 

 

 

xy

x

y

y

yz

z

 

L

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dxdydzdτ .

(3.86)

п

 

 

(τxz wx

+τ yz wy +σz wz )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Работа L внешних сил определяется как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L = Lм + Lп .

 

 

 

 

 

 

 

 

 

(3.87)

 

Количество теплоты, подведенное через к выделенному объему,

перпендикулярно оси x, через грань ABFE:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qτx (x)

= qx dydzdτ ,

 

 

 

 

 

 

 

 

(3.88)

где qx плотность теплового потока, направленная вдоль оси x, Вт/м2.

 

 

Количество теплоты, отведенное через грань DCGH:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

(x + dx)= q

 

+

 

 

dx dydzdτ .

 

 

 

 

(3.89)

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τx

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

Соседние файлы в папке Книги