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

книги / Равновесие анизотропных полых сферических и цилиндрических тел под действием массовых сил

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

щены статьи и монографии [2, 40, 46, 89], а работы [2, 18, 19, 42, 54, 69, 88, 89] посвящены решению задач определения напряженного состояния в окрестности полостей, моделирующих цилиндрические горные выработки в породных массивах.

Решения задач о равновесии полых толстостенных круглых в поперечном сечении труб, изготовленных из материала с цилиндрической анизотропией частного вида, было впервые записано в 1865 г. B. Saint-Venant [113] и W. Vоigt [120]. Кроме того, решению задач о равновесии однородных анизотропных цилиндров при различных видах осесимметричных нагрузок посвящены работы [10, 16, 17, 41, 54, 55, 70, 71] и монография С.Г. Лехницкого [56]. Исследованию равновесия бесконечно протяженных анизотропных полых центрально симметричных тел, материал которых имеет непрерывную неоднородность (например, зависимость деформационных свойств от пространственных координат), посвящены статьи [57, 82, 84, 102, 114]. Различные аспекты решения задач, учитывающие изменение упругих модулей вдоль радиальной или осевой координаты, а также (в простейшем частном случае) условие пропорциональности деформационных постоянных расстоянию, рассматриваются в работах [26, 53, 80, 81]. Решение ряда частных задач для трансверсально-изотропных цилиндров конечной длины пред-

ставлено в статьях [6, 7, 104, 103, 105, 117].

Проблеме исследования напряженно-деформированного состояния находящихся в равновесии под действием внутреннего давления толстостенных тяжелых полых изотропных цилиндрических тел, жестко закрепленных во всей внешней поверхности или опирающихся на ложементы, посвящены работы Л.Л. Кожевниковой, Г.Б. Кузнецова, В.П. Матвеенко, А.А. Рогового и Н.А. Шевелева [44, 45, 51, 67]. Этими авторами с использованием аналитических и полуаналитических методов описано деформирование и взаимодействие внешних поверхностей тяжелых цилиндров с «жесткими» и «податливыми оболочками», а также конечными по длине опорами.

81

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

Несмотря на то, что управление напряженно-деформиро- ванного состоянием в растущих цилиндрических телах, в которых материал вместе с конструкцией и сооружением создается одновременно, позволяет максимально реализовать резервы несущей способности, в настоящее время существует ограниченное число работ, посвященных развитию теории наращиваемых тел. К числу этих работ можно отнести монографии и статьи Н.Х. Арутюняна, А.В. Манжирова и их учеников [3, 4, 5, 65], а также публикации, к которых рассматриваются модели технологических процессов изготовления цилиндрических тел методами аддитивного наращивания. Так, например, в монографии [98] представлено решение линейной осесимметричной задачи о непрерывном наращивании термоупругого ортотропного полого цилиндра в произвольном нестационарном силовом и температурном поле при любом начальном напряженном состоянии материала, присоединяемого к его внутренней или внешней поверхности (при условии упругой реакции на ненаращиваемой границе). При построении этого решения принимались во внимание только внешние воздействия, вызывающие только радиальные смещения частиц, а в качестве возможного приложения

82

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

вались в работах [8, 13, 43, 74, 75, 95, 101].

Рис. 3.3. Тяжелый ортотропный цилиндр, находящийся под действием произвольной нагрузки

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

Рассмотрим линейно-упругий бесконечно протяженный горизонтальный полый толстостенный цилиндр, ограниченный внутренней и внешней боковыми поверхностями, заданными отно-

83

сительно оси симметрии радиусами a и b ( a < b ). Будем считать, что цилиндрическое тело находится в состоянии равновесия под действием массовых сил и произвольно распределенной симметричной относительно вертикальной плоскости, проходящей через ось цилиндра, нагрузки R(θ). Эта нагрузка может быть задана на

еговнешней иливнутреннейбоковой поверхности (рис. 3.3). Будем считать, что материал цилиндрического тела одно-

родный с постоянной плотностью по всему объему, цилиндрически ортотропный (главные оси ортотропии совпадают с цилиндрической ортогональной системой координат r, θ и z, причем z совпадает с осью симметрии цилиндра). В силу бесконечной протяженности рассматриваемого тела и характера действующих сил будем пренебрегать торцевыми эффектами. Поэтому функции радиальных u и окружных v перемещений, радиальных ( σrr и εrr ), окружных ( σθθ и εθθ ), осевых ( σzz ) нормальных напряжения и деформаций, касательных напряжений и сдвиговых деформаций ( σrz и εrz ) не зависят от осевой координаты z.

