Научная статья на тему 'Ударное нагружение полосы с центральной трещиной'

Ударное нагружение полосы с центральной трещиной Текст научной статьи по специальности «Физика»

CC BY
127
16
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ТРЕЩИНА / КОЭФФИЦИЕНТ ИНТЕНСИВНОСТИ НАПРЯЖЕНИЙ / ДИНАМИЧЕСКАЯ МЕХАНИКА РАЗРУШЕНИЯ / КОГЕЗИОННЫЕ КОНЕЧНЫЕ ЭЛЕМЕНТЫ / СИЛЫ СЦЕПЛЕНИЯ / МЕТОД ПРЯМЫХ / СХЕМА КРАНКА-НИКОЛСОН / ЗАДАЧА ЧЕНА / МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ / CRACK / STRESS INTENSITY FACTOR / DYNAMIC FRACTURE MECHANICS / COHESIVE FINITE ELEMENTS / COHESIVE FORCES / METHOD OF LINES / CRANK-NICHOLSON SCHEME / CHEN'S PROBLEM / FINITE ELEMENT METHOD

Аннотация научной статьи по физике, автор научной работы — Малик А. В., Рязанцева И. Э., Лавит И. М.

Рассматривается задача расчета зависимости коэффициента интенсивности напряжений от времени для полосы, находящейся в состоянии плоской деформации и ослабленной неподвижной центральной трещиной нормального разрыва. К основаниям полосы мгновенно прикладывается равномерно распределенная нагрузка, остающаяся далее неизменной. Используется модель трещины с силами сцепления, распределение которых подчиняется постулатам Баренблатта. При этом коэффициент интенсивности напряжений находится в результате вычисления интенсивности высвобожденной энергии, определяемой через силы сцепления. Решение задачи находится новым численным методом, представляющим собой адаптацию метода прямых к решению задач динамической механики разрушения. Для интегрирования по времени используется неявная конечно-разностная схема Кранка-Николсон. Краевые задачи, возникающие на каждом шаге интегрирования по времени, решаются методом конечных элементов. Для того чтобы решение задачи удовлетворяло постулатам Баренблатта, используются специальные когезионные конечные элементы, ранее уже применявшиеся для решения квазистатических задач нелинейной механики разрушения. За счет введения дополнительных степеней свободы в узлах, лежащих на линии трещины, удается обеспечить плавное смыкание кромок трещины в ее кончике, что эквивалентно отсутствию сингулярности полей напряжений и деформаций в ее кончике. При этом силы сцепления вычисляются как реакции связей. Область их действия (зона сцепления) локализована в пределах конечного элемента, прилегающего к кончику трещины. Таким образом, чем мельче сетка конечных элементов, тем точнее удовлетворяется требование теории Баренблатта о малости длины зоны сцепления по сравнению с длиной трещины. Поставленная задача, называемая задачей Чена, решалась ранее разными исследователями, применявшими различные методы. Близость полученных при этом результатов дает основание считать задачу Чена тестовой. Ее решение, полученное разработанным методом, удовлетворительно согласуется с данными других исследователей.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Малик А. В., Рязанцева И. Э., Лавит И. М.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

A shock loading on a bar with a central crack

Linear thermoelasticity is studied of a plane regular truss formed by four families of The paper is concerned with the problem of calculating the time dependence of the stress intensity factor for a plane-strain bar with a stationary central crack caused by the opening mode. A uniformly distributed load has been immediately imposed on the basis of the bar and remained unchanged later. The model of the crack with cohesive forces distributed by Barenblatt’s postulates is used. In this case the stress intensity factor is a result of the calculation of the released energy that is determined by the cohesive forces. The solution is found with a new numerical method that is an adaptation of the method of lines to dynamic fracture mechanics problems. The Crank-Nicolson implicit finite-difference scheme is used for time integration. Boundary problems arising at each step of time integration are solved by the finite element method. The special cohesive finite elements are used, so that the solution of the problem could satisfy Barenblatt’s postulates. Previously these elements were used to solve quasi-static nonlinear fracture mechanics problems. By introducing the additional degrees of freedom of the nodes lying on the crack line, it becomes possible to ensure a smooth closing of the crack edges at its tip; and it is equivalent to the absence of the singularity stress and strain fields at its tip. The cohesive forces are calculated as the constraints. Their field of action (cohesive zone) is localized within the finite element which is adjacent to the tip of the crack. Thus, the smaller the finite element mesh is, the better it satisfies the requirement of Barenblatt’s theory. This requirement concerns the length of the cohesive zone which is small compared to the length of the crack. The stated problem is called Chen’s problem and had earlier been solved by researchers that used different methods. The proximity of the obtained results makes it possible to consider Chen’s problem as a test; and its solution obtained by the developed method agrees well with the data of other researchers.

