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

2982

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

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

n

sin n

 

 

 

 

 

 

 

Vy

 

V00

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

sin

 

 

 

 

 

 

 

 

 

 

 

 

 

n

1

 

 

 

 

 

 

 

Подставляя

Vy

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

sin

An sin n

 

 

 

 

 

 

a sin

,

где

 

 

 

 

 

 

ab

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

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

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

111

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

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

 

 

 

V 2

l / 2

 

 

 

 

Y

 

 

Г z dz .

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

l / 2

 

 

 

 

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

, получаем

 

 

 

 

Y V 2l2

A

n

sin n

sin d

V

2l2A

 

.

 

00

 

 

 

00

1 2

 

0

 

 

 

 

 

 

 

 

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

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

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

 

Y

 

 

l2

Cy

 

 

 

 

A1.

V2

 

 

 

 

S

S

 

2

 

 

 

 

 

 

 

Обозначив

l 2

S

,

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

Cy A1

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

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

Cy a A1 ,

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

где

112

a dCy

d

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

C

y

a ,

 

 

 

 

 

 

C

y

a

00

,

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

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

 

 

 

 

Из этих равенств следует:

 

 

 

 

 

 

 

 

 

 

 

 

Cy

1

1

 

 

Cy

1

 

1

.

00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

a

 

 

 

 

 

a

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

0

 

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

 

1

 

 

1

 

 

1

 

 

, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

a0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

00 C y 1

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

,

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

на

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

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

 

 

Cy

1

1 ;

1

 

 

1

 

 

 

 

 

Cy

 

1

 

2

 

 

 

2

 

 

 

 

2

 

 

(считается, что Cy1

Cy 2

Cy ).

Рис. 8.6

113

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

0

2

Полученная зависимости Cy от

 

0

57,3

Cy

1

1

 

1

2

.

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

формула

позволяет

 

 

пересчитать кривую

для крыла с

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

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

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

 

 

l / 2

X

i

V Г z z dz .

 

00

 

 

l / 2

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

 

X

i

V

2 l2 00

nA

2

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

n

2

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

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

 

 

 

X i

 

 

 

 

 

 

00

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Cxi

 

 

 

 

 

 

 

 

 

nAn .

 

 

 

 

V 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

00

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

00

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

крыла величина

nA 2

положительная,

т.

 

е.

C

x i

, всегда есть

 

n 1

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сопротивление.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Представим значение Cx i , в следующем виде:

 

 

 

 

 

 

00

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 n

nA n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

Cxi

 

 

A1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая, что

 

A1

 

C y

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

nA n

 

 

 

 

 

C xi

 

C y

n 1

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

A1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

114

или

Cxi

C y

2

1

,

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

00

 

2

 

 

 

 

 

nA n

 

1

 

n 1

 

.

 

 

A12

 

 

 

 

 

 

Значения для различных крыльев приведены в рекомендованной литературе [2, 7].

Из значения Cx i следует, что

Cxi

Cxi

 

C y 2 1

1

1

2

.

 

2

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

Эта

формула

позволяет

пересчитать

индуктивное

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

крыла

с удлинением

1

на

индуктивное

 

 

 

 

 

 

 

 

 

 

 

сопротивление крыла с удлинением

2 .

 

 

 

 

 

8.6. Форма крыла в плане с наименьшим индуктивным сопротивлением

Наименьшее индуктивное сопротивление имеет крыло такой формы в плане, у которого 0 или

00

n A2

0

, т.е. когда A

A ... 0 .

n

 

2

3

nz

Вэтом случае циркуляция Г будет выражаться только через

первый член ряда, т.е.

 

 

Γ θ

2 l V00 A1 sin θ .

При θ

π

Г0

2 l V00 A1 ,

 

2

 

 

 

 

 

т.е.

 

A1

 

 

,

 

 

 

 

 

 

 

 

2 l V00

115

где Г0 – циркуляция в среднем сечении крыла.

Так как z

l

cosθ ,

2

 

 

Г 2lV00 A1 1

z2

 

.

l / 2 2

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

Г Г0 1

Z 2

 

.

l / 2

2

 

 

 

Как видно, циркуляция Г z в этом случае меняется по

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

Индуцированная в этом случае скорость равна:

Vy A1 V00 Г2l0 ,

т.е. постоянна по размаху крыла так же, как и угол скоса

Δα

Vy

A

Γ

 

Cy кр

.

 

 

 

 

 

V00

1

2l V00

 

π λ

 

 

 

 

 

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

 

C 2

 

Cxi

y кр

.

πλ

 

 

Циркуляция в любом сечении равна:

1

2 Cy bV00 ,

так что можно записать:

1

C

bV

2l V

A

1

 

 

z2

 

,

 

 

 

 

 

 

2

 

