Атомистические механизмы и кинетика пластической деформации металлов при высокоскоростной деформации

Янилкин, Алексей Витальевич

Атомистические механизмы и кинетика пластической деформации металлов при высокоскоростной деформации тема автореферата и диссертации по физике, 01.04.07 ВАК РФ

Янилкин, Алексей Витальевич АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Долгопрудный МЕСТО ЗАЩИТЫ

2010 ГОД ЗАЩИТЫ

01.04.07 КОД ВАК РФ

Диссертация по физике на тему «Атомистические механизмы и кинетика пластической деформации металлов при высокоскоростной деформации»

Автореферат диссертации на тему "Атомистические механизмы и кинетика пластической деформации металлов при высокоскоростной деформации"

На правах рукописи УДК 538.9+539.3

Янилкин Алексей Витальевич

Атомистические механизмы и кинетика пластической деформации металлов при высокоскоростной деформации

Специальность 01.04.07 - Физика конденсированного состояния

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

2 8 0КТ 2010

Долгопрудный 2010

004611723

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования «Московский физико-технический институт (национальный исследовательский университет)» на кафедре физики высокотемпературных процессов (Базовый институт Объединённый институт высоких температур РАН)

Научный руководитель:

Официальные оппоненты:

доктор физико-математических наук, зав. отд., профессор Норман Генри Эдгарович (Объединённый институт высоких температур РАН)

доктор физико-математических наук, заместитель директора по науке, профессор Рыжов Валентин Николаевич (Институт физики высоких давлений РАН им. Л.Ф. Верещагина, Троицк)

Ведущая организация:

доктор технических наук, зав. отд. профессор Лурье Сергей Альбертович (Институт прикладной механики РАН, Москва)

Институт проблем машиноведения РАН, Санкт Петербург

Защита состоится 27 октября 2010 г. в 16.00 часов на заседании диссертационного совета Д 212.156.06 при Московском физико-техническом институте (национальный исследовательский университет) по адресу: 117393, г. Москва, ул. Профсоюзная, д. 84/32, корпус В-2.

Отзывы направлять по адресу: 141700, г. Долгопрудный Московской обл., Институтский переулок, д. 9., МФТИ, Диссертационный совет Д 212.156.06.

С диссертацией можно ознакомиться в библиотеке МФТИ. Автореферат разослан 27 сентября 2010 г.

Ученый секретарь диссертационного совета кандидат технических наук

Н.П. Чубинский

1 Общая характеристика работы

Диссертация посвящена разработке теоретического описания кинетики и механизмов пластической деформации металлов при высокоскоростной деформации. Построена модель явления и исследованы процесс гомогенного зарождения дислокации, явление фононного торможения, влияние точечных дефектов и препятствий на подвижность дислокаций. Использован метод молекулярной динамики (МД). Полученные результаты обобщены в виде аналитических соотношений.

Актуальность работы. Высокоскоростное нагружение встречается во многих технологически важных приложениях, включая высокоскоростную обработку и формовку материалов, создание ударостойких материалов, пробивание защитных оболочек, ударное повреждение авиационных и космических летательных аппаратов и конструкций. Механический отклик многих материалов (металлов, керамик, стёкол, полимеров) при высокоскоростном деформировании отличается от случая статических нагрузок. Так для многих металлов зависимость напряжения течения от скорости деформирования резко усиливается с превышением скорости деформирования и 103 — 104с-1 [1,2]. Явление можно интерпретировать как следствие изменения механизма движения дислокаций.

При малых скоростях движения дислокации преодолевают препятствия в результате совместного действия приложенного напряжения и тепловых флуктуации. Вследствие этого увеличение температуры сопровождается понижением предела текучести материалов. Для деформирования с высокой скоростью необходимо приложить более высокие напряжения. При скорости деформирования больше некоторой пороговой действующие напряжения оказываются достаточными для обеспечения динамического преодоления препятствий (барьеры Пайерлса), без дополнительного вклада тепловых флуктуации. При этом доминирующим механизмом торможения дислокаций становится перекачка энергии дислокации в колебания кристаллической решетки или, в зависимости от температуры, в электронную подсистему. В отличие от области термофлуктуа-ционной подвижности, скорость дислокаций в динамической области падает с температурой в соответствии с увеличением плотности элементарных возбуждений. Поэтому при высоких скоростях деформирования для некоторых материалов наблюдается (аномальное) возрастание напряжения течения с увеличением температуры [1].

Цель работы. 1) Разработка метода расчёта скорости зарождения дислокаций в бездефектном кристалле. Исследование зависимости скорости зарождения от сдвиговых напряжений и температуры.

2) Развитие моделей для исследования подвижности одиночной дислокаций в монокристалле. Исследование скорости движения в зависимости от сдвиговых напряжений и температуры для ОЦК и ГЦК металлов.

3) Определение механизмов влияния дефектов (вакансий, полостей, включений) на подвижности дислокаций. Исследование их влияния на коэффициент торможения. Определение критических напряжений, необходимых для преодоления дефектов.

4) Развитие многомасштабного подхода для описания скорости пластической деформации на основе результатов МД.

Научная новизна работы. Исследованы атомистические механизмы и кинетика элементарных процессов пластической деформации в условиях высокоскоростной деформации с помощью МД.

1) Исследован механизм гомогенного зарождения дислокаций за счет локального сдвига одной группы атомов относительно другой на вектор Бюргер-са частичной дислокации. Процесс носит стохастический характер. Показано, что зависимость скорости зарождения от сдвиговых напряжений может быть аппроксимирована в рамках формулы Журкова при введении температурной зависимости параметров.

2) Исследована зависимость скорости движения одиночной дислокации в ОЦК и ГЦК кристаллах от сдвиговых напряжений. Для коэффициента фонон-ного трения В, характеризующего наклон линейного участка, получена зависимость от температуры. Показано подобие зависимостей коэффициента трения и периода решетки от температуры.

3) Изучены механизмы влияния дефектов решетки на подвижность дислокации в кристаллах Мо, Си и А1. В ОЦК равномерно распределённые вакансии концентрации с не оказывают существенного влияния на величину критических напряжений <тсг, необходимых для движения дислокации. Наличие вакансий приводит к возрастанию трения В ~ с. В случае ГЦК коэффициент трения практически не возрастает, а увеличивается барьер для движения осг ~ с.

4) Развит подход, позволяющий описывать макроскопическую пластическую деформацию с использованием результатов МД.

Практическая ценность работы. Результаты диссертации могут быть использованы для оценки пластических характеристик материалов при импульсных нагрузках, расчетов кинетики пластической деформации при рассмотрении задач механики (например, в рамках моделей сплошной среды).

Положения, выносимые на защиту.

1. Подход для определения скорости J зарождения дислокаций и зависимость J от сдвиговых напряжений и температуры.

2. Зависимость скорости движения дислокаций от сдвиговых напряжений в динамическом режиме для ОЦК и ГЦК кристаллов на примере Mo, AI и Си. Температурные зависимости коэффициента фононного трения.

3. Характер влияния вакансий на подвижность дислокаций.

4. Механизмы преодоления полостей и включений при движении дислокации в Mo, AI. Зависимость сдвиговых напряжений, необходимых для преодоления препятствий, от расстояния между ними и их размера.

5. Подход, позволяющий рассчитывать кинетику пластической деформации монокристаллов с использованием данных МД расчетов.

Результаты диссертации докладывались на конференциях: МФТИ 2002-08; Уравнения состояния вещества (п. Эльбрус 2005-10); Advanced Problems of Mechanics (Санкт-Петербург 2007-09); Многоуровневые подходы в физической мезо-механике (Томск 2008); Аэрофизика и физическая механика классических и квантовых систем (Москва 2007); Проблемы физики ультракоротких процессов в сильнонеравновесных средах (Новый Афон 2005-09); Физика прочности и пластичности материалов (Самара 2009); 10-ые Забабахинские научные чтения (Снежинск, 2010); Параллельные вычислительные технологии (Уфа 2010); Nucleation theory and applications (Дубна 2009-10); Computational physics (Brussels 2007, Ouro Preto 2008); Shock Compression of Condensed Matter (Hilo 2007); Multiscale Materials Modeling (Tallahassee 2008); Joint U.S.-Russia Conference on Advances in Materials Science (Prague 2009); IUTAM Symposium on Dynamic Fracture and Fragmentation (Austin 2009); New Models and Hydrocodes for Shock Wave Processes (Paris 2010).

Публикации. По материалам диссертации опубликовано 9 работах в реферируемых журналах, 15 работ в сборниках и тезисы конференций.

