ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2015. № 2 (46)
Материаловедение
УДК 681.5
АНАЛИЗ МЕТОДОВ РЕШЕНИЯ ЗАДАЧ СТРУКТУРООБРАЗОВАНИЯ ПРИ АЗОТИРОВАНИИ
Т.А. Бенгина
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
Проведен анализ численных методов решения задачи структурообразования при газовом азотировании. Выделены преимущества и недостатки комбинированного метода, основанного на построении конечно-разностной аппроксимации по времени и применении метода конечных элементов по пространственной координате, и метода с преобразованием координат, трансформирующим изменяющуюся во времени область в стационарную. Установлено, что распределение азота по толщине поверхностного слоя существенно влияет на эксплуатационные характеристики изделия и технико-экономические показатели процесса, поэтому использование численных методов решения задачи помогает прогнозировать структуру азотированного слоя при разных технологических режимах.
Ключевые слова: газовое азотирование, тепломассоперенос, численное решение, метод конечных элементов, конечно-разностная аппроксимация, многофазная диффузия.
Для изменения структуры и свойств поверхностных слоев металлов или сплавов часто используется химико-термическая обработка (ХТО). В результате такой обработки изделий повышается их износостойкость, упрочняется поверхностный слой, увеличивается сопротивление коррозии и предел выносливости. Азотирование как один из видов ХТО заключается в диффузионном насыщении поверхностного слоя металла или сплава азотом путем выдержки при определенной температуре в активно газовой среде. Как следствие этого, изменяются его структура и свойства.
Наиболее существенное влияние на показатели качества обрабатываемого изделия при азотировании оказывает именно распределение азота по глубине упрочняемого слоя. Этот процесс описывается краевой задачей тепломассопере-носа [1, 2].
Глубина азотированного слоя значительно меньше толщины детали, поэтому математическое описание процесса определяется дифференциальным уравнением диффузии для бесконечной пластины толщиной 2R:
дС(х,т) __д_
дт дх1
' (0 дСхт)'
дхг-
(1)
Татьяна Алексеевна Бенгина (к.т.н), доцент кафедры «Высшая математика и прикладная информатика».
где т> 0, 0 < x < R, D¿ (C) - коэффициент диффузии азота в i-той фазе, C(x, г) -концентрация азота в точке х в момент времени т .
Начальное значение концентрации азота в слое C0 (x) и положение границы раздела фаз Si0 (x) произвольно:
C(x,0) = Co( x); (2)
Si (x,T)| т=0 = #ю( x). (3)
Граничное условие на правой границе:
- Di
dC( x, т)
dx
= д (п N -п N),
(4)
x=R
где П n =
P
NH 3
P
- азотный потенциал печной атмосферы; Pnh , Ph - парциаль-
H 2
ные давления аммиака и водорода соответственно; П N - равновесный азотный потенциал соответствующей фазы; Д- - коэффициент массопереноса для --той фазы.
Граничное условие на левой границе:
ас (х, г)
dx
= 0.
(5)
x=0
(6)
Условия слева и справа на границах фаз:
С(х Г х- = С-- ;
С(х ,г)| = Сг+.
Условие сохранения баланса вещества при фазовых переходах определяет положение межфазной границы:
Л
д£ дт
= Di -dC
dx
Su
- Di-, -dC
dx
(7)
Si-
где ji _ С- - С*, г _ 1, N, N - количество фаз [3, 4].
В условиях конкретного производства весьма затруднительно выбрать определенный технологический режим обработки изделия с целью придания ему конкретных эксплуатационных характеристик. Сложность процесса азотирования, его долгая протяженность по времени не позволяют однозначно выбрать необходимый режим. Поэтому прогнозирование процессов формирования структуры и свойств в диффузионной зоне при азотировании является актуальной задачей.
При разработке и освоении технологического процесса азотирования предварительное моделирование и использование численных методов решения задачи структурообразования позволяют оперативно располагать информацией о формировании диффузионного слоя и скорости образования новых фаз при различных технологических режимах. Это в итоге позволяет сокращать сроки проектирования и освоения процесса обработки изделий, повышает качество и стабильность азотирования, увеличивает производительность оборудования.
Известно достаточно много методов решения задач многофазной диффузии. Чаще всего для их решения используются методы с выделением фронтов [3, 4, 5]. Такие расчетные методы имеют высокую точность, но при условии возникно-
вения или рассасывания фаз становятся практически нереализуемыми из-за сложности алгоритма.
Методы сквозного счета позволяют решать такие задачи в теории теплопроводности, однако не могут быть использованы для решения неизотермических задач многофазной диффузии. В связи с этим для решения модели позднее использовался способ, занимающий промежуточное положение между указанными выше методами. Такой подход основан на приведении вариационного соотношения к виду, удобному для построения конечно-разностной аппроксимации по времени, и применении метода конечных элементов по пространственной переменной к возникающим на каждом временном слое краевым задачам [3, 4].
При решении задачи (1)-(7) на первом шаге было получено вариационное соотношение
у Нт) зс й у ^т) д( п ас ^
2 ] д~Лйх = у ] Л—I п Т I' (8)
г=0 £.( т) дТ = й(т) дс V дс У
где л(х) - произвольная функция из Н(0, Ь).
Вместо пространственно-временных конечных элементов применяется МКЭ только по х, а перед разностной аппроксимацией производных по времени используется формула для интегрирования по частям правой части (1) и формула
Чт) зс , й Чт) „ , „ ¿т *к
^ Т = йТ ^ СЛ " ^-0 + (9)
Й(т) Й(т)
N / \ N N Й+1 ,
Vп ддс>=Улп дс -у\ >£*. (10)
гг -1 дх V дх У гг дх гг -1 йх дх
г=0 Йл 1 й 1=0 й
С учетом (9), (10) соотношение (1) перепишется в виде
А N й+1 N йй N йй
—У \сфх-Ус/ й +Уй +
У йл дс . N дс У ] -т-п 1Тёх-Улп -Г-
~=0 й йх дх 1=0 дх
й +1
= 0.
Й
Окончательно получим
й Ь N Й+1 й дс
— |слйх + У | й-П —йх + (Р0Л(0) + Л(Ь) +
Т 0 г=0 й х х (11) N й£. + Ул(Й )-р (с+- с-- ji) = 0.
Пространство пробных функций Н 1(0,1) аппроксимируется конечномерным пространством VN функций, непрерывных на [0, Ь] и линейных в каждом из интервалов (х1;хш). Набор функций Л'-'-'Лу} из VN , в т. ч. лк(х^) = ¿>к], образует базис в VN .
Подставим л = Лк в интегральное тождество (11) и проведем неявную раз-
ностную аппроксимацию по времени:
—^ + Е ] А 17 + ^к|х_° + Фн Л к | х_ь +
-_0 & V (12)
+ £% Е- )(с+- с--- л _
здесь й - шаг разностной схемы по времени; 1к _ |, а ^ означает
о
^/г-й.
Уравнения МКЭ в матричной форме получаются после вычисления и суммирования вкладов конечных элементов :
А(&)С _ Я(£), (13)
где Е_ ), С _ (С1,...,CN)т.
Вычисление концентраций С в узлах сетки (13) и определение положения границ Е путем интегрирования уравнений (7) на каждом шаге по времени осуществляются итерационно.
Если оказывается, что для некоторого - , 1 < - < N +1, имеет место
&-1+1-* > 5+1), считается, что одна из фаз рассосалась. В этом случае одна или несколько границ выбрасываются из рассмотрения, а оставшиеся перенумеровываются. Затем вычисляются Л(&(^+1)), £(Е('5+1)) и система (13) решается относительно С (^г1). После сходимости итераций производится процесс проверки соответствия получившихся значений С1 и CN имеющимся на концах [°;ь] фазам. Если соответствие нарушилось, в точке х _ 0 или х _ Ь устанавливается новая межфазная граница.
Данный метод имеет ряд достоинств: отслеживание появления и рассасывания фаз, фиксирование скачков концентраций на границах, учет переменного числа фаз и границ. Однако при этом коэффициенты диффузии и азотные потенциалы рассматриваются как линейные функции концентрации азота, что приводит к значительному упрощению задачи. Невозможность явно выделить положения границ раздела фаз, а также узнать скорости их перемещения можно также отнести к недостаткам рассмотренного метода. Учесть подобные недостатки помогает метод, основанный на преобразовании координат, которое переводит изменяющуюся во времени область в стационарную. Этот подход впервые был предложен Ландау.
В задаче (1)-(7) рассмотрим изолированную произвольную фазу, заключенную между границами и , и введем безразмерные координаты:
х - Е
х-- ^ >Ъ _ Е Е-1 , 0 <Л< 1. (14)
Е -Е-1
Уравнение (1) в новых координатах примет вид
ас аг
где
(л- -1) %
ах ах
7
ас ал- _ а
аЛ- ах- аЛ-
ас
в- (с )т7
у
чах- у
1
, (15)
(16)
Для крайних фаз уравнение (15) имеет тот же вид, только одна из границ имеет нулевую скорость. Соответственно преобразуются граничные условия: - для последнего уравнения
В
N '
ас а^
ЛN _1
ас 7 7 ах
_Р-1 ^
N
Л- _1
а^ ах
- для первого уравнения
Т У
Скорости у- _ —L находятся на каждом шаге по времени из уравнений (10), йг
_ 0.
71 _0
которые после преобразования имеют вид
В- (с а
ах 7 а
- В-_1 (с)
ас а77i_l
Е-+0
а7i_l
(17)
Е--0
Появление в уравнениях (15) членов, аналогичных конвективным, требует соблюдения консервативности при выборе численной аппроксимации. С учетом преобразования координат расчетная область остается постоянной. А для решения таких задач можно уже использовать численные методы с фиксированным числом узлов в каждой фазе.
В результате конечно-разностной аппроксимации система уравнений (15) примет вид
^с + Л& • А?с • _ £„ (вс). (лхл)2 .
Граничные условия:
Ас АЛЛ _0 _0;
ВN Алс • АхЛ1 _, _ А • - );
ЧН-1
л •АЕ _ в- .Алс 'АхЛ
1 - В--1 А с • АхЛ-1 0>
(18)
(19)
(20) (21)
где Аг, АЛ, Ах - дифференциально-разностный оператор по г, 77, х соответственно.
Для решения системы уравнений (18-21) используется неявный метод с конечным числом итераций, количество которых зависит от порядка системы уравнения.
К преимуществам данного метода следует отнести возможность моделирования процесса с произвольным числом одновременно существующих фаз и их чередования в области [0, ь]. Для этого на очередном временном шаге осуществ-
ляется проверка на появление фазы и выбираются новые граничные условия.
Т У
Расчет концентрации С(х,т) и скоростей уг- =—- производится методом про-
йт
гонки. С полученными значениями скоростей проводятся итерации с целью уточнения скоростей межфазных границ:
- V + V/ V =-Чт~;
£ = 4 + V/ -т, где £, V/ - значения на предыдущем временном слое.
В конце итерационного процесса необходима проверка на исчезновение фаз. Если область существования некоторых фаз оказалась меньше некоторой заранее заданной точности, то из системы уравнений (15) исключаются уравнения, соответствующие таким фазам. Соответственно преобразовываются расчетные области.
Предложенный метод численного решения задач многофазной диффузии позволяет отследить положение межфазной границы, скорость ее передвижения, при этом возможно моделирование переменного числа фаз и их толщины, разрывы концентраций на границах также учитываются. Это является несомненным достоинством при решении задач структурообразования при азотировании, так как учитываются именно те факторы, которые в конечном итоге влияют на свойства поверхностного слоя изделия.
По итогам исследования образцов технического железа были получены данные моделирования. Результаты численных расчетов качественно согласуются с экспериментальными данными, что свидетельствует об адекватности рассмотренных моделей расчетов [3].
В таблице приведены данные расчетов по первому и второму численному методу решения задачи многофазной диффузии (в числителе дроби) и экспериментальные (в знаменателе).
Расчет толщины азотируемого слоя*
№ метода Температура, °С Время процесса, ч Азотный потенциал, ат1/2 Толщина слоя, мкм
а г' е
1 540 9 2,64 300/260 0,7/1-2 1,5/1-2
2 540 9 1,98 600/550 4/4-5 8/7-10
*Материал - техническое железо.
Рассмотренные численные методы решения задачи многофазной диффузии при азотировании позволяют получать различное число фаз, фиксировать их появление и рассасывание, учитывать скачки концентраций на границах. Однако первый из методов затрудняет явное выделение межфазной границы и скорость ее перемещения. Поэтому в случае необходимости точного отслеживания положения межфазных границ целесообразно использовать метод решения, основанный на подстановке Ландау.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Лахтин Ю.М. Регулирование фазового состава и содержания азота в нитридном слое при азотировании стали 38Х2МЮА // МиТОМ. - 1996. - № 1. - С. 5-11.
2. Лахтин Ю.М., Коган Я.Д., Булгач А.А. Азотирование в машиностроении // Сб. научных трудов МАДИ. - 1986. - С.42-49.
3. Бенгина Т.А. Оптимизация технологического процесса газового азотирования: Дисс. ... канд. техн. наук. - Самара: Самар. гос. техн. ун-т, 2008. - 155 с.
4. Коган Я.Д., Булгач А.А. Моделирование на ЭВМ кинетики диффузионного насыщения при газовом азотировании // МиТОМ. - 1984. - № 1. - С. 10-19.
5. Лившиц М.Ю. Теория и алгоритмы оптимального управления термодиффузионными процессами технологической теплофизики по системным критериям качества: Дисс. ... докт. техн. наук. - Самара: Самар. гос. техн. ун-т, 2001. - 425 с.
Статья поступила в редакцию 16 февраля 2015 г.
ANALYSIS METHODS FOR SOLVING PROBLEMS OF STRUCTURE FORMATION DURING NITRIDING
T.A. Bengina
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
The paper carries out the analysis of numerical methods for solving the problem of structure formation during gas nitriding. The paper highlights the advantages and disadvantages of the combined method based on the finite-difference time approximation and the use of finite element method by the spatial coordinate, and the method with coordinate transformation which transforms a time-varying part into stationary. It was found that the distribution of nitrogen throughout of the surface layer significantly affect the performance of the product and the technical and economic parameters of the process, so the use of numerical methods for solving the problem helps predict the structure of the nitrated layer at different technological modes.
Keywords: gas nitriding, heatmass transfer, numeral decision, method of eventual elements, finite difference approximation, multiphase diffusion.
Tatiana A. Bengina (Ph.D. (Techn.)), Associate Professor.