ЭМПИРИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСПРЕДЕЛЕНИЯ ТВЕРДОСТИ α – ТВЕРДОГО РАСТВОРА АЗОТИРОВАННОГО СЛОЯ | Опубликовать статью ВАК, elibrary (НЭБ)

ЭМПИРИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСПРЕДЕЛЕНИЯ ТВЕРДОСТИ
α – ТВЕРДОГО РАСТВОРА АЗОТИРОВАННОГО СЛОЯ

Научная статья

Крукович М.Г.^1, *, Федотова А.Д.²

¹ ORCID: 0000-0001-5563-093Х;

^{1, 2} Российский университет транспорта (МИИТ), Москва, Россия

* Корреспондирующий автор (ya.bormag[at]yandex.ru)

Аннотация

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

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

EMPIRICAL SIMULATION OF THE HARDNESS DISTRIBUTION OF AN α-SOLID SOLUTION
OF A NITRIDED CASE

Research article

Krukovich M.G.^1, *, Fedotova A.D.²

¹ ORCID: 0000-0001-5563-093X;

^{1, 2} Russian University of Transport (MIIT), Moscow, Russia

* Corresponding author (ya.bormag[at]yandex.ru)

Abstract

The study establishes the influence of alloying elements on the kinetics of layer phase growth and hardness based on numerous nitriding experiments in gas, liquid, and solid media of various types of structural steel. The article develops and tests an empirical model for predictive calculation of the distribution of nitrogen content and hardness over the depth of the α – solid solution of the nitrided case. The model provides for the application of the thermodynamic coefficient of nitrogen activity and its dependence on the total effect of the content of alloying elements. The research presents comparative results of numerical experiments for calculating the maximum hardness and its distribution over the depth of the α-layer with real experimental data on various types of steel. The calculation that is based on the presented model shows a higher agreement with the experiments than the calculation methods recommended in the literature.

Keywords: nitrided case, empirical model, hardness, α-phase, alloying elements, structural steel, coefficient of thermodynamic activity of nitrogen, predictive calculation.

Введение

Твердость нитридных фаз азотированного слоя, имеющих ограниченную толщину (0,01 – 0,04 мкм) в зависимости от условий обработки и материала подложки изменяется в достаточно узком интервале и составляет 1000 – 1300 HV [1], [3], [4]. В то же время она может значительно уменьшаться в зависимости от пористости поверхностного слоя. Основную часть азотированного слоя на конструкционных сталях составляет α-твердый раствор с дисперсными включениями нитридов, который и определяет основные эксплуатационные характеристики слоя, сжимающие напряжения на поверхности упрочненных деталей и плавное уменьшение твердости от поверхности к сердцевине детали.

Твердость диффузионной части азотированного слоя представляет собой сумму упрочнения твердого раствора (∆_т.р.) и упрочнения дисперсными частицами (∆_д.ч.). Следовательно, степень упрочнения определяется концентрацией азота в слое и геометрическими параметрами выделяющихся в матрице нитридных фаз. Упрочнение при формировании диффузионной зоны, является результатом торможения дислокаций включениями нитридов, поэтому для расчета целесообразно использовать модели Мотта – Набаppo и Орована. Теория Мотта н Набарро применима для случаев когерентного выделения частиц второй фазы. В этом случае при деформации частицы перерезаются движущейся дислокацией. Для матрицы с некогерентными частицами второй фазы различают три механизма дислокационного взаимодействия: для случая тонких выделений – механизм перерезания частиц, механизм огибания частиц (по Оровану) и механизм переползания (по Хиршу) [2], [5], [9].

Проведенные расчеты [10] показали, что максимальное упрочнение азотированного слоя можно получить как за счет когерентных, так и некогерентных выделений путем соответствующей оптимизации структуры при реализации расчетных показателей: диаметра выделений и дистанции между ними. В то же время расчет суммарной твердости азотированных легированных сталей по вышеупомянутой методике приводит к получению результатов с ошибкой ~ 15% [5].

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

Сущность разработанной модели

Моделирование распределения твердости по толщине слоя α-твердого раствора конструкционных сталей в данной работе проводилось в два этапа. На первом этапе рассчитывалось влияние легирующих элементов (до 6 % по массе) на твердость α-твердого раствора в зоне, примыкающей к зоне нитридов на расстоянии ~ 0,025 мм. А на втором этапе находили формулы влияния легирующих элементов на распределение твердости от поверхности вглубь стальной основы с учетом коэффициента термодинамической активности азота.

Твердость азотированного слоя в начале α – твердого раствора представляется как сумма приращения твердости слоя за счет азота и легирующих элементов, и твердости подложки (термически обработанной упрочняемой стали) [5], [11].

    (1)

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

Влияние легирующих элементов в стали на приращение твердости

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

Суммарное значение твердости α – слоя при воздействии совокупности некоторого количества легирующих элементов определялось по формуле:

    (2)

 

Таблица 1 – Влияние содержания легирующих элементов на приращение твердости α – слоя

	Легирующий элемент	Расчетная формула	Достоверность, R2