Структура и объем диссертации. Диссертация состоит из введения, семи глав и заключения, на 102 страницах, включает 35 рисунков, библиографию из 108 наименований.

2 Содержание работы

В первой главе дается краткое обоснование актуальности, научной новизны и практической значимости задач, рассматриваемых в работе. Сформулированы цели работы и основные положения, выносимые на защиту.

Вторая глава носит обзорный характер. Кратко описываются доступные на сегодняшний день данные о пластическом деформировании материалов при высокоскоростном нагружении. Изложены основные теоретические подходы для описания пластической деформации. Подробнее рассмотрены методы атомистического моделирования, используемые для изучения дефектов, их подвижности и взаимодействия друг с другом.

В разделе 2.1 дается обзор экспериментальных результатов по высокоскоростному деформированию. В качестве основных экспериментов обсуждаются ударно-волновые, в которых достигаются высокие скорости деформирования. Представлен краткий анализ результатов и подчеркнуто отличие от нагружения при обычных условиях.

В разделе 2.2 дается краткое введение в метод МД, который используется для исследования атомистических механизмов. Рассмотрены различные граничные условия. Большое внимание уделено обсуждению потенциалов межатомного взаимодействия и данных, на основе которых была построена их параметризация. Например, в работе [3] отмечено влияние формы потенциала на фазовую диаграмму. Отдельно рассмотрены способы анализа термодинамических и структурных характеристик.

В разделе 2.3 обсуждаются многомасштабные подходы, позволяющие использовать данные об элементарных процессах пластической деформации для макроскопического моделирования. Это позволяет включить информацию о кинетике атомистических механизмов.

Третья глава посвящена исследованию гомогенного зарождения дислокаций под действием сдвиговых напряжений на примере алюминия. Исследованы механизм и стохастические свойства процесса гомогенного зарождения дислокаций. Предложен подход для определения скорости зарождения на основе рассмотрения времён жизни метастабильного состояния. Для получения первичных данных использован метод МД. Определены зависимости скорости нуклеа-ции от значений сдвиговых напряжений для нескольких температур. На основе полученных данных проведено сравнение с теоретическими подходами, сдела-

ны выводы об ограниченности предсказательной силы теоретических подходов. Проведена оценка области сдвиговых напряжений и температур, в которой может реализовываться механизм гомогенного зарождения дислокаций.

В разделе 3.1 обсуждаются различные механизмы зарождения дислокаций (например, при росте кристалла [4]) и роль гетерогенностей. Представлены теоретические подходы, использующиеся для расчета скорости нуклеации 3.

Наиболее полное теоретическое описание основано на выражении для энергии зародыша в рамках приближения линейной теории упругости [5]:

ДС(г) = 2,г.2-" "'2

4 г

1п(—)~ 2 Р

- 1тг2Ьа + 7ГГ27 (2)

2(1 - г/) 4тт

где г - радиус петли, А(3(г) - свободная энергия петли, из - дебаевская частота, Ь - вектор Бюргерса, 7 - энергия дефекта упаковки, V - коэффициент Пуассона, /1 - модуль сдвига, р - характерный размер ядра дислокации.

Другим подходом является формула Журкова [б].

3 = .¡0ехр(-и к°У) (3)

где То - максимальная скорость нуклеации, II - энергия зарождения, V - актива-ционный объем, характеризующий релаксацию упругой энергии при образовании зародыша, и и V не зависят от а и Т. Эта формула зачастую используется для аппроксимации экспериментальных результатов.

В разделе 3.2 проводится сравнение энергии образования петли частичной дислокации, рассчитанной с помощью атомистического моделирования и теории упругости (рис. 1а).

Для теоретического описания были проведены дополнительные расчеты для определения таких величин, как 7 - энергия дефекта упаковки, и - коэффициент Пуассона, /х - модуль сдвига. Что касается р, то его значение выбиралось для каждой величины <тху таким, чтобы обеспечить наилучшее согласие расчета по формуле (2) с результатами МД. Полученная зависимость р{аху) представлена на рис. 16. Из теории [5] следует, что р(0) примерно равно Ъ/4 для винтовой и Ь/8 для краевой. Результаты наших расчетов близки к значению для краевой дислокации, но зародыш петли состоит как из краевой, так и из винтовой дислокаций. Отличие размера ядра от теории, возможно, связано с изменением структуры при больших степенях деформации.

0.00

0.5 1.0 1.5 2.0

а , ГПа

0.0

Рис. 1: а) - Энергия петли дислокации в зависимости от радиуса при Т=0 К. Точки - наш расчет с помощью минимизации энергии при аху\ 1 - 0.0 ГПа, 2 - 0.8 ГПа, 3 - 1.1 ГПа, 4 - 1.4 ГПа, 5 - 1.6 ГПа, сплошные линии - расчет по формуле (2) с использованием зависимости р{<?ху), представленной на рис.б. Для сравнения показаны теоретические значения р(0) [5] для винтовой дислокации (треугольник)

Раздел 3.3 посвящен исследованию механизма зарождения. Показано, что процесс зарождения заключается в локальном сдвиге одной части атомов относительно другой на вектор Бюргерса частичной дислокации. Смещение происходит в плоскости скольжения {111}, наблюдаемый размер зародыша составляет 20-30 атомов (рис.2). Процесс носит термофлуктуационный характер.

Поскольку процесс распада метастабильного состояния носит флуктуацион-ный характер, то моменты времени зарождения дислокаций т случайны. Для определения среднего времени необходимо провести несколько десятков независимых расчётов, подобных рис. За. Все они должны соответствовать одной степени метастабильности макроскопического состояния. Под макросостоянием в данном случае понимается заданная температура Т и степень сдвиговой деформации, которую можно характеризовать значением сдвиговых напряжений аху. На рис. 36 показано распределение времён жизни, процесс зарождения носит Пуассоновский характер. Ансамбль независимых траекторий, соответствующих одному макросостоянию, набирается путём использования различных реализаций распределения Максвелла по скоростям, соответствующих одной и той же температуре.

В разделе 3.4 представлены зависимости от скорости зарождения от сдви-

и краевой (круг).

Рис. 2: Момент зарождения петли дислокации. Т = 300К, аху = 1.76 ГПа. Показаны только атомы с дефектным локальным окружением, (а) - вся ячейка, (б) - вид сверху на область зародыша.

а , ГПа

*У

\ X, ПС

100 150 200 250

Рис. 3: Пример определения скорости зарождения дислокаций из времени жизни метастабильного состояния, (а) - релаксация напряжений после фиксирования степени деформирования, (б) - распределение по временам жизни. Т — 300К. -время динамической памяти.

говых напряжений и температуры (рис. 4). Показано, что полученные зависимости расходятся со значениями, рассчитанными на основе формулы (1). Принимая во внимание хорошее согласие в энергии зарождения, можно предположить, что неправильная форма зародыша вносит существенный вклад в

J,C

/ a

ст , ГПа

1.76 1.78 1.80 1.82 1.50 1.52 1.54 1.56 1.26 1.28 1.30 1.32

Рис. 4: Зависимость скорости иуклеации дислокаций от значений сдвиговых напряжений для трёх температур Т: а) 300К, б) 500К, в) 700К. МД расчёты для двух размеров системы: красные кружки - N = 180000, синие ромбы - 504000 частиц. Сплошные линии - скорости иуклеации, рассчитанные по формуле (2). Пунктир -аппроксимация по формуле (3) всех МД точек с одними и теми же значениями U и V. Штрих-пунктир - аппроксимация по формуле (3) с подбором U и V для каждой температурной зависимости.

энергию зарождения. Полученные результаты могут быть аппроксимированы с помощью формулы Журкова (3) в предположении зависимостей U и У от Т.

В четвертой главе исследуется подвижность одиночной дислокации в кристалле. Рассмотрено поведение дислокации как в недеформированном кристалле, так и под действием сдвиговых напряжений. Расчеты проведены для краевой дислокации в AI, Си и Мо.

В разделе 4.1 представлен подход для создания и моделирования движения одиночной дислокации. Раздел 4.2 посвящен исследованию колебаний дислокационной линии в недеформированном кристалле под действием термических флуктуаций, так называемое "дыхание" дислокаций.

В разделе 4.3 исследуется поведение дислокаций под действием сдвиговых напряжений (рис.5). Показано, что существует три участка на зависимости скорости дислокации от сдвиговых напряжений.

Во-первых, движение дислокации при напряжениях меньше барьера Пай-ерлса. Поскольку преодоление барьера термофлуктуационное, то за времена расчета дислокация практически не смещается.

