Особенности радиационного упрочнения и распределения концентрации точечных дефектов при неравновесной диффузии

УДК 539.2

ОСОБЕННОСТИ РАДИАЦИОННОГО УПРОЧНЕНИЯ И РАСПРЕДЕЛЕНИЯ КОНЦЕНТРАЦИИ ТОЧЕЧНЫХ ДЕФЕКТОВ ПРИ НЕРАВНОВЕСНОЙ ДИФФУЗИИ

© В.В. Красильников, С.Е. Савотченко

Ключевые слова: диффузия; точечные дефекты; радиационное упрочнение; механические свойства; предел текучести; математическая модель; уравнение диффузии; дислокационные трубки.

Работа посвящена анализу особенностей распределения концентрации точечных дефектов, возникающих вследствие неравновесности коэффициента диффузии. Описана эволюция плотности точечных дефектов в облучаемой металлической пластине, неравновесный коэффициент диффузии которых зависит от накопленной дозы. Показано, что предел текучести образца, определяемый распределением точечных дефектов радиационного происхождения, достигает дозового насыщения. Приведено обобщение модели Вуттига-Бирнбаума диффузии по неравновесным дислокационным трубкам. В модели неравновесный коэффициент диффузии точечных дефектов считается зависящим от релаксирующей с течением времени избыточной энергии границ дислокационных трубок. Показано, что увеличение избыточной энергии границ неравновесных дислокационных трубок приводит к росту количества диффузанта, приникающего в образец за одинаковые промежутки времени.

ВВЕДЕНИЕ

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

При моделировании таких процессов использование произвольной зависимости коэффициента диффузии от времени может применяться, например, для управления установлением стационарного распределения примесей в кристалле [1]. Непрерывное изменение коэффициента диффузии с течением времени можно получить, непрерывно меняя температуру кристалла, поскольку коэффициент диффузии зависит от температуры по экспоненциальному закону D(t) = D exp(-E /квтО), где T(t) - температура испытания, Ed - энергия активации диффузии и миграции, D0 = av2, а и v - длина и частота перескоков. Регулируя температуру с течением времени T(t), можно анализировать диффузионно-контролируемые процессы, имеющие большое значение для технологии создания наноструктурных жаропрочных сплавов.

Другим важным аспектом использования диффузи-онно-контролируемых технологий является управление радиационным упрочнением материалов. Неравновесная диффузия точечных дефектов может возникать вследствие накопления дозы облучения. Барьерный механизм радиационного упрочнения непосредственно связывает вклад концентрации различных видов барьеров (точечных дефектов или их кластеров) с изменением предела текучести облучаемого материала [2].

В данной работе приводятся две модели неравновесной диффузии точечных дефектов, в которых анализируется распределение их концентраций.

НЕРАВНОВЕСНАЯ ДИФФУЗИЯ В БАРЬЕРНОМ МЕХАНИЗМЕ РАДИАЦИОННОГО УПРОЧНЕНИЯ

Вклад в прирост предела текучести металлов, связанный с барьерным механизмом торможения дислокаций на определенном типе барьеров, описывается выражением вида [2]:

ЛСТ = Р|0Й

где в - параметр, характеризующий мощность барьера (постоянная величина для данного типа барьера); | -модуль сдвига; Ь - длина вектора Бюргерса; С(х, Г) -объемная плотность точечных дефектов в момент времени Г на расстоянии х от облучаемой поверхности (барьеров данного типа, например, вакансий, или их кластеров); dс - их средний размер. В частности, вакан-сионные и межузельные кластеры имеют средний размер dс < 5 нм, а характеризующий мощность кластерного барьера параметр имеет значение в и 0,2 [2].

Рассмотрим облучаемый образец в виде пластины конечной толщины к. Пусть ось Ох перпендикулярна облучаемой поверхности х = 0, а плоскость уОх совпадает с этой поверхностью. При формулировке модели предполагается следующее [3-4].

