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

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

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

 

i_

JL

(

1

а^

)

1

дБ

dw

 

В

' а р

в

" ' “ а р

А2В

да

да

 

 

 

 

 

 

 

 

 

 

 

 

(4,11,7)

Теперь легко заметить, что члены в уравнениях

(4,11,4), содержа­

щие АХ, AY, AZ, линейны по отношению

к внутренним

силам

Т\, Т%,

а также по отношению к перемещениям и,

v, w и

их

производным.

В случае периодической внешней нагрузки силы Ти

 

также

перио­

дические функции времени; система

(4,11,4)

имеет в этом случае перио­

дические коэффициенты. Полагая

 

 

 

 

 

 

 

 

 

и (а, Р, t) =

V u A(/)cpft(a, р),

 

 

 

 

о (а, р, 0 =

2 о * ( 0 ^ ( « . Р).

 

 

 

 

w(a, Р, t) =

^

wk(t)%k(a, р)

 

 

 

(4,11,8)

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

обыкновенных

дифференциальных, уравне­

. ______ {_________

ний с

периодическими

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

Методы решения таких

систем достаточно

’7 =нр '

хорошо разработаны.

 

 

 

 

§

12. ЦИЛИНДРИЧЕСКАЯ

ОБОЛОЧКА

 

 

Пусть

круговая цилиндрическая

обо­

 

лочка

с радиусом срединной поверхности R

 

и толщиной h загружена, равномерно рас-

ис’

пределенной

радиальной

нагрузкой

q0+

+ qtcosQt и

сжата продольной силой

А> +

 

+ PtcosQt. Систему координат примем согласно рис. (4.1), введя безраз­

мерную продольную координатуа = — . Перемещения обозначим: через

R

и — в направлении образующей, через v — окружное и через w — ра­ диальное.

В уравнениях (4,11,4) в границах принятых допущений, имея в виду цилиндрическую оболочку, следует положить

_

d 2

 

[1 — у

а 2 __ _д* ,

1 — у

а 2

11 _

За*

2

'

ар2

22 ~ ар2 +

2

' да2

 

_

1

-j-

у

а 2

 

(4,12,1)

 

 

 

 

 

 

 

да

Здесь

26 П. М. Огибалов, М. А. Колтунов

Заметив, что внутренние силы, отвечающие исходному безмоментному состоянию, выражаются в виде

Тг

1

(P0 + Ptcos 01),

 

R

 

(4,12,3).

Т2 = R(q0 + <7, cos 60,

что в рассматриваемом случае A = B = R , &i = 0, &2 = — и согласно фор-

Я

мулам (4,11,7)

г, =

да

,

1

/

ди .

\

да

So — — (

---- 4~ Ш) ,

 

2 л

V ар ^

)

 

1

d*w

 

1 /

d2w .

\

1

Я2

да2

2

я 2\

ар2 ^

J

нетрудно выписать

выражения

компонент

приведенной

нагрузки;

АХ, AY, AZ, воспользовавшись формулами (4,11,6):

 

ГДХ = Tl~ Ti

. —

( — + w ) ,

 

 

 

R2

да \

 

ар

)

 

I

ДY =

T J— Г2

 

д2и

 

 

 

я2

 

'

ааар

 

 

 

 

 

 

=

- • ^ [ г . - 5 г

+

т ,

/

a2cy

. \ ]

(4,12,4).

 

 

да2 '

" ч

 

ар2

 

 

Таким образом, принимают вид

d2a

aa2

1 — V*

Eh L

1 + V

2

для цилиндрической оболочки уравнения (4,11,4У

1 — v

 

Фи

1

v

 

a2ti

dw

2

 

ap2 +

 

2

* aa ар

да

+

1

ГЭ

1^

'

dv

,

пг

д*и

Н

R2

da

\

CQ.

 

 

dt2

 

 

 

 

 

д2и

 

d2v

,

1 — v

d2v ,

dw

 

ааар

 

а р 2

^

2

~

а а 2

+

" а р

 

 

1 - V2

Г Т1— Т2

а2^

• /72

а*р

(4,12,5)-

 

 

ЕЛ

I

Я2

ааар

а/2 ]-°-

 

 

 

 

v - £ - +

* r + c W

”’ +

- a r

{^[т‘^

+7'<Ж+1")] +

а« ,

dv

. 2 2 2, 1

1 — V2

 

 

 

 

-f- tn

Решение системы уравнений (4,12,5) ищем в виде

и= U (/)cosraacosAP,

v= V (t) sin па sin /ер,

w— W (0 sin na cos &p,

in R

/

где п = ----- ,

причем i

и k принимают целые положительные значе­

ния. Здесь i показывает число полуволн в меридиональном направле­ нии (/— длина оболочки), k дает число полуволн в окружном направ­ лении. Решение в форме (4,12,6) соответствует тому случаю, когда на

концах оболочки

(z = 0 и z = l) исчезают как радиальные, так

и окруж­

ное перемещения

(иФ 0).

что

уравнения

Легко убедиться непосредственной подстановкой,

(4,12,5) тождественно удовлетворяются, если функции

U(t),

V(t), W(t)

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

 

_ т ( 1 - ^

J

^

+ ( n t + ± J Lk, \ u _

l ±

v nkV - v n W -

 

Eh

 

dt*

 

\

 

2

)

2

 

 

 

 

 

 

1~ *

--1~ Тг n(kV + W) =

0,

 

 

 

 

 

 

Eh

 

R *

K

/

 

 

 

 

- m ( 1

~

v2)

 

E tz

пш

+ —fk2 +

2

 

)

v + m л +Л

 

Eh

 

dt2

 

 

2

 

\

 

 

 

 

 

 

 

1- V

2

JTi-Tg nkU = Q

 

 

(4,12,7)

 

 

 

 

Eh

R

 

 

 

 

v

m

(1 — v2)

 

dPW

— vnU +

kV + c2 {n2 + k2)2W —

 

 

 

 

 

dt*

 

 

 

 

 

 

 

 

 

 

-

 

 

1 7

n « +

7 (jfe*_

1 )]

=

0 .

 

 

 

 

 

Eh

R 2 1

1

v

/J

 

 

 

 

Система уравнений (4,12,7) может быть представлена в матричной

форме

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tn~7T + ^ - T i S 1- T 2S2)f = 0,

 

 

 

 

 

dt2

 

 

 

 

 

 

 

 

где / — вектор с компонентами U, V, tt7,

 

 

 

 

 

 

 

 

п2 +

1~ v k2

— L+-X. nk

 

— vn

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

— v/г

 

k

с2 (n2+

/г2)2

 

 

 

0

 

n k

n

 

 

0

 

n k

n

 

S i= —

n k

0 0

S ,= A n k

 

0

0

 

R 2

0

 

0

 

 

R 2

0

 

0

k 2— 1

 

 

 

 

n 2

 

 

 

 

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

|£ — р2£ | = 0,

а критические параметры продольной сжимающей и радиальной нагру­ зок — из уравнений

= 0 и \R — qRS2\ = 0.

Задача динамической устойчивости приводит в первом приближе­ нии к уравнению

-62£

= 0. (4,12,8)

R ~ ~^R ~ (PO ± T p‘) S l ~ ( 4° ± ' Y q‘) RS2~

 

Если можно пренебречь влиянием тангенциальных сил инерции и тангенциальных составляющих приведенной нагрузки, то задача о ко­ лебаниях цилиндрической оболочки приводится к одному «разрешаю­ щему» уравнению. Для непологой Цилиндрической оболочки оно имеет вид

 

 

д2

 

д2

( у 2 + 1 )2 у 2у 2ф — 0

— v ) ~ £ ^ (

да2

 

ар2 ) у 2ф +

1 — v2

а4Ф _____ z>4=

о

(4,12,9)

 

да4

D

 

 

 

 

Здесь

у ^ 2*!) = до.

Используя (4,12,9) и последнюю из формул (4,12,4), получаем

(V s +

1)« v

V

®

- 0 -

 

V) £

( • £

■—

 

V1® +

 

 

 

 

1 —

V2

аа4 + ■тЯ4D

 

а2

 

 

 

 

 

 

 

 

 

 

 

 

а/2 у 2у 2ф +

 

 

 

 

+

^ -

[

7 ' ‘ ^

 

+

д2

+

1

^

+

1

) ]

л

’ ф “ ° -

( 4 л 2

 

7 '

'

(

Пусть К— длина

полуволны

в

меридиональном

(или

окружном)

 

направлении. Тогда первый член в уравнении

(4,12,10)

будет иметь по-

 

рядок т второй

член — ^

А V

 

 

 

Ш2

. Если длина полу-

 

 

Я. ) третий

 

 

волны мала по сравнению с радиусом, то вторым слагаемым в

(4,12,10)

 

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

 

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

 

 

у 2у 2у 2У2ф +

1 С*

а4Ф

 

R*_ (Т \_

_А_. +

 

 

д2

 

 

а<х*

 

 

D \

R2

да2

^

Я2

* ар2

 

 

 

 

 

 

+

 

 

 

 

0.

 

 

 

 

(4,12,11)

 

Это уравнение полностью соответствует известному уравнению для весьма пологих оболочек.

Итак, возвращаясь к общему уравнению (4,12,10) и полагая в нем

Ф (а, Р, t) f{ t) sin /га cos/ф,

что соответствует случаю (4,12,6), приходим к обыкновенному диффе­ ренциальному уравнению

т Я 4 d2f

где

(я2 + & +

I)2 (/X2 - f А?)2 + (1 — V) л 2 ( л 4 б4) +

- —

с2

л4

g (п, k) =

 

 

 

 

 

(4,12,13)

 

( л 2 +

А2)2

 

 

 

 

 

 

 

 

 

.Введем обозначения:

 

 

 

 

 

 

 

D g (n ,k ) =

P g (n , Л)

= ^ ^

D g (я,

k)

 

(4,12,14)

mR4

 

я2Я2

 

(А2 — I)/?2

 

 

 

 

 

Уравнение (4,12,12) может быть теперь записано в виде

 

 

*L

 

 

 

 

 

(4,12,15)

dt2

 

 

 

 

 

 

 

 

 

 

 

где pt — 2nRTi»,

R

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

задача сведена к извест-

 

 

 

 

 

 

 

ному уравнению.

В заключение отметим, что эти результаты легко обобщаются на случай ортотропной цилиндрической оболочки. В самом деле, пусть Е1, Е2 и V I , V2 — модули упругости и коэффициенты Пуассона в направ­ лениях p= const и a = const соответственно; G — модуль сдвига. Введем следующие дифференциальные операторы:

 

 

 

 

2

 

г

а 2

 

,

 

а 2

 

 

 

 

 

 

V° ~

0 да2

^

2

ар2

 

 

 

 

 

 

 

 

а4

 

+ £4 -

*

 

I

£ —

 

 

 

 

 

 

 

да4

 

 

да2Эр2

 

 

2 Эр4 ’

 

(4,12,16)

 

 

 

 

д*

 

д4

 

+ £

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

да4■+ £3

da2 ар2

 

 

2 Эр4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ 0 — 2G (1 —

 

“Ь ^ 2 v i>

 

 

 

 

 

 

£ 3 =

4G (1 — v xv 2) - f £ XV 2

+

£ 2V

 

(4,12,17)

 

 

 

£4

 

G

£ I V2 —

£ 2VI .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение,

аналогичное

(4,12,11),

для

ортотропной

оболочки прини­

мает

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4у4ф

,

£A (l-viv«)

. *Ф_

^

 

 

 

 

 

V‘V*

+

С2

 

 

 

да4

 

 

 

 

j ^ - ^ z

^

 

^

 

 

 

 

 

 

0

(4 12,18)

 

^

с2л

\ 1

За2

 

ар2

 

 

 

а/2 у

Vl

 

v

где

Ф — функция,

связанная

с радиальным

перемещением

соотноше­

нием (4,12,9).

 

 

 

 

 

 

 

 

 

 

 

 

 

Дальнейшие преобразования вновь приводят к уравнению (4,12,15),

коэффициенты которого для случая опертой оболочки определяются без труда [24].

§ 13. СФЕРИЧЕСКАЯ ОБОЛОЧКА

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

 

 

 

=

(q0 + q tcosQt).

 

 

 

(4,13,1)

 

 

 

 

Обозначим

R — радиус

срединной

по­

 

 

верхности,

h — толщину

И'примем

географи­

 

 

ческие

координаты

ф,

|3

(ф— угол

широты,

 

 

Р — угол долготы; рис.

4.2). Перемещения

то­

 

 

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

 

 

через

и — в направлении

линии ф= сопз1:,

че­

 

 

рез V в направлении линии p= const и w

 

 

радиальное перемещение, положительное в на­

 

 

правлении

внешней нормали.

 

 

 

 

 

 

 

Известно [7], что в случае только радиаль­

 

 

ной нагрузки система уравнений сферической

 

 

оболочки сводится

к одному

разрешающему

 

Рис. 4.2

уравнению

 

 

 

 

 

 

 

 

 

[с2 (v2 +

I)2 - Г

1] (V2 + 2) w =

 

(V2 — 1 — V ) .

(4,13,2)

Здесь

у 2( )— оператор Лапласа на сфере

 

 

 

 

 

 

 

v “( ) =

эшф i

дф

\

 

W — !— 4 - 1 ,

 

 

 

 

 

дф

)

sin2 ф

dp2 J

 

 

 

 

 

с* =

 

/г2

 

 

 

 

 

 

 

 

 

 

12 (1 — v2) R 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В

безмоментном

состоянии

внутренние

силы

оболочки

сводятся

к сжимающим усилиям

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4,13,4)

 

Т г= Т 2 = — R{q0 + qt cos Of).

 

 

 

В случае сферической оболочки

 

 

 

 

 

 

 

 

 

 

kx= k2 = — , А = R,

В = R sin ф,

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

и формулы (4,11,7) принимают вид

 

dPw

 

 

 

 

 

 

 

 

 

= ----

 

 

 

 

 

 

 

 

 

 

 

"аф2"

 

 

 

 

 

 

 

 

 

 

 

R2 Ч

 

 

 

 

 

 

 

 

х2 = ---—( а

sin2 ф

ар2

' «

I

j r

)

 

 

 

2

R2 Ч

 

 

 

Причем, согласно (4,11,6) дополнительная приведенная нагрузка, воз­ никающая при отклонении обрлочки от безмоментного состояния, будет

&Z = - jR ( q 0 + qtcos6/) (xx + x2),

И Л И

AZ = — (q0+ qtcos 0/) (v2 + 2) w.

(4,13,5)

Составляющие

АХ и АУ равны нулю. Кроме внешнего давления (4,13,1)

и приведенной

нагрузки (4,13,5) на

оболочку

действуют

силы инерции

 

т -----д2и , — т -----

d2v , — т

------dPw .

(4,13,6)

 

di2

dt2

dt2

 

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

Z = — (q0 + qtcos 00 — (<7о + qtcos 00 (V2 + 2)w — m

. (4,13,7)

Первое слагаемое дает равномерное сжатие оболочки и может быть отброшено, если под w (г|э, р, t) понимать отклонение от невозмущенного безмоментного Состояния. Тогда уравнение (4,13,2) принимает вид

Jc2(V 2 + I)2 + 1] ( V 3 +

2) w +

J?i±i^£L0Otf_ (V2 + 1 _ v) (v* +

2) а, +

 

mR2

/ о

,

,

ч d2w

п

(4,13,8)

 

Eh

(V

+

1 — V) —— =

о.

 

 

 

 

dt2

 

 

Будем искать решения уравнения

(4,13,8) в классе функций

 

 

ш (ф ,М ) = /(0^(Ф ,Р),

 

(4,13,9)

где f(t) — неизвестная

функция

времени; Р(ф, р)— решения

дифферен­

циального уравнения

 

 

 

 

 

 

 

 

 

V*F + XF = 0,

 

(4,13,10)

удовлетворяющие граничным условиям для до (т. е. условиям непрерыв­ ности и однозначности на сфере). Подстановка в (4,13,8) дает после

•сокращения на /■'(ф, Р)

 

 

[с2 (X — I)2 + 1] — 2) / — - (t?a + q‘cos 90 R

(k — 1 +

v) (A, — 2) / -f-

2Eh

 

 

+ ^ . (X _ 1 + v ) - g -

= °.

(4,13,11)

Введем обозначения:

 

 

* • = R ( t l ? + - i c l <x - 1>a + 4

<4 ’ 1 3 ' 12)

и, внеся их в (4,13,11), получим

----- ?0 + ftyos,6'-. у = о,

(4,13,13)

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

ницы главных областей неустойчивости могут быть найдены по извест­ ным приближенным формулам. В частности, низшая граница

02 = 4р211 — ь + Т *

(4,13,14)

Яо

 

Для практических приложений интересно знать огибающую ниж­ них границ областей неустойчивости. Предположим, что параметр X мо­ жет принимать любые вещественные положительные значения, т. е. до­ пустим, что уравнение (4,13,10) имеет сплошной спектр собственных значений. Для ограниченной области, какой является сфера, спектр соб­ ственных значений дискретный. Но в окрестности интересующих нас значений А, спектр уравнения (4,13,10) все же достаточно «густ», чтобы погрешность, следующая из сделанного допущения, была невелика.

В дальнейшем обозначим

 

 

 

02 =

q(X).

 

 

(4,13,15)

где

 

 

mR* 4 V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q(k) = ( X - 2) с2 (A, — l)2 + 1

 

(<7o+

2

(4,13,16)

 

 

 

 

2Eh

 

 

 

X — 1 4- v

 

 

 

 

Для

определения

огибающей

положим

dq

= 0, отсюда

получим

уравнение для X. Рассмотрим случай

 

 

dX

 

значений

достаточно больших

Л>1. Тогда

 

 

 

 

 

 

 

 

 

q (X )^ c W

+ t

( ‘7в+

2

q ) R

(4,13,17)

 

 

 

 

 

 

2Eh

 

 

 

 

 

dq

 

 

 

 

 

 

 

откуда корень уравнения -2 - = 0:

 

 

 

 

 

 

 

 

dX

 

 

 

 

 

 

 

 

 

 

(яо+ 2

* )

'

 

 

(4.13.18)

 

 

 

 

4Ehe*

 

 

 

и, следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

02:

4£/t

1 —

( . + - Н * R *

(4.13.19)

 

mR2

 

16£^/i2c2

 

Введем обозначения:

 

 

 

 

 

 

 

 

 

 

4Ehc

= Я -

Eh

 

 

 

 

(4.13.20)

 

 

R

mR*

 

 

 

 

 

 

 

 

 

 

 

 

Формула

(4,13,19) принимает вид

 

 

 

 

 

 

 

02=4р2

 

 

 

 

 

 

(4.13.21)

Конечно, следует учитывать, что q** представляет собой приближенное (в смысле сделанных допущений) значение минимального критического давления. Действительно, учтя, что

с2 = — —Л*——

12/?2 (1 — v2) ’

получим хорошо известную формулу

2Eh?

1

(4,13,22)

R*

' у/ 3 (1 — v2)

 

Наконец, займемся определением параметра К в общем случае. Уравнение (4,13,10) на сфере

_ I

_

p

L

дф ,

kF = 0

(4,13,23)

sin “ф

[_

дф

\

sin2 г|) d(i2

 

приводит, как известно, к сферическим функциям. Решение уравнения (4,13,23) находим в виде

К ( ф , Р) = Р ( ф )

s in ftp

(4,13,24)

 

cos ftp.

Условие однозначности Р(ф, р) на сфере требует, чтобы ft было целым числом или нулем (ft = 0, 1,2...). Подстановка в (4,13,23) дает

— —

Г — — ( s i n ф )

] +

( ъ

& ) р = 0.

(4,13,25)

sin гр

iL avdV

dtyагр Jj

\

sin2 тр /

 

Полагая х = созф, приводим уравнение (4,13,25) к виду

г[<1-ч-£Кх-гг5)я“0- <4’13’26)

Это известное уравнение для присоединенных полиномов Лежандра. Уравнение (4,13,26) имеет собственные значения

Х„ = п(п+ 1)

(п = 0 ,1 ,2 ,...) .

(4,13,27)

Каждому собственному значению Я„ соответствует (га+1)

собственных

функций:

 

 

р*(*) = (1 _ * 2)*

~~

Рп{х) (ft = 0, 1,2, ...,п),

(4,13,28)

где

dxk

 

 

 

 

 

dn

 

Рп(*)

 

 

 

2Пп\

dxn

 

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

для

уравнения (4,13,23)

систему его

решений:

 

 

 

 

ft = 0

Р0(ф, Р) = Pn(cosij3),

 

ft = 1

Р_!(ф, Р) = PJf> (cosф) sin Р,

(4,13,29)

 

р1(ф, Р) = Р„ *(cos ф) cos Р,

 

»

»»

k = п F-n(ф, Р) = Я л) (cos ф) sin /гР,

Fn(ф, P) = Р%] (cosy) cosnp.

Известно, что полиномы Лежандра Рп(х) имеют в интервале измене­ ния ф(0, я) ровно п нулей. Присоединенные функции P {nk) (х) имеют со­

ответственно (n—k) нулей.

в нуль на

2k

меридианах,

Так как sin&p и cos&p обращаются

Рп] (х) ввиду только что сказанного — на (пk)

широтах, то

вся

сфера разбивается на «клетки», внутри

которых Е(ф, р)

сохраняет

по­

стоянный знак. Это значит, что число Я определяет вид формулы коле­ баний и, в частности, размеры «полуволн» в меридиональном и широт­ ном направлениях. Чем меньше размеры полуволн, тем больше пара­ метр Я. В этом случае разница между двумя соседними собственными значениями становится малой по сравнению с их величиной, чем оправ­ дывается допущение о непрерывности йзменения Я.

§ 14. ФЛАТТЕР ОБОЛОЧЕК И ПЛАСТИН. ПОСТАНОВКА ЗАДАЧИ

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

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

Теоретическое исследование панельного флаттера в корректной, в физико-математическом смысле, постановке стало возможным после того, как в 1947 г. был установлен закон плоских сечений в аэродина­ мике больших сверхзвуковых скоростей [25].

Анализируя движения тонких твердых тел с большими сверхзвуко­ выми скоростями в различных средах, А. А. Ильюшин обнаружил сле­ дующее общее свойство, названное им законом плоских сечений: «Если вектор скорости какой-нибудь точки тела правильной аэродинамической формы есть V и если поперечные скорости других его точек порядка не более &V, то при установившемся и неустановившемся движениях тело вызывает в окружающей среде только поперечные возмущения, причем давление в любой точке поверхности тела, рассчитанное согласно

этому закону, может отличаться

от

истинного на

величину порядка

не больше

 

 

 

 

1 + /2 =

_L(V + _L')

 

2М2

2

V

М 2 J

 

 

 

 

 

по сравнению с единицей». Здесь М2 = -5—>1 есть число Маха, / ----------

 

 

 

Щ

Vo

параметр Ильюшина, выражающийся через скорость движения тела V, наклон нормали к площадке е и скорость звука в невозмущенной сре­ де и0, т. е. скорость звука в газе на бесконечности; этот параметр имеет фундаментальное значение, поскольку в линеаризированной и нели-