Динамика дислокаций в щелочно-галоидных кристаллах при наложении ультразвука тема автореферата и диссертации по физике, 01.04.07 ВАК РФ
Дегтярёв, Вячеслав Тихонович
АВТОР
|
||||
доктора физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Москва
МЕСТО ЗАЩИТЫ
|
||||
2007
ГОД ЗАЩИТЫ
|
|
01.04.07
КОД ВАК РФ
|
||
|
На правах рукописи
Дегтярев Вячеслав Тихонович
ДИНАМИКА ДИСЛОКАЦИЙ В ЩЕЛОЧНО-ГАЛОИДНЫХ КРИСТАЛЛАХ ПРИ НАЛОЖЕНИИ УЛЬТРАЗВУКА
Специальность 01 04 07 -физика конденсированного состояния
АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора физико-математических наук
Москва - 2007
003061427
Работа выполнена в Калужском филиале государственного образовательного учреждения высшего профессионального образования «Московский государственный технический университет имени Н.Э. Баумана»
Научный консультант-
доктор физико-математических наук,
профессор
Тяпунина Наталья Александровна
Официальные оппоненты.
доктор физико-математических наук,
профессор
доктор физико-математических наук, профессор
доктор технических наук, профессор
Блантер Михаил Соломонович Старостенков Михаил Дмитриевич Абрамов Олег Владимирович
Ведущая организация:
НИИ математики и механики им. акад В И. Смирнова Санкт-Петербургского государственного университета
Защита состоится «¿¿>> e¿j/s»sJp&2007r в i5 час._мин. на заседании
диссертационного совета Д 217.035.01 при Федеральном государственном унитарном предприятии «Центральный научно-исследовательский институт черной металлургии им ИШБардина» по адресу: 105005, Россия, Москва, 2-я Бауманская ул, д.9/23
С диссертацией можно ознакомиться в библиотеке Федерального государственного унитарного предприятия «Центральный научно-исследовательский институт черной металлургии им. И ПБардина»
Автореферат разослан «9 » 2007 г.
Ученый секретарь диссертационного совета
доктор технических наук, . ^
старший научный сотрудник 9л '/¿^¿^¿¿^— Александрова Н.М.
ОБЩАЯ ХАРАСТЕРИСТИКА РАБОТЫ
Актуальность темы работы. Взаимодействие структурных дефектов и полей различной природы относится к фундаментальным проблемам современной физики. Темой данной работы является исследование взаимодействия ультразвукового поля и дефектов кристаллической структуры, установление связи между микропроцессами и макроскопическими свойствами кристаллов в ультразвуковом поле.
Поскольку под действием ультразвука материал может как упрочняться, так и разупрочняться, важно проследить за физическими эффектами, приводящими к таким изменениям Современное состояние теории не позволяет однозначно предсказать, что именно произойдет в процессе воздействия ультразвука упрочнение или разупрочнение кристалла Поэтому выяснение причин и механизмов, вызывающих изменения пластических свойств кристаллов под влиянием высокочастотной вибрации ультразвукового диапазона частот, остается фундаментальной задачей физики конденсированных сред и в тоже время имеет большое прикладное значение
Наиболее эффективным способом исследования взаимодействия полей и дефектов кристаллической структуры является метод компьютерного моделирования. Действительно, при исследовании дислокационных процессов, происходящих за времена порядка ~10'5с в объеме образца, другие методы практически не приемлемы Избирательное травление позволяет установить только начальное и конечное состояния системы дислокаций Просвечивающей электронной микроскопии доступны лишь образцы в виде тонких пленок Процессы в тонких пленках и массивных кристаллах могут существенно отличаться, поэтому закономерности, установленные для пленок, нельзя использовать для массивных кристаллов Кроме того, в реальных условиях воздействие бывает комплексным, и выделить влияние отдельных факторов не представляется возможным Моделирование позволяет выяснить роль отдельных факторов, приводящих к изменению макроскопических свойств материалов, и микромеханизмы, обуславливающие эти изменения Поэтому разработка моделей, алгоритмов и программ для осуществления моделирования дислокационных процессов также является актуальной задачей для современной физики твердого тела и важна для прикладных задач.
Целью настоящей работы явилось- детальное исследование процессов, происходящих при воздействии высокочастотных колебаний ультразвукового диапазона на дислокационные структуры и пластичность материалов; - определение режимов нагружения, позволяющих добиться необходимого изменения пластических свойств образца, будь то упрочнение или разупрочнение
Объекты исследования щелочно-галоидные кристаллы со структурой
хлорида натрия
Для достижения этих целей предстояло решить следующие задачи. 1. Разработать физическую модель и методику моделирования механизмов и процессов, обуславливающих акустопластический эффект, используя строгий динамический подход, основывающийся на решении уравнения движения с учетом поля сил взаимодействия дислокаций и самодействия, а также гибкости скользящих дислокаций и влияния полей, обусловленных внешним нагружением
2 Промоделировать процессы движения и размножения дислокаций в условиях сложнонагруженного состояния кристалла, когда на дислокации леса действует знакопеременная нагрузка, а в плоскости скользящей дислокации действуют или постоянная сила, или сумма постоянной и знакопеременной составляющих поля
3. На основании анализа полученных результатов моделирования выявить зависимости критического напряжения, необходимого для начала пластической деформации, от плотности лесных дислокаций и от параметров ультразвука, то есть установить связь между дислокационными процессами и акустопластическим эффектом
Методы исследования, реализованные для достижения цели работы, основаны на развитии существующих и создании новых, адекватных рассматриваемым физическим явлениям и процессам, математических моделей, связанных с исследованием дислокационных процессов в условиях высокочастотной вибрации
Для анализа напряжения, необходимого для преодоления пробной дислокацией модельной площадки, были разработаны квазистатическая и динамическая модели, с целью исследования особенностей перераспределения дислокаций леса в ультразвуковом поле использована динамическая модель
Достоверность полученных результатов обеспечена детальным рассмотрением физических явлений и процессов, связанных с воздействием ультразвука на кристаллы с дефектами, что позволило создать оригинальные математические модели, реализующие новые подходы к количественному и качественному описанию дислокационных процессов, в присутствии ультразвукового поля и в значительной степени тем, что результаты моделирования совпадают с результатами натурных экспериментов, выполненных другими авторами
Научная новизна данной диссертационной работы заключается в том, что впервые предпринята попытка методом компьютерного моделирования исследовать зарождение и движение дислокаций в ультразвуковом поле в условиях, приближенных к реальной ситуации Впервые предложена математическая модель, описывающая работу источника Франка-Рида, концы которого закреплены лесными дислокациями, совершающими гармонические колебания, и установлены закономерности поведения источника в этой ситуа-
ции
Впервые найдены условия, при которых, при наличии знакопеременной нагрузки, возможно накопление дислокаций в плоскости источника
Предложена модель, позволяющая описывать прохождение пробной дислокацией модельной площадки, в которой движение дислокации описывается в приближении динамической теории
Практическая ценность работы состоит в том, что полученные в ней результаты позволяют предвидеть изменение пластических свойств материалов, подвергающихся высокочастотной вибрации в процессе их эксплуатации, что важно для оценки срока службы и предвидения возможности деградации материалов Кроме того, в современных технологиях обработки материалов все более широкое применение находит ультразвук, причем как в процессе обработки, так и для придания материалам наперед заданных свойств, например, чтобы достичь высокой пластичности материала непосредственно во время технологического воздействия Основание для проведения работ.
План совместных НИР кафедры «Молекулярная физика» физического факультета МГУ имени М.В Ломоносова, кафедры «Компьютерные системы и сети» и кафедры «Системы автоматизированного проектирования» Калужского филиала МГТУ имени Н Э Баумана
Научные положения и научные результаты, выносимые на защиту
1. Установленные впервые закономерности генерации дислокаций в условиях сложнонагруженного состояния, когда точки закрепления источника дислокациями леса совершают вынужденные колебания
2 Обнаруженное влияние соотношения фаз колебаний закрепляющих дислокаций и источника Франка-Рида, в условиях наложения ультразвукового поля, и снижение эффекта упрочнения в случае колеблющегося дислокационного леса по сравнению с неподвижным лесом дислокаций
3 Эффект пластифицирования материала при воздействии на лес дислокаций ультразвукового поля и усиление этого эффекта с увеличением плотности дислокаций леса
4 Явление самоорганизации системы дислокаций леса в ультразвуковом поле, где самоорганизация включает в себя образование динамических дислокационных структур с диполями, мультиполями и стенками дислокаций в качестве компонентов самоорганизованной структуры, причем конечным результатом самоорганизации может стать формирование блочной структуры.
5 Обнаружение устойчивости границ блоков, образовавшихся под действием ультразвука, т е сохранения этих границ после выключения ультразвукового поля. При этом структура границ ультразвукового проис-
хождения отличается от структуры границ термического происхождения (образовавшихся в процессе отжига)
6 Разработанные модель, алгоритм и пакет программ для анализа поведения элементов дислокационной структуры в ультразвуковом поле
7 Результаты моделирования поведения элементов дислокационной структуры в условиях одновременного воздействия на кристалл постоянной и знакопеременной составляющих внешней нагрузки Личный вклад автора в проведенное исследование.
Автором проведен детальный анализ экспериментальных данных, известных из литературы к началу исследования, для щелочно-галоидных кристаллов, на основании которого предложены модели, адекватно описывающие поведение элементов дислокационной структуры в ультразвуковом поле Разработаны алгоритмы и пакет программ, позволяющих моделировать взаимодействие элементов дислокационных структур с ультразвуковым полем и эволюцию дислокационной структуры кристалла под воздействием ультразвука
Автором на кристаллах со структурой NaCl, для которых имеются надежные экспериментальные данные, выполнено компьютерное моделирование поведения элементов дислокационной структуры и ансамблей дислокаций в ультразвуковом поле Проведен анализ полученных компьютерным моделированием результатов и установлены закономерности эволюции дислокационных ансамблей в ультразвуковом поле
Апробация. Результаты диссертационной работы докладывались на 6 Всероссийских и Международных научных конференциях и научно-практических семинарах, в т ч на X Международной конференции «Imperfection interaction and anelasticity phenomena in solids (IIAPS - 10)» (Тула, 2001), Всероссийской конференции, посвященной 100-летию со дня рождения академика Г.В Курдюмова «Дефекты структуры и прочность кристаллов» (Черноголовка, 2002), XLII Международной конференции «Актуальные проблемы прочности» (Калуга, 2004), XLIV Международной конференции «Актуальные проблемы прочности» (Вологда, 2005), третьей Международной конференции «Кристаллофизика 21-го века» (Черноголовка, 2006)
Публикация результатов работы. Материалы исследований, представленных в диссертации, изложены в 30 работах, опубликованных в научных журналах, научно-технических сборниках, материалах, трудах и тезисах докладов Всесоюзных, Российских и Международных конференций, симпозиумов, совещаний и семинаров; ссылки на 30 из них приведены в тексте диссертации
Все основные результаты диссертационной работы опубликованы в 18 статьях, из них 78% (14 статей) опубликовано в ведущих российских журналах, в которых должны быть опубликованы основные научные результаты диссертации на соискание ученой степени доктора наук
Структура и объем работы.
В соответствии с поставленными целями исследования, характером и объемом проведенной работы, диссертация содержит введение, 5 глав, заключение, список литературы Общий объем диссертации составляет 218 страниц машинописного текста, содержащих текст работы, 78 рисунков, 5 таблиц и список использованных источников на 17 страницах, содержащий 142 наименования
КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обосновывается актуальность темы диссертации, ее практическая значимость, определяются цели исследования, приводятся основные положения, выносимые на защиту, изложена структура диссертации
Первая глава посвящена обзору литературы и анализу опубликованных до начала исследования экспериментальных данных, характеризующих поведение материалов в ультразвуковом поле в условиях комплексного на-гружения
Глава вторая содержит описание разработанных и использованных в данной работе моделей и алгоритмов моделирования дислокационных процессов В соответствии с поставленной в работе задачей необходимо было рассмотреть влияние ультразвукового поля на предел текучести. В процессе моделирования были использованы следующие предположения
1 Первоначальные координаты точек пересечения дислокаций леса с плоскостью скольжения пробной дислокации задавались с помощью стандартных программ псевдослучайных чисел
2 Предполагалось, что дислокации леса совершают вынужденные периодические колебания, оставаясь прямолинейными и сохраняя первоначальную ориентацию к линии пересечения плоскости своего скольжения с плоскостью скольжения пробной дислокации
3 При моделировании прохождения пробной дислокации через модельную площадку направление внешней периодической нагрузки, раскачивающей дислокации леса, выбиралось таким образом, чтобы фактор Шмида в системе скольжения пробной дислокации был равен нулю, то есть пробная дислокация оказывалась неподверженной непосредственному воздействию внешней периодической нагрузки
4 Двигающаяся через лес пробная дислокация считалась гибкой, ее форма в зависимости от модели определялась либо с применением квазистатического подхода, либо в процессе решения дифференциального уравнения
5 Самодействие пробной дислокации учитывалось в приближении линейного натяжения, когда каждому элементу дислокационной линии приписывалась погонная энергия, равная где Я - радиус кривизны линии дислокации
6 Считалось, что пробная дислокации открепляется от препятствий в том случае, когда угол между соседними сегментами дислокации становится меньше 10° или когда дислокация огибает группу дислокаций леса так, что две ее различные дуги касаются друг друга Это предположение является прямым следствием континуальной аппроксимации, принятой в современной теории дислокаций, и позволяет при моделировании не принимать во внимание образование порогов
7 Допускалось, что образование зон рекомбинации при пересечении пробной дислокации с дислокациями леса не происходит.
8 Для исключения края модельной площадки учитывалось также действие дислокаций, находящихся вне основной области При этом использовались периодические граничные условия, то есть расположение дислокаций на соседних к основной площадках моделирования задавалось таким же, как на основной модельной площади
9 При моделировании дислокаций леса их движение считалось консервативным
За величину предела текучести при постоянной нагрузке (режим ползучести) принимают напряжение, необходимое для выхода дислокации из объема образца на поверхность В режиме активного нагружения суммарный пластический сдвиг значительно превосходит величину сдвига, который может быть обеспечен дислокациями, содержащимися в исходном состоянии образца В процессе пластической деформации происходит размножение дислокаций Поэтому за предел текучести в режиме активного нагружения принимают напряжение, необходимое для начала размножения дислокаций
В соответствии с вышесказанным при моделировании акустопластиче-ского эффекта были рассмотрены следующие модели
• модель, описывающую влияние ультразвука на процесс прохождения дислокацией модельной площадки,
• модель, в которой учитывается влияние ультразвука на генерацию дислокаций источником Франка-Рида,
• модель, описывающую перераспределение дислокаций в ультразвуковом поле
И в первой, и во второй моделях основным элементом является дислокационный сегмент Скользящая дислокация разбивается лесными дислокациями, служащими центрами сильного закрепления, на совокупность сегментов, а единичный сегмент может служить источником Франка-Рида
Схема, иллюстрирующая использованную модель, представлена на рис 1, где а - плоскость легкого скольжения пробной краевой дислокации Лесные дислокации располагаются в плоскостях, перпендикулярных а, и колеблются под действием ультразвука Направления смещений лесных дислокаций относительно начального положения изображены на рис 1 стрелками Причем в первый полупериод действия ультра-
звука смещение происходит в направлении сплошной стрелки, во второй полупериод - в направлении стрелки, изображенной пунктиром
Причиной, вызывающей движение дислокаций, является нагрузка,
приложенная к образцу Рассматривается условие слож-нонагруженного состояния образца, такое, что на дислокацию, движущуюся в плоскости легкого скольжения, действует постоянная сила, а на дислокации леса действует сила, изменяющаяся во времени по гармоническому закону.
Для определения закономерностей движения дислокаций используются как квазистатическое, так и динамическое приближение теории дислокаций
В квазистатическом приближении в каждой точке скользящей дислокации (хпу,) в каждый момент времени должно выполняться условие равновесия
0,50Ь2ф(х„у„() = Ь[т°"+ Х*/*,,*,*,.^»')], (1)
где ф - локальная кривизна скользящей дислокации, гу- напряжение, обусловленное ./ -ой дислокацией леса, й - область суммирования, удовлетворяющая критерию Предводителева А А - Стратана И В
Моделирование в квазистатическом приближении проводилось по следующему алгоритму
На первом этапе, при фиксированном уровне внешнего напряжения сдвига г'" находилось семейство равновесных конфигураций скользящей дислокации между двумя соседними препятствиями вдоль скользящей дислокации для моментов времени, отвечающих различным значениям фазы колеблющихся дислокаций леса <р из интервала значений Построение
равновесной конфигурации скользящей дислокации при фиксированных значениях т" и (р проводилось при строгом учете дапьнодействующих полей
Рис 1 Схематическое изображение рассматриваемой модели
внутренних напряжений, создаваемых дислокациями леса с учётом ориентирующего действия, которое дислокации леса оказывают на скользящую дислокацию в окрестностях особых точек их пересечения Нахождение предельного г*" для рассматриваемого участка скользящей дислокации между двумя соседними препятствиями вдоль скользящей дислокации проводилось методом постепенного увеличения уровня внешнего напряжения сдвига путем перехода от (г™,), к (ОД+Аг"" вплоть до предельного значения (г"„),, при котором равновесные конфигурации рассматриваемого / -того участка между двумя препятствиями вдоль скользящей дислокации оказываются все еще возможными для моментов времени, отвечающих различным значениям фазы колеблющихся дислокаций леса из интервала значений Найденные семейства равновесных конфигураций скользящей дислокации, отвечающие значениям т'н из интервала [(г™,),,(г"и),]> помещались в архив решений
На втором этапе, на основании полученных решений уравнения (1) для различных 1-ых участков вдоль скользящей дислокации формировался массив значений напряжений (г™„),, отвечающих наибольшему возможному для г-ого участка вдоль скользящей дислокации уровню внешнего напряжения сдвига На основании анализа массива значений (г"и), находилось значение внешнего напряжения (т"')' - тт{(гп""1х),}, при котором проводилось окончательное построение всей равновесной конфигурации скользящей дислокации Уровню (г'")* отвечает предельное значение внешнего напряжения сдвига, при котором закрепление скользящей дислокации на выбранной последовательности препятствий оказывается возможным Наряду с построением равновесной конфигурации скользящей дислокации, отвечающей уровню внешнего напряжения (г'")', проводилась индикация наиболее слабых к уровню т"и участков вдоль скользящей дислокации, открепление от соответствующих лесных дислокаций, продвижение скользящей дислокации вплоть до встречи с новыми дислокациями леса, формирование нового массива дислокационных узлов вдоль скользящей дислокации, после чего описанная выше процедура возобновлялась при (г*"„) = (г™)* и т д
При использовании динамического подхода, среда, в которой движутся дислокации, считается изотропной, обладающей свойством вязкости Для начала движения дислокации в рассматриваемой среде необходимо преодолеть силу, подобную максимальной силе трения покоя в механике
Уравнение движения берется в виде
= , (2)
где ^ = - сила самодействия в приближении линейного натяжения
(Д- радиус кривизны элемента дислокации в точке, для которой рассчитыва-
ется величина силы),
= Вч = Вди/81 - сила вязкого трения,
а ( х + 1/2 х-1/2
1-^--п-^ _ сила взаимо-
Хх + 112)2+у2) [(х-1/2)2+у2))
действия дислокаций, обусловленная действием лесных дислокаций, х, у -координаты точки сегмента скользящей дислокации, в которой рассчитывается значение силы,
I?СТ" - сила, обусловленная внешней нагрузкой, Ъа %\ncot
= +Рех + Рт) - сила, обусловленная наличием стартово-
го напряжения
_ та ,
Согласно сделанным оценкам, — «1, поэтому инерционным членом
Bv
в (2) можно пренебречь
Решение уравнения (2), записанное в относительной системе координат, связанной с скользящей дислокацией, для случая неподвижных центров закрепления концов дислокационного сегмента имеет вид
4
Здесь
5т 1-
-■>/81 > р к
соз(о^) (3)
^ =(-!)* 4СТ° + --5-+ Зт(Х)-5„(УР1) \со$(а„Х)сП.
1 я(2* + 1) I Г ) н 2Ла„,т„) V
я(2Л +1) V
_п(2к + 1) В1 ак ~ , >'~
Я
2Я(\0,т0) ""Г ' )
I Ъосг
Данное решение позволяет по конфигурации сегмента скользящей дислокации в момент времени гд определить смещение точек, образующих сегмент за время с!т, и построить конфигурацию сегмента для момента времени т$+с1т (соответственно конфигурации 1 и 2 на рис 2)
В данной работе предложены два алгоритма, позволяющие моделировать поведение дислокационного сегмента в случае, когда закрепляющие лесные дислокации совершают гармонические колебания по закону
у = уъ$т(ш), (4)
где у - смещение точки закрепления относительно начального положения, уд - амплитуда смещения лесной дислокации, которая находится из уравнения движения винтовой дислокации под действием ультразвука заданной амплитуды, ст - частота ультразвука
Рис. 2. Схема последовательных конфигураций дислокационного сегмента в моменты времени гд (1) и гд +8т (2)
Первый из предложенных алгоритмов учитывает возможность изменения формы сегмента из-за смещения лесных дислокаций, служащих центрами закрепления Согласно этому алгоритму сначала аналитически находится конфигурация дислокационного сегмента в момент времени го + 8т В дальнейшем принимается, что положение центральной точки сегмента не зависит от смещения точек закрепления Смещение всех остальных точек 5, принимается пропорциональным расстоянию вдоль сегмента от соответствующей точки до его центральной точки
где - смещение точки закрепления, Ху ,у} - координаты в лабораторной
системе координаты _)-ой точки сегмента (рис 3)
Это решение может быть применено только для случая, когда прогиб сегмента меньше половины расстояния между точками закрепления Моделирование с использованием данного алгоритма требует большого объема расчетов, а следовательно и машинного времени В связи с этим был разработан другой алгоритм
*1=5П-\ 7ГТ
/=*
к
Рис 3 К расчету смещения точек сегмента при изменении положения точек закрепления (точки закрепления колеблются в противофазе) О, N-1 - точки закрепления, к - центральная точка сегмента, I - произвольная точка сегмента
Второй алгоритм изменение формы сегмента, обусловленной движением лесных дислокаций, не учитывает Он основан на том, что расчет формы сегмента производится в системе координат, привязанной к сегменту, названной выше относительной Как видно из решения (3), смещение произвольной точки сегмента и, зависит от I (см рис 2) Это позволяет реализовать следующий алгоритм также, как в первом случае, аналитически для конфигурации дислокационного сегмента в момент времени т0 рассчитывается смещение точек за время (рис 2) Затем изменяется положение точек закрепления на величину, рассчитанную с использованием (4) за время Это приводит к тому, что относительная система координат поворачивается и растягивается пропорционально изменению / Те аналитические расчеты изменения формы сегмента производятся в относительной системе координат, с учетом меняющегося /, а затем полученная конфигурация
Рис 4 Последовательные конфигурации дислокационного сегмента, полученные в соответствии со вторым алгоритмом
переводится в лабораторную систему координат Последовательные конфигурации дислокационного сегмента, полученные с помощью второго алгоритма, представлены на рис 4 в лабораторной системе координат
Последовательные конфигурации дислокационного сегмента на рис 4 даны с шагом 10£т Как видно из рис 4, в конфигурациях 2, 3, 4, когда точки закрепления сегмента сближаются, стрела прогиба дислокационного сегмента меньше, чем в конфигурациях 6, 7, 8, когда расстояние между точками закрепления увеличивается Такие изменения стрелы прогиба связаны с тем, что сила самодействия обратно пропорциональна радиусу кривизны сегмента, который тем больше, чем больше расстояние между точками закрепления 1
Для сравнения этих двух алгоритмов были исследованы зависимости амплитуды ультразвука от частоты, необходимой для срабатывания источника Франка-Рида, центры закрепления которого двигаются по гармоническому закону Выяснилось, что результаты, полученные с использованием этих алгоритмов, в области частот 60 100 кГц и длин источников менее 5 мкм отличаются менее чем на 1,5% Поэтому в рамках данной диссертации использован второй алгоритм, как менее расчетоемкий
Для анализа поведения скользящей дислокации разработан алгоритм, позволяющий моделировать ее прохождение через лес дислокаций В данном алгоритме скользящая дислокация разбивается закрепляющими ее лесными дислокациями на совокупность гибких сегментов В соответствии с этим на каждом шаге моделирования рассчитывается изменение формы для каждого сегмента По мере продвижения дислокации может произойти открепление сегмента от одних центров и закрепление на других Критерием закрепления дислокационного сегмента на дислокации леса служит попадание сегмента в 5 окрестность лесной дислокации Критерием открепления служит угол между соседними сегментами в точке закрепления - если угол меньше 10°, то считается, что скользящая дислокация освободилась от данного центра закрепления за счет аннигиляции соответствующих участков сегмента
Рассматривается также алгоритм моделирования поведения ансамбля лесных дислокаций в ультразвуковом поле Уравнение движения лесной дислокации берется аналогичным уравнению (1), в котором помимо сил вязкого трения, внешней силы и силы типа сухого трения, учитывается поле сил взаимодействия лесных дислокаций между собой Сила, действующая на дислокацию леса, бралась как сумма сил парного взаимодействия и определялась по формуле-
У
и I X
Здесь /7 =------—— - сила взаимодействия 1-ой и ,|-ой
2 п (Х2 + у^-у-
дислокации леса, где х - расстояние между дислокациями по плоскости «легкого скольжения», А - расстояние между соответствующими плоскостями «легкого скольжения»
Решение уравнения движения лесной дислокации может быть записано в виде
,т . (Увзаимен) У вне ш У старт /С\
<№ =---, (5)
о
где сО/ - смещение лесной дислокации в плоскости легкого скольжения
Зная положение дислокаций на момент времени t, можно рассчитать, используя (5), их смещение за время Ж Как будет показано ниже, ансамбль леса перестраивается под действием ультразвука
В третьей главе рассматривается влияние ультразвука на основные элементы дислокационной структуры
В данной главе показывается, что воздействие гармонической нагрузки на единичную дислокацию может вызывать ее ангармоническое движение
При воздействии ультразвукового поля на дислокационный диполь было установлено, что положения, относительно которых дислокации диполя колеблются, зависят как от поля сил взаимодействия дислокаций, так и от параметров ультразвука Стационарное состояние динамического диполя можно характеризовать углом ср* между плоскостью скольжения дислокаций и плоскостью, содержащей линии, относительно которых дислокации диполя колеблются
В отсутствие ультразвука, когда дислокации в стационарном состоянии неподвижны, ф* = ср, т е углу между плоскостью скольжения и габитусной плоскостью дислокаций диполя Равновесным состояниям краевого диполя, обусловленным силами взаимодействия, соответствуют углы ср=45° (положение устойчивого равновесия) и ф=90° (положение неустойчивого равновесия)
В присутствии внешнего воздействия угол ф изменяется При постоянной во времени нагрузке ф = ср (а), т е изменяется в зависимости от о При знакопеременной нагрузке ф=ф (а0), где а0 - амплитуда ультразвука Рис 5 иллюстрирует изменение конфигурации диполя в зависимости от амплитуды ультразвука а0 при различных значениях расстояния между плоскостями
1, 2) 3 41
О 0,1 0,3 0,4 Ов 0,7 0,9 1,0 1,2 1 3 1 5 1 6 18 1 9 2,1 2,2 2:4
Рис 5 Изменение структуры динамического диполя в зависимости от амплитуды для диполей, имеющих различное плечо й=20,30,40, 50 мкм
скольжения диполя Ь Из рис 5 видно, что можно выделить несколько характерных для структуры диполя диапазонов амплитуд ультразвука Первый диапазон - это те амплитуды, при которых дислокации диполя покоятся, движение не происходит из-за наличия стартовых напряжений Второй диапазон амплитуд - это когда диполь находится в состоянии динамического равновесия, причем угол (р* = 45° Подчеркнем, что дислокации диполя в этом диапазоне амплитуд колеблются около положений равновесия, характерных для устойчивого состояния в отсутствие внешнего возбуждения Однако данная конфигурация диполя сохраняется лишь в сравнительно узком интервале амплитуд В третьем диапазоне амплитуд конфигурация диполя переходит из одного стационарного состояния в другое Обращает на себя внимание тот факт, что в определенных интервалах амплитуд ультразвука ф*=ср*(ст) изменяется плавно, а при достижении некоторой критической амплитуды 5° {И) функция ф*(о°) претерпевает скачок
Дислокационные сегменты в ультразвуковом поле при малых амплитудах совершают периодические движения, которые можно отнести к ангармоническим из-за существования стартового напряжения
По мере увеличения амплитуды ультразвука или при добавлении постоянной составляющей поля дислокационный сегмент может терять устойчивость и начинается размножение дислокаций по механизму Франка-Рида.
Поведение источника Франка-Рида рассматривали при следующих режимах нагружения
1 Лесные дислокации неподвижны, а в плоскости источника действу-
ют постоянная и знакопеременная составляющие нагрузки
2 На скользящую дислокацию действует постоянная сила, а на дислокации леса действует ультразвук, в результате чего последние совершают колебания
3. На скользящую дислокацию действует постоянная и знакопеременная составляющие силы, а лесные дислокации колеблются под действием ультразвука
При первом из приведенных выше режимов нагружения ультразвук играет роль спускового механизма Напряжение ег в плоскости дислокационного сегмента, служащего источником, имеет постоянную crconsl и меняющуюся по гармоническому закону <jq компоненты
cT = a-const+a0sm(wi) (6)
Началом работы источника Франка-Рида служит момент, когда дислокационный сегмент теряет устойчивость В случае действия только постоянной составляющей напряжение потери устойчивости дислокационного сег-
- ; * Gbl мента длинои / определяется по формуле акр - —-— = & const
Пусть условия нагружения таковы, что о const меньше акр При совместном действии постоянной и знакопеременной составляющих напряжения,
Рис 6 Зависимость егсот, (МПа), необходимого для срабатывания источника Франка-Рида, от длины источника при амплитудах ультразвука а - а* = О, Ь - а^ = 0,4 МПа, с - а" = 0,8 МПа (частота ультразвука 60 кГц)
*
начиная с пороговой амплитуды ультразвука сгд , в определенный момент * * * *
времени / достигается напряжение а = оС0т1 + а"о 31п(гй ), при котором дислокационный сегмент теряет устойчивость и начинается размножение дислокаций по механизму Франка-Рида Согласно экспериментальным данным, сгсот1 может принимать даже нулевое значение, т е источник может генерировать дислокации под действием только знакопеременной составляющей нагрузки (рис 6)
Таким образом, ультразвук, действующий в плоскости источника, способствует срабатыванию источников, длина которых меньше / = —^—
Поскольку при своем движении дислокация испытывает сопротивление типа вязкого трения, появляется характерное время генерации источником замкнутой петли Эффект понижения напряжения аС0П5[, необходимого для
Рис 7. Увеличение количества источников, срабатывающих при заданной оС(,„,, когда в плоскости источника действует знакопеременная составляющая
срабатывания источника, обнаруживает зависимость от частоты ультразвука Если концы источника закреплены лесными дислокациями разного знака, то под действием ультразвука они будут двигаться в противофазе. В
(рис 7)
/
ж
/
Ж
I
У
Ж*)
нулевом приближении эту ситуацию можно описать, как изменение эффективной длины источника 1эф = /о + сНыпаЛ В
определенные доли периода
эф
въ
может 2
превосходить
и источник начнет
' кр
Рис 8 Схема изменения эффективной
длины источника за счет колебания лесных дислокаций в противофазе.
а — длина источника уменьшается,
б - длина источника увеличивается
Вероятность появления источников, концы которых колеблются в противофазе, зависит от соотношения между положительными и отрицательными дислокациями в ансамбле леса Эта вероятность максимальна при соотношении положительных и отрицательных дислокаций 1 1
Для случая, когда ультразвук действует и на лесные, и на скользящую дислокации, возможно как упрочнение, так и разупрочнение, а какой из эффектов будет наблюдаться, зависит от соотношения фаз (см рис 9) В один полупериод источник удлиняется; сила, необходимая для его срабатывания, уменьшается, в следующий полупериод длина источника уменьшается, а сила, необходимая для его срабатывания, возрастает. Если изменение длины источника и силы, действующей в его плоскости, в фазе (рис 9а), то будет наблюдаться эффект пластификации, а в противном случае (рис 96) будет происходить упрочнение.
генерировать дислокации при напряжениях и, меньших акр, характерного для состояния источника с неподвижными точками закрепления Это можно видеть из
рис. 8
1
Г
\ ! /
■ > _____я
'1 к *
Рис 9 Возможные соотношения напряжения, действующего в плоскости источника, и расстояния между закрепляющими дислокациями в зависимости от времени
Как показало моделирование эффект упрочнения обнаруживает зависимость не только от параметров ультразвука, но и от длины источника так для коротких источников (длиной менее 3,5 мкм) с ростом амплитуды ультразвука растет эффект упрочнения Для длинных источников, у которых относительное изменение эффективной длины меньше, зависимость эффекта упрочнения от амплитуды ультразвука не монотонна
На основании вышесказанного можно сделать вывод, что для кристаллов, в которых количество лесных дислокаций разного знака одинаково, с ростом плотности дислокаций эффект изменения пластических свойств кристаллов, под действием ультразвука, связанный с размножением дислокаций, будет уменьшаться.
Глава четвертая содержит описание особенностей полигонизации в ультразвуковом поле Как было описано в предыдущей главе 3, в ультразвуковом поле образуются динамические дислокационные диполи Дальнейшее моделирование показало, что в ультразвуковом поле формируются и более сложные упорядоченные дислокационные ансамбли, в том числе мультипо-ли и полигональные границы блоков При формировании мультиполей и границ блоков «затравками» служат дислокационные диполи Число и структура образующихся в ультразвуковом поле полигональных границ зависят от частоты, амплитуды и времени действия ультразвука, а также от начальной плотности и распределения дислокаций, в том числе от соотношения числа дислокаций разного знака в дислокационном ансамбле
В рамках данной диссертационной работы впервые была предпринята попытка произвести моделирование и дальнейший анализ поведения леса краевых дислокаций в ультразвуковом поле Одним из основных условий при моделировании было соответствие количества дислокаций на модельной площадке реальным плотностям дислокаций в кристаллах В связи с этим моделирование проводилось с ансамблем дислокаций от 200 до 300 шт на модельную площадку, что соответствовало средней плотности дислокаций р = 8 106 см'2 На рис 10 представлено начальное расположение дислокаций на модель- _ , „ тт
ной площадке Координаты 10 Начальное расположение дислока-
начального расположения ции на сдельной площадке
дислокаций на модельной площадке задавались по закону случайных чисел
Как показало моделирование, конечный результат существенно зависит от начального распределения дислокаций на модельной площадке и от порядка проведения ЭВМ экспериментов (чередование воздействия ультразвука и релаксации)
Если непосредственно после начальной расстановки дислокаций на модельной площадке подвергнуть систему дислокаций действию ультразвука, то в результате получим картину, представленную на рис. 11 Как видно
из рис И, образовалась ячеистая структура, у которой границы блоков приобретают тонкую структуру Если теперь полученную в результате воздействия ультразвука дислокационную структуру подвергнуть релаксации, то характер расположения ячеек сохраняется, а структура границ блоков меняется То есть, в процессе релаксации дислокационной структуры, созданной ультразвуком, изменяются расположения дис-Рис 11 Распределение дислокаций на модель- локаций, принадлежа-ной площадке, полученное в результате дейст- щих границам ячеек вия ультразвука (ст° = 1,5 МПа, £ = 100 кГц, Результаты, полу-
г = юот) ченные в процессе моде-
лирования, совпадают с
экспериментальными данными, полученными на кафедре молекулярной физики МГУ имени МВ Ломоносова Результат электронной микроскопии кристалла Сс18, подвергшегося облучению ультразвуком, представлен на рис 12
Число сформировавшихся дислокационных стенок и их структура зависят от параметров ультразвука Это видно из графиков зависимости процента дислокаций, вошедших в упорядоченные структуры типа дислокационных стенок от времени действия ультразвукового поля (рис 13) Кривые 14 рис 13 соответствуют проценту дислокаций, вошедших в дислокационные стенки различной длины Кривая 1 - относится к структурам, состоящим из числа дислокаций от 3 до 5 Кривая 2 описывает процент дислокаций, во-
шедших в стенки от б до 9 штук. Кривая 3 показывает изменение во времени процента дислокаций, формирующих стенки от 10 и более штук. Обший процент дислокаций, участвующих в формирование границ ячеистой структуры, характеризует кривая 4.
Рис. 12. Изменение плотности дислокаций после обработки ультразвуком
На основании графиков рис. 13 можно заключить, что зародыши дислокационных стенок, состоящих из 3 ~ 5 дислокаций, формируются за время порядка 10 периодов ультразвука. В дальнейшем, с ростом времени действия ультразвука, процент вошедших дислокаций в комплексы, содержащие от 3 до 5 дислокаций, убывает.
Формирование стенок, состоящих из 6 и более дислокаций (кривые 2 и 3 на рис. 13), начинается по прошествии 20 периодов. Число дислокаций, вошедших в такие стенки, с течением времени возрастает.
Обращает на себя внимание, что по прошествии времени действия ультразвука до 150Т процент дислокаций, вошедших в структуры 3-5, 6-9, 10 - более не достигает стационарного значения. Процесс перестройки дислокационной структуры продолжается в течение всего времени действия ультразвука. Однако, при этом общий суммарный процент дислокаций, вошедших в дислокационные стенки, с течением времени достигает «насыщения» - колеблется в интервале от 70% до 80%. Как видно из рис. 13, характерное вре-
мя, в течение которого «насыщение» достигается, порядка t = 50Т Также видно из графиков рис 13, что в первые 40Т преимущественно организуются дислокационные комплексы, которые содержат от 3 до 5 дислокаций Судя по данным, приведенным на рис 13, для стенок, содержащих 6-9 дислокаций, оптимальное время их формирования имеет порядка 80Т
Рис 13 Зависимость количества дислокаций (в %), вошедших в упорядоченные структуры, от времени действия ультразвука (</ = 0,5 МПа, / = 90 кГц, /=150 Т, N = 200 дисл)
Глава пятая содержит результаты моделирования процесса прохождения пробной скользящей дислокации через колеблющийся дислокационный лес
При моделировании прохождения пробной дислокацией модельной площадки были рассмотрены следующие ситуации
1 Прохождение пробной дислокации через неподвижный лес, дислокации которого расположены по случайному закону
2 Прохождение пробной дислокации через лес, дислокации которого совершают гармонические колебания
3 Прохождение пробной дислокации через ансамбль лесных дислокаций, структура которого получена следующим образом в исходном состоянии дислокации леса располагаются по случайному закону, затем моделируется процесс перераспределения лесных дислокаций под действием ультразвука с учетом поля сил их взаимодействия и релаксации структуры ансамбля после прекращения воздействия ультразвука
4 Прохождение пробной дислокации через колеблющийся лес, структура которого получена по способу, описанному в пункте 3
Для первой и второй ситуации моделирование проводилось как в квазистатическом, так и в динамическом приближении
Для выяснения влияния величины амплитуды колеблющихся дислокаций леса на характеристики процесса движения скользящих дислокации было проведено 8 вариантов моделирования при следующих значениях амплитуды А = 0; 0,05 Л, 0,10 Л, 0,15 Л, 0,20 Л, 0,25 Л, 0,30 Л, 0,35 Я; где Л - среднее расстояние между дислокациями леса, причем для каждого варианта задачи моделирование проводилось для одного и того же начального размещения дислокаций леса Плотность дислокаций леса принималась равной р = 4 Ю10 м"2. Численные значения параметров выбирались применительно к кристаллам цинка Пример прохождения краевой дислокации через дислокационный лес при амплитуде колебания лесных дислокаций А = 0,15 X представлен на рис 14
Рис 14 Последовательные положения скользящей дислокации, двигающейся через колеблющийся дислокационный лес Стрелками обозначены траектории колеблющихся дислокаций леса А = 0,15 X
Результаты проведенного моделирования позволяют установить ряд особенностей в процессах движения скользящих дислокаций через колеблющийся с разными значениями амплитуды дислокационный лес
Было установлено, что возможность дислокаций леса совершать периодические колебания приводит к понижению критического напряжения прохождения пробной дислокацией плошадки моделирования, причем величины относительного разупрочнения при различных значениях амплитуды колебаний дислокации леса сильно отличаются друг от друга Анализ показал, что наилучшими приближениями зависимостей критического напряжения прохождения и относительного разупрочнения от величины амплитуды колебаний являются линейные зависимости с изломом в окрестности значения амплитуды 0,17 Я (рис 15) Весьма примечательно, что излом лежит внутри интервала амплитуд 0,15 Я - 0,20 Я, при прохождении которого также наблюдаются качественные изменения в характере процесса прохождения скользящих дислокаций через колеблющийся лес В самом деле, увеличение амплитуды колебаний в интервале амплитуд 0<А<0,15 Я сопровождается, в основном, уменьшением предельного напряжения сдвига, при котором данная последовательность колеблющихся дислокаций леса может удерживать скользящую дислокацию При этом конфигурации скользящих дислокаций при различных значениях амплитуды в значительной степени повторяют друг друга
При движении скользящих дислокаций через лес дислокаций, колеблющихся с амплитудой А>0,20 Я, указанных повторений в последовательностях колеблющихся дислокаций леса, удерживающих скользящую дислокацию, не наблюдается
Применение квазистатических моделей при рассмотрении вопросов, связанных с влиянием ультразвука, может дать только качественный результат Это связано с тем, что квазистатические модели не учитывают вязкость среды. При воздействии ультразвука на дислокационный лес важную роль
0,35 А/£
Рис 15 График зависимости критического напряжения от амплитуды ультразвука
играет не только амплитуда ультразвука, но и время, когда дислокация подходит к центру закрепления, т к в один период времени расстояния между соседними лесными дислокациями могут увеличиваться как в фазе, что будет способствовать понижению напряжения, необходимого для преодоления модельной площадки, так и в противофазе, а в этом случае трудно предугадать как поведет себя скользящая дислокация. В квазистатической модели считается, что дислокация подходит в наиболее удобный момент Данное утверждение подтверждается и результатами проведенных модельных экспериментов Как уже было сказано, конфигурации скользящей дислокации, полученные для разных амплитуд колебания лесных дислокаций, совпадают, только уменьшается напряжение, необходимое для преодоления модельной площадки Такие совпадения конфигураций связаны с тем, что увеличение амплитуды колебаний приводит к увеличению расстояний между лесными дислокациями; уменьшается сила самодействия и не более Поэтому для исследования влияния ультразвука на процесс движения скользящей дислокации предпочтительнее динамический подход
Установлено, что когда дислокации леса неподвижны, напряжение, необходимое для преодоления площадки, пропорционально а-у[р, что согласуется с ранее полученными, хорошо известными результатами Зависимость напряжения, необходимого для преодоления модельной площадки пробной дислокацией, от плотности дислокаций леса можно видеть на рис 16
е^МПа
0 50ИЕ*03 1С0Е-Р 1 5С£*0» 2ЮЕ-В4 1 ЗЗЕ-М ЗГКЗЕ'М 3,506 - 01
Рис 16 Зависимость напряжения, необходимого для преодоления площадки от корня плотности дислокаций
При колебании лесных дислокаций по гармоническому закону (ситуация 2), напряжение, необходимое для преодоления модельной площадки, уменьшается по сравнению с таковым для неподвижного леса той же плотности Величина эффекта пластификации при воздействии ультразвука зависит от плотности дислокаций, с ростом плотности эффект пластификации увеличивается При плотности дислокаций 109 см'2 эффект составляет 60%, при плотности порядка 106 см"2 - 25%. Также появляется зависимость критического напряжения от соотношения положительных и отрицательных дислокаций в ансамбле леса Максимальная пластификация имеет место, когда число положительных дислокаций равно числу отрицательных. С изменением этого соотношения наблюдается уменьшение эффекта пластификации (рис 17)
Лет
------ — - - - - - 1 I
* ! М '
А к
! I
* . " !
л
I
ООСЕ + ОО 5 ООЕ + ОЭ 1 00Е»04 1,50Е«04 2 ОСЕ+04 2 50Е>04 300Е-Ю4 3 5 0Е»04 I- -1
Рис 17 Отношение изменения напряжения, необходимого для преодоления модельной площадки к напряжению преодоления модельной площадки, для разных соотношений положительных и отрицательных дислокаций (1 - отношение 50 50,2-40 60,3-20 80)
Для случая, когда ансамбль дислокаций леса предварительно подвергся ультразвуковой обработке, наблюдается понижение напряжения, необходимого для преодоления модельной площадки, на величину порядка 10-15% по сравнению с тем, что наблюдается в случае, когда дислокации леса расположены по случайному закону Если при прохождении пробной дислокации через указанный ансамбль включить ультразвук, то напряжение преодоления модельной площадки практически не отличается от представленного на рис 17
На основании результатов, полученных в 5 главе, можно резюмировать
- напряжение, необходимое для преодоления пробной дислокацией модельной площадки, пропорционально а->[р, где р - плотность дислокаций в кристалле,
- коэффициент прозрачности леса а составляет от 0,8йЬ до 0,95СА, а для случая предварительно обработанного ультразвуком леса - от 0,7йЬ до 0,906,
- напряжение, необходимое для преодоления модельной площадки, понижается в случае колебания дислокаций леса и может достигать 40% от напряжения, необходимого в случае неподвижного леса,
- эффект пластификации зависит от плотности дислокаций леса и увеличивается с ее ростом
Заключение В работе рассмотрены возможные механизмы изменения пластичности кристаллов в присутствии ультразвукового поля Основное внимание обращено на роль дислокаций леса в процессе пластификации кристаллов под действием ультразвука
Методом компьютерного моделирования дислокационных процессов изучены эффекты, связанные с дислокациями леса, в условиях сложнона-груженного состояния кристаллов, при котором на скользящие дислокации действует постоянная нагрузка, а на дислокации леса - знакопеременная ультразвукового диапазона частот
Проведено моделирование процессов скольжения дислокаций в поле колеблющихся дислокаций леса и генерации дислокаций источниками Франка-Рида, центрами закрепления в которых служат колеблющиеся дислокации леса
Показано, что в отличие от хорошо известного эффекта упрочнения кристаллов за счет взаимодействия скользящих дислокаций с дислокациями леса может наблюдаться и обратный эффект пластификации кристаллов при условии, если дислокации леса колеблются
ОБЩИЕ РЕЗУЛЬТАТЫ РАБОТЫ И ВЫВОДЫ
Методом математического моделирования детально изучены механизмы влияния ультразвука на эволюцию дислокационной структуры и вызванные ею изменения физических свойств материалов Рассмотрено влияние ультразвука на процесс взаимодействия дислокаций различных систем Впервые исследованы дислокационные взаимодействия в ультразвуковом поле с использованием уравнений динамики дислокаций Проведено сравнение результатов, полученных в приближении квазистатической и динамической моделей Применение динамического подхода позволило установить новые явления и получить более точные количественные оценки эффектам,
ранее рассматривавшимся в квазистатическом приближении Так, например, впервые рассмотрена генерация дислокаций источниками Франка-Рида, концы которых совершают вынужденные колебания, что потребовало учета фазовых соотношений Основное внимание в диссертации уделено двум механизмам, ответственным за изменение пластических свойств материалов в ультразвуковом поле, а именно размножению дислокаций по механизму Франка-Рида и прохождению дислокации через колеблющийся лес, при этом впервые учитывается возможность влияния ультразвука и на скользящую дислокацию, и на дислокации леса Полученные в работе данные моделирования совпали с экспериментальными данными других авторов и объяснили многие ранее известные экспериментальные результаты, что дает основание использовать их для прогнозирования поведения реальных материалов в условиях их эксплуатации
На основании полученных результатов и их анализа могут быть сделаны следующие выводы
1 Методами компьютерного моделирования проведен анализ и выявлены основные закономерности поведения элементов дислокационной структуры в ультразвуковом поле Анализ проведен с помощью специально разработанных модели, алгоритма и пакета программ для проведения компьютерных экспериментов
2 Ультразвуковое воздействие на неупорядоченные дислокационные ансамбли приводит к перераспределению дислокаций, в результате которого формируются упорядоченные динамические дислокационные структуры — диполи, мультиполи и дислокационные стенки В совокупности наблюдавшиеся элементарные процессы интерпретируются как самоорганизация дислокационной подсистемы (полигонизация, образование ячеистой структуры)
3 Самоорганизация дислокационной подсистемы приводит к достижению состояния, в котором суммарная доля дислокаций, вошедших в упорядоченные дислокационные структуры, не изменяется Время достижения этого состояния имеет порядок 102 периодов и зависит и от частоты, и от амплитуды ультразвука
4 Впервые показано, что в ультразвуковом поле формируется система субграниц, каждая из которых образована дислокациями одного знака, но при этом любые две соседние субграницы образованы дислокациями другого знака
5 В границах блоков, сформировавшихся в ультразвуковом поле, дислокации не лежат в одной плоскости, и границы имеют тонкую структуру, т е содержат дислокационные диполи и мультиполи Этим они отличаются от дислокационных стенок, формирующихся при отжиге в отсутствие ультразвукового поля
6 Наложение ультразвукового поля запускает процесс размножения дислокаций по двум каналам Одним из каналов является активация источников Франка-Рида, но при этом благодаря действию ультразвука активация каждого данного источника происходит при меньшем значении постоянного напряжения сдвига (эффект ультразвука сводится к действию спускового механизма) Вторым каналом является увеличение числа источников Франка-Рида
7 Прохождение скользящей дислокации сквозь дислокации леса в условиях наложения ультразвукового поля может давать в результате либо упрочнение, либо пластифицирование материала Знак эффекта определяется соотношением фаз колебаний скользящей дислокации и дислокаций леса
8 В случае колебания дислокаций леса по гармоническому закону, напряжение, необходимое для продвижения скользящей дислокации через единичную площадку, понижается и может составить лишь 60% от величины напряжения, необходимого для прохождения дислокации сквозь лес неподвижных дислокаций
9 Совокупность данных, полученных в результате проведения компьютерных экспериментов, показывает, что при динамическом подходе удается учитывать временные зависимости (эффекты вязкости) и частоту ультразвука Это позволит использовать эти результаты на предприятиях, применяющих ультразвук при получении лент и проволок из непластичных материалов, например, ОАО «Биметалл» и ЗАО «Аметист-Стан» При квазистатическом подходе указанные эффекты обеспечить невозможно.
Основное содержание диссертации отражено в следующих работах
1 Дегтярев В Т, Логинов Б М Исследование влияния тонкой структуры дальнодействующих полей напряжений ансамбля колеблющихся дислокаций на характеристики процесса прохождения скользящих дислокаций //ВИНИТИ -1985 -№7043-В -48с
2 Дегтярев В Т, Логинов Б М, Тяпунина Н А Исследование зависимости характеристик процесса движения скользящих дислокаций через ансамбль колеблющихся дислокаций от величины амплитуды колебаний //ВИНИТИ - 1986 -№1217-Д.-55 с
3 Дегтярев В Т, Логинов Б М, Еремеев А В Исследование зависимости разупрочнения дислокационных ансамблей от амплитуды периодической нагрузки // Физика прочности и пластичности металлов и сплавов- Материалы XI Всесоюзной конференции -Куйбышев, 1986 - С.201
4 Дегтярев В Т, Логинов Б М, Тяпунина Н А Моделирование скольжения дислокаций через дислокационный лес колеблющихся дислокаций в
кристаллах с ГПУ структурой // Кристаллография - 1987 - Т32, №4 -С 967-971
5 Дегтярев В Т, Логинов Б M Моделирование движения дислокаций через колеблющийся дислокационный лес в кристаллах NaCl // Структура дислокаций и механические свойства металлов и сплавов1 Материалы IV Всесоюзного семинара - Свердловск, 1987 - С 52-53
6 Дегтярев В Т., Логинов Б M Моделирование движения дислокаций через колеблющийся лес дислокаций с учетом дальнодействующих полей напряжений в кристаллах с гексагональной плотноупакованной решеткой //Физика металлов и металловедение - 1987 - Т 64, №3 - С 608-610
7. Дегтярев В Т, Логинов Б M, Еремеев А В. Зависимость статистических характеристик процессов движения скользящих дислокаций через дислокационные ансамбли от параметров площадки моделирования // Моделирование на ЭВМ дефектов в кристаллах - Л ФТИ АН СССР, 1988 - С.50-55
8 Математическое моделирование процессов движения линейных дефектов в ГПУ монокристаллах в условиях действия периодической нагрузки /ВТ Дегтярев, Б M Логинов, А В Еремеев, В Р Туганов // Техническая кибернетика, математика, вычислительная техника Материалы 33 Международного коллоквиума - Ильменау (ГДР) РС1, 1988 - С 63-66
9 Дегтярев В Т., Логинов Б M, Тяпунина H А Моделирование движения дислокаций через лес колеблющихся дислокаций с учетом дальнодействующих полей напряжений в кристаллах с решеткой NaCl // Кристаллография -1988 -ТЗЗ,№1 -С 163-166
10 Моделирование процессов движения дислокаций через ансамбли колеблющихся дислокаций /ВТ Дегтярев, Б M Логинов, А В Еремеев, H А Тяпунина // Пластическая деформация материалов в условиях внешних энергетических воздействий. Материалы IX семинара - Новокузнецк, 1988 - С 27-28
11 Дегтярев В Т, Логинов Б M, Еремеев А В Исследование влияния параметров площадки моделирования на статистические характеристики процессов движения скользящих дислокаций через хаотические дислокационные ансамбли//ВИНИТИ -1989 - № 2358-В89. -25с
12 Дегтярев В.Т., Логинов Б M Исследование процессов прохождения скользящей дислокации через ансамбль хаотически распределенных призматических дислокационных петель, подверженных воздействию периодической нагрузки//ВИНИТИ - 1989 -№620-В89 -20с
13 Degtyarev V Т, Loginov В M, Shvedov RN Computer simulation of the glide dislocation movement through chaotic ensembles of prismatic dislocation loops under complex loading conditions // Modelling, measurement & control -1994 - Vol 56, N 3 - P 9-20
14 Дегтярев B.T, Логинов Б M, Ямпольский В И Анализ дефектной структуры кристаллов на основе моделирования процессов комплексного
нагружения I I 165 лет МГТУ им НЭ Баумана Материалы НТК - Калуга, 1995 -Ч 1 -С.91
15.Акустопластический эффект, обусловленный взаимодействием дислокаций / В.Т. Дегтярев, АЮ Лосев, Ф.А Плотников, H.A. Тяпунина // Imperfection interaction and anelasticity phenomena m solids (HAPS - 10) Материалы X Международной научно-практической конференции - Тула, 2001. -С. 34
16 Влияние колебаний лесных дислокаций на движение скользящей дислокации /ВТ Дегтярев, А Ю Лосев, Ф А Плотников, H А. Тяпунина //Дефекты структуры и прочность кристаллов Материалы Всероссийской конференции, посвященной 100-летию со дня рождения академика Г.В. Кур-дюмова - Черноголовка,2002 -С 235
17.Дегтярев В Т Самоорганизация дислокационных ансамблей в ультразвуковом поле // Материалы электронной техники Известия ВУЗов - 2004 -№1 -С 34-37
18 Дегтярев В Т, Лосев А Ю Динамические дислокационные структуры в ультразвуковом поле // Актуальные проблемы прочности. Сборник трудов XLII Международной конференции - Калуга, 2004. - С. 135
19 Пластификация кристаллов ультразвуком, обусловленная взаимодействием дислокаций / ВТ. Дегтярев, АЮ Лосев, ФА Плотников, H А. Тяпунина // Известия Тульского государственного университета Сер Физика.-Тула, 2003.-Вып 3 -С. 3-10
20 Полигонизация в ультразвуковом поле /ВТ Дегтярев, А.Ю. Лосев, Ф А Плотников, H А Тяпунина // Известия РАН Сер Физическая - 2004 -Т. 68, №10 - С. 1516-1517.
21 Дегтярев В Т, Лосев А Ю , Плотников Ф А Динамические дислокационные структуры в ультразвуковом поле, диполи и Триполи // Материаловедение - 2004 - №7. - С. 8-12
22. Дегтярев В Т., Лосев А Ю., Тяпунина H А Влияние ультразвука на процесс генерации дислокаций источником Франка-Рида // Деформация и разрушение материалов. - 2005 - №3. - с 18-21
23 Дегтярев В.Т, Лосев А.Ю , Плотников Ф А. Перераспределение неупорядоченных дислокационных ансамблей в ультразвуковом поле // Наукоемкие технологии. - 2005 - № 3-4, т 6. - С 5-8
24 Акустопластический эффект и факторы на него влияющие /ВТ Дегтярев, А Ю Лосев, Ф А. Плотников, H А Тяпунина // Деформация и разрушение материалов -2005 -№7.-С23-27
25 Дегтярев В.Т., Лосев А Ю Моделирование генерации дислокаций источником Франка-Рида в условиях акустопластического эффекта // Актуальные проблемы прочности. Материалы XLIV международной конференции. - Вологда, 2005 -С 61
26 Дегтярев В Т, Лосев А Ю Процессы, сопровождающие акустопла-
стический эффект // Актуальные проблемы прочности Материалы ХЫУ международной конференции - Вологда, 2005 - С 62
27 Самоорганизация дислокаций в ультразвуковом поле / Г В Бушуева, Г М Зиненкова, Н А Тяпунина, В Т Дегтярев // Кристаллофизика 21-го века Материалы докладов Третьей международной конференции - Черноголовка, 2006 - С 123
28 Дегтярев В Т Изменение пластических свойств материалов в ультразвуковом поле // Наукоемкие технологии. - 2006 - Т 6, № 9 - С 51-52
29 Дегтярев В Т. Моделирование процессов прохождения скользящих дислокаций через ансамбли колеблющихся дислокаций // Наукоемкие технологии -2007 -Т6, № 1 -С 34-37.
30 Дегтярев ВТО возможных механизмах акустопластического эффекта//Доклады РАН -2007 -Т 414, №1 -С28-29
Дегтярев Вячеслав Тихонович
Динамика дислокаций в щелочно-галоидных кристаллах при наложении ультразвука
Автореферат диссертации на соискание ученой степени доктора физико-математических наук
Подписано в печать 11 07.2007г Формат бумаги 60x84 1/16 Бумага типографская № 2. Печать офсетная Уел печ л 2 0 Уч -изд л. 2.0 Тираж 100 экз Заказ № 148
Калужский филиал ГОУ ВПО «Московский государственный технический университет имени Н Э Баумана» 248600, г. Калуга, ул Баженова, 2
Отпечатано в Редакционно-издательском отделе
КФМГТУим Н.Э Баумана 248600, г Калуга, ул Баженова, 2, тел 57-31-87
ВВЕДЕНИЕ.
ГЛАВА 1. ИЗМЕНЕНИЕ МЕХАНИЧЕСКИХ СВОЙСТВ И ДИСЛОКАЦИОННОЙ СТРУКТУРЫ КРИСТАЛЛОВ ПОД ВЛИЯНИЕМ ВЫСОКОЧАСТОТНОЙ ВИБРАЦИИ: СОСТОЯНИЕ ВОПРОСА И ПОСТАНОВКА ЗАДАЧИ.
1.1. Реальные кристаллы и пластическая деформация.
1.2. Изменение механических свойств кристаллов под влиянием высокочастотной вибрации.
1.3. Поведение дислокаций в поле осциллирующего напряжения.
1.4. Влияние ультразвукового поля на дислокационную структуру кристаллов.
1.5. Данные об изменении дислокационной структуры под влиянием . ультразвука.
1.6. Акустопластический эффект.
1.7. Некоторые подходы к моделированию процессов зарождения и прохождения дислокаций.
ГЛАВА 2. МОДЕЛИ, АЛГОРИТМ И ПАКЕТ ПРОГРАММ.
2.1. Модели, использованные для анализа акустопластического эффекта.
2.2. Моделирование процесса прохождения краевой дислокацией модельной площадки в квазистатическом приближении
2.3. Моделирование процесса прохождения краевой дислокацией модельной площадки с использованием уравнения движения.
2.4. Модель и алгоритм для исследования поведения ансамбля дислокаций в ультразвуковом поле.
ГЛАВА 3. ОСНОВНЫЕ ЭЛЕМЕНТЫ ДИСЛОКАЦИОННОЙ СТРУКТУРЫ И ОСОБЕННОСТИ ИХ ДВИЖЕНИЯ В УЛЬТРАЗВУКОВОМ
ПОЛЕ.
3 Л. Поведение единичной дислокации в ультразвуковом поле.
3.2. Поведение дислокационного диполя в ультразвуковом поле.
3.3. Дислокационные сегменты в ультразвуковом поле.
ГЛАВА 4. САМООРГАНИЗАЦИЯ ДИСЛОКАЦИЙ В УЛЬТРАЗВУКОВОМ ПОЛЕ.
4.1. Процесс формирования упорядоченных дислокационных структур в ультразвуковом поле.
4.2. Характерное время формирования ячеистой структуры в ультразвуковом поле.
ГЛАВА 5. ПРОХОЖДЕНИЕ СКОЛЬЗЯЩЕЙ ДИСЛОКАЦИИ ЧЕРЕЗ *
АНСАМБЛЬ КОЛЕБЛЮЩИХСЯ ЛЕСНЫХ ДИСЛОКАЦИЙ.
5.1. Моделирование в квазистатическом приближении.
5.2. Моделирование процесса прохождения в квазидинамическом приближении.
К числу важных проблем физики твердого тела относится создание новых материалов с заданными свойствами. Столь же важным является и вопрос о том, как могут изменяться свойства конструкционных материалов и твердотельных элементов приборов в условиях их эксплуатации. Яркими примерами могут служить проблемы пластичности при обработке трудно-деформируемых материалов, радиационных повреждений, наблюдаемых в ряде элементов ядерных реакторов, и изменения прочностных свойств конструкционных материалов, происходящих при возникновении вибраций.
В этих областях широкое техническое применение получили реальные кристаллы, характеризуемые наличием в них разного рода дефектов - поликристаллы, свойства которых контролируются как межзеренными, так и внутризеренными процессами. Какие из этих процессов станут определяющими, зависит от условий нагружения. Так, при высокочастотных вибрациях основной вклад будут вносить процессы, происходящие внутри зерен. Зерна поликристаллов (кристаллиты), являются, как правило, монокристаллами, размеры которых могут быть существенно различными: от долей микрометра (нанокристаллы) - до сантиметров (бикристаллы).
В трех вышеприведенных проблемах особый интерес представляют исследования, связанные с изменением механических свойств монокристаллов («реальных кристаллов») под действием высокочастотных вибраций ультразвукового диапазона, и прежде всего, изучение микропроцессов, происходящих в объеме образца за времена порядка периода ультразвука 10"5 с, и обусловленных ими изменений свойств материалов.
Основа таких изысканий заключается в исследовании дефектов, определяющих структурно-чувствительные свойства твердых тел, например, дислокаций. Натурные методы изучения поведения дислокаций (избирательное травление и электронная микроскопия) позволяют получить сведения лишь о том, что было до начала ультразвукового воздействия (в исходном состоянии образца) и после его окончания (в финальном состоянии образца).
Результаты электронной микроскопии, полученные на тонких пленках, даже если наблюдения ведутся в процессе ультразвукового воздействия, для описания изменений свойств массивных образцов использовать не корректно, так как свойства тонких пленок контролируются поверхностными эффектами.
Наиболее информативным, и, по-видимому, в настоящее время реально возможным, является метод компьютерного моделирования, поскольку экспериментально выявить вклады различных типов дислокационных ансамблей и факторы, влияющие на движение дислокаций в их связи с изменением пластичности и упрочнения твердых тел, затруднительно.
В качестве объекта исследования целесообразно было взять широко используемые в качестве модельных материалов кристаллы каменной соли №С1. Это разумно сделать по следующим соображениям: №С1 имеет кубическую структуру, структурные дефекты в этих кристаллах хорошо изучены, имеются экспериментальные данные об эволюции дислокационной структуры под действием ультразвука и об акустопластическом эффекте, который заключается в изменении пластических свойств материалов в присутствии ультразвукового поля.
Такие исследования уже в течение значительного времени ведутся зарубежными и отечественными авторами. При моделировании различных процессов пластической деформации они рекомендовали рассматривать процесс эволюции дислокационной структуры при изменяющихся внешних условиях в квазистатическом приближении. По их мнению, при рассмотрении таких процессов необходимо обращать внимание на гибкость дислокаций, силы сопротивления решетки движению дислокаций и на тонкую структуру поля внутренних напряжений, в котором происходит перемещение дислокаций. Показано при моделировании прохождения скользящей дислокации через плоские дислокационные скопления, что напряжение прохождения через них меняется в широких пределах в зависимости от структуры скопления.
Предлагалось при моделировании дислокационные ансамбли разделять на два вида: когда составляющие поля напряжения дислокационных ансамблей на плоскости скольжения гибкой дислокации имеют сингулярные точки и когда поле напряжений ансамбля таких ограничений не имеет. Ансамбль дислокаций, пересекающих плоскость скольжения исследуемых дислокаций, (ансамбль дислокаций леса) и является примером ансамблей первого типа и представляется более интересным случаем. Здесь сингулярные точки отвечают точкам выхода дислокаций леса на плоскость скольжения пробной дислокации.
Наличие особых точек затрудняет поиск решения в окрестности этих сингулярных точек, однако решение этой задачи обогатит новыми данными характеристики процессов движения дислокаций для дальнейшей разработки физической теории пластичности и изыскания возможных способов целенаправленного изменения механических свойств кристаллических материалов. С другой стороны, существование сингулярных точек на плоскости скольжения, например гибкой дислокации, в значительной мере облегчает поиск решения искомой задачи.
В целом же моделирование дает возможность выяснить роль отдельных факторов, приводящих к изменению макроскопических свойств материалов, и механизмов, обуславливающих эти изменения. Поэтому разработка моделей, алгоритмов и программ для осуществления моделирования дислокационных процессов по-прежнему является актуальной.
Посчитали целесообразным изучить теоретические и экспериментальные данные по понижению предела текучести, упрочнению структуры хорошо изученных кристаллов, например ИаО, а также известные модели по изучению условий возникновения, взаимодействия и движения дислокаций в кристаллах при наложении ультразвука. Затем сформулировать цели и задачи работы. После такого анализа цели настоящей работы выглядят следующим образом:
- детальное исследование процессов, происходящих при воздействии высокочастотных колебаний ультразвукового диапазона на дислокационные структуры и пластичность материалов;
- определение режимов нагружения, позволяющих добиться необходимого изменения пластических свойств образца, будь то упрочнение или разупрочнение.
Объекты исследования: щелочно-галоидные кристаллы со структурой хлорида натрия.
При решении поставленной задачи были применены два подхода: квазистатический и динамический. В первом, с использованием уравнения равновесия, при заданных частоте и амплитуде ультразвука находили равновесные дислокационные структуры, во втором - анализировали закономерности движения дислокаций в ультразвуковом поле с использованием уравнения динамики.
Для достижения этих целей предстояло решить следующие задачи:
1. Разработать физическую модель и методику моделирования механизмов и процессов, обуславливающих акустопластический эффект, используя строгий динамический подход, основывающийся на решении уравнения движения с учетом поля сил взаимодействия дислокаций и самодействия, а также гибкости скользящих дислокаций и влияния полей, обусловленных внешним нагружением.
2. Промоделировать процессы движения и размножения дислокаций в условиях сложнонагруженного состояния кристалла, когда на дислокации леса действует знакопеременная нагрузка, а в плоскости скользящей дислокации действуют или постоянная сила, или сумма постоянной и знакопеременной составляющих поля.
3. На основании анализа полученных результатов моделирования выявить зависимости критического напряжения, необходимого для начала пластической деформации, от плотности лесных дислокаций и от параметров ультразвука, то есть установить связь между дислокационными процессами и акустопластическим эффектом.
Методы исследования, реализованные для достижения цели работы, основаны на развитии существующих и создании новых, адекватных рассматриваемым физическим явлениям и процессам, математических моделей исследования дислокационных процессов в условиях высокочастотной вибрации.
С целью оптимизации величин напряжения, необходимого для преодоления пробной дислокацией модельной площадки, были разработаны квазистатическая и динамическая модели, с целью исследования особенностей перераспределения дислокаций леса в ультразвуковом поле использована динамическая модель.
Достоверность полученных результатов обеспечена детальным рассмотрением физических явлений и процессов, связанных с воздействием ультразвука на кристаллы с дефектами, сравнением результатов моделирования с данными других авторов, что позволило создать математические модели, реализующие новые подходы к количественному и качественному описанию дислокационных процессов, в присутствии ультразвукового поля и в значительной степени тем, что основные результаты моделирования совпадают с результатами натурных экспериментов, выполненных другими авторами.
Научная новизна данной диссертационной работы заключается в том, что впервые предпринята попытка методом компьютерного моделирования исследовать влияние ультразвукового поля на предел текучести, характеризующий величину акустопластического эффекта, в результате чего раскрыты различные механизмы его возникновения и оценены их вклады. Впервые предложена математическая модель, описывающая работу источника Франка-Рида, концы которого закреплены лесными дислокациями, совершающими гармонические колебания, и установлены закономерности поведения источника в этой ситуации.
Впервые найдены условия, при которых, при наличии знакопеременной нагрузки, возможно накопление дислокаций в плоскости источника.
Предложена модель, позволяющая описывать прохождение пробной дислокацией модельной площадки, в которой движение дислокации описывается в приближении динамической теории.
Практическая ценность работы состоит в том, что полученные в ней результаты позволяют предвидеть изменение пластических свойств материалов, подвергающихся высокочастотной вибрации в процессе их эксплуатации, что важно для оценки срока службы и предвидения возможности деградации материалов. Кроме того, в современных технологиях обработки материалов все более широкое применение находит ультразвук, причем как в процессе обработки, так и для придания материалам наперед заданных свойств, например, чтобы достичь высокой пластичности материала непосредственно во время технологического воздействия.
Основание для проееденш работ.
План совместных НИР кафедры «Молекулярная физика» физического факультета МГУ имени М.В. Ломоносова, кафедры «Компьютерные системы и сети» и кафедры «Системы автоматизированного проектирования» Калужского филиала МГТУ имени Н.Э. Баумана.
Научные положения и научные результаты, выносимые на защиту:
1. Установленные впервые закономерности генерации дислокаций в условиях сложнонагруженного состояния, когда точки закрепления источника дислокациями леса совершают вынужденные колебания.
2. Обнаруженное влияние соотношения фаз колебаний закрепляющих дислокаций и источника Франка-Рида, в условиях наложения ультразвукового поля, и снижение эффекта упрочнения в случае колеблющегося дислокационного леса по сравнению с неподвижным лесом дислокаций.
3. Эффект пластифицирования материала при воздействии на лес дислокаций ультразвукового поля и усиление этого эффекта с увеличением плотности дислокаций леса.
4. Явление самоорганизации системы дислокаций леса в ультразвуковом поле, где самоорганизация включает в себя образование динамических дислокационных структур с диполями, мультиполями и стенками дислокаций в качестве компонентов самоорганизованной структуры, причем конечным результатом самоорганизации может стать формирование блочной структуры.
5. Обнаружение устойчивости границ блоков, образовавшихся под действием ультразвука, т.е. сохранения этих границ после выключения ультразвукового поля. При этом структура границ ультразвукового происхождения отличается от структуры границ термического происхождения (образовавшихся в процессе отжига).
6. Разработанные модель, алгоритм и пакет программ для анализа поведения элементов дислокационной структуры в ультразвуковом поле.
7. Результаты моделирования поведения элементов дислокационной структуры в условиях одновременного воздействия на кристалл постоянной и знакопеременной составляющих внешней нагрузки.
Личный вклад автора в проведенное исследование.
Автором проведен детальный анализ экспериментальных данных, известных из литературы к началу исследования, для щелочно-галоидных кристаллов, на основании которого предложены модели, адекватно описывающие поведение элементов дислокационной структуры в ультразвуковом поле.
Разработаны алгоритмы и пакет программ, позволяющих моделировать взаимодействие элементов дислокационных структур с ультразвуковым полем и эволюцию дислокационной структуры кристалла под воздействием ультразвука.
Автором на кристаллах со структурой КаС1, для которых имеются надежные экспериментальные данные, выполнено компьютерное моделирование поведения элементов дислокационной структуры и ансамблей дислокаций в ультразвуковом поле. Проведен анализ полученных компьютерным моделированием результатов и установлены закономерности эволюции дислокационных ансамблей в ультразвуковом поле.
Апробация. Результаты диссертационной работы докладывались на 6 Всероссийских и международных научных конференциях и научнопрактических семинарах, в т.ч. на X международной конференции «Imperfection interaction and anelasticity phenomena in solids (IIAPS - 10)» (Тула, 2001); Всероссийской конференции, посвященной 100-летию со дня рождения академика Г.В. Курдюмова «Дефекты структуры и прочность кристаллов» (Черноголовка, 2002); XLII международной конференции «Актуальные проблемы прочности» (Калуга, 2004); XLIV международной конференции «Актуальные проблемы прочности» (Вологда, 2005); третьей международной конференции «Кристаллофизика 21-го века» (Черноголовка, 2006).
Публикация результатов работы. Материалы исследований, представленных в диссертации, изложены в 30 работах, опубликованной в научных журналах, научно-технических сборниках, материалах, трудах и тезисах докладов Всесоюзных, Российских и международных конференций, симпозиумов, совещаний и семинаров; ссылки на 30 из них приведены в тексте диссертации. Все основные результаты диссертационной работы опубликованы в 18 статьях, из них 72% (13 статей) опубликовано в ведущих российских журналах, в которых должны быть опубликованы основные научные результаты диссертации на соискание ученой степени доктора наук.
Структура и объем работы. В соответствии с поставленными целями исследования, характером и объемом проведенной работы, диссертация содержит введение, 5 глав, заключение, список литературы. Общий объем диссертации составляет 218 страниц машинописного текста, содержащих текст работы, 78 рисунков, 5 таблиц и список использованных источников на 17 страницах, содержащий 142 наименования.
ЗАКЛЮЧЕНИЕ И ВЫВОДЫ
В данной работе методом математического моделирования исследованы механизмы влияния ультразвука на эволюцию дислокационной структуры и вызванные ею изменения физических свойств материалов. Рассмотрено влияние ультразвука на процесс взаимодействия дислокаций различных систем. Впервые дислокационные процессы в ультразвуковом поле исследованы с использованием уравнений динамики дислокаций. Приводится сравнение результатов моделирования, полученных с использованием квазистатической и динамической моделей. Применение динамического подхода позволило выявить новые эффекты и получить более точные количественные оценки для процессов, ранее рассматривавшихся в квазистатическом приближении. Так, например, впервые рассмотрена генерация дислокаций источниками Франка-Рида, концы которых совершают вынужденные колебания.
Основное внимание уделено двум механизмам, ответственным за изменение пластических свойств материалов в ультразвуковом поле, а именно: генерации дислокаций источником Франка-Рида и прохождению дислокации через колеблющийся лес дислокаций. Существенно, что при этом учитывается возможность влияния ультразвука и на скользящую дислокацию, и на дислокации леса.
На некоторых моментах выполненной работы необходимо акцентировать внимание. Это прежде всего касается генерации дислокаций источником Франка-Рида.
Если центрами закрепления источника служат лесные дислокации разного знака, то под действием ультразвука концы сегмента будут двигаться в противофазе - в один полупериод действия ультразвука эффективная длина источника уменьшается, во второй - она увеличивается. Физически это можно трактовать как увеличение эффективной длины источника 1эф, а это приводит к тому, что генерация дислокаций может начаться при напряжениях въ въ
7 =-, меньших <7 = —, характерных для состояния источника с непод
Кф к вижными точками закрепления. Т.е. колебание лесных дислокаций приводит к понижению напряжения, необходимого для срабатывания источника. Следовательно, наблюдается эффект пластификации образца в ультразвуковом поле. По результатам моделирования с ростом амплитуды ультразвука эффект пластификации для источников длиной менее 4 мкм в значительной мере возрастает, а для источников больших 4 мкм такой эффект выражен в меньшей степени. Понижение напряжения срабатывания источника Франка-Рида для коротких источников по результатам моделирования может достигать 70%, в отличие от источников длиной более 4 мкм. Хотя в работе четко установлено, что для источников длиной 5 мкм при амплитуде ультразвука 0,4 МПа он еще существует - порядка 20%.
Анализ результатов моделирования действия ультразвука на лесные дислокации показывает, что наблюдается и в этом случае увеличение интервала длин источников, срабатывающих при данном постоянном напряжении: при постоянном напряжении в плоскости источника, равном 4 МПа, срабатывают источники Франка-Рида с длиной, превышающей 2 мкм, а при воздействии на лесные дислокации ультразвука амплитудой 0,1 МПа начинают срабатывать источники с длинами от 1,5 мкм.
Использование этой модели при рассмотрении колебания лесных дислокаций, являющихся центрами закрепления концов источника Франка-Рида, эффект изменения напряжения начала работы источника не обнаружен.
Следовательно, понижение напряжения срабатывания источника Франка-Рида будет обнаруживать зависимость от структуры леса, так как вероятность появления источника, концы которого колеблются в противофазе, зависит от соотношения между положительными и отрицательными дислокациями в ансамбле леса. Следовательно, этот эффект будет зависеть от соотношения положительно и отрицательно ориентированных дислокаций, а не только от плотности дислокаций леса. В этом и состоит отличие данной ситуации от эффекта упрочнения, обусловленного неподвижными дислокациями леса.
Когда на сегмент действует как знакопеременная, так и постоянная нагрузка (в случае сложнонагруженного состояния образца), а на дислокации леса - только знакопеременная, на напряжение срабатывания источника будут влиять оба механизма.
В реальных кристаллах, в которых наиболее вероятен вариант комплексного нагружения, вероятность срабатывания источника Франка-Рида будет зависеть от соотношения отрицательно и положительно ориентированных лесных дислокаций и от того, как в текущем сегменте колеблются закрепляющие дислокации относительно самого сегмента.
Результаты выполненных исследований, в отличие от известных подходов [130-142] к изучению процессов, наблюдаемых в кристаллах, показали, что наряду с образованием динамических дислокационных диполей при наложении ультразвука формируются и более сложные упорядоченные дислокационные ансамбли - мультиполи и полигональные границы блоков, причем «затравками» в этом случае являются именно дислокационные диполи.
Причем число и структура образующихся полигональных границ зависят от частоты, амплитуды и времени действия ультразвука, а также от начальной плотности и распределения дислокаций, в том числе от соотношения числа дислокаций разного знака в дислокационном ансамбле. Это подтверждено при моделировании поведения в ультразвуковом поле ансамбля дислокаций одной системы скольжения.
Поскольку при анализе движения скользящей дислокации следует учитывать и поле ансамбля лесных дислокаций, впервые было произведено моделирование и анализ поведения леса краевых дислокаций в ультразвуковом поле. Основным условием моделирования являлось соответствие количества дислокаций на модельной площадке реальным плотностям дислокаций в кристалле.
Моделирование проводилось с ансамблем дислокаций от 200 до 300 штук, что соответствовало плотности дислокаций р = 8-106 см-2. Рассматривались ансамбли, состоящие из прямолинейных бесконечных лесных дислокаций, и ансамбли, состоящие из дислокаций одного знака и из дислокаций разного знака.
Интересные результаты были получены в процессе моделирования, когда непосредственно после начальной расстановки дислокаций на модельной площадке эту систему дислокаций «подвергли» действию ультразвука, то получили ячеистую структуру, границы блоков которой приобрели тонкую структуру. Подвергнув релаксации полученную в результате воздействия ультразвука эту дислокационную структуру, наблюдали в сравнении, что характер расположения ячеек сохранился, а структура границ блоков изменилась. В процессе релаксации дислокационной структуры, созданной ультразвуком, изменилось расположение дислокаций, принадлежащих границам ячеек. К числу наблюдаемых характерных элементов тонкой структуры границ блоков следует отнести пары дислокаций - дислокационные диполи.
Предложенное моделирование процессов показало, что в неупорядоченном ансамбле дислокаций в ультразвуковом поле происходят процессы самоорганизации дислокаций и формирования ячеистой структуры. Тонкая структура границ ячеек оказывается различной, если процесс полигонизации происходил в результате релаксации в ультразвуковом поле (или в отсутствие ультразвука).
Картина, полученная в результате моделирования, оказалась полностью аналогичной и совпала с экспериментальными результатами, полученными Г.В. Бушуевой с сотрудниками при исследовании образцов СсШ методами электронной микроскопии. Это позволяет заключить, что предложенные методы компьютерного моделирования, программы и алгоритмы на сегодня оптимальны и позволяют с достаточной достоверностью проследить за процессами влияния ультразвука на характер размножения и движения дислокаций и по режимам воздействия ультразвука можно наблюдать как уменьшение предела текучести, так и его увеличение. Стали понятными многие экспериментальные результаты, проанализированные в гл .1.
На основании полученных результатов моделирования и их анализа могут быть сделаны следующие выводы:
1. Ультразвуковое воздействие на неупорядоченные дислокационные ансамбли приводит к перераспределению дислокаций, в результате которого формируются упорядоченные динамические дислокационные структуры - диполи, мультиполи и дислокационные стенки. В совокупности наблюдавшиеся элементарные процессы интерпретируются как самоорганизация дислокационной подсистемы (полигонизация, образование ячеистой структуры).
2. Самоорганизация дислокационной подсистемы приводит к достижению состояния, в котором суммарная доля дислокаций, вошедших в упорядоченные дислокационные структуры, не изменяется. Время достижения этого состояния имеет порядок 102 периодов и зависит и от частоты, и от амплитуды ультразвука.
3. Впервые показано, что в ультразвуковом поле формируется система субграниц, каждая из которых образована дислокациями одного знака, но при этом любые две соседние субграницы образованы дислокациями другого знака.
4. В границах блоков, сформировавшихся в ультразвуковом поле, дислокации не лежат в одной плоскости, и границы имеют тонкую структуру, т.е. содержат дислокационные диполи и мультиполи. Этим они отличаются от дислокационных стенок, формирующихся при отжиге в отсутствие ультразвукового поля.
5. Наложение ультразвукового поля запускает процесс размножения дислокаций по двум каналам. Одним из каналов является активация источников Франка-Рида, но при этом благодаря действию ультразвука активация каждого данного источника происходит при меньшем значении постоянного напряжения сдвига (эффект ультразвука сводится к действию спускового ме
201 ханизма). Вторым каналом является увеличение числа источников Франка-Рида.
6. Прохождение скользящей дислокации сквозь дислокации леса в условиях наложения ультразвукового поля может давать в результате либо упрочнение, либо пластифицирование материала. Знак эффекта определяется соотношением фаз колебаний скользящей дислокации и дислокаций леса.
7. В случае колебания дислокаций леса по гармоническому закону, напряжение, необходимое для продвижения скользящей дислокации через единичную площадку, понижается и может составить лишь 60% от величины напряжения, необходимого для прохождения дислокации сквозь лес неподвижных дислокаций.
8. Совокупность данных, полученных в результате проведения компьютерных экспериментов, показывает, что при динамическом подходе удается учитывать временные зависимости (эффекты вязкости) и частоту ультразвука. Это позволит использовать эти результаты на предприятиях, применяющих ультразвук при получении лент и проволок из непластичных материалов, например, ОАО «Биметалл» и ЗАО «Аметист-Стан». При квазистатическом подходе указанные эффекты обеспечить невозможно.
1. Физика кристаллов с дефектами / A.A. Предводителев, H.A. Тяпунина, Г.М. Зиненкова, Г.В. Бушуева. М.: Изд-во МГУ им. М.В. Ломоносова, 1986.-240 с.
2. Смирнов B.C. Теория обработки металлов давлением. М.: Металлургия, 1978.-496 с.
3. Теория пластических деформаций металлов / Е.П. Унксов, У. Джонс, В.Л. Колмогоров и др. / Под ред. Е.П. Унксова, Г.В. Овчинникова. М.: Машиностроение, 1983. - 598 с.
4. Ивлев Д.Д. Теория идеальной пластичности. М.: Наука, 1966. - 231 с.
5. Соколовский В.В. Теория пластичности. М.: Высшая школа, 1969. -608 с.
6. Сегал В.М. Технологические задачи теории пластичности (методы исследования). Минск: Наука и техника, 1977. - 257 с.
7. Ильюшин A.A. Пластичность. Основы общей математической теории.- М.: Изд-во АН СССР, 1963.-271 с.
8. Ивлев Д.Д., Быковцев Г.И. Теория упрочняющегося пластического тела. -М.: Наука, 1971.-231 с.
9. Колмогоров В.Л. Напряжение, деформация, разрушение. М.: Металлургия, 1970.-229 с.
10. Ю.Работнов Ю.Н. Ползучесть элементов конструкций. М.: Наука, 1966.- 752 с.1 l.Schmid E.G. Plasticity of insonated metals // Japan Jnst. Metals. 1968. -№9.-P. 797.
11. Blacha F., Langenecker B. Plastizitats Intersuchungen von Metallristallen in Ultrashalfeld // Acta Met. 1960. - Bd. 5. - S. 425-449.
12. Кацнельсон A.A., Олемской А.И. Микроскопическая теория неоднородных структур. М.: Йзд-во Моск. ун-та, 1987. - 250 с.
13. Langenecker В. Work hardening of zinc crystals by high-amplitude ultrasonic waves // Proc. Amer. Soc. Test. Mater. 1962. - №82. - P. 602-609.
14. Балалаев Ю.Ф., Холодный В.И. Ультразвуковые разрушения жаропрочной стали Х12Н2-2ТЗМР // Вопросы физики твердого тела. 1969. -№1. - С. 206-210.
15. Балалаев Ю.Ф. Выявление дислокаций ультразвуком // Заводская лаборатория. 1969. - Т.27, № 11. - С. 1362.
16. Классен-Неклюдова М.В., Капустин А.П. Влияние ультразвука на поля упругих напряжений в кристаллах // ДАН СССР. 1951. - №77. -С.1019.
17. Дислокационная структура никеля после ультразвукового нагружения с низкой и ультразвуковой частотами / И.А. Гиндин, И.М. Неклюдов, М.П. Старолат и др. // Физика твердого тела. 1970. - №8. - С. 24562458.
18. Пластическая деформация алюминия под действием ультразвуковых импульсов / И.А. Гиндин, И.М. Неклюдов, Г.Н. Малик, О.И. Волчок // Физика твердого тела. 1969. - Т. 11, №11. - С.3236-3241.
19. Гарбер Р.И., Солошенко И.И., Халдей O.A. Размножение и взаимодействие дислокаций при многократном нагружении // Динамика дислокаций. Харьков: ФТИНТ АН УССР, 1968. - 873 с.
20. Гарбер Р.И., Солошенко И.И., Халдей O.A. Понижение критических напряжений размножения дислокаций при многократном нагружении. // Межвузовская научная конференция по проблемам прочности и пластичности. 1967. - С. 51-52.
21. Полоцкий И.Г., Базелюк Г.Я., Ковш C.B. Действие ультразвука на дислокационную структуру монокристаллов молибдена // ФТТ. 1969. -Т.11, №3. - С.755.
22. Коновалов Е.Г., Ремизовский Э.И. Металловедение и термическая обработка металлов. Минск: Наука и техника, 1965. - 190 с.
23. Пинес Б.Я., Омельяненко И.Ф., Сиренко А.Ф. Влияние ультразвуковых колебаний на кинетику гетеродиффузии в образцах Fe-Al, Ni-Cu и Ni-Si // Физика металлов и металловедение. 1969. - Т.27, №6. - С. 157161.
24. Платков В.Я., Ефименко В.П., Старцев В.И. Изучение некоторых динамических характеристик дислокаций в кристаллах бористого калия методом внутреннего трения // Физика твердого тела. 1968. - Т.9, №10.-С. 2799-2803.
25. Криштал М.А., Головин С.А. Внутреннее трение упорядоченной стали // Физика металлов и металловедение. 1962. - Т. 14, №6. - С. 913-916.
26. Тяпунина H.A., Наими Е.К., Зиненкова Г.М. Действие ультразвука на кристаллы с дефектами. М.: МГУ, 1999. - 215 с.
27. Косевич A.M. Динамическая теория дислокаций // УФН. 1964. -№84.-С. 579-609.
28. Дубнова Г.Н., Инденбом В.Л., Штольберг A.A. О прогибании дислокационного сегмента и источники Франка-Рида // Физика твердого тела. -1968. Т. 10, №6. - С.1760-1768.
29. Дубнова Г.Н., Штольберг А.А. Строгое решение задачи о прогибании сегмента // Динамика дислокаций. Харьков: ФТИНТ АН УССР, 1968.- 873 с.
30. Kochler J.S. Imperfections in perfect crystals // Appl. Phys. 1952. -№ 157.-P. 157-160.
31. Granato A.V., Lucke K., Stern R.M. La resonance des dislocations et son application aux mesures des interactions entre defaults ponctuels et dislocations //Metaux. 1961. -№433. - P. 300-319.
32. Granato A.V., Lucke K., Teutonico L. Jump frequency of a pinned dislocation segment // Bull. Amer. Phys. Soc. 1962. - №3. - P. 223.
33. Brailsford A.D. Abrupt-kink model of dislocation motion // Phys. Rev. -1962. -№3. -P. 1033-1038.
34. Алефельд Г. Сопоставление дислокационных моделей струны и цепочки перегибов // Актуальные вопросы теории дислокаций. М.: Мир, 1968.-227 с.
35. Bode Е. Durch Versetzungen bedingte innere Reibung in Silber // Z. Naturforsch. 1959. - Bd. 14a, №8. - S.762-763.
36. Турков C.K., Шермегор Т.Д. О внутреннем трении при взаимодействии примесных атомов с краевыми дислокациями // Физика твердого тела.- 1964. Т.6, №12. - С.3502-3508.
37. Bauer C.L., Gordon R.B. Mechanism for dislocation pining in the alkali halides // J.Appl.Phys. 1962. - V.33, № 2. - P. 672-682.
38. Whitworth R.W. Some effects of vibration on the internal friction of sodium chloride // Phil. Mag. 1960. - V.5, №53. - P. 425-440.
39. Швидковский Е.Г., Тяпунина Н.А, Белозерова Э.П. Рождение дислокаций при вибрации кристаллов фтористого лития и хлористого натрия // Кристаллография. 1962. - Т.7, № 3. - С. 473-474.
40. Швидковский Е.Г., Тяпунина H.A., Белозерова Э.П. О зарождении дислокаций в кристаллах фтористого лития под влиянием высокочастотной вибрации // Кристаллография. 1966. - Т.11, №4. - С. 651-655.
41. Пинес Б.Я., Омельяненко И.Ф. Размножение дислокаций под действием ультразвуковых колебаний в кристаллических образцах Си, Ni, AI и в монокристаллах LiF и NaCl // Динамика дислокаций. Харьков: ФТИНТ АН УССР, 1968. - 873 с.
42. Швидковский Е.Г., Тяпунина H.A., Белозерова Э.П. Влияние электрического поля на поведение заряженных дислокаций // Кристаллография. 1962. - Т.7, №3. - С. 471-472.
43. Криштал М.А., Выбойщик М.А. Взаимодействия между дислокациями и атомами примесей в металлах и сплавах. Тула: Изд-во ТПИ, 1969. -110 с.
44. Нацик В.Д., Чишко К.А. Физика конденсированного состояния. -Харьков: ФТИНТ АН УССР, 1974. Вып. 33. - С. 44.
45. Кульган И.Н., Тяпунина H.A. Поведение дислокационных петель в ультразвуковом поле // ЖТФ. 1994. - Т.64, №2. - С. 105-113.
46. Blaha F., Langenecker В., Ollshlagel D. Zum plastischen Verhalten von Metallen under Schalleiwirkung // Metallkunde. 1960. - Bd. 51, №11. -S.636-638.
47. Влияние ультразвука на степень совершенства монокристаллов алюминия, выращенных из расплава / И.Г. Полоцкий, Д.Е. Овсиенко, З.А. Ходов и др. // Физика металлов и металловедение. 1966. - Т.21, №5.-С. 727-731.
48. Klim D., Paufler P. Propagation of Cracks in CaP as Rerealed by Measurement of Ultrasonic Damping // Crist. Res. Technol. 1988. - V.23, №8. -P. 121.
49. Благовещенский B.B., Тяпунина H.A. Особенности работы источника Франка-Рида под действием ультразвука // ДАН СССР. 1980. - Т.254, №4. - С.869-872.
50. Никаноров С.П., Кардашев Б.Г. Упругость и дислокационная неупругость кристаллов. М.: Наука, 1985. - 360 с.
51. Лебедев А.Б., Кустов С.Б., Кардашев Б.К. Амплитуднозависимое поглощение ультразвука и акустопластический эффект в процессе активной деформации монокристаллов хлористого натрия // ФТТ. 1982. -Т.24.-С. 3169.
52. Дегтярев В.Т., Логинов Б.М. Исследование влияния тонкой структуры дальнодействующих полей напряжений ансамбля колеблющихся дислокаций на характеристики процесса прохождения скользящих дислокаций //ВИНИТИ. 1985. - №7043-В. -48с.
53. Дегтярев В.Т., Логинов Б.М., Тяпунина H.A. Исследование зависимости характеристик процесса движения скользящих дислокаций через ансамбль колеблющихся дислокаций от величины амплитуды колебаний // ВИНИТИ. 1986. -№2024-Д. -55с.
54. Дегтярев В.Т., Логинов Б.М. Моделирование движения дислокаций через колеблющийся дислокационный лес в кристаллах // Структура дислокаций и механические свойства металлов и сплавов: Материалы IV Всесоюзного семинара. Свердловск, 1987. - С.52-53.
55. Дегтярев В.Т., Логинов Б.М. Моделирование движения дислокаций через колеблющийся лес дислокаций с учетом дальнодействующих полей напряжений в кристаллах с гексагональной плотноупакованной решеткой//Физика металлов и металловедение. 1987. - С.608-610.
56. Дегтярев В.Т., Логинов Б.М., Тяпунина H.A. Моделирование движения дислокаций через лес колеблющихся дислокаций с учетом дальнодействующих полей напряжений в кристаллах с решеткой // Кристаллография. 1988. -Т.ЗЗ, №1. - С.163-166.
57. Дегтярев В.Т., Логинов Б.М., Еремеев A.B. Исследование влияния параметров площадки моделирования на характеристики процессов движения скользящих дислокаций через хаотические дислокационные ансамбли // ВИНИТИ. 1989. - № 2358-В89. -25с.
58. Дегтярев В.Т., Логинов Б.М. Исследование процессов прохождения скользящей дислокации через ансамбль хаотически распределенных призматических дислокационных петель, подверженных воздействию периодической нагрузки // ВИНИТИ. 1989. - №342-В89. - 20с.
59. Дегтярев В.Т., Логинов Б.М. Исследование процессов прохождения скользящей дислокации через ансамбль хаотически распределенных призматических дислокационных петель, подверженных воздействию периодической нагрузки // ВИНИТИ. 1989. - № 620-В89. - 29с.
60. Дегтярев B.T., Логинов Б.М., Ямпольский В.И. Анализ дефектной структуры кристаллов на основе моделирования процессов комплексного нагружения //165 лет МГТУ им. Н.Э. Баумана: Материалы НТК. -Калуга, 1995.-С.91.
61. Дегтярев В.Т., Лосев А.Ю., Плотников Ф.А. Перераспределение неупорядоченных дислокационных ансамблей в ультразвуковом поле // Наукоемкие технологии. 2005. - № 3-4, т. 6. - С. 5-8.
62. Дегтярев В.Т. Самоорганизация дислокационных ансамблей в ультразвуковом поле // Материалы электронной техники. Известия ВУЗов. -2004. -№1. С. 34-37.
63. Дегтярев В.Т., Лосев А.Ю. Динамические дислокационные структуры в ультразвуковом поле // Актуальные проблемы прочности: Сборник трудов XLII Международной конференции. Калуга, 2004. - С. 135.
64. Пластификация кристаллов ультразвуком, обусловленная взаимодействием дислокаций / В.Т. Дегтярев, А.Ю. Лосев, Ф.А. Плотников, H.A. Тяпунина // Труды ТулГУ. Сер. Физика. Тула, 2003. - Вып. 3. -С.3-8.
65. Полигонизация в ультразвуковом поле / В.Т. Дегтярев, А.Ю. Лосев, Ф.А. Плотников, H.A. Тяпунина // Известия РАН. Сер. Физическая. -2004.-Т. 68, №10.- С. 1516-1517.
66. Дегтярев В.Т., Лосев А.Ю., Плотников Ф.А. Динамические дислокационные структуры в ультразвуковом поле: диполи и Триполи // Материаловедение. 2004. - №7. - С. 8-12.
67. Дегтярев В.Т., Лосев А.Ю., Тяпунина H.A. Влияние ультразвука на процесс генерации дислокаций источником Франка-Рида // Деформация и разрушение материалов. 2005. - №3. - С 18-21.
68. Акустопластический эффект и факторы, на него влияющие / В.Т. Дегтярев, А.Ю. Лосев, Ф.А. Плотников, H.A. Тяпунина // Деформация и разрушение материалов. 2005. - №7. - С 23-27.
69. Дегтярев В.Т., Лосев А.Ю. Моделирование генерации дислокаций источником Франка-Рида в условиях акустопластического эффекта // Актуальные проблемы прочности: Материалы XLIV международной конференции. Вологда, 2005. - С. 61.
70. Дегтярев В.Т., Лосев А.Ю. Процессы, сопровождающие акустопласти-ческий эффект // Актуальные проблемы прочности: Материалы XLIV международной конференции. Вологда, 2005. - С.62.
71. Самоорганизация дислокаций в ультразвуковом поле / Г.В. Бушуева, Г.М. Зиненкова, H.A. Тяпунина, В.Т. Дегтярев // Кристалофизика 21-го века: Материалы докладов Третьей международной конференции. -Черноголовка, 2006. С. 123.
72. Дегтярев В.Т. Изменение пластических свойств материалов в ультразвуковом поле // Наукоемкие технологии. 2006. - № 9 - С. 51-52.
73. Дегтярев В.Т. Моделирование процессов прохождения скользящих дислокаций через ансамбли колеблющихся дислокаций // Наукоемкие технологии. 2007. - № 1 - С. 34-37.
74. Предводите л ев A.A., Ничуговский Г.И. Моделирование движения дислокаций через дислокационный лес // Кристаллография. 1972. - Т. 17, №1. - С. 166-171.
75. Предводителев A.A., Ничуговский Г.И., Веселов В.И. Моделирование движения дислокаций через дислокационный лес // Материаловедение. Воронеж: Изд. Воронежского политехнического института, 1975. -Вып. 2. - С.33-48.
76. Логинов Б.М., Предводителев A.A. Моделирование движения дислокаций через лес гибких дислокаций // Моделирование на ЭВМ дефектов в кристаллах. 1981. - Вып. 2. - С.53-55.
77. Логинов Б.М., Дегтярев В.Т., Тяпунина H.A. Моделирование скольжения дислокаций через дислокационный лес колеблющихся дислокаций в кристаллах с ГПУ структурой // Кристаллография. 1987. - Т.32, вып.4. - С.967-971.
78. Дегтярев В.Т. Моделирование процессов прохождения скользящих дислокаций через хаотические ансамбли колеблющихся дислокаций: Автореферат дис. . канд. физ.-мат. наук. -М.: МГУ, 1990. 16 с.
79. Благовещенский В.В. Особенности размножения дислокаций и образования полос скольжения под действием ультразвука: Автореф. дис. . канд. физ.-мат. наук. -М.: МГУ, 1982. 16 с.
80. Предводителев A.A. Подвижность, гибкость дислокаций и влияние этих факторов на их взаимодействие и прохождение через препятствия // Динамика дислокаций. Харьков: ФТИНТ АН УССР, 1968. - 873 с.
81. Бушуева Г.В. Взаимодействие дислокаций в кристаллах цинка и его особенности, выявляемые при электронно-микроскопическом исследовании тонких фольг: Автореф.дис. . канд. физ.-мат. наук. М.: МГУ, 1970.- 16 с.
82. ЮЗ.Стратан И.В., Предводителев A.A., Степанова В.М. Движение отдельных дислокаций в дислокационном ансамбле // ФТТ. 1970. -Т. 12, №3. - С.767.
83. Ничуговский Г.И. Моделирование процесса прохождения скользящих дислокаций через дислокационный лес и полосы скольжения: Дис. . канд. физ.-мат. наук. -М.: МГУ, 1976. 157с.
84. Foreman A.J.E., Makin M.J. Dislocation movement through random array of object// Philosophical magazine. 1966. - V.14, №131. -P.911.
85. Предводителев A.A., Логинов Б.М. Закономерности процесса прохождения дислокаций через гибкий и реагирующий лес дислокаций в области критической плотности дислокаций леса // ФТТ. 1986. - Т.28, №6. - С. 1896.
86. Predvoditelev A.A., Nichugovskii G.I., Veselov V.l. Simulation of dislocation Motion througth a dislocation forest // Physic Status Solid. 1981. -V.65. -P.469.
87. Логинов Б.М., Ничуговский Г.И., Предводителев A.A. Моделирование движения цуга дислокаций через дислокационный лес // Известия ВУЗов. Серия физическая. 1979. - №11. - С.97.
88. Логинов Б.М., Ничуговский Г.И., Предводителев A.A. Моделирование движения цуга дислокаций через дислокационный лес. Л.: ФТИ АН СССР, 1979.-142 с.
89. Логинов Б.М. Моделирование процесса прохождения скользящих дислокаций через дислокационный лес: Дис. . канд. физ.-мат. наук. -М.: МГУ. 1981.-304с.
90. Ш.Еремеев A.B. Моделирование процесса прохождения скользящих дислокаций через ансамбли дислокаций леса: Дис. канд. физ.-мат. наук. М.: МГУ, 1988. - 194с.
91. ПЗ.Нацик В.Д., Чишко К.А. Динамика и звуковое излучение дислокационного источника Франка-Рида. Харьков: ФТИНТ АН УССР, 1976. -7 с.
92. Кулемин A.B. Ультразвук и диффузия в металлах. М.: Металлургия, 1978.- 199 с.
93. Ландау А.И., Выдашенко В.Н. Термоактивированное движение дислокаций через хаотическую сетку точечных препятствий // Металлофизика. 1982. - Т.4, №4. - С.З.
94. Ландау А.И. Распределение углов огибания и длин дислокационных сегментов при статическом зависании дислокационной линии на сетке случайно расположенных локальных препятствий // Динамика дислокаций. Киев: Наукова думка, 1975. - 121 с.
95. Исследование взаимодействия базисных дислокаций с дислокационными кольцами в кристаллах цинка / Г.В. Бушуева, A.A. Предводителев, Р.Д. Фролова, В. Шарфф // ФММ. 1979. - Т.48, №3. -С. 647-653.
96. Предводителев A.A., Игонин С.И. Формирование полос скольжения при пластической деформации кристаллов // Элементарные процессы пластической деформации кристаллов. Киев: Наукова Думка, 1978. -190 с.
97. Предводителев A.A., Игонин С.И. Моделирование на ЭВМ процесса расширения полос скольжения // ФТТ. 1977. - Т.19, №9. - С. 17741777.
98. Игонин С.И., Предводителев A.A. Моделирование движения дислокаций вблизи полос скольжения // Вестн. МГУ. Физика. 1975. - №5. -С. 588-593.
99. Игонин С.И., Предводителев A.A. Особенности работы источников дислокаций в нерегулярном поле напряжений, создаваемом полосой скольжения // Вестн. МГУ. Физика. 1976. - №3. - С. 338-343.
100. Ничуговский Г.И. Моделирование процесса прохождения скользящих дислокаций через дислокационный лес и полосы скольжения: Дис. канд. физ.-мат. наук. -М.: МГУ,1976. 157с.
101. Веселов В.И., Ничуговский Г.И., Предводителев A.A. Моделирование процесса образования полосы скольжения // Известия ВУЗов. Физика.- 1981. Т.24, №9. - С. 82-86.
102. Акустические кристаллы: Справочник / A.A. Блистанов, B.C. Бондаренко, Н.В. Переломова и др. М.: Наука, 1982. - 235с.
103. Предводителев A.A. Подвижность, гибкость дислокаций и влияние этих факторов на их взаимодействие и прохождение через препятствие // Динамика дислокаций. Харьков: ФТИНТ АН УССР, 1968. - 873 с.
104. Дислокационные конфигурации, формирующиеся в кристаллах CdS под действием ультразвука / Г.В. Бушуева, Г.М. Зиненкова, В.И. Решетов, A.A. Хромов // Физика твердого тела. 1994. - Т. 36, №3.-С. 712-718.
105. Кансельсон М.И., Трефилов A.B. Динамика и термодинамика кристаллической решетки. М.: Изд. AT, 2002. - 310 с.
106. Чупрунов Е.Ф., Хохлов А.Ф., Фаддеев М.А. Кристаллография. М.: ФМЛ, 2000. - 285 с.
107. Никифоров К.Г. Физика твердого тела. Часть 1. Теория кристаллической решетки. Зонная теория твердого тела. Калуга: Изд-вао КГПУ им. К.Э. Циолковского, 2005. - 141 с.
108. Sarrak V.l., Mishin V.M., Spector A.Y. Delayed Fracture of High Strength Steels and its Criteria // Mat. Int. Conf. of Eng. Materials and Structures. -Singapore, 1991. P. 844-849.
109. Саррак В.И., Спектор А.Я., Мишин B.M. Сопротивление высокопрочной стали замедленному разрушению // Заводская лаборатория. 1991. -№ 2. - С. 57-58.
110. Мишин В.М. Фундаментальные основы методов диагностики замедленного разрушения высокопрочных сталей // Фундаментальные и прикладные проблемы приборостроения: Научн. труды V Международной науч.-практ. конф. М., 2002. - С.133-135.
111. Островский И.В., Стебленко Л.П., Надточный А.Б. Влияние ультразвуковой обработки на подвижность коротких дислокаций в кристаллах кремния // ФТТ. 2000. - Т. 42, вып. 3. - С. 478-481.
112. Красильников О.М., Доронин И.В., Трунин В.Ф. Математическое моделирование межатомного взаимодействия в жаропрочных сплавах на никелевой основе // Сборник научных трудов ЭПИ МИСИС. Электросталь, 2002 - С. 92-95.
113. Красильников О.М. Устойчивость кристаллической решетки щелочных металлов при отрицательных давлениях // Деформация и разрушение материалов: Сб. статей по материалам первой Международной