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

2982

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

 

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 будет сбегать элементарный вихрь с циркуляцией .

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

MN = dz, – координатой z1 , то

Рис. 8.3

скорость, индуцированная этим полубесконечным вихрем в точке А, очевидно, равна:

dVy 4 (z1 z) .

Скорость в точке А, индуцированная всей вихревой пеленой, будет равна:

106

 

 

 

 

 

 

 

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

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