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

книги / Потенциальная теория неустановившихся сверхзвуковых течений

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
7.57 Mб
Скачать

Обозначим символом

решение (8.2.6), а

соответствующую ему функцию

давления,

 

= Ch

{Ц='=} + SliХф{,’- [V*} + М’Чг* {«*}•

(8.3.4)

По определению

это решение

удовлетворяет условиям

(8.2.3) и (8.2.5),

а также условию (8.3.3). Но функция

не может удовлетворить условию (8.3.2)

при

у* = 0 ,

поскольку производная ф{/*{^'*}

имеет здесь

особенность

вида y*-4i (как на дозвуковых передних кромках в линей­

ном приближении).

 

 

Рассмотрим

теперь функцию

 

f

= Sh X

т - ф1/.’ {V*}]

(8.3.5)

вместе с соответствующим ей потенциалом ф‘“’ (см. ниже (8.3.8)). Поскольку и ф‘” и ф”.’ удовлетворяют уравнению (8.2.2), ему должна удовлетворять и функция ф'-’. Диффе­ ренцируя выражение (8.3.5) по z и полагая Z = О, получим

Iz=O = Sh X [Ф^“ {P j-} -

Ф1У. {v*)h=o

(8.3.6а)

 

Iz=O = - Sh XК * -

V*.] = 0.

(8.3.6Ь)

Кроме того,

поскольку

ф‘^’ тождественно

обращается

в нуль при JC* <

0 и равно

нулю

при Z = 0 и у* < 0 , то

и потенциал ф‘*’,

который

нам пет необходимости опреде­

лять в явном виде, удовлетворяет условиям (8.2.5) и (8.3.2,3)

повсюду, за исключением у* =

0.

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

позицией функций

ф‘^

и ф‘®’:

 

Ф= Ch ХфУ.’

{V*} +

Sh ХФ"’

-1- М -Ч ”’ [ v % (8.3.7)

вместе с соответствующим ей потенциалом удовлетворяет системе соотношений (8.2.2,3,5) и (8.3.2,3), т. е. является

искомым решением (Майлс

1-“’ )). И так как ф‘^’ стремится

к нулю при у* = 0 как

то так же ведет себя и функ­

ция ф‘-’, как это и должно быть на дозвуковой задней кромке (условие Кутта — Жуковского). Потенциал ф^®’ в то же время не должен обращаться в нуль в районе, указан­ ном в соотношении (8.3.2), вследствие наличия следа (пелены).

Зная функцию ф, соответствующий ей потенциал полу­ чим интегрированием выражения (8.3.1) в виде

Ф (л *, t/*, Z, Г ) =

= SChx J Ф f/* - (Jf* — 5) th X, 2,

i)sch x^ | .

b

(8.3.8)

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

Г Л А В А 9

МЕТОДЫ ТЕОРИИ КРЫЛА МАЛОГО УДЛИНЕНИЯ

§ 9.1. Введение

Различные приближенные методы решения граничной задачи о крыле, основанные на предположении, что удли­ нение крыла очень мало, уже рассматривались нами в § 1.3 и были резюмированы в таблице 1 (случаи B^). Основной характерной чертой этих приближенных методов было предположение, что поток в любом поперечном сечении можно рассматривать как двумерный. Последнее предполо­ жение в свою очередь приводило к тому, что продольная координата входила в граничные условия лишь в качестве параметра, а следовательно дифференциальное уравнение, которому, должен удовлетворять потенциал, превращалось в уравнение Лапласа или волновое уравнение в двумерном пространстве. Полученная в результате граничная задача является классической и существует целый ряд разнообраз­ ных методов ее решения.

В этой главе мы рассмотрим применение теории крыла

малого

удлинения к

задаче о тонком заостренном крыле

малого

удлинения,

симметричном

относительно

оси х

(рис. 9.1), C

нестреловидяой задней

кромкой (последнее

ограничение,

как это

будет видно из дальнейшего,

может

быть в некоторой степени ослаблено). В качестве характер­ ной длины, которую в § 3.2 мы обозначали через /, можно принять максимальную хорду крыла с, измеренную от вер­ шины крыла {х = 0) до его задней кромки = 1). Из того, что крыло, по предположению, симметрично, его передние

кромки

могут быть

заданы функцией ± Ь (х),

подчиняю­

щейся

следующим

ограничениям:

 

 

 

Ь {х )= 0 {х )

X-^O

(9.1.1)

-Ь(х)

b'{x)> Q .

(9.1.2)

Первое из этих ограничений, (9.1.1), указывает,

что форма

крыла в плане является заостренной. Второе, (9.1.2),— ис­

ключает

наличие

стре­

ловидных

задних

кро-

мок. В

дополнение

к

»

условиям

(9.1.1,2)

оста­

ются еще II ограничения,

y=-b(x) /

указанные в таблице

1,

 

причем параметр а

(про­

 

порциональный удлине­

нию крыла X) определен так, как это показано на рис. 9.1, а. В част­ ности, следует помнить, что, для того чтобы ли­ нейная теория крыла оставалась справедли­ вой, как толщина крыла, так и безразмерный скос потока V должны быть малыми по сравнению с

Z

+Ь(х)

tV

удлинением X . В

то же

п

 

 

время в случае крыльев

 

 

 

 

 

малого удлинения мож­

Ф-0 Ф^- v

Ф-0 . ^

но ожидать более замет­

-/

 

^

ного проявления

нели­

 

 

 

нейных эффектов.

 

 

 

 

 

Как видно из рис. 9.1,

в)

 

 

принятая

в

этой

главе

 

 

Рис. 9.1. а) Вид в плане крыла ма­

система координат та же,

что и в § 3.3. Поскольку

лого удлинения; б) поперечное сече­

ние крыла малого удлинения; в)

координата

х

входит в

двумерная граничная задача в без­

граничную задачу лишь

размерной

форме.

D

качестве

параметра,

 

 

 

применение

модифицированного преобразования

Лоренца

не

сулит

нам

никаких

преимуществ. Области

S n R

определены теперь просто как \у\<Ь(х)

и |t/|>6(x)

соответственно. Область

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

рассматривается вследствие ограничения (9.1.2), наложен­ ного на форму задней кромки крыла.

Если не ограничиваться рассмотрением форм в плане типа, показанного на рис. 9.1, а, то поперечные сечения, пересекающие задние кромки, должны, вообще говоря, включать вихревую пелену (области Vl?'), где должно выпол­ няться граничное условие (3.3.7). Но тогда течение уже

Рис. 9 .2. Возможные формы задних кромок и метод их эквива­ лентного видоизменения.

не может далее приближенно рассматриваться как двумер­ ное, поскольку интегрирование названного граничного условия по X связывает потенциал в любой точке следа W C потенциалом в соответствующей точке задней кромки, лежа­ щей впереди по потоку ^). Так, на рис. 9.2 а, б потенциал

1) Крылья

малого удлинения с дозвуковыми задними

кромками

рассматривались

Уордом

Мэнглером

Майрлсом

[*^1], Леи-

делом

и Эйхельбреннером

 

 

 

в точке P J зависит от значения потенциала в точке P^. Если задняя кромка рассматриваемого крыла сверхзвуко­ вая, то мы приходим к парадоксу, поскольку ранее мы уже установили, что область не оказывает влияния на значе­ ние потенциала па крыле, если только задняя кромка послед­ него не имеет дозвуковых участков. Однако в рамках теории крыла малого удлинения сверхзвуковая задняя кромка

эквивалентна прямой (перпендикулярной к л) задней кромке. Иными словами, подсчитанный по этой теории эффект от вихревой пелены крыла, в случае сверхзвуковой задней кромки, будет порядка тех членов, которыми мы вначале преиебре1\'1и, а следовательно, и этим эффектом в рамках сделанных допущений мы также должны пренебречь. В практических приложениях теории крыла малого удли­ нения эта трудность может быть обойдена путем замены действительной формы крыла в плане, имеющей сверхзву­ ковые скошенные задние кромки, таким крылом, у которого область следа, расположенная впереди точек исходного К1)ыла, лежащих ниже всех других по потоку, заменена фиктивной поверхностью крыла (см. рис. 9.2 в, г, где гра­ ницы фиктивных поверхностей указаны пунктиром).

Для того чтобы нормировать двумерную граничную задачу, вытекающую из теории крыла малого удлинения,

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

координаты

I = - J -

(9.1.3а)

^ Ь[х) '

 

Ь{х)

(9.1.3Ь)

 

t

(9.1.3с)

X =

Параметры отнесения остаются теми же, что и в § 3.2, только величину Ь следует рассматривать уже как функцию от Xt а поскольку х является теперь лишь параметром, то геометрически задача сведется к задаче, показанной на рис. 9.1, б. Если теперь введем новый потенциал ф, свя­ занный C Cpсоотношением

Ф (х, I/, z,i) = b (х) ф (I, Tl, т; х),

(9.1.4)

ТО граничные условия (3.3.6,8) примут вид

Ij),! In=O=

-W

(1^1 < 1)

(9.1.5)

г1)1п=о =

0

(|^|>1),

(9.1.6)

а уравнение Лапласаа и волновое уравнение; в двумерном пространстве запишут/тся следующним образом:

+ "Фоп = О

(9.1.7)

 

+

=

(9.1.8)

Если имеет место гармоническая зависимость от времени ( § 3.10), единственным изменением (за исключением интер­ претации потенциала -ф как комплексной амплитуды или,

что эквивалентно, замене его на ф) будет изменение в урав­ нении (9.1.8),, которое перейдет в уравнение Гель.мгольца

 

+

+

=

(9.1.9)

в котором

 

 

 

 

 

l = kMb{x).

(9.1.10)

Отметим, что

из условия

& = 0(Х) в случаях В[,

В'[ и В”

(из таблицы

1) вытекают

требования X = O (I),

= P ( I )

и Я » 1 соответственно.

В дальнейшем следует помнить, что операция ди()х}1еренцирования по х не остается инвариантной к преобразова­ нию (9.1.3). Поэтому перепад давления (3.3.5) преобразуется

в этом

случае к

следующему виду:

I = 4

+ 6' (ф -

— !!фп — тфх) + M‘ 4r]z=o*. (9.1.11b)

При гармонической зависимости от времени выражение

(9.1.1 Ib) переходит

в выражение

 

 

/ = 4 [6ф, +

6' (1> - U i - Т1фп) +

(9.1.12)

Теперь мы приступим к рассмотрению возможных методоврешения граничной задачи, формулируемой уравнениями (9.1.7,8,9) C граничными условиями (9.1.5,6). Здесь не будет рассмотрен приближенный метод, именуемый теорией поршня (поскольку, строго говоря, он находится вне рамок теории крыла малого удлинения), однако полезно привести весьма простое решение, которое дает этот метод;

Ф =

V.

(9.1.13)

§9.2. Решение уравнения Лапласа

Вэтом разделе будет рассмотрена граничная задача, формулируемая соотношениями (9.1.5—7), что соответ­

ствует случаям и B.j таблицы 1. Это классическая задача теории потенциала н существует ряд методов ее решения. Из них мы рассмотрим метод конформных отображений (А), и метод, основанный на введении эллиптических коорди­ нат (В). Выбор их обусловлен тем, что в дальнейшем оба метода будут использованы нами в случаях соответственно пространственного тела и высокой частоты колебаний (см. главу 12, § 9.3 и 9.7 и главу 11), В дополнение к этому будет получено и решено соответствующее задаче интегральное уравнение (С), так как этот метод может быть распростра­ нен как на случай высоких частот колебаний крыла малого удлинения (§ 9.4), так и на случай крыла квазималого удли­ нения (§ 9.8).

А. Разрез Tj = Oib, |1|<1

в комплексной плоскости

^ = |4-/Tj отображается на

круг единичного радиуса

в комплексной плоскости Ci с помощью преобразования Жуковского

£ = • !■ « ,-I-C:*)-

(9.2.1)

В полярных координатах

=(9.2.2)

граничные условия (9.1.5) и (9.1.6) примут вид

Решение уравнения Лапласа, удовлетворяющее усло­ виям (9.2.3,4), после разделения переменных получим в сле­ дующем виде:

 

 

 

 

 

 

 

 

 

 

(9.2.5)

 

 

 

 

^ tJ (cos 0) sin (rtO) sin 0 dO,

 

(9.2.6)

a функция

V определена уже

как

функция

переменной

^ = COsG,

а

не у.-

 

 

 

 

 

 

 

Если

рассматривать условие

(9.2.3)

как справедливое

на всем единичном круге (—я, -Ья), то интегрированием

формулы Пуассона (Батеман,

стр. 236 и след.) получим

решение в

виде

 

 

 

 

 

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(9.2.7)

Тождество выражений (9.2.5) и (9.2.7) можно установить,

разлагая (9.2.7) в ряд Фурье,

или

просуммировав

ряд

(9.2.5).

 

 

 

 

 

 

 

 

 

 

В .

 

Переход к эллиптическим

координатам р,

v осуще­

ствляется

преобразованием

 

 

 

 

 

 

 

 

 

 

 

C = ch (p -ftv ),

 

 

 

 

(9.2.8)

которое отображает разрез Tl = О,

[| | < 1 на область

р = 0,

— я < v < + я .

Граничное условие

(9.1.5), рассматривае-

т) Граничное условие (9.2.4), введенное

первоначально

исходя

из соображений симметрии, теперь, по сути дела,

является

излиш­

ним,

поскольку условие (9.2.3) можно рассматривать

как

справед­

ливое

в области ( — я , + я ) , где интервалы (О, я)

и

( — я, 0)

соот­

ветствуют

верхней

( T i = O - I - - U K l )

« нижней

(ii =

0 — ,

| Ц < 1 )

сторонам

крыла соответственно.

 

 

 

 

 

 

 

мое как

справедливое на

= О ± , переходит в следующее:

 

^ H IH= O = — г'sin V

 

( — n < v < + я ),

(9.2.9)

и тогда искомое

решение имеет вид

 

 

 

 

 

IjJ=

2

rt"^a,^e“ ”»^sin(rtv),

 

(9.2.10)

где а,^задано выражением (9.2.6).

 

 

 

С.

Решением

уравнения

(9.1.7),

которое

обращается

в 1|'о(^)

ври ^ l= O + .

является

выражение

 

 

 

 

11)- — И

 

 

 

 

(9.2.11)

Что (9.2.11) удовлетворяет уравнению Лапласа,

легко

убедиться непосредственной подстановкой его в (9.1.7);'

решение вида я"Ч1[(^— |')* +

TI-J"1 является

потенциалом

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

(^', 0). То обстоятельство,

что решение при т| = 0 +

перехо­

дит в ijjd , о -1-) ='фо(1). можно проверить,

если ввести заме­

ну переменной

= ^ + Titg О и,

положив ii =

0-f-,

выпол­

нить интегрирование по 0, или учитывая соотношение

Пределы интегрирования по

определены условием (9.1.6).

Теперь, подставляя выражение (9.2.11) в граничное усло­ вие (9.2.3), составим интегральное уравнение для to. Резуль­

тат

можно

представить

в форме

 

 

 

 

=

|||<1,

(9.2.13)

где

L есть

интегральный оператор

следующего вида:

 

 

 

{ I T ^

S (I-

P)"--ьТ1^}

Интегрируя по частям, с учет9м требования, чтобы to (± 1) равнялось нулю, и выполняя указанные операции, получим

i i o i r )

d^'

91'

(S-^ ')

Соседние файлы в папке книги