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

2982

.pdf
Скачиваний:
0
Добавлен:
15.11.2022
Размер:
2.73 Mб
Скачать

При

известных

значениях

 

радиуса

r0

и

скорости

невозмущенного потока V0

момент диполя находится из равенства

 

 

 

M

2 V r

2 .

 

 

 

 

 

 

 

0

0

 

 

 

 

Подставляя значение момента диполя в формулы для

потенциала

и функции

тока

и учитывая,

что

x

r cos

и

y r sin

, получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

r 2

 

 

 

V0

1

 

0

 

 

r cos

 

и

 

V0

1

 

0

 

 

r sin .

 

 

 

r2

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проекции составляющих скорости в полярных координатах

равны:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vr

 

 

d

и

VS

 

d

1 d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

dS

 

r d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

2

 

Vr

V0

1

0

 

cos

;

VS

 

 

V0

1

0

 

sin .

 

 

r

2

 

 

 

 

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На поверхности цилиндра, т.е. при значении r r0 , имеем:

 

 

 

Vr

 

 

0,

 

VS

2V0 sin .

 

 

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

Точки А и В на цилиндре, для которых и 0 , называются критическими, так как скорость в этих точках равна нулю.

Определим циркуляцию скорости вокруг цилиндра:

2

 

2

2

 

r0

Г VSdS VSrd

V0 r 1

sin d 0.

r 2

0

 

0

Таким образом, имеем бесциркуляционное обтекание кругового цилиндра.

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

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

91

поверхности цилиндра единичной длины на площадку dS , равна

PdS .

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

dX

P cos dS;

dY

P sin

dS.

Интегрируя dX и dY , находим, что

 

2

 

2

 

 

X

P cos r0d ;

Y

P sin

r0d .

0

 

0

 

 

Подставляя

значение

давления

согласно уравнению

Бернулли и значение скорости на поверхности цилиндра в предыдущие равенства, последние представим в виде:

 

 

2

 

 

 

V 2 sin 2

 

 

 

X

r

C

 

 

 

cos

d

0;

 

 

 

 

0

 

2

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2 sin 2

 

 

 

X

r

C

 

 

 

V

sin

d

0.

 

 

 

 

0

 

2

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

Противоречие между экспериментом и теоретическим выводом получило название парадокса Даламбера-Эйлера.

7.2. Циркуляционное обтекание кругового цилиндра

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

В этом случае потенциал скорости и функция тока равны:

 

 

 

r

2

 

 

Г

 

V0

1

 

0

 

r cos

 

 

 

;

 

r 2

2

 

 

 

 

 

 

 

 

 

r

2

 

 

Г

 

 

V0

1

0

 

 

r sin

 

 

ln r.

 

r 2

2

 

 

 

 

 

 

92

Проекции скоростей Vr и VS равны:

 

 

 

 

 

d

 

 

 

r

2

 

 

 

 

 

 

 

Vr

 

V0

1

0

cos

;

 

 

 

 

dr

 

r 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 d

 

 

 

r 2

 

 

 

Г 1

 

VS

 

 

 

 

 

V0 1

 

0

 

sin

 

 

 

 

.

 

2 dS

 

 

r 2

 

 

2 r

 

 

 

 

 

 

 

 

 

Убедимся, что в рассматриваемом случае циркуляция скорости жидкости вокруг цилиндра равна циркуляции вихря:

 

 

2

 

 

2

 

 

 

 

2

 

Г

 

 

Г

VSdS

V0

1

 

r0

 

sin

rd

 

 

 

rd

 

r2

0 2

 

 

 

0

 

 

 

 

 

 

r

 

 

 

 

r 2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V0r 1

 

0

cos

 

Г

Г.

 

 

 

 

 

 

r2

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На поверхности цилиндра, т.е. при r

r0 , имеем V

 

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

V

 

2V sin

 

Г

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

0

 

 

 

2 r0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким

образом,

наложение

циркуляционного потока

нарушает симметрию.

Картина имеет вид, изображенный на рис. 7.2.

Положение критических точек А и В найдем из условия

V 0 , откуда

 

 

 

 

 

S

 

 

 

 

 

sin

 

Г

1

.

kp

 

 

 

4

 

V0 r0

 

 

 

 

