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

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

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

271

Р‘ (т) = р,Рг +р2( р ' - \ Р г)+ р ,(Р '+ Р ')+ р ,Р >+ptP‘ ,

N

й = ( Е ( * |- * .) 0 + ^ )Й г-^-,) • “i = * , (5-12.7)

к1= 1

Л = U f (и, -M„)[(l+21?+21? -S i, )+*>’ f t +Tt X ^ - a J ,)]j ,

A = ( i ( ', - 0 ( i +i ;X ^ - ^ ,) ) . р > = ( ^ 1 ж < № - ё ? А

i 1=1

P« = ( £ [ » , -n . +2(1,- l M ] ( s l - a l,fj,

где постоянные Yj определяются из решения задачи для оди­ ночного цилиндрического слоистого включения, алгоритм по­ строения этих постоянных описан в § 2.11.

Тензор С* (5.12.2) в том же базисе Р’(т) определяется выражениями

С* = КР2 + 2т.(р' - i p 2) + /,(p 3 + P 4) + 4/U>5 +п.Р6,

fc =*e+_ J B L _

п и = т Л ^

L=L+- РРъ

1-4РРА

2(1-РР2Ь2У

* ° I"4Р Р А ’

Р Р

ъ _____ _

, „2

4 /foi

 

, п.=по+рр6+р

1-4РРА

Р*=Ро + 4 - рр5ьь

 

1 -ж .

 

1

(5.12.8)

*,=

 

2 я

4/^„

 

 

В случае однородных волокон (JV = 1) эффективные коэф­ фициенты упругости композита принимают вид

272

р(т-т0)

,Р.=Я

 

2р (м - я )

 

пк=та 1+

2я+(1 -р)(м~М.)

 

 

 

 

, , . „ ( к к , , .

 

 

^

 

^ ( к ) 2

1=1+PI— , к*=К+Р--- т--- ,n.=n+p(n-nj-p(\-p)—— ,

А =ркА '-р)к+р ., Д =

^

 

 

(5.12.9)

 

 

2по

 

 

 

 

 

 

На рис. 5.22 представлено

 

сравнение

рассчитанных

по

 

(5.12.8) значений

эффектив­

 

ных

параметров

упругости

 

волокнистого

композита

с

 

результатами

экспериментов

 

[167].

В [167]

исследовались

 

композиты

с

изотропной

 

матрицей

из

эпоксидной

 

смолы

(Е0 = 5.27ГПа, vo=0.3)

 

и трансверсально-изотроп­

 

ных

углеродных

волокон

 

(Е, = 8

ГПа, Е3 = 410.6 ГПа,

 

Vn = 0.568,

V13 = 0.273, //з=10.2

 

ГПа). Кривая 1 на рис.5.22

 

соответствует

расчетной

за-

висимости^параметра Я,+2т, от объемной концентрации волокон (светлые точки о _ экспериментальные значения),

2-m .(* ),3 -//.(x ),4 -w .(.).

§5.13. Термоупругая деформация композитов со сферическими и цилиндрическими слоистыми включениями

Рассмотрим термическую деформацию композитного ма­ териала в постоянном температурном поле. Деформацию сре-

273

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

родной среде с тензорами модулей упругости С° и коэффици­ ентов линейного расширения а° при действии постоянного температурного поля Т и поля эффективных внешних напря­

жений а , которое наводится окружающими неоднородностя­ ми. В рамках основной гипотезы метода, температурные на­

пряжения в области Vi, занятой / -м включением, определя­ ются из решения уравнения (2.1.17):

о(х) - JS{x - х’)В) c’)cs(x')dx’ = JS(x -

x’)a){x’)ctx' + <т,

У,

У,

 

B]=Bt-B °, а} = а , - а \

Т= 1,

(5.13.1)

где В], at - возмущение тензора упругой податливости и коэф­

фициентов линейного расширения в области Vi. Решение уравнения (5.13.1) можно представить в виде

о (х )= crr (x) + c f(x ),

(5.13.2)

где функции 0х(х) и с /(х ) удовлетворяют уравнению (5.13.1) с правой частью IS {x - х')а](х')с1х' и , соответственно.

Задача о температурных напряжениях а7 (х) в среде с изоли­ рованным сферически слоистым включением решена в §2.9. Выражение для тензора с/(х) в сферической системе коор­ динат (г, и), начало которой находится в центре включения, имеет вид

cF jr,n ) = 3к,(}?, - а))80ф+ 2 M 'r - 3k , -

,

ai_l < r< a i, i = l,2,...,N + l, = 0, aN+l= oo.

(5.13.3)

274

Здесь предполагается, что включение состоит из N слоев с

