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

книги / Прочность, устойчивость, колебания. Т. 3

.pdf
Скачиваний:
25
Добавлен:
13.11.2023
Размер:
24.79 Mб
Скачать

Свободные колебания цилиндрических оболочек

441

в перемещениях в упрощенной форме с учетом параметрических членов имеют вид (Д^,2 = 0)

ен

\ ( & , 1—у

а2 \

и*

1 + у ш

1— V2

[ \а*3

*"

2

*

дз2

)

2

а* а5

 

 

.

V

 

дш 1 .

д-их

 

 

 

 

+

Я

'

дх

1 ~

 

 

 

 

Е1г

Г

1+ V

 

а *иг

/ 1

 

V а2

а2 \

1 — V2

ь

2

* дхдза*а$+

\

2

а*2 +

д р 7 из +

 

 

 

1

 

диз

 

 

д2иг

 

 

 

 

+ Я 'И Г Ь ' Ф ;

 

 

ей

/V

а%

,

1

а«2

,

а» , г0По*а \

т=^г Ы

-а г+ т г -эг + ж

+6'Л'ДАш) ■

 

 

а2и>

" 22

а2ау

=-рЛ

д2ю

*

 

 

дх2

а52

^

а/2

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

 

 

Л

 

тпх

71ПЗ

 

 

1*1

=

С 0 8 —

^

С 05 — ^ —

 

 

и2

,

лтх

лпз

 

 

= А2 з!п — -— 81П —

 

 

 

.

.

птх

тшз

 

 

>= А3

51П ----^----С05---- .

 

Если пренебречь по сравнению с

единицей величинами

порядка

ЬЩ}

дз/гзроз

 

к*р<йЩ/

 

~ВГ ’ 3 также ----- Е

и -----11

, то приближенное значение

частоты можно определить по формуле

 

 

“ а = ' ^

' (Я‘ +

бг(^

+л!!)4(1~ ' ' 2 Г ,+

(26)

+

(А3 Н-я3)2 [А3Фп +

(я3 -

1) Фм]) [(А.3 + я3)3 + п3]-‘

 

 

