
книги / 393
.pdf
метральной плоскости. Эта зависимость является следствием появления сдвиговых деформаций, которые отсутствуют в случае, когда и нагрузка, и само тело обладают одной и той же осью симметрии.
Обратим внимание на то, что в выражениях (16) для параметров ξi и ϕi под знак квадратного корня входят комплексы, состоящие из упру-
гих модулей материала. Поскольку вышеперечисленные параметры должны принимать только действительные значения, можно записать ограничения в форме неравенств на возможные величины упругих модулей материала иусловия существования решения при i >1:
K22 |
+i |
2 |
K22 |
− |
K12 |
|
K12 |
−2 |
|
+1 ≥ ±2 |
(i |
4 |
−2i |
2 |
+1) |
K22 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
K |
|
G |
K |
G |
|
|
K |
|||||||||||
11 |
|
|
rθ |
|
11 |
|
rθ |
|
|
|
|
|
|
|
|
11 |
|
i |
2 |
K22 |
− |
K12 |
|
K12 |
−2 |
|
+1 > 0. |
|
|
|
|
|
|
|
|
||||
|
G |
K |
G |
|||||||
|
|
rθ |
|
11 |
|
rθ |
|
|
|
В частном случае Er = Eθ = Ez = E , µrθ = µθz = µzr = µ и Grθ = G = E[2(1+ ν)], который соответствует изотропному тяжелому
горизонтальному цилиндру, новое аналитическое решение совпадает с результатом, полученным автором [2] с точностью до коэффициентов любого частного решения. Учет анизотропии материала позволяет установить связь между коэффициентами, входящими в равенства (12):
Θ |
= |
ω4Θ −1 |
= |
ρ |
1−µ − 2 |
(µ)2 |
. |
G |
|
||||||
|
ω G |
|
2E |
|
|
||
rθ |
|
3 rθ |
|
|
|
|
|
Покажем одно из возможных приложений найденного аналитического решения на примере задачи о равновесии тяжелого железобетонного цилиндра, наполовину закопанного в грунт. Предположим, что площадь контактной поверхности и реакция основания известны и неизменны во времени, а сама реакция распределена по боковой поверхности цилиндра по параболическому закону
R (θ)= Aθ2 + B θ+C . |
(17) |
41

Тогда, воспользовавшись условиями (1), можно определить коэффициенты уравнения (17) и записать частный вид реакции следующим образом:
R (θ)= |
|
1 |
(b2 −a2 )(8πθ−3π2 −4θ2 )ρ. |
(18) |
|
16 |
|||||
|
|
|
Представим решение задачи для бесконечно протяженного тяжелого толстостенного ортотропного цилиндрического тела со свободной от нагрузок внутренней поверхностью, которое находится на грунтовом основании, в виде суперпозиции решений: «классической» задачи Ламе (i = 0 ), задач о равновесии тяжелого цилиндра с жестко закрепленной внешней и свободной от напряжений внутренней боковой поверхностью (i =1) и, наконец, с жестко закрепленной внутренней поверхностью и заданным неравномерным давлением на внешней поверхности (i >1).
Закрепление цилиндров по внутренней или внешней поверхности приводит к появлению на этих участках дополнительных реакций, суммарный вклад которых
2π
Ra = ∫σ(rri>1)
0
cos |
|
θ + |
π dθ , R |
= |
2πσ(1) |
|
cosθdθ . |
(19) |
|||
|
|||||||||||
|
|
|
2 |
|
b |
|
∫ |
rr |
|
|
|
r =a |
|
|
|
|
|
|
|
r =a |
|
||
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
Поскольку исходная задача не предполагает каких-либо условий закрепления, необходимо скомпенсировать реакции (19). Для этого предположим, что неоднородная реакция на закрепленной поверхности может быть заменена эквивалентной равномерно распределенной нагрузкой интенсивностью
qa = Rπa и qb = Rπb .
Тогда граничные условия для рассматриваемой суперпозиции краевых задач можно представить следующей группой равенств:
σ(i=0) |
|
|
= −p |
− q |
, σ(i=0) |
|
|
|
= −p |
− q , |
(20) |
|||
|
|
|
|
|
||||||||||
rr |
r=a |
a |
|
|
a |
rr |
r |
=b |
|
b |
b |
|
|
|
u(i=1) |
|
|
= v(i=1) |
|
|
= 0 , σ(rri=1) |
|
= 0 , τ(riθ=1) |
|
= 0, |
||||
|
|
|
|
|
|
|||||||||
|
|
r=b |
|
|
|
r=b |
|
|
|
|
r=a |
|
|
r=a |
|
|
|
|
|
|
|
|
|
|
|
42

