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

книги из ГПНТБ / Поляхов Н.Н. Теория нестационарных движений несущей поверхности

.pdf
Скачиваний:
6
Добавлен:
29.10.2023
Размер:
3.78 Mб
Скачать

Скорость Vi вызывается в фиксированной точке следа вих­ рями этого следа и, кроме того, системой присоединенных вихрей профиля. Будем считать, что вихри эти распределены по контуру профиля с линейной плотностью ~(в. Обозначая ско­

рость, вызванную вихрями следа, через Vs, а вихрями профи­ л я — через Vn, получим

V*i =V*s + V * a.

У

Рис. 9.

Как видно из рис. 9, скорость, вызываемая в точке N следа, будет выражаться формулой

Vs f Г Х г ds',

2тс

где f' есть плотность вихрей следа в точке N', а г есть рас­ стояние N N . Скорость же V* представится в виде

1 I ТпХп dsx.

2тс

ЛВС

где гх есть расстояние от точки профиля, в которой поме­ щается вихрь "(п, до точки N, а интегрирование совершается по контуру профиля. Сила, действующая на профиль, на осно­ вании формулы (16) будет равна

й = £ f

+ £ j J Щ ^ -d sd s,

CD CD

CD Л ВС

1

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

Таким образом, сила будет выражаться формулой

7ХУ ~" dsdSv

CD ABC 1

29

Так

как

 

 

 

1 = l k , •(„ = i„ k ,

 

где k

есть орт оси z, то можно

написать,

что

 

й = р Я й

X Тп) dsu

(17)

 

L

 

 

где

 

 

 

 

j

ds

(17а)

 

 

r i

 

 

CD

' 1

 

 

 

 

есть индуктивная скорость, вызываемая вихрями следа на кон­ туре профиля. Формула (17) совпадает с обычной формулой, даваемой теорией крыла конечного размаха для сил индуктив­

ного сопротивления. Для вычисления силы R/ необходимо

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

Заметим, что формулу (17) можно еще установить и дру­ гим путем. Начнем, для простоты, со случая бесконечно тонкой дуги. На основании формул (14) и (14, а) будем иметь:

R* = f I" (и — i v f dz =

у I* (vt ivn)2

dz,

L

L

 

причем контур интегрирования

будет составляться из верхней

и нижней сторон обтекаемой дуги, которая

является поверх­

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

R* =

\

J [(Д2В-

v)H) - 2ivn (vtB -

vm)\ e

dz.

 

LB

 

 

 

 

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

 

 

 

 

Чв ~

v tH= Tn;

v* = J- (v tB+ v tH);

v*n = v nB= v nH,

получим

 

 

 

 

 

 

R* =

ip j

In «

-

*<) e-^dSi = j Tn (u* -

iv*)dsr

 

LB

 

L B

 

 

Скорости и* и v* представляют собой скорости возмуще­ ния, взятые на обтекаемой дуге. Это будут скорости, индуци­ руемые в точках этой дуги свободными вихрями следа и при­ соединенными вихрями fn профиля. Из полученной формулы вытекает, что

я,- = р J ( и Г х т п ) ^ ,

LB

30

причем

—>0. —> —> *

v . = V i + V П

есть индуктивная скорость, вызываемая на профиле свобод­ ными и присоединенными вихрями.

Легко убедиться, что интеграл

К ^ п Х Т п )^ ! = О,

 

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

 

K i— p J ( Vi X Тп) dsiy

(17в)

LB

 

где V/ выражается формулой (17а).

на случай замк­

Для распространения последней формулы

нутого контура следует рассматривать этот контур, как пре­ дельный случай дуги, начало и конец которой стремится к совпадению. Такого рода рассмотрение приведет полученную формулу (17в) к виду (17).

Переходим

теперь к рассмотрению величин,

зависящих от

 

ГФ4 +

дФ*

 

 

dt

 

Из формул (7)

и (6, а) имеем:

 

 

 

X s -\-iYs— /p j^ ro ^ -)— jjpj dz.

(18)

L

Так как циркуляционный поток и поток, создаваемый вихре­ вым следом, обтекают профиль, то функции тока Wr и Ws на контуре профиля имеют постоянное значение. Отсюда следует, что формула (18) может быть записана в виде

X + i Ys = -

ip j" (Г/ч +

dz,

(19)

 

L

 

 

где F4 и /^ — комплексные

потенциалы.

 

 

Отобразим конформно внешность профиля на внешность

круга и предположим, что вихревой след

профиля

таков, что

при отображении он переходит в действительную

полуось—•?

(рис. 10). Общее представление о форме такого следа профиля можно себе составить на основе того, что нулевая линия тока, которая получается при стационарном обтекании профиля поступательным потоком при угле атаки, дающем нулевую циркуляцию, переходит при отображении в контур круга и действительную полуось. Ясно, что след указанного вида вряд ли может существовать в действительности. Однако при дви­ жениях, мало уклоняющихся от прямолинейных и происходя-

