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

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

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

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

361

Таким образом, заменяя Ъ выражением Z* (2.30), из уравнения

(1.8.35) получим уравнение устойчивости рассматриваемой сфе­ рической оболочки:

2(Д- f I)2- f 1 — (Д - f 2) w - f

где

+ ^ - (1- й*А)(А + 1-

у)(А + 2)ц; = 0' (2-31)

 

 

 

сз -

ь*

Eh2

 

1 2 ( 1 — V2 ) f l 2 *

1 0 ( 1 -

v 2 ) № G '

Это уравнение должно решаться при соответствующих гранич­ ных условиях для w. В данном случае граничными условиями для w являются условия непрерывности и однозначности решений на

сфере.

Полагая в уравнении (2.31)

 

 

Дw = — \w,

 

 

 

 

(2.32)

получим для искомой постоянной величины Xследующее характери­

стическое уравнение:

 

 

 

 

 

 

 

[с* (1 - * )2+

1 + h\ + ^

(1 - f Л*Х) (1 -

v -

X)] (2 -

X) =

0.

(2.33)

Отбрасывая тривиальное решение Х=2, имеем

 

 

 

с2 (! — X)2- f 1

А*Х — ^§-(1 - f А*Х) (X— 1

v) =

0.

(2.34)

Отсюда для определения критической силы получим

 

 

2 Eh С2 ( Х - 1 ) 2 + 1 + А *Х

 

 

(2.35)

 

q ~

R (1 + А * Х ) ( Х — 1 - f V) *

 

 

 

 

 

 

Рассматривая q как функцию параметра X, легко получить ми­

нимальное значение выражения (2.35),

которое и будет представ­

лять искомое значение критической силы:

 

 

 

 

*

2Eh 2 с — 2 v c 2 +

Ж

 

 

 

 

 

q

~~ R

1 + 2 h*/c

 

 

 

 

 

или с принятой здесь точностью

 

 

 

 

 

 

q*__

(2с — 2Vc2— й*).

 

 

(2.37)

Полагая в этих формулах h *= 0, получим известное выраже­

ние критической силы

для

изотропной

сферической

оболочки:

*

4 E h ,

 

2 E h

 

V2 h

v A s " !

 

/ о o o \

q = - i r ( c - v C2) = ^ r = ^ y l l / - g - д -

6F j -

 

<2-38)

Рассматривая расчетные формулы (2.36)—(2.38), замечаем, что критическая сила, найденная с учетом явлений, связанных с поперечными сдвигами, при некоторых значениях отношений

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

h/R, E/G' может существенно отличаться от критической силы,

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

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

Пусть оболочка загружена по торцам (s=0, s=l/R ) равномерно распределенной нагрузкой с интенсивностью р. Оболочка шар­

нирно оперта по концам. Тогда для основного безмоментного со­

стояния оболочки

будем иметь

 

 

 

 

 

Т\ = — р,

Ц =

0, S0 — 0.

 

(2.39)

В силу этого из (2. 1) получим

 

 

 

 

 

 

 

 

 

 

