книги / Метод эффективного поля в механике композитных материалов
..pdf71
относительно b\ придем к следующим формулам для сим
метричной Ь\ф и антисимметричной Ь\ар\ частей этого тен зора:
fyap) —h-apXfPxn> \ар\ ~ Ь[ар\> |
|
|
(2.5.20) |
||
Л '= — Е2 + |
5С" |
Е1 |
1 Е2Л , С |
= С с + C j . |
|
3с |
2С + З С |
3 |
у |
|
|
Отсюда, в частности, |
следует, |
что |
поле |
КрХр = £^(х) не |
возмущается сферической неоднородностью, если Ь"ар - анти
симметричный тензор.
Если среда изотропна, а источники поля отсутствуют, то
внешнее поле £а(х) удовлетворяет уравнению V a£„(x)=0. В
случае линейного внешнего поля отсюда следует, что tr b°=
=Ь°аа=0. С учетом этого равенства соотношение для симмет
ричной части тензора Ьхупрощается и принимает вид,
5С0 |
h° |
(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 |
|||
7С |
||||||
|
|
|
' |
зс, у 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 - плоский, эллиптический в плане разрез в уп ругой среде. В этом случае справедлив следующий аналог теоремы о полиномиальной консервативности.
Если внепшее поле а (х) является полиномом степени т, то вектор Ь(х) скачка перемещений на эллиптической тре щине имеет вид