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

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

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

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

где

Далее, из (2.75) имеем

гг- 2 ,

db

1 + M)J +

625 VA/mP4

Vp = 2R±

(i — P т, т р, т -\ - р ) .

И наконец, из (2.70), согласно (2.66), а также приведенным выше численным результатам, для f { найдем

f j qh \ cos ™ da,

о

где q = A i2R*q*/h — безразмерное внешнее давление.

Подставляя значения Vp и Wt в детерминант матрицы (2.76)

и приравняв его нулю, получим для критического значения без­ размерного давления <7=0,223, что имеет место при ге=13.

Теперь, если ту же самую задачу решать, формально полагая, что начальное напряженное состояние оболочки является безмоментным, то для безразмерного критического давления получим /7=0,251 при п=13.

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

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

нормально приложенным к координатной поверхности

у = 0 .

Пусть оболочка

перекрывает прямоугольный

план

(a X b )

и шарнирно оперта

по всему контуру (а = 0, и=а;

р = 0 ,

(3=Ь)

24*

372

 

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

 

 

[ГЛ. Ш

(рис. 64), т. е. имеют место следующие граничные условия:

 

при

я =

0,

я =

а

М 1 — 0, Т х — 0,

w— 0,

г;=

0;

)

(2.77)

при

р=

0,

р=

Ь

М 2=

0,

Т2 =

 

0,

w=

0,

гг =

0.

J

 

 

Система

разрешающих дифференциальных

уравнений (2.45)

в развернутом

виде

перепишется

следующим

образом:

 

F 1= Л22й

+ (Лв +

2Лг) йгзр*

 

^ 4 -

 

 

 

 

 

л п

 

+

 

 

 

 

 

 

 

diw

 

 

 

 

diw

 

d*w

,

 

 

 

+

^21 ^ Г + (^11 +

^22 ~

2d66) daft d№

■°»дрГ"Г

 

 

 

 

 

 

i

7,

d^w , ,

d2w , d 2w d2w

f d2w \2

 

F2S

M

)

д*Ы 7

XT +

 

 

 

da* T

+

22

•Л°я)

 

 

—d

£f.

 

 

■12<Jf|4

2[(ДИ-

a?,)+2(7)6e-

JD«e)l

d i w

(2.78)

а с, г / 7~,

 

п

л v

i

r wT \ nn , ,

 

u

12

 

1 2/

^ \^66

 

ев/1

dfi'2- *"

 

^

Й - К

г +

^22 -

 

2d66) - * 2 -----

 

 

dai

 

 

 

 

 

da2

 

,

d2w

*

d2w

d2w

 

d2f

 

 

 

 

 

 

da2

da2 dp2

 

 

 

 

d2w d2<p

+ 2

d2tf

d2w

■ 9 = 0.

 

 

 

dp2 da2 ~i~ da dp da dfi

 

 

На примере рассматриваемой оболочки осветим некоторые осо­ бенности выпучивания несимметрично собранных слоистых гиб­ ких оболочек.

При граничных условиях (2.77) система уравнений (2.78) может быть проинтегрирована методом Бубнова—Галеркина.

Решение системы

(2.78)

представим следующим образом:

