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

книги / Оболочки и пластины

..pdf
Скачиваний:
21
Добавлен:
19.11.2023
Размер:
71.66 Mб
Скачать

 

 

 

 

Jkn

«1

 

Pi

*

 

(4,19,2)'

где Ckn (As = 1, 2,

 

 

 

 

 

 

n= 1, 2, ;..)— искомые коэффициенты.

 

Подставляя (4,19,2) в уравнение (4,14,15),

потребуем, чтобы полу­

чившаяся функция была ортогональна всем функциям

 

 

 

xVsm(при s — 1,2,

 

ш = 1,2, ...).

 

Если в уравнениях (4,14,15) и (4,19,2) перейти к новой

перемен­

ной | по формуле (4,17,1) и подставить (4,19,2) в (4,14,15),

то после

ряда простых преобразований получим

 

 

 

 

 

 

С2, [(а2 +

п\у + а* -

К(а2 4- п?)2] £

Ck sin kn% -

 

 

 

 

 

 

 

 

k=\

 

 

 

 

 

— Akzi (a1+

/г*)2

Ck cos

 

= 0,

 

(4,19,3)

где

 

 

 

 

k=\

 

 

 

 

 

 

 

 

 

 

 

 

BVRa

 

 

 

 

a = *-n,

 

1

Pi

A =

 

(4,19,4)

 

 

 

/

 

 

 

Ehl

 

 

В вариационной форме Галеркина уравнение (4,19,3) запишется в

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

 

 

 

 

 

 

 

 

 

£

£

6Cfts {[(fl2+ rti)4 +

 

 

+ nf)*] | sin

sin siigdl -

*=1

s= l

 

 

 

 

 

 

0

 

 

 

 

А/гл (a2 +

rt2)21 cos /гл£ sin

 

=

0.

(4,19,5)

 

 

 

 

 

о

 

 

 

 

 

 

Интегралы,

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

(4,19,5),

имеют значения:

 

 

 

«

 

 

 

[ — при s = k,

 

 

 

\ sin kn% sin sn^dl = I

2

 

 

 

 

 

 

о

 

 

 

l

0 при s=jbk,

 

 

 

 

 

 

2s

при k -|- s — нечетном,

 

 

sin sn%dl =

(s2 — fe2) it

 

 

 

 

 

 

 

 

t

 

 

 

 

0

при k 4- s — четном.

 

 

 

 

 

 

 

решений

системы

линейных

Для определения отличных от нуля

однородных алгебраических

уравнений

(4,19,5)

необходимо

и доста­

точно приравнять нулю определитель системы:

 

 

 

т<р - ъ

— — А

 

0

 

15

 

 

3

 

 

 

 

 

 

 

±3А

T

(L-

X)

 

5

 

 

0

 

 

 

0

 

™-А

 

 

 

 

- 2 - А

.. =

0. (4,19,6)

 

 

 

5

 

 

 

 

 

7

 

 

- * - А

0

 

 

» - А

- у ( М - Х ) ..

 

15

 

 

 

 

7

 

 

 

 

 

Каждому собственному значению уравнения (4,19,6) соответствуют отличное от нуля значение Chs системы (4,19,5) и приближенное реше­ ние (4,19,2) рассматриваемой краевой задачи. В уравнении (4,19,6) приведены следующие обозначения:

F = С; (а2 + n\f +

а*

L-=C\ (4а2 +

а?)2 +

2*а*

(а2 + п\?

(4а2 + «I)2

 

 

 

 

К =

cl (9а2 +

п?)2 +

3

 

 

(9а2 -ф- п2)2 ’

 

 

 

 

 

М =

С? (16а2 +

п\)2+

44а4

 

(4,19,7)

(16а2 +

п2)2

 

 

 

 

При расчетах по первому приближению уравнение (4,19,6) прини­

мает вид

 

 

 

 

 

 

(F X) = 0.

 

 

(4,19,8)

Отсюда следует, что все собственные значения рассматриваемой крае­ вой задачи положительны и вещественны, следовательно, независимо от скорости потока невозмущенное движение панели в*классе решений

(4,19,14) устойчиво, а критическая скорость флаттера

цилиндрической

панели равна бесконечности.

(4,19,6), огра­

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

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

 

I 2 — X(F + L ) + F L + -^ 2 = 0.

(4,19,9)

Решение уравнения (4,19,9) даст два корня:

 

 

(4, .9,10)

Из формулы (4,19,10) следует, что собственные значения при рас­ четах по второму приближению существенным образом зависят от ско­

рости потока и при значениях скорости потока

 

 

 

 

А = BVR2

^

3

( F - L )

 

 

 

 

 

Ehl

^

16

 

 

 

 

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

 

 

Re^ =

X1= i - ( F + L),

1тХ =

Х2= [ - | - Л 2- - £ = ^

р

(4,19,11)

Подставляя значения

и %2

 

в уравнение

параболы

устойчивости

(4,14,20)

и учитывая, что

BVR2

после

ряда простых преобра-

А = ——— ,

Ehl

зований получим формулу для определения критической скорости во втором приближении, т. е. скорости потока, при превышении которой появляются комплексные собственные значения, лежащие за пределами параболы устойчивости:

Eh2

U

(4,19,12)

‘'■>тг-Ит[*(,'+ч+ ERR2 Р ( Г - Ц ] }

 

 

 

При расчетах по третьему приближению из уравнения (4,19,6) для собственных значений будем иметь

— Х» + /1Х*- /,* , + /з = 0.

(4,19,13)

Здесь

 

I ^ F + L + K,

 

I, = FL + FK + KL + A * [ ( ^ - y + ( ± J ] ,

(4,19,14)

Полагая X=Xi + iiX2 и отделяя в уравнении (4,19,13) действитель­ ную часть комплексных собственных значений от мнимой, получим си­ стему двух уравнений:

— X] + 3XjX2 -}- / х (X] — Х2) — / 2Хх —|—/3 = 0,

(4,19,15)

3Xl -{- Х2 “I- 2/хХх /2 = 0.

Для дальнейших выкладок удобнее представить уравнение пара­ болы устойчивости (4,14,20) в виде

Л >= ^ { г ~ ~ т ~ У

(4' l9•16,

Подставляя (4,19,16) в уравнения (4,19,15), получим систему

ЗХ?-Х1(2/1+ г) + /2= 0,

(4,19,17)

X? (/х + 8г) — Хх (2/, + Зг/Х) + 3/3= 0.

Для определения общего корня двух многочленов необходимо при­ равнять нулю результат этих уравнений:

3

— (2/ 1т

г)

h

0

 

 

0

 

3

 

- < 2/i f г)

^2

 

 

/ х +

8г - ( 2/, + ЗгЛ)

з/,

0

 

 

0

/х +

 

— (2/2+ Зг/Х)

31,

 

 

или

 

 

 

 

 

 

 

 

Wi + 8г) {21х+

г) -

3 (2/, + Зг/Х)] [3/3(2/х +

г) -

 

-

/2 (2/2 +

Зг/0] + [/, (/х + 8г) - 9/3]2 -

0.

(4,19,19)

Уравнение

(4,19,19)

позволяет

исследовать

характер

изменения

комплексных собственных

значений

рассматриваемой' краевой задачи

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

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

Здесь

I[ = F + L + K + M ,

I2 = MF + KT + FL + MK + ML + KL +

l'z = FKM + MLF+ MKL + KLF + 4А* { [ ( - f - ) 2 + ( - у - ) * ] F +

+[(■£)’+(i)-]-+К™У

 

(4,19,21>

’]

+(тТ]*}•

/, = PMKL +

4Л* [ ( - f ) ’ MF +

( ^ - ) * и . + ( | - у км + ( J L y KLj

 

+

16Л4 tf-x-f № )(!)]■

Подставляя в

уравнение (4,19,20)

значение Я=Я1-ИЯг и отделяя

действительную часть от мнимой, получим

 

Я) — / 1Я1

6Я1Я2 -(- 3/1 ЛГХЛ,2 “Ь Яг “Ь ^2^'1 I 2^>2 — /з ^1 “Ь 1а О,

 

 

 

 

 

 

(4,19,22)

4Я?— 4Л.1Л| — 3/|Я? + l'{kl + 2 /^! — /3 = о.

Заменяя в уравнениях

(4,19,22)

Яг через г%\,

согласно (4,19,16) полу­

чим систему уравнений:

 

 

 

 

Я? (/', + 20г) - Я? (11 r l\

+ 4г2 +

2/ 2) + К

(4г/а +

3/ 3) - М \ = 0,

 

 

 

 

 

 

(4,19,23)

4Я? — Я] (3 /, -)- 4г) -J- Я^ {I\f -f- 2/ 2) — / 3 ="• 0.

Составляя результант системы (4,19,23)

и раскрывая его, получим урав­

нение, позволяющее определить критическую скорость потока в четвер­ том приближении:

(ах0 — /id)3 + {axk — q/i)2 (^О — kd) + 2 (^d — ах0) (cfi —

 

kd) (af — 6/j) + (/d — 60) (каг — c-J^) (ax0 — lxd)

 

_ (fd — bO)2(6/i — a j) + (cfi — kd) (bk — ctf) (bl±aj) = 0.

(4,19,24)

Здесь введены следующие обозначения:

 

 

 

 

аг = l\ + 20r, b =

4г2 + 11г/; +

2/2,

 

сi = 4г/2 -Ь З/3 ,

d =

4/4,

/i =

4,

(4,19,25)

/ = 3/; + 4г, k =

rl\ +

2/2,

0 -

/3.

 

Формулы (4,19,12), (4,19,19), (4,19,24) позволяют вычислить по1 методу Галеркина (во втором, третьем и четвертом приближениях) критические скорости потоков для цилиндрических панелей, шарнирно опертых по всем краям и движущихся в газе со сверхзвуковой скоро­ стью, если известны геометрические размеры панели и постоянные, ха­ рактеризующие газовую среду.

§ 20. ПОЛОГАЯ СФЕРИЧЕСКАЯ ОБОЛОЧКА

Исследуем автоколебания пологой сферической оболочки в случае шарнирного опирания по всем кромкам и в случае жесткого защемлеиия по всему контуру

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

мает вид

 

 

 

 

 

 

 

 

 

 

 

 

 

У4Ф +

1—V2 Ф

ftp со2Ф — BV

дф

+ -|-(йФ = 0.

(4,20,1)

 

 

h.2R2

 

D

 

 

D

~дх

 

 

 

 

 

Вводя

обозначение

 

со2

со = — л(1),

получим

уравнение па-

раболы устойчивости:

 

 

 

 

 

 

 

 

 

 

 

 

 

hpD

 

 

 

 

 

-f-fl'i).

 

 

(4,20,2)

 

 

 

в\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В системе новых безразмерных

координат

х = а% и y=br\

уравнение

(4,20,1)

можно записать в форме

 

 

 

 

 

 

 

 

 

ft4

д*Ф

9 Ji?

А*

д*Ф

.

h*

д*Ф .

12(1- у 2)h 2

ф

 

а*

д|4 +

а2

Ъ2

~д12дц2

+

“ ------Г Т " "Г

R2

 

 

 

 

Ь4

<3т]4

 

 

 

 

 

 

 

BVh*

дФ — Ш>= 0,

 

 

 

 

(4,20,3)

где новое

 

 

aD

5?

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bth*

 

 

 

 

 

Х = Х % \

X ^ - ^ - q 2, Х2

Я-

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Условия шарнирного опирания оболочки по краям будут тождест­

венно удовлетворены,

если искать решение уравнения

(4,20,3)

в виде

 

 

 

 

Ф = sin kn%sin шщ.

 

 

 

 

(4,20,4)

Подставив (4,20,4) в уравнение (4,20,3), получим

 

 

 

 

 

[* ( - f * * + £ • "О ~

* ] 81,1 ^

+ - "‘‘- g *

*

х

 

 

 

 

х

sin knl

 

kn cos kn\ = 0.

 

 

 

(4,20,5)

Применяя к уравнению (4,20,5) метод Бубнова—Галеркина для опре­ деления собственных значений и отличных от нуля решений системы алгебраических уравнений, получим определитель

— (F — X )

—А

0

 

2

3

 

 

т *

 

~ Г

*

0

- Т А

T < L -

X>

1 Решение задачи принадлежит Р. Д. Степанову.

где

R 2

/l2_ i = '4‘ [ 9 - 5 - + ^ - ] , + 12(1- v2> R2

Не воспроизводя здесь выкладок, аналогичных тем, которые были сде­ ланы для цилиндрической панели, получим формулу для определения критических скоростей потока во втором приближении по методу Буб­ нова—Галеркина:

--------\2(F + K) +

------^ -----

(F—/С)2] - (4,20,6)

'

12(1 — V*) L

12В2 (1 — V») v

' ] 4

Критическая скорость в третьем приближении определяется из уравне­ ния

[/. ( / i

+ 8г)- 9 / 3] 2+ [{I, + 8г) (2/, +

г) -

 

- 3 (2/, +

3г/,)] [3/, {21, + г) — / 2 (2/, +

3r lj\ = 0,

(4,20,7)

где

I1 = K + F + L,

I2= F K + KL + F L + A * \ ( l f f + (j - y \ ,

л - т а + ^ [ ? ( ^ - ) , + 4 ( 4 . ) * ] .

а г — параметр параболы устойчивости: rXi = xl-

Отметим, что из решения, полученного для сферической панели, нетрудно получить решение для пластины, шарнирно опертой по всем кромкам; для этого достаточно в (4,20,6) положить R->-oо.

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

Здесь, как и ранее, применяется метод Бубнова—Галеркина, при­ чем в качестве аппроксимирующих функций используются фундамен­ тальные балочные функции. Известно, что балочные функции, ортого­ нальность которых хорошо изучена, не сохраняют это свойство в отно­ шении своих производных первого, второго и третьего порядков, и пото­ му некоторые авторы вводят понятие о квазиортогональности этих функ­ ций, т. е. считают интегралы от произведения вторых производных балочных функций на ту же функцию величиной пренебрегаемой ■. От­ метим необходимость проявлять осторожность при постулировании свойства квазиортогональности фундаментальных балочных функций.

Для решения задачи будем исходить из уравнения (4,20,3), в кото­ ром A,=Ai + tX2— комплексные собственные значения. Координаты пара­ болы устойчивости в этом случае будут

К

h

B,h*

D Я-

1 Ряд новых квадратур вычислен Р. Д. Степановым [48].

Для удобства запишГ введем обозначения:

12 (1 — у2) h2 _

1 — у»

к1 д

_ ВУк*

(4,20,8)

R? ~

с2

' R* '

~ ~ а 6 ~

 

Разрешающую функцию Ф (|, г))

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

балочных функций

 

 

 

(4,20,9)

ф(£,л) = х т ® к я(л)

 

каждая из которых удовлетворяет как дифференциальным уравнениям

X™(I) = %Um( I ) ,

YV(л) = t i v n(л).

 

 

так и условиям защемления сферической панели на краях

 

 

 

5 = 0,

5 = 1

и л = 0, Т)=1.

 

 

Подставив (4,20,9) в (4,20,3), получим

 

 

 

~

С

а д , +

2 ^ - х ту; + ^ Г х тг х

+

 

а4

 

 

 

а 2Ь2

 

Ь4

 

 

+

Г 1 —у2

к4

- А,] ХтУ„ -

=

0.

(4,20,10)

L

2 с

4 Я

 

 

 

 

Если умножить все члены (4,20,10) на

Уп и проинтегрировать по л

в пределах от 0 до 1, получим

 

 

 

 

 

J ^ K U m+ 2 ^

X

m’ + ^ X t X m +

 

^ j

АХ'т= О,

где

 

 

 

 

 

 

 

(4,20,11)

 

 

 

 

 

 

 

 

 

 

 

 

р =

J'У'пУ«*\

 

 

(4,20,12)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если теперь умножить

(4,20,11) на Xs и проинтегрировать по | от

0 до 1, то для определения отличных от нуля решений алгебраических уравнений

 

 

 

ХтХ, +

/71=1

 

 

 

+ 2

$ ^ - X mX s~ A X mXs\ = 0

(4,20,13)

 

a 2b2

J

 

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

[F — X,] un +

2(i

Vu

Awn

Лш12 2Р ^ 2 1^хз

 

-Ае>21

[/С -Ь ]и » +

2 Р -^ -У и - Л ю 1г -Л ш 23

0. (4,20,14)

— Лш32 [L — ^Иад + гр -^ -У зв —^ 'зз

Здесь

F =

я.1-

h\

Q4

 

 

к = Я-2

/I4- +

 

 

a4

L =

 

/14 +

 

 

a4

+

1

fc4

C2

x t 4

1 - *

64

c2

.

1— V2

64

c2

/14

tf4

/14

(4,20,15)

tf4

/14

Д4

ии = | Х Д Д ,

Vij = J X]X,4l,

wi} = j Х',Х&.

(4,20,16)

0

0

о

 

В исследуемом случае эти интегралы имеют следующее значение:

ип =

1,0359,

ш12 = — 3,399,

и22 = 0,9984,

 

ш23 = — 5,512,

vu = — 12,775,

ш13 = 9,9065,

2Q j -

о22 =

—45,977,

w31= — 9,9065,

ши = 0,

 

ш21 = 3,399,

ау22 = 0,

ш23 =

5,512.

 

 

Каждому собственному значению- К уравнения

(4,20,14)

соответ­

ствует отличное от нуля значение Ст системы

(4,20,13) и приближен­

ное решение рассматриваемой краевой задачи.

 

(4,20,14)

прини­

ч При расчетах во втором приближении

уравнение

мает вид

 

 

 

 

 

 

 

 

[F — к]иг1 -1-20-^- vu — Awn

Awn

 

 

 

 

 

 

 

 

 

- 0 .

(4,20,18)

Awn

[К — Ц и 22 +2P

 

v22 — AW22

 

 

 

 

 

a 2b2

 

 

 

 

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

квадратное уравнение

 

 

 

 

 

 

 

где

 

*,2- / ^ + / 2 =

0,

 

 

(4,20,19)

 

 

 

 

 

 

 

 

 

Г

vll

v22 I

__Г Ц*11

/ 1 = F + K 2р-£гг \

«п

+ "!!!

J

 

a2b2

 

U 22

L

 

L «11

 

/>= F/C + 2P -^ - r/C -^ + F - ^ l — 2Pi4 /l4

a262 L

«11

«22 J

a2b2 un a 22

^22 1

«22 J

[wlvv22 + ontt/22] —

- Л Г*

+ F ^ - 1

+ A * - L - [wnw22-

Ш21ау12] + 4p2

L ttU

u22 J

ullu22

a Ь Ui]U22

 

 

 

(4,20,20)

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

Я,3 — %Ч[ + Я,/2 — /з = 0,

(4,20,21)

где

/; = К + F + L + 2Р-^-Г-2и- +

 

^33

«22

 

ааЬа L Иц

« 3 3

"I _

^ Г Wll

^ 2 2

^ 3 3

J

L MU

« 2 2

«33

1 .

J

I2 = FL + KL + F K +

_h*_

[(K + L)

Уп

+ (^ + ^)

^22

+

а2А2

«и

« 2 2

+ (F + /o-^-l-^f(/c + L)-^ + (F + D ^

+ (F + F:)-^1 +

«33 J

 

 

 

«33

J

 

 

 

L

 

 

 

 

 

«Ц

 

 

 

 

«22

 

 

 

 

 

 

 

 

/I8

 

Г

^11«22

 

 

22^33

 

^11^33

^13U31

 

 

 

 

 

 

+ 4P

a464

 

 

 

 

+

-

«

 

 

+

 

 

 

 

 

 

L

«11«22

22«33

«11«33

«11«33

]

 

 

 

 

 

 

 

^11^22

 

 

 

 

 

 

 

«

 

 

 

^п^зз

 

 

 

— 20Л

aW

 

 

|

^11^22 ■+

 

^11^33 ■+

^33^22 +

^22^33

 

 

+

 

 

 

L

«11«22

 

 

 

 

«11«22

 

«11«33

 

«33«22

«22«33

 

«11«33]

 

 

 

 

[

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_|_

^2

I

Ш11Ш22

 

 

^11^33

_J_

 

t^22^33

 

^21^12

1^23^32

 

 

 

 

 

 

 

 

[■

« 11«22

 

 

«11«33

 

 

 

 

«22«33

 

«11«22

«22«33 ]

 

 

 

 

 

 

 

/: = F 7 a — л Г/сь-^- +

F L - ^ _ + F:F - ^ - 1

+

 

 

 

 

 

 

 

% %

 

 

 

 

 

L

 

«11

 

 

 

«22

«33

J

 

 

 

 

 

-f Л»

I

 

+

A:

tdulilss

 

 

 

^22^33

_

^21^12 _

 

^23^32

+

 

 

 

«11«22

 

Г

«11«33+

F

 

«22«33

 

«U«22

 

«22«33 ■']

 

 

 

 

 

+

 

2 р

А4

 

гг

t>ffl

 

K F - ^ -

+ / a

^11

 

 

 

 

 

 

 

 

 

 

- ^

 

ГFL

«22

 

 

+ «11

 

 

 

 

 

 

 

A8

^

a2A2

L

 

 

 

 

 

 

«33

 

 

 

 

 

+

4p*

 

 

 

 

 

 

 

UUU33

I

С* У22У33

 

^13^31

~j __

 

 

 

a4A4

 

 

u llu 22

+ K

 

«11«33

+ F

«22«33

 

«11«33

J

 

 

 

— 4pM

A8

 

'^11U22W33

 

«11^33^22

^33^11^22

^22^31^13

+

 

 

 

a464

|_

«nU22«33 +

«ц«22«зЗ +

 

«11«22«33

 

«11«22«33

 

 

 

+ 2рЛ2 A4

 

 

 

1

 

[ШцОУг^зз +

ииоу22ш33 + а>цШ33у22 +

o>2iou32t;13 —

 

 

 

a2A2

 

un u„2u33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— г>33ВУ21ш 12 +

^12^23^31 —

^

з ^

п

]

+

 

ЭД3 - ^ 7 - [Уи и 22У33— У31я 130221

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

а8А8

.....................................................................Uu«22«3S

 

 

 

 

 

 

 

 

 

[йУиш22ау33 — ш33ау210У12 — шиау23а)32] —

 

 

 

 

 

 

— Л3 «И«22«331

 

 

 

 

— 2рЛ —

L

4^11^22

_|_

р

 

1^11^22

|

р

^22^33

| р

^22^33

+

 

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

 

«22«33

 

«22«33

 

 

 

 

к

 

а2А2 [

 

« 11«22

 

 

 

« 11«22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

я

Юзз«п

 

+

я

U33^11 I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

«11«33

 

 

 

 

«11«33

 

 

 

 

 

 

 

Решение уравнения

(4,20,19) дает два корня

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 ,2 -= А

+ Г -

 

 

 

 

 

(4,20,23)

 

 

 

 

 

 

 

 

 

 

 

 

2

_

 

[

4

 

 

 

 

 

 

 

 

 

Из (4,20,23) можно определить, при'каком значении скорости собствен-

 

ные значения

исследуемой

краевой

задачи становятся

комплексными

 

и движение оболочки в потоке газа становится неустойчивым:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re X — % 1 =

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ImX =

^2 =

[ / 2- - j ] ,/*

 

 

 

 

 

 

 

Подставляя A,i и %2 в уравнение параболы устойчивости

*1 =

B\hрР3 ч

и заменяя А его значением

 

А =

BW

 

a D

после ряда простых преобразований получим формулу для определе­ ния критической скорости потока во втором приближении.

Критическую скорость потока во втором приближении можно по­

лучить и другим способом.

действитель­

Подставив в уравнение (4,20,19) h= X\ + iX2 и отделив

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

 

М — ^2 — ^i^i ~т~ ^2 —

(4,20,24)

--/i^»2 ~ 0.

 

Уравнение параболы устойчивости удобнее представить в форме

Ц =

 

где

(4,20,25)

_ 12В2 (1 — у2)

 

РЕ

 

Подставив (4,20,25) в (4,20,24), получим следующую систему из двух уравнений:

 

 

Ц - г Х 1— 11,к1+

/3 =

0 ,

(4,20,26)

 

 

 

2Хх — = 0.

 

 

 

 

 

 

Для определения общего корня двух многочленов необходимо со­

ставить и приравнять нулю результант этих уравнений:

 

1

~ ( г + 1г)

К

0.

 

 

 

0

— 1

~ ( г + Л)

/*

1

 

 

0

2

~ h

0

2

~ I i

o

0

0

2

- h

0

2

- / ,

После ряда несложных преобразований получим уравнение, из ко­

торого можно определить А:

 

 

 

 

 

 

4/s — 2 / / — /? = 0.

 

(4,20,27)

Скорость V легко теперь определить по формуле

 

вл4