u(i>1) |
|
|
= v(i>1) |
|
|
= 0 , σ(i>1) |
|
|
= a cosi θ, τ(i>1) |
|
= 0 , |
|
|
|
|
|
|
|
|
||||||
|
|
r=a |
|
|
r=a |
rr |
|
r=b |
i |
rθ |
|
r=b |
|
|
|
|
|
|
|
|
|
||||
где ai = (−1)i (a2 −b2 )ρ i2 |
– коэффициенты разложения реакции ос- |
нования в тригонометрический ряд по окружной координате. Эти условия позволяют определить постоянные интегрирования общих ре-
шений (10), (14) и (15).
Обратим внимание на то, что граничные условия (20) позволяют скомпенсировать дополнительные реакции только в радиальном направлении. Несмотря на то что при закреплении внутренней или внешней боковой поверхности возникают усилия в окружном и осевом направлениях, попытка компенсации этих реакций привела бы к получению решения с точностью до константы, отвечающей либо за вращение тяжелого цилиндра как жесткого целого вокруг оси симметрии, либо за перемещение в вертикальном направлении.
На рис. 2 представлено распределение перемещений точек поперечного сечения монолитного железобетонного цилиндра на внутренней, внешней и серединной поверхности по окружной θ и обезразме-
ренной радиальной |
ρ = (r −a) (b −a) координатам. Были выбраны |
|
% |
следующие параметры геометрии тела и упругие постоянные железо-
бетона: a = 3,0 м и |
b = 5,5 м; E |
rr |
= 4 105 МПа, E |
= 8 105 МПа, |
|||||||||
|
|
|
|
|
|
|
|
|
θθ |
|
|
|
|
E |
zz |
= 5,33 105 |
МПа, |
G |
= 5,65 104 |
МПа, |
µ |
θr |
= 0,15 , |
µ |
rz |
= 0,075 |
|
|
|
|
rθ |
|
|
|
|
|
|
|
и µzθ = 0,375 .
Распределения перемещений по окружной координате носят немонотонный характер. Окружные перемещения являются нулевыми в точках, принадлежащих вертикальной диаметральной плоскости, в точках на внутренней, серединной и внешней поверхностях верхнего свода положительны, принимают свои максимальные значения при
θ = 90o , θ = 60o и θ = 45o. На нижнем своде в точках θ =125o , θ =115o
и θ =125o происходит смена знака, а θ =125o , θ =115o и θ =125o соответствуют точкам минимума для этих перемещений v на внутренней, серединной и внешней поверхности.
Радиальные перемещения на верхнем своде минимальны, возрастают по величине при увеличении окружной координаты. Эти перемещения на внутренней поверхности цилиндра всюду положительны,
43

на серединной поверхности принимают положительные значения
до θ =168o и имеют максимумы при θ =180o и θ =120o соответственно. На внешней поверхности зависимость u от окружной координаты
такова, что отрицательные значения имеют место при 0o ≤ θ ≤15o и 160o ≤ θ ≤180o, а максимальные значения достигаются при θ =115o .
Рис. 2. Распределение радиальных ( u,10−6 м) и окружных ( v,10−6 м)
перемещений в монолитном железобетонном цилиндре на внутренней, внешней и серединной поверхностях
Характер зависимостей перемещений точек, принадлежащих верхнему своду тяжелого цилиндра в направлении обезразмеренной радиальной координаты существенно изменяется при увеличении угла θ. Перемещения u монотонно убывают при θ = 0, монотонно возрастают при
θ = π2 и имеют точку локального минимума при θ = π
4 . Зависимости ν~ ρ% имеют локальные максимумы, положение которых смещается ко внутренней поверхности цилиндра при увеличении окружной координаты.
Обратим внимание на еще одну закономерность распределения радиальных перемещений. Производные u по окружной координате
44