Второй участок характеризуется линейной зависимостью скорости движения от сдвиговых напряжений. На основе этой зависимости вводится коэффициент

30 п

20-

10-

v,, А/пс

ГПа

о.о

0.2

0.4

В(Т)/В(0)

0.6 о

А 1 -2 • 3

— 4

— 5

. 'л.

т/е

Рис. 5: а) - Зависимость скорост!г движения дислокации гу в AI от сдвигового напряжения аху. Сплошные линии - аппроксимации по формуле (4). б) - Температурная зависимость коэффициента динамического торможения дислокаций В(Т). Значения температуры и коэффициента торможения для алюминия из МД расчетов обезразмерены на температуру 0 = 230К и величину В(0). 1 - молекулярная динамика; 2 - аппроксимация данных МД по формуле (5); 3 - эксперимент [7]; 4 -эксперимент [8]; 5 - теоретическая зависимость [2].

трения В = аЬ/У4, где - скорость дислокации, а - сдвиговые напряжения. Торможение на этом участке определяется взаимодействием движущейся дислокации с фононной подсистемой.

На третьем участке зависимости наблюдается насыщение с приближением к скорости звука. Проведено сравнение момента выхода на насыщение и скорости распространения поперечных волн. Насыщение наступает при меньших значениях скоростей и в расчетах не удалось значительно приблизиться к скорости звука.

Предложена общая аппроксимационная формула для трёх участков, включая зависимость от температуры:

= 0 если а < <тсг

(а - сгсг) • 6 = В ■ 1^/(1 - 0,*/г>тах)4) если а > асг

где сгсг - барьер Пайерлса (25 МПа), утах - асимптотическая скорость дислокации (2700 м/с). Скорость звука 2800-3100 для Т = 900-100К.

Определена зависимость коэффициента фононного трения от температуры (рис. 56). Для сравнения приведены экспериментальные данные, взятые из работ [7,8] и теоретическая зависимость из [2]. С ростом температуры сказываются отклонения от линейного закона. Отклонение от линейности хорошо описываются функцией четвертой степени:

Абсолютное значение коэффициента торможения дислокации при Т = 300А', определенное из обработки МД расчетов, составило приблизительно 0.9 х 10~5 Па-с, что близко к экспериментальным значениям. По оценкам различных экспериментов оно варьируется от 0.5 х Ю-5 Па-с (затухание ультразвука, 1970 [8]) и 2.6 х Ю-5 Па-с (подвижность дислокаций, 1969 [7]) до 17 х 10~5 Па-с (затухание ультразвука, 1966, см. обзор [2]).

Проведено сравнение коэффициентов трения для различных металлов (рис. 6а). На рис. 66 проводится сопоставление зависимости коэффициента трения от температуры и зависимости периода решетки для двух потенциалов, описывающих Mo (FS - [9], Dud - [10]). Эти потенциалы существенно отличаются от друга поведением коэффициента теплового расширения. В случае потенциала Dud в области Т = 300-900К равновесное значение периода решетки уменьшается, что соответствует отрицательному значению коэффициента теплового расширения. Аналогично ведет себя зависимость коэффициента трения, причем, зависимости совпадают в приведенных координатах. Такое поведение, возможно, определяется общей природой теплового расширения и фононного торможения - энгармонизмом колебания атомов в потенциальной яме.

Раздел 4.4 посвящен исследованию влияния конечных размеров расчетной ячейки на полученные результаты по подвижности дислокаций. Определены геометрические параметры системы, необходимые для моделирования.

В пятой главе рассмотрено влияние точечных дефектов и атомов различных примесей, распределенных в твердом растворе. Такие дефекты создают поля препятствий движению дислокаций, в результате чего внешняя нагрузка, требуемая для такого движения значительно возрастает по сравнению с величиной для чистых металлов. Применительно к высокоскоростной деформации, когда скорость движения дислокаций велика и движение происходит надба-рьерно, возникают эффекты, сходные по механизмам с фононным трением,

15n

B,1Ö5 Па-с

0.5

Т. К

400

800

1200

1600

400

800

1200

1600

Рис. б: а) - температурная зависимость коэффициента фоионного трения AI, Си и Мо. б) - сопоставление приведенных зависимостей коэффициента трения и периода

решетки.

рассмотренным в четвертой главе.

В разделе 5.1 обсуждается статическое взаимодействие краевой дислокации в Мо с вакансией. Показано, что разница энергии образования вакансии вдали от дислокации и на расстоянии г хорошо описывается теорией упругости (Д£у(г) = Р{г) ■ ÖQv, где Р(г) - давление, - активационный объем) на расстояниях 10 А от ядра, а дальше значения с атомистическим моделированием значительно расходятся.

Раздел 5.2 посвящен влиянию вакансий на подвижность дислокаций. Рассмотрены кристаллы ОЦК на примере Мо и ГЦК - Си. Показано, что в Мо введение равномерно распределенных вакансий приводит к увеличению динамического коэффициента торможения В, тогда как минимальные напряжения <тсг, необходимые для движения дислокации в динамическом режиме, практически не меняются (рис. 7а). исг определялась путем линейной аппроксимации зависимости Vd от <т до нулевого значения скоростей. Рост динамического торможения связан со значительным увеличением осцилляции, возникающих в результате отрыва дислокационной линии от вакансии. Такой вывод подтверждается прямо пропорциональной зависимостью коэффициента торможения от концентрации вакансий (рис. 76), то есть пропорционально числу вакансий, с которой взаимодействует дислокация.

В Си наблюдается противоположное влияние вакансий (рис. 8). Наблюда-

v, А/пс

а , ГПа

ху

15.0

12.5

10.0-

5.0

В, 10"5 Па*с б

vac, %

0.0 0.2 0.4 0.6 0.8 1.0

0.0 0.5 1.0 1.5 2.0

2.5

Рис. 7: а) - зависимость скорости дислокации в Мо от сдвиговых напряжений для различных концентраций вакансий: 1 - 0 %, 2 - 0.5 %, 3 - 1 %, 4 - 2 %. б) — Зависимость коэффициента торможения от концентрации вакансий. Пунктирная линия - линейная аппроксимация.

ется значительное увеличение напряжений ггсг, тогда как коэффициент трения практические остается постоянным. Аналогичное поведение наблюдалось при моделировании движения дислокации в сплаве Al-Mg [11]. В работе [12] в сплаве Fe-Cu наблюдается похожее поведение, хотя Fe является ОЦК металлом.

15-1

Vd, А/пс

и»/ к> Ч*-1 ь*)''

■ 1

О 2

♦ 3

0.06

0.04-

0.02

аху, ГПа о.оо

о , МПа

o-J

0.00 0.05 0.10 0.15 0.20 0.0

vac, %

о.з

0.6

0.9

1.2

Рис. 8: а) - зависимость скорости дислокации в Си от сдвиговых напряжений для различных концентраций вакансий: 1 - 0 %, 2 - 0.5 %, 3 - 1 %. Стрелками отмечено значение сгсг. б) - Зависимость критических напряжений от концентрации вакансий. Пунктирная линия - линейная аппроксимация.

Зависимость асг от концентрации вакансий имеет линейный вид (рис. 86). Теоретические модели дают зависимости acr ~ с1/3 или <тсг ~ с2/3 [13]. Возможно, что такое отличие связано с тем, что в теории рассматривается квазистатическое преодоление препятствий, тогда как на рис. 8 представлены значения в случае динамического преодоления.

В шестой главе рассмотрено влияние наноразмерных препятствий (полостей, включений) на пластические свойства. Исследуется роль препятствий как стопоров для дислокаций, в отличии от работы [14], где исследовалось влияния включений на упругие свойства веществ. Исследуются механизм отрыва дислокации от препятствия, критические напряжения и влияние температуры. Рассмотрены полости в Мо, полости и включения Си в AI.

В раздел 6.1 описывается модель для исследования взаимодействия дислокации и препятствий. Модель аналогична модели для подвижности одиночной дислокации, обсуждаемой в разделе 4.1. Отличие заключается во включении препятствия в плоскости скольжения дислокации.

Раздел 6.2 посвящен исследованию механизмов взаимодействия. Под действием сдвиговых напряжений дислокация начинает двигаться и на некотором расстоянии захватывается препятствием. В результате этого энергия в системе падает, также могут изменить знак сдвиговые напряжения. По мере дальнейшей деформации незакрепленная часть линии дислокации продолжает движение, в результате чего происходит ее изгиб. В рассматриваемом случае краевая дислокация по мере поворота линии становится винтовой (рис. 9а). Дислокации, относящиеся к разным частям линии исходной краевой дислокации, имеют противоположный знак, происходит их аннигиляция и отрыв от препятствия. Более детальное рассмотрение обнаружило, что по мере сближения двух частей дислокации они меняют плоскость скольжения в разных направлениях (рис. 96).

