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

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

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

531

При этом Р -базис состоит из шести четырехвалентных тензоров, которые определяются выражениями

P<ifiXft = @а(Х @(1)Р > ^сфХц ~ ^ctfi^Xft > ^afiXft = ^ а / ПХГПц > ( П 1 .1 . 8 )

Р^Х„ = та™Ам > Ком =

= тат^Хт„

Элементы Р и Е -базисов (П1.1.8) и (П 1.1.1) связаны со­

отношениями

 

PJ= а*Е*, Ej = -1 )* Р*,

(П1.1.9)

где матрицы а и а~х имеют следующую структуру

1 0

0

0

-2

1

 

1 0

0

0

2 1

0

1 -1 -1

0

1

 

0

1 1 1

0

1

0

0

1 0

0

-1

(

-1\9 0 0 1 0

0

1

0

0

0

1 0

-1 .

( «

) = 0

0

0

1 0

1

0

0

0

0

1 -1

 

0

0

0

0

1

1

0

0

0

0

0

1

 

0

0

0

0

0

1

(П1.1.10)

Шесть тензоров Р' также образуют замкнутую алгебру от­ носительно операции умножения (П1.1.2). Таблица умноже­ ния этих тензоров имеет вид

 

р 1

р 2

р 3

р 4

р 5

рб

 

 

р 1

р 1

р 2

р 3

0

0

0

р 2

р 2

2Р 2 2Р 3

0

0

0

р 3

0

0

0

р 2

0

р 3

р 4

Р 4

2Р4

2Р б

0

0

0

р 5

0

0

0

0

Р 5/ 2

0

р 6

0

0

0

р 4

0

рб

(П1.1.11)

532

Заметим, что тензор Р х- Р 2/2 ортогонален всем элемен­ там базиса Р', за исключением Рх, т.е.

Рх(р х- \ Р 2) = {р х- \ Р 2)р х = (р х- \ Р 2), Р‘ (РХ- \ Р 2) =

= (р х- \ р 2)р = 0, i * i , (р х- \ р 2)(р х- \ р 2)= (р х- \ Р 2).

(П1.1.12) Очевидно, что линейные оболочки Е и Р -базисов совпа­

дают.

Пусть тензор А в Р -базисе имеет вид

А = а,Р2 +а2(р х- \ Р 2)+ а3Ръ +а4Р4 +а5Р5 +а6Р\

(П1.1.13)

где а( (/ =1,2,...,6) - скалярные коэффициенты. Тогда обрат­

ный тензор А~х (АА~Х=А~ХА х) определяется соотношени­ ем

^ -Р * + — Р5+ ^ - Р 6,

Л

а5

Л

А = 2{а,а6- а 3а4).

 

(П1.1.14)

Пусть тензоры А и В представлены в Р -базисе аналогич­ но (П 1.1.13). Произведение этих тензоров АВ (свертка по двум индексам) определяется соотношением, которое следует из (П1.1.11) и (П1.1.12)

АВ = {2а,Ъ, +а3Ь4)Р2 +a2b2(p x- \ Р2)+ {2а,Ъ3 +а3Ь6)Р3+

+(2 а4Ъ, + а6Ь4)Р* + ^а3Ь5Р^ +{а6Ь6+2 a4b3)b6.

(П1.1.15)

3°. в-базис. Пусть 6^ - проектор на плоскость, ортогональ­ ную вектору т (П1.1.7), а п - единичный вектор, лежащий в этой плоскости (т-п=0). Определим элементы 0-базиса соот­ ношениями

^сфЛ/i =

@ti)P > ^сфХц =

=

> ( П 1. 1. 16)

533

Фф* = П*ПА » > &фАм= ва)(х”м)” (Л»

^ = W /*x"M■

Приведем таблицу умножения элементов 0-базиса:

 

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

20

20

0

0

0

0

0

0

0

0

0

0

0

0

20

20

0

0

0

0

0

0

0

0

(0+ 0)/2

0

0

0

0

0

0

0

0

 

 

 

 

 

 

(П1.1.17)

4°. R -базис. Этот базис составлен из следующих пяти тензо­ ров:

Кр* = W * A ^ , К2фхм= тат^ПхтИ> Крх» = ™а™РхП» ,

K f iX M = n j 'f P x ™ » , K fS X n = n a )n (X m M) m (j}>

(П1.1.18)

где п и т - ортогональные единичные векторы. Свертка тен­ зоров этого базиса по двум индексам характеризуется следую­ щей таблицей

 

R'

R2

R3

R4

R5

R'

R1

0

0

R2

0

R2

0

R2

R3

0

0

R3

R3

0

0

R2

0

R4

0

R4

0

0

0

R5

0

0

0

0

R5/2

(П1.1.19)

534

5* Формулы осреднениядля элементов Е,Р,$ и R -базисов.

При решении задачи осреднения стохастически неоднородной среды возникает необходимость вычисления средних по еди­ ничной сфере или единичной окружности от элементов вве­ денных выше базисов. Приведем значения соответствующих

интегралов по единичной сфере Q ,:

(Е?{п)) = -±-\е(п)(Х1я, / = 1,2,....6,

(П 1.1.20)

(Е') = Е\ (Ег)= Е \ (Е'(п)) = (Е‘ (п)) = \Е\

( £ ’ (»)) = 1 £ ‘ > ( £ '(» ) ) = п ( 2 £ '+ £ 2),

(P*(l»)) = j-Jf'tm Jd n ., / = 1,2, ,6, (П1.1.21)

{P ' M ) = M e * +1E')- (^ W ) = - ft ( £l +3£ ! ).

(Р г(т)) = (Р*{т)) = £ (2 Е 2-Е '),

(/•!(ш)) = ^ £ ' - $ Е г), (р‘ (т)) = ^(2Е ' + £ ! ).

Элементы в и R-базисов осредним по вектору И, лежаще­ му в плоскости, ортогональной вектору т. Указанные сред­

ние, которые равны интегралам по единичной окружности /, в этой плоскости, следующим образом выражаются через эле­ менты Р -базиса:

(в{т,п)) = ^-\в{т,п)с11п, i = 1,2,...,6, (П1.1.22)

I h

(&{т,п)) = Рх{т), (&{т,п)) = Р2{т), (&(т,п)) = \ Р 2{т),

(&(т,п))=^Р2(т), (0(m,n))=ipiW , (&(п))=\{Р2(т)+2Рх(т)\

535

(Rl{m,n)) = ^ - jR i(m,n)dln, / = 1,2,...,5, (П1.1.23) n

{ R ' M ) = $ (p 2{m) + 2Pl{m)), ( R2{myn)) = P6{m),

(JR3(w ,/j))=|JP4(w), {R\m,n))=\p\m), (R5(m,n))=^P\m).

6* Двумерное пространство. Рассмотрим аналог Е и Р -ба­ зисов в двумерном случае. Начнем с Е -базиса, который имеет вид (П1.1.1), но греческие индексы принимают значения 1,2. Существенно, что в двумерном случае не все из шести тензо­

ров Е'(П) являются линейно независимыми и связаны сле­ дующим линейным соотношением

Ех- Е 2 + Е 3 +Е 4 -2 Е 5 = 0.

(П 1.1.24)

Поэтому линейно независимыми в плоском случае оказы­

ваются пять из шести тензоров Е\п) вида

(П 1.1.1) (любые

четыре из первых пяти и Е6{П)).

Для Р -базиса линейная зависимость тензоров (П1.1.8) оче­

видна, поскольку тензоры Р' и Р 2в двумерном случае совпа­ дают. Таблица умножения пяти линейно независимых элемен­ тов Р -базиса имеет вид

 

р х

р 3

Р4

р 5

рб

 

 

р 1

р х

р 3

0

0

0

р 3

0

0

р 1

0

р 3

Р4

Р4

рб

0

0

0

р 5

0

0

0

Р5/2

0

рб

0

0

р 4

0

рб

(П 1.1.25) Произведение двух элементов линейной оболочки Р -бази­

са и В)

536

А = ахР' +а3Р3 +а4Р4 +а5Р5 + а бР б,

(П1.1.26)

В = &.Р1+63Р 3 +Ь4Р4 +65Р5 + г>6р б,

где , Z). - скалярные коэффициенты, определяется соотношением

АВ = (ахЬх +а3Ъ4)Р1+(аг,63 + a 3Z>6)P 3 +

(П1.1.27)

+(«4*i +а6Ь4)Р* +\аъЬъРь+{а4Ь3 +а6Ь6)Р6.

Отсюда следует выражение для обратного тензора А

А = аха6- а 3а4.

(П1.1.28)

Связь между элементами Е и Р -базисов в двумерном слу­ чае определяется следующим образом:

1е '=Р1+2Р5+Р6, е 2=р '+р 3+р 4+р \

\е 3=р 3+р 6, е *=р *+р \ E S=P 5+P 6} E 6=P \ (П11-29)

1>! = Е '-2 Е 5+Е\ Р3 = Е3- Е \

} Р 4 = Е * -Е \

Р5= Е 5- Е 6, Р б = £ б.

<П1Л-3°)

Наконец, формулы осреднения элементов Е и Р -базисов

в двумерном случае имеют вид

 

(Е‘(п)) = 2^ J Е‘{п)ап, / = 1,2,...,6,

(П1.1.31)

(Е') = Е\

(Е3) = Е\ (Е 3(П)) = (Е*(п)) = ±Е \

(Е5{п)) = $Е\ {Е6{п)) = $(2Е' +Е 2),

(Р<Ц = A J pi(” )dln, i = 1,3,4,5 ,6,

(П 1.1.32)

537

(/>'W) = i ( £ 2+2£'), {Р >(п)) = {Р‘ (п)) = $(ЗЕ '-E '),

(P s(n)) = b (2 E '-E 2), (Р‘ (К)) = Ц е ! +2Е'У

П1.2. Изотропный шестивалентный тензор

При решении задачи со скалярным потенциалом в П.2.5 в случае квадратичного внешнего поля возникла необходимость в операциях с шестивалентным изотропным тензором, сим­ метричным по двум парам тензоров. Тензоры такой структуры представляются в виде линейной комбинации следующих тен­ зоров:

Q aglpfipi

^afi^Xfi^p t 5 QaXpftpr

>

QaXpfipx ~ ^■^a(X^p)P^pz->

QaXpfipx ~ ^^a(p^T)fi^Xp ’ (П 1.2.1)

QlCaXpfipt *

^^a^X^p)(j>^t)p >

QiaXpfipT

T)a '

Эти тензоры образуют замкнутую алгебру относительно операции умножения - свертки по трем индексам

Правило умножения тензоров 0-базиса (П1.2.1) представ­ лено в таблице.

 

Q1

Q2

е 3

 

Q4

Q5

Q6

 

0*

30 1

20 1

2 0 1

 

3Q4

2Q4

2Q4

 

Q2

20'

2Q2

2 0 3

 

2Q4

2Q5

2Q6

Q3

З03

20 3

2 0 3

 

30 5

2Q5

2Q5

Q4

2 0 1

20 4

80 1

 

2 0 4

80 4

4Ql +2Q4

Q5

2 0 3

20 5

80 3

 

2Q5

80 5

4 0 3+ 2

0 5

Q6

2 0 3

2Q6

40 1+ 2

0 3

20 5

4 0 4 + 2 0 5

4Q2+2Q6

538

Допустим, что два тензора А и В принадлежат линейной

оболочке -базиса, т.е.

A = Z a iQ ',

B = Z biQ‘ ,

(П1.2.3)

i=i

<=i

 

где at и br скалярные коэффициенты. Произведение этих тензоров может быть найдено с помощью приведенной табли­ цы и имеет вид

АВ = (ахЬх+2афх+2а4Ь2+4 a6b3)Ql +2(а2Ь2 +a6b6)Q2+

“Ьci3bx 2а^Ь2+ 2eg(6| b3) +

(П1.2.4)

+(яД + 2ci2bA+2а4Ь4+4 a6b^)QA+

+[2а2Ь5 +а3Ь3+2а5Ь4 +2ал{Ъ4 +bb)]Q5 + 2[a2b6+a6{b2 +b6)]Q6.

Здесь обозначено

 

Ьх= ЗЬХ-¥2Ъ2+ 2 b3, b2 —bx+ 4 Ь3+ 2 Ь^,

(П1.2.5)

Ь3= Щ + 2bs + 6, Ьл= Ъ2 + Ъл+ 4Ь5+ Ь6.

Используя эти выражения, можно найти шестивалентный

тензор, обратный заданному ( В 1В = jQ 2):

В-' = - З - Й в , - ь а )о ' + ^ 7 - Q 1+ Х - Й 5 . - b.B,)Q' +

Д2 4A J А2

+ 2 - ( » ,в , - * ;в г) е 4+ 3 - ( а д - * ;в 4 ) е 5- 1 7 - е ‘ (П1.2.б)

Д2 Д2 4 Д,

Здесь обозначено

Д1= Ь2(Ь2 + Ь6) —2Z>6 , Д2 —bxb4Ъ2Ь3,

В\ —4д^[^1 (Ь2 2Ь3Ь6], В3= 4д^ {b2b3—bxb6),

539

A

+ *«)-2 *А ], в , = Л - ( * А - » А ) -

П.1.3. Матричные представления элементов Е ,Р ,в и Л-базисов

Для конкретных вычислений удобно воспользоваться изо­ морфизмом между четырехвалентными тензорами, симмет­ ричными по первой и второй паре индексов, и матрицами размерностью (6x6) [23]. При этом двухвалентному симмет­

ричному тензору qap соответствует шестимерный вектор qt

(/=1,2,...,6) с компонентами

Я\ = Я\\ > Яг = Ягг > Яг = Ягг > Яа = <7i2 > Яг = ((в > Яб = Ягг ■

(П1.3.1) Матричные представления элементов £ -базиса ПЫЛ

имеют следующий вид

2

0

0

0

0

0

1 1 1 0 0 0

0

2

0

0

0

0

1 1 1 0 0 0

0

0

2

0

0

0

1 1 1 0 0 0

0

0

0

1 0

0

0 0 0 0 0 0

0

0

0

0

1 0

0 0 0 0 0 0

0

0

0

0

0

1

0 0 0 0 0 0

« I 2

«22

« 3

« 1 « 2

» 1 « 3

« 2 « 3

п\

П\

«з2

0

0

0

« I 2

п \

«з2

« 1 « 2

« 1 « 3

» 2 » 3

п\

«22

«з2

0

0

0

« I 2

п \

«з2

щ п г

« 1 « 3

« 2 « 3

п\

А

п\

0

0

0

0

0

0

0

0

0

п\пг

«1 « 2

« 1 » 2

0

0

0

0

0

0

0

0

0

« 1 « 3 « 1 « 3 « 1 « 3 0 0 0

0

0

0

0

0

0

« 2 « 3 « 2 « 3 « 2 » 3 0 0 0

540

 

4п\

0

0

2пхп2

2и,и3

0

1

0

4« 2

0

2пхп2

 

0

2и2п3

0

0

4«з

0

2«,«3

2и2и3

Е5 =

2пхп2 1пхп2

0

и,2 +и2

 

 

 

4

« 2« 3

« 1 « 3

 

2пхпг

0

2пхпг

п2пг

2

2

w,w2

 

«1

+ « 3

 

0

2п2пг 2«2W3

и,«з

« 1 « 2

« 2 + « 3 2

«I4

п]п\

п]п\

п\п2

п\пг

« , 2« 2 « 3

»l4 2

 

п\п\

пхп\

щп\щ

п\пг

Е6 = « > з 2

п\п\

 

пхп2п\

щп\

щп\

 

П\П\

пхп2п\

2

2

п\п2пг

п\пгпг

 

пхп2

п \п г

щп\щ

пхп\

п]п2п2

п\п\

пхп2п\

п]п2Щ

 

п2п\

щп\щ

пхп2п\

п\п\

Здесь пх,п2,пъ - компоненты вектора п в фиксированной декартовой системе координат.

При записи матричных аналогов элементов Р -базиса (П1.1.8) будем считать, что вектор т направлен вдоль оси х3 декартовой системы координат хх,х2,хг. При этом элементы

Р -базиса принимают вид

 

1

0

0

0

0

0

 

 

1

1

0

0

0

0

 

0

1

0

0

0

0

 

 

1

1

0

0

0

0

Р' =

0

 

0

 

0

0

0

0

7)2

 

0

 

0

 

0

0

0

0

0

0

0

1

0

0

, р

=

0

0

0

0

0

0

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

0

0

0

0

 

 

0

 

0

 

0

0

0

0

 

0

 

0

 

0

0

0

0

 

 

0

 

0

 

0

0

0

0

Соседние файлы в папке книги