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

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

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

§ 6. Окончательное выражение для сил

Из формулы (29)

00 00 _______

j

r(* i t ) V ^ r a dx^ T ''

а

а

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

со

Хг

 

s = рИс Т /1

d x x

] / х\ а2

а

 

 

 

или же

Ys = риса Г

-- — рГ'ис.

J

V ^i- « 2

а

 

Замечая, что первое слагаемое в общей формуле для силы Ys, т. е. слагаемое, дающее силу Н. Е. Жуковского, выражается формулой

 

 

У) = рГив = рГА« с + рГ,яе,

(3°)

где

Г£ — квазистационарная циркуляция,

получим, что

 

 

 

f

^dXx

 

 

 

Y ; -f- Ys — рГfcUc -)- рticd J V x \ — a1

 

 

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

окончательно

выразятся для

пластинки слагаемыми

 

 

 

* в = - р Г Ч ;

Y,

я2

 

 

 

 

 

 

 

Для случая тонкого, мало изогнутого профиля, для которого

приближенно

выполняются

равенства

(28в) и

(29), будем

иметь:

 

 

 

 

 

 

X s =

X's - 9Tfvc,

 

(30а)

 

 

r s = p«ea

+

 

(ЗОв)

 

 

J

У х \— сР

 

 

где

и A Ys даются формулами (25) и

(25а).

 

30

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

X Xjk Д Х,„ + Х т Д Xi Д X s

(30c)

У = yjk + У*>+ Гт + Yi Д r s

причем Xjk и Yjk определяются, согласно теореме Н. Е. Ж у­ ковского, для квазистационарного значения циркуляции.

Для остальных составляющих выше нами были найдены соответствующие выражения. Относительно Х а и Уш следует заметить, что в них содержится влияние следа. Действительно, так как

то

d-Ф=

-f- аг’ ,

 

 

 

 

 

 

 

 

Х ш— рю J X (d<Y>k Д ^ФS Д ^Фг')] == X w k Д X m s

Д Х т Т ',

к „ = .

рш Jy; {йФь Д <1Ф$ Д dФг’) —

Ymk Д

Д УшГ'

В комплексной форме будем иметь:

 

 

 

Х шД

i У'ш— — р<о J 2

 

d z - Po>j dzFdszdz. (30d)

 

L

 

 

L

 

 

Пользуясь

разложениями, полученными выше,

для силы

X s — iYs, будем иметь выражение

 

 

 

 

 

 

оо

 

 

 

 

Д iYsm=

 

+

e~inK)

 

(30e>

 

a

 

n=i

 

 

 

В частности,

для пластинки,

для

которой

Сп — 0

и

 

 

Г = y (jc,

Д У X I — а2),

 

 

получим

 

 

 

 

 

 

х вш= ро)J (хг — v х \ —a2) d x x;

Ysw —о.

 

Заметим, что все силы, зависящие от следа, требуют знания плотности вихрей следа т(Д . 0 -

§ 7. Распределение давления по поверхности профиля

Так как конформное отображение внешности круга на

внешность профиля известно, то определение скорости V'(u', v')t которая получается на профиле при обтекании его потоком, создаваемым вихревым следом на профиле, будет производить­ ся по формуле

40

V' = V°3F>

Rd%

где v 0 дается формулой (28). Скорости возмущенного потока около профиля всегда можно представить в виде

u = uk + u'-, v = vk + v',.

где первые слагаемые соответствуют движению с „квазипостоянной“ циркуляцией, а вторые — влиянию следа и добавоч­ ной циркуляции.

Подставляя эти выражения в формулу (10), получим

P = P k +

Р1(вс—иу — Щ )«'+(v c +сох—v k) v'] —

 

u'2 + v '2 ■ (дФ

d r \

 

_ р ---------------- - р у dt +

dt ) ■>

гдеPk есть давление, соответствующее движению с квазипостоян­ ной циркуляцией.

Так как нормальные составляющие Vn равны нулю, то

I

\7>

/<5Ф . dTk . d V \

/01Ч

p — P k + p W k V — p ^ -----Р \-Ш + Ч Г + ~ м ) ’

(31>

где Wk — мгновенная относительная скорость, которой соот­

ветствует давление рь- Эта скорость равна

Щ =у { Щ — «У—щ У +(Vc +сох-v ky .

Определение величин, отмеченных индексом k, известно..

Величину (?Ф лучше

всего определять по формуле

 

 

о

 

причем Ф5(0, t) равно нулю.

всегда

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

если известна величина т (5)> то

можно найти V0(0, £)

и затем построить распределение

давле­

ния по профилю.

 

 

§ 8. Вычисление моментов

Вычисление моментов относительно начала координат, дей­ ствующих на профиль, следует производить по формуле

M = \ { x d Y - y d X ) ,

1

где

d X = — pdy\ dY — pdx

41,

и р дается формулой (5), так что

dX = — px dy — р [{ис — °>у) и -f- (vc + шх) v\ dy +

d Y = + pndx +

р [(ис — шу) и-f

(^с +

<*>x)v] dx

причем через

агФ

 

й

 

выражение

 

 

обозначено

 

 

ЯФ

+

.

.

.

+ “Ф»-

 

 

ГФ4 + исФ1+

Использовав,

как и

при вычислении сил,

формулы (6а), мы

-получим

d X === — px dy — р(ус + шх) 8ф +

 

+ -j [(v2 — и2) dy +2tivdx] dy + р д- ^ - dy.

 

dy = + P «dx +

p («с — <»y)8Ф +

+ y

[(®2 ~ и2) dx 2uvdy\ dx — p

dx.

Отмечая, что выражения

 

 

 

 

 

dP = 2иvdx + (v 2 — и2) dy,

 

 

dQ =

(v2— и2) dx — 2uvdy

 

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

dP — idQ = i ( u — iv )2 (dx + idy) = i (‘^ X d z ,

для элементарного момента dM получим

dM = x d Y —y d X = рх (xdx + ydy) + p {ucx + vcy) -f-

+

( x d Q - y d P ) - ? d- ^ - { x d x + y d y ) .

Ho x d Q —ydP

есть не

что иное, как действительная часть

выражения

 

 

.и потому

 

 

dM = ръ {xdx +

ydy) + p (ucx + vcy )

 

 

I ^ ( x d x + ydy),

42

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

М = рJ(исх + vcy) 8Ф — ^ Re J^ J z d z

L

f

£

 

-

Р [

+ У^У)-

(32)

Так как 8ф можно представить в виде

 

8Ф = 8ФЛ+ 8Ф5 + 8ФГ|

 

‘где ЗФ^ соответствует

квазистадионарному

течению, а

ЗФЯ и ЗФг' — течению

от наличия вихревого следа, то первое

слагаемое формулы (32) представится в виде

 

р J (исх + vcy) 8Ф =

р j* (исх + v cy) о (фА+

Фг )+

£ J(ы-сх +^сУ)

£ =7М ^ +Мт> +M s*.

£

 

 

 

Величина

 

 

 

Щ = Р j

(«с* + vcy) ЬФк

(33)

i

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

5Ф& = мс3ф, + дс8Ф2 + со8ф3 -(- Г^8Ф4.

Мы умеем определять на контуре круга Фи Ф2, Ф, и Ф4 как функции 6, поэтому определение Мь не представляет собой труда, если мы знаем д:(9) и у (б) согласно формулам (8а).

Так как

 

8 Ф г = [г > ,

ТО

 

 

M r = р £ [ К * ^

^ 88 = Р (мс«о + Vcbо) Г'• (34)

о

 

Это выражение есть момент силы рУсГ', проходящей через точку с координатами а0 и Ьй, относительно начала координат.

Момент M s*, выраженной формулой

M s* = p j ( u ex + Vcy) §Ф5,

L

может быть несколько преобразован. Действительно, так как: функция 4^5 на контуре профиля постоянна, то можно написать

М s* = 1 j (исх

\ dFs dz = pRe

где

Ve* = uc — ivc-

Переходя в плоскость С, будем иметь:

M s* =pRe

dFs(Q

z (С) dC

 

Л

 

где на основании формулы (20)

^ ©

1 (*„ (1

1

Л

2wJ ‘0,Д

С—С'

 

Б

 

причем, как и раньше, о' есть длина дуги вихревого следа, проведенная из начала этого следа в точку, где помещается вихрь с циркуляцией -р, do' (рис. 11); С'— комплексная коор­ дината точки N ' . Так как

 

z(C) = C+ ^

+ 2

^ -

 

 

 

 

 

 

п

 

 

 

ТО

 

 

 

 

 

 

 

 

 

 

 

 

 

V

- T

с

- b l x

 

I s

Lfe

с - с '

 

 

 

 

 

 

 

 

 

* ( ‘ + т + 2

Н Л da

 

 

 

 

 

л = 0

 

 

 

 

 

Производя

интегрирование

по

переменному С с учетом:

формул (23),

получим

 

 

 

 

 

 

 

M s*> - - R e ( p V t*jTr,

f

+

 

 

 

do'

 

l

 

L

 

n~0

 

 

Или, так как

 

 

 

D2

 

/?2л2/6'

 

c' = pV"',

 

=

 

 

 

 

 

с'

 

 

 

 

 

 

 

 

 

 

 

4 2 +

q =

£ ( i +

<?*'),

 

 

44

то, опуская значок штриха, будем иметь:

 

M s*= — Re

f ( 1 + ^ ' )

+

V £л

 

(35)

 

 

 

 

 

 

 

я—О

 

 

Если вихревой след профиля таков, что он при конформном

отображении внешности профиля на внешность

круга

пере­

ходит

в отрицательную,

действительную

 

полуось

?,, то

г =

= о ==

И

 

 

 

 

 

где ^ = |С|.

В частности, для пластинки сп = О, Vc* = ис и

 

Ms* = P«cj \ o f - 2^ , .

(37)

R

 

Вэтом случае, переходя в плоскость z, имеем:

ипотому для точек следа

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

00

(38)

M s* = pucj f { x i — V x 12— a2) d x v

а

Переходим ко второму интегралу формулы (32), а именно интегралу

L

Рассмотрим сначала бесконечно тонкий профиль, т. е. дугу. В этом случае можем написать

Mt = - -f Re j \{vt - ivn)\ - (vt - ivn)2H\ ze~m dz,

L в

причем интеграл берется по верхнему разрезу. Замечая, как и раньше, что

= Тп. «/* = 4

v“* = VnB = ®я"»

45

получим

Mi = — р Re J 'и, (vt* ivn*) e~ibzds1 =

LB

= — p R e .^ t„ (u * — iv *)(x + iy )d sl,

LB

откуда

Mi = —p| f„ (m** + v*y)

LB

или векторно

Mi = p j [г X ( V"/* X Tn)J dsu

LB

где Vi* = Vi + Vn есть индуктивная скорость, вызываемая на профиле системой свободных и присоединенных вихрей. Можно убедиться, что

Г X (УпХТп)] dsl =°>

LB

имы окончательно получим

Mi = p J

X ( Vt X

То)] dsi = j

(г X dRi) ds,

(39)

L B

 

 

L B

 

 

т. e. момент М{

есть главный

момент

элементарных

индук­

тивных сил

 

 

х Tnlrfsi-

 

 

d R i =

p (V{

 

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

Переходим, наконец, к третьему интегралу формулы (32). Этот интеграл имеет вид

-(x d x + y dy) =

L

= — Р I ( « А + ъсф2 + 0)Ф3 + ГФ4 +

(xdx + ydy). (40)

Из этой формулы видно, что будет существовать момент М т, независящий от влияния следа

з

М т = — Р ^ Uk J Фк (xdx + УаУ)

k —1

L

46

где

Wj — UCi U 2 V c , U 3 — о).

Так как Ф^, как функции полярного угла б, мы вычислять умеем, то вычисление момента М т не представляет труда, если мы знаем л;(б) и у (6) в результате конформного отобра­ жения. В частности, для эллипса имеем:

x = a c o s6, y = 6 sin6,

x d x + ydy = ai ^ bi sin 26^6

и, кроме того,

Ф! = — bcos6, Ф2 = — a sin 6, Ф3 = — д2 . 62 sin 20

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

Mm = f ( a 2 - & 2)2<“ .

Наконец, из формулы (40) вытекает выражение для момента,, связанного с наличием следа,

М г' = — Р f (ГФг + Фв) (xdx + ydy) =

s L

a:

1

da',

(41>

 

 

где

 

 

 

 

2 — x — гу,

 

 

F(C ')= j

In — U

d ( z z )

 

 

 

C-Ci

«

 

 

 

 

1

2niz0z0.

 

 

 

C-Cj

 

 

 

 

 

Так как

 

 

 

 

2 (С) - С+ ^ + ^

,

 

 

 

л—0

 

 

 

 

я - 0

 

 

47

то

гг = СС+ 4 2С+ V ??- С+ ~ С+ 4 ' +

 

С

Сл

с

сс

 

 

Л — 1

 

 

 

 

 

\

1 VI СУС|Х

+ '?!2 й + 2 гг‘с+А“Ё ^ + £ J Z j

Кроме того, К = /?2 и,

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

 

 

» - 2Й* + С! +

 

^ й ? "* ц +

 

S г»С<" ' 1’ +

 

 

Л — 1

 

Л =» 1

Если это

выражение помножить на

(f^

1

 

1

Л rfC и по-

tt

 

 

 

 

 

 

 

^

 

 

t - k ) '

 

 

лученный многочлен

проинтегрировать

вдоль

окружности К,

то используя

формулы,

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

С 1 d l

 

2ш

ГМ

rfC

 

n

 

 

 

 

 

J слс —t' —

С'л ’

J t-c-c;

'

 

 

 

 

при « > 0,

опуская

знак

штриха и принимая

ca = z0,

мы

получим

 

 

 

 

 

 

 

 

 

 

 

™ _ - « ( ч + £ ) + Ц - л—1 -р*+л=2 ^ -

 

 

 

 

 

 

2

 

 

2

 

 

 

 

2 » » («+,+

 

2 2

c"cv- (

1

_ i _

И-

,

(42)

 

R2(1 I С”-11

 

 

 

Л — 1

 

 

V

 

 

 

 

 

 

 

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

при v > jx,

 

 

В результате для М / получим

 

 

 

 

 

 

 

 

М

где через Д обозначено выражение, стоящее в формуле (42) в квадратных скобках. Так как для точек, лежащих на следе,

С = гел , Cj = ™ е ib,

48

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