Добавил:
Выполняю работы студентам ОмГТУ, ФТНГ, авиа-ракетных специальностей Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

диОСК

.pdf
Скачиваний:
19
Добавлен:
22.06.2022
Размер:
882.09 Кб
Скачать

Геометрические соотношения типа в цилиндрических координатах для

осесимметричной задачи преобразуются к виду:

εα =

 

 

 

1

 

 

 

u

 

 

 

 

 

 

 

uβ

 

 

 

 

 

 

A

 

 

 

 

 

uγ

 

 

 

1

 

 

u

 

 

 

 

uβ

 

 

 

 

1

 

 

 

 

uγ

 

 

 

 

u

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

α

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

=

 

 

 

 

 

 

 

 

α +

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

=

 

 

α

 

 

 

A

 

 

α

 

 

A B

β

 

 

 

 

R1

 

1

 

 

 

 

 

1 R

 

β

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εβ =

 

 

 

1

 

 

 

 

uβ

+

 

 

 

 

u

 

 

 

 

 

 

B

 

 

+

 

 

uγ

=

 

1

 

 

uβ

+

 

u

 

 

 

 

R

+

 

uγ

=

 

uγ

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

β

 

 

A B

 

α

 

 

 

 

R2

 

R

 

β

1 R

 

α

 

 

R

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εγ =

 

uγ

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

1

 

 

 

 

uγ

 

 

 

 

u

 

 

 

 

 

 

 

u

 

 

 

 

1

 

 

uγ

 

 

 

 

u

 

 

 

 

u

 

 

 

 

 

 

uγ

 

 

,

 

 

 

γαγ =

 

 

 

 

 

 

 

α +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

=

 

 

α

 

+

1

 

 

 

 

 

 

 

 

α

=

 

 

 

α

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

α

 

 

R

γ

 

 

 

α

 

 

 

 

 

 

 

γ

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γβγ =

 

 

uβ

+

 

1

 

 

uγ

 

 

uβ

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Две последних

 

 

γ

 

B

 

 

β

 

R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

u

 

 

 

 

1

 

 

 

 

uβ

 

 

 

 

 

 

 

u

 

 

 

 

A

 

 

uβ

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

деформации равны

γαβ

=

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

нулю по условию

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

β

 

 

A

 

 

 

α

 

A B

β

A B

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

симметричности

 

То есть, подводя итоги:

 

 

 

 

 

 

 

u

 

uγ

 

 

 

uγ

 

 

 

u

 

uγ

 

Совпадает с зависимостями

εα =

α ,

εβ =

 

,

εγ

=

 

,

γαγ

=

α

+

 

,

Сопротивления материалов

R

γ

γ

α

 

α

 

 

 

 

 

 

 

 

 

 

На основании гипотез Кирхгофа получим частную форму выражений для перемещений для рассматриваемого случая

uγ = w(α), uα = u(α)+γ ϑα , ϑα = −w

где w(α) и и(α) — функция прогиба и осевое перемещение точек срединной поверхности, а ( )' обозначает производную по α.

Подставляя данные выражения в геометрические соотношения имеем

 

uα

= u(α)

+γ ϑα

0

+γ κα ,

 

 

 

 

 

 

 

εα =

α

=εα

0

′′

0

=

w

.

εβ =

uγ

=

w

=εβ0

,

 

где

εα

= u (α),

κα = −w , εβ

R

 

 

 

 

 

 

 

 

R

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Закон Гука (относительно напряжений) принимает вид:

σα =

 

 

 

(εα0 + µ εβ0 +γ κα ),

 

 

 

E

E

 

 

 

 

где E =

 

 

(εβ0 + µ εα0 +γ µ κα ).

 

 

σβ =

 

 

1µ2

E

Подставляем в закон Гука выражения для εα0 , εβ0 , κα

и находим усилия и моменты,

выражая их через напряжения

 

 

 

 

 

Nα =

Nβ =

Mα =

M β =

h 2

 

 

 

 

 

w

 

 

 

 

+ µ

