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

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

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

71

относительно b\ придем к следующим формулам для сим­

метричной Ь\ф и антисимметричной Ь\ар\ частей этого тен­ зора:

fyap) —h-apXfPxn> \ар\ ~ Ь[ар\>

 

 

(2.5.20)

Л '= — Е2 +

5С"

Е1

1 Е2Л , С

= С с + C j .

2С + З С

3

у

 

Отсюда, в частности,

следует,

что

поле

КрХр = £^(х) не

возмущается сферической неоднородностью, если Ь"ар - анти­

симметричный тензор.

Если среда изотропна, а источники поля отсутствуют, то

внешнее поле £а(х) удовлетворяет уравнению V a£„(x)=0. В

случае линейного внешнего поля отсюда следует, что tr b°=

=Ь°аа=0. С учетом этого равенства соотношение для симмет­

ричной части тензора Ьхупрощается и принимает вид,

5С0

(2.5.21)

5С о + 2С ,

 

 

У .Квадратичное внешнее поле. Рассмотрим случай квад­

ратичного внешнего поля

 

£°Ах) = КрХХрХх,

(2.5.22)

где b°apX - постоянный тензор, симметричный по двум пос­

ледним индексам. Как следует из формул (2.3.22), (2.3.24) и (2.3.25) общей теории, поле внутри включения будет квадра­ тичным полиномом, у которого отсутствуют линейные по X слагаемые. Это позволяет искать решение уравнения (2.5.13) в виде

saix) = b°a+b2pXxpxx,

(2.5.23)

72

где b°a и blpX- неизвестные постоянные тензоры. Подстановка

этого выражения в (2.5.13.) дает

К +