∂u = − a rβ + |
a2 |
+ a ln r + a |
4 |
+ Br2 |
sin θ− |
||||||
|
|
||||||||||
∂θ |
|
1 |
|
rβ |
3 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|||
∞ |
|
ξ |
a |
|
ϕ |
a |
|
|
|
||
−∑i a1r i + |
|
2 |
|
+ a3r |
i + |
4 |
sin i |
θ |
|||
|
ξ |
ϕ |
|||||||||
i=2 |
|
|
r i |
|
|
r |
i |
|
|
принимают нулевые значения
∂u |
|
|
= |
∂u |
|
|
= |
∂u |
|
r=(a+b) 2 |
= |
∂u |
|
r=(a+b) 2 |
= |
∂u |
|
|
= |
∂u |
|
|
= 0. |
|
|
|
|
|
|
||||||||||||||||||
∂θ |
|
r=a |
|
∂θ |
|
r=a |
|
∂θ |
|
|
∂θ |
|
|
∂θ |
|
r=b |
|
∂θ |
|
r=b |
|
||
|
|
θ=0 |
|
|
|
θ=π |
|
|
|
θ=0 |
|
|
|
θ=π |
|
|
|
θ=0 |
|
|
|
θ=π |
|
в точках поперечных сечений тяжелых ортотропных цилиндров, принадлежащих вертикальной диаметральной плоскости. Этот результат не зависит от того, на каком расстоянии от оси вращения находится рассматриваемая точка.
Библиографический список
1.Матвеенко В.П., Шевелев Н.А. Равновесие упругого и вязкоупругого горизонтального тяжелого цилиндра // Прикладные задачи теории упругости и вязкоупругости. – Свердловск: УНЦ АН СССР, 1976. – С. 77–83.
2.Кузнецов Г.Б. Упругость, вязкоупругость и длительная прочность цилиндрических и сферических тел. – М.: Наука, 1979. – 112 с.
3.Термопрочность деталей машин / И.А. Биргер, Б.Ф. Шорр, И.В. Демьянушко [и др.]. – М.: Машиностроение, 1975. – 455 с.
Получено 29.11.2010
45
УДК 593.3
А.В. Зайцев, А.А. Фукалов
Пермский государственный технический университет
ЭФФЕКТИВНЫЕ МОДУЛИ ОБЪЕМНОГО СЖАТИЯ ДИСПЕРСНО-УПРОЧНЕННЫХ КОМПОЗИТОВ СО СПЛОШНЫМИ И ПОЛЫМИ АНИЗОТРОПНЫМИ СФЕРИЧЕСКИМИ ВКЛЮЧЕНИЯМИ
В рамках полидисперсных моделей механики композитов получены аналитические выражения для эффективных модулей объемного сжатия дисперсно-упрочненных материалов, изотропные матрицы которых армированы трансверсально-изотропными полыми или сплошными сферами различного диаметра, а на межфазных поверхностях выполняются условия идеального сопряжения. На основе полученных решений спрогнозированы эффективные характеристики полиэфирной смолы «Диэпокс 450» для различных наполнений этого связующего мраморной мукой, а также полимербетонов со сплошными сферическими анизотропными мраморными, диабазовыми и гранитными минеральными включениями произвольного размера.
Ключевые слова: дисперсно-упрочненные композиты, анизотропные полные и сплошные сферические включения, полидисперсная модель, идеальный контакт на межфазной поверхности, эффективные модули объемного сжатия; полиэфирное связующее «Диэпокс 450», наполненное мраморной мукой; полимербетон.
Рассмотрим двухфазный композит, армированный сплошными или полыми сферическими включениями различного диаметра. Придерживаясь гипотез полидисперсных моделей механики [1], будем считать, что каждая сферическая частица произвольного радиуса RA окру-
жена слоем связующего толщиной RB − RA таким образом, что отношения c = RARB для каждого армирующего элемента и h = R0
RА для
каждого полого включения являются постоянными величинами. Кроме того, будем предполагать, что линейно упругая однородная матрица изотропна, а материал, из которого изготовлены включения, – сферически трансверсально-изотропный (ось симметрии бесконечного порядка совпадает с осью r сферической ортогональной системы координат r, θ и ϕ с началом в центре составной частицы.
Рассмотрим произвольное включение с неизменным типом упругой симметрии составных частей (частицы армирующего наполнителя и слоя матрицы), находящееся в состоянии однородного гидростатиче-
46

ского растяжения, которое обеспечивается заданными на внешней поверхности равномерно распределенными радиальным давлением p
σrr |
|
r=R |
= p, |
(1) |
||
|
||||||
|
|
|
|
B |
|
|
или перемещением ξ: |
|
|
||||
ur |
|
r =R |
= ξ. |
(2) |
||
|
||||||
|
||||||
|
|
|
|
B |
|
|
Используя точные аналитические решения задач Ламе о равновесии толстостенных трансверсально-изотропных [2] и изотропных сфер [3], запишем выражения для радиальных перемещений ur и напряже-
ний σrr в полом включении и окружающем слое матрицы:
|
I |
|
|
n−1 2 |
|
|
B |
|
|
= 1 |
|
− |
|
|
|
|
B H |
|
|
|
|
|||
u |
r |
= A r |
|
+ |
|
I |
|
, |
σI |
|
rn |
3 2A H |
− |
|
|
I |
2 |
|
|
, |
||||
|
|
n+1 2 |
|
|
||||||||||||||||||||
|
I |
|
|
|
r |
|
rr |
2 |
|
|
I 1 |
|
|
r |
n+3 2 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
H1 = A11(2n −1)+ 4A12 , |
H2 = A11(2n +1)−4A12 , |
|
|
||||||||||||||||||||
|
|
uM |
= A r + BM , |
|
σM |
= 3K |
A − 4 |
GM BM |
, |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
r |
|
M |
|
|
r2 |
|
rr |
|
|
M M |
r |
3 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3)
(4)
где KM и GM – объемный и сдвиговой упругие модули. Здесь и далее
все величины, относящиеся к включениям и матрице, будут отмечены индексами I и M соответственно.
Для дисперсно-упрочненного композита со сплошными включениями в уравнениях (3) исключаются слагаемые с BI , наличие которых в решении приведет к сингулярности радиальных перемещений и напряжений в центре частиц армирующего наполнителя:
|
|
urI = СIrn−1 2 , |
σrrI = 1 |
СI H1 |
. |
|
|
|
|
|
|
(5) |
||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
2 r3 2−n |
|
|
|
|
|
|
|
||||
Содержащийся в равенствах (3) и (5) показатель анизотропии |
||||||||||||||||||
|
|
n = |
1 |
+ 2 |
A22 + A23 − A12 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
4 |
|
|
|
A |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
11 |
|
|
|
|
|
|
|
|
|
|
|
определяется модулями Юнга EI, |
E%I |
и коэффициентами Пуассона υI |
||||||||||||||||
и υI материала включений: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
% |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
% |
|
|
% |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
EI (1− υI ) |
|
|
|
|
|
|
|
EI |
|
|
2 EI |
|
|||||
A11 = |
, A12 |
= |
EIυI |
|
A22 = |
|
|
|
|
|
|
|||||||
|
|
, |
|
|
1 |
− υI |
E%I |
, |
||||||||||
F |
F |
F (1+ υI ) |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
% |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
47 |

