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

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

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

7.3. Расчёт критического внешнего давления

191

Обычно при потере устойчивости от внешнего давления выполня­ ется условие гг2 1. В этом случае из соотношения (7.15) после ми­ нимизации по параметру п получаются следующие расчётною формулв1 для определения критического колвцевого усилия:

fjnKp _ rjnnc

}не)

(7.16)

К 2 - Т Т = 7 ^ Г

\ / в Л О Г 3

ip.

Коэффициент у в соответствии с [8] вычисляется по формуле:

7 = 3 —2,1-г\

О < — < 0,8.

(7.17)

 

Г2

Г2

 

При малой конусности

> 0,8^

полагают 7 =

1 .

 

у г2

 

 

Полученные расчётною формулв1 применимы, если выполняется

условие

 

 

1,75тг

А/=

(7.18)

К 2 < 0,5 Т2Ж; Т2Ж = 7

В ф \ .

крЛ-х.

Если заполнителе достаточно жёсткий на поперечною сдвиги: К 2 > > 0,5 Тж, то можно использоватв модель прямолинейного элемента (типа сдвиговой модели Тимошенко). Тогда разрешающее соотношение

имеет вид

п

В 2

В\ Af

 

Т2 =

,2 Д2

гР2m ax

(7.19)

1

+ и>рп2

Пи

После минимизации этого выражения по параметру п можно по­ лучить следующую формулу для определения критического кольцевого усилия [72]:

грКр _ грж '

Q

Т * \

 

' -

l e

t -

, 7а д

Зная значения критических усилий, можно найти критическое внешнее давление. Для оболочек со слабым на поперечные сдвиги заполнителем 2 < 0,5 Тж):

 

рКР=рпс + _ К ^ .

рнс = 7 1

^ .

4/^1 { D r ) 3

(721)

 

 

 

Рmax

I p ^ax

 

 

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

( к 2 > \

Т» критическое внешнее давление определяется формулой

ркр = р»

, _ л т ?

 

1,75л

(7.22)

 

 

B XDI

 

 

 

16

к2

^

£Ртах

 

 

 

 

 

 

е 3/2

 

 

При

действии

внешнего

давления

в

расчётные

теоретиче­

ские формулы

вводится

поправочный

эмпирический коэффициент

куст = 0,8 -7 0,9

[23,

27].

 

 

 

 

 

192Гл. 7. Устойчивость многослойных композитных конических оболочек

7.4.Устойчивость при действии крутящего момента

исдвигающих усилий

При исследовании устойчивости от действия крутящего момента Мкр или сдвиговое усилий S ограничимся слабоконическими пологи­ ми трёхслойными оболочками. Тогда можно восполвзоваться зависимо­ стями для цилиндрических оболочек. Особенноствю потери устойчиво­ сти подобных оболочек является то, что в случае слабв1х на поперечные сдвиги заполнителей при потере устойчивости образуется много волн как в кольцевом, так и в осевом направлениях [66]. В то же время до­ статочно жёсткие на поперечные сдвиги оболочки теряют устойчивость по классической схеме, образуя одну полуволну в осевом направлении. Поэтому необходимо рассмотреть отдельно каждый из этих типов потери устойчивости.

В случае трёхслойных оболочек со слабым на поперечные сдви­ ги заполнителем по аналогии с цилиндрическими оболочками можно получить следующие зависимости для расчёта критических усилий сдвига [66]:

S KP= y / ( K i + T - p ) К >; г кр k opTrI Z ^ j 1 - ^

ГГ = - ^

J B 2 [)»>■-,

т г = —

A/ S A ;

ртах = ^ т т . (7.23)

Р т а х

V

Р т а х

V

COS 0

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

[ 66] :

 

rp

 

 

 

 

 

,S'Kp

1 - 0 ,3 4 ^ -

 

 

 

л 2

 

S УК

3,3

1,75тг

(7.24)

 

 

 

^ Р т а х

 

Расчётные зависимости (7.23) для оболочек с малыми жесткостями

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

 

лГ кЖ 2 < 0,33 5 Ж.

(7.25)

Расчётные зависимости (7.24) для оболочек с большой жёсткостью

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

 

/ ^ г ^ > 0 , 3

3 5 ж.

(7.26)

В пограничной области

~

0,33 £ ж) следует провести рас­

