Эволюция пластической деформации у затупленной вершины трещины в кристалле

ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

УДК 539.422.23

ЭВОЛЮЦИЯ ПЛАСТИЧЕСКОМ ДЕФОРМАЦИИ У ЗАТУПЛЕННОЙ ВЕРШИНЫ ТРЕЩИНЫ В КРИСТАЛЛЕ*

Д.Н. КАРПИНСКИЙ, С.В. САННИКОВ

(Южный федеральный университет),

Б.В. СОБОЛЬ

Выполнен расчет эволюции пластической деформации у вершины трещины в кристалле а-Ре в условиях плоской деформации (моды I и II). Получены временные распределения пластической деформации, эффективного сдвигового напряжения, коэффициента интенсивности напряжения. Проведено сравнение результатов расчетов для трещин с затупленной и острой вершинами.

Ключевые слова: вершина трещины, пластическая деформация, кристалл, дислокации, плоскости скола, системы легкого скольжения, коэффициент интенсивности напряжения.

Введение. Оценивание влияния формы вершины трещины на характеристики разрушения является актуальной задачей физики прочности и механики разрушения. Выбор модели формы хрупкой трещины существенно влияет на результаты расчета распределения упругого напряжения и деформации в окрестности ее вершины. В простейшем случае трещины-разреза, когда пренебрегают радиусом кривизны вершины р ^ 0, эти распределения можно представить в виде асимптотического ряда по степеням расстояния г(к_1)/2 (к > 0) от вершины до заданной точки [1, 2]. Пренебрежение членами ряда с к < 0 обусловлено условием ограниченности величин перемещений и энергии деформирования у вершины трещины. Расчеты показывают, что распределения напряжений и деформаций содержат лишь особенности в вершине О (г 12) при к = 0, которые могут

исчезнуть при некотором распределении внешней нагрузки [3]. Кроме детального изучения главного члена асимптотики (к = 0), содержащего особенность, в последние годы внимание исследователей привлекает учет несингулярных слагаемых упругого поля вблизи вершины трещины-разреза (к > 1). В работах [4, 5] подробно исследован второй член асимптотики упругого напряжения (Т-напряжение), который в случае трещины-разреза не зависит от г. В этом случае Т-напряжение входит лишь в компоненту тензора напряжения, соответствующую растяжению или сжатию вдоль линии трещины. Расчеты показали, что учет несингулярной составляющей заметно изменяет условия роста трещины [4, 5].

Впервые оценка распределения упругого поля в полярных координатах г и ф у вершины хрупкой трещины с учетом ненулевой величины радиуса вершины р прямолинейной трещины

получена в исследовании [6]. Однако позже была получена более точная формула для распределения компонент тензора напряжений у вершины [7]:

СТг —

К і

Фф :

4л/2тсг

К і

Тгф

4%І2ш

К і

Ф2 (р, г ) СОЭ[ з(р, г ) соэ [ ^

КII

4у[2жг

К іі

4^2жг

4%І2кг

К

Фі^г) яд 12 )+у 2 ^г)яіп Г у

Ф2 ^ г) 8Іп [ 2) + ^4 ^ г) 8Іп [у

Фз (Р,Г ) 8ІП(Ф] + ^5 г) 8ІП^ + 4К= Фз (Р,г )8ІП(ф) + ^6 (Р,Г )ЯІП

Работа выполнена при поддержке Российского фонда фундаментальных исследований (грант 10-08-00839^).

где ф1 (р,г) = 5 -2р-3^ , ф2(р,г) = 3 + 2р + 3(р) , Фз(р,г) = 1 + 2р-3(р) ,

у,(р,г)=1 -7г+2(р), ^2(р.г)=1+2г+2(р) - у(р.?)=1+2г - 2(?) -

у,(р, ? )=3 - 21 ?+15 (р). т,(р, ? )=3+3 ?+125 (р). ув (р, г)=3+9 ? -15 (р).

В формуле (1), аналогично исследованию [7], к 1 и кп - коэффициент интенсивности

напряжения (КИН) эквивалентной трещины, т.е. острой трещины (р^ 0) с той же длиной и ориентацией, как затупленная трещина. Сравнение результатов [6, 7] с результатами расчета методом конечных элементов показало, что с увеличением доли моды II во внешней нагрузке точность аналитических формул для компонент тензора напряжений [7] выше.