В силу малости деформаций будем предполагать неизменность типа упругой симметрии материала цилиндрического тела в процессе приложения нагрузки.

Сформулированные предположения позволяют существенно упростить систему уравнений (3.1)–(3.3), оставив из шести геометрических соотношений Коши (3.1) лишь три:

εrr

=

u ,

εθθ =

u

+ 1 v

,

εrθ =

1

1 u

+

v

v

 

, (3.17)

 

2

 

 

 

r

 

r r ∂θ

 

 

r ∂θ

 

r r

 

из трех уравнений равновесия (3.2) – два:

∂σrr

 

+ 1 ∂σrθ + σrr

− σθθ + F

= 0

,

(3.18)

 

 

 

r

 

r

∂θ

r

 

r

 

 

 

 

 

 

 

 

 

1

∂σθθ

+

∂σrθ

+ 2

σrθ

+ F

= 0 ,

 

 

 

 

 

 

 

r

∂θ

 

r

r

θ

 

 

 

 

 

 

 

 

 

 

84

а из шести определяющих соотношений (3.3) – только четыре

(3.7). Здесь Fr = −γ cosθ и Fθ = γ sin θ – компоненты

вектора

массовых сил, зависящие от удельного веса материала

γ , а ко-

эффициенты Kij , входящие в (3.7), определяются равенствами (3.4). Обратим внимание на то, что для бесконечно протяжен-

ного цилиндрического тела осевые напряжения σzz

не являют-

ся независимыми, определяются при помощи равенства:

μ

rz σrr +

μ

 

 

 

σzz = Ez

 

 

zθ σθθ .

(3.19)

 

Er

Eθ

 

 

Как и ранее, дадим следующее определение. Будем называть горизонтальное толстостенное цилиндрическое тело тяжелым в случае, если величины Fr и Fθ имеют один порядок

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

Запишем неоднородную систему дифференциальных уравнений Ламе в частных производных второго порядка для горизонтального тяжелого ортотропногоцилиндра [27, 28, 32, 33, 99]:

 

 

 

 

 

 

2u

+

 

m2

 

2u

+

1 u

n2

u

 

 

 

 

 

 

 

 

 

 

r2

 

r2

 

∂θ2

r r

r2

 

 

 