|
EI |
|
2 |
EI |
|
|
2 |
EI |
|
A23 = |
|
|
υI + υI |
E%I |
, |
F =1− υI − 2υI |
E%I . |
||
F (1+ υI ) |
|||||||||
|
|
|
% |
|
|
|
% |
|
|
Рассмотрим наиболее простой случай, когда на межфазной поверхности выполняются условия идеального контакта
urI |
|
= urM |
|
|
, σrrI |
|
|
= σrrM |
, |
(6) |
|
r=RA |
|
|
r = RA |
|
|
r=RA |
|
r =RA |
|
|
|
|
|
|
а внутренняя поверхность полых сферических включений свободна от напряжений:
σrr |
|
r=R = 0 . |
(7) |
|
|||
0 |
|
Тогда из решения системы линейных алгебраических уравнений, которая получается при подстановке (4) и (5) в равенства (1) и условия (6), могут быть однозначно вычислены константы интегрирования:
С |
I |
= 2H |
C |
(4G + 3K |
M |
)R3 2−n , |
A = H |
C |
(H +8G |
M |
), |
|
|
M |
A |
M |
1 |
|
BM = −(H1 − 6KM )HCRA3 ,
HC = p 4GM (H1 − 6KM )c3 + 3KM (H1 +8GM ) −1 ,
если композит армирован сплошными сферическими включениями. В случае, когда включения являются полыми, подстановка (3) и (4) в выражения (1) и условия (6) и (7) позволяет записать систему алгебраических уравнений и решить ее относительно неизвестных постоянных интегрирования:
AI = 2H2HI (4GM + 3KM )RA32−n , BI = 2H1HI (4GM + 3KM )h2nRA3
2+n ,
A = H |
I |
(4H G − H H |
H |
|
), B = H |
I |
(H H |
H |
5 |
+ 3H K |
M |
)R3 |
, |
||||||||||||||||
M |
|
|
|
3 M |
|
1 2 |
|
5 |
|
|
M |
|
1 2 |
|
|
|
3 |
A |
|
||||||||||
H |
3 |
= 2 |
( |
H h2n + H |
2) |
, |
|
H |
I |
= p H |
7 |
− A2 H |
H |
5 ( |
4n2 −1 |
− H |
−1 , |
|
|||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
11 4 |
|
|
|
) |
|
8 |
|
|||||||||||
|
|
|
H |
4 |
= 4c3G + 3K |
M |
, |
H |
5 |
= h2n −1, |
|
H |
6 |
= c3 −1, |
|
|
|
||||||||||||
|
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H7 =16A12H5 (A12H4 −6GM KM H6),
H8 = 8A11{A12H4H5 +3H6KMGM 2n(h2n +1)− H5 }.
Если предположить, что армирующие элементы однородно распределены внутри композита, то рассматриваемая среда квазиоднородна и сферически трансверсально-изотропна. Поставим в соответст-
48