внешними радиусами аг Модули объемного сжатия к(, сдвига

fjf и коэффициенты линейного расширения а{ внутри слоя - постоянные, материал слоев - изотропный. Алгоритм вычис­

ления постоянных

и F/j указан в § 2.9.

Тензор < /(х ) является решением задачи для сферически слоистого включения в однородном внешнем поле напряже­

ния а . Решение этой задачи на основании результатов гла­ вы II можно представить в виде

< /(х ) = C{x)\j + А{х)]В°а ,

(5.13.4)

где выражение для тензора А(х) определено в (2.8.20), (2.9.2).

Если У(х)- характеристическая функция области, занятой включениями, то напряжения и деформации в композите в температурном поле Т= 1 удовлетворяют следующим соотно­ шениям:

о{х) = | 5 '(х -х ')[.б 1(х')<з(х0 + « 1(х')]г(х')<&:,> (5.13.5)

f ( x ) = е + | к (д г ')С °[5 1(х:')а{дс') + « 1( х ,) ] ^ ( х , )й^ '-

(5.13.6)

Выражение для эффективного поля напряжений о следует из (5.13.2), (5.13.4) и (5.13.5) после осреднения выражения для локального внешнего поля, в котором находится типичное включение

CT=j5(x-x')([Ar(^') + Af(x')cr*F(x;x')|x)<*:', (5.13.7)

Лг (х) = В1(х )о7 (х) + а' (х ), Л*(х) = Вх(х)С (х)[/ + Л (х)]/Г .

Предполагая статистическую независимость свойств вклю­ чений от их расположения в пространстве, среднее под зна­ ком интеграла в (5.13.7) представим в виде

( [ А т ( х ') + Л*(х')<7 ]F ( X ; X ')|X ) = ( y f + Л ' а М х - х ' ) ,

(5.13.8)

где V - о&ьем типичного включения, п - объемная концент­ рация включений.

Подставив (5.13.8) в (5.13.7) и вычислив интеграл с учетом регуляризации (5.2.11), найдем

а =-£>°(А г + A V ) , D° = j S{x\l-4f{x)]fix. (5.13.9)

Здесь и далее предполагается, что функция 'Р(х) сфериче­

ски симметрична ('Р(х) = хР(|дс|)), при этом тензор опре­ деляется выражением (5.9.17). Разрешив уравнение (5.13.9)

относительно <т*, найдем

a = - [ / + Z)°Af ] '1JD°Ar .

(5.13.10)

Осредняя теперь выражение для тензора деформаций (5.13.6), получим

(е) = а + 7 ? +Л 5сг.

(5.13.11)

Поскольку это выражение представляет собой среднюю деформацию композита при изменении его температуры на один градус, то < а > совпадает с тензором коэффициентов

линейного расширения композита а ,

который с учетом

(5.13.10)

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

 

а

= a + (l + D°AS) Лг,

(5.13.12)

Л'=(Л,)Е2+(Л2)(Е '4 Е г),

ЭКо =1

276

Отсюда следует, что а - изотропный тензор

= а,^<Ф> а

скалярный коэффициент линейного расширения имеет вид

С помощью формул (5.13.2), (5.13.4) и (5.13.10) можно вы­

числить напряжения сг^(г,«) в окрестности произвольного

включения

(5.13.14)

где слагаемыми, пропорциональными £ , учитывается взаимо­ действие между включениями

(5.13.15)

Заметим, что полученное выражение для cr^(r,w) следует

домножить на истинное значение температуры композита Т. В случае однородных включений выражение (5.13.13) для

а, принимает вид (N = 1, а, = а,кх=к)

Результаты расчетов напряжений в соответствии с соотноше­ нием (5.13.14) в некоторых случаях приведены в [57].

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

277

ной вдоль оси армирования. При этом тензоры коэффициен­

тов линейного расширения среды а и / -го слоя aj опреде­ ляются выражениями

= cc^ej^m) + ccommanip, а1(ф= а^в^т ) + ссштатр,

(5.13.17) где 0^(rri)=Sap-mantp\ а°в, aw - коэффициенты линейного

расширения в плоскости, ортогональной оси армирования,

а°т, ocim- те же величины в направлении оси волокна. Тензор­ ный коэффициент линейного расширения композита имеет

структуру, аналогичную (5.13.12). Входящий туда тензор определяется соотношением

D° = С°А°С° - С° ,

(5.13.18)

где тензор А° имеет вид (5.12.4), а тензоры Л* и Л? представ­ ляются в форме

а« = ( а1>р 2+(а2)[р 14 р 2]+(Аз)[р 3+р4]+(а5)р 5+(л6)р 6,

 

F=F(m ),

(5.13.19)

л , = - /> £ ( * , - * . ) ( i + j ; k . A , = - p | ; « 4 X i + i ! k .

i=l

 

1=1

Л2 = -2 р ^ (щ

- » J ( l + 2 i " +2J?)fi +3(г; +Ч )(а,г + а ,1 ,Ы ,

1=1

 

J

Л5=-2р]Г(м - я )(2+%)£ , Л6= - /> £ [(«, -и. )+2(/, - / )l* ]<£ ,

1=1

АТ ф = Ква/, +

/ < _

7 =И =

1

 

\

1=1

- К = - р £ (te i) ’ Am= - / ?Z ((’Л ) >

1=1 1=1

1

 

/, 1 ,

1

-/

и,-"1 1

 

Jo и .,

2Го-11

^__

а,п,

Z=

J

^5 +

ч—

----

 

1

 

 

 

 

 

 

 

 

 

 

278

В случае однородных волокон (N=1) и изотропных матри­ це и включениях тензор а в базисе ортогональных осей, од­ на из которых (ось х3) направлена вдоль волокон, имеет от­

личными от нуля только диагональные компоненты а*, = а12

и «зз. Выражения для этих компонент имеют вид

«п = а» +р{а~ а .)4 , «зз = «о + /> (« - «о3 , (5.13.20)

4 = ^ r [ EAf4(1-/>Xl+ У.Хv~ К,)], ^2=^-{[/>E+(l+ v)(l-p)E0]x

£

 

/ 1

1

— +-

- ( l - p ) v . E

- 4 ( l - p ) v .( l+ v ) ( у - v j } ,

 

 

Vк

д

А = р + Ь - р ) +

О Т

 

К

к

"* к. '

где Е,Е„; v, vo - модули Юнга и коэффициенты Пуассона

включений и матрицы, коэффициенты к.,1. и п, определены в (5.12.9).

Характер погрешности полученных формул для а обсуж­ дался в [76, 86]. Отмечалось, что для материалов типа стекло­

пластиков относительная ошибка в вычислении а по форму­ лам (5.13.12), (5.13.19) не превышает 10-15% вплоть до кон­ центрации волокон, близкой к плотной упаковке.

279

§ 5.14. Приближение точечных дефектов в механике матричных композитных материалов

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

Рассмотрим изолированную неоднородность с модулями упругости С(х) в однородной среде С° и пусть V- конечная область, занятая неоднородностью. Поля напряжений и де­ формаций в среде с неоднородностью можно представить в виде

о(х) = ст(х) + ^8{х-х')т{х')с1х\ /я(х) = Bx(x)o{x)V(x),

s(x) = е°(х) + ^ ¥ { х - х 9)С0гп(х9)сЬс\ Вх(х) = В{х)-В \

(5.14.1)

где т(х) - плотность дислокационных моментов, индуциро­

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

фекту соответствует замене плотности т(х) первым членом ее разложения в ряд по мультиполям (§1.4)

т(х) = MS{х - £)+..., М - ^m{x)dx, (5.14.2)

где £ - точка приведения, которую удобно выбрать в центре тяжести включения.

В случае постоянного в области V внешнего поля плот­ ность т(х) и тензор М определяются выражениями

т{х) = М{х)а\ М = М°а°, М° = J M(x)dx, (5.14.3)

280

где функция М ( Х ) определяется из решения задачи об изоли­

рованной неоднородности в постоянном внешнем поле сг. В частности, в случае эллипсоидальной неоднородности тензор

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

М° = -vB°C'[l + А{а)С']ХВ°,

(5.14.4)

где v - объем включения, тензор А(а) зависит от формы эл­ липсоида и определяется соотношением (2.4.2). Если неодно­ родность представляет собой эллиптическую трещину, кото­

рая раскрывается в поле сг°, то

М °аРХ„ = J (Г )'т П„) > Р -14-5>

где п - нормаль к плоскости трещины, а},а2 ее полуоси, тен­

зор Т° определен в (2.7.11).

Пусть теперь однородная среда содержит множество точе­

чных дефектов. Обозначим через о\ локальное поле, в кото­ ром находится 1-й дефект. Тогда поля напряжений и деформа­ ций в среде с точечными дефектами на основании (5.14.1),

(5.14.2) представляются в форме

 

о{х) = <т (xr) + £ J S{x - х ')м °а*ё(х' -

£ )dx', (5.14.6)

s(x) = £ (х) + £ J К(х - x')C°M°o)S(x' - £ )dx'. (5.14.7)

i

 

 

Здесь М°сг* - коэффициент при S(x -

) при разложении

плотности mt (х) для i-ro

включения в ряд по мультиполям

(5.14.2) .

 

 

Поле сг*, в котором находится к-й дефект, имеет вид

(& ) + Z j s ( &

-х ')М °а Д х '-^ )(к ', (5.14.8)

i* k

 

 

где ^к{к-1,2,...) - точки, в которых имеются точечные дефекты. Рассмотрим простейший пример. Пусть на плоскости име­ ются две прямолинейные трещины, длина каждой из которых

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