Для хрупкой трещины эллиптической формы при смешанном нагружении условие катастрофического роста трещины зависит не от формы трещины, а лишь от радиуса ее вершины р [8], который можно оценить как

р = 16,21 у уСс.

Для а^е удельная поверхностная энергия у^=1,42 Дж/м2 и ос=31,7 ГПа и оценка (2) дает р=0,726 нм (около 3 постоянных решетки). Отсюда следует, что хрупкую трещину в приближении механики сплошной среды, нужно рассматривать как разрез с р^0 [8].

В последние десятилетия обозначился значительный прогресс в создании дислокационных моделей для пластической деформации у вершины трещины [9, 10]. Расчеты взаимодействия дислокации с трещиной-разрезом обнаружили явление экранирования дислокацией ее упругого поля, которое характеризуется снижением КИН по сравнению с КИН трещины без дислокаций. Позже эти расчеты были обобщены для случая клинообразного выреза [11, 12], хотя в этих работах получена корневая особенность распределения напряжения у вершины трещины, соответствующая особенности трещины-разреза. В работе [13] исследовано взаимодействие краевой дислокации с трещиной в форме сплюснутого эллипса. С использованием перехода к трещине-разрезу согласно исследованию [7] получены зависимости КИН от взаимного расположения трещины и дислокации и обнаружены положения дислокации, соответствующие максимальным эффектам экранирования (мода I) и антиэкранирования (мода II). В работе [14] исследованы зависимости КИН модели трещины, состоящей из двух лунок, от положения винтовой и краевой дислокаций относительно вершины. Здесь также получена корневая особенность в распределении напряжения в вершине трещины в форме двух лунок. Следует упомянуть также об оценках величины взаимодействия дислокации с трещиной, в вершине которой находится зона сцепления [15, 16]. В этих работах показано, что если модель трещины с особенностью по напряжениям (высокая величина сцепления) в вершине предсказывает экранирование или антиэкранирование в зависимости от знака вектора Бюргерса дислокации, то при низкой величине сцепления - всегда экранирование вне зависимости от знака вектора Бюргерса дислокации.

Расчет равновесной формы дислокационной пластической зоны у вершины трещины-разреза (см., например, [17]) показал, что компонента напряжения ауу на линии продолжения

трещины достигает максимума на некотором расстоянии от ее вершины. А дислокации, испущенные из вершины, более эффективно экранируют упругое напряжение трещины, чем дислокационные диполи, расположенные вблизи вершины [18]. Другой особенностью распределения дислокаций в пластической зоне является подтвержденное расчетами существование свободной от дислокаций области, прилегающей непосредственно к вершине трещины.

Пластическая деформация у вершины трещины в нагруженном образце существенно затупляет трещину, а это влияет на условия зарождения дислокационных петель из нее [19]. Форма

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

Цель наших исследований -сравнение эволюции пластической деформации у вершины трещины в ОЦК кристалле при учете затупленности вершины трещины, условий нагружения и систем дислокационного скольжения при пластической деформации. В дальнейшем ограничимся расчетами для кристалла а^в (моды I и II) в условиях плоской деформации.

Постановка задачи и метод решения. Рассмотрим плоские задачи об эволюции пластической деформации у вершины трещины длиной 21, расположенной в плоскостях скола {100} и {110} бесконечного ОЦК кристалла с постоянной решетки а (рис.1).

Рис.1. Кристаллографическая схема плоскостей легкого скольжения у вершины трещины

В первом случае выберем направление линии фронта трещины <001> и систему легкого скольжения <111> {110}. В ней пластическая деформация осуществляется путем термоактивированного движения смешанных дислокаций с вектором Бюргерса Ь = а!2 < 111 > . Для трещины в

плоскости {110} направления линий фронтов <111> и <110>, а системы легкого скольжения <111> {110} и <111> {112} соответственно (линии дислокаций параллельны фронтам трещины). Для этих систем скольжения пластическая деформация происходит с помощью тех же дислокаций, как и выше, но во втором случае дислокации чисто краевые. В расчетах учитываются только краевые компоненты смешанных дислокаций.

