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

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

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

для вычисления давления, задаются выражениями

 

 

1

 

 

 

 

 

 

: ^

= 2

^ 1|,<»-М| = 2 т - - ^ т =

(9.7.10а)

 

 

 

X<c“>= 2 - t .

 

(9.7.10b)

Если

произвести

наложение влияний

обеих

кромок

для 1 <

т < 2

(см. § 7.5 и

рис. 7.5), то выражения (9.7.10)

оказываются

справедливы

во всем интервале 0 < г < 2 .

Граничная задача для потенциала

получающаяся

в результате подстановки выражения (9.7.9) в граничное условие (9.7.7), может быть решена как методом Эвварда (см. § 7.3), так и путем наложения конических течений. Выкладки, ведущие к цели, несколько длинны, но окон­ чательный результат, т. е. интеграл от потенциала по размаху и его первая производная по времени (см. Лагерстром [^^°]), имеют вид

= 4

-

(т^ - о

1

+

 

 

 

 

 

 

 

 

 

 

 

 

 

+ i (

1 + 4 т ) (Т»- 4)*+ (I - т) arch ( 4

т)|

(9.7.11а)

И

 

 

 

 

 

 

 

 

 

*"’= 7

{ ( T

- о

T + т

У

(

т

^)}

 

 

 

 

 

 

 

 

(9.7. Ilb )

Величина

 

обращающаяся в нуль при т =

2, должна

быть сложена с х‘“* при 2 <

т < 4 .

 

 

 

 

Переходя к асимптотическому разложению для случая

больших X,

но сохраняя

предположение,

что зависимость

скоса потока

v от времени

задана

соотношением

(9.7.8),

рассмотрим

следующую

граничную задачу для

преобразо­

вания Лапласа:

 

 

 

 

 

 

 

 

 

1 И | ,_ о = - а - ‘

( I I K l ) ,

 

(9.7.12)

 

 

 

r U - o =

0

(| | | > 1).

 

 

(9.7.13)

 

 

 

1|>й +

'И а -а > * --= 0 .

 

 

(9.7.14)

Очевидно, что решение этой задачи может быть получено из результатов § 9.3 простой заменой параметра X на —is

(посколь1

заменено

на е” ). Подставляя

в

выражение

(9.3.3) V — S"^

и используя

разложение (9.3

.7),

получим

а,^ =

^

^

(п — нечетное),

(9.7.15)

откуда (так как при четных /г коэффициенты

обращаются

в пуль) из (9.3.5) для ф* следует выражение

 

 

 

H)»(cosv) =

s - i 2

 

 

 

 

 

 

 

«=0

 

 

 

 

 

 

х

(

^2.4.1 ( - « ) s e 2 „ , i (^v, - J S ' - * )

 

(9.7.16)

Вычисляя интеграл по размаху крыла от выражения (9.7.16), C помощью разложения (9.3.7) получим

Х» =

4

- я 5 - 1 2 [ s ‘r

“ ’ ( - T S ' “) ] * f , „ . , ( - i s ) . (9.7.17)

 

 

п=0

 

 

 

Влияние изменения знака у параметра

<7= —

рассмот­

рено

в

монографии

Маклеклена

2.18

и 13.41, сле­

дуя результатам которого выражения (9.7.17) и (9.3.6,10) примут вид (см. Майлс 1“®^1)

Х* = 4 - ж ' ‘ 3

(9.7.18)

^ 2 n + l( - « )=

-

(9.7.19а)

 

Gekjn+i ^0. — 4^**^

=

( 2

4 ? f f ' ( | 5 = ) { ( 2 r + i ) +

(9.7.19b)

Коэф4)ициенты Л м о ж н о разложить в степеппой ряд

(^no аргументу ~

, радиус сходимости

которого (см. Мак*

леклен

2.17)

будет не менее ^/g, а

опыт показывает,

что ряд сходится в радиусе, близком к 4. Искомые разложе­ ния, полученные с помощью соотношений (2.151), (2.21),

(3.10) и

(3.33)

цитированной

монографии,

имеют вид

 

 

 

 

 

 

 

(9.7.20)

.+,5-.L , (Д. 0+ (| - - • - ) ( ! = ) + ( ' Ц - ¥ -•) ( т

OVo [ (| .) ‘ |) •

 

 

 

 

 

 

 

(9.7.22)

^2ri+l (

^®) =

 

 

 

 

 

 

 

 

1 + (п»+-)-‘ ( 4 0

- ю

[ ( | г У

]

 

= ( 2 r t + i r

, + з ( „ . + „ ) - . ( | , ) + о [ ( i . y j j

( о > 1 ) ,

 

 

 

 

 

 

 

(9.7.23)

 

 

2

= -^s®,

L = y + ln(^-^s^

(9 .7 .2 4 а ,b)

и у — постоянная

Эйлера.

 

 

 

 

Предполагается, что каждая

из функций

имеет

по крайней мере один полюс, так что интеграл %* имеет их бесчисленное множество. Но приведенные выше разложе­ ния (9.7.20—23) служат лишь для определения с приемлемой точностью наименьшего из этих полюсов, а именно

G e k ; (

о

.

.=0,

 

Ч

 

l3=rSo, во

 

5„= ^ 0 ,6 8 2 +

11,114=

1,306е'2>120^

(9.7.25)

где черточкой сверху обозначена комплексная сопряженная величина. Выделяя из выражения для %*(5) сингулярные части (особенности), связанные с этими полюсами, и, пре­ небрегая особенностями более высокого порядка, мы можем разложить выражение (9.7.18) в ряд по возрастающим сте­ пеням параметра s. Этот ряд имеет вид

X*(s) = 4-«»“‘

3

( v ) ”] +

 

Wi=O

 

 

4

 

+ ^-1

V 2

J ^^

+( ^ - т ) ( 1 0 ' +(^'-Т^ +^)(т^> +О(^1п’4

*

 

(9.7.26)

где символом а _1 обозначен вычет, связанный с полюсом

^),

который равен

 

 

а_1 = 0,760ё-»о-72з_

(9.7.27)

Обратное преобразование выражения (9.7.26) дает асимпто­

тические

разложения

 

 

 

X (T) ~ 4

л {

1 +

2Re t0 .,s; v » 4 + 4

т'^ -

 

 

-

I

[|п (4т)— g - ] х-- +

0(т-«1п «т)}

(9.7.28а)

И

 

 

 

 

 

 

Xt(T)

4 я

{ 2 Re

- 4 т ' “+

 

 

 

+ 3 [ In (4т) -

.g. ] т-“ + О (Т-’ 1п> т)}

(9.7.28Ь)

Хотя члены порядка

при т > 1 ,

входящие в выражение

(9.7.26),

при г > о не влияют на асимптотические разложе­

ния (9.7.28), следует отметить, что они могут влиять на эти разложения, если функция v* содержит члены типа 5 '“,

и т. д. (соответствующие членам тип? t, t- и т. д. в выра­ жении для функции о). Отметим также, что хотя, в принципе,

1) При вычислении а.] учитывались член

МЫ никак не обосновали здесь сохранение экспонен­ циальных членов и в то же время пренебрежение членами порядка t " “ ln*t, но расчеты показывают, что экспоненциаль­ ные члены доминируют в асимптотических приближениях на большей части диапазона, где эти приближения приме­ няются.

Результаты численных расчетов по формулам (9.7.10, 11,28) приведены на рис. 9.6. Мы видим, что приближение

Рис. 9.6. Интеграл от потенциала по размаху для единично­ го скачка скоса потока в момент времени т = 0.

(9.7.28) в пределах точности графического изображения

совпадает с более точным результатом

,+ х "’

во всем

диапазоне 2 < t < 4 .

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

 

необходимо

вычис­

лять лишь функцию

при малых

значениях х

(т < 2).

В работе Майлса

приближение Кирхгофа, соответ­

ствующее преобразованию Лапласа

от выражения (9.5.3),

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

значений т, где более точный результат

проще прибли­

жения Кирхгофа.

 

§9.8. Крылья квазималого удлинения и метод исследования низкочастотных колебаний

при околозвуковых скоростях

Результат теории крыла малого удлинения (§ 9.2) может рассматриваться как первый член разложения решения граничной задачи о крыле в ряд по степеням удлинения (см, примечание на странице 156),

Адамс и Сирс (*], (Ч получили следующий член этого

разложения для

случая

стационарного обтекания ^).

В этом разделе мы дадим

обобщение

результатов

Адамса

и Сирса на случай нестационарного

обтекания,

предпо­

лагая, что ArM=O(I). Еслиже ЛМ > 1

и удлинение крыла

не очень мало, то

задача

будет уже

относиться

к слу­

чаям В[ или В'1 из таблицы

1.

 

 

Наши рассуждения здесь будут весьма сходны с рассуж­

дениями § 9,4, но,

поскольку мы теперь сохраняем

произ­

водную по д: в дифференциальном уравнении в частных производных, вместо (|, Ti) следует воспользоваться коорди­ натами {у, г). Будем исходить из дифференциального урав­ нения (3.3.4), предполагая гармоническую зависимость от времени (3.10.2).

Тогда =)

Tvt, +

Ф .. = Р Ч х

+

2»/гМ“Ф , -

( 9 .8 . 1 )

Преобразование

Лапласа

от

этого

уравнения

имеет вид

 

ф5у4-Ф?г

=

 

(9.8.2)

где в отличие от выраокения (9.1.10)

параметр Я определен

соотношением

 

 

 

 

 

 

X2 = p V = 2 tW s -/ fe W ..

(9.8.3)

Уравнение (9.8.2) формально идентично уравнению (9.1.9), если в последнем Я предполагать мнимым, а переменные (|, г|) заменить на (у, г). Тогда, по аналогии C (9.4.1), решение уравнения (9.8.2), обращающееся в ф* (s, у)

при 2 =

о - f и в нуль на бесконечности, молено представить

в виде

 

 

 

 

 

 

 

Ф* (S, у ,

2) = ^

5 ^ '" ^ 1

 

 

1) Дальнейшее

развитие

этого направленн

в работах Лен-

дала

[**•)

и Эйхельбремнсра (®®].

 

уравнения

-)

За

исключением специально оговоренных случаев,

этого

раздела могут быть

записаны или через <р

и о,

или через

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

где теперь уже

и интегрирование по "п проводится в интервале (— оо, + оо). Последнее обстоятельство обусловлено тем, что потенциал Фо (s. "П) зависит, в общем случае, от фо (^, л) при всех зна­ чениях (Пределы интегрирования в выражении (9.8.4) могут быть заменены на ± бтах, однако пределы ± со не менее удобны.) Если выражение (9.8.4) продифференциро­

вать по 2, а обратное преобразование Лапласа

от

получен­

ного результата приравнять — v при Z = О,

(у)

< Ь{х),

то получится интегральное уравнение, аналогичное урав­ нению (9.4.5).

Однако, как и в уравнении (9.4.7), целесообразно сначала выделить сингулярную часть ядра, и тогда, опуская

нулевой

индекс

при

значении

потенциала

ф на

крыле,

получим

 

 

 

 

 

 

+

J g*i{s,

y — r\)(p*{s, 7\)dr] = v,

\y\< b{x),

 

 

 

 

 

 

(9.8.6)

a L — оператор,

определенный

соотношением

(9.2.15)

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

 

1ф =

% (X, Ч) 4т1

 

(9.8.8)

 

. (У— П)

 

 

 

 

 

Интегральное уравнение (9.8.6) можно решать методом

итерации подобно (9.4.24) ^). Тогда первая

итерация

будет

1) Чтобы получить полное и систематическое

разложение реше­

ния по степеням

полуразмаха

6, можно, вообще

говоря, воспользо­

ваться методом,

описанным

соотношениями (4 .4 .28 — 31),

однако

тогда нужно будет включить наряду с членами

вида

и члены

Iog Ь.

 

 

 

 

где <p‘“’ — решение, полученное по теории крыла малого удлинения ILcp" = v\. Первая итерация (9.8.9) дает решение с точностью до членов порядка о\х-Ь^), если только функция g* аппроксимирована первым членом разло­ жения в ряд выражения (9.8.7), т. е. имеет вид (ср. 9.4.25)

g',(y)=

+

( т »'1!'| )] + 0 [x y in (» .|< /|)].

 

 

(9.8.10)

Подставляя это выражение в формулу (9.8.9) и ввделяя член C InX, после перемены порядка операций интегриро­

вания пот] и преобразования Лапласа получим ^)

-Ах)

Ь(Л)

 

 

X \

1])

(9.8.11)

- Ш

 

 

Пользуясь результатами § 9.2, можно показать, что инте­ грал от по хорде стремится к нулю, при х - ^ 0 , как Ь\х), а следовательно, по условию (9.1.1), как х®. Обратное преобразование первого члена (9.8.11), таким образом, будет иметь вид

Нх)