1. Толщина пластины существенно меньше линейных размеров облучаемой поверхности, но превосходящая длины пробега как налетающих частиц, так и первично выбитых атомов.

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

3. Вблизи образовавшихся радиационных дефектов в объеме V не могут образовываться новые каскады, приводящие к изменению плотности барьеров за счет рекомбинационных эффектов.

4. Распределение объемной плотности точечных дефектов по толщине облучаемого образца в зависимости от дозы х = Фг определяется одномерным уравнением диффузионного типа:

§=^ дт

д 2C дх 2

+ K (1 - VC),

где К - начальная скорость образования точечных дефектов; V - эффективный объем взаимодействия между ними.

5. В начальный момент времени считается известным распределение объемной плотности точечных дефектов:

C(x, 0) = C0.

6. На поверхностях пластины отсутствуют точечные дефекты радиационного происхождения:

c.(т)=

4

-Л. JD(s)ds

ж

Со + K J é

I. J D(s)ds

dt\

Лп =

(2n +1)2 ж2, . =0, 1, h2

Выражение (7) описывает зависимость объемной плотности точечных дефектов радиационного происхождения от дозы х = Фг и учитывает неравновесность их коэффициента диффузии. Для более детального анализа необходимо конкретизировать зависимость неравновесного коэффициента диффузии. Рассмотрим два случая.

1. Линейная зависимость коэффициента диффузии от дозы. Зависимость параметра неравновесности от дозы приближенно считается линейной и может быть оценена по формуле [5]:

С(0, t) = C(h, t) = 0.

7. Коэффициент диффузии Б в облучаемом материале считается зависящим от дозы х. Влияние дозы на коэффициент диффузии атомов титана в подложке из железа при ионном облучении, в частности, изучено в [5-6]. Для дозовой зависимости коэффициента диффузии в [5] было приведено выражение:

D(t) = De • (1 + ad),

где Б = Б ехР(~^ / квТ) - равновесный коэффициент диффузии точечных дефектов; Т - температура испытания; ЕЛ - энергия активации диффузии и миграции; Б0 = av2, а и V - длина перескоков и собственная частота колебаний атомов при миграции в равновесных условиях; d - не зависящий от дозы температурно-энергетический фактор; а - зависящий от дозы параметр неравновесности. Температурно-энергетический фактор определяется выражением [5]:

где параметр а0 = /Е0, а - полное поперечное

сечение рассеяния, Е5 - пороговая энергия смещения, g - среднее число смещенных атомов на один первично выбитый атом. Из (5) и (9) получается:

D(t) = De • (1 + ат) ,

где а = а(Д

Использование линейной зависимости неравновесного коэффициента диффузии от дозы (10) позволяет провести интегрирование в (6) и получить явное выражение для коэффициента Фурье разложения плотности дефектов (8):

Сп (t) =

4

ж

-KVt-ЛnDe (т+ ат2 / 2) ,

п eV '{(c{) +

+ K

ж e 2ЛnDeа [erfi^j) - erfi (b2)]},

2Л„ Da

d = E0 -Em + kBT ck„T

En

л/ж

где егй (г) = —^ ^ еу ф - мнимая функция ошибок;

где Ет - энергия миграции при отсутствии облучения; Е0 - максимальная энергия, передаваемая частицам от радиационного потока. Считается, что эта энергия не превышает пороговую энергию смещения атомов в твердом теле, которая для большинства металлов порядка 15 эВ. Учет микроскопической неравновесности возможен при выполнении условия аехр(£т /квТ) > 1.

Математическая формулировка модели представляет собой первую начально-краевую задачу для уравнения (2) с начальным условием (3) и краевыми условиями (4). Решение начально-краевой задачи (2)-(4) пред-ставимо в виде ряда:

C(x, т) = ^ Сп (T)sin

n=0

где коэффициент Фурье:

(2 п + 1)жпх

h

KV + ЛnDe (1 + ai).