В разделе 6.3 исследуется зависимость напряжений отрыва дислокации от радиуса препятствия и расстояния между ними (рис. 10).

Проводится сопоставление с формулой, полученной на основе теории упругости:

А.

<уст — —log

где А = /х • /27г, L - расстояние между полостями, г - радиус полости, В - свободный параметр. Показано, что формула правильно описывает зависимость,

2г В

1 +

2 г

(б)

Рис. 9: Момент перед отрывом дислокации в молибдене от полости, а) - вид сверху, б) - вид спереди. Показаны только атомы с дефектным локальным окружением.

Температура Т — 300А'.

Рис. 10: Зависимость пороговых напряжений, необходимых дислокации для прохождения препятствия, от радиуса препятствия г и расстояния между ними L. а) - полость в Mo. Т = 300К (синие точки), Т — 700К (красные точки). Сплошные линии - значения, рассчитанные по формуле (6) с А=18 Н/м и В=0.46 им (синие), Л=13 Н/м и В=0.26 нм (красные), б) - кластер меди в Al. Т = 300К (синие точки), Т = 800К (красные точки). Сплошные линии - значения, рассчитанные по формуле (6) с А—2.9 Н/м и В=0.4 нм (синие), А=5 Н/м and В-0.7 нм (красные).

0.4-1

0.30.2 o.i Н

ау,ГПа

• 1

2 3

•■> 4

(¡¡>

(HI)

Т,К

—I-1-1-1

200 400 600 800 1000

Рис. 11: Температурная зависимость динамического предела текучести. 1 - экспериментальные результаты монокристаллов А1 [1]; 2 - экспериментальные результаты сплава А1 £>16Т [15]; 3 - оценка для монокристалла; 4 - критические напряжения

для сплава А1 с Си (4%).

но предложенное выражение А = ц ■ /2-я дает в три раза меньшее значение, по сравнению с аппроксимацией.

В седьмой главе рассмотрены подходы, позволяющие на основе результатов МД определять макроскопические характеристики пластической деформации.

В разделе 7.1 на основе результатов моделирования подвижности дислокации в AI и взаимодействия дислокации с включениями Си сделана оценка динамического предела текучести монокристалла и дисперсионно-упрочнённого сплава (рис. 11).

Величину напряжения течения монокристалла можно рассчитать, если считать, что скорость пластического деформирования ё определяется подвижностью дислокаций: е = pmbvo■ Будем считать, что плотность подвижных дислокаций р,п постоянна в экспериментах, не зависит от температуры и условий деформирования. Тогда величина динамического предела текучести «ту как функция температуры определяется температурной зависимостью коэффициента торможения дислокаций В(Т):

Точные значения плотности дислокаций в условиях ударно-волнового сжатия неизвестны, однако видно, что экспериментально полученная температурная зависимость хорошо описывается зависимостью динамического коэффициента торможения дислокаций от температуры В{Т), если величину плотности дислокаций принять равной pm~l-10® см-2. Это значение близко к экспериментальным, если предположить, что плотность в ударной волне увеличивается на два порядка.

Для сплава Al-Cu оценка была сделана на основе формулы (6) для размеров включений 1.5 нм. Стрелки демонстрируют три эффекта: (i) - увеличение динамического предела текучести для монокристаллов с ростом температуры, (¡i) - упрочнение при введение включений, (iü) - разупрочнение сплава с ростом температуры.

Раздел 7.2 посвящен обсуждению многомасштабного подхода для описания пластической деформации. Многомасштабность заключается в описании на двух пространственно-временных уровнях: 1) Для описания элементарных процессов пластической деформации на атомистическом уровне используется МД моделирование. Исследуется поведение одиночной дислокации под действием сдвиговых напряжений. На основе анализа результатов определяются кинетические закономерности данного процесса. Данный метод позволяет описать движение дислокаций с учетом атомарного строения их ядер и рассмотреть влияние тепловых колебаний атомов. 2) Макроскопический уровень представляет собой континуальную модель механики сплошной среды с дислокациями. В модели формулируется самосогласованная система уравнений для описания упруго-пластической деформации металла и эволюции дислокационной подсистемы. Модель была разработана Красниковым B.C. и Майером А.Е. При рассмотрении динамики дислокаций используются результаты МД моделирования.

Такой подход позволяет значительно выйти за рамки пространственных и временных масштабов, доступных МД. В то же время использование при континуальном описании информации об элементарных процессах на атомистическом уровне, полученной с помощью МД, позволяет сократить число свободных параметров.

На основе предложенного подхода проведено численное исследование эволюции дислокационной подсистемы при ударно-волновом нагружении алюминиевой мишени. Исследовано поведение динамического предела текучести с ростом температуры, результаты хорошо согласуются с экспериментальными данными [16] (рис. 12).

Рис. 12: Скорость тыльной поверхности алюминиевой мишени: эксперименты [16] и расчеты (время отсчитывается с момента выхода на поверхность ударной волны), а) - Г=293К; б) - Г=923К.

Удалось воспроизвести следующие экспериментально наблюдаемые явления: возрастание динамического предела текучести с температурой в монокристаллах алюминия, ускоренная релаксация напряжений при упруго-пластическом переходе в области высоких температурах, затухание упругого предвестника по мере его распространения вглубь мишени.

Помимо этого, подход позволяет проследить за эволюцией дислокационной подсистемы. Были получены результаты об изменении плотности дислокации при прохождении ударной волны и затем волны релаксации. Рассчитаны профили скоростей дислокаций вдоль мишени в разные моменты времени. На фронте ударной волны дислокации достигают скорости 1000 м/с, что соответствует динамическому режиму движения дислокаций.

3 Основные результаты и выводы работы.

В работе проведено атомистическое моделирование механизмов и кинетики элементарных процессов пластической деформации (зарождение и движение дислокаций, взаимодействия с дефектами решетки). Проведено сопоставление с теоретическими формулами, основанными на теории упругости, и с экспериментальными результатами. В частности, получены следующие результаты:

1. Определена структура и энергия зародышей дислокационных петель в А1. Получена зависимость скорости зарождения от сдвиговых напряжений и температуры. Показано, что зависимость может быть аппроксимирована в рамках формулы Журкова при введении зависимостей и и V от Т.

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

3. Изучены механизмы влияния точечных дефектов решетки (вакансий) на подвижность дислокации. В Мо равномерно распределенные вакансии концентрации с приводят к возрастанию фононного трения В ~ с, связанного с увеличением осцилляций дислокационной линии. В Си коэффициент трения практически не меняется, увеличиваются критические напряжения асг ~ с.

4. В процессе преодоления полостей и преципитатов наблюдается локальное изменение типа дислокации и плоскости скольжения вблизи препятствия. Получена зависимость величины напряжений, необходимых для преодоления препятствия, от их размера и расстояния между ними.

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

4 Публикации автора по теме диссертации.

В реферируемых журналах:

Янилкин А.В., Жиляев П.А., Куксин А.Ю., Норман Г.Э., Писарев В.В., Стегайлов В.В. Применение суперкомпьютеров для молекулярно-динамического моделирования процессов в конденсированных средах. Вычислительные методы и программирование. 2010. Т. 11. С. 111116.

Kuksin A., Norman G., Stegailov V., Yanilkin A., Zhilyaev P. Dynamic fracture kinetics, influence of temperature and microstructure in the atomistic model of aluminum. International Journal of Fracture. 2010. V. 162. P. 127-136.

Красников B.C., Куксин А.Ю., Майер A.E., Янилкин A.B. Пластическая деформация при высокоскоростном нагружении алюминия. Многомасштабный подход. Физика твердого тела. 2010. Т. 52. С. 1295-1304.

Жиляев П.А., Куксин А.Ю., Стегайлов В.В., Янилкин А.В. Влияние пластической деформации на разрушение монокристалла алюминия при ударно-волновом нагружении. Физика твердого тела. 2010. Т. 52. С. 1508-1512.

Kuksin A.Yu., Norman G.E., Stegailov V.V., Yanilkin A.V. Molecular Simulation as a Scientific Base of Nanotechnologies in Power Engineering. Journal of Engineering Thermophysics. 2009. V. 18. P. 197-226.

Куксин А.Ю., Стегайлов В.В., Янилкин A.B. Атомистическое моделирование пластичности и разрушения нанокристаллической меди при высокоскоростном растяжении. Физика твердого тела. 2008. Т.50. С. 1984-1990.

