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

Свойства рабочих жидкостей. Гидростатика. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
2 Мб
Скачать
☆
3.6. НОРМАЛЬНЫЕ И КАСАТЕЛЬНЫЕ
х
х
х
X
j

р
р
НАПРЯЖЕНИЯ В ЖИДКОСТИ
Экспериментальный и теоретический анализ полей напряжений и деформаций в жидкости показывает наличие нормальных и касатель­ных напряжений: вдоль оси OY нормальное напряжение
2–(2/3)
p е
yy
ской вязкости, относительная деформация
()/ ()/ ()/
формация ются нормальные напряжения вдоль других осей координат:
 
где
Касательные напряжения
 
Рассматривается динамическое равновесие малого объема ΔxΔyΔz, в котором действуют массовые силы
тельные напряжения. Вдоль оси OY
е V х VyVz    
XX ZZ
()/
V
XX
XY X Y
3.7. СИСТЕМА УРАВНЕНИЙ НАВЬЕ–СТОКСА
XYZ
p е p е 
(( )/ ( )/ )
VyV
, где p – давление, μ – коэффициент динамиче-
()/
Vy
yY
. Аналогично формиру-
, объемная де-
2–(2/3), 2–(2/3)
()/
,
Vz
ZZ
(( )/ ( )/ ).
YZ Y Z
.
 
,
XZ X Z
VzVy 
,,
(( )/ ( )/ )
VzV
jj
, нормальные и каса-
YZ
,
,
22
/(1/)(/)(/)()/
dV dt j р yVх
(1 / 3)( / ) ( ) / ( ) / ( ) /
Аналогичные соотношения получены для других направлений. Оконча­тельно система уравнений Навье–Стокса имеет следующий вид:
/ (1/ )( / ) ( )( ) (1/ 3)( ) ( ) / ;
dV dt j
YY Y Y
/ (1/ )( / ) ( )( ) (1/ 3)( ) ( ) / ;
dV dt j
XX X X
 
YY Y
2222

()/ ()/
Vy Vz
    
YY
2222
VyхVy Vyz   
XYZ
yV divVy 
х VdivVх  
41
.
/ (1/ )( / ) ( )( ) (1/ 3)( ) ( ) / ,
р
х y
р
р
р
р
р
р
X
j
Н
dV dt j
ZZ Z Z
где оператор Лапласа
zV divVz 
222222
( )/ ( )/ (...)/
z  
.
Плюс шесть граничных условий и три начальных, а также уравнение неразрывности – задача замкнута. На практике указанная задача Коши решается для слабо сжимаемых жидкостей (например, вода). Конечно­разностные методы показали свою эффективность применительно к си­стеме уравнений Навье–Стокса.
3.8. ЧАСТНЫЕ СЛУЧАИ СИСТЕМЫ УРАВНЕНИЙ НАВЬЕ–СТОКСА
Частным случаем системы уравнений Навье–Стокса является си­стема уравнений Эйлера, которая моделирует течение идеальной жид­кости без учета явлений турбулентности и вязкости:
/(1/)(/);
dV dt j
YY
/(1/)(/);
dV dt j
XX
/(1/)(/).
dV dt j
ZZ
y
х
z
Система уравнений Эйлера преобразуется в уравнение Бернулли. Пусть
VdV j
XX X
VdV j
ZZ Z
равен
,0,
jjg
XY Z
2
0,5 (1/ ) 0dV dр gdz  
(1 / ) ( / )
2
dV V dV V dV V dV
2( )
алов, где
/( ) (1/2 )PgН V Н   
dхV dt dy V dt dz V dt
(1 / ) ( / ) ,
XYYZZ
. После суммирования полученных трех уравнений:
()zz

