CC BY
23
The formation of thick layers in arc welding in argon an experimental model nano-modification special hardfacing flux-cored wire is a new generation multi-layers of ultra-fine and nanomodified songs some blanks for the manufacture of dies for hot forming is considered. The influence of the current, the feed rate of the electrode wire and the temperature of the heated substrate on the formation of the melt cushion and the influence of duration of cooling on the mechanical properties of steel 40ХН is analyzed. Defined thermal cycles on the surface of the substrate at various values of temperature of preliminary heating.
Key words: arc welding, flux-cored wire, thermo-physicalproperties.
Erofeev Vladimir Aleksandrovich, candidate of technical sciences, docent, va erofeevamail.ru, Russia, Tula, Tula State University,
Zakharov Sergey Konstantinovich, candidate of technical sciences, docent, zzzsk19 7l ayandex. ru, Russia, Tula, Tula State University,
Zotova Svetlana Vyacheslavovna, assistant, zsvzot'a yandex. ru, Russia, Tula, Tula State University
УДК 621.785.533; 539.219.3; 548.4
МОДЕЛИРОВАНИЕ ПРОЦЕССА УВЕЛИЧЕНИЯ ДИФФУЗИОННОЙ ПОДВИЖНОСТИ АТОМОВ В ВАКАНСИОННОМ МЕХАНИЗМЕ
ДИФФУЗИИ
П.И. Маленко, С.В. Ильичев, К. Д. Релмасира, А.Ю. Леонов
На основе моделирования методом молекулярной динамики приводится количественная оценка процесса температурного увеличения диффузионной активности атомов. Предлагаемая оценка основана на выявленном в процессе моделирования образовании пустот (дилатонов) и сгущений. В результате в бездефектном кристаллите будет возникать кооперативный механизм диффузии, в то время как в кристаллите с вакансиями будет иметь место сочетание вакансионного и кооперативного механизмов.
Ключевые слова: метод молекулярной динамики, дилатон, сгущение, коэффициент перегрузки, вакансионный механизм, кооперативный механизм диффузии.
Термин «увеличение диффузионной подвижности атомов» является общепринятым в металловедении и качественно объясняет предрасположенность к диффузии, например, при увеличении температуры. Однако для количественной его оценки требуются иные методы, в частности, метод молекулярной динамики (ММД).
124
Сущность вакансионного механизма диффузии состоит в том, что вследствие температурного увеличения диффузионной подвижности атомов некоторая их часть преодолевает потенциальный барьер посредством прыжка (энергию связи Есв в кристаллической решетке). Считается, что преодоление потенциального барьера определяет начало диффузионного процесса вакансионного типа. Таким образом, увеличение диффузионной подвижности атомов постулируется. В предлагаемой статье с использованием ММД делается попытка раскрыть физическую сущность данного процесса.
Материалы и методы исследований. Экспериментальные методы исследования диффузионных процессов в поверхностной зоне трения -сфероионная микроскопия (EIM) и сканирующий туннельный микроскоп (CTM) - используются для проникновения в суть поверхностной диффузии, но доминирующий механизм диффузии с их помощью не может быть идентифицирован.
Металловедческие исследования, основанные на обработке экспериментальных данных и рассматривающих металлы в континуальном приближении, только подтверждают возможность наличия различных механизмов диффузии, но не доказывают их существование.
В этой связи возникла необходимость в анализе механизмов диффузии на атомном уровне, то есть на уровне кристаллических решеток (среда дискретная). Для решения подобных задач используется ММД, основанный на численном интегрировании уравнений движения Ньютона для системы из N частиц (атомов). Число N выбирается обычно от нескольких сотен до десятков тысяч. Важная особенность метода - это задание закона межчастичного взаимодействия - потенциала парного взаимодействия (ППВ). Наибольшую точность дают расчеты ППВ «из первых принципов» (ab initio), в том числе и квантово-химические методы. Однако вследствие больших временных затрат они базируются на ограниченном количестве атомов (n < 1000) и не могут быть перенесены на большие атомные системы. В этой связи в последние годы широко используются полуэмпирические методы, в частности, формализм метода «погруженного атома» (EAM).
Цель данной работы состоит в изучении процесса увеличения диффузионной подвижности атомов путем моделирования вакансионного механизма диффузии с использованием ММД. Объектом исследования послужила теплостойкая сталь 25Х3М3НБЦА с низкотемпературным карбонит-ридным покрытием (процесс никотрирования). Данная сталь используется для изготовления узлов автоматики стрелково-пушечного вооружения, работающих в условиях трения скольжения с ресурсным смазыванием.
Результаты и их обсуждение. Молекулярно-динамические расчеты проводились с помощью стандартного пакета программ XMD, основные характеристики которого описаны в работе [1].
Исследовали увеличение подвижности атомов в a-Fe, объемное содержание которого в никотрированном покрытии составляет более 80 %. Для моделирования процесса был сформирован исходный кристаллит из М = 1458 атомов a-Fe. Сущность эксперимента состояла в исследовании изменения структуры кристаллита в диапазоне температур Т = 300... 900 К за определенный промежуток времени t = NDt, где N = 5-10 - количество шагов, Ат = 10-16 с - шаг моделирования (т = 5-10-13 с). На первом этапе исследовалось изменение количества атомов п в кристаллической решетке в пределах координационного числа к (для a-Fe к = 8): n < к. На рис. 1, а представлены функции распределения в координатах «n - Р(п)» для иссле-
гп\
N
Из графиков сле-
>ViV У
дуемого температурного диапазона P(n)= lim
N
дует, что при комнатной температуре (Т = 300 К) разброс значений n, отличных от к, фактически отсутствует, то есть диффузионная подвижность атомов минимальна. При Т > 600 К имеет место разброс значений n от n = к < 8. При n < 8 происходит разрежение, то есть образование пустот в кристаллите, при n > 8 происходит сгущение (сжатие) атомов в кристаллической решетке, принимающее максимальное значение при Т = 900 К, то есть имеет место максимальная диффузионная подвижность атомов. На рис. 1,б представлено распределение P(n) для решеток, имеющих дефекты (вакансии).
4 5 6 7
4 5 6 7 8 9 Ю 11 12
а
б
Рис. 1. Вероятность расширения-сжатия атомов в кристаллической решетке для различных температур (1 - Т = 300 К, 2 - Т = 600 К, 3 - Т = 900 К): а - без дефектов; б - с дефектами (вакансиями)
На втором этапе полученные результаты трансформировались в двумерное изображение в координатах плоскости «х - у» (рис. 2). Программа моделирования рассчитывала изменения координат атомов под действием температур. На рис. 2 отчетливо видны пустоты и сгущения в расположении атомов. Для оценки влияния дефектов структуры (в данном случае вакансий) в кристаллит была введена вакансия (рис. 2, область 4).
Из полученных данных можно предположить, что увеличение диффузионной подвижности атомов можно количественно оценить по количеству и размерам пустот (и сгущений) в кристаллите. Предлагаемая далее методика количественной оценки увеличения диффузионной подвижности атомов базируется на термофлуктуационной теории механизма разрушения металлов под действием ударов абразивных частиц [2, 3]. Различие предлагаемой методики заключается в том, что она позволяет оценить не процесс абразивного изнашивания при ударе, а механизм диффузии и его параметры.
.4 69 .471 ->.473 ^475 Й77 1479 .481 .4 S3 485
,451 ,453 ■.455 ,457 ¿59 _4 61 .463 .4 65 ^67
,433 -435 ..437 .439 441 .443 -445 ,447 .4 4 9
,415 417 . .419 .4 23 .425 423--
1421 2 9 ,431
Д - ( V 3
399 ,403 -.405 ,407
379 ¿81 ¿83 ^387 ¿389 ¿S91 ,393 395
362 ¿64 £366 4\ а/ (36э\ 371 ч373 375 J77
.344 346 348 С. 3501 --- Д/' 354 35 6 ¿58 >360
,32 6 L 328 3 -J ,330 .334 ~ ~2 зза ¿340 ¿42
г '1
6 of 18
Atom displacements in alfa Fe
Diana Farkas's Fe potential
A = 2.87
box 9 9 9
1458 atoms in box
temp 900K clamp 900
vacancy - atom N368
5000 steps of cmd
RUN 1 STEP 5000 ORIENT Y/Z Thu Jul 21 23:21:25 2016 Page 1 of 1
Рис. 2. Изменение структуры исходного кристаллита в результате увеличения диффузионной подвижности атомов (Т = 900 К) в плоскости <100>: 1 - 30 мкм, 2 - 35 мкм,
3 - 50 мкм, 4 - 30 мкм
Основу методики составляет дилатонный механизм диффузионной подвижности атомов, описывающий образовавшиеся в кристаллите пустоты как совокупность непериодических колеблющихся атомов. Сущность дилатонного механизма отождествляется с переходом через потенциальный барьер не единичных атомов, а с разрушением дилатонов, что приводит к диффузии целой группы разориентированных атомов. Для описания дилатона необходимо определить деформацию межатомных связей [2]
aT, t
e0 =Тln_, (1)
0 3 t0
где a = 1,1710-5 1/град - тепловой коэффициент линейного расширения; т = 5-10-13 c - продолжительность воздействия; т0 = 1/ю0 = 1/(8-1012) = = 1,25-10" - период тепловых колебаний атомов в решетке a-Fe,
£0 = (3,24 ... 4,9)-10-3 - для Т = 600 ... 900 К.
Размер дилатона А определяется по зависимости
D =-0, (2)
s
где а = 2,87-10-10 м - параметр кристаллической решетки; E = 2Д5-1011
2 8 2 Н/м - модуль упругости; а = 6-10 Н/м - реальная прочность стали,
А =(3,3 ... 5)-10-10 м.
Коэффициент перегрузки х = А/a = 1,15 ... 1,75 зависит от наличия дефектов структуры. К примеру, для дилатонов № 1, 2 и 3 (см. рис. 2) размер дилатонов составляет 1,5а; 1,75а; 1,5а. Соответственно х = 1,5; 1,75; 1,5. В то же время для дилатона № 4 (с вакансией) А = 2,5а и х = 2,5. Коэффициент х означает склонность к диффузии атомов. Следовательно, атомы в дилатоне с вакансией диффундируют быстрее, что будет способствовать ускорению диффузионных процессов. Так как количество дилатонов п >> 1, то в данном случае атомы будут совершать не единичные прыжки через потенциальный барьер, характерные для вакансионного механизма, а групповые (п = 30 ... 40), что свойственно кооперативному (эстафетному) механизму [4]. При этом наличие вакансий в исходном кристаллите будет способствовать ускорению диффузионного процесса.
Выводы
1. На основе модельных представлений о процессе температурного увеличения диффузионной активности атомов предложена методика его количественной оценки.
2. Предлагаемая методика базируется на моделировании процесса методом молекулярной динамики на примере кристаллической решетки a-Fe.
3. Моделирование установило, что в результате температурного воздействия изменяется структура исходного кристаллита (плоскость <100>) с образование пустот и сгущений.
128
4. В пустотах (названных дилатонами) происходит деформация межатомных связей, приводящая к их перегрузкам, что в дальнейшем провоцирует начало диффузионного процесса.
5. В бездефектном кристаллите при количестве дилатонов n >> 1 возникает кооперативный механизм диффузии. При наличии в кристаллите вакансий будет иметь место сочетание вакансионного и кооперативного механизмов.
Список литературы
1. Rifkin J. XMD Molecular Dinamics Program [Electronic resource]. University of Connecticut, Center for Materials Simulation, Storrs, CT, 2002. 104 P.
2. Журков С. Н. К вопросу о физической природе прочности //ФТТ. 1980. Т. 22. Вып. 11. С. 3344-3349.
3. Сорокин Г. М., Григорьев С. П., Голова А. Г. О природе ударно-абразивного изнашивания сталей // Проблемы прочности. 1991. № 4. C. 7376.
4. Чудинов В. Г. Кооперативный механизм самодиффузии в металлах // Журнал технической физики. 2000. Т. 70. Вып. 7. С. 133-135.
Маленко Павел Игоревич, канд. техн. наук, доц., malenko a tsii.tula.ru, Россия, Тула, Тульский государственный университет,
Ильичев Сергей Владимирович, инженер-программист, sergil68 a mail.ru, Россия, Тула, ООО «Лаборатория Интеллект»,
Релмасира Капреси Джолан, асп., winnerclassagmail.com, Россия, Тула, Тульский государственный университет,
Леонов Андрей Юрьевич, асп., leonov.emwamail.ru, Россия, Тула, Тульский государственный университет
MODELING OF THE PROCESS OF INCREASING THE DIFFUSION MOBILITY OF ATOMS IN THE VACANCY MECHANISM OF DIFFUSION
P.I.Malenko, S. V.Ilyichev, K.J.Relmasira, A. Yu.Leonov
In article on the basis of modeling by the molecular dynamics method provides a quantitative assessment of the process of temperature increase of diffusion activity of atoms. The proposed assessment is based on identified in the modeling process, the formation of voids (diletantov) and clumps. As a result, in the defect-free crystallite will have the cooperative mechanism of diffusion, whereas in the crystal with vacancies will have a combination of vacancy and cooperative mechanisms.
Key words: molecular dynamics method, dilaton, concentration, overload factor, vacancy mechanism, komparativny the mechanism of diffusion.
129
Malenko Pavel Igorevich, candidate of technical sciences, docent, malen-koatsu. tula. ru, Russia, Tula, Tula State University,
Ilyichev Sergey Vladimirovich, engineer, sergil68 a mail. ru, Russia, Tula, LLC «Laboratory Intelligence»,
Relmasira Kapresi Jolan, postgraduate, winnerclassagmail.com, Russia, Tula, Tula State University,
Leonov Andrei Yurievich, postgraduate, leonov. emw a mail. ru, Russia, Tula, Tula State University
УДК 378.14
ВИРТУАЛЬНАЯ ЛАБОРАТОРНАЯ РАБОТА ПО ИЗУЧЕНИЮ ПРОЦЕССА ОСАДКИ ПРИ СТЫКОВОЙ КОНТАКТНОЙ СВАРКЕ
ОПЛАВЛЕНИЕМ
В. А. Ерофеев, И.Б. Пьянков, А. А. Мещеряков
Рассмотрена проблема организации лабораторного практикума при дистанционном обучении по технологии сварочного производства. В качестве примера рассматривается изучение процесса осадки при стыковой контактной сварке оплавлением. Задачей обучения является формирование навыка определения оптимальных параметров процесса методами компьютерного инженерного анализа, основанного на математическом моделировании объекта. При обучении предлагается построить зависимости критериев качества формирования соединения от параметров сварки при заданных параметрах стыка.
Ключевые слова: дистанционное обучение, виртуальная лабораторная работа, компьютерное моделирование.
Развитие дистанционных форм обучения требует создания виртуального лабораторного практикума. Задачами лабораторного практикума является формирование навыков исследования объектов изучения. При подготовке технологов сварочного производства основной задачей является обучение специалистов по формированию навыков оптимизации технологии, обеспечивающей качество сварных соединений. Решение этой задачи целесообразно выполнить методами компьютерного инженерного анализа, основанного на физико-математическом моделировании сварочного процесса.
В качестве примера рассматривается изучение процесса осадки при стыковой контактной сварке оплавлением.
130
