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

1222

.pdf
Скачиваний:
3
Добавлен:
15.11.2022
Размер:
12.36 Mб
Скачать

равна (1.10)

 

 

 

 

 

 

 

 

(е) = е°,

(е')

= 0.

 

(4.10)

Поэтому

 

 

 

 

 

 

 

(<0 = (Сс : е) = Сс :е 0= о ° ,

(4.11)

и, следовательно, первые два слагаемых в

(4.3) могут быть кон­

кретизированы:

 

 

 

 

 

 

 

=

в °: Сс : е°,

А2=

2 (р) : е°.

(4.12)

Упражнение

4.1.

Доказать,

что

из

формулировки

задачи

(2.4.16), (2.4.17)

следует

 

 

 

 

 

< р :е') = - < е ': С ‘ :е ').С )

(4.13)

Выберем тензор р постоянным внутри каждого компонента ра.

Тогда

 

~

<Р> = £

«аРа-

(4.14);

а = 1

 

 

Так как оценки (4.8) не должны

зависеть от

геометрии тела, то-

мы будем упругое тело сравнения считать неограниченным. То­

гда задачу

теории упругости (2.4.16), (2.4.17) для неограничен­

ной области

можно рассматривать

в силу

(4.14) как задачу о

действии

объемных источников в

упругом

пространстве. Вектор­

ное поле

и'(х), являющееся решением этой задачи, определяется

как сумма решений Кельвина. Поэтому выражение Л4 (4.2) мо­ жет быть легко найдено.

Итак, все величины правой части (4.3) известны:

Л

=

е°: Сс : е°, Л

= 2 £

vapa : е°,

 

 

 

а=1

 

^

Я

~

(е':р )

Я

 

 

At =

= £ va (£a )a :pa. (4.15).

а=1

~

 

а=1

При этом только в выражение Л3 входят характеристики иссле­ дуемого материала. Величина А становится функцией постоянных величин ра. Определив экстремум этой функции, находим оценки

(4.8). ~ Рассмотрим композит, каждый компонент которого изотропен.

Разложим тензоры р, е, е° на шаровую часть

и девиатор:

 

р = -L pi + р . е = — Ы + е, 80 =

-^-е»/ + е»,

(4.16)

Естественно в качестве тела сравнения принять также изотроп­ ный материал

Ccijkl =

+ bifilk)

(4.17)

Предположим, что исследуемый композит является упруго-пласти­ ческим. Тогда из (2.4.50) и (2.4.51) следует, что

оРи{ }

л , = ~ 2

 

 

 

 

E“a,(^ “ ))dp“a)] ’

<4Л8>

°

 

 

(a)

— его интенсивность, т. е.

гтде ра — шаровая часть тензора ра,

Ри

plf1=

(Ра : Ра)1/2, ~Ра =

ра.----- pfj

 

(4.19)

а функция elf1(plf1) является

обратной

к

функции plf1 = oia)(e„a))—

— 2^eia):

 

 

 

^«eS*1]-».

 

(4.20)