0
2
,,
XYZ
х
z
. Дифференциал квадрата полной скорости
, тогда
VdV j
YY Y
(1 / ) ( /
, полагаем: массовые силы
, после интегрирования полных дифференци-
, получается уравнение Бернулли
.
пол
y
,
42
3.9. ТЕЧЕНИЕ КУЭТТА В ПЛОСКОМ КАНАЛЕ
Р
Р
Р
Р
(С ДВИЖУЩЕЙСЯ ВЕРХНЕЙ СТЕНКОЙ)
Рассматривается одномерное стационарное течение несжимаемой вязкой жидкости в плоском канале высотой Н. Верхняя стенка канала
движется со скоростью
V
. Силой тяжести пренебрегаем. Массовых
0
сил нет.
Рис. 13. Схема течения Куэтта
Давление на входе в канал равно
, на выходе –
1
. Градиент дав-
2
ления не изменяется вдоль канала. (ΔР/Δх) – отрицательный, если
Р
. Скорость потока зависит от одной координаты
21
ничные условия: скорость
шение
(1 / ) ( / ) ( / ) ( ) /
х dV dy 

где
*2
/(2 )( / )РН V Рх
0
(0) 0, ( )VVНV
22

/ ) (1 (1 / )VyН V Р y Н
*
0
. Интегрируется соотно-
0
. Тогда
,
– безразмерный градиент давления. Вари-
()VVy
. Гра-
ант А течения потока (рис. 13) – градиент давления отрицательный,
V
– положительная. Вариант В – градиент давления равен нулю,
0
положительная. Вариант С – градиент давления положительный,
V
–
0
V
–
0
положительная. В последнем варианте наблюдается обратное движение жидкости. Течение Куэтта в плоском узком канале моделирует про­цессы в подшипнике скольжения и является первым этапом в расчетах сил трения в соприкасающихся парах при наличии смазки.
43
3.10. ТЕЧЕНИЕ ХАГЕНА–ПУАЗЕЙЛЯ В КРУГЛОЙ ТРУБЕ
Р
Р
Рассматривается стационарное ламинарное течение в круглой трубе вязкой несжимаемой жидкости. Радиус трубы R, плотность жидкости ρ, коэффициент динамической вязкости μ, градиент давления ΔР/Δх – по­стоянный, отрицательный, силой тяжести пренебрегаем, течение осе­симметричное. Система уравнений Навье–Стокса сводится к одному обыкновенному дифференциальному уравнению второго порядка, до­пускающего понижение порядка на единицу:
22
(/) / (1|)/ ,
х dV dr rdV dr
где r – текущий радиус (расстояние от оси трубы), граничные условия: r = 0, dV/dr = 0 (в центре трубы скорость жидкости наибольшая), r = R, V = 0. Понижается порядок уравнения на единицу, полученное соотно-
шение решается методом вариации произвольных постоянных. Выпол­нение первого граничного условия упрощает задачу, так как убирается слагаемое общего решения, содержащее отрицательную степень r. В итоге получается распределение скоростей в следующем виде:
, где
22
SR
сеч
. Средняя скорость потока равна
2
– площадь поперечного сечения трубы,
V
сред
(1 / 4 ) ( / ) ( )V РхRr 
(1 / )SVdS
сеч
2dS rdr – элемент площади, интеграл берется по площади сечения.
Пусть длина трубы равна l, тогда постоянный градиент давления равен
2
//Рх Рl . Средняя скорость в трубе равна
Число Рейнольдса
0,5(64/ Re)( / ) ( )Р ld V
Re /Vd
сред
сред
, где d = 2R, μ = ρν, тогда
2
. Полученное соотношение называется
(/8)(/)VRРl
сред
.
формулой Дарси–Вейсбаха, λ = 64/Re – коэффициент Дарси (коэффици­ент линейного сопротивления). Для турбулентных течений применяется
( 0,25)
коэффициент Блазиуса
/( )h
равны
лин
g

.
0,3164(Re)
. Линейные потери напора
44
3.11. ГИДРАВЛИЧЕСКИЙ УКЛОН.
ДИНАМИЧЕСКАЯ СКОРОСТЬ ПОТОКА
Рассматривается круглая труба длиной l, радиусом R, градиент дав­ления отрицательный, течение происходит с постоянной скоростью. Следовательно, разность давлений уравновешивается трением жидко-
сти о внутреннюю поверхность трубы:
()(2)Р RRl 
, где τ – каса-
2
тельные напряжения, действующие на внутренней поверхности трубы. Гидравлический радиус трубы равен
(/)Р l
трубопровода:
V
пограничного слоя мала и равна
. Гидравлический уклон – потери напора на единицу длины
/hl
, тогда
лин
называется динамической скоростью потока. Пусть толщина
динам
погр
гидр
gR
гидр
(пристенный слой, в котором ско-
рость жидкости возрастает от нулевой до скорости теза Ньютона вязкого трения имеет следующий вид:
. Тогда
0,5RR
. Величина
V
пленки
(/ )R
гидр
гидр
0,5
()gR 
), тогда гипо-
(/).V
пленки погр
Из экспериментальных исследований известно, что отношение
(/)( /)constVV V N
пленки динам динам погр
– постоянная, например,
для некоторых жидкостей N = 11,6. Оценка толщины пограничного слоя
погр динам
/NV
. Следовательно
пропорционально (1/Re), т. е.
погр
чем больше скорость потока, тем меньше толщина пограничного слоя. Указанная зависимость хорошо согласуется с результатами эксперимен­тальных исследований.
3.12. ТУРБУЛЕНТНАЯ ВЯЗКОСТЬ. ФОРМУЛА ПРАНДТЛЯ
Рассматриваются две малые равные коллинеарные площадки ΔS в
V
турбулентном потоке (рис. 14). Скорость нижней площадки
V
ней –
. Поперечная скорость между слоями равна U – скорость пере-
2
, верх-
1
мещения. Расстояние между площадками Δy. Плотность жидкости равна ρ.
Перемещение жидкости между слоями оценивается с помощью со­хранения количества движения. Масса жидкости, ушедшая из верхнего слоя в нижний за время Δt, равна Δm = ρUΔSΔt. Изменение количества
45
Рис. 14. Схема оценки турбулентной вязкости
y
y
L
y
L
движения слоев сти за время Δt: τΔSΔt. Тогда равенство
( ).mV V
21
Импульс сил трения между слоями жидко-
()UV V
преобразуется к
21
следующему виду с учетом экспериментальных данных: скорость пере­мещения
()UVVV 
21
:
()( / )

равна квадрату глубины перемещения
22
(/)
следующий вид:

V
. Объединение полученного соотноше-
22
V
2
L
. Формула Прандтля примет
. Величина
(y)
2
ния с касательными напряжениями при ламинарном течении жидкости записывается как
ной вязкости равен
(/)( )Vy  
2
турб
(/ )
, где коэффициент турбулент-
турб
Vy
. Формула Прандтля имеет
важное значение в области моделирования турбулентных потоков, рас­четов сил трения, величин потерь давления.
3.13. МЕСТНЫЕ СОПРОТИВЛЕНИЯ. ФОРМУЛА БОРДА
Рассматривается течение жидкости плотностью ρ при внезапном увеличении диаметра трубопровода (рис. 15). Для сечений 1-1, 2-2 при­меняется соотношение сохранения количества движения относительно сечения 2-2.
46
рости
,РР
– величины давлений;
12
,VV
в сечениях 1-1, 2-2 соответственно. Импульс сил давления
12
,DD
– диаметры трубопровода, ско-
12
за время Δt равен изменению количества движения жидкости:
( ) 0, 25 0, 25 ( ) .РР DtV DVVt
12 2 2 221

Преобразуем полученное соотношение
22
//РРVVV 
12212
2
к виду
уравнения Бернулли:
/(1/2) /(1/2)(1/2)( )Р ggVР ggV gVV  
112 212
где потери на трение
222
трения 1 2
(1 / 2 )( )hgVV
2
– один из вариантов фор-
,
мулы Борда. Объемный расход жидкости не меняется вдоль трубопро­вода:
0,25 0,25 VDV D
112 2
22
трения 1 1 2
, тогда
22
(1/2 ) 1 ( / )hgVDD

2
.
Рис. 15. Схема внезапного расширения трубопровода
2
Величина

1( / )DD 
2
12
называется коэффициентом местного
сопротивления. В гидравлике любое изменение геометрических пара­метров трубопровода рассматривается как местное сопротивление. Как правило, величины местных сопротивлений определяются опытным пу­тем. Конкретные данные приведены в специальной литературе.
47
3.14. ГИДРАВЛИЧЕСКИЙ РАСЧЕТ
g
g
ТРУБОПРОВОДОВ
3.14.1. ПРЕДВАРИТЕЛЬНЫЕ ПОНЯТИЯ И ОПРЕДЕЛЕНИЯ
В предлагаемом случае расчета трубопроводных систем необходимо учитывать потери на трение и местные сопротивления. Потери напора потока на трение определяется по формуле Дарси–Вейсбаха:
2
– коэффициент гидравлического сопротивления; l – длина трубы;
где
d – диаметр трубы;
– скорость потока жидкости; g – ускорение сво-
w
h
тр

lw
, (3.1)
2
dg
бодного падения. Для расчета местных гидравлических сопротивлений (сужение и расширение, поворот потока, краны, задвижки, клапаны, шайбы и другие местные препятствия течению) используют модифици­рованную формулу Дарси–Вейсбаха:
2
где
– коэффициент местного гидравлического сопротивления.
w
 , (3.2)
h
м
2
Таким образом, суммарные потери равны
Пусть
hh h
тр м
 – объемный вес жидкости;
. (3.3)
,PP
– давления в сече-
12
ниях 1, 2 (рис. 16).
Запишем уравнение Бернулли для сечений 1 и 2:
22
pw p w
11 2 2
zz h
Пьезометрическую высоту
если он неизвестен, и располагаемым напором
  
12
22
gg

р
1
назовем потребным напором
γ
Н
расп
48
. (3.4)
Н
потр
, если он задан.
,
Из уравнения (3.4.) можно найти
p
K
H
K
нимается жидкость. Пусть
р
2
– пьезометрическая высота в конце трубо-
γ
()zz
21
– высоту, на которую под-
провода (давление, которое должен преодолеть поток жидкости на вы­ходе из трубопровода). Сумма этих высот называется статическим напо­ром и представляет собой эквивалентную высоту подъема жидкости:
Рис. 16. Схема простого трубопровода
Сумму потерь
h
можно представить как степенную функцию объ-
емного расхода жидкости
()
Hzz
ст 2 1
m
Q
. Тогда уравнение Бернулли принимает
2
, (3.5)
вид так называемого напорно-расходного уравнения:
m
потр(расп) ст
HKQ
,
(3.6)
где
– сопротивление трубопровода и m – показатель степени (оба па-
раметра зависят от режима движения жидкости в трубопроводе).
Для турбулентного режима

Km

 
l
d
80,316
,,2.
24 4
gd
49
Re
(3.7)
Рассматривается гидравлический расчет разветвленного трубопро-
p
p
Н
вода. Система напорно-расходных уравнений в случае разветвления ос­новного трубопровода на две ветви приведена ниже (3.8):
m
,
m
,
(3.8)
,
Hz
ст1 1
1
, (3.9)
где
коэффициенты
 ,
число Рейнольдса
HH KQ
 
HH KQ
 
QQ Q
 
Hz
ст2 2
,KK
вычисляются по соотношениям (3.7). Вначале
12
Rewd
определяется по данному расходу Q (рис. 17).

нст111

нст222

12
2
,
Можно сразу принять коэффициент
= 0,02 для нулевого приближе-
ния, но это может привести к увеличению числа итераций.
Соотношения (3.8) с учетом (3.9) являются нелинейными относи-
Q
Q
,
тельно неизвестных
. В данном случае систему (3.8) можно пре-
1
2
образовать к одному алгебраическому уравнению второго порядка:
исключить
выразить
относительно
первом трубопроводе находится число Рейнольдса
2320 ламинарный

4000 турбулентный

равлического сопротивления
после вычитания второго уравнения из первого;
м
QQQ
21
Q
. После вычисления скорости потока жидкости в
1
. Далее решается квадратное уравнение
Re

 
. Полученное значение коэффициента гид-
(3.7) сравнивается с найденным ранее. За-
wd
тем вычисляется величина полного расхода в точке ветвления трубопро­вода (3.6). Если относительная погрешность полученных значений коэф­фициентов гидравлического сопротивления превосходит 3 %, тогда реше­ние задачи продолжается методом последовательных приближений. В рам­ках следующего шага необходимо найти относительную погрешность
50
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]