Текст научной работы на тему «Ударное нагружение полосы с центральной трещиной»

Малик А.В., Рязанцева И.Э., Лавит И.М. Ударное нагружение полосы с центральной трещиной // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 2. - С. 125-135. DOI: 10.15593/perm.mech/2017.2.07

Malik A.V., Ryazantseva I.E., Lavit I.M. A shock loading on a bar with a central crack. PNRPUMechanics Bulletin, 2017, no.2, pp. 125-135. DOI: 10.15593/perm.mech/2017.2.07

ВЕСТНИК ПНИПУ. МЕХАНИКА

№ 2,2017 PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/ini7

DOI 10.15593/perm.mech/2017.2.07 УДК 539.375

УДАРНОЕ НАГРУЖЕНИЕ ПОЛОСЫ С ЦЕНТРАЛЬНОЙ ТРЕЩИНОЙ А.В. Малик, И.Э. Рязанцева, И.М. Лавит

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

О СТАТЬЕ

АННОТАЦИЯ

Получена: 17 ноября 2016 г. Принята: 06 мая 2017 г. Опубликована: 30 июня 2017 г.

Ключевые слова:

трещина, коэффициент интенсивности напряжений, динамическая механика разрушения, когезионные конечные элементы, силы сцепления, метод прямых, схема Кранка-Николсон, задача Чена, метод конечных элементов.

Рассматривается задача расчета зависимости коэффициента интенсивности напряжений от времени для полосы, находящейся в состоянии плоской деформации и ослабленной неподвижной центральной трещиной нормального разрыва. К основаниям полосы мгновенно прикладывается равномерно распределенная нагрузка, остающаяся далее неизменной. Используется модель трещины с силами сцепления, распределение которых подчиняется постулатам Баренблатта. При этом коэффициент интенсивности напряжений находится в результате вычисления интенсивности высвобожденной энергии, определяемой через силы сцепления. Решение задачи находится новым численным методом, представляющим собой адаптацию метода прямых к решению задач динамической механики разрушения. Для интегрирования по времени используется неявная конечно-разностная схема Кранка-Николсон. Краевые задачи, возникающие на каждом шаге интегрирования по времени, решаются методом конечных элементов. Для того чтобы решение задачи удовлетворяло постулатам Баренблатта, используются специальные коге-зионные конечные элементы, ранее уже применявшиеся для решения квазистатических задач нелинейной механики разрушения. За счет введения дополнительных степеней свободы в узлах, лежащих на линии трещины, удается обеспечить плавное смыкание кромок трещины в ее кончике, что эквивалентно отсутствию сингулярности полей напряжений и деформаций в ее кончике. При этом силы сцепления вычисляются как реакции связей. Область их действия (зона сцепления) локализована в пределах конечного элемента, прилегающего к кончику трещины. Таким образом, чем мельче сетка конечных элементов, тем точнее удовлетворяется требование теории Баренблатта о малости длины зоны сцепления по сравнению с длиной трещины.

Поставленная задача, называемая задачей Чена, решалась ранее разными исследователями, применявшими различные методы. Близость полученных при этом результатов дает основание считать задачу Чена тестовой. Ее решение, полученное разработанным методом, удовлетворительно согласуется с данными других исследователей.

© ПНИПУ

© Малик Александр Васильевич - аспирант, e-mail: runer10@mail.ru Рязанцева Инна Эдуардовна - магистрант, e-mail: Inna.ori.ru@mail.ru

Лавит Игорь Михайлович - доктор физико-математических наук, профессор, e-mail: IgorLavit@yandex.ru

Alexander V. Malik - PhD Student, e-mail: runer10@mail.ru

Inna E. Ryazantseva - Postgraduate Student, e-mail: Inna.ori.ru@mail.ru

Igor M. Lavit - Doctor of Physical and Mathematical Sciences, Professor, e-mail: IgorLavit@yandex.ru

A SHOCK LOADING ON A BAR WITH A CENTRAL CRACK

A.V. Malik, I.E. Ryazantseva, I.M. Lavit

Tula State University, Tula, Russian Federation

ABSTRACT

Linear thermoelasticity is studied of a plane regular truss formed by four families of The paper is concerned with the problem of calculating the time dependence of the stress intensity factor for a plane-strain bar with a stationary central crack caused by the opening mode. A uniformly distributed load has been immediately imposed on the basis of the bar and remained unchanged later. The model of the crack with cohesive forces distributed by Barenblatt's postulates is used. In this case the stress intensity factor is a result of the calculation of the released energy that is determined by the cohesive forces. The solution is found with a new numerical method that is an adaptation of the method of lines to dynamic fracture mechanics problems. The Crank-Nicolson implicit finite-difference scheme is used for time integration. Boundary problems arising at each step of time integration are solved by the finite element method. The special cohesive finite elements are used, so that the solution of the problem could satisfy Barenblatt's postulates. Previously these elements were used to solve quasi-static nonlinear fracture mechanics problems. By introducing the additional degrees of freedom of the nodes lying on the crack line, it becomes possible to ensure a smooth closing of the crack edges at its tip; and it is equivalent to the absence of the singularity stress and strain fields at its tip. The cohesive forces are calculated as the constraints. Their field of action (cohesive zone) is localized within the finite element which is adjacent to the tip of the crack. Thus, the smaller the finite element mesh is, the better it satisfies the requirement of Barenblatt's theory. This requirement concerns the length of the cohesive zone which is small compared to the length of the crack.

