Моделирование кооперативного механизма диффузии методом молекулярной динамики

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

Моделирование кооперативного механизма диффузии методом молекулярной динамики

П.И. Маленко

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

Ключевые слова: кооперативный механизм диффузии, термический удар, вторичные структуры, структурно-фазовые превращения, трение скольжения со смазочным материалом, метод молекулярной динамики.

диффузия трение молекулярный динамика

Введение

При работе пары трения в поверхностных слоях происходят структурно-фазовые превращения с образованием вторичных структур (ВС), которые определяют их эксплуатационные свойства и, в частности, склонность к схватыванию, обусловленную резким повышением коэффициента трения [1]. Образование ВС связано с диффузионными и полиморфными (бездиффузионными) процессами в поверхностных слоях. При этом, происходит значительное повышение коэффициента диффузии (D = 10-5 ч 10-8 м2/c), что подтверждено экспериментально.

Это объясняется морфологией поверхностей трения, а именно, наличием суб- и микрошероховатостей, вызывающих пульсирующее изменение температуры, то есть термоциклирование [2], которое, в свою очередь, является причиной возникновения термических ударов (ТУ), а следовательно, резкого повышения давления p в поверхностных слоях [3]. Диффузия атомов в твердом теле - это сложный процесс и в этой связи нет единой точки зрения на диффузию, не обладающей внутренними противоречиями [4]. Широко использующийся для простых металлов вакансионный механизм, основанный на “скачкообразной” диффузии, по вакансиям многократно подтверждался, но не доказывался в различных исследованиях. Для других механизмов (обменный, междоузельный, краудионный и другие) “первопричинные” расчеты занижают значения скоростей диффузии в сотни раз и носят эвристический характер.

В настоящее время, при изучении диффузионных процессов применяют методы молекулярной динамики (ММД), позволяющие моделировать исследуемые процессы благодаря используемому механизму, вносимой априорной - физической информации. Процесс диффузии в графическом виде, представляется либо в виде объемного изменения положения атомов, либо в виде траекторий перемещения, диффундирующих атомов во времени и в пространстве.

В работах [5, 6] приведены результаты моделирования ММД процесса перемещения атомов под действием ТУ в поверхностной зоне теплостойкой стали 25Х3М3НБЦА с никотрированным покрытием. Анализ траекторий перемещений атомов позволил определить механизм диффузии как кооперативный (эстафетный) с коэффициентами диффузии D = 10-6 ч 10-8 м2/c.

Материалы и методы исследований

Экспериментальные методы исследования диффузионных процессов и полиморфных превращений в поверхностной зоне трения - сфероионная микроскопия (EIM) и сканирующий туннельный микроскоп (CTM) используются для проникновения в суть поверхностной диффузии, но доминирующий механизм диффузии с их помощью не может быть идентифицирован.

Металловедческие исследования, основанные на обработке экспериментальных данных и рассматривающих металлы в континуальном приближении, только подтверждают возможность наличия различных механизмов диффузии, но не доказывают их существование.

Применяемый для анализа механизмов диффузии на атомном уровне, то есть на уровне кристаллических решеток (среда дискретная) ММД основан на численном интегрировании уравнений движения Ньютона для системы из N частиц (атомов). Число N выбирается обычно от нескольких сотен до десятков тысяч. Важная особенность метода - это задание закона межчастичного взаимодействия - потенциала парного взаимодействия (ППВ). Наибольшую точность дают расчеты ППВ «из первых принципов» (ab initio), в том числе и квантово-химические методы. Однако, вследствие больших временных затрат они базируются на ограниченном количестве атомов (n < 1000) и не могут быть перенесены на большие атомные системы. В этой связи в последние годы широко используются полуэмпирические методы, в частности, формализм метода «погруженного атома» (EAM). Дальнейшая процедура получения наблюдаемых физических величин (в частности, коэффициента диффузии D) состоит в статистической обработке (усреднении) функции диффузии в зависимости от микропеременных - координат атомов r(ф) и скоростей х(ф) атомов системы (кристаллита).

Цель данной работы состоит в изучении кооперативного механизма диффузии под действием ТУ, используя формализм EAM, с последующей проверкой результатов моделирования путем сравнения с экспериментальными данными. Объектом исследования послужила теплостойкая сталь 25Х3М3НБЦА с низкотемпературным карбонитридным покрытием (процесс никотрирования). Данная сталь используется для изготовления узлов автоматики стрелково-пушечного вооружения, работающих в условиях трения скольжения с ресурсным смазыванием.

Результаты и их обсуждение