31

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

тока.

 

 

с циркуляцией

Из гидродинамики известно, что два вихря

ЗГЯ,

вращающиеся в противоположные стороны,

имеют линии

тока

в виде

окружностей радиуса R \ причем (рис. 11)

 

 

 

\ V \ - \ t i \ ~ R * ,

 

где

| Ч I и | Ci | суть расстояния этих вихрей от центра окруж­

ности.

Комплексный потенциал указанной системы вихрей будет

 

 

 

С— С'

 

 

 

 

bFs = п s. ш ------ г

 

 

 

 

2w С—Сц

 

где С =

R ^

есть комплексная ко­

 

ордината точки, лежащей на окруж­

 

ности.

Помещение в центр окружности

Рис. Ю.

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

tr s СК -С ')

bFs

2** с -с ; *

Так как при безотрывном обтекании твердой окружности она является линией тока, то можно утверждать, что комплек­ сный потенциал §FS является потенциалом бесциркуляционного обтекания окружности потоком, создаваемым изолированным вихрем. Если вихри расположены непрерывно вдоль кривой £' так, что

8Г5 = i 0do',

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

(20)

32

Конформное соответствие C= C(z) предполагается извест­ ным, поэтому, зная Fs (С), мы всегда можем найти функцию

Fs {z) = F U z ) ) .

Заметим, кроме того, что

где т и ds' — вихревая плотность и элемент длины вихревого следа в плоскости профиля, есть не что иное, как циркуляция всех вихрей следа, которая в случае, если профиль начал двигаться из состояния покоя, (Г0 = 0) должна быть равна с обратным знаком циркуляции вокруг профиля, т. е.

| ids' — — Г. s

На основании сделанных замечаний формулу (19) можно записать в виде

x + i K = - ^ 4 - J J ’( - l" c + ln-ir^ f-)T * ■ ! « » ' . <21>

к S

где k указывает, что интегрирование в плоскости С совершает­ ся по контуру круга. Но

ибо написанная подстановка дает — 2шг0, где z0 относится

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

 

z(C) =

С+

f -

+ 2

С„

 

^ + Cl

то, замечая, что

 

 

 

л -1

 

 

 

 

 

 

С idt

п

Г

1

rfC

n _ . 1

J с-с -°, J чп С-С' _ 2 с"

кк

(23)

f

- = 2ы1л ;

С- Ci

с -Cl

 

3 Н. Н. Поляхов

33

получим

 

 

 

 

 

h ----- 2тсг ^ Сх

+

p b j

+

2ш (с0 — z0).

 

Интегрирование по С выполнено, поэтому вместо

теперь

можно всюду писать С.

 

 

 

 

Полагая

 

 

 

 

 

с0 =

Z0; С=

ге10; С, =

^ -

е№,

 

получим

 

 

 

 

 

+ г г;

где через f ( s ) обозначено выражение, стоящее в предшест­ вующем интеграле в круглых скобках.

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

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

! / [ » « ) ]

 

 

а, следовательно, в случае

бесконечно длинного следа,

когда

Sj = со, получим

 

 

X's + iY's =

— ip j* yVc - t ds.

(24)

 

S q

 

Особенно просто вести

вычисления для случая пластинки,

с которой сходит вихревой

след, являющийся ее продолже­

нием. В этом случае

все сп

равны нулю, cos6 = — 1 и

:JC=

cos тс =■—

:34

следовательно, полагая

s = — х =

получим

г2 — ЛГ]Г + R* = О,

откуда

 

 

г =

у (*, + J / jc? — а 2) ,

причем точкам, лежащим вне круга, соответствует верхний знак перед радикалом. Тогда

Г ■*! + V х] —,Г* — „2 (Х1 ) ,

где а = 2R.

Формула (24) в этом Случае даст при ис — const

Х = 0,

й " Л 1 т(, _ И Ы Л'"

а

В более общем случае при следе,

совпадающем с отрица­

тельной осью S, получим

 

 

 

 

СО

 

СО

*1

X s i Ys — ip W

йI\ ^

 

Л=1

J

откуда

 

 

 

 

 

СО

00

 

 

* ;= р - з ?

 

 

 

 

Я

л -1

 

1

d%v

Коэффициенты ап и Ьп определяются обычным способом при конформном отображении внешности профиля на внешность круга.

Для тонкого, мало изогнутого профиля можно приближен­ но принять, что

* « 6i + X-

и, следовательно,

Тогда получим

0000

ал —1

3*

35

00

 

^ p V e^ ( - l )

n+i

 

ахл dx v

(25)

 

 

