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

книги / Методы самосогласования механики композитов

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

Таким образом, на основе формул (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

 

 

Отметим, что скалярные поля условных моментов для величины относительного изменения объема KV ] (ξ) и для шаровой состав-

ляющей тензора напряжений

 

 

 

 

рассчитываются

следующим

K[σ] (ξ)

образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ε)

 

 

 

~I

 

 

~I

(3.114)

 

 

 

 

 

 

 

 

 

 

 

 

KV ] (ξ) Kiijj

(ξ) = εii (ξ) =

εV (ξ) ,

 

 

 

 

 

 

1

 

 

 

(σ)

 

1 % I

 

 

 

K[σ] (ξ)

 

 

(ξ) =

(ξ) .

 

 

 

 

Kiijj

 

σii

(3.115)

 

9

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично (3.114) и (3.115) соответствующие смешанные моменты рассчитываются по формулам

172

 

 

 

 

 

 

 

(ε)

(ξ)

~I

 

 

~I

(ξ) ,

 

 

 

 

 

 

 

 

 

 

 

KV γ12 ] (ξ) 2Kii12

= 12

(ξ) = γ12

(3.116)

 

 

 

 

 

 

1

 

 

 

(σ)

 

1

~I

 

 

 

 

 

K[στ12 ] (ξ)

 

 

 

(ξ) =

(ξ) ,

 

 

 

 

 

Kii12

 

σ12

(3.117)

 

 

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где τ12 и γ12 – соответствующие касательные напряжения и угол сдвига, например, в плоскости r1Or2 .

Когда точка ξ лежит в области f-й фазы включений υ( f ) , выполняются соотношения

K[σ] (ξ) = (K( f ) )2 KV ] (ξ) , K[στ12 ] (ξ) = K( f )G( f ) KV γ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

12mn) (ξ) , напряжений

K

12mn) (ξ)

и константы

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ε)

 

 

 

~II

 

 

 

~II

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K[ γ12 ] (ξ) = 4K1212 (ξ) = 12

(ξ) = γ12 (ξ) ,

(3.133)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(σ)

 

 

~II

(ξ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K[ τ12 ] (ξ) = K1212 (ξ)

= σ12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ε)

 

 

 

~II

 

 

 

~II

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K[ γ12 εV ] (ξ) = 2K12ii

(ξ) = εii

(ξ) = εV (ξ) ,

(3.134)

 

 

 

 

 

 

 

 

1

 

 

 

 

(σ)

 

1 ~II

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K[ τ12 σ] (ξ) =

 

K12ii (ξ) =

 

 

σii

(ξ) .

(3.135)

 

 

 

 

 

 

 

 

3

 

3

Когда точка ξ лежит в области f-й фазы включений υ( f ) , тогда

 

 

 

K[ τ12 ] (ξ) = (G( f ) )2

 

 

 

 

 

 

 

 

 

 

 

 

K[ γ12 ] (ξ) .

 

K[ τ12 σ] (ξ) = G( f ) K( f )

K[ γ12 εV ] (ξ) ,

 

Определим связь коэффициентов aDV( f )

и aDD( f ) разложения (3.86)

с решением

~II

(ξ) локально-осредненной

краевой задачи

(3.128),

u

(3.129). Из формулы (3.76) получим

( f )

< εsq εdb

>=

1

~

A12 sq

 

 

u(12)

υ

 

mndb

 

 

 

( m,n)

 

 

 

 

( f ) υ( f

)

или

 

1

~II

 

(ξ)dξ =

 

 

u( m,n)

(ξ)dξ

2υ

 

 

( f ) υ( f )

 

 

 

 

( f )

 

1

~II

 

 

A12 sq

< εsq εdb >=

 

 

εmn

(ξ)dξ

(3.136)

2υ

 

mndb

 

( f ) υ( f )

 

 

 

 

 

 

 

с учетом (3.127) и (3.131). Левая часть формулы (3.136) может быть преобразована к виду

A( f )

< ε

sq

ε

db

>= 2 A( f )

< ε ε

db

>= A( f )

< γ

12

ε

db

>=

12 sq

 

 

1212

12

1212

 

 

 

mndb

 

 

 

 

mndb

 

 

mndb

 

 

 

 

 

= 12 (aDV( f )Vmndb + aDD( f ) Dmndb ) < γ12εdb >,

следовательно, формула (3.136) примет вид

3 aDV δmn < γ12εV

> +aDD Dmndb < γ12εdb

>= υ

εmn (ξ)dξ. (3.137)

1 ( f )

( f )

1

~II

( f ) υ( f )

177

Заменив в (3.137) свободные индексы mn последовательно на mm и на 12 , получим

( f )

 

 

1

 

~II

 

 

 

 

 

 

( f )

 

 

 

 

 

1

~II

 

aDV

< γ12εV

>=

 

 

 

εV

(ξ)dξ ,

aDD

 

< γ12

γ12

>=

 

 

γ12

(ξ)dξ

υ

 

 

 

υ

 

 

 

 

( f ) υ( f )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( f ) υ( f )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~II

(ξ)dξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

εV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a( f ) =

 

 

υ( f

)

 

 

 

 

 

,

 

 

 

 

 

 

(3.138)

 

 

 

 

 

 

υ

 

< γ

 

ε

 

>

 

 

 

 

 

 

 

 

 

 

 

 

DV

 

 

( f )

12

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~II

(ξ)dξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γ12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a( f ) =

 

 

 

υ( f

)

 

 

 

 

 

 

.

 

 

 

 

 

 

(3.139)

 

 

 

 

 

 