Молекулярно-динамические расчеты проводились с помощью стандартного пакета программ XMD, основные характеристики которого описаны в работе [7]. Данный пакет программ позволяет проводить расчеты с циклическими или свободными граничными условиями в приближении многочастичного потенциала EAM.

Сущность кооперативного механизма диффузии состоит в том, что диффузионный процесс осуществляется не за счет “скачков” отдельных атомов, а в результате коллективного перемещения N атомов. Так как энергия связи атомов в решетке Есв>Епот, то атом выходит из кристаллической решетки за n ударов. Значения n рассчитываются из уравнения

. (1)

Последние удары в серии из n ударов можно отнести к термическим флюктуациям, способствующим переходу атома через потенциальный барьер. Известно, что термические флюктуации играют решающую роль на многих этапах эволюции структуры кристаллической решетки, включающей и диффузионные процессы. Диффузионные процессы определяются элементарными атомными перемещениями. Вместе с тем свойства макроскопических тел обусловлены термически активированными перестройками дефектной структуры и кинетикой ансамблей дефектов в структуре [8].

Проанализируем в этой связи механизм выхода атома из кристаллической решетки. Для этих целей был сформирован кристаллит для б-Fe при Т = 500 К из 250 атомов (рис. 1,а). Программное обеспечение позволило осуществить перемещение (диффузию) атома №219 в решетке и зафиксировать мгновенные температуры как в зоне диффундирующего атома, так и в решетке ближайших к нему соседей (атомы №№209, 210, рис. 1,б). Температуры атомов представлены на рис. 1,в.

Рис. 1. Моделирование температур атомов при кооперативном механизме диффузии: а - модельный кристаллит б-Fe;

б - расположение атомов относительно кристаллографического направления <001>; в - изменение мгновенных температур моделируемых атомов

Из рис. 1,в следует, что температура диффундирующего атома №219 в момент выхода из решетки достигает Т = 1600 К, а температуры соседних недиффундирующих атомов не превосходят Т = 800 К. Следовательно, энергия связи Есв в данном месте должна уменьшаться. Для ее определения используем кинетическую теорию прочности [9]. Долговечность кристаллита ф определяется зависимостью

, (2)

а величина прочности на разрыв

, (3)

- где U0 = 4,28 эВ (для б-Fe); k - постоянная Больцмана; - время, соответствующее периоду тепловых колебаний атомов в решетке, ф0 = 2·10-13 с; г ? Vат, для б-Fe г = 8,38·10-30 м3.

Откуда у = удиф = 8,75·109 Па (Т = 600 К).

Расчеты по зависимостям (2) и (3) дают следующие результаты: ф = 2,14·1011 с, у = 1,785·109 Па.

Так как у < удиф, то диффундирующий атом, не совершая скачков, будет покидать решетку, раздвигая соседние атомы.

Непосредственно процесс моделирования кооперативного механизма диффузии производился согласно работе [10]. Траектории частиц, получаемых ММД, подчиняются уравнениям

, , , (4)

где и - воздействие численных ошибок вычислений (погрешностей) численного интегрирования дифференциальных уравнений

, , (5)

- где - сила, действующая на атом.

При наличии внешнего воздействия, например, термического удара с давлением Р

, (6)

- где S - площадь центрального сечения атома.

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

Программа, описанная в работе [7] позволяет вычислять средние значения траекторий атомов , обладающие следующими свойствами погрешностей вычислений:

- они не зависят от точности, с которой рассчитывается траектория;

- могут быть вычислены при определенном разбиении траектории на интервалы фk;

- определяются величиной , то есть полным временным интервалом, для которого рассчитывается траектория, а не числом шагов на траектории.

Особого внимания заслуживает обеспечение термостабилизации системы. При перемещении атомов в результате диффузии может происходить неконтролируемое повышение температуры Т. Для контроля Т к силе F, действующей на каждый атом, на каждом шаге численного интегрирования вводится дополнительное слагаемое - случайная сила, равномерно распределенная в интервале , где G - спектральная плотность «белого» шума, h - шаг интегрирования.

На сформированный из кристаллических решеток кристаллит из N атомов могут накладываться как свободные, так и циклические граничные условия. Свободные условия используются для определения коэффициента диффузии D для конкретного кристаллита. Циклические граничные условия необходимы для кратного увеличения количества атомов. На границе кристаллита фиксируется текущее положение атомов. Затем рассматривается следующий аналогичный кристаллит, для которого начальное положение атомов - это положение атомов на границе предыдущего кристаллита.

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

В качестве примеров далее представлены идеальные кристаллиты нитридов Fe4N и Fe3N. Кристаллиты были сформированы из N = 1000 атомов. На рис. 2,а представлена исходная структура фазы Fe4N, на рис. 2,в - собственно диффузионный процесс.

