Исследование кооперативных явлений, связанных с наличием носителей локализованного свободного объема, на примере кристалла чистого алюминия тема автореферата и диссертации по физике, 01.04.07 ВАК РФ
Тихонова, Татьяна Александровна
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Новокузнецк
МЕСТО ЗАЩИТЫ
|
||||
2010
ГОД ЗАЩИТЫ
|
|
01.04.07
КОД ВАК РФ
|
||
|
На правах рукописи
ТИХОНОВА ТАТЬЯНА АЛЕКСАНДРОВНА
ИССЛЕДОВАНИЕ КООПЕРАТИВНЫХ ЯВЛЕНИЙ,
СВЯЗАННЫХ С НАЛИЧИЕМ НОСИТЕЛЕЙ ЛОКАЛИЗОВАННОГО СВОБОДНОГО ОБЪЕМА, НА ПРИМЕРЕ КРИСТАЛЛА ЧИСТОГО АЛЮМИНИЯ
Специальность 01.04.07 - физика конденсированного состояния
Автореферат диссертации на соискание ученой степени кандидата физико-математических наук
4856053
Барнаул - 2010
4856053
Работа выполнена в ГОУ ВПО «Алтайский государственный технический университет им. И.И. Ползунова и ГОУ ВПО «Кузбасская государственная педагогическая академия»
Научные руководители: доктор физико-математических
Официальные оппоненты:
наук,
профессор, [Неверов Валерий Владимирович
Заслуженный деятель науки РФ, доктор физико-математических наук, профессор Старостенков Михаил Дмитриевич
доктор физико-математических наук, профессор Плотников Владимир Александрович
доктор физико-математических наук, профессор Кащенко Михаил Петрович
Ведущая организация: Кемеровский государственный университет,
г. Кемерово
Защита состоится "23" ноября 2010 года в 14.00 час. на заседании диссертационного совета Д 212.004.04 Алтайского государственного технического университета им. И.И. Ползунова по адресу: 656038, г. Барнаул, пр. Ленина, 46.
С диссертацией можно ознакомиться в библиотеке Алтайского государственного технического университета им. И.И. Ползунова.
Автореферат разослан «22» октября 2010 г.
Ученый секретарь диссертационного совета, кандидат физико-математических наук
Романенко В.В.
Примечание: Отзывы на автореферат, заверенный гербовой печатью организации, просим посылать в 2-х экземплярах на адрес университета.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность проблемы. Теория дефектов кристаллической решетки объясняет механические свойства твердых тел, составляя основу физических и технических теорий пластичности, прочности, разрушения. К простейшим точечным дефектам кристаллической решетки относятся вакансии и межузельные атомы.
При образовании комплексов вакансий и межузельных атомов в материале, структура кристаллической решетки вблизи них искажается за счет смещений атомов в упругих полях дефектов, что вызывает локальное нарушение плотности. Области, окружающие межузельный атом или комплексы межузельных атомов, обладают повышенной локальной плотностью; вакансию или комплексы вакансий - пониженной локальной плотностью.
В процессе структурной релаксации внутри кристалла происходит перераспределение локальной плотности, завершающееся образованием новых дефектов кристаллической решетки. Подобные процессы, связанные с перемещением вещества внутри некой системы, принято называть массопереносом. Локальная перестройка внутренней структуры материала сопровождается излучением материалом механических волн, это явление получило название акустической эмиссии. Несмотря на значительный объем теоретических и экспериментальных исследований, связанных с явлением акустической эмиссии, вопрос о природе акустической эмиссии, сопровождающей пластическую деформацию, остается открытым. Основная сложность при экспериментальных исследованиях заключается в том, что релаксация кристаллической решетки происходит с высокой скоростью. Поэтому в данном случае актуальным является использование метода компьютерного моделирования.
Компьютерное моделирование, являющееся в настоящее время таким же признанным методом исследования как экспериментальный и теоретический метод, начало применяться в физике твердого тела с конца 50х годов 20 века. С его помощью на атомном уровне возможно исследование не только быстропротекающих процессов, как например, аннигиляция точечных дефектов, но и процессов более длительных по времени.
В первой части работы исследован процесс бездиффузионного массопереноса в кристалле алюминия, содержащем комплексы межузельных атомов и вакансий, с помощью компьютерного моделирования по методу молекулярной динамики. Показано, что процесс сопровождается генерацией упругих волн. Во второй части работы решена задача об излучении упругих волн при аннигиляции двух краевых дислокаций, исходя из предположения, что движущие дислокации перемещают массу.
Цель настоящей работы заключается в исследовании методом компьютерного моделирования релаксационных процессов и механизмов
массопереноса в кристалле алюминия, содержащем такие носители свободного и стесненного объема, как вакансии, межузельные атомы и их комплексы; построении аналитического решения для пластического сдвига в бесконечной среде с использованием представлений о сдвиговом переносе массы.
Для достижения поставленной в работе цели были сформулированы следующие задачи:
1. Построить дву- и трехмерную молекулярно-динамическую модель для исследования на атомном уровне релаксационных процессов в кристалле алюминия, содержащем точечные дефекты.
2. Изучить картины атомных смещений вблизи комплексов точечных дефектов на начальной стадии процесса релаксации и итоговые конфигурации кристаллической решетки после релаксации.
3. Исследовать процессы релаксации кристалла алюминия при различных стартовых расстояниях между парами внедренных точечных дефектов, а также влияние деформации сжатия/растяжения расчетного блока, нагревания расчетного блока и влияние стартовых конфигураций комплексов точечных дефектов на процесс релаксации.
4. Выявить механизмы массопереноса, оценить скорости данного процесса, энергетические характеристики.
5. Построить аналитическое решение для пластического сдвига в бесконечной среде, пользуясь представлениями о переносе массы.
6. Решить задачу о разделении упругого поля краевой дислокации на бездилатационную и дилатационную составляющие.
Научная новизна работы заключается в том, что:
1) при помощи метода молекулярной динамики впервые обнаружена стадийность процесса релаксации кристаллической структуры: на начальной стадии зафиксированы максимальные смещения атомов величиной в 0,8 параметра решетки; на последующей стадии наблюдалось увеличение числа направлений кооперативных смещений атомов.
2) показано, что при релаксации атомной структуры кристалла происходит перераспределение локальной избыточной (недостающей) плотности, что приводит к формированию новых комплексов точечных дефектов;
3) установлено, что этот процесс сопровождается генерированием ударного фронта и его последующей трансформацией в упругую продольную волну; показано, что время массопереноса совпадает со временем распространения ударного фронта, а трансформация ударной волны в упругую продольную вызывает эффект акустической эмиссии;
4) наделение структурного дефекта массовой характеристикой позволяет идентифицировать его как механический объект, который обладает импульсом и энергией, зависящими от стартовых условий, таких как радиус
обрезания, расстояние до свободной поверхности, расстояние между дефектами;
5) установлено, что внешний импульс трансформируется в упругие моды разной частоты, по которым распределяется его энергия.
Научная и практическая ценность диссертационной работы заключается в возможности использования развитых представлений в радиационном материаловедении с целью создания новых наноматериалов с заданными свойствами, а также с усовершенствованием свойств уже известных материалов подвергающихся различным экстремальным воздействиям в процессе эксплуатации. Кроме этого, результаты компьютерного моделирования могут быть использованы как демонстрационный материал, отображающий процессы, протекающие в кристаллических структурах, полезный для студентов, осваивающих курс физики и физического материаловедения.
На защиту выносятся следующие положения:
1. Процесс релаксации кристаллической структуры, содержащей такие носители свободного и стесненного объема, как вакансии, межузельные атомы и их комплексы, состоит из нескольких стадий: образование ударного фронта, его последующая трансформация в упругую продольную волну.
2. Зависимость интенсивности процесса релаксации от типа комплексов дефектов.
3. Структурная релаксация при наличии таких носителей свободного и стесненного объема, как вакансии, межузельные атомы и их комплексы в кристалле сопровождается выделением энергии в виде акустической эмиссии.
4. Зависимость параметров упругих волн от вида переносы массы: конвективного или кондуктивного.
Апробация работы. Результаты работы были доложены и обсуждены на следующих научных конференциях:
• 11 Всероссийская конференция студентов-физиков и молодых учёных, Екатеринбург, 2005.
• 12 Всероссийская конференция студентов-физиков и молодых ученых, Новосибирск, 2006.
• Всероссийская научно-практическая конференция «Наука в ВУЗе», Новокузнецк, 2007.
• XVII Петербургские чтения по проблемам прочности, Санкт-Петербург, 2007.
• 13 Всероссийская конференция студентов-физиков и молодых ученых, Ростов-на-Дону, 2007.
• IV Международная конференция «Микромеханизмы пластичности, разрушения и сопутствующих явлений» (МРБР - 2007), Тамбов, 2007.
• Международная научная школа-конференция «Фундаментальное и прикладное материаловедение», Барнаул, 2007.
• 14 Всероссийская конференции студентов-физиков и молодых ученых, Уфа, 2008.
• X Международная школа-семинар "Эволюция дефектных структур в конденсированных средах" (ЭДС - 2008), Бийск, 2008.
• V Международная научно-техническая школа-конференция «Молодые ученые - науке, технологиям и профессиональному образованию в электронике», Москва, 2008.
• Вторая Всероссийская конференция с международным интернет -участием «От наноструктур, наноматериалов и нанотехнологий к Наноиндустрии», Ижевск, 2009.
• XV Международная научно-практическая конференция студентов, аспирантов и молодых ученых, Томск, 2009.
• Международный симпозиум «Перспективные материалы и технологии», Витебск, 2009.
• I региональная научно-практическая конференция «Наноиндустрия Алтая 2009», Бийск, 2009.
• XVII Международная конференции «Физика прочности и пластичности материалов», Самара, 2009.
• Всероссийская конференция с элементами научной школы для молодежи «Новые материалы. Создание, структура, свойства - 2009», Томск, 2009.
• Научно-техническая конференция с международным участием «Ультрадисперсные порошки, наноструктуры, материалы: получение, свойства, применение. V Ставеровские чтения», Красноярск, 2009.
• International Conference «Fundamental Aspects of External Fields Action on Materials», Novokuznetsk, 2010.
• V (XXXVII) Международная научно-практическая конференция «Образование, наука, инновации - вклад молодых исследователей», Кемерово, 2010.
• Всероссийская конференция студентов, аспирантов и молодых ученых, Владивосток, 2010.
• XI Международная школа-семинар "Эволюция дефектных структур в конденсированных средах" (ЭДС - 2010), Бийск, 2010.
Публикации. Результаты работы отражены в 30 публикациях, пять из которых в журналах, включенных в список ВАК для публикации диссертационных работ.
Структура и объем работы. Диссертация состоит из введения, пяти глав, заключения и списка литературы из 131 наименования. Работа изложена на 145 страницах машинописного текста, содержит 1 таблицу и 86 рисунков.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обосновывается актуальность и практическая ценность выбранного направления исследований. Сформулирована цель диссертационной работы, описана научная новизна, практическая значимость. Приведены основные защищаемые положения и краткое содержание работы по главам.
В первой главе дан обзор современных представлений о природе акустической эмиссии, приведены результаты экспериментальных исследований пластической деформации в металлах, сопровождающейся излучением упругих волн.
Во втором разделе уделено внимание работам, в которых проведено решение задач с учетом массопереноса.
Третий раздел главы посвящен описанию современных методов компьютерного моделирования и основных этапов построения компьютерной модели. Отмечено, что значительная часть исследований проблемы образования и последующей эволюции дефектной кристаллической структуры проведена при помощи компьютерного моделирования.
Во второй главе описывается построение компьютерной модели по методу молекулярной динамики. Сущность метода молекулярной динамики заключается в том, что поведение заданной совокупности атомов описывается в рамках классической механики системой обыкновенных дифференциальных уравнений движения Ньютона. Система уравнений решается с помощью численных методов интегрирования дифференциальных уравнений. Температура расчетной ячейки задается через начальные скорости атомов в соответствии с распределением Максвелла. При этом начальные скорости атомов задаются равными по абсолютной величине, но со случайными направлениями. Полная кинетическая энергия соответствует заданной температуре, а суммарный импульс расчетной ячейки должен быть равен нулю.
Определение потенциальной функции, которая бы наиболее корректно описывала межатомное взаимодействие, является важнейшей задачей компьютерного моделирования в физики конденсированных сред. Обычно приходится выбирать между точностью потенциальной функции и ее простотой, так как от этих условий зависит скорость выполнения расчетов и их достоверность. В настоящей работе межатомное взаимодействие задавалось парным центральным потенциалом Морза.
Затем приводится описание результатов компьютерного эксперимента. Вначале в расчетный блок 50 на 50 атомов алюминия импульсно внедрялись две пары межузельных атомов или две пары вакансий на стартовых расстояниях от одного до 12 атомов. После включения процедуры релаксации в каждом из случаев наблюдалось возникновение ударных
стартовых смещений величиной порядка 0,034 нм, затем через 0,5пс времени эксперимента с помощью визуализатора атомных смещений фиксировалась трансформация ударного фронта в упругую продольную волну (рис. 1).
Рис. 1. Картины атомных смещений для расчетного блока с двумя парами межузельных атомов со стартовым расстоянием в 2 межатомных
расстояния в моменты времени 0,01 пс (а) и 0,5 пс (б). Смещения увеличены в 100 раз
б
■V —: Чл
' л ' г. . -
.уф Ъ», . .
ШМ
Для различных расстояний скорости распространения смещений (или скорости массопереноса) были различными, однако в каждом из случаев превышали скорость звука в алюминии от 1,1 до 1,3 раза. После проведения процедуры релаксации расчетный блок подвергался охлаждению, и с помощью визуализатора атомных рядов фиксировались итоговые картины (рис. 2).
Рис. 2. Итоговые структуры после релаксации двух пар межузельных атомов со стартовым расстоянием в два межатомных и
последующего охлаждения (наложение атомных рядов). Выделенные области (а) и (б) содержат ядра дислокаций, образующие дислокационный диполь
Наблюдалось образование дефектов: дислокационных диполей и дислокационных петель. При стартовом расстоянии в 12 межатомных расстояний между парами в итоге полученная структура имеет в одном из направлений (а) две дислокационных петли, разделенных тремя атомными плоскостями. Одна дислокационная петля образована цепочкой из 23 вакансий (рис. За) и определяет область с недостаточной локальной плотностью; вторая дислокационная петля образована цепочкой из 23 межузельных атомов и определяет область с избыточной локальной плотностью (рис. 36). Дислокационные петли образовались в верхней центральной части расчетного блока. В двух других направлениях имеем по одному дислокационному диполю (рис.3).
Рис. 3. Итоговые структуры после релаксации двух пар межузельных атомов со стартовым расстоянием в двенадцать межатомных и последующего охлаждения (наложение атомных рядов). Выделены области: а) дислокационная петля из 23 вакансий б)
дислокационная петля из 23 межузельных атомов
Для исключения влияния свободных границ на внутренние процессы блока были проведены аналогичные эксперименты в блоке с большим числом атомов (80x80). Увеличение числа атомов не повлияло на итоговые конфигурации структуры блока, однако внесло изменения в скорость протекания процесса релаксации: процесс проходил менее интенсивно, итоговые конфигурации образовались лишь спустя 20 пс. О меньшей интенсивности свидетельствовала и более низкая средняя температура в блоке.
Следующий блок аналогичных экспериментов проводился для случая пар из 4 межузельных атомов. Рост избыточной массы привел к увеличению скоростей массопереноса, появлению иных итоговых картин.
В каждом из аналогичных экспериментов для пар вакансий на старте наблюдалось незначительное увеличение температуры до 1-2 К. В области, содержащей вакансии, возникают незначительные смещения группы из 50 атомов вдоль различных направлений кристаллографических плоскостей. В дальнейшем возникают вторичные смещения соседних атомов, и к моменту времени 0,5 пс эксперимента возникает кооперативное движение атомов с небольшими амплитудами, трансформирующееся в упругую продольную волну, что фиксируется по возникновению кооперативных атомных смещений. Оценки скоростей распространения смещений дают значение меньше скорости звука. После проведения процедуры релаксации объединение пар вакансий в единый комплекс не наблюдалось.
Далее приведены результаты компьютерного эксперимента при подвержении расчетного блока различным воздействиям. Наличие в материале внедренных дефектов создает повышенную локальную плотность, деформация растяжения блока способствовала снижению локальной плотности. Блок размером 50 на 50 атомов, содержащий пары межузельных атомов подвергался деформации растяжения с е = 3% и 5 %. Аналогично блок с вакансиями подвергался деформации сжатия в 3 и 5% для увеличения локальной плотности.
В случае, когда расчетный блок с межузельными атомами подвергался деформации сжатия, процесс полной релаксации проходил менее интенсивно, и итоговые картины получались через 22 пс времени эксперимента. В каждом из случаев наблюдался массоперенос с образованием дефектов: дислокационных диполей и петель. Интересно, что схожие структуры получились при различных стартовых условиях. При релаксации блока с парами межузельных атомов на расстоянии 12 межатомных и при релаксации блока с парами межузельных атомов на расстоянии 2 межатомных с дополнительной деформацией сжатия. В обоих случаях конечный дефект - ДП из 23 межузельных атомов и ДП из 23 вакансий в одном из направлений - образовался в верхней части блока.
В отличие от процесса релаксации блока, содержащего пары вакансий, в котором было замечена устойчивость вакансионного комплекса, дополнительное воздействие в виде деформации сжатия оказало влияние на процесс. Процесс релаксации в случае относительной деформации 3% занял 20 пс, 5% - 30 пс, в итоге получены схожие картины.
Далее исследовалось поведение атомов блока размером 50*50, содержащего две пары межузельных атомов, и после процедуры релаксации подвергающегося нагреванию при температуре 500 К в течение 100 пс.
Область, содержащая дефекты, переместилась к правой границе блока, что свидетельствует о том, что сообщение дополнительной энергии способствует переносу массы внутри блока.
В случае, когда в блоке находились две пары вакансий, не удалось проследить влияние нагревания на структурную перестройку материала.
Третья глава посвящена описанию экспериментов, когда в двумерный и трехмерный расчетный блок внедрялись комплексы из 3 и 10 дефектов, расположенных в виде равностороннего треугольника; из 7 и 13 - в виде шестиугольника.
Как и в случае меньшего числа межузельных атомов на старте эксперимента возникают смещения атомов, дискретные по направлению, и их амплитуды большие, поэтому по таким характеристикам подобные смещения можно отнести к смещениям, возникающим при ударной волне.
Со временем число направлений начинает возрастать, и увеличиваются амплитуды смещений. Уже через 0,5 пс в блоке, содержащем три межузельных атома на старте (рис. 4а) можно наблюдать возникновение коллективных смещений атомов, которые образуют упругую продольную волну со сферическим фронтом (рис. 46).
В течение эксперимента динамическая температура колеблется в пределах от 70 К до 80 К. Итоговая конфигурация образуется спустя 5 пс от начала эксперимента, и в дальнейшем не меняется.
ШШШШШ^ШШШШ ШШ^Шх^ЩШЗШ^':
ШШ Ш£:
■. .-л^-л-л-.у ■ -. ■ ■ ■ • - -¿с . . . л
а б
Рис. 4. Картины атомных смещений для расчетного блока, содержащего комплекс из трех межузельных атомов в моменты времени 0,1 пс (а) и 0,65 пс (б) после начала процедуры релаксации. Смещения увеличены в 10 раз
Как видно из рис. 5, происходит изменение стартовой конфигурации: заданный на старте треугольник межузельных атомов трансформируется в дислокационный диполь. Диполь расположен в двух из плотноупакованных направлений типа <110>и <112>. Таким образом, данная конфигурация не является устойчивой.
Рис. 5. Итоговая конфигурация, полученная в результате процесса релаксации расчетного блока, содержащего комплекс межузельных атомов в виде равностороннего треугольника Выделенные области (а) и (б) содержат ядра дислокаций, образующие дислокационный диполь
Заданный на старте треугольник вакансий трансформируется в дислокационную петлю, ядро которой состоит цепочки трех вакансий. Петля расположена в одном из плотноупакованных направлений типа <110>. Таким образом, данная конфигурация не является устойчивой, в процессе релаксации трансформируется в одномерную цепочку. Остальные конфигурации оказались устойчивыми к деформации, например комплекс из семи вакансий, образующий шестиугольник. Данный комплекс может считаться как дислокационный диск, лежащий в плоскости <111> или микропора. Сравнивая начальную и конечную конфигурации комплекса, можно отметить, что в результате релаксации происходит деформация всестороннего сжатия поверхности, ограничивающей комплекс вакансий. В результате этого объем (площадь), занимаемая комплексом, становится равной 0,97 от площади (объема) первоначального комплекса.
В третьей главе приведены результаты эксперимента, который проводился на расчетном блоке, имитирующем трехмерный кристалл
чистого алюминия. Взаимодействие между атомами ограничивалось пятью первыми координационными сферами. В центр расчетного блока внедрялись точечные дефекты. Продольная волна создавалась путем сообщения трем центральным атомам крайней плоскости некоторой скорости вдоль направления <111>. После конструирования расчетного блока включалась процедура релаксации путем разогрева кристалла вблизи О К. Затем при помощи визуализатора распределения потенциальной энергии рассматривались образующиеся после релаксации конфигурации точечных дефектов.
Кроме этого, в работе также использовался визуализатор атомных смещений, представляющий собой линии, соединяющие начальное и конечное положение атомов. При конструировании расчетного блока, дефекты располагались в плоскости {111} так, как показано на рисунке.
Вначале был рассмотрен процесс релаксации кристаллической структуры, содержащей комплекс вакансий, без какого либо воздействия на нее. Так на рис.6а показано распределение энергии атомов в расчетном блоке после процедуры релаксации. Конфигурация атомов в данном случае представляет собой практически идеальный октаэдр, вершины которого ориентированы вдоль направлений <100>, <010> и <001>. При увеличении обрезки энергетических значений атомов, видно, что данный комплекс представляет собой объемную тривакансию (рис.66).
а б
Рис.6. Распределение энергии атомов в расчетном блоке (фрагмент) содержащем комплекс из трех вакансий. Не показаны атомы с энергией связи в интервалах -3.84...-3.28 эВ (а) и -3.39...-3.12 эВ (б) (энергия связи атома в идеальном кристалле -3.34 эВ).
Затем был исследован процесс релаксации структуры, содержащей межузельные атомы. Так при внедрении дефектов в октаэдрическую пустоту, в процессе релаксации межузельные атомы образуют объемный комплекс из трех краудионов, ориентированных вдоль направления <011>.
При внедрении межузельных атомов в энергетически менее выгодное положение - тетраэдрическую пустоту - все трое начинают равномерно удаляться друг от друга, и в результате наблюдается следующая картина.
Продольная волна, возбужденная в кристаллической структуре, вносит изменения в процесс релаксации. Исследование показало, что фронт волны, проходящий через комплекс вакансий, не оказывает на них никакого
влияния. Волна возбуждалась как в начале эксперимента, так и после трансформации вакансионного комплекса.
Совсем иная картина наблюдается в эксперименте с межузельными атомами. Прохождение фронта продольной волны через агрегат межузельных атомов приводит к тому, что в случае расположения их в октаэдрической пустоте образуется комплекс, представленный на рис.7а, а при внедрении их в тетраэдрическую пустоту - объемный краудионный
комплекс, ориентированный вдоль направления < 110 > (рис.76).
а б
Рис.7. Распределение энергии атомов в расчетном блоке (фрагмент) содержащем комплекс из трех межузельных атомов, внедренных в октаэдрическую (а) и тетраэдрическую (б) пустоты после прохождения продольной волны. Не показаны атомы с энергией связи в интервалах -3.82...-2.59 эВ (а) и -3.82. ..-2.64 эВ (б) (энергия связи атома в идеальном кристалле -3.34 эВ).
В четвертой главе описаны теоретические представления о приписывании движущимся дефектам массы, развитые В.В. Неверовым, и на основе решения плоской задачи теории упругости для области с разрезом показано разделение поля дислокации на дилатационную и бездилатационную части.
Пластические сдвиги и трещины перераспределяют массу среды, в которой они развиваются. Распределение массы определяется полями напряжений, связанными с названными объектами. Участки сдвига или трещины представляются линейными разрезами. Решения показывают, что плоскость с разрезом разбивается на области (далее м-области), знак дилатации в которых одинаков (рис. 8). Материал м-областей с отрицательный дилатацией сжат, массовая плотность в этих областях повышена, и, следовательно, они содержат избыточную массу. В областях с положительной дилатацией масса снижена. Масса этих областей частично переместилась в области с отрицательной дилатацией. Движение фронта области пластического сдвига или трещины сопровождается изменением распределения дилатации, а, следовательно, и перемещением массы.
Рис. 8. Линии уровней распределения массовой плотности в поле напряжений пластического незавершенного сдвига. Линии нулевого уровня проходят по осям координат. Области повышенной плотности помечены знаком плюс. Пунктиром показаны границы исходно квадратной части тела после деформации. Смешения границ увеличены.
В микроскопическом описании дефекты, участвующие в сдвиговой пластической деформации, представляют собой дислокации, плоские скопления дислокаций, дисклинации, двойники, прослойки мартенситных фаз, а при мезо- и макроскопическом описании - полосы, пачки, пакеты, стопы скольжения, сдвига или деформации. Во всех этих случаях процесс можно свести к движению дислокаций. Учет переноса массы позволяет присвоить массу дефектам, вызывающим движения. Массу удобно присваивать, например, краю участка сдвига или краю трещины. Разделение полей позволяет установить массу движения.
В общем случае, м-области образуются на неоднородностях функции распределения сдвиговых смещений по плоскости сдвига. У каждой неоднородности формируются две м-области. Они расположены симметрично по разные стороны от плоскости сдвига и имеют разные знаки отклонения массовой плотности (м-диполь). М-диполь сопровождает движение неоднородности, в частности, дислокации и переносит массу. С разных сторон от плоскости сдвига масса перемещается в противоположных направлениях. Масса, которая присваивается дислокации или фронту участка сдвига, равная избыточной (недостающей) массе м-областей, является массой покоя. Поля напряжений описываются вектором смещений u(x,y,z). Известно, что вектор всегда можно представить суммой двух частей:
U =И, +U, j-j-j
одна из которых удовлетворяет условию div й =0(2), а другая - условию
rot и, = 0(3). Из (2) следует, что поле ut описывает движения в несжимаемом
материале. В этом случае масса перемещается вместе с объемом, и такой
перенос является конвективным. Поле (7/ (3) связано с дилатацией. В этом
поле масса перемещается за счет образования и движения м-областей. Масса переносится сквозь объем, поэтому перенос является кондуктивным.
Возможность представления поля деформации в виде двух полей -без дилатации и с дилатацией - следует из функции напряжений Эри. В плоской теории упругости при отсутствии объемных сил напряжения и смещения могут быть выражены через одну вспомогательную функцию U{z) = xp(z) + yq(z) + p](z)^ ^
где р! - гармоническая функция, и 2 — х + ¡у. Первые два слагаемые в (4) представляют гармоническую функцию, обусловленную величиной пропорциональной гидростатическому давлению или дилатации: р(г) + /я(2) = (1/4)][Р(г) + /е(г)]Л, Р{£) = ахх(2) + ауу(2)
здесь СКг) - функция комплексно-сопряженная с Р(г). Сумма этих двух слагаемых
V, =Х1^)+}ф)
является функцией Эри дилатационной части поля деформации. Тогда р1 (г) - является функцией Эри бездилатационного поля.
Так как функции IIх, и р{ гармонические, то каждая из них удовлетворяет условию равновесия и описывает реально возможное поле деформаций в упругой среде. В сумме эти функции удовлетворяют граничным условиям.
Дилатационное и бездилатационное поля найдены для ряда сдвиговых дефектов (рис. 9), а так же для трещины отрыва (рис. 10). Отметим, что бездилатационное поле петли краевой дислокации так же состоит из трех вихрей. Сравнение модулей векторов смещений по размеру точек на рисунке показывает, что у двустороннего скопления дислокаций преобладают движения, обусловленные дилатационной частью. У трещины отрыва движения того и другого поля соизмеримы.
/
¡"г ! ;
<А -: 1 4
Рис. 9. Линии вектора смещения дилатационной части поля двустороннего скопления краевых дислокаций -слева (для малых точек < о 2> Для самых больших | > р ^ ) и для бездилатационной части - справа
(для малых точек |й | < 0 03 Участок сдвига показан утолщенной линией на правом рисунке
Формула (7) записана для конвективного переноса и применима для бездилатационного поля. С дилатационным полем дислокации связаны два движения: первое описывает перемещения элементов объема, а второе -перемещение м-диполя как целого. Массу переносит второе движение. Процесс является кондуктивным. Масса кондуктивного движения, близкая по величине массе экстраплоскостей, для грубых оценок принимается на шесть порядков меньше массы деформируемого блока. Тогда при одинаковых энергетических затратах на движение скорость перемещения м-областей
должна быть на три порядка больше, чем при конвективном движении. В соответствии с этим при прочих равных условиях частота волн, порождаемых кондуктивным переносом, примерно, на три порядка выше, чем при конвективном. Так как у трещин отрыва модули векторов смещения бездилатационного и дилатационного движений близки, то, если дислокационный сдвиг порождает волны высокой частоты, то в спектре трещин отрыва, кроме того, имеется еще и значительная низкочастотная часть.
Кондуктивное движение перемещает массу на длину всего деформируемого блока, что создает большие скорости в правой части второго уравнения (7). В результате масса движения краевой дислокации превосходит массу покоя.
Вид переноса определяет: величину и скорость перемещаемой массы; места, в которые перемещается масса, а, следовательно, высвобождаемая энергия и импульс. Такие места являются источниками упругих волн. Перенос массы определяет развитие динамических процессов пластической деформации: изменение формы зерен, а, следовательно, их взаимодействие при деформации поликристаллов, образование трещин, параметры генерируемых упругих волн и т. д.
Рис. 10. Линии вектора смещения дилатационной части поля трещины отрыва - слева и бездилатационной части - справа (для тонких линий [и, | < 0 15'для самых толстых [ > о 25 ^
В пятой главе, исходя из представления о сдвиговом переносе массы, описан процесс аннигиляции двух краевых дислокаций противоположного знака, расположенных в бесконечной изотропной среде. В среде выделяется квадратный блок, внутри которого находится пара краевых дислокаций противоположного знака, расположенных антисимметрично. Под действием сил взаимного притяжения происходит сближение дислокаций, и затем - их аннигиляция. Дислокации оказывают влияние на упругие поля друг друга, поэтому при их сближении масса, переносимая каждой из них, уменьшается, однако, увеличивается масса, выходящая на границу раздела дислокаций. Для блока, состоящего из материала единичной плотности размером 2000Ь-2000Ь, длина экстраплоскости для каждой из дислокаций составляет ЮООЬ. Вектор Бюргерса принят равным единице (Ь=1). Масса, выносимая на границу аннигиляции, рассчитывается как произведение величины вектора Бюргерса
на длину экстраплоскости. Будем следить за потоками массы (1гор,п - через
горизонтальные стороны блока и .Гвертик - через вертикальные), величиной
массы м-области М и массой, выходящей на границу аннигиляции Мгра11 (8):
и
к(0 ,у,1Му
(8)
Формулы для смещений взяты из решения второй плоской задачи теории упругости. Каждая из этих масс рассчитывается для четверти блока, а затем увеличивается в 4 раза, поскольку массы в соответствующих четвертях имеют разные знаки, и их алгебраическая сумма равна нулю.
График изменения суммарных масс в зависимости от расстояния дислокаций (в векторах Бюргерса) от границы аннигиляции приведен на рис. 11.
м -- Рис. 11. Графики зависимости суммарных
г потоков через горизонтальные стороны блока
"" .. - .[гориз, через вертикальные стороны ./вергик,
_"7"4----- --масса м-области М и масса, выходящая на
границу аннигиляции Мгран от I. (расстояния : краевых дислокаций от границы
■/■•У............ у аннигиляции). Сумма этих масс Ъ в точности
--. равна массе экстраплоскостей КД, вводимых в
^ " материал.
По мере сближения дислокаций вынесенная на границы блока масса, а также масса м-областей перемещается на границу аннигиляции, следовательно, Мгран возрастает. В момент аннигиляции вся масса вынесена на границу. Как видно из рисунка, сумма Е всех масс величина постоянная и в точности равна 2000 (в отн.ед.), т.е. массе экстраплоскостей КД, вводимых в материал.
Полагаем движение дислокаций ускоренным. Считаем, что внешние силы работы не совершают (релаксационный процесс). Уменьшающаяся упругая энергия полностью расходуется на увеличение относительной скорости дислокаций. Затем в теле выбирается поверхность 2 так, чтобы источник Б лежал внутри этой поверхности, а источник Т - вне её. Выбор поверхности 2 произволен, если она разделяет системы внутренних напряжений и лежит в области одновременного существования смещений и® и и,т. Тогда энергия взаимодействия определяется следующим поверхностным интегралом:
2( (9)
Энергия взаимодействия систем напряжений или работа, которую совершают напряжения одной системы напряжений над перемещениями другой, вычисляется для границы, которая разделяет источники систем напряжений. В рассматриваемой модели такой границей является плоскость
аннигиляции (х = 0), а системы напряжений обусловлены наличием поля
мЧ4,ч
/ Л~ м
сдвига (поля напряжений, создаваемого аннигилирующими дислокациями) и полем упругих волн. Эти поля сосуществуют во время аннигиляции и сменяют друг друга.
Слагаемые в подынтегральном выражении формулы (9) равны (работа источника S над перемещениями Т равна работе источника Т над перемещениями S), то можем записать формулу для работы:
А = \axur(<.),y)Jy
(10)
где оху - касательные напряжения, создаваемые дислокациями, иу - смещения точек плоскости аннигиляции.
Тогда для работы, которую совершает этот импульс (S) над генерируемыми волнами (Т), запишем:
А= \dA= \Fsdu = fF\Tdt= \vTdJs
J J J J (11)
По мере возбуждения волна будет распространяться вдоль оси Оу и отводить энергию от места генерации. На время аннигиляции (время сосуществования полей двух источников) будем пренебрегать отводом энергии, тогда за смещения Т можем принять уравнение стоячей волны. Стоячие волны всех частот имеют в точке привязки нулевую амплитуду -узел, это условие учтено использованием функции sin. Тогда выражение для смещений Т имеет вид:
ит(y,t) = и/sinwtsin ку
и для скорости:
vT(y,t) = wu0T cosw/sinAj'
Составим выражение для работы:
л '
г, 0
Связь между периодом и частотой, периодом и длиной волны, длиной волны и волновым числом определяемся формулами:
1ж , 1л
и> = — к- —
с'\ (15),
где - й - скорость упругих продольных волн.
Выбирая период, мы задаем через соотношения (15) все волновые характеристики. Волны максимальной частоты, генерируемые при аннигиляции, лежат в пределах дебаевских частот. На генерацию каждой гармоники тратится определенное количество работы, и, исходя из принципа максимальной работы, можно предсказать, волна какой частоты (длины волны) генерируется (табл. 1). Таким образом, работа, совершаемая над волнами, и параметры генерируемых волн при аннигиляции дислокаций зависят от геометрических свойств образца.
* (14)
ОСНОВНЫЕ ВЫВОДЫ
1. Процесс релаксации кристаллической структуры, содержащей такие носители свободного и стесненного объема, как вакансии, межузельные атомы и их комплексы, протекает с формированием ударного фронта и его последующей трансформацией в упругую продольную волну.
2. Стадия перехода от ударных смещений атомов к возникновению упругих продольных волн увеличивается по времени с ростом числа межузельных атомов в комплексе, т.е. с ростом избыточной локальной массы со знаком «+».
3. Для каждого из комплексов стадия перехода от ударных смещений атомов к возникновению упругих продольных волн уменьшается по времени с ростом числа вакансий в комплексе, т.е. с ростом локальной массы со знаком «-».
4. Время массопереноса совпадает со временем распространения ударного фронта, а его трансформация в упругую продольную волну вызывает эффект акустической эмиссии.
5. Даже при низкоэнергетическом воздействии на материал, когда не образуются точечные дефекты, образующаяся при этом продольная волна способна оказывать влияние как на процесс формирования комплекса дефектов, так и на уже сформировавшийся комплекс.
6. Процесс релаксации кристалла завершается формированием устойчивых структур атомов блока, свидетельствующих о завершении массопереноса.
7. Разработан метод расчета параметров акустических волн, генерируемых при аннигиляции краевых дислокаций.
8. Предложенный способ разделения упругого поля на дилатационную и бездилатационную часть дает возможность выделить в процессе деформации чисто конвективный и кондуктивный перенос.
Таблица 1. Расчет параметров воли, генерируемых при аннигиляции дислокации в алюминии___
Вещество Радиус зерна Значения констант материала и множителей для частных случаев (р, р, Ь) в ед. СИ
Алюминий, 2700,4- Ю10 4,2-10"'0 Алюминий, 2700,4-Ю104,2-10'8
Частота, V, с"1 Длина волны, X, м Частота, V, с"1 Длина нолны, X, м
0,3 мкм 2,96-10' 1,93'Ю"6 - -
0,6 мкм 2,02-10" 3,18'Ю"6 - -
1 мкм 1,34-10'' 4,4610'' - -
3 мкм 4,70-10" 1,27-10"5 - -
6 мкм 2,29-108 2,4810"5 7.80-107 7.6810'5
10 мкм 1,34-10" 4,46-Ю5 4,ЗОЮ7 1,3910А
30 мкм 4,7010' 1,27-Ю-4 2,96 107. 2,0210"*
ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ
1. Тихонова Т.А. Зависимость массы, переносимой краевой дислокацией, от расстояния до свободной поверхности // Материалы научно-практической конференции студентов / Сборник статей под редакцией М.С. Можарова. - Новокузнецк: Изд-во КузГПА, 2004. - С.125-128.
2. Тихонова Т.А. Модель выпрыгивания краевой дислокации на свободную поверхность // Сборник тезисов 11 Всероссийской конференции студентов физиков и молодых учёных. -Екатеринбург, 2005.-С. 118-121.
3. Тихонова Т.А. Модель аннигиляции краевых дислокаций // Сборник тезисов 12 Всероссийской конференции студентов-физиков и молодых ученых. - Новосибирск, 2006. - С. 112-114.
4. Тихонова Т.А. Метод расчета параметров упругих волн, генерируемых при аннигиляции краевых дислокаций // Наука в вузе: современные тенденции. Материалы Всероссийской научно-практической конференции. Часть I. - Новокузнецк: Изд-во КузГПА, 2007. - С. 24-29.
5. Тихонова Т.А., Неверов В В. Модель излучения упругих волн при аннигиляции краевых дислокаций // XVII Петербургские чтения по проблемам прочности. Санкт-Петербург, 10-12 апреля 2007 г.: сборник материалов. - Ч. 2. - СПб., 2007. - С. 94-96.
6. Тихонова Т.А.Генерация упругих волн при аннигиляции краевых дислокаций // Сборник тезисов 13 Всероссийской конференции студентов-физиков и молодых ученых. - Ростов-на-Дону, 2007. - С. 143-144.
7. Тихонова Т.А., Неверов В.В., Маркидонов A.B. Приведенная масса как мера инертности волновых движений // «Микромеханизмы пластичности, разрушения и сопутствующих явлений». Материалы IV Международной школы-конференции. Тамбов, 2007. С. 196-197.
8. Маркидонов A.B., Тихонова Т.А., Неверов В.В. Перемещение массы пластическими сдвигами // Фундаментальные проблемы современного материаловедения. 2007. т.4. №2. С.41-46.
9. Тихонова Т.А., Маркидонов A.B., Неверов В.В. Особенности кондуктивного переноса массы при деформации твердых тел // Фундаментальные проблемы современного материаловедения. 2007. т.4. №4. С.79-82.
10. Маркидонов A.B., Тихонова Т.А. Свойства массы движения пластических сдвигов // «Молодые ученые - 2008». Материалы V Международной научно-технической школы-конференции. Москва, 2008. ч.1. С.64-67.
11. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А. Нелинейный высокоскоростной массоперенос в двумерном кристалле при наличии локальных областей с различной плотностью // Ползуновский альманах. 2008. № 3. С.226-228.
12. Тихонова Т.А. Разделение полей сдвиговых дефектов // Электроника Сибири. - 2008. -№3.-С. 141-143.
13. Маркидонов A.B., Тихонова Т.А., Старостенков М.Д. Влияние начальной температуры на скорость кооперативных атомных смещений // «От наноструктур, наноматериалов и нанотехнологий к Наноиндустрии». Тезисы докладов второй всероссийской конференции с международным интернет - участием. Ижевск, 2009. С. 75.
14. Маркидонов A.B., Тихонова Т.А., Старостенков М.Д. Исследование массопереноса в двумерном нанокристалле чистого AI // «Современные техника и технологии». Труды XV Международной научно-практической конференции студентов, аспирантов и молодых ученых. Томск, 2009. т.З. С.510-512.
15. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А., Медведев H.H. Высокоскоростной перенос массы посредством краудионных столкновений // «Перспективные материалы и технологии». Сборник тезисов международного симпозиума. Витебск, 2009. С.55.
16. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А., Медведев H.H. Примеры высокоскоростного переноса массы в идеальных монокристаллах без разрушения его структуры И «Наноиндустрии Алтая 2009». Тезисы докладов I региональной научно-практической конференции. Бийск, 2009. С.33-35.
17. Старостенков М.Д., Маркидонов A.B., Медведев H.H., Тихонова Т.А. Моделирование переноса массы в виде рядов вакансий и межузельных атомов // «Физика прочности и
пластичности материалов». Тезисы докладов XVII международной конференции. Самара, 2009. С.37.
18. Маркидонов A.B., Тихонова Т.А., Старостенков М.Д. Моделирование массопереноса на примере объемного ГЦК кристалла // «Новые материалы. Создание, структура, свойства - 2009». Труды Всероссийской конференции с элементами научной школы для молодежи. Томск, 2009. С.268-272.
19. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А. Высокоскоростной массоперенос в кристалле при наличии различных конфигураций точечных дефектов // Фундаментальные проблемы современного материаловедения. 2009. т.6. №1. С.12-16.
20. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А., Медведев H.H. Высокоскоростной массоперенос в двумерном кристалле никеля при наличии дислокационных петель различной локальной плотности // Изв. вузов. Черная металлургия. 2009. №6. С.57-60.
21. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А., Потекаев А.И., Кулагина В.В. Высокоскоростной массоперенос в кристаллическом алюминии, содержащем цепочки вакансий и межузельных атомов // Изв. вузов. Физика. 2009. т.52. №9/2. С.139-145.
22. Маркидонов A.B., Тихонова Т.А., Неверова Т.Н., Старостенков М.Д. Моделирование незавершенного массопереноса на примере двумерного нанокристалла // «Ультрадисперсные порошки, наноструктуры, материалы: получение, свойства, применение. V Ставеровские чтения». Труды научно-технической конференции с международным участием / Под ред. В. Е. Редькина. Красноярск, 2009. С.13-15.
23. Starostenkov M.D., Pozidaeva O.V., Tikhonova Т.А., Poletaev G.M., Medvedev N.N. The mechanisms of interaction and annihilation of the aggregates of vacancies and interstitial atoms in a two-dimensional lattice of an ordered alloy of LI 2 superstructure // Fundamental Aspects of External Fields Action on Materials (Book of the International conference articles), 26-28 may, 2010, Advanced Materials Institute Graduate School at Schezhen Tsingua University, Ed. By V. Gromov, Novokuznetsk. - C. 125-135.
24. Неверов B.B., T.A. Тихонова, Маркидонов A.B., Старостенков М.Д. Масса, переносимая при скольжении краевых дислокаций // Известия Алтайского государственного университета. -2010. -№1(65). -С. 146-150.
25. Старостенков М.Д., Пожидаева О.В., Тихонова Т.А., Полетаев Г.М., Медведев H.H. Механизмы взаимодействия и аннигиляции агрегатов вакансий и межузельных атомов в двумерной решетке упорядоченного сплава сверхструктуры L12 // Известия Алтайского государственного университета.-2010.-№ 1(65 ).-С. 173-176.
26. Старостенков М.Д., Маркидонов A.B., Медведев H.H., Тихонова Т.А. Моделирование переноса массы в виде рядов вакансий и межузельных атомов на примере двумерного кристалла //Вестник Самарского гос. техн. ун-та. Сер. Фнз.-мат. науки. -2010. - вып. 1(20). - С. 249-252.
27. Тихонова Т.А., Маркидонов A.B. Наблюдение релаксации комплексов межузельных атомов в металлах с ГЦК-решеткой // Образование, наука, инновации - вклад молодых исследователей: материалы V (XXXVII) Международной научно-практической конференции / Кемеровский госуниверситет. - Кемерово: ООО «ИНТ», 2010. - Вып. 11. - Т. 2. - С. 625-627.
28. Тихонова Т.А., Маркидонов A.B. Релаксационные процессы в кристалле, содержащем различные комплексы вакансий // Тезисы докладов Всероссийской конференции студентов, аспирантов и молодых ученых по физике. Владивосток, 2010. - С. 107-108.
29. Старостенков М.Д., Тихонова Т.А. Исследование механизмов структурной перестройки двумерного кристалла алюминия, содержащего комплексы межузельных атомов // Фундаментальные проблемы современного материаловедения. 2010. т.7. №2. С.55-60.
30. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А. Влияние низкоэнергегических воздействий на релаксационные процессы в ГЦК кристалле (в печати).
Издано в авторской редакции.
Подписано в печать 19.10.10. Формат 60x84 1/16.
Печать - ризография. Усл. п. л. 1,39.
Тираж 100 экз. Заказ 2010-15
Издательство Алтайского государственного технического университета им. И.И. Ползунова, 656038, г. Барнаул, пр-т Ленина, 46. Лицензия на издательскую деятельность ЛР№ 020822 от 21.09.98 г.
ВВЕДЕНИЕ.
ГЛАВА 1. ЛИТЕРАТУРНЫЙ ОБЗОР.
1.1 Результаты исследований пластической деформации, сопровождающейся излучением механических волн.
1.1.1 Понятие акустической эмиссии и ее применение.
1.1.2 История развития представлений о механизмах АЭ.
1.1.3 Направления исследований акустической эмиссии в настоящее время
1.2 Перераспределение массы деформируемых тел неоднородными сдвигами.
1.3 Методы компьютерного моделирования.
ГЛАВА 2. ИССЛЕДОВАНИЕ ПРОЦЕССА РЕЛАКСАЦИИ В КРИСТАЛЛАХ, СОДЕРЖАЩИХ ПАРЫ МЕЖУЗЕЛЬНЫХ АТОМОВ И ВАКАНСИЙ.
2.1 Выбор метода компьютерного моделирования и описание компьютерной модели.
2.2 Изучение процесса релаксации кристалла, содержащего две пары межузельных атомов, в зависимости от расстояния между ними.
2.3 Изучение процесса релаксации кристалла, содержащего две пары вакансий, в зависимости от расстояния между ними.
2.4 Изучение процесса релаксации кристалла, содержащего две пары межузельных атомов и подвергающегося деформации всестороннего растяжения.
2.5 Изучение процесса релаксации кристалла, содержащего две пары вакансий и подвергающегося деформации всестороннего сжатия.
2.6 Изучение процесса релаксации кристалла, содержащего две пары межузельных атомов или вакансий и подвергающегося нагреванию после релаксации.
ГЛАВА 3.ИССЛЕДОВАНИЕ ПРОЦЕССА РЕЛАКСАЦИИ В КРИСТАЛЛАХ, СОДЕРЖАЩИХ КОМПЛЕКСЫ ВАКАНСИЙ И МЕЖУЗЕЛЬНЫХ АТОМОВ.
3.1 Изучение процесса релаксации кристалла, содержащего комплексы межузельных атомов.
3.2 Изучение процесса релаксации кристалла, содержащего комплексы вакансий.
3.3 Изучение процесса релаксации кристалла, содержащего комплексы точечных дефектов, на 3-мерной модели.
ГЛАВА 4. МОДЕЛИ МАССОПЕРЕНОСА ПРИ СКОЛЬЖЕНИИ ДИСЛОКАЦИЙ В КОНТИНУАЛЬНОЙ СРЕДЕ.
4.1 Метод пробного блока для краевой дислокации.
4.2 Присвоенная масса покоя и масса движения дислокации.
4.2.1 Масса покоя пластических сдвигов, трещин и метод ее расчета
4.2.2 Расчет характеристик кондуктивного переноса массы.
4.3 Разделение полей на поле с дилатацией и без дилатации. Разделение поля краевой дислокации.
ГЛАВА 5.МОДЕЛИ МАССОПЕРЕНОСА ПРИ СКОЛЬЖЕНИИ ДИСЛОКАЦИЙ В КОНТИНУАЛЬНОЙ СРЕДЕ.
5.1 Построение модели задачи.
5.2 Расчет работы, затраченной на генерацию звуковых волн.
Актуальность проблемы. Теория дефектов кристаллической решетки объясняет механические свойства твердых тел, составляя основу физических и технических теорий пластичности, прочности, разрушения. К простейшим точечным дефектам кристаллической решетки относятся вакансии и межузельные атомы.
При образовании комплексов вакансий и межузельных атомов в материале, структура кристаллической решетки вблизи них искажается за счет смещений атомов в упругих полях дефектов, что вызывает локальное нарушение плотности. Области, окружающие межузельный атом или комплексы межузельных атомов, обладают повышенной локальной плотностью; вакансию или комплексы вакансий - пониженной локальной плотностью.
В процессе структурной релаксации внутри кристалла происходит перераспределение локальной плотности, завершающееся образованием областей «идеальной» плотности и новых дефектов кристаллической решетки. Подобные процессы, связанные с перемещением вещества внутри некой системы, принято называть массопереносом. Локальная перестройка внутренней структуры материала сопровождается излучением материалом механических волн, это явление получило название акустической эмиссии. Несмотря на значительный объем теоретических и экспериментальных исследований, связанных с явлением акустической эмиссии, вопрос о природе акустической эмиссии, сопровождающей пластическую деформацию, остается открытым. Основная сложность при экспериментальных исследованиях заключается в том, что релаксация кристаллической решетки происходит с высокой скоростью. Поэтому в данном случае актуальным является использование метода компьютерного моделирования.
Компьютерное моделирование, являющееся в настоящее время таким же признанным методом исследования как экспериментальный и теоретический метод, начало применяться в физике твердого тела с конца 5Ох годов 20 века. С его помощью на атомном уровне возможно исследование не только быстропротекающих процессов, как например, аннигиляция точечных дефектов, но и процессов более длительных по времени.
В первой части работы исследован процесс бездиффузионного массопереноса в кристалле алюминия, содержащем комплексы межузельных атомов и вакансий, с помощью компьютерного моделирования по методу молекулярной динамики. Показано, что процесс сопровождается генерацией упругих волн. Во второй части работы решена задача об излучении упругих волн при аннигиляции двух краевых дислокаций, исходя из предположения, что движущие дислокации перемещают массу.
Цель настоящей работы заключается в исследовании методом компьютерного моделирования релаксационных процессов и механизмов массопереноса в кристалле алюминия, содержащем такие носители свободного объема, как межузельные атомы, вакансии и их комплексы; построении аналитического решения для пластического сдвига в бесконечной среде с использованием представлений о сдвиговом переносе массы.
Для достижения поставленной в работе цели были сформулированы следующие задачи:
1. Построить дву- и трехмерную молекулярно-динамическую модель для исследования на атомном уровне релаксационных процессов в кристалле алюминия, содержащем точечные дефекты.
2. Изучить картины атомных смещений вблизи комплексов точечных дефектов на начальной стадии процесса релаксации и итоговые конфигурации кристаллической решетки после релаксации.
3. Исследовать процессы релаксации кристалла алюминия при различных стартовых расстояниях между парами внедренных точечных дефектов, а также влияние деформации сжатия/растяжения расчетного блока, нагревания расчетного блока и влияние стартовых конфигураций комплексов точечных дефектов на процесс релаксации.
4. Выявить механизмы массопереноса, оценить скорости данного процесса, энергетические характеристики.
5. Построить аналитическое решение для пластического сдвига в бесконечной среде, пользуясь представлениями о переносе массы.
6. Решить задачу о разделении упругого поля краевой дислокации на бездилатационную и дилатационную составляющие.
Научная новизна работы заключается в том, что:
1) при помощи метода молекулярной динамики впервые обнаружена стадийность процесса релаксации кристаллической структуры: на начальной стадии зафиксированы максимальные смещения атомов величиной в 0,8 параметра решетки; на последующей стадии наблюдалось увеличения числа направлений кооперативных смещений атомов;
2) показано, что при релаксации атомной структуры кристалла происходит перераспределение локальной избыточной (недостающей) плотности, что приводит к формированию новых комплексов точечных дефектов;
3) замечено, что этот процесс сопровождается генерированием ударной волны и ее последующей трансформацией в продольную волну; показано, что время массопереноса совпадает со временем распространения ударной волны, а трансформация ударной волны в продольную представляет собой эффект акустической эмиссии;
4) наделение структурного дефекта массовой характеристикой позволяет идентифицировать его как механический объект, который обладает импульсом и энергией, зависящими от стартовых условий, таких как радиус обрезания, расстояние до свободной поверхности, расстояние между дефектами;
5) установлено, что внешний импульс трансформируется в упругие моды разной частоты, по которым распределяется его энергия.
Научная и практическая ценность диссертационной работы заключается в возможности использования развитых представлений в радиационном материаловедении с целью создания новых материалов с заданными свойствами, а также с усовершенствованием свойств уже известных материалов подвергающихся различным экстремальным воздействиям в процессе эксплуатации. Кроме этого, результаты компьютерного моделирования могут быть использованы как демонстрационный материал, отображающий процессы, протекающие в кристаллических структурах, полезный для студентов, осваивающих курс физики.
На защиту выносятся следующие положения: I
1. Процесс релаксации кристаллической структуры, содержащей такие носители свободного объема, как межузельные атомы, вакансии и их комплексы, состоит из нескольких стадий: образование ударной волны, ее последующая трансформация в продольную волну.
2. Зависимость интенсивности процесса релаксации от типа комплексов дефектов.
3. Структурная релаксация при наличии таких носителей свободного объема, как межузельные атомы, вакансии и их комплексы в кристалле сопровождается выделением энергии в виде акустической эмиссии.
4. Зависимость параметров упругих волн от вида переносы массы: конвективного или кондуктивного.
Апробация работы. Результаты работы были доложены и обсуждены на следующих научных конференциях:
• 11 Всероссийская конференция студентов-физиков и молодых учёных, Екатеринбург, 2005. • 12 Всероссийская конференция студентов-физиков и молодых ученых, Новосибирск, 2006.
• Всероссийская научно-практическая конференция «Наука в ВУЗе», I
Новокузнецк, 2007.
• XVII Петербургские чтения по проблемам прочности, Санкт-Петербург, 2007.
• 13 Всероссийская конференция студентов-физиков и молодых ученых, Ростов-на-Дону, 2007.
• IV Международная конференция «Микромеханизмы пластичности, разрушения и сопутствующих явлений» (МРБР - 2007), Тамбов, 2007.
• Международная научная школа-конференция «Фундаментальное и прикладное материаловедение», Барнаул, 2007.
• 14 Всероссийская конференции студентов-физиков и молодых ученых, Уфа, 2008.
• X Международная школа-семинар "Эволюция дефектных структур в конденсированных средах" (ЭДС - 2008), Бийск, 2008.
• V Международная научно-техническая школа-конференция «Молодые ученые - науке, технологиям и профессиональному образованию в электронике», Москва, 2008.
• Вторая Всероссийская конференция с международным интернет -участием «От наноструктур, наноматериалов и нанотехнологий к Наноиндустрии», Ижевск, 2009.
• XV Международная научно-практическая конференция студентов, аспирантов и молодых ученых, Томск, 2009.
• Международный симпозиум «Перспективные материалы и технологии», Витебск, 2009. 1
• I региональная научно-практическая конференция «Наноиндустрия Алтая 2009», Бийск, 2009.
• XVII Международная конференции «Физика прочности и пластичности материалов», Самара, 2009.
• Всероссийская конференция с элементами научной школы для молодежи «Новые материалы. Создание, структура, свойства - 2009», Томск, 2009.
• Научно-техническая конференция с международным участием «Ультрадисперсные порошки, наноструктуры, материалы: получение, свойства, применение. V Ставеровские чтения», Красноярск, 2009.
• International Conference «Fundamental Aspects of External Fields Action on Materials», Novokuznetsk, 2010.
• V (XXXVII) Международная научно-практическая конференция «Образование, наука, инновации - вклад молодых исследователей», Кемерово, 2010.
• Всероссийская конференция студентов, аспирантов и молодых ученых, Владивосток, 2010.
• XI Международная школа-семинар "Эволюция дефектных структур в конденсированных средах" (ЭДС - 2010), Бийск, 2010.
Публикации. Результаты работы отражены в 30 публикациях, пять из которых в журналах, включенных в список ВАК для публикации диссертационных работ.
Структура и объем работы. Диссертация состоит из введения, пяти глав, заключения и списка литературы из 131 наименования. Работа изложена на 145 страницах машинописного текста, содержит 1 таблицу и 86 рисунков.
ЗАКЛЮЧЕНИЕ
В диссертационной работе с помощью метода молекулярной динамики проведены исследования основных механизмов бездиффузионного массопереноса при релаксации кристалла алюминия, содержащего локальные области свободного объема. Эксперименты проведены на 2-мерной и 3-мерной компьютерных моделях. Рассмотрено протекание процесса релаксации кристалла при импульсном внедрении в него комплексов межузельных атомов и вакансий. Исследовано влияние деформации кристалла, температурного воздействия на блок на итоговую структуру кристалла. В континуальном приближении решена задача о разделении поля смещений на две составляющие: бездилатационную и дилатационную, показано, что величина смещений в поле с дилатацией на порядок больше, чем для поля без дилатации, и их вклад в массоперенос различен. Для аннигиляции двух краевых дислокаций разработан метод расчета энергии, которая трансформируется в упругие моды разной частоты.
В работе сделаны следующие выводы:
1. Процесс релаксации кристаллической структуры, содержащей такие носители свободного объема, как межузельные атомы, вакансии и их комплексы, протекает с формированием ударной волны и ее последующей трансформацией в продольную.
2. Стадия перехода от ударных смещений атомов к возникновению продольных волн увеличивается по времени с ростом числа межузельных атомов в комплексе, т.е. с ростом избыточной локальной массы со знаком «+».
3. Для каждого из комплексов стадия перехода от ударных смещений атомов к возникновению продольных волн уменьшается по времени с ростом числа вакансий в комплексе, т.е. с ростом локальной массы со знаком «-».
4. Время массопереноса совпадает со временем распространения ударной волны, а трансформация ударной волны в продольную представляет собой эффект акустической эмиссии.
5. Даже при низкоэнергетическом воздействии на материал, когда не образуются точечные дефекты, образующаяся при этом продольная волна способна оказывать влияние как на процесс формирования комплекса дефектов, так и на уже сформировавшийся комплекс.
6. Процесс релаксации кристалла завершается формированием устойчивых структур атомов блока, свидетельствующих о завершении массопереноса.
7. Разработан метод расчета параметров акустических волн, генерируемых при аннигиляции краевых дислокаций.
8. Предложенный способ разделения упругого поля на дилатационную и бездилатационную часть дает возможность выделить в процессе деформации чисто конвективный и кондуктивный перенос.
1. ГОСТ. 27655-81. Акустическая эмиссия. Термины, определения и обозначения. М.: Издательство стандартов, 1988.
2. Гусев О. В. Акустическая эмиссия при деформировании монокристаллов тугоплавких металлов. М.: Наука, 1982. — 106 с. 1
3. Бовенко В.Н. Закономерности акустической эмиссии при деформировании кристалла. Изв.АН СССР. Металлы.1984.№1.С. 129-137.
4. Грешников В.А., Дробот, Ю.Б. Акустическая эмиссия / В.А. Грешников, Ю.Б. Дробот. М.: Издательство стандартов, 1976. - 243 с.
5. Иванов В.И., Белов В.М. Акустико-эмиссионный контроль сварки и сварных соединений. — М.: Машиностроение, 1981. 184 с.
6. Судхзуки Т., Есинага X., Такуети С. Динамика дислокаций и пластичность. / Т. Судхзуки, Х.Есинага, С.Такуети. М.: Мир, 1989. - 296 с.
7. Schofield В.Н. Research on the sources and characteristics of acoustic emission.1. : Acoustic emission, ASTM STP 505, Baltimore, 1972, p. 11-19.
8. Минц Р.И., Кортов B.C., Мелехин В.П. Влияние механизмов пластической деформации на акустическую и экзоэлектронную эмиссию. — В кн.: Металлофизика. Киев, 1973, вып. 44, с. 93-95.
9. Бойко B.C., Гарбер Р.И., Кривенко Л.Ф., Кривуля С.С. Звуковое излучение двойникующих дислокаций при их выходе из кристалла /B.C. Бойко, Р.И. Гарбер, Л.Ф Кривенко, С.С. Кривуля // ФТТ. 1969. - т. 12. - С. 3624-3626.
10. Бойко B.C., Гарбер Р.И., Кривенко Л.Ф., Кривуля С.С. Звуковое излучение двойникующих дислокаций / B.C. Бойко, Р.И. Гарбер, Л.Ф Кривенко, С.С. Кривуля // ФТТ. 1970. - т. 12, №6. - С. 1753-1755.
11. Keiser I. Erkenntnisse und Folgerungen aus der Messung von Gerauschen bei Zugbean-spruchung von metallischen Werkstoffen. — Arch. Eisenhuttenwesen, 1953, Bd. 24, H. 1/2, S. 43-45.
12. Gillis P.P. Dislocation mechanisms as possible sources of acoustic emission.
13. MTRSA, 1971,vol. 11, N3,p. 11-13.
14. Мерсон Д. JI. Физическая природа акустической эмиссии при деформационных процессах в металлах и сплавах: автореферат диссертации на соискание ученой степени д-ра физ.-мат. наук. Барнаул, 2001. 39 с.
15. Нацик В'.Д. Излучение звука дислокацией, выходящей на поверхность кристалла / В.Д. Нацик // Письма в ЖЭТФ. 1968. - №8. - С.324-328.
16. Нацик В.Д., Бурканов А.Н. Излучение рэлеевских волн краевой дислокацией, выходящей на поверхность кристалла / В.Д. Нацик, А.Н. Бурканов // ФТТ. 1972. - т. 14, №5. - С. 1289-1296.
17. Нацик В.Д., Чишко К.А. Звуковое излучение при аннигиляции дислокаций / В.Д. Нацик, К.А. Чишко // ФТТ. 1972. - т. 14, №11. - С. 31263132.
18. Косевич A.M., Маргвелашвили И.Г. Излучение электромагнитных и звуковых волн дислокацией, равномерно движущейся в ионном кристалле. — Изв. АН СССР, Сер. физ, 1967, т. 31, № 5, с. 848-850.
19. Eshelby J.D. Dislocations as a cause of mechanical damping in metals. — Proc. Roy. Soc. L., 1949, vol. A197, N 1050, p. 396-416.
20. Laub T., Eshelby J.D. The velocity of a wave along a dislocation. — Phil. Mag., 1966, vol. 14, N 132, p. 1285-1293.
21. Чишко K.A. Дислокационные механизмы акустической эмиссии пластически деформируемых кристаллов // Автореф. на соискание ученойстепени канд. ф.-м. наук. Харьков: ФТИНТ АН УССР, 1977. 22 с.
22. Нацик В.Д., Чишко К.А. Динамика и звуковое излучение дислокационного источника Франка—Рида. I. Начальная стадия работы источника. — В кн.: физика конденсированного состояния. Харьков: ФТИНТ АН УССР, 1974, вып. 33, с. 44-56.
23. Нацик В. Д., Чишко К.А. Динамика и звуковое излучение дислокационного источника Франка—Рида. II. Формирование дислокационного скопления: Препринт ФТИНТ АН УССР. Харьков, 1976,I26 с.
24. Косевич A.M. Динамическая теория дислокаций. — УФН, 1964, т. 84, № 4, с. 579-609.
25. Бойко B.C., Гарбер Р.И., Кривенко Л.Ф., Кривуля С.С. Переходное излучение звука дислокациями / B.C. Бойко, Р.И. Гарбер, Л.Ф Кривенко, С.С. Кривуля // ФТТ. 1973. -Т.15 , №1. - С.321-323.
26. Косевич A.M. Дислокации в теории упругости / A.M. Косевич Киев: Наук, думка, 1978. - 234 с.
27. Косевич A.M., Нацик В.Д. Уравнение движения дислокации в анизотропных средах / A.M. Косевич, В.Д. Нацик // ФТТ. 1965. - т. 7, №1. -С. 33-41.
28. Скворцов A.A., Орлов A.M., Соловьев A.A. Акустоэмиссионное зондирование линейных дефектов в кремнии / A.A. Скворцов, A.M. Орлов, A.A. Соловьев // ФТТ. 2001. - т.43, №4. - С. 616-618.
29. Скворцов A.A., Литвиненко О.В. Звуковое излучение, вызванное срывом и остановкой краевых дислокаций в изотропной среде / A.A. Соловьев // ФТТ. 2002. - т. 44, №7. - С. 1236-1242.
30. Чишко К.А. Переходное излучение звука винтовой дислокацией, выходящей на поверхность изотропной пластины / К.А. Чишко // ФТТ. -1989. т. 31, №1. - С. 223-229.
31. Чишко К.А. Излучение рэлеевских волн при выходе краевой дислокации на поверхность пластины / К.А. Чишко // ФТТ. 1992. - т. 34, №7. - С. 22362243.
32. Орлов A.M., Скворцов A.A., Фролов В.А. Акустическая эмиссия в дислокационном кремнии при токовых и тепловых воздействиях / A.M. Орлов, A.A. Скворцов, В.А. Фролов // Письма в ЖЭТФ. 1999. - т.25, №3. -С.28-32.
33. Орлов A.M., Скворцов A.A., Фролов В.А. Изменение спектра акустической эмиссии дислокационного кремния при токовых и тепловых воздействиях / A.M. Орлов, A.A. Скворцов, В.А. Фролов // Письма в ЖЭТФ. 1999. - т.25, №12. - С. 52-58.
34. Орлов A.M., Скворцов A.A., Насибов A.C. и др. Акустическая эмиссия в сульфиде кадмия при токовых и тепловых воздействиях / А.М.Орлов, A.A. Скворцов, A.C. Насибов, О.В. Литвиненко // Письма в ЖЭТФ. 2000. - т.26, №22. - С. 36-43.
35. Орлов A.M., Скворцов A.A., Фролов В.А. Магнитостимулированное изменение подвижности дислокаций в пластически деформированном кремнии п-типа / A.M. Орлов, A.A. Скворцов, В.А. Фролов // ФТТ. 2001. -т.43, №7. - С. 1207-1210.
36. Бойко B.C. Дислокационное описание динамического поведения полосы скольжения / B.C. Бойко // ФТТ. 1973. - т. 15, №11. - С. 3426-3427.
37. Бойко B.C., Гарбер, Р.И., Кривенко, Л.Ф. Звуковая эмиссия при аннигиляции дислокационного скопления / B.C. Бойко, Р.И. Гарбер, Л.Ф Кривенко // ФТТ. 1974. - т. 19, №.4. - С. 1233-1255.
38. Бойко B.C., Гарбер Р.И., Кившик В.Ф. Кривенко Л.Ф. Синхронная регистрация перемещения дислокаций и генерируемого ими звукового излучения / B.C. Бойко, Р.И. Гарбер, В.Ф. Кившик, Л.Ф Кривенко // ФТТ. — 1975. т. 17, №.5. - С.1541-1543.
39. Бойко B.C., Гарбер Р.И., Кившик В.Ф. Кривенко Л.Ф. Экспериментальное исследование переходного излучения звука дислокациями при их выходе на поверхность / B.C. Бойко, Р.И. Гарбер, В.Ф. Кившик, Л.Ф Кривенко // ЖЭТФ. 1976. - т.71, № 2(8). - С.708-713.
40. Бойко B.C., Кившик В.Ф. Кривенко Л.Ф. Экспериментальное исследование звукового излучения при аннигиляции дислокаций, в кристалле / B.C. Бойко, В.Ф. Кившик, Л.Ф Кривенко // ЖЭТФ. 1980. - т. 78 , №2. - С. 797-801.
41. Бойко B.C., Кривенко Л.Ф: Исследование пространственного распределения звукового излучения при пересечении поверхности скоплением дислокаций / B.C. Бойко, Л.Ф. Кривенко // ЖЭТФ. 1981. - т.80, №1. - С. 255-261.
42. Бойко B.C., Кившик В.Ф. Кривенко Л.Ф. Условия регистрации импульсов акустической эмиссии, генерируемых при выходе на поверхность отдельных дислокаций / B.C. Бойко, В.Ф. Кившик, Л.Ф Кривенко // ЖЭТФ. 1982. - т. 82 , №2. - С.504-508.
43. Бойко B.C., Кившик В.Ф. Кривенко Л.Ф. Распределение дислокаций в движущемся скоплении и генерируемая им акустическая эмиссия / B.C. Бойко, В.Ф. Кившик, Л.Ф Кривенко // ФТТ. 1984. - т.26, №6. - С. 16251629.
44. Бойко B.C., Иванченко Л.Г., Кривенко Л.Ф. Регистрация сигналов акустической эмиссии при образовании отдельного полумикроскопического двойника в цирконии / B.C. Бойко, Л.Г. Иванченко, Л.Ф. Кривенко // ФТТ. -1984. т.26 , №7. - С. 2207-2209.
45. Бойко B.C., Кривенко Л.Ф. Наблюдение поверхностных акустических волн, генерируемых при выходе скопления дислокаций на поверхность кристалла / B.C. Бойко, Л.Ф Кривенко // Акустический журнал. 1987. - т. 23, №5. - С.821-825.
46. Иосилевский Я.А. / Я.А. Иосилевский // Кристаллография. 1971. - №16.- С.279-282.
47. Нацик В.Д., Чишко К.А. Звуковое излучение дислокаций, движущихся у поверхности кристалла / В.Д. Нацик, К.А. Чишко // ФТТ. 1978. - т. 20, №2.- С. 457-465.
48. Нацик В.Д., Чишко К.А. Акустическая эмиссия дислокаций, выходящих на поверхность кристалла / В.Д. Нацик, К.А. Чишко // Акустический журнал.- 1982. т. 28, №3. - С. 381-389.
49. Кузнецов Н.С. Теория и практика неразрушающего контроля изделий с помощью акустической эмиссии: Мет. пособие. М.: Машиностроение, 1998. -96 с.
50. Криштал М.А., Мерсон Д.Л., Алехин В.П., Зайцев В.А. Распространение пластической деформации по сечению образца и акустическая эмиссия приодноосном растяжении меди // Физика металлов и металловедение. 1987. т.бЗ.вын.5. С. 1011-1016.
51. Березин A.B., Козинкина А.И., Рыбакова JI.M. Акустическая эмиссия и деструкция пластически деформированного кристалла // Дефектоскопия. — 2004.-№3.-С. 9-14.I
52. Бовенко В.Н. Связь акустической эмиссии с предразрушающим состоянием металла. Изв. АН СССР. 1983. т.271. №5. С. 1086-1090.
53. Бовенко В.Н. ФХММ, 1980, № 6, с. 13-16.
54. Осташев В.Е. Распространение звука в движущихся средах. Гл.ред.физ.-мат.лит., 1992.-230 с.
55. Бовенко В.Н. Закономерности акустической эмиссии при деформировании кристалла. Изв.АН СССР. Металлы. 1984.№1.С. 129-137.
56. Бовенко В.Н. Синергетические эффекты при пластической деформации иразрушении кристалла. Изв. АН, серия физическая. 1986. т.50. №3. С.509-513.
57. Бовенко В.Н. Автоколебания акустической эмиссии при деформации и разрушении твердых тел. Сб. Акустическая эмиссия материалов и конструкций. Ростов-на-Дону. 1984.
58. Бибик З.И., Нацик В.Д. Акустическая эмиссия при пластической деформации поликристаллов алюминия высокой чистоты // Металлофизика. 1982. - т.4. - С. 92-99.
59. Криштал М.А., Кацман A.B., Выбойщик М.А. Взаимодействие дислокаций с межузельными атомами и вакансиями // Физика металлов и металловедение. 1984. т. 57. вып. 2. С.374-379.
60. Криштал М.А., Выбойщик М.А., Мерсон Д.Л., Кацман A.B. Влияние микролегирования на модуль упругости меди // Физика металлов и металловедение. 1986. т. 62. Вып. 4. С: 813-819.
61. Криштал М.А. Взаимодействие дислокаций с примесными атомами и свойства металлов // Физика и химия обработки металлов. 1975. № 1. С. 62! 71.
62. Криштал М.А., Троицкий И.В. Равновесное распределение примесныхIатомов вокруг дислокации / В кн. Взаимодействие между дислокациями и атомами примесей в металлах и сплавах. Тула, 1969. С. 218-304.
63. Мерсон Д.Л. , Бергардт A.A., Шендерей П.Э. Исследование неоднородности пластической деформации при растяжении стали // Тезисы докладов XII Всесоюзной конференции по физике прочности и пластичности металлов и сплавов, 1989. Куйбышев.
64. Мерсон Д.Л., Богапов М.А. Влияние площади и состояния поверхности на акустическую эмиссию при деформировании меди // Тезисы докладов XIV Международной конференции «Физика прочности и пластичности материалов», Самара 1995. С. 368-369.
65. Дробот Ю.Б., Корчевский В.А. Исследование связи акустической эмиссии с образованием полос скольжения . при пластическом • деформировании аустенитной стали // Дефектоскопия, 1985. № 6. - С. 38-42.
66. Плотников В.А., Макаров C.B. Акустическая эмиссия при высокотемпературной деформации меди // Вестник ТГУ. Бюллетень оперативной научной информации. 2005. - №44. - С. 58-64.
67. Плотников В.А. Акустическая эмиссия при нагреве деформированного алюминия / В.А. Плотников // Письма в ЖЭТФ. 2001. - т. 27, №15. - С. 2732.
68. Плотников В.А. Акустическая эмиссия при нагреве деформированного алюминия / В.А. Плотников // Письма в ЖЭТФ. 2001. - т. 27, №15. - С. 2732.
69. Плотников В.А. Моделирование взрывной акустической эмиссии при мартенситных превращениях в сплавах / В.А. Плотников // Письма в ЖЭТФ.- 1998. т. 24, №1. - С. 31-38.
70. Плотников В.А. Акустическая эмиссия и динамическая релаксация нехимической энергии при мартенситных превращениях / В.А. Плотников // Письма в ЖЭТФ. 1999. - т. 25, №13. - С. 16-22.
71. Плотников В.А. Инверсия асимметрии акустической эмиссии при мартенситных превращениях в сплавах на основе никелида титана / В.А. Плотников // Письма в ЖЭТФ. 1999. - т. 25, №24. - С. 41-49.
72. Плотников В.А., Макаров С.В. Акустическая эмиссия при высокотемпературной деформации алюминия // Деформация и разрушение материалов. 2005. - №3. - С. 27-31.
73. Потекаев А.И., Плотников В.А. Акустическая диссипация энергии при термоупругих мартенситных превращениях / А.И. Потекаев, В.А. Плотников.- Томск: Изд-во НТЛ, 2004. 196 с.
74. Плотников В.А., Пачин И.М., Акустическая эмиссия при изотермической деформации никелида титана // Тезисы докладов XVII Международной конференции «Петербургские чтения по проблемам прочности». Санкт-Петербург 2007. - С. 275-277.
75. Неверов В.В., Суппес В.Г. Механизм замуровывания поверхностных пленок под действием больших нормальных нагрузок и сдвиговой деформации. // Трение и износ, 1991, т. 12, №3, С. 87 — 90.
76. Неверов В.В., Житников П.П. Поворотные движения материала при сдвиговой пластической деформации тонких слоев. // Изв. вузов СССР. Физика, 1989, № 2, С. 78 82.
77. Неверов В.В., Житников П.П., Горкавенко В.В. Образование областей высокого давления при пластической деформации. // Изв. вузов СССР, Физика, 1990, № 6, С. 85 90.
78. Неверов В.В., Житников П.П. Пластические движения и перемешивание в деформируемых смесях металлов. II Изв. ВУЗов. Физика, 1994, № 2, С. 10 — 15.
79. Неверов В.В. Диссипативная мезоскопическая структура пластической деформации. // ПМТФ, 1993, № 4, С. 128 139.
80. Неверов В.В. Деформационные процессы в тонких слоях при сжатии и сдвиге. Дисс. на соиск. ст. док. физ.-мат. наук. М., МИСИС. 1996. 347 с.
81. Неверов В.В., Антоненко А.И. Теория пластических сдвигов. Перенос массы. Скачки / В.В. Неверов, А.И. Антоненко. Новокузнецк: Изд-во КузГПА, 2005.- 194 с.
82. Неверов В.В., Молотков С.Г. Масса покоя сдвигов и трещин / В.В. Неверов, С.Г. Молотков // Известия ВУЗов. Физика. 2006. - №1. - С. 70-78.
83. Неверов В.В. Перенос массы пластическими сдвигами /В.В. Неверов // Деформация, локализация, разрушение: Сб. статей / Под ред. докт. физ-мат. наук, проф. Л.Б. Зуева. Томск: Изд-во НТЛ, 2005. - 148 с. (с. 57-68).
84. Баренблатт Г.И., Черепанов Г.П. О конечности напряжений на краю произвольной трещины // ПММ. 1961. - №4. - С. 752-753.
85. Партон В.З., Борисковский В.Г. Динамика хрупкого разрушения / В.З. Партон, В.Г. Борисковский. М.: Машиностроение, 1988. - 240 с.
86. Антоненко А.И. Модели скачкообразного развития сдвигов. Дисс. на соиск. ст. канд. физ.-мат. наук. Барнаул, 2004. 128 с.
87. Шишкин Ю.М. Методы машинного моделирования в теории дефектов кристаллов В кн.: Дефекты в кристаллах и их моделирование на ЭВМ.-Л.: Наука, 1980. - С. 77-99.
88. Лихачев В.А., Шудегов В.Е. Принципы организации аморфных ' структур. СПб.: Изд-во С.-Петербургского университета, 1999. — 228 с.
89. Хеерман Д.В. Методы компьютерного эксперимента в теоретической физике: Пер. с англ./ Под ред. С.А. Ахманова. М.: Наука, 1990. - 176 с.
90. Займан Дж. Модели беспорядка. Теоретическая физика однородно неупорядоченных систем: Пер. с англ. М.: Мир, 1982. - 592 с.
91. Андрухова О.В. Компьютерное моделирование атомного упорядочения и фазового перехода порядок-беспорядок в бинарных сплавах стехиометрического состава. Диссертация на соискание ученой степени кандидата физико-математических наук. Барнаул, 1997. - 225 с.
92. Гурова Н.М. Компьютерное моделирование термоактивируемых превращений, протекающих на антифазных и межфазных границах. Диссертация на соискание ученой степени кандидата физико-математических наук. Барнаул, 2000. - 171 с.
93. Nose S. A unified formulation of the constant temperature molecular dynamics methods //J. Chem. Phys. 1984. - V.81, № 1. - P. 511-519.
94. Лихачев В.А., Шудегов B.E. Принципы организации аморфных структур. СПб.: Изд-во С.-Петербургского университета, 1999. - 228 с.
95. Полухин В.А., Ватолин H.A. Моделирование аморфных металлов. М.: Наука, 1985.-288 с.
96. Полетаев Г.М., Краснов В.Ю., Старостенков М.Д. Исследование структуры аморфных металлов // Труды 9-й междунар. научн.-техн. конференции "Композиты в народное хозяйство" (Композит - 2005). -Барнаул: изд-во АлтГТУ, 2005. - С. 129-133.
97. Валуев A.A., Норманн Г.Э., ПодлипчукВ.Ю. Метод молекулярной динамики: теория и приложения / В кн.: Математическое моделирование: Физико-химические свойства вещества. М.: Наука, 1989. - С. 5-40.
98. ПацеваЮ.В. Исследование особенностей самодиффузии в двумерных металлах. Диссертация на соискание ученой степени кандидата физ.-мат. наук, Барнаул, 2005. 136 с.
99. Чирков А.Г., Понаморев А.Г., Чудинов В.Г. Динамические свойства Ni, Cu, Fe в конденсированном состоянии (метод молекулярной динамики) // ЖТФ. 2004. - Т.74, №2. - С. 62-65.
100. Старостенков Д.М., Старостенков М.Д., Демьянов Б.Ф., Полетаев Г.М. Самоорганизация дефектных структур в металлах при нагреве // Фундаментальные проблемы современного материаловедения. 2005. -Т.2, №3. - С. 93-97.
101. Gumbsch P., Zhou S.J. and Holian B.L. Molecular dynamics investigation of dynamic crack stability // The American Physical Society. 1997. - V.55, №6. -P. 3445-3455.
102. Полетаев Г.М. "Моделирование методом молекулярной динамики структурно-энергетических превращений в двумерных металлах и сплавах (MD2)".Свидетельство о гос. регистрации программы для ЭВМ № 2008610486 от 25.01.2008.
103. Царегородцев А.И., Горлов Н.В., Демьянов Б.Ф., Старостенков М.Д. Атомная структура АФГ и ее влияние, на состояние решетки вблизи дислокации в упорядоченных сплавах со сверхструктурой L12 //ФММ, 1984, т.58, вып.2. С.336-343.i
104. Маркидонов A.B. Бездиффузионный механизм массопереноса в кристаллах, содержащих агрегаты вакансий и межузельных атомов // Автореф. на соискание ученой степени канд. ф.-м. наук. Барнаул, 2009. 22 с.
105. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А. Нелинейный высокоскоростной массоперенос в двумерном кристалле при наличиилокальных областей с различной плотностью // Ползуновский альманах.2008. №З.С.226-228.
106. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А., Медведев H.H. Высокоскоростной массоперенос в двумерном кристалле никеля при наличии дислокационных петель различной локальной плотности // Изв. вузов. Черная металлургия. 2009. №6. С.57-60.
107. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А., Потекаев А.И., Кулагина В.В. Высокоскоростной массоперенос в кристаллическом алюминии, содержащем цепочки вакансий и межузельных атомов // Изв. вузов. Физика. 2009. т.52. №9/2. С. 139-145.
108. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А. Высокоскоростной массоперенос в кристалле при наличии различных конфигураций точечных дефектов // Фундаментальные проблемы современного материаловедения.2009. т.6. №1. С.12-16.
109. Исследование характеристик акустико-эмиссионных сигналов при разрушении пористого железа // Известия Алтайского университета. 2009. №1 (61).-С. 124-128.
110. Штремель М.А. Прочность сплавов. Часть I. Дефекты решетки: Учебник для вузов. М.: МИСИС, 1999. - 384 с.
111. Белоус В.А., Лапшин В.И., Марченко И.Г., Неклюдов ИМ. Радиационные технологии модификации поверхности. I. Ионная очистка ивысокодозовая имплантация // Физическая инженерия поверхности. 2003. т.1. № 1.С.40-48.
112. Полетаев Г.М. Атомные механизмы структурно-энергетических превращений в объеме кристаллов и вблизи границ зерен наклона в ГЦК металлах / Диссертация на соискание ученой степени доктора физ.-мат. наук. -Барнаул, 2008. 412 с.
113. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А., Медведев H.H. Высокоскоростной массоперенос в двумерном кристалле никеля при наличии дислокационных петель различной локальной плотности // Изв. вузов. Черная металлургия. 2009. №6. С.57-60.
114. Старостенков М.Д., Маркидонов A.B., Тихонова Т.А., Потекаев А.И., Кулагина В.В. Высокоскоростной массоперенос в кристаллическом алюминии, содержащем цепочки вакансий и межузельных атомов // Изв. вузов. Физика. 2009. т.52. №9/2. С.139-145.
115. Ракитин Р.Ю. Границы зерен в металлах (GB3D) / РОСПАТЕНТ. Свидетельство № 2009610714 от 30 января 2009.
116. Аксенов М.С., Полетаев Г.М., Ракитин Р.Ю., Краснов В.Ю., Старостенков М.Д. Стабильность вакансионных кластеров в ГЦК металлах // Фундаментальные проблемы современного материаловедения. 2005. т.2. №4. С.24-31.
117. Хирт Дж., Лоте И. Теория дислокаций / Дж. Хирт, И. Лоте. М.: Атомиздат, 1972. - 600 с.
118. Мусхелишвили Н.И. Некоторые основные задачи математической теории упругости / Н.И. Мусхелишвили. М.: Наука, 1966. — 708 с.
119. Неверов В.В., Антоненко А.И. Превращение энергии при скачкообразном развитии пластического сдвига / В.В. Неверов, А.И. Антоненко // Физическая мезомеханика. 2004. - т.7, №3. - С. 43-52.
120. Неверов В.В. Масса покоя сдвигов и трещин // Изв. Вузов. Физика. -2006. №1. - С.70-79.
121. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. В 10-ти томах. Т.7. Теория упругости / Л.Д.Ландау, Е.М.Лифшиц. М.: Главная редакция физ.-мат. литературы изд-ва «Наука», 1987. - 248 с.
122. Эшелби Дж. Континуальная теория дислокаций / Дж. Эшелби. М.: Изд-во иностр. лит., 1963. - 248 с.