МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ Тольяттинский государственный университет Инженерно-строительный институт

advertisement
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ
Тольяттинский государственный университет
Инженерно-строительный институт
Кафедра «Водоснабжение и водоотведение»
КУРС ЛЕКЦИЙ
по дисциплине «Гидравлика»
для студентов
строительных специальностей
очной формы обучения
(технология 30/70)
Преподаватель: доцент Калинин А.В.
Тольятти 2007
УДК 532.5
Курс лекций по дисциплине «Гидравлика» для студентов строительных специальностей очной формы обучения (технология 30/70) / Составитель Калинин А.В., Лушкин И.А.
– Тольятти: ТГУ, 2007.
Рассмотрены основные законы покоя и движения жидкостей и газов, гидравлические сопротивления, основы моделирования гидромеханических явлений. Приведены практические задачи
и примеры их решения. Представлены вопросы для самостоятельной подготовки.
Ил. 116; табл. 22; библиогр.: 10 наим.
Составители: Калинин А.В., Лушкин И.А.
Научный редактор: Вдовин Ю.И.
© Тольяттинский государственный университет, 2007
2
Оглавление
Оглавление.........................................................................................................................................................................3
Требования к обязательному минимуму содержания основной образовательной программы по направлению
подготовки дипломированного специалиста «Строительство»....................................................................................5
Лекция 1. Введение в предмет «Гидравлика». Основные свойства жидкостей и газов .............................................6
1.1. Предмет гидравлики.............................................................................................................................................8
1.2. История предмета.................................................................................................................................................8
1.3. Капельные и некапельные жидкости..................................................................................................................9
1.4. Силы, действующие в жидкости .........................................................................................................................9
1.5. Давление и его свойства ....................................................................................................................................10
1.6. Основные физические свойства жидкостей .....................................................................................................12
1.7. Вязкость. Идеальная жидкость..........................................................................................................................15
Лекция 2. Основы гидростатики, динамики и кинематики жидкости .......................................................................17
2.1. Тема 1. Равновесие жидкости............................................................................................................................17
2.1.1. Дифференциальное уравнение равновесия жидкости. Поверхность равного давления .....................20
2.1.2. Основное уравнение гидростатики ..........................................................................................................21
2.1.3. Закон Паскаля.............................................................................................................................................23
2.1.4. Абсолютное, манометрическое и вакуумметрическое давление ..........................................................23
2.1.5. Сила давления на плоские и криволинейные поверхности....................................................................24
2.1.6. Относительный покой жидкости..............................................................................................................29
2.1.7. Закон Архимеда .........................................................................................................................................30
2.1.8. Основное уравнение гидростатики для сжимаемой жидкости..............................................................31
2.1.9. Изотермическая атмосфера.......................................................................................................................32
2.1.10. Неизотермическая атмосфера .................................................................................................................32
2.2. Тема 2. Основы кинематики и динамики жидкости и газа .............................................................................34
2.2.1. Основные понятия кинематики жидкости...............................................................................................37
2.2.2. Уравнение неразрывности ........................................................................................................................39
2.2.3. Виды движения жидкости.........................................................................................................................40
2.2.4. Интегральная формула количества движения.........................................................................................41
2.2.5. Дифференциальное уравнение движения невязкой жидкости (уравнение Эйлера)............................42
2.2.6. Общее уравнение энергии в интегральной форме (Уравнение Бернулли для потока реальной
жидкости) .............................................................................................................................................................43
2.2.7. Три формы представления уравнения Бернулли для потока реальной жидкости ...............................45
2.2.9. Особенности турбулентного и ламинарного течения жидкости. Число Рейнольдса ..........................47
2.2.10. Уравнение Бернулли для элементарной струйки невязкой сжимаемой жидкости............................49
2.2.11. Уравнение Бернулли для потока вязкой сжимаемой жидкости ..........................................................50
Лекция 3. Основы моделирования гидромеханических процессов............................................................................52
3.1. Основы моделирования .....................................................................................................................................53
3.2. Виды подобия. Масштабы моделирования ......................................................................................................54
3.3. Критерии подобия ..............................................................................................................................................55
3.4. Конечно-разностная форма уравнения Навье-Стокса.....................................................................................56
3.5. Общая схема применения численных методов и их реализация на ЭВМ.....................................................58
3.6. Измерительные приборы, используемые при проведении экспериментальных работ................................60
Лекция 4. Гидравлические сопротивления ...................................................................................................................65
4.1. Виды гидравлических сопротивлений..............................................................................................................69
4.2 Сопротивление по длине при движении в цилиндрической трубе при ламинарном течении .....................70
4.3. Формула Дарси-Вейсбаха ..................................................................................................................................73
4.4. Турбулентное движение в гидравлически гладких и шероховатых трубах..................................................73
4.5. Движение жидкости в трубах некруглого сечения..........................................................................................76
4.6. Местные гидравлические сопротивления ........................................................................................................76
3
4.7. Зависимость коэффициентов местных сопротивлений от числа Рейнольдса. Эквивалентная длина ........80
4.8. Кавитация............................................................................................................................................................81
4.9. Истечение жидкостей из отверстия в тонкой стенке.......................................................................................81
4.10. Зависимость коэффициентов истечения от числа Рейнольдса.....................................................................83
4.11. Истечение из насадков .....................................................................................................................................83
4.12. Виды насадков ..................................................................................................................................................84
4.13. Истечение при переменном напоре и под уровень жидкости ......................................................................85
Лекция 5. Практическое применение законов гидравлики .........................................................................................87
5.1. Расчет короткого трубопровода ........................................................................................................................89
5.2. Расчет длинных трубопроводов ........................................................................................................................92
5.2.1. Понятие о простом и сложном напорных трубопроводах .....................................................................92
5.2.2. Расчет трубопроводов, соединенных последовательно и параллельно ................................................93
5.2.3. Гидравлический удар.................................................................................................................................95
5.2.4. Гидравлический таран ...............................................................................................................................96
Рекомендуемая литература ............................................................................................................................................99
4
Требования к обязательному минимуму
содержания основной образовательной
программы по направлению подготовки
дипломированного специалиста
«Строительство»
ОПД.ФС 2.02 – ГИДРАВЛИКА
(90 часов)
Вводные сведения, основные физические свойства жидкостей и газов, основы кинематики,
общие законы и управления статики и динамики жидкостей и газов, силы, действующие в жидкостях, абсолютный и относительный покой (равновесие) жидких сред, модель идеальной (невязкой) жидкости, общая интегральная форма уравнения количества движения и момента количества движения, подобие гидромеханических процессов, общее уравнение энергии в интегральной и дифференциальной формах, турбулентность и ее основные статистические характеристики, конечно-разностные формы уравнений Навье-Стокса и Рейнольдса, общая схема применения численных методов и их реализация на ЭВМ, одномерные потоки жидкостей и газов.
5
Содержание курса лекций «Гидравлика»
Лекция 1.
Введение в предмет «Гидравлика».
Основные свойства жидкостей и газов
Основные понятия: предмет гидравлики; гидромеханическое понятие жидкости; капельные
и некапельные жидкости; силы, действующие в жидкости; давление и его свойства; основные
физические свойства жидкостей и газов; идеальная жидкость.
Вопросы, на которые необходимо найти ответ в ходе изучения темы:
1. Что является объектом изучения в дисциплине «Гидравлика»?
2. Гидравлика, как научная дисциплина.
3. Жидкость – что это?
4. В чем состоит отличие жидкостей от твердых тел и газов?
5. Что в гидравлике является «сплошной средой»?
6. Какими свойствами обладает изотропная жидкость?
7. Какие дисциплины являются базовыми для изучения законов гидравлики?
8. Почему изучение законов движения жидкости считается более сложным, чем изучение
законов движения твердых тел?
9. Почему в гидравлике придается большое значение экспериментальным исследованиям?
10. В каких дисциплинах используются знания, полученные в гидравлике?
11. Какой вклад в понимание законов движения жидкости внесли следующие ученые: Архимед, Леонардо да Винчи, Ньютон, Эйлер, Бернулли, Дарси, Шези, Вейсбах, Прандтль,
Жуковский?
12. Чем капельные жидкости отличаются от некапельных?
13. При каких условиях капельные жидкости становятся некапельными?
14. Какие силы называются массовыми, а какие поверхностными? Почему?
15. От чего зависит величина силы тяжести, силы давления, силы инерции, силы трения?
16. Что называют напряжениями и почему они возникают в жидкости?
17. Какие бывают напряжения в жидкости, в чем их отличие?
18. Как можно измерить давление в точке, находящейся в жидкости?
19. Какими свойствами обладает давление в жидкости? Как это доказать?
20. Как изменение давления dp записать в дифференциальной форме?
21. Как изменение давления dp записать в векторной форме?
22. Как записать вектор перемещения от одной точки к другой?
23. Что называют вектором градиента давления?
24. Какие свойства имеет вектор градиента давления? Как это доказать?
25. Какие существуют единицы давления? Каково их соотношение?
26. Перечислите основные физические свойства жидкостей.
27. Что называют плотностью? Какова связь между плотностью жидкости и удельным весом?
28. Что называют коэффициентом объемного сжатия жидкости, модулем упругости?
29. От чего зависит коэффициент температурного расширения?
30. Как определить плотность жидкости, если ее температура изменилась?
6
31. Что называется вязкостью жидкости?
32. Чем вязкие жидкости отличаются от невязких?
33. Как вязкость влияет на величину силы трения в жидкости?
34. В чем суть гипотезы Ньютона?
35. Что называют градиентом скорости?
36. Как влияет изменение температуры на вязкость жидкостей и газов? Почему?
37. Что показывает эпюра скорости?
38. Какая существует связь между динамической и кинематической вязкостью?
39. Какие существуют единицы измерения вязкости?
40. Как можно определить вязкость жидкости?
41. Что называют текучестью жидкости?
42. Что понимают под идеальной жидкостью?
43. Как определить градусы Энглера?
44. Почему возникает поверхностное натяжение в жидкости?
45. От чего зависит поверхностное натяжение? Почему?
46. Как определить силу поверхностного натяжения?
47. Что показывает краевой угол? От чего зависит его величина?
48. Когда необходимо учитывать величину поверхностного натяжения? Почему?
49. Что влияет на пенообразование в жидкости?
50. От чего зависит способность жидкости растворять газы?
51. Как определить объем газа, который может раствориться в жидкости или выделиться из
нее при изменении давления?
52. Что называют давлением насыщенного пара?
53. Что влияет на изменение агрегатного состояния жидкости?
54. Что называют кипением жидкости?
55. Для чего необходимо знать давление насыщенного пара?
56. Что понимают под идеальным газом?
57. Как зависит вязкость газов от температуры и давления? Почему?
58. В чем суть уравнения состояния идеального газа?
Практические задачи, решение которых может быть найдено после изучения теоретического материала:
Задача 1-1,6.
Резервуар объемом V1 = ___м3 наполнен водой. При увеличении давления на свободной поверхности на величину Δр = 2·106 Па, объем воды уменьшился. Определить объем при увеличении давления.
Задача 1-2,7.
Сколько килограмм мазута необходимо приобрести, чтобы заполнить резервуар объемом V1
=___м3. Максимальная температура хранения мазута t1 = 50°С.
Задача 1-3,8.
Определить наименьший объем расширительного резервуара системы водяного отопления,
чтобы он полностью опорожнялся. Допустимое колебание температуры воды во время перерывов в отоплении Δt° = 25°С. Объем воды V1 =___м3, коэффициент температурного расширения
βt = 0,0006 град-1.
Задача 1-4,9.
7
Давление воды в закрытом сосуде p1 = 5 атм. При повышении температуры давление повысилось на Δр = 0,04 МПа. Изменением плотности и деформацией стенок пренебречь. Определить изменение температуры Δt, если коэффициент температурного расширения βt = 0,2·10–3 °С.
Задача 1-5,10.
При испытании резервуара, давление на начало испытаний было p1 = 5500000 Па, Через некоторое время давление уменьшилось на величину Δр = 35000 Па. Определить объем воздуха,
вышедшего из резервуара через не плотности, если объем резервуара V1 =___м3. Температура
воздуха не изменилась за время испытаний, t = 20°С, а коэффициент объемного сжатия βw =
0,49·10 –9 Па–1.
1.1. Предмет гидравлики
Гидравлика – наука о движении и покое воды и других жидкостей. Жидкостью в гидравлике
представляют как сплошную среду, легко изменяющую форму под действием внешних сил.
Сплошная среда – это масса, физические и механические параметры которой являются функциями координат в выбранной системе отсчета. Молекулярное строение жидкостей заменяется
сплошной средой той же массы.
В данном курсе гидравлики мы считаем, что жидкость имеет одинаковые свойства по всем
направлениям, то есть является изотропной. Законы движения и покоя жидкостей основываются на законах механики сплошной среды и физики.
Вместе с тем, можно сказать, что явления, изучаемые в гидравлике сложнее явлений, которые являются объектом исследования в механике твердого тела, т.к. жидкости легко изменяют
свою форму под действием небольших внешних сил, сжимаемые жидкости (газы) изменяют и
свой объем. Эти свойства связаны с молекулярной структурой строения жидкости. Из-за того,
что движение жидкости очень часто не удается точно математически описать, в гидравлике
приходиться использовать упрощенные математические модели, которые затем уточняются и
дополняются в ходе экспериментального исследования.
Знание законов гидравлики необходимо для решения практических задач теплогазоснабжения: расчета систем водоснабжения, тепловых сетей, теплообменных аппаратов, насосов и т.д.
1.2. История предмета
250 лет до н.э. Архимед установил принципы гидростатики.
XV в. – Леонардо да Винчи положил начало экспериментальной гидравлике (движение воды в каналах, через отверстия).
XVII в. – Торичелли предложил формулу для определения скорости жидкости, вытекающей
из отверстия.
Ньютон высказал основные положения о внутреннем трении в движущихся жидкостях и газах.
XVIII в. – Даниил Бернулли и Леонард Эйлер разработали общие уравнения движения
идеальной жидкости (газа).
XIX в. – Шези, Дарси, Базен, Вейсбах начали опытное изучение движения воды в различных частных случаях и получили множество эмпирических формул. Полученные выводы расходились с теоретическими выводами, полученными ранее.
В конце XIX в. – Петров исследовал трение при ламинарном режиме.
Рейнольдс изучил переход от ламинарного режима к турбулентному, начал изучать гидравлические сопротивления.
XX в. – Жуковский и Прандтль положили начало изучению турбулентных потоков, неустановившегося движения жидкости.
В связи с развитием авиации, космонавтики, теплоэнергетики, машиностроения, автомобилестроения происходит бурное развитие науки как экспериментальное (лазеры, датчики), так и
теоретическое (использование современного математического аппарата, ЭВМ).
8
1.3. Капельные и некапельные жидкости
В гидравлике, в основном, считается, что жидкость практически не изменяет свой объем под
действием внешних сил, т.е. является несжимаемой. К несжимаемым жидкостям относятся все
капельные жидкости: вода, нефть, мазут. В отличие от капельных жидкостей, газы (воздух,
пропан, бутан и т.д.) легко изменяют объем под действием внешних сил, сжимаются, поэтому
их называют сжимаемыми. Любая капельная жидкость может переходить в газообразное состояние при определенной температуре и давлении. Соответственно, газы при понижении температуры и повышении давления могут переходить в жидкое состояние.
1.4. Силы, действующие в жидкости
На произвольно выделенный объем жидкости действуют два вида сил:
Поверхностные:
Р – сила давления
Т – сила трения
Массовые:
G – сила тяжести
I – сила инерции
Массовые силы действуют по всему выделенному объему и пропорциональны его массе
G = mg , I = − ma . Поверхностные силы действуют по поверхности и пропорциональны площади поверхности.
Рассмотрим подробно поверхностные силы. Под влиянием внешних сил, действующих на
выделенный объем возникают соответствующие внутренние силы. Проведем внутри объема
поверхность S, разделяющую его на две части I и II (см. рис. 1.1). Отбросим часть II и для сохранения равновесия введем силы такие же, как и силы с которыми часть II действовала на
часть I. На элементарную площадку Δs разделяющей поверхности действует сила Δf. Площадь
Δs может быть стянута в точку М с координатами x, y, z. В этом случае площадь поверхности
Δs, так и сила Δf стремится к нулю. Отношение силы df к площади поверхности ds стремится к
df
пределу σ =
, который называют напряжением.
ds
dfT
I
df
II
ds
dfN
Рис. 1.1. Напряжение в жидкости
Силу df, действующую на площадь ds можно разложить на две составляющие: тангенциальную и нормальную. Соответственно, напряжение в жидкости может быть тангенциальным
(τ) и нормальным (p). Тангенциальное напряжение, действующее вдоль поверхности ds, называют напряжением трения.
df
τ= T ,
(1.1)
ds
где dfT – сила трения площади ds.
Нормальное напряжение, действующее по нормали к поверхности ds, называют напряжением давления или давлением
df
p= N ,
(1.2)
ds
где df N – сила давления площади ds.
9
Для площади S можно записать
τ=
T
,
S
(1.3)
p=
P
,
S
(1.4)
где Т – сила трения площади S.
где Р – сила давления площади S.
В покоящейся жидкости имеется только нормальное напряжение, тангенциальное напряжение отсутствует.
Сила трения действует вдоль поверхности:
T = τS .
(1.5)
Сила давления направлена по нормали к поверхности:
P = pS .
(1.6)
1.5. Давление и его свойства
В любой точке жидкости имеется давление и его можно измерить, опустив в жидкость стеклянную трубочку с запаянным концом из которой выкачен воздух. Рассмотрим точку М в жидкости, проведем через эту точку поверхность ds (рис. 1.2). Результирующая сила воздействия
всех молекул, находящихся в постоянном движении, на эту поверхность перпендикулярна ds.
r
r
r
Можно записать в векторной форме df N = pdsn , где n – единичный вектор, направленный по
нормали к поверхности ds.
ds
r
n
M
r
fN
Рис. 1.2. Давление в точке
r
df N зависит от величины поверхности ds, но из формулы видно, что давление в точке не зависит от ds.
Свойства давления:
1. Давление в точке в любом напрвлении одинаково и не зависит от ориентации ds. Через
точку М можно провести бесконечное множество поверхностей и сила df N будет зависеть
только от величины ds.
2. Гидростатическое давление является непрерывной функцией координат пространства
p = f ( x, y , z ) .
(1.7)
Понятие о градиенте давления. Рассмотрим точку М, имеющую координаты (x, y, z) и находящуюся в жидкости (рис. 1.3). Давление в точке М – pM . Это давление зависит только от координат точки М. Можно записать pМ = р ( x, y, z ) .
На небольшом расстоянии от точки М находится точка М1 с координатами (x+dx, y+dy,
z+dz).
Давление в М1 отличается от давления рМ на некоторую величину dp: pM1 = pM + dp .
Давление pM1 зависит от координат точки М1: рМ1 = р ( x + dx, y + dy, z + dz) .
Тогда dp = р ( x + dx, y + dy, z + dz ) − р ( x, y, z ) .
10
z
M1
r
k
y
r
i
r
j
x
M
Рис. 1.3. Градиент давления
Т.к. р является функцией координат x, y, z, то величину dp можно записать в дифференциальной форме
∂p
∂p
∂p
dp =
dx + dy + dz.
∂x
∂y
∂z
Вектор перемещения от точки М к точке М1 записывается в форме
uuuur r
r
r
MM 1 = idx + jdy + kdz ,
r r r
где i , j , k – единичные векторы, направленные вдоль осей координат.
Определение: В физике для обозначения изменения некоторой скалярной величины G (температуры, давления) от одной точки к другой используется понятие вектора
uuuuur
∂G r ∂G r ∂G r
gradG =
i+
j+
k.
(1.8)
dx
dy
dz
Значит, вектор градиента давления величин
∂p ∂p ∂p
, , .
∂x ∂y dz
Можно записать в виде
uuuuur
∂p r ∂p r ∂p r
grad p = i +
j+ k.
dx dy
dz
Произведение двух векторов
uuuuur uuuur ∂p
∂p
∂p
grad pMM 1 =
dx + dy + dz
dx
dy
dz
или
(1.9)
(1.10)
uuuuur uuuur
dp = grad pMM 1 .
Вывод: изменение давления dp является скалярным произведением двух векторов grad p и
ММ1.
uuuuur
Свойства вектора grad p :
1. Если точки М и М1 принадлежат поверхности в которой все точки испытывают одинаковые давление, то можно записать
рМ = рМ 1 ,
тогда dp = 0.
uuuuur
Вывод: grad p расположен по нормали к поверхности равного давления, проходящей через
точку М.
2. Предположим, что М1 расположена по нормали к поверхности равного давления, прохоuuuuur
дящей через точку М, тогда dp > 0. Значит скалярное произведение grad p MМ1 имеет положиuuuuur
тельное значение и grad p имеет то же направление, что и ММ1.
11
uuuuur
Вывод: grad p направлен в сторону увеличения давления.
uuuuur
dp
.
3. grad p =
MM 1
uuuuur
Вывод: величина grad p определяется отношением разности давлений в двух точках к расстоянию между этими точками.
Единицы давления.
При измерении атмосферного давления используют единицу давления – бар
5
• 1бар = 10 Па
При измерении при помощи пьезометрических трубок используют единицы длины
• для воды 1 мм.в.ст = 9,8 Па
• тогда 1бар = 9,8 м.в.ст
• для ртути 1 мм.рт.ст = 133,3Па
9810
Можно подсчитать, что: 1 бар =
мм.рт.ст.
133,3
Старая система измерения. Давление, создаваемое телом массой 1 кг на 1 см2:
кгс
1 2 = 0,98 бар .
см
1.6. Основные физические свойства жидкостей
Плотность – масса жидкости m, заключенная в единице объема V
m
ρ= .
V
Плотность меняется при изменении температуры и давления.
Удельный вес – вес жидкости G в единице объема V
G
γ= .
V
Плотность и удельный вес связаны между собой соотношением
γ = gρ .
(1.11)
(1.12)
(1.13)
Сжимаемость – свойство жидкостей изменять объем при изменении давления. Сжимаемость капельных жидкостей характеризуется коэффициентом объемного сжатия, который
представляет собой относительное изменение объема жидкости на единицу изменения давления
1 ΔV
βV = − ⋅
,
(1.14)
V Δp
где V – первоначальный объем жидкости;
Знак «минус» в формуле обусловлен тем, что положительному приращению давления р соответствует отрицательное приращение (уменьшение) объема.
Величина, обратная коэффициенту объемного сжатия называется модулем упругости
1
E0 =
.
(1.15)
βV
В качестве примера, приведем значения модуля упругости стали и воды:
Еводы = 2 × 109 Па, Естали = 2 × 1011 Па . Таким образом, мы видим, что упругость воды всего только
в 100 раз меньше упругости стали, значит, воду можно рассматривать как несжимаемое вещество.
12
Коэффициент объемного сжатия и модуль упругости капельных жидкостей практически не
изменяется при изменении давления, и на практике очень часто их считают неизменными. Част
при расчетах коэффициент объемного сжатия βV воды принимают постоянным и равным
0,49×10–9 Па–1.
Сжимаемость характеризуется также отношением изменения давления к изменению плотности, равным квадрату скорости распространения звука в среде:
dp
a2 =
.
(1.16)
dρ
Очевидно, для малосжимаемой среды при больших изменениях давления изменение плотности незначительно и скорость звука получается большой, и наоборот, при большой сжимаемости скорость звука оказывается малой (для воздуха – 330 м/с).
Для оценки сжимаемости среды при ее движении важно не абсолютное значение скорости
звука а, а относительное, которое называется числом Маха:
Ма = u/a.
Если скорость движения воздуха мала по сравнению со скоростью движения звука в ней,
число Маха мало по сравнению с единицей и движущуюся среду можно рассматривать как несжимаемую жидкость. Скорость воздуха в воздуховодах, газа в газопроводах низкого давления
и газоходах котельных установок не превышает 12 м/с. Следовательно, в практике теплоснабжения и вентиляции газ(воздух) можно рассматривать как несжимаемую жидкость. При движении газов со скоростью более 70 м/с влияние сжимаемости следует учитывать.
В отличие от капельных жидкостей газы характеризуются значительной сжимаемостью и
высокими значениями коэффициента температурного расширения. Зависимость плотности газов от давления и температуры устанавливается уравнением состояния. Для совершенных газов (нет взаимодействия между молекулами) справедливо уравнение Клапейрона, позволяющее
определить плотность газа при известных давлении и температуре
p
,
ρ=
(1.17)
RT
где р – абсолютное давление; R – удельная газовая постоянная (для воздуха R= 283 Дж/кг·К);
Т – абсолютная температура.
В технических расчетах плотность газа, приводят к нормальным физическим условиям
(t = 0ºC; p = 101325 Па) или стандартным условиям (t = 20ºC; p = 101325 Па)
Можно подсчитать, что в стандартных условиях плотность воздуха ρ = 1,21 кг/м3.
При других условиях плотность воздуха можно определить по формуле
р Т0
ρ = ρ0
,
(1.18)
р0 Т
где ρ0, Т0 и р0 – плотность, температура и давление при известных стандартных условиях соответственно.
Сжимаемость газа зависит от характера процесса изменения состояния. Для изотермического процесса сжимаемость воздуха составляет примерно 9,8×104 Па, что превышает в 20000 раз
сжимаемость воды.
Температурное расширение – увеличение объема капельных жидкостей, при увеличении
температуры, характеризуется коэффициентом температурного расширения βt , выражающим относительное увеличение объема жидкости при увеличении температуры на 1 град.
1 ΔV
βt = ⋅
,
(1.19)
V Δt
где ΔV – изменение объема при повышении температуры на величину Δt .
13
Если считать, что плотность не меняется при изменении давления, а только от температуры,
то для расчета изменения плотности капельных жидкостей с изменением температуры можно
использовать формулу
1
ρ = ρ0
,
(1.20)
1 + βt (t − t0 )
где ρ0 – плотность при известной температуре t0 .
Поверхностное натяжение. Капиллярные явления.
Молекулы жидкости, находящиеся у поверхности контакта с другой жидкостью, газом или
твердым телом имеют другую энергию, чем молекулы, находящиеся внутри объема жидкости.
Эта энергия пропорциональна площади поверхности раздела S и характеризуется величиной
коэффициента поверхностного натяжения σ, который зависит от материала соприкасающихся сред, чистоты поверхности и температуры.
На поверхности раздела трех фаз (рис. 1.4): твердой стенки, жидкости и газа образуется
краевой угол θ. Величина угла зависит только от природы соприкасающихся сред, и не зависит
от формы сосуда и силы тяжести. Чем хуже смачивающая способность, тем больше краевой
угол. От явления смачивания зависит поведение жидкости в тонких (капиллярных) трубках, погруженных в жидкость. При плохом смачивании жидкость в трубке поднимается над уровнем
свободной поверхности, при хорошем – опускается.
θ
θ
Рис. 1.4. Капиллярные явления
Влияние сил поверхностного натяжения приходится учитывать при работе с жидкостными
приборами для измерения давления, при истечении жидкости из малых отверстий, при фильтрации, в других случаях, когда силы, действующие в жидкости меньше сил капиллярного натяжения.
Пенообразование. Для пенообразования необходимо, чтобы в жидкости находились смачивающие вещества, которые уменьшают поверхностное натяжение. Смачивающие вещества состоит из двух групп: гидрофильной и гидрофобной. Они создают пену, которая представляет
множество пузырьков воздуха. Пенообразующие добавки используются при изготовлении
ячеистого бетона.
Растворимость газов в жидкостях – способность жидкостей растворять в своем объеме
газы. Количество растворенного газа в единице объема жидкости различно для разных жидкостей и изменяется с изменением давления. Относительный объем газа, растворимого в жидкости до ее полного насыщения прямо пропорционален давлению:
Vг
р
=к 2 ,
(1.21)
Vж
р1
гдеVг – объем растворенного газа при нормальных условиях; Vж – объем жидкости; к – коэффициент растворимости; р1 и р2 – начальное и конечное давление.
При понижении давления в жидкости происходит выделение растворенного в ней газа. Выделение происходит интенсивнее, чем поглощение. Растворимость необходимо учитывать при
расчете работы машин и систем высокого давления, при расчете кавитации.
Давление насыщенного пара. При определенных условиях капельные жидкости превращаются в пар и наоборот. Изменение агрегатного состояния зависит от давления паров жидкости,
14
насыщающих пространство над ней при данной температуре. Интенсивное выделение пара по
всему объему жидкости называется кипением. Температура кипения зависит от давления на
поверхности жидкости.
Таким образом, интенсивное выделение пара (кипение) может происходить при низких температурах, если давление на поверхности пониженное. Это необходимо учитывать при анализе
работы водопроводных систем на участках пониженного давления.
Таблица 1.1
Температура кипения воды, ºC
10
40
80
100
Давление на поверхности, Па
1175
7350
19800
101325
1.7. Вязкость. Идеальная жидкость
Вязкость – свойство жидкостей оказывать сопротивление сдвигу соседних слоев при движении жидкости. Все реальные жидкости обладают вязкостью, которая проявляется в виде
внутреннего трения при относительном перемещении смежных частей жидкости. Свойство, обратное вязкости – текучесть. Текучесть характеризует степень подвижности частиц жидкости.
На рис. 1.5 представлена эпюра скорости вязкой жидкости, движущейся в цилиндрической трубе. Вследствие тормозящего влияния стенки слои жидкости будут двигаться с разными скоростями, значения которых возрастают по мере отдаления от стенки. Рассмотрим два слоя жидкости, движущиеся на расстоянии Δy друг от друга. Слой А движется со скоростью u, а слой B о
скоростью u + Δu . Вследствие разности скоростей слой В сдвигается относительно слоя А на
Δu
– относительный
величину Δu , которая является абсолютным сдвигом слоя А по слою В, а
Δy
сдвиг или градиент скорости. Если расстояние между слоями будет мало, то градиент скорости
du
. Можно также сказать, что градиент скорости показывает интенсивможно записать как
dy
ность изменения скорости в данном сечении. В результате сдвига соседних слоев появляется
касательное напряжение трения.
Δy
A
B
Δu
Рис.1.5. Распределение скоростей в сечении трубы при ламинарном движении
Согласно гипотезе Ньютона, касательное напряжение, возникающее при движении жидкости пропорционально скорости деформации объема жидкости
du
τ=μ ,
(1.22)
dy
где μ – коэффициент пропорциональности или динамический коэффициент вязкости.
Единица измерения динамического коэффициента вязкости – Пуаз (П):
кгс ⋅ с
.
1П = 0,1 Па·с = 0,0102
м2
Вязкость капельных жидкостей с увеличением температуры уменьшается, а газов увеличивается. Это связано с различным молекулярным строением жидкостей и газов. Для определения
вязкости при различных температурах используются эмпирические формулы, значения вязко15
сти для различных жидкостей приводятся в справочниках. Например, вязкость воды при температуре 20 ºC равна 0,01 П.
Наряду с понятием динамической вязкости в гидравлике применяется кинематическая вязкость ν, которая представляет собой отношение динамической вязкости жидкости к ее плотности
μ
ν= .
(1.23)
ρ
см 2
.
с
Для определения кинематической вязкости при различных температурах также используются эмпирические формулы, и ее значения для различных жидкостей приводятся в справочниках.
см 2
.
Кинематическая вязкость воды при температуре 20ºC равняется 0,01
с
Вязкость капельных жидкостей мало зависит от давления в диапазоне до 200 атм. Кинематическая вязкость воздуха (газов) зависит и от давления, и от температуры. Для нормальных
см 2
, т.е. почти в 15 раз больше воды, что свяусловий (t = 20 ºC, p = 100000 Па) νвоздуха = 0,157
с
зано с меньшей плотностью воздуха.
Измерение вязкости проводится при помощи специальных приборов, называемых вискозиметрами. Измерение вязкости вискозиметром Энглера студенты выполняют во время лабораторных работ.
Идеальная жидкость – жидкость, в которой отсутствует вязкость. Представляет собой модель реальной жидкости, и это понятие используется для облегчения решения некоторых задач
гидравлики. Выводы, полученные исходя из свойств невязкой жидкости, приходиться корректировать, вводя поправочные коэффициенты, получаемые в результате экспериментальных исследований.
Практическое применение теоретических знаний
Пример 1-1
Давление в баллоне с кислородом для газовой сварки, расположенного на улице при температуре T1 = – 10 °С, равно pl = 107 Па. Каково будет давление при внесении его в помещение
при температуре T2 = – 20 °С?
Единицы измерения – Стокс (Ст): 1 Ст = 1
Ответ: Давление в баллоне будет равно 1,11⋅107 Па.
Пример 1-2
Насколько увеличится давление в системе водяного отопления, если температура теплоносителя увеличилась с 60 до 80 °С. Коэффициент температурного расширения, при давлении
5,9·105 Па и температуре 70 °С, можно принять равным 0,00056 (1/град).
Ответ: Давление в системе водяного отопления увеличится на 225 атмосфер.
16
Лекция 2.
Основы гидростатики, динамики и
кинематики жидкости
2.1. Тема 1. Равновесие жидкости
Основные понятия: дифференциальное уравнение равновесия жидкости; три формы записи
основного уравнения гидростатики; поверхность равного давления; закон Паскаля; абсолютное,
манометрическое и вакуумметрическое давление; сила давления на плоские и криволинейные
поверхности; относительный покой жидкости; закон Архимеда; относительный покой жидкости.
Вопросы, на которые необходимо найти ответ в ходе изучения темы:
uuuuur
ur
1. Что обозначает выражение grad p = ρ g ?
uuuuur
ur
2. Какие выводы можно сделать, если grad p = ρ g ?
3. Как можно получить основное уравнение гидростатики в дифференциальной форме?
4. В каком случае уравнение Бернулли превращается в основное уравнение гидростатики?
5. Запишите основное уравнение гидростатики в трех его формах. Объясните значение каждого входящего в него слагаемого.
6. Какая плоскость называется плоскостью равного давления?
7. В чем разница между напором и давлением?
8. Что называют удельной энергией положения?
9. Объясните закон Паскаля.
10. Какие механизмы действуют на основе закона Паскаля?
11. Что называют абсолютным давлением, манометрическим давлением, вакуумом?
12. Какое давление будет в напорном потоке, если пьезометрическая линия проходит ниже
его геометрической оси?
13. Как определить давление в любой точке жидкости, находящейся в закрытом резервуаре?
От чего оно зависит?
14. Чему равна результирующая сила давления на стенки резервуара произвольной формы?
Как это доказать?
15. Как определить силу давления жидкости на плоскую горизонтальную поверхность?
16. В чем заключается гидростатический парадокс?
17. Как определить силу давления жидкости на плоскую вертикальную поверхность?
18. Как определить точку приложения результирующей силы давления на плоскую вертикальную поверхность?
19. Как определяется величина силы давления на криволинейные поверхности?
20. Как рассчитать допустимое давление в трубе, чтобы не допустить ее разрыва?
21. Сформулируйте закон Архимеда.
22. Какова природа действия Архимедовой силы?
23. Каковы условия плавания тел, их равновесия?
24. Объясните принцип действия ареометра (прибора для измерения плотности жидкости).
25. Что такое относительный покой жидкости?
26. Какие силы действуют при относительном покое жидкости? Как их записать?
27. Как определить положение свободной поверхности жидкости при ее вращении в цилиндрическом сосуде?
17
28. Как распределяется давление в любой фиксированной круглоцилиндрической поверхности при относительном покое жидкости?
H
d
Практические задачи, решение которых может быть найдено после изучения теоретического материала:
Задача 2.1-1
В сообщающихся сосудах находятся вода и масло. Определить плотность масла, если высота столба воды Н = 150 мм, а разность уровней жидкости в сосудах d = ___мм.
Рис. 2.1. К задаче 2.1-1
Задача 2.1-2
Для измерения высоты налива мазута установлена вертикальная труба. В трубу с малой скоростью подают воздух. Определить высоту Н налива мазута удельным весом γ = 8700 Н/м3, если давление воздуха, поступающего в резервуар, эквивалентно высоте ртути в ртутном манометре hрт = dн мм (ρрт = 13600 кг/м3). Трением при движении воздуха пренебречь.
hрт
H
Рис. 2.2. К задаче 2.1-2
Задача 2.1-3
Для повышения гидростатического давления применяется мультипликатор давления, давление на входе которого р1 = 20 кПа, а диаметры поршней d = dн мм и D = 400 мм. Определить
давление жидкости р2 на выходе из мультипликатора.
p2
d
D
p1
Рис. 2.3. К задаче 2.1-3
Задача 2.1-4
Два сообщающихся цилиндра наполнены жидкостью. В меньший цилиндр диаметром d = dн
мм заключен поршень весом G = 100 Н. На какой высоте H установится уровень жидкости в
большом цилиндре, диаметром D = 800 мм, когда вся система придет в равновесие? Удельный
вес жидкости γ = 9,81 кН/м3. Трением пренебречь.
18
G
H
D
d
Рис. 2.4. К задаче 2.1-4
Задача 2.1-5
Диаметры поршней дифференциального предохранительного клапана равны D = 2dн мм и
d=dн мм. Пренебрегая весом поршней и силой трения, определить давление, при котором клапан откроется, если жесткость пружины с = 50 Н/мм, а ее предварительный натяг: x0 = 12 мм.
d
P
D
Рис. 2.5. К задаче 2.1-5
H
Примечание. Сила пружины F = с x0
Задача 2.1-6
Плотность жидкости определяется погружением в нее поплавка. Вес поплавка в воздухе
равняется 0,72 кН. Вес поплавка, погруженного в испытуемую жидкость G1 = 0,54 кН, вес поплавка, погруженного в воду G2 = 0,56 кН. Определить плотность жидкости.
Задача 2.1-7
Бак водонапорной башни сварен из стальных полос высотой. Определить необходимую
толщину стенки нижней полосы, если допускаемое напряжение на разрыв [σ] = 1×108 Па, диаметр бака D = 100dн м, глубина воды в баке Н = 10 м.
Задача 2.1-8
Круглое отверстие в дне резервуара с жидкостью закрыто пластмассовым шариком, вес которого G = 2,45 Н и радиус r = 0,1 м. Диаметр отверстия d = dн мм. Определить давление р, действующее на шар снизу, при котором отверстие откроется, если глубина воды Н = 2 м.
d
p
Рис. 2.6. К задаче 2.1-8
Задача 2.1-9
Плотность жидкости измеряется при помощи ареометра. Размеры аэрометра: d = dн мм, D =
1,5dн мм, Н = 100 мм, h = 50 мм, масса m = 0,054 кг. Определить плотность жидкости.
19
h
d
H
D
Рис. 2.7. К задаче 2.1-9
Задача 2.1-10
Определить абсолютное давление на поверхности жидкости в сосуде и высоту h, если атмосферное давление равняется 740 мм.рт.ст. Поддерживающая сила F = 10 H, вес сосуда G = 2 Н,
его диаметр d = dн мм. Толщиной стенки сосуда пренебречь.
h
F
d
Рис. 2.8. К задаче 2.1-10
2.1.1. Дифференциальное уравнение равновесия жидкости.
Поверхность равного давления
Предположим, что в точке М находится объем жидкости dV (см. рис. 2.9). На него воздействуют силы давления соседних объемов. Определим результирующую силу давления на объем
dV. dV расположен параллельно осям координат, da, db, dc – его стороны. В точке М давление
обозначим как p. В точках M 1 и M 2 , принадлежащих сторонам параллельным плоскости x0y
давление будет соответственно p1 и p2 . Если рассматривать одну из сторон параллелепипеда,
то результирующая сила давления на эту сторону действует по нормали к ней и ориентирована
внутрь объема dV.
z
M2
da
dc
M
db
M1
k
j
x
i
y
Рис. 2.9. Объем жидкости, находящийся в равновесии
Для результирующей силы сторон объема dV, параллельных плоскости x0y можно записать
uuur r
uuur r
dFz k = p1dadb − p2 dadb или dFz k = ( p1 − p2 )dadb ,
uuur
dFz параллельна оси 0z.
Разность ( p1 − p2 ) можно записать в виде p1 − p2 = ( p1 − p ) − ( p2 − p ) , но в соответствии со
свойством градиента давления можно написать
20
uuuuur uuuuur
uuuuur uuuuur
( p1 − p ) = grad pM MM 1 , ( p2 − p ) = grad pM MM 2 ,
uuuuur
uuuuur uuuuur
откуда ( p1 − p2 ) = grad pM ( MM 1 − MM 2 ) .
uuuuur uuuuur uuuuuur uuuuur uuuuuur uuuuuur
r
Так как ( MM 1 − MM 2 ) = M 2 M + MM 1 = M 2 M 1 и M 2 M 1 = (− k )dc , то
uuuuur
uur
( p1 − p2 ) = grad pM (− k )dc .
uuur r uuuuur
uur
Таким образом, результирующая сила dFz k = grad pM (− k )dcdadb , но dcdadb = dV, oткуда
uuur r uuuuur
uur
dFz k = grad pM (− k )dV .
uuur uuur
Аналогичные результаты мы получим для сил dFx и dFy .
Результирующая всех сил, действующих на объем dV будет соответственно
uuur r
uuuuur
dF k = − grad pM dV .
(2.1)
Выводы:
uuur
uuuuur
1. Результирующая сила dF направлена в противоположную сторону, чем grad pM
uuur
2. dF перпендикулярна плоскости, проходящей через точку М, на которой давления одинаковы и ориентирована в сторону уменьшения давления.
В жидкости, находящейся в покое, действуют:
– сила тяжести
ur
ur
dmg = ρdV g ,
направленная вертикально вниз;
– равнодействующая сила давления
uuuuur
ur
d F = − grad pdV ,
ur uuuuur
ρdV g − grad pdV = 0
uuuuur
ur
или grad p = ρ g .
(2.2)
Выводы:
r
1. Вектор градиента давления направлен вертикально вниз, как и вектор g .
2. В жидкости, находящейся в равновесии давление увеличивается сверху вниз.
3. В покоящейся жидкости плоскости равного давления горизонтальны.
4. В покоящейся жидкости давление в точке зависит только от ординаты z.
uuuuuuruuuuur
uruuuur
Т.к. dp = gradpMM 1 , то с учетом полученного уравнения, можно записать dp = ρ g MM 1 . Т.к.
ur
r
uuuur r
g = − k g и MM 1 = k g , то
dp = −ρgdz .
(2.3)
Нами получено основное уравнение гидростатики в дифференциальной форме.
2.1.2. Основное уравнение гидростатики
В случае несжимаемой жидкости плотность жидкости не зависит от давления, а если принять температуру постоянной, то можно записать
ρ = const.
Для высот в несколько метров ускорение силы тяжести можно считать неизменным. Таким
образом, можно подсчитать разность давления между точками М 1 и М 2 . Проинтегрировав
предыдущее выражение можно получить разность давлений между двумя точками:
21
р2
z2
р1
z1
∫ dp = −ρg ∫ dz ,
p2 − p1` = −ρg ( z2 − z1 )
или p1 + ρgz1 = p2 + ρgz2 = const .
(2.4)
Нами получено основное уравнение гидростатики в поле силы тяжести.
Если принять z1 − z2 = h , то p2 − p1 = ρgh .
Выводы:
1. В покоящейся жидкости давление увеличивается с увеличением глубины.
2. В покоящейся жидкости любая горизонтальная плоскость представляет собой поверхность, на которой в любой точке давление будет неизменным. Такая поверхность называется
поверхностью равного давления.
Три формы записи основного уравнения гидростатики.
Нами было получено основное уравнение гидростатики в форме давлений, т.к. каждый
член уравнения представляет собой давление:
р1 и р2 – статическое давление в точках 1 и 2;
ρgz1 и ρgz2 – давление, создаваемое силой тяжести.
Если разделим основное уравнение гидростатики в форме давлений на ρg, то получим основное уравнение гидростатики в форме напоров (см. рис. 2.10)
p1
p
+ z1 = 2 + z2 = const ,
(2.5)
ρg
ρg
где
p1
p
и 2 – пьезометрические напоры; z1 и z2 – геометрические напоры.
ρg
ρg
p1
ρg
p2
ρg
1
2
z1
z2
0
0
Рис. 2.10. Геометрический и пьезометрический напоры
Если первое уравнение разделить на ρ, то получим основное уравнение гидростатики в
форме удельной энергии
p1
p
+ gz1 = 2 + gz2 = const ,
(2.6)
ρ
ρ
где
p1
p
и 2 – удельная энергия давления; gz1 и gz2 – удельная энергия положения.
ρ
ρ
22
2.1.3. Закон Паскаля
Перепишем основное уравнение гидростатики в форме давлений в следующем виде
p1 + ρgz1 = p0 + ρgz0 = const ,
(2.7)
где р0 – давление на свободной поверхности; z0 – расстояние от свободной поверхности до
плоскости сравнения.
Можно записать это уравнение в другом виде
p1 = p0 + ρgz0 − ρgz1 , или p1 = p0 + ρg ( z0 − z1 ) , или p1 = p0 + ρgh ,
(2.8)
где h – глубина, на которой находится точка 1.
Из этого уравнения следует, что изменение давления на свободной поверхности на величину ±Δр приведет к увеличению давления в точке на ту же величину. Этот вывод и есть
закон Паскаля.
Закон Паскаля используется в гидропрессах и гидроусилителях. Схематично гидропресс
πd 2
представлен на рисунке 2.11. При воздействии на малый поршень площадью s =
с силой F
4
F
создается давление p = . Согласно закона Паскаля это давление передается во все точки жидs
F
4 F πD 2
D2
=F 2 .
кости и поршень площадью S создает усилие P = S = 2 ×
4
s
d
πd
F
D
P
d
Рис. 2.11. Схема гидропресса
Таким образом, можно сделать следующий вывод:
1. Сила P будет больше силы F во столько же раз, во сколько площадь S больше площади s.
2. В отличие от твердых тел жидкость передает не силу, а давление.
2.1.4. Абсолютное, манометрическое и вакуумметрическое давление
Давление можно измерить двумя способами:
1) если принять за начало отсчета атмосферное давление;
2) если принять за начало отсчета абсолютный вакуум, когда давление в объеме отсутствует.
В первом случае давление называется избыточным или манометрическим, во втором –
абсолютным.
Избыточное давление в точке 1 (см. рис. 2.10) p1 = ρgh , где h – глубина, на которой находится точка 1.
Абсолютное давление для этой точки будет p1 = p0 + ρgh , но т.к. давление на свободной поверхности в данном случае равняется атмосферному, то p1 = pатм + ρgh . Это давление будет со23
p1
. Манометрическое давление является разностью
ρg
между абсолютным и атмосферным давлением (см. рис. 2.12). В общем случае абсолютное давление рабс может быть больше или меньше атмосферного, если рабс < ратм , то разность между
атмосферным давлением и абсолютным называется вакуумом
рвак = ратм − рабс .
(2.9)
ответствовать пьезометрическому напору
pабс
pман
pатм
0
pвак
0
Рис. 2.12. Шкалы абсолютного, манометрического и вакуумметрического давлений
2.1.5. Сила давления на плоские и криволинейные поверхности
1. Сила давления на отдельный элемент поверхности.
Точка М (см. рис. 2.13) принадлежит площадке ds, являющейся частью некоторой поверхности. Давление на площадку ds – рМ . Сила давления на площадку ds будет равна
uuur
r
df N = pM dsn ,
(2.10)
r
где n – единичный вектор, ориентированный по нормали к площадке ds.
Предположим, что в резервуаре находится жидкость и газ (см. рис. 2.13).
ds
r
n
M
r
fN
Рис. 2.13. Сила давления на отдельный элемент
Давление газа – р0 . Точка М находится на стенке резервуара, на глубине h. Давление в точке
pМ = p0 + ρgh .
Сила давления на элемент стенки ds
uur
r
df = ( p0 + ρgh)dsn
(2.11)
2. Результирующая сила давления на стенку.
ur
uur
Результирующая сила давления P = ∑ df . Если резервуар имеет произвольную форму (см.
рис. 2.14), то подсчитать результирующую силу довольно сложно, т.к. единичные векторы каждого элемента поверхности направлены в разные стороны. Для определения результирующей
силы прибегнем к следующим рассуждениям. Элемент стенки резервуара ds будет находиться в
uur
неподвижном состоянии, если сила давления жидкости df , будет равна силе реакции материала
uur
uur uur
df R , т.е. df + df R = 0.
24
p0
h
M
ds
r
n
r
df
Рис. 2.14. Результирующая сила давления на стенку поверхности
Весь резервуар испытывает воздействие двух результирующих сил:
1. Равнодействующей сил давления
ur
uur
P = ∑ df ;
2. Реакции материала стенки резервуара
ur
uuur
R = ∑ df R
Исходя из предыдущего уравнения, можно записать
ur ur
ur
ur
P + R = 0 или P = − R .
Сама жидкость в резервуаре находится в равновесии под воздействием двух сил:
• силы реакции материала стенки R;
• силы тяжести, направленной вертикально вниз G
ur
ur
R = −G .
ur
ur
Cсравнивая это уравнение с уравнением P = − R , можно сделать вывод, что
ur ur
(2.12)
P=G.
3. Сила давления жидкости на дно резервуара.
В связи с тем, что в резервуаре произвольной формы очень трудно подсчитать результирующую силу на дно резервуара из-за разнонаправленности векторов сил элементарных площадок, ограничимся только случаем, когда дно резервуара плоское и горизонтальное. На рисунке 2.15 изображен открытый резервуар. Давление на поверхности жидкости р0 , плотность –
ρ, глубина наполнения жидкости – h.
p0
r
df
h
r
df ′
Рис. 2.15. Сила давления на горизонтальное дно резервуаров
Так как дно резервуара плоское и горизонтальное, то каждый элемент поверхности дна будет испытывать давление
р = p0 + ρgh ,
и на него воздействует элементарная сила давления со стороны жидкости
uur
r
df = p dsn
и сила давления со стороны наружного воздуха
uuur
r
df ' = p0 dsn ' .
25
uur uuur
Все элементарные силы df и df ' параллельны между собой.
r
ur
uur
Равнодействующая сила давления воды P1 = ∑ df = ∫ p dsn .
S
r
r
ur
Так как p = const P = pS n = ( p0 + ρgh) S n .
Аналогично равнодействующая сила давления воздуха
ur
ur
P 0 = p0 S n '
Эти две силы вертикальны и действуют в разных направлениях.
Результирующая сила давления на дно резервуара
ur ur uur
P = P1 + P0 , или P = P1 + P0 , или P = ρghS .
(2.13)
Сила Р – вертикальная, направлена вниз и приложена по центру дна резервуара (из соображения симметрии).
Гидростатический парадокс. Независимо от формы резервуара сила давления на дно зависит только от площади S, глубины заполнения h и плотности ρ и не зависит от количества жидкости, находящейся в резервуаре (см. рис. 2.16).
P
P
P
h
Рис. 2.16. Гидростатический парадокс
Опыт Паскаля. Резервуар рассчитан на определенное давление жидкости. В него добавляют небольшое количество воды. Ничего не происходит. Вставляют тонкую трубочку и добавляют гораздо меньшее количество воды – резервуар разрушается.
4. Сила давления на вертикальную прямоугольную стенку.
Пусть прямоугольная стенка длиной l и высотой h сдерживает напор воды (жидкости) плотностью ρ (см. рис. 2.17).
B
O
A
z
h
df
C′
z + dz
C
P
O
l
D
z
Рис. 2.17. Сила давления на вертикальную стенку
Рассмотрим элемент стенки, находящейся на глубине z длиной l и шириной dz. Элемент испытывает давление
uur
r
df = pdsn .
df направлена вертикально к поверхности и приложена в центре элемента на оси О (из соображения симметрии).
Давление на глубине z:
р = ρgz.
26
uur
r
Площадь ds = ldz. Тогда df = ρgzdzln .
ur
Сила давления P uuна
стенку равняется сумме сил, действующие на элементарные площадки
r
ur
uur
P = ∑ df . Все силы df горизонтальные, действуют в одном направлении и приложены на одur
ной вертикальной оси О. Сила P также будет горизонтальна, направлена от жидкости, точка
приложения находится на оси О. Можно посчитать силу давления
r
z =h
ρglh 2
(2.14)
P = ∫ ρgzldzn или P =
.
z =0
2
Т.к. lh = S, то
h
P = ρgS .
(2.15)
2
Определим точку приложения силы Р.
O
2/3h
h
M
C
r
P
z
Рис. 2.18. Определение точки приложения силы Р
Стенка испытывает воздействие всех сил df (см. рис. 2.18). Точка приложения С должна
быть расположена таким образом, чтобы воздействие силы Р в этой точке равнялось воздейстur
uur
ur
ur
ur
вию всех сил df на площадку ds Т.е. P = ∑ df и M O ( P ) = ∑ M O (d f ) , где M O ( P) – момент сиur
ur
ur
лы P относительно точки О; ∑ M O (d f ) – сумма моментов сил d f относительно точки О.
ur
ur
ur
Для P момент силы M O ( P) = P × OC . Для силы d f , приложенной в точке М на глубине z:
ur
M O (d f ) = df × OM , где ОМ = z.
Таким образом P × ОС = ∫
z =h
z =0
ρglz 2 dz :
ρglh 2
ρglh 2
ρglh 2
× ОС =
или
, откуда
3
2
3
2
ОС = h
(2.16)
3
5. Сила давления на криволинейную поверхность.
Рассмотрим поверхность S, на которую с внешней стороны воздействует жидкость, создавая
давление р2 , воздействие жидкости с внутренней стороны – р1 (см. рис. 2.19). Каждый элемент
поверхности площадью ds испытывает воздействие силы давления с внешней стороны
uur
r
df1 = p1dsn ,
r
где n – единичный вектор, направленный по нормали к ds, ориентированный в сторону внешней жидкости.
P × ОС =
27
Δ
p2
r
j
r
n
θ
dS
p1
Рис. 2.19. Определение силы давления на криволинейную поверхность
Равнодействующая от суммы всех элементарных сил, действующих на поверхность изнутри
ur
uur
P1 = ∑ df 1 :
uur
r
Р1 = p1 ∫ dsn.
S
r
Эту силу трудно подсчитать, т.к. векторы n не параллельны между собой.
r
ur
Проекция силы P1 , по направлению единичного вектора j будет
uurr
rr
Р1 j = p1 ∫ n jds .
S
rr
Если θ – угол между направлением Δ и нормалью к поверхности ds, то n j = cos θ , тогда
uur r
р1 Р1 j = p1 ∫ ds × cos θ .
S
uurr
uur r
В то же время ds × cos θ = ds', тогда Р1 j = p1 ∫ ds ' , или Р1 j = p1dS ' , где dS' – проекция поS
верхности S на поверхность, перпендикулярную выбранному направлению Δ.
uurr
Таким образом, можно записать Р1 j = Р1 , как произведение двух векторов и
Р1 = p1dS ' .
(2.17)
Такой же результат мы получим для равнодействующей силы Р2
Р2 = p2 dS ' .
(2.18)
6. Сила давления на цилиндрическую поверхность.
Представим трубу длиной l, с внутренним радиусом R1 и внешним радиусом R2 (см. рис.
2.20). Толщина трубы e = R2 – R1 . В трубе находится жидкость с давление. Внешнее давление –
р2 . Необходимо установить минимальную толщину e, при которой труба не разорвется.
R1
R2
Δ
l
Рис. 2.20. Определение силы давления в цилиндрической трубе
28
Выберем направление Δ, совпадающее с одним из радиусов трубы.
Равнодействующая сил внутреннего давления
Р1 = p1dS ' ,
где S' – проекция поверхности S на плоскость, перпендикулярную направлению Δ. Поверхность
S' является прямоугольником, площадь которого 2 R1 l, тогда
Р1 = p1 2 R1 l .
Аналогичный результат можно получить для силы внешнего давления
Р2 = p2 2 R2 l .
Для второй части трубы мы получим тот же результат.
Если внутреннее давление будет больше внешнего Р1 > Р2 (сила Р1 будет стремиться разорвать трубу), направление Δ было выбрано произвольно, поэтому можно сделать вывод, что
разрыв может произойти по любому направлению. Материал трубы в силу своих физических
свойств будет сопротивляться разрыву. Это сопротивление будет тем больше, чем толще будет
труба. Величина, характеризующая способность материала сопротивляться его разрыву обозначается σ. Сила сопротивления материала
F = σSсопр ,
где Sсопр – площадь сопротивления.
Таким образом, разрыв произойдет в случае, если Р1 − Р2 ≥ F , или p1 2 R1 l − p2 2 R2 l ≥ σSсопр ,
или p1 2 R1 l − p2 2( R1 + e)l ≥ σSсопр , или p1 2 R1 l − p2 2 R1l − p2 2el ≥ σ2el , или p1 R1 − p2 R1 ≥ σe + р2 е ,
откуда следует, что при
R ( p − p2 )
e≤ 1 1
(2.19)
σ + p2
произойдет разрыв трубы.
2.1.6. Относительный покой жидкости
Относительным покоем жидкости называется такое ее состояние, при котором каждая ее
частица сохраняет свое положение относительно твердой стенки движущегося резервуара, в котором находится жидкость (см. рис. 2.21).
z
M
z
z0
dFg
dFr
dF
ω
Рис. 2.21. Относительное равновесие жидкости во вращающемся сосуде
При относительном покое рассматриваются две задачи: определяется форма поверхности
уровня или равного давления и выясняется характер распределения давления. В данном случае
необходимо учитывать силы инерции, дополняющих систему массовых сил, действующих в
покоящейся жидкости.
Рассмотрим случай, когда сосуд с жидкостью вращается вокруг своей оси с постоянной скоростью. Для определения формы свободной поверхности и закона распределения давления выберем вблизи свободной поверхности частицу жидкости массой dm. На эту частицу действует
29
массовая сила dF, направленная по нормали к поверхности. Разложим эту силу на две составляющие: горизонтальную dFr = dmω2 r и вертикальную dFg = −dmg .
Разделив действующие силы на dm, получим дифференциальное уравнение поверхности
уровня
ω2 xdx + ω2 ydy − gdz = 0 или ω2 rdr − gdz = 0 .
Проинтегрировав, получаем
ω2 r 2
(2.20)
− gz = const .
2
Вывод: При вращении резервуара с постоянной скоростью вокруг вертикальной оси поверхностями равного давления будет семейство параболоидов вращения.
Для точки М, находящейся на свободной поверхности жидкости
ω2 r 2
z = z0 +
.
2
Закон распределения давления найдем из дифференциального уравнения гидростатики, которое в данном случае примет вид
dp = ρ(ω2 rdr − gdz ) .
После интегрирования с учетом граничных условий ( r = 0, z = z0 , p = p0 ), получаем:
ω2 r 2
p = p0 + ρ
+ ρg ( z0 − z ) .
2
Если представить, что
p0i = p0 + ρ
ω2 r 2
,
2
то получим уравнение
p = p0i + ρg ( z0 − z ) .
(2.21)
Вывод: Распределение давления подчиняется линейному закону для любой фиксированной
цилиндрической поверхности.
2.1.7. Закон Архимеда
1. Равновесие твердого тела в жидкости.
На тело, находящееся в жидкости (см. рис. 2.22) действуют:
1) сила тяжести;
2) сила давления воды.
S
r
df
M
ds
G
Рис. 2.22. Равновесие твердого тела в жидкости
На каждый элемент поверхности площадью ds действует сила
ur
r
d f = pnds ,
r
где n – единичный вектор, направленный по нормали внутрь тела.
30
Давление p зависит только от положения точки M, соответственно, каждая элементарная сиur
ла d f будет зависеть от расположения элемента и от его формы. Таким образом, равнодейстur
вующая сил d f будет зависеть только от места нахождения тела в жидкости и от величины
внешней поверхности тела S.
2. Равновесие жидкости.
Рассмотрим равновесие жидкости, когда тело извлечено, и жидкость заняла его объем (см.
рис. 2.23). Контур S в жидкости соответствует очертанию тела. Жидкость, находящаяся внутри
этого контура находится в равновесии под действием: 1)собственного веса G'. 2) результирующей силы внешнего давления воды P.
G′
Рис. 2.23. Равновесие объема жидкости в жидкости
Равновесие жидкости, находящейся внутри контура S можно записать в следующем виде:
G'+P = 0 или G' = – P.
Возвращаясь к случаю плавающего в жидкости тела, мы можем сделать вывод, что, т.к. контуры тела S и жидкости, замещающей объем тела равны, то и силы давления жидкости в первом
и во втором случаях должны быть равны. Тогда для тела плавающего в жидкости можно записать
PA = G ' = ρж gVт ,
(2.22)
где PA – сила, действующая на погруженное в жидкость тело; ρ ж – плотность жидкости; Vт –
объем тела.
Закон Архимеда: На твердое тело, погруженное в покоящуюся жидкость, действует сила
гидростатического давления, равная весу жидкости в объеме тела, направленная вертикально
вверх и проходящая через центр тяжести тела.
3. Условия равновесия плавающих тел.
Т.к. вес тела G = ρ т gVт , а архимедова сила PA = ρ ж gVт , то в случае неравенства плотностей
тела жидкости, эти силы будут приложены к различным точкам. Для того чтобы плавающее тело было устойчивым, необходимо, чтобы центр приложения архимедовой силы (центр водоизмещения) находился выше центра приложения силы тяжести.
2.1.8. Основное уравнение гидростатики для сжимаемой жидкости
Выведем основное уравнение гидростатики, учитывая влияние сжимаемости газа, т.е. изменение плотности газа под действием давления. Дифференциальное уравнение равновесия для
среды переменной плотности после интегрирования примет вид
dp
(2.23)
∫ ρ + gz = const.
Для вычисления интеграла необходимо задать закон изменения состояния газа.
31
2.1.9. Изотермическая атмосфера
Если принять, что температура газа не изменяется с изменением высоты, а уравнение состояния записать в следующем виде
p
p
ρ=
= ,
RT C
p
откуда C = – постоянная величина.
ρ
Подставив последнее соотношение в предыдущее уравнение, получим
p
C ln p + gz = const или gz + ln p = const.
ρ
(2.24)
Полученное уравнение отличается от основного уравнения гидростатики жидкости тем, что
давление газа по высоте с учетом его сжимаемости в изотермических условиях распределяется
не по линейному, а по логарифмическому закону.
Запишем предыдущее уравнение для двух высот: на поверхности земли, где z0 = 0, давление
p
p p
p = p0, на высоте z давление будет р. Обозначим z – z0 = h и 0 = H , учитывая, что = 0 ,
ρ0 g
ρ ρ0
p
получим h = H ln , откуда
p0
h
−
p
=e H .
p0
(2.25)
Вывод: Давление уменьшается по высоте по экспоненциальному закону.
Примечание: Полученной формулой можно пользоваться, если высота изменяется на небольшую величину, в пределах нескольких сот метров, т.к. в других случаях формула дает погрешность более 5%.
2.1.10. Неизотермическая атмосфера
Обычно температура воздуха с увеличением высоты уменьшается. Температура дымовых
газов также быстро уменьшается в трубе. Для того, чтобы проинтегрировать основное уравнение гидростатики в этом случае необходимо знать закон изменения температуры с изменением
высоты. Чаще всего для атмосферного воздуха принимают, что температура уменьшается по
линейному закону
T − T0
a=
,
(2.26)
h
где Т0 – температура на поверхности земли; a – градиент температуры.
Можно записать дифференциальное уравнение равновесия в следующем виде
dp
dz
= −g
,
(2.27)
p
R (T0 − az )
g
⎛ T − az ⎞ Ra
p
g T0 − az
=
ln
, или p = p0 ⎜ 0
тогда ln
⎟ , или
p0 Ra
T0
⎝ T0 ⎠
g
Ra
(2.28)
⎛T ⎞
p = p0 ⎜ ⎟ .
⎝ T0 ⎠
Пример: Предположим, что температура понижается на 6ºС с увеличением высоты на 1000
м. Необходимо определить давление на вершине горы Монблан (z = 4800м), если на уровне моря температура Т0 = 30ºС.
32
Решение:
На вершине горы температура Т = 30 – 4800·0,006. Т = 274,2 К, тогда
9,81
p ⎛ 274,3 ⎞ 287⋅4800
=⎜
= 0,557.
⎟
p0 ⎝ 303 ⎠
Вывод: Давление на вершине горы Монблан меньше почти в два раза.
Примеры решения практических задач гидростатики
Пример 2.1-1
К закрытому резервуару для определения давления на свободной поверхности р0 присоединена стеклянная трубка. Спрашивается, какое давление в резервуаре р0, если вода в трубке поднялась на высоту Н = 3 м? Трубка присоединена на глубине h1 = 2 м.
pатм
p0
A
H
h2 A
h1
Рис. 2.24. К примеру 2.1-1
Ответ: Давление на поверхности воды в резервуаре p0 = 107910 Па.
Пример 2.1-2
Определить разность давления в резервуарах А и В, заполненных водой, если разность уровней ртути в U-образном манометре h = 15 см.
h
h2
h1
B
A
Рис. 2.25. К примеру 2.1-2
Ответ: Разность давлений в резервуарах – 18394 Па.
Пример 2.1-3
Определить суммарное усилие, воспринимаемое болтами смотрового люка диаметром d =
0,5 м, расположенного на глубине h = 3 м от свободной поверхности, на которой давление p0 =
0,5 атм.
Ответ: Усилие, воспринимаемое болтами смотрового люка Р = 15700 Н.
Пример 2.1-4
При бурении скважины необходимо определить вес труб, опущенных в скважину, заполненную глинистым раствором плотностью ρр = 2800 кг/м3, длина труб l = 70 м. Один метр таких
труб с муфтами в воздухе весит 300Н. Плотность стали ρст = 7500 кг/м3.
Ответ: Полный вес труб будет равен 902776 Н.
33
p0
pатм
h
d
P
Рис. 2.26. К примеру 2.1-3
2.2. Тема 2. Основы кинематики и динамики жидкости и газа
Основные понятия: расход; мгновенная и средняя скорость; линия тока; трубка тока; уравнение неразрывности; установившееся и не установившееся движение жидкости; равномерное
и неравномерное движение; сплошная среда; количество движения; момент количества движения; дифференциальное уравнение движения невязкой жидкости, общее уравнение энергии в
интегральной форме; три формы представления уравнения Бернулли для потока реальной жидкости; турбулентный и ламинарный режимы течения жидкости; число Рейнольдса; осредненная
и мгновенная скорости; пульсации; турбулентные касательные напряжения.
Вопросы, на которые необходимо найти ответ в ходе изучения темы:
1. Что называют линией тока, трубкой тока, элементарной струйкой?
2. Чем траектория движения частицы отличается от линии тока?
3. Какие свойства трубки тока вы знаете?
4. Зачем в гидравлике вводятся понятия: трубка тока, элементарная струйка, линия тока?
5. В чем суть математической модели Эйлера, которая применяется для описания движения жидкости?
6. Что такое расход жидкости, скорость, средняя скорость, живое сечение потока? Назовите единицы измерения.
7. Какова связь между массовым расходом и объемным?
8. Чем отличается установившееся движение жидкости от неустановившегося, равномерное от неравномерного? Приведите примеры.
9. Как записать вид движения жидкости в математической форме?
10. Привести примеры напорного и безнапорного движения жидкости.
11. Что называют сплошной средой?
12. Объяснить суть уравнения неразрывности.
13. Как уравнение неразрывности можно записать в математической форме?
14. Какие выводы можно сделать из уравнения неразрывности?
15. Для каких условий применима струйная расчетная модель движения реальной жидкости?
16. Как математически записать суммарную проекцию внешних сил, действующих на изолированную массу жидкости в наклонном канале?
17. Что называют коэффициентом Буссинеска?
18. Как записать изменение количества движения жидкости в единицу времени при установившемся движении жидкости?
19. Чему должен быть равен импульс действующих сил при движении жидкости?
20. Какие силы действуют на выделенный объем идеальной жидкости, движущийся равномерно?
34
21. Как записать уравнение равновесия движущегося равномерно объема идеальной жидкости?
22. Написать систему уравнений Эйлера и объяснить ее смысл.
23. Как из системы уравнений Эйлера получить уравнение Бернулли для струйки невязкой
жидкости?
24. Как из уравнения Бернулли для струйки невязкой жидкости получить уравнение Бернулли для двух сечений потока реальной жидкости?
25. В чем отличие уравнения Бернулли для элементарной струйки идеальной жидкости от
уравнения Бернулли потока реальной жидкости?
26. Для чего в уравнение Бернулли вводится коэффициент Кориолиса?
27. Какие выводы можно сделать из уравнения Бернулли?
28. Какие ограничения существуют в применении уравнения Бернулли?
29. Как построить линию полного напора, пьезометрическую линию?
30. В каком случае линия полного напора будет параллельна пьезометрической линии?
31. Записать уравнение Бернулли в форме удельной энергии и объяснить его суть.
32. Записать уравнение Бернулли в форме напоров и объяснить его суть.
33. Записать уравнение Бернулли в форме давления и объяснить его суть.
34. Для каких инженерных расчетов применяется та или иная форма записи уравнения Бернулли? Привести примеры.
35. Чем турбулентное движение жидкости отличается от ламинарного?
36. В чем физический смысл числа Рейнольдса, его практическое значение?
37. Чем характеризуется переходная зона от одного режима к другому?
38. Чему равно значение коэффициента Кориолиса при турбулентном и ламинарном режимах?
39. При каком режиме движения жидкости потери напора больше? Почему?
40. От каких параметров зависят гидравлические потери в ламинарном потоке?
41. Что такое пульсации скорости и давления, чем они вызваны, когда возникают?
42. В чем заключается теория Прандтля?
43. Что называют мгновенной, осредненной скоростью?
44. Чем осредненная скорость отличается от средней?
Практические задачи, решение которых может быть найдено после изучения теоретического материала:
Задача 2.2-1, 6
В закрытом резервуаре поддерживается постоянное манометрическое давление pм = 2 атм,
под действием которого по трубе диаметром d1 = dн мм и общей длине l = 25 м (расстояние от
начала до колена l1 = 5 м) вытекает жидкость при температуре t. Определить расход при напоре
H1 = 2 м.
l1
H1
Рис 2.27. К задаче 2.2-1, 6
Примечание. Потерями напора пренебречь.
35
Задача 2.2-2, 7
Истечение происходит из открытого резервуара при постоянном напоре H1 = 5 м по короткому трубопроводу переменного поперечного сечения с диаметрами d1 = 10dн мм и d2 = 3dн мм
в атмосферу. На втором участке трубопровода имеются два колена с плавным поворотом. Длина первого участка l1 = 0,8 м, длина второго участка l2 = 2 м. Определить скорость истечения υ2
и расход Q2 через трубопровод.
d1
H1
d2
H2
υ1
υ2
l2
l1
Рис 2.28. К задаче 2.2-2, 7
Примечание. Потерями напора пренебречь.
Задача 2.2-3, 8
Из резервуара А, заполненного жидкостью на высоту H1 = 5 м и находящегося под манометрическим давлением рм = 20 кПа, жидкость подается по трубопроводу длиной l = 6 м и диаметром d = 10dн мм в резервуар В на высоту Н = 2 м. Определить расход Q и скорость протекающей по трубопроводу жидкости.
Примечание. Потерями напора пренебречь.
M
H
H1
B
A
Рис. 2.29. К задаче 2.2-3, 8
Задача 2.2-4, 9
К закрытому резервуару, на свободной поверхности которого действует манометрическое
давление рм = 10 кПа, подсоединен трубопровод переменного сечения с диаметрами d1 = 2dн мм
и d2 = 5dн мм, заканчивающийся соплом диаметром dc = d1. Трубопровод подсоединен на глубине Н1 = м. На первом участке длиной l1 = 10 м установлен вентиль. Длина второго участка l2 = 5
м. Определить скорость истечения υ и расход Q вытекающей из сопла жидкости при температуре t и постоянном напоре H1.
36
М
d1
d2
dс
H1
υ
l2
l1
Рис. 2.30. К задаче 2.2-4, 9
Примечание. Потерями напора пренебречь.
Задача 2.2-5, 10
Центробежный насос подает воду с температурой t = 15 °С по трубе диаметром d = dн мм и
длиной l = 27 м. Насос создает давление на выходе рн = 2,65 атм и подачу Q = л/с. Определить
напор H, создаваемый насосом.
H
Рис. 2.31. К задаче 2.2-5, 10
Примечание. Потерями напора пренебречь.
2.2.1. Основные понятия кинематики жидкости
Для описания движения жидкости используется математическая модель. В гидравлике наибольшее распространение получила модель Эйлера, суть которой можно объяснить следующим
образом. Предположим, что точка М движется по некоторой траектории в системе неподвижных координат. Мгновенное значение составляющих скорости вдоль осей координат будет зависеть от положения точки, т.е. от величины координат x, y, z и времени t. Для составляющих
скоростей течения жидкости в рассматриваемой точке u x , u y , u z (см. рис. 2.32), можно записать
функциональные зависимости:
z
uz
uy
M
ux
x
y
Рис. 2.32. Скорость в точке
37
u x = f1 ( x, y, z, t )
u y = f 2 ( x, y , z , t )
(2.29)
u z = f 3 ( x, y , z , t )
Зная для конкретного случая течения значения этих функций, можно для любого момента
времени получить распределение скоростей течения жидкости.
Расход – количество жидкости, проходящей в единицу времени через данное сечение трубопровода. Различают объемный и массовый расходы.
Объемный расход – объем жидкости, проходящий в единицу времени через данное сечение
трубопровода:
Q=
V ⎛ м3 ⎞
⎜ ⎟,
t ⎝ с ⎠
(2.30)
где V – объем жидкости.
Массовый расход – масса жидкости, проходящая в единицу времени через данное сечение:
Qm =
m ⎛ кг ⎞
⎜ ⎟.
t ⎝ с ⎠
(2.31)
Соответственно, Qm = ρQ , где ρ – плотность жидкости.
Траектория – кривая, вдоль которой происходит перемещение частицы жидкости.
Линия тока – кривая, в каждой точке которой вектор скорости движения частицы направлен по касательной к ней (см. рис. 2.33).
Рис.2.33. Линия тока
Трубка тока – поверхность, очерченная вдоль небольшого контура внутри которой вдоль
линии тока перемещаются частицы жидкости. Стенки трубки тока непроницаемы. Площадь поперечного сечения трубки тока мала, поэтому скорости движения в каждой точке равны (см.
рис. 2.34).
Рис. 2.34. Трубка тока
Элементарная струйка – поток жидкости, протекающий в трубке тока Элементарную
струйку можно представить также как совокупность линий тока, проходящих через бесконечно
малое сечение ds, а разность скоростей соседних линий тока бесконечно мала. Расход элементарной струйки dq = uds. Поток жидкости можно представить как совокупность трубок тока, в
которых движутся элементарные струйки.
Q = ∫ dq .
S
Средняя скорость потока – скорость, одинаковая в каждой точке потока в данном сечении,
соответствует реальному расходу
υср =
∑u
i
,
i
где ui – скорость в точке в данном сечении; i – количество точек.
Для потока жидкости, состоящего из нескольких трубок тока можно записать
38
υср =
∑ dq
∑ ds
или υср =
Q
,
S
(2.32)
где S – площадь сечения потока жидкости.
2.2.2. Уравнение неразрывности
Основным условием, которое должно соблюдаться при течении жидкости, является непрерывность изменения параметров потока в зависимости от координат и времени, т.е. при течении
жидкости должны быть соблюдены условия при, которых жидкость должна двигаться в канале
как сплошная среда, без разрывов.
Выделим внутри пространства с движущейся капельной жидкостью неподвижный контур в
форме элементарного параллелепипеда с ребрами dx, dy, dz (см. рис. 2.35). Обозначим скорость
жидкости, которая втекает в левую грань параллелепипеда, через ux . Скорость жидкости, вытекающей из правой грани, вследствие неразрывности поля скоростей равна
z
ux +
ux
r
k
r
j
dx
r
i
dz
∂u x
dx
∂x
dy
x
y
Рис. 2.35. Движение жидкости через контур
∂u x
dx .
∂x
Поскольку рассматриваемый элементарный объем неподвижен, изменение скорости не зависит от времени. В направлении оси х через левую грань втечет за 1 с жидкость массой
( ux dydzρ ) , а вытекает через правую грань
ux +
⎡ ⎛
∂u x ⎞
⎤
⎢ +ρ ⎜ u x + ∂x dx ⎟ dydz ⎥ .
⎠
⎣ ⎝
⎦
Значит, за 1 с из параллелепипеда вытекает в направление оси х жидкости больше, чем втекает, на
∂u
ρ x × dxdydz.
∂x
Аналогичные выражения получаются и для направлений x, y, z. Закон сохранения массы
требует, чтобы сумма трех полученных приращений была равна нулю:
∂u y
∂u
∂u
dxdydz + ρ z dxdydz = 0 .
ρ x dxdydz + ρ
(2.33)
∂x
∂y
∂z
Это уравнение называют уравнением неразрывности, т.к. оно предполагает, что жидкость
является сплошной средой.
Рассмотрим уравнение неразрывности для случая течения струйки при установившемся
движении. Масса жидкости течет в трубке тока (см. рис. 2.34). Пусть левое входное сечение
трубки тока имеет площадь s1 и в этом сечении скорость жидкости u1 , а ее плотность ρ1 . Площадь сечения на выходе из трубки тока s2 , скорость течения жидкости u2 , и ее плотность ρ2 .
Скорости струйки направлены по касательной к стенкам трубки тока, поэтому через стенки об39
мен массой с окружающей жидкостью отсутствует. Через левое сечение втекает в единицу времени масса жидкости Qm1 = ρ1s1u1 . Через правое сечение вытекает в единицу времени масса
жидкости Qm 2 = ρ2 s2u2 . В трубке тока масса жидкости, находящаяся между левым и правым сечениями, остается постоянной, следовательно, условие сплошности потока в трубке тока будет:
Qm1 = Qm 2 = const.
(2.34)
Если плотность жидкости по длине трубки тока не изменяется, т.е. ρ1 = ρ2 , то можно записать для левого и правого сечений:
s1u1 = s2u2 = const или Q1 = Q2 = const.
(2.35)
Полученное уравнение является уравнением неразрывности для трубки тока.
Для потока реальной жидкости уравнение неразрывности записывается в следующем виде:
S1υср1 = S 2 υср 2 = const ,
(2.36)
где S1 и S2 – площади сечения потока в сечениях на входе и на выходе; υср1 и υср2 – средние
скорости потока в этих сечениях.
Можно сделать два важных вывода:
1. При установившемся движении жидкости объемный расход не меняется;
2. При увеличении площади сечения потока жидкости средняя скорость уменьшается, и,
наоборот, при уменьшении сечения - скорость увеличивается.
2.2.3. Виды движения жидкости
Все случаи течения жидкости можно разделить на виды, представленные на рисунке 2.36.
Рис. 2.36. Виды движения жидкости
Установившееся движение жидкости – движение жидкости, при котором все параметры
жидкости (давление, температура, скорость и др.) не изменяются по времени.
u = f ( x, y , z ) .
(2.37)
Для неустановившегося движения:
u = f ( x, y , z , t ) .
(2.38)
Равномерное движение – установившееся движение, при котором скорость по всей длине
потока не изменяется:
u = const .
(2.39)
40
Напорное движение устанавливается в закрытых гидравлических системах, в которых жидкость течет в, основном, под действием силы давления, безнапорное движение наблюдается в
открытых системах, в которых движение жидкости происходит под действием силы тяжести.
2.2.4. Интегральная формула количества движения
В теоретической механике изучается теорема о количестве движения материальных точек,
которую можно применить в гидравлике. Будем считать движение жидкости в канале (см. рис
2.37) установившемся. Распределение давления в сечениях 1-1 и 2-2 гидростатическое. В сечении 1-1 действует сила давления, направленная внутрь выделенного объема жидкости P1 = p1s1 ,
в сечении 2-2 – сила P2 = p2 s2 ( p1 , p2 – давление в сечениях 1-1 и 2-2, s1 , s2 – площади сечений). На элементарной площадке стенки канала действуют: сила реакции стенки, равная силе
давления pds и сила вязкого трения τds. Проекция этих сил на ось движения равна pds sin α +
τds cos α (α – угол наклона элементарной стенки какала к его оси). Суммарная проекция внешних сил, действующих на изолированную массу жидкости, равна
∫ p sin αds + ∫ τ cos αds + p s − p s
1 1
S
2 2
.
S
2
pds
τds
1
θ
P2
P1
α
2
1
Рис. 2.37. Силы, действующие на поток жидкости в канале
Примем, что ось канала наклонена под углом θ к горизонту. Проекция веса жидкости, заключенной между выбранными сечениями равна gρVsin θ (V – объем выделенной массы жидкости).
Проекция всех сил, действующих на изолированную массу, равна
∫ p sin αds + ∫ τ cos αds + p s − p s
1 1
S
2 2
– gρVsin θ.
S
(2.40)
Распределение скорости в контрольных сечениях может быть неравномерным. Через элементарную площадку ds контрольной поверхности в единицу времени переносится количество
движения, равное ρu 2 ds (u – местная скорость течения). Суммарное количество движения равs
но ∫ ρu ds . Средняя скорость в контрольном сечении υср = ∫ υds / s . Количество движения, под2
S
0
считанное по средней скорости равно ρυ s.
2
ср
Отношение количества движения, действительно перенесенного потоком, к количеству
движения, определенного по средней скорости, называется коэффициентом Буссине2
ска α 0 = ∫ υ2 s ( υср
s) .
s
В единицу времени при установившемся движении изменение количества движения составит
2
2
ρ ( α 02 υср2
s2 − α 01υср1
s1 ) = ρQ ( α 02 υср2 − α 01υср1 ) ,
где Q – расход жидкости.
41
Импульс действующих сил должен равняться изменению количества движения массы, на
которую данный импульс действует. Следовательно, при течении жидкости в канале с учетом
принятых условий соблюдается равенство
∫ p sin αds + ∫ τ cos αds + p s − p s
1 1
S
2 2
– gρVsin θ =
S
ρQα 02 υср2 − ρQα 01υср1 .
(2.41)
Нами получено гидравлическое уравнение количества движения.
Вывод: При переходе от сечения 1-1 к сечению 2-2 проекция секундного количества движения потока изменяется на величину, равную сумме проекций всех внешних сил, действующих
на объем потока, заключенный между сечениями 1-1 и 2-2.
Применяя уравнение количества движения, можно решить ряд задач гидравлики, например:
теоретически определить коэффициент местного сопротивления при внезапном расширении
трубы; изменение давления в трубе при равномерном отборе жидкости на одном ее участке; определение расхода газа при всасывании его через цилиндрическую трубу.
2.2.5. Дифференциальное уравнение движения невязкой жидкости (уравнение Эйлера)
Предположим, что в жидкости движется элементарный объем в форме параллелепипеда со
сторонами dx, dy, dz (см. рис. 2.38). На параллелепипед действуют поверхностные силы давления и массовые силы с проекциями X, Y, Z, отнесенными к единице массы. При движении объема возникают силы инерции. Проекции этих сил на оси координат, отнесенные к единице масdu y du z
du
,−
.
сы, равны соответственно: − x , −
dt
dt
dt
z
pdydz
dz
dy
( p + Δp )dydz
dx
x
0
y
Рис. 2.38. Схема равномерного движения объема жидкости
Рассмотрим условие равновесия сил, в проекции на ось х. Сила давления на левую грань –
pdydz, на правую грань
∂p
– (p + dx )dydz,
∂x
∂p
dx =Δp – изменение давления вдоль оси x.
где
∂x
Массовая сила равна Xρdxdydz. Уравнение равновесия запишется в виде
du
∂p
pdydz – (p + dx ) dydz + Xρdxdydz − x ρdxdydz = 0,
∂x
dt
или
dv
∂p
–
dxdydz + Xρdxdydz − x ρdxdydz = 0.
∂x
dt
Разделив каждый член уравнения на ρdxdydz, получим
1 ∂p du x
=
.
X–
ρ ∂x
dt
42
Соответственно для осей и уравнение равновесия будет выглядеть следующим образом
1 ∂p du y
=
Y–
,
ρ ∂y
dt
Z–
1 ∂p du z
=
.
dt
ρ ∂z
Объединив полученные уравнения, получим систему уравнений Эйлера:
1 ∂p du x
=
X–
,
ρ ∂x dt
Y–
1 ∂p du y
=
,
ρ ∂y
dt
Z–
(2.42)
1 ∂p du z
=
.
dt
ρ ∂z
Можно получить полный дифференциал уравнений Эйлера для установившегося движения,
если рассматривать перемещение частиц жидкости вдоль линии тока. Для этого надо умножить
каждое из уравнений системы на соответствующую проекцию элементарного перемещения
частиц dx, dy, dz, и сложить их:
du y
⎞
1 ⎛ ∂p
∂p
∂p ⎞ ⎛ du x
du
dx +
dy + z dz ⎟ = 0.
⎜ dx + dy + dz ⎟ – ⎜
ρ ⎝ ∂x
∂y
∂z ⎠ ⎝ dt
dt
dt
⎠
Т.к. для установившегося течения линии тока совпадают с траекториями движения частиц,
(Xdx + Ydy + Zdz) –
то
dx
dy
dz
= ux , = u y , = uz .
dt
dt
dt
Тогда
du
du x
du
1
1
dx + y dy + z dz = u x du x + u y du x + u z du z = d (u x2 + u y2 + u z2 ) = du 2
2
2
dt
dt
dt
Для установившегося движения давление зависит только от координат, поэтому второй член
уравнения есть полный дифференциал давления dp. Получим
1
1
dp – du 2 = 0
(Xdx + Ydy + Zdz) –
(2.43)
2
ρ
Мы получили дифференциальное уравнение движения невязкой жидкости.
В поле силы тяжести
X = 0, Y = 0, Z = – g,
тогда уравнение запишется в следующем виде
– gdz –
dp du 2
–
= 0.
ρ
2
После интегрирования этого уравнения получаем (при ρ = const) уравнение
p u2
+
= const,
gz +
ρ
2
(2.44)
которое называется уравнением Бернулли для струйки невязкой жидкости.
2.2.6. Общее уравнение энергии в интегральной форме
(Уравнение Бернулли для потока реальной жидкости)
Для двух сечений струйки невязкой жидкости это уравнение будет выглядеть следующим
образом
43
gz1 +
p1 u12
p u2
+ = gz2 + 2 + 2 .
ρ 2
ρ
2
(2.45)
Сумма слева представляет полную удельную энергию струйки в сечении 1-1, сумма справа –
полную удельную энергию струйки в сечении 2-2. Можно записать, что
E1 = E2 .
На практике энергия струйки в начале больше энергии струйки в конце, т.к. часть энергии
теряется на преодолении сил вязкости. В процессе движения вязкой жидкости запас ее механической энергии уменьшается, и на самом деле
E1 > Е2 .
Обозначим энергию, затрачиваемую на преодоление сил сопротивления Eпот. Eпот – это та
часть механической энергии, которая, вследствие вязкости, переходит в тепловую энергию.
Другими словами можно сказать, что Eпот – это часть энергии, которая израсходована на преодоление гидравлических сопротивлений.
Е1 = Е2 + Eпот.
(2.46)
При выводе уравнения Бернулли для элементарной струйки можно было пренебречь изменением скорости и давления в пределах нормальных сечений благодаря их малым величинам. В
потоке жидкости скорости и давления в пределах живых сечений различны, и это необходимо
учитывать. Согласно гипотезе Ньютона, жидкость как бы прилипает к стенкам канала, по которому она течет и ее скорость равна нулю. Но с увеличением расстояния от стенки, скорость
струек увеличивается. Так называемая мощность потока складывается из энергии отдельных
струек
N = ∫ dN ,
S
где N – мощность потока; dN – мощность струйки; S – площадь живого сечения потока.
Для мощности струйки можно записать:
p u2
dN = Ed Qm = (gz +
+
) ρuds,
2
ρ
где ds – площадь живого сечения струйки.
Величина удельной энергии потока равна частному от деления мощности потока на массовый расход Qm
⎛
p u2 ⎞
gz
+
∫S ⎜⎝ ρ + 2 ⎟⎠ρuds
E=
.
Qm
Это уравнение можно разбить на два интеграла
⎛ u2 ⎞
⎛
p⎞
gz
+
ρ
uds
∫S ⎜⎝ ρ ⎟⎠
∫S ⎜⎝ 2 ⎟⎠ρuds
+
E = Eп + Екин =
,
Qm
Qm
где Eп – удельная потенциальная энергия потока относительно выбранной плоскости сравнения; Екин – удельная кинетическая энергия потока.
Для вычисления Eп надо знать закон изменения давления по живому сечению. Для плавноизменяющихся течений ускорения и силы инерции незначительны, поэтому ими можно пренебречь. Экспериментально доказано, что в плавноизменяющемся потоке давления распредеp
ляются по закону гидростатистики gz + = const.
ρ
44
⎛
p⎞
ρ ⎜ gz + ⎟ ∫ υds
ρ⎠S
⎝
=
Eп =
Qm
⎛
p⎞
⎜ gz + ρ ⎟ ρQ
p
⎝
⎠
= gz + .
ρ
Qm
(2.47)
⎛ u2 ⎞
∫S ⎜⎝ 2 ⎟⎠ρuds
Для вычисления интеграла
нужно знать закон распределения скоростей по сечеQm
нию. Умножим и поделим это выражение на υ3ср S .
⎛ u2 ⎞
u3
ρ
uds
2
∫S ⎜⎝ 2 ⎟⎠
∫S 2 ds υср2 Qm
υ3ср S υср
ρυср S 3
=
=
α,
2
Qm
υср SQm
2 Qm
υ3ср S
где α – коэффициент, который учитывает неравномерность распределения скоростей в сечении,
называется коэффициент Кориолиса. Получаем выражение для удельной кинетической энергии потока:
Екин = α
2
υср
.
2
Запишем уравнение Бернулли для двух сечений потока реальной жидкости
gz1 +
2
2
p1 α1υ1 ср
p α 2 υ2ср
+
= gz2 + 2 +
+ Eпот .
ρ
ρ
2
2
(2.48)
(2.49)
Полученное уравнение позволяет сделать следующие выводы:
1. При увеличении кинетической энергии потока от одного сечения к другому потенциальная энергия уменьшается, и, наоборот, с увеличением потенциальной энергии, кинетическая
уменьшается.
2. Коэффициент α тем больше, чем больше скорости отдельных струек отличаются от величины средней скорости. Если скорости всех элементарных струек будут равны средней скорости, то α = 1.
2.2.7. Три формы представления уравнения Бернулли для потока реальной жидкости
Представленное выше уравнение Бернулли является уравнением Бернулли, записанное в
форме удельной энергии, где:
gz1 и gz2 – удельная энергия положения в сечениях 1-1 и 2-2 потока;
p1
p
и 2 – удельная энергия давления в этих сечениях;
ρ
ρ
α1υ12 ср
и
α 2 υ22ср
– удельная кинетическая энергия.
2
2
На практике при выполнении инженерных расчетов обычно применяют две другие формы
представления уравнения Бернулли.
Уравнение Бернулли в форме напоров можно получить, если разделить уравнение в форме
E
удельной энергии на g, обозначив пот = hw .
g
2
2
p1 α1υ1 ср
p2 α 2 υ2ср
z1 +
+
= z2 +
+
+ hw .
ρg
ρg
2g
2g
(2.50)
Каждый член этого выражения имеет размерность длины и может быть представлен некоторой высотой (напором) (см. рис. 2.39).
45
Здесь z1 и z2 – высота положения или геометрический напор (расстояние от плоскости сравнения до центра сечения потока);
p1 p1
p
p
=
и 2 = 2 – пьезометрическая высота или пьезометрический напор (высота подняρg γ
ρg
γ
тия жидкости в пьезометрической трубке под действием давления);
α1υ12 ср
2g
трубки).
и
α 2 υ22ср
2g
– скоростной напор (разность показаний трубки Пито и пьезометрической
2
2
p1 α1υ1 ср
p2 α 2 υ2ср
+
= H 2 – полный напор потока жид+
= H1 и z 2 +
Сумма трех высот z1 +
ρg
ρg
2g
2g
кости.
υ2
2g
p
γ
H
h′
h′′
2
3
1
1
0
d1
d2
z′
d1
3
2
l
l
0
Рис. 2.39. Опытная демонстрация уравнения Бернулли в форме напоров
Разность полных напоров двух живых сечений потока – потеря напора между этими сечениями
hw = H1 − H 2 .
(2.51)
С помощью уравнения Бернулли в форме напоров можно найти высотные отметки жидкостей, которые могут быть достигнуты в данной трубопроводной системе. Это уравнение широко используется при проектировании и гидравлических расчетах водопроводов.
Уравнение Бернулли в форме давлений получаем, если уравнение Бернулли в форме удельной энергии умножим на плотность ρ.
ρgz1 + p1 + ρ
α1υ12 ср
α1υ22ср
= ρgz2 + p2 + ρ
+ ρghw .
2
2
Здесь каждый член имеет размерность давления:
ρgz1 и ρgz2 – гравитационное давление, т. е. давление, создаваемое силой тяжести;
p1 и p2 – статическое давление;
ρ
α1υ12 ср
иρ
α 2 υ22ср
– динамическое давление;
2
ρghw – потери давления на преодоление сил трения и местные сопротивления.
2
46
(2.52)
Вывод: при увеличении скорости движения потока давление на этом участке падает и, наоборот – при уменьшении скорости давление увеличивается.
Уравнение Бернулли в форме давлений применяется для расчета систем вентиляции, газовых стояков внутри зданий и т.д.
2.2.9. Особенности турбулентного и ламинарного течения жидкости. Число Рейнольдса
Наблюдения показывают, что в жидкости возможны две формы движения: ламинарное движение и турбулентное. Проведем следующий опыт. Через стеклянную трубку будем подавать
воду. В начале трубки устанавливаем тонкую трубку, через которую подаем краску. Когда скорость движения воды в стеклянной трубке небольшая, струйка краски, вытекающая из тонкой
трубки, принимает форму нити. Это говорит о том, что отдельные частицы жидкости перемещаются прямолинейно. Жидкость в круглой трубе движется как бы концентрическими кольцевыми слоями, которые не перемешиваются между собой. Такое движение называется ламинарным (слоистым) (см. рис 2.40).
Рис. 2.40. Движение окрашенной жидкости при ламинарном и турбулентном режимах
С увеличением скорости движения в стеклянной трубке струйка краски будет размываться,
терять свою устойчивость и, при больших скоростях, краска будет равномерно окрашивать всю
массу жидкости, что указывает на интенсивное перемешивание всех слоев. Отдельные частицы
жидкости и ее небольшие объемы пребывают в состоянии хаотического и беспорядочного движения. Наряду с общими поступательными движениями имеется поперечное перемещение частиц. Такое движение называется турбулентным (см. рис. 2.40).
Эти два режима движения резко отличаются один от другого, что видно из нижеследующей
таблицы.
Таблица 2.1
Характеристика
Ламинарный режим
Турбулентный режим
Движение
Только продольное
Продольное и поперечное
Потери энергии
ghw ≡ u1
ghw ≡ u 2
Передача тепла
Теплообмен за счет теплопроводности
Теплообмен за счет теплопроводности и конвекции
Эпюра скорости
Параболическая функция
Логарифмическая функция
Коэффициент α
α=2
α ≈1
Условия перехода от ламинарного течения капельной жидкости к турбулентному в круглых
трубках впервые изучил О. Рейнольдс. Он установил, что режим зависит от трех параметров:
средней скорости υср , диаметра d и кинематической вязкости ν. Рейнальдс пришел к выводу,
что существует некоторое критическое значение соотношения этих параметров, являющееся
границей между ламинарными и турбулентными режимами течения, и нашел его:
υср ⋅ d
(2.53)
Re =
= 2320.
ν
Более точные исследования показали, что в интервале чисел Рейнальда от 2000 до 4000 происходит периодическая смена турбулентного и ламинарного режимов. Поэтому можно точно
сказать, что при Re ≤ 2000 режим движения – ламинарный, а при Re ≥ 4000 устанавливается
47
турбулентный режим. В диапазоне чисел Рейнольдса от 2000 до 4000 режим нестабильный, т.е.
может быть и ламинарным, и турбулентным.
При изучении сопротивлений, теплопередачи, явлений, связанных с переносом тепла,
транспортом твердых частиц число Рейнальда является исходным для построения расчетных
зависимостей
Подавляющее число движений жидкости в технике – турбулентные, а не ламинарные. Турбулентные течения значительно сложнее ламинарных, и для их изучения нужны другие методы.
Беспорядочный характер движения отдельных частиц жидкости в турбулентном потоке требует
применения методов статистической механики.
Хаотичность турбулентного движения с кинематической точки зрения означает, что скорость движения в отдельных точках пространства непрерывно изменяется как по величине (см.
рис. 2.41), так и по направлению. Скорость в данной точке турбулентного потока, измеренную в
данный момент времени, называют мгновенной и обозначают u, Экспериментальные исследования показывают, что изменения мгновенной скорости носит случайный характер.
u
Δu
u
u
t
t
Рис. 2.41. График изменения мгновенной скорости
Для описания турбулентного потока вводят понятия осредненной скорости, которой называют среднюю за некоторый промежуток времени скорость в данной точке
t
u=
1
udt ,
t ∫0
где t – достаточно длинный интервал времени.
При равномерном течении жидкости в трубе с постоянным расходом мгновенную скорость,
измеренную в данной точке можно разложить на три составляющие u x , u y , u z .
Каждая из составляющих скоростей изменяется со временем, но для установившегося движения за определенный промежуток времени, определенные во времени значения поперечных
составляющих равны нулю. Если ось х совпадает с осью трубы, то u x = 0, u y = 0, u z = 0 .
Если подобным способом определить осредненные скорости нескольких точек по поперек
трубы, получим эпюру осредненных скоростей по сечению трубы. Осреднение определенных
скоростей дает среднюю скорость потока υср .
Таким образом, осредненную скорость получаем после осреднения по времени мгновенных
скоростей, среднюю скорость получаем после осреднения осредненных скоростей по сечению.
Осредненную скорость можно рассматривать как скорость струйки. При неизменном расходе жидкости эпюра осредненных продольных скоростей в данном живом сечении не изменяется
с течением времени, что и является признаком установившего течения.
С помощью понятия осредненной скорости турбулентный поток с его беспорядочно движущимися массами жидкости заменяют воображаемой моделью потока, представляющей совокупность элементарных струек, скорости которых равны осредненным скоростям по величине и
по направлению. Это означает, что к турбулентному потоку можно применить представление
одномерной гидравлики.
48
Отклонение мгновенной скорости от ее осредненного значения Δu = u − u называют пульсационной скоростью или пульсацией. Замена действительных беспорядочных движений жидких комков на фиктивное струйное движение требует введения некоторых фиктивных сил
взаимодействия между воображаемыми струйками.
Благодаря этому Прандтлем был введен новый вид поверхностных сил и соответствующих
касательных напряжений
τ xy = −ρΔu x Δu y ,
которые называются турбулентными касательными напряжениями. Эти напряжения обусловлены пульсациями или обменом количества движения между соседними слоями жидкости.
Слой, движущийся с большей скоростью, подтягивает за собой отстающий и наоборот, слой,
который движется медленно, тормозит опережающий. Знак «минус» подчеркивает, что сила сопротивления имеет направление, противоположное продольной пульсации. Индексы x и y показывают направление движения слоя и поперечных пульсаций.
Осредненные касательные напряжения называются турбулентными
τ турб = −ρΔu x Δu y .
(2.54)
В схематизированном турбулентном потоке, кроме сил турбулентного обмена, вследствие
пульсации еще проявляются силы внутреннего трения. Полное касательное напряжение турбулентного потока
τ = τ турб + τвязк = −ρΔu x Δu y + μ
dux
.
dy
(2.55)
2.2.10. Уравнение Бернулли для элементарной струйки невязкой сжимаемой жидкости
Плотность сжимаемой жидкости изменяется в процессе движения. Проинтегрировав дифференциальное уравнение одномерного движения жидкости для струйки, получаем:
dp u 2
gz + ∫ + = const .
ρ 2
(2.56)
При установившемся движении влияние сжимаемости практически проявляется только в газах, при анализе течения которых удельной энергией положения можно пренебречь. Тогда:
dp u 2
∫ ρ + 2 = const .
(2.57)
Это уравнение можно назвать уравнением Бернулли для сжимаемой жидкости. Член
dp
∫ ρ характеризует потенциальную энергию газа с учетом преобразования его внутренней энергии.
Вывод: При установившемся течении невязкого газа сумма удельной потенциальной, внутренней и кинетической энергии есть величина постоянная.
dp
Для вычисления интеграла ∫
необходимо знать процесс изменения состояния газа при
ρ
этом течении. Если считать что течение происходит без теплообмена, то для невязкого газа это
будет отвечать изоэнтропическому изменению состояния. Такой процесс описывается уравнением
p
=C,
(2.58)
ρk
где k – показатель адиабаты процесса; С – постоянная.
49
1
dp
pk
вместо ρ его значение из предыдущего уравнения ρ = 1 , тогда
Подставим в интеграл ∫
ρ
Ck
dp
k p
∫ ρ = k −1 ⋅ ρ .
Теперь уравнение Бернулли для сжимаемой жидкости принимает вид
k p u2
⋅ + = const .
k −1 ρ 2
Поскольку отношение
p
= RT
ρ
k
u2
(2.59)
⋅ RT + = const .
k −1
2
Из этого уравнения следует, что изменение скорости вдоль трубки тока сжимаемого газа
связано с изменением температуры. При увеличении скорости температура падает и наоборот.
2.2.11. Уравнение Бернулли для потока вязкой сжимаемой жидкости
При составлении уравнений движения сжимаемой жидкости следует учитывать, что не
только скорости, но и плотности, температуры и давления отдельных струек в пределах живых
сечений неодинаковы, что значительно усложняет исследование. Поэтому поток конечных размеров рассматривают как одну струйку. Заменив в уравнении для струйки скорость струйки u
на среднюю скорость потока υср, можно сразу написать уравнение Бернулли сжимаемой невязкой жидкости:
2
υср
k
(2.60)
⋅ RT +
= const .
k −1
2
Теперь составим уравнение Бернулли для вязкой сжимаемой жидкости, для чего запишем
дифференциальное уравнение движения
dp
+ νd ν + dEп = 0,
ρ
интегрирование которого для сжимаемой жидкости зависит от конкретных условий движения и
закона изменения состояния газа.
При адиабатическом течении, где отсутствует обмен тепла со средой вне границ потока,
можно получить уравнение движения в конечном виде, для чего необходимо применить понятие энтальпии
dp
di = dq + ,
(2.61)
ρ
где q – количество тепла, передаваемое 1 кг газа.
Подставив уравнение энтальпии в уравнение Бернулли, получим
di − dq + νd ν + dEп + 0.
При адиабатическом течении энергия, потерянная на трение, переходит во внутреннее тепло
(dEn = dq), тогда
di + νd ν = 0.
Проинтегрировав, получим
i=
2
υср
= const.
2
Мы получили основное уравнение адиабатического течения газа.
50
(2.62)
Вывод: Сумма удельной кинетической энергии и энтальпии остается неизменной в процессе
движения газа.
Можно доказать, что для воздуха сжимаемостью можно пренебречь, если скорость течения
не превышает 70 м/с, для природного газа – 90 м/с. В системах вентиляции и газопроводов низкого давления скорости течения не превышают указанных пределов, поэтому расчет в этих системах ведется как для несжимаемой жидкости. В этих системах расчет можно вести по уравнению Бернулли в форме давлений
ρgz1 + p1 + ρ
α1υ12ср
= ρgz2 + p2 + ρ
2
α1υ2ср
+ Δpпот .
2
2
Пример применения уравнения Бернулли для расчета коротких трубопроводов
Вода перетекает из резервуара А в резервуар В по трубопроводу с диаметрами d1 = 100 мм и
d2 = 60 мм и длиной l1 = 15 м и l2 = 10 м. Необходимо определить расход воды при разности
уровней в бассейнах H = 300см. Трубопровод стальной сварной, умеренно заржавевший.
H
2
hB
d2
d1
hA
1
0
0
A
B
1
2
Рис. 2.42. К примеру расчета коротких трубопроводов
Примечание. Потерями напора пренебречь.
Ответ: Искомый расход в трубопроводе Q = 0,45 м3/с.
51
Лекция 3.
Основы моделирования гидромеханических процессов
Основные понятия: физическое и математическое моделирование; подобие; масштаб; критерии подобия; автомодельность; напряженное состояние жидкости; уравнение Навье – Стокса;
численные методы; неустранимая и вычисляемая погрешность; погрешности метода; схема
применения численных методов; одномерная модель жидкости; пьезометр; трубка Пито; датчики измерения давления; вертушка; расходомеры; манометр; вакуумметр; барометр.
Вопросы, на которые необходимо найти ответ в ходе изучения темы:
1. Почему применяется моделирование при изучении гидравлических явлений?
2. Какие виды моделирования Вы знаете?
3. Какие явления называются подобными?
4. Какие виды подобия Вы знаете?
5. Что называют масштабом модели?
6. Что называют механическим подобием?
7. Какие два потока являются геометрически подобными?
8. Какие потоки являются кинематически подобными?
9. Какие потоки являются динамически подобными?
10. Приведите пример преобразования одного масштаба в другой.
11. Какие основные силы необходимо учитывать при моделировании гидравлических явлений?
12. Что называется критерием гидравлического подобия?
13. В чем физический смысл критериев подобия Эйлера, Рейнольдса, Фруда, Архимеда?
14. Когда применяется критерий подобия Эйлера?
15. Когда применяется критерий подобия Рейнольдса?
16. Когда применяется критерий подобия Фруда?
17. Когда применяется критерий подобия Архимеда?
18. Почему при моделировании по Фруду невозможно соблюсти критерий Рейнольдса?
19. На какие вопросы должен ответить экспериментатор перед началом исследований?
20. Что делать, если при моделировании по критерию Рейнольдса на модели получают
слишком большие скорости потока?
21. Что называют автомодельностью при моделировании?
22. Объясните принцип действия приборов для измерения давления.
23. Объясните принцип действия приборов для измерения скорости.
24. Объясните принцип действия приборов для измерения расхода.
25. Почему в жидкости возникают напряжения?
26. Какие напряжения возникают на поверхности граней параллелепипеда?
27. Чему равно изменение количества движения жидкости, протекающей через неподвижный объем?
28. Что называют одномерной моделью жидкости?
29. В чем смысл уравнения Навье – Стокса?
30. Что называют численными методами?
31. Что является источником погрешности при решении задач численными методами?
32. Что называют неустранимой погрешностью?
52
33. Что называют погрешностями метода?
34. Почему не учитываются погрешности при решении задач на ЭВМ?
35. От чего зависит выбор численного метола?
36. Какие дополнительные требования предъявляются при выборе численного метода?
37. Что предшествует математическому исследованию?
38. Какие этапы метода вычислительного эксперимента вы знаете?
h1
d1
Практические задачи, решение которых может быть найдено после изучения теоретического материала:
Задача 3-1, 5
Какими будут потери напора на 1 км длины бетонного напорного трубопровода диаметром
5 × __мм, если потери на его воздушной модели (l = 1 м) при скорости движения воздуха 30 м/с
составили 1 м?
Задача 3-2, 6
При испытании на воде модели задвижки в трубе квадратного сечения ( а1 × а1 = 100 × 100
мм) перепад давления при открытии h1 = 30 мм и расходе Q1 = 8 л/с составил Δp1 = 6,4 КПа, а
сила действия потока на задвижку R1 = 48 Н. Определить Δp2 и R2 на натуре при Q2 = 10000
× Qн л/с, если h2 = 0,3, а а2 × а2 = 1 м.
Рис. 3.1. К задаче 3-2, 6
Задача 3-3, 7, 9
Найти отношение кинематических вязкостей жидкостей на натуре и на модели при одновременном соблюдении вязкостного (Reм = Reн) и гравитационного (Frм = Frн) подобия потоков,
если геометрический масштаб моделирования kl .
Задача 3-4, 8, 10
Протекание нефти (вязкость νн = 0,25 СТ) по стальному трубопроводу диаметром 500 мм исследуется на его воздушной модели. Определить скорость движения воздуха на модели ( kl =_),
если ее диаметр равен 50 мм, а скорость течения нефти в натуре составляет 1 м/с.
3.1. Основы моделирования
При изучении гидроаэродинамических явлений необходимо широкое применение эксперимента. Например, все вопросы, касающиеся турбулентного движения жидкости, не
имеют точного теоретического решения, поэтому экспериментальные решения дополняют теоретические. Все виды гидравлических сопротивлений и соответствующие им коэффициенты определяются экспериментальным путем.
Перед постановкой исследования экспериментатор должен знать: каким требованиям должна удовлетворять модель, какие величины надо измерять в опытах, какими приборами надо
пользоваться, на какие полученные величины, прежде всего, необходимо обращать внимание.
Кроме того, необходимо четко знать, что полученные результаты соответствуют явлениям, которые будут иметь место в действительности.
53
Различают математическое моделирование (на ЭВМ) и физическое (на физических моделях). Физическое моделирование проводят на моделях натурных объектов, которые просты в
изготовлении и их размеры позволяют осуществлять в лабораторных условиях эксперименты,
задаваясь различными параметрами модели и исследуемого явления, и выявлять искомые закономерности.
Обоснование моделирования и использование в натуре результатов экспериментов на модели связано с подобием движения в натуре и на модели. Подобными называют явления, происходящие в геометрически подобных системах одинаковой физической природы, когда одинаковые величины (например скорости или силы), действующие в подобных точках, имеют между
собой постоянные отношения, которые называются масштабами.
3.2. Виды подобия. Масштабы моделирования
Для установления подобия гидроаэродинамических явлений между натурой и моделью следует использовать правила механического подобия. Механическое подобие подразумевает выполнение геометрического, кинематического и динамического подобия.
Геометрически подобными являются два потока, если между их соответствующими линейными размерами существует соотношение
lн
= kl ,
(3.1)
lм
где kl – линейный масштаб, показывающий во сколько раз размеры модели изменены по сравнению с размерами натуры; lн и l м – геометрические размеры натуры и модели соответственно
(длина, ширина или высота).
S
Производными от линейного масштаба являются масштаб площадей н = kS и масштаб
Sм
V
объемов н = kV .
Vм
Кинематическими подобными являются два потока, если поля скоростей на модели и в натуре в подобных точках пространства связаны масштабом
υн
= kυ ,
(3.2)
υм
соответственно масштаб ускорений можно выразить отношением
aн
= ka .
aм
Для динамического подобия необходимо, чтобы все силы одинаковой природы, действующие в подобных точках модели и натуры на частицы жидкости, отличались между собой только
постоянными масштабами
Fн
= kF .
Fм
Любой масштаб может быть выражен через другие масштабы, например, масштаб сил можно представить следующими выражениями
k
k
k F = ka km = ν kρ kV = ν kρ kl3 = kν2 kl2 kρ .
kt
kt
Геометрически подобные системы не обязательно будут кинематически и динамически подобными. В то же время динамическое подобие подразумевает автоматически кинематическое
и геометрическое подобие.
54
3.3. Критерии подобия
Условие механического подобия требует равенства на модели и в натуре отношения всех
сил, действующих в данных системах. Однако практически невозможно создать условия подобия сил, определяющих явление. Поэтому на практике стараются соблюдать подобие основных
сил, остальными силами пренебрегают. Устанавливаемые частные условия подобия называются
критериями подобия. При установлении критериев подобия определяют условия, обеспечивающие пропорциональность силам инерции тех действующих сил, которые считаются главными в данном явлении.
Критерий Эйлера. При рассмотрении гидроаэродинамических явлений мы, прежде всего,
сталкиваемся с силами давления. Поэтому первый критерий динамического давления получают,
сравнивая масштабы сил давления и сил инерции
pн
p
= м2 .
(3.3)
2
ρ н uн ρ м u м
p
= idem , где термин idem означает, что условия,
ρ u2
определяемые данным соотношением, должны быть одинаковы на модели и на натуре.
Можно сделать вывод, что критерию Эйлера соответствует равенство отношений сил давления к силам инерции на модели и на натуре.
Критерий Рейнольдса определяет отношение сил внутреннего трения к силам инерции
uн lн ρн u м lмρ м
ulρ
= idem .
=
, Re =
(3.4)
μ
μн
μм
Обозначим это соотношение через Eu =
Критерий Рейнольдса является важнейшей характеристикой исследуемого явления, т.к. от
соотношения между силами инерции и вязкости зависят основные свойства движущейся жидкости. При соблюдении критерия Рейнольдса критерий Эйлера выполняется автоматически.
Если жидкости на модели и в натуре одинаковы ( kρ = 1 ), то при равенстве Re н = Re м
uн lн = u м l м или kl =
1
.
ku
Вывод: При моделировании по Рейнальдсу уменьшение размеров модели в kl раз требует
увеличения скорости движения жидкости на модели в такое же количество раз.
Моделирование со строгим соблюдением подобия сил вязкости встречается редко. Для многих изучаемых явлений при больших числах Re характер движения потока не зависит от изменения числа Re . Например, величина сопротивления трубы в квадратичной области сопротивления не зависит от числа Re . Этот и другие примеры доказывают, что для больших чисел Re
для ряда явлений изменение числа Рейнольдса не влияет на характер явления. Это свойство называется автомодельностью и значительно облегчает исследование. Условие автомодельности
по Рейнольдсу эквивалентно соблюдению критерия Эйлера. Исходя из вышеизложенного, моделирование по критерию Рейнольдса очень часто проводится при изучении напорных потоков.
Критерий Фруда определяет отношение сил инерции к силам тяжести
Fr =
u2
u2
u2
= idem или н = м .
gl
glн gl м
(3.5)
При моделировании по Фруду ku = kl .
Вывод: одновременно нельзя удовлетворить равенство критериев Рейнольдса и Фруда для
одной и той же жидкости на модели и в натуре.
Моделирование по Фруду проводится при изучении открытых потоков, т.к. движение открытых потоков происходит под действием силы тяжести, которая бывает значительно больше
сил вязкости.
55
Критерий Архимеда определяет отношение выталкивающей силы Архимеда к силам инерции
gl ρ − ρн 0 glм ρ м − ρ м 0
gl ρ − ρ0
Ar = 2
= idem или 2н н
= 2
,
(3.6)
u
ρ
uн
uм
ρм
ρн
где
ρ − ρ0
– отношение разности плотностей среды и струи к плотности среды.
ρ
Критерий Архимеда используется при исследовании гравитационных систем воздушного
отопления и вентиляции, где действуют архимедовы силы, возникающие вследствие разности
плотностей двух сред (течение струи в среде другой плотности).
3.4. Конечно-разностная форма уравнения Навье-Стокса
Напряженное состояние жидкости. Законы движения и покоя жидкостей и газов основываются на законах механики сплошной среды, что позволяет рассматривать равновесие и течение жидкости в целом без учета механизма молекулярного движения.
В массе жидкости, которая рассматривается как сплошная среда, под влиянием внешних сил
возникают соответствующие внутренние силы. Оценим порядок значений сил, действующих на
элементарный изолированный объем, имеющий форму параллелепипеда со сторонами dx, dy, dz
(см. рис. 3.2). Вся система движущейся массы отнесена к координатам x, y, z. На плоскостях
граней изолированного параллелепипеда возникают напряжения. Составляющие напряжений,
направленные перпендикулярно грани, называются нормальными напряжениями. Составляющие, находящиеся в плоскости граней, называются касательными напряжениями.
z′
z
σ1
τ1
τ3
τ1
σ2
τ2
τ1
τ2 τ
3
τ3
σ3
y′
τ2
x′
x
y
Рис. 3.2. Силы, действующие на элементарный изолированный объем
На поверхности граней элементарного параллелепипеда возникают три различных по величине касательных напряжения и три нормальных составляющих напряжений: вдоль осей x и y –
напряжение τ1, вдоль осей x и z – τ2 и вдоль осей y и z – τ3. Вдоль оси x действует нормальная
составляющая напряжения σ1, вдоль оси y – σ2 и вдоль оси z – σ3.
На массу жидкости, находящуюся в изолированном объеме, действуют массовые силы, которые пропорциональны третьей степени размера выделенного объема. При прохождении жидкости через изолированный элементарный объем происходит изменение количества движения
массы жидкости. Это изменение вызывает соответствующий импульс сил, действующий на
массу жидкости в изолированном объеме.
Изменение количества движения жидкости, протекающей через рассматриваемый неподвижный объем, пропорционально массе, заключенной в этом объеме и третьей степени его линейного размера. Силы, действующие на поверхности граней и равные возникающим напряжениям, умноженным на соответствующие площади, пропорциональны квадрату характерного
линейного размера. При стягивании рассматриваемого элементарного объема в точку остаются
только силы, связанные с возникающими в этой точке напряжениями.
56
Рассматривая элементарный объем, можно считать, что он находится в равновесии только
под действием сил, возникающих за счет напряжений на его поверхности.
При течении реальных жидкостей в потоке возникают напряжения, которые раскладываются на нормальные и касательные составляющие к площадкам, на которых они действуют. В таком потоке можно рассматривать две системы напряжений: нормальные напряжения (давление), определяемые в любой точке потока; дополнительные напряжения, состоявшие из трех
нормальных и трех касательных составляющих; эта система напряжений зависит в каждой точке потока от ориентации площадки, на которой возникают напряжения.
Выберем в точке, находящейся внутри потока и определяемой координатами x, y, z, систему
прямоугольных координат x′, y′, z′. В плоскостях координат возникнут кроме давления еще три
нормальные σ1, σ2, σ3 и три касательные τ1, τ2, τ3 составляющие дополнительного напряжения.
Значения дополнительного напряжения зависит от физических свойств и характера течения
жидкости. Изолированная элементарно малая масса жидкости находится в момент времени t в
начале координат. Элементарная масса имеет форму параллелепипеда с гранями, параллельными плоскостям координат. Стороны параллелепипеда имеют размеры dx, dy, dz. Рассмотрим
проекцию сил, возникающих на гранях изолированного параллелепипеда под действием дополнительных напряжений, на ось x. В момент, когда он вместе с потоком движется мимо центра
координат, на гранях, нормальных к оси x. действуют силы: σ1dydz и ⎡⎣ σ1 + ( ∂σ1 ∂x ) dx ⎤⎦ dydz .
Суммарная проекция силы, определяемая нормальными составляющими напряжения,
∂σ
⎡⎣ σ1 + ( ∂σ1 ∂x ) dx ⎤⎦ dydz − σ1dydz = 1 dxdydz.
(3.7)
∂x
На гранях, параллельных плоскости координат x0z, действуют напряжения, а, следовательно, и силы вдоль оси x: τ3 dydz и ⎡⎣ τ3 + ( ∂τ3 ∂y ) dy ⎤⎦ dxdz . Эти напряжения приводят к составляющей вдоль оси x, равной
∂τ3
dxdydz .
∂y
На гранях, параллельных плоскости координат x0y, на ось 0x проектируются силы τ 2 dxdy и
∂τ
⎡⎣ τ2 + ( ∂τ2 ∂y ) dz ⎤⎦ dxdy . Вдоль оси эти напряжения дадут составляющую 2 dxdydz .
∂y
Суммарная составляющая сил, возникающих на гранях изолированного элемента жидкости
за счет дополнительного напряжения в проекции на ось x, равна
⎛ ∂σ1 ∂τ3 ∂τ2 ⎞
+
+
(3.8)
⎜
⎟ dxdydz .
⎝ ∂x ∂y ∂z ⎠
Аналогично составляющие дополнительного напряжения, действующие на остальных гранях, в проекциях на оси 0y и 0z дадут составляющие сил:
⎛ ∂σ2 ∂τ3 ∂τ1 ⎞
⎛ ∂σ3 ∂τ1 ∂τ2 ⎞
+
+
+
+
(3.9)
⎜
⎟ dxdydz и ⎜
⎟ dxdydz .
∂x ∂z ⎠
⎝ ∂y
⎝ ∂z ∂y ∂x ⎠
Прибавляя эти силы, отнесенные к единице массы жидкости к правой части уравнений Эйлера, получим условия динамического равновесия в точке потока при течении реальной жидкости.
⎫
∂u x
∂u
∂u
∂u
1 ∂p
1 ⎛ ∂σ ∂τ ∂τ ⎞
+ ux x + u y x + uz x = −
+ X + ⎜ 1 + 3 + 2 ⎟ dxdydz;⎪
∂t
∂t
∂y
∂z
ρ ∂x
ρ ⎝ ∂x ∂y ∂z ⎠
⎪
⎪⎪
∂u y
∂u y
∂u y
∂u y
1 ∂p
1 ⎛ ∂σ ∂τ ∂τ ⎞
+ ux
+ uy
+ uz
=−
+ Y + ⎜ 2 + 3 + 1 ⎟ dxdydz; ⎬
∂t
∂t
∂y
∂z
ρ ∂y
ρ ⎝ ∂y
∂x ∂z ⎠
⎪
⎪
∂u z
∂u
∂u
∂u
1 ∂p
1 ⎛ ∂σ ∂τ ∂τ ⎞
+ ux z + u y z + uz z = −
+ Z + ⎜ 3 + 1 + 2 ⎟ dxdydz. ⎪
∂t
∂t
∂y
∂z
ρ ∂y
ρ ⎝ ∂z ∂y ∂x ⎠
⎭⎪
57
(3.10)
Определим теперь силы, возникающие в точке потока за счет вязкости. Проекция на ось 0x
сил вязкости, отнесенных к единице объема и действующих в точке, определяемой в потоке координатами x, y, z. Проекция на ось 0x сил вязкости, отнесенных к единице объема и действующих в точке, определяемой в потоке координатами x, y, z:
⎡ ∂ 2u
∂σ1 ∂τ3 ∂τ2
∂ ⎛ ∂u y ∂u x ⎞ ∂ ⎛ ∂u x ∂u z
+
+
+ μ ⎢ 2 2x + ⎜
+
+
⎟+ ⎜
∂x ∂y ∂z
∂
∂
∂
∂
∂
∂
∂y
x
y
x
y
z
z
⎝
⎝
⎠
⎣
Дифференцируя уравнение неразрывности по x, получим
∂ ⎛ ∂ux ∂u z ∂u z ⎞
∂ ⎛ ∂u x ⎞
∂ ⎛ ∂u y ∂u z
+
+
=− ⎜
+
⎜
⎟ = 0 или
⎜
⎟
∂x ⎝ ∂x ∂y ∂z ⎠
∂x ⎝ ∂x ⎠
∂x ⎝ ∂y
∂z
Тогда уравнение сил, возникающих за счет вязкости жидкости, равно
⎞⎤
⎟⎥ .
⎠⎦
⎞
⎟.
⎠
(3.11)
(3.12)
⎛ ∂ 2u ∂ 2u ∂ 2 u ⎞
μ ⎜ 2x + 2x + 2x ⎟ .
(3.13)
∂y
∂z ⎠
⎝ ∂x
Аналогично в проекции на оси 0y и 0z дополнительная проекция сил, которые следует учитывать при течении вязких жидкостей, составит:
⎛ ∂ 2u y ∂ 2u y ∂ 2u y ⎞
⎛ ∂ 2u z ∂ 2u z ∂ 2u z ⎞
μ ⎜ 2 + 2 + 2 ⎟ и μ ⎜ 2 + 2 + 2 ⎟.
⎜ ∂x
∂y
∂z ⎟⎠
∂y
∂z ⎠
⎝ ∂x
⎝
Уравнения Эйлера с учетом этих дополнительных сил примут вид:
∂u x
∂u x
∂u x
∂u x
1 ∂p
μ ⎛ ∂ 2 u x ∂ 2 u x ∂ 2u x ⎞ ⎫
+ ux
+ uy
+ uz
=−
+ X + ⎜ 2 + 2 + 2 ⎟ ;⎪
∂t
∂t
∂y
∂z
ρ ∂x
ρ ⎝ ∂x
∂y
∂z ⎠ ⎪
⎪
2
2
2
∂u y
∂u y
∂u y
∂u y
1 ∂p
μ ⎛ ∂ uy ∂ uy ∂ uy ⎞ ⎪
+ ux
+ uy
+ uz
=−
+ Y + ⎜ 2 + 2 + 2 ⎟ ;⎬
∂t
∂t
∂y
∂z
ρ ∂y
ρ ⎜⎝ ∂x
∂y
∂z ⎟⎠ ⎪
⎪
∂u z
∂u
∂u
∂u
1 ∂p
μ ⎛ ∂ 2 u ∂ 2u ∂ 2 u ⎞
+ ux z + u y z + uz z = −
+ Z + ⎜ 2z + 2z + 2z ⎟ . ⎪
∂t
∂t
∂y
∂z
ρ ∂y
ρ ⎝ ∂x
∂y
∂z ⎠ ⎭⎪
(3.14)
(3.15)
Полученная система уравнений называется уравнениями Навье-Стокса.
Уравнение Навье-Стокса, как и уравнение Зйлера, интегрируются только для некоторых частных случаев, но в последние годы, в связи с развитием различных методов решения подобных
задач, полученные уравнения используются все чаще и являются идеальным инструментом, позволяющим получить хороший результат. Наиболее часто для решения уравнений НавьеСтокса применяются численные методы.
3.5. Общая схема применения численных методов и их реализация на ЭВМ
На практике найти точное решение данного уравнения довольно сложно, так как искомое
решение обычно не выражается в привычных для нас элементарных функциях. С появлением
ЭВМ у исследователей появилась возможность решать подобные задачи численными методами.
Под численными методами подразумеваются методы решения задач, сводящиеся к арифметическим и логическим действиям над числами, т.е. к тем действиям, которые выполняет ЭВМ. В
зависимости от сложности задачи, заданной точности и т.д. может потребоваться выполнить от
нескольких десятков до многих миллиардов действий. В первом случае для получения решения
достаточно иметь калькулятор, во втором – потребуется мощная ЭВМ, особенно, если необходимо получить решение в сжатые сроки.
Решение, полученное численным методом, обычно является приближенным, т.е. содержит
некоторую погрешность. Источниками погрешности являются: 1) несоответствие математической задачи изучаемому реальному явлению; 2) погрешность исходных данных: 3) погрешность
метода решения; 4) ошибки округлений в арифметических и других действиях над числами.
58
Погрешность в решении, обусловленная первыми двумя источниками называется неустранимой. Эта погрешность может присутствовать, даже если решение поставленной математической задачи найдено точно. Вопрос о том, насколько хорошо описывает математическая модель
исследуемое явление, проверяется путем сравнения результатов. Влияние погрешности исходных данных часто удается оценить, применяя различные методы: наименьших квадратов, метод
Лагранжа и др.
Численные методы в большинстве случаев сами по себе являются приближенными. Такие
погрешности называются погрешностями метода. Это происходит потому, что численным
методом решается более простая задача, аппроксимирующая исходную задачу. В ряде случаев
используемый численный метод строится на базе бесконечного процесса, который в пределе
приводит к искомому решению. Однако процесс вычисления всегда прерывается на некотором
шаге, что дает приближенное решение.
При решении задач на ЭВМ чаще всего встречаются две ситуации:
1) если количество выполняемых арифметических действий невелико, то, обычно, ошибки
округления не проявляются, так как в ЭВМ числа представляются с 10 и более десятичными
значащими цифрами, а окончательный результат редко бывает нужен более чем с 5 десятичными значащими цифрами.
2) если задача сложная (уравнения с частными производными), то в этом случае погрешности округления в каждом действии не учитываются, так как они взаимокомпенсируются.
Для решения одной и той же задачи могут применяться различные приближенные методы, в
зависимости от требуемой точности вычисления. Численный метод может считаться удачно
выбранным, если его погрешность в несколько раз меньше неустранимой погрешности, а погрешность, возникающая за счет округлений, называемая вычисленной погрешностью, в несколько раз меньше погрешности метода. Если неустранимая погрешность отсутствует, то погрешность метода должна быть несколько меньше заданной неточности.
К численному методу предъявляется еще ряд других требований. Предпочтение отдается
методу, который реализуется с помощью меньшего числа действий, требует меньшей памяти
ЭВМ, является логически более простым, что способствует более быстрой его реализации на
ЭВМ.
Большинство численных методов основано на замене более сложных объектов, уравнений
более простыми. Наиболее удобной в обращении на практике функцией является алгебраический многочлен. Чтобы задать многочлен, нужно задать только конечное число его коэффициентов. Значения многочлена просто вычисляются, его легко продифференцировать, проинтегрировать и т.д. Алгебраические многочлены нашли широкое применение для приближения
функций. Применяются многочлены четырех видов: Тейлора, интерполяционные, равномерного приближения, наилучшего среднеквадратичного приближения (метод наименьших квадратов).
Для интегрирования дифференциальных уравнений, которые не выражаются элементарными функциями, используются различные методы: Рунге-Кутта, Монте-Карло, Эйлера, Гауса и
др. Нелинейные уравнения решаются методом итераций, деления отрезка пополам и др.
Существует большое число задач, где есть хорошо отработанные численные методы и созданные на их основе стандартные программы решения задач. Существует библиотека таких
программ. Исследователь, которому впервые встретилась единичная задача, как правило, вначале ищет стандартную похожую программу, а затем пытается внести в нее изменения, исходя
из имеющихся условий.
На первоначальном этапе исследования обычно используют более простую модель явления,
которая позволяет воспользоваться более простыми методами решения с применением стандартных программ. Затем постепенно переходят к более сложным методам и моделям, добиваясь положительного конечного результата.
Наиболее широко численные методы используются в вычислительных экспериментах – исследовании естественнонаучных проблем, средствами вычислительной математики.
59
Математическому исследованию предшествует выбор физического приближения, т.е. решение вопросов о том, какие факторы надо учесть, а какими можно пренебречь. Далее проводится
исследование проблемы методом вычислительного эксперимента, в котором выделяют несколько этапов:
1. Формируется задача, выбирается физическая модель процесса, решается вопрос о том,
какие физические величины надо учитывать. Проводится описание физической модели
математическим способом (дифференциальные, интегральные и другие уравнения). Полученную математическую модель исследуют методами математической физики, чтобы
установить, правильно ли поставлена задача, хватает ли исходных данных, не противоречат ли они друг другу, существует ли решение поставленной задачи и единственно ли
оно.
2. Построение приближенного численного метода решения задачи, т.е. выбора вычислительного алгоритма. Вычислительный алгоритм – последовательность арифметических и
логических операций, при помощи которых находится приближенное численное решение математической задачи, сформулированной на первом этапе.
3. Программирование вычислительного алгоритма на ЭВМ.
4. Проведение расчетов на ЭВМ.
5. Анализ полученных численных результатов и последующее уточнение математической
модели.
3.6. Измерительные приборы, используемые при
проведении экспериментальных работ
1. Жидкостные манометры прямого действия.
Чувствительность манометра (см. рис. 3.3) определяется по следующей формуле
Δh 1
=
s=
,
Δp gρ
т.е. чувствительность тем больше, чем меньше плотность жидкости. Давление определяется по
разности уровней жидкости в трубках А и В
p = ρg Δ h .
(3.16)
B
ΔH
A
Рис. 3.3. Жидкостный манометр
2. Механические манометры.
Манометр (см. рис. 3.4) состоит из согнутой металлической трубки Т, один конец которой
соединен с резервуаром, в котором измеряется давление. Конец трубки В соединен с рычагом
ВС, который поворачивает стрелку. При повороте стрелки она указывает величину давления.
При избыточном давлении в трубке Т свободный ее конец В начинает распрямляться и приводит в движение стрелку, которая показывает величину давления. Такие манометры отличаются
прочностью.
60
T
B
C
p
Рис. 3.4. Механический манометр
3. Барометры.
Барометры (см. рис. 3.5) используются для измерения атмосферного давления. В лабораторных условиях используется барометр Фортина, позволяющий довольно точно измерить атмосферное давление.
R
H
V
Рис. 3.5. Жидкостный барометр
Принцип действия: барометрическая трубка R опрокинута открытым концом в чашу площадью S. Прибор заполнен ртутью. Чаша прикрыта от попадания пыли тканью, что не мешает
прохождение воздуха сквозь нее. Рядом с трубкой расположена шкала, проградуированная в
мм. При помощи винта V чаша перемещается по вертикали для того, чтобы совместить свободную поверхность ртути с нулем шкалы. Над уровнем ртути в барометрической трубке сохраняется слабое давление, определяемое давлением насыщенных паров ртути. Давление столба ртути высотой H в барометрической трубке соответствует атмосферному давлению. По шкале определяем величину атмосферного давления в мм.рт.ст.
4. Вакуумметры.
Принцип действия механического и жидкостного вакуумметров аналогичен принципу действия механического манометра и жидкостного пьезометра.
5. Трубка Пито–Прандтля.
Трубка Пито–Прандтля (см. рис. 3.6) позволяет одновременно определить величину динамического и статического давления в определенной точке потока.
Через отверстие А происходит измерение динамического давления. Через отверстия М измеряется статическое давление жидкости. Жидкость под действием давления поднимается по соответствующим пьезометрическим трубкам до точек А′ и М′.
61
A′
Δh
M′
M
A
M
Рис. 3.6. Трубка Пито–Прандтля
Так как плотность газа (воздуха) значительно меньше плотности жидкости, то давлением
воздуха можно пренебречь. Разность давления в точках А′ и М′ будет Δp = ρgh . Разность давления Δp зависит от динамического давления на входе в трубку Пито–Прандтля, что следует из
уравнения Бернулли для точек А и М:
ρυ2
,
2
где υ – скорость потока на входе в трубку Пито–Прандтля. Таким образом,
pM = p A +
Δp =
ρυ2
,
2
откуда получаем
υ = 2 gh .
(3.17)
Примечания: 1. Трубка Пито–Прандтля измеряет местную скорость в данной точке, поэтому для определения расхода по формуле Q = υS необходимо измерить местную скорость в нескольких точках сечения для нахождения υ.
2. Для того, чтобы учесть потери на трение в формулу вводится коэффициент φ. Коэффициент φ определяется экспериментально для каждой трубки и вносится в паспорт измерительного
прибора. Обычно φ = 0,97-0,98.
3. Трубка Пито–Прандтля позволяет измерить довольно большие значения скоростей, при
малых скоростях увеличивается погрешность измерения из-за погрешностей манометра (± 0,1
см).
6. Расходомер Вентури.
Расходомер (см. рис. 3.7) служит для измерения расхода жидкости и представляет собой
плавную сходящуюся – расходящуюся вставку, к которой подключается дифферциальный манометр. Для вывода расчетной формулы применим уравнение Бернулли для сечения 1-1 перед
сужением и сечения 2-2 в сужении (α1 = α2 = 1).
1
2
D
d
2
1
h
Рис. 3.7. Расходомер Вентури
62
⎛ υ22 − υ12 ⎞
v22 ⎛ υ12 ⎞
Δp = ρ ⎜
⎟ = ρ ⎜1 − 2 ⎟ .
2 ⎝ υ2 ⎠
⎝ 2 ⎠
(3.18)
С учетом уравнения неразрывности υ1S1 = υ2 S 2
Δp = ρ
Откуда υ2 = 2
υ22 ⎛ d 4 ⎞
⎜1 −
⎟.
2 ⎝ D4 ⎠
(3.19)
πd 2
Δp
Δp
1
1
,
.
=
Q
2
4
4
ρ 1 − (d 4 / D 4 )
ρ 1 − (d / D )
4
Зная перепад давления по дифференциальному манометру, можно для данного диаметра
вставки и трубы определить расход жидкости, протекающей через трубу.
Для ртутного манометра
Δp = (ρрт − ρж ) gh
Обычно расходомеры выпускаются для определенных диаметров труб и его диаметр d и D
известны.
В этом случае формула упрощается
Q = Sμ
2 Δp
,
ρ
где μ – коэффициент, учитывающий конструктивные особенности прибора, в частности
(3.20)
d4
,и
D4
вносится в паспорт прибора.
7. Ротаметры.
Ротаметры (см. рис. 3.8) используются для измерения расхода жидкостей, имеющих слабые
коррозийные свойства. Ротаметр состоит из сужающейся стеклянной трубки и металлического
конусообразного измерителя. На измеритель действуют следующие силы: сила тяжести G, архимедова сила FA, сила динамического давления жидкости F. Для измерителя, находящегося в
покое, можно записать
ρ ж υ2
(3.21)
S,
2
где ρ – плотность жидкости; V – объем измерителя; υ – скорость течения жидкости; S – площадь
сечения измерителя.
Сила тяжести и архимедова сила – величины постоянные, поэтому сила динамического давления жидкости при любом расходе Q = υS будет также постоянной.
G = F + FA или mg = ρжVg +
F
FA
G
Рис. 3.8. Ротаметр
63
При уменьшении расхода измеритель опускается, уменьшается сечение для прохождения
жидкости, соответственно скорость увеличивается, сохраняя постоянным значение силы динамического давления. Если мы проградуируем стеклянную поверхность ротаметра в единицах
измерения расхода, то в зависимости от высоты поднятия измерителя можно определять расход.
64
Лекция 4.
Гидравлические сопротивления
Основные понятия: гидравлические сопротивления; сопротивление по длине; местные сопротивления; коэффициент гидравлического трения; коэффициент эквивалентной шероховатости; график Мурина; гидравлически гладкие и шероховатые трубы; движение жидкости в трубах некруглого сечения; эквивалентная длина; кавитация; истечение жидкости в атмосферу; отверстие в тонкой стенке; насадки; коэффициенты истечения; критическое истечение; истечение
при переменном напоре; истечение под уровень.
Вопросы, на которые необходимо найти ответ в ходе изучения темы:
1. Чем сопротивление по длине отличается от местного сопротивления?
2. Что называют общими потерями?
3. Привести примеры сопротивлений по длине и местных сопротивлений.
4. Почему происходит уменьшение давление при движении жидкости в трубе? Докажите.
5. Когда применяется формула Дарси-Вейсбаха?
6. От каких параметров зависит величина потерь напора по длине при ламинарном режиме?
7. Напишите уравнение равновесия сил, действующих на выделенный цилиндрический
объем жидкости, движущийся равномерно.
8. Как распределяется напряжение трения в трубе при ламинарном режиме? Докажите.
9. Как распределяются скорости в трубе при ламинарном режиме? Как это доказать?
10. Как построить эпюру скорости, если известна величина скорости на оси потока?
11. Как можно определить расход при ламинарном движении жидкости?
12. Как из формулы Пуазейля получить формулу Дарси-Вейсбаха?
13. От чего зависит коэффициент гидравлического трения при ламинарном режиме?
14. От каких факторов зависит коэффициент гидравлического трения при турбулентном течении в трубах?
15. Как выглядит график Мурина? Какие выводы можно сделать из этого графика?
16. Что называют коэффициентом эквивалентной шероховатости?
17. Чем гидравлически гладкие трубы отличаются от гидравлически шероховатых?
18. Сколько зон сопротивления имеется при турбулентном движении? Как это доказать?
19. Почему область шероховатого трения часто называют областью квадратичного сопротивления?
20. Каков порядок определения потерь напора по длине?
21. Каковы особенности определения потерь в трубопроводах некруглого сечения?
22. Что называют гидравлическим радиусом?
23. Что называется простым трубопроводом?
24. Почему при расчете длинных трубопроводов пренебрегают потерями на местных сопротивлениях?
25. Как определяют потери напора по длине в длинных трубопроводах?
26. Что называют модулем расхода?
27. Что называют требуемым напором, как его определить?
28. Как определить величину статического напора?
29. В чем отличие расчета трубопроводов при их параллельном и последовательном соединении?
65
30. Каковы особенности расчета трубопроводов при непрерывной раздаче расхода по пути?
31. Какова особенность расчета разветвления?
32. Почему возникает гидравлический удар в трубопроводах?
33. От чего зависит повышение давления при гидравлическом ударе?
34. От каких параметров зависит скорость распространения гидравлического удара?
35. Чем прямой гидравлический удар отличается от непрямого?
36. Как рассчитать повышение давления при прямом и непрямом гидравлическом ударе?
37. Как рассчитать время закрытия задвижки в магистральном трубопроводе?
38. Каков принцип действия гидравлического тарана?
39. Как избежать гидравлический удар в водопроводных сетях внутри здания?
40. От чего зависят потери напора на местных сопротивлениях?
41. Почему коэффициенты местных сопротивлений определяется почти всегда экспериментально?
42. Что называют эквивалентной длиной?
43. Когда для определения потерь напора на местном сопротивлении используется эквивалентная длина?
44. Каковы причины возникновения потерь напора при резком расширении и сужении трубы?
45. В каком случае коэффициент местного сопротивления при постепенном изменении диаметра трубы будет минимальным?
46. От каких факторов зависит коэффициент местного сопротивления при плавном повороте
трубы?
47. От чего зависит коэффициент местного сопротивления вентиля?
48. Что влияет на характер изменения коэффициента местного сопротивления от числа Re?
49. Как определить потери напора на местных сопротивлениях при ламинарном режиме?
50. Что называют кавитацией?
51. Как избежать кавитации?
52. Что называют отверстием в тонкой стенке, а что насадком?
53. Как определить скорость струи при истечении через отверстие? Как получена эта формула?
54. Как определить расход при истечении через отверстие?
55. Какие коэффициенты истечения вы знаете, как они связаны между собой?
56. Почему происходит сжатие струи при истечении?
57. Какова зависимость коэффициентов истечения от числа Re? Почему?
58. Почему длина насадка должна быть не менее трех диаметров, но не более пяти?
59. Почему при истечении через насадок расход увеличивается?
60. Почему при напорах воды около 14 м коэффициент расхода цилиндрического насадка
резко снижается?
61. Какое истечение через насадок называется безотрывным?
62. Почему расход коноидального насадка больше, чем цилиндрического?
63. Почему не рекомендуется повышать напор воды при истечении через комбинированный
насадок более 4 м?
64. В каких случаях, какие насадки применяются?
65. Почему при определении времени истечения из резервуара нельзя применить уравнение
Бернулли?
66
66. При помощи какого метода определяется время истечения из резервуара, если резервуар
имеет произвольную форму?
67. Как определить время опорожнения цилиндрического резервуара?
68. Что называют истечением при переменном напоре?
69. От чего зависит расход жидкости при истечении под уровень?
Практические задачи, решение которых может быть найдено после изучения теоретического материала:
Задача 4-1
В закрытом резервуаре поддерживается постоянное манометрическое давление pм = 0,82
атм, под действием которого по трубе диаметром d1 = d н мм и общей длине l = 25 м (расстояние
от начала до колена l1 = 5 м) вытекает жидкость при температуре t. Определить расход при напоре H1 = 2 м.
l1
H1
Рис. 4.1. К задаче 4-1
Задача 4-2
Истечение происходит из открытого резервуара при постоянном напоре H1 = 5 м по короткому трубопроводу переменного поперечного сечения с диаметрами d1 = 10 × d н мм и d2 = 100
мм в атмосферу. На втором участке трубопровода имеются два колена с плавным поворотом.
Длина первого участка l1 = 0,8 м, длина второго участка l2 = 2 м. Определить скорость истечения υ2 и расход Q2 через трубопровод.
d1
H1
d2
H2
υ1
υ2
l1
l2
Рис. 4.2. К задаче 4-2
Задача 4-3
Из резервуара А, заполненного жидкостью на высоту H1 = 5 м и находящегося под манометрическим давлением рм = 150 кПа, жидкость подается по трубопроводу длиной l = 6 м и диаметром d = 4 × d н мм в резервуар В на высоту Н = 2 м. Определить расход Q и скорость протекающей по трубопроводу жидкости.
67
M
H
H1
B
A
Рис. 4.3. К задаче 4-3
Задача 4-4
К закрытому резервуару, на свободной поверхности которого действует манометрическое
давление рм = 500 кПа, подсоединен трубопровод переменного сечения с диаметрами d1 = 2 × d н
мм и d2 = 10 × d н мм, заканчивающийся соплом диаметром dc = d1. Трубопровод подсоединен на
глубине Н1 = 8 м. На первом участке длиной l1 = 10 м установлен вентиль. Длина второго участка l2 = 5 м. Определить скорость истечения υ и расход Q вытекающей из сопла жидкости при
температуре t.
М
d2
d1
dс
H1
υ
l2
l1
Рис. 4.4. К задаче 4-4
Задача 4-5
Центробежный насос подает воду с температурой t = 15 °С по трубе диаметром d = 4 × d н
мм и длиной l = 27 м. Насос создает давление на выходе рн = 2,65 атм и подачу Q = 16 л/с. Определить напор H, создаваемый насосом.
H
Рис. 4.5. К задаче 4-5
Задача 4-6
Какую высоту h должна иметь воронка для налива жидкости в бочку из цистерны, чтобы мазут не переливался через край при полном открытии крана. Диаметр крана d1 = d н мм, напор в
цистерне Н = 3 м. Диаметр выходного отверстия воронки d2 = 2d1. Длина крана l = 1 м.
68
d1
H
h
d2
Рис.4.6. К задаче 4-6
Задача 4-7
Призматический резервуар разделен на две части перегородкой. В левом отсеке поддерживается постоянный уровень жидкости. В перегородке имеется круглое отверстие диаметром d1 =
d н мм, расположенное на глубине Н = 3,2 м. Во внешней стенке расположено отверстие диаметром d2 = 100 мм. Определить расход Q при Δh = 1,5 м.
Δh
H
d1
d
Рис.4.7. К задаче 4-7
Задача 4-8
На поршень, диаметром D = 100 мм действует сила Р = 1 kH. Определить скорость движения
поршня при диаметре отверстия d = 0, 2 × d н мм и толщине поршня а = 8 мм.
P
a
d
Рис. 4.8. К задаче 4-8
Задача 4-9
В резервуаре находится 1,1 м воды и 6,2 м нефти плотностью 900 кг/м2. Диаметр резервуара
d = 100 × d н м. Определить время слива воды через короткий патрубок диаметром 100 мм.
Задача 4-10
При исследовании истечения через круглое отверстие диаметром d0 = d н мм, получено: диаметр струи dc = 0, 8d н мм, напор Н = 2 м, время заполнения объема V = 10 л t = 32,8 с. Определить коэффициенты истечения μ, ϕ, ε. Распределение скоростей по сечению струи принять равномерным.
4.1. Виды гидравлических сопротивлений
Получение конкретных зависимостей для расчета потерь энергии при движении жидкости в
трубках и каналах является основным содержанием внутренней задачи гидравлики.
Различают два вида сопротивлений отличающиеся одно от другого по свое структуре: сопротивление по длине и местные сопротивления.
69
Потери по длине. Рассмотрим движение жидкости в горизонтальном трубопроводе постоянного сечения с неизменной эпюрой скоростей, т.е. равномерное движение. Запишем для двух
сечений уравнение Бернулли в форме давлений
ρgz1 + p1 + ρ
Здесь: z1 = z2 , ρ
α1υ12 ср
α1υ12ср
2
= ρgz2 + p2 + ρ
2
α1υ2ср
2
+ ρghw .
α 2 υ22ср
=ρ
, тогда уравнение Бернулли можно записать в следующем
2
2
виде: p1 − p2 = ρghw = ρEпот = pw = рдл , где pдл – давление, теряемое на преодоление сил трения.
Вывод: Работа сил давления расходуется на преодоление сил трения, что и обуславливает
потери механической энергии, которые прямо пропорциональны длине пути движения.
В зависимости от формы записи уравнения Бернулли эти потери называются:
• потерями давления по длине pдл ;
•
потерями удельной энергии по длине Едл ;
•
потерями напора по длине hдл .
Местные потери образуются в результате изменения структуры потока по пути движения
жидкости.
Рассмотрим движение жидкости через частично открытую задвижку в трубопроводе (см.
рис. 4.9).
1
C
2
1
C
2
Рис. 4.9. Задвижка
В отверстии (сечение С-С) скорости увеличиваются, а давление уменьшается. В сечении 2-2,
на некотором расстоянии после задвижки, скорости принимают значения, равные скоростям в
сечении 1-1 перед задвижкой.
На участках 1-С и С-2 наряду с основным течением возникает область вихревого движения.
Скорости движения частиц в этой зоне значительно меньше, чем в основном потоке, что обуславливает возникновение больших напряжений трения из-за большого градиента скорости.
Более подробно о местных гидравлических сопротивлениях будет рассмотрено далее. В зависимости от формы записи местные потери записываются как рм , Ем , или hм .
Общие потери в трубопроводе складываются из потерь по длине и потерь на местных сопротивлениях hw = hдл + hм .
4.2 Сопротивление по длине при движении в цилиндрической трубе при
ламинарном течении
При ламинарном режиме жидкость движется концентрическими слоями. Воспользуемся
формулой Ньютона для напряжений трения, приняв dy = dr
τ = −μ
70
du
.
dr
Знак «минус» указывает на то, что скорость уменьшается в направлении оси r (от центра к
стенке трубы).
Составим уравнение равномерного движения жидкости для выделенного объема длиной l и
радиусом r (см. рис 4.10).
hдл
1
2
T
P1
P2
2r0
r
T
l
2
1
Рис.4.10. Движение жидкости в прямой трубе
На выделенный объем действуют внешние силы: нормальные к живым сечениям: силы давления P1 = p1S , P2 = p2 S и касательные силы сопротивления Т, приложенные к боковой поверхности
T = S τ = 2πrl τ
Уравнение равновесия этих сил относительно направления движения:
P1 = P2 + T или p1πr 2 = p2 πr 2 + 2πrl τ ,
p1 − p2 2τ
( p − p2 )r
=
или τ = 1
.
(4.1)
l
r
2l
Вывод: При ламинарном движении в круглой трубе напряжение трения максимально у
стенки и равно 0 на оси трубы (см. рисунок).
Закон распределения скоростей по сечению трубы можно получить из следующего уравнения
( p1 − p2 )r
du
= −μ
,
dr
2l
( p − p2 )rdr
du = 1
.
2μl
( p1 − p2 )r 2
+ const .
После интегрирования, получаем: u = −
4μl
Константу находим из граничных условий: r = r0 , u = 0 .
0=−
( p1 − p2 )r0 2
( p − p2 )r0 2
+ const откуда 1
= const ,
4μl
4μl
u=−
( p1 − p2 )r 2 ( p1 − p2 )r0 2
+
.
4μl
4μl
Окончательно получаем
( p1 − p2 )(r02 − r 2 )
u=
.
4μl
(4.2)
Вывод: При ламинарном течении скорости в сечении трубки распределяются по параболическому закону (см. рис. 4.10).
71
Максимальная скорость на оси трубы будет при r = 0
umax
⎛ r2 ⎞
( p1 − p2 )r0 2
=
, или u = umax ⎜ 1 − 2 ⎟ .
4μl
⎝ r0 ⎠
(4.3)
Определим величину расхода жидкости через определенное сечение.
Расход элементарной струйки dQ = udS , где dS – площадь сечения трубки тока, dS = 2πrdr
⎛ r2 ⎞
dQ = umax ⎜1 − 2 ⎟ 2πrdr .
⎝ r0 ⎠
Полный расход
r = r0
r = r0
⎛ r2 ⎞
πr02
.
u
rdr
u
1
−
2
π
=
(4.4)
max
∫r =0
∫r =0 max ⎜⎝ r02 ⎟⎠
2
Вывод: Для того чтобы определить расход при ламинарном режиме достаточно измерить
скорость на оси потока и умножить ее на половину площади живого сечения.
Q
Определим среднюю скорость. Согласно определению υср = .
S
dQ =
dQ =
S = πr0 , Q = umax
2
πr02
.
2
Получаем
umax
.
(4.5)
2
Вывод: Средняя скорость при ламинарном режиме в два раза меньше скорости на оси потоυср =
ка.
∫ u dS
3
Коэффициент Кориолиса вычисляется из выражения α =
S
υ3ср S
.
Подставляем значения u и dS, интегрируя, получаем
α лам = 2 .
(4.6)
Для получения закона сопротивления при ламинарном режиме вернемся к формуле расхода
πr02
Q = umax
.
2
Подставим значение umax
( p1 − p2 )r0 2
=
.
4μl
Q=
( p1 − p2 )πr0 4
8μlQ
, откуда ( p1 − p2 ) =
.
8μl
πr04
Разделим на ρg, в результате получаем
( p1 − p2 ) 8μlQ 8νlQ
= 4 = 4 ,
ρg
πr0 ρg πr0 g
где ( p1 − p2 ) – потери давления;
( p1 − p2 )
= hдл – потери напора.
ρg
Получаем
hдл =
8νlQ
.
πr04 g
72
(4.7)
Вывод: Потери напора на преодоление сил сопротивления по длине при ламинарном режиме
прямопропорциональны расходу и длине трубопровода и обратнопропорциональны радиусу
трубы в четверной степени.
4.3. Формула Дарси-Вейсбаха
Заменим в полученной формуле hдл =
8νlQ
πr04 g
Q = υср
Получаем формулу Пуазейля hдл =
разом hдл = 32
32νl υср
gd 2
πd 2
d4
и r04 =
.
4
16
. Перепишем формулу Пуазейля следующим об-
υср
ν l υср
, умножим и разделим на
dd g
2
υср
2
ν l υср 2
64ν l υср
, hдл =
.
hдл = 32
υср d d 2 g
d d g υср
2
2
64 l υср
1
=
Если
, то получаем hдл =
.
Re d 2 g
υср d Re
ν
64
= λ – коэффициент гидравлического трения, то получаем формулу ДарсиRe
Вейсбаха для определения потерь напора по длине
Если принять
hдл = λ
2
l υср
.
d 2g
(4.8)
4.4. Турбулентное движение в гидравлически гладких и
шероховатых трубах
Для ламинарного течения многочисленные экспериментальные исследования подтвердили
справедливость вывода о том, что потери напора на гидравлические сопротивления зависят
только от величины скорости движения потока в первой степени. Соответственно, гидравлический коэффициент трения
64
λ=
.
(4.9)
Re
Опыты, прежде всего Г.А. Мурина с техническими трубопроводами, показали, что для турбулентного режима λ изменяется не только с изменением числа Re, но на величину λ влияет
также техническое состояние трубы. Мурин исследовал 49 труб из различных материалов,
бывших в эксплуатации, с различными диаметрами, при различных скоростях движения жидкости. Результаты опытов были получены в виде нескольких кривых (см. рис. 4.11).
Здесь четко различаются три области сопротивления при турбулентном режиме. Линия I соответствует области гидравлически гладких труб, когда величина λ зависит только от числа Re
и не зависит от материала трубы. Математическая обработка данных показывает, что для этой
области закономерна зависимость
0,3164
λ=
.
(4.10)
Re0,25
d
Эту зависимость можно использовать в диапазоне чисел Re ≤ 20 .
kэ
73
λ
d
kэ
III
II
I
Re
Рис. 4.11. График Мурина
Область II на графике является переходной областью от гидравлически гладких к шероховаd
. Для определения λ в этой области
тым трубам. Величина λ зависит как от числа Re, так и от
kэ
лучше всего подходит формула Альтшуля
⎛ 68 kэ ⎞
λ = 0,11 ⋅ ⎜
+ ⎟
⎝ Re d ⎠
0,25
(4.11)
.
Эту формулу можно использовать в диапазоне чисел 20
d
d
< Re < 500 .
kэ
kэ
Область III является областью гидравлически шероховатых труб. На графике в этой области кривые зависимости λ от Re параллельны между собой, т.е. λ не зависит от числа Re, а опреd
. В этой области λ определяется по формуле
деляется только величиной
kэ
⎛k ⎞
λ = 0,11 ⋅ ⎜ э ⎟
⎝d⎠
d
которую можно использовать при Re > 500 .
kэ
0,25
,
(4.12)
Анализ возможных значений коэффициентов гидравлического трения для различных условий показывает, что трубопроводы для систем теплогазоснабжения и вентиляции работают
преимущественно в переходной области сопротивления. Водопроводные линии чаще всего относятся к области шероховатых труб. Как гидравлически гладкие работают пластмассовые,
алюминиевые, латунные трубы.
⎛
d ⎞
Характер кривых зависимости λ = f ⎜ Re,
⎟ определяется характером обтекания потоком
kэ ⎠
⎝
жидкости в пристеночном слое выступов шероховатости, которые всегда имеются на поверхности трубы (рис. 4.12).
74
Δ
Δ
δ
δ
Рис. 4.12. Движение жидкости в гидравлически гладких и шероховатых трубах
При небольших скоростях движения жидкости частицы обтекают выступы без образования
вихрей, что объясняется малыми силами инерции. Такое обтекание потоком выступов характерно для области гидравлически гладких труб. С увеличением скорости движения, силы инерции частиц жидкости возрастают и возникают отдельные вихри за некоторыми выступами шероховатости. Количество вихрей и их величина возрастает с увеличением скорости движения
жидкости. Такая картина обтекания характерна для переходной области. При дальнейшем увеличении скорости протекании жидкости, вихри располагаются за всеми выступами, их размер
не изменяется, что характерно для области гидравлически шероховатых труб. Размер вихрей
зависит, как мы видим, от размера выступов шероховатостей, их формы, частоты их распределения по поверхности.
В качестве интегральной характеристики состояние внутренней поверхности трубы используется эквивалентная шероховатость kэ , которая определяется экспериментально на основе
гидравлических испытаний различных трубопроводов и приводится в справочниках. Приведем
некоторые значения kэ для труб из различных материалов:
Таблица 3-1
Материал и состояние трубы
kэ , мм
Из цветных металлов и стекла
≈ 0,001
Стальные бесшовные:
– новые
– бывшие в эксплуатации
≈ 0,014
≈ 0,2
Стальные сварные:
– новые
– бывшие в эксплуатации
≈ 0,05
≈ 0,5
Порядок определения потерь напора по длине
1. Определяем число Re.
2. Определяем режим движения.
При ламинарном режиме:
3. Определяем λ.
4. Подсчитываем hдл .
При турбулентном режиме:
3. Находим в справочнике kэ .
4. По значениям
d
и Re определяем область сопротивления
kэ
75
5. По соответствующей формуле определяем λ.
6. Подсчитываем hдл .
4.5. Движение жидкости в трубах некруглого сечения
Для транспорта капельных жидкостей и газов иногда используют трубопроводы некруглого
сечения: овальной, прямоугольной формы. В таких трубах возникают так называемые вторичные течения (рис. 4.13), которые можно наблюдать при подкрашивании потока.
Рис. 4.13. Движение жидкости в трубопроводах некруглого сечения
Вторичные течения происходят в плоскости поперечного сечения трубы: частицы жидкости
движутся при этом от центра трубы к углам. Накладываясь на продольные движения, вторичные течения непрерывно переносят частицы жидкости в угловые области, где наблюдаются
сравнительно высокие продольные скорости. Гидравлическое сопротивление таких труб выше,
чем сопротивление аналогичных круглых труб одинакового поперечного сечения.
При турбулентном движении жидкости коэффициенты гидравлического трения соответствуют коэффициентам для круглых труб, а увеличение сопротивления учитывается тем, что труба некруглого сечения приводится к трубе круглой соответствующего диаметра. Для этого
применяется эквивалентный диаметр d э
d э = 4 Rг ,
(4.13)
S
,
χ
(4.14)
где Rг – гидравлический радиус.
Rг =
где S – площадь сечения трубопровода; χ – смоченный периметр.
Например, для прямоугольного трубопровода со сторонами a и b:
a×b
2ab
Rг =
, dэ =
,
2(a + b)
a+b
для квадратного сечения:
a2
, dэ = a ,
4a
для равнобедренного треугольного сечения со сторонами a и высотой h
a×h
2h
Rг =
, dэ =
.
2 × 3a
3
Rг =
4.6. Местные гидравлические сопротивления
Местные сопротивления вызываются фасонными частями, арматурой, другим оборудованием трубопроводных сетей, которые изменяют величину или направление скорости движения
жидкости на отдельных участках, что всегда связано с появлением дополнительных потерь напора.
Потери напора на местных сопротивлениях определяются по формуле Вейсбаха
76
hм = ξм
υ12
,
2g
(4.15)
где ξм – коэффициент местного сопротивления, который зависит от вида сопротивления и определяется опытным путем.
Основные виды местных потерь напора можно условно разделить на следующие группы:
• потери, связанные с изменением живого сечения потока (резкое или постепенное расширение и сужение потока);
• потери, вызванные изменением направления потока, его поворотом (поворот трубы);
• потери, связанные с протеканием жидкости через арматуру различного типа (вентили,
краны, клапаны, сетки);
• потери, возникшие вследствие отделения одной части потока от другой или слияния
двух потоков (тройники, крестовины и т.д.).
Рассмотрим некоторые виды местных сопротивлений.
Резкое расширение трубопровода.
Как показывают наблюдения, поток, выходящий из узкой трубы, отрывается от стенок и
дальше движется в виде струи, отделенной от остальной жидкости поверхностью раздела (см.
рис. 4.14). На поверхности раздела возникают вихри, которые отрываются и переносятся далее
транзитным потоком. Между транзитным потоком и водоворотной зоной происходит массообмен, но он незначителен. Струя постепенно расширяется и на некотором расстоянии от начала
расширения заполняет все сечение трубы. Вследствие отрыва потока и связанного с этим вихреобразования на участке трубы между сечениями 1-1 и 2-2 наблюдаются значительные потери
напора.
1
2
d2
d1
1
2
Рис. 4.14. Резкое расширение трубопровода
Если принять ряд допущений, то теоретически можно доказать, что потери напора при резком расширении
hрр =
(υ1 − υ2 )2
– формула Борда,
2g
где υ1 и υ2 – средние скорости в трубе до расширения и после. Эту формулу можно привести к
другому виду:
2
2
⎛ υ ⎞ υ2 ⎛ S ⎞ υ2
hрр = ⎜1 − 2 ⎟ 1 = ⎜ 1 − 1 ⎟ 1 .
⎝ υ1 ⎠ 2 g ⎝ S 2 ⎠ 2 g
Если принять
77
2
⎛ S ⎞
ξ рр = ⎜1 − 1 ⎟ ,
(4.16)
⎝ S2 ⎠
коэффициент местного сопротивления при резком расширении, то формула Борда принимает
следующий вид:
hрр = ξрр
υ12
.
2g
В случае присоединения трубы к резервуару можно принять S2 = ∞ , тогда ξ рр = 1 .
Постепенное расширение.
υ1
υ2
α
d1
d2
Рис. 4.15. Постепенное расширение трубопровода
Если расширение происходит постепенно (см. рис. 4.15), то потери напора значительно
уменьшаются. При течении жидкости в диффузоре скорость потока постепенно уменьшается,
уменьшается кинетическая энергия частиц, но увеличивается градиент давления. При некоторых значениях угла расширения α частицы у стенки не могут преодолеть увеличивающееся
давление и останавливаются. При дальнейшем увеличении угла частицы жидкости могут двигаться против основного потока, как при резком расширении. Происходит отрыв основного потока от стенок и вихреобразование. Интенсивность этих явлений возрастает с увеличением угла
S
α и степенью расширения n = 2 .
S1
Потерю напора в диффузоре можно условно рассматривать как сумму потерь на трение и
расширение. При небольших углах α возрастают потери по длине, а сопротивление на расширение становится минимальным. При больших углах α наоборот возрастает сопротивление на
расширение. Коэффициент сопротивления диффузора можно определить по следующей формуле
2
λ
1 ⎞
⎛
⎛ 1⎞
1 − 2 ⎟ + k ⎜1 − ⎟ ,
⎜
(4.17)
n ⎠
⎛α⎞
⎝ n⎠
8sin ⎜ ⎟ ⎝
⎝2⎠
где k – коэффициент смягчения, который зависит от угла α, и его значения приводятся в справочниках
ξ диф =
hдиф = ξдиф
υ12
.
2g
Внезапное сужение.
При внезапном сужении потока (см. рис. 4.16) также образуются водоворотные зоны в результате отрыва от стенок основного потока, но они значительно меньше, чем при резком расширении трубы, поэтому и потери напора значительно меньше. Коэффициент местного сопротивления на внезапное сужение потока можно определить по формуле
78
υ1
υ2
d2
d1
Рис. 4.16. Внезапное сужение трубопровода
ξ рс = 0,5(1 − n) , где n =
S2
.
S1
(4.18)
В случае присоединения трубы к резервуару можно принять S1 = ∞ , тогда ξрс = 0,5 .
Постепенное сужение (конфузор).
Величина сопротивления конфузора будет зависеть от угла конусности конфузора θ. Коэффициент сопротивления можно определить по формуле
ξконф = kконф (1 − n) ,
где kконф = f (θ) , приводится в справочниках.
Поворот трубы (колено).
В результате искривления потока на вогнутой стороне внутренней поверхности трубы давление больше, чем на выпуклой. В связи с этим жидкость движется с различной скоростью, что
способствует отрыву от стенок пограничного слоя и потерям напора (см. рис. 4.17). Величина
коэффициента местного сопротивления ξкол зависит от угла поворота θ, радиуса поворота R,
формы поперечного сечения и приводится в справочниках. Для круглого сечения трубы при θ =
90º. коэффициент сопротивления можно определить по формуле
θ
R
Рис. 4.17. Плавный поворот трубопровода
ξкол = 0, 05 + 0, 2
d
.
R
(4.18)
Другие виды местных сопротивлений.
Коэффициенты местных сопротивлений для большинства сопротивлений приводятся в справочниках, их величина зависит от конструкции. Для ориентировочных расчетов можно пользоваться следующими коэффициентами местного сопротивления:
• задвижка при полном открытии – 0,15;
• вход в трубу при острых кромках – 0,5;
• вентиль с косым затвором при полном открытии (рис. 4.18) – 3;
• симметричный тройник – 1,5.
79
Рис. 4.18. Вентиль
4.7. Зависимость коэффициентов местных сопротивлений от
числа Рейнольдса. Эквивалентная длина
Приведенные данные о коэффициентах местных сопротивлений относятся к турбулентному
режиму движения с большими числами Рейнольдса, где влияние молекулярной вязкости проявляет себя незначительно. При ламинарном или близком к нему течении коэффициенты местных
сопротивлений зависят от числа Рейнольдса (см. рис. 4.19). При малых значениях Re эффект
сопротивления вызван силами вязкости и пропорционален первой степени скорости. Коэффициент сопротивления в этом случае изменяется обратно пропорционально числу Рейнольдса,
А
, где А – постоянная, зависящая от вида местного сопротивления.
т.е. ξ =
Re
ξ
Re
Рис. 4.19. Графики и зависимости коэффициента сопротивления от числа Рейнольдса
При достаточно больших числах Рейнольдса формируются отрывные течения, которые и
являются основной причиной местных сопротивлений при больших значениях Re. В этом случае ξ = const . В первом приближении можно сказать, что при резких переходах в местных сопротивлениях коэффициент ξ не зависит от Re при Re > 3000, а при плавных очертаниях – при
Re > 10000.
В общем виде для коэффициента ξ можно записать
А
ξ=
+ ξкв ,
Re
где ξкв – коэффициент сопротивления при больших числах Re, когда ξ = const. В случае линейного закона сопротивления (наклонная прямая на графике) потери напора можно определить по
эквивалентной длине.
Эквивалентная длина – такая длина прямого участка трубопровода данного диаметра, на
которой потери на трение по длине эквивалентны потери напора, вызываемой данным местным
сопротивлением, т.е. hдл = hм . Таким образом, для определения потери напора на местном сопротивлении, мы мысленно заменяем местное сопротивление прямой трубой эквивалентной
80
длины. Это позволяет нам применить формулу Дарси-Вейсбаха для определения потерь напора
на местном сопротивлении и учесть изменение числа Re.
l υср
hм = λ ⋅ э ⋅
,
d 2g
2
(4.19)
где lэ – эквивалентная длина, приводится в справочниках, зависит от диаметра трубопровода и
вида сопротивления.
4.8. Кавитация
На участках многих местных сопротивлений скорости потока резко возрастают, в результате
чего давление уменьшается. Если давление становится ниже давления насыщенных паров жидкости, возникает кавитация. Источником кавитации являются пузырьки газа и пара, которые
выделяются в сечении с пониженным давлением. Попадая в сечение с нормальным давлением,
пузырьки мгновенно исчезают под действием повышенного давления. В месте исчезновения
пузырьков давление резко увеличивается, повышается температура. Кавитация неблагоприятно
отражается на работе оборудования, т.к. возникает вибрация, шум, эрозия металла. Кавитационные свойства местных сопротивлений оцениваются по критическому значению числа кавитации.
2( p1 − pкр )
k=
,
ρυм2
где pкр = pн.п. – давление насыщенных паров жидкости.
Значение чикла кавитации для различных видов местных сопротивлений определяется экспериментально и приводится в справочниках. Предельно допустимая скорость в трубопроводе
перед местным сопротивлением определяют по формуле
υкр ≤
2 p1 − pн.п.
.
ρkкр
(4.20)
Для проверочного гидравлического расчета трубопровода на бескавитационную работу необходимо, чтобы выполнялось условие υм ≤ υкр .
4.9. Истечение жидкостей из отверстия в тонкой стенке
При истечении жидкости из какого-либо резервуара происходит процесс превращения запаса потенциальной энергии в кинетическую энергию свободной струи. Основной вопрос исследования – определение скорости истечения и расхода жидкости для различных форм отверстий
и насадков.
Рассмотрим истечение жидкости через круглое отверстие в тонкой стенке (рис. 4.20). Отверстием в тонкой стенке называется отверстие, толщина стенок которого составляет не более 1/4 диаметра. Рассмотрим случай, когда жидкость вытекает из резервуара в атмосферу (см.
рис. 4.20). Напишем уравнение энергии в форме напоров для сечения 1-1 и С-С. Обозначим
давление на поверхности жидкости в резервуаре p0, H 0 – заглубление центра отверстия под
уровнем жидкости в сосуде. Струя, вытекающая под давлением столба жидкости высотой H 0 и
разностью давлений ( p0 − pатм ) , при выходе из отверстия сжимается до сечения С-С. В сечении
С-С струйки приблизительно параллельны и движение можно считать плавноизменяющимся,
поэтому для этого сечения можно применить уравнение Бернулли. Такое сжатие обуславливается инерцией частиц жидкости, движущихся при подходе к отверстию по криволинейным траекториям. Степень сжатия струи оценивается коэффициентом
ε = Sс Sо ,
(4.21)
81
p0
1
1
const
pатм
H0
z1
dc
d0
C
0
0
C
Рис. 4.20. Истечение жидкости из отверстия в тонкой стенке
где Sс Sо – отношение площади струи в сжатом сечении к площади отверстия.
Для плоскости сравнения, проведенной относительно оси отверстия, запишем уравнение
Бернулли для движения жидкости от свободной поверхности, где скорость можно принять равной нулю, до сечения С-С
p0 pc αυc2
υc2
H0 +
=
+
+ ξo
,
ρg ρg 2 g
2g
где ξo – коэффициент сопротивления отверстия.
Выполним небольшие преобразования
υc2
p0 − pc υc2
(α + ξo ) или H =
(α + ξ o ) ,
H0 +
=
2g
2g
ρg
где H – расчетный напор
H = H0 +
p0 − pc
.
ρg
(4.22)
Скорость истечения
υc =
1
α + ς0
2 gH = ϕ 2 gH ,
(4.23)
где φ – коэффициент скорости.
Нами получена формула Торричелли для определения скорости струи, вытекающей из резервуара.
Если резервуар открыт, то ( p0 = pатм ) , тогда
H = H 0 и υc = ϕ 2 gH 0 .
(4.24)
Расход через сжатое сечение
Q = υc Sc = Sc ϕ 2 gH = S0
Sc
ϕ 2 gH = S0εϕ 2 gH ,
S0
где ε – коэффициент сжатия струи; S0 – площадь сечения отверстия.
Окончательно получаем:
Q = μS0 2 gH ,
где μ – коэффициент расхода.
82
(4.25)
4.10. Зависимость коэффициентов истечения
от числа Рейнольдса
При истечении вязких жидкостей (например, дизельного топлива через форсунки) или при
истечении с небольшими скоростями маловязких жидкостей (при малых числах Рейнольдса)
будет проявляться зависимость величин коэффициентов истечения μ, φ, ε от Re (см. рис. 4.21).
μ, φ, ε
φ
ε
μ
Re
Рис. 4.21. Зависимость коэффициентов истечения от числа Рейнольдса
При истечении через малое отверстие в тонкой стенке коэффициент скорости φ с увеличением Re возрастает, что связано с уменьшение сил вязкости, что в свою очередь сказывается на
уменьшении коэффициента сопротивления ξ. Коэффициент сжатия уменьшается вследствие
увеличения радиусов кривизны поверхности струи на её участке от кромки до сжатого сечения
С-С (см. рис. 4.22). При Re → ∞ значения коэффициентов φ и ε приближаются к значениям, соответствующим истечению идеальной жидкости (φ = 1, ε = 0,6).
H0
R
C
R
H0
C
C
C
R
R
Рис. 4.22. Истечение жидкости из отверстия
Изменение величины коэффициента расхода μ определяется его зависимостью от коэффициентов φ и ε:
μ = ϕ⋅ε .
(4.26)
Зная характер изменения коэффициентов μ, φ, ε от числа Re при истечении через отверстия
и насадки, можно с большей точностью определить скорость υc = ϕ 2 gH 0 , расход Q и другие
параметры потока.
При больших числах Re (турбулентный режим) коэффициенты истечения постоянны, зависят только от вида отверстия, определяются опытным путем и приводятся в справочниках.
4.11. Истечение из насадков
Насадком называется короткая труба длиной от 3 до 5 его диаметров, присоединенная к отверстию. При расчете насадков потерями напора по длине обычно пренебрегают.
83
Рассмотрим процесс истечения жидкости на примере внешнего цилиндрического насадка
(см. рис. 4.23). На входе в насадок образуется водоворотная зона, которая является источником
потерь напора, поэтому коэффициент скорости насадка меньше, чем круглого отверстия.
H0
C
B
C
B
Рис. 4.23. Истечение жидкости из насадка
После сжатия в сечении С-С струя расширяется до сечения насадка и из насадка выходит
полным сечением. В сжатом сечении С-С скорость потока выше, чем на выходе, а значит, давление в этом сечении меньше, чем при истечении из отверстия, когда давление равно атмосферному. Таким образом, в насадке создается вакуум и эффект «подсасывания», что увеличивает расход через насадок. Если резервуар открыт, а истечение происходит в атмосферу, то расход можно определить по той же формуле, что и для отверстия.
Q = μS н 2 gH 0 ,
где μ – коэффициент расхода насадка; Sн – площадь сечения насадка.
Истечение через насадок, при котором струя полностью заполняет все пространство насадка, называется истечением без отрыва.
Найдем минимальное давление внутри насадка. Истечение происходит в воздух с давлением
p
p
ратм . Давление на свободной поверхности жидкости – р 0 . Расчетный напор – H = 0 − с + H 0 .
ρg ρg
В сжатом сечении струи рс < ратм . Разность давления ратм – рс растет пропорционально напору
Н, т. к. скорость в сжатом сечении увеличивается с увеличением Н. Составив уравнение Бернулли для сечений В-В и С-С, можно доказать, что максимальное значение вакуума для воды
наступает при ратм − 0, 75 Н ρg . При некотором значении Н, которое называется критическим,
давление внутри насадка становится равным 0. Т.к. при дальнейшем увеличении напора давление не может уменьшаться далее, то происходит отрыв струи от стенок из-за того, что наружный воздух прорывается внутрь насадка, и насадок начинает работать как отверстие в тонкой
стенке. Такое истечение называется истечением с отрывом.
При истечении с отрывом расход насадка резко уменьшается, т.к. исчезает эффект «подсасывания».
4.12. Виды насадков
Другие виды насадков применяются для того, чтобы увеличить скорость вытекающей струи
или расход (см. рис. 4.24).
Конический сходящийся насадок применяется для увеличения скорости истечения струи,
ее кинетической энергии. Коэффициенты истечения насадка зависят от угла сужения. При угле
сужения β = 14º, μ = φ = 0,95.
Конический расходящийся насадок применяется для увеличения расхода жидкости, т.к. отверстие на выходе из насадка больше, чем на входе. Коэффициенты истечения насадка также
зависят от угла расширения. При угле расширения α = 6º, μ = φ = 0,47. Эти насадки работают
при небольших напорах, т.к. при увеличении напора больше 3 м может быть отрыв струи от
стенок насадка.
84
d
d
d
d
Рис. 4.24. Виды насадков
В коноидальном насадке вход изготавливается по форме естественно сжимаемой струи, что
обеспечивает уменьшение зоны отрыва и уменьшение сопротивления насадка. Коэффициент
расхода μ = φ = 0,98.
Комбинированный насадок представляет собой комбинацию коноидального и конического
расходящегося насадка. Давление в насадке искусственно снижается, что увеличивает расход
через насадок. Использовать такой насадок можно лишь при небольших напорах от 1 до 4 метров, т.к. при больших напорах в насадке возникает кавитация, в результате чего увеличивается
сопротивление насадка. Длительная кавитация приводит к разрушению насадка.
4.13. Истечение при переменном напоре и под уровень жидкости
Истечение при переменном напоре. При истечении жидкости из резервуаров, бассейнов
очень важно знать время его полного опорожнения. Движение жидкости в этом случае неустановившееся (см. рис. 4.25). За бесконечно малый промежуток времени dt, в течение которого
уровень в сосуде опустится на величину dh, течение можно считать установившимся. За это
время из отверстия вытекает объем жидкости
S
dh
H
h
Sо
Q
Рис. 4.25. Истечение при переменном напоре
dV = −Qdt .
С другой стороны, этот объем можно представить в виде
dV = Sdh ,
где S – площадь свободной поверхности жидкости в резервуаре в момент времени dt.
Приравняв два выражения и подставив значение Q = So μ 2 ghdt , получим
Sdh = So μ 2 gh dt ,
где S – площадь резервуара на уровне h; So – площадь отверстия, через которое сливается жидкость.
85
Время опорожнения резервуара высотой Н
t=−
h=0
1
μSo 2 g
∫
S
h= H
dh
.
h
Для определения времени полного опорожнения резервуара необходимо проинтегрировать
это выражение от h = H до h = 0. Для резервуара с переменной площадью сечения S ≠ const это
сделать трудно, необходимо использовать метод конечных разностей.
Для цилиндрического резервуара, когда S = const
t=−
1
μSo
0
S
2g ∫
H
dh
2 SH
.
=
h μSo 2 gH
(4.27)
Здесь числитель равен удвоенному объему резервуара, а знаменатель представляет расход в
начальный момент опорожнения, т.е. при напоре Н. Таким образом, время полного опорожнения резервуара в два раза больше времени истечения того же объема жидкости при постоянном
напоре.
Истечение под уровень – истечение в пространство, заполненное такой же жидкостью (см.
рис. 4.26). Структура потока при таком истечении не изменяется. Расчетный напор при истечении под уровень представляет собой разность гидростатических напоров по обе стороны стенки, т.е. скорость и расход не зависят от высоты расположения отверстия. Если оба резервуара
открыты, то давление слева от отверстия будет p1 = ρgz1 , справа от отверстия – p2 = ρgz2 . Расход через отверстие будет
Δz
z1
z2
d
Рис. 4.26. Истечение под уровень
Q = S o μ 2 g Δz ,
(4.28)
где Δz – перепад уровней жидкости в резервуарах.
Пример применения уравнения Бернулли для расчета расхода через диафрагму
Нормальная диафрагма – отверстие с достаточно острой кромкой для измерения расхода
жидкости или газа в трубе (рис. 4.27). Измерение проводится при помощи определения разности давлений до диафрагмы и в сечении с минимальным диаметром. Необходимо определить
расход газа через нормальную диафрагму, если ее диаметр d = 5 мм, диаметр трубы D = 100 мм,
перепад давления по дифференциальному манометру: h = 10 мм. рт.ст.
D
d
h
Рис. 4.27. К примеру определения расхода при истечении
Ответ: Расход через диафрагму Q = 0, 0062 м3/с.
86
Лекция 5.
Практическое применение законов гидравлики
Основные понятия: короткий трубопровод; полный напор; пьезометрический напор; магистральный трубопровод; последовательное и параллельное соединение простых трубопроводов;
транзитный расход; непрерывная раздача расхода по пути; гидравлический удар; гидравлический таран.
Вопросы, на которые необходимо найти ответ в ходе изучения темы:
1. Что является целью расчета короткого трубопровода?
2. Какой трубопровод считается коротким?
3. Какие типовые задачи расчета короткого трубопровода Вы знаете?
4. Какое уравнение используется при расчете коротких трубопроводов?
5. Какие требования необходимо выполнить при выборе расчетных сечений?
6. Что называют плоскостью сравнения?
7. Что необходимо определить в начале расчета короткого трубопровода?
8. Какой режим движения жидкости принимается, если неизвестен расход?
9. Как определяются потери напора если трубопровод имеет участки с разными диаметрами?
10. Приведите пример короткого трубопровода?
11. Как определить коэффициент гидравлического трения если расход неизвестен?
12. Какую скорость используют при расчете местного сопротивления?
13. Как определяются полный и пьезометрический напоры?
14. Какие сечения необходимо выбрать при построении линий полного и пьезометрического
напоров?
15. Каков порядок построения линий полного и пьезометрического напоров?
16. Каков порядок определения потерь напора по длине?
17. Что называется простым трубопроводом?
18. Почему при расчете длинных трубопроводов пренебрегают потерями на местных сопротивлениях?
19. Как определяют потери напора по длине в длинных трубопроводах?
20. Что называют модулем расхода?
21. Что называют требуемым напором, как его определить?
22. Как определить величину статического напора?
23. В чем отличие расчета трубопроводов при их параллельном и последовательном соединении?
24. Каковы особенности расчета трубопроводов при непрерывной раздаче расхода по пути?
25. Какова особенность расчета разветвления?
26. Почему возникает гидравлический удар в трубопроводах?
27. От чего зависит повышение давления при гидравлическом ударе?
28. От каких параметров зависит скорость распространения гидравлического удара?
29. Чем прямой гидравлический удар отличается от непрямого?
30. Как рассчитать повышение давления при прямом и непрямом гидравлическом ударе?
31. Как рассчитать время закрытия задвижки в магистральном трубопроводе?
32. Каков принцип действия гидравлического тарана?
87
33. Как избежать гидравлический удар в водопроводных сетях внутри здания?
34. Почему при расчете длинных трубопроводов пренебрегают потерями на местных сопротивлениях?
Практические задачи:
Задача 5-1, 6
Определить расходы в ветвях 1, 2, 3 и повышение давления на участке ВС при гидравлическом ударе в случае мгновенного закрытия задвижки в точке С. Параметры участков: l1 = ll, d1 =
d; l2 = 3/4l, d2 = 3/2d; l3 = 4/3l, d3 = d/2; lBC = 2l, dBC = 2d. Напоры: в точке А НА = 40 м, в точке В
НВ = 5 м. Трубы стальные, новые, толщиной δ = 6 мм. Длина труб 1 = 100 × lн м, диаметр d = 5
× d н мм.
1
Q
2
B
3
A
C
Рис. 5.1. К задаче 5-1, 6
Задача 5-2, 7
От напорного бака идет магистральный трубопровод длиной l1 = 1,5l. В точке В магистраль
разветвляется на две ветви: ВС с расходом QС = 1,5Q и BD с расходом QD = Q = 50 × Qн л/с.
Длины участков: l2 = ll, l3 = 1,1l (l = 100 × lн м, d = 5× d н мм). Отметки земли: C = D. Определить
напор в напорном баке, величину давления на участке АВ при мгновенном закрытии задвижки в
точке В, если толщина стенки трубы δ = 6 мм.
l1
A
C
l2
l3
B
D
Рис. 5.2. К задаче 5-2, 7
Задача 5-3,8
Какое давление должен создать насос в точке А для обеспечения расхода непрерывной раздачи по пути q = 0,05 л/с на 1 п.м, если l AB = 120 м, d AB = 5 × d н мм, а длина и диаметр на участке ВС в два раза меньше, чем на участке ВС? Давление в точке С не должно быть меньше 2
атм. Насколько повысится давление при гидравлическом ударе на участке АВ, если толщина
стенки трубы δ = 4 мм?
A
B
C
q
Рис. 5.3. К задаче 5-3, 8
88
Задача 5-4, 9
Определить отметку водонапорной башни H A , если сеть состоит из новых чугунных труб.
Параметры участков: l1 = 2l, d1 = 1,25d; l2 = l, d 2 = d; l3 = 1,5l, d3 = d; длина магистрального
трубопровода l = 2,5l, диаметр d = 20. Свободный напор Нсв = 10 м. Расход Q = 50 × Qн л/с. Насколько повысится давление на участке А-2 при мгновенном закрытии задвижки в точке 2,если
l A2 = 1/2l Труба толщиной δ = 4 мм (l = 200 м, d =5 × d н мм).
HA
1
Q
2
A
B
3
C
Рис. 5.4. К задаче 5-4, 9
Задача 5-5,10
Из резервуара А вода подается в разветвленную сеть. На одном из участков имеется путевой
объемный расход воды q = 0,03 л/с. Трубы чугунные. Определить распределение расхода в ветвях трубопровода с расходом Q1 = 1,6 Q2 л/с, на параллельных участках. На какой высоте должен находиться резервуар А для самотечной подачи воды в пункты В и С, если расход Q2 = 50
× Qн л/с (l = 400 м, d = 5 × d н мм).
A
l, d
l, d/2
l, d/2
Q1
q
Q2
Рис.5.5. к задаче 5-5, 10
5.1. Расчет короткого трубопровода
Коротким считается трубопровод, в котором потери на местных сопротивлениях превышают
5% от общих потерь. К коротким относятся системы водоснабжения, водоотведения, горячего
водоснабжения внутри зданий (см. рис. 5.6 и 5.7).
Целью расчета короткого трубопровода может быть определение напора или давления в начале трубопровода, потерь напора или потерь давления, а также определение расхода или диаметра трубопровода при известном напоре в его начале.
При расчете напора в начале короткого трубопровода должны быть заданы:
• напор или давление в конце трубопровода;
• расход
• диаметр и длина трубопровода.
При расчете напора в конце короткого трубопровода должны быть заданы:
• напор или давление в начале трубопровода;
• расход;
• диаметр и длина трубопровода.
89
Рис. 5.6. Схема водоснабжения этажа
Рис. 5.7. Схема водоснабжения коттеджей
Если целью расчета является определение расхода трубопровода, то должны быть заданы:
• напор или давление в начале и в конце трубопровода;
• диаметр и длина трубопровода.
В тех случаях, когда необходимо подобрать диаметр трубопровода должны быть известны
напоры (давление), расход и длина трубопровода. В этом случае задача решается подбором
диаметра, соответствующего установленным условиям подачи жидкости.
При расчете коротких трубопроводов применяется уравнение Бернулли для двух выбранных
сечений и уравнение неразрывности. Обычно расчетные сечения выбираются в начале и в конце трубопровода. Обязательным условием является то, что движение на расчетном участке трубопровода должно быть установившемся, а в расчетных сечениях в начале и в конце расчетных
участков трубопровода – равномерным.
Уравнение Бернулли составляется относительно плоскости сравнения, которой может быть
любая горизонтальная плоскость, в том числе линия горизонта, поверхность земли и т.д. Плоскость сравнения рекомендуется проводить с таким расчетом, чтобы упростить решение уравнения Бернулли. Обычно плоскость сравнения проводится или по центру трубы в начале или в
конце трубопровода, или по уровню жидкости в резервуаре, из которого происходит истечение
жидкости или куда она подается.
90
Расчет необходимо начать с определения формы записи уравнения Бернулли. Для расчета
систем водоснабжения и водоотведения обычно применяется уравнение Бернулли в форме напоров:
2
2
p1 α1υ1 ср
p2 α 2 υ2ср
z1 +
+
= z2 +
+
+ hw .
ρg
ρg
2g
2g
(5.1)
Если известны расход и диаметр трубопровода, то вначале определяется режим движения
жидкости для всех участков с различными средними скоростями движения жидкости. Режим
движения жидкости считаем турбулентным, если число Рейнольдса больше или равно 2320, и
ламинарным, если Re<2320. Если расход неизвестен, то в первом приближении для маловязких
жидкостей (вода), режим движения принимается турбулентным, для вязких (густое индустриальное масло) – ламинарным. После определения расхода в первом приближении определяется
число Рейнольдса и по нему уточняется режим движения.
Расчет трубопровода фактически сводится к расчету потерь – местных и по длине. В коротком трубопроводе потери напора по длине определяются по формуле:
2
l υср
hдл = λ
.
d 2g
(5.2)
Если трубопровод имеет участки с разными диаметрами, то потери по длине определяются
для каждого участка отдельно, а затем складываются. Коэффициент гидравлического трения λ
определяется по формулам, которые выбираются в зависимости от режима движения и области
сопротивления. При неизвестном расходе:
• для ламинарного режима коэффициент гидравлического трения определяется по формуле:
64
λ=
,
(5.3)
Re
а затем уточняется;
• для турбулентного режима коэффициент гидравлического трения рассчитывается по
формуле:
0,25
⎛k ⎞
(5.4)
λ = 0,11 ⋅ ⎜ э ⎟ .
⎝d⎠
а затем уточняется. Коэффициент эквивалентной шероховатости трубы определяется по таблицам в зависимости от материала трубы и времени ее эксплуатации.
Потери напора на местном сопротивлении при турбулентном режиме определяются по формуле:
hм = ξ
2
υср
2g
,
(5.5)
где υср – средняя скорость на участке после сопротивления.
Если режим движения ламинарный, то потери напора на местном сопротивлении определяются по формуле:
2
l υср
hм = λ
,
d 2g
(5.6)
где l – эквивалентная длина, определяемая по таблицам в зависимости от вида местного сопротивления.
После определения потерь напора на каждом местном сопротивлении, они складываются.
Считается, что эти потери происходят в данном сечении трубопровода.
После расчета потерь напора по длине и на местных сопротивлениях строятся линии полного и пьезометрического напоров. Для этого трубопровод разбивается на сечения. Количество
91
сечений зависит от количества прямолинейных участков трубопровода и местных сопротивлений. Полный напор в начальном сечении трубопровода определяется по формуле:
H1 = z1 +
2
p1 α1υ1 ср
+
.
ρg
2g
(5.7)
В сечении с местным сопротивлением на графике указывается падение напора на величину
местного сопротивления hм . Обычно в этом месте на графике наблюдается резкое падение величины напора. Величина пьезометрического напора будет меньше полного напора на величину скоростного напора:
H p1 = H1 −
α1υ12ср
2g
.
(5.8)
Следующее сечение выбирается в конце прямого участка трубы перед или в месте расположения следующего сопротивления. Полный напор в этом сечении будет меньше полного напора
в первом сечении на величину потерь напора по длине:
H 2 = H1 − hдл1− 2 ,
(5.9)
где hдл1− 2 – потери по длине между первым и вторым сечениями.
Пьезометрический напор во втором сечении будет меньше полного напора на величину скоростного напора:
H p 2 = H 2 − α2
2
α 2 υ2ср
2g
.
(5.10)
Для упрощения расчета необходимо помнить, что при равномерном движении жидкости линия пьезометрического напора параллельна линии полного напора. Таким образом, последовательно строятся линии полного и пьезометрического напоров по всей длине трубопровода. Благодаря построенным линиям напоров очень легко установить сечение трубопровода, в котором
происходят наибольшие потери напора, что позволяет в случае необходимости предусмотреть
мероприятия по их минимизации. На этом же графике в принятом масштабе наносится линия
геометрического напора, которая проходит по оси трубопровода.
5.2. Расчет длинных трубопроводов
5.2.1. Понятие о простом и сложном напорных трубопроводах
Трубопроводы служат для транспортирования различных жидкостей на различные расстояния. Гидравлический расчет трубопроводов базируется на основных уравнениях гидравлики.
При расчете длинных трубопроводов пренебрегают потерями напора на местных сопротивлениях, которые малы и обычно не превышают 5 % от общих потерь.
Преобразуем формулу Дарси, заменив скорость расходом, поделенным на площадь поперечного сечения трубы
2
l υср
l
hдл = λ
=λ
d 2g
d
где A =
Q2
⎛ πd 2 ⎞
2g ⎜
⎟
⎝ 4 ⎠
2
= AQ 2l ,
(5.11)
8λ
– удельное сопротивление трубопровода.
g π2 d 5
Для квадратичной области сопротивления А зависит только от диаметра трубопровода и от
его шероховатости. Следовательно, значения А можно определить опытным путем для трубопроводов с различной степенью шероховатости и с разными диаметрами.
Величина К =
1
называется модулем расхода и также приводится в справочниках. Тогда
А
92
2
Q
(5.12)
hдл = 2 l .
K
Введение понятия удельного сопротивления трубопровода упрощает расчет, т.к. значения А
приводятся в справочниках в зависимости от диаметра трубы и ее шероховатости.
Запишем уравнение Бернулли для двух сечений трубопровода в его начале и в конце:
2
2
p1 α1υ1 ср
p2 α 2 υ2ср
z1 +
+
= z2 +
+
+ hдл и обозначим:
ρg
ρg
2g
2g
p1
= H тр – требуемый напор, т.е. напор, который должен создать насос в начале трубопроρg
вода;
z2 − z1 = Δz – разность отметок земли в конце и в начале трубопровода;
p2
= H п – пьезометрический напор, т.е. напор в конце трубопровода, который задается при
ρg
проектировании;
H п + Δz = H ст – статический напор.
Принимая во внимание, что в трубопроводе постоянного диаметра
α1υ12 ср
2g
=
α 2 υ22ср
2g
,
тогда уравнение Бернулли примет вид
Q2
(5.13)
l.
K2
Нами получена формула для гидравлического расчета простых, длинных трубопроводов.
Простым называется трубопровод, не имеющий ответвлений. Всякие другие трубопроводы
относят к категории сложных.
H тр = Н ст +
5.2.2. Расчет трубопроводов, соединенных последовательно и параллельно
При последовательном соединении простых трубопроводов разной длины и с различными
диаметрами стык в стык, трубопровод представляет собой простой трубопровод, который можно разделить на несколько участков (см. рис. 5.8). Расчет такого трубопровода не представляет
труда
d1
d2
d3
l1
l2
Q
l3
Рис. 5.8. Последовательное соединение простых трубопроводов
⎛ l
l ⎞
l
H тр = Н ст + Q 2 ⎜ 12 + 22 + 32 ⎟ .
(5.14)
⎝ K1 K 2 K 3 ⎠
При параллельном соединении (см. рис. 5.9) пьезометрический напор в узловых точках А и
В одинаков для всех участков. Расход Q основного трубопровода до деления и после объединения труб один и тот же. Задача расчета состоит в том, чтобы определить расходы в отдельных
ветвях системы Q1 , Q2 , Q3 ,... и также потери напора hдл между точками А и В. Общий расход,
диаметры и длины труб предполагаются известными.
93
Q1 l1 d1
Q
Q
Q2 l2 d2
A
B
Q3 l3 d3
Рис. 5.9. Параллельное соединение простых трубопроводов
Потери напора в любой трубе ответвления одинаковы, так как в обеих общих точках разветвления имеется один и тот же напор Н А и Н В , т.е.
hдл = hдл1 = hдл 2 = hдл3 ,
hдл1 = A1Q12l1 , hдл2 = A2Q22l2 , hдл3 = A3Q32l3 .
Отсюда
A2l2 Q1
,
=
A1l1 Q3
Q1
=
Q2
A3l3 Q2
,
=
A1l1 Q3
A3l3
,
A2l2
(5.15)
т.е. расходы на участках распределяются обратно пропорционально корню квадратному из их
сопротивлений. Кроме того Q = Q1 + Q2 + Q3 . Совместное решение этих уравнений дает возможность найти расходы на участках при заданных их размерах и общем расходе.
Простое разветвление представляет собой схему так называемой «вилки» (см. рис. 5.10.).
В отличие от параллельного соединения напоры в конечных точках С и D могут быть не одинаковы. Расход до разветвления и после соединения равен сумме расходов отдельных участков
Q = Q1 + Q2 .
Q1 l1 d1
Q
A
B
С
Q2 l2 d2
D
Рис. 5.10. Простое разветвление трубопроводов
Для определения расходов в ветвях при заданном общем расходе и размерах труб ветвей необходимо составить уравнения hдлАВ = АQ 2l , h1 = А1Q12l1 , h2 = А2Q2 2l2 и решить их совместно с
уравнением общего расхода.
Непрерывная раздача расхода по пути. На схеме (рис. 5.11) представлен участок трубопровода, на котором расход равномерно разбирается по пути. К такой схеме приближается работа магистрали водопроводной сети вдоль улицы города при правильной планировке поперечных улиц. Возьмем элементарный участок dx на расстоянии х от начала трубопровода. Разбор
расхода на 1 м длины – q. Расход через элементарный участок
Qx = Q − qx .
l
x
dx
Рис. 5.11. Непрерывная раздача расхода по пути
94
Потеря напора на этом участке
dhдл = АQx 2 dx = A(Q − qx) 2 dx .
Потеря напора по всей длине трубопровода
l
l
⎛
⎛ 2
q 2l 3 ⎞
Q 2l ⎞
Q 2l
2
.
=
A
Q
l
−
Q
l
+
=
A
hдл = ∫ dhдл = А∫ (Q − qx) 2 dx = A ⎜ Q 2l − Qql 2 +
=
⎜
⎟
⎟
3
3
3
⎝
⎠
⎝
⎠
0
0
Вывод: При равномерной раздаче потеря энергии в трубопроводе в три раза меньше, чем в
случае транзитного расхода.
На практике обычно встречается смешанный случай, когда часть расхода Qтр проходит тран-
зитом, а другая Qр отбирается вдоль пути (см. рис. 5.12). В этом случае потеря напора опреде2
ляется по формуле hдл = АQрасч
l , где расчетный расход Qрасч определяется по формуле
Qрасч = Qт2 + Qp Qт +
l
Qp2
3
(5.16)
.
Qт
Qт+Qр
Qр
Рис. 5.12. Транзитный расход
5.2.3. Гидравлический удар
Под гидравлическим ударом понимают резкое повышение давления жидкости в трубопроводе, вызванное внезапным изменением скорости течения. Явление гидравлического удара
свойственно только капельным жидкостям, которые почти не деформируются. Гидравлический
удар в водопроводных линиях возникает при быстром закрытии или открытии запорной арматуры. Повышение давления при гидравлическом ударе иногда приводит к разрыву стенок трубы. Физически явление объясняется инерционными усилиями, возникающими в жидкости при
резком изменении скорости движения.
Рассмотрим гидравлический удар на примере простейшей схемы (см. рис. 5.13). Пусть в резервуаре напор воды будет постоянным независимо от изменения скорости течения в трубе.
При полностью открытом кране В в трубопроводе устанавливается скорость движения жидкости υ0 . При быстром закрытии крана жидкость в непосредственной близости от него остановится. Под действием напора движущихся по инерции частиц, давление в этой части трубопровода повысится, что приведет к расширению стенок трубопровода. Переход от движения к покою и повышение давления происходит по всей длине жидкости не мгновенно, а через некоторый промежуток времени, что объясняется тем, что жидкость не является абсолютно несжимаемой, а стенки трубы немного, но деформируются. Движение в трубопроводе при гидравлическом ударе относится к категории неустановившегося, поэтому уравнение Бернулли в данном
случае неприменимо.
υ0
А
В
l
Рис. 5.13. Движение жидкости в магистральном трубопроводе до гидравлического удара
95
Теоретическое обоснование явления гидравлического удара и метод его расчета впервые дал
Н.Е. Жуковский в 1898 г. Жуковский предложил формулу для определения повышения давления, применив закон сохранения количества движения
Δpуд = ρcυ0 ,
(5.18)
где с – скорость распространения гидравлического удара вдоль трубы от крана до резервуара.
Скорость распространения гидравлического удара можно найти, применив закон сохранения
массы с учетом уравнений механики упругих тел
1
с=
,
1
2ρr
(5.19)
+
Eж δEтр
где Еж – модуль упругости жидкости; Етр – модуль упругости трубы; r – радиус трубы; δ –
толщина стенки трубы.
Повышение давления в трубопроводе будет гораздо меньше, если задвижку закрывать не
мгновенно, а постепенно. В этом случае ударная волна успевает достигнуть резервуара, отразиться от него и вернуться к не полностью закрытому крану. Такой гидравлический удар называют непрямым. Повышение давления при непрямом гидравлическом ударе может быть оценено приблизительно, если считать, что его сила уменьшается пропорционально увеличению вре2l
мени закрытия крана В по сравнению с фазой удара, которая рассчитывается по формуле t = .
c
Повышение давления при непрямом ударе
t
Δpнепр = Δруд ,
(5.20)
T
где Т – время закрытия задвижки.
Самым простым методом, позволяющим избежать прямого гидравлического удара, является
медленное закрытие задвижки. Этому требованию вполне удовлетворяют вентили различных
конструкций и задвижки, менее всего – краны и клапаны. На насосных станциях, где имеется
опасность возникновения гидравлического удара при отключении насосного агрегата в связи с
аварией электросети, необходимы самостоятельные мероприятия по борьбе с гидравлическим
ударом.
В водопроводах внутри зданий, где длины участков невелики и фаза удара незначительна,
но есть быстродействующие запорные приспособления (краны), возможно образование непрямого гидравлического удара. Поскольку сила его прямо пропорциональна скорости течения до
удара, то скорость течения воды в сети ограничивают 2,5 м/с.
5.2.4. Гидравлический таран
Гидравлический таран (см. рис. 5.14) служит для подачи воды на небольшую высоту без
использования насосов.
Подача воды осуществляется только с использованием энергии гидравлического удара. Для
пуска тарана в действие открывают ударный клапан 3. Вода начинает поступать по трубопроводу 1 в рабочую камеру 2 и вытекает через клапан 3 наружу. С увеличением скорости движения вытекающей из клапана жидкости клапан поднимается и мгновенно закрывает отверстие в
рабочей камере. Происходит гидравлический удар. Давление в рабочей камере резко повышается, и под действием этого давления открывается клапан 4. Часть воды после открытия клапана 4
поступает в воздушный колпак 5, давление в нем повышается, и часть жидкости поступает в
трубопровод 6. Так как часть жидкости из рабочей камеры поступила в воздушный колпак, то
давление в рабочей камере понижается, в результате под действием силы тяжести открывается
клапан 3, и вода вновь по трубопроводу 1 начинает поступать в рабочую камеру и вытекать через клапан 3 наружу. Затем клапан 3 закрывается, происходит гидравлический удар, процесс
повторяется, и новая порция воды вновь поступает в воздушный колпак и в трубопровод 6.
96
5
H2
H1
6
1
4
3
2
Рис. 5.14. Гидравлический таран
Таким образом, гидравлический таран затрачивает часть расхода на подъем воды по нагнетательному трубопроводу на высоту Н 2 , которая может достигать 10 Н1 . Обычно теряется половина расхода.
Пример расчета короткого трубопровода
Вода перетекает из резервуара А в резервуар В по трубопроводу с диаметрами d1 = 100 мм и
d2 = 60 мм и длиной l1 = 15 м и l2 = 10 м (рис. 5.15). Необходимо определить расход воды при
разности уровней в бассейнах H = 300 см, и построить линии полного, пьезометрического и
геометрического напоров. Трубопровод стальной сварной, умеренно заржавевший.
H
2
2
hB
d2
hA
d1
1
1
0
0
A
B
Рис. 5.15. К примеру расчета короткого трубопровода
Ответ: 1. После уточнения расход в трубопроводе Q = 0,0056 м3/с$
2. Линии полного, пьезометрического и геометрического напоров (рис. 5.17).
Пример расчета магистрального трубопровода
Определить необходимую высоту водонапорной башни h в точке А, если dAB = 200 мм, dAC =
250 мм, dCD = 125 мм. На участке ВС непрерывный расход q = 0,05 л/с, на п.м., а в точках С и D
сосредоточены расходы QC = 10 л/c и QD = 12 л/c. Свободный напор НСB = 10 м, длины новых
трубопроводов: lАВ = 400 м, lBC = 300 м, lCD = 200 м.
Ответ: Высота водонапорной башни равняется 18,65 м.
h
q
A
B
QС
QD
С
D
Рис. 5.16. К расчету потерь напора в магистральном трубопроводе
97
Рис. 5.17. Построение линий полного и пьезометрического напоров
98
Рекомендуемая литература
1. Альтшуль А.Д. и др. Гидравлика и аэродинамика. – М.: Стройиздат, 1987, с. 15-19, 39-40,
99-106, 151-156, 158-175, 196-219,301-316,317-319.
2. Башта Т.М. и др. Гидравлика и гидравлические машины. – М.: Машиностроение, 1981.
3. Вакина В.В. Машиностроительная гидравлика. Примеры расчетов. – Киев: Вища школа,
1987.
4. Лойцянский Д.Г. Механика жидкости и газа. – М.: Наука, 1973, с. 12, 111-114, 132-136.
5. Астров Б.В. Сборник задач по гидравлике. – М.: УДН, 1986, с. 4-101.
6. Смыслов В.В. Гидравлика и аэродинамика. – Киев.: Вища школа, 1979, с. 11-13, 25-40,
95-96,107-124, 137-142,156-176,180-219, 229-233, 240-246, 303-306.
7. Угинчус А.А. Гидравлика и гидравлические машины., 1980.
8. Калинин А. В., Козлов Г. С. Методические указания для студентов очно-заочной формы
обучения по дисциплинам: «Гидравлика», «Гидравлика и гидравлические машины»,
ТГУ, 2003 г.,66 стр.
9. Калинин А. В. Лабораторный практикум по дисциплине «Гидравлика», ТГУ, 2005 г.,45
стр.
Дополнительная литература
10. M. Hanauer Mecanique des fluides. – VONTREUIL, Breal, 1991.
Материально техническое обеспечение дисциплины
На лекционных занятиях преподаватель использует аудиовизуальное оборудование. Необходимое количество экземпляров курса лекций в электронном виде находится в научной библиотеке университета.
Требования к уровню знаний студентов
После завершения курса студент должен знать теоретические основы кинематики и динамики течения жидкостей, законы сопротивления движения жидкостей и газов, основы моделирования гидромеханических явлений, расчётным путём определять характеристики реального потока жидкости и его физической модели.
99
Download