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

Механика композиционных материалов

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

Для определения локальных функций М(0)(£) нужно решить за­ дачу Ж б(0) :

A ijm n № p m & q n ^ mnpq) +

jm nk ^ m n p q \k + D f Jmnk i ^ \n n p q \k l = 0 , (6 .2 0 )

где величины A, В, D находятся по формулам (3.16), (3.17), (3.6). Если локальные функции M ^ mnpq найдены, то находятся ве­

личины К простым осреднением

{ D ijm n k l W m n p q S k r $ ls + &Рт ^ q n ^ rk fys)) = K - ljp q rs »

(6.21)

а по ним — и эффективный тензор упругих податливостей Н:

и итл =

DUmMbH +

( l , 7 - 4 - в / ; ) .p« - " V 3g.

+

 

 

 

 

 

V

У

2 — 6

 

 

Н

(6// Rmn

Rln &inRjm

6/mRin

6/л/?jm).

 

(6.22)

Пусть

теперь

требуется

решить

задачу

В

теории

упругости

в напряжениях

 

 

 

 

 

 

 

 

 

 

OV/./ +

X t = 0, ть7 (£) = 0,

 

 

(6.23)

 

 

 

 

* ,/* /b = S °c .

 

 

 

(6.24)

Если объемные

силы

обладают потенциалом

(3.30),

то

можно

ввести тензор функций напряжений <р (3.31)

и для него

сформу­

лировать

задачу

(3.32) — (3.35). Решение этой

задачи

по

теории

нулевого приближения имеет вид

 

 

 

 

 

 

ф£/ =

Уи М + [a [Mj/iimn (I) Ь и т п

(X) + X ,/],

 

(6.25)

где тензор функций напряжений ф(я) определяется по теории эффективного модуля:

€гри

jmnHpqti tybnjmut = ^rs*

(6.26)

 

jmn

lmnj Is — <5^.

(6.27)

Э тензор % удовлетворяет условиям

[ H p q ii %kn,lmut — 0

(6.27')

e ik l Gjmn t y k t ijm + M knpqrs typq.rslm) K j |s — 0»

где правые части в (6.26) и (6.27) — заданные величины (3.33), (3.35).

Локальные функции М(2)(|) представимы в виде

Mijhmn (£) — ^pln^qmk ^Hpq (£)>

(6.28)

причем М(£) должны находиться из решения уравнений

e ikp G jrl [J k lu v (£ )

+

€ k'p 's e l'm r' X

 

 

 

 

 

X J k lk 'V

M .p rr'uv\m s]\pr 0*

 

 

 

(6.29)

Локальные функции M также должны удовлетворять условиям

(6.7) и (6.8). Как только

уравнения (6.29)

решены, можно

найти

эффективный тензор податливостей

Н:

 

 

 

 

 

 

Н е т =

iiu v (& uk 6/у +

£ kps e im r M

pruv\ms) )•

 

 

(6 .3 0 )

Тензор напряжений

т по

теории

эффективного

модуля

опреде­

ляется по формуле

 

 

 

 

 

 

 

 

 

 

 

 

x i i =

eimn ty k n jm

+ ф

6*7»

 

 

 

(6 -3 1 )

а тензор деформаций еэ — по формуле

(6.18).

 

 

 

 

По теории нулевого приближения тензор напряжений о выра­

жается по формуле

(3.31),

где

тензор <р

находится

из

(6.25),

а тензор деформаций е — по формуле

 

 

 

 

 

 

 

eS?' = Д Ц | )т4;,

 

 

 

 

(6.32)

где тензор J(0> называется

тензором

упругих податливостей

нуле­

вого приближения:

 

 

 

 

 

 

 

 

 

 

 

(£ ) = J

&

v

l

+

^ kps

е Im r

M p r uv\ms ( ^ )l’

 

(6 .3 3 )

При ЭТОМ

 

 

 

 

 

 

 

 

 

 

 

 

 

Н ,ш = (4 % ).

 

 

 

 

(6.34)

Заметим, что тензоры О 0)

(6.12)

и J(0> (6.34) взаимно-обратными

не являются.

 

 

 

 

 

 

 

 

 

 

 

Упражнение 6.1. Доказать, что тензоры

С(°) и J<°) связаны меж­

ду собой соотношениями

 

 

 

 

 

 

 

 

 

 

C\jh (£)=• Cijmn (£) Jmnpq (£)

 

 

 

 

(6.35)

 

(? ) =

J i i« n (? ) c iiS U

(? ) V

 

 

 

<6 -3 6 )

Если при вычислении напряжений по формуле

(6.25)

положить

а = 0, то микронапряжения учитываться не будут.

для

их

решения

Заметим, что во

всех

трех задачах

А,

Б, В

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

Поэтому, найдя для данного композита указанные выше ло­ кальные функции, можно найти для него микроперемещения и

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

Точно так же строится теория нулевого приближения для зада­ чи теплопроводности (4.3), (4.4). А именно:

Т =

Ъ(х) + a

[Pt ) Ь (х) + 0 (*)],

(6.37)

где локальная функция P(Q

удовлетворяет

уравнениям

(4.12):

 

{kmn Pj\т“Ь ^/л)|п = 0.

 

(6.38)

После решения уравнений (6.38) находится

эффективный тензор

теплопроводности Л

(4.11)

 

 

 

 

л , / = < t f* />/.*+ $ )

 

(6.39)

и по теории эффективного модуля определяется Ф(я):

(р Ср)

= Л,7 *,(/ +

Р<?,

(6.40)

Ai/d ./ «,|s = t| (*| n -7 ’c).

(6.41)

Если известны локальные функции N(%)

как решение уравнений

 

 

 

(6.42)

то вычисляется эффективный

тензор р*,

связанный

стензором

теплового расширения (4.32):

 

 

 

( £ /= < Р //- с „« М | 7 > .

(6.43)

Теория нулевого приближения для непериодических регуляр­ ных структур отличается от теории периодических структур тем, что локальные функции (и, может быть, тензор модулей упруго­ сти) зависят от медленных координат. Пусть, например, требуется решить задачу (5.11) — (5.14). Согласно теории нулевого прибли­ жения решение этой задачи ищется в виде

И т= Vm(x) + Ct[W{o)m (£, x )v ,(x ) + N(t)m (|, X) V,.k (*)] + <lWm, (6 .44)

где индекс q—1 у локальных функций, введенных в (5.16), опу­ щен.

Вектор v определяется из решения задачи по теории эффектив­

ного модуля:

 

 

 

 

hilvt +

hiklV{i,k) + hik*lk^vitk,kt + X1=

0;

(6.45)

vt Jv, =

w?,

+ hi,lkvi,k) щ Is, =

So,

(6.46)

а вектор w удовлетворяет однородным уравнениям (6.45)

и

w,\z,= - [ N ml

C,x)vlx) + Л #ц (~ .3»л * W ] k -

(6.47)

Величины h можно определить по формулам (5.28), (5.29):

ft" (х)= ( — rL,;C"m + 2 Г/г С'>'тп х

X (AfJo)ml[r> — r L ) + г L

С‘! Г -

(C

"

(6. 48)

ft'*' (x) = <C'*-""1V{0)m,„ -

r L C 'ta"

+

2ТУгСПшпМ‘кт п +

 

 

+ 2Г<•£'>"* + (C "mnN[kW)mtn).i),

(6.49)

 

ft"'* (x) = (C""»wtf )m|„ +

C "'*),

(6.50)

 

ft"' (x) =

(C """ yVjo,m,n -

r l C " - ) .

(6.51)

а локальные

функции

N(oj и

N(i)

находятся из решения

задач

Ж а (— 1. —2)

и Ж а (— 1. — 1) ((5.25) и

(5.26)):

 

 

 

-

r L

C "m!’)l, = 0,

(6.52)

 

(C""»W{?)ш|„ +

С"'*),, =

0.

(6.53)

Если рассматривается квазипериодическая среда, то после определения эффективного тензора модулей упругости hijlk осталь­ ные величины h можно найти по формулам теории эффективного модуля:

ft" =

-

ft"ra0 rL,/ -

2Г<‘ ft'1™" rL.

 

(6.54)

ft'« =

—ft'*”"TL +

2Г/^ A/w \

 

(6.55)

 

 

f t " '= - f t " m" rL.

 

 

(6.56)

Напряжения по теории

нулевого

приближения

подсчитываются

по формуле

 

 

 

 

 

 

 

 

0(0) = С(0) (£, х) vk(х) Н- C(7of (£, х) у*,; (я),

(6.57)

где

 

 

 

 

 

 

 

 

СЖ =

С"'"”

- rL C "™ ,

 

(6.58)

С{4? = С

« « * ^ ф

+ С "*'.

 

(6.59)

Напряжения по теории эффективного модуля имеют вид

alJ =

Hik(х) vk(х) +

hiikl (л) Vk,i (?).

 

(6.60)

Если при вычислении перемещений по формуле (6.44) поло­

жить а= 0 , то перемещения

будут

удовлетворять

кинематическим

граничным условиям (6.46)

и совпадать

с перемещениями, вычис­

ленными по теории эффективного модуля

(т. е. микроперемещения

учитываться не будут).

 

 

 

 

 

 

 

Упражнение 6.2. Показать, что если при рассмотрении регуляр­

ной структуры можно

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

систему

координат (Г*тп= 0 ),

то теория нулевого приближения

для регу­

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

В заключение заметим, что напряжения по теории нулевого приближения а»/0) можно выразить через напряжения по теории эффективного модуля ацэ= хц по формуле

йц* ( S ) =

Ci}kl ( £ ) Jklmn ( £ ) т т л = Ан'тпШ Тт л >

( 6 . 6 1 )

ГДе Aijmn (£) — так

называемый тензор концентрации

напряже­

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

Аналогично, деформации по теории нулевого приближения е*/<°> выражаются через деформации по теории эффективного модуля еif>==V(iyj) по формуле

