книги / Методы самосогласования механики композитов
..pdfТаким образом, на основе формул (3.103), |
(3.65) приходим |
|||||
к постановке относительно поля перемещений |
|
|
||||
~I |
~ |
~ |
~ |
~ |
(ξ) |
(3.104) |
ud |
(ξ) ≡ uii |
(ξ) = u11 |
(ξ) +u22 (ξ) |
+u33 |
||
|
d |
d |
d |
d |
|
|
первой локально-осредненной краевой задачи
∂ |
|
I |
|
∂ ~ I |
|
|
|
|
amndb |
(ξ) |
|
ud |
(ξ) |
= 0 |
|
∂ξn |
|
||||||
|
|
|
∂ξb |
|
|
с граничными условиями при ξ → ∞
ud |
=< εV |
εdb > ξb , |
~I |
/ |
/ |
где пульсации относительного изменения объема
εV/ (r) ≡ εii/ (r) = εV (r) − εV*
(3.105)
(3.106)
(3.107)
в произвольной точке r композита; εV* – относительное изменение макрообъема композита.
Аналогично (3.101)–(3.107) для расчета EV (ξ) (3.69) можем записать осредненную краевую задачу
∂ |
|
I |
|
∂ |
~I |
|
|
|
amndb |
(ξ) |
|
Uαβ |
(ξ) |
= 0 |
|
|
∂ξb |
||||||
∂ξn |
|
|
dϕψ |
|
|
с граничными условиями при ξ → ∞ вида
U~αβI = δαβIdbϕψξb dϕψ
относительно поля
~I |
~ |
~ |
~ |
|
~ |
|
Uαβ |
(ξ) ≡Uiiαβ |
(ξ) =U11αβ (ξ) +U |
22 |
αβ (ξ) +U |
33αβ (ξ) . |
|
dϕψ |
dϕψ |
|
dϕψ |
dϕψ |
dϕψ |
Поле упругих свойств aI (ξ) в (3.105) и (3.108)
amndbI (ξ) ≡ a11ii |
(ξ) = a1111 (ξ) + 2a1122 (ξ) |
|
mndb |
mndb |
mndb |
(3.108)
(3.109)
(3.110)
171
или
aI (ξ) = 9 |
|
VV (ξ)V + 6 |
|
VD (ξ)D |
(3.111) |
a |
a |
с учетом разложения (3.87), так как выполняются равенства (3.96)
V11ii =1 , D11ii = 0 .
Искомые поля условных моментов деформаций KV(ε) (ξ) (3.101)
рассчитываются через решение |
~I |
(ξ) |
локально-осредненной краевой |
|||||||||
u |
||||||||||||
задачи (3.105) и (3.106): |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
(ε ) |
|
%I |
|
%I |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
KVmn |
(ξ) ≡ εmn |
(ξ) = u(m,n) (ξ), |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
(3.112) |
|
|
(σ) |
|
|
|
% I |
|
I |
|
∂ |
%I |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
||
Kiimn |
(ξ) ≡ σmn (ξ) = amndb (ξ) ∂ξb |
ud |
(ξ), |
с учетом зависимостей (3.72), (3.74), (3.104) и (3.111). Аналогично
(3.112) вычисляются поля (3.102):
|
|
~I |
|
|
|
|
I |
~ I |
|
|
EVαβ |
|
(ξ), Siiαβ |
(ξ) . |
(3.113) |
||||||
(ξ) =Uαβ |
|
(ξ) = amndb |
(ξ)Uαβ |
|||||||
mnϕψ |
m) |
ϕψ,( n |
|
mnϕψ |
|
mϕψ,n |
|
|
Отметим, что скалярные поля условных моментов для величины относительного изменения объема K[εV ] (ξ) и для шаровой состав-
ляющей тензора напряжений |
|
|
|
|
рассчитываются |
следующим |
|||||||||||
K[σ] (ξ) |
|||||||||||||||||
образом: |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
(ε) |
|
|
|
~I |
|
|
~I |
(3.114) |
||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
K[εV ] (ξ) ≡ Kiijj |
(ξ) = εii (ξ) = |
εV (ξ) , |
|||||||||||||||
|
|
|
|
|
|
1 |
|
|
|
(σ) |
|
1 % I |
|
|
|||
|
K[σ] (ξ) ≡ |
|
|
(ξ) = |
(ξ) . |
|
|||||||||||
|
|
|
Kiijj |
|
σii |
(3.115) |
|||||||||||
|
9 |
9 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Аналогично (3.114) и (3.115) соответствующие смешанные моменты рассчитываются по формулам
172
|
|
|
|
|
|
|
(ε) |
(ξ) |
~I |
|
|
~I |
(ξ) , |
|
|||
|
|
|
|
|
|
|
|
|
|
||||||||
K[εV γ12 ] (ξ) ≡ 2Kii12 |
= 2ε12 |
(ξ) = γ12 |
(3.116) |
||||||||||||||
|
|
|
|
|
|
1 |
|
|
|
(σ) |
|
1 |
~I |
|
|
|
|
|
|
K[στ12 ] (ξ) ≡ |
|
|
|
(ξ) = |
(ξ) , |
|
|||||||||
|
|
|
|
Kii12 |
|
σ12 |
(3.117) |
||||||||||
|
|
3 |
3 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где τ12 и γ12 – соответствующие касательные напряжения и угол сдвига, например, в плоскости r1Or2 .
Когда точка ξ лежит в области f-й фазы включений υ( f ) , выполняются соотношения
K[σ] (ξ) = (K( f ) )2 K[εV ] (ξ) , K[στ12 ] (ξ) = K( f )G( f ) K[εV γ12 ] (ξ) ,
где K( f ) и G( f ) |
– модули упругости f-й фазы композита (3.84). |
||||||||||||||||||
Определим связь коэффициентов aVV( f ) |
и aVD( f ) разложения (3.86) |
||||||||||||||||||
~I |
(ξ) локально-осредненной краевой задачи (3.105), |
||||||||||||||||||
с решением u |
|||||||||||||||||||
(3.106). Из формулы (3.76) получим |
|
|
|
|
|
|
|
|
|||||||||||
( f ) |
|
|
|
|
1 |
|
~ |
|
|
|
|
|
1 |
|
~I |
||||
Aiisq |
< εsq εdb |
>= |
|
|
|
|
∫ |
uii |
|
(ξ)dξ = |
|
|
|
|
∫u( m,n) (ξ)dξ |
||||
υ |
|
|
|
υ |
|
|
|
||||||||||||
mndb |
|
|
|
( |
f ) |
υ( f ) |
( m,n) |
|
|
|
( |
f ) |
υ( f ) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( f ) |
|
|
|
|
|
|
|
1 |
|
~I |
|
|
|
|
|
|
|
|
|
Aiisq |
< εsq εdb >= |
|
|
|
∫ εmn |
(ξ)dξ |
(3.118) |
||||||||||
|
|
υ |
|
|
|||||||||||||||
|
|
mndb |
|
|
|
|
|
|
|
( f |
) υ( f ) |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
с учетом (3.104), (3.112). В (3.118) компоненты тензора |
|||||||||||||||||||
|
|
A( f ) |
= δ |
sq |
(a( f )V |
|
|
+ a( f ) D |
|
|
), |
(3.119) |
|||||||
|
|
iisq |
|
|
|
VV |
mndb |
|
VD |
mndb |
|
|
|
||||||
|
|
mndb |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
так как в разложении (3.86) имеем равенства |
|
|
|
|
|||||||||||||||
|
|
|
|
Viisq |
= δsq , Diisq = 0 |
|
|
|
|
|
для тензоров V и D (2.146).
173
Таким образом, левая часть равенства (3.118) может быть пре-
образована к виду |
|
|
|
|
>= (a( f )V |
|
)< ε |
|
|
|
|
|||||||||
|
|
A( f ) |
< ε |
sq |
ε |
db |
+ a( f ) D |
ε |
db |
> , |
(3.120) |
|||||||||
|
|
|
|
iisq |
|
|
VV mndb |
VD mndb |
|
|
ii |
|
|
|
||||||
|
|
|
|
mndb |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
следовательно, формула (3.118) примет вид |
|
|
|
∫ |
|
|
|
|||||||||||||
|
1 |
aVV |
) |
δmn |
< εV εV |
> +aVD |
) |
Dmndb < εV εdb >= |
1 |
|
εmn (ξ)dξ |
(3.121) |
||||||||
( f |
|
|
|
|
|
|
( f |
|
|
|
|
~I |
|
|
||||||
3 |
|
|
|
|
|
|
|
|
|
|
|
|
υ |
( f ) υ( f |
) |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с учетом (3.96), (3.119) и (3.120). Заменив в (3.121) свободные индексы mn последовательно на mm и, например, на 12 , получим
( f |
) |
|
|
|
1 |
~I |
|
|
|
|
1 |
|
|
( f |
) |
|
|
|
|
|
1 |
~I |
|
|||||
aVV |
|
< εV |
εV |
>= |
|
|
∫ εV |
(ξ)dξ |
, |
|
|
|
aVD |
|
< εV |
γ12 |
>= |
|
|
∫ ε12 |
(ξ)dξ |
|||||||
|
υ |
|
|
|
|
|
υ |
|
||||||||||||||||||||
|
|
|
|
|
( f ) υ( f ) |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
( f ) υ( f ) |
|
||||||
или в итоге придем к выражениям |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
∫ |
~I |
(ξ)dξ |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
εV |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
a( f ) |
= |
|
|
υ( f |
) |
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
(3.122) |
|
|
|
|
|
|
|
|
υ |
|
< ε |
|
ε |
|
|
> |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
VV |
|
|
( f ) |
V |
V |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
∫ |
~I |
(ξ)dξ |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
γ12 |
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
a( f ) |
= |
|
|
υ( f |
) |
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
(3.123) |
|
|
|
|
|
|
|
|
υ( f ) < εV γ12 > |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
VD |
|
|
|
|
|
|
|
|
|
|
где ~εVI ≡ ~εiiI – величина относительного изменения объема в точке ξ
расчетной области локально-осредненной краевой задачи (3.105), (3.106). Формулы (3.122) и (3.123) с учетом (3.114) и (3.116) можно записать в виде
|
|
∫ |
|
|
|
|
|
|||||||
|
K[εV ] (ξ)dξ |
|
|
|
||||||||||
a( f ) = |
|
υ( f ) |
|
|
|
|
|
|
, |
(3.124) |
||||
|
υ |
|
|
|
< ε |
|
ε |
|
> |
|||||
VV |
|
( f ) |
V |
V |
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|||||
|
|
∫ |
|
|
|
|
|
|||||||
|
|
K[εV γ12 ] (ξ)dξ |
|
|
|
|||||||||
a( f ) = |
υ( f ) |
|
|
|
|
|
|
|
. |
(3.125) |
||||
|
υ( f ) < εV γ12 > |
|
||||||||||||
VD |
|
|
|
|
174
Таким образом, коэффициенты aVV( f ) и aVD( f ) в качестве параметров самосогласования входят в скалярные поля соответственно aVV (ξ) и aVD (ξ) (3.88), (3.93), (3.94), (3.97) и (3.98), формирующих поле упругих свойств aI (ξ) (3.111) первой локально-осредненной
краевой задачи (3.105), (3.106); при этом коэффициенты aVV( f ) и aVD( f )
вычисляются по формулам (3.124) и (3.125) через решение этой же краевой задачи, что приводит к ее физической нелинейности.
Аналогично (3.101)–(3.125) могут быть определены поля услов-
ных моментов деформаций |
K |
12(εmn) (ξ) , напряжений |
K |
12(σmn) (ξ) |
и константы |
||||||||||
aDV( f ) и aDD( f ) разложения (3.86). Уравнение (3.64) представим в виде |
|||||||||||||||
|
|
|
∂ |
|
|
|
|
|
∂ |
~ |
|
|
|||
|
|
|
|
a12 sq (ξ) |
|
usq (ξ) = 0 |
|
||||||||
|
|
|
|
∂ξb |
|
||||||||||
|
|
|
∂ξn |
|
mndb |
|
d |
|
|
||||||
или |
|
|
|
|
|
|
|
|
|
|
|
||||
|
∂ |
|
∂ ~ |
|
~ |
|
|
||||||||
|
|
a1212 (ξ) |
|
|
|
u12 (ξ) + u21 |
(ξ) = 0 |
(3.126) |
|||||||
|
|
|
|
|
|||||||||||
|
∂ξn mndb |
∂ξb |
d |
|
d |
|
|
||||||||
с учетом симметрии (3.87) |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
a1212 |
= a1221 . |
|
|
|
|
||||
|
|
|
|
|
|
|
mndb |
|
mndb |
|
|
|
|
От (3.64), (3.65) перейдем к постановке относительно поля перемещений
|
~II |
|
~ |
~ |
(ξ) |
(3.127) |
||
|
ud |
(ξ) ≡ u12 (ξ) + u21 |
||||||
|
|
|
|
d |
|
d |
|
|
второй локально-осредненной краевой задачи |
|
|||||||
|
∂ |
|
II |
|
∂ ~II |
|
|
|
|
|
amndb |
(ξ) |
|
ud (ξ) = 0 |
(3.128) |
||
|
∂ξn |
|
||||||
|
|
|
|
∂ξb |
|
|
175
с граничными условиями при ξ → ∞
~II |
(3.129) |
||||
ud =< γ12εdb > ξb . |
|||||
Поле упругих свойств aII (ξ) в (3.128) |
|
||||
amndbII (ξ) ≡ a1212 (ξ) |
|
||||
|
|
mndb |
|
||
или |
|
||||
aII (ξ) = 3 |
|
DV (ξ)V + 2 |
|
DD (ξ)D , |
(3.130) |
a |
a |
с учетом разложения (3.87), так как выполняются равенства (2.146)
|
|
|
V |
|
|
= 0 , |
D |
|
|
= |
1 |
|
. |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
1212 |
|
|
1212 |
|
|
2 |
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Поля деформаций |
~II |
|
|
|
|
|
|
|
|
|
|
|
|
|
~II |
(ξ) , |
соответствую- |
||||
ε |
|
(ξ) и напряжений σ |
|||||||||||||||||||
щие решению для поля перемещений |
|
~II |
(ξ) |
(3.128) и (3.129), можно |
|||||||||||||||||
u |
|
||||||||||||||||||||
рассчитать по формулам |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
%II |
(ξ) |
|
% II |
(ξ) , |
|
|
|
|
|
|
(3.131) |
||||||||
|
|
ε |
= defu |
|
|
|
|
|
|
||||||||||||
|
|
~II |
|
|
|
|
|
II |
|
|
~II |
|
|
|
|
|
|
|
|
||
|
|
σmn |
(ξ) = amndb |
(ξ)εdb (ξ) . |
|
|
|
(3.132) |
|||||||||||||
Искомые поля условных моментов деформаций |
|
|
|||||||||||||||||||
|
|
|
|
|
|
(ε) |
|
|
|
1 ~II |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
K12mn (ξ) = |
|
|
εmn |
(ξ) |
|
|
|
|
||||||||||
|
|
|
2 |
|
|
|
|
|
|||||||||||||
и напряжений |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
(σ) |
|
|
~II |
|
|
|
II |
|
|
|
|
|
∂ |
|
~II |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
K12mn (ξ) |
= |
σmn |
(ξ) = amndb (ξ) |
|
ud (ξ) |
|
|||||||||||||||
∂ξb |
|
с учетом зависимостей (3.131) и (3.132).
Аналогично (3.114)–(3.117) моменты рассчитываются по формулам
176