(3.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n2

+ m2 v

 

q2 + m2 2v

 

 

 

 

ρ

 

 

 

 

 

 

 

+

 

=

 

cosθ,

 

 

 

 

 

 

 

r2

 

∂θ

 

 

r

 

 

 

 

r∂θ

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

2v

+

p2 2v

+

1 v

v

+

 

k2 + 1 2u

+

 

p2

+ 1 u

= −

ρ

sin θ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r2

r2

 

∂θ2

r r

r2

 

r r∂θ

 

 

 

r

2

 

∂θ

Grθ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эта система была получена последовательной подстановкой геометрических соотношений Коши (3.17) в определяющие (3.7), а затем полученного результата – в уравнения равновесия

(3.18). Здесь:

85

m =

Grθ

, q =

K12

, p =

K22

, k =

K21

, (3.21)

 

 

 

 

 

K11

K11

 

Grθ

Grθ

а показатель анизотропии n вычисляется по формуле (3.9). Симметрия нагрузки относительно вертикальной диамет-

ральной плоскости и геометрическая симметрия тяжелого цилиндрического тела предопределяют поиск решения системы дифференциальных уравнений (3.20) в виде тригонометрических рядов по окружной координате θ :

u = ui (r ) cosi θ ,

v = vi (r ) sin i θ ,

i=0

i=0

подстановка которых в исходную систему (3.20) позволяет записать:

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{ui′′ +

 

 

 

 

 

2

+ m

2

)vi

 

(3.22)

 

 

 

 

 

 

 

r

ui+ i (q

 

 

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

 

 

cosi θ}=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

2

 

2

 

 

 

 

 

 

2

 

 

2

 

 

 

 

γ

 

(m

i

 

+ n

 

)ui

+ i(n

 

+ m

 

 

)vi

 

 

cosθ,

r2

 

 

 

 

 

 

K11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{vi′′

+

1

 

 

i(k

2

+

1)ui

 

 

 

 

 

 

 

 

 

 

 

r

vi

 

 

 

 

 

 

 

 

 

 

 

 

i=0

 

 

 

 

 

 

 

 

 

 

 

}sin i θ = −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

( p2i2

+ 1)vi

 

+ i

( p2

+ 1)ui

γ

sin θ.

 

 

 

 

 

 

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Grθ

Обратим внимание на то, что система обыкновенных дифференциальных уравнений (3.22) относительно радиальной координаты r неоднородна при i = 1, а для i = 0 и i > 1 однородна. Поэтому для различных i решение системы (3.22) можно рассматривать независимо и представить перемещения и напряжения в тяжелом цилиндре в виде суперпозиции:

86

 

 

 

 

 

 

 

 

 

 

u = u0 + u1cosθ + uicosiθ ,

v = v1cosθ + visin iθ ,

(3.23)

 

 

 

i=2

 

 

 

i=2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σrr = σ(rr0)

+ σ(rr1)cosθ + σ(rri)cosiθ ,

σrθ = σ(r1θ)cos θ + σ(riθ)sin iθ ,

 

 

 

i=2

 

 

 

 

 

i=2

 

 

 

0

1

i

 

 

μ

 

μ

zθ

 

σθθ = σ(θθ) + σ(θθ)cosθ + σ(θθ)cosiθ ,

σzz = Ez

rz

σrr

+

 

σθθ .

 

Eθ

 

 

 

i=2

 

 

Er

 

 

Представление (3.23) позволяет выделить механический смысл каждого из слагаемых. Так, например, слагаемые с индексом i = 0 отражают вклад равномерно распределенной нагрузки на внутренней и/или внешней боковой поверхности цилиндрического тела и соответствуют решению классической задачи Ламе. Слагаемые при i = 1 позволяют учесть гравитационные силы, а при i > 1 – неравномерно распределенную нагрузку, которую перед подстановкой в граничные условия необходимо разложить в тригонометрические ряды по окружной

координате.

 

 

При i = 0 отсутствуют окружные перемещения

v0 = 0

(и, как следствие, – v

v′′ ≡ 0 ) и касательные напряжения τ(0) ,

0

0

rθ

а единственное уравнение системы (3.22) принимает вид, аналогичный второму уравнению (3.8), если в последнем произвести замену ur на u0 . Решением этого уравнения является вторая

функция (3.10), в которой также нужно заменить ur на u0 . Это

решение было получено ранее в работе [94]. Поэтому распределения окружных, радиальных и осевых напряжений описываются соотношениями (3.11) с коэффициентами (3.12) и равенством

(3.19), в которых нужно заменить σθθ на σ(θθ0) , σrr на σ(rr0) и σzz на σ(zz0) . Второе выражение (3.10) и формула (3.11) являются решением задачи Ламе для бесконечно длинного ортотропного

87

цилиндрического тела, находящегося под действием равномерного внутреннего и/или внешнего давления.

Решение системы дифференциальных уравнений (3.22) при i = 1 представим в виде суперпозиции:

u1 = u1 + u1* , v1 = v1 + v1*

общего решения соответствующих однородной u1 , v1 и одного из частных решений неоднородной системы:

u* = Br2

,

v* = Ar2 .

(3.24)

1

 

1

 

Подстановка (3.24) в (3.22) позволяет записать систему алгебраических уравнений:

(2q2 + m2 n2 )A + (4 m2 n2 )B =

γ

,

(3.25)

 

 

K11

 

(3 p2 )A (2k 2 + p2 + 3)B = − γ ,

Grθ

решением которой являются неизвестные коэффициенты A и B. Определив эти коэффициенты, запишем искомое частное решение (3.24) в виде:

*

 

Θ

 

γ r

2

 

*

ω4Θ − 1

γr

2

, Θ =

ω3Grθ + ω2K11

 

u1

=

 

 

 

,

v1 =

 

 

 

,

 

 

 

ω3Grθ

 

K11 (ω1ω3 + ω2ω4 )

 

 

Grθ

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω = 4 m2

n2 , ω = 2q2

+ m2 n2 ,

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

 

 

ω = 3 p2 , ω = 2k 2 + p2 + 3 .

 

 

 

 

 

 

 

3

 

 

4

 

 

 

 

Решение однородной системы будем искать в виде обобщенных степенных рядов:

ui = a(ji)r j, j=0

 

vi = b(ji)r j.

(3.26)

j=0

88

Подставив (3.26) в (3.22) и приравняв нулю коэффициен-

ты при r j+ν−2 , получим систему линейных алгебраических уравнений:

 

( j + ν )2 (m2 + n2 ) a(ji) +

+ (q2

 

(3.27)

+ m2 )( j + ν ) (n2

+ m2 ) b(ji) = 0,

(k 2 + 1)( j + ν ) + ( p2 + 1) a(ji) + ( j + ν )2 ( p2 + 1) b(ji) = 0 ,

из решения которой определим характеристические числа ν и за-

пишем выражения для радиальных u1 и окружных v1

перемеще-

нийпри i = 1 [27, 28, 32, 33, 99]:

 

 

 

 

 

 

 

 

u

= a(1)r

β + a(1)

1

 

 

+ a(1)ln r + a(1) + Br2

cosθ ,

(3.28)

 

1

 

1

2

 

rβ

 

3

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v =

a(1)α r

β + a(1) α2

+ a(1)α

ln r

+ a(1)α

4

+ Ar2 sin θ ,

1

1

1

2

 

rβ

 

3

3

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

Θ

 

γ,

B =

ω4Θ − 1

γ,

 

 

 

(3.29)

 

 

 

G

ω

G

 

 

 

 

 

 

 

 

rθ

 

 

 

3

rθ

 

 

 

 

 

а также для радиальных σ(rr1) ,

окружных σ(θθ1)

и касательных σ(r1θ)

напряжений, действующих в плоскости, перпендикулярной оси горизонтального цилиндра [27, 28, 32, 33, 99]:

σ(1) = a(1)Q rβ−1

+ a(1)

rr

 

1

1

2

 

 

 

 

 

σ(1) = a(1)O rβ−1

+ a(1)

θθ

 

1

1

2

 

 

 

 

 

σ(1) = G a(1) J rβ−1 a(1) rθ rθ 1 1 2

Q2

rβ+1

O2

rβ+1

J2

rβ+1

+ Q r + a(1)

K11

cosθ , (3.30)

 

 

3

3

 

r

 

 

 

 

 

+ O r + a(1)

K21

 

cos θ ,

 

3

3

 

r

 

 

 

 

 

(B A)r + a3(1) α3 sin θ.

r

89

Осевые напряжения σ(zz1) определяются формулой (3.19),

вкоторойнеобходимопроизвестизамену σθθ

на σ(θθ1) и σrr

на σ(rr1) .

Входящий в (3.30) коэффициент

 

 

 

 

 

 

 

 

 

 

1

 

1

 

K

 

 

 

 

K

 

 

 

β =

K22

 

+

 

 

12

 

2

+

12

 

+ 1

(3.31)

K11

 

 

 

 

 

 

Grθ

 

K11

 

 

 

Grθ

 

 

является показателем анизотропии для горизонтального осесимметричного тела, находящегося под действием равномерно распределенной нагрузки, имеющей вертикальную плоскость симметрии [27, 28, 32, 33, 99]. Обратим внимание на то, что в отличие от показателя n, определяемого формулой (3.9), β зависит от мо-

дуля сдвига в диаметральной плоскости. Эта зависимость является следствием появления сдвиговых деформаций, которые отсутствуют в случае, когда и нагрузка, и само тело обладают одной и той же осью симметрии. Для изотропного материала показатели анизотропии n и β принимают значения: n = 1 и β = 2 .

Коэффициенты, входящие в формулы (3.28) и (3.30), определяются следующим образом:

 

Q1 = K11β + K12 (1+ α1 ) , Q2 = K12 (1+ α2 ) K11β ,

 

(3.32)

 

 

 

Q =

 

ρ

(

2K

 

+ K

) ω4Θ − 1

+ K

 

Θ

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

G

 

11

 

12

ω

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rθ

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

O1 = K21β + K22 (1+ α1 ) , O2 = K22 (1+ α2 ) K21β ,

 

 

 

 

O3 =

 

ρ

(2K21

+ K22 )

ω4Θ − 1 + K22Θ

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Grθ

 

 

 

 

 

ω3

 

 

 

 

 

 

 

 

 

 

J1 = α1 (β − 1) 1 , J2 = α2 (β + 1) + 1,

 

 

 

 

 

α =

 

β (k2

+ 1)+ p2 + 1

, α

 

=

p2

+ 1− β (k2

+ 1)

, α

 

 

= α

 

= −1 ,

 

 

 

 

 

 

 

2

 

 

 

 

3

4

 

 

 

 

 

 

 

 

 

 

 

1

 

β2 p2

 

1

 

 

 

 

β2 p2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

90