Научная статья на тему 'Моделирование охлаждения аморфного сплава'

Моделирование охлаждения аморфного сплава Текст научной статьи по специальности «Физика»

CC BY
90
42
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по физике, автор научной работы — Волкова И. Б., Баранов М. А., Баянкин В. Я.

The molecular dynamics technique is used to investigate the cooling process of the amorphous alloy Fe80B20. The computer simulation results are compared to the real experimental data. A mechanism elucidating the resulting regularities is proposed.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Волкова И. Б., Баранов М. А., Баянкин В. Я.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

SIMULATION OF AMORPHOUS ALLOY COOLING

The molecular dynamics technique is used to investigate the cooling process of the amorphous alloy Fe80B20. The computer simulation results are compared to the real experimental data. A mechanism elucidating the resulting regularities is proposed.

Текст научной работы на тему «Моделирование охлаждения аморфного сплава»

и е 2 , например, ^ = 10 е 2 . Исходя из определения предела пропорциональности, с учетом того, что на упругом участке кривые нагружения при различных е совпадают, получим из (5) такое соотношение для скорости упрочнения на пределе пропорциональности:

а0 , , А „

— -М —— (с .1 -ъ0) -

аг е1

-

- М - -Ч0щ.2 - Оо)"

г 2

откуда следует

—1 /■_ \П _ — 1 \Ш

—(ъпц.1 - ъо) --г— (ъпц.2 - ъо) ,

г1 г 2

Ъпц.1 - Ъ0 1 - _г1 - 10

0пц.2 - ъ0 ) г2

°пцЛ - 00 - ^Ю,

0пц.2 - ъ0 илиъпц.1 - П10"0пц.2 - ъ0(п10 - ^

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

4/т «1 ,8, что подтверждается данными рис. 1 (несо-

ответствие экспериментальных и теоретических зависимостей ст0,05(/) около 5 %).

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

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

ЛИТЕРАТУРА

1. Ильинский А.И. Структура и прочность слоистых и дисперсно-упрочненных пленок. М., 1986. 142 с.

2. Ильинский А.И., Фукс М.Я., Аринкин А.В. и др. // ФММ. 1985. Т. 60. Вып. 5. С. 58.

3. Доценко В.И., Ландау А.И., Пустовалов В.В. Современные проблемы низкотемпературной пластичности материалов. Киев, 1987. 162 с.

4. Борздыка А.М., Гецов Л.Б. Релаксация напряжений в металлах и сплавах. М., 1978. 256 с.

5. Хирш П.Б., Хэмпфри Ф.Дж. // Физика прочности и пластичности / Под ред. Л.К. Гордиенко. М., 1972. С. 158.

6. Малинин Н.Н. Ползучесть и обработка металлов. М., 1986. 210 с.

7. Пуарье Ж. П. Высокотемпературная пластичность кристаллических тел. М., 1982. 272 с.

8. Терлецкий А.С., Белова С.В. // Металлофизика. 1992. Т. 14. № 6. С. 68.

УДК 539.124

МОДЕЛИРОВАНИЕ ОХЛАЖДЕНИЯ АМОРФНОГО СПЛАВА

© И.Б. Волкова, М.А. Баранов, В.Я. Баянкин

Россия, Ижевск, Физико-технический институт УрОРАН

Volkova I.B., Baranov M.A., Bayankin V.Y. Simulation of amorphous alloy cooling. The molecular dynamics technique is used to investigate the cooling process of the amorphous alloy Fe8oB20. The computer simulation results are compared to the real experimental data. A mechanism elucidating the resulting regularities is proposed.

Широко изучены процессы структурной релаксации и кристаллизации аморфных сплавов при средних температурах и температурах вблизи температуры кристаллизации. В то же время, влияние низкотемпературных обработок на структурные превращения и сегрегационные эффекты в поверхностных слоях аморфных систем практически не исследовано, поскольку, на первый взгляд, нет явных причин для каких-либо структурных превращений при охлаждении. Между тем данный вопрос, видимо, требует более тщательного внимания. Так, например, в работе [1] приведен це-

лый ряд экспериментальных эффектов, возникающих при низкотемпературном воздействии на аморфную систему (77-133 К).

Данная работа посвящена исследованию методом молекулярной динамики некоторых экспериментальных данных, приведенных в работе [1], в которой исследовалось влияние выдержки при температуре жидкого азота магнитомягкого аморфного сплава Со57№10Ее581пБ17. Время выдержки составляло 10, 60 и 120 мин. При этом было выявлено, что низкотемпературное воздействие не приводит к заметным изменени-