Куксин А.Ю., Стегайлов В.В., Янилкин А.В. Молекулярно-динамическое моделирование динамики краевой дислокации в алюминии. Доклады академии наук. 2008. Т. 420. С. 467-471, представлена ак. Осипьяном Ю.А.

Стегайлов В.В., Янилкин А.В. Структурные превращения в монокристаллическом железе при ударно-волновом сжатии и растяжении. Исследование методом молекулярной динамики. ЖЭТФ 2007. Т.131. С. 1064-1072.

Kuksin A.Yu., Norman G.E., Stegailov V.V., Yanilkin A.V. Atomistic simulations of structure transitions and fracture in Fe and Al single crystals. Сотр. Phys. Comm. 2007. V. 177, P. 48.

В реферируемых сборниках статей:

Kuksin A.Yu., Yanilkin A.V. Atomistic simulation of dislocation nudeation in aluminum. In Shock Compression of Condensed Matter - 2009, Ed. by M.D. Furnish et al. American Institute of Physics, New York. 2009. P. 841-844.

Stegailov V.V., Kuksin A.Yu., Norman G.E., Yanilkin A.V. Molecular dynamic modeling of plasticity of Al and Al-Cu alloys under dynamic loading / ibid. P. 781-784.

Zhilyaev P.A., Kuksin A.Yu., Stegailov V.V., Yanilkin A.V. Influence of temperature on spall strength: atomistic simulation / ibid. P. 785-788.

Янилкин A.B., Куксин А.Ю., Стегайлов В.В. Атомистическое моделирование пластичности и разрушения нанокристаллической меди при высокоскоростном растяжении / Физико-химическая кинетика в газовой динамике. 2008. Т. 7. www.chemphys.edu.ru/pdf/2008-09-01-032.pdf

Куксин А.Ю., Норман Г.Э., Стегайлов В.В., Янилкин А.В. Описание разрушения монокристаллического А| при высокоскоростном деформировании на основе молекулярно-динамического моделирования / Там же. -031.pdf

Норман Г.Э., Куксин А.Ю., Стегайлов В.В., Янилкин А.В. Многоуровневое моделирование пластичности и разрушения металлов при динамическом нагружении / Там же. -009.pdf

Стегайлов В.В., Куксин А.Ю., Норман Г.Э., Янилкин А.В. Моделирование динамики краевых дислокаций и их взаимодействия с нановключениями в условиях высокоскоростного деформирования / Там же. -030.pdf

Yanilkin A.V., Kuksin A.Yu., Norman G.E., Stegailov V.V. Atomistic simulation of high strain rate plasticity, spall damage and fracture of nanocrystalline Си / In Shock Compression of Condensed Matter - 2007, Ed. by M.D. Furnish et al. American Institute of Physics, New York. 2007. V. 955. P. 347-350.

Norman G.E., Kuksin A.Yu., Stegailov V.V., Yanilkin A.V. Atomistic simulation of plasticity and fracture of crystalline and polycrystalline metals under high strain rate / ibid. P. 329-334.

Kuksin A.Yu., Norman G.E., Stegailov V.V., Yanilkin A.V. Modeling of Al crystal fracture under high-rate strain based on atomistic simulations / ibid. P. 317-320.

Stegailov V.V., Kuksin A.Yu., Norman G.E., Yanilkin A.V. Atomistic study of nanoprecipitates influence on plasticity and fracture of crystalline metals / ibid. P. 339-342.

В сборниках статей:

Стегайлов В.В., Куксин А.Ю., Стариков С.В., Янилкин А.В. Деформационные и прочностные свойства наноструктурных конструкционных материалов / Сборник трудов Международного форума по нанотехнологиям, Москва, 2008, с.460-461.

Kuksin A.Yu., Norman G.E., Stegailov V.V., Yanilkin A.V. The influence of temperature on plasticity of Al alloys: molecular dynamics study / In proceedings of Forth International Conference on Multiscale Materials Modeling (MMM-2008). P. 442-445.

Янилкин А.В., Куксин А.Ю. Моделирование движения дислокаций в монокристаллах алюминия./ Сб. 50-ой научной конференции МФТИ "Современные проблемы Фундаментальных прикладных наук", Москва-Долгопрудный, 2007, с. 175-177.

Янилкин А.В. Разрушение и структурные превращения кристаллического железа при ударнс волновом нагружении. Исследование методом молекулярной динамики. Физика экстремальных состояний вещества-2006, Черноголовка, с. 89.

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

[1] Г. И. Канель, В. Е. Фортов, С. В. Разоренов. УФН 177, 809 (2007).

[2] В. И. Альшиц В. /1. Инденбом. УФН 115, 3 (1975).

[3] N. V. Gribova, Yu. D. Fomin, Daan Frenkel, V. N. Ryzhov. Phys. Rev. E 79, 6 (2009).

[4] M. Ю. Гуткин, К. H. Микаелян, И. А. Овидько. ФТТ 43, 42 (2001).

[5] Дж. Хирт, И. Лоте. Теория дислокаций. Атомиздат, Москва (1972). 857 с.

[6] С. Н. Журков, B.C. Куксенко, В.А. Петров. ДАН СССР 259, 1350 (1981).

[7] J. A. Gorman, D. S. Wood, Т. Jr. Vreeland. Phys. Rev. В 40, 833 (1969).

[8] A. Hikata, R. A. Johnson, C. Elbaum. Phys. Rev. В 2, 4856 (1970).

[9] F. Sinclair. Phil. Mag. A 56, 15 (1987).

[10] P. M. Derlet, D. Nguyen-Manh, S. L. Dudarev. Phys. Rev. В 76, 054107 (2007).

[11] D. L. Olmsted, L. G. Hector Jr, W. A. Curtin, R. J. Clifton. MSMSE 13, 371 (2005).

[12] K. Tapasa, D. J Bacon, Yu. N Osetsky. MSMSE 14, 1153 (2006).

[13] Т. Судзуки, X. Ёсинага, С. Такеути. Динамика дислокаций и пластичность. Мир, Москва (1989). 294 с.

[14] С. А. Лурье, Н. П. Тучкова. Композиты и наноструктуры 2, 25 (2009).

[15] Г. В. Гаркушин, С. В. Разоренов, Г. И. Канель. ФТТ 50, 805 (2008).

[16] G. Kanel, S. Razorenov, К. Baumung, J. Singer. J. Appl. Phys. 90, 136 (2001).

Подписано в печать 23 сентября 2010 г. Объем 1,0 п. л. Тираж 100 экз. Заказ № 726 Отпечатано в Центре оперативной полиграфии ООО «Ол Би Принт» Москва, Ленинский пр-т, д.37

Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Янилкин, Алексей Витальевич

1 Введение

2 Обзор литературы

2.1 Основные свойства дислокаций

2.2 Экспериментальные результаты.

2.2.1 Ударно-волновые эксперименты.

2.2.2 Наноиндентирование

2.3 Атомистическое моделирование.

2.4 Мпогомасштабные подходы.

3 Гомогенное зарождение дислокаций

3.1 Энергия дислокационной петли.

3.1.1 Модуль сдвига и энергия дефекта упаковки.

3.1.2 Энергия петли частичной дислокации.

3.2 Моделирование гомогенного зарождения дислокаций

3.2.1 Молекулярно-динамическая модель.

3.2.2 Механизм гомогенного зарождения дислокации

3.2.3 Метод определения скорости зарождения

3.3 Результаты и обсуждение.

Введение диссертация по физике, на тему "Атомистические механизмы и кинетика пластической деформации металлов при высокоскоростной деформации"

Диссертация посвящена разработке теоретического описания механизмов и кинетики пластической деформации металлов при высокоскоростной деформации. Построена модель явления и исследованы процесс гомогенного зарождения дислокации, явление фононного торможения, влияние точечных дефектов ц препятствий на подвижность дислокаций. Использован метод молекулярной динамики (МД). Полученные результаты обобщены в виде аналитических соотношений.

Актуальность работы. Высокоскоростное нагружение встречается во многих технологически важных приложениях, включая высокоскоростную обработку и формовку материалов, создание ударостойких материалов, пробивание защитных оболочек, ударное повреждение авиационных и космических летательных аппаратов и конструкций. Механический отклик многих материалов (металлов, керамик, стёкол, полимеров) при высокоскоростном деформировании отличается от случая статических нагрузок. Так для многих металлов зависимость напряжения течения от скорости деформирования резко усиливается с превышением скорости деформирования ~ К)-3 — Ю4 с-1 [1,2]. Явление можно интерпретировать как следствие изменения механизма движения дислокаций.

