2982
.pdf
|
1 |
b |
(x)dx |
|
dy . |
|
|
2 V00 |
0 |
x x1 |
|
dx |
|
Полученное равенство |
представляет собой интегральное |
уравнение относительно неизвестной функции x .
Решение основано на введении новой независимой
переменной , определяемой равенством |
|
|
||||||||
|
x |
|
b |
(1 - cos |
) . |
|
|
|||
|
2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|||
Очевидно, что при изменении х в пределах 0 < х< b |
||||||||||
меняется от 0 до . |
|
|
|
|
|
|
|
|
|
|
Решение |
ищется в виде ряда Фурье |
|||||||||
|
|
|
|
|
00 |
|
sin n ) . |
|||
( ) |
2V (A |
0 |
ctg |
|
|
A |
n |
|||
|
||||||||||
|
00 |
|
|
2 n |
|
|
||||
|
|
|
|
|
|
1 |
|
|
После преобразований уравнение, соответствующее условию безотрывного обтекания профиля, принимает вид
|
|
|
|
|
|
00 |
|
|
|
|
|
|
|
|
dy |
|
|
|
||
|
|
|
A0 |
|
|
|
|
An cos n |
|
. |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
n |
1 |
|
|
|
|
|
|
|
dx |
|
|
||
|
Коэффициенты A n |
этого уравнения определяются методом, |
||||||||||||||||||
обычным для рядов Фурье. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Умножая |
обе |
|
|
части |
последнего |
уравнения |
на |
||||||||||||
1, cos , cos 2 , cos 3 |
,... и |
|
интегрируя в |
пределах от 0 до |
, |
|||||||||||||||
получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
A |
|
|
|
1 |
|
dy |
d ; |
|
|
|||||||
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
0 dx |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
A |
|
|
2 |
|
|
dy |
cos n |
d , |
|
|
|||||||
|
|
|
n |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
0 dx |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
где |
dy(x) |
выражено через |
|
с помощью подстановки |
|
|||||||||||||||
|
|
|
||||||||||||||||||
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
b |
(1 |
cos ) . |
|
|
|
|
|
||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
101
Найдя значения коэффициентов, определим погонную циркуляцию X и индуцированные скорости, т.е. задача об
обтекании тонкого профиля решена.
Теперь можно найти подъемную силу, действующую на
профиль: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dY |
|
p (x)dxV00 . |
|
|
|
|
|
|
|
|
||||||||||||||||
Переходя к переменной |
и учитывая, |
что |
dx |
b |
sin d , |
||||||||||||||||||||||||||
2 |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dY |
|
|
b |
|
|
|
|
sin |
d |
|
V00 |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
|
00 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Y pV 2 b [A |
0 |
|
(1 cos ) |
|
|
|
A |
n |
sinn |
|
sin |
] d |
|
|
|||||||||||||||||
00 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
pV2 |
|
b |
(A |
|
|
A1 |
), |
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
00 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
Y |
|
|
|
|
|
2 (A |
|
|
|
A1 |
) . |
|
|
|
||||||||||
y |
b1 pV 2 / 2 |
0 |
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Аналогично вычисляется момент относительно передней |
|||||||||||||||||||||||||||||||
кромки профиля: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
M |
pV |
(x) x dx . |
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Переходя к переменной , после вычислений получаем |
|||||||||||||||||||||||||||||||
|
M |
|
|
|
|
|
pV 2 b2 |
( A A1 |
|
|
A2 |
) , |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
4 |
|
|
|
|
00 |
|
|
0 |
|
|
|
2 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
|
|
|
M |
|
|
|
|
|
|
|
( A |
|
A |
|
|
A2 |
) . |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
m b1b pV 2 / 2 |
2 |
0 |
1 |
|
|
2 |
|
|
|
|
Таким образом, значения подъемной силы и момента определяются лишь первыми тремя коэффициентами разложения погонной циркуляции в ряд.
102
Остальные коэффициенты ряда влияют только на распределение местных скоростей (а стало быть, и давлений) по
профилю. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя значения |
A0 , A1 , A2 |
|
в формулы для С |
y |
и C |
m |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Cy |
2 |
( |
|
|
0 ) |
|
|
|
|
|
|
||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
|
2( |
|
|
|
|
) |
|
1 |
C |
|
, |
|
|
|
m |
0 |
4 |
0 |
4 |
y |
|
|
|
|||||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
где
|
|
|
1 |
|
|
dy |
|
( |
) (1 |
cos |
) d , |
|||
|
|
0 |
|
|
0 dx |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
dy |
( |
) (1 |
cos 2 |
|
) d . |
||
|
|
0 |
|
|
|
0 dx |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||
Угол |
0 |
0 |
называется углом атаки нулевой подъемной |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
силы, т.е. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C |
y |
|
2 ( |
|
0 |
) , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
откуда |
|
|
|
dC y |
2 . |
|
d |
|
|
|
Следовательно, для тонких профилей С y линейно зависит от
, а производная не зависит от формы профиля.
Выражение для коэффициента момента можно записать в
виде
|
|
C |
|
C |
|
1 |
C |
|
, |
|
|
|
m |
m0 |
4 |
y |
|||||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||
где Cm0 2( 0 |
|
0 ) – |
коэффициент момента при нулевой |
|||||||
4 |
||||||||||
|
|
|
|
|
|
|
|
|
подъемной силе.
Таким образом, определение аэродинамических коэффициентов для тонкого профиля сведено к вычислению двух
103
простых интегралов для 0 и 0 , что не представляет труда, так как
форма профиля известна.
Формулы для С y и Cm хорошо согласуются с опытными
данными для профилей небольшой относительной кривизны и толщины.
Теория тонкого профиля находит наибольшее применение лишь в задачах, требующих определения суммарных аэродинамических характеристик профиля, так как получаемое с ее помощью распределение местных скоростей (давлений) по профилю значительно отличается от реального распределения.
8. ТЕОРИЯ КРЫЛА КОНЕЧНОГО РАЗМАХА
8.1. Модели крыла конечного размаха
Ранее было показано, что действие вихря с циркуляцией Г, находящегося в потоке, создает такую же подъемную силу, как и крыло, имеющее ту же циркуляцию.
Таким образом, воздействие крыла на жидкость можно заменить действием вихря бесконечной длины.
Если крыло имеет конечную длину, то подъемная сила крыла будет равна подъемной силе, создаваемой отрезком вихря длиной, равной длине крыла.
По теореме Гельмгольца вихрь не может обрываться в жидкости. Поэтому вихрь выходит за концы крыла и располагается вдоль по потоку, образуя, так называемые вихревые усы (рис. 8.1).
Таким образом, в первом приближении крыло конечного размаха можно заменить одним вихрем П – образной формы в плане.
Вихрь, расположенный по размаху крыла, называется присоединенным, а сбегающие с торцов крыла вихревые усы – свободными вихрями.
104
В рассмотренной простейшей вихревой схеме крыла конечного размаха, согласно теореме Гельмгольца и Стокса циркуляция вдоль размаха крыла получается постоянной.
Следовательно, местная подъемная сила в каждом сечении крыла также должна быть постоянной. Это не подтверждается экспериментами.
Рис. 8.1
Вихревую схему обтекания крыла, учитывающую переменность циркуляции вдоль его размаха,
можно получить, если заменить крыло не одним П–образным вихрем, а системой П–образных вихрей (рис. 8.2.).
Рис. 8.2
Вдоль каждого такого вихря циркуляция будет постоянной, но при переходе от одного вихря к другому будет меняться.
В этой схеме можно положить число вихрей равным бесконечности, и тогда циркуляция вдоль крыла будет плавно изменяться от центра к торцам. При этом за крылом свободные вихри сольются и образуют так называемую вихревую пелену, состоящую из бесконечно большого числа сходящих с крыла вихревых шнуров.
105
Исследования вихревой пелены за крылом показали, что она неустойчива и вскоре после сбегания с крыла свертывается в два вихревых шнура.
Однако рассматривать такую модель крыла сложно с математической точки зрения, поэтому ее заменяют вихревой пеленой.
8.2. Скос потока и индуктивное сопротивление крыла конечного размаха
Рассмотрим действие вихревой пелены, сбегающей с крыла, на произвольную точку А (рис. 8.3.). С каждого элемента крыла dz будет сбегать элементарный вихрь с циркуляцией dГ.
Если положение точки А характеризуется координатой z, а положение элементарного вихря, сбегающего с участка крыла
MN = dz, – координатой z1 , то
Рис. 8.3
скорость, индуцированная этим полубесконечным вихрем в точке А, очевидно, равна:
dГ dVy 4 (z1 z) .
Скорость в точке А, индуцированная всей вихревой пеленой, будет равна:
106
|
|
|
|
|
|
|
dГ |
z1 |
|
|
|
|
|
|
1 |
|
l / z |
|
dz1 |
|
|
||
|
Vy |
|
|
|
, |
|
|
||||
|
4 |
|
l / z |
z1 |
z |
|
|
||||
|
|
|
|
|
|||||||
где Г Г (z1 ) , |
а координата z |
|
для данной точки А есть величина |
||||||||
постоянная. |
|
|
|
|
|
|
|
|
|
|
|
В результате скорость Vy также будет меняться по размаху |
|||||||||||
крыла, т.е. V |
V (z) . |
|
|
|
|
|
|
|
|
|
|
y |
y |
|
|
|
|
|
|
|
|
|
|
Полученный интеграл несобственный. При вычислении он |
|||||||||||
разбивается на два, и точка z1 |
|
|
z исключается. |
|
|
||||||
Фактически это означает, что исключен из рассмотрения |
|||||||||||
элементарный вихрь, проходящий через точку А, |
скорость Vy в |
||||||||||
которой требуется вычислить. |
|
|
|
|
|
|
|
|
|||
Рассмотрим сечение крыла в точке A z (рис. 8.4.). |
|||||||||||
|
|
|
|
|
|
|
В произвольной точке хорды |
||||
|
|
|
|
|
|
|
|
скорость частиц жидкости |
|||
|
|
|
|
|
|
|
|
|
будет состоять из двух |
||
|
|
|
|
|
|
|
|
|
составляющих: скорости |
||
|
|
|
|
|
|
|
|
набегающего потока V0 и |
|||
|
|
|
|
|
|
|
скорости V |
y |
,индуцированной |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
вихрями. |
Рис.8.4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
В результате в рассматриваемой точке |
|||||||||
хорды крыло будет обтекаться некоторой скоростью V, направление |
|||||||||||
которой составит с хордой крыла угол |
и |
|
, называемый |
||||||||
|
|
|
|
|
|
|
|
|
|
||
истинным углом атаки. |
|
|
|
|
|
|
|
|
|
|
|
Таким образом, вблизи крыла происходит изменение угла |
|||||||||||
атаки, т.е. происходит скос потока. |
|
|
|
|
|
|
|||||
Величина скоса потока определяется углом |
|
, называемым |
углом скоса потока.
Значение этого угла невелико и составляет несколько градусов.
Как видно из рис. 8.4,
107
tg |
V V. |
|
y |
Знак минус берется потому, что принято считать угол скоса положительным, когда индуцированная скорость направлена вниз,
т.е. когда Vy отрицательно. |
|
|
|
Учитывая малость угла |
, можно записать: |
||
tg |
Vy |
. |
|
|
|
||
|
V0 |
||
|
|
|
Угол скоса потока так же, как и Vy , меняется вдоль по
размаху крыла.
При отсутствии скоса потока вектор результирующей силы воздействия потока на крыло согласно теореме Жуковского нормален
к направлению невозмущенного потока, т.е. к вектору V0 .
При наличии скоса потока вектор результирующей силы должен быть нормальным к направлению истинной скорости V. Следовательно, скос потока приводит к отклонению результирующей силы от направления, перпендикулярного скорости набегающего
потока, на угол |
. |
По этой причине у результирующей силы появится |
|
составляющая Xi |
по направлению невозмущенного потока. |
Эта дополнительная сила X i , появившаяся в результате скоса
потока, носит, название силы индуктивного сопротивления.
Найдем величину силы индуктивного сопротивления. Для этого выделим на крыле элементарный участок шириной dZ, на который действует подъемная сила dY. Очевидно, что индуктивное
сопротивление этого участка равно: |
|
|
||
|
dXi |
dY |
. |
|
Согласно теореме Жуковского |
dY |
V0 Г z dz . Используя |
||
значения dY и |
, формулу для dX i можно представить в виде |
|||
|
dXi |
Г (z)Vy dz . |
|
Полная сила индуктивного сопротивления равна:
l / 2
X i Г (z) Vy dz .
l / 2
108
Учитывая значение Vy , силу индуктивного сопротивления запишем в виде
|
|
|
|
|
l / 2 |
|
l / 2 |
|
dГ (z) |
|
dz |
|||
Xi |
|
|
|
|
|
|
|
(z) dz |
|
|
1 . |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
4 |
|
l / 2 |
|
l / 2 dz1 |
|
z1 z |
||||||
Зная характеристики крыла и |
Г z1 , |
можно вычислить X i . |
||||||||||||
|
|
|
V2 |
|
|
|
|
|
|
|
||||
Разделив Xi на |
S |
|
|
0 |
, |
|
получим |
|
коэффициент индуктивного |
|||||
|
|
|
|
|||||||||||
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
сопротивления Cxi |
|
Xi / S |
|
V2 |
. |
|
|
|
|
|
||||
|
|
2 |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
И коэффициент сопротивления крыла равен |
||||||||||||||
|
|
|
|
|
Cx |
Cxp |
Cxi , |
|
|
|
где Cxp – коэффициент профильного сопротивления, определяемый
сопротивлением трения и разностью давлений в передней и задней частях крыла.
8.3. Основное уравнение крыла конечного размаха
Выведем уравнение, позволяющее определить закон изменения циркуляции по размаху крыла, т.е. функцию Г z .
При наличии вихревой пелены для элемента крыла dz можно написать два следующих равенства:
dY |
Г (z) V dz ; |
dY |
C |
y |
(z) |
|
|
b (z) dz V |
2 . |
|||||
|
|
|||||||||||||
|
00 |
|
|
|
|
|
2 |
|
|
00 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Приравнивая их, получаем уравнение |
||||||||||||
|
|
связи: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Г (z) |
|
|
1 |
C |
|
(z) b (z) V |
. |
|||||
|
|
|
|
|
y |
|||||||||
|
|
|
|
|
2 |
|
|
|
|
00 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Допустим, что угол атаки в данном |
||||||||||||
|
|
сечении |
отсчитывается |
от так |
||||||||||
|
|
называемой |
|
|
|
|
|
|
аэродинамической |
|||||
|
|
хорды, т.е. от направления, для |
||||||||||||
|
|
которого местный С |
y |
0 (рис. 8.5). |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 8.5
109
С y
или
Обозначив этот угол через |
a , а угол атаки, при котором |
||||
0 , через |
0 , можно записать a |
0 . |
|||
Для малых углов атаки |
dCy |
const a , так что Cy a a |
|||
|
|
||||
d |
|||||
|
|
|
|||
Cy a(a |
a0 ) . |
|
В случае обтекания крыла конечного размаха |
ист |
a |
, |
|
|
где Vy /V .
Таким образом, для крыла конечного размаха
Cy a u a( 0 |
) a( a |
Vy |
) . |
|
V0 |
||||
|
|
|
||
Следовательно, |
|
|
|
Г(z) |
1 |
ab(z)V [ |
|
(z) |
Vy |
] . |
|
a |
|
||||
|
2 |
0 |
|
V |
|
|
|
|
|
|
|
Подставляя сюда значение Vy и опуская индекс а, получаем
|
|
|
|
|
|
|
dr |
dZ1 |
|
|
1 |
|
1 |
|
l / 2 |
|
dZ1 |
||
Г(z) |
ab(z)V[ (z) |
|
|
] |
|||||
2 |
4 |
V -l / 2 |
z1 |
z |
Это уравнение носит название основного интегродифференциального уравнения крыла (интегро-дифференциальное
уравнение потому, что искомая функция Г z входит в него и под знаком производной и под знаком интеграла).
Уравнение позволяет определить циркуляцию крыла Г z по известным функциям b(z) и z , заданным конструкцией крыла.
8.4. Приближенный метод расчета распределения циркуляции по размаху крыла
Имеется ряд приближенных методов расчета распределения циркуляции.
Рассмотрим метод Глауэрта-Трефца.
В основном уравнении сделаем замену переменной, положив
110