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

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

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

1.1. Разрешающие соотношения для трёхслойных оболочек

23

При таких предпосылках жёсткости слоёв, из которых состоят

наружные несущие слои, определяются модулями упругости

тДЮ

тлО)

 

, Ь 2 >

G \2\ коэффициентами Пуассона

v^ и толщинами h^ (к — номер

слоя). Если заполнитель выполнен в виде многослойного пакета, то его сдвиговые жесткости определяются модулями Gj3', G ^ поперечного сдвига и толщинами SW составляющих слоёв. Во многих случаях заполнитель однороден, и тогда его жёсткости на сдвиги определяются модулями сдвига G 13, G23 и толщиной 6 (рис. 1.2).

При исследовании НДС, устойчивости и колебаний многослой­ ных оболочечных конструкций из композитов коэффициенты Пуассона в большинстве случаев слабо влияют на конечный результат. Поэтому обычно принимается, что коэффициенты Пуассона одинаковы для всего пакета и равны некоторым приведенным величинам щ, i/2. Формулы для расчёта приведенных коэффициентов гц, v2 даны в Приложении А.

Считаем также, что оболочка пологая, т. е. параметры Ламе её коор­ динатной поверхности мало изменяются и в пределах рассматриваемой области полагаются постоянными.

На основе указанных предпосылок можно получить следующее дифференциальное уравнение для исследования устойчивости трёх­

слойных пологих оболочек вращения [66, 72]:

 

 

D1 Vf

 

д2

д2

) V4*

V|w +

 

 

^2к^79 +

'ду2

 

 

 

дх2

 

 

 

 

 

 

 

 

 

 

+

( §

V4K + D r VfV4 ) V Aw = V4V4Cqw,

 

 

 

 

 

R

2

 

 

 

Vf

дА

дА

a

 

° 4

V74

° 4

94

_ 94

дхА+ах дх2ду2

+ 9 1dy4 ’

v 2 -

dx4 + «2 dx2dy2

+ /?2V ;

 

V4

dA

+ a,.

 

 

 

V4

+ ^

 

 

 

 

 

 

2*

d x A

2 дх2ду2

 

2 dyA’

R Icte2

R\

ду2)

 

 

 

v f = 1 —{ш\ +Щ2к)

д2

 

f t

 

 

 

 

2 —(W2 + VIK) -Q^2 +^1^2к V4*;

 

<*i

а 2

B 2_

2 1/2;

в 12

ZJi

щ .

K f t

 

2 l 2^

1 + , 2

f t

=

£ 2

 

D\ ’

 

 

 

 

 

Хй to

II

 

*

В *2

-

2 v 2 ;

В *2

 

012 ~ T°1F2

(3*2

 

п *

 

 

В \

Г)*

 

 

D *2

и 2 .

^ 1 2 .

й 2к =

 

U 2 =

K

2~ ’

Ш1К =

Ж

К 2

 

r » « i + , s » J L

+ Т °— -

( 1.1)

1дх2

дхду

+ 2 ду2’

 

24 Гл. 1. Разрешающие соотношения для многослойных оболочек

Т,°, Т2°, S 0 докритические безмоментные усилия, действующие в обо­ лочке.

Уравнение (1.1) можно записать также в следующем операторном

виде:

 

 

 

 

 

 

 

 

 

 

 

_

4

 

 

 

 

 

 

 

 

(DT + D*£lxy) Vfu; + | §

w = qw;

 

 

 

 

 

 

 

 

 

,

