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

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

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

ГЛАВА 2 РАВНОВЕСИЕ УПРУГИХ ТЯЖЕЛЫХ

ТРАНСВЕРСАЛЬНО-ИЗОТРОПНЫХ СОСТАВНЫХ ЦЕНТРАЛЬНО-СИММЕТРИЧНЫХ ТЕЛ

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

2.1. АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ЗАДАЧИ О РАВНОВЕСИИ УПРУГОЙ СОСТАВНОЙ ТРАНСВЕРСАЛЬНО-ИЗОТРОПНОЙ СФЕРЫ ПРИ ДЕЙСТВИИ МАССОВЫХ СИЛ

Рассмотрим равновесие составного линейно-упругого анизотропного центрально симметричного тела, представляющего собой конструкцию, состоящую из двух полых толстостенных сфер с общим центром, ограниченную внутренней и внешней поверхностями радиусов ρ1 и ρ2 ( ρ1 < ρ2 ). В центр

41

поместим начало сферической ортогональной системы координат ρ, θ и ϕ . Будем предполагать, что сферическая поверхность

контакта отдельных частей тела находится на расстоянии ρc от

центра и не изменяет в силу малости деформаций свое положение в процессе нагружения. Считаем, что материал обеих сфер однородный, линейно-упругий, имеющий постоянную плотностью по всему объему, сферически трансверсально-изотропный (относительно любого радиус-вектора, проведенного из центра вданную точку). Все константы и функции, относящиеся к внутренней части составной сферы, будем обозначать верхним индексом1, стоящимвскобках, аквнешней– индексом2 соответственно.

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

( ur(i) и uθ(i) ), радиальные ( σ(rri) и ε(rri) ), меридиональные ( σ(θθi) и ε(θθi) ), окружные ( σ(ϕϕi) и ε(ϕϕi) ) напряжения и деформации, касательные напряжения τ(riθ) и сдвиговые деформации ε(riθ) не

зависят (в силу симметрии тела и внешней нагрузки) от окружной координаты ϕ .

Поэтому уравнения равновесия (1.2) записываем для каждой из частей составного центрально симметричного тела:

 

∂σ(i)

∂τ(i)

 

 

(2σρρ(i)

− σ(ϕϕi) − σ(θθi) + τρθ(i) ctg θ)+ Fρ(i) = 0 , (2.1)

 

ρρ

+ 1

 

ρθ

+ 1

 

 

 

∂θ

 

 

∂ρ

 

ρ

ρ

 

 

 

 

 

 

 

 

∂τρθ(i)

 

+ 1

∂σ(θθi)

+ 1 (σ(i) − σ(i)

)ctg θ + 3τ(i)

+ F(i) = 0 .

 

 

 

∂ρ

 

 

 

 

ρ

∂θ

 

ρ

θθ

ϕϕ

ρθ

θ

Здесь

i = 1 или 2 ,

F(i) = −γ(i) cosθ

и F (i) = γ(i) sin θ – компо-

 

 

 

 

 

 

 

 

 

ρ

 

 

θ

 

ненты вектора массовых сил, γ(i)

– удельный вес материалов.

42

Определяющие соотношения (1.3), представленные в виде:

σ(ρρi) = A11(i)ερρ(i) + A12(i) (ε(ϕϕi) + ε(θθi) ),

σ(i)

= A(i)ε(i)

+ A(i)ε(i)

 

+ A(i)ε(i) ,

(2.2)

ϕϕ

12 ρρ

22 ϕϕ

 

23

θθ

 

σ(i) = A(i)ε(i)

+ A(i)ε(i)

+ A(i)ε(i)

,

τ(i)

= 2A(i)ε(i)

,

θθ 12 ρρ

23 ϕϕ

22 θθ

 

ρθ

44 ρθ

 

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

m(i) A11(i) = E(i) (1− ν(i) ),

m(i) A12(i) = E(i)ν(i) ,

A44(i) = G(i) , (2.3)

A(i)

=

 

 

E(i)

1− ν(i)2

E(i)

 

,

 

 

 

 

 

 

 

 

22

 

 

(1+ ν(i) )m(i)

 

 

 

 

 

 

 

 

 

 

 

E(i)

 

 

A(i)

=

 

 

 

 

E(i)

 

 

ν(i) + ν(i)2

E(i)

 

,

(1

 

 

 

 

 

 

23

 

 

 

 

 

 

E(i)

 

 

 

+ ν(i) )m(i)

 

 

 

 

а m(i) = 1− ν(i) 2

(ν(i)

)2

E(i)

E(i)

; E(i) и

E(i)

– модули Юнга

