книги / Строительная механика, динамика и устойчивость композитных конструкций
..pdfсилы на ось z:
|
|
|
|
Q H |
d d |
Q H |
d d |
(4.1.56) |
|||||
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
||||
qH d H d N H d d 2 N H d d 1 0, |
|
||||||||||||
1 |
|
|
|
Q H |
|
Q H |
k1N k2 N q. |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
(4.1.57) |
||
H |
H |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Составим уравнения равновесия для моментов относительно осей α, β, z. Уравнение моментов относительно оси α:
|
|
H M |
d d |
H M |
d d |
(4.1.58) |
||||
|
|
|
|
|||||||
|
|
|
|
|
|
|
||||
H M d d 1 M H d d 2 Q H d H d 0, |
|
|||||||||
|
H M |
|
H M |
M |
H |
Q H H 0. |
(4.1.59) |
|||
|
|
|
|
|
||||||
|
|
|
|
|
|
|
Уравнение моментов относительно оси β аналогично, а относительно оси z
|
S S k1M k2 M 0. |
(4.1.60) |
|
6. Граничные условия: |
|
|
|
1) |
свободный край: |
|
|
|
N S M H Q 0, |
(4.1.61) |
|
2) |
шарнирно-закреплённый край: |
|
|
|
u v w M 2 |
0, |
(4.1.62) |
3) |
заделка: |
|
|
|
u v w 1 2 |
0. |
(4.1.63) |
161
Это однородные граничные условия, для неоднородных условий на краю заданы какие-либо ненулевые функции. Для замкнутых оболочек часть граничных условий может заменяться условиями периодичности.
7. Энергия упругого деформирования оболочки
U |
2 V |
|
|
e |
e |
|
|
|
|
|
|
|
|
|
|
|
dv, (4.1.64) |
1 |
|
|
|
z |
z |
z |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
z |
используем далее закон Гука, геометрические соотношения и интегрирование по координате z.
4.2. Оболочка вращения
Оболочки вращения осесимметричны, их координатная поверхность представляет собой поверхность вращения (рис. 4.7).
Рис. 4.7. Оболочка вращения
1. Геометрические соотношения. Пересечение плоскости,
которой принадлежит ось z, с координатной поверхностью образует меридиональную кривую (меридиан). Пересечение плоскости, ортогональной оси z, с координатной поверхностью образует окружность, называемую параллелью. Набор меридиональных кривых и параллелей определяет два семейства координатных
162
линий осесимметричной оболочки. Геометрия оболочки вращения определяется формой меридиана.
Рассмотрим оболочку в системе цилиндрических координат r, φ, z (рис. 4.8). Точки поверхности вращения определяются значениями г и р, координата z не является независимой z = z(r). Определим параметры Ламе, когда α – меридиональное и β – окружное направления:
2 |
dr |
2 |
dz |
2 |
dr |
2 |
|
2 |
, |
ds |
|
|
|
z dr |
|
||||
ds 1 z 2 |
1/ 2 dr Hr dr, |
(4.2.1) |
ds rd H d .
Рис. 4.8. Оболочка вращения в цилиндрической системе координат
Определим главные кривизны поверхности вращения: R1 – радиус кривизны меридиана,
k |
|
1 |
|
|
z |
|
|
. |
(4.2.2) |
|
|
|
|
|
|||||
1 |
|
R1 |
|
(z )2 |
3/2 |
|
|
||
|
|
1 |
|
|
|
Для определения радиуса кривизны R2 необходимо рассмотреть треугольник O1O M и O2O M , эту геометрическую задачу предлагается решить самостоятельно (рис. 4.9).
163
Рис. 4.9. Центры кривизн оболочки вращения
2. Кинематические соотношения. Будем рассматривать де-
формирование оболочки при условии, что внешняя нагрузка также обладает осевой симметрией. Тогда деформирование оболочки будет иметь вид (модель С.П. Тимошенко):
uz u z 1 |
, |
(4.2.3) |
wz w, |
|
|
|
|
всилу осевой симметрии отсутствуют окружные перемещения
vz 0 , а это возможно лишь при v 0 и 2 0 . Кроме того, не-
известные функции u, w, 1 функции аргумента r. Рассмотрим геометрические соотношения:
|
|
|
|
er |
r zkr , |
|
|
|
|
|
||
|
|
|
|
e zk , |
|
|
|
|
(4.2.4) |
|||
|
|
|
|
|
er 0 |
|
|
|
|
|
|
|
в силу осевой симметрии, следовательно, |
|
|
|
|
|
|||||||
|
|
|
|
r kr 0, |
|
|
|
|
(4.2.5) |
|||
r |
1 u |
|
1 |
Hr v |
w |
|
1 |
du |
w |
, |
(4.2.6) |
|
|
|
|
|
|
|
|||||||
Hr r |
Hr H |
|
|
|
||||||||
|
|
|
R1 |
Hr dr R1 |
|
|||||||
164 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
v |
|
1 |
|
H |
u |
w |
|
u |
|
w |
, |
|
|
|
|
|
|
|
|
|||||||||
|
|
H |
|
Hr H |
|
r |
|
R2 Hr r R2 |
||||||||
|
|
|
|
|
|
k |
r |
|
|
1 |
1 |
|
1 |
Hr |
2 |
|
1 |
d 1 |
, |
||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
H |
r |
r |
H |
H |
|
|
|
|
H |
r |
dr |
|
||||||||||
|
|
|
|
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
k |
|
|
1 |
2 |
|
|
1 |
|
|
H |
|
|
1 |
|
|
|
1 |
. |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
H |
|
Hr H r |
|
|
|
|
|
Hr r |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.2.7)
(4.2.8)
(4.2.9)
Кроме того, в соответствии с моделью С.П. Тимошенко нормаль поворачивается в плоскости rz:
|
e |
|
|
1 |
|
dw |
u |
, |
(4.2.10) |
||||
|
1 |
|
|
|
|
||||||||
|
rz |
|
|
Hr |
|
dr |
|
R1 |
|
||||
|
|
|
|
|
|
|
|
||||||
оставшиеся |
компоненты |
тензора |
деформаций равны |
нулю: |
|||||||||
e z ezz 0 . |
Если ввести гипотезу о постоянстве нормалей, т.е. |
||||||||||||
перейти к модели Кирхгофа Лява, то следует принять |
|
||||||||||||
|
|
|
|
er 0 |
|
|
|
(4.2.11) |
|||||
и тогда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
u |
|
1 |
|
dw |
(4.2.12) |
|||||
|
|
Hr |
dr |
||||||||||
|
|
|
|
R1 |
|
и, следовательно, тангенциальные деформации в геометрических соотношениях (4.2.6)–(4.2.7) останутся без изменений, а изгибные деформации примут вид
kr |
1 d |
|
u |
|
1 dw |
|
, |
||||
|
|
||||||||||
|
|
|
|
|
|
||||||
|
|
|
|
Hr dr |
|||||||
|
Hr dr |
R1 |
|
|
k |
1 |
|
|
u |
|
1 |
dw |
||
|
|
|
|
|
|
dr |
. |
||
H |
r |
R |
H |
r |
|||||
|
r |
|
|
1 |
|
|
|
|
(4.2.13)
(4.2.14)
Тогда тензор обобщенных деформаций будет содержать четыре нулевые компоненты: r , , r , , определяемые через две
функции перемещений u и v.
165
3. Физические соотношения. Обобщенный закон Гука для оболочки вращения запишется, как и в общем случае, с использованием матрицы обобщенных жесткостей, но количество ненулевых компонент тензора обобщенных нагрузок будет меньше; например, в случае технической теории оболочек вращения соотношения примут вид:
Nr |
|
A |
||
|
|
11 |
||
|
|
|||
N |
|
|
||
|
|
|
||
Mr |
|
|
||
|
|
|
|
|
sim. |
||||
M |
A |
B |
B |
|
r |
|
|
12 |
11 |
12 |
|
|
|
|
A22 |
B12 |
B22 |
|
|
||
|
(4.2.15) |
|||||
|
11 |
12 |
|
|
. |
|
|
|
kr |
|
|||
|
|
|
|
|
|
|
|
|
22 |
|
|||
|
|
k |
|
4.Энергия упругого деформирования оболочки вращения.
Врамках модели Кирхгофа – Лява
U |
1 |
A 2 2A A 2 |
|||
|
2 |
S |
r |
12 r |
22 |
|
|
11 |
11kr2 2 12 kr k 22 k2 Hr H drd
или после интегрирования по φ
2 |
2 |
U A11 r |
2A12 r A22 |
l
11kr2 2 12 kr k 22 k 2 Hr rdr,
где Aij , ij – обобщенные жесткости оболочки.
(4.2.16)
(4.2.17)
5. Уравнения равновесия оболочки вращения. Рассмотрим оболочку, напряженное состояние которой осесимметрично. Составим уравнения равновесия для элемента оболочки (рис. 4.10).
Спроецируем усилия на нормаль к площадке:
N |
|
N d |
r |
r |
d |
d d N |
|
R d d |
||
|
|
|
||||||||
|
|
|
|
|
|
1 |
||||
|
|
|
|
|
|
(4.2.18) |
||||
|
|
|
|
|
Q r |
|
|
|||
|
|
qR1rd d |
d d 0, |
|||||||
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
166
где d sin d ; таким образом, после деления на |
R1rd d по- |
||||||||
лучим |
|
N |
|
|
|
|
|
|
|
|
N |
|
|
1 d |
Q r q 0; |
(4.2.19) |
|||
|
|
|
|
|
|
||||
|
R1 |
R2 |
R1r d |
||||||
|
|
|
|
|
в преобразовании выражения было использовано вспомогательное равенство:
sin |
|
1 |
; |
|
r |
R |
|||
|
|
|||
|
|
2 |
|
второе уравнение получим, проецируя усилия на ось z:
|
|
|
|
|
|
N |
|
|
|
|
r |
|
|||||
N rd sin N |
|
|
|
|
d r |
|
|
|
d d sin ( d ) |
||||||||
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
qR1rd d cos Q rd cos |
||||||||||||||
|
Q |
Q d |
|
r |
r |
d |
d cos d 0 |
||||||||||
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 d |
|
|
Q cos N |
|
sin r q cos . |
|||||||||
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
rR d |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.2.20)
(4.2.21)
(4.2.22)
Рис. 4.10. Усилия и моменты, действующие на элемент срединной поверхности оболочки вращения
167
Уравнение для моментов составляем относительно оси, касательной к окружности,
1 |
|
d |
M |
|
r cos M |
|
Q 0. |
(4.2.23) |
|
|
|||||||
rR1 |
|
d |
r |
|
|
|||
|
|
|
|
|
Полученных уравнений достаточно для описания деформирования оболочки вращения. Например, в модели С.П. Тимошенко величины N , N , M , M ,Q выражаются через перемещения u, w
иугол поворота 1 , и решается система относительно этих трех функций. В модели Кирхгофа – Лява величины N , N , M , M вы-
ражаются через u, w, и система разрешается относительно u, w,Q .
Для определения искомого решения система полученных обыкновенных дифференциальных уравнений дополняется соответствующими граничными условиями
4.3. Пологие оболочки
Пологой оболочкой называется оболочка, геометрия которой может быть отождествлена с геометрией ее проекции на плоскость (рис. 4.11), а, b – характерные размеры проекции, f (x, y) –
функция, описывающая координатную поверхность, т.е. стрела подъема оболочки.
Рис. 4.11. Пологая оболочка
168
Пологими считаются оболочки, для которых выполняется неравенство
max f |
|
1 . |
(4.3.1) |
|
min(a,b) |
||||
|
5 |
|
В дальнейшем будем рассматривать прямоугольные в плане оболочки.
1. Геометрия оболочки. Так как оболочка пологая, то система координат α, β совпадает с системой координат х, у (см. рис. 4.11) и выражение для квадрата длины дуги
dS 2 |
H 2 d 2 |
|
H 2 d 2 |
(4.3.2) |
||||||||
может быть записано в виде |
|
|
|
|
|
|
|
|
|
|
||
dS 2 |
dx2 |
dy2 , |
(4.3.3) |
|||||||||
т.е. можно считать, что H |
1, H 1 |
и главные кривизны опре- |
||||||||||
деляются соотношениями: |
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
1 |
|
|
2 f |
, |
(4.3.4) |
||||
|
R |
|
|
|
x2 |
|||||||
1 |
|
|
|
|
|
|
|
|
||||
|
|
|
|
1 |
|
|
|
|
|
|
||
k |
|
|
1 |
|
|
2 |
f |
(4.3.5) |
||||
2 |
|
R |
y |
2 |
||||||||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
и являются малыми величинами. Следовательно, гауссова кривизна
Г k1k2 0, |
(4.3.6) |
поэтому пологие оболочки называют оболочками с нулевой гауссовой кривизной.
2. Геометрические соотношения. Учитывая, что выполняют-
ся равенства H |
H 1 |
и то, что x, y, |
получим выраже- |
|||||||||
ния для деформаций |
|
|
|
|
|
|
|
|
|
|||
x |
1 u |
|
1 |
H v |
w |
|
u |
|
w |
, |
||
|
|
|
|
|
|
|
||||||
H |
|
|
|
|
||||||||
|
|
H H |
R1 |
x R1 |
169
|
|
|
|
y |
|
v |
|
|
w |
, |
|
|
|
|
|
|
||||||||||||
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R2 |
|
|
|
|
|
|
(4.3.7) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
u |
|
|
|
v |
|
|
|
|
|
|
||||||
|
|
|
|
xy |
|
|
|
, |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
x |
|
|
|
|
|
|
|
||||
e |
|
|
|
1 |
|
|
|
w |
|
|
u |
|
|
|
w , |
|
||||||||||||
1 |
|
|
|
|
|
|
|
|
1 |
|
||||||||||||||||||
13 |
|
|
|
H |
|
|
|
R1 |
|
|
|
|
|
|
|
x |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.3.8) |
||||||||||||
|
|
|
|
|
|
1 |
|
w |
|
|
v |
|
|
|
|
|
|
w ; |
||||||||||
e |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
2 |
|
|
|
|
2 |
|
||||||||||||||||||||||
23 |
|
|
|
H |
|
|
|
R2 |
|
|
|
|
|
|
|
y |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
kx |
1 |
1 |
|
|
|
1 |
|
|
|
|
H |
2 |
1 |
, |
||||||||||||||
|
|
H H |
|
|||||||||||||||||||||||||
|
H |
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
||||||||||||
|
|
|
|
|
|
|
ky |
2 , |
|
|
|
|
|
|
|
|
|
(4.3.9) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
kxy |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
2 |
|
|
1 |
|
|
|
|
|
|
2 |
. |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
x |
|
|
|
|
В последних двух уравнениях пренебрегаем слагаемыми u
R2
и |
v |
, считая, что тангенциальное смещение не оказывает на |
|
R |
|||
|
|
||
|
2 |
|
сдвиг существенного влияния.
Компоненты тензора деформации оболочки в произвольной точке оболочки на расстоянии z от координатной поверхности
e11 |
x zkx , |
|
e22 |
y zky , |
(4.3.10) |
e12 |
xy zkxy . |
|
Таким образом, компоненты тензора деформаций eij выра-
жаются через функцию прогиба w, тангенциальные смещения координатной поверхности и, v и углы поворота нормали 1 и 2 .
170