Автоволновая механика пластичности металлов
- Авторы: Зуев ЛБ1, Баранникова СА1
- Учреждения:
- Институт физики прочности и материаловедения СО РАН
- Выпуск: № 1 (2019)
- Страницы: 50-64
- Раздел: Статьи
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/449
- DOI: https://doi.org/10.15593/perm.mech/2019.1.05
- Цитировать
Аннотация
Предложена модель развития локализованного пластического течения твердых тел, основанная на представлениях о взаимодействии носителей пластичности и сигналов акустической эмиссии, возникающих при развитии элементарных актов пластичности. Экспериментально показано, что локализация пластического течения является общим признаком всех деформационных процессов и может наблюдаться на всем протяжении процесса от предела текучести до формирования шейки. Специфические картины макроскопической локализации деформации - паттерны локализованного пластического течения - являются автоволнами, которые генерируются в деформируемой среде за счет упорядочения ее дефектной структуры. Паттерны локализованной пластичности коррелируют со стадиями деформационного упрочнения, выявляемыми при механических испытаниях металлов и сплавов. Предложена и проанализирована двухкомпонентная модель развития локализованной деформации, учитывающая взаимодействие упругой и пластической компонент деформации. В рамках этой модели описана генерация автоволн локализованной пластической деформации в ряде материалов. Установлены общие закономерности развития локализованного пластического течения на разных стадиях деформационного упрочнения. Введен упругопластический инвариант деформации, устанавливающий взаимосвязь между упругими и пластическими свойствами материалов. Показано, что основные свойства деформируемых сред являются следствиями такого инварианта. В частности, следствиями из упругопластического инварианта деформации являются: зависимость скорости распространения автоволн локализованной пластической деформации от коэффициента деформационного упрочнения, дисперсионное уравнение для автоволн, масштабный эффект локализации пластической деформации, зависимость автоволновых параметров от размера зерна, уравнение автоволны локализованной деформации, коэффициент деформационного упрочнения, уравнение дислокационной динамики Тейлора-Орована, соотношение Холла-Петча, зависимость плотности подвижных дислокаций от деформации.
Ключевые слова
Полный текст
Введение Многолетнее изучение природы пластической деформации твердых тел привело к заключению, что пространственно-временнáя неоднородность пластического течения наблюдается на любом его этапе и является наиболее важным атрибутом деформационного процесса. Справедливость этого заключения отчетливо прослеживается на всех масштабных уровнях деформационного процесса [1], начиная с дислокационного и заканчивая макроскопическим. При этом дислокационные механизмы развития локализованной пластичности изучены очень детально, в основном методом трансмиссионной электронной микроскопии тонких фольг [2]. Макроскопические аспекты локализации исследовались в значительно меньшей степени. В наших работах, результаты которых обобщены в монографиях [3, 4], макромасштабные аспекты пластического течения изучались в моно- и поликристаллах чистых металлов и сплавов разного состава с ГЦК, ОЦК, ГПУ и тетрагональной решетками, а также в неметаллических материалах (керамика, щелочно-галоидные монокристаллы, горные породы). В ходе этих экспериментов было показано, что пластическая деформация твердых тел имеет макроскопически локализованный характер на протяжении всего процесса пластического течения, т.е. от предела текучести до разрушения. Возникающие при этом в образцах пространственно-временные структуры локализации макроскопического масштаба спонтанно генерируются при деформации с постоянной скоростью, так что деформируемая среда в ходе процесса пластического течения самопроизвольно расслаивается на чередующиеся друг с другом деформирующиеся и недеформирующиеся в данный момент времени объемы, как это показано на рис. 1. На наблюдаемой поверхности образца деформирующиеся и недеформирующиеся объемы образуют специфическую картину распределения очагов локализованной пластичности. В нелинейных средах разной природы самопроизвольно возникают аналогичные картины, которые в современной научной литературе называются паттернами [5-7]. Этим термином далее будут обозначаться картины локализованной деформации при их интерпретации, качественном и количественном анализе. Паттерны локализации пластической деформации удается визуализировать с помощью методов спекл-фотографии, развитых специально для подобного анализа [3]. Это дает возможность определить пространственный (длину автоволны) l и временнóй T периоды процесса с помощью построения диаграмм «Положение деформационного очага - время» (X-t диаграммы), как это показано на рис. 2, а также вычислить волновое число , частоту и скорость движения очагов . Рассматривая проблему пластического течения, Зегер и Франк [5] указали на важную роль формирования специфических дефектных структур, возникающих в деформируемой среде на макроскопическом масштабном уровне. С общей точки зрения, предложенной Хакеном, структурообразование является процессом самоорганизации среды, а «система называется самоорганизующейся, если она без специфического воздействия извне обретает какую-то пространственную, временнýю или функциональную структуру» [6]. Рис. 1. Последовательные (1→2→3) паттерны локализации пластичности в монокристалле сплава Fe-Ni-Cr на стадии линейного деформационного упрочнения. Темным выделены активные зоны деформации Fig. 1. Sequential (1 → 2 → 3) patterns of localized plasticity in a single crystal of the alloy Fe-Ni-Cr at the stage of linear work hardening. The dark color corresponds to active deformation zones а б Рис. 2. К методике определения пространственного и временнóго периодов локализованной пластической деформации с помощью диаграмм Fe-Ni-Cr монокристаллов, ориентированных вдоль направлений (а) и (б) на стадии линейного деформационного упрочнения Fig.2. To the method determining the spatial and temporal periods of the localized plastic deformation using diagrams; Fe - Ni - Cr single crystals oriented along the axis (a) and (b) at the stage of linear work hardening Здесь обсуждаются основные закономерности рождения и развития паттернов локализации деформации на разных стадиях деформационного процесса при растяжении с постоянной скоростью. Такой подход принципиально отличается от традиционного анализа процессов пластической деформации, в рамках которого максимальное внимание уделяется описанию событий на дислокационном уровне, но практически игнорируется роль кристаллической решетки, которой обычно отводится всего лишь роль резервуара кристаллических дефектов, непосредственно отвечающих за пластическое течение. Работа призвана продемонстрировать существование связи между решеточными характеристиками и закономерностями локализованного пластического течения и подчеркнуть особую роль этой связи в развитии локализованного на макроскопическом масштабном уровне пластического течения. 1. Локализация пластического течения и автоволновая природа пластической деформации Возникающие при пластическом течении деформируемой среды паттерны локализованной деформации имеют характерный масштаб (радиус корреляции) ~10-2 м. Всесторонний анализ эволюции паттернов в ходе деформационных процессов, проведенный на большом числе деформируемых материалов с разными структурами и механическими свойствами, позволил установить следующие общие закономерности процесса развития локализованной пластической деформации [3, 4, 8, 9]: - паттерны макролокализации возникают самопроизвольно при деформации с постоянной скоростью или в режиме ползучести, не требуя для своего появления специального воздействия; - паттерны локализации закономерно меняются по мере развития пластического течения: между их формами и стадийностью пластического течения прослеживается тесная связь; - на некоторых этапах развития деформации паттерны локализации демонстрируют отчетливо выраженную пространственно-временнýю периодичность. Интерпретация физической природы паттернов локализации основана на гипотезе о том, что они являются проекциями автоволновых процессов [4], возникающих в объеме деформируемой среды, на наблюдаемую поверхность образца. Практически одновременно автор работы [10] сформулировал гипотезу о волновом характере пластического течения. Эта точка зрения, очевидно, была определена представлениями о волнах пластичности (волнах Кольского), связанных с распространением фронтов пластической деформации [11, 12]. Эта скорость [11, 12]. В то же время, как было экспериментально установлено в [3], скорость распространения автоволн зависит от коэффициента деформационного упрочнения , как . Различие зависимостей и порождает проблему выбора между возможностями использовать для описания локализованной пластичности волновые или автоволновые процессы. Такой выбор должен также учитывать существенное различие скоростей распространения: 10-4 м/с [3, 4], а 102 м/с [11, 12]. Указанная проблема может быть решена на основе следующих рассуждений. Можно ожидать, что генерации волновых и автоволновых процессов обусловливают разное поведение термодинамических характеристик деформируемых систем, в частности энтропии. Для проверки этого предположения примем, что скорость пластической деформации , и применим для анализа кинетическое уравнение Тейлора-Орована для дислокационной деформации [2] , где b -вектор Бюргерса; - плотность подвижных дислокаций. Скорость термически активированного движения дислокаций определяется соотношением [2] (1) в котором - термодинамический потенциал Гиббса; U - внутренняя энергия; S - энтропия; - работа деформирующего напряжения s; g - активационный объем процесса; - постоянная Больцмана, а T - температура. Из-за близости механических характеристик исследованных материалов при грубой оценке можно считать, что . Тогда из уравнения (1) следует , т.е. зависимости для волновых и автоволновых процессов, представленные на рис. 3 в координатах , качественно соответствуют зависимостям для них. При построении этого графика данные о скорости волн локализации деформации на стадиях легкого скольжения монокристаллов (прямая ) и стадиях линейного деформационного упрочнения моно- и поликристаллов металлов и сплавов (прямая ) определены экспериментально; скорости распространения волн пластичности (прямая ) для этих же металлов и сплавов вычислены с помощью выражения [12] с использованием рассчитанных по кривым нагружения значений и справочных данных о плотности исследуемых материалов. Рис. 3. Изменение энтропии при формировании волн пластичности и автоволн локализации пластического течения на стадиях легкого скольжения и линейного деформационного упрочнения Fig. 3. Entropy change during the formation of plasticity waves and autowaves of localized plastic flow at the stages of easy glide and linear work hardening Анализируя графики зависимостей для обоих типов волн, представленные на рис. 3, можно видеть, что в случае автоволн локализованной пластичности (прямые 1 и 2) энтропия системы линейно убывает с ростом коэффициента деформационного упрочнения, как на стадии легкого скольжения, так и на стадии линейного деформационного упрочнения, т.е. для этих случаев . Естественнно, это справедливо для термодинамически открытой системы, какой является деформируемый образец. Напротив, при генерации волн пластичности (прямая 3) энтропия системы возрастает ( ), что характерно для любых диссипативных процессов. Таким образом, данные рис. 3 демонстрируют принципиальное различие между обнаруженными автоволнами локализованного пластического течения и другими волновыми процессами, возможными при пластическом деформировании твердых тел. Снижение энтропии деформируемой системы при генерации автоволн пластичности показывает, что она является самоорганизующейся [6, 13], т.е. точка зрения на паттерны локализованной пластичности как на автоволны является предпочтительной. Развивая этот тезис, обсудим некоторые важные свойства деформируемой среды. Прежде всего она является нелинейной, на что указывает сложный характер зависимости «напряжение-деформация». Далее, деформируемая среда в процессе деформации необратимо меняет свою структуру и свойства, запоминая механические воздействия. Наконец, деформируемая среда является активной [6, 7, 14], т.е. содержит распределенные по объему источники энергии. Последними служат концентраторы напряжений с их полями упругих напряжений [1, 3, 4]. Согласно принятым представлениям о природе самоорганизации [14] в открытых системах, такая среда способна генерировать автоволновые процессы. 2. Уравнения автоволновых процессов пластического течения Закономерности явлений локализации пластического течения показывают, что паттерны локализации определяются автоволновыми процессами в активных средах. Автоволны описываются решениями параболических дифференциальных уравнений типа [14], которые получаются добавлением нелинейной функции в правую часть уравнения Фурье . В то же время упругие волны описываются гиперболическими дифференциальными уравнениями типа , в которых скорость распространения определяется материальными константами среды. Известно [4, 14], что адекватное описание генерации автоволн в активных средах базируется на учете кинетики двух антагонистических факторов - активирующего (автокаталитического) и ингибирующего (демпфирующего). При пластическом течении в качестве активатора процесса выбирается деформация e, а в качестве ингибитора - напряжение s [3]. Скорости изменения и в таком случае определяются дифференциальными уравнениями (2) (3) Нелинейные функции и суть точечные кинетики деформаций и напряжений соответственно, а члены со вторыми пространственными производными описывают перераспределение деформаций и напряжений в объеме образца за счет случайных процессов диффузионного типа, причем и - коэффициенты с размерностью L2∙T-1. Смысл их двойной индексации будет объяснен ниже. Фактически при таком подходе уравнения (2) и (3) постулируются, однако они могут быть получены строгим образом. Так, уравнение кинетики пластической деформации (активатора) следует из условия неразрывности пластического течения [15]: (4) где - поток деформации в поле ее градиента. Если в общем случае , то (5) где (6) - нелинейная функция деформации и напряжения. Уравнение кинетики релаксации напряжений (ингибитора) следует из уравнения Эйлера для потока вязкой среды П, записанного в форме [16] (7) В вязкой среде - тензор плотности потока импульса; - единичный тензор; p - давление, а и - компоненты скорости потока. Тензор напряжений есть сумма упругих и вязких напряжений. При пластической деформации или Скорость релаксации упругих напряжений [17] , где M - упругий модуль системы «образец-испытательная машина»; 10-5…10-4 Па·с - коэффициент квазивязкого торможения дислокаций, а - скорость их надбарьерного движения. Вязкие напряжения возникают из-за неоднородности поля внутренних напряжений и линейно связаны с изменениями скорости распространения упругих волн в деформируемой среде: [3], где - скорость распространения поперечных упругих волн в отсутствие напряжений, а . Если принять, что , где - динамическая вязкость среды, то . Тогда скорость релаксации вязких напряжений Окончательно (8) где - коэффициент переноса с размерностью L2∙T-1. Очевидно, уравнения (6) и (8) по форме совпадают с постулированными выше уравнениями (2) и (3). Правая часть уравнения (8) может рассматриваться как сумма скоростей релаксации упругих и вязких напряжений. Нелинейная функция g(ε,σ) учитывает перераспределение упругих напряжений между контактирующими микрообъемами на фронте уже существующего очага пластичности, а член ответствен за стохастические процессы перераспределения напряжений в образце на макроскопическом масштабе. Таким образом, деформационный поток оказывается состоящим из «гидродинамической» и «диффузионной» компонент. Первая из них описывается нелинейными функциями и в уравнениях (6) и (8) и связана с непрерывным движением вдоль образца фронта деформации при последовательной активации локальных концентраторов напряжений на нем. Вторая компонента определяется членами и в этих же уравнениях и инициирует пластическую деформацию на макроскопическом расстоянии от существующего фронта за счет релаксации одного из концентраторов в этой области и связанного с этим рождения дислокаций. Подобный эффект «заброски деформации» является причиной порождения автоволны локализованного пластического течения. 3. Локализация и стадийность пластического течения Если принять во внимание, что характерной и хорошо изученной закономерностью любого процесса пластического течения является его многостадийность [1, 2], то паттерны локализации целесообразно сопоставить с действующими на соответствующих стадиях течения законами деформационного упрочнения. Стадии упрочнения на кривой пластического течения формально могут быть выделены с помощью уравнения Людвика [1, 15]. Здесь ; K - модуль упрочнения, а n - показатель упрочнения. Отыскивая участки диаграммы течения , для которых , можно выделять стадии деформационного упрочнения: площадку текучести, стадии легкого скольжения, линейного и параболического деформационного упрочнения и предразрушения. Сопоставление стадий процесса с данными о паттернах локализации пластического течения показало, что каждой стадии пластического течения отвечает определенный паттерн, практически не зависящий от особенностей структуры конкретного материала. Это наблюдение, справедливое для всех исследованных до настоящего момента материалов, позволило установить полное однозначное соответствие форм локализации и стадий деформационного упрочнения. Так, площадке текучести соответствует движение вдоль образца уединенного очага локализованной деформации, на стадии линейного деформационного упрочнения наблюдается движение по образцу с постоянной скоростью группы эквидистантных очагов локализации, на параболической стадии деформационного упрочнения в образце возникает система эквидистантно расположенных неподвижных очагов локализованной деформации, и, наконец, стадии предразрушения соответствует согласованное движение очагов локализованной пластичности, приводящее к их «коллапсу» и последующему вязкому разрушению. Этот порядок схематически представлен на рис. 4. Рис. 4. диаграммы для разных стадий деформационного упрочнения Fig. 4. diagrams for different stages of work hardening Эти закономерности выполняются для всех материалов независимо от их состава, структуры и механизма пластической деформации. Даже достаточно широкие вариации этих характеристик влекут за собой лишь количественные изменения паттернов локализации, не меняя их формы. Можно утверждать, что наблюдаемые паттерны локализации исчерпываются перечисленными выше формами, число которых всегда совпадает с числом наблюдаемых в материале стадий кривой пластического течения. Сравнение наблюдаемых закономерностей с описанными в литературе структурами автоволновых процессов [14] показывает, что на каждой стадии пластического течения возникает специфический тип паттерна локализации, которому в соответствие может быть поставлена определенная автоволновая мода. Это позволяет утверждать, что в рамках обсуждаемой модели пластическое течение следует рассматривать как эволюцию автоволновых мод локализации, протекающую начиная от предела текучести, вплоть до разрушения. Иначе говоря, процесс пластической деформации, включая переход от пластического течения к вязкому разрушению, есть последовательная и закономерная смена типов автоволн локализованной деформации, идущая в следующем порядке: автоволна переключения → фазовая автоволна → стационарная диссипативная структура → коллапс автоволны. В результате такого анализа паттерны локализации пластического течения, характерные для каждой стадии деформационного упрочнения, отождествляются с определенными модами (типами) автоволновых процессов [14]. Отсюда следует Правило соответствия, которое представлено в форме табл. 1 (см. также рис. 4). Согласно Правилу соответствия, многостадийный процесс пластического течения может рассматриваться как эволюция автоволновой картины деформации. При наблюдении это проявляется как закономерная перестройка наблюдаемых паттернов локализованной пластичности. Такая точка зрения, согласно которой процесс пластического течения характеризуется пространственно-временной периодичностью, получила развитие, например, в работах [18-24]. Остановимся еще на одной важной особенности эволюции автоволновых процессов в деформируемой среде. Известно [10], что экспериментальное изучение различных типов автоволн, например в химических, а тем более в биологических системах, требует создания генератора, индивидуального для каждого варианта процесса. Такие генераторы могут отличаться друг от друга типом или кинетикой химических реакций в них, химическим составом, температурным режимом, размерами и другими параметрами. Таблица 1 Правило соответствия стадий деформационного упрочнения и автоволновых мод Table 1 The rule of conformity of the work hardening stages and autowave modes Стадия Обозначение на рис. 4 Характеристика Мода автоволны Площадка текучести I Автоволна переключения Легкое скольжение и линейное деформационное упрочнение II Фазовая автоволна Параболическое упрочнение III Стационарная диссипативная структура Предразрушение IV Коллапс автоволны а б в г Рис. 5. Характеристики фазовой автоволны локализованного пластического течения: зависимость скорости распространения от коэффициента деформационного упрочнения (а); дисперсионная зависимость (б); зависимость длины автоволны от размера зерна (в); зависимость длины автоволны от длины образца (г) Fig. 5. Characteristics of the phase autowave of the localized plastic flow: the dependence of the propagation velocity on the work hardening coefficient (a); dispersion dependence (b); dependence of the autowave length on the grain size (c); dependence of autowave length on the sample length (g) Однако, как следует из изложенного выше, в твердых телах генерация всех типов автоволн локализованной пластичности осуществляется при растяжении с постоянной скоростью [4]. По этой причине процесс генерации различных автоволновых мод может быть реализован при значительно меньших экспериментальных усилиях. По существу, деформируемое твердое тело может рассматриваться как универсальный генератор автоволновых процессов, удобный для их моделирования и экспериментального изучения [4]. Формирующаяся на стадиях легкого скольжения (в монокристаллах) и линейного деформационного упрочнения фазовая автоволна соответствует условию постоянства фазы . Для этого типа автоволн возможно определение длины и скорость распространения. В наших работах [3, 4] были установлены зависимости этих характеристик от условий эксперимента и структуры материала. Так, для скорости распространения фазовых автоволн имеем (рис. 5, а) (9) Закон дисперсии фазовых автоволн имеет квадратичный вид (рис. 5, б): (10) Зависимость длины автоволны от размера зерна задается логистическим соотношением (функцией Ферхюльста, рис. 5, в) (11) а длина автоволны логарифмически зависит от длины образца L (рис. 5, г): (12) Соотношения (9)-(12) и соответствующие графики, показанные на рис. 5, в, достаточно полно характеризуют процессы развития локализованного пластического течения, однако набор констант C и в них требует специального объяснения, которое будет приведено ниже. 4. Инвариант упругой и пластической деформации Признание того факта, что пластическое течение твердых тел представляет собой автоволновой процесс, приводит к очевидной аналогии с упругой деформацией, в основе развития которой лежат процессы распространения упругих волн. При сравнении количественных характеристик упругих волн и автоволн локализованного пластического течения была обнаружена тесная линейная корреляция между парными произведениями макропараметров и , характерных для стадии линейного деформационного упрочнения, и характеристик кристаллической решетки соответствующих металлов: 10-10 м - расстояние между плотно упакованными плоскостями и 103 м/с - скорость распространения упругих поперечных волн. Отношение этих произведений, вычисленное усреднением по 18 металлам, приведенным в табл. 2, имеет вид (13) Соотношение (13) было названо упругопластическим инвариантом деформации. Таблица 2 Сравнение величин и для металлов Table 2 Comparison of quantities and for metals ×107м2/с Стадия линейного деформационного упрочнения Cu Zn Al Zr Ti V Nb α-Fe γ-Fe Ni Co Mo 3,6 3,7 7,9 3,7 2,5 2,8 1,8 2,55 2,2 2,1 3,0 1,2 4,8 11,9 7,5 11,9 7,9 6,2 5,3 4,7 6,5 6,0 6,0 7,4 0,75 0,3 1,1 0,3 0,3 0,45 0,33 0,54 0,34 0,35 0,5 0,2 ×107м2/с Стадия линейного деформационного упрочнения Стадия легкого скольжения Sn Mg Cd In Pb Ta Hf α-Fe γ-Fe Cu Zn Ni Sn 2,4 9,9 0,9 2,6 3,2 1,1 1,0 7,4 2,9 1,9 1,0 1,3 3,3 5,3 15,8 3,5 2,2 2,0 4,7 4,2 6,5 6,0 4,7 5,0 6,0 4,9 0,65 0,63 0,2 1,2 1,6 0,2 0,24 1,1 0,49 0,4 0,2 0,2 0,67 Очевидно, инвариант (13) играет очень важную роль в описании пластичности, количественно связывая характеристики упругих волн с характеристиками автоволн локализации пластического течения. Заданная соотношением (13) закономерность отражает взаимосвязь малых упругих и больших пластических деформаций, а наличие скоростей в инварианте (13) указывает на связь контролирующих деформацию процессов - перераспределения упругих напряжений с характерной скоростью и перераспределения областей локализованной пластичности со скоростью , причем 107…108. Существование инварианта (13) указывает на некорректность традиционного допущения, принимаемого по умолчанию практически в любом подходе к проблеме описания пластичности. Оно состоит в том, что в силу аддитивности вкладов упругой и пластической деформаций в общую при очевидном условии упругим вкладом допустимо пренебречь и полагать . Инвариант (13), однако, указывает на взаимообусловленность упругой и пластической деформаций в процессе пластического формоизменения. Для объяснения природы упругопластического инварианта примем во внимание развиваемую в последние годы точку зрения [3, 4], согласно которой локализация пластической деформации является результатом самоорганизации (структурообразования) в активной нелинейной деформируемой среде. Главным признаком развития процессов самоорганизации в термодинамически открытой системе, какой является деформируемая среда, служит уменьшение ее энтропии [13]. Как было показано выше, такое условие реализуется при формировании автоволн локализованного пластического течения. Поэтому использование энтропии для объяснения природы процессов локализации пластической деформации обоснованно и перспективно. Будем рассматривать пластическое течение как самосогласованную пространственно-временнýю трансформацию полей напряжений и пластических деформаций [3, 4]. При такой трансформации релаксация напряжений вызывает деформацию, а изменение деформации инициирует перестройку поля напряжений. Кинетика этих процессов контролируется скоростями, входящими в инвариант (13): для упругого поля и для поля пластических деформаций. Полагая теперь, что в основе трансформации полей лежат смещения частиц среды, рассмотрим связь упругих (обратимых) и пластических (необратимых) смещений при малом отклонении деформируемой системы от равновесия. В этом случае скорости смещений при трансформациях полей деформаций и напряжений с точностью до малых первого порядка могут быть приняты линейными по градиентам пластических и упругих деформаций и записаны в виде (14 a) (14 б) где из соображений размерности принято, что и . Из-за существенной нелинейности связи деформации и напряжения, следующей из формы кривой течения , необходимо учесть также возникновение дополнительных скоростей (15 a) (15 б) Окончательно система уравнений для пластической и упругой составляющей скоростей смещений приобретает вид (16 a) (16 б) Коэффициенты уравнений этой системы можно представить в виде матрицы [23] (17) в которой, согласно принципу симметрии кинетических коэффициентов Онзагера [7, 18, 23], недиагональные компоненты равны, т.е. Диагональные коэффициенты и , которые служат коэффициентами автоволновых уравнений локализованной пластичности (2) и (3), не обязаны быть равными; так, в [3], например, рассмотрен случай, когда В уравнении (13) длины и суть пространственные масштабы трансформации полей упругой и пластической деформации, а скорости и характеризуют кинетику такой трансформации. Уравнение (13) в виде (18) сводит инвариант (13) к отношению масштабного и кинетического факторов. Первый интерпретируется как число возможных мест зарождения автоволны локализованной пластической деформации, а второй характеризует выбор скорости автоволны из спектра возможных скоростей в твердом теле . Тогда отношения и имеют смысл термодинамических вероятностей [4]. В таком случае возможна численная оценка по уравнению (18) путем расчета изменения энтропии системы при самопроизвольном формировании автоволн локализованного пластического течения. В силу аддитивности энтропии запишем ее изменение при генерации автоволны как сумму масштабного и кинетического вкладов: (19) Условие отражает уменьшение энтропии при формировании автоволн локализации пластического течения. Для его выполнения необходимо, чтобы хотя бы одно слагаемое в уравнении (19) было отрицательным. Используя формулу Больцмана, получим для масштабного вклада (20) так как Здесь постоянная Больцмана. Полагая, как сказано выше, кинетический вклад отрицательным, получаем (21) Знаки величин и в уравнениях (20) и (21) указывают на различие вкладов масштабного и кинетического факторов в характер развития локализованной пластической деформации. Из уравнений (19)-(21) следует (22) и, соответственно, (23) Окончательно (24) В этом случае (25) в расчете на элементарный релаксационный акт [24]. 5. Некоторые следствия из упругопластического инварианта Упругопластический инвариант деформации (13) имеет ряд важных следствий, охватывающих, в частности, некоторые ранее экспериментально установленные закономерности. В частности, перечисленные выше зависимости (9)-(12) могут быть выведены из уравнения упругопластического инварианта (13). Например, дифференцируя уравнение (13) по деформации e, получаем (26) Записав это уравнение относительно , имеем (27) где величина c не зависит от пластической деформации, так что . Тогда (28) Простые преобразования уравнения (28) приводят к уравнению (29) совпадающему с ранее полученной экспериментальным путем зависимостью (9), если положить, что коэффициент деформационного упрочнения есть отношение двух структурных параметров λ и [4], т.е. . Чтобы получить дисперсионное уравнение, запишем уравнение (13) в форме (30) где . Если , то В этом случае (31) Таким образом, если принять, что , то закон дисперсии автоволн локализованного пластического течения приобретает квадратичную форму: (32) которая ранее была установлена экспериментально [4] и соответствует уравнению (10). Запишем далее уравнение инварианта (13) в форме (33) и учтем, что скорости и зависят от размера зерна [4]. В таком случае дифференцирование соотношения (33) по дает (34) Преобразование уравнения (34) приводит к (35) где и , поскольку Решение уравнения (35) есть логистическая функция (36) где , а C - константа интегрирования. Зависимость (36) экспериментально получена в [3] для размеров зерна 5·10-3 15 мм в алюминии и описывается уравнением (11). Наконец, продифференцируем выражение (13) по длине образца L: (37) В этом случае (38) При условии приходим к , т.е. к зависимости (12) , которая также была установлена ранее экспериментально [4]. В [4] были рассмотрены в общей сложности более десяти следствий из упругопластического инварианта (13). Это подчеркивает важную роль обсуждаемого соотношения в развитии представлений о природе пластичности материалов. 6. Упругопластический инвариант как основное уравнение автоволновой механики пластичности Очень важно, что автоволновое уравнение для деформации (2) тесно связано с инвариантом (13). Чтобы показать это, запишем (39) и примем, что - пластическая деформация. Применяя оператор , в котором коэффициент имеет размерность L2·T-1, к левой и правой частям этого уравнения, соответственно получаем (40) Скорость распространения ультразвука слабо зависит от деформации [4], так что Кроме того, в силу (13) . В этом случае (41) что эквивалентно уравнению (6) для скорости деформации, если . Развивая этот результат, можно установить связь автоволнового подхода к проблеме пластичности с теорией дислокаций [8, 17]. Известно, что многие дислокационные теории пластичности базируются на уравнении Тейлора-Орована, определяющего скорость пластической деформации: (42) Сравним теперь уравнения (41) и (42) и покажем, что первый член в правой части уравнения (41) аналогичен выражению в уравнении (42). Действительно, пусть ( - частота Дебая), а . Примем также, что . Тогда если дислокации распределены хаотически, то , где l - пробег дислокаций. Следовательно, (43) Скорость дислокаций связана со скоростью , так что . Окончательно (44) Подобная связь автоволнового уравнения с уравнением дислокационной кинетики была установлена нами ранее другим путем [4]. Очевидно, что уравнение (44) отличается от уравнения Тейлора-Орована (42) только наличием члена в правой части. Это означает, что уравнение Тейлора-Орована является частным случаем уравнения (44), а само уравнение (44) есть расширенный вариант уравнения дислокационной кинетики, который включает как «гидродинамическую» (зависящую от скорости) , так и «диффузионно-подобную» компоненты деформационного потока. Таким образом, автоволновая модель пластического течения и теория дислокаций не противоречат друг другу, а являются взаимно дополнительными. При малых плотностях хаотически распределенных дислокаций использование уравнения (42) допускает получение правильных результатов. Однако при больших пластических деформациях и соответствующих им больших плотностях дефектов теория дислокаций становится нелинейной и возникает необходимость использования уравнения (44). 7. Двухкомпонентная модель пластического течения В рамках развиваемой модели находится объяснение генерации автоволн локализованной пластической деформации. Действительно, если снова использовать уравнение Тейлора-Орована (42), то задаваемое испытательной машиной условие может быть удовлетворено только при , т.е. при достаточной плотности подвижных дислокаций и достаточной скорости их движения. Это условие может не выполняться из-за снижения плотности подвижных дислокаций при росте деформации [27] или за счет падения скорости движения подвижных дислокаций при уменьшении эффективного напряжения, действующего на дислокацию, от до [2] при деформационном упрочнении. В таком случае условие обеспечивается включением в работу диффузионно-подобного механизма, который контролируется членом в уравнении (44) и порождает новые очаги локализованного пластического течения на расстоянии от исходного. Это является причиной «заброски деформации» и приводит к формированию автоволны локализованного пластического течения. Развиваемая модель локализованного пластического течения базируется на инварианте (13), который косвенно указывает на важную роль акустических свойств (фононной подсистемы) кристалла в формировании автоволновых картин локализации пластического течения. Этому сразу же может быть дано объяснение, связанное с тем, что в скорость макроскопической пластической деформации делают вклад только движущиеся между локальными препятствиями дислокации, для которых . Их движение в таком случае является надбарьерным и контролируется торможением фононным и электронным газами [28]. Состояние этих газов отражается на динамике пластического течения, в том числе и на макромасштабном уровне. В развиваемой модели, однако, фононной подсистеме отводится более важная роль в описании развития локализованного пластического течения твердых тел. Процесс деформации включает два типа взаимосвязанных событий (двухкомпонентная модель), согласованно протекающих в деформируемой среде и описывающихся уравнениями (2) и (3). С одной стороны, это скачкообразные релаксационные акты движения дислокаций, их ансамблей и в макроскопическом пределе - автоволн локализованного пластического течения, сопровождающиеся акустической эмиссией [29], т.е. генерацией упругих волн при перераспределении упругих деформаций в ходе релаксации. С другой стороны, известен так называемый акустопластический эффект, состоящий в инициировании пластических сдвигов при действии ультразвуковых импульсов [30]. Эти эффекты изучались ранее независимо друг от друга, но развиваемая модель построена на учете их взаимосвязи. Смысл рассматриваемой модели состоит в объединении этих двух явлений. Если элементарный релаксационный сдвиг порождает импульс акустической эмиссии, то последний, в свою очередь, может инициировать развитие нового сдвига в другой области кристалла. Таким образом, перераспределение существующего упругого поля концентраторов, в свою очередь, инициирует новые сдвиги в подсистеме движущихся дислокаций. Этот подход согласуется с идеей, объясняющей самоорганизацию в сложных системах спонтанным расслоением самоорганизующихся систем на взаимосвязанные информационную и динамическую подсистемы [31]. Конкретный выбор таких подсистем для случая деформируемой среды несложен. В активной деформируемой среде, которая моделируется набором концентраторов напряжений, первую можно связать с импульсами акустической эмиссии при дислокационных сдвигах, а вторую - с собственно сдвигами. Сценарий развития пластического течения, схематически показанный на рис. 6, в этом случае таков. Релаксация концентраторов напряжений генерирует импульсы акустической эмиссии и формирует сигнальную подсистему. Энергия акустических сигналов, поглощаемая другими концентраторами, инициирует их распад, вызывая новые сдвиги в динамической подсистеме образца, которые сопровождаются новыми импульсами акустической эмиссии. Рис. 6. Схема работы двухкомпонентной модели развития автоволновой деформации Fig. 6. The scheme of the two-component model of the autowave deformation development Оценим возможности предложенной модели, сравнив времена ожидания t термоактивированного акта сдвиговой релаксации [32] без акустического импульса 5·10-5 c (45) и при его действии (46) При вычислениях по формулам (45) и (46) принято, что 0,5 эВ, а эВ. При действии на концентратор акустического импульса с амплитудой упругой деформации энергия активации снижается на 0,1 эВ. Сравнение величин и , даже при очевидной грубости оценки, подтверждает справедливость обсуждаемой модели. Развиваемая модель основана на взаимосвязи явлений в акустической и дислокационной подсистемах и способна объяснить крупномасштабную периодичность расположения очагов локализованной пластичности в деформируемой среде, исходно содержащей только дислокации с их характерным масштабом ~b. Для этой цели достаточно оценить расстояние, на котором акт новой релаксации будет иметь наибольшую вероятность. Для этого используем описанный в [32] эффект расщепления поперечной ультразвуковой волны в поле упругих напряжений. Пусть при элементарном сдвиге эмитируется импульс поперечных упругих волн с частотой w. В упругонапряженной области он расщепляется на две взаимно перпендикулярно поляризованные волны, распространяющиеся со скоростями и и имеющие длины волн и . Для разности в этом случае имеем [33] (47) Для расчета по формуле (47) примем, что разность главных нормальных напряжений в напряженной области » 108 Пa, r » 5×103 кг/м3, » 3×103 м/с, а отвечающая максимуму энергии в спектре акустического излучения при пластической деформации металлов частота w » 106 Гц [29]. При таких условиях »10-4 м. Вероятность нового сдвига максимальна там, где совпадают максимумы квадратов напряжений в обеих волнах, т.е. максимальна упругая энергия волны. Это соответствует 10-2 м, что порядку величины близко к наблюдаемой длине автоволны и объясняет дистанцию «заброски деформации», т.е. зарождения очагов пластичности на расстоянии ~l от существующего фронта деформации за счет процессов «диффузионного» типа, за которые ответственны члены и в уравнениях (2) и (3) [34]. Заключение Экспериментальные данные показали, что пластическое течение твердых тел макроскопически локализовано на всем протяжении кривой пластического течения, причем существует ограниченное количество форм локализации, а каждая из них определяется действующим законом деформационного упрочнения. Закономерности формоизменения определяются сменой форм макролокализации, т.е. последовательной генерацией нескольких типов автоволн локализованного пластического течения. На основе автоволновых представлений предложена двухкомпонентная модель развития локализованной пластичности, согласно которой пластическое течение является результатом взаимосогласованной активации дислокационных сдвигов и импульсов акустической эмиссии. Характеристики упругой и пластической деформации в деформируемой среде связаны упругопластическим инвариантом, существование которого определяется связью процессов перестройки полей упругой и пластической деформации. Следствия из инварианта позволяют в достаточной мере точно описывать основные закономерности локализованного на макромасштабном уровне пластического течения.Об авторах
Л Б Зуев
Институт физики прочности и материаловедения СО РАН
С А Баранникова
Институт физики прочности и материаловедения СО РАН
Список литературы
- Pelleg J. Mechanical Properties of Materials. - Dordrecht: Springer, 2013. - 634 p.
- Messerschmidt U. Dislocation Dynamics during Plastic Deformation. - Berlin: Springer, 2010. - 503 p.
- Зуев Л.Б., Данилов В.И., Баранникова С.А. Физика макролокализации пластического течения. - Новосибирск: Наука, 2008. - 327 с.
- Зуев Л.Б. Автоволновая пластичность. Локализация и коллективные моды. - М.: Физматлит, 2018. - 207 с.
- Seeger A., Frank W. Structure formation by dissipative processes in crystals with high defect densities // Non-Linear Phenomena in Materials Science. Eds. L.P. Kubin and G. Martin. - New York: Trans Tech Publications, 1987. - P. 125-138.
- Хакен Г. Информация и самоорганизация. Макроскопический подход к сложным системам. - М.: URSS, 2014. - 317 с.
- Скотт Э. Нелинейная наука. Рождение и развитие когерентных структур. - М.: Физматлит, 2007. - 559 с.
- Зуев Л.Б., Баранникова С.А., Лунев А.Г. Упругопластический инвариант деформации в металлах // Успехи физики металлов. - 2018. - Т. 18, № 4. - С. 379-481.
- Зуев Л.Б., Баранникова С.А., Семухин Б.С. Кинетика развития паттернов макролокализации пластического течения металлов // ФТТ. - 2018. - № 7. - С. 1358-1364.
- Hähner P. Theory of solitary plastic waves // Appl. Phys. - 1994. - Vol. A58. - No. 1. - Р. 41-58.
- Кольский Г. Волны напряжений в твердых телах. - М.: ИЛ, 1955. - 192 c.
- Шестопалов Л.М. Деформация металлов и волны пластичности в них. - М.-Л.: Изд-во АН СССР, 1958. - 268 с.
- Климонтович Ю.Л. Введение в физику открытых систем. - М.: Янус-К, 2002. - 284 с.
- Васильев В.А., Романовский Ю.М., Яхно В.Г. Автоволновые процессы. - М.: Наука, 1987. - 240 с.
- Хилл Р. Математическая теория пластичности. - М.: ГИТТЛ, 1956. - 407 с.
- Ландау Л.Д., Лифшиц Е.М. Гидродинамика. - М.: Наука, 1986. - 736 с.
- Пустовалов В.В. Скачкообразная деформация металлов и сплавов при низких температурах // Физ. низк. темп. - 2008. - Т. 34, № 9. - С. 871-913.
- Упругопластический переход в железе: структурные и термодинамические особенности / О.А. Плехов, О.Б. Наймарк, N. Saintier, T. Palin-Luc // ЖТФ. - 2009. - Т. 7, № 8. - С. 56-61.
- Третьякова Т.В., Вильдеман В.Э. Пространственно-временная неоднородность процессов неупругого деформирования металлов. - М.: Физматлит, 2016. - 118 с.
- Природа упругопластического инварианта деформации / Л.Б. Зуев, В.И. Данилов, С.А. Баранникова, Н.А. Плосков // ЖТФ. - 2018. - Т. 88, № 6. - С. 858-862.
- Reyne B., Manach P.-Y., Moes N. Macroscpoic consequences of Poibert-Luders and Portevin-Le Chatelier bands during tensile Deformation in Al-Mg alloys // Mat. Sci. Eng. A. - 2019. - Vol. 746. - P. 187-196.
- Kobelev N.P., Lebyodkin M.A., Lebedkina T.A. Role of self-organization of dislocations in the onset and kinetics of macroscopic plastic instability // Met. Mat. Trans. A: Phys. Met. Mat. Sci. - 2017. - Vol. 48(3). - P. 965-974.
- Taupin V., Chevy J., Fressengeas C. Effects of grain-to-grain interactions on shear strain localization in Al-Cu-Li rolled sheets // Int. J. Sol. Str. - 2016. - Vol. 99. - P. 71-81.
- Tretyakova T., Wildemann V. Study of spatial-time inhomogeneity of inelastic deformation and failure in bodies with concentrators by using the digital image correlation and infrared analysis // Proc. Str. Integ. - 2017. - Vol. 5. - P. 318-324.
- Румер Ю.Б., Рывкин М.Ш. Термодинамика, статистическая физика и кинетика. - М.: Наука, 1977. - 552 с.
- Слуцкер А.И., Гиляров В.Л., Лукьяненко А.С. Особенности энергетики адиабатически нагружаемого ангармонического осциллятора // ФТТ. - 2006. - Т. 48, № 10. - С. 1832-1837.
- Gilman J.J. Micromechanics of plastic flow at a constant stress // J. Appl. Phys. - 1965. - Vol. 36. - No. 9. - Р. 2772-2777.
- Al’shits V.I., Indenbom V.L. Mechanism of dislocation drag // Dislocations in Solids. - Vol. 7. - Amsterdam: Elsevier, 1986. - Р. 43-111.
- Williams R.V. Acoustic Emission. - Bristol: Adam Hilger, 1980. - 412 p.
- Malygin G.A. Acoustic-plastic effect and mechanism of stress superposition // Solids State Phys. - 2000. - Vol. 42. - No. 1. - P. 69-75.
- Кадомцев Б.Б. Динамика и информация. - М.: Ред. УФН, 1997. - 399 с.
- Caillard D., Martin J.L. Thermally Activated Mechanisms in Crystal Plasticity. - Oxford: Elsevier, 2003. - 452 p.
- Tokuoka T., Iwashimizu Yu. Acoustical birefringence of ultrasound waves in deformed isotropic elastic materials // Int. J. Solids Structures. - 1968. - Vol. 4. - No. 2. - Р. 383-389.
- Основные соотношения автоволновой модели пластического течения / Л.Б. Зуев, С.А. Баранникова, А.Г. Лунев, С.В. Колосов, А.М. Жармухамбетова // Известия вузов. Физика. - 2018. - Т. 61, № 9. - С. 142-148.