Очевидно, что при Г 4 V0r0 будут существовать два критических угла в 3 и 4 квадратах, как показано на рис. 7.2.

93

Рис. 7.2

 

При Г 4 V0 r имеем sin kp

1 и kp 2700 .

Это означает, что обе критические точки совмещены в одной. Если принять, что Г 4 V0 r0 , то получим, что sin kp 1,

что невозможно.

Это означает, что критическая точка сходит с цилиндра.

7.3. Теорема Н.Е. Жуковского о подъемной силе крыла бесконечного размаха

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

Этот теоретический вывод справедлив для тела произвольной

формы.

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

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

Н.Е.Жуковский, используя модель идеальной жидкости, предложил искать источник силового воздействия потока на тело в образовании циркуляции.

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

94

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

Жуковский заменил реальное крыло некоторым воображаемым жидким крылом, ограниченным замкнутой линией тока, предположив, что внутри этого жидкого крыла расположен вихрь (или система вихрей) с интенсивностью, равной сумме интенсивностей вихрей, образованных в тонком слое на поверхности тела при обтекании его реальной жидкостью. Такой вихрь Жуковский назвал присоединенным вихрем.

Докажем теорему Жуковского о подъемной силе для частного случая обтекания кругового цилиндра при наличии циркуляции вокруг него.

Определим силовое воздействие потока на круговой цилиндр с радиусом r при его циркуляционном обтекании идеальной жидкостью.

Поскольку жидкость идеальная, то на поверхности цилиндра будут действовать только силы давления.

Выделим на цилиндре участок высотой, равной единице.

В окрестности произвольной точки, находящейся на окружности цилиндра, элементарная сила давления равна pdS , где

dS - элементарная длина дуги. Ее проекции на оси координат равны:

dX

p cos

dS

dY

p sin

dS

Используя уравнение Бернулли, после преобразования и интегрирования, предыдущие равенства представим в виде

 

2

 

pV 2

d ;

X

r

C

 

 

cos

 

 

 

0

 

2

 

 

 

0

 

 

 

 

 

 

 

 

 

 

2

 

 

pV 2

d .

Y

r

C

 

 

sin

 

 

 

0

 

 

 

 

 

0

2

 

В результате имеем

X 0 .

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

После определения проекции силы на ось y получим, что

95

к лобовому

Y pV0 Г .

Эта формула называется формулой Жуковского и является основной формулой аэродинамики.

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

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

Вданном случае в точках, лежащих на одной вертикали, давление на верхней поверхности цилиндра меньше, чем на нижней,

аскорости течения, наоборот, на верхней поверхности больше, чем на нижней.

Врассмотренном выше циркуляционном обтекании цилиндра последний можно заменить одиночным вихрем той же интенсивности, что и цилиндр. Формально можно заменить любое обтекаемое тело соответствующей системой вихрей.

На возможности замены тела эквивалентной системой вихрей впервые указал Н.Е.Жуковский.

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

Метод конформного преобразования изложен в приведенной литературе [2]

7.4. Теория тонкого профиля

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

Отношение подъемной силы крыла Y сопротивлению Q называется качеством крыла K , т.е.

96

 

 

 

CyS

V2

 

 

 

 

Y

 

 

 

Cy

 

2

 

 

 

 

 

K

Q

 

CxS

V2

 

Cx

.

 

 

 

2

 

 

 

 

 

 

 

 

 

 

Как правило, крыло самолета имеет форму, симметричную относительно некоторой плоскости – плоскости симметрии.

Любое сечение крыла плоскостью, параллельной плоскости симметрии крыла, называется профилем крыла.

В разных сечениях профиль крыла может быть различным по форме, размерам и ориентации.

Один из возможных профилей изображен на рис. 7.3.

Точка А - передняя кромка профиля, В - задняя кромка, или точка схода профиля.

Рис. 7.3

Линия АВ, соединяющая две наиболее удаленные точки профиля, т.е. переднюю. и заднюю кромки профиля, называется хордой профиля.

Хорда делит профиль на две части – верхнюю и нижнюю. Угол между хордой профиля b и направлением

невозмущенного потока называется углом атаки , если вектор

скорости невозмущенного потока параллелен плоскости профиля. При изучении геометрических характеристик профиля

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

97

