CC BY
14
DOI: 10.5862/JPM.242.14 УДК: 539.3
В.В. Тихомиров
Санкт-Петербургский политехнический университет Петра Великого
ПОЛУБЕСКОНЕЧНАЯ ТРЕЩИНА МОДЫ III В БИМАТЕРИАЛЬНОМ КЛИНЕ
Получено точное решение антиплоской задачи для полубесконечной интерфейсной трещины в кусочно-однородном клине при самоуравновешенной нагрузке на ее берегах. Рассмотрены три типа граничных условий на сторонах клина: обе стороны свободны от напряжений, обе стороны защемлены и одна сторона свободна, а вторая защемлена. В результате применения метода Винера—Хопфа решение представлено в квадратурах. Для коэффициентов интенсивности напряжений получены функции Грина, а в случае клина геометрически симметричной структуры для них найдены простые формулы. Изучена сингулярность напряжений в вершине клина. В отличие от однородной структуры, асимптотика напряжений вблизи вершины при некоторых значениях параметров композита может иметь два сингулярных члена.
ИНТЕРФЕйСНАЯ ТРЕЩИНА МОДЫ III, БИМАТЕРИАЛЬНЫЙ КЛИН, КОЭФФИЦИЕНТ ИНТЕНСИВНОСТИ НАПРЯЖЕНИЙ.
Введение
Исследованию полей напряжений при антиплоской деформации клиновидных областей, содержащих дефекты типа трещин, посвящен ряд публикаций. Помимо самостоятельного значения, интерес к изучению этой проблемы связан, например, с возникновением и распространением трещин в зоне высокой концентрации напряжений вблизи вершины острого выреза. Как правило, в литературе рассматривались трещины конечной длины, которые ориентированы вдоль интерфейса разнородных материалов композита, причем как изотропных, так и анизотропных.
Одной из первых работ в этом направлении явилась статья Ф. Эрдогана и Г. Гупты [1], в которой рассмотрена интерфейсная трещина в композитном клине, получено сингулярное интегральное уравнение и использован численный метод его решения. В работе [2] для случая нескольких колли-неарных трещин, расположенных на оси симметрии изотропного клина, построено точное решение. Анализ трещины, исходящей из вершины анизотропного клина, был проведен А. Шахани [3]. На основе метода
дислокаций в работе [4] для нескольких произвольно ориентированных трещин в изотропном клине получена система сингулярных интегральных уравнений, для решения которой была предложена численная процедура. К. Ву, Ю. Дзенис и В. Зоу [5] исследовали методом комплексных потенциалов равновесие клина, состоящего из двух трансверсально-изотропных пьезоэлектрических материалов, с интерфейсной трещиной конечной длины.
В работе [6] было установлено, что решение антиплоской задачи для композитного анизотропного клина может быть получено из решения соответствующей задачи, в которой материалы являются изотропными, с помощью преобразования полярных координат. По этой причине решения, полученные для изотропных компонент композита, приобретают особую важность.
эффективный способ построения точных решений задач рассматриваемого направления продемонстрирован в работах Х. Беома и Х. Янга [7, 8] для интерфейсной трещины и трещины, лежащей вне зоны интерфейса. Для случая нагружения граней клина сосредоточенными силами авторы свели задачу к скалярному уравнению Ви-
нера — Хопфа. Однако при решении этого уравнения для факторизации функций использовался метод бесконечных произведений, что несколько снижает эффективность подхода.
В настоящей работе исследовано равновесие кусочно-однородного клина с полубесконечной интерфейсной трещиной моды III, на берегах которой приложена самоуравновешенная нагрузка. Рассмотрены три варианта граничных условий на сторонах клина:
1) обе стороны свободны от напряжений;
2) обе стороны защемлены;
3) одна сторона свободна, а вторая защемлена.
В каждом из вариантов с помощью интегрального преобразования Меллина задача сведена к задаче Римана [9]. На основе техники, развитой в работе [10], построено ее точное решение в квадратурах. При этом факторизация функций выполнена с использованием интегралов типа Коши. В случае геометрически симметричной структуры композита получены простые формулы для коэффициента интенсивности напряжений (КИН) в вершине трещины. Исследованы предельные ситуации, когда относительная жесткость композита стремится к нулю или к бесконечности. Изучена сингулярность напряжений в вершине клина.
оказывается, что в отличие от случая однородного материала, при некоторых значениях параметров композиции асимптотика напряжений вблизи этой особой точки будет содержать два сингулярных члена. Данный факт не был отмечен в предыдущих работах рассматриваемой тематики [11, 12].
Постановка задачи и ее сведение к задаче Римана
Рассмотрим полубесконечную интерфейсную трещину моды III между двумя разнородными клиньями с углами а1 и а2 (рис. 1). Материалы клиньев в областях Ц и 0.2 будем считать изотропными, однородными и имеющими модули сдвига | и |2. Контакт материалов вне трещины предпо-
Рис. 1. Биматериальный клин с интерфейсной полубесконечной трещиной:
□2 — области материала с модулями сдвига | и
а,, а.
12;
р б, г — геометрические параметры; g(r) — самоуравновешенная нагрузка, приложенная к берегам трещины
лагается идеальным. Без ущерба общности будем считать, что вершина трещины находится на расстоянии единичной длины от вершины композитного клина. К берегам трещины приложена самоуравновешенная нагрузка g (г) (г — полярный радиус).
Как известно, перемещения ж в областях Ц и 0.2 в рассматриваемом случае являются гармоническими функциями:
д Ж
1 д ж . 1 дж.
г2 дб2
. + __
г дг
= 0 (7=1,2), (1)
а напряжения в полярных координатах г и 0 определяются формулами
тбд ='
I- Ж
дб
дж -
т • = |. —-
ГЦ дг
(2)
На линии б = 0 напряжения и перемещения удовлетворяют условиям
тб1(г, +0) = т0г2 (г, -0),
ж1(г, +0) = ж2(г, -0) (0 < г < 1) (3)
тб1(г, +0) = тб2(г, -0) = g(r) (1 < г < да).
Рассмотрим три варианта граничных условий на сторонах клина:
(г, а1) = тб1(г, -а2) = 0 (задача1); (4)
м^(г, а1) = ^2(г, -а2) = 0 (задача 2); (5)
т0г1(г, а1) = 0, *2(г, -а2) = 0 (задача 3). (6)
Решение каждой задачи ищем в виде интегралов Меллина (п — номер задачи):
wjn (r, 0) = 2П7 W (P'0)r"Pdp
(j = 1,2; n = 1,2,3), j (r' 0) = ^ Í Tejn (P' 0)r-P-1dp,
(7)
V (r, 9) = ¿J J Tj (P' "
где трансформанты перемещений и напряжений в соответствии с равенствами (l) и (2) определяются следующими формулами:
Wjn (P' 0) = AM (p) sin p0 + B M (p) cos p0; (8) Tj(P,0) = ЦjP[Ajn(p)cosp0 - Bjn(p)sinp9];
Tzn (p, 0) = -M Ajn (p) sin p0 + Bjn (p) cos p0].
Согласно условиям регулярности, контур интегрирования L расположен параллельно мнимой оси в полосе
-s < Rep <8 (s, 8 > 0).
Величины Ajn (p) и Bjn (p) находятся из равенств (3) и соответствующих рассматриваемой задаче условий (4) — (6).
Из условий (3), с учетом представлений (7) и (8) имеем:
MipAn p) = Ц2 pA2n( p) = Tn+( p) + G-( p),
-p[Bin( p) -B2n( p)] =^n- (p),
где
Tn+ (p) = J Vin (r, +0)rpdr,
0
да
G- (p) = J g (r )rpdr,
i
Un-(p) = J ^ [Win(r, +0) - W2n(r, -0)]rpdr.
l dr
функции G (p) и Un- (p) регулярны и не имеют нулей в левой от контура L полуплоскости Q-, а Tn+ (p) — в правой полуплоскости Q+ [13].
Используя граничные условия на сторо-
нах клина (4) — (6), приходим к неоднородным скалярным задачам Римана [9]:
ЗДТП+Cp) + G-(p)] - q^UJp) = 0
(9)
(p e L; n = 1,2,3).
Здесь q = (-1)n, в качестве контура L может быть взята мнимая ось, а функции Fn (p) имеют вид
F1(p) = ц ctg pa1 + ctg pa2;
F2( p) = ц tg pa: + tg pa2; (10)
F3(p) = ц ctg pa1 - tg pa2.
Входящая в эти формулы постоянная ц = ц2/ц является относительной жесткостью композита, которая может изменяться в пределах от нуля до бесконечности.
Решение задач Римана
Ключевым моментом при решении задач (9) является факторизация функций. Представим их коэффициенты (10) в виде
Fn (p) = (1 + ц) pqMn (p), (11)
Mn(p) = X(Р)Фп(p), X(p) = p ctg p«1,
где
Ф1p) = (1 + ц)-1(ц + tg pa ctg pa2);
®2( p) = (1 + Ц)-1(Ц + ctg pa1 tg pa2); (12)
Фз(p) = (1 + ц)-1(ц - tg pax tg pa2).
Факторизация функции X(p) осуществляется элементарно [13]:
x+ (p)
x (p) =
x± (p) =
x- (p)' п Г(1 ± pa^n)
(13)
[а1 Г(1/2 ± ра^п)'
где Г(х) — гамма-функция.
На мнимой оси при р = И функции (12) непрерывны, не имеют нулей и полюсов, их индекс равен нулю, и при t ^ да экспоненциально стремятся к единице. Поэтому справедливы следующие равенства [10, 13]:
Фп+ (р)
фп (p) =
Фп- ( p)'
Фп± (p) = exp
_!_ г dt
2п/ i t - p
(p й £)■
С учетом четности функций (12) аналитические функции в областях и могут быть представлены в форме
р 11п Фя (/
фя± (Р) = exp
(15)
п 0 + р2
Таким образом,
Ип(р) = Мп+( р)Мп-( р),
Мп± (р) = Х±9 (р )ФИ± (р).
В результате использования формул (11), (13) — (15), перегруппировки слагаемых в уравнениях (9) и применения теоремы Лиувилля [13], получаем:
Мп+ (р)Т+ (р) + бп+ (р) = 912(1 + I)-1 р~9 ^ X Мп (р)^ (р) - 0п-(р) = /п(р),
(16)
где
а± (p) =
1- j
М)
2п/ ¿ t - p
dt,
e„(t) = F,(t g (t )•
(17)
(1 + Ц) tq
Из оценок членов в равенстве (16) при p ^ ж вытекает, что единые аналитические функции Jn(p) = const = Cn, причем
Q = 0.
Константы С1 и С3 определим из уравнения (16) при значении р = 0. Учитывая формулы (17), находим:
Сп = Мп+ (0)в_ (0) - 1
2 п/(1 + ц)
xjMn_(t)Fn(t)G (t)dt (n = 1,3).
(18)
(20)
Запишем функции (10) в задачах 1 и 3 в виде
Fn(t) = (1 + ц)Дп(0//п(0; (19) A1(t) = sin(a1 + a2)t + m sin(a1 - a2)t;
Д3(t) = cos(a1 + a2)t - m cos(a1 - a2)t; f1 (t) = 2 sin t a1 sin t a2; f2(t) = 2 sin ta1 cos ta2.
В формулы (20) входит биупругая постоянная m = (1 - ц) / (1 + ц). Для всех сочетаний модулей сдвига материалов эта величина удовлетворяет неравенству Iml <1,
причем значение т = 0 отвечает однородной среде.
Из представлений (19) вытекает, что подынтегральные функции в формулах (18) имеют в полуплоскости полюсы в точках = -пк/а и ¿/2 = -пк/а2 (к = 1,2,...) в задаче 1, а в задаче 3 — в точках и зк2 = - п(2к - 1)/(2а2). Тогда, используя теорему о вычетах, находим:
С1 = ^^Чг^(0) - -+- х
-^Х M1- (tkl)G_ (tkl)-a1 k=1
1
M1- (tk 2 )G_ (tk 2)
a2 k=1
-ц-G (0) - —x
a1(l + ц) 1 + ц
(21)
I I _
-^Х M3- (tk1)G- (tk 1)-
a1 k=1
Х M3- (sk 2 )G_ (Sk 2)
a2 k=1
Из уравнения (16) для п = 1, 3 получаем:
Тп+ (р) = [Сп - бп+ (р)]Х+-1(р)Ф-+ (р). (22)
При р ^ да, согласно формулам (13), (14) и (17), имеют место оценки
Х+ (р) = 0(р2), Фп±(р) = 0(1),
6п+ (р) = 0(р1).
следовательно, представление (22) приводит к асимптотике
Тп+ (р) ~ Спр-12 (р ^ да).
Отсюда по теореме абелева типа [13] заключаем, что асимптотика напряжений при г ^ 1 - 0 имеет вид
С
,0)
(п = 1,3). (23)
л/п(1 - г)
В задаче 2 из уравнения (16) находим 72+(р) = -Х+(р)Ф-+(р)02+(р).
Отсюда получаем, что 72+ (р) ~ -Вр 12 при р ^ да, где
B =
1
i C (t )Л,
2ni(1 + ц) L t X (t) F2(t) = (1 + ^(t V/2(t), (24) A2(t) = sin(a1 + a2)t - m sin(a1 - a2)t, (25) /2(t) = 2 cos ta1 cos ta2.
Тогда, аналогично предыдущим случаям, в задаче 2 будем иметь:
B
8z 21
(г,0) ~ --
(r ^ 1 - 0). (26)
- г)
В результате применения теоремы о вычетах константу В можно представить в виде
B =
2
2 » р^ф-^^ (27) п(1 + ц) £ £ (2к - 1)X- () '
01 = ц Р2 = 1> = -п(2к - 1)/(2а ;) (к = 1,2,...).
Коэффициенты интенсивности напряжений
Определим коэффициент интенсивности напряжений (КИН) в вершине трещины в каждой задаче формулой
= Шп^т^ (г ,0).
г ^1-0
Тогда, используя асимптотики (23) и (26), получаем:
К^а,, а2, ц) = ЛСп (п = 1,3),
(28)
К 111 (а1, а2, ц) = -42.В.
в случае геометрически симметричного клина формулы для КИН в задачах 1 и 2 существенно упрощаются. При а1 = а2 = а из формул (12) вытекает, что ФДр) = Ф2(р) = 1 и, следовательно, согласно (14), Фп± (р) = 1 (п = 1,2). Учитывая равенства (13) , (15), (21) и (27), в итоге из (28) будем иметь:
К<»(а, а, ц) = К<».(а) =
е-(0) + ЕГ--
a
K2)(a, a, ц) = K^(a) =
(29)
= 2 £ (2k_-1)_Н g(_JL (2 k +1)
Vak=0 2kk! 4 2a '
Отсюда вытекает, что в случае клина
геометрически симметричного строения с одинаковыми граничными условиями на его гранях КИН в вершине трещины не зависит от относительной жесткости материалов и совпадает с КИН для однородного клина.
Если к берегам трещины приложены сосредоточенные силы Т0 на расстоянии г0 > 1 от вершины клина, т. е. g(r) =-Т05(г - г) (5(г-г0) — дельта-функция), то б- (р) = -Т0г0р. Тогда ряды в (29) суммируются, и для КИН получаем простые представления
к« = т I2 Го
к(2) = т
■/(2a )
VroV° - 1 ' 2 1 Vr,Va-1.
,/(2a )
(30)
Поскольку КЩ„/К^» = г0" а > 1 при г0 > 1 и 0 < а < п, то КИН в вершине трещины для клина со свободными гранями всегда превышает его значение для случая, когда грани клина защемлены.
В частности, при п = 1 и а = п получаем задачу о взаимодействии двух коллине-арных интерфейсных трещин, к одной из которых приложена точечная нагрузка. В этом случае
= Т^2г0/[п(Г) - 1)].
В задаче 2 при а = п имеем случай взаимодействия полубесконечной интерфейсной трещины с полубесконечным жестким включением. Для КИН в вершине трещины согласно (30) получаем
К% = Т^2/[п(г0 - 1)].
Возвращаясь к общему случаю а ^ а2, найдем асимптотики КИН для предельных значений относительной жесткости композита. Устремляя ц к бесконечности, используя формулы (21), (27) и (28), а также обозначения (29), получаем:
К ^ К«»(а), К(2 ^ К(2?.(а1>,
(31)
К ^К«»(а) (ц^да).
Таким образом, в каждой из задач асимптотики КИН определяются соответствующими выражениями для случая геометри-
На рис. 2 представлены зависимости нормализованных КИН в вершине трещины, нагруженной сосредоточенными силами, для а1 = п/ 2, а2 = 3 п/ 2 и г0 = 4 от параметра относительной жесткости клина. При таком выборе углов задача 1 соответствует случаю взаимодействия двух ортогональных трещин, задача 2 — трещины и жесткого включения, а задача 3 — трещины и отслоившегося жесткого включения. Приведенные данные показывают, что неоднородность композита может вызывать как эффект увеличения КИН по сравнению со случаем однородной среды, так и эффект его уменьшения. Наиболее сильно эти эффекты проявляются при малых и больших значениях относительной жесткости в задаче 3.
Поля напряжений вблизи вершины клина
Исследуем поля напряжений в вершине клина при г ^ 0. На линии б = 0 касательные напряжения на основании формул (7) и условий (3) имеют вид
№
0.6
0.5 -\-1-1-1-1-1-1-1-1-
0 2 4 6 8 ц
Рис. 2. Зависимости нормализованных КИН от модуля сдвига ц в вершине трещины, расположенной в композитной плоскости при а: = п/2, а2 = 3п/2 и г0 = 4; номера кривых
соответствуют номерам задач
чески симметричного однородного клина с углом раствора, равным 2 а:. При этом для больших значений относительной жесткости предельные значения КИН в задачах 1 и 3 совпадают.
Когда параметр ц близок к нулю, будем иметь:
К ^ К«,^), КЦ ^ К®^),
(32)
Кщ ^ К^Ю (| ^ 0).
Отсюда вытекает, что и в случае малой относительной жесткости композита предельные значения КИН выражаются формулами (29) для однородного клина при а = а2. Заметим, что асимптотики (31) и (32) также могут быть получены путем решения предельных для (9) задач, возникающих при | ^ да и | ^ 0.
Введем в рассмотрение нормализованные КИН:
N(п) = К«(а1, а2,|)/ К«(а1, а2,1),
отражающие влияние на КИН неоднородности композиционной структуры.
хб- (г, 0) = ¿7 {{Тп+ (р) + (р)]г"-Ыр.(33)
С помощью равенств (15) и (16) выражение для подынтегральной функции можно записать следующим образом:
Тп+ (р) + О- (р) = (1 + ц)р" М + {?\
мп -(р)тп(р)
(п = 1,2,3).
Тогда с учетом представлений (19) и (24) формула (33) приобретает вид
^(г -0) = Ъ [ гп( р) г'р^
где У, (р) = [Сп + (р)]р'М;-( р).
В результате применения теоремы о вычетах в полюсах, определяемых отрицательными корнями -р(кп) уравнений (см. (20 ) и (25 ))
Дп (р, ап а2, т) = 0, (34)
получаем:
да (п)
Тб^п (г, 0) = X Ак (аl, а2, Ц)г
к=0
(Г« = 1 - р«).
(35)
Коэффициенты этого ряда находятся по формулам
Ак(а1, а2, ц) = Уп(-р<п))/п(рГУДп(рГ),
где Дп(р) = d Дп( р)/ ^р.
Из (35) вытекает, что сингулярные поля напряжений в окрестности вершины клина порождаются корнями уравнений (34), расположенными в полосе 0 < Ке р(кп) < 1. По аналогии с выкладками, представленными в работе [10], можно показать, что все числа р(п) в этой полосе являются вещественными и однократными. Функции (20) и (25) обладают следующими свойствами:
Д1 (р, а1, а2, т) = Д1 (р, а2, а1, -т); (36)
Д2 (р, а1, а2, т) = Д1 (р, а1, а2, -т); (37)
Д3(р, а1, а2, т) = Д3( р, а2, а1, т). (38)
Из равенств (36) и (38) следует, что в задачах 1 и 3 достаточно рассмотреть, например, корни уравнений (34) при а1 < а2 для -1 < т < 1. В силу свойства (37) корни в задаче 2 совпадают с корнями в задаче 1
при замене параметра т на —т.
Не касаясь в данной статье проведенного детального анализа корней уравнений
(34), отметим только некоторые свойства корней, расположенных в интервале (0, 1). В зависимости от параметров композиционной структуры уравнения (34) могут не иметь корней, иметь один корень или два корня в данном интервале. В вершине клина эти корни могут порождать как слабые (0 < Х[п) < 0,5), так и сильные (0,5 < Х<п) < 1) сингулярности. При наличии двух корней в отмеченном интервале могут быть только две ситуации: 0,5 < р1п) < р2п) < 1 или 0 < р! < 0,5 < рр] < 1. Иными словами, такие распределения корней внутри интервала определяют либо две слабые сингулярности поля напряжений при г ^ 0, либо одну сильную и одну слабую сингулярности. существование двух сильных сингулярностей в асимптотике напряжений невозможно.
В задачах 1 и 2 отсутствие сингулярных членов в представлениях напряжений
(35) имеет место, например, при любом т е (-1 , 1) для углов раствора клина, не превышающих п/2. В задаче 3 подобная ситуация для углов того же диапазона возникает при т < 0.
Для геометрически симметричных композитных клиньев (а1 = а2 = а) в задачах 1 и 2 уравнения (34) имеют в интервале (0, 1) только один корень, если а > п/2. Этот корень не зависит от биупругой постоянной и совпадает с корнем для однородного клина р1(п) = п/(2а). В задаче 3 в этом случае возможны три варианта. Во-первых, при углах а < а» = 0,5 агееоз т корни, меньшие единицы, отсутствуют. Во-вторых, для углов, удовлетворяющих неравенству а» < а < п - а», в интересующем нас интервале существует один корень. Наконец, для углов п - а» < а < п напряжения в угловой точке клина будут содержать два сингулярных слагаемых, порождаемых корнями
р13) = агееоз ти/(2а),
р23) = ^а - агееоз т/(2а).
На рис. 3, а приведены зависимости первых двух корней уравнения (34) в задаче 1 от угла а1 при различных значениях относи-
тельной жесткости, когда композитная среда занимает всю плоскость (а1 + а2 = 2п) и содержит две интерфейсные трещины. При
а)
0 < ц < 1 (биупругая постоянная т > 0) уравнение (34) имеет в интервале (0, 1) два корня: 0 < р! < 0,5 < < 1, если
Рис. 3. Зависимости первых двух корней уравнений (34) от угла а: при а: + а2 = 2п в задаче 1 (а)
и задаче 3 (Ь) для различных значений относительной жесткости 2 (номера кривых 1 и 1'), 10 (2 и 2'), 100 (3 и 3'), 0,5 (4 и 4'), 0,1 (5 и 5'), 0,01 (6 и 6') (кривые, соответствующих второму
корню, отмечены штрихом)
0 < а1 < п/2, а при п/2 < а1 < п — один корень 0 < р^1' <0,5. В случае ц > 1 (т < 0) во всем диапазоне изменения угла а1 корни превосходят 0,5. Это значит, что асимптотика напряжений в вершине трещины, совпадающей с угловой точкой линии раздела материалов, может содержать один или два сингулярных члена. При этом показатель сингулярности отличен от классического значения, равного 0,5. Для двух коллинеар-ных трещин (а = п) классический результат А-11' = 0,5 сохраняется.
Задача 3 при а1 + а2 = 2п представляет задачу для составной упругой плоскости, содержащей полубесконечную трещину и отслоившееся полубесконечное жесткое включение. В этом случае при любых значениях относительной жесткости композита ц е (0, да) во всем диапазоне изменения угла а1 уравнение (34) имеет два корня (см. рис. 3, Ь), порождающие сильную и слабую особенности у напряжений в вершине клина.
Заключение
В данной работе с помощью метода Винера — Хопфа получено точное решение задачи о равновесии биматериального клина, ослабленного полубесконечной межфазной трещиной продольного сдвига, при трех типах граничных условий на его сторонах. В случае композита геометрически симметричной структуры и приложения к берегам трещины сосредоточенных сил, найдены простые алгебраические представления функций Грина для коэффициентов интенсивности напряжений. Исследована сингулярность напряжений в вершине клина. В отличие от однородного клина, асимптотика напряжений в вершине клина при некоторых значениях параметров неоднородной структуры может содержать два сингулярных члена.
Данный подход позволяет получить точные решения как для случая приложения нагрузки на сторонах клина, так и случая клина с надрезом, исходящим из его вершины.
список литературы
[1] Erdogan F., Gupta G.D. Bonded wedges with an interface crack under anti-plane shear loading // Int. J. Fract. 1975. Vol. 11, No. 4. Pp. 583-593.
[2] Mkhitaryan S.M., Melkoumian N., Lin B.B. Stress-strain state of a cracked elastic wedge under anti-plane deformation with mixed boundary conditions on its faces // Int. J. Fract. 2001. Vol. 108. No. 4. Pp. 291-315.
[3] Shahani A.R. Mode III stress intensity factors for edge-cracked circular shafts, bonded wedges, bonded half planes and DCB's // Int. J. Solids Struct. 2003. Vol. 40. No. 24. Pp. 6567-6576.
[4] Faal R.T., Fotuhi A.R., Fariborz S.J., Daghyani H.R. Antiplane stress analysis of an isotropic wedge with multiple cracks // Int. J. Solids Struct. 2004. Vol. 41. No. 16-17. Pp. 4535-4550.
[5] Wu X., Dzenis Y., Zou W. Screw dislocation interacting with an interfacial edge crack between two bonded piezoelectric wedges // Int . J. Fract. 2002. Vol. 117. No. 3. Pp. 9-14.
[6] Beom H.G., Lee J.W., Cui C.B. Analysis of a kinked crack in an anisotropic material under antiplane deformation // J. Mech. Sci. Technol. 2012. Vol. 26. No. 2. Pp. 411-416.
[7] Beom H.G., Jang H.S. Interfacial wedge cracks in dissimilar anisotropic materials under antiplane shear // Int. J. Eng. Sci. 2012. Vol. 56. No. 1. Pp. 49-62.
[8] Beom H.G., Jang H.S. A crack emanating from a wedge in dissimilar anisotropic materials under antiplane shear // Int. J. Fract. 2012. Vol. 177. No. 2. Pp. 191-198.
[9] Gakhov F.D. Boundary Value Problems. London: Pergamon Press, 1966.
[10] Nuller B., Ryvkin M., Chudnovsky A.A. Closed form solution for a crack approaching an interface // J. Mech. Mat. Struct. 2006. Vol. 1. No. 8. Pp. 1405-1423.
[11] Ma C.C., Hour B.L. Analysis of dissimilar anisotropic wedges subjected to antiplane shear deformation // Int. J. Solids Struct. 1989. Vol. 25. No. 11. Pp. 1295-1309.
[12] Shahani A.R. Mode III stress intensity factors in an interfacial crack in dissimilar bounded materials // Arch. Appl. Mech. 2006.Vol. 75. No. 6-7, Pp. 405-411.
[13] Noble B. Method based on the WienerHopf technique for solution of partial differential equations. Oxford: Pergamon Press, 1958.
сведения об авторе
тИхоМИРоВ Виктор Васильевич — кандидат физико-математических наук, заместитель директора по образовательной деятельности Санкт-Петербургского политехнического университета Петра Великого.
195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29 уюйкЬ@таП.ги
Tikhomirov V.V. A SEMI-INFINITE CRACK OF MODE III IN THE BIMATERIAL WEDGE.
An exact solution of the antiplane problem for a semi-infinite interface crack in a piecewise-homogeneous wedge under a self-balanced load on its sides has been obtained. Three types of boundary conditions on the wedge sides were examined: the both sides being stress-free; both sides being clamped, and one side being stress-free with the second one clamped. As a result of using the Wiener-Hopf method, the solution was represented in quadratures. The Green's functions were obtained for stress intensity factors; in the case of a geometrically symmetrical wedge structure simple formulae were found for these functions. The stress singularity in the apex of the wedge was studied. In contrast to the homogeneous wedge structure the asymptotic of the stresses near the apex was established to have sometimes two singular terms for some values of the composite parameters.
INTERFACIAL MODE III CRACK, BIMATERIAL WEDGE, STRESS INTENSITY FACTOR, STRESS SINGULARITY.
REFERENCES
[1] F. Erdogan, G.D. Gupta, Bonded wedges with an interface crack under anti-plane shear loading, Int. J. Fract. 11 (4) (1975) 583-593.
[2] S.M. Mkhitaryan, N. Melkoumian, B.B.
Lin, Stress-strain state of a cracked elastic wedge under anti-plane deformation with mixed boundary conditions on its faces, Int. J. Fract. 108 (4) (2001) 291-315.
[3] A.R. Shahani, Mode III stress intensity factors for edge-cracked circular shafts, bonded wedges, bonded half planes and DCB's, Int. J. Solids Struct. 40 (24) (2003) 6567-6576.
[4] R.T. Faal, A.R. Fotuhi, S.J. Fariborz, H.R. Daghyani, Antiplane stress analysis of an isotropic wedge with multiple cracks, Int. J. Solids Struct. 41(16-17) (2004) 4535-4550.
[5] X. Wu, Y. Dzenis, W. Zou, Screw dislocation interacting with an interfacial edge crack between two bonded piezoelectric wedges, Int. J. Fract. 117(3) (2002) 9-14.
[6] H.G. Beom, J.W. Lee, C.B. Cui, Analysis of a kinked crack in an anisotropic material under antiplane deformation, J. Mech. Sci. Technol. 26(2) (2012) 411-416.
[7] H.G. Beom, H.S. Jang, Interfacial wedge cracks in dissimilar anisotropic materials under antiplane shear, Int. J. Eng. Sci. 56 (1) (2012) 49-62.
[8] H.G. Beom, H.S. Jang, A crack emanating from a wedge in dissimilar anisotropic materials under antiplane shear, Int. J. Fract. 177 (2) (2012) 191-198.
[9] F.D. Gakhov, Boundary Value Problems, London, Pergamon Press, 1966.
[10] B. Nuller, M. Ryvkin, A.A. Chudnovsky,
Closed form solution for a crack approaching an interface, J. Mech. Mat. Struct. 1(8) (2006) 1405-1423.
[11] C.C. Ma, B.L. Hour, Analysis of dissimilar anisotropic wedges subjected to antiplane shear deformation, Int. J. Solids Struct. 25(11) (1989) 1295 -1309.
[12] A.R. Shahani, Mode III stress intensity factors in an interfacial crack in dissimilar bounded materials, Arch. Appl. Mech. 75 (6-7) (2006) 405-411.
[13] B. Noble, Method based on the WienerHopf technique for solution of partial differential equations, Oxford, Pergamon Press, 1958.
THE AuTHOR
TIKHOMIROV Victor V.
Peter the Great St. Petersburg Polytechnic University
29 Politechnicheskaya St., St. Petersburg, 195251, Russian Federation
victikh@mail.ru
© Санкт-Петербургский политехнический университет Петра Великого, 2016
