Математическое моделирование взаимодействия лазерного излучения с плазменным потоком при лазерно-плазменном напылении

Математическое моделирование взаимодействия лазерного излучения с плазменным потоком при лазерно-плазменном напылении

Чащин Е.А., Шилов И.В., кандидаты техн. наук, Балашова С.А., инж.

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

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

Mathematical modeling of interacting laser light with plasma stream at laser-plasma spraying

Chashchin E.A., Shilov I.V., Cand.Tech.Sci., Balashova S.A., a management engineer

In work, by a synthesis of model of a laser-plasma spraying, on the basis of the analysis of absorption of laser light in plasma stream, optimal conditions of interacting of two concentrated energy flows are defined at plotting ceramic coatings AhO3.

Keywords: modulated laser radiations, plasma stream, spraying, an adhesion strength, mathematical model, thermal processes.

Исследование процесса лазерно-

плазменного напыления поверхностных керамических покрытий на основе А12Оз, выполненное нами ранее [1], показало, что существует возможность повышения в 1,5-2 раза адгезионной прочности напыляемых покрытий при введении модулированного лазерного излучения (ЛИ) в плазменный поток (ПП). Совместное с нагревом ПП воздействие быстропротекаю-щих процессов нагрева импульсами ЛИ нано-секундного диапазона ограничивает возможность использования известных аналитических зависимостей [2, 3]. Это делает актуальной разработку математической модели взаимодействия лазерного излучения с плазменным потоком, содержащим мелкодисперсную фазу, позволяющую определить оптимальные условия взаимодействия двух концентрированных потоков энергии. Решение этой задачи выполним путем синтеза модели лазерно-

плазменного напыления и оптимизации дополнительного теплового воздействия ЛИ на частицы плазмонапыляемого порошка А12О3 на основе анализа поглощения ЛИ в ПП.

Рассмотрим поглощение ЛИ в ПП. Принимаем, что для обеспечения равномерного прогрева ПП при дополнительном воздействии поглощение ЛИ должно иметь объемный характер. Также принимаем, что эффекты, связанные с резонансным поглощением ПП энергии ЛИ, пренебрежимо малы. Тогда степень поглощения ПП дополнительной энергии ЛИ удобно определить по формуле Унзольда -Крамерса [4]:

1бп2 в6г2кГп ( е*-ЬуУ, (

-ехр|----- II 1 - ехр| —

3л/3 с • h4v3

kT

kT

где с - скорость света; V - частота излучения; Z - заряд иона; е* - потенциал ионизации атома; Т - температура ПП; п - концентрация

атомов в ПП; к - постоянная Больцмана; Л -постоянная Планка.

Большинство известных технологических установок обеспечивают генерацию низкотемпературной плазмы с температурой до 20 кК. Результаты расчета по формуле Унзольда -Крамерса приведены на рис. 1. Анализ полученных зависимостей (рис. 1) показывает, что при температуре менее 10 кК ПП практически «прозрачен» для ЛИ, поэтому ЛИ целесообразно вводить в ПП в области максимальной температуры потока, а именно, на срезе сопла плазмотрона, где температура ПП может достигать значений в интервале 10-20 кК. При воздействии ЛИ на ПП в этой области толщина

§и

*=1

Кис

скин-слоя, в котором происходит ос-

лабление излучения в е раз, нелинейно уменьшается с ростом температуры ПП. Однако даже при температуре ПП 20 кК толщина скин-слоя достигает 4 см, что значительно превышает сечение ПП на срезе сопла плазмотрона и, как следствие, делает малоэффективным поглощение ПП энергии ЛИ.

Т, кК

Рис. 1. Зависимость коэффициента объемного поглощения аис и толщины скин-слоя Зис от температуры ПП

Известно, что модулированное ЛИ с малой длительностью импульсов и высокой

плотностью мощности приводит к возникновению оптического пробоя, что позволяет перевести газовую среду в плазменное состояние. При этом в области оптического пробоя поглощается до 95 % энергии падающего лазерного излучения [5]. Принимаем, что совмещение ПП в зоне его выхода из плазмотрона с областью оптического пробоя обеспечивает повышение температуры плазмы, а следовательно, и напыляемого материала за счет практически полного поглощения энергии ЛИ.

Известно [6], что обязательным условием возникновения пробоя является нарастание концентрации возбужденных электронов при превышении частоты ионизации V,- над частотой диффузионных потерь у0:

1

Е = ю/(а • п -а • е)|^(тее*)/3^2 .

Тогда интенсивность ЛИ, обеспечивающего возникновение пробоя, удобно определить как

I = сее0Е2,

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

При формировании ЛИ оптического пробоя распространение фронта ионизации может проходить в разных режимах [6]. Наиболее предпочтительным является режим све-то-детонационной волны (скорость 5-20 км/с, интенсивность ЛИ 20-1000 МВт/см2), который характеризуется хорошей устойчивостью и средними значениями пороговой интенсивности пробоя.

В настоящее время достигнуты результаты генерации лазеров модулированного ЛИ на базе технологических ИАГ:Ш с длительностью импульса 60-350 нс и плотностью мощности до 2-5 ГВт/см2. Учитывая, что интенсивность ЛИ может быть уменьшена изменением условий фокусировки, а также, что геометрические размеры ПП на срезе сопла плазмотрона, как правило, имеют значения в интервале 1-10 мм, принимаем, что при формировании оптического пробоя в ПП фронт ионизации распространяется в режиме свето-

детонационной волны (СДВ) со скоростью

VСДВ =(у2 -1)) I^,

1 ’ 1Р0)

