книги / Устойчивость и колебания трехслойных оболочек
..pdfперечных |
сдвигов <p=ai,1+ 02,2 (где хь х2— ортогональные ко |
|
ординаты |
в срединной плоскости заполнителя; ai, |
аг — попереч |
ные сдвиги в плоскости * 0z и yOz соответственно; |
z — расстоя |
ние по нормали срединной поверхности заполнителя).
Э. И. Григолюк (10] вариационным методом установил гра ничные условия и получил систему разрешающих уравнений ко нечных прогибав, произвольно нагретых по толщине и по по верхности упругих пологих трехслойных оболочек с легким за полнителем, когда несущие слои ортотропны в механическом и термическом смысле, а оси их ортотропии совпадают.
Система четырех уравнений, содержащая до, F, |
<*1, аг для изо |
тропных несущих слоев сведена последовательно |
к трем (до, F, |
cf) и Двум (до, F) нелинейным уравнениям. Здесь |
впервые в те |
ории слоистых оболочек была сформулирована гипотеза о ли нейном распределении касательных перемещений по высоте пакета, позволившая методологически строить эту теорию в ду хе теории однослойных оболочек. Принималось, что несущие слои, передающие изгиб, и кручение, испытывают конечные про гибы, а заполнитель воспринимает только малый поперечный сдвиг. Гипотеза Кирхгоффа—Лява о прямой и нерастяжимой нормали несущих слоев и предположение о прямолинейности нормали в заполнителе удовлетворяют принятому линейному закону распределения касательных перемещений по толщине оболочки. Одновременно для случая изотропных несущих слоев дана система двух нелинейных уравнений (до, F ), найденных при условии, что срединные поверхности несущих слоев присоеди нены к крайним поверхностям заполнителя.
Развивая прежнюю работу (10] Э. И. Григолюк выдвинул [11] общую теорию произвольно нагруженных и нагретых поло гих оболочек с ортотропными жесткими несущими слоями и ортотропным жестким заполнителем, сопротивляющимся нагруз кам в плоскости, параллельной срединной поверхности. Несжи маемый в поперечном направлении заполнитель и несущие слои испытывают конечные прогибы; механические характеристики несущих слоев и заполнителя различны, но оси ортотропии их параллельны. Полученная в общем случае система четырех не линейных уравнений (до, F, ai, аг) для изотропных слоев с оди наковым коэффициентом Пуассона сводится к системе трех уравнений (до, F,<р). В [10, 11] дается также уравнение прогибов до, описывающее устойчивость пологих оболочек при малых прогибах.
Случай трехслойной оболочки несимметричной структуры может быть изучен на основе компактной системы уравнений.
Для случая оболочки с изотропными несущими слоями и не сжимаемым трансверсально изотропным жестким заполните лем 1в постановке [11] систему пяти нелинейных уравнений мож но свести [15— 14] к трем дифференциальным уравнениям, со держащим функцию прогиба %, силовую функцию F и функцию
71
углов поворота <р, выражая прогиб w и углы поворота нормалЕ к срединной поверхности заполнителя (xi, хг) в виде
W= (1_T уа) х:
(2. 122
где
“i = * i + «\i; a, = va+ w .„ |
(2. 123, |
V2 = d2{)/dx12-{-dfi()/dx22 — оператор Лапласа; h — суммарная тол щина пакета; р — коэффициент поперечного сдвига; у, Ф — коэф фициенты, зависящие от толщин и механических характеристи оболочки. При этом одно уравнение оказывается независимы! (1 —fl)h2V 2cp=2pcp (v — коэффициенты Пуассона).
Система дифференциальных уравнений с %, F, <р имеет общи> порядок, равный двенадцати, одно (ф) из уравнений этой cnq темы не связано с остальными и при решении частных задач он» может не приниматься во внимание. Тогда задача сводится Е решению системы двух нелинейных уравнений (%, F ), обща; структура которой весьма напоминает соответствующие урав нения теории конечных прогибов однослойных оболочек Мар герра {27]. Более того, эта система является разрешающей — че рез основные функции %> F выражаются все перемещения и уси лия и, следовательно, граничные условия. Последним достоин сивом многие редуцированные системы, отмеченные выше, ш обладают, и для решения на их основе конкретных задач при ходится вновь возвращаться к более громоздкой системе диф ференциальных уравнений.
Г л а в а 3
УСТОЙЧИВОСТЬ И КОЛЕБАНИЯ КРУГОВЫХ ЦИЛИНДРИЧЕСКИ) ОБОЛОЧЕК
1. ЛИНЕАРИЗИРОВАННЫЕ УРАВНЕНИЯ УСТОЙЧИВОСТИ
Относя круговую цилиндрическую оболочку‘радиуса R, дли ны I к координатам х\— х, хг= s, измеренным соответственно па образующей и дуге круга нормального поперечного сечения, со гласно (2.99)— (2.100) приходим к следующим линеаризирован^
.72
ным уравнениям местной потери устойчивости цилиндрических трехслойных оболочек (рис. 13):
|
v w |
= |
~ |
— (1 - — v2) z; |
(3.1) |
||||
|
|
|
R |
dx2 |
\ |
p |
уЛ |
|
|
D I 1 — — у ) v2v2x + — |
^ — |
[ № n ~ + |
|
||||||
|
|
) |
|
R |
dx2 |
\ |
и дх?~ |
|
|
+ |
2ЛГ° — |
— I-JVS2— 1 (1 |
8 |
V2) x = 0 . |
(3.2) |
||||
~ |
12 dxds ~ |
22dsi! |
I |
J |
Л |
|
|||
Здесь |
|
|
|
|
|
|
|
|
|
|
|
|
|
dx2 |
1 ds2 |
|
|
|
Рис. 13. Круговая цилиндрическая оболочка
Для свободно опертой замкнутой цилиндрической оболочки граничные условия имеют вид
77= V 2^;’= x = V 2x := V 2V2x — 0 при JC= 0 и х = 1 . (3.3)
Систему уравнений (3.1) — (3.2) с помощью введения разреша ющей функции xi
Х = v ^ , |
£ ( 1 - у V2) Хх |
(3.4) |
сведем к одному уравнению устойчивости
° ( 1 - у v‘ ) ^ ^ + 5 ? й ( 1 - f '»■) *‘ ~
- ( П £ + 2 Л ^ + Л ^ ) ( 1 - f v j v v n - o . (3 .5 ,
7 3
Граничные условия для свободно опертых краев относительно функции Xt запишутся следующим образом:
X i = V * X i = V * V * X i = V aV2V2X i = V 2V2V *X x = 0 при х — 0 и х = 1 .
(3 .6)
Перейдем к решению конкретных задач по определению крити ческой нагрузки свободно опертой круговой тонкой упругой трехслойной цилиндрической оболочки при различных внешних воздействиях.
2. КРИТИЧЕСКАЯ НАГРУЗКА ПРИ РАВНОМЕРНОМ ОСЕВОМ СЖАТИИ
Рассмотрим [14] местную потерю устойчивости - опертой по торцам круговой трехслойной цилиндрической оболочки, под верженной действию равномерно распределенной по контуру сжимающей силы N. Под действие этого усилия в оболочке при докритическом безмоментном состоянии возникнут удельные усилия, равные
№ = - |
N -; ^ = о ; N O 2= о. |
(3 .7 ) |
|
2яЛ |
|
Задача сводится к |
решению уравнения устойчивости |
|
|
|
|
(3 .8) |
которое ищем в виде |
|
|
|
|
. тпх |
ns |
(3 .9) |
X i= X osln —— cos — |
|||
где m — число полуволн |
по образующей |
цилиндра, л — число |
|
волн по окружности, образующихся |
в момент потери устойчи |
||
вости оболочек; хо— постоянная. |
|
|
|
Вводя (3. 9) в (3. 8), |
придем к формуле для критического |
||
осевого усилия |
|
|
|
74
Подбирая целоисчисленные значения т, п, соответствую щие минимуму выражения (3. 10), найдем наименьшее значение критической осевой силы
|
|
|
-^Kpmin= |
2я/?(ЛГ°1)крт1п= |
— |
РтЫ- |
(3.11) |
|||
|
Проведенные вычисления показывают, что при k^ O , |
Х .^ 0,4, |
||||||||
р-* > |
16 /^практически не зависит от параметра % (см. таблицу). |
|||||||||
|
Таблица значений pmin* при |
р2= 100, |
0 = 0,05, |
/г = 1, п —0 |
||||||
|
X |
1,00 |
1,1 |
1,2 |
1,3 |
|
1,4 |
|
1,5 |
1,6 |
. |
• |
5,4218 |
5,4212 |
5,4013 |
5,4134 |
5,4136 |
|
5,4013 |
5,4093 |
|
Рт\п |
|
|||||||||
т |
7 |
7 |
8 |
9 |
|
9 |
|
10 |
11 |
|
На рис. 14 приведены графики |
зависимости |
наименьшего |
||||||||
критического |
значения |
параметра |
Р* = |
Р*т1в, |
определяемого из |
|||||
формулы |
|
|
I |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( 3 .12) |
от параметра |
сдвига k |
для значений р2= 64, |
1000, 5000, |
10000, |
||||||
15000,30000 при 0 = 0, 1. |
|
|
|
|
|
|
С достаточно хорошим приближением эту зависимость мож но получить, считая т21№+п21п2 и т 2Д 2 непрерывными аргумен тами и проводя минимизацию выражения для р*.
Полагая |
|
|
|
перепишем формулу (3. 12) в виде |
|
||
Р* |
1 -j- ftkmi |
(3.13) |
|
1 Н- km\ |
|||
|
|
||
(^ -< m 1< o °, m1< m 8< ooj . |
|
||
Легко видеть, что минимум выражения (3.13) |
достигается |
||
на границе изменения переменных, определяемой |
равенством |
||
|
гп{—т2, |
(3.14) |
75
что возможно, когда форма потери устойчивости оболочки осесемметрична. Выражение для р* приобретает вид
|
|
Р* |
1 -}- Ыгтпл |
. |
и.2 |
|
|
|
|
(3.15) |
|
|
|
---------- Lml + |
— • |
|
|
|
|
||||
|
|
|
1 4 * ktn.\ |
|
Ш\ |
|
|
|
|
|
|
Рассмотрим два возможных случая. |
|
|
|
|
|
|
|||||
1. Пусть |
/гт! — 1, |
где |
т , |
соответствует |
/>*1п. |
Здесь |
член |
||||
1 и им можно |
пренебречь, |
вследствие |
чего формула |
||||||||
|
|
|
(3. 15) запишется так: |
|
|
|
|||||
|
|
|
|
р *_ ^ |
ml |
I |
!*2 |
_ |
|
(3.16) |
|
|
|
|
|
|
1 + k mi |
mi |
|
|
|
||
|
|
|
Из условия d p*/dm l= 0 |
имеем |
|
||||||
|
|
|
m |
t* |
; |
/СП= М 2 - ^ ) . |
|
(3.17) |
|||
|
|
|
1 — |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 14. Зависимость критического значения |
||||||||
|
|
|
продольно |
сжимающего |
напряжения р* = |
||||||
|
|
|
= р min круговой цилиндрической |
оболочки |
|||||||
|
|
|
от |
параметра |
сдвига |
k |
при 0= 0,1 |
для ря |
|||
|
|
|
|
|
|
да значений |
р,2 |
|
|
||
Формула |
(3.17) дает |
достаточно |
точные |
значения |
P*mln> |
||||||
если выполняется условие |
|
|
|
|
|
|
|
|
|
||
|
|
|
1 -р А > & . |
|
|
|
|
|
|
(3.18) |
|
В частности, при k = 0 |
получим известную формулу для кри |
тического значения параметра продольного сжимающего усилия однородной цилиндрической оболочки.
2 . Когда •bkm.i~\, где rrii соответствует |
Ф°Р~ |
|||
мула для р* приобретает вид |
|
|
|
|
|
1 + ИЙ/И] |
|
(3. |
19) |
|
mi |
|
||
|
|
|
|
|
Определяя корень уравнения d p*ldm i= 0 , |
найдем |
|
||
т г |
_м_ |
|
(3 .20) |
|
p*m\n— v-V* |
k |
|||
|
|
|
|
Уравнение (3. 19) с достаточной для практических расчетов, точностью может быть использовано при ^> 0,2.
Согласно (3.20) и (3. 11) при &>0,2 имеем
11/кр min" |
~ ~ |
\ |
f — ~~----- b Oha— |
(3.21) |
|
2R |
V |
3 (1 — v2) 1 , 3 eT |
|
76
Отсюда для оболочек с легким заполнителем получим
£iAi2 |
£ 2Л22 |
кр щ1п' |
2R / 3 ( 1 — V22) |
2R ^ 3 (1 — Vl2) |
|
|
•‘ч |
и |
т |
т |
(3.22) |
|
|
|
|
|
|
|
|
|
Пренебрегая |
изгибной жесткостью несущих слоев, найдем |
||||||
|
(^И )крга.л~ -0Л8[ 1 + 4 _ ( _!4 Г ^ ) 2] • |
(3-23> |
|||||
Отбрасывая |
в скобке |
второй |
член, |
|
|||
придем к общеизвестной формуле |
|
|
|
||||
|
(-^ц)кртт~ |
Oh3, |
|
(3.24) |
|
||
полученной в предположении, что средин |
|
||||||
ные поверхности несущих слоев прикреп |
|
||||||
лены к поверхностям заполнителя. |
|
|
|
||||
На рис. 15 представлены зависимости |
|
||||||
P* = Pmia |
°т десятичного логарифма па |
|
|||||
раметра |
ц2 при |
k = 0 , |
0,01, |
0,05, 0,1, |
0,2, |
|
|
0,5, 1000, |
когда 0= 0,05.3* |
|
|
|
|
Рис. 15. Зависимость критического значения про дольного сжимающего напряжения р* круговой цилиндрической оболочки от десятичного лога рифма параметра |х2 при 0= 0,05 и ряда значений параметра сдвига k
3.НЕРАВНОМЕРНОЕ ОСЕВОЕ СЖАТИЕ
Вслучае, когда осевая сжимающая сила N приложена экс центрично к торцу цилиндрической оболочки, удельное усилие
при докритическом безмоментном состоянии оболочки вычисля ется по формуле
№ = - ( — + - ^ - c o s — ) , |
(3.25) |
гдеМдзг— изгибающий момент, равный
Mw3r= N e —NRs |
(3.26) |
(е — эксцентриситет силы N ).
77
Решение уравнения устойчивости
° ( 1 - |
т |
i |
7w |
, |
’z‘ + U |
s r ( |
' - |
f Л |
«*+ |
|
|
+ i ^ ( 1+ 2 i c° s i ) S ( 1 - T vV ' ' * ^ 0 |
(3.27) |
||||||||||
|
|||||||||||
ищем в форме |
|
|
|
|
со |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
®=sin |
|
^ |
a ne ins'R ; |
|
|
|
|||
|
|
X i — sin |
^ |
a ~ l e inslR , |
|
|
(3.28) |
||||
где |
|
|
|
|
|
|
|
|
|
|
|
Вводя обозначения |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
h W . 8_ 12/?2(1 — y2) . |
|
(3.29) |
||||
|
|
2Dix3 |
|
P/?2 ’ ^ ~ |
Я4Л20 |
|
|||||
|
|
|
|
|
|||||||
|
|
|
|
i |
l |
|
e |
|
|
|
|
|
|
|
|
к — — : s = |
— |
|
|
|
|
||
и подставляя |
(3.28) |
в уравнение (3 .27), |
найдем |
связь |
между |
||||||
коэффициентами ряда (3. 28) |
|
|
|
|
|
|
|
||||
a n+iP*&— (Фя— P *)a n+ a n- iP * s = 0 |
при п ф |
О |
(3.30) |
||||||||
|
|
ф0— ^1 -]-2е — j |
р * |
При П = |
0. |
|
(3.31) |
||||
Здесь |
|
|
|
|
|
|
|
|
|
|
|
, |
|
„ / /к2 |
п2 \ |
Infl |
л2 \ |
|
/и2 |
|
|
||
t . ~ |
|
|
|
|
|
|
I |
г» |
X2 |
|
(3. 32) |
, |
( т2 |
|
n2\ |
|
/к2 |
//и2 |
л2 ’ |
||||
|
1 |
+ k l ----+ ---- |
|
X2 |
|
|
|
|
|
||
|
|
\ Х2 ^ |
я2} |
|
|
|
|
|
|
||
Преобразуя (3. 30) |
|
|
|
|
|
|
|
|
|
||
|
|
«л- 1 |
|
Фл — Р*— ЕР* «я *-1 |
|
|
|
|
|||
получим, что для достаточно больших п |
|
|
|
|
|||||||
|
|
|
|
«л |
_ |
Р*Е |
|
|
|
(3.33) |
|
|
|
|
|
«л- 1 |
|
фл— Р* |
|
|
|
||
|
|
|
|
|
|
|
|
|
78
Возвращаясь к Оо, будем иметь бесконечную |
цепную дробь |
||||
£ L |
|
р2 |
|
|
(3.34) |
«о |
|
р*е |
|
|
|
<h — Р* |
|
л |
|
|
|
|
|
р*£ |
|
|
|
|
h —p * — Фз—Р* |
|
|
|
|
Теперь на основании (3.31) |
|
имеет уравнение I B |
виде |
бесконеч |
|
ной цепной дроби |
|
|
|
|
|
Фо — Р* |
|
р*£2 |
|
|
(3.35) |
2р* |
|
р*е? |
|
|
|
|
|
|
|
||
h — р * |
|
р*г |
|
|
|
|
|
Ф2— Р* |
е2 |
|
|
|
|
Фз~Р* |
|
|
|
|
|
|
|
|
|
из которого можно определить как р * при заданном |
е, так и в |
||||
при заданном р*. |
|
|
|
|
|
Расчеты показывают, что р* с высокой точностью и всегда с приближением снизу может быть определено по формуле
NR |
Фо |
2лЮ Р |
1 + 2е ’ |
т. е. критическая сила, приложенная вне центра, будет равна или больше критической центральной силы, деленной на вели чину 1 + 2 е, где е — эксцентриситет приложения силы.
4.КРИТИЧЕСКАЯ НАГРУЗКА ПРИ ВНЕШНЕМ РАВНОМЕРНОМ
ПОПЕРЕЧНОМ ДАВЛЕНИИ
Исследуем потерю устойчивости опертой по торцам тонкой упругой круговой пологой трехслойной цилиндрической оболоч ки, подверженной действию внешнего равномерного поперечно го давления q [14, 15]. Под действием этой нагрузки в оболочке при докритическом безмоментном состоянии .возникнут удель ные усилия, равные
№и = 0 ; |
N ^ - q R |
; №п = 0 . |
(3.36) |
Задача сводится к решению уравнения устойчивости |
|
||
° ( ' |
|
£ 0 — 1г v j b + |
|
|
|
|
(3.37) |
Отыскиваем решение уравнения (3. 37) в виде |
|
||
Ъ. = |
. тех |
ns |
/О 00\ |
Хо sin — |
COS — . |
(3.38) |
Здесь п — число волн в окружном направлении; хо— постоянная.
79
Подставляя |
(3. 38) |
в |
(3.37), |
приходим к уравнению |
|
|
|||||||
|
1 -м » |
|
|
|
|
п2_у |
|
|
|
|
|
|
|
д* = |
|
|
|
|
Я 2 ) |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
rfi |
Х« |
* ( a L |
+ |
± |
) 2 |
* |
||
|
1 - м |
|
|
|
|||||||||
|
|
|
|
Я2 |
|
я \я2 |
|
X2 ) |
|
||||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
(3.39) |
|
где |
|
д* |
12^ |
3(1 - У 2) |
. |
а. |
12^2(1 —у2) . |
|
|
|
|
||
|
|
|
|
(3.40) |
|||||||||
|
|
4 |
ЯЛЗДя2 |
|
|
Я40Л2 |
|
|
|
|
|
||
|
|
Я2А2 _ |
I |
|
|
|
|
|
|
|
|
|
|
|
k — - |
’ |
|
|
|
|
|
|
|
|
|
|
|
|
|
т |
R |
|
|
|
|
|
|
|
|
|
|
На рис. 16 приведены графики |
|
|
|
|
|
|
|||||||
минимального значения критическо |
|
|
|
|
|
|
|||||||
го давления |
q |
— 9,min |
в зависимо |
|
|
|
|
|
|
||||
сти от коэффициента %для значений |
|
|
|
|
|
|
|||||||
параметров при ц2= 15000 и ■0=0,05, |
|
|
|
|
|
|
|||||||
когда 6 = 0 ; 0,04; 0,08; 0,2; 0,8; 2; 1000. |
|
|
|
|
|
|
|||||||
Рис. 16. Зависимость критического значения |
|
|
|
|
|
|
|||||||
внешнего |
равномерного поперечного давле |
|
|
|
|
|
|
||||||
ния |
|
круговой |
цилиндрической |
|
|
|
|
|
|
||||
оболочки |
от |
отношения X=l!R (/ — длина, |
0 1 2 |
3 9 5 |
6 |
7 8 |
9 |
7 . |
|||||
R — радиус) |
при |
р2= Г5 000, |
6= 0,05 |
k |
для |
||||||||
ряда |
значений параметра |
сдвига |
|
|
|
|
|
|
|
5. КРИТИЧЕСКАЯ НАГРУЗКА ПРИ ВНЕШНЕМ РАВНОМЕРНОМ ВСЕСТОРОННЕМ ДАВЛЕНИИ
Перейдем к изучению потери устойчивости опертой по тор цам тонкой упругой круговой пологой трехслойной цилиндри ческой оболочки, подверженной действию внешнего равномерно го всестороннего давления (14, 15]. Вследствие этой нагрузки в оболочке при докритическом безмоментном состоянии возникнут удельные усилия, равные
N\x= - \ q R \ №22= - q R ; |
0. |
(3.41) |
|
Задача сводится к решению уравнения устойчивости |
|
||
° (1 — т - v ‘ ) v ’ v W * - + - £ £ |
( 1—f ^ * + |
|
|
+**(т2 £dx2+1£ds2)(■ - f v ) V |
V x , = . 0, |
(3. 42) |
|
которое имеет решение в форме |
|
|
|
пх |
ns |
|
(3.43) |
Хх=Хо sin — cos — |
|
80