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

книги / Оболочки и пластины

..pdf
Скачиваний:
21
Добавлен:
19.11.2023
Размер:
71.66 Mб
Скачать

границе области

неустойчивости не

 

возникают

пластические

деформа­

CM

ID

ции, оболочку следует называть по­

Ef

логой.

 

 

 

 

 

 

справед­

 

Приведенное решение

Ч

ливо для оболочек с постоянной га­

СЗ

уссовой

кривизной

(цилиндриче­ Н

ская, сферическая и др.).

 

 

 

 

Из решения видно, что первое

 

приближение дает завышенное зна­

 

чение «верхних» и заниженное зна­

 

чение «нижних»

критических нагру­

 

зок. Видно также, что с увеличением

 

отношения сторон панели жесткость

 

оболочки

уменьшается.

 

 

 

 

 

 

Таким образом, первые два при­

 

ближения дают

практически точное

 

решение.

Анализ расчетов показал,

 

что последующие третье и четвертое

 

приближения

не

вносят

существен­

 

ных поправок ни в величину проги­

 

бов, ни в величину критических на­

 

грузок.

 

 

 

теперь

результаты

 

 

Рассмотрим

 

решения

об

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

панели

 

цилиндрической

жесткой

пологой

 

оболочки,

 

подвергнутой

 

сжатию

 

вдоль образующих [10]. Аппроксими­

 

рующие функции выбирали

в виде

 

суммы двух, трех и четырех членов

 

двойного тригонометрического ряда,

 

удовлетворяющие условиям шарнир­

 

ного

закрепления кромок, применя­

 

ли процедуру Бубнова—Галеркина и

 

полученные

 

системы

нелинейных

 

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

 

вали методом итераций на элек­

 

тронной

машине

«Стрела»

в

Вы­

 

числительном

центре МГУ.

 

 

 

Результаты вычислений получе­

 

ны в виде таблиц, которые перенесе­

 

ны на графики (нагрузка — прогиб),

 

часть из которых приводим здесь.

 

Рассмотрим, например,

график

 

на рис. 5.19,

построенный

для

ци-^

 

линдрических панелей с отношением"

 

сторон у —1

и суммой главных кри-

 

 

 

 

 