1	Ni	5,1095x3 – 54,795x2 + 192,96x + 92,079	0,991
2	Si	47,643x + 150	0,9967
3	Mn	72,558x + 153,41	0,9991
4	V	387,68x0,4038	0,9962
5	Cr	8,1988x3 – 94,883x2 + 385,87x + 105,29	0,9968
6	Al	4,193x3 – 65,062x2 + 375,31x + 101,12	0,9975
7	Mo	19,038x3 – 213,94x2 + 715,2x – 2,4302	0,9815
8	W	3,2338x3 – 35,937x2 + 128,76x + 97,625	0,9909
9	Ti	8,2836x3 – 97,2x2 + 399,83x + 101,85	0,9968

Примечание: х – содержание соответствующего легирующего элемента

 

Расчет коэффициента термодинамической активности азота

Коэффициент термодинамической активности азота (γ^N) представляет собой меру отклонения реального раствора от равновесного состояния. Он тесно связан с термодинамической активностью азота и с распределением его по толщине слоя. При этом по мере увеличения расстояния от поверхности наблюдается снижение концентрации азота и, как следствие, повышение коэффициента активности и приближение его к единице. Следовательно, он характеризует способность к поглощению азота конкретной сталью или сплавом и уменьшается от поверхности до сердцевины детали.

      (3)

где С_N – максимальное содержание азота в слое; С₀ – равновесное содержание азота в α – железе;  – термодинамический коэффициент активности азота в железе;  – термодинамический коэффициент активности азота в легированной стали;  – максимальное содержание азота в легированной стали.

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

Например, для содержания хрома в интервале до 5% по массе формула имеет вид [3]:

      (4)

В этом случае термодинамический коэффициент активности азота от содержания хрома составит:

 при 3% Cr, а концентрация азота в этом случае составит:

    (5)

 

Таблица 2 – Влияние легирующих элементов (% по массе) на термодинамический коэффициент
активности азота в стали

	Легирующий элемент	Логарифм активности азота: lgγ(N)	Концентрационный интервал, %
1	Cr	– 0,54*(% Cr)0,5	0 – 5
		– 0,0067(137+%Cr)	5 – 13
2	Si	0,5* (% Si)	0 – 3
3	Mo	– 0,29* (% Mo)	0 – 5
4	Mn	– 0,21* (% Mn)	0 – 3
		– 0,28* (% Mn)	0 – 5
5	V	– 0,895* (% V)0,4	0 – 2,5
6	C	0,41* (% C)	0 – 1,0
7	Al	0	Теоретическое влияние
		– 0,5*(% Al)	Расчетное влияние
8	Ni	0,019* (% Ni)	0 – 5
9	Ti	– 0,588* (% Ti)	0 – 2
10	Zr	– 0,22* (% Zr)	0 – 0,15
11	Nb	– 0,59* (% Mo)	0 – 0,35
12	W	– 0,118* (% W)	0 – 8,8
13	S	0,8*(%S)	S > 0,1
14	Pb	0,7*(%Pb)0,5	Pb > 0,1

 

Термодинамический коэффициент активности азота стали от ряда легирующих элементов определяли по формуле:      (6)

В соответствии с полученным значением определяется максимальное содержание азота в α – твердом растворе

     (7)

Концентрация азота по толщине α – слоя в конкретной точке (y_i) определяется по формуле, полученной из решения уравнений Фика и кинетическим закономерностям роста фаз слоя:

      (8)

где:  – минимальная концентрация азота в стали, не вызывающая изменение структуры и свойств;

 y_α – толщина слоя α – твердого раствора, определенная металлографическим методом;

 y_нитр – толщина слоя нитридов.

Сравнение рассчитанных значений твердости в начале α – слоя с экспериментальными данными (Таблица 3), полученными при обработке в расплавах солей, показал хорошее соответствие.

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

       (9)

Тогда твердость слоя в точке (i) составит:

       (10)

 

Таблица 3 – Сравнительные данные экспериментов и расчета твердости

Марка стали	Твердость на расстоянии 25мкм от поверхности, кгс/мм2		Погрешность, %
	Эксперимент	Расчет
10	280	297	6,07
35	375	350	6,9
80	405	438	8,2
35ХГМ	650	660	1,5
30Х2МНЮА	840	821	2,2
22ГФС	610	605	0,8
40Х5МФС	1190	1169	1,9
Х	700	649	7,2
30Х3М	880	865	2,9
35Х2Н4М	685	678	1,2
35ХГН2	630	621	1,2
А45Г2	450	446	0,9
АС15Г	315	318	0,9

Примечание: условия обработки: 570 ^оС; 1,5 ч; 38 % CNO^–

 

Полученные расчетные результаты численного эксперимента (Рис. 1) показывают хорошее совпадение с экспериментальными данными, а представленная эмпирическая методика позволяет легко разработать алгоритм расчета.

 

Рис. 1 – Сравнение результатов реальных и численных экспериментов:  Рассчитанные значения;
 Экспериментальные значения

Примечание: условия обработки: 570 ^оС; 1,5 ч; 38 % CNO^–

 

Заключение

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

Конфликт интересов Не указан.

Conflict of Interest None declared.