-Ь(*)

Ь(х)

= ^ I [ - Т + 1 п (у | у - 1 1 | )]ф "’Ч ^ Ч )‘'Л- (9.8-12)

х )_______

Заметим, что при получении решения в ряде случаев бывает целесообразно сохранить функцию в форме (9.8.11) и выполнить операции преобразования Лапласа после проведения необходимых интегрирований по хорде крыла с помощ,ью преобразования Лапласа.

где символом Л обозначен оператор

Рассмотрим обратное преобразование последнего члена выражения (9.8.11). Обозначим

Ь(.х)

 

 

X" (■«) = \

(^1 л)

(9.8.14)

-« *)

 

 

Принимая, что P"i = 0 ( l) и пользуясь соотношением (9.8.3), разложим InX в ряд

1пХ = 1п(М + у 1 п [

] +

 

+ 4 l n [ ^ ± i ^ M

± l T L ]

(9.8.15)

Теперь, применяя теорему Фальтунга и ссылаясь на монографию Кемпбелла и Фостера [^®] (893, 894,2), можно показать, что

i?-4?^My + ln(Ps)]^X<“)} =

 

 

= AiS-^ {

], 5д а о )| _

(9 8 .16а)

X

 

 

= Л [(1пР)х<">-| 5 In(X-I)Xhhdldl]

(9.8.16b)

Jg-I |хЧ п

j;x i» )} =

 

=

( 1 - е х р [ - « Ш ( М ± 1 ) - ‘ ( * - | ) 1 } - ~ (9.8.17)

откуда

-Н х)

^

^

 

X

 

X Ф<»'(X,Ч) йлIr S (*- 1)х<« <И+

О

+ J [1 -

V-6) COS (Л мр'" {X - 1))1

I .

 

 

(9.8.18)

Если в выражение (9.8.18) подставить разложение для функции ф<®) из (9.2.17) и выполнить интегрирование по ij, получим

оС О ..»+ 1 л { [| п (I f » ) - у

] ft*(х )< ’ (X)-

 

о

 

 

JC

 

 

+ J l l - e-i'‘M23-8(,:_6)cos (ШР*2 (^x - 1))]

^

m=l

 

 

X COS

arccos-^^J I

(9.8.19)

Решение интегрального уравнения (9.8.9) с правой частью заданной выражением (9.8.18) или (9.8.19), получается теперь непосредственно по аналогии с решением, данным в § 9.2, в виде

ф(‘)(д?, у)

(х) 2

(л) Sin [ п arccos

(9.8.20)

n=5sl

где

п

=

J *)(;:, ^COS G)sin (п0) SinGлЮ. (9.8.21)

 

о

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