где I - интенсивность излучения; р0 - плотность газовой среды; у - показатель адиабаты газа.

Рассмотрим зависимость интенсивности ЛИ от температуры ПП, при которой СДВ распространяется на расстояние, перекрывающее ПП. Результаты расчета, полученные для различных условий фокусировки ЛИ приведены на рис. 2.

Рис. 2. Зависимости пороговой интенсивности пробоя I: (1-3) в азоте от начальной температуры Т: 1 - диаметр лазерного луча 1 мм; 2 - 5 мм; 3 - 10 мм. Зависимость расстояния I (4-6), пройденного фронтом ионизации в режиме СДВ, от Т: 4 - при длительности импульса 300 нс; 5 -200 нс; 6 - 100 нс

Анализ полученных зависимостей (рис. 2) показывает, что с увеличением температуры ПП интенсивность ЛИ, при которой возникает пробой, нелинейно снижается. Также снижается расстояние, пройденное фронтом ионизации, при уменьшении интенсивности ЛИ.

Следует отметить, что при воздействии импульсов ЛИ с длительностью до 300 нс время формирования фронта ионизации в зоне оптического пробоя составляет 1-5 нс, что много меньше длительности импульсов. Поэтому, для уменьшения затрат машинного времени, расчеты выполним с итерационно изменяемым шагом по времени. Изменение температуры ПП удобно найти решением одномерной задачи нагрева среды с объемным поглощением:

ДТ = аис(Тр1 ) • 1 • ^

рр1 • Ср1

где t - время воздействия излучения; Рр1 -плотность ПП; Ср - теплоемкость плазменного потока, принимаемая при температуре ПП более 10000 К равной 13 кДж/(кг-К).

В качестве критерия изменения шага по времени принимаем изменение прироста температуры в 2 раза.

Таким образом, проведенный анализ определил: объемный и равномерный прогрев ПП с начальной температурой 6-20 кК обеспечивается созданием ионизационного фронта в режиме СДВ на удаленной от ИАГ :Ш лазера стороне ПП путем фокусировки луча ЛИ со следующими параметрами излучения: длительность импульсов 100-300 нс; энергия импульсов 0,1-0,5 Дж; интенсивность 107-109 Вт/см2.

При выполнении операции синтеза устройства лазерно-плазменного напыления [7] полагали, что наблюдаемый эффект повышения адгезионной прочности керамических покрытий вызван повышением температуры в контактной области за счет интенсификации

теплового воздействия на частицы А12О3. В этом приближении представляется достаточным с помощью соответствующего выбора параметров ЛИ добиться повышения температуры частиц до температуры 3000-3500 К [8]. Таким образом, задачей этапа является математическое описание влияния ЛИ на температуру частиц Л120з при совместном воздействии ЛИ и ПП. Расчет выполняли в сферическом приближении частиц А12О3, подвод теплоты полагали равномерным по всей поверхности (рис. 3).

С учетом сделанных приближений температуру частиц напыляемого материала удобно определить из решения единого нестационарного уравнения теплопроводности: дТ/дТ + (/ дгадТ) = а(Т )ДТ, где а(Т) = Хт(Т)/ с(Т)р - температуропроводность; Хт(Т) - теплопроводность; р - плотность материала; V- вектор скорости фронта испарения; с(Т) - удельная теплоемкость.