y 00

 

00

1

 

 

l / 2

2

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4l A

 

 

 

 

z2

 

 

 

 

 

 

 

 

b

1

1

 

 

 

 

.

 

 

 

 

 

Cy

 

 

l / 2 2

 

 

116

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

9.OCHOВЫ ГАЗОВОЙ ДИНАМИКИ

9.1.Основные уравнения газовой динамики

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

Считается, что f P,T , где Р – давление,

Т – температура.

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

dVx

Vx

Vx

Vy

Vx

Vz

Vx

1

 

P

;

dt

x

y

z

 

ρ

 

x

 

 

 

 

 

 

dVy

V

Vy

V

Vy

V

Vy

1

 

P

;

 

 

 

 

 

 

 

 

dt

x

x

y

y

z

z

 

ρ

 

y

 

 

 

 

 

 

 

 

dVz

 

Vx

Vz

 

Vy

 

Vz

 

Vz

 

Vz

1

 

P

,

 

dt

 

x

 

 

y

 

 

z

 

ρ

 

z

 

 

 

 

 

 

 

 

 

 

 

 

уравнение неразрывности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρVx

 

 

ρVz

 

 

ρVz

 

0.

 

 

 

 

 

 

 

 

x

 

 

y

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Неизвестные Vx ,Vy ,Vz ,

, P – функции координат и времени.

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

состояния газа:

 

 

P

g RT

 

 

 

ρ

 

 

Вэто уравнение вошла шестая неизвестная - температура Т.

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

117

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

Вэтом случае dQ = 0 – шестое уравнение, где Q – количество тепла в газе.

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

Для изоэнтропического процесса справедливо уравнение Пуансона:

Pconst ,

ρk

где k

Cρ

отношение теплоемкостей.

Cv

 

 

9.2. Уравнение Бернулли для адиабатических процессов

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

Это соотношение называется уравнением Бернулли и имеет

вид

V2

 

dP

C,

2

 

ρ

 

 

причем С есть функция линий тока.

Дифференцируя уравнение Пуассона P k C , получаем dP C K ρk 1.

После интегрирования приходим к соотношению

 

 

 

dP

C

 

k

ρk 1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из уравнения Пуассона

C

 

P

k

следует что

 

 

 

dP

 

 

P

 

 

k

 

ρk 1

 

k

 

 

P

.

 

 

 

 

ρk k

 

 

 

 

 

 

 

ρ

 

1

 

 

k

1 ρ

118

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

 

 

 

 

 

k

 

 

 

P

 

V

2

 

C .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

1 ρ

2

 

 

 

 

 

 

 

 

 

 

Используя уравнение состояния газа, приходим к следующей

форме уравнения Бернулли:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

g RT

V 2

C .

 

 

 

k

1

2

 

 

Как известно из физики, скорость распространения звука в

среде определяется формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

dP

 

или

a2

 

dP

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

представим в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

dP

 

C K

k 1

K

P

 

 

K g RT ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

g R T

 

a2 / K .

 

 

 

Подставляя значение

g RT в уравнение Бернулли, приведем

его к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a2

 

V 2

 

C .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из этого уравнения следует, что с ростом скорости движения газа скорость распространения звука в нем уменьшается, а с уменьшением скорости движения газа – увеличивается.

Используя понятие теплосодержания, или энтальпии i , где

 

i

 

C

T и

i

C

 

 

P

,

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

ρR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и учитывая, что gR

C

p

 

C

V

и

k

C

p

 

/ C

V

, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

Cp

 

 

P

 

 

k

 

P

 

 

 

k

 

g RT

Cp CV

 

 

k 1

 

 

 

 

k 1

Подставляя значение энтальпии в уравнение Бернулли, представим его в виде

119

 

V 2

i

 

 

C .

2

 

 

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

9.3. Основные соотношения для одномерных газовых потоков

Рассмотрим изоэнтропическое течение газа в трубе переменной плотности (рис. 9.1)

Рис. 9.1

В сечении A1 все параметры газа будем обозначать индексом 1, а эти же параметры в сечении А будем записывать без индекса.

Составим уравнение Бернулли для двух сечений A1 и А в

виде:

 

 

 

 

V 2

 

 

V 2

 

 

 

 

 

 

 

i

 

1

 

i

 

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если предположить,

что

V1

 

 

0 и

обозначить

параметры

покоящегося газа индексом "о", то получим

 

 

 

 

 

 

 

 

i0

i

V 2

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

g RT0

k

 

 

g RT

V 2

 

 

 

 

 

 

 

 

 

.

 

 

k

 

k

1

2

 

 

1

 

 

 

 

 

 

 

 

Из полученных

соотношений видно,

 

что при i0

const c

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

В несжимаемой жидкости температура меняется только при

120

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