CC BY
32
Устинов К.Б. О расслоении полосы по границе раздела упругих свойств. Часть 2. Случай сдвиговой трещины // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2016. - №2. -С. 131-142. DOI: 10.15593/perm.mech/ 2016.1.09
Ustinov K.B. On delamination of a strip along the boundary between two elastic layers. Part 2, Сase of shear crack. PNRPUMechanics Bulletin. 2016. No. 2. Рр. 131-142. DOI: 10.15593/perm.mech/2016.2.09
ВЕСТНИК ПНИПУ. МЕХАНИКА № 2, 2016 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
Б01 10.15593/регт.тесЬ/2016.2.09 УДК 593.3
О РАССЛОЕНИИ ПОЛОСЫ ПО ГРАНИЦЕ РАЗДЕЛА УПРУГИХ СВОЙСТВ ЧАСТЬ 2. СЛУЧАЙ СДВИГОВОЙ ТРЕЩИНЫ
К. Б. Устинов
Институт проблем механики им. А.Ю. Ишлинского РАН. Москва, Россия Московский авиационный институт, Москва, Россия
О СТАТЬЕ АННОТАЦИЯ
Рассматривается задача о полосе, образованной двумя изотропными упругими полосами меньшей и в общем случае различной толщины, обладающими различными упругими свойствами, с полубесконечной трещиной, проходящей по границе раздела упругих свойств и нагружаемой на бесконечности системой усилий. В первой части работы [1] представлена математическая формулировка, и с помощью преобразования Лапласа задача сведена к однородной задаче Римана-Гильберта с матричным коэффициентом. В предположении возможности пренебрежения влиянием нормальных напряжений на сдвиговые смещения и сдвиговых напряжений на нормальные смещения задача сведена к двум скалярным задачам Римана-Гильберта. Данная постановка может рассматриваться как приближенная для общего случая (данное приближение при этом заведомо не хуже общепринятого приближения, заключающегося в рассмотрении узкого слоя в рамках теории балок либо стержней) и как точная - для случая, когда прилегающие слои могут скользить друг относительно друга, но удерживаются силами адгезии от нормального отрыва.
Путем факторизации получено точное аналитическое решение одной из сформулированных в [1] задач, а именно задачи о сдвиге. Получены асимптотические выражения для смещений берегов трещины вдали от ее вершины. Показано, что ведущие члены асимптотики смещений берегов трещины вдали от вершины соответствуют смещению стержня при граничных условиях типа упругой заделки, т.е. условиях пропорциональности смещения в точке заделки приложенному усилию. Для данного коэффициента пропорциональности получено аналитическое выражение. Также получены асимптотические выражения для поля напряжений вблизи вершины трещины (коэффициент интенсивности напряжений и скорость высвобождения энергии).
©ПНИПУ
Получена: 22 марта 2016 г. Принята: 25 мая 2016 г. Опубликована: 30 июня 2016 г.
Ключевые слова: отслоение, интерфейсная трещина, факторизация, упругая заделка
© Устинов Константин Борисович - доктор физико-математических наук, доцент, старший научный сотрудник, e-mail: ustinov@ipmnet.ru
Konstantin B. Ustinov - Doctor of Physical and Mathematical Sciences, Associate Professor, Senior Researcher, e-mail: ustinov@ipmnet.ru
ON DELAMINATION OF A STRIP ALONG THE BOUNDARY BETWEEN TWO ELASTIC LAYERS. PART 2, СASE OF SHEAR CRACK
K.B. Ustinov
A. Yu. Ishlinsky Institute for Problems in Mechanics RAS. Moscow, Russian Federation Moscow Aviation Institute, Moscow, Russian Federation
ARTICLE INFO ABSTRACT
Received: 22 March 2016 The problem of a strip, composed by two isotropic elastic layers of different elastic
Accepted: 25 May 2016 properties and thicknesses, separated by a semi-infinite crack located along the line
Published: 30 June 2016 between the layers is considered. The mechanical load is supposed to be applied at
- infinity. In the first part of the study [1] the mathematical formulation of the problem and
Keywords. its reduction to a homogeneous Riemann-Hilbert problem by application Laplace trans-
delamination, form was presented. Under the assumption of possibility to neglect the cross-terms relat-
interface crack, ed to the influence of the normal stresses to the shear displacements and the shear
factorization, stresses to the normal displacements, the problem is reduced to two scalar Riemann-
elastic clamping Hilbert 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, where the two layers may slide but may not separate due to cohesion.
By means of factorization procedure the exact analytical solution has been obtained for one of the formulated in [1] scalar problems, namely, the problem of a shear crack. The asymptotical expression has been derived for the relative displacements of the crack faces far from the crack tip. It is shown, that the leading asymptotic terms of these relative displacements correspond to a rode under the boundary condition of the type of elastic clamping. i.e. the proportionality of the displacement of the clamping point to the applied force. The analytical expression for this coefficient has been obtained under the accepted assumptions. The asymptotical expression for the stress field near the crack tip (stress intensity factor and energy release rate) is also derived.
©PNRPU
1. Постановка задачи. Исходная конфигурация
В первой части работы [1] была сформулирована задача теории упругости в постановке плоской деформации о полубесконечной трещине, проходящей вдоль части границы, разделяющей два изотропных упругих слоя с произвольными соотношениями толщин и упругих свойств. Самоуравновешенная нагрузка, характеризуемая главным вектором и изгибающим моментом, предполагалась приложенной на бесконечности со стороны трещины. Путем применения преобразования Лапласа задача была сведена к задаче Ри-мана-Гильберта с матричным коэффициентом. Там же в предположении возможности пренебрежения влиянием нормальных напряжений на сдвиговые смещения и сдвиговых напряжений на нормальные смещения задача сведена к двум скалярным задачам Римана-Гильберта. Одна из задач - задача о нормальном отрыве - была решена.
Далее рассматривается вторая из намеченных задач - задача о сдвиговой трещине. Для решения используется та же постановка [2], а именно - предположение о возможности пренебречь указанным перекрестным влиянием нормальных напряжений на касательные смещения и касательных напряжений на нормальные смещения. Данный подход может рассматриваться как приближенный, причем рассматриваемое приближение заведомо не хуже общепринятого приближения, заключающегося в рассмотрении узкого слоя в рамках теории балок либо стержней [3-7]. Однако существует прототип данной математической формулировки, для которого получаемое решение является точным, а именно
для случая, когда прилегающие слои могут скользить друг относительно друга, но удерживаются силами адгезии от нормального отрыва. В любом случае для граничных условий, получаемых исходя из указанной постановки, решение получено строго. Исходная конфигурация задачи представлена на рис. 1.
>'=1 ■ 7
ГГ,,( V) -)—Г>
У = т
Рис. 1. Расслоение составного слоя сдвиговой трещиной. Геометрия задачи
В условиях плоской деформации рассматривается неоднородная упругая полоса -к < у < 1, составленная из двух изотропных частей с различными упругими свойствами. Вдоль части границы у = 0, х < 0 имеется полный контакт, вдоль оставшейся части границы - у = 0, х > 0 контакт отсутствует. Все величины, относящиеся к нижней полосе 0 < у <-к, обозначаются индексом 1, относящиеся к верхней 0 < у < 1 - индексом 2. Предполагается, что все поверхности свободны от напряжений,
ауу = аху = 0, при у = 1, у = -к, и при у = 0, х > 0, (1)
а нагрузка с эквивалентными главным вектором Т приложена на бесконечности, так что
0
Т = {а ^ (х,0)А. (2)
—ад
Здесь и, V - компоненты вектора смещения; ахх, ауу, аху - компоненты тензора напряжений. Условие (2), безусловно, не охватывает весь класс однородных задач для рассматриваемой конфигурации, оставляя в стороне рассмотрение, например, однородного на-гружения на бесконечности. Однако для большинства подобных задач решения могут быть получены элементарными методами [8]. Условия сопряжения на границе имеют вид
ауу =а®, а® =а£, и(1) = и(2), V® = ^2) при у = 0, х < 0. (3)
Модифицированные для условий плоской деформации модули Юнга и коэффициенты Пуассона материалов обозначим Е(|),у(г), I = 1,2 соответственно. С обычными модулями Юнга Е0(|) и коэффициентами Пуассона у0(г) они связаны соотношениями
Е (1) V (1)
Е (I) =_Е_ V (I) = __(4)
^ Л (I )\Л (I )\> у „ а)
(1 ^0(|) )(1 ^)' "1 — V
(I)
0
Предположение о возможности пренебрежения влиянием нормальных напряжений, действующих на линии продолжения трещины, на сдвиговые смещения равносильно пренебрежению самими нормальными напряжениями. Это приводит к дополнительному условию
о« = g(J = 0 при y = 0, 3x. (5)
Данное условие неявно используется при моделировании узкого слоя стержнем при контакте с другим телом [3-7].
Для каждой из полос выполняются основные уравнения теории упругости, применение преобразования Лапласа к которым приводит к задаче Римана-Гильберта на мнимой оси ( p е L ) (см. Приложение):
F+(p) = K (p) F_(p) , p е L, (6)
77(2) ш p¡
F+ ( p) = Er ¡J- [u(2) (x,0) - u(1) (x, 0)] e" pxdx, (7)
F_(p) = ¡ оxy (x,0) e"pxdx, (8)
—ад
K ( p) = (sin p cos p — p) / d2 + ^ (sin h p cos hp — hp) / d1, (9)
d1 = sin2 (hp) — (hp)2, d2 = sin2 p — p2, (10)
л = E(2)/ E(1). (11)
Кроме того, должны выполняться условия в опорных точках (в рассматриваемом случае - в нуле и на бесконечности), характеризующие возможность наличия, тип и порядок сингулярностей в данных точках. Так, для р, стремящегося к нулю, из уравнения (2) следует
F— (p) = T + o (1), Re p ^ 0 —. (12)
Для р, стремящегося к бесконечности, исходя из наличия корневой особенности в нуле оригинала в напряжениях, на основании теоремы абелевого типа [9] следует
F+(p)=-T=Kii +O(p232), Rep ^—ад. (13)
V-2 p
Кроме того, из элементарных соображений следует, что горизонтальные смещения верхней полосы при x ^ ад будут расти как 4Tx / E(2), а нижней - как —4Tx / hE(1), что в Лаплас-образах согласно [10], (7) даст
p 2(1 + Vh)T + o(1), Rep ^ 0 +. (14)
p
Для решения поставленной задачи необходимо представить ее коэффициент K (p)
в виде произведения двух функций, аналитичных в правой и левой полуплоскостях комплексного переменного р:
K ( p ) = Л—1 ( p )Л+( p ). (15)
(16)
После нахождения Л+ (р) окончательное решение дается с помощью теоремы Лиу-вилля
^+(р) = Л+(р)П(р), Яе(р)> 0, (р ) = Л—(р )П(р), Яе (р )< 0.
Здесь П(р) - функция, подлежащая определению. Поскольку индекс (в определении
[11]) коэффициента краевой задачи, определяемой уравнением (10), равен минус единице, данная функция может содержать лишь одну константу и имеет вид
П( р ) = П 0 рп. (17)
Здесь п - целая постоянная, определяемая характером поведения решения в нуле либо на бесконечности, П0 - постоянная, определяемая из граничных условий.
2. Решение задачи Римана-Гильберта для сдвиговой трещины
Уравнения (10), (15) могут быть записаны в виде
Л:1 (p)Л+ (p) = -(1 + ц)ctg p Gl (p),
p e L,
G1 (p )="-tg p
1 + ц
sin p cos p - p + sin hp cos hp - hp sin2 p -p2 sin2 hp -(hp)2
p eL.
(18) (19)
Здесь Ь - мнимая ось.
Применением стандартного приема факторизации арктангенса через гамма-функцию решение (18) может быть выражено следующим образом:
Г 1 ш J±(p) = expГ-— Jln
1 + ц
th p
sh s ch s - s sh hs ch hs - hs
sh2 s - s2
- + ц-
sh2 hs - (hs)2
ds
is - p
Асимптотические разложения вблизи р = 0 имеют вид
Л- (p ) = i
1
p
yj1 + ц/h V2-
О (p2).
Л+ (p ) = üj 1 + ц / h
V—
1 I 4ln2
1 + 1-+ S1 I p
—
О (p2).
Здесь
S1 (ц) = dlJpl
dp
1 r d ,n
J /7c
— 0 ds
p=0 0
th s
f
1 + ц
sh s ch s - s sh hs ch hs - hs
Y
sh2 s - s2
- + ц-
sh2 hs - (hs)2
ds s
(20) (21)
(22)
(23)
(24)
Интегрирование последнего выражения по частям с учетом свойств интеграла Коши [12, 13] дает
1 ш
S1 (ц) = -J ln
— 0
th s
f
2 (1 + ц / h)
sh s ch s - s sh hs ch hs - hs
\
sh2 s - s2
• + ц-
sh2 hs - (hs)2
ds ~~2'
(25)
Обе формы (24), (25) не представляют сложности для численного интегрирования. Однако для анализа последнее выражение полезно преобразовать, разбив на два интеграла, первый из которых вычисляется в замкнутой форме
51 (л) = 52 Ы + Зз (л) /4,
- •
" Ш £ " ds - ОТ -1 1" " 1 1" ds _2 ^ (- -1)п 21п2
£ £ 2 -ад бЬ £ £ £ - П=1 п -
4
5з (л, И ) = -{ 1п
— г,
(
бЬ s еЬ s - s БЬ еЬ --+л-
Л
бЬ2 s - s2
бЬ2 hs - (Ию )2
2 (1 + л / h)
Из (22) на основании (16) находим
П (р ) = -^42-у1+ц1ИТр 1. Подстановка (29) и (23) в (16) и разложение для р ^ 0 дает
ds
„2 '
(26)
(27)
(28)
(29)
^+(р ) = 2Т (1 + л / h)
- + 53
Р
+ О (р)
(30)
Выражение (30) есть асимптотика образа Лапласа от производной скачка смещения. На основании свойств преобразования Лапласа (например, [10]) образ самого смещения получается делением его на р . Члены разложения высшего порядка могут быть получены дальнейшим умножением на р произвольное число раз, применением обратного преобразования Лапласа с последующим дифференцированием оригинала такое же число раз. Ограничиваясь двумя членами, с учетом (7), окончательно получаем для х ^ ад
Е (2)и = 4Т (1 + л / И ) х + КиТТ, Кит =(1 + л/И)5з (л,и).
(31)
(32)
В формуле (31) первый член соответствует смещению, посчитанному в стержневом приближении. Коэффициент 4 появляется за счет того, что сила приложена не по центрам верхней и нижней полосы (как при расчетах в стержневом приближении), вызывая изгибы, сопровождающиеся дополнительными смещениями нижнего края верхней полосы на величину 3Т / Е(2) и верхнего края нижней полосы на величину 3Т / Е(1). Оставшийся член, деленный на Т , может рассматриваться как коэффициент упругой заделки, т. е. как коэффициент пропорциональности между смещением в точке заделки и приложенным усилием. Численные значения представлены на рис. 2, 3. Для больших л / И соотношение приближенно описывается формулой
КиТ = 0,3И + 0,84л.
(33)
Знание коэффициента упругой заделки полезно при решении задач, связанных с отслоением покрытий [14-19].
0
Рис. 2. Значение коэффициента упругой заделки в зависимости от отношения модулей: сплошная линия - И = 1; пунктирная линия - И = 5; точечная линия - И = 20
Рис. 3. Значение коэффициента упругой заделки в зависимости от И: сплошная линия - л = 0,5; пунктирная - л = 1; точечная - л = 5; штрих-пунктирная - л = 10
3. Определение параметров разрушения
При р ^ -ад для (20) справедливо следующее асимптотическое разложение:
Л-(р )=/'тттЧ+о ((~р).
У-(1 + л)
(34)
Подстановка (34) и (29) в (16) и разложение для больших по абсолютной величине отрицательных р ^ -ад дает
F_(p ) =
T
v^p V
2 (1 + л/ h )
1 + л
+ о I-
(-p1/2 )■
(35)
Последнее выражение немедленно дает асимптотику напряжений при х ^-0
V2T 1l + л/ h
° XV =
V^HV 1
+ о\\х
-1/2
(36)
Отсюда коэффициент интенсивности напряжений (КИН)
(37)
л
С другой стороны, из первого члена (31) находим изменение упругой энергии системы при продвижении трещины:
^=±
ёх
г 1 т 4Т 1 4Т Л 2Т2
/ё^^ФЁ^+л/А). (38)
Сравнение (37) и (38) показывает, что связь между параметрами разрушения имеет традиционный вид [20, 21]
Г = К2^Л= К2 Е(1) +Е(2) (39)
2Е(2) 2Е(1)Е(2) ' (39)
Полученные выражения для КИН (37) и скорости высвобождения энергии (38) для Л = 1, " = 1 совпадает с решением, полученным элементарным методом (например [22]). Для л = 1, " , что соответствует полубесконечной трещине параллельной границе полуплоскости, значение выражения для скорости высвобождения энергии (38) совпадает с выражением, полученным в [23, 24]. Однако в общем случае для тел с классическим условием контакта (3) в точном решении будут присутствовать обе моды, поэтому приведенный анализ для определения параметров разрушения может рассматриваться только как приближенный.
Согласно (36) поле напряжений вблизи вершины трещины описывается простой корневой особенностью: осцилляции не наблюдается. Очевидно, что отсутствие осциллирующих членов следует из принятых в постановке задачи условий.
Заключение
Получено решение однородной задачи о полубесконечной сдвиговой трещине, проходящей вдоль интерфейса, разделяющей два упругих слоя с отличающимися свойствами. В предположении возможности пренебрежения влияния нормальных напряжений на сдвиг с помощью применения преобразования Лапласа задача сведена к задаче Римана-Гильберта. Путем факторизации получены асимптотические выражения для смещений берегов трещины вблизи и вдали от ее вершины. Показано, что ведущие члены асимптотики смещений берегов трещины вдали от вершины трещины соответствуют смещению стержня при граничных условиях типа упругой заделки, т. е. условиях пропорциональности смещения в точке заделки продольному усилию. Получена асимптотика поля напряжений вблизи вершины трещины, имеющего корневую особенность, соответствующая коэффициенту интенсивности напряжений. Показано, что принятые при постановке задачи допущения приводят к отсутствию осциллирующих членов в поле напряжений и смещений при приближении к вершине трещины. Вычислена скорость высвобождения энергии при росте трещины.
Работа выполнена при финансовой поддержке РФФИ, проект 14-01-00855.
Приложение
Вывод основного уравнения рассматриваемой задачи Римана-Гильберта
В декартовой системе координат рассматривается плоскодеформированное состояние полосы у е [0,1], верхняя граница у = 1 которой свободна от напряжений,
ауу ( х,1) = а^ ( х,1) = 0, (40)
а к нижней границе у = 0 приложена система нагрузок
а^ (х,0) = дх (х), (41)
а уу ( х,0) = ду ( х). (42)
Здесь дх (х), ду (х) - заданные функции. Изложение следует в основном в [11].
Полная система уравнений теории упругости включает:
1) уравнения равновесия при отсутствие объемных сил
^ + = 0, (43)
дх ду
= 0. (44)
дх ду
Здесь ахх, а аху - компоненты тензора напряжений;
2) уравнения Закона Гука, выраженные через компоненты смещений и, V,
Е(2) | = ( +ауу )-(1 + у(2)) (45)
Е(2)
^ди + дv ^ ду дх
2 (1 + у (2) )а ху. (46)
Здесь Е(2) , у(2) - модуль Юнга и коэффициент Пуассона полосы;
3) уравнение совместности, выраженное через компоненты тензора напряжений
Лд2 д2
чдх2 ду
+ Т-2 |(ахх +а,у) = 0. (47)
Применим к уравнениям (43)-(47) и граничным условиям (40)-(42) двустороннее преобразование Лапласа, определяющее образ / (р) функции / (х, у) как
ад
I(р,у)= | I(х,у)рхсы. (48)
-ад
Обратное преобразование при этом определяется как
I(х,у) = --1- ГI(р,у)ерЧр (49)
2га \
где направление обхода по контуру Ь, соответствующему мнимой оси, осуществляется сверху вниз.
После этого находим связь между действующими на нижней границе полосы нагруз-
ди dv
ками и производными от компонент смещений —(x,0), — (x,0) (подробнее [1]):
dx дх
^ gf (Р )=^Mfm ^ (p)+d % (p)+,, (p),
d = sin2 p - p2. (51)
Здесь
W си /-Ч
gf (p)=^dd~xu (x,0)e-pxdx, gf (p)= ¡^v (x,0)e-pxdx. (52)
В силу условия (5), которое согласно (48) приводит к обнулению qy (р) = 0, система
(50)-(51) распадается, первое из уравнений (50) преобразуется к скалярному уравнению, связывающему образы касательных напряжений и производной тангенциальных смещений
Е-(р) = ™РС°;Р-Р q, (р). (53)
Проделывая аналогичную процедуру для полосы у е [-А, 0], вычисляем связь Лаплас-образов от величин напряжения и разности смещения (6)—(10). Здесь использовано свойство равенства нулю скачка смещения для отрицательной полуоси и равенства нулю напряжений для положительной полуоси.
Более подробный вывод представлен в [1], где выведено матричное уравнение, которое можно свести к скалярному при выполнении условия (5), т.е. положив в (7) из [1]
К12 = К21 = 0 .
Библиографический список
1. Устинов К.Б. О расслоении полосы по границе раздела упругих свойств. Ч. 1. Постановка задачи, случай нормального отрыва // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. - № 4. - С. 226-245.
2. Устинов К.Б. О сдвиговом отслоении тонкой полосы от полуплоскости // Изв. РАН. МТТ. - 2014. - № 6. - С. 141-152.
3. 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.
4. Попов Г.Я. Изгиб полубесконечной плиты, лежащей на линейно-деформируемом основании // ПММ. - 1961. - № 25, Вып. 2. - С. 342-355.
5. Alblas J.B., Kuypers W.J.J. On the diffusion of load from a stiffener into an infinite wedge-shaped plate // Appl. Sci. Research, Sec. A. - 1965-1966. - Vol. 15. - Iss. 1. - P. 429-439.
6. Каландия А.И. Математические методы двумерной упругости. - М.: Наука, 1973. - 303 с.
7. Салганик Р.Л., Устинов К.Б. Задача об упругозаделанной пластине, моделирующей частично отслоившееся от подложки покрытие (плоская деформация) // Изв. РАН МТТ. - 2012. -№ 4.- C. 50-62.
8. Suo Z., Hutchinson J.W. Interface Crack Between Two Elastic Layers // Int. J. Fract. - 1990. -Vol. 43. - P. 1-18.
9. Нобл Б. Применение метода Винера-Хопфа для решения дифференциальных уравнений в частных производных. - М.: Изд-во иностр. лит., 1962. - 279 с.
10. Дёч Г. Руководство к практическому применению преобразования Лапласа. - М.: Гос. изд-во физ.-мат. лит. 1958. - 207 с.
11. Khrapkov A.A. Wiener-Hopf method in mixed elasticity theory problems. - Sait Petersburg,
2001.
12. Мусхелишвили Н.И. Сингулярные интегральные уравнения. - 2-е изд. - М.: Наука, 1962. -600 с.
13. Гахов Ф.Д. Краевые задачи. - М.: Наука, 1977. - 640 с.
14. 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.
15. Cotterell B., Chen Z. Buckling and cracking of thin film on compliant substrates under compression // Int. J. Fract. - 2000. - Vol. 104. - No. 2. - P. 169-179.
16. 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. Cimetière, J. Grilhé // Acta materialia. - 2005. - Vol. 53. - P. 441-447.
17. Гольдштейн Р.В., Устинов К.Б., Ченцов А.В. Оценка влияния податливости подложки на напряжения потери устойчивости отслоившегося покрытия // Вычисл. мех. сплошных сред. -2011. - Т. 4, № 3. - С. 48-57.
18. Ustinov K.B. On influence of substrate compliance on delamination and buckling of coatings // Eng. Fail. Analys. - 2015. - Vol. 48B - P. 338-344.
19. Устинов К.Б., Каспарова Е.А. Оценка влияния кривизны и податливости основания на параметры отслоения покрытия // Деформация и разрушение материалов. - 2015. - № 3. - С. 28-35.
20. Салганик Р. Л. О хрупком разрушении склеенных тел // ПММ. - 1963. - Т. 27, № 5. -C. 957-962.
21. Malyshev B.M., Salganik R.L. The strength of adhesive joints using the theory of crack // Int. J. Fracture Mechanics. - 1965. - Vol. 1. - No. 2. - P. 114-128.
22. Hutchinson J.W., Suo Z. Mixed Mode Cracking in Layered Materials // Advances in Applied Mechanics. Eds. J.W. Hutchinson, T.Y. Wu. - 1992. - Vol. 29. - P. 63-191.
23. Златин А.Н., Храпков A.A. Полубесконечная трещина, параллельная границе упругой полуплоскости // Докл. АН СССР. - 1986. - Т. 31. - С. 810-813.
24. Златин А.Н., Храпков A.A. Упругая полуплоскость, ослабленная трещиной, параллельной ее границе // ЛГУ. Исследования по упругости и пластичности. - 1990. - T. 16. Проблемы современной механики разрушения. - С. 68-75.
References
1. Ustinov K.B. O rassloenii polosy po granice razdela uprugikh svoistv. Chast' 1. Postanovka zadachi, sluchai normal'nogo otryva [On delamination of a stipe along the boundary between two elastic layers. Part 1. Problem formulation, the case of normal crack] VestnikPNRPUMechanics Bulletin, 2015. no. 4, pp. 226-245.
2. Ustinov K.B., On Shear Separation of a Thin Strip from the Half-Plane. Mech. Solids., 2014, vol. 49, no. 6, pp. 713-724.
3. Koiter W.T. On the diffusion of load from a stiffener into a sheet. Quart. J. Mech. Appl. Math, 1955, vol. 8, pt 2, pp. 164-178.
4. Popov G.Ya. Bending of a semi-infinite plate resting on a linearly deformable foundation J. Appl. Math. Mech., 1961, vol. 25, iss. 2, pp. 502-520.
5. Alblas J.B., Kuypers W.J.J. On the diffusion of load from a stiffener into an infinite wedge-shaped plate. Appl. Sci. Research, Sec. A., 1965-1966, vol. 15, iss. 1, pp. 429-439.
6. Kalandia A.I. Matematicheskie metody dvumernoi uprugosti [Mathematical methods of 2-D elasticity]. Moscow: Nauka. 1973. 303 p.
7. 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.
8. Suo Z., Hutchinson J.W. Interface Crack Between Two Elastic Layers. Int. J. Fract. 1990, vol. 43, pp. 1-18.
9. Noble B., Methods Based on the Wiener-Hopf Technique for the Solution of Partial Differential Equations (Pergamon Press, London etc., 1958; Mir, Moscow, 1962).
10. Doetsch G. Handbook der Laplace-Transformation (Birkhäuser, Basel, 1946; Fizmatlit, Moscow, 1958).
11. Khrapkov A.A. Wiener-Hopf method in mixed elasticity theory problems. Sait Petersburg, 2001.
12. Muskhelishvili N.I. Singular integral equations; boundary problems of function theory and their application to mathematical physics. Groningen: P. Noordhoff, 1953. 600 p.
13. 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.
14. 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.
15. Cotterell B., Chen Z. Buckling and cracking of thin film on compliant substrates under compression. Int. J. Fract., 2000, vol. 104, no. 2, pp. 169-179.
16. Parry G., Colin J., Coupeau C., Foucher F., Cimetière A., Grilhé 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.
17. Goldstein R.V., Ustinov K.B., Chentsov A.V. Otsenka vliianiia podatlivosti podlozhki na napriazheniia poteri ustoichivosti otsloivshegosia pokrytiya [Influence of substrate compliance on stresses generating stability loss of a delaminated coating. Comput. continuum mech]. Vychislitelnaya mechanika sploshnyh sred, 2011, vol. 4, no. 3, pp. 48-57.
18. Ustinov K.B. On influence of substrate compliance on delamination and buckling of coatings. Eng. Fail. Analys, 2015, vol. 48B, pp. 338-344.
19. 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.
20. Salganik R.L. The Brittle Fracture of Cemented Bodies. J. Appl. Math. Mech., 1963, vol. 27, no.5, pp.1468-1478.
21. Malyshev B.M., Salganik R.L. The strength of adhesive joints using the theory of crack. Int. J. Fracture Mechanics, 1965, vol. 1, no. 2, pp. 114-128.
22. Hutchinson J.W., Suo Z. Mixed Mode Cracking in Layered Materials. Advances in Applied Mechanics. Eds. J.W. Hutchinson, T.Y. Wu, 1992, vol. 29, pp. 63-191.
23. 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.
24. Zlatin A.N., Khrapkov A.A. Uprugaia poluploskost, oslablennaia treschinoi, parallelnoi eio granitse [Elastic Half-Plane Weakened by a Crack Parallel to Its Boundary]. Issledovaniya po uprugosti i plastichnosti, vol. 16. Problemy sovremennoi mechaniki razrusheniya. Leningradskii gosudarstvennyi universitet, 1990, pp. 68-75.