Первоначально атомы N перемещаются более активно, чем атомы Fe (рис. 2,б). При этом происходит как выход атомов N из кристаллита в окружающую среду, так и их скопление непосредственно у поверхности трения. В дальнейшем атомы N полностью покидают кристаллит и активизируется процесс диффузии г-Fe (рис. 2,в). Атомы г-Fe распределяются по поверхности кристалла в виде аморфной структуры. Учитывая высокую степень ионизации кислорода (сродство к электрону), происходит процесс образования вторичных структур - Fe2О3.

Рис. 2. Расположение атомов в кристалле ГЦК решетки фазы г-Fe4N (Т = 600 К, удиф = 8,9 ГПа, n = 1000 атомов):

а - исходное состояние; б - диффузия N; в - диффузия г-Fe

Аналогичным образом происходит процесс диффузии в фазе защитного никотрированного слоя е-Fe3N. На рис. 3 представлены исходная структура фазы е-Fe3N (рис. 3,а) и собственно диффузионные процессы (рис. 3,б) под действием ТУ.

Рис. 3. Расположение атомов в кристалле ОЦК решетки фазы е-Fe3N (Т = 800 К, удиф = 8,9 ГПа, n = 2400 атомов):

а - исходное состояние; б - диффузия N

На рис. 4,а,б представлен кристаллит с дефектом - межзеренной границей “б-Fe-FeS2”. В данном случае атомы серы, как более лёгкие, диффундируют из межзеренной границы гораздо быстрее, чем атомы б-Fe.

Рис. 4. Расположение атомов в кристалле решетки фазы FeS2

(Т = 600 К, удиф = 8,9 ГПа, n = 2620 атомов):

а - исходное состояние; б - диффузия S2

Определенный интерес представляет анализ кооперативного механизма диффузии для фазы Mo2C. Диффузия через межзеренную границу “г-Fe-Mo2C” представлена на рис. 5.

Рис. 5. Расположение атомов на межзеренной границе “г-Fe-Mo2C” (Т = 900 К, удиф = 14 ГПа, n = 1000 атомов):

а - исходное состояние; б - диффузия Mo и C

Атомы C диффундируют с большой скоростью, выходя из межзеренной границы и фазы г-Fe. В то же время атомы Mo “протискиваются” между атомами Fe по извилистой траектории.

По аналогичному сценарию происходит и процесс образования других вторичных структур. Из представленных схем очевидно, что механизм диффузии является кооперативным (эстафетным), когда одновременно диффундирует большое количество атомов (n > 30) и возможны их столкновения. Следовательно, можно утверждать, что отсутствует доминирующее влияние вакансионного механизма.

Оценку коэффициента диффузии D для г-Fe и N (фаза г-Fe4N) произведем по формуле диффузии Эйнштейна

, (7)

- где ф - время диффузии; ?r2=x2+y2+z2 - квадрат перемещения атомов относительно осей координат.

На рис. 6 показана зависимость среднего квадрата смещений атомов от времени. Полученные зависимости нелинейные, причем перемещение атомов N по отношению к атомам г-Fe происходит ускоренно.

Рис. 6. Зависимость среднего квадрата смещений атомов

от времени для фазы Fe4N (Т = 600 К, удиф = 8,9 ГПа)

В таблице 1 приведены коэффициенты D для различных временных шагов, рассчитанные по формуле (7).

Полученные значения коэффициентов Dг-Fe несколько завышены от экспериментальных Dэксп ? 10-6 … 10-7 м2/c. Расчетные значения D были получены для идеального кристалла без учета влияния дефектов решетки.

Для сравнения дадим оценку коэффициента диффузии D для б-Fe и S (фаза FeS2), то есть в качестве дефекта используется межзеренная граница.

Таблица 1

Значения коэффициентов диффузии D для фазы г-Fe4N

при Т = 600 К и давлении удиф = 8,9 ГПа

На рис. 7 и 8 показана зависимость среднего квадрата смещений атомов по времени для Fe и S соответственно для межзеренной границы “б-Fe-FeS2”.

Рис. 7. Зависимость среднего квадрата смещений атомов

от времени для Fe

Для сбора и последующего анализа данных, необходимых при вычислении коэффициента диффузии, взяты по пять частиц Fe и S, находящихся около центра сферы приложения внешней силы. Среднее значение на рис. 7 и 8 показано красной линией.

В таблице 2 приведены коэффициенты D для фазы FeS2, рассчитанные по формуле (7).

Рис. 8. Зависимость среднего квадрата смещений атомов

