МОДЕЛИРОВАНИЕ ЭВОЛЮЦИИ ЗЕРЕННОЙ СТРУКТУРЫ В ПРОЦЕССЕ КРИСТАЛЛИЗАЦИИ МЕТОДОМ КЛЕТОЧНЫХ АВТОМАТОВ

Механические свойства металлических изделий в значительной степени определяются зеренной структурой, ключевыми характеристиками которой являются размер, форма и кристаллографическая ориентация кристаллических решеток зерен. Особое внимание при изготовлении ответственных деталей необходимо уделять росту столбчатых зерен, наличие которых характерно для аддитивного производства и литья. Их рост приводит к анизотропии, а следовательно, к снижению прочностных свойств изделия. Для прогнозирования и целенаправленного управления этими характеристиками необходимо глубокое понимание процессов кристаллизации, в ходе которой формируется зеренная структура. Целью исследования является разработка и верификация математической модели кристаллизации на основе метода клеточных автоматов, позволяющей прогнозировать эволюцию зеренной структуры с учетом стохастического зарождения и анизотропного роста. Разработанная модель включает стохастические алгоритмы гетерогенного и гомогенного зарождения зерен на поверхности и в объеме расплава. Для сохранения первоначальной кристаллографической ориентации растущего зерна был применен метод децентрированного квадрата. Проведена серия численных экспериментов по затвердеванию сплава Al7Si. Исследовано влияние скорости охлаждения, критической температуры зарождения и максимальной плотности зародышей на формирование зеренной структуры. Результаты показали, что модель корректно воспроизводит ключевые физические механизмы: рост столбчатых зерен от границ с последующим переходом к равноосной структуре в центре, а также конкурентный рост зерен с предпочтительной ориентацией. Установлено, что увеличение скорости охлаждения и плотности зародышей приводит к уменьшению среднего размера зерна, а рост критической температуры зарождения способствует расширению зоны столбчатых зерен. Реализованная и верифицированная двумерная модель процесса кристаллизации представляет собой начальный этап разрабатываемого программного комплекса, объединяющего метод клеточных автоматов и метод гидродинамики сглаженных частиц для моделирования процессов проволочной наплавки металлов.

Modeling the grain structure evolution during crystallization by cellular automata

The mechanical properties of metallic components are largely governed by the grain structure, whose key characteristics are grain size, morphology, and the crystallographic orientation of grain lattices. When manufacturing critical parts, particular attention must be paid to the growth of columnar grains, which are characteristic of additive manufacturing and casting processes. Their growth induces anisotropy and, consequently, reduces strength. Predictive and targeted control of these characteristics requires a deep understanding of the solidification processes by which the grain structure forms. The aim of this study is to develop and verify a cellular-automaton-based mathematical model of solidification that predicts the evolution of grain structure while accounting for stochastic nucleation and anisotropic growth. The model incorporates stochastic algorithms for heterogeneous and homogeneous nucleation on the surface and in the bulk of the melt. To preserve the initial crystallographic orientation of a growing grain, we employ the decentered square algorithm. A series of numerical experiments was performed for the solidification of an Al7Si alloy. We investigated the influence of cooling rate, critical nucleation temperature, and maximum nucleation density on the resulting grain structure. The results demonstrate that the model accurately reproduces key physical mechanisms, including the growth of columnar grains from boundaries, followed by a transition to an equiaxed structure in the center (columnar-to-equiaxed transition, CET), as well as the competitive growth of grains with preferential orientation. It is established that increasing the cooling rate and nucleation density decreases the mean grain size, whereas increasing the critical nucleation temperature broadens the columnar zone. The developed and validated two-dimensional solidification model represents the initial stage of the developing software framework, combining the cellular automaton method with the smoothed particle hydrodynamics method for modeling metal wire-arc additive manufacturing processes.