Energy product in a crack-like defect model under loading of mode II type

Abstract


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.

Full Text

Введение Классическая постановка задачи о трещине моды II рассматривает бесконечное линейно-упругое пространство с дефектом в виде математического разреза, поверхности которого нагружены противоположными векторами напряжений [1, 2]. В этом случае на продолжении математического разреза в слое нулевой толщины реализуются только сдвиговые напряжения. Для образцов конечных размеров предлагается ряд схем нагружения [3-6], в которых на продолжении трещины реализуется напряженное состояние, соответствующее поперечному сдвигу. Данные схемы также принято рассматривать как трещину моды II. Отметим, что для образцов конечных размеров общие решения строятся, как правило, на основе численных процедур, основанных на методе конечных элементов [7-12]. Аналитические решения получаются, как правило, в рамках упрощающих гипотез [13-19]. Так, в работе [6] для консоли с центральной трещиной на основе гипотез теории пластин Тимошенко получено решение задачи и определено значение J-интеграла. В этом случае в консолях учитывались деформации растяжения и сдвига. Влияние обжатия консолей в результате внешнего воздействия не рассматривалось. Альтернативной моделью трещины в виде математического разреза является разрез физический с характерной толщиной . В данном представлении особую роль играет описание взаимодействия материального объема, лежащего на продолжении физического разреза, со смежной средой. Важно, чтобы та или иная модель не приводила к сингулярности в области окончания физического разреза. При нагружении трещиноподобного дефекта нормальным отрывом отметим модель в виде связей Прандтля, предложенную в статье [20] и используемую в работах [21-23]. В модели Ф. Макклинтока [24] предложено использовать для материала слоя усредненные характеристики напряженно-деформированного состояния (НДС). Постановки конкретных задач исходя из представления Ф. Макклинтока были реализованы на основе теории слоя взаимодействия в работах [18-20], [26-28]. В частности, в статье [27] введено понятие энергетического произведения (ЭП) в виде произведения приращения удельной свободной энергии в элементе слоя и толщины слоя. В статье [18] о нагружении двухконсольной балки (ДКБ-образца) нормальным отрывом (нагружение по моде I) показано, что при уменьшении отношения толщины слоя к высоте консоли имеет место сходимость ЭП к поверхностной энергии, рассчитанной через вязкость разрушения . В данной работе предложена схема нагружения ДКБ-образца в которой реализуется состояние, близкое к моде II. На основе общей вариационной постановки задачи [26-28] получены упрощенные постановки с учетом и без учета обжатия консолей. Из аналитических решений показано, что при уменьшении отношения толщины слоя к толщине консоли имеет место сходимость ЭП в слое для данного вида нагружения. 1. Постановка задачи Рассмотрим случай нагружения ДКБ-образца согласно схеме, показанной на рис. 1. Правые торцы консолей 1 и 2 образца свободны от внешней нагрузки. Верхняя граница консоли 1 и нижняя граница консоли 2 закреплены от вертикальных перемещений. На левых торцах консолей действует горизонтальная антисимметричная нагрузка интенсивностью . Рис. 1. Схема сдвигового воздействия на тонкий слой Fig. 1. Shear effect on a thin layer Для описания взаимодействия слоя 3 с телами 1 и 2 применим концепцию «слоя взаимодействия», развитую в работах [26-28]. В этом случае условия равновесия тел 1 и 2 запишем в вариационной форме для тела 1: (1) и тела 2 (2) где - вектор перемещений; - вектора внешней нагрузки тел 1 и 2; , - площади тел 1 и 2; , - тензоры напряжений и деформаций; , - тензоры средних напряжений и деформаций слоя с компонентами: (3) (4) где - компоненты векторов перемещений верхней и нижней границ слоя соответственно; ; - граница приложения внешней нагрузки для тела 1 и 2. Постулируется жесткое сцепление между границами области 3 и областями 1, 2: (5) Примем определяющие соотношения в форме закона Гука: (6) где , - модуль упругости и коэффициент Пуассона тела; - объемное расширение, ; - символ Кронекера; . Для материала слоя взаимодействия 3 определяющие соотношения считаем справедливыми для средних компонент тензоров напряжений и деформаций: (7) Распределение перемещений в консолях при данном нагружении принимаем в следующем виде: (8) (9) Вертикальные перемещения тождественно удовлетворяют граничным условиям на верхней границе консоли 1 и нижней границе консоли 2: В этом случае деформацию консолей будут определять компоненты тензора деформаций (10) (11) (12) где - горизонтальное перемещение верхней (+) или нижней (-) границы слоя; символ преобразуется в знак (-) для , а для преобразуется в знак (+). При распределении деформаций в консолях в виде (10)-(12) работа внутренних напряжений в теле 1 рассчитывается следующим образом: (13) в теле 2: (14) Введем в рассмотрение обобщенные силы: (15) (16) (17) (18) и обобщенные моменты: (19) (20) Интегрируем по частям ряд слагаемых в (13), (14) с учетом (15)-(20) и четвертое и пятое слагаемые в левой части (1) и (2) с учетом (8) и (9): (21) (22) (23) (24) (25) (26) Рассмотрим правые части уравнений (1) и (2). При положительном обходе контура работа внешних напряжений рассчитывается так: (27) (28) где Подставив (13),(15), (16) (21)-(28) в (1) и приравняв слагаемые при равных вариациях, приходим к двум системам дифференциальных уравнений: для участка (29) для участка (30) с условиями сопряжения (31) (32) (33) (34) и с естественными граничными условиями: на левом торце (35) (36) на правом торце (37) (38) Аналогичным образом из (2) приходим к системам дифференциальных уравнений: для участка (39) для участка (40) с условиями сопряжения (41) (42) (43) (44) и естественными граничными условиями: на левом торце (45) (46) правом торце (47) (48) С учетом выражений (10)-(12) и условием плоской деформации ( ) определяющие соотношения (6) запишем в виде (49) (50) (51) где Из (7)-(9), (3) и (4) получим связь средних напряжений в слое взаимодействия с его граничными перемещениями: (52) (53) (54) Таким образом, система дифференциальных уравнений (29), (30), (39), (40) с учетом связей (10)-(12), (15)-(20) и определяющих соотношений (49)-(54) является замкнутой относительно четырех неизвестных функций . Каждая искомая функция входит в систему со второй производной включительно. Удовлетворение восьми граничных условий (35)-(38), (45)-(48) с учетом сопряжения решений (31)-(34), (41)-(44) даст решение поставленной задачи. 2. Решение задачи Рассмотрим решение поставленной задачи в виде для полубесконечного образца ( ). В этом случае движение консолей противоположно, с одинаковым модулем обжатия. Для решения задачи достаточно рассмотреть только систему уравнений (29), (30) при следующих выражениях обобщенных сил: (55) (56) обобщенного момента (57) и средних напряжений в слое (58) (59) (60) Проинтегрируем второе уравнение системы (29) с учетом граничного условия (35) и связи (55): (61) Запишем первое уравнение системы (29) с учетом (56), (57) в виде (62) Из (61), (62) приходим к следующему неоднородному дифференциальному уравнению второго порядка: (63) где Запишем общее решение (63): (64) Из граничного условия (36) и выражения (57) находим связь постоянных интегрирования (65) Таким образом, из (64), (65) имеем (66) а из (61) и (66) получаем распределение горизонтального перемещения на участке : (67) Рассмотрим решение на участке . Запишем систему (30) с учетом выражений (55)-(60): (68) где ; ; ; . Систему линейных дифференциальных уравнений второго порядка (68) сведем к системе линейных дифференциальных уравнений первого порядка: (69) Характеристическое уравнение системы (69) запишем в виде (70) где Коэффициенты биквадратного уравнения (70) представим в виде Дискриминант уравнения (70) имеет вид (71) где На рис. 2 приведем зависимость от отношения для различных значений коэффициента Пуассона. Из приведенных зависимостей видно, что для выбранного коэффициента Пуассона существует пороговое значение отношения для которого становится положительной величиной. При дальнейшем изложении будем рассматривать только такие толщины слоя и консоли, при которых . Из (70) в этом случае получаем собственные значения: где Рис. 2. Зависимость от отношения : (1) ; (2) ; (3) Fig. 2. Dependence of СС on ratio : (1) ; (2) ; (3) Запишем решение (69) в следующем виде: (72) (73) где При удовлетворении на бесконечности граничных условий (37), (38) получаем . Найдем постоянные в (66), (67) и в (72), (73) из условий сопряжения решений (31)-(34): (74) Запишем решение системы (74) в следующем виде: где ; ; . Рассмотрим упрощенный вариант решения рассмотренной задачи, положив в задании поля перемещений (8), (9) . В этом случае из (29)-(38) с учетом (58) приходим к следующей постановке задачи: для участка (75) для участка (76) с условиями сопряжения (77) (78) и граничным условием в левом торце (79) и правом торце (80) где . В постановке (75)-(80) связь (76) соответствует постановке «дифференциального сдвига», предложенной в [29] и обсуждаемой в [30]. Решение (75)-(80) запишем в виде (81) Введем в рассмотрение энергетическое произведение (ЭП) [27] среднее на элементе слоя размером : (82) где ; - средняя удельная свободная энергия по толщине слоя. Из (81) находим (83) При выражение (83) принимает постоянное значение (84) Рассмотрим консоль со следующими геометрическими и механическими характеристиками: м, м, , Пa. На рис. 3 показано поведение функции в зависимости от десятичного логарифма отношения для единичной внешней нагрузки. Рис. 3. Зависимость величины ЭП, отнесенной к значению (84), от отношения : 1 - ЭП рассчитано по формуле (83); 2 - ЭП рассчитано на основе формул (72) и (73) Fig. 3. Dependence of the EP magnitude, referred to the value (84), on the ratio : 1 - EP calculated by the formula (83); 2 - EP calculated on the basis of formulas (72) and (73) На рис. 3 график 1 определяет соотношение (83), а график 2 построен на основе решений (72), (73). Из полученных зависимостей видно, что имеет место сходимость ЭП. При отношении значения ЭП практически не меняются. Однако учет обжатия консоли дает снижение значения ЭП на 20 %. Объяснение данному факту можно дать, рассмотрев зависимости деформации в слое с учетом обжатия и без него; , где - горизонтальное перемещение в модели с обжатием. На рис. 4 построены деформации в слое с учетом и без учета обжатия консоли, а также слагаемые деформации с учетом обжатия на участке при толщине слоя, равной . График 1 определяет деформацию , график 2 - деформацию , график 3 - слагаемое деформации , график 4 - слагаемое деформации . Значения деформаций отнесены к значению деформации на торце слоя. Рис. 4. Зависимость деформаций в слое: 1 - 2 - ; 3 - ; 4 - Fig. 4. Dependence of strains in the layer: 1 - ; 2 - ; 3 - ; 4 - Из рис. 4 видно, что деформации имеют противоположный знак и составляют порядка 5 % от деформации на торце слоя. В итоге деформация в слое с учетом обжатия консоли составляет около 93 % от деформации . При толщине слоя характер решения меняется. Деформация составляет порядка 0,03 % от деформации , и ее можно не принимать во внимание. Однако в этом случае деформация будет равна 89 % от деформации , что приводит к результату, показанному на рис. 3 для ЭП. Заключение Рассмотрено НДС ДКБ-образца, близкое к сдвиговому, в рамках модели слоя взаимодействия. Гипотеза об однородном распределении главных деформаций по толщине консолей позволила получить аналитическое решение задачи. Установлено, что ЭП при стремлении толщины слоя к нулю на некотором отрезке, начиная с порогового значения, практически не зависит от толщины слоя. Определено, что учет обжатия консолей на 20 % уменьшает пороговое значение ЭП. Если в результате эксперимента определить критическое значение внешней нагрузки при заданной толщине консоли, то можно найти критическое пороговое значение ЭП. Универсальность данного значения можно установить, сравнивая величины критических нагрузок, определенных теоретически из условия постоянства ЭП при различных толщинах консолей с результатами экспериментов.

About the authors

V V Glagolev

Tula State University

L V Glagolev

Tula State University

A A Markin

Tula State University

References

  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-Lagrangian 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 с.
  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. doi: 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. DOI.org/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. - 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. 319-330. https://doi.org/10.1016/j.ijadhadh.2008.06.005

Statistics

Views

Abstract - 365

PDF (Russian) - 214

Cited-By


PlumX


Copyright (c) 2019 Glagolev V.V., Glagolev L.V., Markin A.A.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies