CC BY
130
УДК 669.14.018.291:621.891:539.219.3
АНОМАЛЬНОЕ УСКОРЕНИЕ ДИФФУЗИИ В ПОВЕРХНОСТНЫХ СЛОЯХ СТАЛЕЙ ПРИ ТРЕНИИ СКОЛЬЖЕНИЯ СО СМАЗОЧНЫМ МАТЕРИАЛОМ
П.И. Маленко, А.Ю. Леонов
Показано, что аномальное ускорение диффузии в поверхностных слоях сталей при трении скольжения со смазочным материалом вызвано высокочастотными температурными воздействиями в виде термических ударов, определяемых морфологией контактирующих поверхностей. В результате воздействия термических ударов атом (атом отдачи) покидает кристаллическую решетку, приобретая кинетическую энергию Екин с образованием пары “междоузельный атом-вакансия”. Число атомов отдачи превышает число вакансий в 102 ... 106раз, поэтому механизм диффузии будет носить преимущественно междоузельный характер. Количество одновременно появившихся атомов отдачи п>>1 в совокупности с столкновениями атомов эстафетного типа позволяет говорить о кооперативном механизме диффузии, кинетические характеристики которого предполагается определить при моделировании диффузионного процесса методом молекулярной динамики.
Ключевые слова: аномальное ускорение диффузии, трение со смазкой, термический удар, термоупругие напряжения, механизмы диффузии, коэффициент диффузии.
Экспериментальные исследования фазового состава вторичных структур, образующихся в поверхностных слоях сталей при трении скольжения со смазочным материалом, выявили аномально высокую скорость диффузионного массопереноса в направлении поверхности трения [1, 2]. Подобная аномалия определяется морфологией контактирующих поверхностей в виде суб- и микрошероховатостей и силой трения, в результате чего в зоне трения возникают механические и термические удары.
Динамическая реакция среды на механические удары заключается в возникновении вязких контактных напряжений т £. Динамическая реакция поверхностных слоев на термические удары отражается в виде термоупругих напряжений: диффузионной (тепловой) волны и продольной упругой волны расширения.
Цель работы - исследовать влияние динамической реакции среды на диффузионные процессы в поверхностном слое стали при трении скольжения со смазочным материалом.
Методика проведения исследований. Исследования производились расчетно-аналитическими методами на основе теории термоупругости, теории теплопроводности и теплопередачи, динамической теории кристаллической решетки с последующим сопоставлением полученных результатов с экспериментальными данными.
Результаты исследований и их обсуждение. Рассмотрим структурный состав поверхностного слоя стали на уровне кристаллической решетки. Контактные и упругие напряжения, проникая вглубь поверхностного слоя, воздействуют на диффузионную подвижность атомов в решетке. Оценка влияния упруго-волновых воздействий на динамику поведения атомов может быть осуществлена с помощью потенциалов парных взаимодействий (ППВ) п(г) [3]. На рис. 1 схематически представлен ППВ атомов в равновесном состоянии и при волновом нагружении механическими и термическими ударами.
Рис. 1. Потенциалы парного взаимодействия атомов n(r): а - в равновесном состоянии; б - волновом нагружении механическим и тепловым ударом
При механическом ударе атом смещается из положения rmin, соответствующего минимальной потенциальной энергии (глубине потенциальной ямы) Есв, на расстояние x. При термическом ударе смещение Х2, пропорциональное энергии удара Еуд, уменьшает энергию Есв. Если
Еуд < Есв, то атом смещается в решетке на расстояния x и Х2. При Еуд > Есв атом выбивается из решетки, приобретая некоторую скорость V, то есть кинетическую энергию Екин. В зависимости от величины Екин либо образуется пара “вакансия-междоузельный атом”, либо при Екин > 2Есв может произойти каскад столкновений. При подобном сценарии для описания диффузионного процесса необходимо проводить анализ перемещения отдельных атомов и эволюции исследуемого ансамбля атомов в целом, что возможно при моделировании процесса посредством использования метода молекулярной динамики (ММД) при известных ППВ.
Оценим диффузионную подвижность атомов в кристаллической решетке. Следуя работе [4] рассмотрим случай, когда силы взаимодейст-
вия между атомами будут центральными, то есть потенциальная энергия взаимодействия Епот пропорциональна квадрату расстояния между атомами. При расширении-сжатии изменение Епот на одну ячейку (атом) равняется работе, совершаемой над ячейкой силами давления р
АЕпот = 3е^ар, (1)
где е = — - относительная линейная деформация ячейки; Е - модуль уп-
Е
ругости, Уа - объем ячейки.
Результаты расчетов по зависимости (1) для кристаллической решетки а - Ев, приведенные в таблице, показывают, что для исследуемого диапазона давлений Епот < Есв = 4,5 эВ. Следовательно, атом не покинет решетку за один термический удар, однако, в течение времени релаксации 1р = 10-4...10-6 с он будет оставаться в положениях Х1 и х^. Так как
—6 —8
время между ударами т = 10 ...10 с [2], то атом выйдет из решетки за п
ударов с определенной скоростью V. Подобные атомы называются атомами отдачи.
Средняя скорость V выхода атома из решетки определяется из следующих соображений. Затухание диффузионной (тепловой) волны происходит на глубине поверхностного слоя И »1 мкм [5], в то время как для звуковой волны И »1 мм [6]. Помимо этого вязкие напряжения т^ значительно меньше давления Р, следовательно, влиянием диффузионной волны и вязких напряжений на среднюю скорость V можно пренебречь.
Тогда плотность мощности звуковой волны (поток энергии в единице объема) Ж можно определить по зависимости [6]
Ж = Р • Vзв, (2)
где V - скорость звуковой волны.
Мощность, приходящаяся на одну решетку Жреш = Ж • ^яч.
В свою очередь
Жат = Еупр • V, (3)
ЭФ( х)
Эх
парного взаимодействия.
Для определения кинетических характеристик атомов отдачи необходимо выбрать конкретный тип ППВ. Методика расчета ППВ сама по себе является сложной задачей. В этой связи обычно используют стандартные ППВ, хорошо зарекомендовавшие себя для определенных классов задач. Для рассматриваемого случая целесообразно представить потенциал как результат “стыковки” потенциала Борна-Майера (ПБМ) с потенциалом Леннарда-Джонса (ПЛД) (рис. 1).
где Еупр =---- упругая сила, действующая на атом; Ф( х) - потенциал
Известия ТулГУ. Технические науки. 2013. Вып. 10 Расчетные значения параметров (сталь 30ХН2МФА)
Параметр Т, к
470 570 670 770
P, ГПа 4,8 8,0 15,6 21,0
%k, ГПа 1,8
Vа, м3 0,8410-29
E 0,023 0,038 0,074 0,100
АБпот, эВ 0,05 0,14 0,54 0,94
N 90 32 9 5
ш-10-12, 1/с 3,20 3,65 4,15 4,72
F Н 1 упр> Н 0,410-9
1 • 09, м 1,8 1,6 1,4 1,24
W • 10-14, Вт/м2 0,12 0,20 0,39 0,52
£яч (а-Fe), м 4,1310-20
Wат • 106, Вт/ат 0,5 0,83 1,6 2,15
Vат • 10-3, м/с 1,25 2,08 4,0 5,4
Д • 106, м2/с 2,25 3,3 5,6 6,7
Eкин, Дж/эВ 0,76 10 - 20 0,03 2,2 10 - 20 0,14 7,85 10 - 20 0,49 1,43 •10-19 0,89
Потенциал Борна-Майера является потенциалом отталкивающего типа и имеет вид
п(г) = A • e-r/b, (4)
где A - предэкспоненциальный множитель, Ь = —- коэффициент,
(^2 )
имеющий размерность длины [м], а§ - боровский радиус водорода, z\ и Z2 - порядковые номера элементов в таблице Менделеева.
В рассматриваемом “склеенном” потенциале (рис.1) потенциал Борна-Майера носит вспомогательный характер и служит для нахождения параметра о в основном потенциале Леннарда-Джонса (5). При г >о
энергия межатомного взаимодеиствия становится характеристикой данного вещества. К примеру, для пары Ев - Ев • о = 2,5 • 10-10 м.
п(г ) = 4Е(
св
о
(5)
г • = 21/6 •о
'Ш1И ^ ^ •
В результате дифференцирования зависимости (5) получим следующее соотношение
Еупр = -4ЕСв 0,0о35, (6)
где Есв = 4,5...5 эВ для пары Ев - Ев [7]. Значение Еупр приведено в таблице. Там же даются искомые значения скорости V атомов отдачи.
Оценочные значения коэффициента диффузии находятся по формуле Д = V • 1, где 1 - длина звуковой (тепловой) волны, воздействующей на кристаллическую решетку.
Исходя из известных модельных представлений, основанных на использовании коэффициента теплопроводности 1 механизм теплопроводности в веществах, находящихся в различных агрегатных состояниях, различен. Соответственно и величины длин волн сильно разнятся.
С точки зрения авторов работы [8] для определения 1 вместо употребляемого коэффициента теплопроводности следует использовать значение коэффициента температуропроводности “а”. В этом случае физическое содержание понятия коэффициента “а” эквивалентно скорости тепловой (звуковой) волны на расстоянии, соответствующем ее длине. Поэтому используем для расчетов следующие зависимости
V2 а
ю = ^ и 1 = —, (7)
а Vзв
где ю - частота колебаний решетки (см. таблицу).
Там же приводятся значения коэффициентов диффузии Д и кинетической энергии, приобретаемыми атомами отдачи, рассчитанные по формуле екин=.
Рассчитанные значения коэффициентов диффузии соответствуют коэффициентам диффузии, вызванным механическими ударными взаимодействиями, лазерной обработкой и другими процессами, приведенными в работе [9].
Что же касается преимущественного направления диффузионного потока в сторону поверхности трения (направление оси 2), то оно связано с переносом атомов отдачи волной расширения с отрицательным давлени-
ем и напряжениями расширения оz > <
Выводы
1. Выход атомов из кристаллической решетки возможен за несколько ударов суб- и микрошероховатостей в пределах времени релаксации
^рел .
2. Кинетической энергии атомов отдачи достаточно для выбивания других атомов из решетки, поэтому следует ожидать появления большого количества междоузельных атомов.
3. Концентрация диффузионно возбужденных атомов в 102 ... 106 раз превышает концентрацию вакансий. Следовательно, механизм диффузии будет определяться междоузельными атомами.
4. Процесс диффузии в поверхностной зоне будет носить кооперативный характер, в связи, с чем имеет смысл говорить о кооперативном механизме диффузии. Термин “кооперативный” означает, во-первых, что практически одновременно в результате термических ударов образуется большое количество атомов отдачи. Во-вторых, первичные атомы отдачи, имея достаточную энергию, будут сталкиваться с другими атомами, то есть диффузия будет носить эстафетный характер. Особенности кооперативного механизма диффузии предполагается в дальнейшем установить при моделировании диффузионных процессов с помощью метода молекулярной динамики.
Работа представлена на второй Международной Интернет-конференции по металлургии и металлообработке, проведенной в ТулГУ 1 мая - 30 июня 2013 г.
Список литературы
1. Любарский И.М., Палатник Л.С. Металлофизика трения. М.: Металлургия, 1976. 176 с.
2. Маленко П.И., Зеленко В.К., Левин Д.М. Температурные поля и эксплуатационные свойства пар трения скольжения со смазочным материалом / Под ред. Ю.Н. Дроздова. М.: Машиностроение, 2011. 239 с.
3. Особенности процесса динамического разрушения металлов при воздействии теплового удара, вызываемого импульсами проникающих излучений и мощных импульсов лазерного излучения / Учаев А. Я. [и др.]. // VII Забабахинские научные чтения. Снежинск. 2003. С. 1-5.
4. Борн М., Хуан Кунь Динамическая теория кристаллических решеток. М.: Изд-во иностр. лит., 1958. 488 с.
5. Карташов Э.М., Ремизова О.И. Модельные представления термического удара при импульсных и пульсирующих тепловых нагрузках на
основе обобщенного уравнения энергии // Математическое моделирование. 2005. Т. 17. № 4. С. 81-95.
6. Ландау Л.Д., Лифшиц Е.М. Теория упругости. М.: Физматлит, 2001. 259 с.
7. Шульце Г. Металлофизика. М.: Мир, 1971. 503 с.
8. Коломейцев В.В., Коломейцева Е.Ф., Суворов С.А. Феноменологическая теория температуропроводности // Огнеупоры и техническая керамика. 2001. № 4. С. 35-36.
9. Bekrenev A.N., Kamashev A.V. Features of phase transformations passing abd mass transport in metals under intensive external reactions // J. of Physics and Chemistry of Solids. 2001. Vol. 62. P. 647-651.
Маленко Павел Игоревич, канд. техн. наук, доц., malenkoatsn. tula.ru, Россия, Тула, Тульский государственный университет,
Леонов Андрей Юрьевич, аспирант, [email protected], Россия, Тула, Тульский государственный университет
ABNORMAL ACCELERATION OF DIFFUSION IN THE SURFACE LAYERS OF STEELS
AT SLIDING FRICTION LUBRICANT
P.I.Malenko, A. Yu.Leonov
It is shown that the anomalous acceleration of diffusion in the surface layers of steels at sliding friction lubricant caused by high-temperature effects in the form of thermal shock, determined by the morphology of the contacting surfaces. As a result of exposure to thermal shock atom (atom of bestowal) leaves the crystal lattice, gaining kinetic energy Ekin with the formation of pairs of “interstitial atom-vacancy ”. The number of the recoil atoms exceeds the number of vacancies in 102 ... 106 times, so the mechanism of diffusion will wear mostly interstitial character. Number of simultaneously appeared recoil atoms n > > 1 together with the collisions of atoms relay type allows to speak about cooperative mechanism of diffusion, kinetic characteristics which is intended to set when modeling the diffusion process by the method of molecular dynamics.
Key words: anomalous acceleration of diffusion, friction with lubrication, thermal shock, thermoelastic stresses, mechanisms of diffusion, diffusion coefficient.
Malenko Pavel Igorevich, candidate of technical science, docent, malenko atsii. tula. ru, Russia, Tula, Tula State University,
Leonov Andrey Yuryevich, postgraduate, leon13-23@mail. ru, Russia, Tula, Tula State University