К плоскостям кристалла у = приложены однородные напряжения растяжения

а уу(0 = а а ($) (мода I) и сдвига % уу(Г) = %а () (мода II), монотонно возрастающие до некоторых значений атах, %тах, достаточных для пластического деформирования кристалла, но недостаточных для роста трещины. После достижения нагрузками аа ($) и %а ($) своих максимальных значений расчет проводился в режиме релаксации напряжения у вершины трещины, которому соответствуют постоянные внешние напряжения. В данном расчете принята пропорциональность двух типов внешних нагрузок аа($) = к%а($), постоянная к является параметром задачи.

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

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

Предположим, что пластические зоны, образующиеся в кристалле, имеют в плоскости 0ху линейные размеры настолько малые по сравнению с I, что правомерно описание пластической деформации у вершины полубесконечной трещины. Плоскости {110} и {112}, пересекаясь с плоскостью 0 ху, образуют на ней два семейства линий скольжения, направления которых заданы

векторами Ъ1к (к=1, 2).

Скорость пластического деформирования, обусловленного движением дислокаций у вершины трещины, дается формулой [20]:

и 0{1 - [ав(г, $)/%0] }

dвk

---------— в 0exp

signe e

квТ 0

где и0 - энергия активации; кв - постоянная Больцмана; т0 - температура; е0 и %0 - постоянные; аке - эффективное напряжение сдвига в к-й плоскости легкого скольжения [20]

аке(г,$) = ак(г,$) -а.Хг,t)signаk(г,$) при |ак(г,$)| >|аДг,$)|,

11,,

ак(г,$) = 0 при |ак(г,$)| <|аДг,()\.

В формуле (4) сдвиговое напряжение в плоскостях легкого скольжения

ак

а

напряжение, препятствующее пластическому сдвигу за счет трения решетки а0 и локального упрочнения материала а f, вычисляемого по формуле:

_ 2

I

к =1

e f~ el

где ei и m - постоянные.

¿ Ів(r,t )|

В формуле (5) a£(r, t) - сдвиговая компонента тензоров напряжения (1) для k-й плоскости легкого скольжения, которую в ряде расчетов дополняли вторым членом асимптотики (Crr = Cr+Tcos2Ф, сфф = сф + Tsin2ф, сгф = хгф-Tsinф^ф), а также силами сцепления берегов трещины у ее вершины; d (r, t) - дальнодействующее упругое напряжение, создаваемое дисло-

кациями одного знака в пластическом зоне

elk

k —1 Dk

где Dk - часть пластической зоны, образовавшейся в верхней полуплоскости в результате скольжения дислокаций вдоль ^.

В формуле (8) Арk(z',t) - плотность эффективных дислокаций, связанных с деформацией вk (r, t) соотношением:

APk^, t) = - вk

bd Sk

a с k (z') определяет напряжение, создаваемое дислокацией в упругой плоскости с полубесконеч-ным разрезом, и задается выражением:

с k

где c ' k (r, r') - компоненты тензора, которые характеризуют собственное поле пары дислокаций, а компоненты тензора с " k (r, r') описывают поле напряжений полубесконечной трещины-разреза, нагруженной на ее берегах усилиями, равными по величине и противоположными по знаку усилиями, создаваемыми в сплошном кристалле напряжениями с ' k (r, r') на месте трещины.

Расчеты показали, что учет формы трещины не вносит существенных изменений в clk (r, t) по сравнению с трещиной-разрезом на расстояниях, больших радиуса вершины (r >р) или размера зоны сцепления.

Уравнения (9)-(17) образуют систему, из которой при начальных

в k

и граничных условиях

с

определяются вk(r,t),ck(r,t).

Данную систему уравнений решали численно с изменяющимся в зависимости от вычислительной ситуации шагом интегрирования At,-. Автоматический выбор шага по времени учитывал

ограничение на максимальную скорость деформации в формуле (3) max (в(r, t))< 0,1 с-1. При

этом ограничении остается справедливой формула (3), которая предполагает термоактивированный механизм движения дислокаций. Метод решения системы (3)-(12) аналогичен методу в работе [20]. В расчетах задавали не Cmax и хтах, а соответствующее значение максимального КИН

x^max x^max

K i и K п .

