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

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

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

Г Л А В А 11

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

§11.1. Введение

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

эффективное удлинение которого (Л = = 2^Ь, см. (7.5.15)) меньше единицы ^), и тем самым дополним резуль­ таты главы 7, где рассматривались крылья с эс|хрективным удлинением, превышающим единицу. Как указывалось в §7 .5, две боковые кромки крыла всегда взаийодействуют при Л < 1. Хотя взаимодействие последовательных отра­ жений от кромок можно рассчитать по методу Шварцшильда (см. § 9.7), более удобным является метод, основан­ ный на введении эллиптических координат, за исключением начальных стадий переходных движений (ср. § 9.7 и § 11.4). Определенными достоинствами обладает и метод сведения к интегральному уравнению (§ 9.4), использованный Стыо-

артсоном

при решении стационарной

задачи.

 

Формулируя

рассматриваемую задачу,

удобно

ввести

независимые переменные

т], которыми мы пользовались

в теории крыла малого удлинения (см. (9.1.3а, Ь)), поместив начало коор; инат в середине передней кромки крыла

(а не в угловой точке). Тогда соотношения (7.1.1—3) при­ мут вид

 

=

( 11. 1. 1)

_____________ Ф;|т,=о= - 6 У *

( ! I K 1)

(11.1.2)

9 В принципе, результаты § 11,2 справедливы для крыльев про­ извольного удлинения, но их практическое применение ограничено значениями А. Пример из § 11.3 показывает, что выражение

(11.2.7) справедливо при А < I, а теория крыла малого удлинения— при А < 1 / 2 .

где \ вновь определено выражением (3.11.6),

а не (9.1.10).

Граничная задача,

формулируемая соотношениями

(11.1.1—3), формально

отличается от задачи

(9.7.12— 14)

лишь более общим видом граничного условия на крыле. Но имеется и весьма важное отличие, состоящее в том, что система (11.1.1—3) является результатом преобразования Фурье по модифицированному времени i' и преобразова­ ния Лапласа по координате х', направленной по потоку, тогда как система (9.7.12— 14) получена преобразованием Лапласа по истинному времени г, а производными потен­ циала скоростей по координате JC пренебрегается. То обстоя­ тельство, что формально возможно исследовать задачу о прямоугольном крыле малого удлинения в рамках теории крыла очень малого удлинения, сохраняя, однако, про­ изводные по координате, направленной вдоль потока, является следствием постоянства величины Ь.

C другой стороны, следует подчеркнуть, что функция Ь, характеризующая размах прямоугольного крыла, имеет разрыв при А' = о и тем самым оказывается нарушенным условие (9.1.1), налагаемое в теории крыла малого удли­ нения.

§ 11.2. Решение задачи в эллиптических координатах Установив аналогию между задачами (11.1.1—3) и

и(9.7.12— 14), искомое решение у поверхности крыла можно получить модификацией выражений (9.7.15,16) или (9.3.2,3),

ионо будет иметь вид (см. Майлс 1^“®1)

ф* = 6 S

( — *^^)se„ (V, — q),

(11.2.1)

где теперь

 

 

 

= ^ 5 V^*(6cosv)se,^(v,

— ^)sinvdv,

(11.2.2)

P

_ Oefen (0.

—ч)

2 з\

 

Gefe; (О,

—<7)

(Jl.^ .d )

Если, как в § 9,7, функция V* не зависит от у, вы ние (11.2.2) упростится и примет вид

{ — д)

(/г— нечетное),

(11.2.5)

а преобразование Лапласа

от

интеграла по размаху (здесь

этот интеграл отличается

множителем Ь~ от аналогичного

выражения из главы 9) имеет вид (ср. с (9.7.18))

 

X*= \ 0*dy-- i

S

[ ( ( / ) ] = =

я , . . , ( -

1Щ .

- Ь

Ti=I

 

 

( 11.2.6)

 

 

 

 

Следуя (9.7.20—26), выражения (11.2.1)

и (11.2.6) могут

быть разложены в ряд по степеням и логарифмам Xb. Однако здесь нам нужно произвести обратное преобразование по s,

тогда как и SQ, согласно (9.7.25), представляют собой зна­ чения Xb, при которых функция F имеет полюса. Искомое разложение функции х* (см. (9.7.26)) имеет вид

X* = Y ^ bW * ja.iSor-» [ (6s — г„)-1 +

TO=O

+ 5 _ iV o ' [(^ S -^ o )'M -^ ;' S

+

+ 0{X^^bЧn^Xb)^ ,

(11.2.7)

где теперь

 

( 11. 2. 8)

a величины SQ ” заданы равенствами (9.7.25) и (9.7.27) соответственно ^).

§ 11.3. Производные устойчивости

Выражение (11.2.7)' можно использовать для расчета производных продольной устойчивости. Мы приведем здесь расчеты производной коэффициента подъемной силы по углу атаки (Майлс, [^“Ч) и границы области устойчивости при движении C одной степенью свободы ^).

Koa(Jxj)Hциепт подъемной силы прямоугольного крыла

вфизических координатах х, у имеет вид

1ь

+(И.3.1)

О- ь

где перепад давлений на крыле подставлен из соотношения (3.3.5), а производная потенциала по времени (р, сохранена для целей § 11.4. Поскольку потенциал ф на кромках обра­ щается в нуль, последовательность операций дис})ференцирования по X или ^ и интегрирования по у может быть изме­ нена, и выражение (11.3.1) примет вид

1

ь

о

1

 

 

= T

+

(П.3.2Ь)

1) Значение a_i нз

(9.7.27) является

более точным, чем данное

в работе Майлса [®®б], где членами порядка О (s®) в расчетах пренебрега.аось.

®) Выражения для остальных производных устойчивости были получены Майлсом [®®®]. Численные результаты см. в работе Вино­ града и Майлса [®®*], а для случая очень малых удлинений н пысо- Kiix частот— Мазельского [*®®J.

Обращаясь теперь к случаю установившегося обтек

ния (A = H = 0), из выражения (И .3.2Ь) получим

 

 

 

2 .

(11

.3.3)

' Ь

\х=1<

6Ь ^ Ii|.х'=р-1’

 

 

Вычислим значение Cr, (11.3.3) для плоского

крыла

под

углом атаки а, для которого

 

 

 

U= а,

y * = as“^.

(11

.3.4)

Подставляя эти значения в выражение (11.2.7)

и полагая

H = O (т. е. Го = SQ и Я, =

s), после подстановки полученного

выражения в формулу (11.3.3) получим

 

 

={ o j [ ( 6 s - s „ ) '4 -

+S ( т Т ] + “-[№ »-'«)“ +

TTi=о

m=0

где, в функции L, Я, = s. Поскольку члены вида s™ при m > 0 не влияют на окончательный результат (а члены вида s™ In s оказывают влияние), выражение (11.3.5) можно переписать в виде

+ 1 + а т » ' п ( 4 ) + ( А ) % »

Производя обратное преобразование от выражения (11.3.6), полагая х' = и исключая Ь с помощью соотношения (7.5.15), после дифференцирования по а получим

- l , 1 6 0 e x p ( = ^ ^ ) c o s [ ( ^ 2 ^ ^ - 1 - 0 , 2 9 9 ] }

( 1 1 . 3 . 7 )

График этой величины, поделенной на значение для крыла бесконечного размаха Сх,о = 4Р"\ представлен на рис. 11.1,

Рис. 11.1. Производная коэффициента подъемной силы, вычисленная по формуле (11.3.7); для сравнения приведены значения, полученные по теории крыла малого удлинения.

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

нения

Точное (в линейной постановке) выражение для производ­ ной коэффициента подъемной силы по углу атаки при Л >1/2 было получено Лагерстромом [“ ®] в виде

CL« =

('’ > ! )

(11.3.8а))*

*) См. работу Р.

Джонса [^б]. Влияние нелинейности (по а ) ,

связанной со сходом

вихрей, и^ледовалось Ченом

^ который

показал его важность при малых X.

CLO = J

+

 

+ 1 ц + Л - Ч > ^ Т ^ “- ^ 1 п Ш ' " ^ }

1< 1.

 

 

(11.3.8b)

При всех значениях Л < 1 приближенное выражение (11.3.7) дает численные результаты, совпадающие с точными

Рис. 11.2. Границы устойчивости прямоугольного крыла относительно низкочастотных продольных колебаний; кривые для А > 1 перестроены с рис. 7 .7.

(11.3.8) ДО третьей значащей цифры. Аналогичное сравнение значений производной коэффициента продольного момента Сма показало расхождение на 3% при A = I.

Границы области устойчивости для случая движений C одной степенью свободы, представляющие собой, по-види­ мому, наиболее строгую проверку точности приближенного выражения (11.2.7), приведены на рис. 7.8 и рис. 11.2 (Виноград и Майлс 1®°®]). Сравнение с соответствующими результатами § 7.6 при А = 1 (см. точки, обозначенные кружками, на рис. 11.2) указывает на прекрасное совпаде­ ние обоих решений. Из этого можно заключить, что резуль­ таты, полученные в §11.2, дают необходимое развитие тео­ рии неустановившегося обтекания на случай прямоуголь­

ных крыльев C

1 и K = O (1), тогда как при А < 1/2

во многих случаях

применима теория крыла малого удли­

нения. Более того, если к велико, а A = О (1), как при око­ лозвуковом обтекании, получающиеся результаты после преобразования к физическим координатам (Майлс [-“®]) очень сходны с (10.6.19—21) и, по-видимому, обеспечивают

такую же точность (см. § 10.6). Если, с другой

сторонЬт,

А >

1, как это имеет место в случае переходных движений

(см.

§ 11.4), и эффективное удлинение мало <

1), метод

§9.7 обеспечивает необходимую степень приближения. Расчетные данные, приведенные на рис. 11.2, указы­

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

удлинения неустойчивость не может иметь

места при

А < 0,65, а следовательно, вообще при M = 1.

Этот вывод

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

рассматривая треугольное крыло (§ 10.6). В этой связи

отметим, что результаты § 9.8, использованные

в

§ 10.6,

не применимы в задаче о прямоугольном крыле.

 

 

§ 11.4. Нагрузки при порыве ветра

 

 

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

случай входа

крыла

в единичный ступенчатый порыв,

для которого

функция

V имеет вид

 

 

 

 

 

 

(11.4.1)

где / и T безразмерные времена из § 3.3 и § 9.1 соответствен­ но. При малых значениях аргумента (т — х/МЬ) производ­ ные по времени в дифференциальном уравнении велики по

сравнению с производными по координате в направлении потока. Поэтому последними можно пренебречь, так что оказывается применимой теория крыла малого удлинения (§ 9.7). При больших значениях аргумента (т — х/М6) можно воспользоваться обратными преобразованиями соот­ ношений, полученных в § 11.2, однако сами эти обратные преобразования весьма сложны. Для практических целей, видимо, достаточно воспользоваться теорией крыла малого

удлинения

при

Л

< 1/2 (см. § 11.3) и результатами § 7 .8

при i4 >

1, а

в

промежуточной области воспользоваться

интерполяцией.

 

 

Выражение для коэффициента подъемной силы, соответ­ ствующее закону изменения скоса потока па крыле (11.4.1), получится, если в выражение (11.3.2) вместо % подставить

величину Ь^х(х — х/ЬЛЬ),

полученную в § 9 .7 по

теории

крыла малого удлинения, где множитель Ь-

связан здесь

C иным определением функции %. После подстановки получим

CL = I [ Ь^Х ( — i i r )

+ М-‘ S 4 ^ ^ Ь--х ( ^

- 4 )

d x ]

 

 

 

(11.4.2)

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

 

 

CL = 2 6х (T) = 26х ( 4 )

 

(11 .4.3)

Если в выражении (11.4.3) выделить асимптотическое (по теории крыла малого удлинения) значение производной

коэффициента подъемной силы по углу атаки (^nb = у itX j , получим

Cbit)

Ix(4)

 

2

(11.4.4)

 

 

Правая часть этого выражения представлена графи­ чески на рис. 9.6. Полезно сравнить этот результат с соот­ ветствующими данными для треугольного крыла малого удлинения (рис. 10.7,6).

Сравнивая результаты (11.4.4) и (7.8.8), заметим, что

Принципиальные различия между двумя этими резуль­ татами состоят в следующем: (а) коэффициент, вычислен­ ный по формуле (7.8.8), стремится к своему конечному зна­ чению монотонно, тогда как величина (11.4.4) приближается

к

своему предельному значению, колеблясь около него;

(Ь)

величина

(7.8.8) достигает своего конечного значения

за конечный

промежуток времени ^ = (М — 1)“^, а величи­

на

(11.4.4) за

бесконечно большое время (хотя при т > 4

отличие от предельного значения очень мало). Первое из отличий вызвано интерференцией боковых кромок при Л < 1, а второе является просто следствием использования приближенной теории крыла малого удлинения, поскольку (см. рис. 5.5 и относящиеся к нему рассуждения) предель­ ное значение коэффициента подъемной силы при всех зна­ чениях А должно достигаться к моменту времени ^= (М — 1)"^.

Однако для крыльев с эффективным

удлинением Л < 1/2

величина т согласно (11.4.5) при (М — 1) t = 1 всегда боль­

ше 4 (если Л = 1/2 и M = V 2, г = 9 ,7

при (М — 1) / = 1),

и потому погрешности очень малы. В

заключение отметим,

что если выполнены оба неравенства / < (М + 1)'^ и / < 26

(т. е. рассматривается начальная стадия входа в порыв), выражение (11.4.4) тождественно выражению (7.8.8а).

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