№nDe,

, = KV + Л nDe

л/2Л nDea

Подстановка (7) и (11) в (1) дает значения предела текучести в различных точках пластины в зависимости от дозы [3-4]. Анализ (7) приводит к выводу, что в любой фиксированной точке пластины зависимость плотности дефектов и предела текучести от дозы не является монотонной. При малых дозах происходит рост предела текучести по степенному закону, согласующийся качественно с экспериментальными данными [7-8]. Далее при определенной дозе наблюдается максимум предела текучести, после которой величина предела текучести начинает медленно уменьшаться. Ограниченность модели линейного нестационарного коэффициента диффузии обусловливает ее применимость только при малых дозах. С увеличением дозы,

é

a0t

a

z

0

а

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

2. Зависимость коэффициента диффузии от дозы с учетом эффекта насыщения. Для параметра неравновесности в [5] была предложена формула

а = оФ? /(1 + оФ?),

которая учитывает эффект насыщения коэффициента диффузии с накоплением дозы. С учетом (12) коэффициент диффузии (5) примет вид:

Б(Г) = Бе (1 + ЛоФГ /(1+стФ?)).

С ростом дозы неравновесный коэффициент диффузии (13) монотонно возрастает и достигает значения

Пе • (1 + а),

при котором наступает насыщение.

Воспользовавшись дозовой зависимостью коэффициента диффузии (13), выражение (8) можно представить в виде:

4е

-(КУ „0„ )т

с" О = ,0 п

%(2п +1) х (1 + от п <|со + К|

I е(кУ пО„ )í

(1 + ст/)м

где |и = ХД^/о. При произвольных значениях параметров интеграл в (15) может быть вычислен только численно.

Подстановка (15) в (7) позволяет получить распределение плотности дефектов в облучаемой пластине, коэффициент диффузии которых учитывает дозовое насыщение:

этих зависимостей выражение (17) удобно представить в виде Лст = ВС , где . Для численно-

го анализа предела текучести (17) использовались модельные значения параметров (они указаны в подписях к рисункам).

Из выражения (17) следует, что предел текучести достаточно быстро достигает насыщения (кривая 3 на рис. 1). Это обусловлено тем, что неравновесный коэффициент диффузии выходит на значение насыщения (14).

Также на форму дозовой зависимости (17) существенно влияет значение температурно-энергетического фактора (6).

При определенных значениях параметров форма кривой дозовой зависимости предела текучести не является монотонной. При накоплении определенной дозы предел текучести достигает максимума, затем немного уменьшается до значения насыщения и далее с ростом дозы уже не меняется (кривые 1 и 2 на рис. 1). Такие формы дозовых зависимостей (17) определяются в первую очередь величиной полного сечения рассеяния, с ростом которого происходит более быстрое достижение величины насыщения предела текучести.

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

С( х, т) =

4е~КУт^ е~'кп°»т (1 + от/п

п=0

т е(КУ+ХпО„ )г

+ К | ---&} БШ

0

1п

(2п + 1)%х И '

