книги / Оболочки и пластины
..pdfПодставляя теперь в выражения (5,1,9) для |
Т и ..., |
N ь N 2 значения |
||
Ох, Оу,..., тХу из (5,1,7), учитывая (5,1,8) и (5,1,1), |
имеем |
|
||
Л_ |
_л |
|
|
|
2 |
2 |
|
|
|
7*1 = |
£г |
|
d z = |
|
1 — V2 V дхг |
дф |
|||
|
) |
|
|
|
Eh |
(ef0 |
+ |
ve° ). |
|
||
|
|
1— v2 |
|
||||||
|
|
|
x x |
‘ |
|
y y } |
|
||
Внося сюда значения |
e \x и & |
из (5,1,1), получим |
|
||||||
Т г = - £ h |
(е° |
+ ve° ) = |
Eh |
|
Г |
|
д* |
+ — ( — |
Y — kjW + |
1 — v2 |
хх |
уу |
1— v2 |
L |
|
2 \ дх |
) |
Аналогично найдем Т2 |
и S = S \ = S 2. |
|
|
|
|||
Таким образом, получаем важные формулы |
|
||||||
• Гди, , |
v dva |
-L |\ |
|
У |
dw у |
||
v2 [ дх |
|
ду |
|
2 V |
|||
|
|
|
f dw |
V _L |
v |
/ ^ |
|
'iS . 4- v * ft+ _L 1 |
|
||||||
, д у |
|
дх |
2 '^ ду J + Т [ дх |
|
|||
s = |
|
£ Л |
Г-ди0 |
dv0 |
, |
dw |
4 = л |
|
2 ( l « * v ду) |
'V дх |
|
дх |
ду ) |
которые дают связь между тангенциальными усилиями, действующими в срединном слое, и его перемещениями. Срезывающие усилия будут
определены ниже.
Исключая из равенств (5,1,12) неизвестные перемещения срединно го слоя и0 и V 0 , можно выразить усилия Т ь Г2, S через прогибы w .
Предварительно заметим, что из (5,1,12) можно получить такие со отношения:
Пользуясь этими |равенствами, составим уравнение
(Тг - vT2) + |
(Г2 - v7\) - 2 (1 + v) |
д хд у
d2w |
i |
d2w |
^ |
~дф |
2 |
дх2 |
Г |
Но согласно правилам дифференцирования
д2 |
\ |
± |
( |
dw |
- V |
I |
J |
д |
Г |
dw |
' |
d2w |
1\ |
|
- |
( |
' |
d2w |
'\ 2 |
, |
dw |
d3w |
||
ду2 |
L |
2 |
ч |
) |
|
ду |
1 |
^ |
дхду |
Г |
|
\ |
. дхду |
.) |
+ |
~дх |
'' дхду2 ' |
|||||||
|
|
|
|
~дх |
|
|
|
|
|
|
J |
|
|
|
||||||||||
д2 |
г ± { . dw • V |
I |
J |
д |
Г |
dw |
|
& w |
'} |
|
- |
( |
’ |
d*w |
'\ 2 |
, |
dw |
d3w |
||||||
дх2 |
L 2 |
|
\ |
ду |
) |
|
дх |
L |
I T |
|
дхду |
И |
|
|
. дхду |
,) |
‘ |
ду |
дх2ду ’ |
|||||
|
|
д2 |
( . dw |
|
dw Л _ |
|
|
f |
dbw |
|
|
dw |
+ |
dw |
|
dPw |
\ |
_ |
||||||
|
|
дхду |
\ |
дх |
|
|
ду Г |
|
дх |
'\ |
дхду |
|
|
ду |
"& Г |
ду2 |
) |
|
||||||
|
|
|
|
d3w |
|
|
|
|
|
|
|
|
|
|
|
dPw |
dw |
|
d3w |
|
||||
|
|
|
дх2ду |
|
|
|
|
|
|
|
|
|
|
|
|
ду2 |
дх |
|
дхду2 ' |
|||||
Внося это |
в предыдущее уравнение, |
получим |
|
|
|
|
|
|
||||||||||||||||
|
|
д2 |
( Т |
г - у Т |
г) |
|
д2 |
(T t |
- v |
T |
J |
- |
2 (1 +v ) |
|
d2S |
|
|
|||||||
|
|
ду2 |
|
|
|
|
|
дх2 |
|
|
|
|
|
|
|
|
|
|
дхду |
|
||||
|
|
= |
|
E h |
|
|
|
|
|
d*w |
|
d2w |
|
|
^ |
|
d2w __^ |
|
d2w *| |
(5,1,14) |
||||
|
|
|
|
|
|
|
|
~дх? |
|
~ду2 |
|
|
|
'~ д у 2 |
|
2 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дх2 J ' |
Уравнение (5,1,14) дает искомую связь между усилиями, действующи
ми в срединной плоскости оболочки, и ее прогибом |
до. |
|
|||||
Введя функцию напряжений ф по формулам |
|
|
|||||
T x = |
дг<р |
|
|
s |
= - h |
- ^ Ц |
(5,1,15) |
h |
|
|
|||||
|
ду2 ’ |
|
|
|
|
дхду |
|
уравнение (5,1,14) после подстановки значений (5,1,15) примет вид |
|||||||
- j V2V2<P+ К |
d2w |
d2w |
|
d2w |
d2w |
= 0, |
(5,1,16) |
~дф + |
^2 ~д& |
|
"дх2 |
~ду* |
|||
|
|
|
|||||
где |
|
|
|
|
|
|
|
|
V2V2(P = |
-тт- + |
2 |
о |
+ |
|
|
|
|
дх4 |
|
дх2ду2 |
|
|
|
V2V2 — бигармонический оператор |
Лапласа. |
|
|
Теперь перейдем к построению уравнения равновесия.
Пусть на оболочку действуют нагрузки: поперечная q = q ( x , у ), дей ствующая нормально к срединной поверхности, сжимающие или растя
гивающие усилия h p ( x ) |
и h r ( у ) , |
приложенные |
нормально к кромкам |
||||||||||
оболочки, и сдвигающие усилия Их, действующие вдоль кромок. |
|||||||||||||
Выделим |
|
из |
оболочки |
двумя |
|
|
|
|
|
|
|
||
парами взаимно |
перпендикулярных |
|
|
/ |
|
|
|
|
|||||
плоскостей, параллельных до де- |
|
|
dN . ) |
. |
Ni dy r d |
||||||||
формации оболочки плоскостям х о г |
|
|
+ |
|
|
|
У |
||||||
и y o z , элемент |
h d x d y . |
|
грани |
|
|
|
|
|
|
|
|||
Усилия, |
действующие на |
(T^ ^ I i d;] dx |
|
|
|
|
|||||||
этого элемента, показаны на рис. 5.4 |
' г |
ду |
|
|
|
|
|
||||||
При изгибе |
срединного |
слоя |
|
d^ d y \ d i ^ ^ ^ r |
|
|
|||||||
эти усилия |
поворачиваются |
в про-' |
|
у |
* |
|
| |
|
|
||||
странстве. Проектируя их на под |
|
|
|
|
|
|
|
||||||
вижные координатные оси и отбра |
|
(T<+ f r dxh |
f a |
l x |
dxh |
||||||||
сывая малые величины более высо |
|
|
|
|
|
|
|
||||||
кого порядка, получим |
|
|
|
|
|
|
|
Рис. 5.4 |
|
||||
|
+ |
dx)dy~Tldy+(Sa+^ d |
i d y ) |
|
d x ~ |
S id x ‘ |
|||||||
— N j d y — . — d x — ( N 1 + - ^ d x ) d y — |
|
d*w |
d x = 0, |
||||||||||
|
|
||||||||||||
|
1 |
2 |
dx2 |
|
\ |
1 |
dx |
J |
2 |
|
dx2 |
|
|
— T 2d x + ^ T 2 + - ^ - d y ^ d x + ( s t + ^ - d x ^ d y - S xd y - |
|||||||||||||
— (^2 + |
~ r d y ^ d x 4 - |
• -45- d y — № |
4 r |
• - Й - dy = °> |
|||||||||
|
|
|
dy |
J |
|
dy* |
|
|
|
|
dy* |
|
+ ( г ’ + ^ |
+ г 0 Ч т § Г + Т ^ ) + |
||||
+ |
H—~ |
d x -f- |
d y — |
d2w |
-I- |
|
|
|
dxdy |
|
|
+ ( s * + ^ d y + |
S2) |
~ S k d y + ( Ni + |
ddx |
+ |
|
*** N l ) dy |
+ ( N 2 + ^ - d y — N ^ d x + q ( x , y ) d x d y = 0.
Проводя в первых двух уравнениях сокращения и отбрасывая чле ны третьего порядка малости, имеем
dTi |
|
dS2 |
= А \ |
|
дх |
~Г |
ду |
||
дх" |
||||
дТ2 |
^ |
|
хг |
|
ду |
|
дх |
— -А2 ду* |
ное состояние оболочки. Полагая в (5,1.16) и .(5,1,24) начальные кривиз ны k\ = 0 и k 2= 0 , получим уравнения для пластин большого прогиба:
v V < p = £ r r — Y - — -— 1. |
|
|||||
V V f |
[v дхду ) |
дх' |
дф |
J' |
|
|
|
|
|
dx°~ |
dy2 |
|
|
Dy2y2tiy — q = h Ldf |
d2w |
, д2ф |
(pW |
|
^2ф |
d2w |
дх2 |
dx2 |
dy2 |
|
дхду |
дхду |
Задача об изгибе и устойчивости пластин и оболочек, как видно из предыдущего, сводится к интегрированию системы совместных нелиней ных дифференциальных уравнений в частных производных четвертого порядка:
1 |
2 |
2 |
\ |
|
U |
|
I |
и |
|
|
I |
дх2 |
|
|
|
|
|
|
|
|
|
|
|
|
Е |
V |
V Y |
-г |
1 |
|
г |
2 |
дхг |
^ |
|
дф |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D y V ^ |
|
h ( k г J |
■ + ь * |
|
) - |
|
|
|
|
|
|
|
|||||||
-ч |
|
|
|
|
|
|
|
|
|
|
|
|
1 ду2 |
|
|
|
|
|
|
|
||||
д2ф |
d^w |
|
. |
d2y |
d2w |
|
2 |
- |
^ . |
д2до )-< 7 = 0. (5,1,25) |
||||||||||||||
|
|
|
% 2 |
дх2 |
|
‘ |
дх2 |
‘ ду2 |
|
|
дхду |
|
дхду |
|
|
|
|
|
|
|||||
Связь между усилиями, действующими в срединном |
слое, и переме |
|||||||||||||||||||||||
щениями определяется формулами (5,4,12): |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
Eh |
Г |
ди0 |
|
до0 |
|
|
/ |
dw ' |
|
V |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
■4- V |
|
|
|
|
|
|
.)‘ + |
|
|
|
|
|
|
|
|
|
|
||
1 — V2 1. |
дх |
|
ду |
|
|
\ |
дх |
т |
|
|
|
|
|
|
|
|
|
|||||||
|
Eh |
Г dvQ |
|
ди0 |
1 |
1 |
/ |
dw ' |
|
V |
|
|
|
|
|
|
|
|
|
|||||
1 — v2 [ |
дх ■+ |
V |
ду |
+ |
т |
V ду |
.)’ + |
~2 |
|
|
|
|
|
|
(5, |
,26) |
||||||||
|
|
|
|
S = |
|
Eh |
|
/ |
ди0 |
|
до, |
|
dw |
|
dw |
\ |
|
|
||||||
|
|
|
|
|
|
|
\ |
ду |
|
|
—2— |
1------- |
|
ду |
) ' |
|
|
|
|
|
||||
|
|
|
|
|
2 (1 ■Ф- v) |
|
|
дх |
|
дх |
|
|
|
|
|
|
||||||||
В случае |
|
сферической |
оболочки |
(k i = k 2= k ) |
уравнения |
(5,1,16) |
и |
|||||||||||||||||
(5,1,25) |
примут |
вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
— v 2V2(P + k\/*w + |
^ |
|
d2w |
/ |
|
|
** |
о , |
|
|
|
|
|||||||||
|
|
|
Е |
|
|
|
|
|
|
|
дх2 |
|
ду2 |
\ |
дхду I |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dPw |
|
о |
д2^ |
|
d2w |
\ |
|
|
|
|
|
|
|
|
|
\ |
ду2 |
|
дх2 |
дх2 |
|
— 2 |
дхду |
дхду |
J - < 7 = |
0. |
|||||||||
|
|
|
|
|
|
ду2 |
|
|
||||||||||||||||
Для |
цилиндрической |
оболочки |
(fti=0, |
fc2# 0 ) |
уравнения |
будут |
||||||||||||||||||
такими: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
- j W |
< P + f c 2 |
дх2 |
\ |
dh) |
d2w — |
( ' |
|
у |
|
|
|
|
|
|||||||||
|
|
|
~д& |
|
|
|
\ . |
дхду |
) |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
V |
|
dPw |
|
ааФ |
|
dhs) |
|
|
|
d2w |
\ |
|
|
|
||
|
|
|
|
|
|
|
|
2 |
|
дх2 ■+ |
■дх* |
|
ду2 |
|
дхду |
|
- 7 |
= |
0 . |
|||||
Если знаки |
кривизн k { |
и k 2 |
разные, |
например, &i>0, k 2< 0 , |
то эти |
|||||||||||||||||||
уравнения |
могут |
быть |
использованы |
в |
качестве |
основных |
уравнений |
|||||||||||||||||
теории гибких оболочек отрицательной гауссовой |
кривизны. |
|
|
|
|
Следовательно, граничные условия для закрепления кромок по схе ме (б) запишутся так:
d2w = |
о, и = |
v = 0. |
дх2 |
|
|
Эти условия говорят о том, что геометрия кромок оболочки при ее |
||
изгибе остается прежней и они не |
смещаются. |
|
В третьем случае (в) закрепления |
кромок граничные условия |
|
будут |
|
|
и- const, V = const,
что выражает факт смещения кромок панели при'деформации оболочки параллельно самим себе. Возможны также комбинации этих способов шарнирного закрепления.
Рис. 5.7
Аналогично можно рассуждать и в случае защемленных кромок (рис. 5.7).
Граничные условия,- соответствующие схеме (с), можно записать
так:
|
w = |
dw |
Л |
|
|
|
---- = |
0, |
|
||
агд> _ |
дх |
|
|
||
т = |
- |
0, |
|||
OV = |
= 0, |
||||
ду2 |
|
|
д)еду |
||
а граничные условия схемы |
(d) |
|
|
|
|
w |
dw |
u = |
v = |
0. |
|
= |
дх
В первом случае (с) при деформации пластин кромки искривляются, во втором (d) — кромки остаются прямолинейными и не смещаются.
Влинейной теории в таком различии необходимости нет, и схемы
(с)и (d) можно считать эквивалентными.
Вслучае, когда на кромках пластины заданы перемещения и и V,
после интегрирования основных уравнений (5,1,16) и (5,1,25) по найден ной функции сриш следует составить общие выражения перемещений и и о, которые следует затем подчинить заданным граничным условиям на кромках. Другие условия здесь не рассматриваем.
§ 2. ПАНЕЛЬ ГИБКОЙ ПОЛОГОЙ ОБОЛОЧКИ
Рассмотрим задачи об изгибе и устойчивости панелей гибких поло гих оболочек (геометрическая нелинейность).
Гибкими называются оболочки, у которых прогиб сравним с их толщиной. Для этого учитывают геометрическую нелинейность, выражая деформации через перемещения в срединной поверхности зависимостя ми (5,1,1), т. е. в ряду для компонентов деформации удерживаются
члены второго порядка малости относительно прогибов w = i ( x , у ) . Зада ча может быть сведена к интегрированию системы двух нелинейных
уравнений |
|
(5,1,16) |
и (5,1,25) |
|
d2w| . |
|
|
|
|
|
|
|
|||
1 |
V |
о |
« Л |
* |
д2Ш . |
, |
d2w |
02W |
|
/ |
dPw |
|
|||
Ф = — |
У |
Ф + |
« 1 |
-------------Ь |
^ 2 -------- |
дх2 |
ду* |
|
\ дхду |
|
|||||
Е |
v |
v |
т |
1 |
ду* |
^ |
2 |
дх* |
|
|
|||||
|
|
|
|
|
|
|
|
( к |
^ |
+ |
k 1 - |
°2Ф |
\ |
|
|
|
|
|
|
|
|
|
|
\ |
дх* |
|
ду* |
) |
|
|
|
|
|
|
<Э2ф |
|
д2w |
, |
о2ф |
0*W |
|
с, |
О2Ф |
|
d*w |
\ |
(5,2,2) |
|
|
|
ду2 |
~~ д х + |
дх2 |
ду* |
|
дхду |
дхду |
= 0 |
|||||
|
|
|
|
) |
|
||||||||||
при соответствующих граничных условиях. |
упругих |
оболочек решают |
|||||||||||||
Большинство задач |
нелинейной теории |
приближенными методами, причем из-за сложности обычно ограничи ваются решением в рядах в первом приближении.
Здесь приведем решения задач об изгибе и устойчивости гибкой пологой оболочки, построенных в первом [3] и в более высоких прибли жениях [4] методом Бубнова — Галеркина.
Решения таких задач представляют интерес с точки зрения возмож ности установить области неустойчивости оболочек. Линейная теория оболочек дает возможность установить только верхнюю границу этой области. Эксперименты показывают, что значения расчетных критиче ских напряжений не совпадают с экспериментальными и значительно их превосходят.
Рассмотрение задачи об изгибе и устойчивости оболочек в свете нелинейной теории дает возможность предвидеть поведение оболочки после потери устойчивости и установить не только верхнюю границу области неустойчивости, соответствующую критическим напряжениям, полученным по линейной теории, но и ее нижнюю границу.
Приведем некоторые соображения общего характера о решении задач этого типа. Пусть к оболочке, каким-то образом закрепленной на жестком кусочно-гладком кон туре, приложены произвольная поперечная нагрузка q(x, y) и нормальные к кромкам сжимающие или растягивающие напряжения, компоненты которых вдоль осей оу и ох будут р(х) и г (у). Требуется определить зависимость между действующими на обо лочку внешними усилиями и ее прогибами, не считая последние малыми. Для решения этой задачи необходимо проинтегрировать уравнения (5,2,1) и (5,2,2). Ввиду того что методы точного интегрирования этих уравнений еще не найдены, будем искать их при ближенные решения в виде рядов
ф = |
А,пп [ип (х) vn (у) - |
е (*) - х т . |
||
т |
п |
fmnXmЕ |
(*) |
|
|
' - Е |
Уп (у), |
тп
где Л шп и fmn — неизвестные постоянные, а функции Um (x), Vn (y)y Xni(*), Yn {y) выбирают заранее таким образом, чтобы удовлетворялись все статические и геометри
ческие контурные условия.
Применяя метод Бубнова — Галеркина, внесем эти выражения для ф и w в урав нения (5,2,1) и (5,2,2), затем, умножая в соответствии с их физическим смыслом пер вое на -вариацию функции ф, а второе на вариацию функции w и учитывая независи мость вариаций параметров 6Amn между собой и вариаций параметров 6fmn между собой, проинтегрируем полученные выражения по области, ограниченной контуром обо лочки. Из полученной системы нелинейных алгебраических уравнений найдем неизвест
ные параметры А тип И /тп-
Возможность применения метода Бубнова — Галеркина к операторам Ф и w изу чалась в работах И. И. Во-ровича [5] и других авторов. Об этом же рассказано в [6], где приводится доказательство, предложенное А. Р. Ржаницыным.