При малых скоростях движения дислокации преодолевают препятствия в результате совместного действия приложенного напряжения и тепловых флуктуацпй. Вследствие этого увеличение температуры сопровождается понижением предела текучести материалов. Для деформирования с высокой скоростью необходимо приложить более высокие напряжения. При скорости деформирования больше некоторой пороговой действующие напряжения оказываются достаточными для обеспечения динамического преодоления препятствий (барьеры Пайерлса), без дополнительного вклада тепловых флуктуаций. При этом доминирующим механизмом торможения дислокаций становится перекачка энергии дислокации в колебания кристаллической решетки или, в зависимости от температуры, в электронную подсистему. В отличие от области термофлуктуационной подвижности, скорость дислокаций в динамической области падает с температурой в соответствии с увеличением гоготиости элементарных возбуждений. Поэтому при высоких скоростях деформирования для некоторых материалов наблюдается (аномальное) возрастание напряжения течения с увеличением температуры [1].

Применительно к высокоскоростному деформированию сплавов интерес представляет конкуренция механизмов тсрмофлуктуациониого преодоления упрочняющих препятствий и фононного трения дислокаций. При этом равномерно распределенные дефекты могут приводить к увеличению коэффициента трения. В результате этого зависимость динамического предела текучести о г температуры зависит как I от скорости деформирования, так и от микроструктуры. Влияние температуры на предел текучести валено в некоторых технологических процессах, где, вследствие быстротечности, деформация и разрушение обрабатываемого материала протекают практически в адиабатических условиях и сопровождаются возрастанием температуры. Рост температуры может приводить к падению предела текучести и образованию адиабатических сдвигов.

2) Развитие моделей для исследования подвижности одиночной дислокаций в монокристалле. Исследование скорости движения в зависимости от сдвиговых напряжений п температуры для ОЦК и ГЦК металлов.

1) Исследован механизм гомогенного зарождения дислокаций за счег локального сдвига одной группы атомов относительно другой иа вектор Бюргерса частичной дислокации. Процесс носит стохастический характер. Показано, что зависимость скорости зарождения от сдвиговых напряжений может быть аппроксимирована в рамках формулы Аррениуса при введении температурной зависимости параметров.

2) Исследована зависимость скорости движения одиночной дислокации в ОЦК и ГЦК кристаллах от сдвиговых напряжений. Для коэффициента фононного трения В, характеризующего наклон линейного участка, получена зависимость от температуры. Показано подобие зависимостей коэффициента трения и периода решетки от температуры.

3) Изучены механизмы влияния дефектов решетки на подвижность дислокации в кристаллах Мо, Си и А1. В ОЦК равномерно распределённые вакансии концентрации с не оказывают существенного влияния на величину критических напряжений <тсг. необходимых для движения дислокации. Наличие вакансий приводит к возрастанию трения В ~ с. В случае ГЦК коэффициент трения практически не возрастает, а увеличивается барьер для движения acr ~ с.

Положения, выносимые на защиту.

2. Зависимость скорости движения дислокаций от сдвиговых напряжений в динамическом режиме для ОЦК и ГЦК кристаллов иа примере Mo, А1 и Си. Температурные зависимости коэффициента фононного трения.

3. Характер влияния вакансий и атомов примесей на подвижность дислокаций.

4. Механизмы преодоления полостей и включений при движении дислокации в Mo, А1. Зависимость сдвиговых напряжений, необходимых для преодоления препятствий, от расстояния между ними и их размера.

Результаты диссертации докладывались на конференциях: МФТИ 200208; Уравнения состояния вещества (п. Эльбрус 2005-10); Advanced Problems of Mechanics (Санкт-Петербург 2007-09); Многоуровневые подходы в физической мезомеха-нике (Томск 2008); Аэрофизика и физическая механика классических и квантовых систем (Москва 2007); Проблемы физики ультракоротких процессов в сильнонеравновесных средах (Новый Афон 2005-09); Физика прочности и пластичности материалов (Самара 2009); 10-ые Забабахинские научные чтения (Снежинск, 2010); Параллельные вычислительные технологии (Уфа 2010); Nucleation theory and applications (Дубна 2009-10); Computational physics (Brussels 2007, Ouro Preto 2008); Shock Compression of Condensed Matter (Hilo 2007); Multiscale Materials Modeling (Tallahassee 2008); Joint U.S.-Russia Conference on Advances in Materials Science (Prague 2009); IUTAM Symposium on Dynamic Fracture and Fragmentation (Austin 2009); New Models and Hydrocodes for Shock Wave Processes (Paris 2010).

В главе 2 проводится краткий обзор состояния исследований пластического деформирования при высокоскоростном нагружении. Приводятся основные методы исследования и основные результаты. Рассмотрены методы моделирования пластической деформации, в том числе многомасштабный подход. Глава 3 посвящена гомогенному зарождению дислокаций. Исследуются стохастические свойства процесса зарождения и зависимость скорости зарождения от величины сдвиговых напряжений и температуры. В главе 4 исследуется подвижность одиночной дислокаций под действием сдвиговых напряжений в зависимости ог температуры. В главах 5, 6 рассматривается влияние дефектов па подвижность дислокаций. Глава 7 проводится расчет кинетики пластической деформации па основе полученпых результатов МД и проводится сравнение с экспериментальными результатами.

Глава 2

Обзор литературы

В данной главе будут кратко рассмотрены основные экспериментальные результаты по исследованию пластической деформации, более подробные обзоры будут сделаны для каждой последующих глав по-отдельности. Также в этой главе будет подробно обсужден метод молекулярной динамики, используемый в данной работе, будут кратко рассмотрены другие методы моделирования и многомасштабные подходы для описания пластической деформации. Прежде чем перейти к описанию экспериментальных результатов и теоретических подходов рассмотрим основные свойства и механизм перемещения дислокаций - основных переносчиков пластической деформации.

Заключение диссертации по теме "Физика конденсированного состояния"

Основные результаты представленной работы следующие:

1. Определена структура и энергия зародышей дислокационных петель в А1. Получена зависимость скорости зарождения от сдвиговых напряжений и температуры. Показано, что зависимость может быть аппроксимирована в рамках формулы Ар-репиуса при введении зависимостей энергии активации и активациониого объема от температуры.

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

3. Изучены механизмы влияния точечных дефектов решетки (вакансий) на подвижность дислокации. В Мо равномерно распределенные вакансии концентрации приводят к возрастанию фононного трения пропорционально концентрации, связанного с увеличением осцилляций дислокационной линии. В Си коэффициент трения практически не меняется, увеличиваются критические напряжения пропорционально концентрации.

5. Предложен подход, позволяющий использовать результаты МД расчетов в механике сплошной среды: скорость гомогенного зарождения и коэффициент фонон-ного трения.

Список источников диссертации и автореферата по физике, кандидата физико-математических наук, Янилкин, Алексей Витальевич, Долгопрудный

1. Капель Г. И., Фортов В. Е., Разоренов С. В. Ударные волны в физике конденсированного состояния // Успехи физических наук. — 2007. — Т. 177. — С. 809.

2. Алъшиц В. И., Пнденбом В. Л. Динамическое торможение дислокаций // Успехи физических наук. — 1975.— Т. 115,— С. 3.

3. Fienkel Y. I. // Z. Physics. — 1926. — Vol. 37. — P. 572.

4. Хирт Д., Лоте И. Теория дислокаций. — Москва: Атомиздат, 1972.

5. Kanel G. I., Razorenov S. V., Fortov V. E. Shock-wave phenomena and the properties of condensed matter. — New York: Springer, 2004.— P. 321.

6. Armstrong R. W., Arnold W., Zerilh F. J. Dislocation mechanics of copper and iron in high rate deformation tests // Journal of Applied Physics.— 2009.- Vol. 105, no. 2.— P. 023511. http://link.aip.org/lmk/JAPIAU/vl05/i2/p02351 l/sl&Agg=doi.

7. Schuh C. A., Mason J. K., Land A. C. Quantitative insight into dislocation nucleation from high-temperature nanoindentation experiments / / Nature materials. — 2005, — Vol. 4, no. 8. — Pp. 617-621.

8. Головин Ю. И. Наноиндентпрование и механические свойства твердых тел в субмикрообъемах, тонких приповерхностных слоях и пленках // Физика Твердого Тела. 2008. — Т. 50, № 12. - С. 2113-2142.

9. Wo Р. С., Zuo L. Ngan А. Я. W. Time-dependent incipient plasticity in ni3al as observed in nanoindentation // Journal of Materials Research.— 2005. —Vol. 20. no. 2, —Pp. 489-495.