Рис. 3. Модель состояния частицы А12О3

Поскольку частица А12О3 тугоплавкая и имеет малый радиус 10-100 мкм, то справедливо предполагать, что сначала она вся полностью расплавляется и только потом начинается процесс испарения. С учетом этого, для уменьшения затрат машинного времени область нагрева моделировали без разделения на твердое тело и расплав, а изменение свойств материала учитывали в эффективной сглаженной функции удельной теплоемкости с(Т), которая, кроме сглаживания перепада в значении при фазовом переходе «твердое тело - расплав», включает в себя удельную энергию плавления. С учетом сделанных упрощений, граничные условия на поверхности поглощения-испарения имеют вид д

-Хт —Т = Ят а), дп

Ят ^) = (1 - R(Ts)) (г, t) - V (Т(г)) + яv, где ЯТ - интенсивность тепловыделения на облучаемой поверхности; д/дп - дифференци-

рование по нормали к поверхности канала; 1±(г,0 - распределение интенсивности падающего излучения в плоскости, нормальной оси луча; 1/ = Vs р 1/ - потери плотности мощности на испарение; - лучистый обмен меж-

ду поверхностью частицы и ПП; Я - отражательная способность поверхности, учитывающая ее температуру ^. Скорость фронта испарения ^ задается также температурой поверхности ^.

Перепишем уравнения теплопроводности в сферической системе координат:

1 д

АТ = АТ -^—\г2— I,

г2 дг

дТ

дг

(V, дгабТ) = V

дТ

дг '

Поскольку вблизи оси, где 1/г , воз-

никает неопределенность в операторе АгТ, чтобы этого избежать, применим удобную подстановку по зависимости

Р (г) = гт (г).

С учетом сделанных допущений, уравнение теплопроводности удобно представить в следующем виде:

дТ=1дР дР = дТ Р дТ=1 (дР - Р

ді г ді ’ дг дг г ’ дг г [ дг г

Умножим уравнение теплопроводности на г и с учетом подстановки получим его аналог для Р в виде

дР гО

—=а(Т)А„г (р)+—,

где А„ (Р) =

дР

дг

д2Р

дг

2

Здесь деление на г остается, однако на оси сферы отпадает необходимость вычислять Р, так как Р(г = 0) = 0 , и при любом малом шаге по координате г неопределенности около оси не возникнет.

Граничные условия для расчетной области в виде центрально-симметричной сферы размера гсф запишутся в виде Р (г = 0) = 0,

д ч —xt ~д)£Р (г = гсф) = гсф Ят '

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

Для каждого внутреннего узла расчетной сетки по г переход от момента времени (S-1) к моменту S описывается уравнением

г є - г є-1 1 І 'і

А

гЛ

где

гі+і гі-і

V'г' = VrNг— СК0Р0СТЬ сдвига і-го узла сетки при

движении внешней границы узла сетки со скоростью фронта испарения V,-.

При поверхностном поглощении ЛИ член уравнения, учитывающий объемное тепловыделение О\/, можно опустить для всех узлов сетки, кроме крайнего узла (с индексом і = М), соответствующего облучаемой поверхности сферической частицы. Тогда объемное тепловыделение для поверхностного узла удобно получить из интенсивности теплового потока дт:

д

о.. =Л ат\ ^ ІТ_

V дгЧт'г=гсф 2Аг

'сф

Таким образом, при перенесении теплового потока дт из граничного условия непосредственно в дифференциальное уравнение (в ОУ для поверхностного узла) для облучаемой поверхности граничное условие упростится до выражения

зт

дг

= 0,

2'сф

которое для подстановочной переменной Г запишется в виде

зр

дг

Г№+1 = Г№-1

2Аг

№

г=Гсф 'сф 'сф

С учетом сделанных упрощений, уравнения, соответствующие шагу по времени S вдоль одной координаты, для каждого из внутренних узлов расчетного стержня приводятся к линейным уравнениям относительно 3-х узловых значений параметра Р и могут быть записаны в виде

Лр*-1 - с I Р^ + В,.р/?1 = Ъ,,

где А =Аат) Гі

В =

Аг

Аіа

Аг 2

Аг Vг і

1 -

2 а

Аг Vгi

ап

2Аіа

1 -

Аг2 Vг¡

г-1

2г, а(Т) у

Данное уравнение фактически уже приведено к каноническому линейному виду для метода прогонки относительно РД , Р¡s и Р+ .

Для крайнего узла с индексом / = Мг исключим соответствующее мнимому узлу зна-

чение Г,

используя граничное условие крайнего узла, и запишем выражение

А РМг-1 - Сі РМг = Ь і,

где

А = 2^; Сі = 1+2^Г1 О;

Аг2 Аг2 ^ Мг)