The stated problem is called Chen's problem and had earlier been solved by researchers that used different methods. The proximity of the obtained results makes it possible to consider Chen's problem as a test; and its solution obtained by the developed method agrees well with the data of other researchers.

© PNRPU

Введение

Формулировка условий старта трещины при динамическом нагружении - сложная проблема, до сих пор окончательно не решенная [1-3]. Хотя предлагаемые теории различны (см., например, работы [4-6]), все они включают в себя как составную часть зависимость коэффициента интенсивности напряжений (КИН) от времени. Ниже везде речь идет только о трещинах нормального разрыва, и поэтому под КИН понимается величина KI. Определение функции ), где ? - время, представляет собой самостоятельную задачу, аналитическое решение которой возможно только для бесконечных областей [1-3]. Существующие численные методы (см. обзоры [7-10]) обладают теми или иными недостатками, поэтому разработка новых эффективных методов расчета представляет интерес.

По-видимому, наиболее перспективным направлением вычислительной динамики разрушения является создание методов, основанных на модели когезионной трещины, кромки которой притягиваются друг к другу силами сцепления [11]. Главное достоинство этих методов - алгоритмическая простота, позволяющая решать ими сложные задачи. Первый такой метод был предложен в работе [12] и применен к исследованию ветвления трещин. Дальнейшее развитие - это, как правило, объединение подхода работы [12] с тем или иным вариантом метода конечных элементов, а также использование различных конституционных соотношений для сил сцепления (см., например, работы [11, 13-16]).

ARTICLE INFO

Received: 17 November 2016 Accepted: 06 May 2017 Published: 30 June 2017

Keywords:

crack; stress intensity factor, dynamic fracture mechanics, cohesive finite elements, cohesive forces, method of lines, Crank-Nicholson scheme, Chen's problem, finite element method.

Для того чтобы вычислять КИН, используя модель когезионной трещины, необходимо, чтобы распределение сил сцепления подчинялось постулатам Баренблатта [17, 18, 2]. Большинство методов, в том числе и упомянутые выше, этому требованию не удовлетворяет. Одним из исключений является метод, изложенный в работе [19] и предназначенный для решения квазистатических задач. Его обобщение на динамические задачи неизвестно и, по-видимому, невозможно, так как в нем существенным этапом является вычисление J -интеграла, а последний перестает быть инвариантным при необходимости учитывать инерционные силы. Есть и другой метод, который также был применен ранее для решения квазистатических задач [20]. Этот метод допускает обобщение на динамические задачи. Оно представлено ниже. В нем, по сравнению с работой [20], учтены силы инерции и зависимость всех полевых величин от времени. Некоторые аспекты получившегося метода решения динамических задач приведены в статье [21].

Он представляет собой модификацию метода прямых (метода Роте) [22-24]. С использованием конечно-разностной дискретизации по времени решение задачи сводится к совокупности последовательно решаемых краевых задач теории упругости. Выбор конечно-разностной схемы интегрирования по времени неоднозначен: возможно применение любой надежной схемы, как явной, так и неявной. В настоящей работе использована широко применяемая, в том числе и совместно с методом конечных элементов, неявная схема Кранка-Николсон [25].

Для решения краевых задач используется метод конечных элементов [25]. Проблема согласования распределения сил сцепления с постулатами Баренблатта решается с помощью включения в конечно-элементную сетку специальных когезионных элементов [20].

Оценка эффективности численного метода основывается на решении тестовых задач. Обычно эти задачи имеют аналитическое решение. Но в механике разрушения таких задач для конечных областей нет, и поэтому тестовые задачи приходится выбирать из задач, решаемых численно. И если подходящих статических задач можно найти немало [26], то для случая внезапно прикладываемой нагрузки - только одну. Это задача о полосе с неподвижной центральной трещиной; к основаниям полосы мгновенно прикладывается равномерно распределенная нагрузка, остающаяся далее неизменной. Впервые эту задачу решил Чен [27]; она была предметом ряда последующих исследований ([1, 28-33] и др.). В работе [30] эта задача названа задачей Чена, и это название закрепилось в более поздних публикациях.