чёты по формулам (7.23), (7.24) и выбрать минимальный результат.

7.4. Действие крутящего момента и сдвигающих усилий

193

Критические

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

момента М кр и эквивалентной

перерезывающей силы QKp определяются зависимостями:

 

y/Ki K2 < 0,33 ^

:

 

 

 

 

 

 

 

№> = 2np2mi^ /( K i + T-P) K2;

 

 

 

QKp = npminV(Ki + T r ) K2 ;

 

(7.27)

y/Ki K2 > 0,33 ^

:

 

 

 

 

 

 

 

П^Ж.\

 

 

6,6п pmmn

VbD 5 ;

 

М кр = М ж

1 - 0,34 Т2

\ ; М ж=

 

 

K2) ;

 

 

71/2

 

 

QKp

T2?\

;

Q

=

3,3п pm/4

^ в D .

(7.28)

K2

;

 

 

 

71/2

 

 

Соотношения (7.27), (7.28) получены с помощью формул (7.23), (7.24) и с учетом зависимостей

М = 2nR2S ; Q = nRS.

7 С.Н. Сухинин

Г л а в а 8

УСТОЙЧИВОСТЬ многослойных

СФЕРИЧЕСКИХ ОБОЛОЧЕК ПРИ ДЕЙСТВИИ ВНЕШНЕГО ДАВЛЕНИЯ

8.1. Разрешающие соотношения

Рассмотрим устойчивость сферического сегмента с углом полураствора $. Значению д = 7г/2 соответствует полусфера, значению = д — полная сфера. Разрешающие соотношения на основе модели лома­ ной линии применительно к исследованию устойчивости трёхслойных и многослойных сферических оболочек получаются из зависимостей (1.1)-т-(1.5), если положить в них вместо радиусов кривизн R \, R 2 радиус сферической оболочки R = R\ = R 2. При S = 0 получаются разрешающие соотношения для анализа устойчивости сферического сегмента в случае действия наружного давления [71]:

 

T = ( D f +

nz?f)

 

 

A2

B 2

1 + Ф2

 

 

R2

 

1+ ф2 + F2 (ф)

A2

 

 

 

 

 

 

F\

(ф) = 1 + а-1Ф2 + Ф\ФА\

 

F2 (ф) = 1 + сх2ф2 + /32фА\

 

 

Ш7Г

 

 

 

pR

£ 2

 

 

В2

ф = п / A; A

T = T\ T2

~ Y ; A

D \ ’ ih P

Вф

 

 

^ A 2 , 0

\

B 2

0

 

 

/0 N

 

ai = ---- b 2 V2

a2 =—-----2 г/2-

 

 

(8.1)

 

 

 

 

 

B \2

 

 

 

 

Здесь D™c

суммарная изгибная

жёсткость

несущих

слоёв,

D\ =

= 77“с + D* — изгибная жёсткость трёхслойного пакета; D"c + П D* —

эффективная

изгибная

жёсткость трёхслойного пакета; В 2 = В2х

х (1 —^1^2);

В2, i/\, щ

жёсткость пакета с лёгким заполнителем

в широтном направлении и соответствующие коэффициенты попереч­ ной деформации; р — наружное давление; од, а2, фи Ф2 параметры анизотропии трёхслойной конструкции. Величина П(А,ф), определяю­ щая влияние поперечных сдвигов (0 < Q ф 1), имеет вид (1.5).

Критическое значение усилия Ткр и соответствующие значения Акр, фкр параметров волнообразования получаются из (8.1) минимизацией по А и ф.

Из разрешающего уравнения (8.1), построенного на основе общей модели ломаной линии, как частные случаи получаются зависимости,

8.2. Устойчивость изотропных оболочек (классическая модель)

195

соответствующие моделям более низкого уровня. Так, если Q.D*

3> D™, то применима модель

прямой линии,

и

соотношение

(8.1)

примет вид:

 

 

\+ф2

 

А

А2ВД)

_

 

1 + Л2 (ип +си2Ф2) (1 + Ф2) R2