ям структуры сплава, образцы остались рентгеноаморфными. В то же время, низкотемпературная обработка вызвала концентрационные изменения в поверх- &

ностных слоях сплава. Оценку поверхностного пере- ю

распределения компонентов аморфных сплавов произ- 1

водили методом оже-спектроскопии после очистки го

поверхности («0,5 нм), когда загрязнение и адсорбиро- ^

ванные примеси на поверхности практически отсутствуют. Показано, что атомы металлоидов обогащают о

поверхность, а атомы Со и N1 обедняют. В работе [1] к

также указывалось, что подобное поведение компонентов наблюдалось при охлаждении аморфной системы Ее-№-Р-Б. §

Для выявления механизма данного эффекта использовалось молекулярно-динамическое моделирование. Методом молекулярной динамики исследовался кристаллит сферической формы, состоящий из 769 атомов железа и 193 атомов бора, что соответствует составу Ее80Б20. Температурный режим в компьютерном эксперименте был выбран таким образом, чтобы на поверхности кристаллита поддерживался постоянный градиент температур, равный 100 К/нм. Причем на поверхности образца поддерживалась постоянная температура 100 К, а в его глубине - 300 К. Глубина расчета составляла 60000 шагов, при временном шаге 2-10-15 с. При таких временах невозможно полностью проследить за диффузионными процессами, но есть надежда, что удастся уловить некоторые тенденции.

Результатом компьютерного эксперимента явилось увеличение количества атомов бора на поверхности кристаллита. При этом число атомов железа на поверхности практически не менялось. Однако относительная концентрация и того, и другого компонента, естественно, изменилась соответствующим образом. На рис. 1 демонстрируется расчет относительной концентрации для 7 атомных слоев, причем отсчет ведется из глубины образца, поверхностный слой имеет порядковый номер 6. Из рисунка 1 следует, что обогащение поверхностных слоев идет, прежде всего, за счет более глубоких слоев (так, увеличение атомов в 6-м слое привело к уменьшению в 4-м). Аналогично, число атомов в 3-м слое увеличивается за счет атомов 1-го. В общем, возникает тенденция к увеличению относительной концентрации атомов бора на поверхности.

Кроме того, в ходе компьютерного эксперимента было замечено, что миграция бора происходит в составе железо-борных кластеров. Об этом свидетельствует и анализ функции радиального распределения атомов. Подобные функции были получены при исследовании других воздействий на аналогичную систему и подробно изучены в работах [2, 3].

Анализируя результаты вычислительного эксперимента, можно сделать вывод о том, что, по крайней мере, на первом этапе эксперимента, связанного с охлаждением аморфной системы, перераспределение компонент происходит путем миграции атомов бора в составе железо-борных кластеров. Подобная миграция может быть описана в рамках теории движения макроскопических включений, разработанной Гегузиным и Кривоглазом и подробно изложенной в монографии [4]. Ранее в работах [2, 3] элементы этой теории применялись для описания перераспределения компонентов аморфизующихся систем в условиях деформации изгибом и ионного облучения. Расчетные данные с использованием данной теории хорошо коррелировали с экспериментом. По всей вероятности, достижения

слои

Рис. 1. Зависимость концентрации бора по глубине (6 атомный слой - поверхность) для различных времен охлаждения (1 шаг -2-10-14 с)

О

0.08---------1------1------1------1-------1------1

0 2 4 6

номер слоя

Рис. 2. Зависимость концентрации бора по глубине (6 атомный слой - поверхность) для различных времен нагрева (1 шаг -2-10-14 с)

указанной теории возможно применить и для описания поведения компонент аморфного сплава при охлаждении. Так, следуя работе [4], включение будет двигаться в поле градиента температур относительно поверхности со скоростью:

V = 2 С а Г1 В , (1)

где

