Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2021. Т. 25, № 1. С. 67—82 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1816
Механика деформируемого твёрдого тела
УДК 517.958:531-133
Плоско-деформированное состояние равномерно кусочно-однородной плоскости с периодической системой полубесконечных межфазных трещин
© В. Н. Акопян, А. А. Григорян
Институт механики НАН Республики Армения,
Республика Армения, 0019, Ереван, пр. Маршала Баграмяна, 24Б.
Аннотация
Рассмотрено плоско-деформированное состояние равномерно кусочно-однородной плоскости, полученной при помощи поочередного соединения двух разнородных полос, которая по линиям стыков разнородных полос расслаблена периодической системой двух полубесконечных межфазных трещин и деформируется под воздействием нормальных нагрузок, приложенных к берегам трещин. Выделена базовая ячейка задачи в виде двухкомпонентной полосы и при помощи обобщенного преобразования Фурье получена определяющая система уравнений задачи в виде одного сингулярного интегрального уравнения второго рода относительно комплексной комбинации контактных напряжений в зоне стыка полос.
В частном случае путем устремления ширины полос к бесконечности получено определяющее уравнение задачи для двухкомпонентной плоскости из двух разнородных полуплоскостей с двумя полубесконечными межфазными трещинами и построено его точное решение. Получено также определяющее уравнение поставленной задачи в виде одного сингулярного интегрального уравнения первого рода относительно нормальных контактных напряжений еще в одном частном случае, когда все полосы изготовлены из одного и того же материала, т.е. в случае однородной плоскости, расслабленной периодической системой параллельных полубесконечных трещин.
Научная статья
Q ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Акопян В. Н., Григорян А. А. Плоско-деформированное состояние равномерно кусочно-однородной плоскости с периодической системой полубесконечных межфазных трещин // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2021. Т. 25, № 1. С. 67-82. https://doi.org/10.14498/vsgtu1816. Сведения об авторах
Ваграм Наслетникович Акопян& https://orcid.org/0000-0003-3684-9471
доктор физико-математических наук, профессор; главный научный сотрудник; отд. теории
упругости и вязкоупругости; e-mail: [email protected]
Арам Арутюнович Григорян © https://orcid.org/0000-0001-7582-1960 аспирант; младший научный сотрудник; отд. теории упругости и вязкоупругости; e-mail: [email protected]
© Самарский государственный технический университет
67
В общем же случае определено поведение искомой функции в концевых точках интервала интегрирования и решение задачи численно-аналитическим методом механических квадратур сведено к решению системы алгебраических уравнений. Получены простые формулы для определения коэффициентов интенсивности напряжений, ^интеграла Черепанова-Райса и раскрытия трещин. Проведен численный расчет. Выявлены закономерности изменения контактных напряжений и интеграла Черепанова-Райса в концевых точках трещин в зависимости от упругих характеристик разнородных полос и геометрических параметров задачи.
Ключевые слова: периодическая задача, смешанная задача, кусочно-однородная плоскость, межфазные трещины.
Получение: 10 августа 2020 г. / Исправление: 29 января 2021 г. / Принятие: 12 февраля 2021 г. / Публикация онлайн: 24 февраля 2021 г.
Введение. Определение локальных полей напряжений вокруг одинарно расположенных или периодически повторяющихся дефектов в однородных и составных массивных телах всегда было и остается одним из приоритетных направлений развития механики разрушения и смешанных задач теории упругости. Изучению этой проблемы и решению различных задач для упругих однородных и составных массивных тел с различными конечными и полубесконечными концентраторами напряжений типа трещин, штампов и включений посвящено много работ. Полученные в них многие основополагающие результаты приведены в монографиях [1-7] и в избранных трудах Г.Я. Попова [8]. Для кусочно-однородных равномерно слоистых тел с периодическими и двоякопериодическими внутренними или межфазными дефектами аналогичные исследования начали проводиться совсем недавно.
В работах [9-11] поставлены и методами разрывных решений и сингулярных интегральных уравнений построены замкнутые или эффективные решения ряда периодических и двоякопериодических антиплоских и плоских задач для кусочно-однородного, равномерно слоистого пространства с межфазными дефектами различного типа. Особо отметим работу [9], которая более тесно связана с рассматриваемой здесь задачей. Здесь построены разрывные решения уравнений теории упругости для кусочно-однородной равномерно слоистой плоскости с межфазными двоякопериодическими конечными дефектами и получены решения двух конкретных задач, когда дефекты представляют собой трещину и абсолютно жесткое включение. Укажем также на работу [11], где построены замкнутые решения двух антиплоских задач для кусочно-однородного равномерно слоистого пространства с периодической системой межфазных туннельных полубесконечных трещин.
1. Постановка задачи и вывод определяющих уравнений. Пусть кусочно-однородная плоскость, изготовленная при помощи поочередного соединения двух разнородных полос равной толщины 2к с коэффициентами Ламэ А1 и ^2, и находящаяся в плоско-деформированном состоянии, на линиях соединения полос у = 2пк, п € Ъ, по бесконечным отрезкам Ь = (-те, —а) и (а, ж) расслаблена периодической системой параллельных полубесконечных трещин (рис. 1).
Рис. 1. Кусочно-однородная плоскость [Figure 1. The piecewise-homogeneous plane with a periodic system of parallel semi-infinite cracks]
Будем полагать, что плоскость деформируется под воздействием одинаковых, самоуравновешенных распределенных нагрузок р(х), действующих на берегах трещин и имеющих конечную результирующую Р.
Необходимо построить решение поставленной задачи и изучить закономерности изменения контактных напряжений и коэффициентов интенсивности напряжений в концевых точках трещин в зависимости от соотношения упругих постоянных и геометрических характеристик разнородных полос.
При такой постановке задачи линии у = (2п + 1)Л,, п € Ъ, являются линиями симметрии, вследствие чего напряженное состояние в составных полосах, находящихся между линиями симметрии у = (2п — 1)Л, и у = (2п + 1)Л,, будет одинаковым. Это означает, что для определения напряженно-деформированного состояния кусочно-однородной плоскости достаточно рассмотреть напряженно-деформированное состояние только двухкомпонентной полосы (базовой ячейки) между линиями симметрии у = ±Л, (рис. 2).
Поставленная задача математически эквивалентна следующей граничной задаче для базовой ячейки:
V] (х, (-1У+1к) = 0, -те <х < те; 2 = 1,2 т^ (ж, (-1У+1К) =0, -те < ж < те; ] = 1, 2
аУ\х, 0) = -'р(х),
у
т^(х, 0) =0,
тг
(1) ху
(1) У
(х, 0)
= т(2)
ху
(х, 0),
аХ' (х, 0) = а(у2)(х, 0),
х
М > а; ] = 1, 2;
|ж| > а; 2 = 1, 2; |ж| < а; |ж| < а; |ж| < а; Ы < а.
(1)
^(ж, 0) = У2(х, 0), „и1(х, 0) = и2(х, 0),
Здесь из (х,у) и V^ (х,у), ] = 1, 2 — соответственно горизонтальные и вертикальные составляющие смещений точек разнородных полос, удовлетворяющие, каждая в своей области определения, уравнениям Ламэ и связанные
с компонентами нормальных ау) (х,у) и касательных тХ3у(х,у) напряжений известными соотношениями [2].
Для решения поставленной задачи мысленно разделим базовую ячейку на две однородные полосы и введем в рассмотрение неизвестные нормальные и касательные контактные напряжения д(х) и т(ж), действующие в зонах контакта разнородных полос:
'тХ3у)(х, 0) = т(х), |ж| < а; ] = 1, 2; аЦ\х, 0) = д(х),
|ж| < а; ] = 1, 2.
(2)
Далее решим вспомогательные задачи для каждой из разнородных полос, занимающих соответственно области {-те < х < те;0 ^ у ^ К} и {-те < х < те; -К ^ у ^ 0}, когда на лицевой линии первой полосы у = К и на лицевой линии второй полосы у = -К заданы условия симметрии, а на линии у = 0 заданы напряжения, и определим смещения точек зоны контакта обеих полос через введенные неизвестные контактные напряжения д(х) и т(ж), т. е. решим граничную задачу (1), заменив последние четыре условия условиями (2) при каждом из значений ] = 1, 2.
С этой целью, следуя [7], решения уравнений Ламэ представим в виде интегралов Фурье:
из (х,У) = ^ Ъ (х,У) = ^
—те
гоо
[(А, + 1уВ**) сЪ(зу) + (В, + гуА**) ъЬ(8у)]е—г°Х<18, [(С, - у А**) сЦзу) + (^ - у В*) вЦ^)] е—г°хй8,
(3)
где Aj, Bj, С^ и И] — неизвестные постоянные, подлежащие определению, а коэффициенты А* и В* определяются формулами
е е
А** = ТГ - ^), В* = - (С, - Ш,-);
к,
к
к
+
Ху + ^
3 = 1, 2.
При этом напряжения записываются следующим образом:
оо
Vy)(%, У) = Г
к
8 Б* + в А* — у В*] сЪ^у) + + С* + эВ* — У А*]
}(*Ь
Ч8,
Т~ху(х, У) = 1
В* — 8С* + ¡*гзуА*] сЪ(зу) + + з А* — зБ* + ¡*гз у В**]
Чв; (4)
4Л =
$(*) = (х* + 2^*
3 = 1, 2.
X* + 3ц*' 2 X* + 3ц*
Используя представления (3) и (4), удовлетворим условиям указанных вспомогательных граничных задач и выразим неизвестные коэффициенты, входящие в эти представления, через трансформанты Фурье функций О(х) и Т(х). Для постоянных А* и С* получим следующие выражения:
А1 = —
С1 = —
Х10
(1)
2г 8 ¡1
сЪР 8ЪР — цР/к^ сЪР 8Ъ Р + Р
ОМ- ^ сЪ Р Т(а)
°($) 28^1 (сЪР +Р)Т (8) (1) -сър 8ЪР — ¡ф/к^
28ц2(сЬР 8ИР + Р)
2г 8¡¡1
сЪР 8Ъ Р + Р
П*),
А =
к 2$
(2)
2г в ¡2
сЪР 8ИР — ^2Р/к2$1]
сЬР 8Ъ Р + Р
^ сЪ2Р ,, ад,
С2 =
О(*) +
к2$
(2)
2г в ¡2
28(сЪР 8ЪР + Р) сЪР 8ЪР — ¡2Р/к2$1)
сЪР 8ЪР + Р
Остальные же постоянные даются формулами
г$(*) „ 1 ^ .. г8 л
В* = А* =— К(1+^ +
2к*$2*)
Т(з).
О(з),
Б* =--К А* +
2в $
(*)
в = к, (1+т»)С*+
2к*$2*)
О(8),
где
О(х) =
О(8) = 2К
д(х), |ж| < а, —р(х), |ж| > а,
Т (х) =
О(х) егзх<1х, Т(з) = — 2 ■
\ т(х), |ж| < а, 1 0, |ж| > а,
Т (х) егзх<1х, Р = Ы.
(5)
Используя полученные выражения для коэффициентов А* и С*, = 1, 2, и значение интеграла [12]
Г™ 2 г
81§п(,§ ) ег8х(1з = —,
и —С© %
оо
эо
1
1
1
1
2
2
оо
эо
определим производные от смещений точек зоны контакта разнородных полос (-те < х < те, j = 1, 2):
(х, 0) = к^г((Я!) + (-1)' к^ У—те Т^
(X, 0) = -к^т(Х) + (-1у-^ /те
(з) ]и\
' Х ^ 2тт^ 7—те 5 - Ж
( -Х)(((8)<18 + [ к!^^ -Х)Т(з)(18.
— оо J—оо
(6)
Здесь
'(^аЛ = К°')(Ж) = 1 - ^ _Я_(
(Ж) = К22 (Ж)= ,1—те еИ? +
К$(х) = ¿(-1)^+11-Ч [те ¿"Ь,
К21 (Х) = г( 1) 2тт^ ]—те еЪр8ЪР + ? 6
Теперь, используя соотношения (6), удовлетворим последним двум соотношениям (1). В итоге, учитывая (5), после некоторых несложных выкладок для определения контактных напряжений д(х) и т(х) получим следующую систему определяющих сингулярных интегральных уравнений (-а < х < а):
. , , В [а ф)^ Г п .
Ад(х)--—^--К11(,в -х)а(8)й8 +
К }—а 8-Х ]—а
Г
+ ^12 (^ - ж) = Р1(х),
а
— а
а а
- Ат (х) - -Г ^^ - / К21(8 - х)д(8)с18 +
(7)
а
+ ^22(5 - ж)г(«)^5 = ^(ж),
а
где
С ^ Я
Е11(,) = Е22(,) = -/_те е^?^е
0 . . г Б /те е-'^' еИ | - | 4
Е12(Ж) = У—те еЬ^Л Я 16** е21Ы = * Ге—'" 8Ь |?| +
Л21(Х) 2тт /-со еИЯ + Я
8Ь |Я| + |Я| ;
оо
Л(х) = /(2>(х) - ^(х). Р2(х) = ¡!?\х) - !); А = хЩ!! - в = оф + с = А +1 - 2.)
(-1У ( Г—а Гр(з)(18
,2"(х) = +}. ) Т5 +
А! а + I ^К?^ -Х)Р(8)С13.
— те ■)а /
$\х) =([ а + [ ""Ы^ -х)р(8)с18. 3 = 1. 2.
\./—те ■)а /
Систему определяющих уравнений (7) нужно рассматривать вместе с условием равновесия верхней или нижней полуплоскости, которое при самоуравновешенных нагрузках можно представить в виде
I д(в =(! + / ) р( в )(1в = Р. I = 0. (8)
■)—а \./—те '¡а / ■)—а
Таким образом, решение поставленной задачи свелось к решению определяющей системы (7) при условиях (8). Отметим, что когда все полосы изготовлены из одного и того же материала, т. е. в случае однородной плоскости с периодической системой параллельных полубесконечных трещин А = с = = Л (х) =0 и из первого уравнения (7) для определения касательных контактных напряжений приходим к однородному сингулярному интегральному уравнению первого рода при втором однородном условии (8). Откуда следует, что т(х) = 0. Для определения же нормальных контактных напряжений из второго уравнения (7) получим следующее сингулярное интегральное уравнение первого рода
В Га [а „ . . . . , ^ . .
— /--+ Я21(8 - х)д(8)а8 = —г2(.х). -а < х < а.
К ]—а 8-х ]—а
которое нужно рассматривать при первом условии (8). Отметим, что при помощи несложных математических выкладок из полученного уравнения нетрудно получить определяющее интегральное уравнение этой же задачи, приведенной в [1]. В общем случае, чтобы построить решение определяющей системы интегральных уравнений, из первого уравнения (7) вычтем второе уравнение, первоначально умножив его на мнимую единицу. В результате, введя комплексную комбинацию контактных напряжений х(х) = д(х) -гт(х), приходим к следующему сингулярному интегральному уравнению второго рода:
В [а Х(8)(8
Ах(х) + В Г + Г Яп(8 - х)хШ8+ ■г .]—а 8 х ,]—а
а
Г
+ Я12(8 - х)х(в)й8 = Л(х). -а < х < а. (9)
■) —а
а
а
При этом условия (7) примут вид
/ х(^ = Р- (10)
■)—а
Здесь
дП(ж) = -Пц(х) + 21 [Д21 (х) - Д12(ж)] ,
^(х) = 21 [^21(х) + Д12(х)] , Г (х) = Г\(х) + г*Ь(ж),
а черточка над знаком функций здесь и в дальнейшем будет обозначать комплексно-сопряженную величину этих функций.
2. Составная плоскость с двумя полубесконечными межфазными трещинами. Прежде чем перейти к решению определяющей системы (9) в общем случае, рассмотрим еще один частный случай поставленной задачи, когда ширина полос стремится к бесконечности, т. е. рассмотрим задачу о плоско-деформированном состоянии составной упругой плоскости из двух разнородных полуплоскостей, содержащей симметрично расположенные две межфазные полубесконечные трещины. Заметим, что указанная задача рассматривается впервые и представляет самостоятельный интерес. В этом случае регулярные части в уравнении (9) исчезают и мы приходим к следующему определяющему сингулярному интегральному уравнению:
Ах(х) + —Г = Г (х), -а<х<а, (11)
К% 3—а $ х
которое нужно решать при условиях (10).
Решение уравнения (11) при условиях (10) будет даваться формулой [11]
( )_ 1 ¡г( N аш(х) Г Р\
Х(х) _ А(1 - а2) Г (х) ' Ла Ф)(* -х)Г
Р ) , .
+--—ш(х), -а<х<а,
'
где
ф) _ (х + а)—(1/2)—^(а -х)—(1/2)+^, а _ 7 _ ^ 1п ^ + ^ -
Вычислим коэффициенты интенсивности разрушающих напряжений и интеграл Черепанова—Райса в концевых точках трещины х _ ±а:
К [(-а)-7] _ К7 [(-а)7] - гКп [(-а)7] _ /2к 11ш ш—1(х)х(х) _
х^( — 1)] (а+0)
Л/2к (а [а Г 81ПК7„1 . _ „
-Р }, 3 _ 1, 2-
■I— Г
А ' —
--< --I - —
(1 - а2)\Акг .)—аш(,з)(,з - (-а)7) к
Тогда интеграл Черепанова—Райса в концевых точках трещины будет даваться формулой [13]
3(±а) _ /¿К(±а)К(±а),
где
1 /1 - 1 - \ в
2\ Е1 Е2 ) 4
3. Решение определяющего уравнения в общем случае. Решение определяющего сингулярного интегрального уравнения (9) в общем случае будем строить методом механических квадратур. Для этого при помощи замены переменных {х. з} = [а]. а£} запишем уравнение (9) на интервале (-1.1):
((]) + " (1 + /'1 Q* (Z - г,Ыт+
■ ]_1 С-] ]-1
<■1
1*
+ -т)Р(№ = (]]). (12)
где
а а
р(] = рх(а]]). р*(г]) = рал(а]]).
а- . Л 1 Г Ве "2|/?| + 2С/ - 2В\/3\
(т) = АQll(аT]) = 4К1й100-СЪ+*Ъ(3 + 13 еФ,/(р
/ Л а , . В Г™ ¿МЧр , , .
1(т]) = А^2(ат]) = ША ]_ж Ж/зЖ/зГ/з. 1 = н/а.
Условие (10) при этом принимает вид
(0^ = 1. (13)
Учитывая регулярность функции Q*j(т), несложно установить, что искомая функция р(т]) в концевых точках интервала интегрирования имеет показательную особенность и ее можно представить в виде
р(т) = р*(т)(1 + т)—(1/2)—^(1 - т)—(1/2)+^.
где р* (]) — непрерывная функция, ограниченная вплоть до концов отрезка [-1.1]. Подставляя значение функции р*(]) в (12) и (13) и используя соотношения, приведенные в [14], по стандартной процедуре, придем к системе алгебраических уравнений относительно значений р* (где ^ — корни многочлена Якоби р(—(1/2)+^-1/2)-1) (£) = о.
После определения функции р* (^) нетрудно восстановить функцию р(]), -1 < ] < 1, тем самым определив комплексную комбинацию безразмерных контактных напряжений, действующих в зоне стыка разнородных полос. После чего довольно просто определить комплексные коэффициенты интенсивности разрушающих напряжений и, следовательно, известные 7-интегралы Черепанова—Райса по формулам
К (±1) = К7 (±1) - г Кц (±1) =
= Т ]){1/2)тг1 [аЦ)(ат) - гт$(ат)] =
= ^L lim (1 ^ ^(V2^^) = ^^2-(1/2)*V(±1), а а
J (±1) = ßK (±1)К(±1) = ^И±1)|2.
Выпишем также формулу для определения раскрытия трещины. С этой целью, используя второе из соотношений (6), составим нормальную составляющую дислокаций смещений точек берегов полубесконечных трещин и, перейдя на интервал (-1,1), представим их при помощи функции ^(rj) в следующем виде:
d
v'(arj) = v[(arj) - v'2(ar]) = -1 +
+ J1 Qi(£ - V) Re^OdÇ + Ц1 Q2(Ç - V) Im^(£)d£} - F2{ar}),
где
-aD , Л i Г е-1 ßl sh |ß | + Iß | ivß/,a
Qi(n) = -^iM = -Chßstiß + ß1 ^ß/ldß,
ße ißl/ldß
fcM = ^амщ ) = ^ £
Б у " 2В1 у-те е^ + £
Тогда безразмерное раскрытие трещин, находящихся соответственно на интервалах ( а, те) и (-те, -а), можно определить по формулам:
v( arj) a
v( a tï ) f1 V*(rf) = = v'(a£)d£, -те < r] < -1.
a 1
,u( a il) f1
v*(rj) = = v'(a£)d£, 1 < r] < те,
a Ji
v( a r] ) I"11
a J-i
4. Численные расчеты. Проведен численный расчет и изучены закономерности изменения контактных напряжений, действующих в зонах контактов разнородных полос, и безразмерных интегралов Черепанова—Райса J*(±a) = 102 ■ (а/Р) ■ J(±1) в концевых точках трещин в зависимости от соотношения у = в случае фиксированных значений коэффициентов Пуассона Vj, j = 1, 2, а также от изменения параметра I = h/а. При этом считается, что кусочно-однородная плоскость деформируется под воздействием симметрично расположенных относительно оси Оу сосредоточенных нагрузок величины Р/2, т. е. принято (а/Р) ■ р(аг]) = [ô(r] — b) + ô(r] + b)]/2, b = щ/а, где щ — расстояние сосредоточенных нагрузок от начала координат и (Р/а) ■ = 0.1.
Результаты численных расчетов приведены на рис. 3-6.
На рис. 3 приведены графики безразмерных контактных напряжений в зависимости от параметра у в случае, когда v\ = 0.3, v2 = 0.25, 6 = 2 и I = 1. Они показывают, что нормальные контактные напряжения мало зависят от параметра у (графики при у = 1 и у = 0.5 практически совпадают), в то время как касательные контактные напряжения при приближении параметра у к единице стремятся к нулю.
Г1 Г]
Рис. 3. Безразмерные контактные напряжения (при различных значениях параметра слева — нормальные контактные напряжения, справа — тангенциальные контактные напряжения
[Figure 3. Dimensionless contact stresses when = 0.3, = 0.25, b = 2, and I = 1 for different values of the parameter on the left — normal contact stresses; on the right —
tangential contact stresses]
На рис. 4 приведены графики изменения значения интеграла Черепанова— Райса ,]*(а) = ,]*(-а) в концевых точках трещин в зависимости от параметра Ь, описывающего степень удаленности точек приложения внешней нагрузки от вершин трещин, в случае, когда и\ = 0.3, = 0.25 и I = 1 при некоторых значениях Из них видно, что при удалении точек приложения внешних нагрузок от концевых точек трещин интеграл Чере-панова—Райса, монотонно убывая, стремится к определенному пределу. Причем чем модуль сдвига первой из полос больше, тем меньше значение интеграла Черепанова—Райса, т. е. меньше вероятность распространения трещины.
На рис. 5 приведены графики безразмерных контактных напряжений в зависимости от параметра I в случае, когда ц, = 2, Ь = 2, и\ = 0.3 и и2 = 0.25. Они показывают, что при возрастании I, т. е. при увеличении ширины полосы, когда длина контактной зоны остается неизменной, нормальные контактные напряжения в средней части контактной зоны уменьшаются, а у концевых точек зоны контакта увеличиваются. Касательные же контактные напряжения при этом по абсолютной величине уменьшаются.
На рис. 6 приведены графики интеграла Черепанова—Райса ,]*(а) = ,]*(-а) в вершинах трещин в случае, когда = 0.3 р2 = 0.25, Ь = 2: слева — в зависимости от параметра I при значениях параметра ^ = 0.5, 1 и 5 и справа — в зависимости от параметра ^ при значениях параметра I = 0.5, 1 и 5.
Из рис. 6 слева следует, что при увеличении ширины полос интеграл ■1*(±а) возрастает, стремясь к определенному пределу, соответствующему значению интеграла Черепанова—Райса в случае двухкомпонентной плоскости с двумя симметричными полубесконечными трещинами. При этом, как и выше, чем жестче первая из полос, тем меньше значение интеграла Черепа-нова—Райса.
Рис. 4. Изменение интеграла J*(a) [Figure 4. Variation of J-integral depending on the parameter b when V! = 0.3, V2 = 0.25, I = 1]
-0.2
Рис. 5. Безразмерные контактные напряжения (при различных значениях параметра I): слева — нормальные контактные напряжения, справа — тангенциальные контактные напряжения
[Figure 5. Dimensionless contact stresses when ^ = 2, b = 2, v-^ = 0.3, and = 0.25 for different values of the parameter I: on the left — normal contact stresses; on the right —
tangential contact stresses]
J* (a) J* (a)
Рис. 6. Изменение интеграла J* (a) [Figure 6. Variation of J-integral when v-^ = 0.3 = 0.25, b = 2: on the left — depending on the parameter I; on the right — depending on the parameter
Из рис. 6 справа следует, что при возрастании параметра ^ значение интеграла, монотонно убывая, стремится к постоянному значению, которое представляет из себя значение интеграла Черепанова—Райса в случае, когда первая из полос абсолютно жесткая.
5. Заключение. Таким образом, выведено определяющее уравнение плоской задачи теории упругости для равномерно кусочно-однородного пространства с периодической системой двух полубесконечных параллельных межфазных трещин в виде сингулярного интегрального уравнения второго рода относительно комплексной комбинации контактных напряжений, действующих в зонах стыка разнородных полос.
Выяснено поведение искомой функции в концевых точках интервалов интегрирования и построено решение определяющего уравнения численно-аналитическим методом механических квадратур. Получены простые формулы для определения комплексного коэффициента интенсивности разрушающих напряжений и .]-интеграла Черепанова—Райса в концевых точках трещин.
При помощи численных расчетов изучено поведение контактных напряжений и 7-интеграла в зависимости от физико-механических и геометрических характеристик задачи. Показано, что при выбранных параметрах увеличение жесткости одной из полос, когда жесткость второй полосы остается неизменной, значение ■]-интеграла стремится к определенной постоянной, со-
ответствующей значению интеграла Черепанова—Райса в случае, когда первая из полос абсолютно жесткая.
Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем.
Авторский вклад и ответственность. В.Н. Акопян —идея исследования, поиск методов решения, руководство и консультирование, чистовик рукописи. А.А. Григорян — анализ литературы, реализация алгоритмов на компьютере, расчеты и визуализация результатов, черновик рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.
Финансирование. Исследование выполнялось без финансирования.
Библиографический список
1. Панасюк В. В., Саврук М. П., Дацышин А. П. Распределение напряжений около трещин в пластинах и оболочках. Киев: Наукова думка, 1976. 443 с.
2. Мусхелишвили Н. И. Некоторые основные задачи математической теории упругости. М.: Наука, 1966. 708 с.
3. Развитие теории контактных задач в СССР/ ред. Л. А. Галин. М.: Наука, 1976. 493 с.
4. Попов Г. Я. Концентрация упругих напряжений возле штампов, разрезов, тонких включений и подкреплений. М.: Наука, 1982. 344 с.
5. Бережницкий Л. Т., Панасюк В. В., Стащук Н. Г. Взаимодействие жестких линейных включений и трещин в деформируемом теле. Киев: Наукова думка, 1983. 288 с.
6. Барзокас Д. И., Фильштинский Л. А., Фильштинский М. Л. Актуальные проблемы связанных физических полей в деформируемых телах. Т. 1. М., Ижевск, 2010. 864 с.
7. Акопян В. Н. Смешанные граничные задачи о взаимодействии сплошных деформируемых тел с концентраторами напряжений различных типов. Ереван: Гитутюн, 2014. 322 с.
8. Попов Г. Я. Избранные труды. В 2-х томах. Одесса: ВМВ, 2007.
9. Акопян В. Н., Даштоян Л. Л. Разрывные решения двоякопериодической задачи для кусочно-однородной плоскости с межфазными дефектами // Мех. композ. матер., 2017. Т. 53, №5. С. 863-879.
10. Hakobyan V. N., Sahakyan A. V., Aghayan K. L. Periodic problem for a plane composed of two-layer strips with a system of longitudinal internal inclusions and cracks / Wave Dynamics, Mechanics and Physics of Microstructured Metamaterials. vol. 109. Cham: Springer, 2019. pp. 11-22. https://doi.org/10.1007/978-3-030-17470-5_2.
11. Hakobyan V. N., Grigoryan A. H. Anti-plane stressed state uniformly piece-homogeneous space with a periodic system of parallel semi-infinite interfacial cracks // J. Phys.: Conf. Ser., 2020. vol. 1474, 012017. https://doi.org/10.1088/1742-6596/1474/1Z012017.
12. Брычков Ю. А., Прудников А. П. Интегральные преобразования обобщенных функций / Итоги науки и техн. Сер. Мат. анал., Т. 20. М.: ВИНИТИ, 1982. С. 78-115.
13. Stress Intensity Factors Handbook. vol. 1 / ed. Y. Murakami. Oxford: Pergamon, 1987.
14. Sahakyan A. V., Amirjanyan H. A. Method of mechanical quadratures for solving singular integral equations of various types// J. Phys.: Conf. Ser., 2018. vol. 991, 012070. https:// doi.org/10.1088/1742-6596/991/1/012070.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2021, vol. 25, no. 1, pp. 67-82_
ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1816
MSC: 74A45
Plane stress state of a uniformly piece-wise homogeneous plane with a periodic system of semi-infinite interphase cracks
© V. N. Hakobyan, A. A. Grigoryan
Institute of Mechanics, National Academy of Sciences of the Republic of Armenia,
24B, Marshal Baghramyan ave., Yerevan, 0019, Republic of Armenia.
Abstract
The plane stress state of a uniformly piecewise-homogeneous plane obtained by alternately joining two dissimilar strips is considered, which along the lines of joints of dissimilar strips is weakened by a periodic system of two semi-infinite interfacial cracks and is deformed using normal loads applied to the crack banks. The basic cell of the problem in the form of a two-component strip is considered and, using the generalized Fourier transform, a governing system of equations for the problem is obtained in the form of one singular integral equation of the second kind for a complex combination of contact stresses in the junction zone of the strips.
As a special case, tending the height of the strips to infinity, the governing equation of the problem for a two-component plane of two dissimilar halfplanes with two semi-infinite interfacial cracks is obtained and its exact solution is constructed. The governing equation for the stated problem is also obtained in the form of one singular integral equation of the first kind with respect to normal contact stresses in another particular case, when all strips are made of the same material, i.e. in the case of a homogeneous plane, a weakened periodic system of parallel, two semi-infinite cracks.
In the general case, the behavior of the unknown function at the end points of the integration interval is determined and the solution of the problem by the numerical-analytical method of mechanical quadratures is reduced to solving a system of algebraic equations. Simple formulas are obtained to determine the intensity factors, the Cherepanov-Rice integral and crack opening. A numerical calculation has been performed. Regularities of changes in contact stresses and the Cherepanov-Rice integral at the endpoints of cracks are revealed, depending on the elastic characteristics of heterogeneous strips and the geometric parameters of the problem.
Research Article
3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Hakobyan V. N., Grigoryan A. A. Plane stress state of a uniformly piece-wise homogeneous plane with a periodic system of semi-infinite interphase cracks, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2021, vol. 25, no. 1, pp. 67-82. https://doi.org/10.14498/vsgtu1816 (In Russian). Authors' Details:
Vahram N. Hakobyan https://orcid.org/0000-0003-3684-9471
Dr. Phys. & Math. Sci., Professor; Chief Researcher; Dept. of Elasticity and Viscoelasticity;
e-mail: [email protected]
Aram A. Grigoryan © https://orcid.org/0000-0001-7582-1960 Postgraduate Student; Junior Researcher; Dept. of Elasticity and Viscoelasticity; e-mail: [email protected]
80
© Samara State Technical University
Keywords: periodic problem, mixed boundary value problem, piece-wise homogeneous plane, interface cracks.
Received: 10th August, 2020 / Revised: 29th January, 2021 / Accepted: 12th February, 2021 / First online: 24th February, 2021
Competing interests. We declare that we have no conflicts of interest in the authorship and publication of this article.
Authors' contributions and responsibilities. V.N. Hakobyan: Idea of study; Search for methods of solution; Supervision and consulting; Writing — review & editing. A.A. Grigoryan: Literature review; Implementation of the computer algorithms; Performing calculations and visualizing results; Writing — original draft. The authors are absolutely responsible for submitting the final manuscript in print. Each author has approved the final version of manuscript.
Funding. This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.
References
1. Panasyuk V. V., Savruk M. P., Datsyshin A. P. Raspredelenie napriazhenii okolo treshchin v plastinakh i obolochkakh [Stress Distribution Around Cracks in Plates and Shells]. Kiev, Naukova Dumka, 1976, 443 pp. (In Russian)
2. Muskhelishvili N. I. Some Basic Problems of the Mathematical Theory of Elasticity. Netherlands, Springer, 1977, xxxi+732 pp. https://doi.org/10.1007/978-94-017-3034-1
3. Razvitie teorii kontaktnykh zadach v SSSR [The Development of Theory of Contact Problems in USSR], ed. L. A. Galin. Moscow, Nauka, 1976, 493 pp. (In Russian)
4. Popov G. Ya. Kontsentratsiia uprugikh napriazhenii vozle shtampov, razrezov, tonkikh vkli-uchenii i podkreplenii [Concentration of Elastic Stress around Stamps, Cuts, Thin Inclusions, and Reinforcements]. Moscow, Nauka, 1982, 344 pp. (In Russian)
5. Berezhnitskii L. T., Panasyuk V. V., Stashchuk N. G. Vzaimodeistvie zhestkikh lineinykh vkliuchenii i treshchin v deformiruemom tele [The Interaction of Rigid Linear Inclusions and Cracks in a Deformable Body]. Kiev, Naukova Dumka, 1983, 288 pp. (In Russian)
6. Barzokas D. I., Fil'shtinskii L. A., Fil'shtinskii M. L. Aktual'nye problemy sviazannykh fizich-eskikh polei v deformiruemykh telakh [Actual Problems of Coupled Physical Fields in Deformable Bodies], vol. 1. Moscow, Izhevsk, 2010, 864 pp. (In Russian)
7. Hakobyan V. N. Smeshannye granichnye zadachi o vzaimodeistvii sploshnykh deformirue-mykh tel s kontsentratorami napriazhenii razlichnykh tipov [Mixed Boundary-Value Problems on the Interaction of Continuous Deformable Bodies with Stress Concentrators of Various Types]. Yerevan, Gitutiun, 2014, 322 pp. (In Russian)
8. Popov G. Ya. Izbrannye trudy [Selected Works]. In 2 volumes. Odessa, VMV, 2007 (In Russian).
9. Hakobyan V. N., Dashtoyan L. L. Discontinuous solutions of a doubly periodic problem for a piecewise homogeneous plate with interphase defects, Mech. Compos. Mater., 2017, vol.53, no. 5, pp. 601-612. https://doi.org/10.1007/s11029-017-9690-8.
10. Hakobyan V. N., Sahakyan A. V., Aghayan K. L. Periodic problem for a plane composed of two-layer strips with a system of longitudinal internal inclusions and cracks, In: Wave Dynamics, Mechanics and Physics of Microstructured Metamaterials, vol. 109. Cham, Springer,
2019, pp. 11-22. https://doi.org/10.1007/978-3-030-17470-5_2.
11. Hakobyan V. N., Grigoryan A. H. Anti-plane stressed state uniformly piece-homogeneous space with a periodic system of parallel semi-infinite interfacial cracks, J. Phys.: Conf. Ser.,
2020, vol. 1474, 012017. https://doi.org/10.1088/1742-6596/1474/1Z012017.
12. Brychkov Yu. A., Prudnikov A. P. Integral transforms of generalized functions, J. Soviet Math., 1986, vol.34, no. 3, pp. 1630-1655. https://doi.org/10.1007/BF01262407.
13. Stress Intensity Factors Handbook, vol. 1, ed. Y. Murakami. Oxford, Pergamon, 1987.
14. Sahakyan A. V., Amirjanyan H. A. Method of mechanical quadratures for solving singular integral equations of various types, J. Phys.: Conf. Ser., 2018, vol.991, 012070. https:// doi.org/10.1088/1742-6596/991/1/012070.