Считать задачу Чена тестовой позволяет то, что авторы упомянутых исследований получили близкие результаты, причем они применяли различные методы. Так, в работах [27, 30] использовался метод конечных разностей для интегрирования и по времени, и по пространственным координатам с применением явной трехслойной схемы. В работах [1, 28] для интегрирования по пространству использовался метод конечных элементов, причем в работе [1] для моделирования особенности поля напряжений в кончике трещины применялись специальные сингулярные конечные элементы. В результате получалась система обыкновенных дифференциальных уравнений, которая интегрировалась численно по явной трехслойной схеме Ньюмарка. В работах [29, 31-33] применялся метод граничных элементов. В работе [29] для моделирования особенности поля напряжений использовались сингулярные граничные элементы; были найдены зависящие от времени функции Грина, благодаря чему интегрирование по времени свелось к вычислению определенных интегралов. В работах [31, 32] применялось преобразование Лапласа по времени, а в работе [33] - преобразование Фурье.

1. Постановка задачи

Рассматриваются малые деформации изотропного однородного линейно-упругого тела. Тело содержит трещину, длина которой неизменна; форма тела и распределение нагрузок таковы, что эта трещина является трещиной нормального отрыва, то есть

К, * 0; К11 = Кш = 0,

где К1, К11, Кш - коэффициенты интенсивности напряжений [1].

Пусть рассматриваемое тело - полоса, находящаяся в состоянии плоской деформации. Поперечное сечение полосы изображено на рис. 1.

В начальный момент времени к двум противоположным сторонам прямоугольника, параллельным трещине, мгновенно прикладывается равномерно распределенная нагрузка величиной q, далее остающаяся неизменной. Ставится задача определить зависимость К1 (7).

2H

П A I a i I а n

H

2a

W - a

2W

Задача имеет две оси симметрии, проходящие через середины противоположных сторон прямоугольника. Поэтому расчетная схема - это четверть сечения, вырезанная осями симметрии (рис. 2).

q

Рис. 1. Поперечное сечение полосы -прямоугольник с симметрично расположенной

прямолинейной центральной трещиной Fig. 1. The cross-section of the strip. This is the rectangle with a symmetrical rectilinear central crack

Рис. 2. Расчетная схема задачи Чена Fig. 2. The design model of Chen's problem

Граничные условия задачи определяются условиями нагружения и условиями симметрии. Пусть ui - вектор перемещений в декартовой системе координат х (см. рис. 2);

- внешняя распределенная нагрузка; г = 1,2. Граничные условия записываются следующим образом. На участке контура АВ : и1 = 0, р2 = 0 - условия симметрии; на участке ВС : р1 = 0, р2 = q; на участке СБ : р1 = 0 , р2 = 0; на участке БЕ : и2 = 0, р1 = 0 -условия симметрии; на участке ЕА (кромка трещины): р1 = 0, р2 = 0 .

Деформирование полосы определяется следующими соотношениями [1, 2]:

q

=1 (дкит + дтик ); акт = 3КгЬкт + Ю ( - гЬкт );

8= 3 Вкт5кт ; дт°кт = Рд t\

дк = д/дхк; д( = д/дг; к, т = 1,2,

где вкт - тензор деформаций; 5кт - символ Кронекера; в - средняя деформация, К, G -объемный модуль и модуль сдвига; скт - тензор напряжений; р - плотность материала. В начальный момент времени (I = 0) полоса покоится.

2. Метод решения задачи

Численное решение задачи получается на основе принципа возможных перемещений

\(рд^к 5ик +°кт 5в кт ) =\рк Ьщё1;

(2)

где - площадь области (см. рис. 1); I - ее граничный контур; 5 - символ вариации; Ук - поле скоростей.

В соответствии с методом прямых [22-24] перейдем в уравнении (2) к конечно-разностному представлению, используя неявную схему Кранка-Николсон [25]. Пусть Ы -величина шага интегрирования по времени, п =1, 2... - номер шага. Конечно-разностное представление производной по времени на п-м шаге имеет вид

д

п п-1

У - У

ы

где уп, уп - значения у на границах п -го временного интеграла (шага). Величины, не содержащие производных по времени, представляются на п-м шаге интегрирования по времени в виде

У - 2(уп+Уп-1).

С учетом этих формул конечно-разностное представление уравнения (2) записывается в виде

Р( -VI 1)

М

5ик + 2 ((т т )5вкт

& = 2 к+рп-1)),

(3)

п п-1

к -ик =-1уп + Уп

М

2 (+*1).

Величины с индексами п -1 известны из решения для предыдущего шага. Из системы (3) определяются величины с индексами п .

Выразим из второго уравнения системы (3) величину уЩ и подставим в первое уравнение. Получим

5

