The deformation model of the DCB-sample with elastoplastic properties
- Authors: Glagolev V.V1, Markin A.A1
- Affiliations:
- Tula State University
- Issue: No 2 (2021)
- Pages: 55-63
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/1773
- DOI: https://doi.org/10.15593/perm.mech/2021.2.06
- Cite item
Abstract
The loading of a strip with a crack-like defect according to mode I is considered. In contrast to the classical representation of a crack in the form of a mathematical section, the proposed model defines a crack as a physical cut with a characteristic linear size. The mental continuation of a physical cut in a solid forms an interaction layer (IL). It is important that the stress-strain state of the layer at a finite value of the linear parameter does not introduce a singularity into the crack model. The process of elastoplastic deformation with a constant layer length is considered. We obtained a simplified analytical solution to the problem of deformation of two elastic bodies connected by a thin layer with elastoplastic properties. The dependence of the displacement and stress fields on the length and thickness of the interaction layer has been found. It is shown that, under the classical plasticity condition, the range of variation of the external load leading to a purely elastic behavior is possible only for a finite layer thickness. As the layer thickness tends to zero, as in the Dugdale model, the plasticity region is formed at an arbitrarily small external load. For small layer thicknesses, a local plasticity criterion is proposed, by using which it is possible to distinguish the intervals of the external load variations associated with elastic and plastic deformations. The local plasticity condition, determined by the critical value of the energy product, makes it possible to reflect the stage of elastic deformation at an arbitrarily small finite thickness of the interaction layer. An asymptotic dependence of the external load on the IL thickness and the reduced length of the plastic zone is obtained. At the same time, the separation of the external load into elastic and plastic components is preserved. From the analysis of the experimental data, an estimate of the elastic limit of the energy product for the AV138 adhesive was obtained.
Full Text
Введение Экспериментальное исследование трещиностойкости материалов использует в качестве образца для изучения трещины нормального отрыва двухконсольную балку (ДКБ-образец). На образцы и порядок проведения испытаний действует соответствующий регламент, определяемый по ГОСТу. В процессе испытаний определяется критическое значение коэффициента интенсивности напряжений (КИН) - вязкость разрушения [1-3]. При этом используется определение КИН в рамках классической модели трещины, рассматриваемой в виде математического разреза в линейно-упругой среде. Непосредственно измеряемой характеристикой данного эксперимента является зависимость расклинивающего усилия от перемещения, в которой выделяются линейный и нелинейный участки. Для выяснения причины нелинейности необходимо проводить промежуточные разгрузки с целью выявления наличия пластических деформаций либо их отсутствия. Учет возможной физической нелинейности, обусловленной пластичностью, находит отражение в моделях трещины нулевой толщины путем введения поправки Ирвина - Орована [4, 5], сил сцепления [6-12], тонкой пластической зоны [13, 14], пластического коэффициента интенсивности напряжений [15], поперечника пластической области [16, 17]. Как правило, образование пластической зоны в этих моделях происходит с начала нагружения, а стадия чисто упругого деформирования не отражается. Классическая постановка упругопластической задачи [18-21] подразумевает наличие условия перехода материальной области из упругого состояния в пластическое. Если область не содержит источников сингулярности напряжений, то из решения упругой задачи определяется внешняя нагрузка, при которой достигается предел упругости. В этом случае в рассмотрение вводится модель трещины в виде физического разреза [22]. Модель упругого деформирования трещины в виде физического разреза может быть рассмотрена на основе связей Прандтля [23-26], однако вопрос о выборе характерного размера и жесткости связей открыт. В работах [27-31] рассмотрена модель, основанная на концепции слоя взаимодействия, для которой материальные свойства слоя соответствуют поврежденному телу, а характерный размер определяется из решения обратной задачи [30]. В данной работе на основе общей вариационной постановки [30] получено упрощенное аналитическое решение задачи деформирования двух упругих тел, связанных тонким слоем с упругопластическими свойствами. Рассматривается процесс упругопластического деформирования при неизменной длине слоя. Полученное решение в упругой области не содержит осцилляций, в отличие от решения в балочном приближении [28]. Показано, что при классическом условии текучести диапазон внешней нагрузки, приводящей к чисто упругому поведению, возможен только при конечной толщине слоя. При стремлении толщины слоя к нулю, как и в задаче Дагдейла, область пластичности образуется при сколь угодно малой внешней нагрузке. Для малых толщин слоя, вплоть до нулевых, предложен локальный критерий пластичности, при использовании которого возможно выделение внешних нагрузок, связанных с упругим и пластическим деформированием. Из анализа экспериментальных данных [32] приведена оценка зоны пластического деформирования при использовании локального критерия. Постановка задачи нормального отрыва На рисунке показана двухконсольная балка (ДКБ-образец) , состоящая из трех тел. Пластины 1 и 2 с одинаковыми толщинами по длине связаны слоем взаимодействия 3 толщиной . Материал пластин принимаем линейно-упругим, а материал слоя взаимодействия - идеально упругопластическим. Пластические деформации слоя ограничены длиной . Правый торец образца жестко закреплен от горизонтальных и вертикальных перемещений, на левых торцах консолей действует вертикальная симметричная нагрузка интенсивностью . Вся остальная поверхность образца свободна от внешней нагрузки. Рис. Модель ДКБ-образца Fig. Model of the DCB-sample Для описания взаимодействия слоя 3 с телами 1 и 2 применим концепцию «слоя взаимодействия», развитую в работах [24-29]. В этом случае условия равновесие тел 1 и 2 запишем в вариационной форме: - для тела 1: (1) - тела 2: (2) где , - площади тел 1 и 2; , - тензоры напряжений и деформаций; , - тензоры средних напряжений и деформаций слоя с компонентами: , , , (3) (4) где - компоненты векторов перемещений верхней и нижней границ слоя соответственно; ; - граница приложения внешней нагрузки для тела 1 и 2. Постулируется жесткое сцепление между границами области 3 и областями 1, 2: Для материала пластин примем определяющие соотношения в форме закона Гука: (5) где , - модуль упругости и коэффициент Пуассона тела; - объемная деформация; - символ Кронекера; . Предполагаем, что напряженное состояние слоя для данного вида нагружения определятся одной компонентой тензора средних напряжений. Для материала слоя взаимодействия 3 определяющие соотношения на стадии обратимого деформирования принимаем в виде (6) а на стадии пластического течения - (7) где - предел текучести материала слоя; - модуль упругости материала слоя. Соотношения (6) соответствуют представлению типа упругих связей Прандтля [23], а (7) - модели Дагдейла [14]. В силу симметрии задачи проекции поля перемещений удовлетворяют условиям , , а вектор распределенной внешней нагрузки - . Таким образом, достаточно ограничиться рассмотрением тела 1. С учетом того, что в слое взаимодействия для данного нагружения имеет место (6) и (7), система уравнений (1)-(2) преобразуется в вариационное уравнение: (8) Решение системы (5)-(8) сводится к определению поля перемещений в теле 1 (см. рисунок) c учетом граничных условий на его торцах: (9) (10) (11) (12) Для упрощения задачи принимаем, что поле перемещений в теле 1 определено следующим образом: (13) (14) Входящий в представление (13), (14) параметр имеет геометрический смысл малого угла поворота материальной нормали к плоскости в теле 1. Согласно распределению (13)-(14) деформации в консоли будут определяться в виде (15) (16) (17) Выражения (15), (16), как теория Тимошенко [33] и работы [34-36], учитывают сдвиговые деформации и повороты нормалей в теле. Рассмотрим работу внутренних напряжений для тела 1 с учетом заданных полей деформаций (15), (16): (18) Введем в рассмотрение обобщенные силы: (19) (20) и обобщенный момент: (21) Интегрируем по частям ряд слагаемых в правой части (18) с учетом (19)-(21): (22) (23) где m,n - соответствующие пределы интегрирования по координате . Рассмотрим правую часть (8). Найдем вектор напряжений на левом торце консоли . Работа напряжений составит: (24) Подставив (18), (22)-(24) в (8) и приравняв слагаемые при равных вариациях, приходим к двум системам дифференциальных уравнений: - для участка : (25) - для участков : (26) с условиями сопряжения: (27) (28) (29) (30) (31) (32) и естественным граничным условиям на левом торце: (33) (34) (35) На правом торце из (9), (10) с учетом (13), (14) рассматриваем граничные условия: (36) (37) (38) С учетом (15)-(17) и условия плоской деформации ( ) определяющие соотношения (5) запишем в виде (39) (40) где Запишем выражения обобщенных сил (19), (20) с учетом (39), (40): (41) (42) а также момента (21): (43) В результате задача (25)-(38), (6), (7) с учетом (41)-(43) становится замкнутой относительно трех неизвестных функций: , , . Решение задачи. Рассмотрим частное решение задачи при . Запишем решение (25) при удовлетворении граничных условий (33)-(35): (44) где - постоянные интегрирования, . Общее решение (26) для участка запишем в виде (45) а для участка : (46) где - постоянные интегрирования; ; ; ; ; . Удовлетворение решений (44)-(46) двенадцати условиям сопряжения (27)-(32) и трем граничным условиям (36)-(38) дает решение поставленной задачи, из которого, при , приходим к следующему значению напряжения в слое: (47) При из решения задачи приходим к асимптотическому выражению вертикального перемещения на границе упругой и пластической зон: из последнего выражения, с учетом (3), (6) и (7), - к связи: (48) Из выражения (48) находим представление внешней нагрузки с условием конечности напряжения в виде (49) где Если , то из (49) приходим к выражению Из последнего выражения следует, что при переходе к математическому разрезу стадия упругого деформирования отсутствует и при сколь угодно малой величине внешней нагрузки возникает пластическая деформация. Такая же ситуация имеет место в модели Дагдейла [14] и в моделях с использованием когезионных сил сцепления [6-12]. Введение поправки Ирвина приводит к увеличению длины трещины за счет зоны пластичности на величину . Поскольку коэффициент интенсивности напряжений пропорционален внешней нагрузке, то и в модели квазихрупкого разрушения пластическая область возникает в начале нагружения. Таким образом, в случае математического разреза классические условия пластичности приводят к невозможности описания стадии чисто упругого деформирования. Локальное условие пластичности. Сформулируем условие пластичности, позволяющее отразить стадию упругого деформирования при вырождении физического разреза в математический. Предполагаем, что в диапазоне толщины слоя условие пластичности имеет вид (50) где - энергетическое произведение [30]. В условии (50) используется сингулярное представление напряжения в виде (47). Из (49) выразим положительное значение длины пластической области: (51) Преобразуем выражение (51) с учетом (50) при : (52) где - приведенная безразмерная длина пластической зоны. В отличие от предыдущих зависимостей внешней нагрузки от длины пластической зоны, выражение (52) учитывает стадию упругого деформирования и при переходе к математическому разрезу, когда . Здесь - упругая составляющая внешней нагрузки. Определяя в эксперименте значение упругой составляющей, находим предел упругости энергетического произведения в виде (53) В статье [32] приведены результаты экспериментов по определению трещиностойкости адгезива . Образец состоит из двух пластин толщиной B в направлении, ортогональном плоскости (см. рисунок), соединенных тонким слоем адгезива на части пластин. Образец нагружается расклинивающими силами . Определяется зависимость силы N от перемещения точки ее приложения - u. Приведены линейные размеры образца, мм: и значение модуля Юнга . Из приведенных зависимостей внешней силы от перемещения точки ее приложения можно оценить величину упругой составляющей - 30 Н. Подставляя приведенные значения в формулу (53), получим значение Н/мм. В статье [32] приведено среднее значение трещиностойкости Н/мм, полученное в результате испытания шести образцов. Выделим нагрузку, связанную с пластической деформацией , тогда и с ростом пластической зоны растет параметр . Используя экспериментальные зависимости внешней нагрузки от перемещения, исходя из выражения , можно установить закон изменения параметра от u. Заключение В результате упрощенного аналитического решения задачи упругопластического деформирования ДКБ-образца получены следующие результаты: 1. Установлена зависимость полей перемещений и напряжений от длины и толщины слоя взаимодействия. Напряжения обратно пропорциональны корню квадратному от толщины слоя и убывают по экспоненте вдоль длины слоя. 2. Использование классического условия перехода в пластическое состояние отражает упругую стадию деформирования только при конечной толщине слоя взаимодействия. При переходе к математическому разрезу пластическая составляющая деформации, как и в моделях Ирвина - Орована и Дагдейла, возникает с начала внешнего нагружения. 3. Локальное условие пластичности, определяемое критическим значением энергетического произведения, позволяет отразить стадию упругого деформирования при сколь угодно малой конечной толщине слоя взаимодействия. 4. Для экспериментальной конкретизации локального условия пластичности достаточно стандартного эксперимента над ДКБ-образцом с промежуточными разгрузками, в результате которого определяется упругая составляющая внешней нагрузки.About the authors
V. V Glagolev
Tula State University
A. A Markin
Tula State University
References
- Broberg K.B. Cracks and fracture. - Academic Press, London, UK, 1999. - 752 p.
- Черепанов Г.П. Механика хрупкого разрушения. - М.: Наука, 1974. - 640 с.
- Мураками Ю. (ред.) Справочник по коэффициентам интенсивности напряжений. Т. 1. - М.: Мир, 1990. - 448 с.
- Irvin G.R. Linear fracture mechanics, fracture transition, and fracture control // Engineering Fracture Mechanics. - 1968. - Vol. 1, no. 2. - P. 241-257. doi: 10.1016/0013-7944(68)90001-5
- Orowan E.О. In: Proc. Symposium on internal stresses in metals and alloys // London: Institute of Metals. - 1948. - P. 451.
- Баренблатт Г.И. О равновесных трещинах, образующихся при хрупком разрушении // Прикладная математика и механика. - 1959. - Т. 23, № 3. - С. 434-444.
- Kumar N., Rajagopal A., Pandey M. A rate independent cohesive zone model for modeling failure in quasi-brittle materials // Mechanics of Advanced Materials and Structures. - 2015. - Vol. 22, no. 8. - P. 681-696. doi: 10.1080/15376494.2013.855852
- Experimental characterization of cohesive zone models for thin adhesive layers loaded in mode I mode II, and mixed-mode I/II by the use of a direct method / G. Lélias, E. Paroissien, F. Lachaud, J. Morlier // International Journal of Solids and Structures. - 2019. - Vol. 158. - P. 90-115. doi: 10.1016/j.ijsolstr.2018.09.005
- Thanh L.T., Belaya L.A., Lavit I.M. A solution to the problem of elastic half-plane with a cohesive edge crack // Journal of Physics: Conference Series. - 2018. - Vol. 973, no. 1. - id. 12020. doi: 10.1088/1742-6596/973/1/012020.
- Лавит И.М. Об устойчивом росте трещины в упругопластическом материале // Проблемы прочности. - 1988. - № 7. - С. 18-23.
- The cohesive zone model: advantages, limitations and challenges / M. Elices, G.V. Guinea, J. Gómez, J. Planas // Engineering Fracture Mechanics. - 2002. - Vol. 69, no. 2. - P. 137-163. doi: 10.1016/S0013-7944(01)00083-2
- Перельмутер М.Н. Критерий роста трещин со связями в концевой области // Прикладная математика и механика. - 2007. - Т. 71, № 1. - С. 152-171.
- Леонов М. Я., Панасюк В. В. Развитие мельчайших трещин в твердом теле // Прикл. механика. - 1959. - Т. 5, № 4. - С. 391-401.
- Dugdale D.S. Yielding of steel sheets containing slits // J. Mech. Phys. Solids. - 1960. - Vol. 8, no. 2. - P. 100-104. http://dx.doi.org/10.1016/0022-5096(60)90013-2
- Захаров А.П., Шлянников В.Н., Иштыряков И.С. Пластический коэффициент интенсивности напряжений в задачах механики разрушения // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2019. - № 2. - С. 100-115. doi: 10.15593/perm.mech/2019.2.08
- Корнев В.М. Обобщенный достаточный критерий прочности. Описание зоны предразрушения // Прикладная механика и техническая физика. - 2002. - Т. 43, № 5. - С. 153-161.
- Kornev V.M., Kurguzov V.D. Multiparametric sufficient criterion of quasi-brittle fracture for complicated stress state // Engineering Fracture Mechanics. - 2008. - Vol. 75, no. 5. - P. 1099-1113. doi: 10.1016/j.engfracmech.2007.04.023
- Быковцев Г.И., Ивлев Д.Д. Теория пластичности. - Владивосток: Дальнаука, 1998. - 528 с.
- Ильюшин А.А. Пластичность. Основы общей математической теории. - М.: Изд-во АН СССР, 1963. - 272 с.
- Ишлинский А.Ю., Ивлев Д.Д. Математическая теория пластичности. - М.: ФИЗМАТЛИТ, 2001. - 701 с.
- Маркин А.А., Соколова М.Ю. Термомеханика упругопластического деформирования. - М.: ФИЗМАТЛИТ, 2013. - 320 с.
- Макклинток Ф. Пластические аспекты разрушения // Разрушение. Т. 3 - М.: Мир, 1976. - С. 67-262.
- Prandtl L. Ein Gedankenmodell für den Zerreißvorgang spröder Körper // Zeitschrift für Angewandte Mathematik und Mechanik. - 1933. - Vol. 13. - P. 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. doi: 10.1007/s10704-011-9637-3)
- Ентов В.М., Салганик Р.Л. К модели хрупкого разрушения Прандтля // Изв. АН СССР. МТТ. - 1968. - № 6. - С. 87-99.
- Салганик Р.Л., Мищенко А.А., Федотов А.А. Модель трещины Прандтля и ее применение для решения задачи механики контактного взаимодействия // К 75-летию со дня рождения профессора Владимира Марковича Ентова. - Ижевск: Ин-т компьютерных исследований, 2012. - 180 с.
- Салганик Р.Л., Мищенко А.А., Федотов А.А. Напряженное состояние в окрестности выработки, пройденной в глубокозалегающем горизонтальном пласте // Физико-технические проблемы разработки полезных ископаемых. - 2015. - № 2. - С. 24-33.
- Shear loading of a cantilever with a crack-like defect explicitly including linear parameters / F. Berto, V.V. Glagolev, L.V. Glagolev, A.A. Markin // International Journal of Solids and Structures. - 2020. - Vol. 193-194. - P. 447-454. doi: 10.1016/j.ijsolstr.2020.02.039
- 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
- Глаголев В.В., Глаголев Л.В., Маркин А.А. Определение напряженно-деформированного состояния упругопластических тел с боковым трещиноподобным дефектом с использованием модели с линейным размером // Прикладная механика и техническая физика. - 2018. - Т. 59, № 6. - С. 143-154. doi: 10.15372/PMTF20180610
- 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
- Berto F., Glagolev V.V., Markin A.A. Relationship between Jc and the dissipation energy in the adhesive layer of a layered composite // International Journal of Fracture. - 2020. - Vol. 224, no. 2. - P. 277-284. doi: 10.1007/s10704-020-00464-0
- Santos M.A.S., Campilho R.D.S.G. Mixed-mode fracture analysis of composite bonded joints considering adhesives of different ductility // International Journal of Fracture. - 2017. - Vol. 207, no. 1. - P. 55-71. doi: 10.1007/s10704-017-0219-x
- Тимошенко С.П., Войновский-Кригер С. Пластины и оболочки. - М.: Физматгиз, 1963. - 636 с.
- Mattei O., Bardella L. A structural model for plane sandwich beams including transverse core deformability and arbitrary boundary conditions // Eur. J. Mech. A-Solid. - 2016. - Vol. 58. - P. 172-186. doi: 10.1016/j.euromechsol.2016.01.015
- Panettieri E., Fanteria D., Danzi F. Delaminations growth in compression after impact test simulations: Influence of cohesive elements parameters on numerical results // Composite Structures. - 2016. - Vol. 137. - P. 140-147. doi: 10.1016/J.COMPSTRUCT.2015.11.018
- Panteghini A., Bardella L. Structural theory and finite element modelling of linear elastic sandwich beams subject to severe boundary conditions // Eur. J. Mech. A-Solid. - 2017. - Vol. 61. - P. 393-407. doi: 10.1016/j.euromechsol.2016.10.012