книги / Потенциальная теория неустановившихся сверхзвуковых течений
..pdfГ Л А В А 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а).