10. Хеермаи Д. В. Методы компьютерного эксперимента в теоретической физике.—Москва: Наука, 1990.— Р. 176.

11. Frenkel D., Sm.it В. Understanding Molecular Simulation. From Algorithms to Applications. — San Diego: Academic Press, 2002. — P. 638.

12. Schlick T. Molecular Modeling and Simulation: An Interdisciplinary Guide. — New York: Springei-Verlag, 2002. — P. 620.

13. Baidakov V. G., Protsenko S. P. Singular point of a system of lennard-jones particles at negative pressures // Phys. Rev. Lett. — 2005. — Vol. 95. — P. 015701.

14. Байдаков В. Г., Проценко С. П. Спинодаль и линия плавления леппард-джонсовского кристалла при отрицательном давлении // ДАН.— 2005.— Т. 402, № 6. С. 754-758.

15. Жуховицкий Д. И. // Коллоидный журнал. — 2003. — Т. 65. — С. 480.

16. Formation of high density amorphous ice by decompression of ice'vii and ice viii at 135 k / C. McBride, C. Vega, E. Sanz, J. L. F. Abascal // ,/. Chem. Phys.— 2004,- Vol. 121, no. 23, —Pp. 11907-11911.

17. The range of meta stability of ice-water melting for two simple models of water / C. McBride, C. Vega, E. Sanz et al. // Molecidar physics. — 2005.— Vol. 103, no. 1,- Pp. 1-5.

18. Zvbov V. I., Caparica A. A. A statistical mechanical study of thermidynamic properties of solid sodium under pressures based on an effective interatomic potential // Int. J. of Modern Phys. B. — 2004. — Vol. 18, no. 15,- P. 2185.

19. Разрушение твердой пленки в результате действия ультракороткого лазерного импульса / С. И. Аписимов, В. В. Жаховский, Н. А. Иногамов и др. // Письма в ЖЭТФ. — 2003. — Т. 77, № 11.- С. 731-736.

20. Ivanov D. S., Zhigilei L. V. Combined atomistic-continuum modeling of short-pulse laser melting and disintegration of metal films // Phys. Rev. B.— 2003.— Vol. 68. P. 064114.

21. Daw M. S., Baskes M. I. Semiempirical, quantum mechanical calculation of hydrogen embrittlement in metals // Phys. Rev. Lett. — 1983. —Apr. — Vol. 50, no. 17, —Pp. 1285-1288.

22. Daw M. S., Baskes M. I. Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals // Phys. Rev. B. — 1984. — Jun. — Vol. 29, no. 12. Pp. 6443-6453.

23. Foiles S. M., Baskes M. I., Daw M. S. Embedded-atom-method functions for the fee metals cu, ag, au, ni, pd, pt, and their alloys // Phys. Rev. B. — 1986. — Jun. — Vol. 33, no. 12, — Pp. 7983-7991.

24. Interatomic potentials for monoatomic metals from experimental data a,nd ab initio calculations / Y. Mishin, D. Farkas, M. J. Mehl, D. A. Papaconstantopoulos // Phys. Rev. В.— 1999,- Vol. 59, no. 5. — Pp. 3393-3407.

25. Atomistic studies of segregation and diffusion in al-cu grain boundaries / X.-Y. Liu, X. Wei, S. M. Foiles, J. B. Adams // Appl. Phys. Lett. — 1998. Vol. 72, no. 13.— P. 1578.

26. Liu X.-Y., Ercolessi F., Adams J. Aluminium interatomic potential from density functional theory calculations with improved stacking fault energy // Modell. Simul. Mater. Sci. Eng. 2004. — Vol. 12, no. 4. — Pp. 665-670.

27. Analysis of semi-empirical interatomic potentials appropriate for simulation of crystalline and liquid al and cu / M. I. Mendelev, M. J. Kramer, C. A. Becker, M. Asta // Philosophical Magazine.— 2008, — Vol. 88, no. 12,— P. 1723 — 1750.

28. Belashchenko D. K. // Russian Journal of Physical Chemistry. — 2006. —Vol. 80, no. 5. — P. 758.

29. Белащепко Д. К. Компьютерное моделирование жидких и аморфных веществ. — Москва: МИСИС, 2005. — С. 408.

30. Germann Т., Kadau К. Trillion-atom molecular dynamics becomes a reality // International Journal of Modern Physics C.— 2008. — Vol. 19, no. 9.— Pp. 13151319.

31. General formulation of pressure and stiess tensor for arbitrary many-body interaction potentials under periodic boundary conditions. // The Journal of chemical ■physics.— 2009,— Vol. 131, no. 15.— P. 154107. http://www.ncbi.nlm.nih.gov/pubmed/205G8817.

32. Lebovitz J. L., Percus J. K., Verlet L. Ensemble dependence of fluctuations with application to machine computations // Phys. Rev. — L967. — Vol. 153, no. 1.— Pp. 250-254.

33. Kelchner C. L., Plimpton S. J. Hamilton J. C. Dislocation nucleation and defect structure during suiface indentation // Phys. Rev. В.— 1998.— Vol. 58.— P. 11085.

34. Faken D., II. J. Systematic analysis of local atomic structure combined with 3d computei graphics // Computational Materials Science.— 1994. — Vol. 2, no. 2.— Pp. 279-286.

35. Кикып A. Y., Yamlkm A. V". Kinetic model of fracture at high strain rates in the example of crystalline aluminum // Doklady Physics. — 2007. — Vol. 52, no. 4. — Pp. 186-190.

36. Судзуки Т., Ёсинага X., Такеути С. Динамика дислокаций и пластичность. — Москва: Мир, 1989. — С. 296.

37. Strength differences arising from homogeneous versus heterogeneous dislocation nucleation / H. Bei, Y. F. Gao, S. Shim et al. // Phys. Rev. B. 2008. - Vol. 77. -P. 060103.

38. Pop-in effect as homogeneous nucleation of dislocations during nanoindentation / D. Lorenz, A. Zeckzer, U. Hilpert et al. // Phys. Rev. В. 2003.- Vol. 67.— P. 172101.

39. Gvtkin M. Y., Ovid'ko I. A. Special mechanism for dislocation nucleation in nanomaterials // Applied Physics Letters.— 2006. — Vol. 88, no. 21, — P. 211901. http://link.aip.org/link/APPLAB/v88/i21/p211901/sl&Agg=doi.

40. Гуткин M. Ю., Овидько И. А. Гомогенное зарождение дислокационных петель скольжения в нанокерамиках // Физика Твердого Тела. — 2008. — Vol. 50, по. 4. — Pp. 630-638.

41. Nanodisturbances in deformed Gum Metal / M. Y. Gutkin, T. Ishizaki, S. Kuramoto, I. A. Ovidko // Acta Material™.— 2006,— Vol. 54, no. 9.— Pp. 2489-2499. http://linkinghub.elsevier.coni/retrieve/pii/S1359645406001029.

42. Gvtkin M. Y., Ovidko I. A. Homogeneous nucleation of dislocation loops in nanocrystalline metals and ceramics // Acta Materialia. — 2008. — Vol. 56, no. 7. — Pp. 1642-1649. http://liiikinghub.elsevier.com/retrieve/pii/Sl359645407008476.

43. Гуткин M. Ю., Микаелян К. H., О виды; о И. А. Зарождение и развитие частичных дислокаций несоответствия и дефектов упаковки в тонкопленочных гетероструктурах // Физика твердого тела. — 2001, — Т. 43, 1.— С. 42.

44. Бобылев С. В., Овидъко И. А. Зарождение дислокаций па аморфных межзе-реипых границах в деформируемых наиокерамиках // Физика Твердого Тела. — 2008. Т. 50, № 4. - С. 617-623.

45. Dislocation nucleation induced by a shock wave in a perfect crystal: Molecular dynamics simulations and elastic calculations / D. Tanguy, M. Mareschal, P. S. Lomdahl et al. // Phys. Rev. B. — 2003. — Vol. 68. P. 144111.

46. Cao В. Brmga E., Meyers M. A. Shock compression of monocrystalline copper: Atomistic simulations // Metallurgical and Materials Transactions A. — 2007. — Vol. 38A, no. 11.— Pp. 2681-2688.

47. Predictive modeling of nanoindentation-induced homogeneous dislocation nucleation in copper / T. Zhu, J. Li, K. J. Van Vliet et al. // Journal of the Mechanics and Physics of Solids.— 2004. — Vol. 52, — Pp. 691-724.

48. Quantifying the eaily stages of plasticity through nanoscale experiments and simulations / K. Van Vliet, J. Li, T. Zhu et al. // Physical Review B.— 2003.— Vol. 67, no. 10,—P. 104105. http://link.aps.org/doi/10.1103/PhysRevB.67.104105.

