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

книги / Общая теория анизотропных оболочек

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
21.14 Mб
Скачать

3 11

СВО БО ДН Ы Е К О Л ЕБА Н И Я

351

Затруднительно также отыскание общего решения уравнения (1. 23). Поэтому ограничиваемся классом решений, являющихся решениями следующего дифференциального уравнения:

+

=

(1.24)

и удовлетворяющих граничным условиям для и>, т. е. условиям непрерывности и однозначности w на сфере.

Из (1.23), согласно (1.24), для частот собственных колебаний трансверсально изотропной сферической оболочки получим

2 _ gE п __р\ с'г (X— I )2+

1 + х*х

 

(1.25)

<0

ТоД 2 {

( l + A * X ) (

X - l + v )

 

 

 

 

Формулой (1.25) собственные частоты оболочки представлены

с помощью пока еще неизвестного параметра X.

 

 

Значения параметра X определим из уравнения (1.24),

которое

подробно изучено.

Решение

уравнения

(1.24)

ищется

в виде

 

 

 

sin к В,

 

 

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

функции W (а., (3) на сфере

требует,

чтобы к было нулем

или целым числом (& = 0 , 1 , 2 ,.. .).

 

Подстановка значения W в (1.24) с учетом обозначения cos а =

= х приводит к известному уравнению

 

 

 

 

- * ’ > - £ ] + ( * - А

) ^

0

 

 

(к =

0 ,

1, 2, ...),

 

 

 

собственными функциями которого являются присоединенные функции Лежандра.

Нетривиальное ограниченное решение этого уравнения суще­

ствует лишь при собственных значениях

 

Х„ =

п(п +

1)

(ге =

0 , 1, 2, ...),

(1.26)

которым соответствуют следующие собственные функции:

 

^

(*) =

( ! - * T

* £ EP. (* )

(1-27)

 

(* =

0 ,

1, 2,

...),

 

где полином Лежандра

 

 

 

 

Ря(х) = 2яп ! ^ К * 2— !)*]■

352

К О Л ЕБА Н И Я И УСТОЙЧИВОСТЬ

О БО ЛО ЧЕК

[ГЛ . I l l

 

Возвращаясь

к

уравнению

(1.24),

находим

систему

его

решений:

 

 

 

 

 

 

 

 

к— 0:

W0(а,

(3) =

Рпcos а;

 

 

 

к=

1:

W_! (а,

Р) =

P^l) (cos а) sin |3,

 

 

 

 

 

Wj (а,

Р) =

(cos а) cos Р;

 

 

 

к — п:

1¥_я (a, ft) =

PM (cos а) sin п р,

 

 

 

 

 

W„ (а,

Р) =

Р№ (cos а) cos п р.

 

 

 

Напомним, что полиномы Лежандра Р п(а:) в интервале

( — 1,

+ 1 ), т. е. в интервале изменения а (0, я), имеют п нулей. Присоеди­ ненные функции Лежандра P ik) (х) имеют п—к нулей.

Как известно, sin п р и cos пр обращаются в нуль на 2п мери­ дианах, а Р (пк) (х) (в интервале 0 ^ а ^ тс) — на п—к широтах. Следовательно, этими меридианами и широтами (узловыми линиями) оболочка разбивается на участки, внутри которых Wn не меняет свой знак. Это значит, что параметр \ определяет вид формы коле­ баний оболочки. Чем больше параметр X, тем меньше размеры по­ луволн.

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

д*ш

I о

 

.

| Д 4

т0 ,

(1.28)

ds

4- 2т

-— г -\-nw -\- -г.-----— h

'

'

1

Du

g

 

где

§ 1]

СВО БО ДН Ы Е КО Л ЕБА Н И Я

353

Решение уравнения (1.28) ищем в виде

 

w{s,

t) = Акs

i n s

i n <s>kt,

(1.30)

где A k — неизвестные

постоянные,

<ой — частоты

колебаний.