для растяжения вдоль ρ и в ортогональном к нему направлении;

G(i) – модуль сдвига для диаметральной плоскости; ν(i) и ν(i) – коэффициенты Пуассона, характеризующие поперечную деформацию в направлениях θ и ϕ при растяжении вдоль радиальной ко-

ординаты ρ и в направлении ρ при растяжении в направлениях θ и ϕ соответственно.

Последовательная подстановка геометрических соотношений Коши, такжезаписанныхдлякаждойизчастейсоставноготела:

 

 

(

)

 

 

 

(i)

 

 

 

 

i

 

 

 

 

uθ

 

 

ε(i)

=

uρ

 

,

ε(i) = 1

 

+ u(i) ,

(2.4)

∂ρ

 

∂θ

ρρ

 

 

 

θθ ρ

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(i)

 

 

(i)

ε(i)

= 1

(u(i) ctg θ + u(i) ) ,

ε(i)

= 1

 

uθ

+ 1

 

uρ

∂ρ

∂θ

ϕϕ

ρ

θ

ρ

ρθ

2

 

ρ

 

 

 

 

 

 

 

 

 

 

 

( ) uθi

43

в определяющие (2.2), а затем полученного результата – в уравнения равновесия (2.1) позволяет записать неоднородные системы дифференциальных уравнений Ламе в частных производных:

 

 

 

 

 

 

 

 

 

2

 

(i)

1

 

 

 

 

 

 

 

(i)

 

 

 

 

 

 

2

 

(i)

 

 

 

 

 

 

 

 

 

 

 

A44(i)

 

uθ

 

 

+

2A44(i)

uθ

+ (A12(i) + A44(i) )

 

uρ

 

+

 

 

 

 

 

 

 

 

 

 

∂ρ2

 

 

ρ

 

∂ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ρ∂θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

(i)

 

 

 

(i)

 

 

 

 

 

 

 

 

 

 

 

 

 

(i)

 

 

 

+

 

 

A22(i)

 

 

uθ

+

 

uθ

 

ctg θ

+ (A22(i) + A23(i)

+ 2A44(i) )

uρ

(2.5)

ρ2

∂θ2

∂θ

 

∂θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(A(i) + 2A(i)

+

A(i)

ctg2 θ)u(i) = −γ(i) sin θ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

 

 

 

44

 

22

 

 

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

(i)

 

 

1

 

 

 

 

 

 

 

(i)

 

 

 

 

 

2

 

(i)

 

 

 

 

(i)

 

 

 

A11(i)

 

 

uρ

+

 

2A11(i)

 

uρ

 

+

(A12(i) +

A44(i) )

 

uθ

 

+

uθ

ctg

θ

+

 

 

 

 

 

 

 

 

∂ρ

 

∂ρ∂θ

 

 

 

 

∂ρ2

 

 

ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ρ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2u(i)

 

 

u(i)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

A44(i)

 

 

 

 

 

ρ

+

 

ρ

 

ctg

θ + 2(A12(i) A22(i)

A23(i) )uρ(i) +

 

ρ2

 

 

∂θ2

∂θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(i)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ (A(i)

A(i)

A(i)

A(i) )

uθ

 

+ u(i)

ctg θ

= γ(i) cosθ,

 

 

 

 

 

 

 

 

 

 

 

 

12

 

 

 

 

22

 

 

23

 

44

 

∂θ

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

лено в виде рядов [63, 72, 106, 107, 118, 119]:

uρ(i) = uρ(in) (ρ) Pn (cosθ) , uθ(i) n=0

dPn (cosθ)

 

= uθ(in) (ρ)

, (2.6)

dθ

n=0

 

 

 

по полиномам Лежандра, подстановка которых в выражения (2.5) и последующее приравнивание коэффициентов левых и правых частей уравнений при одинаковых функциях аргумента θ позволяет получить бесконечное число систем обыкновенных дифференциальных уравнений [34, 35, 38, 39]:

 

(i)

 

(i)

 

(i) 1

(i)

 

 

(i)

 

1

 

(i)

 

 

(i)

1

(i)

 

 

(i)

1

 

 

(i)

 

(i)

 

 

 

a

 

u′′

 

+ a

2n ρ

u

 

+ a

 

 

 

u

 

 

+ a

 

ρ

u

 

+ a

 

 

 

u

 

 

=

A

 