Щ{ф)'

 

А

А

 

 

,0 г,\

W1 -

Ш2

 

 

(8-2)

Здесь А> А — жёсткости ортотропной оболочки на изгибы в двух

направлениях; шь

^

-

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

перечные

сдвиги;

К\,

А

жёсткости

пакета на

поперечные сдви­

ги [41].

 

 

 

 

 

 

В том

случае,

когда

жёсткость на

поперечные

сдвиги заполни­

теля велика, т. е. податливости сщ, малы и выполняется условие А2(сщ Т и>2'ф2) -С 1, справедлива классическая модель оболочек с неиз­ менной нормалью. Тогда из соотношений (8.2) получается зависимость для классических сферических оболочек:

= А

х*Щ)

1

R 2

1 + А

Щ{ф)

8.2. Устойчивость изотропных оболочек (классическая модель)

Рассмотрим сначала решение задачи об устойчивости изотропных сферических оболочек, построенное на основе классической модели оболочек с неизменной нормалью. Количественный критерий приме­ нимости классической модели будет дан ниже. В случае изотропных оболочек соотношение (8.3) приводится к виду:

Т = - ^ х + — ; х = А2(1 + ф 2) = А2 + п2;

 

Rz

х

 

 

В = В { 1 - г /2).

(8.4)

Минимизируя выражение (8.4) по параметру х, получим классиче­ ский результат:

TKp = P ^ R = 2 v ^ ; Л2р + п2р = А2, =

(8.5)

Как видно, в случае изотропных сферических оболочек форма волнообразования при потере устойчивости (параметры Акр и пкр) однозначно не определена: параметры А и п могут быть любыми, лишь бы выполнялось условие А2р + п2р = А^.. Таким образом, одному значению критической нагрузки Ткр соответствует множество форм волнообразования. Это позволяет предположить [4, 69], что результаты

7*

196 Гл. 8. Устойчивость многослойных сферических оболочек

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

В случае однородных по толщине h оболочек с характеристиками упругости Е н г жёсткости определяются соотношениями

В = Е h , D = 12(1 _ v2y

а формулы для расчёта критических параметров примут вид [8, 81]:

Ркр

Eh2

1,21 Е

R 2

 

Д2’

 

/ M l - * ' 2)

 

К Р + пкР = >L = 3,3

R

(8

.6)

 

8.3.Изотропные трёхслойные оболочки

слёгким заполнителем

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

Т =

D*

х

В

 

Dnc +

д 2 + Т :

 

 

1 + U!X

 

D*

х = А2

(1 + ф2) = А2 + п2.

(8.7)

LU

Д Д 2’

Здесь б?3, S — соответственно модуль упругости заполнителя при поперечных сдвигах и его толщина; Н — расстояние между нейтраль­ ными поверхностями несущих слоёв.

Если заполнитель достаточно слабый на поперечные сдвиги (шх з> 1, К < 0,5 Т*), то из соотношений (8.7) прямой минимизацией

можно получить:

 

Ткр = T0HV l - Ь 2 + К;

 

 

Л2р + гг2р = А2С\ / 1 - Ь 2; Ь = Щ ^<

1;

 

 

Анс2

B R 2

(8.8)

£ )чс ’

При выводе формул (8.8) принималось во внимание характерное для трёхслойных конструкций соотношение D* к D.

Если жёсткость К заполнителя на поперечные сдвиги достаточно ве­ лика (UJX <с 1, К ф 0,5Т^), то анализ устойчивости можно проводить

8.3. Изотропные трёхслойные оболочки с лёгким заполнителем

197

на основе модели прямой линии, В этом случае соотношение (8.7) примет вид:

D

В

D

D

Ъ 2 Х + — ’

х = А2(1 + ф2) = А2 + п2. (8.9)

Rz

х

 

1 + и>х’

Здесь величина шх <^i 1, и, как показал численный анализ, она слабо влияет на форму волнообразования. Учитывая это, проминимизируем выражение (8.9) по параметру х , считая величину D постоянной. Тогда

ГГ

/уж

(8.10)

Т кр= /-g— ;

xKP = \ i V Y T T 0 - £0 = 4 , .