w — f sin Хяsin (ф, X—

,

 

(2.79)

ср= ф sin Хаsin fip, fi.=

.

Представление (2.79)

удов­

летворяет условиям шарнирного опирания (2.77), а искомые функции найдутся из системы урав­ нений (2.78).

Согласно методу Бубнова — Галеркина необходимо выписать следующие уравнения:

а Ь

 

J J Fxsin Ха sin da dp = 0,

 

о

о

(2.80)

а

Ь

 

ij | F2sin Хя sin р.р da dp = 0.

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

373

С помощью (2.78)—(2.79) осуществляя преобразования (2.80)

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

инагрузкой:

 

 

 

Г/3- * / » +

tf = 9,

 

(2.81)

где

 

 

 

 

 

 

 

t =

(Du - Don) X*- f 2 [(D„ -

Z>?s) +

2(D U -

Dg,)J X y +

 

_ l_ /n

n n \ ,4 I |fe]jj-2 - b k y k 2 —

^21 ^

— 2 d 8g)X2[i.2 — d 12^ .*|2

 

22

2

 

+ (^66 + 2J412) X2p.2+ Ацр*

s __ ^вХ[1 fejjJ.2 -[_ fcg\2 —

— (^11 ~b ^22 — 2dgg) X2[l2 —

 

 

o b

Л 2 2 Х4 +

(^ 6 0

+ 2 4 J 2) X2(i2 +

(2.82)

r _ _ 512X2(12_________________ 1_________________

 

 

 

 

a b

 

 

 

 

 

q =

^ j

j q sin Xa sin fip da dp.

 

 

 

 

о

0

 

 

 

 

 

 

Существенное влияние на характер зависимости / — q оказы­

вает коэффициент s при /2.

 

 

 

 

Рассматривая (2.82),

замечаем, что

коэффициент

s сущест­

венным образом зависит не только от кривизн координатной по­ верхности 7 = 0 (что имеет место в однородных и симметрично

собранных слоистых оболочках, когда координатная поверхность 7 = 0 является срединной поверхностью оболочки), но и от коэф­ фициентов dik и тем самым от жесткостей взаимного влияния К {к,

которыми характеризуется несимметричное строение слоев по тол­ щине оболочки.

Несимметричность строения оболочки по толщине при задан­ ных значениях кривизн кх и к2 может привести к увеличению или уменьшению коэффициента s, что в свою очередь приведет к уве­

личению или уменьшению склонности оболочки к хлопку. Интересно заметить также, что при некоторых значениях жест­

костей K ik даже пластинка (fc^O, к2= 0) может терять устойчи­

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

Таким образом, при заданных значениях кривизн многослой­ ная оболочка в зависимости от характера слоистости может пре­ терпевать деформации различных типов. В случае анизотропных слоистых оболочек исследование характера деформирования сле­ дует вести, рассматривая совместно как геометрические, так и ме­ ханические параметры оболочки, т. е. исходя из поведения коэф­ фициента s уравнения (2.81).

В случае, когда s достаточно мало (очевидно, для этого нет необходимости, чтобы оболочка была весьма пологой), или s < 0 ,

374

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

[ГЛ. III

с увеличением стрелы прогиба / параметр нагрузки монотонно возрастает (рис. 65, а). С увеличением s кривая q—/ приобретает

предельную точку и оболочка может деформироваться хлопком (рис. 65, 6). Дальнейшее увеличение параметра приводит к диаг­

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

Рис. 65.

начальной ветви (рис. 65, в). Эта диаграмма, как известно, отно­

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

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

доразумениям.

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

3.

чивания длинной слоистой ортотропной цилиндрической панели, шарнирно закрепленной по длинным сторонам (р = 0, р= Ь ) и загруженной нормально прило­ женной к координатной поверхности у = 0 , равномерно распределен­ ной нагрузкой с интенсивностью q (рис. 66).

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

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

375

Очевидно, для рассматриваемой панели

R 1 = о о , R 2 = R , Z ~ q .

Имея в виду рассмотрение задачи цилиндрического изгиба, из

системы уравнений (2.45),

согласно

(2.46), (2.47),

(2.49) и

(2.56), получим следующее

уравнение

равновесия для

панели:

А,

А- 52\

d^W

Та

I

(2.83)

"c^7 W

W

^

+ q -

 

 

Здесь и в последующем принимаем, что усилие Т2 по ширине панели остается постоянным. Вообще говоря, Т2— Т2 ({3), однако

точные исследования показывают, что в достаточно широком диа­ пазоне изменения прогиба центральной линии (Р=5/2) панели (О <1 w/h ^ 2,5) изменяемостью Тй по ширине панели можно

пренебрегать.

Пусть сила Т2сжимающая. Тогда, полагая Т2 = — Т, можно пе­

реписать уравнение (2.83) следующим образом:

 

 

d*w ,

. 2Й2Ш

X2 I

9^22

 

(2.84)

dp

dp

R 1 DnCn — K*

 

 

 

 

X2—

тс22

 

 

(2.85)

 

^22^*22

 

 

 

 

 

 

 

Общее решение уравнения (2.84) представится так:

 

 

w = C1sin Хр-|- C2cos Xp-f- C3p-f-C4-f

у *

(2.86)

Постоянные интегрирования Ci определим из условий шарнир­

ного закрепления:

 

 

 

 

 

при р = 0, р= Ь

w=

0, М 2 — 0.

 

 

(2.87)

С помощью соотношений упругости и условия 8Х=0 удовлет­ воряя граничным условиям (2.87), получим для нормального перемещения оболочки выражение

i - ) ^ l l p (2.88)

при условии, что sin ХЬ^=0 .

В случае же, когда sin Х 5 = 0 , т. е. ХЬ=2тш (ге=1,2,3,. . .),

имеем

» = С. М т - у ) ] , ™ » хе - ‘ > -

- ( т - т ) ^ С -

<2-89>

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

что соответствует несимметричному выпучиванию панели относи­ тельно сечения р = 6/2-

Из приведенных формул (2.88) и (2.89) нетрудно получить

формулу для определения прогиба центрального сечения панели

р=6/2:

f ~ [Й + А ("7"“ i ) ] ( ' ~ 1) + ( х т ) Г ■

(2.90)

cosy

 

В окончательных представлениях (2.89)—(2.90) пока еще оста­ ется неопределенным внутреннее усилие Т.

Значение Т определим из условия несмещаемости опор оболоч­ ки. Полное взаимное сближение опор 8 должно быть равно нулю,

т. е.

*

5 =

- 5

^

р = ° -

<2 -9 1 )

Из формул (2.52)—(2.58) для dv/dfi в случае цилиндрического

изгиба получим

d v

w

1 /dig\2

. K ^ d - w

Т_

(t\ qr>.

dp

R

2 \dp)

"Г C22 dpa

C22‘

*

>

Подставляя значение dv/dfi из (2.92) в (2.91) и учитывая

(2.88),

будем иметь

г ^ = М ё + р ( т - ? Ж ‘ ^ - т - ¥ ^ т ) +

+ я ( £ - £ ) - ‘ & ( т - т ) -

<2-м >

Для дальнейших рассуждений целесообразно ввести безраз­ мерные параметры

9 =

дЫ

 

к - Л

 

 

 

* ~ R h ’ *

A ’

Г

т

 

(2.94)

е=

* 2 2

 

 

 

 

с22а *

¥° ™ - с 2а

характеризующие соответственно нагрузку, кривизну, прогиб, усилие и эксцентриситет закрепления.

Учитывая (2.94), из (2.90) и (2.93) получим уравнения, ко­ торые устанавливают связь между расчетными величинами

S г]

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

377

(q*,

р, 5, е) панели:

 

 

 

[q* -J- 4'f2(4р2е— Л)]2{2р3—(- 3[5 (ср — tg <р)-f- ? tg2<р]} -|-

 

 

-f. 16ф2 [д* -f- 4ср2(4ср2е— Л)] [(/с— 4ср2е— 12е) р3-f-

 

 

- f

3[(/с — 4<р2е) (р — tg cp)- f 4<p2etg <p]> - f

 

 

Г

КЪ

1

 

£= « ( г г у О + з У 2 ( s b - ' Ь Ф ' - w <2-9l?>

В случае же, когда нагрузка приложена со стороны вогнутости панели и сила Т%является растягивающей, взамен (2.95) и (2.96)

получим

[д* -f- 4ср2(4р2е А)]2{2ср3— 3 [5 (<р — th. <р) — р th2р]} —

— 16<р2 [д* -f- 4р2(4-f2e-f- &)] ((k -f- 4р2е— 12е) р3

— 3 [(& 4р2е) (р — th р) — 4р2еth р]} —

s= ‘ (a r?-1 ) + w [ 2 (5Ef-1)+ f,^

+

w

<2-98>

Для иллюстрации

на примере панели

с

параметрами

12 (/>22С22—Х|2)/С|2Л2= 1,

/с=15 покажем

характер изменения

закономерности д*— £ в зависимости от величины эксцентриситета закрепления е и от знака нагрузки q.

Задаваясь значением р , из (2.95) или (2.97) (в зависимости от знака силы Т г) определим соответствующие значения q*. Далее по найденным значениям q* из (2.96) или (2.98) находим соответ­

ствующее значение относительного прогиба.

На рис. 67 приведены кривые зависимости «нагрузка — про­ гиб» в случае, когда тангенциальная сила Тг сжимающая, т. е.

когда равномерное давление приложено со стороны выпуклости оболочки. Графики построены для двух значений эксцентриситета закрепления е. Значение е=0 отвечает центральному шарнирному

закреплению оболочки (оболочка

однородная, с толщиной А/2),

а значение е = —0,25 относится к

случаю закрепления нижних

кромок оболочки.

 

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

378

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

[ГЛ. III

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

панели прогибы являются монотонно возрастающими функциями давления (рис. 68). Что же касается вопроса влияния эксцентри­

ситета закрепления, то и в этом случае перемещение опор к центру кривизны приводит к увеличению прогибов оболочки.

§ 3. Некоторые задачи анизотропных оболочек, подверженных действию динамически приложенных нагрузок

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

Из этой серии задач будут рассмотрены досконально изу­ ченные задачи динамической устойчивости, а также задачи удара.

3 ]

ДИНАМИЧЕСКИ ПРИЛОЖЕННЫЕ НАГРУЗКИ

379

 

1. Задача динамической устойчивости многослойной

орто-

тронной пологой оболочки. Рассмотрим нелинейную задачу дина­ мической устойчивости слоистой ортотропной пологой оболочки, собранной из нечетного числа однородных ортотропных слоев, симметрично расположенных относительно срединной поверхности оболочки (см. гл. I, § 10).

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

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

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

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

пряженное состояние является

безмоментным и характеризуется

тангенциальными силами Т\ (а,

р, t), Т\ ( а, р,

t).

 

Ha основании принятых предположений из

(1.5.46),

согласно

(1.5), (1.6), (2.1) и (2.54)—(2.58), получим следующую

систему

дифференциальных уравнений динамической устойчивости поло­

гой гибкой

оболочки (.4=1,

В = 1,

кг= Щ 1, кг= Щ *) (рис. 69):

 

- L 1

 

\

I

л

^

4 -

 

 

 

 

e-h~Ai2) да2

 

 

 

 

 

 

 

 

 

,

,

d2w .

,

d2w

, d2w d2w

 

 

 

Da +

2 D j » L - 9+

D

*

-

 

 

(3.1)

 

 

 

__,

<?2<p_____д-w d

д- w d*w .

 

 

 

— fc2<fo2

 

 

 

 

dp-

da°-

 

 

4

2 d2?

d*w

rn ^ w

i To

I Г

d- w

= 0

,

"'f

firt riQ.

/)*г /)f«

 

 

 

 

 

 

 

да d& da d'ji

 

 

 

 

 

 

 

 

380

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

[ГЛ. 1U

где для коэффициентов уравнений при осесимметрично собранной ортотропной оболочке имеем

__ £22 д _ . £ п д

Qn , ^22— Q0 . ^ 0

 

I ___£12

о — г

г

 

 

 

 

 

*12----- Q

»

-- °П°2

 

 

 

 

 

 

 

m

 

 

— ^ * + l ) j»

 

 

s r * + 1 / w

+ 2 s ; * ( v

 

 

 

 

* = 1

 

 

 

 

 

 

 

 

 

 

“ 1

 

_2

 

 

 

 

 

 

 

D

 

 

 

 

 

-

^

+ 1 ) J ,

3

 

 

 

 

 

 

L

 

s = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г

 

m

 

 

 

 

1------------1 +

 

СР=

L

 

8Z=1

*

1

Si“-

 

 

L

A

+ I + 2

t

 

 

 

(3.2)

(3.3)

(3.4)

здесь у,. — удельный вес материала i-ro слоя оболочки, g — уско­ рение силы тяжести, hi — расстояние г-го слоя от срединной по­ верхности оболочки, которая одновременно является срединной поверхностью среднего (те+1) слоя (рис. 33).

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

Исходные уравнения рассматриваемой задачи в вариацион­ ной записи принимают вид

 

|Y2(4 )!p -f-V ^ -f-y£ (u >, n>)j8(f>dadp =

0,

 

\$ [^ i

w Vfltp L (w, с р ) + 7 ^ + г ф

+

(3.5)

 

+ Cp S - ]5l/7dadP =

0’

 

где, как всегда,

L2{А) =

А22^

+

(Лвв -(- 2А12) да2^

,

 

Ll (D) =

Dn ; g

+

2(Z>12-f-2Z>ee)

______D

JL

 

*

22

*

 

 

 

 

 

да*

* ‘ = т г , & + ж £ ' k ' = R ' 1'

 

 

г /

\

о Р 2"

d'w ( </iwV I

 

 

L{w,

ш) =

2[_—

 

 

 

 

j .

v

d2w d%f

. d^w дЪр

0 <?2sp (Pw

 

L \w’ V — -^2 +