CC BY
362
УДК 621.1+536.2
А.Г.Князева
Институт физики прочности и материаловедения СО РАН (Томск)
МОДЕЛИРОВАНИЕ ФИЗИЧЕСКИХ И ХИМИЧЕСКИХ ЯВЛЕНИЙ В ПРОЦЕССАХ ОБРАБОТКИ ПОВЕРХНОСТЕЙ МАТЕРИАЛОВ ВЫСОКОЭНЕРГЕТИЧЕСКИМИ ИСТОЧНИКАМИ
Abstract
In this paper the traditional conceptions and some effects concerned with irreversible processes models are described. The possibilities of the surface treatment model taking into account the previously applied powdered layer shrinkage and the model of carbon diffusion between covering and substrate in the cooling down condition are presented in a short way.
Под высокотемпературными технологическими процессами зачастую понимают такие технологические процессы, которые реализуются при использовании концентрированных источников энергии. Типичными примерами высокотемпературных технологических процессов являются лазерная и электронная технологии сварки, резки, термической обработки; дуговая сварка и другие способы соединения материалов; плазменные технологии нанесения покрытий и поверхностной обработки; ионные технологии; кислородная резка; совмещенные технологии резки, сварки, наплавки; процессы получения тонких пленок и выращивание монокристаллов; процессы химической и диффузионной обработки поверхностей материалов и др. К высокотемпературным технологиям можно отнести и многие технологии получения новых материалов в химической промышленности; переработки и сжигания природного топлива; различные металлургические процессы. Каждый из названных технологических процессов, в свою очередь, включает множество частных технологий, зависящих от конкретного технического решения, условий и материалов, для которых процесс предназначен. Современные методы обработки поверхностей и нанесения покрытий связаны с преобразованием различных видов энергии в тепловую, с процессами нагрева и охлаждения и разнообразными сопутствующими явлениями, которые и приводят к изменению состава, структуры и свойств обрабатываемых слоев. Преобразования структуры в высокотемпературных технологических процессах происходят на разных стадиях, включают различные уровни (мезо-, микро- и нано-) и зависят от разнообразных физических явлений, которые характеризуются разными масштабами и протекают совместно. Традиционно математическое моделирование в этой области сводится к расчету формы эффективных источников энергии, процессов плавления и кристаллизации в условиях высокоэнергетических воздействий, явления распыления поверхностного слоя или гидродинамики расплавов. Но формирование структуры и свойств поверхностных слоев происходит непосредственно в процессе обработки и является следствием разнообразных необратимых физико-химических процессов, не учитываемых в традиционных моделях. На традиционных представлениях, некоторых эффектах, связанных с учетом в моделях необратимых процессов, а также на описании возможностей этих моделей мы далее и остановимся.
Эффективный источник тепла
Число способов модификации поверхностей велико, каждый из них имеет свои особенности (рис. 1). Более того, под общим названием скрываются различные технологические модификации. Тем не менее технологические процессы поверхностной обработки высокоэнергетическими источниками имеют ряд общих черт, что и послужило в свое время базой для развития теплофизических основ этих процессов. Для всех высокотемпературных технологий характерны существенная неравновесность процессов, связанная с неоднородным распределением температуры и ее изменением во времени; высокие скорости нагрева и охлаждения различных элементов системы; наличие сложного теплообмена; существование нескольких различных фаз; соотношение между которыми изменяется; разнообразные физикохимические явления, сопутствующие нагреву и охлаждению или лежащие в основе технологии. Но, в первую очередь, во всех высокотемпературных технологиях происходит преобразование различных видов энергии в тепловую энергию и (или) ее непосредственное использование для получения, переработки и модификации материалов и их поверхностей.
Рис. 1. Упрощенная классификация методов обработки поверхностей
Использование разных источников, связанных с видами воздействия, приводит к разнообразным задачам теплопроводности с различными начальными и граничными условиями. При анализе тепловых процессов реальные источники заменяют идеализированными, форма которых в той или иной степени приближена к фактической. В зависимости от характера воздействия выделяют источники объемные и поверхностные, неподвижные и движущиеся, распределенные и точечные. Такая классификация является достаточно условной и берет свое начало в задачах сварки и резки [1,2], но может быть использована и для задач поверхностной обработки.
„ Один и тот же источник энергии может быть
отнесен к разным типам. Например, характер распределения энергии в электронном луче может быть описан кривой Гаусса (рис.2) с максимумом в объеме [3]:
q(x, y, г) =
г) = ^ехр
- (х 2 + у 2 )к - (г - И)2 к2
где ^ максимальная плотность мощности источника тепла (определяемая силой тока и разностью потенциалов), ^ -
коэффициент сосредоточенности потока электронов (аналогично лазерному излучению), к - положение максимума тепловыделения в объеме материала. Можно определить эффективный радиус электронного луча
*1=14^1 и эффективную область рассеяния энергии или
проникновения электронов в объем *2 = 1 ■'[^2 . Тогда при условии Хт = ^К1* >> *2, где
к = Х/(ер) - коэффициент температуропроводности, г* - характерное для
Рис. 2. К определению вида источника
технологического процесса время (например, время наблюдения), хт - характерный тепловой масштаб, источник энергии, связанный с электронным лучом, можно считать поверхностным. Это вполне подходит для электронов с энергиями 10-500 Эв, когда
плотности мощности не превышают 105 ^ 106 Вт/см2. Но это же совершенно не годится для быстрых электронов, а тем более для релятивистских, энергия которых превышает 1 Мэв.
Если эффективный радиус электронного луча много меньше характерного размера Г0 поверхности, на которую он действует, можно говорить о точечном источнике. Если же *1 ~ *0, где *0 - эффективный размер обрабатываемой
поверхности, источник тепла становится распределенным, т.е. зависящим от координат.
В технологиях обработки поверхностей материалов, так же, как в процессах сварки и резки, большую роль играют характер перемещения детали и режим сканирования электронного луча [3,4], что приводит к понятию подвижного источника энергии.
Например, при обработке поверхности детали, перемещающейся вдоль оси Ох со скоростью V, сканирующим электронным лучом (рис.3) часто можно использовать функцию вида (в случае пилообразных колебаний луча)
/(X У, г) = ехР
у 2 + (х - Уі)2 (г - И)2
Я?
Я2
х < Их
Рис. 3. Движущийся вдоль
поверхности источник тела
где кх - полуширина сканирования. Для поверхностного источника к = 0 .
С увеличением сложности характера тепловыделения возрастает и число параметров, необходимых для описания характера воздействия.
Иного вида источники тепловыделения, связанные с внешним воздействием, имеют место в процессах ионной имплантации и плазменного напыления [5].
Геометрия детали
Наряду с характером внешнего воздействия при постановке задач поверхностной обработки материалов большую роль играет геометрия детали [2,6,7]. При формулировке моделей процессов обработки в декартовой системе координат среди тел классических форм выделяют стержни, пластины (параллелепипеды, у которых толщина много меньше ширины и длины) и трехмерные тела. Если речь идет о термической обработке тел вращения (сканирующим электронным лучом, излучением лазера и т.п.), то целесообразно перейти к цилиндрической системе координат. В этом случае уравнение теплопроводности принимает вид
дТ
ср = -У-3 + д(г,ф,г,г), 3 = -ХУТ , (1)
ді
где
-пг дТ і ді ді ^ т і д/г\ і д д
^Т = ^Т“ег +-1-ЄФ + ^“ег, V' ^ = —:т(г^г )_|-1-^Ф + Т",
1 дТ
дТ
ід
ід
д
дг г дф дг г дг г дф дг
Jr,Jф- компоненты вектора плотности теплового потока, ег ,Єф,ег - единичные векторы цилиндрической системы координат; с, р, X - теплофизические характеристики (теплоемкость, плотность и теплопроводность), зависящие от температуры и координат. В общем случае мы приходим к трехмерным задачам теплопроводности. Учет конкретных особенностей технологий позволяет существенно упростить задачу или, наоборот, учесть в ней более «тонкие» вещи.
Например, при термической обработке поверхностным источником энергии боковой поверхности полого цилиндра (рис. 4), внутренняя поверхность которого теплоизолирована, уравнение теплопроводности (1) дополняется граничными условиями (в (1) источник тепла в этом случае отсутствует),
г = Я1: — = 0, (2)
дг
г =
^2 : Х-^ = д(р,г.г)-а(Т -Те)-с в ГГ
дг '
г = 0, Ь :
4 7^4
- Тк
).
дТ
дг
= 0
(3)
(4)
(иногда требуется учесть теплообмен излучением на внутренней поверхности цилиндра).
В случае цилиндра с покрытием толщиной к мы приходим к сопряженной задаче теплопроводности: уравнение (1) без источника, в котором все величины имеют индекс «1», имеет место для основной детали, а с индексом «2» - для покрытия. Условие (3) с Х = Х2 записывается при г = Я + к, а при г = Я имеет место условие идеального теплового контакта,
цилиндра
Х
Г
дх
= Х
Я2 -0
дТ_
' дх
Я2 +0
% -0=Г
Я2 + 0
(5)
Если покрытие можно считать термически тонким по сравнению с основой,
проинтегрируем уравнение (1) по координате г в пределах от Я до Я + к с учетом
условий (3), (5). Тогда придем к задаче, включающей уравнение теплопроводности (1) без источников со свойствами основы ( с,Р1, Х1), условия (2), (4) и
дТ / ч _ \ _4 , дТ
Х1 = д(т, г, г) - а(Т - Те ) - С в Т4 - С2р2к 8(
к
д
. дТ
Хо ——
д
дТ
+ к— I Х 2 — дг 1 2 дг
(2 + к) др^др д I • (6)
Для металлов, обладающих высокой теплопроводностью, представляет интерес еще одно приближение, которое корректно для тонкостенных труб с покрытиями. Интегрируя (1) по всей толщине трубы от Я до ^2 + к, придем к двумерному уравнению
дТ
[1 (2 - Я1) + С2Р2к] = дг
д (, дТ л Хр^р
др
-^(хг —
дг 1 г дг
+ )/(ф, г, г) - с в Т4 -а(Т - Те ), (7)
где
Х =Х К2 -К1 +______________
Ф 1 Я (2 + к)2
Хг = Х1(я2 -Я1 )+Х2к ,
2 1 2 2
которое, вообще говоря, описывает анизотропный процесс теплопроводности.
Если источник равномерно распределен по толщине цилиндра, что может представлять интерес для плазменных технологий, то от (7) легко перейдем к одномерному уравнению.
Задачам теплопроводности с подвижными и неподвижными источниками тепла, посвящалось и посвящается множество публикаций в журналах ИФЖ, «Теплофизика и аэромеханика», ТВТ, «Сварочное производство», «Heat- and Mass Transfer» и другие, поскольку они имеют широкие приложения в различных областях физики.
Не секрет, что процессы поверхностной обработки материалов высокоэнергетическими источниками сопровождаются, в первую очередь, плавлением материала, а остывание деталей за движущимся источником и после обработки непосредственно связано с его кристаллизацией. Экспериментатору в подобных исследованиях интересны не только величина максимальной температуры, но и размер и форма ванны расплава, размер зоны термического влияния (ЗТВ), достигаемые в тех или иных режимах обработки. С одной стороны, размер ванны расплава можно оценить условно по величине температуры плавления, используя решения задач теплопроводности, в том числе аналитические. Но поглощение тепла в результате плавления может сказаться на всем балансе тепла, на величине максимальной температуры, поэтому необходим учет явлений плавления и кристаллизации в динамике.
Из физики и термодинамики известно, что в окрестности фазовых переходов первого рода теплоемкость вещества резко возрастает, что отражает зависимость
где Тт - температура плавления (кристаллизации), Ьт - скрытая теплота плавления,
дельта-функция Дирака. При численном решении задач теплопроводности с фазовыми переходами дельта-функция заменяется дельта-образной функцией при неизменном энергетическом балансе. Удобной в этом отношении является функция
где 5 - полуинтервал сглаживания. Экспериментальные данные по теплофизическим характеристикам веществ вдали от температуры плавления аппроксимируются полиномами второй или третьей степени.
В выборе параметров 5 -образной функции существуют свои тонкости, которые необходимо принимать во внимание. Если все сделано корректно, то на температурной кривой в окрестности температуры плавления должно появиться плато, причем оно будет заметно и на кривых температура-время в фиксированных точках, и на пространственных распределениях в различные моменты времени.
В качестве простого примера рассмотрим задачу о термической обработке материала с покрытием при условии малости их общей толщины движущимся источником энергии. Простейшая математическая постановка задачи в этом случае может быть представлена в виде
Плавление и кристаллизация
CLPL,Т^^
ср =
CsPs,Т < ^
(8)
индекс «Ь» относится к жидкости (к расплаву), индекс «5 » - к твердой фазе; 5(х)
x = О, Hx : — = О, у = О,Hy : — = О ,
cx у ду
t = О: Т = ТО,
(9)
где
/(х, У, г) = ехр[- кх (у2 + (х - Vt)2 )] |х| < Их,
кх - ширина сканирования, Нх, Ну - размеры пластины в направлениях х и у, -
толщина основы, эффективные теплофизические свойства вычисляются аналогично предыдущему (осредняются по толщине покрытия с основой), а теплоемкости основы С\Р\ и покрытия С2Р2 меняются в соответствии с (8).
Конечно, эта же модель справедлива, если покрытия нет. Примем в этом случае, что пластина изготовлена из железа (сь = 0,825, с1Ь = 0,608 Дж/(г-К), р^ = 7,61, рц = 7,87 г/см3, А,1 = 0,45 Дж/(см-с-К), Тт1 = 1810 К, Ьт1 = 247 Дж/г). В любом случае через некоторое время после начала движения источника значение максимальной температуры перестает изменяться, что видно из рис. 5, а. Плато на температурных кривых соответствует температуре плавления. Плато имеется на кривой и перед источником, но вследствие большого градиента температуры оно слабо выражено. Кроме максимальной температуры можно построить пространственные распределения, что иллюстрирует рис. 6. Анализируя данные расчетов Т(х,у,t), можно исследовать эволюцию ванны расплава и зоны термического влияния (ЗТВ), которую определяют по условно заданной температуре, ниже которой в материале не наблюдается видимых изменений в структуре. Соответствующие рис. 6 форма ванны расплава и ЗТВ показаны на рис. 7 темным и светло-серым цветом. Первая из изолиний на всех рисунках отвечает температуре Т = 350 К. Внешний источник движется вдоль оси Ох (как это показано на рис. 3).
а
б
Рис. 5. Распределение температуры вдоль оси движения источника при обработке
4 2
железной пластины дд = 7 -10 Вт/см с учетом плавления (а) и без учета плавления t = 1. - 0,02;2. - 0,1;3. - 0,5;4. - 0,9;5. -1,7;6. - 2,5;7. - 3,3;8. - 4,1;9. - 4,9;10. - 5,7;11. - 6,5
Если плавление явно не учитывается, то можно, конечно, построить границу зоны термического влияния и предположительной ванны расплава, но вследствие того, что величина максимальной температуры теперь стала иной (рис. 5,б), результаты такого исследования также изменятся. Состав фаз, определенный на основе диаграмм состояния по рассчитанной в тепловой задаче температуре, в любом случае будет весьма приближенным. Дополнительное усложнение тепловых моделей связано с учетом зависимостей теплофизических свойств от температуры, что также приводит к численным изменениям температуры и размеров ЗТВ и ванны расплава.
42
Рис. 6.Эволюция температуры во времени дд = 7 -10 Вт/см при нагреве железной пластины с учетом плавления: t = (а) -1,3; (б) -1,7; (б) - 2,5; (г) - 3,3; (д) - 4,9; (е) - 6,5 с
а б в
где Рис. 6. Эволюция во времени ванны расплава (темно серый цвет) и ЗТВ (светло серый цвет) при нагреве железной пластины с учетом плавления. В незакрашенной области температура
изменяется от 900 до 350 К: д0 = 7 -104 Вт/см2 t = (а) -1,3; (б) -1,7;
(б) - 2,5; (г) - 3,3; (д) - 4,9; (е) - 6,5 с
Если осуществляется нагрев подвижным источником материала с покрытием, то вследствие дополнительных затрат тепла численные значения температуры изменяются. Цель подобных исследований заключается в анализе зависимости максимальной температуры, размеров ванны расплава от толщины и типа покрытия при варьировании технологических параметров.
Это, пожалуй, все, что могут дать чисто теплофизические задачи.
Качественно результаты не изменяются и при переходе к трехмерным моделям. Дополнительные возможности имеются лишь для определения глубины проплавления и распределения температуры вследствие объемного тепловыделения.
Замечу, что при исследовании процессов плавления-кристаллизации появляется класс задач теплопроводности с подвижными границами, который представляет интерес для различных областей физики. Выделение подвижных границ раздела фаз в задачах поверхностной модификации не всегда целесообразно, так как возникающие при этом математические и алгоритмические трудности [8] не стоят получаемых при этом результатов.
В случае неподвижных внешних источников (например, импульсная электронно-лучевая и лазерная обработка) для процессов нагрева и охлаждения тел, в том числе многослойных, в первом приближении подходят классические аналитические методы. Целый ряд полезных для технологий обработки задач теории теплопроводности с решениями содержится, например, в [9-11]. Некоторые из задач воспроизведены в [12,13].
Физико-химические превращения в зоне обработки
Усложнение (усовершенствование) моделей поверхностной обработки может идти по пути учета разнообразных физико-химических явлений, в результате которых и происходит формирование свойств покрытий. В технологиях различного типа явления, лимитирующие формирование свойств, - различны.
Например, в некоторых ситуациях для придания нужных свойств предварительно нанесенному порошковому покрытию требуется термическая обработка, приводящая к уплотнению этого слоя или к усадке. С физической точки зрения усадка связана с процессами плавления и массопереноса в различных фазах, с появлением сил сцепления между отдельными частицами (адгезионной или химической природы). Все эти процессы - термически активируемые и происходят на микро- и мезоуровнях, т.е. в масштабе отдельно взятых частиц: в объеме частиц и на границах их раздела. Макрокинетическое описание процессов, приводящих к усадке, в первом приближении можно свести к поиску кинетического закона, параметры которого зависят от процессов, протекающих на «низших» структурных уровнях.
В рамках представлений, развитых в [14-16], кинетику усадки описывает уравнение
йк2
(
0,
ъ
V 2 шт
-1
ехр
( Е.
Ъ2 > Ъшт;
(10)
Ъ2 >Ъшт,
1: = 0 : ъ2 = ъ20 ,
где Ъ2шп - минимально возможная толщина покрытия, к2 - константа скорости усадки, Е2 - энергия активации процесса усадки (ее величина определяется процессом, лимитирующим усадку на микроуровне), р - показатель степени в законе усадки, который может быть как целой, так и дробной величиной (зависит, например, от характера поверхностей частиц). Все эти параметры, с одной стороны, -
экспериментально определяемые, а с другой - кинетика процессов спекания является очень развитой областью физики, в рамках которой для нахождения кинетических констант можно формулировать и исследовать специальные частные подмодели.
Для иллюстрации дополнительных возможностей, которые появляются в модели процесса поверхностной обработки, предположим, что на поверхность материала, изготовленного из железа (свойства которого выписаны выше), нанесен слой порошка толщиной ^2 = ^20 = 0,1 см. Примем его свойства следующими: X2 = 0,08 Вт/(см-К), с2з = С2Ь = 1,25 Дж/(г-К), Р52 =РL2 = 3,75 г/см3, Тт2 = 2318 К, Ьт2 = 389 Дж/г.
Тогда модель (9), (10) в первом приближении позволит описать изменение рельефа поверхности вследствие усадки порошкового слоя в процессе поверхностной
4 2
обработки. Так, полагая ^ = 7 -10 Вт/см , получим следующую картину.
Во-первых, температура в этом случае оказывается несколько ниже, чем в предыдущем случае, так как потребовались дополнительные затраты тепла на прогрев покрытия. Во-вторых, температура плавления материала покрытия не достигается, т.е. усадка происходит в «твердой» фазе. В-третьих, иной получается и форма температурных кривых, поскольку эффективные теплофизические характеристики в процессе усадки изменяются:
(ср)е/ = С1Р1 + С2Р2 Н2 (тVН1 , Хв// = Х1 + Х2 Н2 (тУН1 . (11)
Изменение рельефа поверхности со временем показано на рис. 7 и рис. 8 где представлены распределения толщины порошкового слоя в различные моменты времени. Видно, что следствием малой по длительности начальной нестационарной стадии процесса обработки является неравномерная толщина покрытия у торца х = 0. Неравномерная толщина вдоль оси Оу , очевидно, связана с неравномерным нагревом и потерями тепла вследствие теплопроводности. Чтобы выровнять поверхностный слой, требуется многопроходная обработка или большая ширина сканирования, что и используется в технологиях обработки. Влияние ширины сканирования видно из сравнения верхней и нижней строк (рис. 8). Параметры для расчета (Е2 = 100 Дж/моль,
^2шт = 0,03 см, ^2 = 5 -10 см/с) намеренно выбраны такими, чтобы влияние было наглядным. Наблюдаемая экспериментально усадка может достигать 30-40 %.
а б
Рис. 7. Распределение температуры и толщины покрытия вдоль оси движения
4 2
источника в различные моменты времени: ^0 = 7 -10 Вт/см , V = 0,5 см/с, Л1 = 0,03, Ну = 0,05 см
а
б
в
Рис. 8. Толщина покрытия в различные моменты времени. а-в д0 = 7 -10 , Ну = 0,5; г-е до = 7 -104, Ну = 1,5 ; I = а, г - 4,1; б, д - 5,7; в, е - 7,3 с
В трехмерном варианте модели поверхностной обработки с учетом усадки порошкового слоя [15,16] обнаруживаются подобные закономерности. Но для корректного описания наблюдаемых процессов в выбранной системе необходимо учесть перемешивание частиц в расплаве, образующемся под порошковым слоем, и поглощение энергии электронного луча в объеме.
Технология поверхностной обработки сканирующим электронным лучом может
существенно отличаться от описанных простых ситуаций. Например, в [17] для
модификации поверхностей с целью придания им высоких прочностных свойств
используют целенаправленное изменение состава. В формирующуюся ванну расплава
подают модифицирующие частицы, химическая природа и свойства которых могут
быть различны. Если частицы - инертные и тугоплавкие, то в результате обработки
образуется дисперсно-упрочненный поверхностный слой, или фактически
композиционное покрытие. Если частицы могут растворяться в расплаве, то в
зависимости от технологических параметров после обработки получается либо
гомогенное покрытие, либо композиционное, причем состав матрицы также получается
различным, что непосредственно связано с процессом растворения. Математические
модели для различных вариантов такого процесса обработки, позволяющие
предсказывать и описывать эти эффекты, представлены в [18-20] и включают
кинетическое уравнение для доли частиц
р ( \ ( Е„
= Жр - кр Ф(р )ехр
&
__р
ЯГ
у
где пр - объемная доля частиц в расплаве, Жр - источник частиц, а формальнокинетические параметры зависят от конкретной физической ситуации. В теории растворения [21] считается, что энергия активации Ер определяется либо диффузией в
расплаве, либо химическими превращениями между частицами и основой, а константа растворения кр целиком зависит от гидродинамической картины течения расплава (в
моделях растворения эта зависимость считается известной). Кинетическая функция
ф(пр) опять же определяется процессами, протекающими на микроуровне. Для ее
нахождения имеется целый ряд моделей, часть из которых описана в [22, 23].
Частицы оказывают влияние и на эффективные теплофизические свойства, что учитывается в моделях. В случае растворимых частиц изменяется состав матрицы, что приводит к изменению температуры плавления и к необходимости использовать для описания процессов плавления и кристаллизации более сложную теорию (см. раздел 5).
Картина существенно усложняется, если процесс растворения сопровождается химическими реакциями в объеме. Вариант модели для этой ситуации описан в [24].
Существуют технологии поверхностной обработки, сочетающие в себе достоинства энергии электронного луча и синтеза в конденсированной фазе -фактически осуществляется управляемый синтез покрытия на подложке [25]. Организация процесса синтеза также может быть различной. Чисто теплофизические модели здесь более не подходят, и на помощь приходит аппарат макрокинетики [26]. Задача теплопроводности дополняется системой уравнений химической кинетики, а в уравнении теплопроводности появляется слагаемое, учитывающее суммарное тепловыделение в химических реакциях. Модели такого типа позволяют описать эволюцию фазового и химического состава покрытий в ходе обработки, дать состав покрытия к концу процесса, оценить роль химических источников энергии [27-30].
Но даже при таком усложнении (усовершенствовании) моделей ряд реально наблюдаемых явлений остается за кадром.
Перераспределение легирующих элементов
Не секрет, что свойства покрытий и характер их связи с основой зависят от диффузионных процессов различной природы: от перемешивания частиц в расплаве, от диффузии растворенных элементов, от перераспределения примесей в процессах обработки (в ходе нагрева и остывания). Перераспределение элементов приводит к формированию градиентной структуры поверхностного слоя, что важно для придания материалам нужных свойств.
Так, в процессе электронно-лучевой наплавки покрытия с высоким содержанием углерода на подложку из малоуглеродистой стали переходная зона формируется за счет твердофазной диффузии углерода в подложке от поверхности подложка-покрытие. В результате происходит обогащение приграничной области подложки углеродом. При формулировке модели в [31] предполагается, что медленный процесс диффузии углерода из наплавленного покрытия в подложку происходит на стадии остывания покрытия, а перераспределением элементов между слоями непосредственно в ходе наплавки можно пренебречь. Это условие обеспечивается малостью скорости диффузии по сравнению с теплопроводностью в твердых телах.
Пусть после прохождения электронного луча по поверхности обрабатываемого материала (Бе) (в процессе нанесения покрытия) на ней остается «дорожка» расплавленного присадочного материала (Бе + 0,8% С), в процессе остывания которого и происходит образование переходной зоны (рис. 9). Полагаем, что в каждом сечении обрабатываемой поверхности, перпендикулярном направлению движения электронного луча, условия формирования диффузионной зоны одинаковы. Процессы распространения тепла и массы (легирующего элемента) в подложке опишем обычными уравнениями теплопроводности и диффузии
к
Дп >Тг
ср
дТ
X
2 и
7(0
7
Рис. 9. Иллюстрация к математической постановке задачи: 1 - присадочный
материал, 2 - подложка
срИ -дТ1 = И
ді
ді да 2
= 4^2
дх і дх
+ -
= А( о
дх і дх
+ -
ду
д
X
ду
о
да
2
\
(12)
(13)
д1 дх ^ дх ) ду ^ ду )
Полагаем, что толщина слоя присадочного материала к много меньше, чем толщина подложки. Это позволяет не рассматривать распределение температуры и концентрации в нем. Тогда в качестве граничного условия для температуры при у = 0 будет справедливо дифференциальное уравнение
,4 гг 4'
-1 \ дП Т + Х ГдТ2 ^
дх ч дх у=-0 1 ду у у = +0
-Сов
Ті4 - Те 4 ) ,
(14)
позволяющее определить поток тепла, поступающего в подложку, при этом очевидно
Т1 = Т2»у = °.
Аналогичные условие имеют место для концентрации легирующего элемента:
и =иіА о ®0і
ді і дх дх
+ о
у = -0
^да2 ^
ду
(15)
а1 = а2.
Ось оу является осью симметрии, поэтому здесь справедливо условие
дТі = 0 даі дх ’ дх
На бесконечном удалении от нагретой области источники и стоки тепла и массы отсутствуют,
дТі „ даі
■ = 0, і = 1,2, х = 0.
(16)
X,-
дх
■ = 0 , О,——— = 0, і = 1,2, х ^ да ;
дх
X2 ^ = 0, О2а'2
(17)
, . , „ (18)
ду ду
Если вдоль оси Оу обрабатываемый образец имеет конечные размеры, то на его поверхности, противоположной нагретой, для температуры ставится условие типа (14) для к=0, а для концентрации легирующего элемента - условие (16).
В начальный момент времени 1 = 0 имеем
а1 = а°, Т1 = тз; а2 = ^ т2 = /2 (х> у) . (19)
В пользу такого приближения говорит и тот факт, что вследствие конвективного перемешивания в процессе наплавки температура расплавленного слоя и концентрация в нем выравниваются.
В (12) - (19) а( - концентрация легирующего элемента, индекс / = 1 относится к слою присадочного материала, индекс / = 2 - к подложке; Б - коэффициент диффузии, зависящий в общем случае от температуры и концентрации легирующего элемента, Те - температура окружающей среды, а0 - начальное значение концентрации
легирующего элемента в расплаве, Т - его начальная температура; к = к°, если |х| < Я , и к = 0 при |х| > Я, Я - половина ширины сканирования электронного луча. Остальные
обозначения - аналогичны предыдущему.
Начальное значение температуры расплава и начальное распределение температуры в подложке определяются параметрами внешнего источника и могут быть оценены
д
независимо либо на основе данных эксперимента, либо с помощью простейших аналитических зависимостей [3,2].
Эта модель позволяет определить фазовую структуру переходной зоны, распределение концентраций в ней и форму зоны термического влияния, формирующиеся в процессе остывания покрытия и подложки.
В отличие от простейшей модели, описанной выше (в разделе 3 этой статьи), здесь при моделировании процесса плавления использованы представления теории двухфазной зоны [32].
Для теплоемкости, плотности и коэффициента теплопроводности приняты зависимости
с ^ •(!-п) + сь-П + б-Г-П,
Р = Р 5-(1 -П)+Р ь-П, (20)
X = X 5 -(1 -п) +А,ь 'п, где п - доля жидкой фазы, определяемая из соотношения
П = 1 —
n
(21)
1 liq 1 sol
Tliq, Tsol - температуры ликвидуса и солидуса на диаграмме состояния Fe - C [33], n -варьируется для разных сплавов и зависит от содержания примесей; Q - скрытая теплота плавления, индекс "s" - относится к твердому состоянию, ’L - к жидкости. Температуры ликвидуса Tliq и солидуса Tsol можно аппроксимировать линейными
зависимостями
Tsol = as • a + bs , Tliq = ai • a + bl , (22)
T — 0 T — 0
где bs = bs =0, as = —----, aL = -L0---, Ts0 , Tlq - температуры солидуса и ликвидуса
a0 a0
при содержании углерода в стали ао, 0 - температура плавления (кристаллизации) чистого железа (т.е. для а = 0 ).
В случае плавления чистого железа изменение теплоемкости в окрестности температуры плавления описывается зависимостью, аналогичной (8). Для теплоемкости a-фазы справедливо такое же выражение, включающее свои параметры.
Величины полуинтервала сглаживания и теплоту структурного фазового перехода у о а находим, аппроксимируя экспериментальные данные [34] по теплоемкостям в окрестности температуры 0j. Имеем:
са = 0,447, cY = 0,608, cL = 0,825 Дж/(г.К); р s = 7,87, р L = 7,61 г/см3; Ts0 = 1643,
Tl0 = 1763, 0 = 1810, 01 = 1042 К; Q = 247,1 Дж/г, a0 = 0,8%.
В некоторых случаях вместо (20) - (22) удобнее использовать аппроксимацию данных эксперимента [34,35] по теплоемкости чистого железа в виде
с = -0,422 + 0,005•T-8,651-10-6 • T2 + 5,159-10-9T3, T <01; (23)
с = 115,636-0,312• T + 3,158-10-4 • T2 -1,409-10"7 • T3 + 2,342-10"11 • T4, T >01.
Зависимость теплоемкости от температуры без учета процессов в двухфазной зоне (включая окрестность точки плавления железа), основанная на аппроксимации экспериментальных данных и использовании 5 -образной функции в окрестности точки плавления, показана на рис. 10.
С
р
12
0 500 1000 1500 2000
Рис. 10. Зависимость теплоемкости от температуры
Рис. 11. Доля жидкой фазы в двухфазной зоне в зависимости от температуры
Расчеты показали, что в области фазового перехода (плавления) теплоемкость в приближении двухфазной зоны меняется слабо. Зависимость доли жидкой фазы от температуры представлена на рис. 11
Экспериментальные данные по теплопроводности сталей с содержанием углерода до 0,8 % (рис. 12) можно описать единой зависимостью [31] (Вт/см2)
X = 1,222 - 0,0016 • Т - 8,3617 • 10-7Т2 -1,1399 • 10-10Т3.
(24)
Величина коэффициента диффузии углерода, естественно зависит от температуры и концентрации. В расчетах использованы следующие зависимости
где
В* =
[Ах, Т <е1;
[А , е1 ^ Т < ТНд • Коэффициенты диффузии в твердых фазах зависят от температуры по закону Аррениуса
Рис. 12. Теплопроводность чистого железа
Аа = 0,02 • ехр(- 84219/ЯТ),
Ау = [0,04 + 0,08а] • ехр(-131356/ЯТ) см2/с. (19)
Коэффициент диффузии углерода в жидком железе А, входящий в (25), в общем случае не описывается аррениусовской зависимостью. В [36] наиболее надежной считают зависимость вида
А =(0,84 + 0,26 • а )-10 -4 см2/с,
справедливую для содержания углерода до 2,4 %.
Начальное распределение температуры аппроксимировано функцией вида
/( у) = ТтахехР(- а • х 2 )• ехр(- у • у).
Приведем лишь некоторые результаты исследования модели.
(20)
Очевидно, что температура наплавленной области постепенно уменьшается (рис. 13, нижняя кривая) вследствие прогрева подложки по механизму теплопроводности.
Распределение концентрации углерода имеет аналогичный вид, но глубина его проникновения незначительна (составляет несколько микрон). Это проиллюстрировано на рис. 14, где представлено распределение концентрации вдоль оси симметрии.
Рис. 13. Зависимость температуры поверхности от времени
Рис. 14. Распределение концентрации углерода в поверхностном слое ^=1-0,00025; 2-0,0015; 3-0,00525; 4-0,07525; 5-5,0 с
Вследствие уменьшения температуры к моменту времени 0,007 с диффузия существенно замедляется и в дальнейшем происходит весьма слабое изменение размеров диффузионной зоны за время остывания подложки.
Интересно отметить, что вследствие зависимости температур ликвидуса и солидуса от концентрации углерода кристаллизация раньше начинается с поверхности, чему способствует теплоотдача, хотя и в меньшей степени.
Увеличение температуры расплавленной области в начальный момент времени приводит к более значительному прогреву подложки (верхняя кривая на рис. 14), более широкой диффузионной зоне, следствием чего, очевидно, будет более плавное изменение свойств в окрестности границы раздела материалов в результате обработки.
Вообще говоря, для описания процессов поверхностной обработки существует целый ряд более простых моделей, включающих в себя много качественных эффектов. В общем случае диффузионные явления, реализующиеся в технологических процессах в неравновесных условиях, следует изучать с учетом разнообразных перекрестных явлений.
Механические напряжения в зоне обработки
Вопрос об оценке механических напряжений в зоне обработки является весьма актуальным. Известные в механике подходы к оценке остаточных напряжений, вообще говоря, не дают представления о том, с какими необратимыми процессами они связаны. На практике же важно знать, как изменяются напряжения в зоне обработки при изменении технологических параметров, состава обрабатываемых слоев и сопутствующих физико-химических превращений. В этом случае плодотворным оказывается путь, основанный на термодинамике необратимых процессов, вернее, ее обобщении на деформируемые среды с химическими превращениями и диффузией. Подробно обобщение модели механики с использованием идей термодинамики описано, например, в [37, 38]. Здесь ограничимся простым примером.
Механическое поведение плоского образца (с покрытием или без него), который подвергается обработке источником энергии, движущимся вдоль поверхности в направлении одной из осей, можно описать в рамках представления об обобщенном плоском напряженном состоянии [39, 40].
Тогда для пластины, изображенной на рис. 3, при условии, что ширина сканирования равна ширине образца, имеют место соотношения:
С гг = с гх = С гу = 0 , с хх = с хх (х,у ), с уу = с уу (х,у).
Учитывая, что скорость распространения тепла при теплопроводности и
скорость движения источника по поверхности много меньше скорости распространения механических возмущений, для оценки компонент тензора напряжений в первом приближении можем использовать идеи, описанные в [40] для задач термоупругости. Тогда для того, чтобы найти с хх и с уу, нам потребуется решить задачу о равновесии
при заданном распределении температуры Т(х, ^). Поскольку поле температуры при условии кх = Нх не зависит от координаты у, можно предположить, что от координаты у не зависят и компоненты тензоров напряжений и деформаций.
Тогда для решения задачи о равновесии можно воспользоваться условиями совместности и условиями равенства нулю суммарных сил и суммарных моментов по контуру пластины. В чистом виде этот путь приемлем, когда не достигается температура плавления.
Именно такой способ оценки напряжений в диффузионной зоне использован в [41,42].
Более строгая постановка задачи для условий, проиллюстрированных на рис. 3, может быть сформулирована в рамках теории термоупругости для пластин с учетом возникающих прогибов и моментов. Но это выходит за рамки настоящей работы.
Заключение
Заметим, что все описанные в статье модели исследованы численно. Методы численного решения имеют свои особенности, описать которые в рамках этой статьи не представляется возможным Основная проблема - в разномасштабности физических процессов, учитываемых в моделях. Так, коэффициенты диффузии и температуропроводности для твердых сред отличаются на несколько порядков. Это приводит к тому, что характерные размеры зоны прогрева и диффузионной зоны существенно различны, особенно если принять во внимание аррениусовскую зависимость коэффициента диффузии от температуры. На порядки могут отличаться и парциальные коэффициенты диффузии в многокомпонентных системах. В свою очередь, области влияния чисто температурных напряжений и напряжений, возникающих в диффузионной зоне, также различны. Все это приводит к необходимости разработки специальных численных алгоритмов. Вариант алгоритма для неизотермических задач диффузии предложен, например, в [43]. Другая трудность в численном исследовании связана с тем, что в случае воздействия концентрированных источников энергии области, занимаемые разными фазами и обладающие различными теплофизическими и реологическими свойствами, существенно различаются по размерам. Основные физико-химические изменения происходят в области, много меньшей размеров всего образца. Это опять приводит к поиску специальных методов решения, что представляет собой отдельную нетривиальную задачу.
Библиографический список
1. Рыкалин Н.Н. Расчеты тепловых процессов при сварке. - М.: Машгиз, 1951.—
296 с.
2. Теория сварочных процессов: учебное пособие / под ред. В.В. Фролова. — М.: Высшая школа, 1988. — 559 с.
3. Рыкалин Н.Н., Зуев И.В., Углов А.А. Основы электронно-лучевой обработки материалов. — М.: Машиностроение, 1978. — 239 с.
4. Шипко А.А., Поболь И.Л., Урбан И.Г. Упрочнение сталей и сплавов с использованием электронно-лучевого нагрева. — Минск: Наука и техника, 1995. — 280 с.
5. Модифицирование и легирование поверхности лазерными, ионными и электронными пучками / под ред. Дж. Поута. — М.: Машиностроение, 1987. — 424 с.
6. Демидов В.Н., Князева А.Г. Расчет термических режимов электроннолучевой обработки поверхностей деталей, имеющих форму фигур вращения //Proceeding 12th International conference on radiation physics and chemistry of inorganic materials, 23—27 September, 2003. — Томск: Изд-во ТПУ. — С.200—204.
7. Демидов В.Н., Князева А.Г. Численное исследование режимов термической
обработки поверхностей деталей, имеющих форму фигур вращения // Физическая мезомеханика. — 2004. — Т. 7, Спец. вып. — Ч. 1. — С.54—58. ( Материалы
международной конференции по физической мезомеханике, компьютерному конструированию и разработке новых материалов. 23—28 августа, Томск, 2004)
8. Численное моделирование теплофизических процессов при лазерно-лучевой сварке с образованием парового канала / А.Н. Черепанов, В.П. Шапеев, В.М. Фомин [и др] // ПМТФ. — 2006. — Т.47, № 5. — С. 88—96.
9. Лыков А.В. Теория теплопроводности. — М.: Высшая школа, 1967. — 599 с.
10. Егерев В.К. Диффузионная кинетика в неподвижных средах. — М.: Наука, 1970. — 239 с.
11. Карслоу Г., Егер Д. Теплопроводность твердых тел. — М.: Наука, 1964. —
487 с.
12. Солоненко О.П., Головин А.А. Нестационарный сопряженный теплообмен и фазовые превращения при высокоэнергетической обработке поверхности. Часть 2: Моделирование технологических процессов // Теплофизика и аэромеханика. — 2007. — № 4. — Т.14. — С. 623—638.
13. Головин А.А., Солоненко О.П. Нестационарный сопряженный теплообмен и фазовые превращения при высокоэнергетической обработке поверхности. Часть 1: Вычислительный метод и его реализация // Теплофизика и аэромеханика. — 2007. — № 3.
— С. 413—427.
14. Князева А.Г., Букрина Н.В., Поболь И.Л. Моделирование электроннолучевой обработки поверхностей с учетом усадки порошкового слоя // Радиационная
физика твердого тела: матер. XIV Международ. совещания, Севастополь, 5—10 июля 2004 г. — Севастополь, 2004. — C. 434—438.
15. Букрина Н.В., Князева А.Г. Трехмерная модель термической обработки поверхностей внешним источником энергии с учетом усадки порошкового слоя //Физико-химические процессы в неорганических материалах: матер. IX международ. конф., 22—25 сентября, Кемерово, 2004. — Т.2. — С. 376—380.
16. Букрина Н.В., Князева А.Г., Поболь И.Л. Влияние кинетики усадки
порошкового слоя на режимы формирования покрытия в процессе ЭЛН // Физическая мезомеханика: матер. международ. конф. по физической мезомеханике,
компьютерному конструированию и разработке новых материалов. 23—28 августа, Томск, 2004.—Т. 7, Спец. выпуск. — Ч. 2. — С. 193—196.
17. Электронно-лучевая наплавка износостойких композиционных покрытий на основе карбида титана / В.Е. Панин, В.Г. Дураков, Г.А. Прибытков, С.И. Белюк [и др.] // ФХОМ. — 1997. — №2. — С. 54—58.
18. Крюкова О.Н., Князева А. Г. Влияние динамики поступления частиц в расплав на фазовую структуру и свойства покрытия, формирующегося в процессе электронно-лучевой наплавки // Физическая мезомеханика — 2004.— Т.7, Спец. выпуск.
— Ч. 2. — С. 205—208.
19. Крюкова О.Н., Князева А. Г. Физико-химические явления в ванне расплава в процессе электронно-лучевой обработки поверхностей с модифицирующими частицами // Физическая мезомеханика— 2006. — Т. 9, Спец. выпуск. — С. 65—68.
20. Крюкова О.Н., Князева А. Г. Критические явления при растворении частиц в ванне расплава в процессе электронно-лучевой наплавки покрытий // ПМТФ. — 2007. — Т. 48. — № 1. — С. 131—142.
21. Вигдорчик Е. М., Шейнин А. Б. Математическое моделирование
непрерывных процессов растворения. — Л.: Химия, 1971. — 248 с.
22. Романков П.Г., Рашковская Н.Б., Фролов В.Ф. Массообменные процессы химической технологии. — Л.: Химия,1975. — 336 с.
23. Химическая гидродинамика: справоч. пособие / А.М. Кутепов, А.Д. Полянин, З.Д. Запрянов, А.В. Вязьмин, Д. А. Казенин. — М.: Бюро Квантум, 1996. — 336 с.
24. Крюкова О.Н. Численное исследование модели электронно-лучевой наплавки покрытий с модифицирующими частицами с учетом физико-химических превращений // Известия ТПУ. — 2006. — № 6. — С. 120—125.
25. Поболь И.Л. Научные и технологические основы обработки конструкционных и инструментальных материалов и получения изделий с использованием электронно-лучевого воздействия: дис. д-ра техн. наук. — Минск: ФТИ НАН, 2007. — 340 с.
26. Самораспространяющийся высокотемпературный синтез: теория и практика / отв. ред. А.Е. Сычев. — Черноголовка: Территория, 2001. — 432 с.
27. Knyazeva A.G., Pobol I.L., Gordienko A.I. Coating formation in SHS-regime during thermal treatment of material by moving energy source // 7-th International Conference on modification of Materials with Particle Beams and Plasma Flows, Tomsk, Russia, 25—30 July 2004. — P. 178—183.
28. Knyazeva A.G., Pobol I.L. Composition Coating Formation on the Base of Intermetallide Phase Using Electron-Beam Treatment // Proceedings. 8th International conference on modification of materials with particle beams and plasma flows. Tomsk. Russia. 10—15 September. 2006. — Известия вузов, сер. Физика. — 2006. — Т.49. — № 8 (Приложение) Изд-во Томск: ТГУ. — С. 396—399.
29. Сорокова С.Н., Князева А.Г. Моделирование формирования фазовой структуры покрытия в процессе электронно-лучевой обработки с использованием синтеза в твердой фазе // Теоретические основы химической технологий. — 2008. — Т.42, №4. — С. 457—465.
30. Sorokova S.N., Knyazeva A.G. Modeling of coating modification using condensed phase synthesis // Book of reports of the International Workshop ”Nonequilibrium Processes in Combustion and Plasma Based Technologies” (Minsk, August 23—28, 2008). — P .65—69.
31. Князева А.Г., Дураков В.Г., Прибытков Г. А. Формирование диффузионной зоны в подложке в процессе электронно-лучевой наплавки покрытий // in Proceeding of 6th Int. Conference on Modification of Materials with Particle Beams and Plasma Flows. Proceeding. 23—28 September, Russia, Tomsk, 2002. — С. 613—616.
32. Борисов В.Т. Теория двухфазной зоны металлического слитка. - М.: Машиностроение, 1987. - 223 с.
33. Диаграммы состояния систем на основе железа: справочник / под ред. О. А. Банных, М.Е. Дрица. - М., 1986. - 440 с.
34. Физические свойства сталей и сплавов, применяемых в энергетике: справочник / под ред. Б.Е. Неймарка, - М.-Л., 1967. - 240 с.
35. Зиновьев В. А. Теплофизические свойства металлов при высоких температурах. - М.: Машиностроение, 1989. - 384 с.
36. Ершов Г.С., Майборода В.П. Диффузия в металлических расплавах. - Киев: Наукова думка, 1990. - 224 с.
37. Князева А. Г. О моделировании необратимых процессов в материалах с
большим числом внутренних поверхностей // Физическая мезомеханика. - 2003. - Т.6.
- № 5. - С. 11-27.
38. Князева А.Г. Диффузия и реология в локально-равновесной термодинамике //Математическое моделирование систем и процессов / под ред. П.В. Трусова. Пермь: Изд-во Перм. гос. техн. ун-та, 2005. - С. 45-60.
39. Безухов Н.И. Основы теории упругости, пластичности и ползучести. - М.: Высшая школа, 1951. - 537 с.
40. Боли Б., Уэйнер Дж. Теория температурных напряжений. - М.: Мир, 1964. -510 с.
41. Bukrina N.V., Knyazeva A.G. Estimation of mechanical stress in diffusion zone during the cooling of built-up coating // Proceeding 9th International Conference on Modification of Materials with Particle Beams and Plasma Flows. Tomsk. Russia, 21-26 September 2008.- Р. 451-455.
42. Букрина Н. В., Князева А. Г.Об оценке механических напряжений в материале с покрытием в условиях ионной бомбардировки // Физическая мезомеханика.- 2008.- Т.11.- №5.- С. 95-102.
43. Букрина Н.В., Князева А.Г. Алгоритм численного решения задач неизотермической диффузии, встречающихся в процессах поверхностной обработки // Физическая мезомеханика. - Т. 9. - № 2. - 2006. - С. 55-62.
Получено 20.06.2009