( 1 -

 

 

 

 

 

Формулу (26) можно представить также в виде ©2 = (4(1+в1<Рп .^ Я,ф22).

442 Колебания оболочек

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

можно записать

 

 

 

 

 

 

©2 =

(1 +

+

§2^22 ~~

| ^ 121);

 

здесь

 

кп (А,2 + л2)2

 

 

 

5 -

 

 

 

 

X4 +

63 (1 — V»)-! (X3 +

п3)1 •

 

Минимальное значение частоты

три заданных фп , фаг, ф12, К,

Л, I

и т достигается при значении п — лоПри фп = Фаз = 0, т =

1 ве­

 

 

 

личину

л о

определяют по фор­

 

 

 

муле

 

1

 

 

 

 

 

 

 

 

 

 

л0 =

4,2(1 — V3) 8 X

 

 

 

 

* (

)+ (

)+

 

 

 

Это значение4

получено4

в пред­

 

 

 

положении,

что

 

 

 

 

 

У*икр

 

а.2 « л|).

 

 

 

Рис. 14

 

 

(й/со0

от отно­

 

 

Зависимость

ражена

на рис. 14 (эта

кривая

шения фха/ф12дгр

при л = л 0 изоб­

практически

не зависит от геометри­

ческих

размеров оболочки); здесь

 

 

 

 

 

 

«Риде» =

0,74 (1 -

V3) 8

( М

- )

4

 

 

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

127, 371.

Влияние поперечных сдвигов и инерции вращения. Учет деформации сдвига и инерции вращения приводит к некоторому уточнению частот, найденных в рамках гипотез Кирхгофа:Лява. Это уточнение тем суще­ ственней, чем.меньше длина полуволны форм колебаний. Кроме того, появляются новые формы собственных колебаний, соответствующие более высоким частотам. Для получения уточненных дифференциаль­ ных уравнений вводим триортогональную систему координат: х — вдоль

образующей; 5 — в окружном

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

средней поверхности;

г — перпендикулярно средней

поверхности. Предполагаем, что пере­

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

виде

 

 

 

 

 

«1 (г) =

иг (х ,

«;

I) +

гр, (х,

«

();

«2(г) = и 2 (*.

 

0 +

гр, (х,

г.

<);

 

«з(г) =

ш (х ,

*; <)•

 

 

Колебания сферических оболочек

443

Система уравнений имеет

вид

 

 

^ 11«1 -I- ^ 12^2 +

^13^ + ^14^1 +

^ 1*02 =

Ф

^21М 4~ ^Й2и2 +

^23^ "Ь ^2401

^2»Рг =

0;

^31^1 ~Ь

+ ЦэМ + ^34^1 + ^'35р2 = 0;

^•41^1 +

^>42^2 Ч~ ^13^ “Ь ^44Р» Н“ ^4Г>Рз —0,

1*Ы И | Н “

^ 5 2 ^ 2

^ 5 3 ^ “ 1“ ^ - 5 4 ^ 1 "1 “

^ В б Р г =

0 *

Результаты

[38] приведенной

частоты

“ »= м ( — )Г~) У ъ

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

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

КОЛЕБАНИЯ СФЕРИЧЕСКИХ ОБОЛОЧЕК

Основные соотношения. Уравнения колебаний. Пусть срединная поверхность (рис. 16) сферической оболочки отнесена к географической системе координат ^ = 0,

*2 = Ф (в — угол широты; <р — угол дол­ готы). Параметры Ламе выразятся фор­ мулами

# 1 = /? ; н г — К 5Н1 0,

где К — радиус кривизны оболочки.

444

 

 

 

 

 

К о л еб а н и я

оболочек

 

 

 

Дифференциальные уравнения

колебаний будут

 

дN

+ " Ш

 

дЯ

 

{Ыи ~ *1п) с‘ е ® _

 

Ж

Т ' ж

+

"

 

 

 

 

 

= рН К ^ - К Ч х ,

 

 

 

ал?,»

,

Ж ё - ^ + 2^ ^ е - о » =

 

ав

+ '

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

Ра д ^ т - -/?</»;

 

 

 

&0\

,

1

 

дфа +

<?1с18Э + (^и + ^ 22) =

<30

 

5Н1 0

 

дгш

 

 

 

 

 

 

 

 

 

=

 

 

/?<73;

 

 

 

 

 

 

 

 

— рнк

 

+

 

 

 

дМп

 

1

 

^

+ <м„ - М»»)с1е 0-

<?1 Й= 0:

50

5Н1 0

 

5Л4

 

 

 

1

 

дМгг

+

2Л 1 „ с» ев -0 » Л = 0 .

_ Л _|____________

50

8Ш0

 

 

 

 

 

 

 

Условия и моменты определяют по формулам

 

 

Рк

 

 

+

' ’е22);

 

 

Рк

 

Vе11);

Ми = у—

^ (еи

ЛГ22 = у^ - у2 (е23 +

Ни = и 21 =

Ек

( в1а

1 2 /? Щг ) ;

 

1 + 7

 

М п

= Э (хи

+ 7х2г);

Л422 = Е>(х22 +

7Хц);

 

 

Л422 = М21 = Э (1 — V) х18.

 

 

Компоненты деформации еа^ и изменения кривизны:

 

 

 

 

е11= х ( ж + “’) ;

 

 

 

 

^=тН^-ж+“,с‘е9+“');

 

е11= 2 ^

( ш

т - ж

+ ж

-

“»с1е е ) ;

 

 

 

 

 

* 11

 

1

/ 52ш

 

5% \

 

 

 

 

 

 

 

:

/?2

\502

 

59 )

'

 

 

 

___1

Г _ 1 _

5 /

1

\

 

*22 — ^2

[ 5ш 0 ’ 5ф \ 81П 0 ’ 5ф ~ ~ и%)

 

 

___1__

 

+

с*е Э (

ж

 

- и1) ]

;

 

 

 

Г

1

 

52ш

 

 

д / _ 1

дои

 

%1Ъ

 

К2

 

1 зт

0 ’ 50 5ф

 

С8 ° \ зт 0

’ 5ф

 

 

 

 

 

 

\

____ 1_

 

 

50 Г

 

 

 

 

 

 

 

V

 

5Ш 0 * 5ф

 

 

(27)

(28)

(29)

Колебания сферических оболочек

445

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

(30)

Соотношения упругости также упрощаются:

^12 — Мц = '

Е12-

Для отыскания преимущественно изгибных форм колебаний могут быть использованы другие упрощенные дифференциальные уравнения. Уравнения упрощают введением двух предположений: 1) влиянием тан­ генциальных сил инерции можно пренебречь; 2) вкладом тангенциаль­ ных смещений в изменение кривизн можно пренебречь. Последователь­ ное проведение этих гипотез и введение функции усилий х приводит к дифференциальным уравнениям

