1222
.pdfравна (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 и]) и т. д.