, (2.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1n

 

ρn

 

ρn

 

3n ρ2

ρn

 

 

4n

θn

 

5n ρ2

 

 

θn

 

n

 

 

 

 

 

 

 

(i)

 

(i)

 

(i) 1

 

(i)

 

 

 

(i)

1

 

(i)

 

 

(i) 1

 

(i)

 

 

(i)

1

 

(i)

 

 

 

(i)

 

 

b

 

u′′

+ b

 

 

u

+ b

 

 

 

u

 

+ b

 

 

u

+ b

 

 

 

 

u

 

=

B

 

,

 

 

 

 

 

ρ2

 

 

 

 

ρ2

 

 

 

 

1n

 

θn

2n ρ

 

θn

 

 

3n

 

θn

 

4n ρ

 

ρn

 

 

5n

 

ρn

 

 

n

 

44

где

a1(ni) = A11(i) , a2(in) = 2A11(i) , a3(in) = 2(A12(i) A22(i) A23(i) )A44(i)n(n + 1) , a4(in) = (A12(i) + A44(i) )n(n + 1) , a5(in) = (A12(i) A22(i) A23(i) A44(i) )n(n + 1) ,

b(i) = A(i)

, b(i) = 2A(i) , b(i) = − A(i) 2A(i)

A(i) (n2

+ n 1),

1n

44

 

2n

 

 

44

3n

23

 

44

 

 

22

 

 

 

b(i) = − A(i) A(i) ,

b(i) = − A(i)

A(i)

2A(i) ,

 

 

 

 

4n

12

44

5n

 

22

 

23

 

44

 

 

An(i)

 

(i)

 

 

 

Bn(i)

 

 

(i)

 

 

 

 

= γ

 

,

n = 1 ,

=

−γ

 

 

,

n = 1 ,

 

 

 

0, n = 0 , n > 1 ;

 

0, n = 0 , n > 1 .

Определение перемещений в точках тяжелой толстостенной трансверсально-изотропной составной сферы сводится к решению n самостоятельных задач (2.7), каждая из которых должна быть дополнена граничными условиями. Для этого необходимо разложить в ряды по полиномам Лежандра заданные на границах тела перемещения, напряжения или их комбинации, отражающие условия закрепления и приложенные нагрузки.

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

 

 

 

 

 

u

(i)

= C

(i)

 

1 2k(i)

+ C

(i)

1 2+k(i)

+

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

ρ

 

 

 

 

 

 

 

 

 

 

ρ

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

(i)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ x(i)C(i) + x(i)

C21

+ x(i)C(i)ρ1 2+t(i)

+ x(i)C(i)ρ1 2t(i)

+ H

(i)ρ2

cosθ +

(2.8)

 

 

1

1

2

 

ρ

 

3

 

3

 

 

 

4

4

 

 

 

ρ

 

 

1

1

1

 

 

1

 

1

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ(i)

 

 

 

 

λ(i)

 

λ(i)

 

λ(i)

 

 

+ x1(i)C1(i)ρ

1n + x2(i)C2(i)ρ

2n + x3(i)C3(i)ρ

3n + x4(i)C4(i)

ρ

4n Pn (cosθ) ,

 

 

 

n

n

 

 

 

n n

 

 

n

n

 

n

n

 

 

 

 

 

n=2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

45

 

 

 

 

(i)

 

 

 

 

 

 

 

 

 

(i)

 

 

u(i) = C(i) +

C21

 

+ C(i)ρ1 2+t

(i) + C(i)ρ1 2t(i) + H

ρ2 sin θ +

ρ

 

θ

θ

 

11

31

 

 

 

 

41

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ(i)

 

 

λ(i)

 

 

 

λ(i)

λ(i)

dP

(cosθ)

 

+ C1(i)ρ

1n + C2(i)ρ

2n + C3(i)ρ

3n + C4(i)ρ

4n

 

 

n

 

.

 

 

 

 

 

 

n

 

 

n

 

 

 

n

n

 

 

 

 

dθ

 

 

n=2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь, как и прежде,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k(i) =

1 + 2

A22(i) + A23(i) A12(i)

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

A11(i)

 

 

 

 

 

 

 

t(i) =

9 +

 

2(A22(i) + A23(i) A12(i) )

+

A11(i) (A22(i) + A23(i) )2A12(i) (A12(i) + 2A44(i) )

 

 

A11(i)

 

 

 

A11(i) A44(i)

 

 

 

 

4

 

 

 

 

 

 

 

 

 

являются показателями анизотропии для центрально-симмет- ричного тела, находящегося под действием равномерно распределенных центрально симметричной и вертикальной осесимметричной нагрузки [34, 35, 38, 39], а

x3(1i)

x4(1i)

 

x(i) = −1 ,

x(i)

=

A22(i) + A23(i) + 2A44(i)

 

,

 

 

 

 

 

(i)

(i)

 

(i)

(i)

 

 

 

 

 

 

 

11

 

21

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

A

A

A

 

 

 

 

 

 

 

 

 

 

 

 

 

12

22

 

23

44

 

 

 

 

 

 

 

(i)

(i)

(i)

 

(i)

 

(i)

 

(i)

+

(i)

)t

