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

2864

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

Поэтому вихрь выходит за концы крыла и располагается вдоль по потоку, образуя, так называемые вихревые усы (рис.

8.1).

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

свободными вихрями. Рис. 8.1 В рассмотренной простейшей вихревой схеме крыла

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

Следовательно, местная подъемная сила в каждом сечении крыла также должна быть постоянной. Это не подтверждается экспериментами.

Вихревую схему обтекания крыла, учитывающую переменность циркуляции вдоль его размаха, можно получить, если заменить крыло не одним П – образным вихрем, а системой П – образных вихрей (рис. 8.2.).

Рис. 8.2

101

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

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

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

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

8.2. Скос потока и индуктивное сопротивление крыла конечного размаха

Рассмотрим действие вихревой пелены, сбегающей с крыла, на произвольную точку А (рис. 8.3.). С каждого элемента крыла dz будет сбегать элементарный вихрь с циркуляцией .

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

= dz, – координатой z1 , то Рис. 8.3

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

dVy 4 (z1 z)

102

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

1

Vy 4

 

 

z

 

 

 

 

 

l / z

 

dz

1

 

 

 

,

 

1

 

l / z

z1 z

где Г Г (z1 ) , а координата z для данной точки А есть величина, постоянная.

В результате скорость Vy также будет меняться по

размаху крыла, т.е. Vy Vy (z) .

Полученный интеграл несобственный. При вычислении он разбивается на два, и точка z1 z исключается.

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

Рассмотрим сечение крыла в точке A z

(рис. 6.4.).

 

В

произвольной

точке

 

хорды

скорость

частиц

 

жидкости

будет

состоять

 

из

двух

составляющих:

 

скорости

 

набегающего

 

потока

V0

и

скорости Vy

Рис. 8.4

индуцированной вихрями.

В результате в рассматриваемой точке хорды крыло

будет обтекаться некоторой скоростью V, направление которой

составит с хордой крыла

угол

и

 

,

называемый

 

 

 

 

 

 

истинным углом атаки.

Таким образом, вблизи крыла происходит изменение

угла атаки, т.е. происходит скос потока.

 

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

,

называемым углом скоса потока.

 

Значение этого угла невелико и составляет несколько градусов.

103

Как видно из рис. 8.4. tg

Vy V . Знак минус

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

Учитывая малость угла , можно записать:

tg

Vy

V0

 

Угол скоса потока так же, как и Vy , меняется вдоль по размаху

крыла.

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

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

скорости набегающего потока, на угол

.

По этой причине у результирующей силы появится

составляющая X i по направлению невозмущенного потока.

Эта дополнительная сила X i , появившаяся в результате

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

Найдем величину силы индуктивного сопротивления. Для этого выделим на крыле элементарный участок шириной dZ, на который действует подъемная сила dY. Очевидно, что индуктивное сопротивление этого участка равно:

dXi

dY

.

 

 

Согласно теореме

Жуковского

dY

V0 Г z dz .

Используя значения dY и

,

формулу

для

dX i можно

представить в виде:

 

 

 

 

dXi

Г (z)Vy dz .

 

 

104

Полная сила индуктивного сопротивления равна:

l / 2

X i Г (z) Vy dz .

l / 2

Учитывая

значение

Vy ,

силу

 

индуктивного

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

 

 

 

 

 

 

 

X i

 

 

l / 2

l / 2 (z)

 

dz

 

 

 

 

(z) dz

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

4 l / 2

l / 2

dz1

 

z1

z

Зная характеристики крыла и

Г z1

, можно

вычислить X i .

 

 

 

 

V

2

 

 

 

 

 

Разделив

X i

на S

0

, получим коэффициент индуктивного

2

 

 

 

 

 

 

 

 

 

 

 

сопротивления C

 

X

 

/ S

 

V2

.

xi

i

 

 

 

 

 

 

 

 

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

.

 

 

2

y

00

 

 

 

 

 

Допустим, что угол атаки в данном сечении отсчитывается от так называемой аэродинамической хорды, т.е. от направления, для которого местный Сy 0 (рис. 8.5.).

105

Обозначив этот угол через

 

a , а угол атаки, при котором

Сy 0 , через 0 , можно записать,

 

a

0 .

 

 

 

Для

 

малых

 

углов

 

атаки

 

dCy

const a ,

так что Cy

a a ,

 

 

 

 

d

 

 

 

 

 

 

 

или Cy

a(a a0 ) .

 

 

 

 

 

 

 

Рис. 8.5

 

В

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

размаха

ист

a

,

где

 

 

 

 

 

 

Vy /V .

Рис 8.5