( _

d2

_

d2

 

 

 

 

 

 

 

x y

 

1

[ Ш 2к

дх2 +Шыд ^

 

( 1.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В формулах (1.1), (1.2) введены следующие обозначения:

В\,

В2у

 

В \2

жёсткости

трёхслойного

пакета

на растяжение-

сжатие и сдвиг В ПЛОСКОСТИ X , у;

 

 

 

 

 

 

В 2 = В 2( \ - щ и 2)\

В *2 = B2(l -

v\v2);

 

 

 

 

4В

в{2)

 

 

а в ^{

в {2)

 

 

4^12 Д [2

 

_

в

 

 

 

Bt

 

 

 

 

 

в12

4д(') # (2)

 

 

в ,

 

 

 

в 2

;

В 12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в \2 , в[2\

 

В ^ \

?(2)

жёсткости на растяжение-сжатие

в\1 в:

'

 

В

\2

 

J2

^12 >^1

>^2

>^12

 

 

 

 

 

 

и сдвиги соответствующих несущих слоев;

 

 

 

D 1,

D2, D 12 —

минимальные

изгибные

и крутильная жёсткости

трёхслойного пакета;

 

 

 

 

 

 

 

 

 

 

D\ = В Г + D*;

D2 = D r + D*2-

D n = D\\ + D*l2;

D r,

D2C, D r — суммарные минимальные изгибные и крутильная

жёсткости несущих слоёв;

 

 

 

 

 

 

 

 

D*, D2, D*2 — минимальные изгибные и крутильная жёсткости

трёхслойного пакета с безмоментными (В"с = D2C=

 

= 0) несущи­

ми слоями;

 

жёсткости

трёхслойного

пакета на поперечные сдвиги

К\,

К 2

 

в плоскостях xz и yz соответственно.

Методы определения указанных жесткостей изложены в Приложе­ нии А.

Решение дифференциального уравнения (1.1) будем искать в виде

[79]:

 

 

 

 

с,у)

. \ х

пу

w(x,y)

Xx + ny

w0sm — cos— ;

Wo sin ■

 

R

R

 

~R

(1.3)

rmrR

A =

£

m, n — параметры волнообразования; R, £ — характерные размеры оболочки.

1.1. Разрешающие соотношения для трёхслойных оболочек

25

Подставляя (1.3) в (1.1), с помощью соотношения (1.2) получим разрешающее уравнение для отыскания критических параметров:

q{X,n) [£>Г + D ffi (А, п)} Ф,(А,п)

B 2R2

Фд(А, п) _

 

 

 

R 2

+

R \

Ф2(А,п) ’

 

 

 

 

Ф|(А,п)

 

 

1 + (W2K^2 + и ып2) Ф!(А,п) .

f2(A, п)

 

ФС(А, п)

 

 

 

 

 

 

 

 

 

 

 

Ф] (А, п) = А4 + «1 А2п2 + /?in4;

Ф2(А, п) =

А4 + а 2А2п2 + /32п4;

 

 

 

 

 

 

 

2

Ф2 (А, п) = А4 + а|А 2п2 + Р2п ;

фд (А, п) = ( ^ 2 + J ^ П2^ ;

Фс(А, п) = 1 + (сЩ+ W2K)A2 ф 2 +

к)п2 + Ш1о;2к Ф|(А, п);

_

_

D*

 

_

 

Щ

 

Д2 “

K \R 2 ;

~ Д2 ~

K r f 2'’

_ WIK _

^t2

_ W2K _

Р 12

WlK

Д2

^ Д 2 ’

W2K

W

Д2Д2:

g(A,n) = А2Т! + 2 5 Агг + гг2Т2;

(1.4)

Ti, Т2, S — сжимающие усилия, действующие до потери устойчивости. Критические значения усилий Т), Т2, S определяются минимизаци­

ей выражения (1.4) по параметрам А и п волнообразования.

При решении конкретных задач часто вместо параметров волнообразования А и и удобно использовать параметры А и ф = —. Производя эту замену, из соотношений (1.4) можно получить:

q(Ki>) = [D^ + D*ln(X,i;)]X2Р \(Ф)

B 2R 2

FR (VQ

R 2 +

R \

X2F2(V>)’

Fi(Ф) = Х4 ф 1(Л>те) = 1 + а\ф2 + AV’4;

р2{Ф) = ^ Ф 2(А,п) = 1+ а 2ф + м> ;

2

р и(Ф) = ^ Фд(А,п) =

1 + т /

26 Гл. 1. Разрешающие соотношения для многослойных оболочек

f2i(A, ф) + А2(сщ + ш2ф2) 0.2(\,ф ) ’

1 + Х2(ш2к + Ш1кф2)

«1 (А,ф)

 

 

 

L '

WIW2K\ 2F2 (Ф)

 

,

/9

 

П2 (А, ф)

+ и!2ф2

 

 

 

 

Я( \ , ф ) = Т 1 + 2 ф 8 + ф2Т2;

ф = ^ .

(1.5)

Уравнения устойчивости (1.4), (1.5) трёхслойных оболочек, постро­ енные на основе модели ломаной линии, совпадают с аналогичными соотношениями для классических оболочек с неизменной нормалью, если величину D"C+ QD* заменить полной изгибной жёсткостью D\ = D"c + D*. Следовательно, величину

у/;фф n r | п D \

( 1.6)

можно назвать эффективной изгибной жёсткостью трёхслойной обо­ лочки. Она представляет собой изгибную жёсткость D\ оболочки, уменьшенную за счёт влияния поперечных сдвигов. Влияние попе­ речных сдвигов определяется функцией ЩХ,ф), которая зависит как от податливостей wi, ш2, ui\K, ш2к заполнителя на поперечные сдвиги, так и от параметров А и ф волнообразования. При отсутствии запол­ нителя податливости на сдвиги стремятся к бесконечности, а функ­ ция Q(A, ф) становится равной нулю, т. е. Dfc + Q, D* = D*c, и несущие слои сопротивляются независимо друг от друга. Если заполнитель абсолютно жёсткий на сдвиги, то податливости равны нулю, Q = 1, и оболочка сопротивляется в соответствии с классической моделью неизменной нормали. Назовём функцию Q (А, ф) функцией влияния поперечных сдвигов. Эта функция изменяется в зависимости от сдвиго­ вой жёсткости заполнителя в пределах от нуля до единицы: 0 ^ Q ^ 1.

В общем случае можно записать

= D?c + ttD*i = Q D i + (\ - Q ) D ^ C.

Для достаточно жёстких на поперечные сдвиги заполнителей величина Q ~ 1, а поскольку в трёхслойных оболочках D*c «С D ь получим

Д]ФФ = П Д .

В этом случае оболочка сопротивляется в соответствии с моделью прямой линии [37].

1.2. Трёхслойные с жёстким заполнителем и композитные оболочки 27

Если заполнитель имеет малую жёсткость на поперечные сдвиги, то £2 1, и D|*c ~ £2D*, что соответствует модели ломаной линии.

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

Рассмотрим слабые заполнители с высокой податливостью на попе­

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

 

 

сщА2.Р2 » 1 .

 

 

(1.7)

Тогда функцию £1 влияния поперечных

сдвигов можно

представить

в виде:

 

 

 

 

 

 

к , + к 2ф2 2 /

 

_ К , + К 2ф2 2\ .

£2 (А,ф)

V

 

 

D*X2F*

)

D\X2Fx

 

 

кх+к2ф2

2

«

1.

 

( 1.8)

Щ Х Щ

R

 

 

 

 

 

 

 

 

Соотношения (1.7), (1.8) характерны для модели ломаной линии.

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

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

 

W,A2F2 « 1.

(1.9)

Тогда функцию £2 влияния поперечных сдвигов можно представить в виде (£2i и £22 и 1):

________ 1________

к 12

1

ЩХ,ф)

к12 + Х2’ к 12

> А2.

1 + А2(сщ + и2ф2)

W1+ и 2ф2

(1.Ю)

Зависимости (1.8), (1.10) используем в дальнейшем для получения расчётных формул.

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

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

В этом случае несущие слои и заполнитель по существу равноправны,

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

1.2. Трёхслойные с жёстким заполнителем и композитные оболочки 29

А

( _

д2

д2 \

| f v | j V > = v |v 4 qw.

V f

к дх2

V 42W -

 

Р

+UJlKdy2)

 

( 1. 12)

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

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

^ V lV2 + | f

V |V 4fi

V|V§. qw;

2 ,

_ д2

д2

 

 

(М 3)

Порядок этого уравнения равен восьми, оно учитывает поперечные сдвиги параметрически и соответствует модели прямолинейного эле­ мента. Заметим, что входящие в податливости ш\, Ш2 жёсткости К\, К 2 на поперечные сдвиги должны вычисляться с учётом неоднородности пакета по толщине. Формулы для расчёта величин К\, К 2 приведены в Приложении А.

Уравнение (1.13) по аналогии с (1.2) можно записать в операторном

виде:

 

 

 

 

Ч

Щ

+ EI

3 L

(1.14)

 

R 2

V |

'

Входящие в коэффициенты и операторы дифференциальных урав­ нений (1.12)—(1.14) жёсткости на изгиб и кручение являются соответ­ ствующими минимальными жесткостями многослойного пакета. Они вычисляются в соответствии с Приложением А.

Представляя решение в виде (1.3), найдём соотношение для расчёта критических параметров, аналогичное зависимости (1.4):

q(X,n)

D 1

Oi(A,n)

B 2R 2 Фд (А, п)

1 + А2сщ + n2w2

R 2

 

R 2

 

(1.15)

 

 

Ф2(А, n)

Переходя к параметрам А и ф п

получим

 

 

 

А’

 

 

 

 

 

A

A2F ,(A

B 2R 2

FR (ф)

 

1 + А2(сщ + ш2ф2)

R 2

 

+

R2

(1.16)

 

 

АЩ(фУ

30 Гл. 1. Разрешающие соотношения для многослойных оболочек

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

П Т = 1 +А2(и Ц + С ^ 2)'

(1Л7)

Сравнивая эту величину с функцией (1.10) влияния поперечных сдвигов в трёхслойных оболочках с достаточно жёстким на сдвиги заполнителем, видим, что они полностью совпадают. Отсюда следует, что соотношение (1.9) может служить критерием применимости модели прямой линии к оболочкам: если выполняется условие

А2(^2к + ^ 1к^2) - с 1; где Ш1к = т

^ =

( Ш )

то можно применять модель прямой линии.

1.3. Разрешающие соотношения для многослойных композитных оболочек (классическая модель)

Во многих практически важных случаях при расчёте многослойных оболочек на устойчивость хорошие результаты даёт теория, основанная на классической гипотезе о неизменной нормали. Из соотношений (1.16), (1.17) следует, что модель оболочки с неизменной нормалью применима, если

A2(CL>I + си2ф2) < 1.

(1-19)

В этом случае дифференциальные соотношения (1.13), (1.14) примут вид:

( L >IVIV 2 +

v !i) w = V 2qw; q = D lV \ + | | ^ .

(1.20)

Зависимости (1.15), (1.16) для расчёта критических параметров запишутся следующим образом:

<?(А,п) E i ® i(A ,n) +

B 2R2 Фд(А, п)

 

R 2

R\

Ф2(А, п) ’

 

А А2

B 2R 2

Дя(А

( 1.21)

д {\,ф )

Д |

а Щ(ФУ

д 2

 

Входящие в зависимости (1.20), (1.21) жёсткости представляют собой минимальные жёсткости многослойного пакета.

1.4. Полубезмоментная модель и непологие оболочки

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

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