П = 1

 

 

a

 

 

 

 

oo

 

 

•*1

 

 

 

Ys = pVt aИ

\ 1

 

"I-

 

 

1 Л ? — ■

 

+ pVe ' 2 l { - l ) n+1nan j

-

H - 4

^ S - i j : i -

<r;). - + 4 r ; '

,25'a)

 

Л - 1

«

1

 

 

 

 

где (Ув)

соответствует

случаю неизогнутой пластины, для

которой

все а„

и Ьп равны нулю.

 

 

 

§ 5. Определение циркуляции

Скорость на профиле в проекции на касательное напрарление будет выражаться формулой

d®d0 _

d®,

d®2

d®3

J_

ds

d0 ds

d0

v,с d0

+ © d0 ^ 2tc^ d0

ds

 

 

 

 

 

М

Для того чтобы в особой

точке

преобразования,

для которой

Ц - равно нулю, скорость vs не

обращалась

в бесконечность,

необходимо найти Г из условия

d®.

d®t

 

Г = — 2тс \и,

dФ^

а?Ф2

(26)

с

dO

d0

+ “ ~ d r +

d0

 

где 6B соответствует заднему острому концу профиля.

Выше мы выбрали направление осей I и так,

что бв ока­

залось равным тс.

 

 

 

 

 

Если бы вихревого склада не было, то циркуляция опре­ делилась бы из условия

_ /

d®a

, d<b.о ,

(1Ф*

(26а)

Гk

2тс I

^0

“Ь ’Ос Ж~ ^

т~Ж

 

 

 

 

 

'“ «В

Таким образом, циркуляция Г при наличии вихревого следа складывается из циркуляции Г& и дополнительной циркуля­ ции Г', которую следует наложить на профиль для того, чтобы удовлетворить постулату С. А. Чаплыгина при обтекании про­ филя потоком, создаваемым вихревым следом. Как видно из формулы (26), эта циркуляция равна

Г = — 2тс

dd>s

db

 

36

Найдем сначала скорость ^ | - , получающуюся на окруж­

ности при обтекании ее потоком, который создается одним изолированным вихрем с циркуляцией §Г5. На основании фор­ мулы (20) имеем:

?£?1 Л I _ J______ Z _ \

где

Л

+

С - С ' , ] ’

 

 

 

 

r = Rel\ C' =

r V ',

 

Таким образом,

 

 

 

 

оГ^. 1

С-------- -------------- Г =

rfC

1 +

2я/ С

С2 — С(£' +

q) + C'Cj J

_ 8£,1 1 2га С 1

г' + р- — 2/? cos (0 — б7)

Вособой точке 9 равно тс и потому скорость там будет

равна

 

*£s I

 

R2

V' (тс)

,

/?2

 

R

 

/•

+

р + 2/? cos 0',

Для того чтобы скорость в этой точке была равна нулю, до­ статочно на окружность наложить циркуляцию

(рис. 12)

ЗП = - STs X

 

г'

Л2

X I

г’

R 2

Рис. 12.

г ' +

—,

-\-2R cos 6'

Тогда общая скорость на окружности будет

В случае вихревого следа с плотностью вихрей 70 получим

 

,

R*

 

 

 

 

 

r ~ F

 

 

 

 

 

г' +—г

1R cos(

1

da,

(27)

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

R2

 

 

 

 

 

г' + jr — 2/?cos91

 

r 4"

R12

 

 

da.

 

(28)

4- 2R cos ^0 —0^

 

 

 

причем интеграл берется вдоль линии вихревого следа в пло­

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

Г'2 = Г2 + V2,

Г' COS 0i =

 

da = d? | / l

+

(g ,)2.

Особенно простой вид написанные формулы примут в случае,

когда

мало

по

сравнению

с единицей.

Действительно,

тогда

 

 

daz^dt;,

cos 0j ps 1.

 

 

 

 

 

 

 

Замечая, что приближенно

для

вихрей в следе в плоскости z

можно

принять (опуская значки штриха у г')

 

 

 

 

 

Г +

^

 

. Я2

 

 

 

 

 

COS 0!: ■г А----

 

 

 

 

 

 

 

 

 

1 #*

 

 

и, следовательно,

 

_

х1 + угх1 —а2

 

 

 

 

 

 

(28а)

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

?

= V x i - a f ,

 

(28в)

получим

 

 

 

 

 

 

 

 

 

 

г

=

 

 

 

d x u

 

(29)

причем

интеграл

берется

вдоль

действительной

оси х г в пло­

скости

профиля.

Точно так же

 

 

 

 

 

®0(Я. е) =

_ 1_

 

 

 

-yfх{* аг

d x v .

 

iR

 

 

 

хг+ acos 0

38

Соседние файлы в папке книги из ГПНТБ