134 |
Устойчивость оболочек |
Дополняя соотношения (2) нелинейными членами, получим сле дующие выражения для деформаций в срединной поверхности:
|
|
_ |
ди |
до |
дао |
дао |
|
|
|
|
|
^~~ |
ду |
' дх |
дх |
|
ду |
|
|
Уравнение совместности деформаций принимает вид |
|||||||||
дЧх |
, дЧу |
|
д2у |
|
1 |
|
|
1 |
д2ао |
ду2 |
дх2 |
дх ду |
|
~ г |
М в * |
|
(31) |
||
|
|
|
|||||||
где 1 (и;, ао) — оператор, |
|
|
|
|
|
|
|
||
|
С {ао, но) — 2 Г д2ао |
д2ао |
/ |
д2ю |
\ 2~| |
(32) |
|||
|
|
|
[ |
дх2 * |
ду1 |
\ |
дхду |
) ] ‘ |
|
Изменения кривизн и кручение срединной поверхности определяют по формулам (3).
Соотношения Гука (8) остаются прежними. Первые два уравнения равновесия отвечают соотношениям (5), третье уравнение равновесия
принимает форму |
|
|
|
|
|
дОх , |
д()у . |
, д2ао |
{ 1 |
д2ао \ |
|
~дГ + |
+ |
Ь ё + |
а* \~ К |
+ ~ д ^ ) + |
|
|
+ |
2 т й - Ц |
г + |
‘' = ° - |
(33) |
|
|
Для поперечных сил справедливы соотношения (10). Подставим эти выражения в уравнение (33), а соотношения (14) — в уравнение (31),
тогда, вводя функцию напряжений по формулам |
(16), уравнения (33) |
и (31) приводим к следующему виду: |
|
= ^ (® . |
(34> |
= - 4 - 1. |
(35) |
Оператор 1 в применении к функциям ао, Ф:
|
Устойчивость оболочек в пределах упругости |
135 |
|
||||||||||
Если оболочка до нагружения имеет начальные прогибы щ (х, у), |
|
||||||||||||
то выражения для деформаций получают вид |
|
|
|
|
|
|
|||||||
|
|
|
ди . \ ( дш у |
1_ / дщ \ 2 |
|
|
|||||||
|
*■* |
_ |
дх + |
2 \ |
дх |
) |
2 \ |
дх |
) |
• |
|
|
|
|
до |
1 |
|
. |
1 |
/ дш \2 |
|
1 |
/ |
дщ \2 |
( |
||
|
= - — |
|
|
ди) |
дни |
дщ |
|
|
|
(37) |
|||
|
_ |
|
|
|
дщ |
|
|
|
|||||
|
^ |
ду "*™дх |
*’ "дх™ ’ |
ду |
дх |
|
|
ду |
|
|
|
||
здесь |
ш — полный |
прогиб. |
Повторяя |
вывод |
|
основных |
уравнен |
|
|||||
вместо |
выражений |
(34) и |
(35) |
получим |
|
|
|
|
|
|
|||
|
— |
(ге — а>о)=Ц |
|
. |
1 |
д2Ф |
(38) |
|
|||||
|
|
ф) +!«?■ |
; |
|
|||||||||
-^Г V 4® = — - у - (*• (“ >, ю) — |
(а»«, “ >„)] — |
|
|
|
(ш — |
а»0). (39) |
|
||||||
При интегрировании приведенных выше линейных или нелинейных уравнении необходимо удовлетворить граничным условиям. Для тор цовых сечений замкнутой оболочки эти условия формулируются так: при шарнирном опирании оболочки по краям х = 0, х = /, (направле ние координатных осей соответствует рис. 3) для точек краев должно выполняться условие
при защемлении оболочки по краям
ш = 0, ^дх= 0.
Приведем условия, касающиеся перемещений и, о, а также усилий в срединной поверхности. Если точки краев свободно смещаются вдоль образующей и по дуге, то в этих точках должно быть
ох = 0; т = 0.
В случае несмещающнхся кромок следует положить и = 0; о = 0.
Замкнутые круговые оболочки
Сжатие замкнутой оболочки вдоль образующей. Рассмотрим замкну тую круговую цилиндрическую оболочку длиной I , шарнирно опертую * по торцам, подвергающуюся сжатию вдоль образующей усилиями р, равномерно распределенными вдоль дуговых кромок (рис. 4).
Об о з н а ч е н и я : К — радиус срединной поверхности оболочки;
к— толщина оболочки. Исследуя устойчивость оболочки в малом, определяем верхнее критическое напряжение; при этом исходным
• Обозначения на схемах закрепления краев см. рнс. 2, гл. 2.
136 |
Устойчивость оболочек |
является дифференциальное уравнение (21) Для рассматриваемого случая оно принимает вид
О а . Е |
д*ха |
. / |
Я3 |
дхл |
(40) |
|
Приведем первый вариант решения, в котором предполагаем, что поверхность оболочки после выпучивания является осесимметричной, т. е. что поперечные сечения остаются круговыми. В этом случае про гиб ш будет зависеть только от х\ уравнение (40) переходит в следующее:
|
|
Л Лх» + Р |
йх* |
|
Е |
~дхг |
= |
0. |
(41) |
|
|
+ |
Ж |
||||||
|
|
В соответствии |
с граничными условия |
||||||
|
|
ми принимаем выражение для |
прогиба |
||||||
|
|
|
, |
. |
тпх |
, |
|
|
... . |
|
|
ш = / $ т — |
|
|
(41а) |
||||
|
|
где т —число полуволн изогнутой поверх |
|||||||
|
|
ности по образующей оболочки. |
в урав |
||||||
|
|
Подставляя |
выражение (41а) |
||||||
|
|
нение (41), находим |
[1] |
|
|
|
|
||
|
|
Р |
2 - |
|
I |
А |
|
|
(42) |
|
|
Л |
’ |
№ + |
X2 |
’ |
|
||
|
|
|
|
|
|||||
* |
тпК |
я/? „ |
|
|
|
|
|
|
. |
здесь А = |
— |
—. Приравниваем нулю производную от р |
по Л; |
||||||
а1.Х
при этом считаем т > 1. Получаем следующее выражение для
К = У \ 2 (1 - V 3) |
(42а) |
Подстановка выражения (42а) в формулу (42) приводит к следующему значению верхнего критического напряжения рв:
|
1 |
с |
Л |
(43) |
|
Ра |
]Л }(1— V2) |
Я |
|||
|
|||||
При V = 0,3 имеем |
|
|
|
|
|
|
р, = 0,605Е А . |
|
(44) |
||
Из выражения (42а) |
получаем |
длину |
полуволны |
||
1х = |
- 5 ------- --------- У М . |
(45) |
|||
|
V 12(1 — V3) |
|
|
||
В случае весьма короткой оболочки, если |
С 1) , нужно в фор- |
||||
муле (42) принять т = 1 и пренебречь вторым членом. Тогда
лЮ
Устойчивость оболочек в пределах упругости |
137 |
Рассмотрим вариант решения задачи в допущении, что поверхность оболочки после выпучивания не является осесимметричной. Для про гиба и) принимают следующее выражение, удовлетворяющее гранич ным условиям:
. . тях |
. |
пи |
(47) |
1В)= [ 51П — |
8Ш |
: |
где т — число полуволн по образующей оболочки; п — число полных волн вдоль окружности. Исходим из дифференциального уравнения (40); подставляя в это уравнение выражение (47), получим (1]
|
И ( |
т ая 2 |
п2 у |
|
|
|
Е |
гп*л4 |
|
|
|||
~ П \ |
~ |
+ 1 ? ) |
|
+ |
|
|
|
|
|
|
|||
|
|
/ т2л2 |
|
п2 \ - |
|
т2л2 |
|
(48) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Введем безразмерные параметры по формулам |
|
|
|||||||||||
|
РК . |
|
|
|
|
|
|
|
|
Л2Л |
(49) |
||
Р = ЕК ’ |
|
|
|
л/, |
|
» |
|
Я |
|||||
|
|
|
|
|
|
||||||||
Параметры *0* |
и г| |
выражаются |
через |
длины полуволн |
изогнутой |
||||||||
|
дуги |
^ |
|
лЯ \ |
|
и |
|
по |
образующей |
| |
|||
по формулам |
|
|
> - |
п |
) |
|
|
|
|
|
|
||
|
|
1ц . |
|
|
|
|
|
|
|
|
|
||
|
|
■&= |
|
Л |
|
|
|
|
|
(50) |
|||
|
|
|
1х ' |
|
|
~ |
|
Ч |
' |
|
|
||
|
|
|
|
|
|
|
|
|
|
||||
Параметр б1 характеризует |
очертания |
вмятины; |
т| — длину полу |
||||||||||
волны 1У. Тогда вместо формул (48) можно записать |
|
|
|||||||||||
Р ~ |
|
1 |
|
(1 + |
Ф2)2 |
|
Л + |
б 2 |
(51) |
||||
|
|
|
(1 + |
■&*)* V |
|||||||||
12(1— -V2) |
|
№ |
|
|
|
|
|||||||
Считая числа т и п достаточно |
большими, находим |
минимум р |
|||||||||||
из условия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
др_ |
|
0, |
|
|
|
|
|
|
|
|
|
|
|
др |
|
|
|
|
|
|
|
||
где |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
(1 + |
б 2)3 |
|
|
|
|
|
||||
|
|
|
Р = |
Л- |
|
|
|
||||||
|
|
|
|
№ |
|
|
|
|
|
||||
Из этого условия находим значение р: |
|
|
|
|
|
||||||||
|
|
|
р = |
К 1 2 (1 — V2) |
|
|
|
||||||
и верхнее критическое значение |
параметра р |
|
|
||||||||||
|
Ре = |
|
|
1— |
■■ |
|
0,605. |
|
(52) |
||||
|
^3 (1 — V2) |
|
|
|
|
|
|
||||||
138 |
Устойчивость оболочек |
Приведенным решением не устанавливается однозначно форма волнообразования оболочки; вытекает лишь следующее условие, которому должны удовлетворять величины т| и д :
(53)
В случае выпучивания оболочки с образованием квадратных волн (д = 1) получаем г\ = 0,825; тогда
(54)
Верхнее критическое напряжение рв, определяемое в соответствии с формулой (52), в точности совпадает с формулой (43). Следовательно, потеря устойчивости оболочки в малом с образованием вмятин, распо ложенных в шахматном порядке, происходит при том же напряжении, что и в случае осесимметричного выпучивания.
Все приведенные выше формулы относятся к случаю п ^ 4. Прак тически это соответствует выпучиванию оболочки средней длины. Границы применимости теории оболочек средней длины в случае сжатия определяют, исходя из следующего условия, уточненного по сравнению с выражением (1):
(55)
Выражением (55) определяется область значений I для которых
можно считать справедливой формулу (43).
Рассмотрим случай |
слабо выраженного волнообразования (л = 2, |
п = 3), когда оболочка |
выпучивается с образованием длинных волн. |
При получении расчетных формул для этого случая используют урав нение (29). В выражении (28) следует положить
Принимая
где
из уравнения (29) получим [1]
Устойчивость оболочек в пределах упругости |
139 |
|
здесь обозначено е = |
и введен прежний параметр к = |
. в |
Из условия минимизации р по к находим к2 — е1/2л2 (п2— 1)/(1 — у2)1/2 и величину
|
Л _______ I_____ |
|
п8 — 1 |
(57) |
|
Р" ~ УЗ (1 — у‘) ' |
"2 + 1 |
||
|
|
|||
При V = 0,3 |
будет |
|
|
|
|
Р .* 0,605 |
~ |
| . |
(58) |
В случае п = |
2 находим рв = 0,363, что составляет 0,6 |
от «класси |
||
ческого» значения, определяемого по формуле (43). При /1 = 3 верхнее критическое напряжение составляет 0,8 от «классического» значения.
При |
п2 > 1 получаем |
значение ра 0,605, что совпадает с результа |
том |
по формуле (43). |
|
Следовательно, для весьма длинных оболочек получаются понижен ные, по сравнению с формулой (43), значения критического напряжения.
Характер выпучивания реальных оболочек средней длины не со ответствует ни одному из этих вариантов. В действительности, вместо вмятин прямоугольного очертания, расположенных в шахматном по рядке и обращенных к центру и от центра кривизны, образуются ромбо видные вмятины, глубина которых уже в первоначальный момент
сравнима с |
толщиной оболочки. Эти вмятины появляются обычно |
в процессе |
резко выраженного хлопка оболочки. Отсюда вытекает |
необходимость решения задачи с позиции нелинейной теории. |
|
В книге |
[1] изложено несколько вариантов решения нелинейной |
задачи. Решение по методу Ритца состоит в выборе аппроксимирующего выражения для прогиба до, содержащего несколько варьируемых пара метров, и подстановке этого выражения в уравнение (39). В результате
интегрирования этого уравнения определяют функцию |
напряжений |
в срединной поверхности Ф. Находят полную энергию системы |
|
Э = 11е + 1)и - Ш , |
(59) |
где 1!с — потенциальная энергия деформации срединной поверхности;
Уи — потенциальная |
энергия изгиба; 117— работа внешних |
сил. Ве |
|
личины Ци и Ус определяют по формулам |
|
||
11с = |
| |
[(V2®)2 - (1 + V) Ц Ф , Ф)1 йх Лу, |
(60) |
|
| |
|(у 2о03 — (1 — V) Ц ш . ш)| (1х <1у, |
(61) |
|
Р |
|
|
где р = 2я/?/, — площадь поверхности оболочки. Далее полную энер гию системы варьируют по параметрам прогиба и находят диаграмму равновесных форм оболочки. В процессе выпучивания оболочки число н размеры вмятин являются переменными, поэтому диаграмма равно весных форм представляет собой огибающую серии кривых, отвечающих