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

книги / Прикладные задачи устойчивости многослойных композитных оболочек

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

3.2. Многослойные цилиндрические оболочки (классическая модель) 51

3 .2 .1 . О севое сж атие. В этом случае, полагая S = Тг = 0, найдём

Критическое усилие Т*р находится отсюда минимизацией по обоб­ щённому параметру волнообразования А(1 + ф2):

ТГР = Т10 = |\ / Ш З ; В = В{ 1 - г /2);

(ЗЛО)

В случае однородное изотропных оболочек, полагая

Eh3

В = Eh-, D = 12(1 - 1/2)

(Е, v, h — соответственно модуль упругости, коэффициент Пуассо­ на и толщина оболочки), из зависимостей (ЗЛО) получим известные формулы для расчёта критических напряжений и параметров волнооб­ разования:

тк р

Eh

0 ,6 0 5 ^ ;

л/3( 1 - г/2) R

 

i t

Обратим внимание на то, что при расчёте на устойчивость изотроп­ ных цилиндрических оболочек при осевом сжатии форма волнообра­ зования не определена однозначно: параметры волнообразования Акр и пкр могут принимать произвольные значения, лишь бы между ними выполнялись соотношения (ЗЛО) или (3.11). Таким образом, одной критической силе соответствует множество вариантов форм волнооб­ разования. Это многообразие форм [36] является одной из причин большого различия между теоретическими и экспериментальными зна­ чениями критических напряжений [4, 10, 16, 45, 64, 81]. В связи с этим различием необходимо вводить поправочный коэффициент, умножая на него значения критических усилий и напряжений в формулах (ЗЛО), (3.11). Поправочный коэффициент куст можно принимать в соответ­ ствии с рекомендациями [8]:

R /h

250

500

750

1000

1500

h

0,3

0,23

0,20

0,165

0,15

fVycr

52 Гл. 3. Общая устойчивость изотропных цилиндрических оболочек

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

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

4 р Т ^ t ^AIB R2

 

V 2 V В Д2 < R < 2 V D

 