При расчете эволюции пластической деформации в кристалле a-Fe выбраны следующие значения постоянных: 2l = 103 м, Kmax = 0,2 МПа-м1/2, T0 = 300 К; ц = 86,4 ГПа, v = 0,29 , a = 2,867 -10 10 м, в0 = 101 с1, U0 = 0,9 эВ, xci = 330 МПа, с0 = 18 МПа, ci = 2 ГПа, m = 1 [20].

На основе результатов расчетов эволюции пластической деформации далее была вычислена временная эволюция КИН. В расчете предполагалось, что для КИН трещины имеет место представление [21]:

Ki

В формуле (13) KV(t), KVI(t) - см. формулу (1), а Kp(t) - поправка, учитывающая влияние пластической деформации на КИН:

Kp

k=1 Dk

В формуле (14) кр - КИН трещины, соответствующий действию одной дислокации [21]. В некоторых расчетах выражения (13) дополняли поправкой за счет влияния области сцепления

для моды I к&) = кГ(0 + Re[кp(t)] + КС, где КС = “а0^ [15].

Результаты расчетов и их обсуждение. Проанализируем вначале результаты для модели с закругленной вершиной трещины. Ранее было получено [22], что распределения пластической деформации с учетом Т-напряжений и без них заметно различаются для всех видов ориентации плоскости скола, направления линии фронта трещины и плоскости легкого скольжения в кристалле. Влияние Т-напряжений при всех вариантах расчета приводит к укорочению времени достижения равновесного распределения пластической деформации у вершины трещины в кристалле. В работе [22] для расчета а£ (г, t) использованы компоненты тензора напряжений [см. формулу (1)] при 0 в смещенных узлах г + г0. В данной работе в расчетах учитывался конечный радиус вершины трещины (р> 0) [см. формулу (1)]. Результаты расчетов показали, что форма пластической зоны и распределение пластической деформации в ней заметно зависят от радиуса вершины р, но основные особенности распределения пластической деформации (с учетом

Т-напряжений и без них), полученные в [22], сохранились. На рис.2 показаны равновесные распределения пластической деформации у вершины трещины моды II. Для сравнения на рис.3 приведено аналогичное распределение при р ^ 0. Результаты расчетов также подтвердили выводы [22] о том, что наиболее наглядное представление о влиянии Т-напряжений дают временные зависимости КИН и угла направления роста трещины.

У, МКМ --------------------------------------

3-

-2,4 0 х МКМ 2,4 -2’4 ° .V, МКМ 2,4

а) б)

Рис.2. Равновесное распределение пластической деформации без учета Т-напряжений и радиусом затупления вершины р = 0,12 мкм по сопряженным плоскостям скольжения. Изолинии соответствуют величине сдвиговой пластической де-

1 2 формации (%): а- 8 = 0 - граница пластической области; 1 - 0,036; 2 —0,017; 3 - 0,09; 4 —0,07; б - 8 = 0 - грани-

ца пластической области; 1 —0,035; 2 —0,07; 3 —0,1; 4 —0,17. Плоскости скольжения дислокаций {112}, плоскость трещины {100}, направление линии фронта трещины <001>, мода II

а) б)

Рис.3. Равновесное распределение пластической деформации без учета Т-напряжений и радиусом затупления вершины

р = 0 по сопряженным плоскостям скольжения. Изолинии соответствуют величине сдвиговой пластической деформации

1 2 (%): а - 8 =0 - граница пластической области; 1 - 0,18; б - 8 =0 - граница пластической области; 1 —0,05; 2 —0,1;

3 —0,2. Кристаллографические параметры те же, что и на рис.2

Исследования показали, что эти величины чувствительны также к величине радиуса вершины. Кривые 1 и 2 на рис.4 демонстрируют эффект экранирования вершины трещины для моды I.

а) 6)

Рис.4. Временная зависимость коэффициентов интенсивности напряжения моды I без учета Т-напряжений и радиусом затупления вершины: а- р = 0,12 мкм; б - р = 0 . Кривые 1 и 2 определяют эволюцию к г(1) соответственно без учета пластической деформации и с ее учетом (экранирование вершины дислокациями). Кристаллографические параметры те же, что и на рис.2

а) б)