ь = г5-1 А£ гсфЯт

і=~ Мг - 2Аг с(т) .

Таким образом, для «теплового стержня», содержащего узлы расчетной сетки, задача нахождения неизвестных узловых значений сводится к системе линейных уравнений с 3-диагональной матрицей, а значит, для ее решения применим экономичный метод прогонки.

Для внутренних узлов расчетного «теплового стержня» вдоль координаты г рекуррентные формулы прогоночных коэффициентов X и п имеют вид

Хі =

В

п = Апі-1- Ь Сі- АхУ 1 сі - А-Хі-1

и связаны с узловыми значениями основным рекуррентным выражением

Р? = х,р,+1 + П/.

Для 0-го узла расчетного «теплового стержня» прогоночные коэффициенты определяются с использованием рекуррентного соотношения и граничного условия 1-го рода и записываются как

Р0S =%0Р^ +П0 = 0 ° П0 = 0.

Тогда распределение Р(г ) находится из обратной прогонки, начиная с конечного узла, по формуле

Р№ = ^ + Х№ П№-1 = П№, так как х№ = 0,

1 - 1ыг Ум-1

и далее по рекуррентному соотношению.

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

Характерные результаты расчета приведены на рис. 5. Анализ полученных зависимостей показывает, что в результате воздействия импульсного ЛИ температура поверхности сферической частицы существенно превышает температуру плавления напыляемого материала. В течение импульса ЛИ температура повышается до максимального значения и 6000 К, а после окончания действия импульса ЛИ температура поверхности интенсивно снижается до 3000-3500 К и стабилизируется, поскольку охлаждение поверхности частицы через испарение вещества прекращается. Кроме того, воздействие ЛИ сопровождается уменьшением начальных размеров частицы, кото-

рые после окончания импульса остаются неизменными.

При воздействии на частицу керамики ЛИ с длительностью импульса от 1-50 нс и энергией импульса от 0,1 до 0,25 Дж происходит нагрев частицы до 7000-9000 К, и при этом она испаряется на 60-80 %. А увеличение длительности импульсов свыше 400 нс приводит к полному испарению напыляемых частиц.

Рис. 4. Алгоритм определения температуры частицы А1203

Таким образом, для лазерноплазменного напыления А1203 целесообразным представляется воздействие импульсами модулированного ЛИ с длительностью в диапазоне от 100 до 300 нс и энергией импульса от 0,1 до 0,4 Дж во всем диапазоне температур ПП. При этом частица испаряется на 30-50%,

а температура стабилизируется на значениях 3000-3500 К, и частица напыляемого материала попадает на поверхность напыляемой детали полностью расплавленной.

Выполним оптимизацию дополнительного теплового воздействия ЛИ на частицы мелкодисперсного порошка А12О3, напыляемого плазменным методом. В качестве критерия оптимизации будем считать достижение максимума адгезионной прочности сцепления напыляемого покрытия при максимальной энергетической эффективности дополнительного воздействия ЛИ.

а)

0 1 ш’7 2 ю’’3 ю-’4 10 ’ 5 ю’’6 ю’? I, с

б)

О 1-Ю7 2-Ю1 3 10 7 4-Ю 7 5-10 7 6 10 1 1с

в)

Рис. 5. Зависимости изменения от времени 1: а - температуры поверхности частицы Т; б - скорости испарения частицы У„с„; в - радиуса частицы Р

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

N(t) / N0 = 1 - exp[-(v t)/exp(Ea/k • Tk)],

где v - частота собственных колебаний атомов, находящихся в контакте; t - время затвердевания частицы порошка; Еа - энергия активации поверхности подложки.

При условии идеального теплового контакта подложки с расплавленной частицей, что возможно при физико-химическом взаимодействии с привариванием покрытия к подложке, температура в зоне контакта определяется из выражения

Тк =[Ке(Тч - Т0)/(КЕ+ ФК))] + Т0, где Тч и Т0 - температура частицы и подложки в момент осаждения покрытия; Ф(ак) - функция интеграла вероятности; Ке = (^/А, 2) л/ a2/ a1

