Полный текст

Введение Способность металлических материалов в достаточно высокой степени деформироваться пластически лежит в основе всей промышленной цивилизации. В силу большой сложности процессов пластического деформирования и многообразия физических механизмов реализации пластических деформаций первые модели их теоретического описания появились значительно позднее моделей упругого деформирования. В настоящее время сложились три основных подхода к проблеме пластического деформирования поликристаллических материалов, связанные в первую очередь с масштабными уровнями рассмотрения структуры: макроскопический (однородное тело), мезоскопический (учет размеров и структуры зерен), микроскопический (уровень отдельных дислокаций). На макроскопическом уровне применяются чисто феноменологические континуальные теории, базирующиеся на разложении полной деформации на упругую и пластическую части. Пластическая часть оценивается из термодинамических соображений касательно диссипативных процессов в материале (принцип максимальной пластической работы Р. Хилла и т.д.). Для учета влияния зеренной структуры вводятся некоторые внутренние переменные, подлежащие идентификации. Феноменологический подход использует классические модели пластичности. При мезоскопическом рассмотрении пластические свойства поликристаллического агрегата выводятся из свойств и закономерностей пластического деформирования монокристаллов, соответствующих зернам материала. Свойства отдельного зерна и соответствующего (макроскопического) монокристалла считаются идентичными. Физическим механизмом деформирования зерна является сдвиг в определенных кристаллографических системах скольжения. Решается краевая задача деформирования поликристалла как составного тела. Большое распространение получили методы компьютерного моделирования структуры поликристаллических тел с последующим решением для них краевых задач численными методами. Для зерен-монокристаллов используются различные модели пластического деформирования. В англоязычной литературе это наиболее часто применяемый подход в исследовании пластичности поликристаллов. Соответствующее направление получило название вычислительной кристаллической пластичности (Computational Crystal Plasticity). Разработано много вариантов этого метода, опубликованы многочисленные обзоры [1-5]. Метод вычислительной кристаллической пластичности позволяет проследить эволюцию деформирования индивидуальных зерен, что особенно важно в проблемах высокоцикловой и гиперцикловой усталости, где зарождение усталостных повреждений происходит в результате накопления малых пластических деформаций в отдельных зернах. Очаги повреждений возникают в специфически ориентированных зернах, окруженных характерной микроструктурой. Эти зерна получили название «мягкие зерна», а микроструктуры - название микроструктурных экстримов или экстремальных микроструктур [6, 7, 8]. Для поиска таких микроструктур применяют различные вычислительные и экспериментальные методы [9, 10, 11]. Вероятность образования экстремальных микроструктур весьма мала, но они являются причиной большого разброса усталостной долговечности [12] в режимах сверхвысокоцикловой усталости. Вычислительный поиск случайно образующихся экстремальных структур методом компьютерного моделирования практически невозможен ввиду необходимости генерировать колоссальное количество моделей. На микроскопическом уровне используются дискретные или континуальные модели отдельных дислокаций и их взаимодействий [13]. Для связи микроскопических параметров, применяемых на этом уровне, с макроскопическим пластическим откликом вводятся дополнительные модельные соображения. Все три вышеупомянутых подхода активно развиваются многочисленными научными школами в теории пластичности. В настоящей работе используется мезоскопический уровень рассмотрения проблемы пластического деформирования поликристаллов. Для решения краевых задач упругопластического деформирования моно- и поликристаллов развивается новый метод, основанный на теоретико-полевом подходе к проблеме деформирования неоднородных тел. 1. Краевая задача упругопластического деформирования поликристаллических тел в теоретико-полевой формулировке Рассмотрим произвольное неоднородное упругопластическое тело. Пусть - тензор полной глобальной неоднородной деформации в теле; - полный неоднородный вектор перемещений. Будут рассматриваться только малые упругопластические деформаций, так что полные деформации и перемещения связаны соотношением , (1) где - градиент перемещений. Для деформаций и векторов перемещений считается справедливым аддитивное разложение на упругие и пластические части: , (2) . Будем использовать деформационную теорию пластичности. Определяющие соотношения для всего тела запишем в виде , где - глобальный тензор секущих модулей упругости. Секущий модуль зависит не только от координат, но и от текущего значения тензора деформаций. Однако для краткости эту зависимость явно не будем указывать, но всегда подразумевать. Считаем тензор секущих модулей обладающим обычными свойствами симметрии по перестановкам индексов: . (3) Будем рассматривать только макроскопически изотропные тела. Это означает, что осредненный по объему тела тензор (4) однороден и изотропен, то есть представим в виде , (5) где и - шаровая и девиаторная части единичного тензора 4-го ранга , где - символ Кронекера, а и - осредненные по объему тела модули объемного сжатия и сдвига. Для получения искомых уравнений теоретико-полевого подхода поставим две краевые упругопластические задачи в традиционной дифференциальной форме. Первая задача для неоднородного упругопластического тела, не обязательно поликристаллического, с кинематическими граничными условиями имеет вид (6) . Здесь - внешняя поверхность тела; - заданные перемещения на поверхности тела. Разлагаем тензор секущих модулей на осредненную по объему тела и флуктуирующую составляющие: . (7) Осредненный тензор не зависит от координат, но зависит от осредненной макродеформации тела. Осредненная макродеформация является в данном контексте внешним параметром и не затрагивает структуру уравнений (не участвует в дифференцировании по координатам). Уравнения равновесия в (6) примут вид . (8) Введем функцию Грина для однородного упругопластического тела с осредненным тензором секущих модулей упругости посредством уравнения . (9) Для большого тела за решение (9) можно взять функцию Грина для неограниченной среды, обращающуюся на бесконечности в ноль. Для нахождения ее вида поступим так же, как это делается для упругой задачи [14]. Перейдем в (9) к образам Фурье-переменных по формулам , . Однородный тензор не зависит от координат и не подвергается Фурье-преобразованию. Для Фурье-образа функции Грина получается уравнение , по форме совпадающее с аналогичным уравнением для функции Грина упругой задачи [14]. Поэтому для среды с однородным изотропным тензором типа (5) искомая функция Грина , (10) где - единичный вектор вдоль направления радиуса-вектора , а - упругопластический материальный параметр. Выражение (10) внешне совпадает с упругой функцией Грина, если в последней формально заменить средние упругие модули на секущие , . Упругопластическая функция Грина зависит от осредненного тензора деформации как от внешнего параметра. Производные функции Грина обладают очевидными свойствами: (11) Умножим уравнение (8) на и просуммируем по i, одновременно переобозначив аргумент : Уравнение (9) умножаем на и также суммируем по i: . В первых членах последних двух уравнений распространим производные по на все сомножители, вычтя соответствующие члены: (12) (13) В уравнениях (12) и (13) свободный индекс m, немые индексы (по которым производится суммирование) можно переобозначать. Сделав во вторых членах (12) и (13) замену , можно видеть, что они одинаковы. Вычитая из (12) уравнение (13) и применяя теорему Гаусса для объемного интеграла от первых членов (с дивергенциями), получим (14) Уравнение (14) есть интегродифференциальное уравнение для глобального вектора перемещений в рамках деформационной модели пластичности неоднородного тела. Вторая краевая задача ставится для однородного упругопластического тела с осредненным тензором секущих модулей. Вектор перемещений для этого тела отметим звездочкой . Уравнения равновесия в перемещениях и граничные условия для этого тела , (15) где - заданные перемещения на поверхности тела для второй задачи. Вводим для этой задачи функцию Грина, совпадающую с функцией Грина в уравнении (9): . (16) Умножаем (15) на , а (16) на : (17) Переписываем систему (17), комбинируя производные: (18) (19) Снова замечаем, что вторые члены в уравнениях (18) и (19) одинаковы. Вычитаем из уравнения (18) уравнение (19), интегрируем по по объему тела, применяем теорему Гаусса для членов с дивергенциями и получаем интегродифференциальное уравнение для перемещений в однородном теле (20) Если выбрать граничные условия для уравнений равновесия (8) и (15) обеих задач одинаковыми, то поверхностные интегралы в интегродиффреренциальных уравнениях (14) и (20) совпадают. Вычитая одно уравнение из другого получаем (21) Это интегродифференциальное уравнение для полного вектора перемещений в неоднородном упругопластическом теле. Влияние граничных условий на решение перенесено в задачу для перемещений в однородном теле. Решение однородной задачи можно считать известным (это классическая задача деформационной теории пластичности однородного тела). При этом само решение может быть неоднородным. Таким образом, как и в упругой задаче [15], уравнение (21) разделяет влияние на решение граничных условий и неоднородностей тела. Дифференцированием (21) по переходим к градиентам перемещений: (22) Уравнение (22) неудобно тем, что содержит (в случае микронеоднородных тел, например поликристаллов) производную от быстроосциллирующих или даже разрывных функций . Чтобы уйти от этих производных, применим прием, впервые использованный И.М. Лифшицем в статье 1946 г. [16]. Заменяем функции пробными, совпадающими с исходными везде, кроме узких областей вблизи границ зерен, где пробные функции непрерывны и дифференцируемы, что устраняет неопределенности с производными. Вносим под знак первой производной (чтобы получить полную дивергенцию) и одновременно вычитаем соответствующее выражение: (23) В последней производной в (23) замечаем, что . Преобразуя первый объемный интеграл в поверхностный, получаем еще одну форму интегрального уравнения: (24) Уравнение (24) не содержит производных от разрывных функций, и в нем можно вернуться к исходным быстроосциллирующим и разрывным функциям . Уравнение (24) предпочтительнее для микронеоднородных тел. Для них в поверхностном интеграле стоит произведение быстроосциллирующей около нуля функции на гладкие функции. Поэтому по любому достаточно большому участку поверхности тела (для поликристалла много больше площади сечения зерна), а значит, и для всего тела интеграл обращается в ноль. Получаем искомое точное для больших поликристаллических (или микронеоднородных) тел уравнение для градиентов деформаций . (25) Уравнение (25) отличается от соответствующего уравнения для упругого случая [15] формальной заменой тензора модулей упругости на тензор секущих модулей. В этом удобство использования секущих, а не касательных модулей, применяемых в инкрементальной формулировке задач пластичности, в контексте развиваемого метода. Инкрементальная формулировка позволяет решить более широкий круг задач пластичности (хотя не все задачи пластичности решаемы этим подходом; например, задачи пластической устойчивости точнее решаются в деформационной формулировке). Однако инкрементальные уравнения сложнее, а методы их решения трудоемки. Производя в (25) симметризацию по индексам m и n, переходим к уравнению для тензора полных деформаций , (26) где - модифицированный тензор Грина упругопластической задачи для однородного тела, . (27) Уравнение (26) является исходным для последующего анализа. Для исследования деформаций в отдельных зернах представим глобальный тензор секущих модулей и глобальный тензор деформаций в виде суперпозиции локальных тензоров с помощью индикаторных функций (равных единице внутри зерна и нулю вне его): , (28) , (29) где - однородный тензор секущих модулей ξ-го зерна (рассматриваемого как монокристалл) в глобальной системе координат; - полная мезодеформация ξ-го зерна; M - полное число зерен в поликристалле (макроскопически большое число). Подстановка разложений (28) и (29) в (26) приводит к системе интегральных уравнений для мезодеформаций в зернах: (30) В уравнениях (30) интегралы под знаком суммы отражают неупругое взаимодействие деформаций. В нем удобно выделить в отдельное слагаемое внутризеренное взаимодействие в рассматриваемом ξ-м зерне: (31) Система (31) содержит макроскопически большое число уравнений для деформаций во всех зернах. В теоретическом отношении интерес представляют как решения для каждого зерна (мезодеформации и мезонапряжения), так и осредненные по объему тела макронапряжения. Последнее, будучи найденным, дает макроскопические определяющие соотношения. Ключевым в уравнениях (31) является концепция неупругого взаимодействия деформаций. Для решения уравнений (31) с учетом взаимодействий применяем теорию возмущений по межзеренному взаимодействию [17], отнеся внутризеренное взаимодействие деформаций к нулевому приближению. Разложение решения в бесконечный ряд суммы поправок всех порядков (32) дает бесконечную (точную) систему уравнений для поправок: , (33) ............................................................. и т.д., где - отклонение секущего модуля ξ-го зерна от макроскопического осредненного модуля. Для приближенного решения системы уравнений (33) применим кусочно-постоянную аппроксимацию для мезодеформаций. Таким образом, осуществляется дискретизация континуальной краевой задачи. Ограничимся пренебрежением неоднородностями мезодеформаций внутри отдельных зерен. Тогда все поправки в (33) становятся постоянными числами, и для них получается система алгебраических уравнений , , , (34) , .......................................... и т.д., где - (35) тензорные коэффициенты, зависящие от конфигурации структуры поликристалла. Микроструктура поликристалла задается совокупностью коэффициентов . Каждая строка в (34) есть в действительности подсистема из шести уравнений. Введем новый тензор, обратный тензорному выражению в квадратных скобках в (34), . (36) Вычислительные эксперименты с зернами-монокристаллами различных сингоний и всевозможных случайных ориентаций кристаллографических осей зерен в поликристалле (от этого зависят коэффициенты ) показали, что такое обращение всегда возможно. Для практических вычислений система (34), тензор (36) и аналогичные выражения ниже записываются в матричном виде. Например, тензоры и заменяются на квадратные матрицы и размерностью 9х9. Используя симметрию введенных тензоров по перестановке индексов первой и второй пары, можно сократить ранг матриц до 6х6. Матричные операции выполняются численно с помощью стандартных современных пакетов программ для решения задач линейной алгебры. С введением тензоров (36) система (34) перепишется так: , , , (37) , ........................ и т.д., где введен еще один тензор . (38) Сложив все уравнения в (37), получим уравнение, связывающее полные деформации во всех зернах, . (39) Из последнего уравнения видно, что тензор описывает взаимодействие упругопластических деформаций в ξ-м и η-м зернах, а - внутризеренное взаимодействие. Система (34) нелинейная, так как коэффициенты зависят от искомой полной деформации и мезодеформаций. Таким образом, каждая подсистема в (34) и система в целом являются нелинейными алгебраическими. Они решаются методом последовательных итераций. Применен следующий алгоритм вычислений. Задается фиксированное значение макродеформации . На первой итерации все секущие модули полагаются равными упругим модулям. Численно находятся коэффициенты в левых частях уравнений с этими модулями, и решается линеаризованная система (34). Полученное решение принимается за нулевое приближение. На следующем шаге по найденному нулевому приближению для мезодеформаций уточняются (перевычисляются) значения секущих модулей по формулам из подразд. 2 настоящей статьи, уточняются все коэффициенты в левых частях уравнений, снова решается линеаризованная система. Процесс продолжается до достижения стационарного состояния с требуемой точностью. Вычислительные эксперименты показали, что такой алгоритм дает устойчивое решение с точностью не хуже 1 % не более чем за десять итераций. Алгоритм реализован в среде Java для высокопроизводительной настольной платформы. Значение мезодеформации в фиксированном по форме и ориентации зерне зависит от взаимодействия с окружающими зернами. Интенсивность этого взаимодействия (при фиксированных формах и размерах зерен) зависит от взаимной ориентации кристаллографических осей взаимодействующих зерен и расстояния между ними. В определенных конфигурациях это взаимодействие максимально. При анализе локализованных критических явлений эти конфигурации микроструктуры представляют значительный интерес. Им соответствуют экстремальные флуктуации в фиксированном зерне, которые ведут к инициации разрушения в таких зернах. В приложении теоретико-полевого подхода к упругопластическому деформированию поликристаллов ограничимся анализом таких экстремальных состояний на начальных этапах пластического макродеформирования при малых макроскопических пластических деформациях. Структура левых частей всех уравнений в системе (34) одинакова. Нужен единственный код решения системы нелинейных уравнений размерностью 9х9 или 6х6. Последовательным применением кода сначала для первого уравнения (поправки нулевого порядка), а затем для всех остальных, подставляя в правые части решения для предыдущих поправок, можно получить решения формально с любой точностью. Для практического решения уравнений (34) необходимо задаться моделью деформирования одиночного зерна-монокристалла и исследовать сходимость процедуры последовательного нахождения поправок. Кроме того, нужно оценить минимально необходимое число взаимодействующих зерен, входящих в правые части уравнений. Оно связано со скоростью уменьшения интенсивности взаимодействия при увеличении расстояния между зернами. 2. Модель упругопластического деформирования монокристалла зерна В работе принимается, что единственным физическим механизмом пластической деформации является скольжение по кристаллографическим плоскостям с наибольшими плотностями упаковки атомов в решетках монокристаллов. Скольжение происходит в соответствии с законом Шмида: скольжение наступает при превышении разрешающих напряжений критических значений в каждой системе скольжения. При этом предполагается, что при наличии упрочнения рост (инкремент) разрешающих напряжений превышает скорость (инкремент) роста критических разрешающих напряжений в активных системах скольжения. В этом параграфе все величины относятся к кристаллографической системе координат ξ-го зерна. Пусть зерно-монокристалл подвергнут однородной пластической деформации, так, что в некоторых системах скольжения произошли пластические сдвиги. Будем нумеровать системы скольжения в монокристалле индексом g. Пусть - нормаль к плоскости скольжения системы g, а - направление скольжения в этой системе. Ориентационный тензор Шмида системы g имеет вид . (40) Если обозначить величину скольжения в системе g как , то пластическая часть полной деформации в ξ-м зерне представима в виде [18] , (41) где N - количество систем скольжения. Уравнение (41) можно разрешить относительно . Умножим обе части на и просуммируем по индексам i и j: . (42) Выражение (42) есть система линейных уравнений относительно N неизвестных . Введем квадратную матрицу в пространстве индексов систем скольжения, которые пишем в верхнем положении (в отличие от «координатных» индексов - индексов ортов систем координат, записываемых внизу): . (43) Тогда имеем и решение , (44) где обратная матрица определяется соотношением , а - символ Кронекера по индексам систем скольжения. Можно выразить сдвиги через полные деформации с помощью условий скольжения (N условий), выраженных через деформации . Для этого упругопластические определяющие соотношения деформационного типа запишем в двух вариантах: (45) и , (46) где и - тензоры секущих модулей и модулей упругости соответственно. Разрешающее напряжение сдвига в системе скольжения g связано с тензором напряжений в зерне соотношением . (47) Подставляя в (47) формулы (46) и (44), имеем . (48) Если система скольжения активна, то для нее должно быть (49) или , (50) если , находимое по , отрицательно. Здесь - критическое разрешающее напряжение Шмида в системе скольжения g. Кроме того, при наличии упрочнения должно выполняться дополнительное условие для инкрементов , (51) иначе деформирование будет продолжаться упруго. Ниже для краткости временно опускаем индекс зерна ξ, но все выражения относятся к конкретному зерну. Предполагаем линейный закон упрочнения, который запишем в виде , (52) где - начальное значение критического напряжения; - матрица упрочнения. Подставляем (48) и (52) в условие скольжения (49): (53) Переписываем (53) так: . (54) Выражение (54) есть система NxN нелинейных (так как коэффициенты - квадратные скобки - зависят от ) уравнений относительно . Формально преобразуем ее с помощью матричных операций. В квадратных скобках в (54) стоит квадратная матрица. Обозначим эту матрицу , (55) и перепишем (54) в кратком виде, используя общепринятое правило суммирования по повторяющимся индексам: . (56) Если ввести матрицу , обратную , , (57) то получим формальное решение системы (56) для в виде . (58) Уравнение (58) выражает пластические сдвиги в активных системах через полную деформацию, содержащую упругую и пластическую части. На самом деле (58) не решение, а новое уравнение, так как правая часть зависит от через матрицу . С помощью (58) получим связь пластической части деформации с полной деформацией: (59) Возвращаемся к закону Гука (46) и подставим в него (59): . (60) Сравнивая (60) с определяющим соотношением упругопластического деформирования (45), получим выражение для тензора секущих модулей зерна в кристаллографической системе координат кристалла . (61) Уравнением (61) завершается построение модели упругопластического деформирования монокристалла-зерна в рамках деформационного подхода. Уравнение (61) для секущих модулей по структуре аналогично уравнениям для касательных модулей в инкрементальном подходе, полученным в ряде работ [19-22]. Однако инкрементальные уравнения для деформаций более сложны в решении и в упомянутых работах практически не решаются. Вместо этого применяют различные схемы самосогласования [23, 24] . Зависимость секущего модуля от полной деформации сосредоточена в матрице через и уравнение (58). Эта зависимость, несмотря на внешнюю простоту (61), весьма сложная и нелинейная. Формула (61) справедлива, если хотя бы одна отлична от нуля. Иначе . Это следует из того, что при всех становится несправедливым исходное соотношение (54), в котором левая часть обращается в ноль, а правая нет. В уравнении (61) содержатся тензоры с обычными пространственными индексами (ijkl) и матрицы по индексам систем скольжения (gh). Уравнению (61) можно придать характерный для уравнений квантовой теории взаимодействующих полей «перенормировочный» вид [25]: , (62) где - (63) обычный тензор. Пластические деформации эквивалентны наличию дополнительных фиктивных объемных сил [26, 27]. Можно сказать, что эти силы трансформируют («перенормируют») упругие модули в секущие, а есть перенормирующий множитель, зависящий от полной деформации. В квантовой теории взаимодействие полей описывается величинами типа в лагранжиане, а система взаимодействующих полей сводится к невзаимодействующим полям, но при этом перенормируются свойства полей по формулам типа (62). Идентификация матриц упрочнения для конкретных материалов по экспериментальным данным весьма проблематична. Представляется интересным рассмотреть случай пренебрежимо малого упрочнения и независимости систем скольжения. В этом случае сдвиги в плоскостях скольжения упруго-идеально-пластические. Тогда условие для инкрементов выполняется тривиально, а условия скольжения (49) в системах запишутся так: , g = 1, 2, ... N. (64) В (64) отдельные могут быть положительными или отрицательными, соответственно, они могут давать уменьшающий или увеличивающий вклад (с учетом умножения на тензор Шмида) в разрешающее напряжение . Переписываем (64) в виде системы линейных уравнений для : . (65) Решение для записываем по аналогии с (54) в виде , (66) в котором матрица имеет простую в сравнении с (55) форму: . (67) Геометрический смысл (66) проясняется, если рассмотреть формально случай всего одной плоскости скольжения без упрочнения и в ней одного направления скольжения. Тогда матрица упрощается: (по индексу g нет суммирования), и обратная матриц легко находится: , а для сдвигов получается простое алгебраическое выражение . Легко показать, что - сдвиговая компонента тензора модулей упругости в системе координат с ортами , , в плоскости , - формально упругое разрешающее сдвиговое напряжение в плоскости ; и получаем . Это - закон Гука для сдвига в плоскости . Каждое слагаемое в круглых скобках (66) имеет смысл, если в данной системе h , (68) тогда конкретный вклад в сдвиг положителен, или если выполняется условие , (69) то соответствующий вклад в отрицателен и (66) надо писать в виде . (70) В диапазоне полных деформаций (для каждой системы скольжения) (71) надо полагать соответствующий вклад в равным нулю. Выражения (66) и (70) можно объединить в единую формулу (72) При значениях полной деформации в диапазоне условий (71) выражения (66), (70), (72) для пластических сдвигов неприменимы, они обращаются в ноль. Для связи пластической деформации с полной получается выражение , (73) а закон Гука принимает вид (ниже снова вводим индекс зерна ξ) (74) Его можно записать иначе, тождественно введя во втором члене полную деформацию следующим образом: (75) где - второй инвариант. Тогда для тензора секущих модулей зерна получается следующее выражение: (76) Из экспериментально идентифицируемых параметров в выражение для секущих модулей (76) входят только критические разрешающие напряжения в системах скольжения, в настоящее время измеренные для многих материалов. Дальнейшие упрощения возможны, если ограничиться учетом лишь одной системы скольжения. Этот случай достаточно реалистичен. Например, в монокристалле чистого титана имеются три плоскости скольжения [28]: призматическая, базовая и пирамидальная, критические напряжения для которых соответственно равны 60, 120 и 120-180 МПа [29]. Так что при малых деформациях активируется только одна призматическая система скольжения. Кроме того, в поликристаллическом титане до макродеформаций 3 % при одноосном растяжении упрочнение практически отсутствует [30]. Наличие одной активной системы скольжения характерно для многих кристаллов с гексагональной плотно упакованной решеткой [31]. Ограничение учетом одной системы скольжения в таких анизотропных материалах представляется адекватным при рассмотрении критических явлений типа зарождения пластических сдвигов в отдельных зернах, имеющих место при малых макронапряжениях и макродеформациях поликристаллов. В одной системе скольжения g при отсутствии упрочнения разрешающее сдвиговое напряжение по мере роста деформации линейно увеличивается до критического значения и далее остается постоянным: , при . (77) Условие скольжения (77) принимает вид , (78) откуда (79) и . (80) Закон Гука принимает вид (81) Учтем, что в (81) произведение есть скаляр (зависящий от деформации), а также скаляром является постоянный знаменатель, который обозначим через : . (82) Поэтому (81) можно переписать в виде (83) Отсюда получаем выражения для тензора секущих модулей в рассматриваемом случае: . (84) Выражение (84) справедливо для деформаций , в противном случае . 3. Особенности упругопластического взаимодействия зерен и экстремальные микроструктуры Упругопластическое взаимодействие деформаций разных зерен описывается коэффициентами тензора . Для каждой пары зерен тензор имеет 36 компонентов. Свойства симметрии по перестановкам такие же, как для случая упругого деформирования [17]: тензор симметричен по перестановке индексов первой и второй пары, но несимметричен по перестановке пар индексов, что делает взаимодействие несимметричным относительно взаимодействующих зерен. Величина коэффициентов взаимодействия зависит от расстояния между зернами и ориентации их кристаллографических осей. Вычисления для различных типов зерен показывают, что зависимость от расстояния и ориентации осей зерен в упругопластическом случае подобна такой же зависимости для упругого взаимодействия [17]. Величина коэффициентов сильно зависит от упругой анизотропии зерен и достаточно медленно убывает с расстоянием между зернами. Коэффициенты упругопластических межзеренных взаимодействий зависят от мезодеформаций в зернах, в отличие от чисто упругих взаимодействий, которые не зависят от деформированного состояния зерен [17]. Эту зависимость исследуем на примере двух соседних зерен титана кубической формы. Зерна расположим вдоль оси координат Х3. Оси гексагональной упругой симметрии обоих зерен направим вдоль оси Х3. В обоих зернах активны только призматические системы скольжения с вектором нормали плоскости скольжения и вектором направления скольжения в выбранной системе координат. Зерна подвергнуты деформации простого сдвига в плоскости Х1-Х2, так что единственная ненулевая компонента деформаций в зернах . Для вычислений используем зависимости секущих модулей зерен от мезодеформаций по формулам (84), соответствующим одной активной системе скольжения. По вычисленным секущим модулям для заданной мезодеформации численно находятся коэффициенты по формулам (35), а затем по формулам (38). В таблице приведены вычисленные значения характерных компонент тензоров коэффициентов межзеренного взаимодействия - нормальных с нормальными ( ) и сдвиговых со сдвиговыми ( ) деформациями в зависимости от полных деформаций в зернах. Рассмотрены четыре возможных комбинации взаимодействий: упруго-упругое, упруго-упругопластическое, упругопластическое-упругое и упругопластическое-упругопластическое взаимодействия. Зависимость коэффициентов взаимодействия от полной деформации в зернах Interaction coefficients dependence on full strain in grains 0 0,1 % 0,2 % 0,3 % 0,4 % 0,5 % Упруго-упругое -3,596 -3,596 -3,596 -3,596 -3,596 -3,596 6,754 6,754 6,754 6,754 6,754 6,754 Упруго-упругопласт. -3,596 -3,663 -3,670 -3,628 -3,624 -3,610 6,754 10,708 22,159 25,915 27,855 28,955 Упругопласт.-упругое -3,596 -3,589 -3,645 -3,703 -3,658 -3,686 6,754 7,372 9,561 10,440 10,956 11,219 Упругопласт.-упругопласт. -3,596 -3,595 -3,598 -3,626 -3,646 -3,656 6,754 11,723 31,442 40,127 44,988 47,849 Из таблицы видно, что интенсивность взаимодействия двух пластически деформированных зерен кратно выше, чем тех же двух упругодеформированных (при полной деформации 0,5 % более чем в 7 раз). Взаимодействие двух по-разному деформированных зерен несимметрично. Действие пластически деформированного зерна на упругое в среднем в 2,5 раза интенсивнее действия упругодеформированного зерна на пластически деформированное. Это означает, что пластически деформированное зерно при росте макродеформации оказывает возрастающее влияние на соседние упругие зерна, способствуя развитию в них пластических деформаций. Таким образом, локально возникшая в одном зерне пластическая деформация быстро распространяется на соседние зерна, способствуя делокализации очагов пластического деформирования и гомогенизации мезодеформаций. Очаги пластических деформаций возникают в специфических кластерах зерен, в которых флуктуирующие поля деформаций принимают экстремальные значения. Разработанный формализм позволяет находить такие кластеры в случайной микроструктуре поликристалла при различных типах макронагружения, одноосных и многоосных, и исследовать эволюцию деформированного состояния зерен в кластере при увеличении макродеформации всего поликристалла. Для вычислительного поиска нужно сгенерировать большой ансамбль случайных микроструктур поликристалла, решить для каждой реализации микроструктуры из ансамбля краевую задачу деформирования и выбрать микроструктуру с экстремальными деформациями в каком-либо выбранном зерне. Рассмотрим кластер зерна размером 3х3х3, погруженный в осредненную среду остальных зерен поликристаллического титана (рисунок, слева), и найдем конфигурацию микроструктуры, приводящую к максимальной деформации в центральном зерне кластера. Хотя в таком кластере взаимодействуют только ближайшие зерна-соседи (для кубических зерен одинакового размера число ближайших соседей центрального зерна равно 26), то все равно ясно, что нужно сгенерировать неимоверно большое число случайных конфигураций микроструктур кластера и решить для каждой такое же число краевых задач, из решений которых выбрать экстремальное. Аддитивность предлагаемого подхода позволяет сократить объем необходимой вычислительной работы на много порядков. Достаточно найти экстремальную ориентацию каждого зерна-соседа по отдельности. Решить для него «парную» краевую задачу, а затем сложить результаты решений «парных» задач для всех соседей. Рассмотрим для примера еще раз макроскопический сдвиг поликристаллического титана в плоскости Х1Х2 с учетом только одной призматической системы скольжения и без упрочнения. Без учета межзеренных взаимодействий максимальные разрешающие сдвиговые напряжения генерируются в зернах, в которых плоскость пирамидального скольжения параллельна плоскости Х1Х3. Это зерна предпочтительной в отношении пластического деформирования ориентации, мягкие зерна [6, 8]. В нашем подходе этому соответствует нулевое решение по межзеренному взаимодействию - решение первого уравнения в (34). Фиксируем в качестве такого «мягкого зерна» центральное зерно кластера. X3 X2 X1 Рис. Ориентации зерен в экстремальном кластере Fig. Grains orientations in extreme cluster Взаимодействие с соседями увеличивает деформацию в мягком зерне. Поправки от взаимодействия даются правыми частями уравнений в (34). Величины добавок зависят от ориентации кристаллографических осей взаимодействующих зерен. Разработанный алгоритм для заданной макродеформации находит ориентации соседних зерен, дающие максимальные вклады, на большом массиве ориентаций и визуализирует эти ориентации. Найденная структура экстремального кластера, приводящая к максимальному значению сдвиговой деформации в центральном зерне при макродеформации сдвига , показана в правой части рисунка. Кластер образует симметричный паттерн. Для наглядности зерна разделены на три слоя, параллельные плоскости сдвига, и условно разнесены. Обоюдоострые стрелки изображают оси гексагональной упругой симметрии зерен, а их наклон имитирует углы, образуемые этими осями с осью Х3 системы координат в экстремальной конфигурации. Проследим эволюцию деформирования двух исследуемых взаимодействующих зерен - центрального зерна и зерна в зените кластера. Вычисления по разработанному алгоритму решения систем уравнений (34) для мезодеформаций в первом порядке теории возмущений (чему соответствуют первые два уравнения в системе) с учетом взаимодействий всех зерен в кластере дают следующие результаты. До макродеформации все зерна кластера деформируются упруго. При этой макродеформации мезодеформация в центральном зерне кластера (на 33 % больше макродеформации), а разрешающее сдвиговое напряжение в призматической системе скольжения сравнивается с критическим: Мезодеформация зерна в зените , а разрешающее сдвиговое напряжение в нем меньше критического. В остальных зернах мезодеформации значительно меньше (кроме зерна в надире, эволюционирующего так же, как зерно в зените). Мезодеформации в рассматриваемых зернах отличаются в раз. При макродеформации начинается пластическое скольжение в центральном зерне. Полная мезодеформация в нем увеличивается до Полная мезодеформация в зенитном зерне увеличивается до Различие деформаций в паре зерен уменьшается до раз. При макродеформации полные мезодеформации в рассматриваемых зернах становятся соответственно: и Межзеренная неоднородность деформаций становится еще меньше: Таким образом, упругопластическое взаимодействие зерен приводит к существенному уравниванию деформаций в рассмотренной паре зерен. С такой же скоростью происходит выравнивание деформаций в зерне в надире кластера. Деформации в остальных зернах кластера меньше, но пластическое деформирование начинается и в них и ведет к снижению неоднородностей деформаций. В целом происходит гомогенизация деформаций в зернах при увеличении макродеформации, качественно описанная выше, даже при относительно малых макродеформациях в несколько десятых долей процента. При больших макродеформациях активизируются другие системы скольжения, гомогенизация усиливается и распространяется на остальные зерна кластера. Эффект гомогенизации проявляется и на уровне одного зерна, если его разбить на малые субзерна и применить описанную в работе [32] методику вычисления неоднородных мезодеформаций. Все субзерна в зерне имеют одинаковую ориентацию кристаллографических осей и систем скольжения. Вследствие этого коэффициенты взаимодействия субзерен в одном зерне выше, чем межзеренные коэффициенты, и эффект гомогенизации проявляется сильнее. Экспериментальные измерения деформаций в индивидуальных зернах с помощью современных источников синхротронного рентгеновского излучения, методов микрорентгенотомографии и дифракции нейтронов [33, 34] подтверждают высокую степень однородности пластических мезодеформаций. Заметная неоднородность наблюдается только вблизи стыков границ зерен. Заключение Развитый ранее авторами теоретико-полевой подход к проблеме упругого деформирования поликристаллических материалов распространен на упругопластические задачи. В постановке задачи используются определяющие соотношения деформационного типа. Дифференциальная форма краевой задачи упругопластического деформирования трансформируется в систему интегральных уравнений для мезодеформаций в зернах. Посредством кусочно-постоянной аппроксимации мезодеформаций система интегральных уравнений преобразуется в систему нелинейных алгебраических уравнений. Последние решаются методом теории возмущений по межзеренному взаимодействию и итераций по нелинейностям. Деформация в каждой точке поликристалла представляется в виде вкладов от парных взаимодействий со всеми зернами. Вследствие аддитивности этих вкладов решение краевой задачи для поликристаллического тела, содержащего макроскопически большое число зерен, разбивается на сумму решений двухзеренных задач. Численное решение парных задач значительно проще, чем для сплошного поликристаллического тела, составленного из множества зерен. Последнее обстоятельство позволяет кардинально, на много порядков, сократить вычислительную трудоемкость решения некоторых задач, таких как поиск экстремальных микроструктур, в которых генерируются экстремально большие мезодеформации. На примере поликристаллического титана исследованы особенности упругопластического взаимодействия зерен. Показано, что его интенсивность зависит от деформированного состояния зерен и существенно выше интенсивности упругого взаимодействия. При полных мезодеформациях в паре зерен 0,5 % это различие составляет 7 раз. Рост интенсивности упругопластического взаимодействия зерен при увеличении макродеформации поликристаллического тела приводит к гомогенизации мезодеформаций, что может быть одной из причин успеха простых моделей типа Тэйлора-Линя в описании пластического деформирования поликристаллов, полностью игнорирующих различие мезодеформаций в зернах.

Об авторах

А А Ташкинов

Пермский национальный исследовательский политехнический университет

В Е Шавшуков

Пермский национальный исследовательский политехнический университет

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

Статистика

Просмотры

Аннотация - 330

PDF (Russian) - 246