CC BY
37
УДК 621.921.34:621.7.044.2
СТАТИСТИЧЕСКОЕ МОДЕЛИРОВАНИЕ КОРРЕЛЯЦИИ ЛОКАЛЬНОЙ ПЛОТНОСТИ МАКРОСТРУКТУРЫ ПРИ ОДНОСТОРОННЕМ ПРЕССОВАНИИ ПОРОШКОВЫХ
МАТЕРИАЛОВ
© 2011 Г Г. Винокуров1, О Н. Попов2
1 Институт физико-технических проблем Севера им. В.П. Ларионова СО РАН,
г. Якутск
2 НИИ математики, Институт математики и информатики Северо-Восточного
федерального университета им. М.К. Амосова, г. Якутск
Поступила в редакцию 23.01.2011
Методом статистического моделирования изучена макроструктура порошковых материалов, полученных односторонним прессованием. Использована вероятностно-геометрическая модель макроструктуры порошковой среды с программной реализацией на языках Paskal и Delphi. Проведены расчеты корреляции локальной плотности прессованных порошковых материалов.
Ключевые слова: порошковые материалы, прессование, макроструктура, статистическое моделирование
В настоящее время существуют многочисленные работы по математическому моделированию процесса прессования порошков [1-6]. С развитием информационных технологий характеристики макроструктуры порошковой среды также могут удовлетворительно оцениваться из результатов статистического моделирования по вероятностно-геометрическим системам моно- и полидисперсных частиц [7-10]. Как показывают работы в этом направлении, изменение характеристик качественно отражает основные статистические закономерности формирования макроструктуры порошковых материалов в технологических процессах спекания, прессования и напыления. Одной из актуальных задач данных подходов является установление корреляционных характеристик локальной плотности, описывающих взаимосвязь элементов макроструктуры прессованных порошковых материалов и, безусловно, связанных с их свойствами переноса и прочности.
Цель настоящей работы - выявление корреляции локальной плотности макроструктуры порошковых материалов при одностороннем прессовании на основе статистического моделирования вероятностно-геометрической системой порошковой среды.
Винокуров Геннадий Георгиевич, кандидат технических наук, ведущий научный сотрудник Отдела материаловедения. Е-шаИ: [email protected] Попов Олег Николаевич, старший преподаватель. Е-шаП: [email protected]
Статистическое моделирование корреляционных характеристик макроструктуры порошковых материалов при одностороннем прессовании. В данной работе для моделирования корреляции локальной плотности порошковой среды использована элементарная двумерная модель Монте-Карло формирования макроструктуры порошков при одностороннем прессовании. Модель построена на основе статистического подхода, предложенного в работе [7], где в качестве количественных характеристик макроструктуры выбраны продольные, поперечные (к направлению прессования) координационные числа частиц и локальная пористость при каждом шаге прессования. Экспериментальное определение данных характеристик макроструктуры в процессе прессования представляется практически невозможным из-за большого количества частиц порошкового материала. Статистический подход моделирует коллективное поведение частиц порошкового материала, на основе которого можно качественно описать процесс уплотнения порошковых материалов при прессовании.
Начальное состояние прессуемого порошкового материала задается прямоугольной матрицей макроструктуры с размерностью (Ь*М), которая случайным образом заполняется частицами порошка. Следует отметить, что начальная пористость в двумерных моделях задается с учетом заполнения частицами сечения трехмерного порошкового тела. Состояние
макроструктуры порошкового материала на определенном шаге прессования получается из предыдущего состояния переносом частиц первого столбца по горизонтали вправо в ближайшую незаполненную ячейку; после переноса всех заполненных элементов первый столбец удаляется. Таким образом, процессу одностороннего прессования порошкового материала соответствует последовательное сокращение одного крайнего столбца матрицы. Если все ячейки в строке окажутся заполненными (возникновение перколяции), процесс прессования считается законченным [7]. В работе по данной двумерной модели Монте-Карло при каждом шаге прессования вычислялись локальная плотность столбцов матрицы и соответствующие коэффициенты корреляции между ними. Для модели разработаны программы на языках Pascal и Delphi, которые позволяют описать изменение макроструктуры в матрицах размером до 1000x1000 с усреднением по 10000 реализациям.
Актуальной задачей статистического моделирования является установление влияния состояния макроструктуры порошкового материала на данном этапе прессования на ее последующие состояния. Необходимо установить уровень этого влияния, и выяснить расстояние, на которое распространяется влияние в процессе прессования. Данные закономерности процесса прессования математически представляются поведением корреляционной функции локальной плотности столбцов x вероятностно-геометрической системы. Далее и на рисунках используется значение относительной локальной плотности р, введены дискретные переменные для шага прессования и столбцов матрицы макроструктуры; используются обозначения: pjf - локальная плотность сечения m при шаге прессования tk; г(р{, -коэффициент корреляции между локальной плотностью сечений (столбцов) l и m, соответственно, при шагах прессования tj и tk.
Рассмотрим объемные сечения, соответствующие столбцам x-1, x, x+1 вероятностно-геометрической системы макроструктуры порошкового материала, перпендикулярные направлению прессования (рис. 1). Пусть при шаге прессования ti состояния сечений x-1, x, x+1 характеризуются своими локальными плотностями px-11 рх, px+11 соответственно. Обозначим через At шаг прессования, при котором частицы материала проходят путь длиной, равный характерному размеру частицы. Тогда при шаге прессования t^=t1 + At состояние сечений x и x+1 будут описываться некоторыми значениями плотности px2, px+12 соответственно (рис. 1). Следует предположить,
что в основном на величину р2 влияют значения рх.\1 и рх\ поэтому в работе исследованы корреляционные зависимости данных случайных величин.
Коэффициент корреляции локальной плотности макроструктуры порошкового материала определялся по формуле:
r (px , Pk)=
(P'x -<P'x >)(Pk~ < Pk >)
где х, х' — порядковые номера столбцов матрицы; <Рх'>, <Рхк> - средние значения, ах', -среднеквадратические отклонения плотности столбцов матрицы при шагах прессования ^ и соответственно. Для обеспечения достоверности результатов и сглаживания колебаний коэффициента корреляции усреднение в формуле проводилось по 10000 реализациям модели Монте-Карло.
с>
1 Рх-1 1 Рх 1 Рх+1 h
х-\ X i + 1
2 рх 2 рх+1 t2
X J + 1
Н аправление прессования
Рис. 1. Схема сечений порошкового материала при шагах прессования ^ и (2.
Обсуждение результатов. Для анализа закономерностей корреляции локальной плотности целесообразным является также рассмотрение случайной упаковки частиц с малой плотностью. Такие порошковые материалы с существенно пористой макроструктурой, т.н. трехмерные капиллярно-пористые покрытия недавно разработаны для медицинских приложений [10]. На рис. 2 приведены коэффициенты корреляции г локальной плотности слоев в начальной стадии прессования, показывающие влияние величин рх-11 и рх на уровень рх2; также показано изменение относительной плотности р по сечению порошкового материала. Как видно из графиков на рис. 2,а, при высокой начальной плотности уровень коэффициента корреляции г(рх.\, рх2) даже вблизи границы прессования, где наблюдается уплотнение порошкового материала, является низким, корреляционная связь несущественна («0,4-0,2). По мере удаления от границы прессования с понижением плотности корреляция г(рх.\, рх2) вовсе стремится к нулю, так как макроструктура материала приближается к начальной равновесной упаковке. Коэффициент корреляции г(рх, рх2), наоборот, стремится к 1, вследствие уменьшения влияния соседнего сечения х-1.
№
14
16
18
30
б)
Рис. 2. Изменение плотности и корреляции по сечению порошкового материала в начальной стадии прессования: а) начальная плотность 0,75, 1\ =1, ¿2=2; б) начальная плотность 0,25, 1\ =12, ¿2=13; 1-р; 2 - г(рхЛ\ рх2); 3 - г(рх\ рх2); матрица размера 300x1200
Как видно из графиков на рис. 2б при низкой начальной плотности характер изменения коэффициента корреляции г(рх-1 , рх2) соответствует поведению локальной плотности вблизи границы прессования, где наблюдается уплотнения порошкового материала; причём корреляционная связь достаточно существенна («0,8-0,6). По мере удаления от границы прессования с понижением плотности корреляция г(рх.\, рх2) стремится к нулю, так как макроструктура материала также приближается к начальной равновесной упаковке. Коэффициент корреляции г(рх, рх2), наоборот, стремится к 1, также вследствие снижения влияния соседнего сечения х-1. Такие же закономерности, но более с четким проявлением можно наблюдать при больших значениях степени прессования порошковых материалов. На рис. 3 приведены коэффициенты корреляции локальной плотности слоев при шагах прессования 100-
101, показывающие влияние величин рх-\ и рх на уровень рх2; также показано изменение относительной плотности р по сечению порошкового материала.
127
137
147
б)
Рис. 3. Изменение плотности и корреляции по сечению порошкового материала при больших степенях прессования: а) начальная плотность 0,75, ¿1 =100, ¿2=101; б) начальная плотность 0,25, ¿1 =100, ¿2=101; 1-р; 2 - г(рхЛ\ рх2); 3 - г(рх, рх2); матрица размера 300x1200
Как видно из графиков на рис 3а при высокой начальной плотности вблизи границы прессования наблюдается несущественная корреляция г(рх-\, рх2). Отсюда можно сделать вывод, что при прессовании реальных порошковых материалов, плотность которых всегда не меньше плотности исходной засыпки (стадия нестабильных и стабильных пространственных структур [4] с плотностью не более «0,7-0,8), значение локальной плотности сечения, в основном, определяется только состоянием этого же сечения на предыдущем шаге
прессования. Как видно из графиков на рис. 3б, для низкой начальной плотности характер изменения коэффициента корреляции г(рхл , рх2) также соответствует поведению локальной плотности вблизи границы прессования; корреляционная связь достаточно существенна («0,8-0,6). По мере удаления от границы прессования с понижением плотности корреляция г(Рхл, Рх) стремится к нулю; в то же время коэффициент корреляции г(рх, рх) также возрастает до 1.
Экспериментальными исследованиями авторов ранее установлено, что существует корреляция локальной плотности различных слоев газотермического покрытия, которая исчезает на расстоянии порядка диаметра напыляемых частиц [11]. Актуальной становится задача выявления аналогичной корреляции между локальными плотностями слоёв прессованного порошкового материала, перпендикулярных направлению прессования. Для этого следует установить корреляцию локальной плотности соседних и близких слоев. Более четкие закономерности наблюдаются для случайной упаковки частиц с малой плотностью вблизи границы прессования; на рис. 4 показана корреляция г плотности столбца х матрицы, полученной на сотом шаге прессования, с плотностью столбца х+1 на /-м шаге прессования вблизи границы прессования, где 1=101,102,103,104,105. Как видно из рисунка, влияние состояния сечения х матрицы на соседнее сечение практически исчезает после дальнейших пяти шагов прессования - падает с «0,8 до «0,2. Вдали от границы прессования корреляция несущественна с самого начала. В случае высокой начальной плотности (0,75) такая корреляция практически отсутствует уже на первом шаге (рис. 3,а график 2).
. ше
Рис. 4. Коэффициент корреляции г(рх.\100, рх') локальной плотности соседних сечений порошкового материала вблизи границы прессования; матрица размера 300^1200; начальная плотность 0,25; 1=101 - 105
Выводы:
1. В начальной стадии прессования при высокой начальной плотности порошковых материалов практически отсутствует влияние локальной плотности столбца матрицы макроструктуры на локальную плотность соседнего столбца при последующих шагах прессования. Это влияние существенно только для низкой начальной плотности вблизи границы прессования, коэффициент корреляции ~ 0,8-0,6.
2. Влияние текущего состояния столбца на его локальную плотность при следующем шаге прессования является существенным независимо от начальной плотности, коэффициент корреляции близок к 1.
Таким образом, локальная плотность столбца матрицы макроструктуры прессованных порошковых материалов практически определяется его состоянием на предыдущем шаге прессования.
СПИСОК ЛИТЕРАТУРЫ:
1. Друянов, Б.А. Математическое моделирование процессов обработки давлением металлических порошков и пористых тел / Б.А. Друянов, И.Д. Ра-домысельский, М.Б. Штерн // Порошковая металлургия. 1981. № 3. С. 6-11.
2. Новофастовский, А.Л. Математическое моделирование процесса прессования порошков методом крупных частиц // Порошковая металлургия. 1989. № 12. С. 4-7.
3. Феноменологические теории прессования порошков / М.Б. Штерн и др. - Киев.: Наук. думка, 1982. 140 с.
4. Радомысельский, И.Д. Некоторые особенности уплотнения порошков на разных стадиях прессования / И.Д. Радомысельский, Н.И. Щербань // Порошковая металлургия. 1980. № 11. С. 12-19.
5. Николаев, А.Н. Связь между давлением и плотностью прессовок из металлических порошков // Порошковая металлургия. 1975. № 6. С. 32-42.
6. Рогозин, В.Д. Уравнение прессования порошков // Порошковая металлургия 1981. № 6. С. 28-31.
7. Каминский, В.М. Двумерная стохастическая модель уплотнения порошковых материалов / В.М. Каминский, А.Н. Николенко, И.Я. Сидоренко // Порошковая металлургия. 1982. №2. С. 29-31.
8. Николенко, А.Н. Анализ случайной упаковки идентичных частиц. Общая теория / А.Н. Николенко, М. С. Ковальченко // Порошковая металлургия. 1985. № 11. С. 38-41.
9. Кадушников, Р.М. Геометрическое моделирование структуры полидисперсных материалов / Р.М. Кадушников, А.Р. Бекетов // Порошковая металлургия. 1989. № 10. С. 69-74.
10. Гнедовец, А.Г. Модель формирования макроструктуры покрытий при плазменном напылении / А . Г. Гнедовец, В. И. Калита // Физика и химия обработки материалов. 2007. № 1. С. 30-39.
11. Винокуров, Г. Г. Статистические подходы для описания формирования и изнашивания макроструктуры порошковых покрытий и материалов, полученных высокоэнергетическими методами / Г.Г. Винокуров, О.Н. Попов. - М.: Academia, 2009. 184 с.
STATISTICAL MODELING OF LOCAL DENSITY MACROSTRUCTURE CORRELATION AT ONE-SIDED PRESSING OF POWDER MATERIALS
© 2011 G.G. Vinokurov1, O.N. Popov2
1 Institute of Physicotechnical Problems of the North named after V.P. Larionov SB RAS,
Yakutsk
2 Scientific Research Institute of Mathematics, Institute of Mathematics and Computer Science of Northeast Federal University named after M.K. Amosov, Yakutsk
By method of statistical modeling it is studied macrostructure of powder materials received by one-sided pressing. The stochastic geometrical model of powder macrostructure with program realization in Paskal and Delphi languages is used. Calculations of correlation of pressed powder materials local density are lead.
Key words: powder materials, pressing, macrostructure, statistical simulation
Gennadiy Vinokurov, Candidate of Technical Sciences, Leading Fellow at the Materials Technology Department. E-mail: [email protected] Oleg Popov, Senior Teacher. E-mail: [email protected]