49. Mechanical instabilities of homogeneous crystals / J. Wang, J. Li, S. Yip et al. // Phys. Rev. В. 1995.-Vol. 52.- Pp. 12627-12635.

50. Лившиц E. M., Питасвский Л. П. Теоретическая физика. Т. V. Статистическая физика. — Москва: Наука, 1979. — С. 528.

51. Хи G. S. A. A. Homogeneous nucleation of dislocation loops under stress in perfect crystals // Philisophical magazine letters. — 2000. — Vol. 80, no. 9. — Pp. 605-611.

52. Tschopp M. ASpearot D. E.; McDowell D. L. Atomistic simulations of homogeneous dislocation nucleation in single crystal copper // Modelling Simul. Mater. Sci. Eng. 2007. - Vol. 15. - Pp. 693-709.

53. Куксин А. Ю. Стегайлов В. В., Янилкип А. В. Молекулярно-динамическое моделирование динамики краевой дислокации в алюминии // Доклады Академии Наук. 2008. - Т. 420, № 4. - С. 467-471.

54. Пластическая деформация при высокоскоростном иагружепии. Многомасштабный подход. / В. С. Красников, А. Ю. Куксин, А. Е. Манер, А. В. Япил-кин // Физика Твердого Тела. — 2010. — Т. 7.—С. 1295-1304.

55. Psakhie S. G., Zolnikov К. P., Kryzhevich D. S. Elementary atomistic mechanism of crystal plasticity // Physics Letters A. — 2007. — Vol. 367, no. 3. — Pp. 250-253.

56. Norman G. E., Stegailov V. V. Simulation of ideal crystal superheating and decay // Mol. Simul. — 2004. Vol. 30, no. 9. — Pp. 397-406.

57. Норман Г. Э., Стегайлов В. В. Гомогенная пуклеация в перегретом кристалле. Молекулярно-динамический расчет // Доклады Академии Наук. — 2002. — Т. 386, № 3,- С. 328.

58. Standards for molecular dynamics modelling and simulation of relaxation / A. Y. Kuksin, I. V. Morozov, G. E. Norman et al. // Mol. Simul.— 2005.— Vol. 31, no. 14-15.—Pp. 1005-1017.

59. Hikata A., Johnson R. A., Elbaum C. Interaction of dislocations with elections and with phonons // Phys. Rev. В. — 1970. —Vol. 2. — P. 4856.

60. Gorman J. A., Wood D. S., Vreeland T. J. Mobility of dislocations in aluminum // Phys. Rev. B. 1969. - Vol. 40. - P. 833.

61. Судзуки Т., Есинага X., Такеути С Динамика дислокаций и пластичность.— Москва: Мир, 1989. — С. 294.

62. Nmomiya Т. // J. Phys. Soc. Jpn.- 1974. — Vol. 36.- P. 399.

63. Daw M. S.; Foiles S. M., Bashes M. I. // Mater. Sci. Rep. 1992,- Vol. 9,-P. 251.

64. Osetsky Y. N., Bacon D. J. An atomic-level model foi studying the dynamics of edge dislocations in metals // Model Simul. Mater. Sci.Eng. — 2003. —Vol. 11.— P. 427.

65. Cai J., Ye Y. Y. Simple analytical embedded-atom-potential model including a long-range force for fee metals and their alloys // Phys. Rev. В. — 1996.— Vol. 54. P. 8398.

66. Compact and dissociated dislocations in aluminum: Implications for deformation / S. G. Siinivasan, X. Z. Liao, M. I. Baskes et al. // Phys. Rev. Lett.— 2005.— Vol. 94. P. 125502.

67. Dynamic properties of screw dislocations in cu: A molecular dynamics study / D. Mordehai, Y. Ashkenazy, I. Kelson, G. Makov // Phys. Пси. В.- 2003.— Vol. 67.-P. 024112.

68. Kresse G., Furfhmuller J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. // Physical review. B. — 1996. — Vol. 54, no. 16, — Pp. 11169-11186. http://www.ncbi.nlra.nih.gov/pubmed/9984901.

69. Kresse G., Joubert D. From ultrasoft pscudopotentials to the projector augmented-wave method // Physical Review B.— 1999.— Vol. 59, no. 3. — Pp. 11-19.

70. Finnis, Sinclair // Phil. Mag. A. — 1987,-Vol. 56. — Pp. 15-30.

71. Derlet P. M., Nguycn-Manh D., Dudarev S. L. Multiscale modeling of crowdion and vacancy defects in body-centered-cubic transition metals // Phys. Rev. B.— 2007. Aug. - Vol. 76, no. 5.- P. 054107.

72. Clouet E. The vacancy-edge dislocation interaction in fee metals: a comparison between atomic simulations and elasticity theory // Acta Materialia. — 2006,- Vol. 54, no. 13.- Pp. 3543-3552. http://linkinghub.elsevier.com/retrieve/pii/S1359645406002424.

73. RUSSELL K., BROWN L. A dispeision strengthening model based on differing elastic moduli applied to the iron-copper system // Acta Metallurgies— 1972,- Vol. 20, no. 7.— Pp. 969-974. http://linkinghub.elseviei-.coin/retrieve/pii/0001616072900910.

74. Wirth B. D. Bulatov V. V. de la Rubia T. D. Dislocation-stacking fault tetrahedron interactions in cu // Journal of Engineering Materials and Technology.— 2002,— Vol. 124, no. 3.— Pp. 329-334. http://link.aip.org/link/? JYT/124/329/1.

75. Bacon J., Kocks U. F., Scattergood R. 0. // Philos. Mag.— 2005.— Vol. 28.— P. 1241.

76. Hafano Т., Matsvi H. Molecular dynamics investigation of dislocation pinning by a nanovoid in copper ,// Phys. Rev. B. — 2005. — Vol. 72. — P. 094105.

77. Garkushin G. V., Razorenov S. V., Kanel G. I. Submicrosecond strength of the dl6t aluminum alloy at room and elevated temperatures // Physics of the solid state. — 2008. — Vol. 50. — P. 805.

78. Dynamic yield and tensile strength of aluminum single crystals at temperatures up to the melting point / G. I. Kanel, S. V. Razorenov, K. Baumung, J. Singer // J. Appl. Phys.— 2001.- Vol. 90. no. 1, — Pp. 136-143.

79. Косевич С. В. Динамическая теория дислокаций // Успехи Физических Наук. — 1964. Т. LXXXIV, № 4. - С. 579.

80. Колгатиь С. Н., Хачатуръянец А. В. // Теплофизика Высоких Температур. — 1982, — Т. 20, № 3, — С. 90-94.

81. Уилкинс M. JI. Расчет упруго-пластических течений // Вычислительные методы в гидродинамике / Под ред. Б. Олдера, С. Ферпбаха, М. Ротенберга.— Мир, Москва, 2007.— С. 212-263.

82. Talion J. L., Wolf enden A. Temperature dependence of the elastic constants of aluminum // J. Phys. Chem. Solids. — 1979. — Vol. 40. — P. 831.

83. Horstemeyer M. F. Dashes M. /., Plimpton S. J. // Acta mater.— 2001.— Vol. 49. Pp. 4363-4374.

84. Ananthakrishna G. // Acta mater. — 2007. — Vol. 440. — Pp. 113-259.

85. Ландау Л. Д., Лифшиц Е. M. Теоретическая физика. T. VII. Теория упругости.— Москва: Наука, 1987.— С. 248.

86. Свенссон Т. Образование дислокаций в чистом алюминии при квазистатическом и ударном нагружении // Ударные волны и явления высокоскоростной деформации металлов / Под ред. М. А. Мейерса, JI. Е. Мурра. — Металлургия, Москва, 1984,- С. 164-176.

87. Малыгин Г. А. // Успехи Физических Наук, — 1999. — Т. 169, № 9.— С. 979985.

88. Малыгин Г. А. // Физика Твердого Тела. — 2007. — Т. 49, № 6.— С. 961-982.

89. Малыгин Г. А. // Физика Твердого Тела. — 2005. — Т. 47, № 2. — С. 236-241.

90. Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. T.V. Статистическая физика, — Москва: Наука, 1976,— С. 584.

91. Яловец А. П. Расчет течений среды при воздействии интенсивных потоков заряженных частиц // Прикладная Механика и Техническая Физика. — 1997.— T. 1.-С. 151-166.

92. Asay J. R. 11 Int. J. Impact Engng.— 1997, —Vol. 20. — Pp. 21-61.