Принимая (1.30), полностью

удовлетворяем граничным усло­

виям, а удовлетворяя

уравнению (1.28), получим

для частот

колебаний следующую

формулу:

 

 

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

сматриваемого частного случая. Изме­

 

няя материал оболочки или, что очень

 

важно, ориентацию главных направ­

 

лений упругости материала оболочки,

 

можно

проектировать

оболочки

 

с требуемыми

частотами

свободных

 

колебаний.

 

 

 

 

 

Пусть,

например, оболочка изго­

 

товлена из ортотропного материала,

 

главные направления упругости ко­

 

торого в

каждой

точке

оболочки

 

с главными геометрическими направ­

 

лениями оболочки

(s и 0) могут со­

 

ставлять произвольный угол.

В ча­

 

стности,

пусть

в

каждой

точке

 

оболочки главное направление упру­

Рис. 61.

гости составляет с главным геоме­

 

трическим

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

s

угол ?

 

(рис. 61). Пусть, далее, оболочка составлена из однородного материала (Са:=hBik); тогда для коэффициентов /тайга, входящих

в расчетную формулу

(1.31), будем иметь

 

т __о ^16(Д11^26 —

 

 

Вп{ВпВвв-Щв) '

(1.32)

__12Л- (ВцВ%2 В \г) 8 еs ~Ь ^.В^В^ф ^ъ — Дц^?б — ^22^16

 

/г-

В п (ВцВвд — В j’e)

 

Имея значения упругих постоянных материала оболочки в глав­ ных направлениях упругости, с помощью формул (1.1.42) илй (II.9.38) легко определить значения коэффициентов Bik, входя­ щих в (1.32).

Имея значения всех В(к, с помощью формул (1.32) определим значения /тайга при различных значениях угла <ри далее с помощью

23 С. А. Амбарцумян

354

КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК

[ГЛ. III

формулы (1.31) найдем значения частот колебаний оболочки при различной ориентации ортотропного материала в теле оболочки.

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

§ 2. Некоторые вопросы статической устойчивости анизотропных оболочек

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

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

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

Упрощенная система дифференциальных уравнений устойчи­ вости анизотропных оболочек может быть получена на основании уравнений теории пологих оболочек (см. гл. I, §§ 5 и 14).

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

известные представления:

 

 

 

 

 

Z* — — Т\ч — Т\жг S \

X* = 0,

У* = 0.

(2.1)

Заменяя в соответствующих уравнениях

(например,

(1.5.14),

(1.5.21),

(1.5.27), (1.5.46),

(1.14.12),

(1.14.30), (1.14.41)

и др.)

Z , X , Y

значениями Z * ,

X * , Y *,

получим уравнения локаль­

ной устойчивости соответствующей анизотропной слоистой обо­ лочки.

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

§ 2J СТАТИЧЕСКАЯ УСТОЙЧИВОСТЬ АНИЗОТРОПНЫХ ОБОЛОЧЕК

355

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

Здесь, сознавая незавершенность исследований в этой области, при рассмотрении задач устйчивости оболочек с учетом попереч­ ных сдвигов все-таки будем пользоваться представлениями (2.1),

которые в первом приближении позволяют выяснить влияние по­

перечных сдвигов на устой­

 

чивость

анизотропной обо­

 

лочки.

 

 

 

1. Устойчивость

пологой

 

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

 

панели. Рассмотрим

задачу

 

статической устойчивости по­

 

логой ортотропной

круго­

 

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

 

(i?1=co,

R 2= R ),

сжатой

Рис. 62.

вдоль образующих равномер­

 

но распределенной нагрузкой

 

с интенсивностью р, приложенной по торцам оболочки (рис. 62).

Пусть оболочка изготовлена так, что главные направления упругости материала ее совпадают с координатными направле­ ниями а, [3(А = \ , 2?=1). Пусть, далее, оболочка шарнирно оперта по всему контуру (а = а , а=0; $=Ь, [3=0).

Очевидно, что в рассматриваемой оболочке в безмоментном

состоянии появляются следующие внутренние силы:

 

 

Ц = — Р, Л>= 0, S° = 0 .

 

 

(2.2)

Тогда, согласно (2. 1), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.3)

Разрешающую

систему

получим

из

системы

уравнений

(1.7.56). Принимая

(1.7) и

вместо Z подставляя

 

Z*,

получим

1 d*F , АЗ г <?2

<?2

 

|

,

 

№ (

д* .

 

i T ^ '+ T 2' L ^ z'11+

 

 

 

 

 

То

 

 

 

+ «55 | г)

- т ] * +

р t 1 -

-JJ (% A i +

 

 

-f- аиЬ22) -(- а44а55щ

(L nL 22

£f2)J

 

= О,

(2.4)

(L aL B - Щ

 

 

(amLn +

aHL22) +

 

 

 

 

hi

,,

,

LJll)J 0a2

0,

 

 

~ T aU «55 100 '^11^22

 

где для операторов L ik имеем формулы

(1.9).

 

 

 

23*

356

КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК

[ГЛ. III

Решение системы (2.4) представим в виде

 

F(n,

Р) =

sin

sin

К =

 

f

(a,

p) =

f mnsinXIasinX2p,

Х2 = - ^ ,

 

где m a n — целые числа, которые показывают числа полуволн

по направлениям а и р соответственно.

Принимая (2.5), удовлетворим условиям шарнирного опира­ ю т по всему контуру.

Подставляя значения F и Ч*- из (2.5) в систему уравнений (2.4), получим относительно коэффициентов Fmn и Чгтп однород­

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

£ - = т я г l b [ЦЬи+ 2X1X2612+ х^ 2+

 

+

(аиЦ +

orJ ‘D (Ьпь2, - ft?2) ] +

^

, (2.6)

где,

как и раньше,

 

 

 

 

ди„=

1. +

Tf7 (as6fcn +

аМ

+ flaOssт (ьн6и — ьи).

 

Ьц — Х»Яц -f- Х|566, fe22 =

Х|В22 -{- Х?506, Ь12 =

Х(Х2 (В 12 -f- Soe)-

Формула (2.6) получена на основании уточненной теории, учитывающей явления поперечных сдвигов. В частном случае классической теории, т. е. когда пренебрегают явлениями, свя­ занными с поперечными сдвигами, полагая a44=l/G 23= 0 , я55= = 1/G13= 0, получим

р ;„= - щ - (цьп +

2х м 2+ чъп)

 

^Рввыо______Х(

(2.8)

 

 

 

Д2

ЬпЬТ1Щг

 

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

дой точке плоскость изотропии материала

 

оболочки

параллельна

срединной поверхности оболочки, в силу (1.15) получим

 

.

D

(Xf +

Xj)2

,

E h

X?

 

(2.9)

~

X2

1 + h*a? (Х| + Ц )

'

Д2

(Х| +

XI)2 *

 

 

 

где, как всегда,

 

 

...

l b 2

 

Е

 

 

(2.10)

D

Eh?

 

 

 

 

 

 

 

 

 

 

 

 

 

12 (1 — V2) »

П ~ ' 10 о2

(1 — v2) G ' "

 

 

Рассматривая формулы (2.6)— (2.10),

 

замечаем,

что критиче­

ская сила вообще и величины, входящие в зти формулы и представ-

§ 2] СТАТИЧЕСКАЯ УСТОЙЧИВОСТЬ АНИЗОТРОПНЫХ ОВОЛОЧЕК

357

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

срединной поверхности панели при потери устойчивости. С уве­ личением чисел та=1, 2, 3,. . . и п = 1, 2, 3,. . . влияние членов,

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

Далее, следует указать, что с увеличением параметра h*

увеличивается расхождение между значениями критических сил, найденными по классической теории (k * = 0 ) и по уточненной тео­

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

Наконец, укажем, что, минимизируя р*пп по т и п, из (2.6), (2.8) и (2.9) получим значения т и п, при которых р*тп получат

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

координат выбрана так, что i?1=co,

R 2= R ,

А ~ 1, 5 =1 . Пусть,

далее,

рассматриваемая оболочка шарнирно

оперта по

торцам

(а = 0,

а=1).

 

 

 

 

а. О б о л о ч к а н а х о д и т с я п о д с о в м е с т н ы м

д е й с т в и е м о с е в о г о

и

р а в н о м е р н о

р а с ­

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

н о г о д а в л е н и я

д.

Для начального безмоментного состоя­

ния имеем

 

 

 

Т\ =

— р, T% = — gR, S° = 0,

(2.11)

тогда для Z* получим

 

— p d?w

 

Z* =

( 2. 12)

 

 

da2

 

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

Д2|,_____ V

^

I

1

d3w

*

 

R

доз Т

Л

дадр2

 

2

+ ч

d3w

1

<?3ц>

(2.13)

 

Л

да^д$

~ Л " Тр* ’

 

 

DA*w + § - (1 - Л *Д )^ + (1 - h * A )l* (p ^

358 КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК [ГЛ. III

где

D :

Eh3

 

Eh2

 

Д _

02

(2.14)

12(1 — м2)»

**= ■ 10 (1 — v2) G' * — 0а2"1"<$2*

Решение системы (2.13)

представим следующим образом:

 

 

и =

.

тка .

иЙ

,

 

 

 

A

cos —

sin

 

 

 

 

V =

г,

. mm

cos

ий

 

(2.15)

 

В

sin —

,

 

 

 

10=

г,

. «пса .

иЗ

,

 

 

 

t

Sin —— Sin

-£■

 

 

где А , В ,

С — некоторые постоянные, I — длина

оболочки, т,

п — целые

числа,

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

средин­

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

оболочки.

 

 

 

Принимая (2.15), удовлетворяем условиям шарнирного опира-

ния по торцам.

 

 

 

 

Подставляя значения и (а , (3),

v (а, В), w ( а,

В)

из (2.15)

в систему

уравнений (2.13), после

некоторых преобразований

получим следующую однородную систему алгебраических урав­

нений относительно неизвестных постоянных А ,

В, С:

А ( X 2

п2) =

С Х (v X 2 п 2) ,

 

 

В (X 2 +

и2) =

— С [(2 +

v ) X 2 +

п21п,

(2.16)

(Х2 + и2)2

(

Х+

рХ2 +

 

9Д п 2 _ ^ _ ^

С \С21 + h* (X2+

И«)

1

2-f « 2)2

E h

 

где

 

 

hi

 

 

 

m n R

 

 

 

E №

 

Cz = 1 2

( 1

— v * ) Д 2 *

h* = 10(1 — v2) RiG " (2.17)

Очевидно, что тривиальное решение системы

(2.16) А = В =

—С= 0 нас не будет интересовать.

 

 

 

Из условия существования нетривиального решения, из треть­ его уравнения однородной системы (2.16), получим соотношение для определения критической комбинации нагрузок р и q:

рХ2-f-q R n i

_ .9. (X2-Н n i)i

,____ X*

(2 18)

E h

~

1 -f h* (Х2 + «2)

(Х2 -f « 2)2

 

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

(h *= 0) формула (2.18)

пере­

пишется следующим образом:

 

 

 

рХ2-|-дД«2

о /*, tt I

«КО I

Х4

 

E h

 

: С2(X2 - f П2)2 ■

(Х2 -f. «2)2

(2.19)

Уравнение (2.18) отличается от соответствующего уравнения классической теории (2.19) лишь наличием члена с множителем/г*,

§ 2] СТАТИ ЧЕСКАЯ УСТОЙЧИВОСТЬ АН И ЗО ТРО П Н Ы Х ОБО ЛО ЧЕК

359

т. е. члена, учитывающего явления поперечных сдвигов. Поэтому дальнейшее исследование уравнения (2.18) ничем не отличается от анализа классического уравнения (2.19).

Не вдаваясь в подробности, укажем, что с помощью уравнения (2.18) нетрудно определить те комбинации чисел т и п , при ко­ торых нагрузки р и q получат свои критические комбинации в слу­

чае учета поперечных сдвигов.

 

б. О б о л о ч к а н а х о д и т с я п о д

д е й с т в и е м

к р у т я щ и х м о м е н т о в М 0, п р и л о ж е н н ы х п о

к о н ц а м о б о л о ч к и .

Такая

нагрузка в оболочке вызывает

следующие внутренние силы безмоментного состояния:

 

 

Z* = 0,

Г» = 0, S° = MJ2*Д2.

 

(2.20)

тогда для Z* получим

MQ d2w

 

 

 

 

 

(2. 21)

 

z*= ui?2 dctdp •

 

В этом случае первые два уравнения системы

(1.8.23)

остаются

неизменными,

а третье уравнение

принимает

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

ДД*ш +

г * , .

d *w

 

 

 

§ .(1

- 5 - ( 1 - ^ д) д гт а - = ° -

<2-22>

 

 

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

приближенно,

а именно:

 

 

 

 

при а = 0,

a= l

W — 0,

 

 

(2.23)

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

Такая форма удовлетворения граничных условий позволяет решить задачу с помощью лишь уравнения (2.22), без учета пер­ вых двух уравнений системы (1.8.23).

Решение уравнения (2.22), удовлетворяющее условиям (2.23), можно представить в виде

w = ^ C { sin

(2.24)

»=i

 

тогда условия (2.23) сводятся к следующим уравнениям:

<?1 + <?2 =

0,

Ci sin (r*x - п | - )+ С2sin ( ц2£- — п |-) =

(2.25)

0.

Условия (2.25) имеют место при любом р и отличных от нуля значениях С{, лишь если

sin (ft, — рх) 2 if = 0,

360 КОЛЕБАНИЯ И УСТОЙЧИВОСТЬ ОБОЛОЧЕК [ГЛ. Ill

т.

е. когда

2mnR

 

 

 

 

 

1*2 1*1

(т =

1,

2, ...).

(2.26)

 

I

 

 

 

 

 

Подставляя (2.24) в (2.22), получим

 

 

 

2

с -( щ г ^ г ,яж ^ + ^ rt [i + £

W ■+" * ) ] -

 

 

М,

 

 

 

 

 

 

s - [ i

 

 

 

" Ю = ° -

< 2 -2 7 >

Отсюда для определения критического значения крутящего мо­

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

 

 

 

 

 

((4 +

п2)2

 

IА

(2.28)

 

«ft- [1 + h* (|4 + пЩ

П(ft2+ » 2)2!•

 

 

где с* и h* определяются по формулам (2.17). Здесь рД£=1, 2) —

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

(2.26).

Кроме того, параметры ft должны быть такими, чтобы выражение

(2.28) для критического значения крутящего момента имело одно

и то же значение при £=1

и £—2.

 

 

 

 

На основании сказанного выше нетрудно каким-либо извест­

ным способом определить критическое

значение крутящего мо­

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

3. Об устойчивости трансверсально изотропной сферической оболочки. Пусть замкнутая сферическая оболочка находится под действием равномерно распределенного по поверхности нормаль­ ного внешнего давления g=const.

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

щие внутренние силы:

 

 

=

Т\ = — Ц -, S0 = 0.

(2.29)

Тогда в силу (1.8.26) и (2.29) из (2.1) получим

 

Z* =

- £ - ( * + 2 ) w,

(2.30)

где А — безразмерный (А = В , B = R sin а) оператор Лапласа,

который в принятых здесь географических координатах имеет вид

Д = ж т й ( 81па1 ) + Ж т й -