Multilevel model of polycrystals: application to assessing the effect of texture and grains misorientation on the critical deformation of the dynamic recrystallization initiation
- Authors: Kondratev N.S1, Trusov P.V1, Podsedertsev A.N1
- Affiliations:
- Perm National Research Polytechnic University
- Issue: No 4 (2021)
- Pages: 83-97
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/1955
- DOI: https://doi.org/10.15593/perm.mech/2021.4.09
- Cite item
Abstract
The current state of materials constitutive models and the computer technology development make it possible to numerically implement complex multilevel models that allow describing the material structure evolution. In this regard, it is possible to formulate optimal control problem for metal forming processes in order to create the required performance characteristics of finished products and their ingots. To solve this problem in this study, the effective multilevel modeling approach is used to describe the thermomechanical treatment of polycrystalline materials. The model is based on this approach with the introduction of internal variables, in which the carriers and physical mechanisms of the processes of hot intense plastic deformation are explicitly considered. At deformation temperatures order of 0.5 homologous and above, recrystallization process have a special effect on the formation and change of the grain and defect material structure. The paper considers the problem of determining the critical deformation of dynamic recrystallization initiation, that depending on the material texture and the mutual misorientation of neighboring grains. Numerical experiments of the multilevel model are used to simulate two stages of inelastic deformation for this purpose. At the first stage, cold inelastic deformation by simple shear and compression is considered, that leading to the formation of a corresponding texture. At the second stage, uniaxial hot tension deformation is considered. The initial distribution of crystallographic grain orientation is assumed to be uniform. Two variants of the grains mutual misorientation with the prescribed increased and decreased values of the average misorientation angles are considered. The recrystallization process is not explicitly modeled. The current model is intended to assess the recrystallization critical deformation. It is shown that the mutual misorientation of grains, rather than texture, has the most influence on the critical deformation. An increase in the angle of grains mutual misorientation contributes to an earlier start of the dynamic recrystallization process. The formation of a deformation texture leads to a decrease in the angle of mutual misorientation, and, accordingly, to a decrease in dynamic recrystallization intensity. Despite this, with an increase of deformation, the driving force of recrystallization (the average value of the difference of stored energy between neighbor grains) is increases, which leads to the implementation of dynamic recrystallization.
Full Text
Введение Применение традиционных методов обработки поликристаллических материалов давлением, таких как прокатка, прессование, волочение, ковка, экструзия и др., остается актуальным на сегодняшний день. Интенсивная пластическая деформация приводит к значительному увеличению плотности дислокаций в объеме и на границах зерен поликристаллов [1-4]. Это, в свою очередь, является причиной перехода рассматриваемой системы (поликристалла) в термодинамически нестабильное состояние [5-6], увеличения напряжений, необходимых для продолжения деформации, и возможного разрушения [7-8]. Как правило, отмеченные выше процессы обработки металлов реализуются при повышенных гомологических температурах [9-12] или являются многостадийными, состоящими из последовательности операций формования с промежуточными стадиями отжига [6, 13]. Обработка металлов при повышенных температурах сопровождается соответствующими процессами на мезо- и микроуровне, которые стремятся снизить плотность дефектов и уменьшить свободную энергию материала [6, 14-15]. Основными процессами изменения микроструктуры при горячем деформировании, ведущими в том числе к разупрочнению, являются динамический возврат и динамическая рекристаллизации [6, 14-16]. Динамический возврат проходит при относительно низких температурах, вплоть до комнатной, и отвечает за процессы формирования малоугловых границ и их эволюции. В процессе рекристаллизации происходит формирование и эволюция высокоугловых границ, их миграция [15, 17]. Наиболее значимыми с точки зрения изменения структуры и механических макросвойств являются процессы динамической рекристаллизации [6, 13-14, 16]. Основная движущая сила первичной рекристаллизации - разность запасенной энергии между соседними зернами [6, 13-17]. На начальной стадии пластического деформирования основным фактором, определяющим запасенную на дефектах энергию в зерне, является увеличение плотности дислокаций по активным системам и их взаимодействие [18-19]. Активация систем скольжения зависит от ориентации кристалла относительно осей обработки материала. Взаимодействие дислокаций определяется как внешними параметрами воздействий (температурой, скоростью деформирования и др.), так и характеристиками материала (типом решетки, энергией дефекта упаковки) и начальным состоянием его структуры [6, 18]. Поскольку основной причиной первичной рекристаллизации является разность накопленных энергий соседствующих кристаллитов, то помимо распределения ориентаций зерен другим существенным фактором является взаимная ориентация соседних зерен. Прямая зависимость между текущим состоянием структуры и макрохарактеристиками поликристаллов подтверждена множеством экспериментальных данных [4, 6, 13, 20]. Хорошо известно, что в результате интенсивной горячей пластической деформации происходит существенное изменение (обычно измельчение) зеренной структуры, которое во многом определяется процессом динамической рекристаллизации [6, 14]. Для многих эффективных механических свойств поликристаллического материала (упругих, пластических, ползучести, твердости, усталости и ряда других) определяющее значение имеет зеренная структура (форма и размеры зерен) и текстура материала [21-25]. Формирование равноосных зерен и их измельчение в процессе рекристаллизации способствует реализации межзеренного скольжения и переходу материала в режим структурной сверхпластичности [26-28], которое обеспечивает достижение аномально больших деформаций при почти постоянных значениях напряжений. Отмеченное выше подчеркивает актуальность исследования динамической рекристаллизации в качестве как отдельного, так и аккомодационного процесса (например, при реализации деформирования в режиме структурной сверхпластичности). Следовательно, актуальной проблемой является формулировка конститутивной модели поликристаллического материала, в основу которой положено явное физическое описание процессов различных структурно-масштабных уровней (включая рекристаллизацию). Для решений этой задачи эффективным инструментом являются многоуровневые модели с внутренними переменными [19]. Многоуровневые модели можно разделить на два основных типа - прямые и статистические [29-30]. Для прямых моделей характерным является пространственное рассмотрение полевых величин (напряжений, деформаций и других внутренних переменных) для каждого кристаллита (зерна, субзерна, фрагмента), входящего в рассматриваемый представительный макрообъем, с учетом его ориентации, формы и границ. Статистические модели рассматривают представительный макрообъем материала («макроточку») как совокупность отдельных зерен, объединяемых в поликристаллический агрегат гипотезой связи (Фойгта, Рейса, Крёнера). Прямые модели обладают большей точностью, при этом являются более ресурсоемкими по сравнению со статистическими. В данной работе в рамках многоуровневого подхода рассматривается расширенная статистическая модель неупругого деформирования, где явно учитывается состояние соседних зерен посредством учета их контакта по сопрягающимся границам. Целью работы является исследование влияния формируемой текстуры в процессе пластического деформирования и взаимной разориентации соседних зерен на величину критической деформации начала динамической рекристаллизации с использованием статистической модели неупругого деформирования. 1. Прерывистая динамическая рекристаллизация В процессе пластической деформации при повышенных температурах формируются практически бездефектные (малодефектные) зерна, которые в дальнейшем увеличивают свой объем за счет окружающих более дефектных зерен [6, 13-17]. В результате этого процесса, носящего название динамической рекристаллизации, образуется новая зеренная структура с низкой плотностью дефектов (дислокаций) при полной или частичной замене исходных зерен. В случае предварительной холодной деформации и последующего нагрева при отсутствии силовых воздействий также происходит процесс изменения зеренной и дефектной структуры материала, названный статической рекристаллизацией [6, 13-17]. Между процессами статической и динамической рекристаллизации есть фундаментальное отличие, которое заключается в том, что при статической рекристаллизации после деформационного упрочнения происходит резкий переход из метастабильного состояния в относительно стабильное состояние, соответствующее температуре проведенного отжига [5-6]. Такой переход аналогичен фазовому переходу. При динамической рекристаллизации и высоких степенях деформации происходит переход к установившемуся динамическому равновесию микроструктуры при сохраняющемся среднем уровне свободной энергии системы. Динамическую рекристаллизацию обычно разделяют на два типа - прерывистая и непрерывная динамическая рекристаллизация [6, 13-14, 16-17]. Некоторые исследователи выделяют в отдельный тип геометрическую рекристаллизацию [15]. В материалах с низкой энергией дефекта упаковки, где затруднены процессы возврата, происходит образование зародышей гетерогенным путем. В этом случае могут быть четко разделены стадии образования зародышей рекристаллизации и их последующий рост. Эта рекристаллизация носит название прерывистой или классической [5-6]. В материалах с высокой энергией дефекта происходит образование ячеистой структуры (с формированием соответствующих малоугловых границ) в результате активного действия процесса динамического возврата. В дальнейшем малоугловые границы переходят в высокоугловые эволюционным путем. В этом случае процесс рекристаллизации происходит гомогенно (однородно) таким образом, что явно не выделяются стадии образования и роста новых зерен, в силу чего для такой рекристаллизации применяется термин непрерывная. При этом не всегда удается четко провести раздел между данными типами прерывистой или непрерывной рекристаллизации [5, 15]. Оба типа рекристаллизации имеют сходные черты и могут происходить одновременно, дополняя друг друга, или последовательно. Более того, в зависимости от параметров воздействия для рассматриваемого материала может наблюдаться переход между различными типами рекристаллизации [31]. Для выявления реализуемого типа рекристаллизации в экспериментальных исследованиях применяется методика дифракции отраженных электронов (англ. сокращ. EBSD) [6, 17]. Предметом исследования в данной работе является прерывистая динамическая рекристаллизация. Данный тип рекристаллизации является наиболее характерным для широкого класса поликристаллических материалов: меди, никеля, аустенитных сталей, некоторых сплавов алюминия и др. [5-6, 15]. В этом случае на ранней стадии горячей пластической деформации образуются бездефектные зародыши, которые поглощают более дефектные соседние зерна. Такие зародыши рекристаллизации не образуются в результате случайных флуктуаций, как в случае фазовых переходов, и имеют характерный размер порядка микрометров [6]. Основной и широко принятый для большинства металлов механизм образования зародышей для прерывистой динамической рекристаллизации основан на движении («выгибании», «вздутии») участков границ зерен, ассоциированных с субзернами, которые исходно существовали в поликристалле [32-35]. Это явление впервые было исследовано Беком и Сперри [36]. Критерий образования новых рекристаллизованных зерен сформулирован Бейли и Хиршем [37]. Согласно этому критерию уменьшение локальной объемной энергии за счет устранения дефектов должно быть больше увеличения межзеренной энергии в результате увеличения площади границы при ее перемещении. В этом случае преимуществом к росту будут обладать зародыши рекристаллизации (субзерна) большего размера. При высокотемпературном неупругом деформировании прерывистая динамическая рекристаллизация отражается на диаграмме нагружения [5-6, 14, 17]. Ее характерный вид с одиночными или множественными пиками показан на рис. 1. Характерными параметрами диаграммы являются: 1) критическая деформация εс, до достижения которой динамическая рекристаллизация не инициируется; 2) максимальное значение напряжений σmax и деформация εmax (соответствует локальному, неустойчивому динамическому равновесию между процессами упрочнения и разупрочнения); 3) стационарное значение напряжений σs и деформация εs (соответствует началу стационарного участка на диаграмме нагружения). Основными параметрами воздействий, определяющих эволюцию структуры, являются температура испытания θ и скорость деформирования d. Как правило, для большинства материалов характерно, что при повышении температуры и уменьшении скорости деформации число пиков на диаграмме нагружения увеличивается (рис. 1, а) [5-6, 14, 17]. Основным параметром структуры, влияющим на вид диаграммы, является начальный средний размер зерен d0. Значения d0 выше критического (при фиксированных параметрах θ и d из диапазона реализации рекристаллизации) будут приводить к одиночному пику на диаграмме, ниже - к множественному (рис. 1, б) [38-39]. Рис. 1. Схема влияния температуры θ, скорости деформирования d (а) и начального среднего размера зерен d0 (б) на характерный вид диаграммы нагружения при реализации прерывистой динамической рекристаллизации Fig. 1. Scheme of the influence of temperature θ, strain rate d (a) and initial average grain size d0 (б) on the loading diagram characteristic form in discontinuous dynamic recrystallization В литературе широко освещена проблема влияния θ, d и d0 на напряженно-деформированное состояние или вид диаграммы нагружения [6, 14]. С другой стороны, более детального обсуждения требует проблема определения критической величины интенсивности деформаций εс и влияния на нее структурно значимых параметров, описывающих состояние поликристаллического материала. Проведенный предварительный анализ имеющихся данных показывает, что на εс основное влияние должны оказывать углы взаимной разориентации между соседними зернами и текстура поликристалла. В связи с этим с применением многоуровневого моделирования неупругого деформирования предлагается более детальное рассмотрение этой проблемы. 2. Подходы и методы моделирования динамической рекристаллизации, постановка численного эксперимента Для моделирования динамической рекристаллизации существует множество подходов и методов. Наиболее упрощенными для исследования рекристаллизации являются модели макрофеноменологического типа. В таких моделях рассматривается эмпирическое описание доли рекристаллизованного материала, в основу которого положено кинетические соотношения типа Джонсона - Мела - Аврами - Колмогорова (англ. - Johnson - Mehl - Avrami - Kolmogorov, сокращ. JMAK) [40-41]. Для моделирования напряженно-деформированного состояния применяются зависимости (как правило, степенного типа) напряжений (или скорости напряжений) от параметров воздействия (температуры, скорости деформации). В феноменологических соотношениях часто применяется параметр Зинера - Холломона Z, включающий в себя скорость деформирования и температуру [40, 42]. В моделях данного класса в зависимости от сложности нагружения может потребоваться введение в рассмотрение нескольких экспериментально установленных функций, а при моделировании многоэтапных процессов также могут потребоваться данные об изменении значений параметров модели. Развитием подхода JMAK являются модели, которые учитывают дополнительных характеристики зеренной структуры - размеры и форму зерен, положение межзеренной границы [43-44]. Существуют модификации JMAK-моделей, основанные на вероятностном подходе, как правило, с использованием метода Монте-Карло [45-46]. Другим наиболее развитым является физический подход, включающий явное рассмотрение эволюции структуры материала [6, 19]. В этом случае применяются физически обоснованные модели прерывистой динамической рекристаллизации. В рамках данного класса моделей вводятся переменные и параметры, характеризующие структуру материала, причины образования зародышей рекристаллизации и миграции межзеренных границ. Вследствие этого такие модели могут быть использованы для описания формирования требуемой структуры материала и являются актуальными для металлургической и обрабатывающей отраслей. В рамках физического подхода существуют несколько типов моделей. Созданы самосогласованные [47-48] и прямые [49-50] физические модели, включающие описание динамической рекристаллизации. Прямые модели, для численной реализации ориентированные на использование метода конечных элементов (МКЭ), можно разделить на два типа. В моделях первого типа рассматриваются в явном виде форма и размеры зеренной структуры, в силу ограниченности вычислительных ресурсов они не применимы для описания конструкций и деталей, но для анализа эволюции внутренней структуры являются весьма эффективным инструментом. При описании рекристаллизации в моделях первого типа одной из основных проблем является описание зеренной структуры, которая существенно эволюционирует в ходе горячей пластической деформации. Для этой цели используются несколько подходов: 1) метод Монте - Карло; 2) метод клеточных автоматов; 3) метод фазового поля; 4) метод заданного уровня (англ. - level set method); 5) метод, основанный на теории графов [6]. Прямые модели второго типа, при вычислении отклика в точках интегрирования использующие статистические подмодели, применяются для определения отклика материала на макроуровне (уровне реальных конструкций), при этом учитывая в явном виде структуру материала. В рамках настоящего исследования рассматривается расширенная физическая статистическая модель для анализа неупругого поведения представительного объема поликристалла. Как сказано выше, исследуется неупругое деформирование с позиций анализа начальной стадии образования зародышей рекристаллизации, влияния на нее взаимной ориентации соседних зерен и текстуры. Проблема явного моделирования рекристаллизации является актуальной, но выходит за рамки рассматриваемого исследования. Полную постановку задачи описания деформирования представительного макрообъема поликристаллического металла на основе двухуровневой статистической модели можно найти в [19]. Математическая формулировка модели для отдельного зерна, рассматриваемого как элемент мезоуровня, содержит следующие соотношения: (1) , (2) , (3) , (4) , (5) , (6) (7) , (8) , (9) , (10) . (11) Соответствующие начальные условия имеют следующий вид: (12) Связи переменных мезоуровня с переменными макроуровня (передача сверху вниз) осуществляется с применением расширенной гипотезы Фойгта: . (13) Выше использованы следующее обозначения: - взвешенный тензор напряжений Кирхгоффа мезоуровня, - его коротационная производная, σ - тензор напряжений Коши мезоуровня; , - транспонированный градиент скорости перемещений мезоуровня и макроуровня; - оператор Гамильтона в текущей лагранжевой системе координат; V, - вектор скорости перемещений макроуровня и мезоуровня; - мера скорости деформации, удовлетворяющая требованию объективности; , - упругая и неупругая составляющие меры скорости деформации мезоуровня; - плотность материала в отсчетной и текущей конфигурациях; - скорость сдвига по k-й системе внутризеренного скольжения; - скорость сдвига по системе скольжения при достижении касательным напряжением критического напряжения сдвига; m - показатель скоростной чувствительности материала; - единичные векторы направления скольжения и нормали к плоскости k-й системы скольжения; H - функция Хевисайда; - сдвиговые и критические напряжения k-й системы; - составляющие критических напряжения k-й системы, обусловленные соответственно взаимодействием дислокаций с дислокационной структурой и границами зерен; - матрица, описывающая упрочнение кристалла за счет взаимодействия дислокаций с дислокациями леса; qlat - параметр латентного упрочнения; - напряжения насыщения; h0, a - параметры, описывающие упрочнение материала; Ns, Nf - число систем скольжения и число плоских участков границы (фасеток) рассматриваемого зерна; S, Sj - площадь зерна и площадь j-й фасетки границы зерна; - параметр, характеризующий взаимодействие дислокаций s-й системы скольжения рассматриваемого зерна c j-й фасеткой (определение этого параметра приведено в работах [51-52]); - компонента касательных барьерных напряжений, действующих на k-ю систему скольжений со стороны дислокаций ориентационного несоответствия, образованных в фасетке в результате скольжения дислокаций по s-й системе; η - безразмерный параметр, определяемый в ходе процедуры идентификации модели; - тензор упругих свойств рассматриваемого кристаллита, компоненты которого определены и постоянны в актуальном базисе ki (в отсчетной конфигурации векторы базиса обозначаются как ) жесткой подвижной системы координат (ПСК), вращающейся со спином w. Жесткая ортонормированная система координат связана с одним кристаллографическим направлением и кристаллографической плоскостью, содержащей это направление. В отсчетной конфигурации базис ПСК совпадает с базисом кристаллографической системы координат (КСК). Базис КСК при деформировании может искажаться, базис ПСК является квазитвердым. Движение относительно ПСК является собственно деформационным. С определением ПСК и спина w можно ознакомиться в работах [53-54]); - тензор, совмещающий ПСК с лабораторной системой координат (ЛСК); - дельта Кронекера; точка над соответствующими переменными мезоуровня означает производную по времени t, верхний индекс «T» - операцию транспонирования соответствующей тензорной величины 2-го ранга. Переход к переменным макроуровня осуществляет в результате осреднения по объему соответствующих переменных мезоуровня. Для описания упрочнения за счет взаимодействия подвижных дислокаций с дислокациями леса (соотношение для ) принимался известный закон упрочнения [55-56], закон зернограничного упрочнения (соотношение для ) описан в работах [51-52]. При высокотемпературном деформировании рассматривается эффект, связанный с началом прерывистой динамической рекристаллизации. Полагается справедливым критерий Бейли - Хирша образования новых рекристаллизованных зерен, согласно которому уменьшение локальной объемной энергии за счет устранения дефектов должно быть больше или равно увеличению зернограничной энергии в результате увеличения протяженности границы при образовании рекристаллизованного зерна [37]: , (14) где - увеличение площади границы при изменении объема зерна на величину в результате рекристаллизации, - разность удельных на единицу объема запасаемых энергией в зернах i и j, - удельная на единицу площади энергия межзеренной границы между i-м и j-м зернами. Критерий Бейли - Хирша сформулирован для описания образования новых зерен из зародышей рекристаллизации (субзерен). Вопрос моделирования эволюции размеров субзерен r выходит за рамки данной работы. Принимается, что исходный размер субзерен подчиняется закону Рэлея [34-35] и в процессе деформирования не изменяется. На рис. 2 показана гистограмма распределения размеров субзерен по закону Рэлея со средним размером субзерен r0 равным 0,25 мкм [34]. Зародыши рекристаллизации расположены вблизи границы, которая представлена совокупностью фасеток. При численной реализации каждой фасетке рассматриваемого зерна приписывалась выборка размеров субзерен (зародышей рекристаллизации) согласно закону Рэлея. Критерий рекристаллизации (14) проверялся для каждого субзерна. Рис. 2. Гистограмма распределения размеров субзерен по закону Рэлея Fig. 2. Histogram of subgrain size distribution according to Rayleigh’s law В научной литературе отсутствует четкое определение критической деформации εс. Обычно величина εс определяется как деформация, при которой начинается процесс рекристаллизации [6, 15]. Следует обратить внимание, что локально выполнение критерия рекристаллизации возможно в некоторых областях металлических образцов, особенно если речь идет о поликристаллах. В рамках данной работы исследуется представительный объем поликристалла. Между каждой парой соседних зерен, одно из которых ориентировано на легкое скольжение, а во втором скольжение затруднено, может выполняться критерий (14). При этом в целом в объеме поликристалла между другими зернами этот критерий не будет выполняться. Вследствие этого для определения величины критической деформации εс вводится в рассмотрение объемная доля зародышей рекристаллизации (субзерен) (по отношению к представительному объему) fr, для которых выполняется критерий рекристаллизации (14). Характерная зависимость fr от интенсивности деформаций εu показана на рис. 3. Полагается, что значению критической деформации εс соответствует начальный этап активного перехода зародышей рекристаллизации в новые рекристаллизованные зерна согласно критерию (14). В этом случае мерой критической Рис. 3. Характерная зависимость доли активных зародышей fr от интенсивности деформации εu Fig. 3. Characteristic dependence of the active nuclei fraction fr on the deformation intensity εu деформации εс может служить угол наклона φ касательной к графику функции . В настоящей работе принято, что достижение угла критического значения φс, равного 15°, соответствует значению критической деформации εс инициации рекристаллизации. 3. Результаты моделирования и их анализ В работе исследуется неупругое деформирование поликристалла меди, используемые параметры многоуровневой модели приведены в табл. 1. Данные по доле запасенной энергии a в поликристалле меди сильно отличаются в различных источниках, приемлемым можно считать значение 10 % согласно обзору, представленному в работе [57]. Таблица 1 Параметры материала многоуровневой модели неупругого деформирования Table 1 Material parameters of the inelastic deformation multilevel model Параметр Значение Источник 170 ГПа [58] 124 ГПа [58] 75 ГПа [58] 16 МПа [56] 148 МПа Процедура идентификации h0 300 МПа Процедура идентификации a 1,5 Процедура идентификации qlat 1,4 [56] 0,001 с-1 [55] m 83 [55] η 1,4 Процедура идентификации a 0,10 [57] egb 0,625 Дж/м2 [6] r0 0,25 мкм [34] Для реализации расширенной статистической модели необходимо задать зеренную структуру поликристалла в отсчетной конфигурации. С процедурой ее формирования можно ознакомиться в работе [59]. Для этого на основе экспериментальных данных устанавливаются статические законы для распределения среднего размера зерна и сферичности (отношение площади поверхности сферы (с объемом, равным объему рассматриваемого зерна) к площади поверхности зерна) зерен поликристалла меди [60]. С применением свободно распространяемого программного обеспечения Neper [61] построена зеренная структура в трехмерном пространстве, которая далее передается в расчетный модуль статистической модели неупругого деформирования. В приведенных ниже численных экспериментах топология зеренной структуры задавалась упрощенно - в виде кубов с ребром R, соответствующего характерному размеру зерна 3 мкм [60]. В статистической модели в качестве одного из параметров в каждый момент деформирования устанавливается ориентация подвижной системы координат рассматриваемого зерна относительно условно неподвижной лабораторной системы координат, описываемая тензором o. Тогда взаимная разориентация соседних i-го и j-го зерен будет определена тензором . Далее по этому тензору определялся угол взаимной разориентации φr. Следует обратить внимание, что для одной и той же текстуры (или ее отсутствия) взаимная ориентация между соседними зернами может отличаться. Этого можно добиться изменением «назначения» конкретной ориентации зернам из сгенерированной фиксированной выборки ориентаций. В отсчетной конфигурации рассмотрено равномерное распределение ориентаций ПСК зерен, полюсная фигура для выборки этого распределения плоскости {001} показана рис. 4 (а). Для сгенерированной выборки ориентаций зерен был реализован алгоритм перебора взаимной разориентации соседних зерен с целью получения повышенных и пониженных средних значений взаимных разориентации; далее эти варианты обозначаются соответственно как и . На рис. 4 (б-в) показаны гистограммы распределения углов взаимной разориентации φr для вариантов и . Средние значения углов φr получились равными 53° и 27° для вариантов и соответственно. В численных экспериментах рассматриваются два этапа неупругого деформирования. На первом этапе анализируется холодное неупругое деформирование представительного макрообъема поликристаллической меди. Законы кинематического нагружения заданы предписанным градиентом скорости перемещений, рассматриваются два варианта холодной деформации - простой сдвиг и сжатие : (15) где pi - ортонормированный базис ЛСК. В расчетах приняты следующие значения параметров кинематического нагружения: . Многоуровневая модель неупругого деформирования (1)-(13) построена в рамках геометрически нелинейной постановки, учитывающей взаимодействие соседних кристаллитов; различные взаимные разориентировки зерен φr при фиксированном распределении ориентаций будут приводить к активации различных систем скольжения, следовательно, к различным поворотам зерен [54]. На рис. 5 показаны полюсные фигуры плоскостей {100}, {110}, {111} для кинематического воздействия при и величине накопленной деформации εu, равной 0,45. При воздействии для варианта получается близкая полюсная фигура. На втором этапе рассматривалось горячее одноосное растяжение представительного макрообъема при повышенных гомологических температурах: . (16) При этом, поскольку рассматривается медь (материал с низким значением энергии дефекта упаковки), полагается, что процессами возврата, ведущими к уменьшению критических напряжений, можно пренебречь. Единственный аспект высокотемпературного деформирования, который учитывается в расчетах, - это выход на начальную стадию рекристаллизации, определяемую критерием (14), и вычисление величины критической деформации εс. Были рассмотрены различные значения величины интенсивности предварительной холодной деформации εpr до величины 0,15 с шагом 0,05, после которой осуществлялся этап горячего одноосного деформирования. Результаты проведенных вычислительных экспериментов сведены в табл. 2. а б в Рис. 4. Полюсная фигура {001} для выборки равномерного распределения ориентаций подвижной системы координат (а), гистограммы углов взаимной разориентации в случае (б) и (в) Fig. 4. The pole figure {001} of uniform distribution orientations of the moving coordinate system (a), histograms of mutual misorientation angles of φr in case (б) and (в) Рис. 5. Расчетные полюсные фигуры {100}, {110}, {111} при εu = 0,45 для заданного и углов взаимной разориентации Fig. 5. Calculated pole figures {100}, {110}, {111} at εu = 0.45 for and mutual misorientation angles Таблица 2 Изменения параметров структуры в результате двухэтапного нагружения Table 2 Evolution structure parameters as a result of two-stage loading № εиp φr0 φ1, ° φ2, ° φ3, ° fr εс 1 0 - 53 53 50,5 0,070 0,194 2 0 - 27 27 28,4 0,015 0,235 3 0,05 53 52,4 50,2 0,013 0,183 4 0,10 53 52,4 50,4 0,016 0.129 5 0,15 53 52,3 50,5 0,018 0,065 6 0,05 27 27.1 28,6 0,023 0,233 7 0,10 27 27.2 28,5 0,022 0,182 8 0,15 27 27,3 28,3 0,018 0,129 9 0,05 53 52,5 50,2 0,013 0,188 10 0,10 53 52,4 50,4 0,016 0.135 11 0,15 53 52,3 50,5 0,018 0,072 12 0,05 27 27.2 28,6 0,024 0,235 13 0,10 27 27,4 28,5 0,022 0,185 14 0,15 27 27.5 28,3 0,019 0,135 В табл. 2 использованы следующие обозначения: εup - интенсивность предварительной деформации, φr0 - вариант задания углов взаимной разориентации в отсчетной конфигурации, - предписанный градиент скорости перемещений предварительной деформации, φ1, φ2, φ3 - средний угол взаимной разориентации между соседними зернами на начало предварительного этапа пластической деформации (в конфигурации К01), начало этапа горячего деформирования (К02), окончания горячего деформирования при достижении величины критической деформации εс, отчитываемой от конфигурации К02, fr - доля активных зародышей рекристаллизации. Физические модели позволяют получать более детальную информацию об эволюции материальной структуры. Приведенные в табл. 2 результаты вычислительных экспериментов № 1, 2 соответствуют случаю отсутствия предварительной холодной деформации. Значение εс для всех проведенных вычислительных экспериментов находится в интервале от 0,072 до 0,235. Полученные результаты свидетельствуют о том, что на величину εс большее влияние оказывает угол взаимной разориентации между соседними зернами φr, а не формируемая текстура. При прочих равных условиях за счет изменения φr удается сместить εс на величину 0,064 (эксперименты № 5 и 8) и 0,063 (эксперименты 11 и 14). На рис. 6 показана эволюция среднего значения взаимного угла разориентации φr для двух вариантов предварительного нагружения и величины предварительной холодной деформации εup = 0,1 с последующим горячим деформировании до величины εс. Смена вида нагружения (при переходе от первого ко второму этапу) хорошо видна на полученных зависимостях рис. 6 (резкое изменение значения производной ). При назначении случайным образом ориентаций соседних зерен из сгенерированной выборки (без применения алгоритма по получению , ), угол взаимной разориентации составляет примерно 40°. В этом случае, как и для , формирование текстуры ведет к уменьшению угла φr. Для варианта практически любые развороты ПСК приводят к увеличению угла φr. Текстура, формируемая в момент окончания отмеченных выше двухэтапных процессов, показана на рис. 7. Следует отметить, что несмотря на наличие предварительной деформации εup = 0,1, финальная текстура в рассматриваемых случаях определяется наложенной деформацией второго этапа . Рис. 6. Эволюция среднего значения угла взаимной разориентации φr при предварительной деформации εup = 0,1 и заданных (а, б) и (в, г) при различных значениях углов взаимной разориентации (а, в) и (б, г) Fig. 6. Evolution of the mean value mutual misorientation angle φr under preliminary deformation εup = 0.1 and prescribed (a, б) and (в, г) for different values of the mutual misorientation angles (a, в) and (б, г) а б в г Рис. 7. Текстура, сформированная в момент окончания процесса предварительной деформации εup = 0,1 и заданных (а, б) и (в, г), при различных значениях углов взаимной разориентации (а, в) и (б, г) Fig. 7. Simulated texture formed at the end of the process with preliminary deformation εup = 0.1 and prescribed (a, б) and (в, г) for different values of the mutual misorientation angles (a, в) and (б, г) Рис. 8. Гистограммы распределения модуля разности запасенной энергии между соседними зернами в момент окончания процесса деформирования с предварительной деформации εup = 0,1 и заданных (а, б) и (в, г), при различных значениях углов взаимной разориентации (а, в) и (б, г) Fig. 8. Histograms of the modulus stored energy difference between grains at the end of the process with preliminary deformation εup = 0.1 and prescribed (a, б) and (в, г) for different values of the mutual misorientation angles (a, в) and (б, г) На рис. 8 для исследуемых процессов показана гистограмма распределения модуля разности запасенной энергии между соседними зернами. Полученные результаты свидетельствуют, что важнейшее влияние на величину критической деформации оказывает угол взаимной разориентации соседних зерен - большие значения последнего способствуют более раннему началу динамической рекристаллизации. Меньшее влияние на εс оказывает текстура материала. Полученные результаты расчетов демонстрируют, что формирование деформационной текстуры приводит к уменьшению угла взаимной разориентации зерен для случаев или задания разориентации соседних зерен по равномерному закону и, следовательно, уменьшению интенсивности процесса динамической рекристаллизации. Заключение Интенсивная пластическая деформация сопровождается глубокой перестройкой материальной структуры, как дефектной, так и зеренной. Поликристаллический материал, подвергаемый горячей деформацией, испытывает большое разнообразие процессов, определяющих его структуру, наиболее значимым из которых является рекристаллизация. В работе исследована проблема определения влияния текстуры и угла взаимной разориентации между соседними зернами на критическое значение интенсивности деформаций инициации динамической рекристаллизации. Для этого использован аппарат многоуровневого моделирования с введением внутренних переменных, позволяющий в явном виде исследовать формирование и эволюцию структуры материал. Показано, что определяющее влияние на процесс рекристаллизации оказывает угол взаимной разориентации между соседними зернами. Обычно формирование деформационной текстуры приводит к уменьшению угла взаимной разориентации и, как следствие, к уменьшению интенсивности динамической рекристаллизации. Реальный процесс динамической рекристаллизации сопровождается появлением новых рекристаллизованных зерен, которые окружены большеугловыми границами с соседними зернами. Другими словами, формирование новых зерен будет способствовать увеличению угла взаимной разориентации и, соответственно, дальнейшей реализации динамической рекристаллизации. В данной работе этот эффект не описан, поскольку представленная модель ориентирована на исследование начальной стадии рекристаллизации и определение значения критической деформации εс. Для дальнейшего исследования актуальными представляются рассмотрение эволюции распределения взаимной ориентации зерен в ходе процессов рекристаллизации различных металлов и сплавов, образование новых высокоугловых границ в результате образования новых рекристаллизованных зерен.About the authors
N. S Kondratev
Perm National Research Polytechnic University
P. V Trusov
Perm National Research Polytechnic University
A. N Podsedertsev
Perm National Research Polytechnic University
References
- Hirth J.P. The influence of grain boundaries on mechanical properties // Metallurgical Transactions. - 1972. - Vol. 3 (12). - P. 3047-3067. doi: 10.1007/BF02661312
- Valiev R.Z., Islamgaliev R.K., Alexandrov I.V. Bulk nanostructured materials from severe plastic deformation // Progress in materials science. - 2000. - Vol. 45 (2). - P. 103-189. doi: 10.1016/S0965-9773(99)00061-6
- Accumulation of Defects During Plastic Deformation of Polycrystals of the Meso-and Microscale Grain Size / N. Koneva, L. Trishkina, N. Popova, E. Kozlov // Russian Physics Journal. - 2014. - Vol. 57 (2). - P. 187-196. doi: 10.1007/s11182-014-0224-8
- Plastic deformation of nanostructured materials / A.M. Glezer, E.V. Kozlov, N.A. Koneva, N.A. Popova, I.A. Kurzina. - Boca Raton: CRC Press, 2017. - 334 p.
- Montheillet F. Moving grain boundaries during hot deformation of metals: dynamic recrystallization // Moving Interfaces in Crystalline Solids. - 2004. - Vol. 453. - P. 203-256. doi: 10.1007/3-211-27404-9_5
- Rollett A., Rohrer G.S., Humphreys J. Recrystallization and related annealing phenomena. - Oxford: Elsevier, 2017. - 734 p.
- Panin V.E. Overview on mesomechanics of plastic deformation and fracture of solids // Theoretical and applied fracture mechanics. - 1998. - Vol. 30 (1). - P. 1-11. doi: 10.1016/S0167-8442 (98) 00038-X
- Hertzberg R.W., Vinci R.P., Hertzberg J.L. Deformation and fracture mechanics of engineering materials. - New York: John Wiley & Sons, 2020. - 784 p.
- Mechanical properties of aluminium-graphene composite synthesized by powder metallurgy and hot extrusion / S.N. Kumar, R. Keshavamurthy, M.R. Haseebuddin, P.G. Koppad // Transactions of the Indian Institute of Metals. - 2017. - Vol. 70 (3). - P. 605-613. doi: 10.1007/s12666-017-1070-5
- Nayak K.C., Date P.P. Hot deformation flow behavior of powder metallurgy based Al-SiC and Al-Al2O3 composite in a single step and two-step uni-axial compression // Materials Characterization. - 2019. - Vol. 151. - P. 563-581. doi: 10.1016/j.matchar.2019.03.047
- Comparison of hot deformation behaviour and microstructural evolution for Ti-5Al-5V-5Mo-3Cr alloys prepared by powder metallurgy and ingot metallurgy approaches / Q. Zhao, F. Yang, R. Torrens, L. Bolzoni // Materials & Design. - 2019. - Vol. 169. - P. 107682. doi: 10.1016/j.matdes.2019.107682
- Hot Deformation Characteristics and Processing Map of 1Cr12Ni2Mo2WVNb Martensitic Stainless Steel / C. Zhao, J. Zhang, B. Yang, Y.F. Li, J.F. Huang, Y. Lian // Steel Research International. - 2020. - Vol. 91 (7). - P. 2000020. doi: 10.1002/srin.202000020
- Thermo-mechanical processing of metallic materials / B. Verlinden, J. Driver, I. Samajdar, R.D. Doherty. - Oxford: Elsevier, 2007. - 560 p.
- Горелик С.С., Добаткин С.В., Капуткина Л.М. Рекристаллизация металлов и сплавов. - М.: МИСИС, 2005. - 432 c.
- Huang K., Logé R.E. A review of dynamic recrystallization phenomena in metallic materials // Materials & Design. - 2016. - Vol. 111. - P. 548-574. doi: 10.1016/j.matdes.2016.09.012
- Thermal Mechanisms of Grain Refinement in Steels: A Review / Z. Nasiri, S. Ghaemifar, M. Naghizadeh, H. Mirzadeh // Met. Mater.Int. - 2021. - Vol. 27. - P. 2078-2094. doi: 10.1007/s12540-020-00700-1
- Current issues in recrystallization: a review / R.D. Doherty, D.A. Hughes, F.J. Humpheys, J.J. Jonas, D. Juul Jensen, M.E. Kassner, W.E. King, T.R. McNelley, H.J. McQueen, A.D. Rollett // Materials Science and Engineering: A. - 1997. - Vol. 238 (2). - P. 219-274. doi: 10.1016/S0921-5093 (97) 00424-3
- Honeycombe R.W.K. The plastic deformation of metals. - London: Hodder Arnold, 1968. - 483 p.
- Трусов П.В., Швейкин А.И. Многоуровневые модели моно- и поликристаллических материалов: теория, алгоритмы, примеры применения. - Новосибирск: Изд-во СО РАН, 2019. - 605 с. doi: 10.15372/MULTILEVEL2019TPV
- Valiev R.Z., Langdon T.G. Principles of equal-channel angular pressing as a processing tool for grain refinement // Progress in materials science. - 2006. - Vol. 51 (7). - P. 881-981. doi: 10.1016/j.pmatsci.2006.02.003
- Armstrong R.W. The influence of polycrystal grain size on several mechanical properties of materials // Metallurgical and Materials Transactions B. - 1970. - Vol. 1 (5). - P. 1169-1176. doi: 10.1007/BF02900227
- Grain size effects on dislocation and twinning mediated plasticity in magnesium / H. Fan, S. Aubry, A. Arsenlis, J.A. El-Awady // Scripta Materialia. - 2016. - Vol. 112. - P. 50-53. doi: 10.1016/j.scriptamat.2015.09.008
- Crystal plasticity simulation of the effect of grain size on the fatigue behavior of polycrystalline Inconel 718 / A. Cruzado, S. Lucarini, J. LLorca, J. Segurado // International Journal of Fatigue. - 2018. - Vol. 113. - P. 236-245. doi: 10.1016/j.ijfatigue.2018.04.018
- Knezevic M., Beyerlein I.J. Multiscale modeling of microstructure - property relationships of polycrystalline metals during thermo-mechanical deformation // Advanced Engineering Materials. - 2018. - Vol. 20 (4). - P. 1700956. doi: 10.1002/adem.201700956
- Predicting elastic anisotropy of dual-phase steels based on crystal mechanics and microstructure / A.M. Cantara, M. Zecevic, A. Eghtesad, C.M. Poulin, M. Knezevic // International Journal of Mechanical Sciences. - 2019. - Vol. 151. - P. 639-649. doi: 10.1016/j.ijmecsci.2018.12.021
- Superplasticity induced by the competitive DRX between BCC beta and HCP alpha in Ti-4Al-3V-2Mo-2Fe alloy /j. Shen, Y. Sun, Y. Ning, H. Yu, Z. Yao, L. Hu // Materials Characterization. - 2019. - Vol. 153. - P. 304-317. doi: 10.1016/j.matchar.2019.05.014
- High strain rate superplasticity in Al-Zn-Mg-based alloy: Microstructural design, deformation behavior, and modeling / O. Yakovtseva, M. Sitkina, A.O. Mosleh, A. Mikhaylovskaya // Materials. - 2020. - Vol. 13 (9). - P. 2098. doi: 10.3390/ma13092098
- Shveykin A., Trusov P., Sharifullina E. Statistical crystal plasticity model advanced for grain boundary sliding description // Crystals. - 2020. - Vol. 10 (9). - P. 822. doi: 10.3390/cryst10090822
- Trusov P.V., Shveykin A.I. Multilevel crystal plasticity models of single-and polycrystals. Direct models // Physical mesomechanics. - 2013. - Vol. 16 (2). - P. 99-124. doi: 10.1134/S1029959913020021
- Trusov P.V., Shveykin A.I. Multilevel crystal plasticity models of single-and polycrystals. Statistical models // Physical mesomechanics. - 2013. - Vol. 16 (1). - P. 23-33. doi: 10.1134/S1029959913010037
- Dehghan-Manshadi A., Hodgson P.D. Dependency of recrystallization mechanism to the initial grain size // Metallurgical and Materials Transactions A. - 2008. - Vol. 39 (12). - P. 2830-2840. doi: 10.1007/s11661-008-9656-5
- Roberts W., Boden H., Ahlblom B. Dynamic recrystallization kinetics // Metal Science. - 1979. - Vol. 13 (3-4). - P. 195-205. doi: 10.1179/msc.1979.13.3-4.195
- Ponge D., Gottstein G. Necklace formation during dynamic recrystallization: mechanisms and impact on flow behavior // ActaMaterialia. - 1998. - Vol. 46 (1). - P. 69-80. doi: 10.1016/S1359-6454 (97) 00233-4
- Zurob H.S., Brechet Y., Dunlop J. Quantitative criterion for recrystallization nucleation in single-phase alloys: Prediction of critical strains and incubation times // Acta Materialia. - 2006. - Vol. 54 (15). - P. 3983-3990. doi: 10.1016/j.actamat.2006.04.028
- Modelling discontinuous dynamic recrystallization using a physically based model for nucleation / D.G. Cram, H.S. Zurob, Y.J.M. Brechet, C.R. Hutchinson // ActaMaterialia. - 2009. - Vol. 57 (17). - P. 5218-5228. doi: 10.4028/www.scientific.net/MSF.715-716.492
- Beck P.A., Sperry P.R. Strain induced grain boundary migration in high purity aluminum // Journal of applied physics. - 1950. - Vol. 21 (2). - P. 150-152. doi: 10.1063/1.1699614
- Bailey J.E., Hirsch P.B. The recrystallization process in some polycrystalline metals // Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. - The Royal Society, 1962. - Vol. 267 (1328). - P. 11-30. doi: 10.1098/rspa.1962.0080
- Sakai T., Akben M.G., Jonas J.J. Dynamic recrystallization during the transient deformation of a vanadium microalloyed steel // Acta Metallurgica. - 1983. - Vol. 31 (4). - P. 631-641. doi: 10.1016/0001-6160(83)90053-6
- Sakai T., Jonas J.J. Overview no. 35 dynamic recrystallization: mechanical and microstructural considerations // Acta Metallurgica. - 1984. - Vol. 32 (2). - P. 189-209. doi: 10.1016/0001-6160(84)90049-X
- The identification of dynamic recrystallization and constitutive modeling during hot deformation of Ti55511 titanium alloy / K. Tan, J. Li, Z. Guan, J. Yang, J. Shu // Materials & Design. - 2015. - Vol. 84. - P. 204-211. doi: 10.1016/j.matdes.2015.06.093
- Model-ling for the dynamic recrystallization evolution of Ti-6Al-4V alloy in two-phase temperature range and a wide strain rate range / G.Z. Quan, G.C. Luo, J.T. Liang, D.S. Wu, A. Mao, Q. Liu // Computational Materials Science. - 2015. - Vol. 97. - P. 136-147. doi: 10.1016/j.commatsci.2014.10.009
- Sellars C.M. Recrystallization of metals during hot deformation // Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences. - 1978. - Vol. 288 (1350). - P. 147-158. doi: 10.1098/rsta.1978.0010
- Effects of het-erogeneity on recrystallization kinetics of nanocrystalline copper pre-pared by dynamic plastic deformation / F. Lin, Y. Zhang, N. Tao, W. Pantleon, D.J. Jensen // Acta Materialia. - 2014. - Vol. 72. - P. 252-261. doi: 10.1016/j.actamat.2014.03.036
- Micro-structure-based modeling of residual yield strength and strain hardening after fire exposure of aluminum alloy 5083-H116 / P.T. Summers, A.P. Mouritz, S.W. Case, B.Y. Lattimer // Materials Science and Engineering: A. - 2015. - Vol. 632. - P. 14-28. doi: 10.1016/j.msea.2015.02.026
- Radhakrishnan B., Sarma G.B., Zacharia T. Modeling the kinetics and microstructural evolution during static recrystallization-Monte Carlo simulation of recrystallization // Acta materialia. - 1998. - Vol. 46. (12). - P. 4415-4433. doi: 10.1016/S1359-6454 (98) 00077-9
- Peczak P. A Monte Carlo study of influence of deformation temperature on dynamic recrystallization // Acta Metallurgica et Materialia. - 1995. - Vol. 43 (3). - P. 1279-1291. doi: 10.1016/0956-7151 (94) 00280-U
- A polycrystal plasticity based discontinuous dynamic recrystallization simulation method and its application to copper / G. Zhou, Z. Li, D. Li, Y. Peng, H.S. Zurob, P. Wu // International Journal of Plasticity. - 2017. - Vol. 91. - P. 48-76. doi: 10.1016/j.ijplas.2017.01.001
- Beyerlein I.J., Knezevic M. Review of microstructure and micromechanism-based constitutive modeling of polycrystals with a low-symmetry crystal structure // Journal of Materials Research. - 2018. - Vol. 33 (22). - P. 3711-3738. doi: 10.1557/jmr.2018.333
- Zhao P., Wang Y., Niezgoda S.R. Microstructural and micromechanical evolution during dynamic recrystallization // International Journal of Plasticity. - 2018. - Vol. 100. - P. 52-68. doi: 10.1016/j.ijplas.2017.09.009
- Sarrazola D.A.R., Muñoz D.P., Bernacki M. A new numerical framework for the full field modeling of dynamic recrystallization in a CPFEM context // Computational Materials Science. - 2020. - Vol. 179. - 109645. doi: 10.1016/j.commatsci.2020.109645
- Кондратьев Н.С., Трусов П.В. О мере разориентации систем скольжения соседних кристаллитов в поликристаллическом агрегате // Вестник ПНИПУ. Механика. - 2012. - № 2. - С. 112-127.
- Кондратьев Н.С., Трусов П.В. Описание упрочнения систем дислокационного скольжения за счет границ кристаллитов в поликристаллическом агрегате // Вестник ПНИПУ. Механика. - Пермь: Изд-во ПНИПУ, 2012. - № 3. - С. 78-97.
- Trusov P.V., Shveykin A.I. On motion decomposition and constitutive relations in geometrically nonlinear elastoviscoplasticity of crystallites // Physical Mesomechanics. - 2017. - Vol. 20 (4). - P. 377-391. doi: 10.1134/S1029959917040026
- Shveykin A.I., Trusov P.V., Kondratev N.S. Multiplicative Representation of the Deformation Gradient Tensor in Geometrically Nonlinear Multilevel Constitutive Models // Lobachevskii Journal of Mathematics. - 2021. - Vol. 42 (8). - P. 2047-2055. doi: 10.1134/S1995080221080291
- Bronkhorst C.A., Kalidindi S.R., Anand L. Polycrystalline plasticity and the evolution of crystallographic texture in FCC metals // Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences. - 1992. - Vol. 341 (1662). - P. 443-477. doi: 10.1098/rsta.1992.0111
- Anand L. Single-crystal elasto-viscoplasticity: application to texture evolution in polycrystalline metals at large strains // Computer methods in applied mechanics and engineering. - 2004. - Vol. 193 (48-51). - P. 5359-5383. doi: 10.1016/j.cma.2003.12.068
- Bever M.B., Holt D.L., Titchener A.L. The stored energy of cold work // Progress in materials science. - 1973. - Vol. 17. - P. 5-177. doi: 10.1016/0079-6425 (73) 90001-7
- Simmons G., Wang H. Single Crystal Elastic Properties and Calculated Aggregate Properties: A HANDBOOK // Cambridge (MA): MITPress, 1971. - 370 p.
- Kondratev N.S., Podsedertsev A.N., Trusov P.V. The grain structure construction of polycrystals for modified two-level crystal plasticity statistical models // Procedia Structural Integrity. - 2021. - (in print)
- Suresh K.S., Rollett A.D., Suwas S. Evolution of microstructure and texture during deformation and recrystallization of heavily rolled Cu-Cu multilayer // Metallurgical and Materials Transactions A. - 2013. - Vol. 44 (8). - P. 3866-3881.
- Quey R., Renversade L. Optimal polyhedral description of 3D polycrystals: Method and application to statistical and synchrotron X-ray diffraction data // Computer Methods in Applied Mechanics and Engineering. - 2018. - Vol. 330. - P. 308-333. doi: 10.1016/j.cma.2017.10.029