вие составным частицам, находящимся в условиях однородного гидростатического растяжения, эквивалентные включения, упругие модули которых являются эффективными деформационными характеристиками композита. Из равенства радиальных перемещений на внешних границах эквивалентных однородных включений
pRB ,
3Kθφ*
и составных частиц
uM |
|
r =RB |
= A |
R |
+ |
BM |
|
||||||
r |
|
M |
B |
|
2 |
|
|
|
|
|
|
|
RB |
|
|
|
|
|
|
могут быть найдены оценки сверху для аналогов эффективного модуля объемного сжатия для сферически трансверсально-изотропного дис- персно-упрочненного материала со сплошными
Kθϕ* = |
|
4G c3 |
(H − 6K |
M |
)+ 3K |
M |
(H +8G |
M |
) |
||||
|
M |
1 |
|
|
1 |
|
|
|
|||||
|
|
|
|
3 |
|
|
|
|
|
|
|
||
|
|
3 H1 +8GM |
− c |
(H1 − 6KM ) |
|
|
|
|
|||||
и полыми |
2H3 ( A12H4 − 6H6KMGM )− A11H4H9 |
|
|
||||||||||
Kθϕ* = |
, |
||||||||||||
|
|||||||||||||
|
|
|
3H1H2H6H5 + 3H3H10 |
|
|
|
|
(8)
(9)
H |
9 |
= H |
2 |
(1− 2n)+ h2nH |
(1+ 2n), |
H = 4G +3c3K |
M |
|
|
|
1 |
|
10 |
M |
включениями. Обратим внимание на то, что замена условий (1) на кинематические (2), обеспечивающие однородное перемещение внешней поверхности составной сферы, позволяет определить константы интегрирования для сплошных
С |
I |
= 2ξ(4G + 3K |
M |
)H |
C |
R3 2−n , |
A = ξH |
C |
(H +8G |
), |
|
M |
|
A |
M |
1 M |
|
BM = −ξ(H1 − 6KM )HCRA3 , HC = {RB H1 +8GM −(H1 −6KM )c3 }−1
и полых
AI = −2ξH2HI (4GM + 3KM )RA32−n , BI = −2ξH1HI (4GM + 3KM )h2nRA3
2+n , AM = ξHI [H1H2H5 − 4H3GM ], BM = −ξHI [H1H2H5 + 3H3KM ]RA3 ,
H11 = A11 −4A12 ,
HI = {RB H5H11(2H8 + H6H11)−4A11n(H10 + H10h2n + A11H6H5n) }−1
49