Рис.5. Временная зависимость коэффициентов интенсивности напряжения моды II без учета Т-напряжений и радиусом затупления вершины: а- р = 0,12 мкм; б - р = 0 . Кривые 1 и 2 определяют эволюцию к г(1) соответственно без учета пластической деформации и с ее учетом (антиэкранирование вершины дислокациями). Кристаллографические параметры те же, что и на рис.2

Кривые 1 и 2 на рис. 5 подтверждают существование «аномальной» эволюции к п, обнаруженное в [22]: при трех ориентациях плоскостей скола и линий фронта трещины при скольжении дислокаций в плоскостях легкого скольжения, образующих угол 60° с плоскостью скола, с учетом Т-напряжений и без них, эволюция к п с учетом пластической релаксации дает более вы-

сокие значения к п по сравнению с эволюцией к п для хрупкой трещины. Данное явление соответствует антиэкранированию вершины трещины.

Кроме полученных результатов расчетов, отметим, что для данной модели (мода I) трещины со сцеплением берегов вблизи вершины качественно повторяют особенности эволюции КИН за счет влияния дислокаций из пластической зоны на КИН, полученные выше. Расчеты в [15] выполнены в предположении «замороженности» дислокаций вблизи вершины, а в нашем случае происходят значительные перераспределения дислокаций у вершины трещины.

Заключение. Затупление трещины и Т-напряжение заметно влияют на равновесное распределение пластической деформации у вершины трещины и на временные зависимости коэффициентов интенсивности напряжения, а также приводят к уменьшению времени достижения равновесия. Подтверждено существование явления антиэкранирования дислокациями вершины затупленной трещины, которая соответствует большим величинам к п при учете пластической релаксации по сравнению с к п для хрупкой трещины. Эта «аномалия» определяется плоскостью скольжения дислокаций и не исчезает при учете Т-напряжения. Результаты расчетов, учитывающих наличие зоны сцепления, не подтверждают оценки [15], которая получена для «замороженных» дислокаций.

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

1. Sinclair G.B. Stress singularities in classical elasticity-I: removal, interpretation and analysis /

G.B. Sinclair // Applied Mechanics Review. - 2004. - Vol.57. - P.251-297.

2. Sinclair G.B. Stress singularities in classical elasticity-II: asymptotic identification / G.B. Sinclair // Applied Mechanics Review. - 2004. - Vol.57. - P.385-439.

3. Chan Y.-S. Dependence of crack singularity on loading functions / Y.-S. Chan, G.H. Paulino,

B.-F. Feng, A. Sutradhar // Mechanics Research Communication. - 2010. - Vol.37. - P.191-197.

4. Fett T. T-stress and stress intensity factor solutions for 2-dimensional cracks / T. Fett. -VDI-Verlag, 2002. - 270 p.

5. Chen Y.Z. A rigorous derivation for T-stress in line crack problem / Y.Z. Chen, X.Y. Lin,

Z.X. Wang // Engineering Fracture Mechanics. - 2010. - Vol.77. - P.753-757.

6. Creager M. Elastic field equations for blunt cracks with reference to stress corrosion cracking

/ M. Creager, P.C. Paris // Int. J. Fract. Mech. - 1967. - Vol.3. - P.247-252.

7. Kulmer G. Influence of the root radius of crack-like notches on the fracture load of brittle components / G. Kulmer, H.A. Richard // Arch. Appl. Mech. - 2006. - Vol.76. - P.711-723.

8. Li Z. The brittle fracture criterion based on the maximum tensile stress on the surface of

blunt crack tip / Z. Li, C. Ji, Y. Li, L.R. Xu // Mechanics Research Communication. - 2007. - Vol.34. -P.472-479.

9. Lin I.-H. Cleavage, dislocation emission and shielding for cracks under general loading /

I.-H. Lin, R. Thomson // Acta metallurgica. - 1986. - Vol.34. - P.187-206.

10. Weertman J. Dislocation based fracture mechanics / J. Weertman. - Singapore: World Scientific, 1996. - 290 p.

11. Ohr S.M. Elastic interaction of a wedge crack with a screw dislocation / S.M. Ohr,

S.-J. Chang, R. Thompson // J. Appl. Phys. - 1985. - Vol.57, N6. - P.1839-1843.