(2-40>

Будем изучать лишь осесимметричные1формы потери устойчи­

вости рассматриваемой оболочки.

Согласно

(1.3.36) и

(2.40)

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

 

d*w

■2m d^w

■nw •

12Д2

d^w

n

(2.41)

dsi

 

 

h?B

Р Ш

= 0’

 

где, как и раньше, для коэффициентов т и п имеем (1.29).

Решение уравнения, удовлетворяющее условиям шарнирного опирания, ищем в форме

 

. kr.R

(2.42)

w

= 2 > * Sin —т— s,

где l — длина оболочки, R — радиус кривизны, fk — неизвестные постоянные, к — целые числа, характеризующие формы потери

устойчивости.

Подставляя значение w из (2.42) в уравнение (2.41), получим следующую формулу для определения значений р , при которых

возможны осесимметричные формы потери устойчивости:

[7 Д У

/ к к Д \ * ,

1

__ АЗЯ ц [_ (~ )

— 2т\ ~ г ) +

ге

Pic

 

(2.43)

12Д2 fkiz.

 

Отсюда легко определить наименьшее значение рк, т. е. кри­

тическое давление

* _ _ №Вit

{\Jnт),

(2.44)

6Д2

 

 

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

363

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

 

Рассматривая формулы (2.43), (2.44), (1.29), замечаем,

что

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

Bik легко определяются с помощью формул

(1.1.34), (1.1.42).

Ясно также, что каждой группе коэффи­ циентов B ik будут соответствовать свои ко­

эффициенты т и п .

 

 

Рассматривая приведенные

результаты,

видим, что,

меняя

ориентацию

материала

в оболочке,

можно

существенно изменить

характеристики устойчивости оболочки.

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

главные направления упругости совпадают с главными геометри­ ческими направлениями оболочки а, [3, у (см. гл. I).

Очевидно, для слоистых оболочек в общем случае безмоментное напряженное состояние невозможно. В этом легко убедиться, рассматривая соотношения упругости многослойной оболочки (1.11.2). Поэтому при решении задач устойчивости многослойных оболочек, вообще говоря, следует докритическое состояние обо­ лочки считать моментным.

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

Выберем координатную систему так, чтобы коэффициенты первой квадратичной формы были равны единице (Я = 1, В = 1), а координатная поверхность f = 0 проходила внутри какого-либо

слоя оболочки (в частном случае она может совпадать с поверх­ ностью контакта двух каких-либо слоев, т. е. А = 8<) (рис. 63).

Пусть начальное моментное состояние оболочки характеризу­ ется нормальным перемещением w ( а, р) и функцией напряжения

364

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

[ГД, III

<р (а, (3). Тогда для определения начального напряженного состоя­ ния мы можем использовать уравнения (1.14.64)*):

 

(D D°) wL 3 (d) 9 Уд? L(w ,

f ) — Z ,

'

 

£2(Л)<р +

L3(d)w +

VRw-\--^ L(w ,

w) =

0,

(2.45)

 

j

где, как и раньше,

 

 

 

 

 

 

 

 

 

L, (D -

DO) = (D n -

D»u) £

+

(D22 -

DM »L +

 

 

 

 

 

+

2 [(Z>12— Z»?2) + 2 (Z>66— Z>ge)J

d*

 

 

 

 

 

 

 

 

 

 

da* dpz

 

L2 (A) — A22jfcl +

(^66 + 2^1г) 0a‘>d p +

 

<?4

 

 

(2.46)

All dpi »

 

 

№) =

^21^

 

+

d22

2d6C) ^ 5-555+

^12 ^pT.

 

 

VA= — — + — —

 

 

 

 

 

 

 

 

Л

Д2

R x

<?p ’

 

 

 

 

 

 

 

 

, ,

.

d-Шf)°9

.

Г?-ф

0

f)2W

 

 

 

^

 

da- d p

1

</p J

r ~ 2 ;d a 0 $

да dfi

 

 

 

 

 

 

d a

2

 

 

 

(2- « )

 

 

 

 

 

 

 

 

 

 

 

а для жесткостей, входящих в линейные операторы (2.46), в случае ортотропной оболочки имеем

3

=

2 я$.,(з.- V , ) >

 

 

«=1

3

 

 

 

 

п

 

К Р< =

т 2

к8? ~ 8?-i) - 2Д (8« - V . )]-

(2.48)

в= 1

п

Лу* = 4 2

К8*’ - 8W

-

ЗА М -

+