V 1 + £о

 

 

Если пренебречь квадратом малой величины £о, то расчётные фор­ мулы (8.10) примут вид:

Т кр =

7 ;>.

_о_

 

 

 

Т'ж.

B R 2

 

l + ifi-

А кр + « к р

(8.П)

 

 

Как видно из (8.8), (8.10), (8.11), критические усилия определяются тремя величинами: Тцс — критическое усилие, соответствующее обо­ лочке с раздельно работающими несущими слоями; Тж критическое усилие трёхслойной оболочки с абсолютно жёстким —*• оо) на по­ перечные сдвиги заполнителем; К — жёсткость заполнителя на попе­ речные сдвиги — критическое усилие при так называемой сдвиговой форме потери устойчивости (b = 1). Этими тремя величинами (Т™, Тж, К) определяются обобщённые жёсткости трёхслойной изотропной сферической оболочки в задаче об устойчивости.

Соотношение между обобщёнными жесткостями оболочки опреде­ ляет две характерные области потери устойчивости. Так , в области слабых на поперечные сдвиги заполнителей < 0,5Т^, b < 1), как следует из соотношений (8.8), критическое усилие пропорционально жёсткости К заполнителя и слабо зависит от жёсткости Т™ несу­ щих слоёв, поскольку для силовых трёхслойных конструкций харак­ терна зависимость К Т™. Если заполнитель достаточно жёсткий > 0,5 Т*, b > 1), то, как следует из соотношений (8.10), (8.11), критическое усилие определяется обобщённой жёсткостью Тж, а попе­ речные сдвиги играют поправочную роль и слабо влияют на критиче­ ское усилие. Поправка от сдвигов определяется величиной

 

Т0Ж

_ В

А с

Ж

( 8 . 12)

 

~ Ж Ж '

Если Дс <с 1, то можно применять классическую модель оболочек с неизменной нормалью и пользоваться формулами (8.5).

198 Гл. 8. Устойчивость многослойных сферических оболочек

Зависимости (8.10), (8.11) пригодны также для расчёта на устойчи­ вость многослойных сферических оболочек, сопротивляющихся по мо­ дели прямой линии либо по классической модели. Для этого доста­ точно при вычислении величины Г * подставлять соответствующие жёсткости оболочки на растяжение и изгиб; кроме того, необходимо принимать во внимание жёсткость К многослойного пакета на попе­ речные сдвиги [32, 41].

8.4. Устойчивость многослойных ортотропных оболочек (классическая модель)

Рассмотрим ортотропные оболочки сначала на основе классической модели. Критические характеристики таких оболочек определяются минимизацией выражения (8.3) по параметрам волнообразования А и 1ф. Проводя минимизацию по параметру А, найдём [71]:

Т = ^ = к ( ф ) ^ ^ Т Ш = к ( ф ) Т 0; То = | 0 ^ А ;

к(ф)

Р\{Ф) Л2 _ \ 2 ! + V > 2

^2 _

2Д 2

(8.13)

В Д Г

ж

V А

 

 

Величина Го соответствует критическому наружному давлению при осесимметричной форме (ф = 0) потери устойчивости. При этом кри­ тическое значение параметра А становится равным Аж (Акр = Аж).

В общем случае волнообразование неосесимметрично (ф Ф 0), и для определения критических параметров Ткр, Акр, фкр необходимо прове­ сти минимизацию функции к(ф) по ф. Проводя минимизацию, получим формулы для расчёта критической нагрузки и параметров волнообра­ зования при потере устойчивости многослойных ортотропных сфериче­ ских оболочек при действии наружного давления (од < а 2):

Т кр = Р к р И

кортГ),

к,орт

В Д кр) _

/ а 1+ 2 i/3 ~ ^

 

2

 

 

Р 2 ( Ф к р )

у «2 + 2\ f f i

 

 

2

 

I + ф^р

 

А2

2

 

 

Акр

 

 

 

кр'

 

 

 

 

 

^кр

 

 

 

В Д к р ) Г 2 ( ^ к р )

 

 

 

 

/2

\ J (A-/?2)2+ ( a 2- a i) { P \ o t 2 - l k o L i ) - (А -

/32)

