СВ-СИНТЕЗ
ДИНАМИКА ТЕПЛОВЫХ ПРОЦЕССОВ ПРИ ФОРМИРОВАНИИ ЗАЩИТНЫХ ПОКРЫТИЙ В ТЕХНОЛОГИЯХ ДЕТОНАЦИОННО -
ГАЗОВОГО НАПЫЛЕНИЯ
В.Ю.Филимонов, К.Б.Кошелев, В.И.Яковлев
Применение технологий детонационно -газового напыления (ДГН) позволяет эффективно решить проблему получения защитных покрытий из композиционных материалов деталей и узлов машин ввиду того, что процесс ДГН характеризуется высокой скоростью дисперсного потока, до 1000 м/с, и температурами достигающими 3000 К [1]. Однако указанный метод является импульсным. В процессе нанесения покрытий происходит кристаллизация капель напыляемого компонента за
время порядка 10-3 ^ 10^с, толщина наносимого за один цикл слоя, составляет несколько микрометров, в объеме формирующегося слоя возможны фазовые превращения. Все это создает большие трудности экспериментального изучения динамики процессов образования покрытия в процессе ДГН. Следовательно, создание теоретической прогнозной модели, описывающей разогрев поверхности контакта основы для напыления, с напыляемым материалом, является актуальной задачей, поскольку температура контакта является главным динамическим параметром, определяющим адгезионные характеристики покрытия.
Развитые на сегодняшний день теоретические модели, относятся в основном к описанию одиночных актов взаимодействия расплавленных частиц с поверхностью основы [1,2], однако эти процессы не допускают прямой экспериментальной проверки, по этой причине, при моделировании, необходим иной подход к описанию указанных процессов, с точки зрения измеряемых макроскопических величин.
В настоящей работе проведено теоретическое исследование макрокинетики процессов разогрева поверхности основы для напыления с учетом возможных фазовых превращений.
При построении теоретической модели использовались следующие упрощающие предположения:
1. Основа для напыления является плоской, конечной толщины. Площадь торцевой поверхности значительно меньше площади поверхности напыления.
2. Поток напыляемых частиц ламинарный, распределение скорости в потоке отсутствует как по абсолютной величине, так и по направлению.
3. Скорость потока постоянна по времени и направлена нормально к поверхности основы. Рост напыляемого слоя происходит с постоянной скоростью.
4. Рост конвективного потока и его спад происходят мгновенно, через конечный интервал времени, циклично.
С учетом возможных фазовых превращений, распределение температуры в области напыляемой основы следует разбить на три характерных зоны (рис.1). Зона 1 ( )< г < к()) - напыляемый слой в жидкой фазе. Здесь )- граница зоны плавление -кристаллизация, к(1) - внешняя граница поверхности напыления, скорость движения которой за интервал времени напыления постоянна, и определяется через скорость конвективного потока с использованием уравнения неразрывности (р1У1 = р"У*). Зона 2
(г < г < I
('))
напыляемый слой в конден-
сированной фазе). Здесь г0 - толщина основы. В случае I = г0 - напыляемый слой полностью в жидкой фазе. В случае ) = к({) напыляемый слой полностью в конденсированной фазе. На границе г0 - накладываются
граничные условия четвертого рода для уравнения теплопроводности. На границе ) накладываются граничные условия для поверхности плавления в режиме Стефана.
На границе к(1:) накладывается граничное условие второго рода для конвективного потока. На границе г = 0 накладывается граничное условие третьего рода (закон Ньютона - Рихмана). Система уравнений теплопроводности для областей в обезразмеренной форме запишется в виде (1):
Го г© ВД
Рис.1. Схематическое распределение температуры в областях: 1 - жидкая фаза, 2 - напыляемый слой в конденсированной фазе, 3 - основа
5©
а©2
дт
д2© x2 д2 ©
д1=ь д© *<х < у(Т
дт дx
_ Ь2^ ; 1 < x <*(т)
дx
(1)
д©3 д 2©3
0 < x < 1
дт дx2 с набором граничных условий:
x = АЛ ^ = в (1 -© у ) ^ = uн,
дx 1 dт
©1 *_© 2 I * =©,
x = 1> © 2 =©31,
д©2\ , д©3 21 = d1 3
дx
дx 1
■ = *{т), (2) (3)
* dт
0, при© и, при©
ю
12
д©„
д*
*-0 ы2
1
у ,
д©
(4)
д*
*+0
д© 3
x _ 0, 3
дx
0 = Вт ^0 (5)
где ©
© =
T - Т
ал
т' - т0
0. т
и^ x __* _
к
а
а,,
у —, Ь _-*-, Ь2 _ и
-УН
а
а
а,
4 d2 _ 4 ,
42 42
При начальных условиях:
1, при x > 1
[0, при 0 < x < 1 с'р'у *г0
т _ 0, ©_
*
с
в
в12 _
в2 _
4
С рУ *г0
42
, при *< у(т) , при * _ у(т)
4 (т' - Т0)
а3р1Ь
4 (- Т0)
, при * <
*< у(т)
а3р2 Ь
, при*_ у(т).
Здесь Г0, Т* - температура окружающей среды и дисперсного потока соответственно,
а
1-3
, Т1-3 ,4-
температуропроводность,
температура и теплопроводность в соответствующих областях (рис.1), г - координата, t - время, г0 - толщина основы, У - скорость роста внешней границы слоя к,
* *
ср, р - удельная теплоемкость и плотность
газодисперсного потока, р1, р2 - плотности жидкой и твердой фаз соответственно, Ь - удельная теплота плавления напыляемого порошка, Н - обобщенная функция, вид которой представлен на рис.2:
н
Жг
1-1
Рис.2. Н - функция, используемая в модели: здесь 11 - время взаимодействия струи с основой, 12 - промежуток между циклами напыления
Двойной индекс обозначает то, что внешняя граница зоны 1 может находиться либо в конденсированной, либо в жидкой фазе.
Модель (1 - 6), аналогична модели, представленной в работе [2]. Однако автора-
т
3
г
г
г
0
0
0
г
0
г
0
В.Ю.ФИЛИМОНОВ, К.Б.КОШЕЛЕВ, В.И.ЯКОВЛЕВ
ми [2] изучался процесс непрерывного (а не циклического) напыления для полубесконечной основы, без учета кинематических и теп-лофизических параметров газодисперсной струи.
Параметры Ъ{, di выбирались характерными для напыляемых на сталь металлов: Ъх = 1,5, Ъ2 = 2, d1 = 1,2, d2 = 1,3.
Предполагалось р1 « р2 (расчет показывает, что отличие в плотностях на 20 - 30% не оказывает существенного влияния на значимые
зависимости задачи). Параметры т1 и т2 фиксировались и выбирались характерными для процессов детонационного напыления:
т1 = 10-3, т2 = 10-4 . Параметр оп оказывает влияние лишь на скорость плавления в напыляемом слое, но не на закономерности разогрева поверхности за времена существенно большие длительности одного цикла. Варьируемыми являлись параметры в12, 01, Bi (01 - безразмерная температура
плавления напыляемого материала). Расчет проводился методами численного моделирования.
На рис.3 представлена характерная зависимость температуры напыляемой поверхности и координаты фронта плавления от времени при наличии фазового превращения в объеме слоя. Из графика следует, что термограмма носит осциллирующий характер. На амплитуду осцилляций оказывает влияние скважность импульсов напыления (время т2). В точке фазового превращения зависимость претерпевает излом, что связано с ощутимым
различием параметров Рх, в2. Если значения указанных параметров близки, излом менее рельефный.
Рис.3. Зависимость обезразмеренной температуры 0 и координаты фронта плавления 1(т) от времени для значений: 01 = 0,25; Д = 0,16; в2 = 0,05; Bi = 0,05
На протяжении конечного промежутка времени, координата фронта плавления осциллирует с периодом, равным времени скважности. Длительность этого промежутка времени определяется временем нагрева поверхности до температуры плавления 01,
другими словами колебания фронта имеют место тогда, когда внешняя поверхность напыляемого слоя имеет температуру выше температуры плавления, а внутренняя ниже. После прогрева колебания исчезают, и далее напыляемый слой находится в жидкой фазе. Однако колебательный режим может иметь устойчивый характер при определенных значениях теплофизических параметров детонационной струи и основы. На рис.4 представлена диаграмма плавления на плоскости параметров в-Bi. Области ниже затемненных зон соответствуют твердофазному состоянию напыляемого слоя, область выше зон - жидкофазному. Промежуточные зоны -устойчивый колебательный режим плавления.
0 О 1 0.2 0 3 В1
Рис.4. Диаграмма плавления на плоскости параметров в- Bi: 1 - 01 = 0,25, 2 - 01 = 0,3, 3 - 0, = 0,3
В случае если система находится в области параметров плавления (выше затемненных зон), закономерно встает вопрос о временах достижения точки плавления. На рис.5 представлены зависимости времени плавления от критерия Bi для различных значений ©1 и в. На рисунке вертикальные линии иллюстрируют границы зон плавления (значения Bi, в при которых плавления не происходит, математически время плавления становится бесконечным). Монотонный рост зависимости т(В^) объясняется ростом времени достижения заданной температуры ликвидус с ростом теплоотвода. С другой стороны рост параметра в (кривые 1 - 4), при заданном Bi, уменьшает время плавления, т.к. при этом возрастает плотность теплового потока в момент соударения струи с основой, что ускоряет прогрев до температуры плавления и дает возможность достигнуть границы зоны плавления при большем теплоотво-де.
с)
Рис.5. Зависимость обезразмеренного времени плавления от критерия Ы: а) © _ 0,25, Ь)
©1 _ 0,3, с) © 1 _ 0,35; 1 - в _ 0,05, 2 - в _ 0,1, 3 -в_ 0,15,4 - в_ 0,2
Из приведенного анализа следует, что изменяя значения параметров дисперсной струи, или теплофизические характеристики основы можно добиться плавления напыляемого материала непосредственно на основе напыления, либо по заданным значениям указанных параметров прогнозировать возможность ликвации, при выполнении последней, прогнозировать время напыления до начала плавления, соответствующее количество циклов и толщину слоя нанесенного покрытия, при этом основными значимыми параметрами задачи для данного напыляемого
материала являются: в_ с*рр*У*г0/А (отношение конвективного теплопереноса на поверхность к кондуктивному теплопереносу в пластине) и Bi _ог0/4 (отношение внешнего теплового сопротивления пластины основы, к внутреннему).
ЛИТЕРАТУРА
1.Порошковая металлургия и напыленные покрытия. Под.ред. Б.С.Митина. М.: Металлургия, 1987.
2. Гавриленко Т.П., Николаев Ю.А., Прохоров Е.С., Ульяницкий В.Ю. О механизмах образования покрытий при газотермическом напылении.// Физ. гор. и взрыва, 1990, т.26, №2. - С.110 - 123.