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

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

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

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

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

чтобы удовлетворить уравнениям

движения.

 

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

чине с деформациями изгиба и

практически эти

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

узкой полосой у краев.

 

характера колебаний, возникающие

Если представить себе постепенные изменения

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

§ 5. ПОЛОГАЯ ОБОЛОЧКА

Дифференциальное уравнение собственных колебаний пологих обо­ лочек и слегка искривленных пластин в наиболее общем случае можно представить в виде [7]

(4,5,1)

Здесь

a F = F ( x ,

у , / ) — функция

напряжений; у — объемный вес материала

оболочки;

С — коэффициент

упругого основания; g — ускорение силы

тяжести.

Дифференциальное уравнение (4,5,1) .с соответствующими гранич­ ными условиями, присоединенными к нему, полностью описывают коле­ бания пологих оболочек как обыкновенных (с = 0), так и лежащих на упругом основании, находящихся под действием заданной внешней по­

верхностной нагрузки Z = Z ( x y у , t ) и испытывающих

заданные внут­

ренние осевые (нормальные и сдвигающие) силы № \у N%t

S°.

Эти силы могут быть как статические (не изменяющиеся во време­

ни), так и динамические (изменяющиеся во времени по заданному, на­ пример, периодическому закону).

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

лениям заданные сжимающие силы N°i = const, N% = const, будем иметь дифференциальное уравнение

Частные интегралы этого уравнения при указанных выше гранич­

ных условиях и при k% N°i, N 2, С, h, рассматриваемых как по­ стоянные величины, могут быть определены в такой форме:

Внося (4,5,4) в (4,5,3) и совершая согласно (4,5,2) необходимые операции дифференцирования, получим после простых преобразований формулу для всех частот собственных колебаний пологой, шарнирно опертой на прямоугольном контуре оболочки:

yh_ a 2 = ------------------

1----------------

I

------ е * ------

r / ^ L V + Г — Т Г +

2 тп

 

 

12 (1 — v2)