+J К^ (х - x ')C i {Ъ\ + blMvx'Mx'v)ctx’ = b°afsxxf3xx .

 

v

 

(2.5.24)

 

 

 

Продифференцировав это уравнение дважды по коорди­

натам X и положив затем х = 0, получим

 

Ь<фХ+ ActftPXvpCftsbsvp —Ьф, •

 

(2.5.25)

При этом учтено, что

 

 

VjK^Or-*')<&' =о

при x eV ,

(2.5.26)

так как

v

 

 

результат интегрирования - постоянный тензор, а

тензор А2 определяется выражением (см. Приложение П2.2)

, (2.5.27)

^1 Если среда изотропна, а включение имеет форму шара, то

А2- изотропный шестивалентный тензор. Для его представле­ ния удобно ввести следующий базис из шестивалентных тен­ зоров, составленных из символов Кронекера

QaXpfipi ~ ^affixp^pz »

QaXpfipz ~ ^afi^X{p^z)p у

QaXpfipz ~ ^^а{р

рт> Qita*, = W . < A A .. Р-5-28)

QaXpfipz ~~ ^^а{Х^р)(р^т)Р> QaXpppv= 4др(хдр)(р8г)а ■

В этом базисе имеем

( Q 4 Q 3) + T AT (Q2+Q4 +Q5+Q6)- (2.5.29) 14С0

 

73

Разрешая теперь уравнение (2.5.25) относительно

с по­

мощью приведенных в Приложении П1.2 формул обращения шестивалентных тензоров, заданных в базисе (2.5.28), получим

 

L2

_

д2

то

 

(2.5.30)

 

^CL^jj. ~~

'‘а.Рцкрт®крт

 

Л2 =

О Д

^(Q' + Q3) + 1+ #

24 [l+^ f (Q4 Q5)4 Q6

 

 

 

'

зс, у 1

,

.-1

 

 

 

1+ £

 

 

 

dx=

1 + ^ М ,

^2 =

 

 

 

 

7С„

 

с..

Вернувшись к уравнению (2.5.24), положим в нем х=0. В

результате найдем

 

 

 

 

ba+ AaxCXpbp - JappzCpvbvpt,

(2.5.31)

 

J appz = - J

 

 

(2.5.32)

 

 

 

к

 

 

 

где тензор /1°^ определен в (2.5.10). Заметим, что интеграл

(2.5.32) имеет слабую особенность и может быть вычислен непосредственно. Для изотропной среды и сферического включения радиуса а имеем

^=^г(е2-зе'),

v-1

л-=( i + A

 

3с.

(2.5.33) Подставляя сюда выражение (2.5.30), получаем окончатель­

но

К = с? С,Л°

ж

1+

Cj_

 

(2.5.34)

15С JaXX

 

зс

ЛЛа

 

74

Если источники поля отсутствуют, то b°XXflxfl= 0, откуда

следует, что b°UfJ = 0 . В этом случае формула (2.5.34) для век­

тора Ъ°а принимает вид

(2.5.35)

§ 2.6. Трещина в однородной упругой среде

Рассмотрим однородную среду, в которой имеется сплюс­ нутая полость, занимающая конечную односвязную область V с гладкой границей. Интегральные уравнения для среды с по­ лостью имеют вид (2.1.24) и (2.1.26):

£•(*)- jK (x -x ')C 0f +(x')d5c' = £°(*)>

(2.6.1)

(2.6.2)

Здесь под £+(х)= е(х) V(x) понимается предельное значе­

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

напряжений вне полости через поле £+(х).

Рассмотрим предельный переход от сплюснутой полости к трещине - разрезу по расположенной внутри Vгладкой повер­ хности Q . Выберем в точкех e Q локальную систему коорди­ нат с осью z , направленной вдоль нормали п(х) к поверхнос­ ти Q . Пусть И(х) - поперечный размер полости, zx(x,h),

z2(x,h) - координаты пересечения оси z с границей области

V, причем zx,z2—>0 при h —>0. Для фиксированной точки

х V ядро S(x) в уравнении (2.6.2) - гладкая ограниченная функция. Следовательно, главный член асимптотики поля на­ пряжений вдали от полости имеет вид:

75

о(х) = сг°(х) +

(2.6.3)

п

 

z2(x,h)

 

~e(x,h) = J e+(x + zh(x))dz, x e Q .

(2.6.4)

Z\{x,h)

 

Величина ~ё(х,И)может быть интерпретирована как коэф­ фициент при главном члене разложения функции е+(х) в ряд по мультиполям, сосредоточенным на поверхности Q (§ 1.4):

£+(х) = f(x,/;)Q (x)-h..

(2.6.5)

Здесь П(х) - дельта-функция, сосредоточенная на по­ верхности Q и определенная соотношением (1.1.22).

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

па(х)&ар(х) и перемещений иа(х) удовлетворяют условиям

па{х)(тар(х) = 0, К(х)] = 6а(х), х е П, (2.6.6)

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

ничения на внешнее поле <т (х) .

Поскольку поле перемещений разрывно на Q , то тензор деформаций в среде с трещиной есть сумма регулярной d и

сингулярной ё1составляющих:

£(х) = ег(х) + £*(х), esap(x) = n(a(x)bp)(x)П(х). (2.6.7)

Из

сравнения этого равенства с (2.6.5) следует,

что

при

h —>0

невыписанные слагаемые в (2.6.5) исчезают,

а

пре­

дельное значение ~ё(х) функции ~ё+{х,И) имеет вид

 

 

£а/3(дг) = Нш?;/,(х,/;) = /?(а(х )^ (х ), х е О . (2.6.8)

76

Отсюда и из (2.6.4) видно, что при И—>0 "нормальная" компонента тензора е+(х) должна стремиться к бесконечнос­

ти, тогда как функция s(x) конечна.

Таким образом, поле напряжений вне трещины при учете соотношений (2.6.3), (2.6.8) и симметрии тензора S(x) имеет вид

° 4 (x) = o fafi(x ) + jS & e (x - x ,)nx(x')bM(x')dn'. (2.6.9) Q

Отсюда и из первого условия (2.6.6) следует уравнение для векторного поля Ь(х) на Q:

(Tapbp){x) = J 7^(х,х')6Д х'>Ю ' = nfi(x)cfafi(x), (2.6.10)

Q

 

ТаР(х,х') = -nx(x)SXafiM(x - х')пм(х').

(2.6.11)

Заметим, что оператор Т в уравнении (2.6.10) может быть

записан в форме интегрального оператора с

ядром Т(х,х')

лишь условно, поскольку соответствующий интеграл фор­ мально расходится при х для сколь угодно гладкой

функции 6(х) (7’(х,х')~|х-х'| 3 при X—»х').

Нетрудно видеть, что представление (2.6.9) совпадает с полем напряжений в упругой среде, содержащей дислокацию Сомильяны (1.2.34), а уравнение (2.6.10) эквивалентно урав­ нению (1.2.36), к которому приводит дислокационная модель

трещины (§ 1.2).

Рассмотрим подробнее обобщенную функцию Т(х,х') - ядро интеграла (2.6.10). Пусть для простоты Q есть область на

плоскости х3 = 0 в трехмерном пространстве, а функция Ь(х),

х = х(х,,х2) продолжается нулем вне этой области. Тогда Т

есть оператор свертки с обобщенной функцией Т(х), которая порождается обобщенной функцией вида (1.2.9). Из (2.6.11) имеем

Tafiixi,x2) = *2,х3)1*=о • (2.6.12)

77

Отсюда следует, что преобразование Фурье Т*р{кх,к2)

функции Тсф(х1,х2) имеет вид

 

Т'а1}(к„к2) = — J S;afB(k^k2,k3)dk2,

(2.6.13)

где S*(к) определено в (1.2.11). Нетрудно убедиться, что этот интеграл сходится абсолютно и определяет четную однород­ ную функцию первой степени по к . Поэтому, если, напри­

мер, Ь(х) - функция из пространства S(R2) (затухающая на бесконечности быстрее любой отрицательной степени |лг|), то действие на нее оператора Т определяется соотношением

(ТЬ)(х) = -^\Г(к)Ь*(к)е-**с1к, к = к(кик2), (2.6.14) (2я) '

где интеграл является абсолютно сходящимся. Этой фор­ мулой дается регулярное представление оператора Т на

функциях из S(R2)в случае плоской области Q .

Если Q - произвольная поверхность Ляпунова, ограни­ ченная гладким контуром Г, то регулярное представление оператора в (2 6.10) на достаточно гладких функциях постро­ ено в Приложении П3.2 и определяется соотношением

( W ( x ) = |

( * ') - * , W * * 1' -

 

Q

 

 

-n v(x)j rotxpZvaXfJ(x -

x')dTfiM(x ),

(2.6.15)

Г

в котором для существования интегралов достаточно, чтобы функция 6(х) была непрерывно дифференцируема на Q и обращалась в ноль на контуре Г области Q .

Вслучае произвольной гладкой поверхности Q оператор

Тявляется обобщенным псевдодифференциальным операто­

ром с главным однородным символом Т*(к) вида (2.6.13) -

78

однородной функцией степени 1. Общая теория разрешимос­ ти уравнений типа (2.6.10) изложена в [148]. В частности, для бесконечно дифференцируемой правой части в (2.6.10) это уравнение имеет единственное решение, причем его асимп­ тотика вблизи гладкой границы поверхности Q (контура Г) имеет вид

Ь(х)=0(хо) ^ + О(гУ2),

(2.6.16)

где Г - расстояние от точки х e Q до хоеГ

по нормали к Г,

/?(хо) - бесконечно дифференцируемая вдоль Г функция.

Рассмотрим асимптотику поля напряжений о(х) вне тре­

щины в окрестности ее кромки

Г. Пусть у 1,у 2,у 3- локальная

декартова система координат в

точке хо е Г , причем ось у3

направлена по предельной к Q

нормали в точке хо, ось у2 -

по касательной к Г, тогда ось _у, лежит в касательной к Q

плоскости в точке хо. Учитывая асимптотику (2,6.16), имеем следующее выражение для асимптотики вектора скачка пере­

мещений на трещине в окрестности точки хо е Г

а д = Д ^ )л Я + 0ОГ2)-

 

(2.6.17)

Используя уравнение (2.6.8), запишем выражение для

тензора напряжений в точке у

с координатами ух= —г cos в,

у2= 0, у3 = —тsin 0, где г - расстояние от точки у

до начала

координат, в- полярный угол в плоскости (у,,_у3) (£

- у { /г)

о(у)=а(у)+-^ J S{c os#^, , £ , s

i n

) A/ ^’^ + 0 ( d ).

Л' ГОСг)

 

 

 

(2.6.18) Здесь учтено, что S(x)- четная однородная функция сте­

пени (-3). Интегральный множитель при г~ш в этом выра­ жении при г —>0 имеет конечное значение

J(d, хо) = lim -Jro(y) , r - » 0 .

(2.6.19)

 

 

 

79

Учитывая, что при г —» О поверхность Q(r)

в координатах

<£. переходит в полуплоскость (£3 = 0, £ ,> ()),

получим

сле­

дующее выражение для компонент тензора J(0,xo):

 

Jap(e>xо) = за1в*(О)пЛх.)0м(хо')>

(2.6.20)

 

* ( 0 =

J S ( c o s 0 + £ , & , s i n 0 ) ^ 2, (2.6.21)

 

0

- о о

 

 

где п(хо) - предельное значение нормали к Q в точке хо е Г .

Таким образом,

функция J(0,xo) представляется в

виде

двух сомножителей, первый из которых S(ff)n(xa) не зависит от формы поверхности Q и действующего в среде внешнего поля и определяется только локальной ориентацией края Q в точке хо е Г . Второй сомножитель (вектор /?(хо)) есть функ­ ционал всей поверхности Q и внешнего поля <т (х ).

Из (2.6.19) следует, что

 

Ыу) = -U J(0, хо) + 0 (а )

(2.6.22)

у/г

 

и поэтому J(6,xo) можно назвать тензорным коэффициентом интенсивности напряжений.

Функция J(0,xo) допускает наглядную интерпретацию, ес­ ли учесть, что интеграл

оо

(2.6.23)

-00

по существу есть аналог тензора S(x) в плоской задаче о де­ формации и сложном сдвиге (в безразмерных координатах <£) однородной среды с модулями упругости С°, причем нормаль к плоскости деформирования (£,,£3) направлена вдоль оси

. При этом тензор J(9,хо) совпадает с тензором напряже-

80

 

 

ний в точке

£, = - cos в,

4з= ~ s*n когда вдоль положи­

тельной оси

задан скачок вектора перемещений, уменьша­

ющийся по

закону

Соответствующая плотность

дислокационных моментов т имеет вид

 

. £з) = А* (*■И » (Х.

> 0 ,

^ (^ .^ 2 .^ з) = 0» 6 < ° -

(2.6.24)

Тензор J(0,xc) можно представить в

виде суммы трех

тензоров, соответствующих трем компонентам вектора Р(хо) в

осях УХ,У2,УЪ(не суммировать по / !)

J = J '+ J 2+J\ Jiafi(e,xo) = safaj(e)nx(xo)/3i(xo). (2.6.25)

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

ны характеризуют коэффициентами напряжений К ^К ^К щ

[144]. Связь этих

коэффициентов с компонентами тензоров

/ ( 0 , х ) дается соотношениями

к,(х.)=

кпоо=

К ш( х ) =

(2.6.26)

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

тора /?0О -

§ 2.7. Эллиптическая трещина

Пусть Q - плоский, эллиптический в плане разрез в уп­ ругой среде. В этом случае справедлив следующий аналог теоремы о полиномиальной консервативности.

Если внепшее поле а (х) является полиномом степени т, то вектор Ь(х) скачка перемещений на эллиптической тре­ щине имеет вид

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