elf1(plf1) = [of’(elf1) -

 

alf1(elf1) =

2pa [1 -

coa (elf1)] e ? 1,

 

(4.21)

где (oa (elf1) — функция

пластичности А. А. Ильюшина. В

частно-

сти, если композит является упругим, то в (4.21)

следует поло­

жить (Оа = 0, и поэтому

 

 

 

 

 

 

 

c !f1= 2ра elf1,

elf

>) =

-2(^ а1 ю

■,

(4.22)

и величина Л3 (4.18) будет равна

 

 

 

 

 

 

 

°2

 

 

Ra. : Ra

~1

 

Л3 -------S

[

Р а

 

 

(4.23)

9(Ка- К с)

 

2(Ра—Н-О

-1’

 

 

а = 1

 

 

 

 

 

 

Величины Л! и Л2 для изотропного композита можно записать в виде

Ai = К < + 2р.‘е» : е»,

А, =

»<* (-f- Р А +

2pa : *» ). (4.24)

 

a = l

 

 

Найдем величину Л4 для

изотропного материала. Для

этого к

уравнениям (2.4.16)

 

 

 

 

(кс + р,с) grad div и' +

рсДи = — Div р

(4.25)

применим преобразование

Фурье

(приложение

IV). Обозначая

через U(со) образ Фурье

вектора

и'(х)у а через

Р (ш) — образ

Фурье тензора р(х), получим

где |со |— длина вектора ш.

Заметим, что решение уравнений

(4.25) не изменится от того,

что мы к тензору р прибавим посто­

янный тензор. Поэтому вместо тензора р можно ввести в (4.25) «нормированный» тензор р'\

 

р =

£' +

(р>,

<р'> =

0.

 

 

(4.27)

Пусть

Р' (со)— образ

Фурье

тензора р'.

Умножая уравнение

(4.26)

скалярно на вектор со, получим

 

 

 

 

 

(о•U = i-

т -Р' •(й

 

 

 

(4.28)

 

(^+2рс)|со|2

 

 

 

 

 

 

 

 

 

 

Подставляя (4.28) в (4.26), выразим вектор U:

 

 

 

и = -

-[—

Yfk (k-P'-k) +

P'-kl,

 

(4.29)

 

рс |0)|

 

 

 

 

 

 

где к — единичный вектор:

 

 

 

 

 

 

 

 

Xc+\ic

i>3Kc+Pc

=

1

(4.30)

 

 

Хс+2цс

3Kc+4\ic 2( 1—vc)

 

 

 

Теперь

путем обратного

преобразования

Фурье к

(4.29) можно

получить решение и'. Однако нам нужна более скромная инфор­ мация. Пользуясь обобщенной формулой Парсеваля, получим

2^4= ( Р : e’dvx=

j Р (ш): E(Z)dVo,

(4 .31)

где £ ( CD) — образ

Фурье

тензора

е'(х). Пусть

(со)— образ

Фурье Gradtt(jc). Из (4.29) видно

 

 

 

и

=

\(-лс1

P ' k ) k ( g ) k — (P' Д ) ® k]

 

 

( W, =

-+Г l

k,P'bl -

P'ikkkkj) ) .

(4.32)

Подставляя это

выражение

в (4.31)

силу симметрии тензора

Р в (4.31) вместо £(ы ) можно подставить U(со)), получим

fP:[(x'k-F-k)l®%-(P’:k)®k]dVa. (4.33)

Система уравнений для определения

ра будет нелинейной для

упруго-пластического композита:

~

£

ВаеРе=

 

 

 

(4.39)

Э=1

 

 

 

 

 

 

где

 

 

 

 

 

 

Ва — |

— M l — «Ь). если а =

(4.40)

°f 1 М э. если а=£р,

 

 

 

 

 

лричем

 

 

 

 

 

 

7 , Л -

W

l

_ J _

(4.41)

( Р “

 

* £ * >

 

'

Р<“ > ’

 

если материал упруго-пластическии, и

 

 

 

2 (Р“в>) =

2(ца- ц с)

(4.42)

 

если материал является упругим.

 

 

 

проверкой, что

Упражнение 4.2. Убедиться

непосредственной

решением системы уравнений (4.37) является

 

(Р) =

3

\

т

 

 

(4.43)

0

1+аст

 

 

где

(4.44)

Упражнение 4.3. Убедиться непосредственной проверкой, что решением системы уравнений (4.39) для упругого композита яв­ ляется

 

 

<£> =

п

е°,

 

 

 

 

(4.45)

 

 

1-\-Ьсп

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

a o J -iP T - M

" '

 

С

 

(4-46>

0 = 1

 

 

 

 

 

 

 

 

Обозначим корни

систем уравнений (4.37)

и

(4.39),

соответ­

ственно через (ра)0. (£а)о-

(Для

упругого

случая

они

собраны в

комбинации <(р)0>

(см.

(4.43))

и ((р)о)

(см.

(4.45)).

Подстав­

ляя эти корни в выражение & или А (4.1*), получим

 

 

а = - = m l + 2р£ : £ +

(-2- ((р)о> в. +

2 <(р)0> : * )

-

Неравенство (4.4) будет выполняться, сравнения выбрать компонент композита

ми, например Kq,

Обозначим

 

К" =

Kq +

т

2ц" =

9 1 -\-аст

 

если в качестве тела с наименьшими модуля­

2ц, +

П

(4.55)

1 -\-Ьсп

Неравенство

(4.5) будет выполняться, ес­

 

ли в качестве тела сравнения выбрать ком­

 

понент композита с наибольшими модуля­

 

ми, например К\, щ. Обозначим

 

 

К' = Кг

9(1 + а ст)

2р' = 2рх + -

1+М

 

 

 

 

 

 

 

 

 

 

 

(4.56)

 

Сравнение выражений (4.53) и (4.54) по­

 

казывает, что из (4.8) следует

 

 

 

К ' < К * < К р '< р * < р " .

(4.57)

 

Это и есть так

называемая

вилка

Хаши-

 

на— Штрикмана.

показана

зависимость величин

К' и К" от

На рис.

13,

где

объемной

концентрации y = v x для

двухкомпонентного

композита

(/Ci> А^2) ,

эта

вилка

изображается

заштрихованной

областью.

Для сравнения на этом же рисунке изображена вилка Фойгта— Рейсса.

Вилку Хашина— Штрикмана не удается сузить, если не учи­ тывать геометрии композита, хотя для многих композитов ц она

оказывается

достаточно широкой. Для

примера

в

табл. 1

и 2

указаны

значения

К',

К"

заодно и KF, KR) для

двухкомпо­

нентного

композита (/Ci> /С2)

при

различных

объемных

концен­

трациях у. Все величины /С',

К", KF, KR считаются

безразмерны­

ми и отнесены к /С2- В табл.

1 /Ci = 2,

в табл.

2 К\ = 20;

в обеих

таблицах vi = v 2= 1/3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица

1

Nsss\

0

0.1

0.2

0,3

0,4

0,5

0,6

0,7

0,8

0,9

 

 

 

 

 

 

 

 

 

 

 

 

 

KF

1

1,10

1,20

1,30

1,40

1,50

1,60

1,70

1,80

1,90

2

К”

1

1,07

1,14

1,22

1,31

1,40

1,50

1,61

1,73

1,86

2

К'

1

1,06

1,13

1,20

1,29

1,37

1,47

1,58

1,71

1,84

2

к*

1

1,05

1,11

| 1,18

1,25

1,33

1,43 | 1,54

1,67

1,81

2

X

0

0.1

0 ,2

0 ,3

0,4

0,5

0,6

0,7

0 .8

0 ,9

 

 

 

 

 

 

 

 

 

 

 

 

KF

1

2,90

4,80

6,70

8,60

10,50

12,40

14,30

16,20

18,10

20

К”

1

1,74

2,60

3,58

4,73

6,10

7,74

9,76

12,30

15,58

20

К'

1

1,15

1,34

1,58

1,86

2,30

2,88

3,77

5,30

8,54

20

 

1

1,10

1,23

1,40

1,61

1,90

2,32

2,98

4,17

6,90

20

Упражнение 4.4. Показать, что для упругого двухкомпонент­ ного композита (/Ci>/Сг, M-i>

 

 

XlK2+ “

(1--Y)^2^2+

g YP2X1

 

 

К' =

Xi+

 

Рг—Y(-Ki—X2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.58)

 

 

Ki ^X2+

g

YPij +

~ (1—Y)PIX2

 

K" = -----------------1------------------------------

 

 

 

 

/Ci+ —

I-I2—Y(/^i—/^2)

 

W =

И-2 +

1

 

 

6 (X 2+ 2p 2) (1 -Y )

*

 

 

 

 

 

 

 

 

 

Pi—[fi2

5

(3/C2-f-4jx2)|J’2

(4.59)

 

 

 

 

 

 

 

 

 

 

 

P" =

Pi +

 

 

1 -Y

 

 

 

1

,

_6_

(/CI+2HI )Y

 

 

 

 

 

 

 

 

 

 

 

p2—Pi

 

5

(3/Ci+4/ii)/ii

 

Упражнение

4.5.

Показать,

что

при

pi = p2

(тело Хилла из

(4.58)) следует К' = К"-

 

 

 

 

 

 

 

 

§ 5. Некоторые

методы

определения

эффективных характеристик

Для приближенного определения эффективных характеристик существует много методов. Один из самых простых методов — метод вириального разложения — применим в случае, когда кон­ центрация одного из компонентов мала. Метод основан на разло­ жении эффективных тензоров модулей упругости и упругих по­ датливостей в ряд по концентрации компонента (если она доста­ точно мала).

Вметоде самосогласования принимается, что каждый компо­ нент имеет специальную форму (чаще всего форму эллипсоида или шара) и рассматривается как включение, при этом связую­ щим служит материал с искомыми эффективными свойствами. Этот метод приводит к успеху благодаря результату Дж. Эшелби, который показал, что внутри шарового упругого поля возни­ кает однородное напряженное состояние.

Втеориях смесей предполагается, что в каждой точке среды

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

Существует много других методов определения эффективных характеристик среды (а также определения микроперемещений и микронапряжений), однако самым распространенным методом, пожалуй, является метод теории случайных функций.

В этом методе тензор модулей упругости считается случайной функцией, представимой в виде суммы статистически среднего тензора модулей упругости и тензора, описывающего флюктуационные добавки. Как правило, принимается гипотеза эргодичности: среднее по объему совпадает со средним статистическим. (Прав­ да, здесь объем, по которому совершается осреднение, связан с характерным размером неоднородности, и поэтому средние вели­ чины, вообще говоря, зависят от координат).

Запишем

 

Сцы М = {Cijki) + C jki (х),

(5.1)

где тензор С' описывает флюктуационные добавки, а тензор ( С ) от координат не зависит (так называемая статистически одно­ родная среда).

Пусть, например, требуется решить задачу теории упругости с заданными на границе перемещениями

[({Cijki) + Cijki) ((uk) j +

Uk,i)],j =

0,

(5.2)

 

({щ) +

tii) |s =

u°i,

 

 

(5.3)

где {и) — статистическое среднее искомого

поля

перемещений,

и' — соответствующие

флюктуации. Решить

задачу

(5.2), (5.3) со

случайным тензором

модулей

упругости — это

значит по задан­

ным статистическим характеристикам механических свойств сре­ ды: среднего тензора модулей упругости (Ci!kl) y тензора корре­ ляционных функций второго порядка (Cijki Cmnpq) и т. д. найти статистические характеристики поля перемещений: среднего век­ тора перемещений ( щ) , корреляционные функции второго по­ рядка (tii и]) и т. д.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]