( l, 2 y / h /R < | < 3 y/R/h'j .

(3.12)

Если оболочка короткая, т. е. выполняется условие

 

 

D

 

R < V 2

; < 1,2 y/h/R^j ,

 

V Ш ¥

 

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

Ткр 2D

Для однородных по толщине оболочек отсюда получаем известные [8] зависимости:

7T2E h 2

h2

i

1.2у / Щ .

12(1 -1/2)^ -°<9° £

£2 ’

д <

В случае длинных оболочек, когда

 

 

 

| > \ \ j BR?/D (

| > 3

у / Щ

у

расчёт их на устойчивость следует проводить по формулам для сжатых стержней [81]

т кР = 2B R 2 2£2 '

3 .2 .2 . Внеш нее давление. Для исследования критических пара­ метров при действии внешнего давления р примем известное положе­ ние, что при потере устойчивости образуется одна полуволна в осевом направлении:

токр = 1 ;

3.2. Многослойные цилиндрические оболочки (классическая модель) 53

(R, I — соответственно радиус и длина оболочки).

Полагая в зависимости (3.9) Т\ = S = 0, получим соотношение для отыскания критических параметров пкр и Т*р при потере устойчивости цилиндрической изотропной пологой оболочки от действия внешнего

давления:

 

 

.0(1 + А[/п2) п2

(3.13)

Т2 =

+

R 2

 

n6(1 + А2/п 2)

T2 = p R ;

Х\

R

~Т'

 

 

Критическое значение кольцевого усилия Т2 получается минимиза­ цией по дискретному параметру п. Если параметр п достаточно велик, то его можно считать непрерывно изменяющимся. Проведя в этом случае минимизацию при условии, что величина А2/п 2, содержащаяся в скобках, мала, найдём:

ллкр

1.75тг

V B D 3 ;

п:

4/ЗОЙ2

Акр — А).

(3.14)

2

£Д>/2

 

кр

~ с Г

 

 

Если положить

 

 

 

 

 

 

 

В = Eh; D

Eh3

 

(3.15)

 

 

12(1 - и 2) ’

 

 

 

 

 

 

 

то из соотношений (3.14) получим известные формулы П.Ф. Папковича:

ТКР =

0,85

Eh5/2

=

R

R

 

2 ^ 3 / 4

IRI/2 ’

7,7

 

(!

- к р

J

\ h ’

- г/2)

 

 

 

Р кр

0,85

Eh 5/ 2

(3.16)

 

(1 _ ^2)3/4

£R3/ 2

 

 

 

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

 

_ Д(п2 - 1 )

ВХ\

2

R 2

п4(п2 —1)'

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

Как показывают соотношения (3.14), при потере устойчивости ци­ линдрической оболочки от внешнего давления параметры Акр и пкр волнообразования однозначно определены. Это позволяет предпола­ гать, что теоретические и экспериментальные значения критических усилий должны хорошо согласоваться. Опыт показывает, что различие

54 Гл. 3. Общая устойчивость изотропных цилиндрических оболочек

теоретических и экспериментальных результатов для многослойных оболочек находится в пределах 10 ч- 15%.

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

6 I / D / B R 2 < ^ < 2 {/BR?/D

( з у / Щ < ^ < 3 ^ / Щ \ . (3.17)

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

| < б f o r m e

( ^ < 3 у / щ у

критические усилия в них соответствуют критическим усилиям шар­ нирно опёртой прямоугольной пластины [8]:

7T2D

т,7

Длинные оболочки, для которых выполняются условия

£ > 2 P & / D ( | > 3 у / щ у

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

ЗГ)

п о

т кр д_

± 2 — д 2 ’

кр —

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

ется одна полуволна в осевом направлении ^Акр = А) = С учётом этого из соотношений (3.9) можно получить для пологих оболочек:

D n 3(1 +А 2/гг2)2

вх3

Ai =

R

25 =

+

(3.18)

A7R 2

n5(1 + А2/ п 2)

 

~Г'

Это соотношение соответствует решению (1.3), которое в рассмат­ риваемом случае не удовлетворяет граничным условиям шарнирного опирания и даёт заниженный результат. Используя результаты [24], где граничные условия удовлетворены точно, можно найти следующее выражение для расчёта критических параметров (п2 3> 1, п2 з> А[):

5 кр =

3,3

 

D5 ;

B R 2

 

£ 1 /2 Д З /4

у[вг

П к р = 2 -7 Л 1

 

 

0,7

Д/г5/4

^кр

R

[R

(3.19)

(1 -|/2 )5 /8

Р Р

;

15,8 УГ

 

 

£ y h

 

3.2. Многослойные цилиндрические оболочки (классическая модель) 55

Зависимости (3.19) относятся к шарнирно опёртым оболочкам. Если

осуществляется жёсткая заделка, то коэффициент

3,3 в формуле для

5 кр следует заменить на 3,96, а коэффициент 0,7 в

формуле для ткр —

на 0,84 [24, 81].

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

R/h

250

500

750

1000

1500

h

0,80

0,70

0,63

0,58

0,5

fVycT

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

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

 

 

 

6 {/D /BR2 < д < 5 y\ /fBmRP2/D

 

< |

< 7 y /R ih } . (3.20)

Для коротких оболочек, когда

 

 

| < 6 { / D / B R 2

- < 3yfh fR ) ,

следуя [8], найдём

 

 

 

 

Дкр = ks

TT2D

 

2Eh2

~ё~

?кр ks 12(1 - v 2)P '

Здесь ks = 5,34 для

шарнирно

закреплённых

оболочек, ks = 8,97

для жёстко защемлённых.

Если длина оболочки достаточно велика, то в кольцевом направле­ нии образуются две волны (пкр = 2), а величина критического усилия не зависит от длины оболочки и граничных условий. Полагая пкр = 2,

с помощью (3.9) найдём:

3D вх3

AR2 + 48 '

Минимизируя это соотношение по параметру А, придём к формуле,

аналогичной известной формуле Шверина [14, 66]:

 

2

B D 3

48D

(3.21)

W

А:кр

нкр —2.

дз/2

Т Ш ’

 

В частном случае однородных по

толщине конструкций,

полагая

в <321)

 

 

 

В = Eh; D =

 

 

12 (1 - г/2) ’

56 Гл. 3. Общая устойчивость изотропных цилиндрических оболочек

получим формулу Шверина [79]

Е( b _ Y 12

3\/2 (1 - j/2 )3/4 \ R J

3 .2 .4 . О зависимостях м еж ду критическими усилиями в и зо ­ тропны х цилиндрических оболочках. Выпишем полученные выше значения критических усилий для многослойных цилиндрических обо­ лочек при осевом сжатии, внешнем боковом давлении и кручении:

Ткр —— л/ BD ■

Т кр —

\ / B D 3 •

S Kp —

^

\/"R3Г)5

R v v u ,

1 2 -

l R l / 2 w u u ,

Ь -

i l / 2 R 3 / A

У Е и

Нетрудно видеть, что между этими критическими нагрузками су­ ществует зависимость

sKf>=^/т,крт2кр

(3.22)

Точно так же, между

известными критическими напряжениями

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

 

1

Eh

0,85

Eh? / 2

< Р =

~пR ’

(l _ г,2)3/4 I R 1/2’

3(1 - v2)

ткр _

0,7

Eh 5/4

 

 

Т - V 2)5/8 / 2 ДЗ/4 '

 

— выполняется соотношение

 

 

ткр

^.кр ^.кр

(3.23)

 

 

 

Заметим, что соотношения (3.22) и (3.23) получены для случая шарнирного опирания цилиндрических оболочек.

3.3. Устойчивость трёхслойных цилиндрических оболочек с лёгким заполнителем

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

3.3. Трёхслойные цилиндрические оболочки с лёгким заполнителем

57

А2Т, + 2Аn S + п2Т2 = [DHC+ П (А, п) D*]

2 , „ 2\

ВХА

 

 

 

R 2

(А2 +гг2)'

 

 

 

 

 

. 2 /1

/2,\ ^

В

21

^ + 2 ^ + ^ 2Т2 = [Я“ + П (А, ф) D *] — 1

j

~

 

 

 

л

А2 (1 + фА)

П(А,п) = О (А, ф) =

1 + (х>(А2 + гг2)

1 + <^гА2(1 + ф2) ’

 

Ф

п

_ D*

 

 

(3.24)

А ’

 

 

 

 

 

 

 

 

Разрешающие соотношения (3.5) на основе полубезмоментной мо­ дели непологих трёхслойных оболочек в случае изотропных цилиндри­

ческих оболочек примут вид:

 

 

 

 

 

1)

S + T2 (n2 - 1) =

 

 

 

 

 

 

D*

 

("2 - о 2

 

В А4

 

DHCА_________—-------------

+

ТГ

(3.25)

 

1 + ix>A2 + ш (п2 —1)

Д2

3 .3 .1 . О севое

сж атие. В этом

случае,

полагая

в

зависимостях

(3.24) Т2 = S = 0, найдём

 

 

 

 

 

 

D*

А2 (1 + ф2)

 

В

 

Г Г \-

1 + о;А2 (1 + ф2)

R 2

+

 

 

 

 

А2 (1 + ф2)

 

Можно показать, что минимум величины Т\ достигается при ф = О, т. е. трёхслойные изотропные оболочки с лёгким заполнителем теряют устойчивость по осесимметричной форме (пкр = 0). В результате полу­ чим следующее выражение для расчёта критических усилий:

D*

А2

В

 

Т, = DHC+ 1 + о;А2

Д2 +

А2'

(3.26)

Проводя минимизацию этого выражения по параметру волнообра­ зования А, получим выражение для расчёта его критического значения:

_

Ао

 

 

Ai

 

 

 

BR2

Ак р —

1+ —

 

ГГ

+ ■

 

ГГ

|

D*

 

 

1 т

 

(l+wAjp)

 

 

 

 

 

 

( 1 +

 

 

( l +wA^)

 

 

 

 

 

 

 

 

 

 

 

л4

B R 2

л4

B R 2

ВД2

х4

 

(3.27)

 

А°

DHC

А*

D*

D

Аж'

 

 

 

 

58 Гл. 3. Общая устойчивость изотропных цилиндрических оболочек

Из соотношений (3.27) следует, что при общей потере устойчивости от осевого сжатия трёхслойная изотропная оболочка образует меньше осесимметричных волн, чем её раздельно работающие несущие слои (Акр = Ао ) и больше, чем соответствующая трёхслойная оболочка с абсолютно жёстким на сдвиги заполнителем (Акр = Аж). Таким об­ разом, для критического параметра Акр в трёхслойной цилиндрической изотропной оболочке выполняется неравенство

Аж ^ Акп ^ Ап

(3.28)

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

Если заполнитель имеет малую сдвиговую жёсткость, то выполня­ ется условие wA2 1. Используя это условие, найдём из зависимостей

(3.27):

2К 2К

К Г = К (1 - ь 2)'.

■10 ■10

R

Здесь Т \о — критическое усилие в трёхслойной оболочке с абсолютно жёстким на сдвиги заполнителем; К — жёсткость трёхслойного пакета

на поперечные сдвиги. Пренебрегая в соотношении (3.26) единицей 2

по сравнению с величиной (wA2p) , приведём его к виду

DHCА2

,

в ( ' ~ ъ 2)

+ К.

(3.30)

Т\

+

А2

R 2

 

 

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

Т,кр = 7|Н0С\ / 1 - Ъ 2 + К; А^р = AQ(1 —62);

пкр = 0.

(3.31)

Здесь Т[ос = у/ B D nc — критическое усилие

при осевом

сжатии,

соответствующее оболочке с раздельно сопротивляющимися несущими слоями. Зависимости (3.31) можно получить также непосредственной минимизацией по параметру А выражения (3.30).

Как видно, полученные расчётные формулы (3.31) пригодны, если выполняется условие

2К

(3.32)

< 1

чо

3.3. Трёхслойные цилиндрические оболочки с лёгким заполнителем 59

В случае выполнения равенства 2К = наступает так называемая сдвиговая форма потери устойчивости, и расчётные формулы принима­ ют вид:

Т*р = К- Акр = 0.

(3.33)

Если заполнитель имеет достаточно высокую жёсткость на сдвиги, то шА^р ~ 1. В этом случае можно из зависимостей (3.27) получить

\2 ™

Л- = Т ^ Л Г “<

= #

" - - f t

<334>

а соотношение (3.26) представить в виде

 

 

Тх

D

А2

В

(3.35)

1 +wA2 Д2 + Д2-

 

 

Зависимость (3.35) по существу определяет критические параметры при осевом сжатии трёхслойной оболочки, сопротивляющейся в со­ ответствии с моделью прямолинейного элемента. Подставляя в выра­ жение (3.35) значение А = Акр, получим формулу для расчёта кри­ тических усилий в изотропных трёхслойных оболочках с жёстким на сдвиги заполнителем:

т к Р =

Т;

 

ФУК

 

 

40

1 _ f i o

(2К > Т\ю)

(3.36)

 

10

 

1 + ±1°

 

АК

 

 

 

 

 

 

 

 

2К

 

 

 

 

Как следует из формулы (3.36), она пригодна, если выполняется

условие

 

 

 

К

>

\ т

(3.37)

Формулу (3.36) можно получить также непосредственной миними­ зацией зависимости (3.35) по параметру А. Сравнение критических усилий, полученных по формулам (3.31), (3.36), с результатами мини­ мизации по параметру А общего выражения (3.26) показало их прак­ тическое совпадение: погрешность не превосходит одного процента.

Из формул (3.31), (3.36) следует, что критические усилия при сжатии трёхслойных цилиндрических оболочек определяются тремя величинами:

TJQC — критическое усилие, соответствующее оболочке с раздельно работающими несущими слоями (жесткостная характеристика несу­ щих слоёв);

TjQ — критическое усилие при осевом сжатии трёхслойной обо­ лочки с абсолютно жёстким на поперечные сдвиги заполнителем (жесткостная характеристика всего трёхслойного пакета);

К — критическое усилие трёхслойной оболочки при так называемой сдвиговой форме потери устойчивости (жёсткость трёхслойного пакета на поперечные сдвиги).

60 Гл. 3. Общая устойчивость изотропных цилиндрических оболочек

Эти три величины естественно назвать обобщёнными жесткостями трёхслойных цилиндрических оболочек при осевом сжатии. Обобщён­ ные жёсткости определяют основные закономерности потери устойчи­ вости. Так, из соотношений (3.31), (3.36) следует, что для трёхслой­ ных оболочек характерны две области влияния поперечных сдвигов. В области 2К < Т * (слабый на сдвиги заполнитель) критическое усилие практически прямо пропорционально сдвиговой жёсткости К заполнителя (Т*р и К) и слабо зависит от жёсткости несущих слоёв. В этой области критическое усилие рационально увеличивать за счёт повышения жёсткости К оболочки на поперечные сдвиги. В области достаточно жёстких на сдвиги заполнителей (2К > Т^) в соответствии

ссоотношениями (3.36) критическое усилие определяется величиной

асдвиги играют поправочную роль. Так, увеличение жёсткости запол­ нителя на сдвиги в несколько раз приводит к увеличению критического усилия лишь на несколько процентов. Отсюда следует, что в области

2К > Т[* повышения критического усилия рационально добиваться за счёт повышения жёсткости В на растяжение-сжатие или увеличения расстояния Н между нейтральными поверхностями несущих слоёв.

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

оболочки:

 

К > Tffi.

(3.38)

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

Обобщённые жёсткости и соотношение между ними определяют области применения различных расчётных моделей. Известно, что если пренебречь собственной изгибной жёсткостью несущих слоёв (Т^ = 0), то придём от модели ломаной линии к модели прямолинейно­ го элемента. Из соотношения (3.31) следует, что переход к модели пря­ молинейного элемента оправдан, если выполняется условие К Tf^. Поскольку это условие, как только что было показано, является од­ новременно и условием (3.38) эффективности трёхслойных оболочек, то для трёхслойных конструкций во многих случаях при расчётах на устойчивость можно использовать модель прямолинейного элемента. При этом погрешность Дпр применения такой модели составляет

Соседние файлы в папке книги