включений, а также получить оценки снизу для эффективных модулей K*θϕ из равенства радиальных напряжений на внешних границах эквивалентных однородных включений
|
|
|
|
|
|
ξ |
* |
|
||
|
|
|
|
|
|
|||||
|
|
σrr |
r =R |
= 3 |
|
Kθϕ |
||||
RB |
||||||||||
|
|
|
|
B |
|
|
|
|
||
и составных частиц |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
||||
σM |
|
r=RB |
= 3K |
A |
− 4 |
GM BM |
. |
|||
|
||||||||||
|
|
|||||||||
rr |
|
|
|
M M |
3 |
|
||||
|
|
|
|
|
|
|
|
RB |
||
|
|
|
|
|
|
|
Несложно показать, что эти оценки совпадут с полученными выражениями (8) и (9). Следовательно, полученные решения являются точными в рамках ограничений, используемых в полидисперсных моделях. Обратим также внимание на то, что при подстановке R0 = 0 в уравне-
ния (9), последние упрощаются и принимают вид (8).
В частном случае E%I = EI = 9GI KI (3KI + GI ) и υ%I = υI = (3KI − 2GI )
(6KI + 2GI ) из выражений (8) и (9) следуют эффективные модули объемного сжатия двухфазных дисперсноупрочненных композитов, изотропная матрица которых содержит однородно распределенные изотропные ( n =1,5) сплошные
K* = |
4vIGM (KI − KM )+ KM (3KI + 4GM ) = |
|
|
|
|||||
|
3vI (KI − KM )−3KI − 4GM |
|
|
|
(10) |
||||
|
KM + ( KI − KM ) |
vI ( 4GM + |
3KM ) |
|
|||||
= |
|
|
|
||||||
4G + 3K |
I |
+ 3v |
( K |
M |
− K |
) |
|
||
|
|
M |
I |
|
I |
|
|
и полые
K* = 4GI KI (1− h3)(4vIGM + 3K(M )+ G)M KM (1− vI )(4GI + 3h3KI )(11) 3KI h3 + 4GI (1− vI ) 1− h3 + 4GI (4GM + 3vI KM )
сферические включения. Обратим внимание на то, что в формулах (10) и (11) проведена замена c3 = vI (где vI — объемное наполнение),
а равенство (10) в точности совпадает с выражением, впервые записанным З. Хашином [1].
В качестве примера рассмотрим задачу прогнозирования эффективных модулей объемного сжатия полимербетонов – дисперсноупрочненных композитов на основе полиэфирных смол, армированных
50