ОДДш + -1-Д х + р Л ^ = 0 ;

(31)

ДД% — К Дш = 0.

Связь усилий в срединной поверхности с функцией % н выра­ жение оператора Лапласа Д даются формулами

(32)

Иногда удобно исключить из системы (31) функцию усилий

1> ДДДиу -1- Дш + рИ Дщ = 0.

(33)

446

Колебания оболочек

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

до = до (г, ф) еш ,

(34)

где г, ф — координаты точки в полярной системе координат. Подстав­

ляя выражение (34) в уравнение (33),

убеждаемся,

что свободная

форма колебаний должна удовлетворять уравнению

 

ДДДдо — 64 Ддо =

О,

(35)

где

 

 

Но решение уравнения (35) является суммой решений дифферен­

циальных уравнений в частных производных

 

Ддо = 0 ; Дю — Ьгю = 0; Дм» 62до = 0.

(36)

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

Свободные формы колебаний удобно искать в виде

до (г, Ф) = / (г) соз гсф.

Тогда решения уравнений (36) будут иметь, соответственно, следу­ ющий вид:

 

/ =

 

 

 

;

 

 

 

 

/ = В г 1п (Ьт) +

В 2К п (Ьг)\

! =

С У л (Ьг) +

С 2 У п (Ь г).

Так как в вершине сегмента смещения должны

быть ограничены,

то А2 — Вй = С2 = 0.

Окончательно,

функция / (г) имеет вид

 

! (г) = Агя + В1п (Ьг) + а

п (Ьг),

 

(38)

Аналогичные выкладки для функции напряжений (усилий)

приводят к формуле

 

%=

*Ф(г) соз пуеш

 

 

 

 

 

 

 

 

 

 

 

+ (О = Ргп +

 

 

г*+* +

р р (В1п (Ьг) -

С1п (И ].

Для

иллюстрации

удобно

рассмотреть

осесимметричные свобод­

ные колебания сферического сегмента

с защемленным

контуром. Гра­

ничные

условия при

г =

г0

 

 

 

 

 

 

 

 

и 1

= «а =

до =

ддо

0.

 

(39)

 

 

 

 

 

 

~дг

 

 

 

Колебания сферических оболочек

447

Соответствующее характеристическое уравнение (уравнение частот)

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

 

 

Л> (Ьг0) и (*/•„)+Л (Ьг„) /о (»/•„) + 4 (1 + У У ± ^ Га)

= °- (40>

 

- к )

 

Для скользящей заделки (осесимметричные колебания)

 

Л Г п = ш = - ^ = 0

 

(41)

уравнение частот будет

 

 

(ЬГо) /1 (Лто) 4 -Л (^о) Л) (^го)---- НЬГц ^1

^1

 

Эта задача рассмотрена Г. А. Ван-Фо-Фы и В. Н. Буйволом [12]. Если в выражении (38) постоянную А положить равной нулю и сохра­ нить из краевых условий (39) или (41) лишь условия

ДО= тг- = 0,

д г

то получится приведенное в работе [28] неточное уравнение

 

К

ы

1'п ( М -

Л, ( Ч ) 'п ( Ьга) = о-

(“3)

Однако

это

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

определении частот

приводит

лишь

к несущественной

погрешности

[12].

Сопоставление

частот осесим­

метричных

колебаний защемленного и

опертого сегментов проведено

втабл. 5.

5.Частоты колебаний со* 10 3 в рад/сек сферического сегмента

радиуса /*„, высоты 5

=

= 60; ~ = 0.984-107 м*/сек: '

Защемленный край

Неподвижно опертый

Н о м ер

 

 

 

 

$

 

часто­

 

 

 

 

го

 

ты

 

5

10

16

5

 

0

0

 

60

60

60

60

 

 

 

1

4,4

9

16

22

2, 1

9

2

17,16

19

22

29

12,76

15

3

38,39

39

40

43

31,85

32

о

16

60

 

16

21

20

29

32

38

448

Колебания оболочек

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

й2иг

2

<1иг

2

2

Л

 

 

}Г*г + * и г - 0

 

 

 

 

 

 

