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

книги / Устойчивость и колебания трехслойных оболочек

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

перечных

сдвигов <p=ai,1+ 02,2 (где хь х2— ортогональные ко­

ординаты

в срединной плоскости заполнителя; ai,

аг — попереч­

ные сдвиги в плоскости * 0z и yOz соответственно;

z — расстоя­

ние по нормали срединной поверхности заполнителя).

Э. И. Григолюк (10] вариационным методом установил гра­ ничные условия и получил систему разрешающих уравнений ко­ нечных прогибав, произвольно нагретых по толщине и по по­ верхности упругих пологих трехслойных оболочек с легким за­ полнителем, когда несущие слои ортотропны в механическом и термическом смысле, а оси их ортотропии совпадают.

Система четырех уравнений, содержащая до, F,

<*1, аг для изо­

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

к трем (до, F,

cf) и Двум (до, F) нелинейным уравнениям. Здесь

впервые в те­

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

Развивая прежнюю работу (10] Э. И. Григолюк выдвинул [11] общую теорию произвольно нагруженных и нагретых поло­ гих оболочек с ортотропными жесткими несущими слоями и ортотропным жестким заполнителем, сопротивляющимся нагруз­ кам в плоскости, параллельной срединной поверхности. Несжи­ маемый в поперечном направлении заполнитель и несущие слои испытывают конечные прогибы; механические характеристики несущих слоев и заполнителя различны, но оси ортотропии их параллельны. Полученная в общем случае система четырех не­ линейных уравнений (до, F, ai, аг) для изотропных слоев с оди­ наковым коэффициентом Пуассона сводится к системе трех уравнений (до, F,<р). В [10, 11] дается также уравнение прогибов до, описывающее устойчивость пологих оболочек при малых прогибах.

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

Для случая оболочки с изотропными несущими слоями и не­ сжимаемым трансверсально изотропным жестким заполните­ лем 1в постановке [11] систему пяти нелинейных уравнений мож­ но свести [15— 14] к трем дифференциальным уравнениям, со­ держащим функцию прогиба %, силовую функцию F и функцию

71

углов поворота <р, выражая прогиб w и углы поворота нормалЕ к срединной поверхности заполнителя (xi, хг) в виде

W= (1_T уа) х:

