Глаголев В.В., Глаголев Л.В., Маркин А.А. Энергетическое произведение в модели трещиноподобного дефекта при нагружении типа моды II // Вестник Пермского национального исследовательского политехнического университета. Механика. 2019. № 4. С. 48-58. DOI: 10.15593/perm.mech/2019.4.05
Glagolev V.V., Glagolev L.V., Markin A.A. Energy product in a crack-like defect model under loading of mode II type. PNRPU Mechanics Bulletin, 2019, no. 4, pp. 48-58. DOI: 10.15593/perm.mech/2019.4.05
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 4, 2019 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
DOI: 10Л5593/perm.mech/2019A05 УДК 539.375
ЭНЕРГЕТИЧЕСКОЕ ПРОИЗВЕДЕНИЕ В МОДЕЛИ ТРЕЩИНОПОДОБНОГО ДЕФЕКТА ПРИ НАГРУЖЕНИИ ТИПА МОДЫ II
В.В. Глаголев, Л.В. Глаголев, А.А. Маркин
Тульский государственный университет, Тула, Россия
О СТАТЬЕ АННОТАЦИЯ
Рассматривается нагружение трещиноподобного дефекта по моде II. В отличие от классического представления трещины в виде математического разреза предлагаемая модель определяет трещину в виде физического разреза с характерным линейным размером. Мысленное продолжение физического разреза в твердом теле формирует слой взаимодействия. Существенно, что напряженно-деформированное состояние слоя не вносит в модель трещины сингулярность. Произведение приращения удельной свободной энергии в торцевом квадратном элементе слоя на линейный размер определяет его энергетическое произведение (ЭП). Объектом исследования является двухконсольный образец, а предметом исследования - энергетическое произведение в торцевом элементе слоя взаимодействия. Внешняя нагрузка консолей приводит к их горизонтальным антисимметричным перемещениям, которые формируют в слое взаимодействия однородные сдвиговые деформации. Из условий равновесия консолей в вариационной форме с учетом гипотезы однородности осевых деформаций и с учетом их обжатия получена система дифференциальных уравнений, связывающая напряженное состояние в слое и консолях. Исследовано решение характеристического уравнения системы при различных отношениях толщин слоя и консолей. Показано, что при отношениях меньше определенного значения, зависящего от коэффициента Пуассона, имеют место вещественные корни. В рамках вещественных корней характеристического уравнения получено аналитическое решение поставленной задачи. При условии пренебрежения обжатием консолей найдено упрощенное решение. Определены деформации в слое с учетом обжатия консолей и без него. Проведен анализ зависимости ЭП от отношений толщин слоя и консолей. Показано, что при отношении толщин 10-6 и менее ЭП практически не меняет своего значения. Учет обжатия консолей дает различие в значениях ЭП порядка 20 % по отношению к упрощенному решению задачи.
©ПНИПУ
Получена: 11 марта 2019 г. Принята: 14 ноября 2019 г. Опубликована: 30 декабря 2019 г.
Ключевые слова:
энергетическое произведение, линейный размер, трещина моды II.
© Глаголев Вадим Вадимович-д.ф.-м.н., проф., зав. каф., e-mail: vadim@tsu.tula.ru, : 0000-0003-0371-7704
Глаголев Леонид Вадимович - к.ф.-м.н., инж., e-mail: Ien4ic92@gmail.com, : 0000-0003-2313-2084 Маркин Алексей Александрович - д.ф.-м.н., проф., e-mail: markin-nikram@yandex.ru, : 0000-0003-1456-8281
Vadim V. Glagolev - Doctor of Physical and Mathematical Sciences, Professor, Head of Department, e-mail: vadim@tsu.tula.ru, : 0000-0003-0371 -7704
Leonid V. Glagolev - CSc in Physical and Mathematical Sciences, Engineer, e-mail: len4ic92@gmail.com,
: 0000-0003-2313-2084 Alexey A. Markin - Doctor of Physical and Mathematical Sciences, Professor, e-mail: markin-nikram@yandex.ru, : 0000-0003-1456-8281
.-4 ' Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International Qy License (CC BY-NC 4.0)
л; This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
ENERGY PRODUCT IN A CRACK-LIKE DEFECT MODEL UNDER LOADING OF MODE II TYPE
V.V. Glagolev, L.V. Glagolev, A.A. Markin
Tula State University, Tula, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 11 March 2019 Accepted: 14 Novembert 2019 Published: 30 December 2019
Keywords:
energy product, linear size, crack mode II.
The loading of a crack-like defect in mode II is considered. In contrast to the classical representation of a crack in the form of a mathematical cut, the proposed model defines a crack in the form of a physical cut with a characteristic linear size. The mental continuation of a physical cut in a solid forms an interaction layer. It is significant that the stress-strain state of the layer does not introduce a singularity to the crack model. The product of the increment of the specific free energy in the face square element of the layer by the linear size determines its energy product. The object of the study is a double-cantilever sample, and the subject of study is the energy product in the face element of the interaction layer. The external load of the cantilevers leads to their horizontal antisymmetric displacements, which form uniform shear deformations in the interaction layer. From the equilibrium conditions of the cantilevers in the variation form, taking into account the hypothesis of axial deformation homogeneity and their reduction, a system of differential equations is obtained, which relates the stress state in the layer and the cantilevers. The solution of the characteristic equation of the system is investigated for various ratios of layer thickness and cantilevers. It is shown that when the relationship is less than a certain value, depending on the Poisson's ratio, real roots take place. In the framework of the real roots of the characteristic equation, an analytical solution of the problem is obtained. Subject to the neglect of compression cantilevers found a simplified solution. The deformations in the layer are determined taking into account the compression of the consoles and without it. The analysis of the dependence of the energy product on the relationship of the thickness of the layer and cantilevers. It is shown that with a thickness ratio of 10-6 or less, the energy product practically does not change its value. Accounting for the compression of cantilevers gives a difference in the values of the energy product of the order of 20 % in relation to the simplified solution of the problem.
©PNRPU
Введение
Классическая постановка задачи о трещине моды II рассматривает бесконечное линейно-упругое пространство с дефектом в виде математического разреза, поверхности которого нагружены противоположными векторами напряжений [1, 2]. В этом случае на продолжении математического разреза в слое нулевой толщины реализуются только сдвиговые напряжения. Для образцов конечных размеров предлагается ряд схем нагружения [3-6], в которых на продолжении трещины реализуется напряженное состояние, соответствующее поперечному сдвигу. Данные схемы также принято рассматривать как трещину моды II. Отметим, что для образцов конечных размеров общие решения строятся, как правило, на основе численных процедур, основанных на методе конечных элементов [7-12]. Аналитические решения получаются, как правило, в рамках упрощающих гипотез [13-19]. Так, в работе [6] для консоли с центральной трещиной на основе гипотез теории пластин Тимошенко получено решение задачи и определено значение 1-интеграла. В этом случае в консолях учитывались деформации растяжения и сдвига. Влияние обжатия консолей в результате внешнего воздействия не рассматривалось.
Альтернативной моделью трещины в виде математического разреза является разрез физический с характерной толщиной 50. В данном представлении особую роль играет описание взаимодействия материального объема, лежащего на продолжении физического разреза, со смежной средой. Важно, чтобы та или иная модель не приводила к сингулярности в области окончания физического разреза. При нагружении трещинопо-добного дефекта нормальным отрывом отметим модель в виде связей Прандтля, предложенную в статье [20] и используемую в работах [21-23]. В модели Ф. Макк-линтока [24] предложено использовать для материала слоя усредненные характеристики напряженно-деформированного состояния (НДС). Постановки конкретных задач исходя из представления Ф. Макклинтока были реализованы на основе теории слоя взаимодействия в работах [18-20], [26-28]. В частности, в статье [27] введено понятие энергетического произведения (ЭП) в виде произведения приращения удельной свободной энергии в элементе слоя и толщины слоя. В статье [18] о нагружении двухконсольной балки (ДКБ-образца) нормальным отрывом (нагружение по моде I) показано, что при уменьшении отношения толщины слоя к высоте консоли имеет место сходимость ЭП к поверхностной энергии, рассчитанной через вязкость разрушения Кс .
В данной работе предложена схема нагружения ДКБ-образца в которой реализуется состояние, близкое к моде II. На основе общей вариационной постановки задачи [26-28] получены упрощенные постановки с учетом и без учета обжатия консолей. Из аналитических решений показано, что при уменьшении отношения толщины слоя к толщине консоли имеет место сходимость ЭП в слое для данного вида нагружения.
1 0.55„
22 (X )=5" J СТ22 (Х2 ) dX2
50 -0.550 1 0.550
11 ( X1 )=5" J CT11 ( X1' X2 ) dX2'
( X ) =
0 -0.550 iu2>(X1 )- U2-(X1 ) Л
1. Постановка задачи
Рассмотрим случай нагружения ДКБ-образца согласно схеме, показанной на рис. 1. Правые торцы консолей 1 и 2 образца свободны от внешней нагрузки. Верхняя граница консоли 1 и нижняя граница консоли 2 закреплены от вертикальных перемещений. На левых торцах консолей действует горизонтальная антисимметричная нагрузка интенсивностью Р .
a H 1
*— À À h A (7) A
0 ! © I
— ©
Рис. 1. Схема сдвигового воздействия на тонкий слой Fig. 1. Shear effect on a thin layer
Для описания взаимодействия слоя 3 с телами 1 и 2 применим концепцию «слоя взаимодействия», развитую в работах [26-28]. В этом случае условия равновесия тел 1 и 2 запишем в вариационной форме для тела 1 :
J о ••5£ds +Jct22 5u+dx1 +J<512 5u+dx1 +
S1 t t
I d5u+ d5u+ I 1 + 0,5501 fä11-l—dx1 + fä12-2-dx1 I = f P •5udl (1)
I \ dx1 \ dx1 ) L
и тела 2
J о • -5sds - J ct225u2 dx1 - J ст125u1 dx1
+ 0,5501 |"стп ^^ ёх1 + [ст12 1 = [Р2-8ий?/, (2)
и дх1 \ дх1 ) I
где и - вектор перемещений; Р1 ,Р2 - вектора внешней нагрузки тел 1 и 2; ^, - площади тел 1 и 2; с, г - тензоры напряжений и деформаций; с, г - тензоры средних напряжений и деформаций слоя с компонентами:
1 0.55о
°21 ( Х1 ) = ^12 (X )=— | (х2 ) dX2,
d5u,
еп (x1 ) = О,5
du+ ( x1 ) du1 ( x )
ôx1
ôx1
(3)
821 ( x1 ) ^12 ( x1 )
= 0,5
u+ ( x1 )- u-( x1 )
+ 0,5
du2+ (x1 ) du1 (x1 )
ôx1
ôx1
(4)
где щ - компоненты векторов перемещений верхней и нижней границ слоя соответственно; к = 1,2; Ц, Ь2 -граница приложения внешней нагрузки для тела 1 и 2. Постулируется жесткое сцепление между границами области 3 и областями 1, 2:
и+ = и(X,80/2); и- = и(х,- 80/2) х е[0;^]. (5)
Примем определяющие соотношения в форме закона Гука:
E I V s
(1+)^ + T-V8^
(6)
где Е, V - модуль упругости и коэффициент Пуассона тела; в - объемное расширение, в = вп +е22 +е33; 8в -символ Кронекера; I, ] = 1,2,3.
Для материала слоя взаимодействия 3 определяющие соотношения считаем справедливыми для средних компонент тензоров напряжений и деформаций:
1 + V
— E I V СТ.. =-1 8,, +--85,.,. I.
1 - 2v
(7)
Распределение перемещений в консолях при данном нагружении принимаем в следующем виде:
uf (x1, x2 ) = u+ ( x ),
u21)( x1, x2 ) = y 1 ( x1 )( h + 50/2 - x2 ),
u1(2) (X1, x2 ) = u1-(x1 ),
u22) (x, x2 ) = y 2 (x1 ) (h +50 /2 + x2 ).
(8)
(9)
Вертикальные перемещения тождественно удовлетворяют граничным условиям на верхней границе консоли 1 и нижней границе консоли 2:
,W
1x2 =h+50/2
= 0, u
(2)
1x2 =-h-50/2
= 0.
S
2
2
0 -0.55
В этом случае деформацию консолей будут определять компоненты тензора деформаций
su (Xj, х2)
du± (x,)
dx,
( Xj, X2 ) = +Vk ,
sf2 (x,, x2) = 0,5 d^ (Xj) (h + S0/2 + x2)
(10) (11) (12)
где и± - горизонтальное перемещение верхней (+) или нижней (-) границы слоя; символ + преобразуется в знак (-) для к = 1, а для к = 2 преобразуется в знак (+).
При распределении деформаций в консолях в виде (10)-(12) работа внутренних напряжений в теле 1 рассчитывается следующим образом:
J а Seds = {" J
-0 лА+50/2 -а J 50/2
(
ст,
d 5и+ dx,
-ct 22 +ct12 x
d dx1
(h + 5 0/2 - X2)
dx1dx2 +
//
h+50/2
J+0J 5„/2 ( .
+ CT,
JJ
J+0 Ji
d dx,
d 5m,+
CT11 --CT225Vl +
dx1
(h + 5„/2 - x2)
dx, dx2,
в теле 2:
Г ? j С-0 С-5^2 d5U, 5
J ^ -5£dS = J-а j-h-5^2 CTl1 ~~dX +CT-^2 +
(
d 5m,
v dx,
+ CT,,
^ (h + 5„/2 + x2)
v dx,
dx1 dx2 +
//
bi+01
-50/2
+ CT,
+0 J -h-50 /2
' d 5y 2 v dx,
d 5m,
ct,, ——- + ст225у2 +
v 11 dx,
W
(h + 50 /2 + x2)
//
Введем в рассмотрение обобщенные силы:
Ю (Х1 ) = 1^/2СТ11^Х2.
А+50/2
6,2 (xi ) = i5o+/20 CT22 dx2
2 / \ f-5()/2
Öll (xi ) = J-h-50/2 CTUdx2,
-h-5^2
Ö222 (x, ) = J
-5^2 h-5^/2 '
dx^
и обобщенные моменты:
1 / \ Г^+8о/2 / / \
Ми2 (Х1 ) = ]8 /2 °12 (А + 80/2-Х2 ) -Х2
(13)
dx,dx2. (14)
(15)
(16)
(17)
(18)
(19)
М122 (X ) = Х (Л + 8о/2 + *2 ) ^ (20)
Интегрируем по частям ряд слагаемых в (13), (14) с учетом (15)-(20) и четвертое и пятое слагаемые в левой части (1) и (2) с учетом (8) и (9):
|-061 ^=а>±Г=" -Г0°8и±(21)
-х, и1 =-а •,-а ахг
г« , -8и± , , п=1 г1 сЮк, .
¿Ю ^ = бХ* 1X1=+0 - 1+0^5и±ёх1, (22)
'x =+0 J+0 dx,
i-0 , d5ш, , ix, =-0 i*-0 dMMk
Mkn —^dx, = M^* -J 05^1,(23)
-а dX x, =-а J-а r]\-
J Mf2 Äd = M1k,5yk p=' -f -dMll5Vkdx1, (24)
J+0 12 Ay 1 12 Tklx,=+0 J+0 Ay k 1
dx.
dM
dx
J ct,,-— dx = ct,,5m,+ Г -J —115u±dx,. (25)
J+0 11 Ay 1 11 1 Ix, =+0 J+0 Ay 11 v '
x =+0 J+° dx,
d ct,,
J+0 11 dx,
x =+0 J+0 dx,
et _ d5m2 r
I ст12-dx, = I
J+0 12 Ay 1 Jz
_ d5hyk
J+0 12 dx,
dx
dx, =
=Л ,1х;=+0 -|+0 ^Хт*-х1 ]. (26)
Рассмотрим правые части уравнений (1) и (2). При положительном обходе контура работа внешних напряжений рассчитывается так:
, 50/2 J P -5ud/ = - J Р5и+ (-dx2)
5^2
J
L, h+50/2
= - Р05М+ | =-Ö15M1+[
2 5^2
J P -5ud/ = - J P5m- (-dx2)
L, h +50/2
= Ph5 m,- = 6,5m-
(27)
(28)
где Ю1 = РЛ.
Подставив (13),(15), (16) (21)-(28) в (1) и приравняв слагаемые при равных вариациях, приходим к двум системам дифференциальных уравнений: для участка х1 е[-а;0)
dx,
для участка x, е (0; t]
dM,2 + öj2 = 0, d^l = 0,
dM1
dx,
d ct,.
(29)
dx,
12 + 6,2 + 0,550= hCT22,
dx
dÖli dx,
+ 0,55,
d ct,, dx
= ст,
(30)
s
X
+
s
01а§в1ву V. V., 01а§в1ву Ь. V., Маткт А. А. /Р№Р1?Мвскст^ БиПеШ 4 (2019) 48-58
с условиями сопряжения
правом торце
и, = и, ,
1 1х, =-0 1 1х, =+0 '
х. =-0 = ^1 х. =+0,
МЦ |х =-0 =( М1 + 0,580 кстИ )| х 011! |х =-0 =(б1\ + 0,580ап )|
х1 =+0
и с естественными граничными условиями: на левом торце
(31)
(32)
(33)
(34)
= е-
М1
= 0,
на правом торце
(е111 + 0,5^ )| х = = 0,
(М12 + 0,580кст12 )| 0.
(35)
(36)
(37)
(38)
Аналогичным образом из (2) приходим к системам дифференциальных уравнений: для участка х1 е[-с;0)
йМ2и йх,
=е2
для участка х1 е (0; £]
йМ,
йх1
й ст.,
= 0,
(39)
_ ^ + 0,580^^йх1"+ М22 = 0,
е + 0,580 ^ = -ст21,
йх * йх.
с условиями сопряжения
^2 1 х1 =-0 =^2| х1 =+0,
М2 |х=-0 =( М2 + 0,580 ист12 )| х=-0=(е2+0,580стп )|
х1 =+0
х1 =+0
и естественными граничными условиями: на левом торце
а2! =^ = -е,,
(40)
(41)
(42)
(43)
(44)
М,2.
= 0,
(45)
(46)
(б + 0,580стп )| х =^ = 0, (М122 + 0,580кст12 )| 0.
(48)
С учетом выражений (10)-(12) и условием плоской деформации (е33 = 0) определяющие соотношения (6) запишем в виде
к йи ±
стк1 = т к,
йх1
к йи ±
ст22 = + ту к + т1-г~
йх1
стк2 = Ь | ^ (к +8„/2 + х2 )
где т = -
Е (1 -V)
йх1
-; т1 = -
(49)
(50)
(51)
Е V
-; Ь = ■
Е
(1 + v)(l-2vУ ~ (1 + v)(l - 2v)' 2 (1 + V) *
Из (7)-(9), (3) и (4) получим связь средних напряжений в слое взаимодействия с его граничными перемещениями:
тй (< +и-) т>1,, ч
8„
йх.
= Т~ к (^1 ) +
т1
й (и+ + и1 ) йх1
ст12 = Ь
(+ 0,5к Г^ + йу ЛЛ
(52)
(53)
(54)
Таким образом, система дифференциальных уравнений (29), (30), (39), (40) с учетом связей (10)-(12), (15)-(20) и определяющих соотношений (49)-(54) является замкнутой относительно четырех неизвестных функций , м-, у 2. Каждая искомая функция входит в систему со второй производной включительно. Удовлетворение восьми граничных условий (35)-(38), (45)-(48) с учетом сопряжения решений (31)-(34), (41)-(44) даст решение поставленной задачи.
2. Решение задачи
Рассмотрим решение поставленной задачи в виде и+ = -и" , у 1 = у2 для полубесконечного образца (£ ^ да ). В этом случае движение консолей противоположно, с одинаковым модулем обжатия. Для решения задачи достаточно рассмотреть только систему уравнений (29), (30) при следующих выражениях обобщенных сил:
, , , |*к+80/2 I йи,
б (х1) = Ц ст11йх1 = к| т-т1у1
йи +
йх,
(55)
22
2
-т-
-и+ —X,
(56)
1 / \ Г"+00/2
Ю2 (Х1 ) = ^/2 Х22 —Х2 = И | -БХ1 + обобщенного момента < (Х1 ) = ^Х (И + 80/2-Х2) —Х2 = ^—Х (57) и средних напряжений в слое
Стц = 0,
х22 = 0,
Х12 = Т
2и+ , — х, —— + И——
80 —хх
(58)
(59)
(60)
Проинтегрируем второе уравнение системы (29) с учетом граничного условия (35) и связи (55):
,+ 01 б1 и1 =—Х1 +
1 Б 1 ИБ
(61)
Запишем первое уравнение системы (29) с учетом (56), (57) в виде
ТИ2
^х" = ох1 - БЦ+.
(62)
Из (61), (62) приходим к следующему неоднородному дифференциальному уравнению второго порядка:
х1'-| и I Х1 = а,
(63)
3 (Б2 - Б12)
2 - - - / п и 3бю
где я2 = -> 0; А =--Щ-
ТБ ТБИ
Запишем общее решение (63):
^ -ь ~ ь бю
Х1 = с1е И + С2еИ -- ^
И (Б2 - Б12)'
(64)
Из граничного условия (36) и выражения (57) находим связь постоянных интегрирования
- 22ок
С1 = С2е И . Таким образом, из (64), (65) имеем
( Я (2а+Х) я Л
Х1 = С2
+ еИ
БЮ
И ( Б2 - Б12)'
(65)
(66)
а из (61) и (66) получаем распределение горизонтального перемещения на участке х1 е [-а;0):
и1 = С2
Б1И
Я(2а+х) Л
Бя
Ql Х1 (Б2 - 2 Б12) ИБ (Б2 - Б12)
+ С.
(67)
Рассмотрим решение на участке х1 е (0;«]. Запишем систему (30) с учетом выражений (55)-(60):
[х1' = К2х1 - К3щ [<+ = Ки+ + К1х1,
(68)
где К =
21
И80Б '
К1 = -
Т + Б1
Б
К 2 = -
3Б
ТИ21 1 +1.5^° ' И
К 3 =
3 (Т + Б1)
ТИ2
Систему линейных дифференциальных уравнений второго порядка (68) сведем к системе линейных дифференциальных уравнений первого порядка:
(69)
Характеристическое уравнение системы (69) запишем в виде
( и+Л (0 0 0 1 Л ( и+ +
Х1 0 0 1 0 х
w 0 К 2 0 -К 3 w
V У У ,К 0 К1 0 У V
Х4 + ЬХ2 + с = 0,
(70)
где Ь = К1К3 - К2 - К; с = КК2.
Коэффициенты биквадратного уравнения (70) представим в виде
Ь = -
т И8„
й38„1 1 +1,5-
т =
1 - 2 V 1 -V
6 (1 -V)
1,5
(1 -2v)(; +1,5И (1 -2V)(;-V) И
Дискриминант уравнения (70) имеет вид
— =
( и л2
V И80 У
(71)
где и2 = т2 -
24
8„
1 +1,5 ^Л И /г
На рис. 2 приведем зависимость п2 от отношения И/80 для различных значений коэффициента Пуассона.
Из приведенных зависимостей видно, что для выбранного коэффициента Пуассона существует пороговое значение отношения И/80, для которого п2 становится положительной величиной. При дальнейшем из-
6
с =
ложении будем рассматривать только такие толщины слоя и консоли, при которых и2 > 0 . Из (70) в этом случае получаем собственные значения:
X,, = ±.
X 3,4 =±.
где m1 = 0,5 (m + n) > 0; m2 = 0,5 (m - n) > 0.
Рис. 2. Зависимость n2 от отношения h/S0 : (1) v = 0,1; (2) v = 0,2; (5) v = 0,3
Fig. 2. Dependence of СС n2 on ratio h/50 : (1) v = 0,1; (2) v = 0,2; (5) v = 0,3
Запишем решение (69) в следующем виде:
U+=Z S/'\
j=1
4
=Z Sjqje%
(72)
(73)
j=1
где q = 2
X (1-v)-
(1 - 2v)
hS0Xj ,
При удовлетворении на бесконечности граничных условий (37), (38) получаем = = 0.
Найдем постоянные С2, С3 в (66), (67) и Б2, Б4 в (72), (73) из условий сопряжения решений (31)-(34):
,( 2ag Л
а
а
а
1 - e
V
D1h
D
( - 22ag Л 1 + e h v
+ C3 — S2 + S4,
D1Q
h (D2 - D12)
= S2 q2 + S4 q4,
h2 g
2ag Л С 1 - e h 1 = S,
„ . h3 50h q2 X 21 — + -T-
+S„
61
( ( q4X 4
h3 Snh
- + -3 2
v
2 Л Л + h
2\ Л
+ h
j j
+ (74)
h = S2 (DX2 - Dlq2) + S4 (DX4 - D1q4).
Запишем решение системы (74) в следующем виде: RR2 - VR3
S. = -
R4 - V1R3 - R1R2
S2 = V - V1S4, C2 = R + R1S4,
C3 — S2 + S4 C2
D1h
D
(
_ Л
1 - e h v j
где
V =
а
- 2ag
M = 1 + e h :
h(DX2 - D1q2)'
D1Q,
V1 =
B2 =
(DX4 - D1q4) ; (DX2 - D1q2) '
Vq2 + B2
h ( D2 - D12):
R = -
R1 =
q, - V 1q2 _
M
h sh
h2 g ( -— R 2 = £«.| 1 - e h
M R3 = q2 X 2
+ h; R4 = q4X41 — +
h3 5„ h
3 2 ) " \ 3 2
Рассмотрим упрощенный вариант решения рассмотренной задачи, положив в задании поля перемещений (8), (9) у1 = у2 = 0. В этом случае из (29)-(38) с учетом (58) приходим к следующей постановке задачи: для участка х1 е[-с;0)
й8п
dx,
■ = 0,
для участка x1 е (0; i]
dQL= 5
dx,
с условиями сопряжения
QnL„ = 6111
1 =-0 lx1 =+0
и граничным условием в левом торце
QiL = 6
и правом торце
- г I 2u+
где ст„ = L | —L
8„
6111
= 0,
(75)
(76)
(77)
(78)
(79)
(80)
а (x, )=hDdf.
В постановке (75)-(80) связь (76) соответствует постановке «дифференциального сдвига», предложенной в [29] и обсуждаемой в [30].
Решение (75)-(80) запишем в виде
u+ = 61 u+ = - 61
Dh М 2LDh
x1 е [-a;0),
Sq _e^dh50 x1 x
2LDh
1 е[0; i).
j-1
x =i
Введем в рассмотрение энергетическое произведение (ЭП) [27] среднее на элементе слоя размером
80 х80 :
2 у = 80 ф,
1 80
где ф = — |ф(х1) —х1 ; ф(х1) - средняя удельная
80 0
бодная энергия по толщине слоя. Из (81) находим
(82)
сво-
2 у = -
(ß. )2
1
2 LDh5„
i I2LS,
-2,--
, 4 Dh - 1
V
Л
(83)
/
При 5„/h ^ 0 выражение (83) принимает постоянное значение
2% =
м
Dh
(84)
Рассмотрим консоль со следующими геометрическими и механическими характеристиками: а = 0,1 м,
И = 0,05 м, V = 0,3, Е = 2-10" Па. На рис. 3 показано поведение функции у/у„ в зависимости от десятичного логарифма отношения 80/И для единичной внешней нагрузки.
X %
1
0,95 0,9 0,85 0,8 0,75
1 2
-8
-7
-6
-5
3 ig(S„//0
Рис. 3. Зависимость величины ЭП, отнесенной к значению (84), от отношения 80/h : 1 - ЭП рассчитано по формуле (83);
2 - ЭП рассчитано на основе формул (72) и (73)
Fig. 3. Dependence of the EP magnitude, referred to the value (84), on the ratio 80/h : 1 - EP calculated by the formula (83);
2 - EP calculated on the basis of formulas (72) and (73)
На рис. 3 график 1 определяет соотношение (83), а график 2 построен на основе решений (72), (73). Из полученных зависимостей видно, что имеет место сходимость ЭП. При отношении 50/h < 10-6 значения ЭП практически не меняются. Однако учет обжатия консоли дает снижение значения ЭП на 20 %.
Объяснение данному факту можно дать, рассмотрев
, _c u+c h dш. зависимости деформации в слое е12 =--1---
50 2 dx.
с учетом обжатия и без него; е^ = —, где и+с - гори-
80
зонтальное перемещение в модели с обжатием. На рис. 4 построены деформации ёР, р = с, w,;с, 2с в слое с учетом и без учета обжатия консоли, а также слагаемые деформации с учетом обжатия на участке х1 е [0, а/4] при толщине слоя, равной 80 = 10-2 И . График 1 определяет деформацию е^ , график 2 - деформацию ёс2, график 3 - сла-
гаемое деформации е.2с = , график 4 - слагаемое е,
5
-2с 12
деформации е12с = И—Х. Значения деформаций отнесе-
2 —х1
ны к значению деформации ё^ на торце слоя.
Рис. 4. Зависимость деформаций в слое: 1 - е';
Fig. 4. Dependence of strains in the layer: 1 - е.' ;
Из рис. 4 видно, что деформации s.2/ имеют противоположный знак и составляют порядка 5 % от деформации е.' на торце слоя. В итоге деформация в слое ё12 с учетом обжатия консоли составляет около 93 % от деформации е.'.
При толщине слоя 50 = Ю-6 h характер решения меняется. Деформация s22c составляет порядка 0,03 % от деформации е.', и ее можно не принимать во внимание. Однако в этом случае деформация ё^ будет равна 89 % от деформации е.', что приводит к результату, показанному на рис. 3 для ЭП.
Заключение
Рассмотрено НДС ДКБ-образца, близкое к сдвиговому, в рамках модели слоя взаимодействия. Гипотеза об однородном распределении главных деформаций по толщине консолей позволила получить аналитическое решение задачи. Установлено, что ЭП при стремлении толщины слоя к нулю на некотором отрезке, начиная с порогового значения, практически не зависит от тол-
+
+ с
щины слоя. Определено, что учет обжатия консолей на 20 % уменьшает пороговое значение ЭП. Если в результате эксперимента определить критическое значение внешней нагрузки при заданной толщине консоли, то можно найти критическое пороговое значение ЭП. Универсальность данного значения можно установить, сравнивая величины критических нагрузок, определенных теоретически из условия постоянства ЭП при различных толщинах консолей с результатами экспериментов.
Библиографический список
1. Broberg K.B. Cracks and fracture. - London: Academic Press, UK, 1999. - 752 p.
2. Черепанов Г.П. Механика хрупкого разрушения. - М.: Наука, 1974. - 640 с.
3. Thouless M.D., Hutchinson J.W., Liniger E.G. Plane-strain, buckling-driven delamination of thin films: Model experiments and mode-II fracture // Acta Mater. - 1992. - Vol. 40. - No. 10. -Р. 2639-2649. https://doi.org/10.1016/0956-7151(92)90333-A
4. Samudrala O., Huang Y., Rosakis A.J. Subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone // Journal of the Mechanics and Physics of Solids. -2002. - Vol. 50. - No. 6. - Р. 1231-1268. https://doi.org/10.1016/S0022-5096(01)00129-6
5. Fang X., Charalambides P.G. The fracture mechanics of cantilever beams with an embedded sharp crack under end force loading // Engineering Fracture Mechanics. - 2015. - Vol. 149. -Р. 1-17. https://doi.org/10.1016/j.engfracmech.2015.09.039
6. Fang X., Charalambides P.G. A J-integral approach in characterizing the mechanics of a horizontal crack embedded in a cantilever beam under an end transverse force // Engineering Fracture Mechanics. - 2017. - Vol. 169. - Р. 35-53. https://doi.org/10.1016/j.engfracmech.2016.11.010
7. Davidson B.D., Hu H., Schapery R.A. An analytical crack-tip element for layered elastic structures // Journal of Applied Mechanics. - 1995. - Vol. 62. - No. 2. - Р. 294-305. DOI:10.1115/1.2895931
8. Bruno D., Greco F., Lonetti P. Computation of energy release rate and mode separation in delaminated composite plates by using plate and interface variables // Mechanics of Advanced Materials and Structures. - 2005. - Vol. 12. - No. 4. - Р. 285-304. https://doi.org/10.1080/15376490590953563
9. Yarasca J., Mantari J.L., Arciniega R.A. Hermite-Lagran-gian finite element formulation to study functionally graded sandwich beams // Composite Structures. - 2016. - Vol. 140. - No. 15. -Р. 567-581. https://doi.org/10.1016/j.compstruct.2016.01.015
10. Todo M., Jar P.-Y.B., Takahashi K. Initiation of a mode-II interlaminar crack from an insert film in the end-notched flexure composite specimen // Composites Science and Technology. -2000. - Vol. 60. - No. 2. - Р. 263-272. https://doi.org/10.1016/S0266-3538(99)00119-0
11. Shahverdi M., Vassilopoulos A.P., Keller T. Mixed-Mode I/II fracture behavior of asymmetric composite joints // Procedia Structural Integrity. - 2016. - Vol. 2. - Р. 1886-1893. https://doi.org/10.1016/j.prostr.2016.06.237
12. Morais A.B. Analysis of the metal adhesively bonded double cantilever beam specimen // International Journal of Adhesion and Adhesives. - 2015. - Vol. 61. - Р. 8-14. https://doi.org/10.1016/j.ijadhadh.2015.04.010
13. Тимошенко С.П., Войновский-Кригер С. Пластины и оболочки - М.: Физматгиз, 1963. - 636 с.
Благодарность
Работа выполнена при финансовой поддержке РФФИ, проекты № 18-31-20053 и № 19-41-710001 р_а.
Acknowledgment
The reported study was funded by RFBR according to the research projects No 18-31-20053 and No 19-41-710001 г_а.
14. Williams J.G. On the calculation of energy release rates for cracked laminates // International Journal of Fracture. -1988. - Vol. 36. - No. 2. - Р. 101-119. https://doi.org/10.1007/BF00017790
15. Suo Z., Hutchinson J.W. Interface crack between two elastic layers // International Journal of Fracture. - 1990. -Vol. 43. - No. 1. - Р. 1-18. https://doi.org/10.1007/BF00018123
16. Mantari J.L., Yarasca J. A simple and accurate generalized shear deformation theory for beams // Composite Structures. - 2015. - Vol. 134. - Р. 593-601. https://doi.org/10.1016/j.compstruct.2015.08.073
17. Joseph R.P., Wang B.L., Samali B. Size effects on double cantilever beam fracture mechanics specimen based on strain gradient theory // Engineering Fracture Mechanics. - 2017. - Vol. 169. -Р. 309-320. https://doi.org/10.1016/j .engfracmech.2016.10.013
18. Berto F., Glagolev V.V., Markin A.A. A body failure model with a notch based on the scalable linear parameter // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 4. -С. 93-97. D0I:10.15593/perm.mech/2018.4.08
19. Гаврилкина М.В., Глаголев В.В., Маркин А.А. К решению одной задачи механики разрушения // Прикладная механика и техническая физика. - 2007. - Т. 48, № 4. -С. 121-127.
20. Глаголев В.В., Маркин А.А. Модель трещины с масштабируемым линейным параметром // Актуальные проблемы прикладной математики, информатики и механики: сб. тр. междунар. науч.-техн. конф. - Воронеж, 2018. -С. 1082-1087.
21. Prandtl L. Ein Gedankenmodell für den Zerreißvorgang spröder Körper // Zeitschrift für Angewandte Mathematik und Mechanik. - 1933. - Vol. 13. - S. 129-133. = Prandtl L., Knauss W.G. A thought model for the fracture of brittle solids // International Journal of Fracture. - 2011. - Vol. 171. - No. 2. -P. 105-109. https://doi.org/10.1007/s10704-011-9637-3
22. Ентов В.М., Салганик Р.Л. К модели хрупкого разрушения Прандтля // Изв. АН СССР. МТТ. - 1968. - № 6. -С. 87-99.
23. Салганик Р.Л., Мищенко А.А., Федотов А.А. Модель трещины Прандтля и ее применение для решения задачи механики контактного взаимодействия // К 75-летию со дня рождения профессора Владимира Марковича Ентова; Ин-т ком-пьют. исслед. - Ижевск, 2012. - 180 с.
24. Салганик Р.Л., Мищенко А.А., Федотов А.А. Напряженное состояние в окрестности выработки, пройденной в глубокозалегающем горизонтальном пласте // Физико-технические проблемы разработки полезных ископаемых. -2015. - № 2. - С. 24-33.
25. Макклинток Ф. Пластические аспекты разрушения // Разрушение. Т. 3 - М.: Мир, 1976. - С. 67-262.
26. Глаголев В.В., Глаголев Л.В., Маркин А.А. Определение напряженно-деформированного состояния упруго-пластических тел с боковым трещиноподобным дефектом с использованием модели с линейным размером // Прикладная механика и техническая физика. - 2018. - Т. 59, № 6. -С. 143-154.
27. Glagolev V.V., Markin A.A. Fracture models for solid bodies, based on a linear scale parameter // International Journal of Solids and Structures. - 2019. - Vol. 158. - P. 141-149. D0I.org/10.1016/j.ijsolstr.2018.09.002
References
1. Broberg K.B. Cracks and fracture. 1999. London, UK, Academic Press, 752 p.
2. Cherepanov G.P. Mekhanika hrupkogo razrusheniya [Mechanics of brittle failure]. Moscow, Nauka, 1974, 640 p.
3. Thouless M.D., Hutchinson J.W., Liniger E.G. Plane-strain, buckling-driven delamination of thin films: Model experiments and mode-II fracture. Acta Mater., 1992, vol. 40, no. 10, pp. 2639-2649. DOI: 10.1016/0956-7151(92)90333-A
4. Samudrala O., Huang Y., Rosakis A.J. Subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone. Journal of the Mechanics and Physics of Solids, 2002, vol. 50, no. 6, pp. 1231-1268. DOI: 10.1016/S0022-5096(01 )00129-6
5. Fang X., Charalambides P.G. The fracture mechanics of cantilever beams with an embedded sharp crack under end force loading. Engineering Fracture Mechanics, 2015, vol. 149, pp. 1-17. DOI: 10.1016/j.engfracmech.2015.09.039
6. Fang X., Charalambides P.G. A J-integral approach in characterizing the mechanics of a horizontal crack embedded in a cantilever beam under an end transverse force. Engineering Fracture Mechanics, 2017, Vol. 169, pp. 35-53. DOI: 10.1016/j.engfracmech.2016.11.010
7. Davidson B.D., Hu H., and Schapery R.A. An analytical crack-tip element for layered elastic structures. Journal of Applied Mechanics, 1995, vol. 62, no. 2, pp. 294-305. DOI:10.1115/1.2895931
8. Bruno D., Greco F., Lonetti P. Computation of energy release rate and mode separation in delaminated composite plates by using plate and interface variables. Mechanics of Advanced Materials and Structures, 2005, vol. 12, no. 4, pp. 285-304. DOI: 10.1080/15376490590953563
9. Yarasca J., Mantari J.L., Arciniega R.A. Hermite-Lagran-gian finite element formulation to study functionally graded sandwich beams. Composite Structures, 2016, vol. 140, no. 15, pp. 567-581. DOI: 10.1016/j.compstruct.2016.01.015
10. Todo M., Jar P.-Y.B., Takahashi K. Initiation of a mode-II interlaminar crack from an insert film in the end-notched flexure composite specimen. Composites Science and Technology, 2000, vol. 60, no. 2, pp. 263-272. DOI: 10.1016/S0266-3538(99)00119-0
11. Shahverdi M., Vassilopoulos A.P., Keller T. Mixed-Mode I/II fracture behavior of asymmetric composite joints. Procedia Structural Integrity, 2016, vol. 2, pp. 1886-1893. DOI: 10.1016/j.prostr.2016.06.237
12. Morais A.B. Analysis of the metal adhesively bonded double cantilever beam specimen. International Journal of Adhe-
28. Glagolev, V.V., Glagolev, L.V., Markin A.A. Stress-strain state of elastoplastic bodies with crack // Acta Mechanica Solida Sinica. - 2015. - Vol. 28. - No. 4. - P. 375-383. DOI:10.1016/S0894-9166(15)30023-9
29. Volkersen O. Die Nietkraftverteilung in zugbeanspruchten Nietverbindungen mit konstanten Laschenquerschnitten // Luftfahrtforschung. - 1938. - Vol. 15. - P. 41-47.
30. Analytical models of adhesively bonded joints. Part I: Literature survey / L.F.M. Da Silva, P.J.C. das Neves, R.D. Adams, J.K. Spelt // Int. J. Adhes. - 2009. - Vol. 29. - No. 3. - P. 319330. https://doi.org/10.1016/jijadhadh.2008.06.005
sion and Adhesives, 2015, vol. 61, pp. 8-14. DOI: 10.1016/j.ijadhadh.2015.04.010
13. Timoshenko S.P., Vojnovskij-Kriger S. Plastiny i obo-lochki [Plates and shells]. Moscow, Fizmatgiz,1963, 636 p.
14 Williams J.G. On the calculation of energy release rates for cracked laminates. International Journal of Fracture, 1988, vol. 36, no. 2, pp. 101-119. DOI: 10.1007/BF00017790
15. Suo Z., Hutchinson J.W. Interface crack between two elastic layers. International Journal of Fracture, 1990, vol. 43, no. 1, pp. 1-18. DOI: 10.1007/BF00018123
16. Mantari J.L., Yarasca J. A simple and accurate generalized shear deformation theory for beams. Composite Structures, 2015, vol. 134, pp. 593-601. DOI: 10.1016/j .compstruct.2015.08.073
17. Joseph R.P., Wang B.L., Samali B. Size effects on double cantilever beam fracture mechanics specimen based on strain gradient theory. Engineering Fracture Mechanics, 2017, vol. 169, pp. 309-320. DOI: 10.1016/j.engfracmech.2016.10.013
18. Berto F., Glagolev V.V., Markin A.A. A body failure model with a notch based on the scalable linear parameter. PNRPU Mechanics Bulletin, 2018, no. 4, pp. 93-97. DOI:10.15593/perm.mech/2018.4.08
19. Gavrilkina M.V., Glagolev V.V., Markin A.A. Solution of one problem of fracture mechanics. Journal of Applied Mechanics and Technical Physics, 2007, vol. 48, no 4, pp. 571-576. DOI: 10.1007/s10808-007-0072-1
20. Glagolev V.V., Markin A.A. Model' treshchiny s masshtabiruemym linejnym parametrom [Crack model with scalable linear parameter]. Aktual'nye problemy prikladnoj matematiki, informatiki i mekhaniki: sbornik trudov Mezhdunarodnoj nauchno-tekhnicheskoj konferencii (Actual problems of applied mathematics, computer science and mechanics: a collection of works of the International Scientific and Technical Conference). Voronezh, 2018, pp. 1082-1087.
21. Prandtl L., Knauss W.G. A thought model for the fracture of brittle solids. International Journal of Fracture, 2011, vol. 171, no. 2, pp. 105-109. DOI: 10.1007/s10704-011-9637-3
22. Entov V.M., Salganik R.L. K modeli hrupkogo razrusheniya Prandtlya [To the Prandtl model of brittle fracture]. Izv. AN SSSR. MTT, 1968, no. 6, pp. 87-99.
23. Salganik R.L., Mishchenko A.A., Fedotov A.A. Model' treshchiny Prandtlya i ee primenenie dlya resheniya zadachi mekhaniki kontaktnogo vzaimodejstviya [Prandtl crack model and its application for solving the problem of contact interaction mechanics]. Izhevsk: Institut komp'yuternyh issledovanij, 2012, 180 p.
24. Salganik R.L., Mishchenko A.A., Fedotov A.A. Stress state in the vicinity of excavation in deep horizontal bed. Journal of Mining Science, 2015, vol. 51, no. 2, pp. 220-227. DOI: 10.1134/S1062739115020027
25. Makklintok F. Plasticheskie aspekty razrusheniya [Plastic aspects of destruction]. Razrushenie T. 3, Moscow, Mir, 1976, pp. 67-262.
26. Glagolev V.V., Glagolev L.V., Markin A.A. Determining the stress-strain state of elastic-plastic solids with a lateral crack-like defect with the use of a model with a linear size. Journal of Applied Mechanics and Technical Physics, 2018, vol. 59, no. 6, pp. 1085-1094.
27. Glagolev V.V., Markin A.A. Fracture models for solid bodies, based on a linear scale parameter. International Journal of
Solids and Structures, 2019, vol. 158, pp. 141-149. DOI: 10.1016/j.ijsolstr.2018.09.002
28. Glagolev, V.V., Glagolev, L.V., Markin A.A. Stress-strain state of elastoplastic bodies with crack. Acta Mechanica Solida Sinica, 2015, vol. 28, no. 4, pp. 375-383. DOI: 10.1016/S0894-9166( 15)30023-9
29. Volkersen, O. Die Nietkraftverteilung in zugbeanspruchten Nietverbindungen mit konstanten Laschenquerschnitten. Luftfahrtforschung,1938, vol. 15, pp. 41-47.
30. Da Silva L.F.M., das Neves P.J.C., Adams R.D., Spelt J.K. Analytical models of adhesively bonded joints Part I: Literature survey. Int. J. Adhes. Adhes. 2009, vol. 29, no. 3, pp. 319-330. DOI: 10.1016/j.ijadhadh.2008.06.005