Устинов К.Б. О расслоении полосы по границе раздела упругих свойств. Часть 1. Постановка задачи, случай нормального отрыва // Вестник Пермского национального исследовательского политехнического университета. Механика. -2015. - № 4. - С. 226-245. DOI: 10.15593/perm.mech/2015.4.13
Ustinov K.B. On delamination of a stipe along the boundary between two elastic layers. Part 1. Problem formulation, the case of normal crack. PNRPUMechanics Bulletin. 2015. No. 4. Рр. 226-245. DOI: 10.15593/perm.mech/2015.4.13
ВЕСТНИК ПНИПУ. МЕХАНИКА № 4,2015
PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
DOI: 10.15593/perm.mech/2015.4.13 УДК 593.3
О РАССЛОЕНИИ ПОЛОСЫ ПО ГРАНИЦЕ РАЗДЕЛА УПРУГИХ СВОЙСТВ ЧАСТЬ 1. ПОСТАНОВКА ЗАДАЧИ, СЛУЧАЙ НОРМАЛЬНОГО ОТРЫВА
К.Б. Устинов
Институт проблем механики им. А.Ю. Ишлинского РАН. Москва, Россия
О СТАТЬЕ
АННОТАЦИЯ
Получена: 8 ноября 2015 г. Принята: 17 ноября 2015 г. Опубликована: 25 декабря 2015 г.
Ключевые слова: отслоение, интерфейсная трещина, факторизация, упругая заделка
Рассмотрена задача о полосе, состоящей из двух изотропных упругих полос различной толщины, обладающих различными свойствами, разделенных полубесконечной трещиной, проходящей по границе раздела упругих свойств и нагружаемой на бесконечности системой усилий с ненулевым главным вектором и моментом. Путем применения преобразования Лапласа задача сведена к однородной матричной задаче Римана. В предположении возможности пренебрежения влиянием нормальных напряжений на сдвиговые смещения и сдвиговых напряжений на нормальные смещения задача сведена к двум скалярным задачам Римана. Данная постановка может рассматриваться как приближенная для общего случая (данное приближение при этом заведомо не хуже общепринятого приближения, заключающегося в рассмотрении узкого слоя в рамках теории балок либо стержней) и как точная - для случая, когда прилегающие слои могут скользить друг относительно друга, но удерживаются силами адгезии (например, вандерваальсовыми). Путем факторизации получено точное аналитическое решение одной из указанных задач, а именно задачи о нормальном отрыве. Получены асимптотические выражения для смещений берегов трещины вдали от ее вершины. Показано, что ведущие члены асимптотики смещений берегов трещины вдали от ее вершины соответствуют смещению балки при граничных условиях типа упругой заделки, т.е условиям пропорциональности смещения и угла поворота в точке заделки вектору усилия и изгибающему моменту посредством матрицы коэффициентов податливости.Для данных коэффициентов получены аналитические выражения. Также получены асимптотические выражения для поля напряжений вблизи вершины трещины (коэффициент интенсивности напряжений).
©ПНИПУ
© Устинов Константин Борисович - кандидат физико-математических наук, доцент, старший научный сотрудник лаборатории, e-mail: ustinov@ipmnet.ru
Konstantin B. Ustinov - Ph.D. of Physical and Mathematical Sciences, Associate Professor, Senior researcher, e-mail: ustinov@ipmnet.ru
226
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
ON DELAMINATION OF A STIPE ALONG THE BOUNDARY BETWEEN TWO ELASTIC LAYERS
PART 1, PROBLEM FORMULATION, THE CASE OF NORMAL CRACK K.B. Ustinov
A.Yu. Ishlinsky Institute for Problems in Mechanics RAS. Moscow, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 8 November 2015 Accepted: 17 November 2015 Published: 25 December 2015
Keywords: delamination, interface crack, factorization, elastic clamping
The problem of a strip, composed by two isotropic elastic layers of different elastic properties and thicknesses, separated by a semi-infinite crack located along the line between the layers, is considered. The mechanical load with nonzero total force and moment is supposed to be applied at infinity. By means of Laplace transformation the problem is reduced to a homogeneous Riemann problem. Under the assumption of possibility to neglect the cross-terms related to the influence of the normal stresses to the shier displacements and the shier stresses to the normal displacements the problem is reduced to two scalar Riemann problems. Such a formulation may be considered as an approximation for the general case (which is not worse than the traditional beam or rode approximation) and as the exact one for the case, when the two layers may slide but may not separate due to cohesion, e.g. by van-der-Waals forces. By means of factorization procedure the exact analytical solution has been obtained for one of the formulated scalar problems, namely, the problem of the normal separation. The asymptotical expression has been derived for the relative displacements of the crack faces far from its tip. It is shown that the leading asymptotic terms of these relative displacements correspond to a beam deflection under the boundary condition of the type of generalized elastic clamping. i.e. the proportionality of the displacement and angle of rotation of the clamping point to the total vector and bending moment of the applied load by means of the matrix of coefficients of compliance. The analytical expressions for these coefficients have been obtained. The asymptotical expression for the stress field near the crack tip (stress intensity factor) was also derived.
©PNRPU
Введение
Задача продольного расслоения полосы, составленной из двух полос меньшей толщины, обладающих различными упругими свойствами, решалась многими авторами (например, [1-5]), при этом в основном интерес был сосредоточен на поведении решения вблизи вершины трещины. Однако в ряде задач, например при исследовании потери устойчивости отслаивающимися покрытиями, значение имеет поведение поля смещения вдали от вершины трещины [6-18]. Одной из наиболее простых конфигурация является полоса, составленная из двух полос меньшей толщины. Если длина отслоения много больше толщины полосы, для ряда случаев отслоение может рассматриваться как полубесконечное. Если толщина одного из слоев много больше протяженности отслоения, то такой слой может рассматриваться как полуплоскость; эффективное и элегантное решение задачи для такого случая, полученное путем применения преобразования Лапласа и сведения ее к матричной задаче Римана, было дано в работах [19-22] для одинаковых упругих свойств слоя и полуплоскости и в работе [17] для различающихся свойств. Ключевым моментом решения являлась факторизация матричного коэффициента. В указанных работах получены асимптотические представления для трансформант поля напряжений на линии продолжения трещины, и в конечном счете вычислены коэффициенты интенсивности напряжений (КИН) в зависимости от главных вектора и момента напряжений, действующих на продолжении трещины.
227
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
Однако интерес к задаче не ограничивается вычислением КИН. В частности, из указанного решения можно извлечь асимптотику смещений берегов трещины вдали от ее вершины [14]. При этом при наличии ненулевых главных вектора и момента действующих нагрузок ведущие члены асимптотики смещений берегов трещины будут, очевидно, соответствовать смещения балки (пластины), подверженной действию указанных главных вектора и момента при некоторых граничных условиях. Следует подчеркнуть, что эти граничные условия вовсе не являются условиями жесткой заделки. Они представляют собой [7, 10-15, 17] условия обобщенной упругой заделки, т.е. условия пропорциональности смещений и угла поворота в точке заделки действующим компонентам главных вектора и момента нагрузки. В работе [17] данное решение было обобщено на случай отслоения и основания, обладающих различными упругими свойствами, хотя и связанных некоторыми дополнительными ограничениями. В работе [16] рассматривалась возможность пренебрежения влиянием нормальных к линии трещины напряжений на сдвиговые компоненты смещения (и сдвиговых компонент напряжения на нормальные компоненты смещения), в результате было получено приближенное решение о сдвиговой трещине на границе раздела. Там же было показано, что подобное решение по крайней мере не хуже традиционного подхода - замены отслаивающейся части полосы одномерным объектом (стержнем). В настоящей статье рассматривается аналогичная задача для расслаиваемой составной полосы.
1. Общая постановка задачи о расслоении полосы, состоящей из двух упругих слоев по границе контакта
В условиях плоской деформации рассматривается неоднородная упругая полоса 1 < y <—h, составленная из двух изотропных частей с различными упругими свойствами, контакт которых проходит вдоль линии y — 0. Вдоль части границы y — 0, x < 0 имеется полный контакт, вдоль оставшейся части границы y — 0, x > 0 контакт отсутствует (рис. 1). Все величины, относящиеся к нижней полосе 1 < y <—h, будут обозначаться индексом 1, все величины, относящиеся к верхней полосе 0 < y < 1, индексом 2. Модули Юнга и коэффициенты Пуассона (модифицированные для условий плоской деформации) материалов полуплоскости и полосы обозначим E), v(i), i — 1,2 соответственно. С обычными модулями Юнга E0(i) и коэффициентами Пуассона v0(i) они связаны соотношениями
E(i) — ■
E (i)
(1 — v0" ') (1 + v0('1 у
vnl — ■
(i)
1 — v0" >■
(1)
В постановке плоской деформации рассматривается однородная задача: предполагается, что все поверхности свободны от напряжений,
аyy — аxy — 0, при y — 1, y — —h, и при y — 0, x > 0,
xy
(2)
а нагрузка с эквивалентными главным вектором (T, N) и моментом M приложена на бесконечности, так что
M — | xа^ (x,0)dx, N — { а^ (x,0)dx, — T — { аxy (x,0)dx.
(3)
228
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
Условия сопряжения на границе имеют вид
3_ II 3_ II to
yy yy xy xy
u(1) = u(2), v(1) = v(2) при y = °,
Здесь u, v - компоненты вектора смещения; axx, ayy, axy
женин.
x < 0. (4)
- компоненты тензора напря-
A y=\ ' У
crvi(.v), crt,.(y) £i2\vl2>
y = -h L *
X
Рис. 1. Расслоение составного слоя. Геометрия задачи
В окрестности нуля поле напряжений должно быть интегрируемо.
Из механики разрушения известно, что это поле имеет корневую особенность:
(5)
(6)
Однако для построения решения можно довольствоваться условием (5), а условие (6) получится из решения. В случае различных значений упругих постоянных здесь могут присутствовать осциллирующие члены.
Воспользуемся результатами работы [22], в которой предлагается связь между образами двустороннего преобразования Лапласа от производной компонент вектора смещения
ГО ^
g(xl) (p) = I ^(l) (x,°)pxdx
—ГО ГО CS
gf (p ) = I ^v'1 > (x,0)e-pxdx
(7)
и напряжениями
ГО
qx (p )= |a xy (x,°)e ~pxdx,
qy (p )= |a yy (x,°)e_pxdx
(8)
229
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
для полосы 0 < y < 1 (при поставленных граничных условиях свободной внешней границы ayy = axy = 0, при y =1). Вывод данных выражений дан в Приложении 1.
E(2) ~2~
gx2) (p)
E(2) ~2~
g y2)
(p )
sin p cos p - p
dn
n 2 1 -v(2)
qx (p)+J^qy (p)+—2 qy 22
1 -v(2) / 4 p2 / ч sin p cos p + p
qx (p)^ qx (p)+-
d
d-,
qy (p),
p € L
(9)
1 • 2 2 d2 = sin p - p .
(10)
Здесь L — мнимая ось.
Рассмотрим преобразование Лапласа от следующих величин:
F+(p )
F
-
: 0 дх
(p )= }
:>(x,°)-v <1’(x-0)le- (2) (x,0)-u(1) (x,0)J (11)
. ayy (x’0)le-p>x axy (x,0)J (12)
Здесь использовано свойство равенства нулю скачка смещения для отрицательной полуоси и равенства нулю напряжений для положительной полуоси. На основании (4) подынтегральное выражение (11) тождественно равно нулю для x < 0, и, следовательно, F+( p)
аналитична в правой полуплоскости (Re p > 0). Аналогично согласно (2) подынтегральное выражение (12) равно нулю для x > 0 , и, следовательно, F-(p) аналитична в левой полуплоскости (Rep < 0). Здесь аналогично [19-22] предполагается выполнение условия
убывания напряжений на бесконечности при x ^ -да . В окрестности нуля выполняются условия (6). При x ^ +да смещения могут расти как полином третей степени (что соответствует ненулевой поперечной силе на бесконечности), следовательно, в Лаплас-образе F+ (p) возможно появление полюса в нуле до третьего порядка включительно.
Комбинируя (7)—(12), получаем уравнение матричной задачи Римана, аналогичное [19-22]:
F+ (p) = K(p)F (p), (13)
где
K (p ) =
a11 a1l''
V a21 a22
, p € L,
(14)
11 =(sinpcosp + p)/d2 + p(sinhp coshp + hp)/d1,
12 = -a21 = -p2 / d2 + p (h p)2 / d1 -p'/ 2,
(15)
a
22
=(
sin p cos p - p)/d2 + p(sin h p cos h p - hp)/d1,
230
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
(16)
ц = E(2)/ E(1), р' = 1 - v(2) — р(1 - v(1) ).
р' = 1 — v(2)
(17)
(18)
Основная сложность состоит в факторизации матричного коэффициента K (p), т.е. в представлении его в виде
ного р. det X± (p)^ 0 в соответствующих плоскостях вплоть до границы. В настоящее
время общее точное аналитическое решение указанной задачи неизвестно; можно указать несколько путей построения решения:
- применение аппроксимаций для матричного коэффициента (14) задачи методами [23-26];
- построение приближенных решений в предположении возможностью пренебречь перекрестными членами в (14) (положить а12 = а21 = 0) [16];
- построение решения для частных случаев, для которых факторизация может быть осуществлена [27-34].
Далее решение будет получено с использованием второго из указанных методов.
2. Задача о трещине нормального отрыва. Постановка
С помощью предположения о возможности пренебрежения в (14) влиянием перекрестных членов, что соответствует а12 = -а21 = 0, матричная задача распадается на две скалярные. Данное приближение заведомо не хуже традиционной замены слоя однородным объектом - стержнем либо балкой [35-37]. Рассмотрим задачу о нормальном отрыве.
K (p ) = X+ (p )X* (р), (19)
где функции X± (p) аналитичны в правой и левой полуплоскости комплексного перемен-
F+(p) = K(p)F—(p), p e L
(20)
(21)
0
(22)
K (p) = (sin p cos p + p)/ d2 +p(sin h p cos h p + hp)/ d1,
d1 = sin2 (hp} — (hp)2 , d2 = sin2 p — p2.
Аналог (19) становится скалярным уравнением
(23)
(24)
231
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
После нахождения Л± (р) окончательное решение дается с помощью теоремы Лиу-вилля
f+(p)=Л+(р)П(р)’ Re(p)>0’
F-(p)^-(p)n(p), Re(p)< 0. (25)
Здесь П(р) - функция, подлежащая определению.
Второе из уравнений (3) и условия (6) для интересующих нас величин после трансформации принимают соответственно вид
F_(p ) = N + Mp + o (р), Re р ^ 0 -. (26)
F-(p )=Т2р K+ O (р (,‘2}
Re р ^ -да.
(27)
Из элементарных соображений также следует, что производные от нормальных смещения для нижней кромки верхней полосы при x ^ да будут расти как
(\2Mx - 6 Nx2)/ E(2), а нижней - как -{\2Mx - 6 Nx 2)/ (e(1)), что в Лаплас-образах согласно (17) даст
F+(p } = 6 (1 + л / hJ)
N
3
р 2
Re р ^ 0 +.
(28)
Условие отсутствия скачка смещений в нуле в трансформантах записывается следующим образом
F+(p } = O (р ~к}, k < 1, Re р ^-да. (29)
Данные условия будут использованы для определения П( р }.
3. Решение задачи Римана
Уравнения (23) могут быть записаны в виде (здесь L - мнимая ось) Л-1 (р)Л+(р}=(1 + Л)^3 р G2 (р}, Р е L,
G2 (p)=----tg3 Р
1 + п
sin р cos р + р + ^ sin hp cos hp + hp
sin р - р
sin2 hp - (hp)2
, Р е L.
(30)
(31)
С помощью стандартного приема факторизации арктангенса через гамма-функцию решение может быть выражено следующим образом:
, \ I----- Г3 (1 + р / л) , ч
Л+(Р )^/1 + ЛГ3 ( } J2+ (Р}, (32)
Г ^1/2 + р / л)
, ч 1 Г3 (1/2 - р / л) , ч
Л- (Р) = -1Т=т37—f) J2- (Р) V1 + л Г (-р / л)
(33)
232
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
J2±(Р )_ eXP
I №
— J ln
2 к
——th3 s
1 + p
sh s ch s + s sh hs ch hs + hs
sh2 s - s2 + P sh2 hs - (hs)2
ds
is - p
(34)
Для получения решения осталось определить функцию П( p), число независимых параметров в которой определяется индексом задачи и поведением решения в опорных точках. Наличие опорных точек только в нуле и на бесконечности определяет вид П(р)
как сумму слагаемых с положительными и отрицательными степенями р. Более детальное исследование показывает, что наиболее общим видом является
П(Р)
(35)
Здесь D2,D3 - постоянные, определяемые из граничных условий. Наличие членов со степенями р большими минус двух привело бы к нарушению условия (29), а меньшими минус трех - к нарушению условия (28).
Для определения постоянных D2, D3 потребуется разложение функции Л+(р) вблизи нуля (два члена). Однако для определения параметров упругой заделки необходимо четыре члена. Для внеинтегрального члена (32) имеем
Г3 (1 + р / к) Г3 (1/2 + р/к)
к
4К
л 6ln2 (18ln2 2
1 +----Р + ---2“
к V к
\ р2 + r 6 6 ln3 2-Ф)] 3ln2~]
V к3 к J
J
(36)
Здесь ^(х) - дзета-функция Римана, ^(3)« 1,20206.
Для интегрального члена (32) из (34), (97), (101), (103), положив t _ is , имеем
J2+ (р )_ h±1/2 (0 Ш + I1P +
112 +1"
2 2 1
2 + 1 р + 6
13 + 2I1I2 +113
(37)
h (s)
—^th3 s
1 + р
sh s ch s + s sh hs ch hs + hs
sh2 s - s2 Л sh2 hs - (hs)2
h (0)
6 ( + h3) h3 (1 + p)
1 j d ln h (s) ds _ 1 j h'(s) ds _ 1 j h (s) ds
2к J ds s 2к J h (s) s 2к^ h (0) s2’
2
h" (0) 2 (-ph2 + 4h3)
h (0) 5 (p + h3)
(38)
(39)
(40)
(41)
233
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
1 rd3lnh(s) ds _ 1 rf h(s) 3 h'(s)h"(s) . h,3(s)^
2nr ds3 s 2n h (s) h2 (s) h3 (s)
ds
s
(42)
Выражение для интеграла I1 (40) можно преобразовать следующим образом:
.. <х>
1J ln
n
Я 0
из (40), (38), (43) следует
1 bo £ 1 3 ds 3 f — _— 1 ln 1 >> +-> 1 ds _ 3 f ■ 1 1"
_ s _ s2 2n •' —f _ s _ s2 n r —f _ sh s s _
ds_ _ 6 ^ (-1)n _ 6ln2 s n n_1 n n
I1 _ I1A
6ln2
n
I1A _ ^
1 r d ln hA (s) ds _ 1 r hA'(s) ds _ 1 r hA (s) ds
2nГ ds s 2пГ hA (s) s 2пГ hA (0) s2’
hA (s)_ s"
1 + 2
sh s ch s + s sh hs ch hs + hs -----------+ 2---------------
sh2 s - sz
sh2 hs - (hs)2
(43)
(44)
(45)
(46)
Подстановка выражения для функции n(p) (35) и разложения вблизи нуля функции Л+ (p) (результата подстановки (36)-(39) в (32)) в первое из равенств (25) дает разложение Лаплас-образа производной от скачка смещения вблизи нуля. С точностью до двух ведущих членов оно имеет вид
F+(p)_ — ^ ^t*2/ — [ D2 n + D3 (nI1 + 6ln2)] + o (p “2 ) Re P ^ 0 +. (47)
p n n p 4 '
Сравнение (47) с (28) позволяет вычислить значение констант, входящих в П( р):
D2 1 + 2/h3 [nM + (nI1 + 6ln2)],
D3 _-n372^! + 2/h3N.
(48)
Таким образом, задача решена: подстановка значений констант (48) в (35), а затем вместе с (32)-(34) в (25) дает Лаплас-образы напряжений и производной от разности смещений вдоль границы раздела упругих свойств.
4. Определение асимптотики смещения берегов разреза вдали от вершины трещины
Лаплас-образ указанной асимптотики может быть получен путем подстановки (48), (35), (32), (36)-(39) в (25). Результат данной подстановки может быть записан в виде
F+ (р)_ 6 ( + 2 / h3
M N
2 3
V P P )
Ku M M + K^N К,л№ + K.^N
2 P
f
I
V
u N _ + _ uM
u N
- + O(p), Rep ^ 0 +. (49)
KNM _ 12(1 + 2/h3)l I1 + — _ 12(1 + 2/h1 )i,a.
v 4 n ) v !
(50)
234
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
Ku N = 6 О + Л/ h3 )
, I2 (Т 61п2л2
1 —2 + I, +--
2 V п )
= 6 (1 + л / h3)
1 -I +I 2
1 +11Л
(51)
KuM ~
= 6 (1 + л / h3)
л I2 (г 61п2л
-1 + — + L +----
2 V п )
= 6 (1 + л / h3)
-1 + ^ + 2
и
(52)
Ku N =
= 2 (1 + л / h3)
)-13 - 36С3(3)+ 2 (г 61п2 ^ I1 + 3 = 2 (1 + л / h3) ) - 36«(3) + 21 3 13 3 ~|~3Л
п V п ) п
(53)
Первый член в (49) соответствует прогибу балки, нагружаемой силой и моментом. Оставшиеся члены могут интерпретироваться как поправки к элементарной балочной теории (типа модели балки Тимошенко, учитывающей вклад в смещения от поперечной силы) и дополнительные смещения, вызываемые граничными условиями, отличающимися от условий жесткой заделки. Величины, определяемые формулами (50)-(53), являются при этом коэффициентами пропорциональности между поперечной силой и моментом, с одной стороны, и углом поворота и нормальным смещением - с другой.
Поскольку из приведенного решения непосредственно невозможно выделить указанные причины и при этом естественно ожидать, что матрица податливости должна быть симметрична, для модели эквивалентной балки естественно принять граничное условие в симметричном виде:
Kи NN = KuM = 6 (1 + Л / h )/. (54)
Зависимость KuM от h для различных отношений модулей представлена на рис. 2.
K um
Рис. 2. Зависимость коэффициента упругой заделки KuM от относительной толщины основания h для различных отношений модулей: сплошная линия - ц = 0,5 ; пунктирная - ц = 1; точечная - ц = 5 ; штрих-пунктирная - ц = 10
235
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
Для толстых оснований (h ^да) справедлива асимптотика
Ku M _ 6I1A •
hA (s)= -
+ P
sh s ch s + s „
-T2----Г + 2P
sh s - s
(55)
(56)
Здесь I1A по-прежнему определяется (45). Зависимость Ku M от отношения модулей р для больших h (h ^да) представлена на рис. 3. Данная зависимость неплохо аппроксимируется формулой
Ku M = 1 + 7,5р1/3.
(57)
K u'M
Рис. 3. Зависимость коэффициента упругой заделки KuM от отношения модулей р для больших толщин подложек (( ^да)
5. Определение параметров разрушения
При p ^ -да для (33) справедливо следующее асимптотическое разложение:
Л-()=-Щ+р+о(рГ) (58)
Подстановка (58), (35), (48) в (25) и разложение для больших отрицательных p ^ -да
дает
F_(p)^6^V£[M + Ih<N ] + o (-p-1'2). (59)
Последнее выражение немедленно дает асимптотику напряжений при x ^-0
а
yy
7^1 + р/ h3 х/гё]Хл/-+'П
[M + I1AN ] + o (“1/2).
(60)
236
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
Отсюда КИН
(61)
Скорость высвобождения энергии [38, 39] при этом имеет вид
(62)
Полученное выражение для скорости высвобождения энергии для р = 1, h = 1 совпадает с элементарным решением (например, [2]). При рассмотрении задачи о распространении интерфейсной трещины в классической постановке приведенный анализ для определения параметров разрушения может рассматриваться лишь как приближенный, поскольку в общем случае в точном решении будут присутствовать обе моды. Однако задача имеет и самостоятельное значение [13] для случая, когда прилегающие слои могут скользить друг относительно друга, но удерживаются силами адгезии (например, вандер-ваальсовыми).
Заключение
Получено и исследовано решение однородной задачи о полубесконечной трещине, проходящей вдоль границы, разделяющей два упругих слоя с отличающимися свойствами. Путем применения преобразования Лапласа задача сведена к однородной матричной задаче Римана. В предположении возможности пренебрежения влиянием нормальных напряжений на сдвиговые смещения и сдвиговых напряжений на нормальные смещения задача сведена к двум скалярным задачам Римана. Получена асимптотика поля напряжений вблизи вершины трещины, имеющего корневую особенность при отсутствии осциллирующих членов. Посчитан коэффициент интенсивности напряжений. Показано, что ведущие члены асимптотики смещений берегов трещины вдали от вершины трещины соответствуют смещению балки при граничных условиях типа упругой заделки, то есть условиях пропорциональности смещения в точке заделки вектору усилия и изгибающему моменту.
Работа выполнена при финансовой поддержке Программы Президиума РАН № I.3.
Приложение 1. Вывод связи между производной от скачка смещения и напряжениями на границе слоя
Рассмотрим плоскодеформированное состояние полосы у е[0,1] в декартовой системе координат, верхняя граница у = 1 которой свободна от напряжений.
ауу (x,1) = °*y (x,1) = 0
xy
(63)
а к нижней границе у = 0 приложена некоторая система нагрузок
(64)
237
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
yy
(x,0) = qy (x).
(65)
Здесь qx (x), qy (x) - заданные функции. В основном изложение следует работам [19-22]. Найдем связь между действующей на нижней границе полосы системой нагрузок
и производными от компонент смещений — (x,0), — (x,0).
dx dx
Полная система уравнений теории упругости включает:
1) уравнения равновесия в отсутствие объемных сил
da да
xy
dx
да
xy
dy
да
yy
dx dy
- 0,
- 0,
(66)
(67)
где аxx , аyy , аxy
компоненты тензора напряжений;
2) уравнения закона Гука, выраженные через компоненты смещений u, v:
E(2) du -(
dx
-(xx +ayy )-
E
(2)
du + dv Kdy dx ^
)-(1 + V(2) )d yy ,
2 ^ + V(2) ) xy ,
где E
(2) ,.(2)
— модуль Юнга и коэффициент Пуассона полосы;
d
2
f d2
Kdx2 dy2 у
(axx +аyy )- 0.
(68)
(69)
3) уравнение совместности, выраженное через компоненты тензора напряжений
(70)
Для решения поставленной задачи подвергнем уравнение (69) следующим преобра-
du
зованиям. Подставим значение — из уравнения (68) в продифференцированное по x
dx
уравнение (69).
E(2)
d2v / r^\f . da
dx2
(1 + V121)
xy
da
+ -
yy
dx dy
d(
a xx + a yy
dy
(71)
Подстановка сюда уравнений равновесия (67) дает
E <2) dE-(1+V (2) -d(+a7
dx2 v ' dy dy
(72)
dx2 v ! dy dy
Для решения системы уравнений (66), (67), (68), (72), (70) используем двустороннее преобразование Лапласа, определяющее образ f (p) функции f (x, y),
238
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
да f (p> у) = f f (x у Кp%dx- —да (73)
Обратное преобразование при этом определяется как
f (x У) ——2^7 f f (Р> У)epxdp, (74)
где направление обхода по контуру L, соответствующему мнимой оси, осуществляется сверху вниз. Трансформанты от уравнений (66), (67), (68), (72), (70) и граничных условий (63), (64), (65) примут вид
- ( ),d&xy (p, У) А ^ (p’y)+ dy = 0> (75)
а ( ч d а уу tp, У) п p xy(p у )+ dy = 0 (76)
( d2 ^ p + г 2 [axx (p, у ) + ауу (p, у)] — 0, V ау у (77)
Е(2)й'^ у) — <тxx ^ у) + ауу ^ у) — ( + v(2)) <туу ^ у), (78)
pE(2)v'^ у) = p (l + v(2)) axx ^, у)— d [axx ^, у ) + ауу CР, у )], (79)
аxy (p,l) — ауу Cp,l) — 0 (80)
аxy (p,0) — 4x (p) (81)
ауу (Л0) = Чу (p). (82)
Здесь , дй . 8v й =—, v =—, dx dx а (83)
да 4x (p) = f аxy (x,0)e^Р^dx,
да Чу (p)= f аyy (x,0) e~pXdx —да (84)
есть трансформанты от напряжений на нижней границе. Решение системы обыкновенных дифференциальных уравнений (75)-(79) с граничными условиями (80)-(82) дает выражения для трансформант производных смещения через действующие на нижнем крае полосы (у = 0) напряжения. Для поставленной задачи интерес
представляют значения трансформант производных смещения на этой же линии (у = 0):
239
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
E(2) 2
E(2) 2
g™ (p) g? (p)
n2 1 -v(2)
d qx (p)+Ejqy (p)+~2qy 2
sm p cos p - p
1 - v(2) / 4 p2 / ч sin p cos p + p
2 q (p)~qx (p)+‘
d
qy (p ),
(85)
d = sin2 p - p2. (86)
Здесь
gX2) (p )= I (x,0)e_ pxdx>
» я (87)
g У2) (p ) = 1-xv (x,0)e~pxdx.
Коэффициенты в полученных выражениях аналитичны на всей мнимой оси за исключением нуля, где они имеют полюса до третьего порядка включительно. Полюс третьего порядка в нуле трансформанты соответствует квадратичному росту оригинала на бесконечности. Следовательно, из полученных выражений, в частности, следует, что при отсутствии особенностей для напряжений производные от смещения могут расти как квадрат расстояния, а сами смещения как куб расстояния. Для полосы такое поведение является общеизвестным фактом, в частности, в рамках элементарной балочной теории, являющейся при определенных условиях строгой асимптотикой точной теории, смещения вдали от приложенной системы сил с ненулевым главным вектором растут как третья степень расстояния.
Приложение 2. Разложение факторизующих функций вблизи нуля
Пусть существует факторизация функции h (t)
h (t ) = H+ (z)/ H-(z), (88)
Hi(z) = ^i1 В. dt (89)
Пусть функция h(t) четная h(-t) = h(t). Тогда функция lnh(t)тоже четная. Согласно формулам Сохоцкого-Племеля предельные значения H+(t0 ) выражаются следующим
образом:
ln H+ (0)- ln H- (0 ) = ln h (0), (90)
ln H+(t0) + ln H (t0 )= 1 Iln h(t) dt. + v ’ v ’ niJ t-10 (91)
Из последнего равенства с учетом четности ln h (t) для t0 = 0 ln H+ (0) = - ln H- (0). Подставляя полученное равенство в (90), получаем следует
240
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
Н±_ м±( 0)_л11'2 (о).
(92)
Для получения следующих членов разложения вычислим k-ю логарифмическую производную от обеих частей равенства (88):
d -lnH+(z)-ттlnH-(z)_ттlnh(z).
(93)
j k ' \ ' _j k - \ / ik
dz dz dz
На основании формул Сохоцкого-Племеля предельные значения функций
£_
dZ
-ln Н+ (z) в нуле следующие:
Тт ln н±(z )=±i dLln h (z )+ik,
dz 2 dz
(94)
_ 1 . dk ln h (t) dt
k 2ni dtk t ’
(95)
причем если функция h (t) четная, то для всех нечетных k в правой части равенства (94)
будут обнуляться первые члены, а для всех четных k - вторые.
Выпишем несколько первых членов получаемого таким образом разложения. Для первой логарифмической производной
Н +'(°) _ Н-'(°) _ 1 Н+(0) Н-(о) 1
1 ^ d ln h (t) dt _ 1 r h'(t) dt 2niГ dt t 2niГ h (t) t
(96)
(97)
Для второй логарифмической производной
H± "(0) -H’±2(0) _±
H±(0) H±2 (0) ~2h(0).
С учетом (96) данное выражение преобразуется к виду
H±'' (0) _± hM +I j
H±(0) 2h(0) 1'
Для третьей производной получаем
H± ’’’ (0) - 3 H '± (0 )H± ’’ (0) + 2 H ’±3 (0) _ r H±(0) H±2 (0) H±3 (0) 3
I3 _
1 . d3 ln h (t) dt _ 1 . 2ni dt3 t 2ni
h-(t) -3 h’(t)h-(t) 2k3)3 h (t) h2 (t) h3 (t) t
С учетом (96) и (99) выражение (100) преобразуется к виду
H± ’’’ (0) _ j ± 3 h ’’ (0) т + ,3 H±(0) 3 2 h(0) 1 1'
(98)
(99)
(100)
(101)
(102)
241
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
В результате для разложения факторизующих функций вблизи нуля может использоваться следующее представление:
H±(z ) = h±yl (0 ф + Z,z +1
1 гч _ + £! +1 1 2 л 1 z л— 3 h"(0) 3 i3 ± ;/I +13
2h(0) 1 6 3 2 h(0) 1 1
z3 1 +
O(z4). (103)
Данный вывод не представляет сложности, однако в литературе (например, [40, 41]) обычно ограничиваются лишь двумя членами разложения.
Библиографический список
1. Suo Z., Hutchinson J.W. Interface crack between two elastic layers // Int. J. Fract. - 1990. -Vol. 43. - P. 1-18.
2. Hutchinson J.W., Suo Z. Mixed Mode Cracking in Layered Materials // Adv. Appl. Mech. / eds. J.W. Hutchinson, T.Y. Wu. - 1992. - Vol. 29. - P. 63-191.
3. Александров В.М., Мхитарян C.M. Контактные задачи для тел с тонкими покрытиями и прослойками. - М.: Наука. 1983. - 487 с.
4. Григолюк Э.И., Толкачев В.М. Контактные задачи теории пластин и оболочек. - М.: Машиностроение. - 1980. - 415 с.
5. Греков М.А. Сингулярная плоская задача теории упругости. - СПб.: Изд-во C.-Петерб. ун-та, 2001. - 192 с.
6. Ustinov K.B., Dyskin A.V., Germanovich L.N. Asymptotic analysis of extensive crack growth parallel to free boundary // 3rd Int. Conf. Localized Damage - 94. - 1994. - P. 623-630.
7. Cotterell B., Chen Z. Buckling and cracking of thin film on compliant substrates under compression // Int. J. Frac. - 2000. - Vol. 104. - No. 2. - P. 169-179.
8. Dyskin A.V., Germanovich L.N., Ustinov K.B. Asymptotic analysis of crack interaction with free boundary // Int. J. Solids Structures. - 2000. - Vol. 37. - P. 857-886.
9. Yu H.H., He M.Y., Hutchinson J.W. Edge effects in thin film delamination // Acta Mater. - 2001. -Vol. 49. - P. 93-107.
10. Yu H.-H., Hutchinson J.W. Influence of substrate compliance on buckling delamination of thin films // Int. J. Fract. - 2002. - Vol. 113. - P. 39-55.
11. Effect of substrate compliance on the global unilateral post-buckling of coatings: AFM observation and finite element calculations / G. Parry, J. Colin, C. Coupeau, F. Foucher, A. Cimetiere, J. Grilhe // Acta materialia. - 2005. - Vol. 53. - P. 441-447.
12. Гольдштейн P.B., Устинов К.Б., Ченцов А.В. Оценка влияния податливости подложки на напряжения потери устойчивости отслоившегося покрытия // Вычисл. мех. спл. сред. - 2011. -Т. 4, № 3. - С. 48-57.
13. Салганик Р.Л., Устинов К.Б. Задача об упруго заделанной пластине, моделирующей частично отслоившееся от подложки покрытие (плоская деформация) // Известия РАН МТТ. - 2012. -№ 4. - C. 50-62.
14. Устинов К.Б. Еще раз к задаче о полуплоскости, ослабленной полубесконечной трещиной, параллельной границе // Вестник ПНИПУ. Механика. - 2013. - № 4. - C. 138-168.
15. Ustinov K.B. On influence of substrate compliance on delamination and buckling of coatings // Eng. Fail. Analys. - 2015. - Vol. 48B. - P. 338-344.
16. Устинов К.Б.. О сдвиговом отслоении тонкой полосы от полуплоскости // Изв. PAH. МТТ. - 2014. - № 6. - С. 141-152.
17. Устинов К.Б.. Об отслоении слоя от полуплоскости; условия упругой заделки для пластины, эквивалентной слою // Изв. PAH. МТТ. - 2015. - № 1. - С. 75-95.
242
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
18. Устинов К.Б., Каспарова Е.А. Оценка влияния кривизны и податливости основания на параметры отслоения покрытия // Деформация и разрушение материалов. - 2015. - № 3 - С. 28-35.
19. Златин А.Н., Храпков A.A. Полубесконечная трещина, параллельная границе упругой полуплоскости // Докл. АН СССР. - 1986. - Т. 31. - С. 1009-1010.
20. Златин А.Н., Храпков A.A. Упругая полуплоскость, ослабленная трещиной, параллельной ее границе // Исследования по упругости и пластичности. Проблемы современной механики разрушения. - 1990. - T. 16. - С. 68-75.
21. Златин А.Н., Храпков A.A. Векторная задача Римана с ненулевым индексом показателя матрицы-коэффициента // Изв. ВНИИГ им. Б.Е. Веденеева. - 1985. - Т. 181. - С. 12-16.
22. Khrapkov A.A. Wiener-Hopf method in mixed elasticity theory problems. - S.-P., 2001.
23. Ворович И.И., Александров В.М., Бабешко В.А. Неклассические смешанные задачи теории упругости. - М.: Наука, 1974. - 456 с.
24. Ворович И.И., Бабешко В.А. Динамические смешанные задачи теории упругости для неклассических областей. - М., 1979. - 320 с.
25. Бабешко В.А. Обобщенный метод факторизации в пространственных динамических смешанных задачах теории упругости. - М.: Наука, 1984. - 256 с.
26. Abrahams I.D. On the non-commutative of factorization of Wienner-Hopf kernels of Khrapkov type // Proc. Roy. Soc. London A. - 1998. - Vol. 454. - Р. 1719-1743.
27. Heins A.E. System of Wiener-Hopf equations // Proceeding of Symposia in Applied mathematics II. - McGraw-Hill, 1950. - Р. 76-81.
28. Чеботарев Г.Н., К решению в замкнутой форме краевой задачи Римана для систем n пар функций // Учен. зап. Казан. ун-та. - 1956. - T. 116б. - Кн. 4. - С. 31-58.
29. Daniele V.G. On the factorization of Wiener-Hopf matrices in problem solvable with Hurd’s method // Trans ANTENNAS propagate. - 1978. - Vol. 26. - Р. 614-616.
30. Jones D.S. Commutative Wiener-Hopf factorization of a matrix // Proc. R. Soc. A. - 1984. -Vol. 393. - Р. 185-192.
31. Моисеев Н.Г. О факторизации матриц-функций специального вида // Докл. АН СССР. -1989. - Т. 305, № 1. - C. 44-47.
32. Antipov Y.A., Moiseev N.G. Exact solution of the plane problem for a composite plane with a cut across the boundary between two media // J. Appl. Math. Mech. - 1991. - Vol. 55. - Р. 531-539.
33. Antipov D.A., Silvestrov V.V. Factorization on a Riemann surface in scuttering theory // QR. J. Mech. Appl. Math. - 2002. - Vol. 55. - Р. 607-654.
34. Antipov D.A., Silvestrov V.V. Vector functional difference equation in electromagnetic scuttering // IMA J. Appl. Math. - 2004. - Vol. 69. - No. 1. - Р. 27-69.
35. Koiter W.T. On the diffusion of load from a stiffener into a sheet // Quart. J. Mech. Appl. Math. -1955. - Vol. 8. - Part 2. - P. 164-178.
36. Alblas J.B., Kuypers W.J.J. On the diffusion of load from a stiffener into an infinite wedgeshaped plate // Applied Scientific Research, Section A. - 1965-1966. - Vol. 15. - Iss. 1. - P. 429-439.
37. Каландия А.И. Математические методы двумерной упругости. - М. Наука, 1973. - 303 с.
38. Салганик Р.Л. О хрупком разрушении склеенных тел // ПММ. - 1963. - Т. 27. - № 5. -C. 957-962.
39. Malyshev B.M., Salganik R.L. The strength of adhesive joints using the theory of crack // Int. J. Fract. Mech. - 1965. - Vol. 1. - No. 2. - C. 114-128.
40. Мусхелишвили Н.И. Сингулярные интегральные уравнения. - 2-е изд. - М.: Наука, 1962. - 600 с.
41. Гахов Ф.Д. Краевые задачи. - М.: Наука, 1977. - 640 с.
243
Ustinov K.B. /PNRPUMechanics Bulletin 4 (2015) 226-245
References
1. Suo Z., Hutchinson J.W. Interface crack between two elastic layers. Int. J. Fract. 1990, vol. 43,
pp. 1-18.
2. Hutchinson J.W., Suo Z. Mixed Mode Cracking in Layered Materials. Adv. Appl. Mech. eds. J.W. Hutchinson, T.Y. Wu, 1992, vol. 29, pp. 63-191.
3. Alexandrov V.M., Mkhitaryan S.M., Kontaktnye zadachi dlya tel s tonkimi pokrytiyami i prosloikami [Contact Problems for Bodies with Thin Coatings and Interlayers]. Moscow: Nauka, 1983. 487 p.
4. Grigoluk E.I., Tolkachev V.M. Kontaktnye zadachi teorii plastin I obolochek [Contact problems of plate and shell theory]. Moscow: Mashinostroenie, 1980. 415 p.
5. Grekov M.A. Singulyarnaya ploskaya zadacha teorii uprugosti [Singular Plane Problem of Elasticity]. St. Petersburgskii universitet, 2001. 192 p.
6. Ustinov K.B., Dyskin A.V., Germanovich L.N. Asymptotic analysis of extensive crack growth parallel to free boundary. 3rd Int. Conf. Localized Damage 94, 1994, pp. 623-630.
7. Cotterell B., Chen Z. Buckling and cracking of thin film on compliant substrates under compression. Int. J. Frac., 2000, vol. 104, no. 2, pp. 169-179.
8. Dyskin A.V., Germanovich L.N., Ustinov K.B. Asymptotic analysis of crack interaction with free boundary. Int. J. Solids Structures, 2000, vol. 37, pp. 857-886.
9. Yu H.H., He M.Y., Hutchinson J.W. Edge effects in thin film delamination. Acta Mater., 2001, vol. 49, pp. 93-107.
10. Yu H.-H., Hutchinson J.W. Influence of substrate compliance on buckling delamination of thin films. Int. J. Fract., 2002, vol. 113, pp. 39-55.
11. Parry G., Colin J., Coupeau C., Foucher F., Cimetiere A., Grilhe J., Effect of substrate compliance on the global unilateral post-buckling of coatings: AFM observation and finite element calculations. Acta materialia, 2005, vol. 53, pp. 441-447.
12. Goldstein R.V., Ustinov K.B., Chentsov A.V. Influence of substrate compliance on stresses generating stability loss of a delaminated coating. Comput. continuum mech., 2011, vol. 4, no. 3, pp. 48-57.
13. Salganik R.L., Ustinov K.B. Deformation problem for an elastically fixed plate modeling a coating partially delaminated from the substrate (plane strain). Mech. of Solids., 2012, vol. 47, no. 4, pp. 415-425.
14. Ustinov K.B. Once more to the problem on a half-plane, weakened by a semi-infinite crack, parallel to its boundary. PNRPU Mechanics Bulletin, 2013, no. 4, pp. 138-168.
15. Ustinov K.B. On influence of substrate compliance on delamination and buckling of coatings. Eng. Fail. Analys., 2015, vol. 48B, pp. 338-344.
16. Ustinov K.B., On Shear Separation of a Thin Strip from the Half-Plane. Mech. Solids., 2014, vol. 49, no. 6, pp. 713-724.
17. Ustinov K.B., On separation of a layer from the half-plane: elastic fixation conditions for a plate equivalent to the layer. Mech. Solids, 2015, vol. 50, no. 1, pp. 62-80.
18. Ustinov K.B., Kasparova E.A. Otsenka vliianiia krivizny i podatlivosti osnovaniia na parametry otsloeniia pokrytiia [Influence assessment of curvature and yielding of substrate on parameters of coating delamination]. Deformatsiia i Razrushenie Materialov, 2015, no. 3, pp. 28-35.
19. Zlatin A.N., Khrapkov A.A. A Semi-Infinite Crack Parallel to the Boundary of the Elastic HalfPlane. Doklady Akademii nauk SSSR, 1986, vol. 31, pp. 1009-1010.
20. Zlatin A.N., Khrapkov A.A. Uprugaia poluploskost, oslablennaia treshchinoi, parallelnoi eyo granitse [Elastic half-plane weakened by a crack parallel to its boundary], Issledovaniya po uprugosti i plastichnosti. Problemy sovremennoi mekhaniki razrusheniia, 1990, vol. 16, pp. 68-75.
21. Zlatin A.N., Khrapkov A.A. Vektornaia zadacha Rimana s nenulevym indeksom pokazatelia matritsy-koeffitsienta [Vector Riemann Problem with Zero Index of Matrix-Coefficient Exponent]. Izvestiia nauchno-issledovatelskogo instituta gidrotekhniki imeni B.E. Vedeneeva, 1985, vol. 181, pp. 12-16.
244
Устинов К.Б. /Вестник ПНИПУ. Механика 4 (2015) 226-245
22. Khrapkov A.A. Wiener-Hopf method in mixed elasticity theory problems. Saint Petersburg, 2001.
23. Vorovich I.I., И.И., Aleksandrov V.M., Babeshko V.A. Neklassicheskie smeshannye zadachi teorii uprugosti [Nonclassical mixed problems in theory of elasticity]. Moscow: Nauka, 1974. 456 p.
24. Vorovich I.I., Babeshko V.A. Dinamicheskie smeshannye zadachi teorii uprugosti dlia neklassicheskikh oblastei [Dynamical mixed problems of elasticity theory for non-classical domains]. Moscow: Nauka, 1979. 320 p.
25. Babeshko V.A. Obobshchennyi metod faktorizatsii v prostranstvennykh smeshannykh zadachakh teorii uprugosti [Generalized method of factorization in 3-D dynamical mixed problems of theory of elasticity]. Moscow: Nauka, 1984. 256 p.
26. Abrahams I.D. On the non-commutative of factorization of Wienner-Hopf kernels of Khrapkov type. Proc. Roy. Soc. London A., 1998, vol. 454, pp. 1719-1743.
27. Heins A.E. System of Wiener-Hopf equations. Proceeding of Symposia in Applied mathematics II. McGraw-Hill. 1950, pp. 76-81.
28. Chebotarev G.N. K resheniiu v zamknutoi forme kraevoi zadachi Rimana dlia system n par funktsii [To Solution of Riemann Boundary-Value Problem in Closed Form for Systems of n Pairs of Functions]. Uchenye zapiski Kazanskogo universiteta, 1956, vol. 1166, no. 4, pp. 31-58.
29. Daniele V.G. On the factorization of Wiener-Hopf matrices in problem solvable with Hurd’s method. Trans ANTENNAS propagate, 1978, vol. 26, pp. 614-616.
30. Jones D.S. Commutative Wiener-Hopf factorization of a matrix. Proc. R. Soc. A., 1984, vol. 393, pp. 185-192.
31. Moiseev N.G., Factorization of Matrix Functions of Special Form. Sov. Math. Dokl., 1989, vol. 39, pp. 264-267.
32. Antipov Y.A., Moiseev N.G. Exact solution of the plane problem for a composite plane with a cut across the boundary between two media. J. Appl. Math. Mech., 1991, vol. 55, pp. 531-539.
33. Antipov D.A., Silvestrov V.V. Factorization on a Riemann surface in scuttering theory. QR. J. Mech. Appl. Math., 2002, vol. 55, pp. 607-654.
34. Antipov D.A., Silvestrov V.V. Vector functional difference equation in electromagnetic scuttering. IMA Journal of Applied Math., 2004, vol. 69, no. 1, pp. 27-69.
35. Koiter W.T. On the diffusion of load from a stiffener into a sheet. Quart. J. Mech. Appl. Math., 1955, vol. 8, part 2, pp. 164-178.
36. Alblas J.B., Kuypers W.J.J. On the diffusion of load from a stiffener into an infinite wedgeshaped plate. Appl. Sci. Res., Sec. A., 1965-1966, vol. 15, iss. 1, pp. 429-439.
37. Kalandiia A.I. Matematicheskie metody dvumernoi uprugosti [Mathematical Methods of TwoDimensional Elasticity]. Moscow: Nauka, 1973. 303 p.
38. Salganik R. L., The Brittle Fracture of Cemented Bodies. J. Appl. Math. Mech., 1963, vol. 27, no. 5, pp. 1468-1478.
39. Malyshev B.M., Salganik R.L. The strength of adhesive joints using the theory of crack. Int. J. Frac. Mech., 1965, vol. 1, no. 2, pp. 114-128.
40. Muskhelishvili N.I. Singular integral equations; boundary problems of function theory and their application to mathematical physics. Groningen: P. Noordhoff, 1953. 600 p.
41. Gakhov F.D. Boundary Value Problems. International Series of Monographs in Pure and Applied Mathematics, vol. 85. Oxford/London/Edinburgh/New York/Paris/Frankfurt, Pergamon Press, 1966. 564 p.
245