ЗА2 (3 , - 8 .

9 = 1

 

 

 

 

 

 

 

а также следующие комбинации из них:

 

 

А

— £22.

»

А

А

1

 

— Q0

Л 22 — Q0 .

Лвв — ^•66

 

__ ^ 1 2

 

о — С Г

 

 

(2.49)

А

»

___Г2

"

12

о„

~о — l/iit'22

 

12

*) Ради преемственности результатов здесь главные кривизны недеформированной координатной поверхности ■ (= () кг и к2 считаются положитель­

ными, если центр кривизны лежит с выпуклой стороны оболочки.

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

365

(2.50)

Л?, =

1 В Д 2- 2К 1гС12К 12+

К\2Сп\2^,

 

О?, -

[* п К 12С22- (К пК 22+

К\2) С12+ а д , С „ ] 2 ,1,

(2.51)

=

\K\fiu 2К22К 12С12-(- К12С22] 2ц1.

 

Таким образом, системой (2.45) будет определяться началь­ ное моментное состояние оболочки. Однако при определенных со­ отношениях геометрических и механических параметров оболочки это состояние может оказаться неустойчивым, и оболочка под действием малых возмущений 8и> и &р перейдет к новой, устой­

чивой форме равновесия. Тогда во втором состоянии будем иметь

w 1 = w - f - 8 ш ,

<р1 =

tp - f - 8 ? .

(2.52)

Подставляя (2.52) в систему уравнений (2.45), получим следую­ щую систему линейных дифференциальных «уравнений в вариа­ циях» относительно 8ш и 8®:

(D — Z/°)8u>— L 3(d )bf — VflBep = L (8ц;, ?)-f-L(w, Sep),

L 2 (A ) 8x>-f-L3 (d) 8w -f-Vx8u; = — L(bw, w).

(2.53)

 

Таким образом, решение задач устойчивости моментного на­ пряженного состояния многослойной анизотропной оболочки при­ водится к интегрированию двух систем дифференциальных урав­ нений: (2.45) и (2.53). К этим системам уравнений должны быть присоединены граничные условия, которые для системы (2.45) имеют обычный вид, а для системы (2.53) однородны и на основа­ нии (2.52) вытекают из граничных условий начального напряжен­ ного состояния. При этом необходимы будут также полученные ранее представления

366

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

1ГЛ. Ill

Мг= ■

- [(£>п -

 

Щд £

+

(Dlt -

D°12) ^ ] W+

 

~

[( 0 22-

 

ДО,) ^

+

12-

D0J ^ ] ш

+

(2.55)

 

 

 

 

 

 

 

 

+

( d12^2 + d22i^)?>

 

 

66

 

Щв)

д2w

 

^вб

 

 

 

 

 

 

да(J(.

 

с бб дад$

 

 

 

 

т ,

д2<?

тг= S

 

 

д2?

(2.56)

 

 

д$2'

-

* =

dadfj *

 

 

Тг =

Сцб, +

^ 12£2“Ь -^11*1

^ 12* 2’

 

 

 

?2 =

^ 22Е2+

Cl2£5 Н~ К22*2

 

(2.57)

 

 

 

S =

 

 

+

t

 

 

 

 

 

 

 

 

 

 

'^6б'С

 

 

 

 

 

 

 

ди . i [dw\2

W

 

д2W

 

 

£i — да + ¥ w ) ~~1Ь ' х

да2

 

 

 

 

до

,

1 fdw \ 2

W

х-

d2w

(2.58)

 

S2= Ж + У ( ж ) "

 

dfr*

 

(!) =

ди

 

до . dw dw

 

 

О d2w

 

 

Ж

 

 

+

- , X •-= -

~ да др '

 

 

 

 

 

 

да <ЭВ

 

 

 

На основании приведенных уравнений и расчетных формул могут быть решены задачи устойчивости моментного состояния

различных типов анизотропных

слоистых оболочек.

1. Рассмотрим многослойную

замкнутую круговую (Л1=оэ,

R 2= R ) цилиндрическую оболочку, шарнирно опертую по торце­ вым линиям (и = 0, a=l\ I — длина оболочки) координатной по­

верхности Y= 0. Пусть оболочка загружена нормально приложен­ ной к внешней поверхности осесимметричной равномерно рас­ пределенной нагрузкой с интенсивностью q, т. е. Z = q , Х =0 ,

Г=0.

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

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

367

к которой должны быть присоединены следующие граничные

условия:

 

 

 

 

 

при

а ~ О

М х — 0,

w = О,

Т1=

0; |

при

а — 1

М 4 = 0,

w— О

Т1 =

0. j

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

d 2 (5 ш )

. d“’w d 2 ( 8 9 )

(2.61)

^2

d£2 »

 

й 2ш d 2 ( 8ш )

: d^2 (jpa •

К уравнениям (2.61) необходимо присоединить следующие однородные граничные условия:

при

а =

0

8Mj = 0,

bw= 0,

ЬТ1=

0,

оу= 0;

(2.62)

при

а =

1

ЬМ10,

bw 0,

8Г, =

0,

8у== 0.

 

Из системы уравнений (2.59), в силу граничных условий (2.60)

и соотношений

(2.54)—(2.58), получим для

определения w ( а)

следующее

дифференциальное

уравнение:

 

 

 

 

 

 

 

{^21

d^W

 

W

(2.63)

 

 

 

 

Л22#

dot

A ^ m ~ q

 

 

 

 

 

Решая уравнение (2.63) и удовлетворяя граничным условиям (2.60), получим для искомого перемещения и? (а) следующее выра­ жение:

w ~ С1е*“ cos ta -f- С2е~елcos ta -f- C3e,a sin ta -j-

 

 

-f- Cie~m sin ta -f- A22R2q,

(2.64)

где

 

 

 

s — ^22 (On — ^?i) + Щ.\— ^21 V*

 

2Д

+

J ’

(2.65)

^22 ( Д 11 Щ 1 )

~ b ^21 ~H <?21

 

^

 

27? [^22(^11

T5fi) + <i|i|

 

а для постоянных интегрирования с точностью 1 ±е~2‘1да 1 имеем

Cl (a sin tl — И22Яе~*; cos tl) Rq,

C2 — — (a sin tl -f-A^JR) Rq,

(2 .6 6 )

C3= — (a cos tl -f- A 22Re~“ sin tl) Rq, C4— (aesl — A 22Re~sl sin tl) Rq,

368

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

[ГЛ. III

где

а_ - М 22 (Рц ~ Р °и ) + ( « 2 - t * ) R л ,

- 2sIM 22(JD11-B ? 1) + df1]

Впоследующем нам необходимо также знать значение cPy/da.2.

Из (2.59), согласно (2.54)—(2.60), легко получить

d2cp 1 ( д d^w , w\

(2.67)

~A ^V 21de*"“г'д / ’

Точное решение уравнения начального моментного состояния оболочки (2.64) представим в виде бесконечного тригонометриче­ ского ряда по косинусам:

 

 

 

 

u?=

2 /fcCosXfcK’

 

 

 

 

 

 

 

 

 

 

 

к = 0

 

 

 

 

 

 

(2 .68)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 = 2 * * cosX*a-

 

 

 

 

 

 

 

 

 

 

 

fc=0

 

 

 

 

 

 

где согласно второму уравнению системы (2.59)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.69)

Далее,

в силу (2.64)— (2.66)

для коэффициента разложения fk

получим

СкС3) е“

 

 

Скс4) е - 1] (-1 )* cos Я +

 

 

4 =

{[(<4 -

-

(С, +

 

 

+

Г

+

С3) е*' +

(С*С2-

С4) е -'] ( - l ) fc sin Я -

С, +

С2 +

 

 

I

г

 

г \\ —

82+

+

х%______ 1_

 

 

 

 

 

 

 

 

I

(S2 + f2 +

>.2 )2 _

4S2\*

Т

 

 

 

 

 

 

 

 

 

 

 

I

кк

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.70)

 

 

 

 

 

 

 

 

-f- A22R 2q — ^ cos -j- ос da,

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, = 7

( 1 -

, . + ? + x i ) ( * = 0 ,

i.

2, 3.......со).

 

Решение системы уравнений устойчивости (2.61)

будем искать

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

 

 

 

 

 

8m =

COSp„p2

amSinXma>

Xm=

^

 

 

 

 

 

 

 

 

m=l

 

 

 

 

 

(2.71)

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

8?

=

COS fl.p 2

Ь»

Sin Хт й>

t1» =

J .

 

 

m=l

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

369

Представления (2.71) удовлетворяют граничным условиям (2.60) и условиям замкнутости оболочки по координате (3.

Подставляя в (2.61) значения w,d?<f/dcP, bw и 8<р соответственно

из (2.68) и (2.71), получим следующую систему алгебраических уравнений:

«

-

К

-

 

-

 

 

 

+

4

 

 

 

=1

 

 

 

 

 

 

 

 

+ 2 Ф т-А,— 2

 

2

 

 

 

 

 

 

m=r

» i= l

 

я»—1

 

 

 

 

 

 

 

- 2

 

4- 2

х2* Л * А . )

- о,

 

(2.72)

 

 

т=г

 

 

т=1

/

 

 

 

^Зря

]|) Я р “Ь ®2p»^p4~ “2~ ( 2

*r-iJr-mam4~

 

 

 

 

 

 

 

' *и=1

 

 

 

 

 

 

 

 

 

4” 2Kn-rfm-ram

2^m+rfm+ram\~

'

 

 

 

Wl=r

 

 

=1

/

 

 

где для Ф<?„ имеем

 

 

 

 

 

 

 

 

=

(Ди -

Хг +

~ Л"*) ^

+

 

 

 

 

 

 

 

+ 2 [(/ > „ - Д?2) +

2 ( А * -

В Д

х*р «,

 

 

Ф 2рт =

^22^p

4 " И 66 4 - 2^12) ХрРя'

 

 

 

 

ф3рв =

d2l\* +

4- (dn 4- d22-

2dee) X*X*.

 

 

 

 

Исключая b из системы уравнений (2.72) и пренебрегая не­

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

p — i

(V p — 2WU+ W2p) ар— 2 t iw p-ш— w р+т)ап

-

2 {Wm. p- W m+p)a m = 0.

(2.74)

m=p+l

Здесь

 

+ Ф 1 рпФ 2 я »

Vr = 2 R A „ ±

/*2Ф2ря

 

(2.75)

(Z = p — /га, /га — р, m -j-p ).

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

370

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

1ГЛ. Ш

Очевидно, тривиальное нулевое решение системы (2.74) нас не будет интересовать. Тогда, приравнивая нулю детерминант матрицы коэффициентов системы (2.74), получим условие для определения критическоё силы. Для наглядности приведем уко­ роченный вид этой бесконечной матрицы:

I Vi - 2W 0 + W1+1

 

 

+

w 2+1

“ ^3-! +

^3+1 • • •

—^2-1 +

^1+2

v 2 -

2W

0 +

W 2+2

—И'з-г +

^3+2 ■• •

—^3-1 +

^1+3

~

W3~2

+

И^г+З

F3 — 2W o +

W$+3 . ..

Таким образом, поставленная задача, т. е. определение кри­ тического значения внешнего давления Z = q * , привелась к отыска­

нию наименьшего собственного значения матрицы (2.76). Рассмотрим численный пример. Пусть оболочка составлена

из двух ортотропных слоев одинаковой толщины =0,5 h (t= l, 2).

Пусть, далее, для механических характеристик слоев имеем

Е{Ъ =

Ш ,

E W = E ,

vp =

v(i>=0,

G [» =

0,2E,

E[V =

2E,

Е р = 0,5Е, v(*) =

v ^ )= 0 ,

G $ = 0 ,2 E .

Принимаем также, что

1/Я=п/Ъ,

h I R — 1/100 и

A = h , т. е.

закреплена внутренняя поверхность оболочки*).

Тогда, согласно (2.48), (2.49) и неоднократно использованным

обозначениям

В'и =

ЕЦ(1 — vjv‘),

В{2=

ЕЦ( 1 — vjv*), В\2=

= ч\ЕЦ(\ — v‘v‘) =

V} Е\ (1 — V*v|),

В*е =

G{v

для коэффициентов

А а*

D iv

D ik получим

 

 

 

 

 

 

 

A l1 ~ 6E h »

А 22 — 3£А »

Л в — E h »

 

А Ч — 0>

dii =

2

d22 ~

7

^ 6 ==

1

 

 

==

g" h,

j2

~2

 

^ 12

 

Dn =

3Eh\

Dn= ^ E & ,

DM=^E № ,

D a =

0,

 

|b7r,

D\2= ^ E h \ D

l = L Eh\

A °2 =

0.

 

D n-

 

 

 

 

 

 

 

 

 

 

Из

(2.73) для Ф(

получим

 

 

 

 

 

 

 

 

ф.

F.h3 . . .

-

 

Л . . .

 

 

 

 

*40

 

 

 

 

 

 

 

 

фз ^ = - Ш

АХ/ (т ) 2

 

 

 

 

 

*) Иные значения Дздесь не рассматриваются. Однако укажем, что спе­ циальным исследованием, выполненным на основании изложенной здесь теории, установлено влияние величины Д на значение критической силы.