книги / Метод эффективного поля в механике композитных материалов
..pdf141
непрерывно дифференцируемую ограниченную функцию b(t) по формуле
Ь(П |
(3.6.5) |
|
J( t - n 2 |
||
Оt - n |
которая является следствием регуляризации (3.2.26).
Вид асимптотики решения уравнения (3.6.4) при Н О и асимптотики (3.6.1) в окрестности контура Г совпадают [148].
Для построения асимптотики b(t) подействуем на обе час ти (3.6.4) оператором преобразования Меллина (2.8.5). Учиты
вая равенство |
|
|
°° ,s-1 |
(1 -5 ) л: (/'Г , *'>о, |
(3.6.6) |
|
tg(*s— 1) Я"
где интеграл понимается в смысле (3.6.5) и существует при О < R e5< 2, из (3.6.4) получим соотношение
b’ (s - q ) + } S. |
r b \ s - \ ) = f ( s ) , |
(3.6.7) |
tg (5 — 1) Я" |
|
|
Здесь b*(s) и f*(s) |
- преобразования Меллина решения и |
правой части уравнения (3.6.4).
Будем рассматривать непрерывные ограниченные решения (3.6.4), убывающие на бесконечности быстрее любой положи
тельной степени Г 1. Преобразование Меллина b*(s) такого решения является аналитической функцией в правой полу
плоскости (Re5>0) комплексной плоскости S. Функцию f*(s)
в (3.6.7) будем считать аналитической в полосе q < Re5 < 2
(q < 2). Тогда в полосе q< Re5 < 2 интегралы Меллина (2.8.5) от общих частей (3.6.4) имеют смысл.
Далее ограничимся исследованием асимптотики (3.6.4) при t-+0. Множество краев включений, которые описываются функцией h(x) вида (3.6.2), можно разделить на три качест венно разных класса: затупленный край, для которого 0<<7<1,
142
остроконечный край q= 1 и край, имеющий точку возврата
?>1 .
Рассмотрим каждый из этих случаев.
1°. Затупленный край (0 <<7< 1). Введем замену перемен ных (s - q ) —>s и перепишем соотношение (3.6.7) в форме
* > ) = — ( r r ’ t / ’ C s + O - ^ - t f + i ) ] - (3.6.8)
S7C
Для выяснения вида асимптотики функции b(t) при t —>О контур интегрирования в формуле обращения преобразования Меллина (2.8.5) возьмем в виде прямой Re5 = т(0 < т< 1/2), параллельной мнимой оси комплексной плоскости S . Как из вестно [135], поведение функции b(t) при / - > 0 определяется
особыми точками b*(s), лежащими слева от контура интегри рования и имеющими наименьшие по абсолютной величине вещественные части. Ближайшей к контуру интегрирования особой точкой множителя (tgS7t)/S7t в правой части (3.6.8) является простой полюс при 5=-1/2. Нетрудно убедиться, что
в полуплоскости Re s >-1/2 функция b*(s) не имеет особых точек. Действительно, если предположить, например, наличие
полюса у функции b*(s) при s = so, где -l/2 < R eso<0, то из
соотношения (3.6.8) следует, что полюса b*(s) должны быть также в точках
**=■*. ± * (1 -0 ), |
*= 1,2,.... |
(3.6.9) |
Но тогда предположение об аналитичности Ъ*(s) в полу |
||
плоскости R es> 0 будет нарушено. |
|
|
Поскольку f(t) - |
бесконечно |
дифференцируемая функ |
ция, то в полуплоскости Re 5 < 0 ее преобразование Меллина f*(s) имеет простые полюса в точках s - 0 ,-1 ,-2 ,.... Таким образом, полюса функции /*(5 + 1) в (3.6.8) расположены в точках s= -1 ,-2,... и компенсируются нулями функции tg5^.
|
|
|
143 |
Поэтому ближайший к прямой |
Re 5 =-1/2 полюс |
функции |
|
b*(s) находится из условия |
|
|
|
s - q + 1 = —1/2. |
(3.6.10) |
||
Если q * 0, то это полюс первого порядка, а при |
q = 0 - |
||
второго. |
|
|
|
Заметим, что наличие полюса функции b*(s) |
в |
точке |
|
S - - M 2 в силу (3.6.8) порождает систему полюсов этой функ |
|||
ции в точках |
|
|
|
* (!- ? )• |
(3-6.11) |
Поскольку в рассматриваемом случае (1 -д )> 0, то все эти
полюса расположены в полуплоскости Re s < - 1 / 2.
Итак, поскольку ближайшие к контуру интегрирования Res= тполюса b’ (s) находятся в точках s = - 1/2 и s= q -(3/2),
то решение уравнения (3.6.4) при »0 представляется в фор ме
m |
+ * ,('). |
A = R e s » -(i)L _ w , (3.6.12) |
|
|
fO (/3/2ln/) |
при |
q = 0 |
l U |
_ [ 0 ( /3/2"‘') |
при |
q> 0 |
2*. Остроконечный край (q=l). В этом случае соотноше
ния (3.6.8) можно разрешить относительно К (s):
vjr |
|
b-(s) = H-l( s ) f \ s + 1), Hap = 8сф+ ------ rap. (3.6.13) |
|
tg S7T |
|
Следовательно, полюса функции V (s), лежащие слева от |
|
контура инте1рирования и не совпадающие |
с полюсами |
/*(5 + 1), определяются из условия |
|
dettf(5) = 0 . |
(3.6.14) |
Используя определение тензора Т°(3.6.4) и явное выраже |
|
ние для тензора S(x) (1.2.9), можно показать, |
что в случае |
144
изотропных матрицы и включений полюсами функции К (s) в
полуплоскости R e J < 0 являются ненулевые корни независи мых уравнений
* |
(1 + аг„) |
s |
( I - ® » ) |
t g J ^ = - |
- ----- - t - S , |
t g S 7 T = - - , tgS7C = |
------------------S , |
|
$ |
$ |
$ |
hjt. |
£ > 0 , |
0 <агс < 1. |
(3.6.15) |
|
|
|
Первые два корня Sjk)(/ = 1,2) каждого из этих трех (£=1Д 3) уравнений, расположенные в левой полуплоскости, удовлет воряют условиям
|
- 2 < £ ' < - % . |
(3.6.16) |
Таким образом, при q = 1 асимптотика решения уравнения |
||
(3.6.4) |
при t —>0 имеет вид |
|
|
b(t) = f J/3kf s' ) + 0 (Г '* ), |
(3.6.17) |
|
к=\ |
|
где |
- корни трансцендентных уравнений (3.6.15), ле |
жащие в левой полуплоскости и наиболее близкие к мнимой оси, Рк - некоторые коэффициенты.
3\ Края, имеющие точку возврата (\<q<2). Введем за мену ( s - q ) —>s и перепишем (3.6.8) в виде
b\s) = f \ s + q ) - ( S- l + <l)* Г У ( д - 1 + ? ) . (3.6.18)
t g ( s - l + q)7r
Рассмотрим особые точки правой части слева от контура интегрирования Res= г, 0 < г < 1 /2 . Очевидно, полюсом яв ляется точка s - —q, а в силу (3.6.8), полюса имеются только в точках
|
145 |
sk = ~ + k(\-q), к =1,2,... . |
(3.6.19) |
Заметим, что точка s = —1 не является полюсом функции |
|
Ъ*(s), так как при этом полюс выражения |
(s-<7+l) |
в правой части (3.6.8) компенсируется нулем множителя tgsn.
Как и в случае Г, здесь можно предположить, что Ъ*(s)
имеет полюс при 5 = - 1 / 2 . |
Но тогда в силу (3.6.8) |
функция |
|
Ъ* (s) должна иметь полюса в точках |
|
|
|
5, = - \ - k ( \ - q ) , |
к= 1,2,... . |
(3.6.20) |
|
Поскольку в данном случае ( 1 - # ) < 0 , |
то такая |
функция |
не будет аналитической в правой полуплоскости. Поэтому в случае существования непрерывного ограниченного решения
уравнения (3.6.4) полюс функции tgsn при s - - 1 / 2 должен
компенсироваться нулем функции /* (s +1) - Ъ* (s - q +1).
Таким образом, при \<q<2 асимптотика b(t) при t —> 0
определяется равенством |
|
b(t) = /32tq + 0 (t2ql). |
(3.6.21) |
Полученные результаты позволяют утверждать, что асим птотика решения исходного уравнения (3.6.1) в окрестности
точки хо е Г при гладкой ограниченной правой части имеет следующий вид
Ъ(х) = Р(х0)Гг+Ь1(г), |
|
|
(3.6.22) |
|
1/2, |
0<q<\ |
0(rV2-q), |
0 < <7 < 1 |
|
у = <ух, |
q=\ , *,(/■) = |
0(гГг), |
q= 1 |
(3.6.23) |
q, |
\<q<2 |
0(г2ч~х), |
\<q<2 |
|
Здесь параметр q определяет асимптотику формы включе
ния в окрестности края £2 по (3.6.2), У\,у2~ коРни трансцен
10 1937
146
дентных уравнений вида (3.6.15), /?(хо) - гладкая на Г функ ция. Аналогичные значения для частных значений параметра q в (3.6.2) получены в [130].
Рассмотрим теперь асимптотику решения уравнения (3.3.18) у края тонкого жесткого включения, предполагая, что
асимптотика функции h(x) определена в (3.6.8). Модельное уравнение, аналогичное (3.6.4) будет в данном случае иметь вид
+ |
1> 0 , (3.6.24) |
00 |
|
и°<фхм= W a p vpJ © |
r^( l , x 2, 0 ) 0 ^ (n)dx2,(3.6.25) |
- 0 0 |
|
М=К'0(п)ВЩп), /а,(О = (м)а0хм® ^ г Ш рХ(),
где интегральный оператор в (3.6.24) определен правой частью (3.6.5).
Тем же путем, что и выше, можно показать, что для затуп
ленных краев (0 < q < 1) асимптотика <т(/) при t —> О имеет вид
ст(/)=Д|'я+а,(0, Д(О =0(/м1пО, <? =о,
CTl( 0 = O (/m _,) > 0 <<? < 1 |
(3.6.26) |
|
В случае остроконечного включения (^=1) тензор сг(/) |
||
имеет следующую асимптотику: |
|
|
oit) = 'Z P T ? + |
), |
(3.6.27) |
/=1 |
|
|
где 5,(,), s(p - ближайшие к мнимой оси и лежащие в левой по луплоскости R es< 0 корни трансцендентного уравнения
det[©(«)+ |
t/* ] = 0. |
(3.6.28) |
147
В случае изотропных среды и включения это уравнение эквивалентно следующим двум независимым уравнениям
t g j f f = - 7^2 ■8 °^Д, |
tgS7t= |
~ |
•(3.6.29) |
Ч |
|
% |
Км. |
В случае включений, |
имеющих |
точку |
возврата (q> 1), |
асимтотика ~a(t) имеет вид |
|
|
|
ст(0 = P f + 0 (t2q~l). |
|
(3.6.30) |
Перейдем к рассмотрению асимптотики поля напряжений в среде в окрестности краев тонких включений. При этом воспользуемся интегральными представлениями (3.3.11), (3.3.18) полей напряжений в среде с тонким включением. Поскольку указанные представления являются внешним пре
делом решения задачи о тонком включении при дх, д2—>0, то имеется в виду промежуточная асимптотика на расстояниях от контура Г , больших по сравнению с характерным попереч ным размером включения И, но малых по сравнению с ос тальными его линейными размерами.
Начнем с выражения (3.3.11) для напряжения в среде с тонким податливым включением и рассмотрим асимптотику ст(х) вне поверхности Q в окрестности ее кромки Г. Пусть
У\>Уг>Уз - локальные декартовы координаты в точке причем ось уг направлена по предельной к П нормали в точ
ке хе, ось у2- по касательной к Г, тогда ось у, лежит в каса
тельной к Q плоскости в точке хо. Учитывая (3.6.22), имеем асимптотику вектора Ь(Х) в окрестности точки хое Г
6 0 ) = Р(х0К |
+ 0 О Г ) > |
Г1> Г > 0 • |
(3.6.31) |
Используя (3.3.11), (3.6.31), запишем выражение для тен |
|||
зора напряжений в |
точке у с |
координатами |
ух= -г cos q>, |
у2~ 0, уг = -г sin (р, где г - расстояние от точки у до хо, (р - полярный угол в плоскости (у,,Уз):
148
<Ку) = гГ~' J ^ (co s^ + ^ ,^ 2,sin ^ + ^ 3) x (3.6.32)
0 (0
хп(г£)/3(х0)%с1П4+ 0(cro), £ =y,/r.
Здесь учтено, что S(X) - однородная степени (-3) четная функция. При т—> 0 интеграл в (3.6.32) стремится к конечно
му пределу и, следовательно, при у < 1 напряжения на Г име
ют особенность типа гу~х. Если у> 1, то в окрестности края напряжения ограничены.
Аналогично случаю трещины (см. §2.6), введем тензорный коэффициент интенсивности напряжений J(<p,x0) соотноше нием
J(<p,х0) = lim rl~rcr(y) при г —» 0 . |
(3.6.33) |
|
Учитывая, что при г —> 0 поверхность П(г) в координатах |
||
£ переходит в полуплоскость (£3=0, £,>()), |
для компонент |
|
тензора J(<p,xo) получим выражение |
|
|
•/«*(?>>*-) = |
(?>)«л (*. )*„ (* .), |
(3.6.34) |
s(<p) = J J S ( c o s , ^2,sin <p)d%2 , (3.6.35)
О—оо
где п(хо) - предельное значение нормали к Q в точке хо е Г .
Таким образом, функция J((p,xo) представима в виде двух
сомножителей, первый из которых S (<р)П(хо) определяется только локальной формой и ориентацией края тонкого вклю чения в точке хо. Второй сомножитель - вектор /?(хо) есть функционал всей поверхности Q , формы включения и внеш
него поля cr(x).
Аналогичным путем из (3.3.18) можно получить вид асим птотики напряжений в окрестности края тонкого жесткого
включения |
|
а(х) = rr~'J(<p,xo) + 0(<J ), |
(3.6.36) |
149
где у - показатель степени, определяющий вид асимптотики
Ь(х) в окрестности края Q . Тензор интенсивности напряже
ний J(<p,x„) в этом соотношении имеет вид |
|
Jap(<P>X*) = -SafiXM ((PWxuvpPvpix.) , |
(3.6.37) |
где тензор S (<р) определен соотношением (3.6.35), а коэффи
циент Р(хо) определяет асимптотику ст(Х) в окрестности точ
ки е Г (д(х) = Р(хо)гг+...).
§ 3.7. Тонкие включения эллипсоидальной формы
Рассмотрим случай тонких включений эллипсоидальной формы. При этом срединная поверхность включения имеет
форму эллипса с полуосями ах,а2, а функция h(X) - пред ставляется в виде
h(x) = 2hz(x), z(xi,x2) = yj l - { x l/а, f ~{x2/a2f , (3.7.1)
где хх,х2 - декартовы координаты, связанные с главными ося ми эллипса Q , ах,а2 - его полуоси.
Обозначим через Vh(X) - характеристическую функцию эл
липсоидальной области с полуосями занятую включе
нием, и пусть Рт(Х) - есть полином степени т по координа там хх,х2,хг. В силу теоремы о полиномиальной консерватив
ности эллипсоида [82] интегральный оператор К с ядром К(Х) (2.1.10), переводит функцию
Рт,н(х) = РЛх) К ( * № Г 1 |
(3.7.2) |
в полином степени т внутри эллипсоидальной |
области Vh. |
Очевидно, что этим же свойством обладают операторы U и Т |
|
с ядрами |
|
U(x) = 0(и)К (х)0(и), T(x) = nS(x)n, |
(3.7.3) |
150
где п - нормаль к срединной поверхности эллипсоида, 0 (и )- проектор (2.2.7). Операторы Ти U в уравнениях (3.3.15),
(3.3.24) для тонких включений связаны с операторами Т и U соотношениями
(ТЬ)(х) = (ТЬП)(х) = | Г (х -х 'Ж л :')П (х ')^ ', х еП,
(Ua)(x) = (UoQ)(x) = J U (х - x')a(x')Q(x')dx'. (3.7.4) Символы U и Т имеют вид
С/*(£) = 0(л)К *(к)0(п), Г (к ) = -nS*(к)п, (3.7.5)
и являются однородными функциями нулевой степени по к.
Поэтому U и Т ограничены, а следовательно, и замкнуты в пространстве обобщенных функций Соболева-Слободецкого
H(R3) [148]. Отсюда следуют равенства
jim UPmh = U(liraPmh), |
lim TPmh = T(lim Pmh), |
(3.7.6) |
||
Л—>0 |
h->0 |
h->0 |
/J—>0 |
|
предел в правых частях которых имеет вид |
|
|||
|
jim Км (х ) = Pm(x)z(x)Q(x). |
|
(3.7.7) |
Поскольку в области Vhлевые части (3.7.6) являются поли номами степени тп для всех h , то правые части для предель ного образа функции Pmh являются полиномами на П . Отсю
да и из (3.7.4) получаем следующее свойство операторов Ти U в уравнениях (3.3.15), (3.3.24).
Если Q - эллиптическая поверхность, то операторы Ти U переводят полином степени ш, домноженный на функцию
z(x) вида (3.7.1), в полином той же степени на Q .
Отсюда следует, что в случае тонкого эллиптического включения и полиномиального внешнего поля решения урав нений (3.3.15), (3.3.24) определяются выражениями
Ь(х) = В (x)z(x) , а(х) = Q (x)z(x) , |
(3.7.8) |
где В(х) и Q(x) - полиномы той же степени Ш, что и внеш нее поле. В частности, если сг (или е ) - постоянный тензор,