i.n 2 /„ .n _,n-1 \ -.,n-1 Vk (Uk - uk j-vk .

4p

(At j

2 К§Uk + aL5skm

dS = {( + pkn-1 jdl + {

4p

n-1

At

^ + vn 1

At

§Uk -a"kmdSkm

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

dS. (5)

Таким образом, система (3) распадается на два последовательно решаемых уравнения. Вначале из уравнения (5) определяются перемещения ипк, затем из уравнения (4) -

скорости упк.

Уравнение (5) представляет собой вариационное уравнение по пространственным переменным, которое необходимо решать на каждом шаге интегрирования по времени. Его решение находится методом конечных элементов [25].

Специфика, вносимая трещиной в краевую задачу теорию упругости, учитывается с помощью специальных когезионных конечных элементов - так, как это описано в работе [20].

Когезионные элементы составляют горизонтальный ряд, прилегающий к кромке трещины (рис. 3).

В локальных координатах [25] и обычные, и когезионные элементы представляют собой одинаковые квадраты. Глобальные координаты точек элемента определяются как

хт = Ц (г^Ц ^2)Х1; /,/ = 1,2; ^2 е[-1;1];

Ц (*) = 2 (1 - г); Ц (*) = 2 (1+

где Хьт - заданные глобальные декартовы координаты узлов (верхние индексы обозначают но-

Рис. 3. Конечно-элементная сетка (принципиальная схема).

Когезионные элементы заштрихованы мер узла в локальной нумерации); г1з г2 - ло-

Fig. 3. Finite-element mesh (schematic). Cohesive elements are shaded

кальные координаты точек элемента; Ц (г) -

интерполяционные полиномы Лагранжа. Перемещения точек обычного конечного элемента задаются аналогичной формулой [25]

ит = Ц (г )Ц (г2)ит, (6)

где ит - узловые перемещения.

Перемещения точек когезионного конечного элемента определяются формулами [20]

иг = Ц (г)Ц (г2)и*;

1 (7)

и2 = Н (гх )Ц (г 2 )и 21 + Ц (г )Ц (г2 + - Ах1И1+2 (гх)Ц (г2 ^

где Ах1 - размер элемента по оси абсцисс; и'3 - значения производной ди2/дх1 в //-м узле (эти величины определяются только для узлов, лежащих на оси абсцисс - прямой, на которой расположена трещина); Нт (г) - интерполяционные полиномы Эрмита,

H (*) = , ( - 3z + z3 )• H 2 (z ) = 4 (2 + 3z - z3), H (z) = 4(l -z-z2 + z3), H4 (z) = i(-1

I ,2.3

-1 - z + z + z

Представление перемещений в виде (6), (7) обеспечивает межэлементную непрерывность поля перемещений во всей области 5. Необходимость введения дополнительных узловых степеней свободы обусловлена требованием плавного смыкания кромок

трещины в ее кончике [17].

Конечно-элементная сетка строится так, чтобы кончик трещины совпадал с каким-либо узлом. В этом узле полагается не только и2 = 0, но и ди2/ дх1 = 0, что и обеспечивает упомянутую плавность смыкания кромок трещины.

Эти условия представляют собой с механической точки зрения связи, уменьшающие число степеней свободы узла, совпадающего с кончиком трещины. Реакции этих связей, QB и Мв (рис. 4), определяются обычным об-

разом после решения системы линейных алгебраических уравнений относительно узловых неизвестных.

При продвижении трещины на длину ко-

Рис. 4. Реакции связей и узловые перемещения в узлах конечного элемента, прилегающего к кончику трещины (узел В)

at the nodes of the finite element which is adjacent to the crack tip (node B)

нечного элемента AXi эти реакции связей pig. 4. Reactions of ties and nodal displacements уменьшаются (по модулю) до нуля, а соответствующие им перемещения u2 и угол поворота ди2/дхх приобретают конечные значения.

Обычно принимается допущение о неизменяемости конфигурации концевой области при продвижении трещины [17]. Это значит, что перемещения, которые получит узел B при продвижении трещины на величину Ах1, будут равны перемещениям узла A при совпадении кончика трещины с узлом B. При этом величина высвобожденной энергии, отнесенная к приращению длины трещины Ах1, определяется как

J = --

Ax,

B"- 2 A

+ MB

ди

2 А

дх,

(8)

1 J

Здесь учтено, что трещина имеет две кромки.

Можно показать [17, 18], что на расстояниях, значительно превышающих длину ко-гезионной зоны, поле напряжений определяется асимптотическими формулами [1-3], получаемыми при решении задач при отсутствии сил сцепления. При этом связь между КИН и величиной J дается формулой [1, 2]

KI =±,

I EJ

TW

(9)

где Е - модуль Юнга; V - коэффициент Пуассона.

Если из формулы (8) получается, что J < 0, то это значит, что корневая асимптотика [1-3] еще не установилась и, следовательно, КИН равен нулю. Знак минус в формуле (9)

соответствует случаю, когда и2А < 0 . Не следует считать, что при этом трещина закрывается и теория становится неприменимой. Разрез, не имеющий ширины, как модель трещины - это удобная абстракция. В реальных трещинах расстояние между кромками всегда конечно. Отрицательное значение и2А приводит к уменьшению этого расстояния, но не обязательно к смыканию кромок.

3. Решение задачи Чена

Расчеты проводились при следующих исходных данных [27]: модуль Юнга Е = 2 -105 МПа, коэффициент Пуассона у = 0,3, плотность р = 5000 кг/м3, нагрузка q = 400 МПа, размеры полосы: Ж = 10 мм, Н = 20 мм, длина трещины а = 2,4 мм. Рассматриваемый интервал времени tmax = 14 мкс.

Объемный модуль и модуль сдвига выражаются через Е и V известными формулами

Е „ Е

K =

G =■

3 (1 - 2v)' 2 (1 + v)'

Сходимость численного решения исследовалась на различных конечно-элементных сетках при различном количестве шагов по времени N. В таблице приведены результаты вычисления максимальной величины безразмерного КИН

* K7 K*= 1

урШ

для трех вариантов расчета, отличающихся значениями N и числом конечных элементов. Во всех вариантах отношение числа конечных элементов по оси ординат п2 к числу конечных элементов по оси абсцисс п принималось постоянным, п2 /п1 = 2.

Максимальные значения K* для различных вариантов расчета The highest values of K* for different calculation variants

N n K*

450 50 2,639

900 50 2,652

900 100 2,648

На рис. 5 представлены результаты расчета при N = 900, п1 = 100 и п2 = 200. При

этом относительная погрешность расчета, согласно данным таблицы, не превышает 1 %. На рис. 5 по оси абсцисс отложено безразмерное время [30]

т =

c1t 2W'

где c1 - скорость волны расширения [1-3],

Cl 1

3K + 4G

к7 2,5

2

1,5

1

0,5 0

-0,5

Рис. 5. Зависимость КИН от времени: 1 - решение Чена [27];

2 - результаты работы [30]; 3 - результаты расчета разработанным методом Fig. 5. The dependence of stress intensity factor on time: 1 is Chen's solution [27];

2 are the results of [30]; 3 are the calculation results based on the developed method

Как следует из рис. 5, результаты расчетов изложенным методом практически совпадают с результатами работ [27, 30]. Резкие качественные изменения различных участков графиков обусловлены распространением, взаимодействием и отражением от границ области волн расширения, сдвига и волн Рэлея. Обсуждение этих эффектов можно найти в работах [1, 27, 30].

Заключение

Решение задачи Чена изложенным методом свидетельствует о его приемлемой точности. Другие его достоинства: простота и возможность применения к более сложным динамическим задачам, в частности задачам о распространении трещин, в том числе и с учетом пластического деформирования. Примеры решения таких задач в квазистатической постановке даны в работе [20].

Библиографический список

1. Партон В.З., Борисковский В.Г. Динамика хрупкого разрушения. - М.: Машиностроение, 1988. - 240 с.

2. Freund L.B. Dynamic fracture mechanics. - New York: Cambridge university press, 1998. - 563 p.

3. Ravi-Chandar K. Dynamic fracture. - Amsterdam: Elsevier, 2004. - 254 p.

4. Морозов Н.Ф., Петров Ю.В. Проблемы динамики разрушения твердых тел. - СПб.: Изд-во С.-Петерб. у-та, 1997. - 129 с.

5. Zhao Y.-P. Suggestion of a new criterion of dynamic fracture // Int. J. Fract. - 1995. - Vol. 71. -P. R77-R78.

6. Liu C., Knauss W.G., Rosakis A.J. On the modeling of fracture of brittle solids // Int. J. Fract. -1998. - Vol. 90. - P. 103-118.

7. Song J.-H., Wang H., Belytschko T. A comparative study on finite element methods for dynamic fracture // Comput. Mech. - 2008. - Vol. 42. - P. 239-250.

8. Братов В.А. Численные модели динамики разрушения // Вычислительная механика сплошных сред. - 2009. - Т. 2, № 3. - С. 5-16.

9. Rabczuk T. Computational methods for fracture in brittle and quasi-brittle solids: state-of-the-art review and future perspectives // ISRN Applied Mathematics. - 2013. - Vol. 2013. - Art. ID849231. -38 p. - URL: http://dx.doi.org/10.1155/2013/849231.

10. Fineberg J., Bouchbinder E. Recent developments in dynamic fracture: some perspectives // Int. J. Fract. - 2015. - Vol. 196. - P. 33-57.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

11. The cohesive zone model: advantages, limitations and challenges / M. Elices, G.V. Guinea, J. Gomez, J. Planas // Eng. Fract. Mech. - 2002. - Vol. 69. - P. 137-163.

12. Xu X.-P., Needleman A. Numerical simulations of fast crack growth in brittle solids // J. Mech. Phys. Solids. - 1994. - Vol. 42. - P. 1397-1434.

13. De Borst R., Remmers J.J.C., Needleman A. Mesh-independent discrete numerical representations of cohesive-zone models // Eng. Fract. Mech. - 2006. - Vol. 73. - P. 160-177.

14. Bardenhagen S.G., Nairn J.A., Lu H. Simulation of dynamic fracture with the material point method using a mixed J-integral and cohesive law approach // Int. J. Fract. - 2011. - Vol. 170. - P. 49-66.

15. Agwai A., Guven I., Madenci E. Predicting crack propagation with peridynamics: a comparative study // Int. J. Fract. - 2011. - Vol. 171. - P. 65-78.

16. Javidrad F., Mashayekhy M. A cohesive zone model for crack growth simulation in AISI 304 steel // J. Solid Mechanics. - 2014. - Vol. 6. - P. 378-388.

17. Баренблатт Г.И. Математическая теория равновесных трещин, образующихся при хрупком разрушении // Журн. прикл. механики и техн. физики. - 1961. - № 4. - С. 3-56.

18. Willis J.R. A comparison of the fracture criteria of Griffith and Barenblatt // J. Mech. Phys. Solids. - 1967. - Vol. 15. - P. 151-162.

19. Moes N., Belytschko T. Extended finite element method for cohesive crack growth // Eng. Fract. Mech. - 2002. - Vol. 69. - P. 813-833.

20. Лавит И.М. Об устойчивом росте трещины в упругопластическом материале // Проблемы прочности. - 1988. - № 7. - С. 18-23.

21. Малик А.В., Белая Л.А., Лавит И.М. О динамическом нагружении тела с трещиной в условиях плоской деформации // Фундаментальные и прикладные проблемы техники и технологии. -2016. - № 1 (315). - С. 3-10.

22. Михлин С.Г. Прямые методы в математической физике. - М.: ГИТТЛ, 1950. - 428 с.

23. Лисковец О.А. Метод прямых // Дифференциальные уравнения. - 1965. - Т. 1, № 12. -С.1662-1678.

24. Ладыженская О.А. Краевые задачи математической физики. - М.: Наука, 1973. - 408 с.

25. Зенкевич О. Метод конечных элементов в технике. - М.: Мир, 1975. - 541 с.

26. Справочник по коэффициентам интенсивности напряжений. Т. 1 / под ред. Ю. Мурака-ми. - М.: Мир, 1990. - 448 с.

27. Chen Y.M. Numerical computation of dynamic stress intensity factors by a lagrangian finite-difference method (the HEMP code) // Eng. Fract. Mech. - 1975. - Vol. 7. - P. 653-660.

28. Brickstad B. A FEM analysis of crack arrest experiments // Int. J. Fract. - 1983. - Vol. 21. - P. 177-194.

29. Israil A.S.M., Dargush G.F. Dynamic fracture mechanics studies by time-domain BEM // Eng. Fract. Mech. - 1991. - Vol. 39. - P. 315-328.

30. Lin X., Ballmann J. Re-consideration of Chen's problem by finite difference method // Eng. Fract. Mech. - 1993. - Vol. 44. - P. 735-739.

31. Wen P.H., Aliabadi M.H., Rooke D.P. Application of the weight function method to two-dimensional elastodynamics fracture mechanics // Int. J. Fract. - 1996. - Vol. 76. - P. 193-206.

32. Wen P.H., Aliabadi M.H., Rooke D.P. A contour integral method for dynamic stress intensity factors // Theor. and Appl. Fract. Mechanics. - 1997. - Vol. 27. - P. 29-41.

33. Phan A.-V. A non-singular boundary integral formula for frequency domain analysis of the dynamic T-stress // Int. J. Fract. - 2012. - Vol. 173. - P. 37-48.

References

1. Parton V.Z., Boriskovsky V.G. Dynamic Fracture Mechanics. vol. 1: Stationary Cracks. New York, Hemisphere Publishing Corporation, 1989, 322 p.

2. Freund L.B. Dynamic fracture mechanics. New York, Cambridge university press, 1998, 563 p.

3. Ravi-Chandar K. Dynamic fracture. Amsterdam, Elsevier, 2004, 254 p.

4. Morozov N.F., Petrov Y.V. Problemy dinamiki razrusheniia tverdykh tel [ Problems of dynamics fracture of solids]. Sankt-Peterburg, Izdatel'stvo Sankt-Peterburgskogo universiteta, 1997, 129 p.

5. Zhao Y.-P. Suggestion of a new criterion of dynamic fracture. Int. J. Fract., 1995, vol. 71, pp. R77-R78.

6. Liu C., Knauss W.G., Rosakis A.J. On the modeling of fracture of brittle solids. Int. J. Fract., 1998, vol. 90, pp. 103-118.

7. Song J.-H., Wang H., Belytschko T. A comparative study on finite element methods for dynamic fracture. Comput. Mech, 2008, vol. 42, pp. 239-250.

8. Bratov V. A., Numerical models of the dynamics of destruction. Computational continuum mechanics, 2009, vol. 2, no.3, pp. 5-16.

9. Rabczuk T. Computational methods for fracture in brittle and quasi-brittle solids: state-of-the-art review and future perspectives. ISRN Applied Mathematics, 2013, vol. 2013, 38 p. DOI: http://dx.doi.org/10.1155/2013/849231

10. Fineberg J., Bouchbinder E. Recent developments in dynamic fracture: some perspectives. Int. J. Fract., 2015, vol. 196, pp. 33-57.

11. Elices M., Guinea G.V., Gomez J., Planas J. The cohesive zone model: advantages, limitations and challenges. Eng. Fract. Mech., 2002, vol. 69, pp. 137-163.

12. Xu X.-P., Needleman A. Numerical simulations of fast crack growth in brittle solids. J. Mech. Phys. Solids, 1994, vol. 42, pp. 1397-1434.

13. De Borst R., Remmers J.J.C., Needleman A. Mesh-independent discrete numerical representations of cohesive-zone models. Eng. Fract. Mech., 2006, vol. 73, pp. 160-177.

14. Bardenhagen S.G., Nairn J.A., Lu H. Simulation of dynamic fracture with the material point method using a mixed J-integral and cohesive law approach. Int. J. Fract., 2011, vol. 170, pp. 49-66.

15. Agwai A., Guven I., Madenci E. Predicting crack propagation with peridynamics: a comparative study. Int. J. Fract., 2011, vol. 171, pp. 65-78.

16. Javidrad F., Mashayekhy M. A cohesive zone model for crack growth simulation in AISI 304 steel. J. Solid Mechanics, 2014, vol. 6, pp. 378-388.

17. Barenblatt G.I. The mathematical theory of equilibrium cracks in brittle fracture. Adv. Appl. Mech., 1962, vol. 7, pp. 55-129.

18. Willis J.R. A comparison of the fracture criteria of Griffith and Barenblatt. J. Mech. Phys. Solids, 1967, vol. 15, pp. 151-162.

19. Moes N., Belytschko T. Extended finite element method for cohesive crack growth. Eng. Fract. Mech., 2002, vol. 69, pp. 813-833.

20. Lavit I.M. Stable crack growth in an elastoplastic material. Strength of Materials, 1988, vol. 20, no. 7, pp. 854-860.

21. Malik A.V., Belaya L.A., Lavit I.M. About dynamic loading the body with a crack under plane strain. Fundamental and applied problems of engineering and technology, 2016, vol. 315, no. 1, pp. 3-10.

22. Mikhlin S.G. Variational methods in mathematical physics. Oxford, Pergamon Press, 1964, 428 p.

23. Liskovets O.A. The method of lines. Differential equations, 1965, vol. 1, pp. 1308-1323.

24. Ladyzhenskaya O.A. The boundary value problems of mathematical physics. Berlin, Springer-Verlag, 1985, 408 p.

25. Zienkiewicz O.C. The finite element method in engineering science. London, McGraw-Hill, 1977, 541 p.

26. Murakami Y. Stress intensity factors handbook. vol. I. Oxford, Pergamon Press, 1990, 448 p.

27. Chen Y.M. Numerical computation of dynamic stress intensity factors by a lagrangian finite-difference method (the HEMP code). Eng. Fract. Mech., 1975, vol. 7, pp. 653-660.

28. Brickstad B. A FEM analysis of crack arrest experiments. Int. J. Fract., 1983, vol. 21, pp. 177-194.

29. Israil A.S.M., Dargush G.F. Dynamic fracture mechanics studies by time-domain BEM. Eng. Fract. Mech., 1991, vol. 39, pp. 315-328.

30. Lin X., Ballmann J. Re-consideration of Chen's problem by finite difference method. Eng. Fract. Mech., 1993, vol. 44, pp. 735-739.

31. Wen P.H., Aliabadi M.H., Rooke D.P. Application of the weight function method to two-dimensional elastodynamics fracture mechanics. Int. J. Fract, 1996, vol. 76, pp. 193-206.

32. Wen P.H., Aliabadi M.H., Rooke D.P. A contour integral method for dynamic stress intensity factors. Theor. and Appl. Fract. Mechanics, 1997, vol. 27, pp. 29-41.

33. Phan A.-V. A non-singular boundary integral formula for frequency domain analysis of the dynamic T-stress. Int. J. Fract, 2012, vol. 173, pp. 37-48.

i Надоели баннеры? Вы всегда можете отключить рекламу.