υ( f ) < γ12 γ12 >

 

 

 

 

 

 

 

 

 

 

 

 

 

DD

 

 

 

 

 

 

 

 

 

 

Формулы (3.138) и (3.139) с учетом (3.133) можно записать в виде

aDV( f )

aDD( f )

K[ γ12εV ] (ξ)dξ

= υ( f )

 

 

 

 

 

 

,

 

υ( f ) < γ12εV >

 

 

 

 

 

 

 

K[ γ12 ] (ξ)dξ

 

=

υ( f )

 

 

 

 

 

 

.

υ

( f )

< γ

12

γ

12

>

 

 

 

 

 

 

 

 

 

(3.140)

(3.141)

Таким образом, коэффициенты aDV( f ) и aDD( f ) являются параметрами самосогласования второй локально-осредненной краевой задачи (3.128) и (3.129), так как поле упругих свойств aII (ξ) зависит от

коэффициентов aDV( f ) и aDD( f ) (3.88), (3.95), (3.96), (3.99), (3.100),

(3.130), и вычисляются через решение этой же краевой задачи по формулам (3.140) и (3.141). Вторая локально-осредненная краевая задача (3.128), (3.129) так же, как и первая (3.105), (3.106), является физически нелинейной.

178

Отметим, что имеют место равенства для вычисления осредненных по области f-й фазы моментов:

< εV εV >f =

 

1

 

 

 

 

 

 

 

K[εV ] (ξ)dξ = aVV( f )

υ

 

 

 

 

 

 

( f ) υ( f )

 

 

 

 

 

 

< εV γ12

>f

=

 

1

 

 

 

 

 

 

 

 

K[εV γ12 ] (ξ)dξ = aVD( f )

υ

 

 

 

 

 

 

( f ) υ( f )

 

 

 

 

 

< γ12εV

>f

=

 

1

 

 

 

 

 

 

 

 

K[ γ12εV ] (ξ)dξ = aDV( f )

υ

 

 

 

 

 

 

( f ) υ( f )

 

 

 

 

 

< γ12 γ12 >f =

 

1

 

 

 

 

 

 

 

K[ γ12 ] (ξ)dξ = aDD( f )

υ

 

 

 

 

 

 

 

( f ) υ( f )

 

 

 

 

 

 

< εV εV >,

(3.142)

< εV γ12

> ,

(3.143)

< γ12εV

> ,

(3.144)

< γ12 γ12

>

(3.145)

с учетом формул (3.57), (3.124), (3.125) и (3.140), (3.141). Из (3.143)

и (3.144) следует свойство симметрии для коэффициентов

a( f ) = a( f ) ;

(3.146)

VD

DV

 

однако (3.146) будет выполняться лишь приближенно, что обусловлено принятой аппроксимацией (3.86) для тензора A( f ) .

Расчетные схемы первой (3.105), (3.106) и второй (3.128), (3.129) локально-осредненных краевых задач теории упругости – это одиночное включение с неоднородным переходным слоем в однородной неограниченной среде при заданном на бесконечности всестороннем растяжении или сдвиге. Упругие свойства этого одиночного включения, переходного слоя и неограниченной среды для первой и второй краевых задач различны и определяются соответственно полями aI (ξ) (3.111) и aII (ξ) (3.130). Упругие свойства среды не равны

эффективным упругим свойствам композита. Граничные условия, а именно величины всестороннего растяжения или сдвига среды на бесконечности для первой и второй краевых задач также различны и рассчитываются по соответствующим формулам (3.106) или (3.129).

179

Упругие свойства переходного слоя учитывают особенности случайного взаимного расположения включений в композите.

Уточнить решение краевой задачи (3.51), (3.52) для квазипериодических структур возможно, если представить искомое решение K (u ) (ξ1 , ξ) в виде разложения

 

 

 

 

 

 

(u ) (ξ1 ,ξ) =

 

(u ) p (ξ1 ,ξ) +

 

(u )o (ξ1 ,ξ),

(3.147)

K

K

K

где

 

(u ) p (ξ1 , ξ) – считающиеся известным решение задачи

 

K

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a p

 

(ξ)

 

 

 

 

 

(u ) p (ξ , ξ) = 0,

(3.148)

 

 

 

 

 

 

 

K

 

 

 

∂ξn

 

 

∂ξq

∂ξb

 

 

 

 

mndb

 

 

 

sd

1

 

 

 

 

 

 

 

 

 

 

 

ijsq

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sd(u ) p

=< εsq εdb > ξ(q1)ξb

при ξ → ∞

(3.149)

 

 

 

 

 

K

для соответствующей периодической структуры, и перейти к решению задачи

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

(ξ)

 

 

 

 

(u )o (ξ , ξ) + f

 

(ξ) = 0,

(3.150)

 

 

K

 

 

mndb

∂ξq

 

m

∂ξn

 

 

∂ξb

1

 

 

 

 

 

ijsq

 

 

(1)

 

 

sd

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sd(u )o

= 0 при ξ → ∞

 

 

(3.151)

 

 

 

K

 

 

относительно отклонений K (u )o (ξ1 , ξ) (3.147); поле объемных сил

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (1) (ξ)

 

ao

(ξ)

 

 

 

p (u ) (ξ

 

, ξ) ,

 

 

K

1

 

 

∂ξq

 

m

∂ξn

mndb

 

∂ξb

 

ij

 

 

 

ijsq

 

(1)

 

 

sd

 

 

 

 

 

 

 

 

 

 

 

 

поле отклонений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ao (ξ) = a(ξ) a p (ξ) ;

 

 

отметим требование выполнения равенства

 

 

 

 

< εpεp

>=< εε >

 

(3.152)

180