(аа

р^а

\

 

 

Х + 2 Ц / ’

 

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

напряжение

будет

 

стг =

(Я, + 2р)

- у - .

 

Решение дифференциального уравнения может быть записано в об­

щем виде

Л /

Сх 51П аг +

Со со® аг

иг

йг \

 

г

=

 

 

Если гг — внутренний радиус; г2 — внешний радиус полого шара,

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

 

 

1аг1+ (аМ — 0

>8а г 1

Ъ,аг2 + (аУ-2— С) агг

( “ М — ? ) —& Г1

а г ,

( ( ? г \ — 5 ) — ? а л 2 а г г

В случае очень тонкого шара (гг <=* г2 ^ В) приближенное значение частоты может быть найдено в явном виде

гп 1 /

(Х + 2 |1 )С (3 - 0

V

 

 

л Г

2Е

® ”

Г

р/?я (1 — V) *

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

Влияние деформации сдвига и инерции вращения на колебания сфе­ рических оболочек. Учет деформации поперечного сдвига и инерции вращения приводит к уточнению частот и форм колебаний, найденных на основе гипотез Кирхгофа-Лява. Это уточнение тем существенней, чем меньше длина полуволн форм колебаний. Кроме того, появляются новые формы колебаний, соответствующие более высоким частотам. Пусть сферическая оболочка постоянной толщины Ни радиусом средин­ ной поверхности К отнесена к полярной системе координат (г, 0, <р). Решение ищем в форме

“в = Их + гфц иФ = и2 + грх; иг - ш.

Колебания сферических оболочек

449

Уточненные дифференциальные уравнения движения будут

Ж + Ш 9 - 1 & +

с ' е 6 - * =

Ж ~ « к

Ж + Ш о - Ж + 2Ы^

в - ° ^

р Ш ^

- Щг;

Ж + 5ТЙ-0’§ Г + <?«с ‘ 8 0 + <"*• +

^

Ш + * «

ж+ ж е ‘ж +

дМ12

1

•^

+ 2Ми с18 в — Л<?2 == рАЛЧ»*! ф ,

50

"Г 51П 0

где

 

 

 

Л2

 

к = 1

+ 62;

/<! = ! + 1,8б2; 62 =

 

12#2 *

В частном случае свободных колебаний замкнутой сферической

оболочки уравнение частот при к« =

1,2, V =

имеет вид (&2—осред-

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

сдвига)

 

 

 

|

6* [— 3,6М3 -I- (6,8 + 8.2Л0 М3 +

(9,2 — 6,337/ —

 

— 5.6Л/2) М + (4,27 — 7.2Л/ + 0,53/V3 + Щ | +

+

Л42 — М (2 + ЛГ) + -5- N

 

(ЗкМ + 2 — Щ X

 

х ( з М + 2 — ^

 

- * ) = 0 ;

здесь

 

рй)2 (1 — V2) К2

 

 

 

 

М =

 

N = п(п + 1);

 

Е

;

п— номер присоединенной функции Лежандра. Приведенная задача рассмотрена Прасадом [45].

Нелинейные колебания пологих сферических оболочек. Для изуче­ ния осесимметричных нелинейных колебаний пологих сферических оболочек может быть использован метод Бубнова-Галеркина. Пусть сферический сегмент отнесен к системе координат, которая в плане является полярной. Тогда уравнение срединной поверхности будет

Компоненты деформации выразятся формулами

е

ди

5/

дур

1

/ дур \ а

1

11

Аг '

дг +

2

^“д Г ) ;

«22

 

 

450

 

К о л еб а н и я

оболочек

 

Условие совместности деформаций имеет вид

д

1гс ,

. ___V . *и_____}_(Л Е А *

дг

(«и)

«и -

~ЗГ

дг

2 \ дг )

Удобно ввести следующие

обозначения:

 

*

Ю “

т

 

‘ -■

Г1

рЬг\

И

н

ОЯ

в

VI .

-

йН

 

 

ОН

9 4

 

(?— осевая сила, приложенная в полюсе оболочки; ф — силовая функция,

^22г\

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

у д3ш*

д2ю*

 

 

д ?

 

 

д2т*

0 1 .

 

 

2я *

* (Рф

дф

1

 

1 —V2

 

6

+

 

|

^

б2

 

х

Г А{*

дьо*

1

/

дно \ 2 1

= 0.

[ й1 ' дЬ + 2 \

) ]

Удобно исходить из следующих граничных условии при | = 1

• т к м '*

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

а =

0,

р =

0

— жестко защемленный край;

а =

0,

р =

оо — свободно

защемленный край;

а =

оо,

р =

0

— жестко опертый край;

а =

оо,

р =

оо — свободно

опертый край.