Full Text

Введение Неустойчивое пластическое течение (прерывистая текучесть), которое проявляется на диаграммах нагружения в виде ступенек или зубцов различного типа, можно наблюдать во многих сплавах, подвергнутых растяжению (Мак Рейнольдс, 1949), сжатию (Bell, 1973) или кручению (Диллон, 1963) [3] в определенном диапазоне скоростей деформаций, напряжений и температур. Такую пластическую неустойчивость обычно связывают с эффектом Портевена-Ле Шателье (ПЛШ). Проблема появления прерывистой текучести имеет большое практическое значение, поскольку эффект ПЛШ играет отрицательную роль при обработке промышленных сплавов, ухудшая механические свойства изделий и пластичность материалов. Деформирование в режиме прерывистой пластичности, особенно на заключительных стадиях технологических процессов обработки металлов давлением, ведет к существенному ухудшению качества поверхности (шероховатости), что приводит к снижению усталостной прочности, коррозионной стойкости и аэродинамических характеристик изделий. В связи с этим экспериментальным исследованиям данного эффекта посвящено большое число работ как отечественных, так и зарубежных авторов. В первой части обзора [24] был представлен краткий обзор работ, посвященных описанию физических механизмов и экспериментальных исследований прерывистой пластичности, рассмотрены наиболее известные макрофеноменологические модели. Следует отметить, что свойства материала на макроуровне в значительной степени определяются его микроструктурой, которая в большинстве макрофеноменологических моделей не учитывается. Кроме того, сложность физической интерпретации соотношений и параметров, абстрагирование от вопросов эволюции микроструктуры делают макрофеноменологические модели малопригодными для разработки новых материалов и предсказания их свойств. В настоящей статье анализируются некоторые теоретические работы, основанные на физических теориях пластичности, посвященные описанию особенностей деформирования сплавов в температурно-скоростных диапазонах, в которых существенное влияние на поведение материалов оказывают диффузионные процессы. Особое внимание уделено описанию эффекта ПЛШ, возникновение которого большинство авторов связывают с взаимодействием дислокаций с атмосферами примесных атомов. Из-за возросшего в последнее время интереса к эффекту Портевена-Ле Шателье, одной из главных задач в изучении которого является определение критерия потери устойчивости пластического деформирования и установление области существования данного типа неустойчивости, в литературе представлено большое число моделей различных типов для его описания, в частности - феноменологические [71, 97], модели, основанные на физических теориях пластичности [35, 66, 77, 87, 89, 90, 95], модели дислокационной динамики [15, 35, 36, 37]; следует отметить также обзорные работы, содержащие описания различных моделей [57, 82, 93, 96] и физических механизмов, ответственных за появление данного эффекта [9, 10]. 1. Эффект Портевена-Ле Шателье Эффект Портевена-Ле Шателье как проявление неустойчивости пластического течения можно наблюдать для широкого класса сплавов при произвольных видах нагружения в определенных диапазонах скоростей деформации и температур. В экспериментах на одноосное нагружение эффект ПЛШ проявляется в следующем: при низких скоростях деформирования и повышенной температуре диаграмма σ-ε приобретает пилообразную форму («зубчики») при «жестком» нагружении; при «мягком» нагружении диаграмма становится ступенчатой [2, 3]. В ряде работ, посвященных исследованию эффекта ПЛШ, проявление последнего связывают с возникновением полос сдвига, наблюдаемых на боковой поверхности плоских образцов, которые представляют собой области локализованной деформации [17, 28, 30, 59, 60, 70, 95]. Выделяют три основных типа проявления эффекта Портевена-Ле Шателье: 1) тип А - появление и движение вдоль оси образца одиночной (уединенной) деформационной полосы, которое может происходить многократно; 2) тип В - деформационные полосы появляются и исчезают в осциллирующем или перемежающемся режиме, распространяясь вдоль образца (stop-and-go); 3) тип С - полосы возникают (и пропадают) случайным образом по длине образца. Приведенная классификация является условной; в реальных опытах может наблюдаться различное сочетание этих трех типов [4, 12]. Закономерности локализации пластической деформации в образце при его нагружении, как отмечается в ряде работ [5, 6, 13, 27, 28], часто коррелируют со скачкообразным характером деформационной кривой, при этом каждый скачок напряжения на диаграмме σ-ε связывают с возникновением или распространением деформационной полосы по боковой поверхности образца. В работах ряда авторов предлагаются методы подавления эффекта ПЛШ [7, 27, 29, 84]. Сопоставление экспериментальных данных позволяет сделать вывод, что эффект ПЛШ проявляется при отрицательной скоростной чувствительности напряжения течения m [35, 47, 50, 59, 60, 64, 71]. В условиях одноосного нагружения при заданной температуре скоростную чувствительность можно определить следующим соотношением: где при ступенчатом увеличении скорости деформации от до величина напряжения изменяется от до . Краткий обзор известных теоретических и экспериментальных работ описания неустойчивого пластического течения содержится в [82]. Все модели описания эффекта ПЛШ можно разделить на две основные группы: (I) модели, основанные на макроскопическом описании деформационных полос, предложенные Kubin и Estrin [71] и затем модифицированные другими исследователями [97]; (II) модели, основанные на микроскопическом описании динамического деформационного старения с использованием внутренней переменной - времени старения, предложенные MacCormick [79, 80] и модифицированные в работах [59, 83]. В работе [82] представлена модификация модели MacCormick [80] для никелевого сплава, учитывающая динамическое деформационное старение (ДДС). Расчеты проведены для одноосного растяжения при малых деформациях. Параметры для жаропрочного никелевого сплава были получены из испытаний на растяжение при постоянной скорости деформации. В работе [70] для воспроизведения распространяющихся полос деформации используется прямая упруговязкопластическая модель поликристалла [25]. Численные эксперименты на растяжение при постоянной скорости деформации 10-3 с-1 для Al-Mg сплава были проведены для образцов прямоугольного сечения 0,5×2×12 мм (расчеты проведены с шагом по времени t = 0,005 с). В результате моделирования авторам удалось воспроизвести деформационные полосы, наблюдаемые в эксперименте, которые распространялись непрерывно от одного конца образца к другому, где они либо отражались, либо начинали распространяться заново, что характерно для полос деформации типа А. В работе [81] представлены результаты экспериментальных и численных исследований эффекта ПЛШ при простом сдвиге в Al-Mg сплаве при различных скоростях деформации и температурах. Прямые наблюдения поверхности образца с использованием метода корреляции цифровых изображений позволили выявить различные типы полос деформации в зависимости от напряжений и скорости деформации. Для описания поведение алюминиевых сплавов, в которых зачастую при деформировании могут проявляться различные типы неустойчивости пластического течения, в частности эффект Портевена-Ле-Шателье, используются конститутивные модели MacCormick и Johnson-Cookс использованием анизотропного и изотропного закона упрочнения соответственно. 2. Влияние состава и микроструктуры материала на характер деформационных кривых При численном исследовании процессов деформации важную роль играет корректное описание свойств материала деформируемого образца в широком диапазоне величин пластических деформаций, скоростей деформаций и температур. В большинстве случаев уравнения состояния в математических моделях процессов деформирования формируются на основании тщательно проведённых экспериментов и учитывают некоторое количество эмпирических констант. Известно, что многие материалы демонстрируют сходные закономерности, многообразие свойств чистых металлов и сплавов зависит не только от физико-химических взаимодействий и микроструктуры, возникающих на стадии создания материала (например, кристаллизацией из расплавов), но и от процессов, происходящих на стадиях последующей термомеханической обработки заготовок. Большинство используемых в промышленности сплавов образуются на основе некоторого основного металла (например, железа - в сталях), остальные компоненты могут составлять от десятых (и сотых) долей до десятков процентов атомной массы; в дальнейшем атомы последних будут называться «примесными». Для рассматриваемой в предлагаемой работе прерывистой текучести к классу наиболее важных процессов относятся диффузионные процессы, изучению и описанию которых посвящено огромное число теоретических и экспериментальных работ. В последнее время в литературе появляется большое количество статей, в которых представлены исследования влияния на параметры скачкообразной деформации микроструктуры сплавов и концентрационных характеристик примесных атомов [5, 11, 18, 21, 31, 78, 84, 94, 101]. В данном разделе приведены некоторые теоретические и экспериментальные работы, в которых отражено влияние микроструктуры и состава материала на вид деформационных кривых. Пластическая деформация твердых тел обусловлена главным образом движением краевых дислокаций, которые при своем движении преодолевают различные препятствия: точечные (примесные атомы), линейные (барьеры Ломера-Коттрелла и другие), поверхностные (границы зерен, двойников и т.д.) и объемные (выделения вторичных фаз и другие). В общем случае дислокации располагаются и движутся в разных плоскостях скольжения, в том числе и пересекающихся. Подвижность дислокаций в неидеальных кристаллах уменьшается за счет их взаимодействия друг с другом и с другими дефектами. Такое взаимодействие приводит к торможению движущихся дислокаций; к иммобилизованным дислокациям могут диффундировать атомы примесей. Влияние точечных дефектов на свойства материалов весьма велико, атомы примеси образуют вокруг дислокаций зону повышенной концентрации, что затрудняет движение последних [14, 19, 21, 22, 26, 45, 46, 69]. При этом дислокации взаимодействуют с атомами примеси упругими полями напряжений, примесные атомы испытывают со стороны дислокаций силы притяжения или отталкивания. С одной стороны от плоскости скольжения расположена область гидростатического (всестороннего) сжатия, а с другой - гидростатического растяжения. Чем больше фактор размерного несоответствия атомов растворенного элемента и основы, тем сильнее эти атомы тормозят дислокации и тем больше энергия упругого взаимодействия дислокации с примесным атомом. Атомы внедрения притягиваются к области гидростатического растяжения и располагаются не в совершенной области решетки, где создают поля значительных напряжений, а под краем экстраплоскости, где им энергетически выгоднее размещаться. Притяжение атомов примесей, вызванное разными причинами, приводит к «осаждению» этих атомов вдоль края экстраплоскости, такую область повышенной концентрации примесных атомов внедрения вдоль линии дислокации называют атмосферой Коттрелла. Атомы примеси в атмосфере Коттрелла «привязаны» к дислокации силами упругого взаимодействия. Дислокация при скольжении стремится увлечь за собой атмосферу Коттрелла, которая, в отличие от самой дислокации, может перемещаться на расстояния больше межатомного только диффузией, т.е. сравнительно медленно. Поэтому атмосфера из растворенных атомов способна перемещаться вместе с дислокацией, находящейся в центре этой атмосферы, лишь при высоких температурах и очень малых скоростях скольжения дислокации. При увеличении скорости скольжения атмосфера отстает от ядра дислокации, дислокация «освобождается» от дополнительного тормозящего её движение воздействия, происходит падение критического сдвигового напряжения. Аналогичный эффект возникает одновременно для значительного количества дислокаций, что на макроуровне находит отражение в резком уменьшении напряжения течения. В работе [31]приводятся результаты исследования сплавов АМг3 и АМг6, в которых существенное влияние на механическую прочность и характеристики скачкообразной деформации играет микроструктура. Обнаружены границы существования скачкообразной деформации, переход от зубчатой формы кривой σ-ε к гладкой наблюдается при добавлении 0,5% Si в сплавах АМг3, переход от монотонных зависимостей к скачкообразным при нагружении образцов холоднокатаного сплава АМг6 происходит с ростом температуры отжига образцов. Авторы предполагают, что наличие в границах зерен примесных частиц Mg или Al подавляют макроскопическую скачкообразную деформацию, влияют на формирование атмосфер Коттрелла. Авторы статьи [26] предлагают способ повышения прочности металлических материалов на примере Al сплавов. Высокая плотность дислокаций в границах зерен и отсутствие дислокаций в их внутренних областях в ультрамелкозернистых материалах повышает роль границ зерен как источников и стоков дислокаций. В результате создаются благоприятные условия для стекания примесных атомов в область границ зерен. Дислокации оказываются погруженными в твердый раствор с повышенной концентрацией примесных атомов. Показано, что при варьировании фактора размерного несоответствия и подбором примесных элементов с соответствующими радиусами атомов и их концентрации в Al сплавах, подвергнутых интенсивной пластической деформации, можно повысить вклад напряжения закрепления дислокаций в предел текучести ультрамелкозернистых сплавов. В [84] представлены результаты исследования влияния микроструктуры образцов на проявления эффекта Портвена-Ле Шателье в Al-Mg-Sc сплаве с нерекристаллизованной, частично перекристаллизованной и полностью рекристаллизованной зеренной структурой. Авторами было установлено, что измельчение зерна способствует возникновению неустойчивости пластического течения, появление первого скачка критического напряжения происходит при более низких скоростях деформирования, амплитуда скачков, наблюдаемых при комнатных температурах, возрастает в несколько раз, что указывает на значительное увеличение вклада в упрочнение за счет взаимодействия атомов растворенного вещества с дислокациями. При этом измельчение зерна уменьшает локализацию деформации и снижает шероховатость поверхности. На микроскопическом уровне эффект ПЛШ, как предполагается в работах [8, 20, 47, 48, 55, 56, 59, 60, 65, 71, 81], связан с динамическим деформационным старением, обусловленным взаимодействием между подвижными дислокациями, временно остановившимися на препятствиях (дислокациях леса или других дефектах), и диффундирующими атомами примеси, которые дополнительно закрепляют дислокации. MacCormick повышенную диффузию атомов растворенного вещества объяснял увеличением плотности дислокаций и увеличением концентрации вакансий в процессе деформирования. В работах [58, 59, 70, 72] предполагается, что существенную роль в ДДС играет диффузия атомов примеси вдоль мобильных дислокаций и дислокаций леса, а не объемная диффузия и вакансионный перенос вещества. Авторы работ [53, 100] также объясняют высокую диффузию примесных атомов при деформировании сплавов туннельной диффузией. Данный подход предсказывает очень высокие значения концентрации примесных атомов в области гидростатического растяжения и очень низкие - в области гидростатического сжатия. Для описания неустойчивости пластического течения, связанного с эффектом Портевена-Ле Шателье, Balik и Lukac [45] предложили модель трения, в которой примесные атомы внедрения располагаются вдоль линии дислокации, образуя атмосферы, и распространяются за счет диффузии вдоль линии дислокации. В работе анализируется зависимость напряжения течения от скорости движения дислокации, температуры и концентрации растворенного вещества. Показано, что динамическое деформационное старение и отрицательная чувствительность к скорости деформации зависят от наличия примеси. Были получены зависимости напряжения от температуры для трех различных концентраций примеси и четырех значений скорости дислокаций, показано, что вид расчетных кривых соответствует экспериментальным зависимостям. В работе [47] отмечается влияние динамического деформационного старения на эффект ПЛШ. Авторы попытались воспроизвести данный тип неустойчивости при одноосном нагружении, учитывая только динамическое старение, ведущее к отрицательной скоростной чувствительности. Для теоретического описания используется модель, которая позволяет отделить скачкообразное и монотонное поведение на кривой деформации. Модель включает в себя систему дифференциальных уравнений в частных производных и однородные граничные условия, для которой в случае одноосного нагружения получено аналитическое решение. Проявление эффекта ПЛШ в условиях ДДС зависит от ряда факторов: температуры, средней скорости деформации, наличия растворенных веществ и барьеров дислокационной и недислокационной природы. Основой существующих моделей для описания ДДС является понимание того, как атомы растворенного вещества в сплаве взаимодействуют с дислокациями в процессе деформации. Если времена диффузии атомов примеси при некоторых скоростях деформации и температурах соизмеримы с временами активации движения дислокаций, то зависимость σ-ε может приобрести зубчатую форму. Увеличение скорости деформации (при данной температуре) и, соответственно, уменьшение времени остановки (задержки)дислокаций на препятствиях может привести к снижению ДДС, что, в свою очередь, облегчает термическую активацию дислокаций и может привести к отрицательной скоростной чувствительности напряжения течения. В работе [53] рассмотрены некоторые существующие модели для описания механизма ДДС и их основные недостатки. Показано, что механизм ДДС и отрицательная чувствительность к скорости деформирования могут быть объяснены туннельной диффузией атомов растворенного вещества непосредственно через плоскость скольжения от области гидростатического сжатия в область растяжения. Данный механизм диффузии вызывает большую разницу энтальпии между атомами растворенного вещества по обе стороны от плоскости скольжения, что приводит к дополнительному упрочнению при взаимодействий дислокаций и примесных атомов. В работе используется кинетическая модель для описания движения атомов Mg в Al-Mg сплаве, предлагаются эволюционные уравнения для концентрации атомов, решение производится методом Монте-Карло. Представлены карты распределения атомов Mg при 500 К и исходной концентрации Mg - 0,05 %. Результаты, полученные с использованием модели, находятся в хорошем соответствии с экспериментальными данными. Для объяснения влияния ДДС на эффект Портевена-Ле Шателье Zhang и Curtin предложили модель [100], которая описывает сопротивление движению дислокаций примесными атомами при деформировании Al-Mg сплавов, в основе которой лежит механизм туннельной диффузии. Авторы изучали процессы диффузии в Al-Mg сплавах при различных температурах и концентрациях примеси, в особенности - роль решеточной и туннельной диффузии. Дифференциальные уравнения диффузии используются для получения стационарного распределения облаков примесных атомов в зависимости от скорости дислокаций, температуры и концентрации Mg. Моделирование показало, что при скоростях дислокаций, превышающих некоторое критическое значение, соответствующее максимальному напряжению торможения дислокаций, концентрация Mg резко уменьшается в области гидростатического сжатия, атомы притягиваются к области гидростатического растяжения, находящейся под краем экстраплоскости, где их расположение дает выигрыш в энергии. Это приводит к асимметричному распределению атомов Mg вокруг дислокаций. При скоростях дислокаций, значительно превышающих критическую, сопротивление движению дислокаций уменьшается, и облака примесных атомов уже не могут внести существенный вклад в напряжение торможения дислокаций. Результаты численного моделирования показывают хорошее соответствие экспериментальным данным при деформировании в широком диапазоне концентраций и температур; отметим, они согласуются с экспериментальными данными, приведенными в более поздней работе [32]. Влияние примесей на закономерности пластической деформации кристаллов с позиций эффекта динамического старения рассмотрено в статье [21]. Представленная в работе модель состоит из двух подмоделей: первая описывает кинетику иммобилизации отдельных дислокаций за счет увеличения примесей, вторая описывает коллективную динамику дислокационного ансамбля. Объединенная модель позволяет рассчитывать концентрационную и температурную зависимости вклада примесных атомов в деформирующее напряжение; используемая теория объясняет аномальную температурную зависимость предела текучести, экспериментально наблюдаемую в ряде материалов. В обзорной статье [57] рассматриваются различные механизмы взаимодействия между атомами растворенного вещества и движущимися дислокациями. Предложенные для их описания теории зачастую работают в узком диапазоне, в котором динамическое деформационное старение оказывает значительное явление на процессы, происходящие в материале, что ограничивает область применения данных моделей. Во многих из предложенных моделей приводятся только общие формулировки и качественные характеристики, а некоторые детали, например специфические формы уравнений или значения материальных констант, не всегда описываются в публикациях. Таким образом, можно констатировать, что наиболее распространенные модели, опирающиеся на рассмотрение влияния динамического деформационного старения основаны на классических представлениях Коттрелла [52], согласно которым возникновение эффекта ПЛШ объясняют диффузией примесных атомов к дислокациям и «закрепление» последних за счет взаимодействия с атомами примеси. 3. Модели, основанные на физических теориях вязкопластичности В последние 15-20 лет в связи с широким распространением для исследования поведения материалов при неупругом деформировании многоуровневых моделей, базирующихся на различных физических теориях пластичности и вязкопластичности, растет количество работ, в которых указанный класс моделей применяется для анализа прерывистой пластичности. Для многих материалов, проявляющих прерывистую текучесть, наблюдаются два качественно различных режима зависимости критического напряжения от скорости деформации и температуры эксперимента, которые принято обозначать как нормальное и обратное (аномальное или инверсное) поведение. Нормальное прерывистое течение ПЛШ возникает, когда критическая пластическая деформация (деформация появления первого скачка напряжения) возрастает с ростом скорости деформирования и уменьшением температуры. При уменьшении критической деформации с ростом скорости деформирования и уменьшением температуры испытания возникает обратная зависимость, которую называют аномальным течением. Исследованию эффекта Портевена-Ле Шателье как с экспериментальной, так и с теоретической точки зрения посвящена работа [50]. В работе предлагается модель, которая качественно объясняет экспериментально наблюдаемую обратную скоростную и аномальную температурную зависимость напряжения течения в пересыщенных твердых растворах. Рассмотрено влияние гомогенных и гетерогенных включений вторичных фаз (осадков) на характеристики эффекта ПЛШ в концентрированных твердых растворах, условия возникновения прерывистой текучести, проведено сравнение с экспериментальными данными. Выявлено, что в концентрированных твердых растворах критическая деформация, необходимая для возникновения эффекта Портвена-Ле Шателье, должна испытывать «обратную» температурную зависимость по сравнению с разбавленными твердыми растворами. Механизм возникновения аномального эффекта Портвена-Ле Шателье, т.е. увеличения (а не уменьшения с ростом температуры и понижением скорости деформации величины критической деформации появления первых скачков нагрузки на диаграммах деформации, обсуждается в [15]. В работе используется квазистатическая модель динамического деформационного старения, которая учитывает, что движение дислокаций в процессе пластической деформации происходит с остановками их на барьерах, которые дислокации преодолевают с помощью тепловых флуктуаций. Нормальный эффект ПЛШ объясняется блокировкой дислокаций межузельными атомами, обратный - закреплением дислокации атомами замещения. Нормальное поведение связывают с возникновением полос типа А и В, напротив, полосы типа С появляются и распространяются при обратном эффекте Портевена-Ле Шателье при блокировке дислокаций атомами замещения). В работе [49] приведены одномерная и трехмерная версии модели, учитывающие ДДС, которые адекватно воспроизводят как нормальный, так и обратный эффект Портевена-Ле Шателье. Для одномерной модели по аналогии с работой [85] предложено аддитивное разложение напряжения течения. Упрочнение определяется по закону Войса. При расчетах предполагается, что время, необходимое для активации движения дислокаций, равно времени ожидания на локальных препятствиях. Предлагаемая нелинейная упруговязкопластическая модель в трехмерной постановке была реализована в пакете ABAQUS. Данные по одноосному растяжению образцов из алюминиевого сплава 2024 были использованы для идентификации модели. Результаты численных расчетов показали качественное и количественное соответствие экспериментальным данным. Модель для определения сопротивления деформации, представленная в виде системы интегродифференциальных уравнений, описывающих в совокупности упрочнение материала в результате приращения плотности дислокаций и блокирования движения свободных дислокаций границами зерен и субзерен, учитывающая разупрочнение из-за возврата и рекристаллизации, предложена в [91]. Физические процессы упрочнения и разупрочнения происходят при изменении формы деформационной кривой и влияют на физико-механические свойства сплавов. В сплавах с большим количеством включений может наблюдаться барьерный эффект, проявляющийся в остановке подвижных свободных дислокаций на границах зерен, субзерен, примесных атомах и включениях. В модель сопротивления деформации были введены дополнительные внутренние параметры, отвечающие за блокирование движения подвижных дислокаций примесными атомами и включениями. Модель позволяет описать сопротивление деформации сплавов при высокотемпературной деформации с учетом комбинированного влияния упрочнения, восстановления, рекристаллизации и барьерных эффектов (блокировки свободных дислокаций примесными атомами), а также эффект Портвена-Ле Шателье. Параметры модели определялись в процессе ее идентификации по экспериментальным данным, полученным при одноосном растяжении образцов сплава АМг6в диапазоне скоростей деформации от 5 до 25 с-1 и температурах 400 и 500 °С. Предложенная в более ранней работе данных авторов [23] модель, описывающая сопротивление деформации металломатричного композита 15 % SiC-Alпри высокотемпературной пластическойдеформации (при температурах 450 и 500 ºС), для диапазона скоростей деформаций от 0,1 до 10 с-1 с достаточно хорошей точностьюописывает поведение материала. Известно, что движение дислокаций не является непрерывным, а носит скачкообразный характер, при этом свободное движение чередуется с остановками на локализованных препятствиях в течение некоторого времени ожидания [12, 83]. За время диффузии (время, необходимое для ДДС) примесные атомы диффундируют к дислокациям и служат дополнительным препятствием для термоактивируемого открепления дислокаций. Время ожидания уменьшается с ростом скорости деформации. Существует два основных вида взаимодействия дислокаций с неподвижными точечными дефектами: термоактивированный механизм и атермический (силовое взаимодействие возникает при превышении некоторого критического напряжения). Термоактивированный механизм взаимодействия дислокаций с точечными дефектами зависит от скорости движения дислокации и температуры, атермическое взаимодействие возникает при остановке дислокаций на барьерах Пайерлса или точечных включениях и зависит от скалярной плотности дислокаций. В ряде работ [16, 51, 61, 62, 63, 69] показан вклад данных механизмов, которые могут объяснить характер движения дислокаций и взаимодействия с точечными дефектами. В работе [92] для описания деформационного старения используется вязкопластическая модель, учитывающая чувствительность к скорости деформирования. Точный тип основных мигрирующих дефектов неизвестен, это могут быть вакансии, примесные атомы или другие дефекты. Концентрация дефектов в работе описывается следующим выражением: где - время ожидания дислокаций напрепятствиях; - время диффузии; n - коэффициент, изменяющийся от 0 до 1, значение которого связано с энергией взаимодействия между линией дислокации и движущимся дефектом. Адекватность модели проверялась на основе сопоставления с экспериментальными данными по одноосному сжатию монокристаллов оксида урана. В работе [93] приводится следующая формула для описания средней концентрации атомов растворенного вещества вокруг дислокаций: где - значение насыщения концентрации атомов растворенного вещества на дислокациях; - номинальная концентрация растворенного вещества во всем объеме кристалла; - время ожидания дислокаций на препятствиях; - время, необходимое для диффузии примесей к дислокациям; n = 2/3 - коэффициент; W - энергия взаимодействия между атомом растворенного вещества и дислокацией; - постоянная Больцмана; - температура. 4. Исследования пластической деформации твердых тел на основе дислокационной динамики Процессы неупругого деформирования и свойства поликристаллических материалов на макроуровне, как показывают многочисленные экспериментальные и теоретические исследования, существенным образом определяются состоянием эволюционирующей мезо- и микроструктуры материала. Очевидно, что для качественного воспроизведения пластической неустойчивости необходимо создание моделей, способных описывать самоорганизацию микроструктурных процессов, которые, в конечном счете, могут привести к спонтанному появлению локализации деформации. Неоднородная деформация связана с передвижением элементарных носителей пластической деформации - дислокаций. Дислокации и барьеры дислокационной (Ломера-Коттрелла, иммобильные дислокации) и недислокационной природы (облака примесных атомов, твердых частиц, выделений вторичной фазы и т.д.) рассматривают на мезо- и микроуровне. Понимание многоуровневого характера пластической деформации и ее неоднородности, которая может проявляться на различных масштабных уровнях, пришло относительно недавно [39, 66, 68]. Для корректного описания пластической деформации и ее неоднородности, для учета наиболее важных реальных физических механизмов необходимо изучать поведение материала на более низких, чем макроуровень, масштабных уровнях. В данном разделе проведен анализ некоторых моделей, описывающих эволюцию дислокационной субструктуры. В обзорной статье [73] представлены теории упрочнения, основанные на рассмотрении дислокационных субструктур, и описаны некоторые экспериментальные работы. Предполагается, что основным механизмом неупругого деформирования является движение краевых дислокаций, которые при этом взаимодействуют между собой и с дефектами кристаллической решетки: точечными (примесные атомы), линейными (барьеры Ломера-Коттрелла и др.) и объемными (выделения вторичных фаз и другие). Моделирование на основе дислокационной динамики является наиболее универсальным подходом к изучению поведения больших массивов дислокаций, их согласованного движения и взаимодействия; данные модели позволяют получить зависимости плотности дислокаций от параметров воздействия (степени деформации, скорости деформации и температуры), а также от параметров структуры (размеры зерен, субзерен, фрагментов, их ориентации). Достоинством таких моделей является способность описать структуру, прочность и пластичность металлов и сплавов на микроструктурном уровне. Дислокационная модель, развитая для анализа эволюции микроструктуры и механизмов деформации нанокристаллических материалов с размером зерна менее 100 нм, предлагается в [1].В модели рассматриваются совместные действия различных конкурирующих деформационных механизмов в процессе деформации образцов. Учитывается общая плотность дислокаций, плотность дислокаций во внутренних областях возрастает за счет действия источников Франка-Рида в стенках границ, обращенных во внутреннюю область фрагментов. Плотность дислокаций убывает за счет их стока в границы фрагментов. Модель позволила рассчитать эволюцию плотности частичных дислокаций во внутренних областях фрагментов, двойников, плотности дислокаций в границах фрагментов и полной плотности дислокаций, разориентировки между фрагментами, объемную долю двойников в зависимости отстепени деформации. Для идентификации параметров модели используются экспериментальные данные по растяжению медных образцов. На основе модели было оценено влияние отдельных деформационных механизмов на экспериментально наблюдаемые особенности механического поведения наноматериалов. В работе [64] показано, что в условиях отрицательной чувствительности к скорости деформации макроскопическая неустойчивость пластического течения может проявляться при кооперативном движении дислокаций. Переход от микроскопического описания неустойчивости пластического течения к макроскопическому может быть реализован за счет введения пространственной связи, как, например, было предложено в [74]. Отмечается, что выбор связей между микроскопическими и макроскопическими параметрами имеет важное значение в описании эффекта ПЛШ. Современные дислокационные модели позволяют описать поведение линейных дефектов при пластическом деформировании, их использование дает возможность проследить поля распределения и движения краевых дислокаций. Во многих моделях используется разделение дислокаций на различные типы (мобильные - иммобильные [34, 41, 42, 43, 44, 68, 76, 98, 99], статистически накопленные и геометрически необходимые [95]). Для описания плотностей дислокаций с одинаковыми плоскостями залегания и отличающимися векторами Бюргерса в работах [33, 39, 54, 67, 75, 86] используется разделение краевых дислокаций на положительные и отрицательные. Положительный или отрицательный знак отражает направление движения одиночной дислокации при предписанном сдвиговом напряжении. Упругопластическая конститутивная модель для описания эволюции плотности дислокаций в образцах монокристалла алюминия предлагается в [39]. Градиент деформации представляется в виде мультипликативного разложения на упругую и пластическую составляющие по аналогии с работой [40]. В рассмотрение вводятся положительные и отрицательные краевые и винтовые дислокации, учитываются такие дислокационные механизмы, как размножение дислокаций, их аннигиляция, закрепление дислокаций примесными атомами. Скорость сдвига на каждой системе скольжения (СС) определяется из уравнения Орована с учетом эволюции плотности дислокаций, приведены соотношения для определения скорости подвижных краевых и винтовых дислокаций. Плотности винтовых и краевых дислокаций рассчитываются на каждой системе скольжения. Реализация модели осуществлялась с использованием метода конечных элементов (МКЭ) в программном пакете ABAQUS, параметры модели были определены из экспериментальных данных. В работе [43] представлена конститутивная вязкопластическая дислокационная модель для описания поведения поликристаллических металлов и сплавов с ГЦК-решеткой при высоких скоростях деформации. Данная модель основана на использовании изотропного закона упрочнения и учитывает в своем рассмотрении скольжение дислокаций. Скорость сдвига на каждой СС определяется из уравнения Орована с учетом плотностей различных типов дислокаций и средней скорости их движения. В работе приводится разделение всех дислокаций на мобильные и иммобильные с введением механизмов их зарождения, размножения, аннигиляции и остановки на барьерах (границах зерен, дислокациях леса и облаках примесных атомов). Плотности подвижных и неподвижных дислокаций используются в качестве внутренних переменных, для которых сформулированы отдельные эволюционные уравнения. В работе обсуждается влияние включений из примесных атомов (их формы, размера и концентрации) на скольжение дислокаций. В предложенной модели не рассматривается неоднородное зарождение дислокаций от источников. Полученные результаты расчетов с использованием модели для образцов из алюминиевого сплава 6061-T6 показали хорошее соответствие экспериментальным данным. Модификация вязкопластической модели[43], основанная на использовании анизотропного закона упрочнения для исследования ударно-волнового нагружения, предложена в работе [44]. Скорость сдвига в каждой СС определяется плотностями соответствующих типов дислокаций и средними скоростями их движения (по уравнению Орована). Представлены эволюционные уравнения для плотностей различных типов дислокаций, учитывающие их однородное и неоднородное зарождение. Проведено сравнение результатов расчетов с применением модели с экспериментальными данными, полученными при ударно - волновом нагружении трех поликристаллических металлов (Cu, Ni, и Al). Недостатком данной модели является наличие большого числа параметров, используемых для описания эволюции субструктуры, которые трудно однозначно определить из имеющихся экспериментальных данных. Модель дислокационной динамики для описания поведения монокристаллов α-Fe и сплавов Fe-Ni, Fe-Cr, основанная на построении системы эволюционных уравнений для плотностей различных типов дислокаций, представлена в работе [76].Основным источником дислокаций полагаются источники Франка-Рида, представленные в виде конечных дислокационных сегментов; общая плотность дислокаций разделяется на мобильные и иммобильные. В рамках дислокационного подхода подвижность винтовых дислокаций относительно краевых дислокации считается незначительной. Получены графики степенной зависимости подвижности краевых дислокаций от температуры в сплаве Fe-Ni при изменении концентрации никеля от 5 до 20 %. Значения параметров, используемых в модели, определяются в процессе ее идентификации по экспериментальным данным. Анализ полученных результатов показал, что подвижность краевых дислокаций выше в сплавах Fe-Ni по сравнению с α-Fе и возрастает при увеличении концентрации Ni. Подвижность дислокаций в сплавах Fe-Cr сравнима с подвижностью краевых дислокаций в чистом железе. Основным преимуществом модели является возможность предсказать упрочнение материала за счет учета взаимодействия дислокаций и точечных дефектов. В статье [41] рассматривается модификация ранее предложенной авторами модели [42, 76] для описания поведения поликристаллов от микро- до макромасштабов. Скорость сдвига и упрочнение кристалла рассчитываются исходя из скорости движения дислокаций и скоростей их реакций. Скорости сдвига на каждой СС определяются с использованием классического закона Орована, скорости дислокаций зависят от сдвиговых и критических напряжений, а также чувствительности к скорости деформирования. В работе приведены эволюционные уравнения для плотностей дислокаций, общая плотность дислокаций разделяется на плотности положительных и отрицательных дислокаций. В целом модель незначительно отличается от изложенной в цитируемых выше работах. Результаты моделирования обнаруживают хорошее приближение к экспериментальным данным при анализе деформирования образцов из сплава титана и меди. К данному направлению относится упругопластическая модель, изложенная в работе [99]. В модели различаются следующие типы дислокаций: дислокации леса, положительные и отрицательные мобильные дислокации на СС. Поликристалл представляется в виде набора зерен, каждое зерно имеет эллипсоидальную форму и определенную ориентацию. Упрочнение кристалла описывается аддитивной суммой слагаемых: упрочнение кристалла при остановке дислокаций на дислокациях леса, при остановке дислокаций границами зерен, при торможении дислокаций на включениях из примесных атомов. Описана эволюция плотности дислокаций, а также эволюция плотности барьеров. Данная модель удовлетворительно описывает циклическое сжатие при больших пластических деформациях образцов алюминиевого сплава AA6022-T4 и может быть использована для описания деформирования широкого круга алюминиевых сплавов. Одна из первых моделей, для описания эффекта ПЛШ с учетом дислокационной динамики принадлежит Schoeck [88]. В работе рассматривается влияние ДДС и чувствительности к скорости деформирования как возможные причины возникновения эффекта ПЛШ. В модели выделяют только подвижные дислокации, плотность которых зависит от времени и скорости деформирования в явном виде и описывается с помощью дифференциального кинетического уравнения. В силу того что в процессе своего движения дислокации останавливаются на барьерах различной природы, за время ожидания примесные атомы диффундируют к дислокациям и дополнительно закрепляют их на препятствиях. Возможность переноса атомов примеси вместе с дислокациями зависит от скорости деформации и температуры. Согласно теории Коттрелла, дислокации, окруженные примесными атомами, движутся вместе с примесными облаками при низких скоростях деформирования. В случае высоких скоростей деформации дислокации движутся, «не замечая» примесных атомов. На кривых деформации данные интервалы соответствуют монотонному характеру кривой, эффект ПЛШ возникает, если временные масштабы остановок дислокаций и диффузии примесей сопоставимы. В работе учитывается среднее числолокализованных препятствий, которые пересекают дислокации в процессе своего движения до остановки на барьере, приведены выражения для оценки времени ожидания дислокаций на препятствиях и времени диффузии. Было оценено влияние отдельных деформационных механизмов на экспериментально наблюдаемые особенности поведения материала, высказаны предположения о связи наблюдаемых макроскопических неоднородностей возникновение деформационных полос на боковой поверхности образцов при их деформировании) с микроскопическими неустойчивостями. Авторы работы предполагают, что проявление макроскопической неустойчивости пластического течения зависит от характеристик деформационной машины и нагружающего устройства. В работе [56] предлагается вязкоупругая конститутивная модель, учитывающая динамическое деформационное старение, а также процессы динамического возврата и рекристаллизации. В основе модели лежат эволюционные уравнения для плотностей мобильных дислокаций и дислокаций леса, плотности дислокаций различных типов играют роль внутренних переменных. Учет микроструктуры позволил построить простую математическую модель, которая содержит небольшое количество материальных параметров. Для реализации модели используется метод конечных элементов МКЭ (пакет ABAQUS), параметры модели были определены из одноосных испытаний. Результаты затем сравнивались с экспериментальными данными и продемонстрировали хорошее соответствие между экспериментом и прогнозом модели. В работе Кубин и Эстрин [71] представлена математическая модель для описания неоднородного пластического течения. Локальная связь между напряжением и скоростью деформации задается функцией , которая может быть убывающей в определенном диапазоне скоростей деформации: , где h - представляет собой модуль упрочнения. Авторы связывают отрицательную скоростную чувствительность с динамическим деформационным старением и рассматривают изменение плотности подвижных дислокаций, плотности дислокаций леса и диффузионные процессы. Подвижные дислокации задерживаются на стопорах различной природы и дополнительно закрепляются атмосферами из примесных атомов. Описать макроскопическое поведение материала при пластической деформации невозможно без записи уравнения пластического течения, которое устанавливает связь между деформирующим напряжением и скоростью деформации. Скорость макроскопического пластического течения зачастую связывают с микроскопическими величинами, которые характеризуют дислокационную подсистему деформируемого материала, соотношением Орована , где - скорость сдвиговой деформации; b - модуль вектора Бюргерса; - плотность мобильных дислокаций; - средняя скорость их движения. Задача определения скорости пластического течения сводится к нахождению плотности подвижных дислокаций и средней скорости их движения. В настоящее время в литературе представлено большое число работ, посвященных экспериментальным и численным оценкам данных величин. Переход от микроскопического описания неустойчивости пластического течения к макроскопическому с использованием уравнения Орована или соотношений, полученных на основе уравнения Орована, используется в моделях различных авторов [55, 59, 60]. В работе [55] предлагается стохастическая модель для теоретического описания динамического деформационного старения и эффекта Портевена-Ле Шателье. Скорость пластической деформации зависит от плотности мобильных дислокаций и задается уравнением Орована. Движение подвижных дислокаций в работе рассматривается как прыжки петель дислокаций. Дислокации в процессе термической активации за время ожидания и остановки на препятствиях приобретают зигзагообразную форму. Напряжение течения должно быть увеличено для преодоления дислокациями препятствия. Процесс старения в работе рассматривается как стохастический процесс, происходящий между членами статистического ансамбля, которыми являются отдельные дислокационные петли или двойные петли в случае, когда дислокации тормозятся дислокациями леса. Изменение дислокационной плотности в объеме, в котором наряду с процессом размножения дислокаций происходит их блокировка и аннигиляция, описывается нелинейными дифференциальными уравнениями. В результате численных экспериментов на образцах из AlMg сплава были получены однородные и скачкообразные зависимости одноосного напряжения от продольной деформации в зависимости от приложенной скорости деформации. Кривые растяжения σ-ε с учетом ДДС были получены при неизменной структуре, т.е. плотность подвижных дислокаций и плотность дислокаций леса задавалась постоянной. Однако, несмотря на хорошее соответствие экспериментальных и расчетных данных для различных температур, для воспроизведения всей кривой напряжение-деформация необходимо вводить в рассмотрение изменение плотностей мобильных дислокаций и дислокаций леса. Динамическая модель для описания эффекта Портевена-Ле Шателье предложена в [59]. В статье представлена модификация модели [71], которая позволяет воспроизводить наблюдаемые типы полос деформации, отрицательную чувствительность к скорости деформации, критическую деформацию, отделяющую монотонное поведение кривой напряжение-деформация от скачкообразного, а также описать влияние ДДС на сложные динамические процессы. Отличительной особенностью модели является формулировка отдельных эволюционных уравнений для плотностей дислокаций и скоростей дислокаций. Поле напряжений получают с помощью решения краевой задачи, рассматривается одноосное растяжение при постоянной скорости деформации, скорости пластической деформации вводятся через уравнение Орована. В более поздней работе авторов [60] представлена модель, основанная на континуальной теории дислокаций. Поле упругих взаимодействий между дислокациями способно сгенерировать их коллективное поведение и самоорганизацию на локальных участках, например появление полос деформации. Средняя скорость дислокаций в плоскости скольжения описывается следующим соотношением: где - «начальная» скорость дислокаций при , определяемая из экспериментальных данных; - энтальпия, активационный объем, постоянная Больцмана и температура; - сдвиговое напряжение; - критическое напряжение для преодоления препятствий. Для определения изменения плотности и скорости дислокаций записаны феноменологические соотношения. Данная модель позволяет описать распространяющийся фронт пластического течения в однородном материале. Авторы ряда работ для описания дислокационной динамики используют вязкопластические модели [35, 36, 37]. Общая плотность дислокаций представляется как сумма плотностей мобильных и иммобильных дислокаций, для описания взаимодействия дислокаций с облаками примесных атомов, образующих атмосферы Коттрелла, вводятся эволюционные уравнения для плотности барьеров. Динамическая модель, представляющая собой модификацию модели [37] и ориентированная на описание эволюции полос деформации, рассмотрена в [48]. Механизм эффекта прерывистой текучести связывают с динамическим деформационным старением, обусловленным взаимодействием между подвижными дислокациями, временно остановившимися на препятствиях и диффундирующими атомами примеси, которые стремятся дополнительно закрепить дислокации. Модель основана на совокупности дифференциальных эволюционных уравнений для плотностей мобильных, лесовых и «захваченных» (облаками примесных атомов) дислокаций. Модель позволяет описать некоторые особенности эффекта ПЛШ, а именно возникновение эффекта в определенном диапазоне скоростей деформации и температур, возникновение отрицательной скоростной чувствительности. В работе [35] авторы предлагают модель, основанную на совокупности эволюционных уравнений для плотностей мобильных, лесовых и «захваченных» (облаками примесных атомов) дислокаций. Эволюционные соотношения в виде дифференциальных уравнений в частных производных включают члены, отвечающие за генерацию мобильных дислокаций, их захват диффундирующими примесными атомами, аннигиляцию, переход иммобильных дислокаций в мобильные за счет термической и механической активации. Для анализа модели использованы методы нелинейной динамики; показано, что модель качественно описывает наличие отрицательной скоростной чувствительности (снижение напряжения течения при росте скорости деформации) в определенном диапазоне скоростей деформаций, образование известных типов полос скольжения (А, В и С). Модель, содержащая эволюционные уравнения для плотности иммобильных дислокаций и концентрации примесных атомов, представлена в [77]. Напряжения течения полагаются равными сумме сопротивления от дальнодействующих барьеров дислокационной природы и составляющей, зависящей от температуры, скорости деформации и концентрации примесных атомов. Основным механизмом диффузии примесных атомов к дислокациям считается вакансионный, от концентрации вакансий зависит коэффициент диффузии примеси; получено соотношение для определения кинетики изменения концентрации вакансий. Для идентификации модели использованы результаты на сжатие до деформации 0,6 образцов из нержавеющей стали 316L при скоростях деформации (5×10-4÷10) с-1 и температурах от комнатной до 1300 оС. Другая серия экспериментов, реализованных по более сложным программам нагружения, включающим выдержки при различных температурах и резкое изменение скорости деформации, была проведена для верификации модели. Результаты сопоставления теоретических и экспериментальных данных свидетельствуют об адекватности модели, за исключением описания влияния деформационного старения на поведение материала в интервале температур от 600 до 800 оС. Модель не позволяет адекватно описать начало динамического деформационного старения для данных температур. В [95] рассматривается модель монокристалла, основанная на континуальной теории дислокаций мезоуровня и ориентированная на исследование прерывистой пластичности, реализующейся за счет образования полос скольжения. В качестве определяющего соотношения используется закон Гука в скоростной форме, градиент вектора скорости перемещений представляется суммой упругой и пластической составляющих. Пластическая составляющая градиента скорости перемещений учитывает как статистически накопленные дислокации (дислокации, хаотически распределенные в кристалле и не приводящие к дальнодействующим напряжениям), так и геометрически необходимые дислокации, называемые в работе «полярными дислокациями» (дислокации, образующие дислокационные субструктуры типа стенок и др., приводящие к разориентировкам частей кристалла и дальнодействующим полям напряжений). Трансляционная составляющая пластической составляющей градиента скорости перемещений устанавливается уравнением Орована, в котором скорость перемещений (ансамбля) дислокаций по системе скольжения определяется экспоненциальной функцией эффективного напряжения. Эффективное напряжение на системе скольжения определяется разностью действующего касательного напряжения и трех составляющих сопротивления деформации: сопротивления решетки, барьеров дислокационной природы (дислокаций леса) и вклада от деформационного старения. Приведены эволюционные уравнения для плотностей мобильных дислокаций и дислокаций леса, а также феноменологические соотношения для характерного времени старения и концентрации атомов примеси. Для реализации модели используется метод конечных элементов. Рассматривается одноосное растяжение монокристаллического образца размерами 2×0,5×12 мм, аппроксимированного совокупностью 12×6×200 шестигранных элементов; положение оси растяжения выбрано таким образом, что реализуется одиночный сдвиг. Получены и анализируются результаты расчетов для скоростей деформации в диапазоне (3×10-5 ÷ 6×10-5) с-1. Показано, что при всех исследованных скоростях деформации описывается «зуб текучести», при скорости деформации 3×10-5 с-1диаграмма σ-ε приобретает монотонный характер. Для остальных анализируемых скоростей деформации имеют место неоднородные по времени и пространству неупругие деформации (движущаяся полоса сдвига) и зубчатые диаграммы σ-ε. Детально анализируется влияние распределения примесных атомов и дислокаций на поля напряжений в окрестности полосы скольжения. Результаты детальных экспериментальных исследований локализации деформаций в специально препарированных тонких бикристаллических образцах представлены в [87]. В процессе деформирования in-situ с помощью цифровых (оптической и инфракрасной) камер определялись перемещения случайно распределенных точек на поверхности образца (которые затем использовались для расчета деформаций) и локальные температуры. Эксперименты были ориентированы на идентификацию и верификацию разрабатываемой авторами модели. Последняя основана на МКЭ (пакет ABAQUS) и вязкопластической физической теории для расчета скоростей сдвига по системам скольжения; критические напряжения устанавливались с помощью дислокационной подмодели, описывающей анизотропное (деформационное и латентное) упрочнение и разупрочнение за счет динамического возврата. По скоростям сдвига и касательным напряжениям на СС рассчитывалась локальная мощность источников тепловыделения (диссипации). Сравнение экспериментальных и теоретических результатов свидетельствует об удовлетворительном соответствии по полям деформации, несколько хуже согласуются данные по тепловым источникам. Исследованию прерывистой пластичности ГЦК-материалов (медь, медные сплавы, аустенитные нержавеющие стали) при криогенных (близких к 0 К) температурах посвящена статья [89]. Для указанных материалов установлено, что прерывистое течение имеет место при температурах ниже характерных значений (порядка 20-40 К), при которых происходит переход от движения винтовых к активизации краевых дислокаций. В качестве возможных механизмов прерывистой пластичности отмечаются фазовые переходы, динамическое деформационное старение, прорыв барьеров Ломера-Коттрелла скоплениями дислокаций, основное внимание в работе уделено двум последним. Для описания прерывистой пластичности используется дислокационная модель, включающая эволюционные уравнения для плотностей дислокаций и барьеров Ломера-Коттрелла. На основе связи между микропараметрами и макропеременными осуществляется переход от модели микроуровня непосредственно к (одномерной) модели макроуровня. Отдельный «зубец» прерывистой пластичности представляется четырьмя стадиями: 1) упругое деформирование, 2) пластическое монотонное деформирование, 3) резкое падение напряжений, 4) релаксация напряжений за счет движения групп дислокаций. Для идентификации модели использованы данные по одноосному растяжению образцов из нержавеющей стали 316 при температуре жидкого гелия (4,2 К). Результаты численных расчетов обнаруживают удовлетворительное соответствие экспериментальным данным. Вплотную к рассмотренной примыкает работа [90], в которой особое внимание уделяется анализу многоосного деформирования. Для описания деформационного старения в [38] использовалась микроскопическая модель, основанная на совокупности дифференциальных эволюционных уравнений для плотностей мобильных и иммобильных дислокаций, скорость сдвига определялась уравнением Орована. Результаты численного анализа сопоставляются с экспериментальными данными для образцов из пересыщенного твердого раствора закаленных сплавов системы Al-Mg-Si, подвергаемых одноосному растяжению в диапазоне температур от 25 до 235 °С и скоростях деформации от 10-4 до 5∙10-2 с-1. Численная реализация была выполнена в ABAQUS с помощью метода конечных элементов. Авторы проводили исследование влияния концентрации примеси, температуры, деформации и скорости деформации на характер деформационной кривой. Было установлено, что в зависимости от скорости деформации и температуры за возникновение ДДС могут быть ответственны два механизма. При низких температурах доминирующим является механизм остановки дислокаций на атомах растворенного вещества, а при высоких температурах - на выделениях вторичных фаз (осадков), которые образуются в процессе деформирования. Влияние ДДС становится более заметным при низких скоростях деформации для всех температур. Вклад в скорость пластического течения от остановки дислокаций атмосферами примесных атомов пропорционален концентрации атомов растворенного вещества вдоль линии дислокации. Вклад вторичных фаз в напряжение течения зависит от объемной доли осадков, среднего радиус атомов, критического радиуса включений, когда механизм перерезания твердых частиц изменяется на механизм переползания. Показано, что оба механизма, ответственные за проявление ДДС, приводят к отрицательной чувствительности скорости деформации. Начальное напряжение пластического течения возрастает с увеличением температуры деформации из-за наличия осадков в образцах, предварительно отожженных при температуре 520 °С в течение 2 часов. Предлагаемая модель позволила получить адекватные деформационные кривые, а также воспроизвести полосы локализованной деформации. В статье [66] представлена конститутивная пластическая модель для описания поведения поликристаллических жаропрочных сплавов на основе никеля. В рассмотрение вводятся три масштабных уровня: макроуровень, соответствующий представительному объему поликристалла, мезоуровень (или уровень отдельного зерна) и уровень субзерен. С помощью модели исследуется монотонное нагружение. Экспериментальные исследования показали, что на уровне субзеренважен учет таких механизмов, как микродвойникование, размножение, аннигиляция, скольжение дислокаций. Результаты моделирования на уровне субзерна, где вводятся плотности дислокаций, используются для определения конститутивных параметров для модели мезоуровня. На мезоуровне структура поликристаллического жаропрочного сплава на основе никеля моделируется с использованием дополнительной пластической подмодели с помощью метода конечных элементов. Результаты моделирования на каждом уровне сравниваются с экспериментальными данными. Использование трехуровневой модели позволило описать субзеренную структуру поликристаллического материала. Существенным преимуществом этой модели является высокая эффективность, которая позволяет включать в рассмотрение микроструктуру поликристаллов, при одновременном снижении вычислительных затратбез значительной потери точности. В статье [75] предложена новая континуальная дислокационная модель, которая включает в себя всестороннее рассмотрение упругих взаимодействий между дислокациями, а также описание непрерывного движения, размножения дислокаций и их аннигиляцию. Представлены формулировки эволюционных уравнений для плотностей дислокаций, которые описывают упругие взаимодействия между дислокациями, что стало возможным в результате обобщения описания упругих взаимодействий между участками дислокаций. Эффективность численной реализации достигается благодаря использованию упрощенных математических уравнений, которые содержат только элементарные функции. В модели используется динамическая постановка,в которой плотности мобильных и иммобильных дислокаций не вводятся отдельно, проводится расчет всех действующих сил на дислокации каждого типа, а затем рассчитывается эволюция плотности дислокаций с помощью законов дислокационной динамики. Численная реализация была выполнена с помощью метода конечных объемов при использовании периодических граничных условий для модели идеального монокристалла алюминия. Построены кривые σ-ε для образцов различных размеров, показано, что предел прочности при растяжении возрастает с уменьшением размера образца, обнаружено, что образцы малых размеров демонстрируют зуб текучести. Результаты, полученные в модели, показали хорошее соответствие экспериментальным данным. Заключение Предлагается краткий обзор теоретических работ, основанных на физических теориях пластичности, посвященных описанию особенностей деформирования сплавов в температурно-скоростных диапазонах, в которых существенное влияние на поведение материалов оказывают диффузионные процессы. Особое внимание уделено описанию эффекта Портевена-Ле Шателье. В силу того что процессы неупругого деформирования и свойства поликристаллических материалов на макроуровне существенным образом определяются состоянием эволюционирующей мезо- и микроструктуры материала, в статье затронуты также некоторые работы, содержащие описание эволюции дислокационных субструктур, результаты исследования влияния на характеристики скачкообразной деформации микроструктуры сплавов и концентрации примесных атомов. Основной задачей при исследовании эффекта ПЛШ, играющего отрицательную роль при обработке промышленных сплавов (ухудшает механические свойства изделий, увеличивает шероховатость поверхностей изделий, понижает усталостную прочность и коррозионную стойкость), является определение области существования неустойчивости с целью ее исключения из технологических режимов обработки изделий. Математическое моделирование «тонкостей» процессов эволюции дислокационных субструктур и точечных дефектов, особенно во внутренних областях образцов, выступает зачастую единственной возможностью количественного анализа микроструктурных процессов, так как натурный эксперимент сопряжен с огромными принципиальными сложностями и высокими затратами. Анализ макрофеноменологических моделей, представляющих предмет предыдущей публикации [24], и моделей, основанных на физических теориях пластичности, позволяет рекомендовать многоуровневый подход на основе физических теорий пластичности в качестве наиболее предпочтительного для описания прерывистой деформации. Предсказательная способность моделей во многом зависит от задания исходных данных. На основе экспериментальных и теоретических работ, посвященных описанию процессов неупругого деформирования и свойствам поликристаллических материалов, были выделены наиболее значимые механизмы для описания неустойчивого пластического течения: взаимодействие дислокаций с барьерами дислокационной (Ломера-Коттрелла, иммобильные дислокации) и недислокационной природы (облака примесных атомов, твердых частиц, выделений вторичной фазы и т.д.). В разрабатываемой авторами многоуровневой модели поликристаллических металлов планируется учет вышеупомянутых механизмов, верификация параметров модели будет осуществляться по имеющимся в литературе экспериментальным данным.

About the authors

P V Trusov

Perm National Research Polytechnic University

E A Chechulina

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 322

PDF (Russian) - 111