12. Zhang T.-Y. Interaction of an edge dislocation with a wedge crack / T.-Y. Zhang, P. Tong,

H. Quyang, S. Lee // J. Appl. Phys. - 1995. - Vol.78, N8. - P.4873-4881.

13. Li T. The stress intensity factor of an edge dislocation near an elliptically blunted crack tip /

T. Li, Z. Li // Int. J. Fract. - 2007. - Vol.144. - P.45-52.

14. Xie C. Cracking characteristics of mixed mode dislocations near a lip-like mode crack /

C. Xie, Y.W. Liu, Q.H. Fang, M. Deng // Theoretical and Applied Fracture Mechanics. - 2009. - Vol.51. -P.139-143.

15. Bhandakkar T.K. Dislocation shielding of a cohesive crack / T.K. Bhandakkar, A.C. ^ng,

W.A. Curtin, H. Gao // J. Mech. Phys. Sol. - 2010. - Vol.58. - P.530-541.

16. Song J. Dislocation shielding and crack tip decohesion at the atomic scale / J. Song, W.A. Curtin, T.K. Bhandakkar, H.J. Gao // Acta Mater. - 2010. - Vol.58. - P.5933-5940.

17. Qian C.F. Micro-mechanical analysis and TEM study of crack initiation in dislocation free zone / C.F. Qian, W.Y. Chu, L.J. Qiao // Int. J. Fract. - 2002. - Vol.117. - P.313-321.

18. Qian C.F. Simulation of edge dislocation emission from a crack under in-plane shear / C.F. Qian, H.F. Li, Z.J. Jiang, D.C. Fan // Theoretical and Applied Fracture Mechanics. - 2005. - Vol.43. -P.369-375.

19. Fischer L.L. The effect of crack blunting on the competition between dislocation and cleavage / L.L. Fischer, G.E. Beltz // J. Mech. Phys. Sol. - 2001. - Vol.49. - P.635-654.

20. Карпинский Д.Н. Расчет эволюции пластической деформации у вершины трещины и связанные с ней явления / Д.Н. Карпинский, С.В. Санников // Прикладная механика и техническая физика. - 1993. - №3. - С.154-160.

21. Карпинский Д.Н. Расчет КИН с учетом пластической деформации / Д.Н. Карпинский,

С.В. Санников // Письма в ЖТФ. - 1985. - Т.11, №24. - С.1481-1487.

22. Карпинский Д.Н. Влияние Т-напряжений на эволюцию пластической деформации у вершины трещины в кристалле при смешанной моде нагружения / Д.Н. Карпинский, С.В. Санников // Деформация и разрушение материалов. - 2009. - №2. - С.20-24.

Материал поступил в редакцию 18.05.2011.

References

1. Sinclair G.B. Stress singularities in classical elasticity-I: removal, interpretation and analysis / G.B. Sinclair // Applied Mechanics Review. - 2004. - Vol.57. - P.251-297.

2. Sinclair G.B. Stress singularities in classical elasticity-II: asymptotic identification /

G.B. Sinclair // Applied Mechanics Review. - 2004. - Vol.57. - P.385-439.

3. Chan Y.-S. Dependence of crack singularity on loading functions / Y.-S. Chan, G.H. Paulino,

B.-F. Feng, A. Sutradhar // Mechanics Research Communication. - 2010. - Vol.37. - P.191-197.

4. Fett T. T-stress and stress intensity factor solutions for 2-dimensional cracks / T. Fett. -VDI-Verlag, 2002. - 270 p.

5. Chen Y.Z. A rigorous derivation for T-stress in line crack problem / Y.Z. Chen, X.Y. Lin, Z.X. Wang // Engineering Fracture Mechanics. - 2010. - Vol.77. - P.753-757.

6. Creager M. Elastic field equations for blunt cracks with reference to stress corrosion cracking / M. Creager, P.C. Paris // Int. J. Fract. Mech. - 1967. - Vol.3. - P.247-252.

7. Kulmer G. Influence of the root radius of crack-like notches on the fracture load of brittle components / G. Kulmer, H.A. Richard // Arch. Appl. Mech. - 2006. - Vol.76. - P.711-723.