Таким образом, для крыла конечного размаха

 

 

 

 

 

 

 

 

 

 

 

 

 

Vy

 

Cy

 

a u

a( 0

 

) a( a

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V0

 

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г (z)

1

ab(z)V [

 

 

(z)

Vy

 

].

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

2

0

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя

 

сюда

 

значение

 

Vy

и

 

опуская индекс а,

получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

dZ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

l / 2

dZ

1

 

Г (z)

 

 

 

 

 

z

]

 

2 ab(z)V[ (z)

4 V -l / 2

 

 

z1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

Это уравнение носит название основного интегро-

дифференциального

уравнения

крыла

 

(интегро-

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

Г z

входит

в него и под знаком производной и под знаком интеграла).

Уравнение позволяет определить

циркуляцию

крыла

Г z по известным

функциям b(z)

и

z

, заданным

конструкцией крыла.

 

 

 

 

 

106

8.4. Приближенный метод расчета распределения

циркуляции по размаху крыла

Имеется ряд приближенных методов расчета распределения циркуляции.

Рассмотрим метод Глауэрта-Трефца.

В основном уравнении сделаем замену переменной, положив:

z

l

cos

(0

;

 

l

 

z

l

)

2

2

 

2

 

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

l

sin d

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате замены переменной Г z

перейдет в Г .

Представим неизвестную функцию Г

 

в виде ряда Фурье:

 

 

 

00

 

 

 

 

 

 

 

 

 

( )

2lV00

sin

 

n

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

Представим через новую переменную Vy . После преобразования получим:

 

 

00 nA sin n

Vy

V00

 

n

 

 

sin

 

 

n 1

Подставляя Vy в интегро-дифференциальное уравнение,

после преобразования получаем:

 

 

n sin

An sin n

 

a sin

где

 

ab

 

 

 

 

 

 

 

 

4l

 

 

 

 

 

 

Очевидно, что полученное равенство должно выполнятся в любом сечении крыла, т.е. при любом .

107

Это условие дает возможность определить коэффициенты An .

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

В результате получим конечную систему алгебраических уравнений для определения коэффициентов An . Необходимо

отметить, что этот метод эффективен только для крыльев простейших форм.

8.5. Определение аэродинамических сил, действующих на

 

крыло конечного размаха

 

 

 

Согласно теореме Жуковского

 

 

 

 

 

 

 

 

l / 2

 

 

 

 

 

 

Y

V 2

Г z dz

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

l / 2

 

 

 

 

Переходя к переменной

, получаем:

 

 

 

 

Y

V 2l 2

 

A sin n

sin d

V

2l 2 A

 

 

 

 

 

00

0

n

 

00

1 2

 

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

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

Коэффициент подъемной силы крыла равен:

C

 

 

Y

 

 

l

2

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

V 2

S

 

S

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Обозначив

l 2

S

где – удлинение крыла, получим:

108

Cy A1

Найдем формулу для пересчета плоского крыла с одного удлинения на другое.

Для крыла в целом имеем:

Cy a

A1 ,

отсюда для крыла определенной формы в плане a A1 ,

где a

dCy

 

d

 

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

Cy a ,

Cy a0 00 ,

где 00 – угол атаки крыла бесконечного размаха. Из этих равенств следует:

00 Cy

1

1

 

 

Cy

1 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

0

 

Введя обозначение

1

 

1

 

1

, получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

00 Cy 1

Эта формула позволяет пересчитать зависимость Cy по

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

Если имеются два крыла, составленные из одних и тех же профилей, но с различными удлинениями 1 и 2 , то для них можно написать равенства:

109

 

 

Cy

1

 

1

 

 

1

 

 

 

 

1

 

 

 

 

 

Cy

 

1

 

2

 

 

 

2

 

 

 

 

2

(считается, что Cy1 Cy 2 Cy ).

Вычитая из первого уравнения второе, получаем;

 

 

 

 

 

 

 

 

 

 

 

Рис. 8.6

0

0

57,3

Cy

1

1

1

2

 

,

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

57,3

Cy

1

1

1

2

.

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

Полученная формула позволяет пересчитать кривую

зависимости Cy от

для крыла с

 

 

1 , на крыло с 2 (рис. 8.6.)

Значения для различных крыльев приводятся в рекомендованной литературе.

Сила индуктивного сопротивления крыла равна:

 

l / 2

 

X i

V00 Г z

z dz

 

l / 2

 

После преобразования получим:

 

 

 

2

 

00

 

2

 

 

X

 

V

l

2

nA

 

.

 

 

 

 

 

 

i

00

 

 

n 1

n

 

2

 

Коэффициент индуктивного сопротивления равен:

110

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