[ U /

U Л

+ E h [ * ■

(

™У

Ч

т

О -

NТ ] (■f + У -c

N ’ ( T - )

T

<4 ’ w

>

Полагая

в формуле

(4,5,5)

№\ =

N 2 =

С = 0,

получим

соотно­

 

шение для

частоты

собственных

колебаний

простой

пологой

оболочки,

 

т. е. оболочки, не испытывающей заданных начальных напряжений и

 

свободно опертой только на прямоугольном контуре (упругое основа­

 

ние при С = 0 отсутствует):

 

 

 

 

 

 

 

yh

2 _

Eh?n*

~g~

m/l “

12(1 — v2)

f),+(t)‘]+a [ЧтУЧтУУ[(тУЧтУГ

(4,5,6)

Наименьшая частота, соответствующая основному тону колебаний, при

т = п = \ будет

yh

,

£ А*я«

(4,5,7)

®т I =

----------------12(1 — v2)

g

*4

 

Если оболочка очерчена по части сферической поверхности, то при

k x = ^2 = —

получаем более простую формулу:

 

 

R

 

 

 

 

 

 

 

Y*

Ш2

£А3

[ ( тп у

/_яя

\*1 *

Eh_

(4,5,8)

g

т'п

12 (1 —v2)

а ) Ч

Ь

) J

^

'

Формула (4,5,7) показывает, что наименьшая частота колебаний будет иметь наименьшее значение для оболочек отрицательной гаус­

совой кривизны [7].

U fa

Эта частота при —\ = -----L совпадает с частотой собственных

сг

колебаний шарнирно опертой на всем контуре прямоугольной пласти­ ны:

y>Lu2

Ehw /_L + J_y

g

*4

12(1 - va) l a ‘

i’ J

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

§ 6. ПОСТАНОВКА И РЕШЕНИЕ ЗАДАЧИ О СОБСТВЕННЫХ ПОПЕРЕЧНЫХ КОЛЕБАНИЯХ ПЛАСТИНЫ

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

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

+

Г D n — + 4Dle

 

 

+ 2 (D12 + 2 D j

dx*dif

+

dP

yh

L

dx*

16 dx3dy

^ K 12

M'i

 

 

 

i

A n

 

|

p.

*1

л

 

 

 

 

+

4D 2e^

3 -

+

D ’ ’ - ^ r J

= °>

 

 

для ортотропной пластины оно примет вид

^ L + ^ ( D l J ^ L + 2D 3 - ^ - + D 2 - ^ ~ ) = 0. dt* yh \ 1 дх* дх2ду2 ду* I

В случае изотропной пластины, как известно,

Di = D2 = D3= D =

ЕЛ3

12(1 — va)

и (4,6,2) принимает вид

d2w .

gD

4

Л

-------- h -

V w = 0.

dt2

yh

v

 

(4.6.1)

(4.6.2)

(4,6,3)

(4,6,4)

Прогиб w должен удовлетворять граничным условиям, зависящим от способа закрепления контура пластины, и начальным условиям

w = ш0(X, у ),

= v 0 (X, у ) при t = 0,

(4,6,5)

 

01

 

где w 0— заданный начальный прогиб в точке (х, у ) .

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

кой.

Опишем кратко ход решения методом Фурье задачи о свободных поперечных колебаниях пластины.

W m n (x, у ) :

Введем дифференциальный оператор L :

L “ 0 ,1

+ 4 D “

а л Г

+ 2 ( 0 ,! + 2£>*')

+

 

+ 4D " - £ F + d * £ ; '

 

<4 Д 6 >

Пользуясь (4,6,6), перепишем уравнение (4,6,1):

 

 

 

+ - ^ - L w

= 0 .

 

-(4,6,7)

 

dt2

 

yh

 

 

4

Представим решение этого уравнения в виде произведения

 

w = (A cos p t +

В sin pi) W (x,

у).

$ , 6,8)

Здесь

 

 

 

 

 

 

есть частота собственного колебания пластины.

 

 

Внося (4,6,8)

в (4,6,7), получим для

W уравнение

 

 

L W

- p 2 ^ - W

= 0.

 

(4,6,9)

 

 

 

8

 

 

 

Найдем решение этого уравнения,

удовлетворяющее

граничным усло­

виям. Граничные

условия для

функции

W (х, у )

в случае заделанного,

опертого и свободного края не отличаются от известных условий для прогиба w ( x , y ) при равновесии пластины.

Представим решение уравнения (4, 6, 9) в виде ряда Фурье, содер­ жащего произвольные постоянные.

Выполнив требования, что функция W (х, у ) удовлетворяет гранич­ ным условиям и является решением уравнения (4,6,9), получим си­ стему однородных уравнений для неизвестных произвольных постоян­ ных; эта система имеет решения, отличные от нуля, если ее определитель А ( р) равен нулю. Отсюда получаем уравнение частот

Д(Р) = 0.

(4,6,10)

Уравнение (4,6,10) имеет бесконечное множество решений,

кото­

рые составляют спектр частот данной пластины. Частоты зависят от

двух параметров т и п ( т = 1 , 2, 3, ...,

п =

1 , 2, 3,

...). Наименьшая ча­

стота называется частотой основного тона,

остальные — частотами выс­

ших порядков, или частотами обертонов.

функция

W m n { x , y ) , которая

Каждой частоте р тп соответствует

определится с точностью до постоянного множителя. Функции W m n ( x , y ) , называемые собственными функциями, определяют формы изогнутой поверхности, соответствующие частотам р тп.

Определение прогиба w ( x , y , t ) в любой точке и в любой момент времени осуществляют так: разлагают заданные начальный прогиб w 0 и начальную скорость v 0 в ряды по собственным функциям

находят коэффициенты < w

и рт „

и берут решение уравнения

(4,6,7)

в виде суммы всех

решений

вида

(4,6,8).

Причем если коэффициенты

c im n , fW удается определить, то нахождение постоянных А и В

не со­

ставит труда.

 

 

 

 

 

Таким образом, общее решение будет

 

 

w =

cos pmnt +

ЭЮPmnt ) W ma(x, у).

(4,6,12)

 

 

 

Рт п

 

 

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

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

 

yh (\Lw +

d2w

 

d2w

+

cw N ~дх?

N„

<4 д 1 з >

Для определения частот собственных колебаний и собственных функций можно применить тот же метод Фурье, который в этом случае приводит к более сложному, по сравнению с (4,6,9), уравнению для функции W ( x , у ) .

В качестве иллюстрации рассмотрим пример определения частот собственных ко­

лебаний и собственных функций методом

Фурье.

 

 

 

 

 

 

 

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

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

лежала

на оси ох, а сторона b

на оси оу\ пластина со свободно опертым контуром. Следовательно, функция

W (х, у)

должна

удовлетворять условиям:

 

 

 

 

 

 

 

 

 

 

 

при х =

0 и х =

а

 

 

 

d2W

 

d2W

 

 

 

 

 

 

 

 

 

W = О,

 

=

0 ,

 

 

 

(4,6,14)

при у =

0 и у =

Ъ

 

 

 

дх2 + v2 ду2

 

 

 

 

 

 

 

 

 

d2W

 

d2W

 

 

 

 

 

 

 

 

 

W =

0,

+ Vi

=

0 .

 

 

 

(4,6,15)

 

 

 

ду2

дх2

 

 

 

Этим условиям

будет удовлетворять

выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

тпх

 

 

 

 

 

 

(4,6,16)

 

 

 

W = sin

а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где т и п — целые числа.

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение движения

(4Д$)

в данном случае имеет вид

 

 

 

 

 

d*W

 

 

d*W

 

d*W

 

 

vh

W =

0.

(4,6,17)

 

 

Dx — ^ 4 - 2 D 3-------

4- D o ------- — p2 - L~

 

 

дх4

 

 

дх2ду2

 

dy*

 

 

g

 

 

 

Его решение будем искать в виде

(4,6,116). После внесения

i(4,6,16)

в (4,6,17),

получим

 

 

пт V

л

/

тпп2 \ 2

/ пп \*

п

yh

л

,

 

 

(— J

+ 2 D 3( - — j

+ D a( —

J - P 2 —

= 0.

(4,6,18)

Отсюда

находим

частоты ртп:

 

 

 

 

 

 

 

 

 

 

 

 

 

л2 * / g _ f

^ { mb V

 

п ( mb \ 2

^

 

\/

У Dl ( т ) + 2 D s a ( ~ )* °2П

(4,6,19)

 

Частота основного

тона

равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ла т г ~ т .. т а *

 

T v*

 

 

 

 

<4’«'2о>

 

 

p"

~

^

V

 

- ^

V

 

 

 

 

 

 

 

 

В случае

квадратной

пластины

со стороной

а частоты

будут соответственно

 

 

Рт п = - J

]

/

у Г

^ Dim* + 2D3m W +

D&*,

 

(4,6,21)

 

 

 

Pii =

-^ -

j

/

"

V D l + 2D3 + Di .

 

 

(4,6,22)

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

 

 

 

 

Р т п = —

1 /

- з г Ц —

j * п * | ,

 

 

(4,6,23)

 

 

 

Ри - S /f [OfM

 

 

(4,6,24)

 

 

 

 

 

 

В случае изотропной квадратной пластины со стороной а частоты будут соответ­

ственно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рт п =

j t

л / ~

g£_

(т2 + п * ).

 

 

 

(4,6,25)

 

 

 

 

 

 

 

а2

У

уh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2я2

,

f~ g D

 

 

 

 

 

(4,6,26)

 

 

 

 

 

 

P il~

 

а2

V

 

yh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Форму колебаний, соответствующую данным частотам рЖп, определяют собствен­

ные функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V /т п — sin

тпх

 

ппу

 

 

 

 

(4,6,27)

 

 

 

 

-------sin

Ь

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

Полный прогиб пластины

представится

 

как

результат

наложения бесконечного

ряда прогибов

вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ -

 

 

 

,

 

 

 

 

ч

тпх

.

ппу

,

(4,6,28)

 

и>тп (Атп cos Pmnt -t- Втп sin pmnt) sin

 

sin

-

а именно:

 

 

 

 

 

 

 

 

 

 

 

 

a

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

OO

00

 

 

 

 

 

 

Pmn .

 

Д .

mnx

.

nny

 

W =

 

 

a m/i COS Pmnt +

 

(4,6,29)

 

 

------

Sin pmnt. Js i n -------

Sin— — .

 

m=1££n= (1

 

 

 

 

 

pmn

 

 

/ 0

.

 

b

 

Здесь Ctm n, Pmn — коэффициенты разложения

начального

прогиба и начальной скоро­

сти в ряды по собственным

функциям;

коэффициенты

рядов определяются

по фор­

мулам

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

а

Ь

 

 

тпх

 

ппу

,

,

 

 

 

 

 

 

 

г

С

 

 

 

 

 

 

 

 

a m n ~ ~ a b

V

)

^ o s i n —

 

sin — ^ — dxdy,

 

 

 

 

 

 

 

 

 

о

о

 

 

 

 

 

 

 

 

 

 

(4,6,30)

 

 

 

 

 

 

а

Ь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ппу

J

 

 

 

 

 

 

Л

 

4

Г

Г

 

« я *

.

J

 

 

 

 

 

Pmn =

— Г-

\

\

«О S14

-------

Sin— ;—

dxdy.

 

 

 

 

 

 

 

ab

J

.J

 

 

 

a

 

b

 

 

 

 

 

 

 

 

 

 

 

о

о

 

 

 

 

 

 

 

 

 

 

 

Выражение собственной функции W (хуу)

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

две противоположные

стороны,

например,

х= 0

и х = а, свободно

оперты,

а

другие за­

креплены как

угодно или свободны, будет

 

 

 

 

 

 

 

 

 

 

 

 

■+■Сзm sh

mzik^y

 

,

mnkiy

 

mnx

 

 

 

 

+

Сш c h --------

 

a

 

 

 

 

 

 

 

а

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

(4,6,31)

 

 

 

— Д А

4-

(

p2yha4

 

A _

 

 

 

 

 

 

2

m*n4g

 

A

 

(4,6,32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ppyha4

 

A

 

 

 

 

 

 

— А А + D2 m4n4g

 

 

 

 

 

Функция

W (x,y)

вида (4,6,31) удовлетворяет условиям

на

опертых

сторонах. Из

условий на

других

сторонах

получается

система

четырех

однородных

уравнений

с четырьмя неизвестными С!т ,

С2т, С3тп,

С47П; частоты определяются из

уравнения,

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

Каждому значению m= 1, 2, 3, соответствует бесконечный ряд частот Р т п , так как уравнение Д(р)=0 получается трансцендентным и имеет бесконечное число ре­ шений.

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

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

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

деляется одной обобщенной координатой q ( t ) .

Прогиб ее представим

в виде

(4,6,33)

w = q ( t ) W ( x i y ) ,

где W ( x , y ) есть функция, удовлетворяющая граничнымA ‘ условиям и изображающая приближенно форму изогнутой поверхности колеблю­ щейся пластины.

Составим уравнение движения Лагранжа 2-го рода при одной сте­ пени свободы, которое имеет вид

(

дТ

\

дТ

дУ = 0.

(4,6,34)

dt \

dq

)

dq

' dq

 

Здесь и далее точка над буквой обозначает дифференцирование по времени; Т — кинетическая энергия и V — потенциальная энергия си­ стемы.

Для ортотропной и однородной пластины Т к V будут

d*W d*W

(4,6,35)

+

Обозначим во

втором из

(4,6,35) подынтегральное

выражение через

V { W ) ; тогда

 

 

 

 

 

v

= - £ - ^ V { W ) d x d y .

(4,6,36)

Уравнение

Лагранжа

(4,6,34)

после внесения в

него выражений

Т и V принимает вид

 

 

 

 

 

q +

p 2q = О,

(4.6.37)

где р — частота, определяемая формулой

 

 

pi = J g _

tt'V(W)dxdy

(4.6.38)

 

 

УЬ

Я W2dxdy

 

 

 

Выражение для прогиба получается в виде

 

 

w =

(A cos р/ + В sin pt) W.

(4.6.39)

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

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

Определим функцию W (х, у ) так, чтобы она давала минимум вы­ ражения (4,6,38). Для приближенного решения этой вариационной задачи воспользуемся методом Ритца.

Зададим выражение для W в виде

^

=

аЕ

» £.л

 

тп

п

 

где Wmn — непрерывные функции, зависящие от двух параметров

т и п и удовлетворяющие всем

граничным* условиям. Тогда (4,6,38)

представится в виде дроби

 

 

_2 _

М (Атп)

(4,6,41)

Р ~

N (Атп)

 

числитель и знаменатель которой — однородные квадратичные функции коэффициентов А тп. Найдем минимум этой функции. Составим частные производные д р 21дАтп и приравняем их нулю. Очевидно,

ар2

1

дм

М

dN

1

Г дм

2 dN \

dAmn

N

дАтп

N*

дАтп

N

V дАтп

дАщп )

= 1 Г - 1 Г - ( М - Р 2^ -

<4’6’42>

N дАтп

 

Следовательно, задача отыскания минимума дроби (4,6,41) равно­ сильна задаче отыскания минимума выражения

S = M — p2N ,

(4,6,43)

или

5=И[^ F{W)~p2W]dx y■

(4,6>4)

Внося сюда значение W по формуле (4,6,40), получаем 5 в виде одно­ родной квадратичной функции коэффициентов А тп. Отыскание мини­ мума функции S приводит к системе однородных уравнений первой степени относительно коэффициентов Атп\ из условия равенства нулю определителя этой системы получим уравнение частот:

Л(р) = 0.

(4,6,45)

Наименьшее решение этого уравнения, отличное от нуля, и будет при­ ближенным значением частоты основного тона.

Для иллюстрации метода Релея—Ритца рассмотрим пример. Опре­ делим частоту основного тона прямоугольной ортотропной пластины, у которой все четыре стороны жестко защемлены. Точное решение этой задачи нам неизвестно. Приближенное решение получим, полагая

W mn = ( > - ~ ~ У ( у 2 - -£-)* *"Ч/П,

(4,6,46)

где т , п — целые числа.

Функции W mn будут удовлетворять граничным условиям:

 

I

о.

 

 

 

 

 

 

 

 

при х — ± —

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wmn =

=

0,

 

(4,6,47)

при у =

±

Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w mn =

=

0.

 

(4,6,48)

 

 

 

 

 

 

°У

 

 

 

 

В

качестве

первого

приближения

возьмем первый

член

суммы

(4,6,40):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W = А (2 - -^-)2 ( у 2 - - j

 

(4 ,6 ,4 9 )

Внося это в (4 ,6 ,3 8 ) и производя интегрирование, получим

 

 

 

 

f t . =

- Т

Г - 1

+

 

° ’571

 

t 4 ’6 ' 50»

Для квадратной

изотропной пластины с жесткостью

D получим

 

 

 

 

 

 

 

 

 

 

(4Д 51)

Можно указать и иное выражение для W mn:

 

 

 

 

W

= |^1— (— 1)'»cos 2тпх

£l — (— l)"cos —jp -J,

(4,6,52)

 

 

w тп

 

 

 

 

 

 

 

 

которое также удовлетворяет всем граничным условиям. Первое при­ ближение для частоты основного тона получим, полагая

W = A h . +

c o s ^ - ^ 1 + cos

(4,6 53)

после вычислении находим

 

 

 

Рц = J I J L

i /

Z

i /

(4,6,54)

Ьг

У

yh

У a4

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

(4,6,55)

Сравнение выражений (4,6,51) и (4,6,55) дает незначительную разницу; она составляет 3,33% по отношению к меньшей.

§ 7. в ы н у ж д е н н ы е к о л е б а н и я о б о л о ч к и , п о с т а н о в к а з а д а ч и

Выпишем систему дифференциальных уравнений осесимметричных гармонических колебаний тонкой упругой оболочки вращения в линиях главных кривизн, которая, как известно [8], имеет вид

A da \ АВ J da

+(1-v)Gsk

1

dw \

 

 

AR2

da )

 

) +

i r { f )

 

 

+

A

da

[

\ Rl

 

R2 )

J

 

 

 

(1 - -v®)p®

 

,

/la

i

1

j

_

( '

1

 

to

 

H

3

l

ARi

9

da

\ .

AB

 

о

 

 

dw

 

 

 

 

 

dw

,

 

 

 

da ] +

K i s , <

 

 

+ t )

} +

1 — у®

Х =

0,

 

 

 

 

 

 

(4,7,1)

2Eh

 

 

 

 

 

 

 

 

 

 

АВ

f —

+ —

Y dBu

 

1 — V .

j L

/ ^

_

/ _ L

+

_ 2 v _

+ _ L \ w +

4 Ri

R2

J

da

 

AB

d a \ R 2 )

{ R *

 

 

R XR 2

R \ J

 

 

( l - v ® ) p ®

 

h2

( 1

d / В \

d f

1

\

d Г Bu \ ,

 

 

 

 

w ---------- \ ------

 

 

 

 

 

 

 

 

 

 

Rl

 

 

3

(A B

da

U

/

da \

AB

)

da

V Ri ) +

 

 

+

1V

d (

В \ / 1 dw

и \ .

 

 

 

 

AB

d (

 

 

Д A ' da

R i J

 

 

 

 

 

da V\

 

 

 

 

. —

f — ^ ~ ~ (

—■— ^ —— (

') — —1 +

Z = 0

( p

Р<ааЯо

AB da

\ A

J

da \

AB J

da \

A

J da I

2Eh

\

 

Здесь Ro — некоторый характерный радиус кривизны оболочки, р — без­ размерный параметр частоты колебаний; остальные обозначения обще­ приняты; множитель cosco£ при искомых величинах опущен.

Соотношения упругости для тангенциальных усилий имеют вид

2Eh

2Eh

1 — V®