σα dγ = E u

,

h 2

 

 

 

 

 

R

h 2

 

 

w

 

 

 

 

 

 

 

 

 

 

σ β dγ = E

R

+ µ u,

h 2

 

 

 

 

 

h 2

σα γ dγ = −D w′′,

h2

h 2

σ β γ dγ = −µ D w′′,

h2

здесь, как и ранее: E = E(1µ2 ), D = E h3 12

Выведем теперь уравнения равновесия. Выделим из оболочки элемент, показанный на рисунке. б и приравняем нулю сумму всех сил, действующих в осевом направлении α, сумму проекций сил на нормаль к поверхности и сумму моментов относительно элемента параллели (остальные уравнения равновесия удовлетворяются тождественно). Получим:

Nα′ +q = 0, Qα′ −

Nβ

+ p = 0,

Mα′ −Qα = 0

 

 

R

 

Полученные соотношения представляют собой систему основных уравнений для осесимметрично нагруженной цилиндрической оболочки. Семь полученных уравнений включают столько же неизвестных — два перемещения и, w и пять

усилий и моментов Nα, Nβ, Qα, Mα, Mβ. Как отмечалось ранее, эту систему можно свести к уравнениям, включающим в качестве неизвестных перемещения.

Из последнего уравнения и выражений для погонных нагрузок получаем:

Qα = Mα′ = −D w′′′

Подставляем в первые два уравнения из последней группы выражения погоных нагрузок

 

′′

 

µ

 

q

,

1

 

 

w

 

h2

 

IV

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

+

R

w = −

E h

R

 

µ u

+

 

+

12

w

 

=

E h

 

 

 

 

 

 

 

 

R

 

 

 

 

Во многих случаях задача является статически определимой относительно усилий Να — они могут быть получены путем интегрирования уравнения,

Nα′ +q = 0

причем произвольная постоянная определяется по известному значению Να = N на каком-либо торце оболочки. В этих случаях вместо двух уравнений можно получить одно уравнение относительно прогиба w.

Пусть на оболочку действует только

внутреннее давление р. При q=0 из

первого уравнения получим Nα=N (если

р=const, то N = 0,5 p R

). Исключая из

формул и’ найдем

Nβ′ = ERh w + µ N

В этом случае второе уравнение после подстановки выражений для Qα и Νβ

примет вид:

D

w

IV

+

E h

 

w = p µ

 

N

 

 

 

 

 

 

R2

 

 

R

 

 

 

или

 

 

IV

+4

 

 

4

 

 

 

1

 

 

 

N

w

 

 

k

 

w =

 

 

p µ

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

где

k

4

=

 

E h

 

=

3

(1µ2 )

 

 

 

 

 

 

 

4 R2 D

 

 

R2 h2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Краевой эффект и безмоментное состояние

Общее решение дифференциального уравнения запишем в виде:

w = ek α (C1 cos kα +C2 sin kα)+ek α (C3 cos kα +C4 sin kα)+ w0

где w0(α) — частное решение неоднородного уравнения, зависящее от вида правой части. Четыре произвольных постоянных C1234 определяются из четырех граничных условий при α=0 и α=l.

Если оболочка рассматривается как полубесконечная (l→∞), то из условия ограниченности решения при α→∞ следует положить С2 = С3 = 0. Тогда

w = ek α (C1 cos kα +C2 sin kα)+ w0

Здесь константы С1 и С2 определяются из граничных условий на краю α=0; они характеризуют быстро затухающее решение типа краевого эффекта, описывающее местный изгиб оболочки вблизи края α=0. Так как eπ≈0,043, то влиянием краевого эффекта практически можно пренебречь при α≥π·k или, если

положить μ≈0,3, при α 2,5 R h . Например при R/h =100 и R/h =400

ширина

зоны краевого эффекта равна 0,25R и 0,125R соответственно.

 

 

Если длинa оболочки превышает ширину зоны краевого эффекта ( l 2,5

R h

), то