от времени для S

Таблица 2

Значения коэффициентов диффузии D для фазы FeS2

при Т = 600 К и давлении удиф = 8,9 ГПа

Аналогично были рассчитаны значения коэффициентов диффузии для “г-Fe-Mo2C” (табл. 3).

Таблица 3

Значения коэффициентов диффузии D для фазы Mo2C

при Т = 900 К и давлении удиф = 14 ГПа

Рассмотрим влияние кооперативного механизма самодиффузии Fe на структуру б-Fe, содержащую дефект в виде вакансии. Для более наглядного представления на рис. 9 приведены траектории диффундирующих атомов Fe, спроектированных на плоскость <101>.

Рис. 9. Изменение структуры исходного кристаллита б-Fe

с одной вакансией при кооперативном механизме диффузии

(Т = 900 К, n = 104 шагов)

Из приведенного рисунка следует, что в процессе диффузии образуется вакансионный кластер, перерастающий затем в вакансионную пору.

Выводы

1. Моделирование диффузионных процессов в сформированных из кристаллических решеток кристаллитах показало, что как в поверхностных фазах, так и в фазах подслоя под действием ТУ в процессе диффузии одновременно участвует большое количество атомов (n > 30 атомов), то есть речь идет о кооперативной диффузии. Данная диффузия носит эстафетный характер, когда в процессе перемещения атомов происходят их столкновения и передача кинетической энергии.

2. При кооперативной диффузии преимущественная роль вакансионной диффузии, характерная для обычного термического воздействия в виде плавного нагрева значительно уменьшается.

3. Для поверхностных фаз, например, для фазы Fe4N имеет место повышенный коэффициент диффузии, что связано с уменьшенной поверхностной энергией кристаллита по сравнению с объемной энергией, характерной для фаз подслоя.

4. Коэффициенты диффузии для примесей (N, S, С) значительно превосходят коэффициенты самодиффузии Fe как в поверхностных слоях, так и в подслое.

5. Кооперативный механизм способствует порообразованию за счет увеличения вакансионных кластеров.

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

1. Маленко П. И., Зеленко В. К., Левин Д. М. Температурные поля и эксплуатационные свойства пар трения скольжения со смазочным материалом / Под ред. Ю. Н. Дроздова. М.: Машиностроение. 2011. 239 с.

2. Маленко П. И. Исследование температур на дискретных субшероховатых поверхностях при трении скольжения со смазочным материалом // Вестник машиностроения. 2011. № 7. С. 38-42.

3. Дроздов Ю. Н., Маленко П. И. Структурно-фазовые превращения в поверхностных слоях сталей при трении скольжения // Трение и износ. 2014. Т.35. № 1. С. 87-98.

4. Чудинов В. Г. Кооперативный механизм самодиффузии в металлах // Журнал технической физики. 2000. Т. 70. Вып. 7. С. 133-135.

5. Дроздов Ю. Н., Маленко П. И. Применение метода молекулярной динамики для анализа кинетики структурно-фазовых превращений в поверхностных слоях упрочненных теплостойких сталей при трении скольжения с ресурсным смазыванием // Фундаментальные исследования и инновационные технологии в машиностроении: тез. докл. III Международ. науч. конф. М.: Издательский дом Спектр. 2014. С. 90-92.

6. Дроздов Ю. Н., Маленко П. И. Исследование методом молекулярной динамики диффузионных процессов в поверхностных слоях теплостойких сталей с покрытием при трении скольжения с ресурсным смазыванием // УРАЛ-ТРИБО: материалы Международ. науч.-техн. конф. ученых, ведущих специалистов в обл. трибологии и аспирантов. Челябинск: Издательский центр ЮУрГУ. 2014. С. 7-9.

7. Rifkin J. XMD Molecular Dinamics Program [Electronic resource] // University of Connecticut, Center for Materials Simulation, Storrs, CT, 2002. 104 P. URL: http://xmd.SourceForge.net/ (accessed: 18.02.2011).

8. Термически активированные процессы в кристаллах: сб. статей: пер. с англ. / Под ред. А. Н.Орлова. М.: Мир, 1973. 212 с.

9. Физика прочности и пластичности: сб. статей / Под ред. С. Н. Журкова. Л.: Наука, 1986. 152 с.

10. Маленко, П. И. Исследование методом молекулярной динамики диффузионных процессов в поверхностных слоях никотрированных теплостойких сталей в процессе трения скольжения с ресурсным смазыванием // Вестник ЮУрГУ. Серия: Машиностроение. 2014. Т. 14, № 3. С. 22-29.

Размещено на Allbest.ru