ьл(л

k Ъ2

 

 

 

визн

+ и2 =

h— Н—hг— >равных

 

0; 6;

18;

24;

32.

Тонкими

линия­

 

ми нанесены

графики

зависимости

 

параметра

прогиба

х г

=

- j - (f

 

СО СМ 05 СО

СМ

0 0

2

О°0*ч

о*ю о*

Is-см

00 СМ ID

CM CD ОО СО оОО

Tf CD

м»CD СО о00

00 CD

00 СО —

^СО О

00

СО 05 СМ

ЮCD ^

СО ч* 00 CD — СМ

<5 a

смоо о

0,64

3,63

6,5

CD

 

 

00 Ю

 

 

 

 

 

 

0,0

0,0

0,0

 

ю —

0,0

0,018

0,02

о м* см

 

 

 

(NJNCN

0,64

3,64

6,48

— г>-

 

 

 

 

 

оо ю

 

 

 

Ю

CM ^

141,2

93,4

119,4

0 0 —

CM

 

 

 

T f CM

 

 

 

 

 

 

141,2

93,4

119,4

__ID CM

141,2

93,5

119,4

00 CM CO

CM

CM

ID

 

 

 

 

CO CM CD

142.1

90,2

112.1

ID

CM

 

 

 

4

CM

II

 

-о I a

стрела прогиба) от параметра нагрузки г* = - - при решении задачи

Енг

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

рического

ряда. Зависимость нагрузка — прогиб

в этом

случае имеет

 

 

вид

 

 

 

 

 

 

 

(У )х 1

12(1

■ х +

 

16и.

X2 +

 

 

 

 

 

 

 

— V2)

/ т

Т

" ’ О

у

У

 

 

“ ’ О

 

 

512

,

 

 

 

 

 

 

+ ----------------------

X3.

 

 

 

 

 

 

9я2( v +

т У

 

 

 

 

Здесь х2 = —-------параметр кривизны. h

Результаты решения задачи в более высоких Приближениях (два, три и четыре члена ряда) ввиду их малого отличия друг от ДруГа укла-

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

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

вчетвертом приближении от

третьего, третьего

от

второго

не­

^

£hz

значительно, что видно также из

 

 

табл. 5.3, где для некоторых про­

 

 

гибов (*1=0,2; 1,0; 2,0) даны зна­

 

 

чения

параметров

нагрузки

гг-

 

 

(/= 1,2,3,4),

получаемых

в

 

раз­

 

 

ных приближениях,

и

разницы

 

 

между

ними Д,- = Г/

 

100%.

 

 

Видно,

что

если

П

 

 

 

 

 

Д!2 достигает

 

 

иногда значительных величин, по­

 

 

рядка 50%, то Д2з, Дз4, Д24 нигде

 

 

не превосходит 0,5%. Это говорит

 

 

о сходимости процесса Бубнова—

 

 

Галеркина

для подобных

нели­

 

 

нейных -задач.

 

зависимости

 

 

Аналитические

 

 

нагрузка — прогибы для

решений

 

 

задач в более высоких приближе-,

 

 

ниях здесь не приводятся

ввиду

 

 

их громоздкости.

 

 

 

 

 

 

 

Легче пользоваться графика­

 

 

ми и таблицами. В табл. 5.4 при­

 

 

ведены решения для цилиндриче­

 

 

ских панелей с отношением сторон

 

 

у = 2, которые графически изобра­

 

 

жаются подобно

рис. 5.19.

 

 

 

 

Рассмотрим

теперь

результаты решения задачи по определению

нормальных, касательных и интенсивности напряжений в срединной по­ верхности гибкой пологой оболочки при том жр выборе аппроксимирую­ щих функций [И]. Нормальные и касательные напряжения в срединной поверхности гибкой пологой оболочки вычисляются по формулам (5 15)

и имеют вид

 

 

v

 

П*П*

т п х ,

п п у

д у 2

££ Ь2

sin ----- sin

Ь

 

тп

а2ср

т п п 2

л

т пх

п п у

д х д у

А тп cos ----- COS

ab

 

 

 

(т — 1, 3, 5, 7;

п

 

1, 3, 5, 7).

 

Интенсивность напряжений определяется формулой

2

V (ах ауУ* + т2.

 

cti-fa.

Д|2 »%

А».%

*54,%

Лм,%

0

0

1 О iz

10О

0,2

3,672

3,673

3,673

3,673

0,02

0,00

0,00

0,00

1,0

5,056

5,0{>2

5,052

5,052

0,07

0,00

0,00

0,00

2.0

9,379

9,252

9,251

9,251

1,35

0,10

0,00

0,10

0,2

4,097

4,098

4,098

4,098

0,02

0,00

0,00

0,00

1

3,533

3,543

3,543

3,543

0,28

0,00

0,00

0,00

2

5,422

5,4094

5,4096

5,4096

0,24

0,003

0,00

0,03

0,8 к .н .

3,5016

3,5082

3,5082

3,5082

0,18

0,00

0,00

0,00

0,2

6,347

6,346

6,346

3,346

0,015

0,00

0,00

0,00

1

.3,840

3,840

3,840

3,840

0,00

0,00

0,00

0,00

2

3,300

3,329

3,330

3,330

0,87

0,03

0,00

%0,03

к.н.

3,159

3,179

3,179

3,179

0,63

0,00

0,00

0,00

0,2

10,47

10,41

10,41

10,41

0,57

0,00

0,00

0,00

1

6,190

5,941

5,941

5,941

4,02

0,00

0,00

0,00

2

2,997

3,012

3,012

3,012

0,50

0,00

0,00

0,00

к.н.

2,586

2,642

2,642

2,642

2,16

0,00

0,00

0,00

 

0,2

 

27,01

26,98

26,97

26,97

0,11

 

0,037

 

0,00

 

0,037

32

1,0

 

18,01

17,61

17,61

17,61

2,22

 

0,00

 

0,00

 

0,00

2,0

 

9,349

8,835

8,836

8,836

5,50

 

0,01

/

0,00

,

0, 0 J

 

к.н.

1

0,3183 /

0,7 1 6 0 /

0 ,7 1 2 5 /

0 ,7 1 2 5 {

0 ,9 6

/

0 ,5 6

0 ,0 0

1

0 ,5 6

 

/

 

 

 

 

 

 

 

 

 

 

 

 

ai+ a*

с

0

а

0

1 9

Q 9

 

 

 

 

 

 

Aiti%

Дгз»%

Д.4,%

Л*4,%

0 ,2

5 ,643

5,643

5,643

5,6 4 3

0 ,0 0

0,0 0

0 ,0 0

0,0 0

1,0

6,571

6,572

6,572

6,572

0,015

0,00

0,00

0,00

2 ,0

9,348

9,312

9,312

9,312

0,38

0,00

0 ,0 0

0,00

0 ,2

5,958

5,958

5,958

5,958

0,00

0,00

0,00

0,00

1,0

5,958

5,602

5,602

5,602

0,07

0,00

0,00

0,00

2 ,0

6,808

6,813

6,813

6,813

0,07

0,00

0,00

0,00

к.н.

5,576

5,580

5,580

5,580

0 ,07

0,00

0,00

0,00

0 ,2

7,397

7,397

7,397

7,397

0,00

0,00

0 ,0 0

0,00

1,0

5,792

5,795

5,795

5,795

0,05

0,00

0,00

0,00

2 ,0

5,447

5,464

5,464

5,464

0,31

0,00

0,00

0,00

к.н.

5,364

5,375

5,375

5,375

0,22

0,00

0,00

0,00

0 ,2

13,77

13,77

13,77

13,777

0,00

0,00

0,00

0,00

1,0

9,683

9,667

9,667

9,667

1,16

0,00

0,00

0,00

2 ,0

6,225

6,216

6,216

6,216

0,14

0,00

0,00

0,00

к.н.

4,482

4,545

4,545

4,545

1,40

0,00

0 ,00

0,00

0 ,2

20,62

20,62

20,62

20,62

 

0,00

0,00

0,00

0,00

1,0

14,87

14,82

14,82

14,82

1

0,33

0,00

0 ,0 0

0,00

2 ,0

9,337

9,256

9,256

9,256

0,87

0,00

0,00

0,00

к.н.

3,533

3,706

3,706

3,706

 

4 ,9 0

0,00

0,00

0,00

(

Придерживаясь прежних обозначений, имеем

 

 

 

 

Лц = PiEh.2,

 

 

где

 

 

 

 

 

>

 

 

 

 

 

 

 

 

При решении задачи в первом приближении получим

 

 

а('> =

А11

sin

лх

s in - ^

=

 

Eh?

1—

/

 

 

 

, пх

, пи

4

 

з v

sin —

sin —1-

'

1 у

\

 

 

 

Ц^т)

 

 

 

 

 

 

Во втором приближении

 

 

 

 

 

где

 

Азз = РзEh2,

 

 

 

9 (ai +

а 2)

,

4

648

 

Рз =

 

 

Г у

[ 4

 

Хз + Т х' " 1 Г ¥ з + 1^

 

 

81л2(’+Т>

 

 

 

 

 

 

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

 

 

 

Щ)

,

9я2

, Злх . Зли

 

 

 

 

-А33—-

sin----- sin — -

Ц т + у )

 

 

16

 

sin —

sin -^L —

 

% + —у

a

b

 

 

 

 

 

 

— 9я2ф

______ Г_9 (ах+ а2).,

+

J _ A-2.

81я2 ( ,

J _ Y

L

4

 

5

 

 

 

( Y +

Y

/

 

 

 

 

 

648

 

12x^1 sin

sin

 

---------X 1x3+

 

Здесь

35

 

 

J

a

b

 

 

 

 

 

 

 

 

 

 

ip

Eh?

 

 

 

Подобные формулы получим при решении задачи в третьем и чет­ вертом приближениях.

Аналогично строятся зависимости касательных и интенсивности на­ пряжений от прогиба. На рис. 5.20 приведены значения о/ф, т/ф в раз­ ных точках квадратной панели оболочки с параметром кривизны «1 + 02= 24 при прогибе f=3,lh. Здесь же построены кривые равных на­ пряжений а/ф и т/ф.

Из таблиц функций о+>, т<!'\ согласно которым строили эти кривые, видно, что в самых неблагоприятных случаях отличие значений напря­ жений, полученных в первом и втором приближениях, не столь велико,

как для зависимости нагрузка — прогиб. Разница

между

значениями

напряжений, полученных при решении задачи во

втором

и третьем,

третьем и четвертом приближениях, как видно из таблиц и графиков, незначительна. При меньших прогибах эта разница еще меньше.

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

Рис. 5.20

 

 

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

панели

гибкой

пологой оболочки, работающей в упругой области.

 

 

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

четверти панели с у =1,5; ai + a2= 24 и при разных

прогибах

х = 0,6ft;

ft; 1,25ft; 3ft; 3,1ft приведены на рис. 5.21.

 

 

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

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

Рассмотрим зависимости напряжений а, т и а* от прогибов в разных точках (120) четверти панели (в остальных трех четвертях картина

аналогична). На рис. 5.22 приведены графики таких зависимостей. На­ пример, в точке 15 с координатами х = а/8, г/= а/4 построены кривые за­ висимости а, т, Oi от прогиба этой точки.

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

34 П. М. Огибалов, М. А. Колтунов

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

Большое количество работ по динамической устойчивости панелей гибких оболочек решено методом Бубнова—Галеркина. При этом вво­ дится сила инерции и затухания в дополнение к внешней нагрузке и осуществляются, например, такие же операции, которые приведены здесь с анализом решений. Поэтому, не останавливаясь на этих вопросах додробно, приведем одно из решений по динамической устойчивости по­ логой цилиндрической оболочки, данное Г В. Мишенковым (12].

§ 3. ДИНАМИЧЕСКАЯ УСТОЙЧИВОСТЬ ПАНЕЛИ ПОЛОГОЙ ОБОЛОЧКИ

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

материальной

системы (тела) начинают неограниченно возрастать

во

времени, что

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

а затем и к разрушению.

к

Впервые

вопрос о динамической устойчивости применительно

упругим стержням в 1924 г. рассматривал Н. М. Беляев [13]. Четверть века спустя М. А. Лаврентьев и А. Ю. Ишлинский [14] исследовали яв­ ления, обусловленные ударными воздействиями или внезапным прило­ жением на стержень нагрузок. Они показали, что при внезапном прило­ жении нагрузки, превышающей п-ю критическую статическую силу, возможно появление п-й устойчивой формы равновесия, имеющей' п по­ луволн. Этот результат авторы подтвердили экспериментально. Извест­ ны и зарубежные работы [15, 16], посвященные исследованию устойчи­ вости стержней.

В’ 1938—1939 гг. В. Н. Челомей [17] рассмотрел ряд задач о дина­ мической устойчивости авиационных конструкций. Ряд исследований этого рода теоретического и экспериментального характера выполнили другие авторы [18—20].

В 1955 г. А. И. Блохиной [21] решена задача о динамической устой­ чивости цилиндрической оболочки, шарнирно опертой по краям. Извест­ ны исследования других авторов, советских и зарубежных. Довольно полный обзор современных направлений в области динамики пластин и оболочек и большая библиография даны В. В. Болотиным; интересую­ щихся отсылаем к его статье К Изложим решение, принадлежащее Г В. Мишенкову [12].

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

С учетом приведенных предположений деформацию оболочки опи­ сываем системой нелинейных уравнений:1

1 «Труды II всесоюзной конференции по теории пластин и оболочек». Киев, 1962.