Энергетическое произведение в модели трещиноподобного дефекта при нагружении типа моды II

Аннотация


Рассматривается нагружение трещиноподобного дефекта по моде II. В отличие от классического представления трещины в виде математического разреза предлагаемая модель определяет трещину в виде физического разреза с характерным линейным размером. Мысленное продолжение физического разреза в твердом теле формирует слой взаимодействия. Существенно, что напряженно-деформированное состояние слоя не вносит в модель трещины сингулярность. Произведение приращения удельной свободной энергии в торцевом квадратном элементе слоя на линейный размер определяет его энергетическое произведение (ЭП). Объектом исследования является двухконсольный образец, а предметом исследования - энергетическое произведение в торцевом элементе слоя взаимодействия. Внешняя нагрузка консолей приводит к их горизонтальным антисимметричным перемещениям, которые формируют в слое взаимодействия однородные сдвиговые деформации. Из условий равновесия консолей в вариационной форме с учетом гипотезы однородности осевых деформаций и c учетом их обжатия получена система дифференциальных уравнений, связывающая напряженное состояние в слое и консолях. Исследовано решение характеристического уравнения системы при различных отношениях толщин слоя и консолей. Показано, что при отношениях меньше определенного значения, зависящего от коэффициента Пуассона, имеют место вещественные корни. В рамках вещественных корней характеристического уравнения получено аналитическое решение поставленной задачи. При условии пренебрежения обжатием консолей найдено упрощенное решение. Определены деформации в слое с учетом обжатия консолей и без него. Проведен анализ зависимости ЭП от отношений толщин слоя и консолей. Показано, что при отношении толщин 10-6 и менее ЭП практически не меняет своего значения. Учет обжатия консолей дает различие в значениях ЭП порядка 20 % по отношению к упрощенному решению задачи.

Полный текст

Введение Классическая постановка задачи о трещине моды 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 % уменьшает пороговое значение ЭП. Если в результате эксперимента определить критическое значение внешней нагрузки при заданной толщине консоли, то можно найти критическое пороговое значение ЭП. Универсальность данного значения можно установить, сравнивая величины критических нагрузок, определенных теоретически из условия постоянства ЭП при различных толщинах консолей с результатами экспериментов.

Об авторах

В В Глаголев

Тульский государственный университет

Л В Глаголев

Тульский государственный университет

А А Маркин

Тульский государственный университет

Список литературы

  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

Статистика

Просмотры

Аннотация - 362

PDF (Russian) - 214

Cited-By


PlumX


© Глаголев В.В., Глаголев Л.В., Маркин А.А., 2019

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution-NonCommercial 4.0 International License.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах