книги / Равновесие анизотропных полых сферических и цилиндрических тел под действием массовых сил
..pdfГЛАВА 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 2−k(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 2−t(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 2−t(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) ) |
ρ1−3 2−k(1) + |
||||
0 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
(2.11) |
|
+C2(10) A11(1) (k(1) − |
1 )+ |
|
|
|
|
|||
2A12(1) |
ρ1−3 2+k(1) = − p, |
|||||||
|
|
|
2 |
|
|
|
|
|
ρ−1 2−k(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) ) ρc−32−k(1) + +C2(10) A11(1) (k(1) − 12)+ 2A12(1) ρc−32+k(1) = = C1(02) 2A12(2) − A11(2) (12 + k(2) ) ρc−32−k(2) + +C2(20) A11(2) (k(2) − 12)+ 2A12(2) ρc−32+k(2) ,
49
ρc−12−k(1) C1(01) + ρc−12+k(1) C2(10) = ρc−12−k(2) C1(02) + ρc−12+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(α ) (1− 2k(α ) ). |
|||||||||||||||||||||
+ |
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