- критерий тепловой активности частицы по отношению к подложке; Х1 и Х2 - коэффициенты теплопроводности порошка и подложки; ai и а2 - коэффициенты температуропроводности порошка и подложки .

Результаты расчета приведены на рис. 6.

Рис. 6. Зависимость относительной прочности сцепления частиц с подложкой ОД/Ио от температуры в зоне контакта Тк, К: 1 - основа Ре; 2 - основа сталь 12Х18Н9Т

Анализ полученных зависимостей (рис. 6) показывает, что относительная прочность сцепления покрытия с основой при достижении температуры в зоне контакта 1500 К для железа и 1800 К для стали 12Х18Н9Т достигает максимально возможного значения. Температура частиц, которая вызывает нагрев в зоне контакта до таких величин, составляет и 3000 К, а при дальнейшем увеличении температуры частиц А12О3 улучшения прочности покрытия не наблюдается. При изменении температуры ПП энергетические параметры изменяются. В результате последовательности проведенных расчетов была получена область рациональных режимов для лазерноплазменного напыления (рис. 7).

Рис. 7. Рабочие режимы процесса лазерно-плазменного напыления: 1 - длительность импульса 100 нс, 2 - 150 нс, 3 - 200нс, 4 - 250 нс, 5 - 300 нс

Заключение

Предложенная математическая модель описывает тепловые процессы при взаимодействии ЛИ с ПП, содержащим мелкодисперсную фазу. Оптимальные условия взаимодействия двух концентрированных потоков энергии лежат в диапазоне: температура ПП - 60-20 кК; длительность импульсов модулированного ЛИ -100-300 нс; энергия импульса ЛИ - 0,1-0,4 Дж. Результаты моделирования качественно соответствуют известным теоретическим и экспериментальным данным.

Список литературы

1. Повышение эксплуатационных характеристик керамических покрытий / С.А. Балашова, Е.А. Чащин, И.В. Шилов, А.А. Митрофанов // Вестник СГТУ. - 2009. - № 2 (38). - С. 72-79.

2. Кудинов В.В. Плазменные покрытия. - М.: Наука, 1977.

3. Пузряков А.Ф. Теоретические основы технологии плазменного напыления. - М.: МГТУ им. Н.Э. Баумана, 2003.

4. Булгаков А.В., Булгакова М.Н. Тепловая модель импульсной лазерной абляции в условиях образования и нагрева плазмы, поглощающей излучение // Квантовая электроника. - 1999. - Т. 27. - № 2. - С. 154-158.

5. Балашова С.А., Чащин Е.А., Шилов И.В. Лазерно-плазменное напыление керамических покрытий: Междунар. молодежная науч. конф. (XXXIV Гагаринские чтения). Ч. 3. - М., 2008. - С. 126-127.

6. Воздействие лазерного излучения на материалы / Р.В. Арутюнян, В.Ю. Баранов, Л.А. Большов и др. -М.: Наука, 1989.

7. Полезная модель 75391 РФ. Плазмотрон для

лазерно-плазменного нанесения покрытий; заявл.

01.04.08. Балашовой С.А. и др.

8. Свидетельство об отраслевой регистрации разработки. Расчет температуры частицы при воздействии плазменного потока. № 12240 от 29.01.2009г., авторов Балашова С.А., Чащин Е.А., Воротнев П.В.

9. Погребняк А.Д., Тюрин Ю.Н. Модификация свойств материалов и осаждение покрытий с помощью плазменных струй // Успехи физических наук. - 2005. -Т.175. - № 5. - С. 515-545.

Чащин Евгений Анатольевич,

ГОУ ВПО «Ковровская государственная технологическая академия имени В.А. Дегтярёва», кандидат технических наук, доцент, зав. кафедрой электротехники, телефон (49232) 3-13-47, e-mail: kanircha@list.ru

Шилов Игорь Вячеславович,

ГОУ ВПО «Ковровская государственная технологическая академия имени В.А. Дегтярёва», кандидат технических наук, доцент кафедры лазерной физики и технологии, телефон (49232) 3-13-47, e-mail: kanircha@list.ru

Балашова Светлана Александровна,

ФГУП «ВНИИ "Сигнал"», ведущий инженер,

ГОУ ВПО «Ковровская государственная технологическая академия имени В.А. Дегтярёва»,

аспирант кафедры электротехники,

тел. (49232 )3-13-47,

e-mail: BalashovaSveta@list.ru