8. Li Z. The brittle fracture criterion based on the maximum tensile stress on the surface of blunt crack tip / Z. Li, C. Ji, Y. Li, L.R. Xu // Mechanics Research Communication. - 2007. - Vol.34. -P.472-479.

9. Lin I.-H. Cleavage, dislocation emission and shielding for cracks under general loading /

I.-H. Lin, R. Thomson // Acta metallurgica. - 1986. - Vol.34. - P.187-206.

10. Weertman J. Dislocation based fracture mechanics / J. Weertman. - Singapore: World Scientific, 1996. - 290 p.

11. Ohr S.M. Elastic interaction of a wedge crack with a screw dislocation / S.M. Ohr,

S.-J. Chang, R. Thompson // J. Appl. Phys. - 1985. - Vol.57, N6. - P.1839-1843.

12. Zhang T.-Y. Interaction of an edge dislocation with a wedge crack / T.-Y. Zhang, P. Tong,

H. Quyang, S. Lee // J. Appl. Phys. - 1995. - Vol.78, N8. - P.4873-4881.

13. Li T. The stress intensity factor of an edge dislocation near an elliptically blunted crack tip / T. Li, Z. Li // Int. J. Fract. - 2007. - Vol.144. - P.45-52.

14. Xie C. Cracking characteristics of mixed mode dislocations near a lip-like mode crack /

C. Xie, Y.W. Liu, Q.H. Fang, M. Deng // Theoretical and Applied Fracture Mechanics. - 2009. - Vol.51. -P.139-143.

15. Bhandakkar T.K. Dislocation shielding of a cohesive crack / T.K. Bhandakkar, A.C. Chang, W.A. Curtin, H. Gao // J. Mech. Phys. Sol. - 2010. - Vol.58. - P.530-541.

16. Song J. Dislocation shielding and crack tip decohesion at the atomic scale / J. Song, W.A. Curtin, T.K. Bhandakkar, H.J. Gao // Acta Mater. - 2010. - Vol.58. - P.5933-5940.

17. Qian C.F. Micro-mechanical analysis and TEM study of crack initiation in dislocation free

zone / C.F. Qian, W.Y. Chu, L.J. Qiao // Int. J. Fract. - 2002. - Vol.117. - P.313-321.

18. Qian C.F. Simulation of edge dislocation emission from a crack under in-plane shear /

C.F. Qian, H.F. Li, Z.J. Jiang, D.C. Fan // Theoretical and Applied Fracture Mechanics. - 2005. - Vol.43. - P.369-375.

19. Fischer L.L. The effect of crack blunting on the competition between dislocation and cleavage / L.L. Fischer, G.E. Beltz // J. Mech. Phys. Sol. - 2001. - Vol.49. - P.635-654.

20. Karpinskij D.N. Raschyot e'volyucii plasticheskoj deformacii u vershiny' treshhiny' i svyazanny'e s nej yavleniya / D.N. Karpinskij, S.V. Sannikov // Prikladnaya mexanika i texnicheskaya fizika. - 1993. - #3. - S.154-160. - In Russian.

21. Karpinskij D.N. Raschyot KIN s uchyotom plasticheskoj deformacii / D.N. Karpinskij,

S.V. Sannikov // Pis' ma v ZHTF. - 1985. - T.11, #24. - S.1481-1487. - In Russian.

22. Karpinskij D.N. Vliyanie T-napryazhenij na e'volyuciyu plasticheskoj deformacii u vershiny'

treshhiny' v kristalle pri smeshannoj mode nagruzheniya / D.N. Karpinskij, S.V. Sannikov //

Deformaciya i razrushenie materialov. - 2009. - #2. - S.20-24. - In Russian.

PLASTIC DEFORMATION EVOLUTION AT BLUNT CRACK TIP IN CRYSTAL

D.N. KARPINSKIY, S.V. SANNIKOV

(Southern Federal University),

B.V. SOBOL

The plastic deformation evolution at the crack tip in a-Fe crystal under the condition of the plane deformation (modes I II) is calculated. Temporal distributions for the plastic deformation, effective shear stress, stress intensity factor are derived. The calculation results for the blunt and acute crack tips are compared.

Keywords: crack tip, plastic deformation, crystal, dislocations, cleavage faces, easy slip systems, stress intensity factor.