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

книги / Строительная механика, динамика и устойчивость композитных конструкций

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
3.67 Mб
Скачать

силы на ось 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