®i/* (Е)г= JIjkl (£) Cklmn(£) и(т,п)=Вцтп(|) U(m,n),

(6.62)

где Вцтп(1) — так называемый тензор концентрации

деформа­

ций. Тензоры А (Е) и В(|) взаимно-обратными не являются.

Упражнение 6.3. Показать, что тензоры А(£) и В(£) связаны между собой зависимостями

Bljmn (£) = Jiikl (£) ^ klpq (£) hpqmn»

4ц,пЛ) = СиА ) Вы„ (Ё)Ямтя. €

(6.63)

Можно получить оценку точности нулевого приближения в за­ висимости от параметра а:

Нет—а(0)Н*.а<У а-С, C= const.

(6.64)

§ 7. Неидеальный контакт

Рассмотрим поверхность Г0, разделяющую один компонент композита от другого. Согласно (1.2.16) ^ на этой границе осуще­ ствляется идеальный контакт, если на ней выполнены условия

= w<-2> (oj/n = o|/,)) « J = 0,

(7.1)

где индексом (1) помечены величины, относящиеся к компонен­

ту (1), а индексом (2) — к компоненту (2), п — единичный век­ тор нормали, внещней для компонента (1) и внутренней для ком­ понента (2). Другими словами, величины на поверхности Г0, по­ меченные индексом (1), получены предельным переходом со сто­ роны компонента (1), а индексом (2) — со стороны компонен­ та (2). Условия (7.1) можно записать сокращенно:

[К ]]г .= 0. [ К л/]]г .= 0.

(7.2)

В § 2 гл. 1 было показано, что обобщенное решение задачи МДТТ автоматически удовлетворяет условиям (7.1) на каждой поверхности, разделяющей компоненты композита.

В § 2 настоящей главы введены понятия псевдоперемещений и псевдонапряжений, полезные при интерпретации задачи Ж(<7).

Предположим, например, нам нужно решить уравнения

(6.6), ко­

торые, используя обозначения

(2.20), (2.21), можно

записать

в виде

 

 

2*/<й*)|/ (I) =