С = (А - Х>У(2А + ХО, (2)

А и А0, в свою очередь, - теплопроводности матрицы и включения соответственно, а - безразмерное термодиффузионное соотношение (порядок которого составляет от 1 до 10), Т - температура образца в месте нахождения включения, В - коэффициент самодиффузии атомов матрицы, УТ„ - градиент температуры вдали от включения. Для определения порядка и знака ^ можно сделать обоснованное рядом экспериментальных дан-

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

Vпов = -а Т-1Б УТ„ . (3)

Кстати, при построении подобных моделей часто не учитывается этот факт. Если принять вышеприведенные положения, то получается, что металл-металлоидные кластеры будут двигаться в направлении поверхности со скоростью, в три раза превышающей скорость движения самой поверхности. В этом случае будет наблюдаться эффект, сходный с эффектом, описанным в работе [1]. Попробуем оценить величину коэффициента диффузии, необходимую для того, чтобы было возможно объяснить изучаемый эффект в рамках описанной модели. Если предположить, что V ~ 10~10 м/с, С = -1, УТ ~ 1012 К/м, Т ~ 300 К, а ~ 10, то В ~ 10-21, что хорошо совпадает с данными работы [5]. При этом следует понимать, что основное перераспределение компонент будет происходить в области температуры 300 К. Если предложенный подход верен, то в случае разогрева поверхности аморфного сплава должно наблюдаться обратное перераспределение компонент. Подобное исследование было проделано в рамках метода молекулярной динамики. Параметры компьютерного эксперимента совпадали с параметра-

ми, указанными выше, но направление градиента температур было изменено на противоположное (результаты приведены на рис. 2). Видно хорошее качественное согласие с предложенной моделью.

ЗАКЛЮЧЕНИЕ

Таким образом, методом молекулярной динамики исследовано поведение аморфной системы Ее80Б20 при температурных воздействиях в области низких и высоких температур. Результаты компьютерного эксперимента находятся в качественном согласии с результатами реального эксперимента. Анализ полученных результатов позволил предложить механизм перераспределения компонентов (прежде всего, атомов металлоидов) в верхних слоях аморфной системы при температурном воздействии. В основе предложенной модели лежит теория движения макроскопических включений, подробно разработанная Гегузиным и Кривоглазом.

ЛИТЕРАТУРА

1. Баянкин В.Я., Васильев В.Ю., Волкова И.Б. // Металлы. 1997. № 1. С. 164-167.

2. Волкова И.Б., Баранов М.А., Баянкин В.Я. // Материаловедение. 1998. № 6. С. 2-8.

3. Баранов М.А., Баянкин В.Я. // Материаловедение. 1999. № 4. С. 18-25.

4. Гегузин Я.Е., Кривоглаз М.А. Движение макроскопических включений в твердых телах. М.: Металлургия, 1971. 344 с.

5. Бокштейн Б. С., Клингер Л.М., Уварова Е.Н. Аморфные металлические сплавы: Науч. тр. № 147. МИСиС / Под ред. Ю.А. Скакова. М.: Металлургия, 1983. С. 81-86.

УДК 669.71З:669

ОСОБЕННОСТИ ИЗМЕНЕНИЙ СТРУКТУРЫ И МЕХАНИЧЕСКИХ СВОЙСТВ ПРИ ОТЖИГЕ СУБМИКРОКРИСТАЛЛИЧЕСКОГО СПЛАВА А1 - 3 % М§

© Н.К. Ценев, А.М. Шаммазов, В.Д. Суханов, Б.А. Козачук

Россия, Уфа, Государственный нефтяной технический университет

Tsenev N.K., Shammazov A.M., Suhanov V.D., Kozachuk B.A. Features of structure changes and mechanical properties during the annealing of submicrocrystalline of alloy Al - 3 % Mg. Using the pattern Al-based alloy (Al - 3 % Mg) processed by severe plastic deformation techniques, the article looks at the TEM studies of the specific submicrocrystalline (SMC) structure. The annealing of the alloy results in considerable changes in microstructure: grain boundary structure, grain size. It is revealed that abnormal grain growth leads to a sharp increase of lattice dislocation density enhancing the microhardness. Two possible reasons of the unusual behaviour in the SMC alloy during annealing are considered.

В последние годы субмикрокристаллические (СМК) материалы ( < 1мкм), получаемые методами интенсивной пластической деформации, привлекают пристальное внимание специалистов в области физики прочности и пластичности [1-3]. Этот интерес обусловлен необычными механическими свойствами, проявляемыми СМК материалами: при комнатной температуре наблюдается одновременное повышение значений прочности и пластичности, обнаруживается эффект низкотемпературной высокоскоростной сверхпла-

стичности (НВССП) [4-6] и др. Очевидно, что столь необычное поведение СМК сплавов обусловлено протяженностью и состоянием границ зерен, их термостабильностью. Однако структурные изменения, происходящие в процессе роста зерен в таких материалах, практически не изучены. В этой связи данная работа посвящена экспериментальным наблюдениям изменений дислокационной структуры, происходящих при отжиге на примере субмикрокристаллического модельного сплава А1 - 3 % М&

i Надоели баннеры? Вы всегда можете отключить рекламу.