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

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

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

более простую форму, = у), которую мы имели выше. При таком уравнении отраженной линии Маха не удается использовать теорему свертывания, с помощью

которой доказывается

справедливость выражения (7. 4.6).

Приведенные выше результаты первоначально были полу­

чены (Майлс

П^®®1) из представления граничной

задачи в виде интегрального уравнения Винера — Хопфа. Численный метод, основанный на теории конических тече­ ний, был дан в работе Экыома [Ч. Частный случай, когда заданный скос потока не зависит от г/, исследовался также

Гудменом

1®-J, использовавшим метод Гарднера (§4,3),

P OTTOM

 

решавшим задачу методом,

напоминавшим

метод,

примененный Лэмбом [^^“],

к

исследованию

задачи

дифракции

на полуплоскости (§ 7.2)

и, наконец,

Уоткинсом

["”®],

получившим несколько первых

чле­

нов разложения решения в ряд по степеням частоты.

Не

полностью законченное исследование задачи дано в работе Галина [’®], который воспользовался разложением в ряд Фурье в методе Магнарадзе (§4 .2). (См. также работу Стыоартсона [^’®]).

§ 7.5. Прямоугольное крыло конечного размаха

Решение для потенциала в окрестности угла прямо­ угольного крыла можно обобщить на целое прямо­ угольное крыло в том случае, когда отсутствует взаимодей­ ствие боковых, кромок. Рассмотрим крыло, изображенное на рис. 7.5 в системе координат {х\ у). Любая точка зоны I лежит впереди линий Маха, проходящих через передние угловые точки, а следовательно, здесь не сказывается влияние боковых кромок. Точки зоны I I лежат ниже по потоку от линии Маха, проходящей через левую угловую точку, но выше такой же линии, проходящей через правую угловую точку, а следовательно, на них оказывает влияние лишь левая боковая кромка. Точки зоны Г Г, наоборот, подвержены воздействию лишь правой боковой кромки. Наконец, точки зоны I I I лежат ниже по потоку от обеих указанных линий Маха и, таким образом, захвачены влия­ нием обеих боковых кромок, действующих, однако, незави­ симо, при соблюдении условия с ' < 2 6 . Таким образом, результаты, полученные для изолированного угла, позво-

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

Рис. 7.5. Прямоугольное крыло отно­

сительного удлинения Я, = 26/с; пока­ заны области влияния; только передней кромки (У), одной нз боковых кромок (/у H УУ')н обеих боковых кромок (У/У).

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

обтекаются теперь

не независимо.

 

Принципиально влияние этих отражений может быть

найдено методом

Шварцшильда

разработанным им

для исследования

задачи о дифракции

на щели и при­

мененным в работе Гана [®*] к случаю стационарного обтекания. Иначе решение нестационарной задачи для'этого

случая может быть

получено непосредственным примене­

нием метода §4 .2

к

существующим (см. работы Гана или

Лагерстрома

решениям стационарной задачи. Но полу­

ченные таким способом результаты могут оказаться слиш­ ком громоздкими для того, чтобы найти себе практическое применение, поэтому более целесообразным здесь будет, по-видимому, непосредственное решение задачи (7.1.1—3)

вэллиптических цилиндрических координатах (см. главу 11).

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

по размаху, взятых с соответствующим весом. Кроме того, из соображений симметрии (функцию v, а следовательно, и ф всегда можно представить в виде двух компонент — симметричной и антисимметричной относительно у = Ь) можно рассматривать интегрирование не от О,до 2Ь, а от О до Ь. В задачах флаттера и устойчивости обычно можно принять, что распределение скоса имеет вид

V {x’,y ) = ( b - y r V , { x r

(7.5.1)

в соответствии с общепринятой практикой представления отклонений крыла в форме полинома. Исходя из этого, рассмотрим интеграл с весом следующего вида:

О

(х') = б-*"-"'! ^ (6 - уГФ^ {х\ у , о+ )dy, (7.5.2)

где Ф„ — потенциал, соответствующий распределению ско­

са, заданному уравнением (7.5.1), а множитель

вве­

ден для того, чтобы нормировать интеграл.

 

Подставив функцию (7.5.1) и потенциал (7.4.5,6) в выра­

жение

(7.5.2),

запишем

интеграл

в

виде

 

 

 

 

 

X'

 

 

 

(^') = ( т +

п +

1 )-^

W o (^' -

 

1)] Vn И) d l -

 

X'

rt/2

 

 

 

 

 

- 1

(I) d i H J ,

[X [X' _ I) Sin 01 ^

X . . » - 1) “ 5 61 de,

 

 

 

 

 

 

 

(7.5.3)

Хтпп(*^) = 1 A -m -n - 1

2

 

 

 

 

y b - y ^ d y

J

 

(7.5.4)

Первый член в выражении (7.5.3) является, по суще­

ству,

двумерным,

(т. е. вытекающим из теории

полоски),

а оставшийся интеграл может быть представлен в аналогич­

ной форме после выполнения интегрирований по г/, т)

и Q..

Область интегрирования по у

ц для функции

Xnm

показана на рис, 7.6, а. Чтобы упростить пределы, введем

юз

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

Подобласти 5 i и 5г (рис. 7.6, а) преобразуются теперь

Р и с . 7 . 6 .

а) О б л а с т ь и н т е г р и р о в а н и я

в в ы р а ж е н и и

( 7 . 5 . 4 ) ; б) об»

 

л а с т ь и н т е г р и р о в а н и я в в ы р а л с е н н и ( 7 . 5 . 6 ) .

В подобласти,

показанные на

рис. 7.6 б,

и выражение

(7.5.4)

примет

вид

 

 

яъ

5{Ь-уУ^{Ь-уЛ-исоъфУ^<1у.

о^!+COS

(7.5.6)

Этот интеграл можно сразу разложить в ряд по степеням (и/6), для того чтобы получить для функции Xmn(«) выраже­ ние вида

 

m + n + l

 

х«»М =

S

“Г ’ Г г ) "

 

h=0

 

Общее выражение для

ajy^>

в функции параметров т, п

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

Воспользовавшись выражением (7.5.7), можно выпол­ нить интегрирование по 0, учитывая, что имеют место соот­ ношения

я/2

' ' • = ( т У \ ^ ( ’« 's in e ) ^ (c o s e )''d e . (7.5,8а)

О

-------^ ( l ^ + о

(Л > 1).

 

(7.5.8Ь)

где (7.5.8а) вычислялся как интеграл типа СоЕшна (см.

Магнус и

Оберхеттингер

30). В результате, учиты­

вая, что /о = о,

для функции

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

выражение

вида

 

Ч^„,Л^') =

( т +

п + 1 ) - 1 х

 

 

 

771+71+1

X 5 J J x ( X ' -

1)1 ( !) d | + 2

“Г " Т ( I -I-1 )

о^l=l

X*

 

X 5 (X' - |У'«У,„г Ix (X' - 1)1 V„ (i) d l.

(7.5.9)

о

 

,в простейшем случае, п = О, представляет интерес

просто

интеграл

от потенциала

при m =

0.

Полагая

в (7.5.6)

т = п =

О, получим

 

 

 

 

 

X o o (W )= I -- ^

 

(7.5.10)

Подставляя

в выражение (7.5.9) и вводя триго­

нометрический эквивалент функции

Бесселя

7i/„:

 

 

 

 

 

 

(7.5.11)

получим

 

 

 

 

 

Ч'оо=

j { + l x ( x '- e ) l - 5 i i l J ^ ^ = B - ‘}

 

(7.5.12)

В частном случае гармонических колебаний, преобра­ зуя выражение (7.5.12) к физическим переменным с помощью соотношений (3.10.6,7) и х = рх' и вводя соответствующую

переменную интегрирован

будем иметь

ъ

Фо(^> y ,0 + )d y =

о

л*

= $g<,.(*-5)».(6)<'5. (7.5.13)

 

(7.5.14)

Л = 2рб=рХ

(7.5.15)

есть эффективное удлинение крыла.

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

Из (7.5.14) видно, что чем меньше приведенная частота колебаний k, тем большее значение приобретает влияние удлинения крыла, поскольку первый и второй члены в фигур­ ных скобках выражения (7.5.14) с ростом k убывают соответственно как и В пределе, при очень малых k, пренебрегая в выражениях (7.5.13) и (7.5.14) членами по­ рядка kr, получим

?о о (^ )=

[ 1 - ^ ( ^ - 5 ) ] 4 ,(1 ) « .

 

(7.5.16)

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

^) Результаты

численных

расчетов см. в работах Майлса [^®®]

и Нельсона, Райнн

и Уоткинса

Случай элерона был рассчитан

Майлсом [!” ] H Берманом [^®]. Обширные расчеты были проведены

для де(|)ормнруемых крыльев Брандстаттером и Моршки (®*J.

§7.6. Продольное демпфирование

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

(ср, § 5.3), так что

O(X) = а [1 + Ik (X -X o )].

(7.6.1)

Поскольку пдследовательность операций

дифференци­

рования по х и интегрирования по размаху (у) несуществен­ на, безразмерная подъемная сила полосы, занимающей половину размаха (протяженностью, равной Ь), согласно

соотношению (3.3.5), равна

 

 

 

6L

«Лф).

(7.6.2)

 

:46(ф^ +

Подставляя в

эту формулу выражение (7.5.16) и прене­

брегая членами

порядка 0(^*),

получим

 

f

= к { » > ) - j [ 4 +

f - ^

( - » ]

<^5}.

 

 

 

 

(7.6.3)

где

QQ есть производная

коэффициента

подъемной

силы

.по углу атаки для двумерного крыла (профиля). Подста­ новка зависимости (7.6.1) в выражение (7.6.3) и вычисление

коэффициента момента

по формуле

 

Сл1 = 4

\ (Xo-X) 3L dx,

(7.6.4)

дают нам следующие выражения Для вращательных произ­ водных (ср. §5.3):

^Afg= — Ao { [ 12+ ( " 2 “ “ 2А

Границу области устойчивос:ти относительно продоль­ ных колебаний с одной степенью свободы определим из

уравнения {СMq + C Mi) = О, которое сведется к следующему

бр*( 2 А - 1 ) X q = S (2р*- 1 ) Л - 2 ( р 2 _ 1 ) ±

± Y W + ^ / (З Л - 2 ) 2 - р 2 (6 Л 2 _ 2 ) . (7.6.7)

Области устойчивости для крыльев с удлинениями, рав­ ными Я,= 2, Я, = 4 и Л =оо, показаны на рис. 7.7.

В плоскости (к, М) граница области возможной неустой­ чивости, при любых положениях оси, получится, если

Рис. 7.7. Границы устойчивости прямоугольного крыла относительно .низкочастотных продольных колебаний [см. уравнение 7.6.7)]«

потребовать, чтобы два значения Xq, определяемые уравне­ нием (7.6.7),-сливались. Соответствующие значения удли­ нения определятся из уравнения

I — (RP+ 1) ■L '" к ^ — 6---------

J

•Этот результат показан в виде верхней ветви (Л > 1) на рис. 7.8. Соответствующие результаты для Л < 1 (§11.3) !Представлены нижней ветвью. Видно, что неустойчивость

невозможна

как

при у^1линениях Л, < 1,93 для Л > 1 ,

так

и. при Л <

0,65.

Здесь, однако, следует лодчеркнуть,

что

аппроксимация по А ухудшается с.при($лиженйем числа M

к единице, поскольку (кроме ограничения Л > 1) в выра­ жении (7.5.16) мы пренебрегли членом порядка (Ш-/р^)2 (см. также замечания в конце § 11.3).

16.

15

^A=^O

*А~2,0

1.2

S-AHO

0,1 0.2 0.3 0.^_ 0.5

J/JL

Рис. 7 .8. Граница устойчивости относи­ тельно продольных колебании вокруг произвольной оси (см. уравнение (7.6.8)

при А '> 1 и § 11.3 при А < 1).

Влияние конечной толщины прямоугольного крыла на его продольное демпфирование исследовалось Экыомом [®], который использовал результаты двумерной теории, аналогичные рассмотренным в § 13.3.

§ 7.7. Случай произвольной зависимости от времени

Если скос потока v является произвольно заданной функцией времени, то потенциал на верхней поверхности прямоугольного крыла может быть получен как обратное преобразование Фурье от выражения (7.4.7). Прежде всего, найдем обратное преобразование от потенциала, получен­ ного по теории полоски (7.4.5); оно имеет вид (ср. с (5.4.2)):

t'Mx'-l)

v'il. у, Qdtj

(7.7.1)

 

В ТО время как обратным преобразованием от потенциала стационарного обтекания (7.3.6) будет выражение

Фо

1

у' (^. Л. П dr\

(7.7.2)

 

 

 

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

Рассмотрим сначала обратное преобразование выраже­ ния (7.4.7). Учитывая соотношения

 

 

 

 

(7.7.3)

J

Уо(хи)е1«''^х = 0

( U l > « )

и применяя

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

 

 

2

( * '- E ) s ln e

 

 

b

о

- (* '-1)1Sln■ о к(Ж' -

6)= Sln2 е-С2

2 arCBln V M *' - 1) сое О

 

Х Ж

 

S

 

(7.7.4)

Если для получения обратного преобразования выраже­ ния (7.4.6) воспользоваться соотношениями (7.7.3) и затем результат проинтегрировать по частям, потенциал будет

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