I

,01/1.

фкр =

- --------------------

д---------

д---------------------------

 

 

« - 7=-.

(8.14)

 

 

P\0t2-P20t\

 

 

у/Р

 

Как

следует из

(8.14),

параметры

волнообразования Акр

и пкр

в рассматриваемом

случае однозначно определены. Если закономер­

ности, замеченные

для цилиндрических оболочек [4, 69], верны

и для сферических,

то однозначная определённость волнообразования

8.4. Многослойные ортотропные оболочки (классическая модель) 199

при потере устойчивости указывает на возможность удовлетворитель­ ного согласования теоретических и экспериментальных результатов. Как известно, в случае изотропных оболочек, когда форма волнооб­ разования теоретическим путём однозначно не определена, различие между и теоретическими и экспериментальными значениями критиче­ ских усилий достигает Зч-4 раз [8, 18]. Систематизированные экспе­ риментальные результаты по устойчивости ортотропных сферических оболочек нам неизвестны.

При подстановке в формулы (8.13), (8.14) соответствующих жёстко­ стей можно рассчитать критическое внешнее давление для вафельных сферических оболочек [40].

Для однородных по толщине h ортотропных сферических оболочек зависимости (8.14) принимают вид:

укр _ РкрД _ £

Ф',кр

ор т Т0;

Л2

= А2

____1+ \J~fi_____

А1 = 3,46,

I Еч( 1 V\V2)

R

^кр

''я

а 1 + 2^/Д ) («2 + 2\[Ji)

 

ЖI

h'

 

 

 

 

 

 

 

 

 

 

 

_

л/ЁфЁч

h2

_

/2G]2( 1+ yJV\V2 )

(8.15)

 

~

/3 (1 -1 /1 1 ^ )

R '

°рт

V

 

 

 

 

 

Здесь Ei, Е2, G12, v\, щ — характеристики упругости ортотропного материала оболочки.

Если анизотропия оболочек такова, что а\ > «2, то оболочки теряют устойчивость по осесимметричной форме. Это следует из соотношения (8.13), поскольку в этом случае величина к(ф) достигает минимума при фКр = 0 (/сорт = О- Критические параметры при этом определяются соотношениями:

Ткр= ^ = Т о = | У в д ; Акр= А 0; ^ кр= 0 (пкр = 0). (8.16)

Расчёт критических нагрузок в случае а\ = ач также проводится по формуле (8.16). Однако форму волнообразования в этом случае нельзя определить однозначно. Анализ зависимости (8.3) показывает, что критические параметры волнообразования при он = ач определя­ ются формулой

АкР + у/Р-«кР = А^

B 4 R 2

(8.17)

 

“ о Г

Таким образом, форма волнообразования при потере устойчивости от наружного давления в случае а\ = ач может быть любой, лишь бы выполнялось условие (8.17).

200Гл. 8. Устойчивость многослойных сферических оболочек

8.5.Устойчивость трёхслойных ортотропных оболочек

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

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

К i + К 2ф2

R 2

Л2 Fi(tp)

(8.18)

А

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

R 2 К х + К 2ф2 (

R 2 К х+ К 2ф2

(8.19)

А ' АЩ ( ф )

А ~

' AЩ ( ф )

 

Подставляя это значение П в формулу (8.1), получим:

A ^ I A ) . А

1 А + А < А 2.

 

R 2 Х+ F2 ' х + 1 + ф2

 

1 + 1 У

(8.20)

В* В 2 1

1 + ф2

 

 

Проминимизируем это выражение по параметру х и найдём:

'НС

1 -

А А ;

+

А + А А .

Т = к(ф) Тс

 

 

О

 

1 + ф2

 

1 + А

\V

т0™ = д v 'A A

Л2

= Л2

1 +Ф1р

1

( 2А

1 + -^ ^2Р

 

( . )

^кр

''нс

 

 

ук

1

8 21

 

 

^ ( А р ) т2(а р)

^

Г "

+ Ф1р

 

 

 

 

 

 

 

 

 

 

Критическое

значение

параметра ф находится минимизацией по

ф величины Т

в формуле

(8.21). Рассмотрим случай изотропного

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