В этой системе координат осей уравнения верхнего и нижнего контуров профиля соответственно имеют вид:

Yв f1 (X) ,

Yн f2 (X) .

Толщина профиля в произвольной точке хорды выражается как разность ординат точек Yв и Yн .

Наибольшая разность YB YH называется максимальной толщиной, или просто толщиной, профиля и обозначается буквой C .

Отношение максимальной толщины профиля C к хорде b

носит название относительной толщины профиля C , т.е.

 

 

 

C

 

 

 

C

 

 

 

 

 

, или в процентах: C

100 0 0

 

 

C

 

 

 

 

b

 

 

 

 

 

 

 

 

b

 

 

Линия,

соединяющая середины

 

отрезков

YB YH ,

построенных в разных точках хорды, называется средней линией профиля.

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

Наибольшая ордината средней линии называется кривизной профиля f , а ее отношение к хорде b называется относительной кривизной, т.е.

 

 

f

 

 

 

f

 

 

 

, или в процентах, f %

100% .

f

 

 

b

 

 

 

 

b

 

Абсциссы наибольшей толщины профиля и наибольшей

кривизны соответственно обозначаются через Xc и Xf . Отношения этих величин к хорде b носят названия

относительных абсцисс соответственно толщины и кривизны, т.е.

 

 

 

 

 

Xc

 

 

 

 

Xf

 

 

 

Xc

Xf

 

 

b

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значения Xc для

дозвуковых

профилей колеблются в

пределах 25-30 % , у сверхзвуковых - 40-50 % .

Относительные толщины аэродинамических профилей крыльев и лопастей винтов обычно находятся в пределах от 3-4 % до

20-25 % .

Тонкие профили применяются на концах лопастей винтов и на крыльях сверхзвуковых самолетов.

98

Относительная кривизна современных профилей крыльев и лопастей винтов обычно не превышает 2 % .

Кривизна профилей самолетов 20-30-х годов доходила до 6-8 % .

Ограничимся рассмотрением простейшей задачи теории тонкого профиля, основанной на замене профиля непрерывно распределенными по его средней линии вихрями и впервые рассмотренной английским ученым Глауэртом.

Введем понятие о погонной циркуляции в данной точке

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

Выделим на этой кривой элемент S и обозначим через Г циркуляцию скорости от вихрей, расположенных на элементе S .

Погонной циркуляцией в точке называется

lim

Г

 

.

 

 

 

 

S 0

S

 

dS

 

Введя понятие погонной

циркуляции

x , обратимся к

рассмотрению теории тонкого профиля.

Допустим, что плоский установившийся поток идеальной несжимаемой жидкости, движущейся на бесконечности с заданной

скоростью V00 , обтекает тонкий профиль, длина которого ОА = b

(рис. 7.4)

Поместим начало координат в передней кромке профиля, ось

x направим вдоль хорды

к задней кромке

профиля,

ось y -

перпендикулярно ей вверх (рис. 7.4).

 

 

 

Введем следующие упрощения:

 

 

 

– будем считать

профиль

 

достаточно

тонким,

мало

 

отличающимся

от своей

средней

 

линии ОА;

 

 

 

будем

полагать

профиль

 

слабоизогнутым, так что его средняя

 

линия ОА незначительно отличается

Рис. 7.4

от прямой ОА, являющейся хордой

 

профиля;

 

 

99

y(x)

–будем полагать, что обтекание этого тонкого слабоизогнутого профиля происходит под малым углом атаки.

Заменим профиль непрерывно расположенными по его средней линии ОА вихрями, напряжение которых на элементе dx обозначим x dx .

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

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

Учитывая, что профиль слабо изогнут, будем считать, что вихри расположены вдоль хорды профиля, а не на его средней линии.

Тогда скорость, индуцируемая вихрями в точке x1 хорды, будет равна:

b

(x)dx

 

Vi

.

 

2 (x x1 )

 

0

 

Эту индуцированную скорость Vi x1 на хорде профиля

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

Наличие этой скорости изменит угол атаки на величину ,

так что

tg

Vi

V00

 

и местный угол атаки будет равен VVi .

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

Vi dy , V00 dx

где y – заданное уравнение средней линии профиля, т.е. dydx – известная функция.

Подставляя в полученное уравнение значение Vi , получаем

100

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]