(i)

 

 

=

A44

A12 2

(A22

+ A23

)3A44

+ 2(A12

A44

 

,

 

 

2(A12(i) + A44(i) )2 A11(i) (A22(i)

+ A23(i) + 2A44(i) )

 

 

 

 

 

 

 

 

 

 

(i)

(i)

(i)

+

(i)

 

(i)

 

(i)

+

(i)

)t

(i)

 

 

=

A44

A12 2

(A22

A23

)3A44

2(A12

A44

 

 

,

 

 

2(A12(i) + A44(i) )2 A11(i) (A22(i)

+ A23(i) + 2A44(i) )

 

 

 

 

 

 

 

 

 

 

Hρ(i) = γ((ii)) (A22(i) + A23(i) 2A12(i) 2A44(i) ),

H

46

H (i) =

 

γ(i)

(2A(i) A(i)

A(i) 2A(i) ),

 

 

θ

 

H (i)

11

22

23

44

 

 

 

 

 

 

H(i) = 2 A11(i) (A22(i) + A23(i) 4A44(i) )+ 2A44(i) (A22(i) + A23(i) 3A12(i) )2A12(i)2 .

Входящие в уравнения

(2.8)

характеристические числа

λ(ji) ( j = 1, 2, 3 и 4 )

 

должны быть вычислены по отдельности

n

 

 

 

 

 

 

для внутренней (4

числа при

i =1)

и внешней (4 числа при

i = 2 ) составной части сферы. Эти характеристические числа определяются для каждого n > 1 из решения нелинейного уравнения четвертой степени:

(i)

(

λ

(i)

 

(i)

(i)

(

λ

(i)

+1

+ α

(i)

=

(

α

(i)

λ

(i)

+ α

(i)

(i)

λ

(i)

(i)

,

λ

n

 

+1 + β λ

n

 

3n

4n

n

 

 

β

n

+ β

n

 

 

)

3n

n

 

 

)

 

 

 

 

 

 

 

5n )(

4n

 

5n )

 

 

 

 

 

 

 

(i)

=

 

λ(ji)

(λ(ji) + 1− α(4in) )+ β3(in) − α5(in)

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

x jn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ(i)

(

λ(i) + 1− β(i)

)

+ α(i)

− β(i)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

n

 

 

j

 

4n

 

 

 

 

3n

 

5n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где α(i)

= a(i)

a(i)

и β(i)

= b(i)

b(i) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jn

 

jn

1n

 

 

 

jn

 

 

jn

1n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Постоянные интегрирования C1(0i) и C2(i0) определяются из граничных условий для краевой задачи, соответствующей n = 0 ; константы C1(1i) , C2(1i) , C3(1i) и C4(1i) – из граничных условий крае-

вой задачи при n =1, а постоянные C(i) ,

C(i) ,

C

(i)

и C(i) – из

1

2

n

3

4

n

n

 

 

n

 

граничных условий краевых задач, записанных

для каждого

n > 1 . Обратим внимание на то, что при n = 0

и n >1 системы

дифференциальных уравнений однородны и имеют тривиальные решения: uρ(in) = 0 и uθ(in) = 0 , если граничные условия, разложен-

ные в ряды по полиномам Лежандра при n = 0 и n >1, также однородны.

47

2.2. РАВНОВЕСИЕ УПРУГОЙ ТЯЖЕЛОЙ СОСТАВНОЙ ТРАНСВЕРСАЛЬНО-ИЗОТРОПНОЙ СФЕРЫ С ЖЕСТКО ЗАКРЕПЛЕННОЙ ВНЕШНЕЙ ПОВЕРХНОСТЬЮ

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

пряжения:

u

(1)

 

 

 

= u(2)

 

 

 

 

, u

(1)

 

 

 

 

 

= u(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

ρ=ρc

ρ

 

 

ρ=ρc

 

θ

 

ρ=ρc

θ

 

ρ=ρc

τρθ(1)

 

 

= τρθ(2)

 

 

, σρρ(1)

 

 

 

= σρρ(2)

 

 

 

ρ=ρ

 

ρ=ρ

c

 

 

ρ=ρ

c

 

 

ρ=ρ

c

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, (2.9)

.

Рис. 2.1. Тяжелая трансверсально-изотропная составная сфера

48

Как и ранее, будем предполагать, что материал обеих составных частей линейно-упругий однородный, с постоянной плотностью по всему объему, сферически трансверсальноизотропный, относительно любого радиус-вектора, проведенного из центра в рассматриваемую точку.

Поскольку граничные условия на контактной (2.9), внешней и внутренней поверхности (ограничены радиусами ρ1 и ρ2 , проведенными из общего центра):

uρ(2)

 

ρ=ρ

= 0 , uθ(2)

 

ρ=ρ

= 0 ,

τρθ(1)

 

ρ=ρ

= 0 , σρρ(1)

 

ρ=ρ = − p , (2.10)

 

 

 

 

 

 

 

2

 

 

2

 

1

 

1

не нарушают осевой симметрии задачи, то ее решение может быть получено с помощью метода, описанного в параграфе 2.1.

Дополним каждую из неоднородных систем дифференциальных уравнений (2.7) при n = 0, 1, 2, граничными условиями (2.9) и (2.10), разложенными в ряды по меридиональной координате θ . Константы интегрирования для n = 0 , входящие в слагаемые выражений (2.8), записываются из решения системы четырех линейных алгебраических уравнений [34]:

C1(1)

2A12(1)

A11(1) (1

+ k(1) )

ρ13 2k(1) +

0

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

(2.11)

+C2(10) A11(1) (k(1)

1 )+

 

 

 

 

2A12(1)

ρ13 2+k(1) = − p,

 

 

 

2

 

 

 

 

 

ρ1 2k(2) C(2) + ρ1 2+k(2) C

(2) = 0 ,

2

 

1

 

2

 

 

2

0

 

 

0

 

 

 

 

 

C1(01) 2A12(1) A11(1) (12 + k(1) ) ρc32k(1) + +C2(10) A11(1) (k(1) 12)+ 2A12(1) ρc32+k(1) = = C1(02) 2A12(2) A11(2) (12 + k(2) ) ρc32k(2) + +C2(20) A11(2) (k(2) 12)+ 2A12(2) ρc32+k(2) ,

49

ρc12k(1) C1(01) + ρc12+k(1) C2(10) = ρc12k(2) C1(02) + ρc12+k(2) C2(20)

следующим образом [34]:

Z0(α ) = 4A11(β)B(α)k(β)ρβ2k(β) b+(β)c(β) (b(β)c(α ) + B(α ) ), (2.12)

C1(01)Z0(2) = 2 pρc2k(1) ρ132+ k(1) (b(1)c(2) + B(2) ),

 

 

C(1)Z

(2) = −2 pρ3 2+k(1) (b(1)c(2) + B(2) ),

 

 

 

 

 

 

 

2

0

0

 

 

1

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

) =

 

1

 

 

(1)

 

(2)

ρ13 2+k

(1)

,

β =

1, α = 2,

 

C2(

)Z0(

8A11( ) pk(1)ρсk

 

+k

 

 

 

 

α = 1.

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1(

2

)Z0(

) = −8A11( )

pk(1)ρсk

(1)

+k

(2)

ρ13 2+k

(1)

ρ22k

(2)

,

 

 

 

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с(α ) = ρ

2k(α )

− ρ

2k(α ) , B(α ) = b(

α )ρ

2k(α ) b(α )ρ 2k(α )

,

 

 

 

 

с

 

 

α

 

 

 

 

 

 

 

 

+

 

α

 

 

c

 

b(α ) = 4A(α )

A(α ) (1+ 2k

(α ) )

, b(α )

= 4A(α ) A(α ) (12k(α ) ).

+

12

 

 

 

11

 

 

 

 

 

 

 

 

 

 

12

 

11

 

 

 

При n =1 постоянные интегрирования в формулах (2.8) находятся из решения восьми линейных алгебраических уравнений:

2C

(1) (x(1)

+ 2)ρ1 2+t(1)

+ C(1) (2x(1)

2t(1) + 3)ρ1+2t(1)

+

 

2

2

 

1

 

3

3

 

 

 

1

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

(2.13)

 

 

(1) (2x(1)

+ 2t(1)

+ 3)ρ = 2

(H

(1) H

(1) )ρ7 2+t(1)

 

+C

,

 

 

 

4

4

 

 

1

 

θ

 

ρ

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

2C2(11) 2A12(1) (x2(1) +1)A11(1) x2(1) ρ112+t(1) +

+C3(11) 4A12(1) (x3(1) +1)+ A11(1) x3(1) (2t(1) 1) ρ11+2t(1) +

+C4(11) 4A12(1) (x4(1) +1)A11(1) x4(1) (2t(1) +1) ρ1 =

= −4 Hρ(1) (A11(1) + A12(1) )+ Hθ(1) A12(1) ρ172+t(1) ,

50