Саш © »

(7-3)

где компоненты 2*/*/ тензора псевдонапряжений первого уровня,

или просто тензора псевдонапряжений,

связаны

с компонентами

и тк1 вектора

псевдоперемещений по закону Гука

(2.21):

 

= CijmnUrrHkl),ny

(7.4)

где

 

 

 

 

 

Um(kl)=

Nmkl.

 

(7.5)

Решением

задачи (6.6) — (6.8)

для

определения локальных

функций (вектора псевдоперемещений) является непрерывный вектор псевдоперемещений Ощ), для которого на границе Го в си­

лу (7.1)

выполняются условия

 

 

 

 

 

£/$) = Й|1!) (ЛГЙ =

ЛГ$|).

 

(7.6)

а в силу (7.2)

— условия

 

 

 

 

 

[(Sl/V«) + Cj/ii) 2щы) +

Cqli)\ tif =

0.

(7.7)

Условия

(7.6)

и (7.7) также можно записать сокращенно:

 

 

[[£/(/г/)]]г0 = 0* [[(2г/(*о + c ifkl) tij]]Го =

0.

(7.8)

Псевдонапряжения 2(И) при переходе через поверхность Го терпят разрыв, ибо «объемные силы» Сцы,1 в уравнениях равновесия (7.3) являются производными от разрывных на поверхности Го функций координат.

Условия (7.5) и (7.6) и являются условиями идеального кон­ такта на поверхности Го, так как они автоматически выполняют­ ся при решении задачи (6.6) — (6.8).

Однако при деформировании реального композита возможны нарушения идеального контакта: отслаивание, проскальзывание одного компонента относительно другого и т. п.

Например, в случае, если на границе Го, представляющей со­ бой плоскость х3 = const, происходит проскальзывание одного ком­ понента композита относительно другого без отслаивания, условия контакта будут иметь вид

4 " = 4 2), ой> = ОЙ». ОЙ» = а!з2) = ОЙ» = ^ > = 0.

(7.9)

Можно записать и более сложные условия, например условия, когда на границе Го учитывается трение или частичное отслаива­ ние одного компонента от другого. Однако такие граничные усло­ вия будут нелинейными, а наша основная задача сейчас состоит в том, чтобы, сохранив методику осреднения, описанную в преды­ дущих главах, учесть условия некоторой неидеальности контакта на границе Г0.

Вспомним о граничных условиях контактного типа (1.2.7). Они являются линейными, но описывают достаточно широкий класс

граничных условий. С их помощью условия контакта на

грани­

це Го можно записать в следующем виде:

 

 

 

 

 

«У « Г

+ 4Р <4n = eg> 4 ?’ 4 2) + 4 ? и?\

р =

1,.

_.. .6,

(7.Ю)

где [a*1)], [&W],

[а(2)], [6(2)] — некоторые матрицы 6X3,

/г*1)

век­

тор внешней

нормали на Го для компонента

(1),

а

п

— для

компонента (2).

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

позита

(6.1),

(6.2)

при условиях контакта на каждой поверхно­

сти Г0,

разделяющей один компонент от другого в виде (7.7).

Будем

решать

эту

задачу, например, методом нулевого прибли­

жения, т. е. искать решение в виде (6.3), где v — решение задачи (6.4), (6.5) по теории эффективного модуля. Однако сам эффек­ тивный тензор модулей упругости h для этого случая будет, вооб­ ще говоря, отличаться от соответствующего тензора для задачи с идеальным контактом. Ведь для определения эффективных мо­ дулей нужно найти решение задачи (6.6) — (6.8) для определения локальных функций. В случае идеального контакта на поверхно­ сти Го выполняются условия (7.8), а в случае неидеального кон­ такта их нужно в соответствии с (7.10) заменить условиями

[[Яр; (%ijw) + Сцы) tij + frpi£/;(fto]]r0 = 0 •

(7.11)

Упражнение 7.1. Подставляя в (7.10) выражения для переме­ щений (6.3)

“ i = vt (х) + a [Nijk(£) Vj'k + Wt]

(7.12)

и напряжений (6.11)

tfi/= CUml(l)[Vm,l(x) + aNmpq\i (fj Vp,q (x) + aWm.l],

(7.13)

по теории нулевого приближения показать, что они эквивалент­

ны условиям (7.11), так как для вектора v выполняется условие

[Й ]г .= 0.

(7.14)

Упражнение 7.2. Показать, что если для упругого стержня (см. § 1), являющегося двухкомпонентным композитом, условия идеального контакта на границе раздела компонентов £= у

[[Щ] = 0, [[£ + £ Л /]] = 0

(7.15)

заменить на условия неидеального контакта вида

[[N]] = О,

[[Е+ E N -+ аЩ] = О,

(7.16)

то эффективный модуль будет иметь вид

 

 

2

{Е)

 

h =

У(1—У)

+(а(1) - я (2)) Е{1) Е {2)

(7.17)

} +(‘<,,-‘<2,)/£(,,£,2)

^

( т

 

причем если в (7.16)

 

aO>=a<2),

(7.18)

 

 

то эффективный модуль совпадает с эффективным модулем в слу­ чае идеального контакта (7.15):

A= I / ( T )' €

(7Л9)

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

§ 8. Плоская задача теории упругости

Если в теле возможно существование плоско-деформированно­ го, плоско-напряженного или обобщенного плоско-напряженного состояний, то естественно сформулировать двумерную задачу МДТТ. Сделаем это на примере задачи теории упругости.

Для задачи теории упругости в перемещениях плоская задача ничем по существу не отличается от пространственной, разобран­

ной в § 2: нужно только

формально все малые латинские буквы

в индексах § 2 заменить

на большие (см. приложение I). Сле­

дует иметь в виду, что тензоры модулей упругости и упругих по­ датливостей в плоской задаче имеют в общем случае 6 незави­ симых констант (против 21 в пространственной задаче). Эти кон­ станты будем называть приведенными.

Обозначим тензор, обратный

к

тензору CUKL, через

7UKL,

а тензор, обратный к JUKL, — через CuKL\

 

CJJMNJMNKL = 2 IJMNJMNKL =

{&IK§JL + 6 /zA//<).

(8.1)

Так как для плоской деформации

 

 

 

ei3= 0 ( £ = 1 ,2 ,3 ),

(8.2)

то связь между тензорами напряжений и деформаций имеет вид

ff/j = CUKL^KL, &и = JIJKLPKL,

(8.3)

при этом

<*з/ = CZJKLJKLMNGMN ( / = 1 ,2 , 3 ).

(8.4)

Для плоского напряженного состояния

 

<7/3=0 (i= 1 ,2 ,3 ),

(8.5)

и связь между напряжениями и деформациями будет

 

<Уи = C u K L ^ K L ,

Ь и =

jIJKL&KL,

(8.6)

причем

 

 

 

ез/ = J V K L C K L M N

&M N

(i= 1,2, 3).

(8.7)

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

Плоская задача теории упругости в напряжениях заключается в решении двух уравнений равновесия (1.3.32)

o’/y .j +

X /

= 0 ;

X / = — ф (/

(8 .8 )

и уравнения совместности

(1.3.31)

 

 

е ш

GJ N

[JI J K L

(*) 0/а], M N

(8.9)

относительно трех независимых компонент плоского тензора на­

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

на конту­

ре Г, ограничивающем замкнутую односвязную

область

2

(1.3.33):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gutij |г = 5?.

 

 

 

(8.10)

Здесь и далее J*==J

при

плоском напряженном

состоянии

и

J *= J

при плоской

деформации. Точно

так

же С* =

С

при пло­

ской

деформации и

С* = С при плоском

напряженном

состоянии.

Если ввести

функцию напряжений Эри <р (1.3.34)

 

 

 

 

 

 

<7/а= ёкрб1рф,рр + Ф6ях>

 

 

(8.11)

то уравнения

равновесия

(8.8)

удовлетворятся тождественно,

а уравнение совместности

(8.9) даст

 

 

 

 

 

где

6lM 6JN6KP G L R [JI J K L M

ф,рд],млг =

Л (*)»

 

(8.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т ]= G I M CJ N ^ K L

{JI J K L $)),M N

 

(8.13)

Граничные условия

(8.10) приобретут вид

 

 

 

 

 

 

(^IK^JL^.KL^J +

Фл/) |г = 5/.

 

(8.14)

Представляя функцию напряжений <р в виде

 

 

 

 

 

 

Ф = £

и 'Ч '.’..»,(I)'I’.V -* ,W .

 

(8-15)

Я= 0

где

 

М(0)= 1, M W sO для ^<0,

(8.16)

и проделывая обычные процедуры техники осреднения, получим

две рекуррентные

последовательности задач Д в(£), &=0, 1, 2,

и Жв(<7), <7= 0, 1,

которые были разобраны в § 3 при решении

задачи В.

 

Первая из этих последовательностей возникает в результате

представления

 

 

 

 

 

ф(*) =

£ а 'х И Й ,

 

(8.17)

 

 

Р = 0

 

 

причем каждая

задача Д в(р),

р = 0, 1, 2,

заключается

в оты­

скании функции

путем решения уравнения

 

 

cm GJN бцр eLp Н

 

 

(8.18)

и удовлетворяет граничным условиям

 

 

 

eIKeJL X w -4 | r= s “ (4

 

(8.19)

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

Входные

данные задачи Дв(р) определяются из решения

задач

Д в(г), г=

1, 2, .... р— 1:

 

 

Т){р) = - £ H Z ..K r+ A 7 ' \ « ’ р >°>

<8-20)

S / ^ = S ? — Фп/|г,

где

T u lg l.K ^ ^ e/pejRM^^K^PR + («/к?+а х

X eJP + ejR ^P ip) M ^ + "V l |P + е‘ кя^ к ^ гМ.к1. к1=

=eipeJRr PK*)l

.KiitjlR + eipeJK(/tt

(8.22)

Мк,+.'к,+11|Р+ МкУ.к^Ьрк^. (8.23)

Эффективный тензор упругих податливостей Н и величины определяются при решении второй рекуррентной последователь­ ности задач Жв(<7).