{Со +

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

Ло = 2Р|Ае

т

I

КУт/2

( ю „-КО,, т

I

V п=0

е п°»т (1 + от

+ К

Г е(КУ+лп°»)' . (2п + 1)%х41/2

{Со +

(1 + о/)м

Результаты численного анализа предела текучести (17) приведены на рис. 1-2. Для построения графиков

Рис. 1. Зависимости предела текучести (17) от дозы при фиксированных значениях параметров: С0 = 0,01; к = 1; V = 0,7; Де = 1; В = 0,8; к = 10; х = 5; d = 10 и различных значениях сечения рассеяния: (1) - о = 0,5; (2) - о = 1,5; (3) - о = 8 в отн. ед.

Д(Т 0.7^ 0.6 0.5 0.4 0.3 0.2 0.1

0

Рис. 2. Зависимости предела текучести (17) от дозы при фиксированных значениях параметров: С0 = 0,01; к = 1; V = 0,7; Де = 1; В = 0,8; к = 10; х = 5; о = 5 и различных значениях температурно-энергетического фактора: (1) - d = 10;

(2) - d = 15; (3) - d = 25 в отн. ед.

х

0

т

%

т

И

0

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

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

НЕРАВНОВЕСНАЯ ДИФФУЗИЯ ПО ДИСЛОКАЦИОННЫМ ТРУБКАМ

Известно несколько моделей, описывающих диффузию точечных дефектов по изолированным дислокационным трубкам: Р. Смолуховского [9], А. Ле-Клера и А. Рабиновича [10], М. Вуттига и Х. Бирнбаума [11]. В них дислокации рассматриваются в виде цилиндрических трубок радиуса г[, которые проходят через кристалл перпендикулярно его поверхности, на которую нанесен диффузант - тонкий слой примесных атомов, диффундирующих в образце. Диффузия внутри дислокационной трубки подчиняется закону Фика для однородной среды, а также коэффициент диффузии Dd внутри трубки считается много больше коэффициента Б в объеме кристалла.

Рассмотрим обобщение модели Вуттига и Бирн-баума [11-12], в которой неравновесность дислокационных трубок будет учитываться посредством зависимости коэффициента от времени Г [13-17]. При формулировке модели предполагается следующее.

1. Аналогично неравновесным границам зерен [13], избыточная энергия границ неравновесных дислокационных трубок (на один атом): АЕ(/) = АЕ0 exp(-/ / /0), где АЕ0 - избыточная энергия релаксации границ дислокационных трубок; г0 - характерное время релаксации.

2. Неравновесный коэффициент диффузии по дислокационным трубкам связан с избыточной энергией их границ выражением:

5. Уравнение, описывающее распределение концентрации диффузанта в дислокационной трубке при г < г[, при сделанных предположениях записывается в виде:

д2с, 2О дев

^ = О, (/)-- +---

д/ * &2 г* дг

Подстановка (19) в (20) приводит к уравнению при Г < и

дг=о* ^ д/

д 2сг д22

где q2 = 2Б / г21п(гт / г) . При условии Dd >> Б недооценка отвода диффузанта из дислокационной трубки в объем в уравнении (21) при малых временах оказывается несущественной [11-12].

6. В начальный момент времени образец не содержал диффузанта:

Cd(z, 0) = 608(7),

где 8(7) - дельта-функция Дирака (правосторонняя); с0 - постоянная концентрация диффузанта, поддерживаемая на поверхности образца.

7. Так как рассматривается диффузия из постоянного источника, то на входе дислокационной трубки задается граничное условие:

с^0, г) = С0.

8. На выходе дислокационной трубки считается, что диффузант не накапливается, а распределяется по поверхности образца, поэтому следует использовать модель «открытой трубки», в которой на тыльной стороне трубки задается граничное условие:

с£К Г) = 0.

Б* (/) = Бге exp(АE(/)/ квТ),

где - коэффициент диффузии по равновесным дислокационным трубкам; Т - температура отжига.

3. Концентрация внутри трубки с[ считается независящей от расстояния г от оси симметрии дислокационной трубки радиуса г[ [11]. Пусть ось 07 перпендикулярна поверхности кристалла в глубину от поверхности в плоскости 7 = 0, в которой находится слой диффузанта, а плоскость 7 = к является тыльной стороной образца толщины к. Дислокационные трубки проходят через весь материал.

4. Радиальный градиент концентрации диффузан-та в зерне считается не зависящим от времени, концентрация диффузанта в зерне cg определяется выражением при г > г[ [11]:

cg = СГ 1п

где гт - половина расстояния между дислокациями.

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

Решение начально-краевой задачи (21)-(24) представляется в виде разложения в обобщенный ряд Фурье:

(2, /) = С0 Л - 2 - - X Рк (/^1П1

к=1

где рк (/) = 1 е-?2/-(як / Н)21 (/) |1 + ? 21 е?2х+(А / Н)21 (х)*х|, /

I(/) = |О(эГ .

0

Выражение (25) описывает распределение концентрации диффузанта внутри дислокационной трубки. С использованием (18) можно получить:

г

С

г

Н

0,30

0,25 -

0,20 -

0,15

0,10

0,05

0,00

cd 1 c0

Г

1 (J^""

1 /

t / to

0

1 2 3 4 5

Рис. 3. Типичные зависимости (25) при фиксированных h = 10, z = 3, q2 = 0,8 в отн. ед. и при различных значениях ДЕ>: (1) - АЕ0 = 0; (2) - ЛЕ0 = 0,7квТ; (3) - АЕ0 = 1,5ВД (4) -АЕ0 = 2,2квТ

10 8 6 4 2 0

0

1

2

3

Рис. 4. Типичные зависимости (28) при фиксированных h = 10, q2 = 0,8 в отн. ед. и при различных значениях ЛЕ0: (1) - ЛЕ0 = 0; (2) - ЛЕ0 = 0,7квТ; (3) - ЛЕ0 = 1,5квТ

I (t) = Det0 < Ei

ЛЕ 0

V kBT

- Ei

ЛЕ

V kBT

0 g- / io

где Е^) - интегральная экспонента.

Подстановка данного выражения в (25) приводит к распределению концентрации диффузанта в дислокационной трубке с учетом неравновесного коэффициента диффузии (рис. 3). Из анализа (25) с учетом конкретных зависимостей (18) и (26) следует, что с ростом избыточной энергии релаксации неравновесных дислокационных трубок ЛЕ0 происходит увеличение концентрации диффузанта на фиксированной глубине в данный момент времени (рис. 3).

При больших значениях ЛЕ0 зависимость концентрации диффузанта от времени не является монотонной, а в определенный момент времени Гт = Гт(ЛЕ0) достигается максимальное значение концентрации диффузанта на определенной глубине. С дальнейшим течением времени Г > Гт концентрация диффузанта начинает уменьшаться и достигает определенного значения, которое не зависит от ЛЕ0 и представляет собой стационарную концентрацию при диффузии по равновесным дислокационным трубкам.

В рамках модели «открытой трубки» количество диффузанта, прошедшего за время Г через образец и собравшегося на единице площади тыльной поверхности кристалла, выражается через его концентрацию cd [11-12]:

Q(t) = -nrdfPd J

Dd

dcd

dz

dT,

где р, - плотность дислокаций.

Подставка (25) в (27) приводит к выражению для количества диффузанта, прошедшего через образец [16], в явном аналитическом виде (рис. 4):

Q(t) =

= Qо (0 + 2| (т)1 (-1)4 (т>Г* 2т-(Л / И)2 7 (т)Лт

I 0 к=1 j

т

где Ак (т) = 1 + д21 е^^ / И)2 7 ; Qо = го^С) / И.

0

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

ЗАКЛЮЧЕНИЕ

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

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

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

t

h

z

0

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

ЛИТЕРАТУРА

1. Панкратов Е.Л. // ЖТФ. 2004. Т. 74. Вып. 1. С. 115-119.

2. Физика радиационных явлений и радиационное материаловедение / Паршин А.М. и др. Белгород, 1998. 378 с.

3. Колобов Ю.Р., Красильников В.В., Савотченко С.Е., Прозорова М.С. // Актуальные проблемы прочности: материалы 47 Меж-дунар. конф. Н. Новгород, 2008. С. 36-38.

4. Красильников В.В., Савотченко С.Е., Прозорова М.С. // Материаловедение. 2009. № 11. С. 59-63.

5. Степанов В.А. // ЖТФ. 1998. Т. 68. № 8. С. 67-71.

6. Ковалева М.Г. Автореф. дис. ... канд. физ.-мат. наук. Белгород: БелГУ, 2006. 152 с.

7. Неклюдов И.М., Камышанченко Н.В. // Научные ведомости БГПУ. 1996. № 1. С. 18-33.

8. Зеленский В.Ф., Неклюдов И.М. // ВАНТ. Сер.: ФРП. 1984. Вып. 1 (29). С. 46-73.

9. Smoluchowski R. // Phys. Rev. 1952. V. 87. P. 482.

10. Le Clair A.D., Rabinovitch A. // J. Phys. C: Solid State Phys. 1981. V. 14. P. 3863.

11. WuttigM., Birnbaum H.K. // Phys. Rev. 1966. V. 147. P. 495-505.

12. Kaur I., Mishin Y., Gust W. Fundamentals of Grain and Interphase Boundary Diffusion. Chichester: John Wiley & Sons Ltd, 1995. 512 p.

13. Назаров А.А. // ФТТ. 2003. Т. 45. Вып. 6. С. 1112-1114.

14. Красильников В.В., Савотченко С.Е. // 18 Петербургские чтения по проблемам прочности и роста кристаллов: в 2 ч. СПб., 2008. Ч. 2. С. 128-130.

15. Красильников В.В., Савотченко С.Е. // Физико-математическое моделирование систем: материалы 5 Междунар. семинара. Воронеж: ВТГУ, 2008. Ч. 1. С. 131-136.

16. Красильников В.В., Савотченко С.Е. // Известия РАН. Серия Физическая. 2009. Т. 73. № 9. С. 1348-1354.

17. Савотченко С.Е. // Известия ТулГУ. Естественные науки. 2009. № 3. С. 204-218.

Поступила в редакцию 23 декабря 2014 г.

Krasilnikov V.V., Savotchenko S.E. PECULIARITIES OF RADIATION HARDENING AND CONCENTRATION DISTRIBUTION OF POINT DEFECTS UNDER NON-EQUILIBRIUM DIFFUSION

The paper discusses the analysis of features of the distribution of the concentration of point defects arising from the non-equilibrium diffusion coefficient. An evolution of the density of point defects in the irradiated metal plate, a non-equilibrium diffusion coefficient which depends on the cumulative dose was described. It is shown that the yield stress of the sample determined by the distribution of point defects caused by radiation reaches the dose saturation. The generalized Wuttig-Birnbaum model of diffusion on non-equilibrium dislocation pipes is introduced here. In the model the non-equilibrium diffusion coefficient of point defects is considered dependent on relaxing excess energy over time of boundaries of dislocation tubes. It is shown that an increase in the excess energy of boundaries of non-equilibrium dislocation tubes leads to an increase in the number of diffusing, penetrating into the sample at regular intervals.

Key words: diffusion; point defects; radiation hardening; mechanical properties; yield stress; mathematical model; diffusion equation; dislocation tubes.

Красильников Владимир Владимирович, Белгородский государственный национальный исследовательский университет, г. Белгород, Российская Федерация, доктор физико-математических наук, профессор кафедры «Материаловедение и нанотехнологии», старший научный сотрудник, e-mail: kras@bsu.edu.ru

Krasilnikov Vladimir Vladimirovich, Belgorod State National Research University, Belgorod, Russian Federation, Doctor of Physics and Mathematics, Professor of "Materials and Nanotechnology" Department, Senior Scientific Worker, e-mail: kras@bsu.edu.ru

Савотченко Сергей Евгеньевич, Белгородский юридический институт Министерства внутренних дел Российской Федерации, г. Белгород, Российская Федерация, доктор физико-математических наук, доцент, профессор кафедры «Математика», e-mail: savotchenko@hotbox.ru

Savotchenko Sergey Evgenyevich, Belgorod Law University of the Russian Interior Ministry of Russian Federation, Belgorod, Russian Federation, Doctor of Physics and Mathematics, Associate Professor, Professor of "Mathematics" Department, e-mail: savotchenko@hotbox.ru