(2. 122

где

“i = * i + «\i; a, = va+ w .„

(2. 123,

V2 = d2{)/dx12-{-dfi()/dx22 — оператор Лапласа; h — суммарная тол щина пакета; р — коэффициент поперечного сдвига; у, Ф — коэф­ фициенты, зависящие от толщин и механических характеристи оболочки. При этом одно уравнение оказывается независимы! (1 —fl)h2V 2cp=2pcp (v — коэффициенты Пуассона).

Система дифференциальных уравнений с %, F, <р имеет общи> порядок, равный двенадцати, одно (ф) из уравнений этой cnq темы не связано с остальными и при решении частных задач он» может не приниматься во внимание. Тогда задача сводится Е решению системы двух нелинейных уравнений (%, F ), обща; структура которой весьма напоминает соответствующие урав нения теории конечных прогибов однослойных оболочек Мар герра {27]. Более того, эта система является разрешающей — че рез основные функции %> F выражаются все перемещения и уси лия и, следовательно, граничные условия. Последним достоин сивом многие редуцированные системы, отмеченные выше, ш обладают, и для решения на их основе конкретных задач при­ ходится вновь возвращаться к более громоздкой системе диф ференциальных уравнений.

Г л а в а 3

УСТОЙЧИВОСТЬ И КОЛЕБАНИЯ КРУГОВЫХ ЦИЛИНДРИЧЕСКИ) ОБОЛОЧЕК

1. ЛИНЕАРИЗИРОВАННЫЕ УРАВНЕНИЯ УСТОЙЧИВОСТИ

Относя круговую цилиндрическую оболочку‘радиуса R, дли­ ны I к координатам х\— х, хг= s, измеренным соответственно па образующей и дуге круга нормального поперечного сечения, со­ гласно (2.99)— (2.100) приходим к следующим линеаризирован^

.72

ным уравнениям местной потери устойчивости цилиндрических трехслойных оболочек (рис. 13):

 

v w

=

~

(1 - — v2) z;

(3.1)

 

 

 

R

dx2

\

p

уЛ

 

D I 1 — — у ) v2v2x + —

^ —

[ № n ~ +

 

 

 

)

 

R

dx2

\

и дх?~

 

+

2ЛГ° —

— I-JVS2— 1 (1

8

V2) x = 0 .

(3.2)

~

12 dxds ~

22dsi!

I

J

Л

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

1 ds2

 

 

 

Рис. 13. Круговая цилиндрическая оболочка

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

77= V 2^;’= x = V 2x := V 2V2x — 0 при JC= 0 и х = 1 . (3.3)

Систему уравнений (3.1) — (3.2) с помощью введения разреша­ ющей функции xi

Х = v ^ ,

£ ( 1 - у V2) Хх

(3.4)

сведем к одному уравнению устойчивости

° ( 1 - у v‘ ) ^ ^ + 5 ? й ( 1 - f '»■) *‘ ~

- ( П £ + 2 Л ^ + Л ^ ) ( 1 - f v j v v n - o . (3 .5 ,

7 3

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

X i = V * X i = V * V * X i = V aV2V2X i = V 2V2V *X x = 0 при х — 0 и х = 1 .

(3 .6)

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

2. КРИТИЧЕСКАЯ НАГРУЗКА ПРИ РАВНОМЕРНОМ ОСЕВОМ СЖАТИИ

Рассмотрим [14] местную потерю устойчивости - опертой по торцам круговой трехслойной цилиндрической оболочки, под­ верженной действию равномерно распределенной по контуру сжимающей силы N. Под действие этого усилия в оболочке при докритическом безмоментном состоянии возникнут удельные усилия, равные

= -

N -; ^ = о ; N O 2= о.

(3 .7 )

 

2яЛ

 

Задача сводится к

решению уравнения устойчивости

 

 

 

 

(3 .8)

которое ищем в виде

 

 

 

 

. тпх

ns

(3 .9)

X i= X osln —— cos —

где m — число полуволн

по образующей

цилиндра, л — число

волн по окружности, образующихся

в момент потери устойчи­

вости оболочек; хо— постоянная.

 

 

Вводя (3. 9) в (3. 8),

придем к формуле для критического

осевого усилия

 

 

 

74

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

 

 

 

-^Kpmin=

2я/?(ЛГ°1)крт1п=

РтЫ-

(3.11)

 

Проведенные вычисления показывают, что при k^ O ,

Х .^ 0,4,

р-* >

16 /^практически не зависит от параметра % (см. таблицу).

 

Таблица значений pmin* при

р2= 100,

0 = 0,05,

/г = 1, п —0

 

X

1,00

1,1

1,2

1,3

 

1,4

 

1,5

1,6

.

5,4218

5,4212

5,4013

5,4134

5,4136

 

5,4013

5,4093

Рт\п

 

т

7

7

8

9

 

9

 

10

11

На рис. 14 приведены графики

зависимости

наименьшего

критического

значения

параметра

Р* =

Р*т1в,

определяемого из

формулы

 

 

I

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 3 .12)

от параметра

сдвига k

для значений р2= 64,

1000, 5000,

10000,

15000,30000 при 0 = 0, 1.

 

 

 

 

 

 

С достаточно хорошим приближением эту зависимость мож­ но получить, считая т21№+п21п2 и т 2Д 2 непрерывными аргумен­ тами и проводя минимизацию выражения для р*.

Полагая

 

 

перепишем формулу (3. 12) в виде

 

Р*

1 -j- ftkmi

(3.13)

1 Н- km\

 

 

(^ -< m 1< o °, m1< m 8< ooj .

 

Легко видеть, что минимум выражения (3.13)

достигается

на границе изменения переменных, определяемой

равенством

 

гп{—т2,

(3.14)

75

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

 

 

Р*

1 -}- Ыгтпл

.

и.2

 

 

 

 

(3.15)

 

 

---------- Lml +

— •

 

 

 

 

 

 

 

1 4 * ktn.\

 

Ш\

 

 

 

 

 

Рассмотрим два возможных случая.

 

 

 

 

 

 

1. Пусть

/гт! — 1,

где

т ,

соответствует

/>*1п.

Здесь

член

1 и им можно

пренебречь,

вследствие

чего формула

 

 

 

(3. 15) запишется так:

 

 

 

 

 

 

 

р *_ ^

ml

I

!*2

_

 

(3.16)

 

 

 

 

 

1 + k mi

mi

 

 

 

 

 

 

Из условия d p*/dm l= 0

имеем

 

 

 

 

m

t*

;

/СП= М 2 - ^ ) .

 

(3.17)

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 14. Зависимость критического значения

 

 

 

продольно

сжимающего

напряжения р* =

 

 

 

= р min круговой цилиндрической

оболочки

 

 

 

от

параметра

сдвига

k

при 0= 0,1

для ря­

 

 

 

 

 

 

да значений

р,2

 

 

Формула

(3.17) дает

достаточно

точные

значения

P*mln>

если выполняется условие

 

 

 

 

 

 

 

 

 

 

 

 

1 -р А > & .

 

 

 

 

 

 

(3.18)

В частности, при k = 0

получим известную формулу для кри­

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

2 . Когда •bkm.i~\, где rrii соответствует

Ф°Р~

мула для р* приобретает вид

 

 

 

 

1 + ИЙ/И]

 

(3.

19)

 

mi

 

 

 

 

 

Определяя корень уравнения d p*ldm i= 0 ,

найдем

 

т г

_м_

 

(3 .20)

p*m\n— v-V*

k

 

 

 

 

Уравнение (3. 19) с достаточной для практических расчетов, точностью может быть использовано при ^> 0,2.

Согласно (3.20) и (3. 11) при &>0,2 имеем

11/кр min"

~ ~

\

f — ~~----- b Oha—

(3.21)

 

2R

V

3 (1 — v2) 1 , 3 eT

 

76

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

£iAi2

£ 2Л22

кр щ1п'

2R / 3 ( 1 — V22)

2R ^ 3 (1 — Vl2)

 

 

•‘ч

и

т

т

(3.22)

 

 

 

 

 

 

 

Пренебрегая

изгибной жесткостью несущих слоев, найдем

 

(^И )крга.л~ -0Л8[ 1 + 4 _ ( _!4 Г ^ ) 2] •

(3-23>

Отбрасывая

в скобке

второй

член,

 

придем к общеизвестной формуле

 

 

 

 

(-^ц)кртт~

Oh3,

 

(3.24)

 

полученной в предположении, что средин­

 

ные поверхности несущих слоев прикреп­

 

лены к поверхностям заполнителя.

 

 

 

На рис. 15 представлены зависимости

 

P* = Pmia

°т десятичного логарифма па­

 

раметра

ц2 при

k = 0 ,

0,01,

0,05, 0,1,

0,2,

 

0,5, 1000,

когда 0= 0,05.3*

 

 

 

 

Рис. 15. Зависимость критического значения про­ дольного сжимающего напряжения р* круговой цилиндрической оболочки от десятичного лога­ рифма параметра |х2 при 0= 0,05 и ряда значений параметра сдвига k

3.НЕРАВНОМЕРНОЕ ОСЕВОЕ СЖАТИЕ

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

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

№ = - ( — + - ^ - c o s — ) ,

(3.25)

гдеМдзг— изгибающий момент, равный

Mw3r= N e —NRs

(3.26)

— эксцентриситет силы N ).

77

Решение уравнения устойчивости

° ( 1 -

т

i

7w

,

’z‘ + U

s r (

' -

f Л

«*+

 

+ i ^ ( 1+ 2 i c° s i ) S ( 1 - T vV ' ' * ^ 0

(3.27)

 

ищем в форме

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

®=sin

 

^

a ne ins'R ;

 

 

 

 

 

X i — sin

^

a ~ l e inslR ,

 

 

(3.28)

где

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h W . 8_ 12/?2(1 — y2) .

 

(3.29)

 

 

2Dix3

 

P/?2 ’ ^ ~

Я4Л20

 

 

 

 

 

 

 

 

 

 

i

l

 

e

 

 

 

 

 

 

 

 

к — — : s =

 

 

 

 

и подставляя

(3.28)

в уравнение (3 .27),

найдем

связь

между

коэффициентами ряда (3. 28)

 

 

 

 

 

 

 

a n+iP*&— (Фя— P *)a n+ a n- iP * s = 0

при п ф

О

(3.30)

 

 

ф0— ^1 -]-2е — j

р *

При П =

0.

 

(3.31)

Здесь

 

 

 

 

 

 

 

 

 

 

 

,

 

„ / /к2

п2 \

Infl

л2 \

 

/и2

 

 

t . ~

 

 

 

 

 

 

I

г»

X2

 

(3. 32)

,

( т2

 

n2\

 

2

//и2

л2 ’

 

1

+ k l ----+ ----

 

X2

 

 

 

 

 

 

 

\ Х2 ^

я2}

 

 

 

 

 

 

Преобразуя (3. 30)

 

 

 

 

 

 

 

 

 

 

 

«л- 1

 

Фл — Р*— ЕР* «я *-1

 

 

 

 

получим, что для достаточно больших п

 

 

 

 

 

 

 

 

«л

_

Р*Е

 

 

 

(3.33)

 

 

 

 

«л- 1

 

фл— Р*

 

 

 

 

 

 

 

 

 

 

 

 

78

Возвращаясь к Оо, будем иметь бесконечную

цепную дробь

£ L

 

р2

 

 

(3.34)

«о

 

р*е

 

 

 

<h — Р*

 

л

 

 

 

 

 

р*£

 

 

 

 

h —p * — Фз—Р*

 

 

 

Теперь на основании (3.31)

 

имеет уравнение I B

виде

бесконеч­

ной цепной дроби

 

 

 

 

 

Фо Р*

 

р*£2

 

 

(3.35)

2р*

 

р*е?

 

 

 

 

 

 

h — р *

 

р*г

 

 

 

 

Ф2— Р*

е2

 

 

 

 

Фз~Р*

 

 

 

 

 

 

 

из которого можно определить как р * при заданном

е, так и в

при заданном р*.

 

 

 

 

 

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

NR

Фо

2лЮ Р

1 + 2е ’

т. е. критическая сила, приложенная вне центра, будет равна или больше критической центральной силы, деленной на вели­ чину 1 + 2 е, где е — эксцентриситет приложения силы.

4.КРИТИЧЕСКАЯ НАГРУЗКА ПРИ ВНЕШНЕМ РАВНОМЕРНОМ

ПОПЕРЕЧНОМ ДАВЛЕНИИ

Исследуем потерю устойчивости опертой по торцам тонкой упругой круговой пологой трехслойной цилиндрической оболоч­ ки, подверженной действию внешнего равномерного поперечно­ го давления q [14, 15]. Под действием этой нагрузки в оболочке при докритическом безмоментном состоянии .возникнут удель­ ные усилия, равные

№и = 0 ;

N ^ - q R

; №п = 0 .

(3.36)

Задача сводится к решению уравнения устойчивости

 

° ( '

 

£ 0 — 1г v j b +

 

 

 

 

(3.37)

Отыскиваем решение уравнения (3. 37) в виде

 

Ъ. =

. тех

ns

00\

Хо sin —

COS — .

(3.38)

Здесь п — число волн в окружном направлении; хо— постоянная.

79

Подставляя

(3. 38)

в

(3.37),

приходим к уравнению

 

 

 

1 -м »

 

 

 

 

п2_у

 

 

 

 

 

 

д* =

 

 

 

 

Я 2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rfi

Х«

* ( a L

+

±

) 2

*

 

1 - м

 

 

 

 

 

 

 

Я2

 

я \я2

 

X2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.39)

где

 

д*

12^

3(1 - У 2)

.

а.

12^2(1 у2) .

 

 

 

 

 

 

 

 

(3.40)

 

 

4

ЯЛЗДя2

 

 

Я40Л2

 

 

 

 

 

 

 

Я2А2 _

I

 

 

 

 

 

 

 

 

 

 

k — -

 

 

 

 

 

 

 

 

 

 

 

 

т

R

 

 

 

 

 

 

 

 

 

На рис. 16 приведены графики

 

 

 

 

 

 

минимального значения критическо­

 

 

 

 

 

 

го давления

q

— 9,min

в зависимо­

 

 

 

 

 

 

сти от коэффициента %для значений

 

 

 

 

 

 

параметров при ц2= 15000 и ■0=0,05,

 

 

 

 

 

 

когда 6 = 0 ; 0,04; 0,08; 0,2; 0,8; 2; 1000.

 

 

 

 

 

 

Рис. 16. Зависимость критического значения

 

 

 

 

 

 

внешнего

равномерного поперечного давле­

 

 

 

 

 

 

ния

 

круговой

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

 

 

 

 

 

 

оболочки

от

отношения X=l!R (/ — длина,

0 1 2

3 9 5

6

7 8

9

7 .

R — радиус)

при

р2= Г5 000,

6= 0,05

k

для

ряда

значений параметра

сдвига

 

 

 

 

 

 

 

5. КРИТИЧЕСКАЯ НАГРУЗКА ПРИ ВНЕШНЕМ РАВНОМЕРНОМ ВСЕСТОРОННЕМ ДАВЛЕНИИ

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

N\x= - \ q R \ №22= - q R ;

0.

(3.41)

Задача сводится к решению уравнения устойчивости

 

° (1 — т - v ‘ ) v ’ v W * - + - £ £

( 1—f ^ * +

 

+**(т2 £dx2+1£ds2)(■ - f v ) V

V x , = . 0,

(3. 42)

которое имеет решение в форме

 

 

 

пх

ns

 

(3.43)

Хх=Хо sin — cos —

 

80