такую оболочку можно считать длинной и изгиб вблизи каждого из двух краев и можно рассматривать независимо друг от друга, используя решение для полубесконечной оболочки (при этом в каждом случае координата α отсчитывается от рассматриваемого края в направлении другого края).

,

Если правая часть дифференциального уравнения представлена в виде полинома по степеням α с показателями, не превышающими трех, то частное решение имеет

вид:

 

R2

 

N

w0

 

=

 

p µ

 

 

E h

 

 

 

 

R

В более общем случае, если правая часть ДУ меняется достаточно плавно, данное выражение может быть использовано как приближенное частное решение.

Полученное частное может быть также получено как решение ДУ при D=0.

Из равенств и соотношений, определяющих распределение напряжений по толщине оболочки, следует, что при D=0, Mα=Mβ=Qα=0 и σα=Nα/h, σβ=Nβ/h

т. е. напряжения равномерно распределены по толщине. Такое напряженное состояние оболочки называется безмоментным и соответственно частное решение называется безмоментным решением.

Таким образом, можно считать, что общее решение задачи при осесимметричной деформации тонкой достаточно длинной цилиндрической оболочки складывается

из безмоментного решения w0 и решения краевого эффекта wk.

Краевой эффект возникает в том случае, когда на краю оболочка нагружается моментным образом (поперечными силами, моментами), соединяется с другой оболочкой шпангоутом и т. п. или закрепляется в отношении поперечного перемещения w и угла поворота θα. Угол поворота нормали, изгибающие моменты и перерезывающая сила согласно формулам имеют вид:

ϑα = −w′ = k ek α [C1 (cos kα +sin kα)C2 (cos kα sin kα)],

Mα = −D w′′ = −2 D k 2 ek α (C1 sin kα C2 cos kα),

M β = µ Mα ,

Qα = −D w′′′ = −2 D k3 ek α [C1 (cos kα sin kα)C2 (cos kα +sin kα)]

Определим постоянные С1 и С2 для различных типов граничных условий.

Для защемленного края (см. рис. б) при α=0 следует принять:

w =ϑα = 0

при этом С12=-w0

w = w0 [1ek α (cos kα sin kα)]

Если на торце оболочка закреплена шарнирно (рис. в), то

α = 0 w = Mα = 0

получаем

C1 = −w0 , C2 = 0

решение имеет вид:

w = w0 (1ek α cos kα).

Из полученных равенств следует, что при увеличении α прогиб приближается к безмоментной составляющей, т. е. на достаточном удалении от края w=w0

Если оболочка не закреплена в отношении радиального перемещения w и угла поворота ϑα, то при α=0 Mα=Qα=0 и следовательно С12=0 и w=-w0 т. е. при таком закреплении краевой эффект отсутствует и прогиб полностью определяется безмоментным решением.

Таким образом, достаточно длинная цилиндрическая оболочка почти всюду (за исключением быть может узких участков вблизи краев) находится в безмоментном состоянии, а в окрестности краев при соответствующих граничных условиях реализуется локальное изгибное состояние — краевой эффект. Эта особенность напряженного состояния характерна не только для цилиндрической, но и для более широкого класса оболочек и позволяет в ряде случаев значительно упростить расчет.

Безмоментное напряженное состояние

Основная система уравнений безмоментной теории оболочек для случая осесимметричной деформации принимает следующий вид:

(r N

α

 

)

N

β

R cosα + q R r = 0,

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

N

α

 

+

 

N

β

= p,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

R

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

εα

=

 

 

1

 

(u′+ w), εβ

=

1

(u cosα + w sin α),

 

R1

r

 

 

 

 

 

 

 

 

 

 

 

 

Nα = E h (εα + µ εβ ), Nβ = E h (εβ + µ εα )

Радиальное и осевое перемещения определяются

ur = u cosα + w sin α, ux = u sin α w cosα

угол поворота нормали к срединной поверхности

ϑα = 1 (u w)

R1

В уравнениях принято qα = q, qγ = р, ( )' обозначает производную по α.

Соседние файлы в предмете Основы расчета оболочек