Математическое моделирование пластической деформации кристаллов

1.13. Визуализация

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

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

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

Число и тип связей, которые имеет атом, определяют локальную кристаллическую структуру. Например, атомы в совершенном ГЦК кристалле имеет 12 связей типа 421, тогда как в совершенном ГПУ кристалле имеют 6 связей типа 421 и 6 типа 422.

Использование CNA позволяет сделать видимыми при моделировании дислокации, границы зерен и дефекты упаковки. Например, при деформировании кристаллов меди, с помощью CNA классифицируют атомы на три класса: ГЦК, ГПУ и “другие”, где в класс “другие” попадают атомы, которые имеют число связей, отличное от 12, или тип их связей отличен от 421 и 422. Тогда внутренние дефекты упаковки видны как две сопряженные (111) плоскости ГПУ атомов, несвойственные дефекты упаковки видны как две (111) плоскости атомов ГПУ, разделенных (111) плоскостью атомов ГЦК. Границы двойников видны как одиночные (111) плоскости ГПУ атомов. Ядра дислокации и границы зерен состоят из атомов класса “другие”, хотя границы зерен содержат также немного ГПУ атомов.

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

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

2. Моделирования пластической деформации ГПУ-кристаллов

Автором была создана программа для изучения пластичности в двумерных кристаллах. Двумерные системы были выбраны, чтобы обойти проблемы, связанные с высокими требованиями к вычислительным ресурсам в случае трехмерных систем и визуализацией результатов вычислений. Для решения уравнений движения использовался алгоритм Верле со скоростью. В качестве потенциала взаимодействия между атомами был выбран потенциал Леннарда-Джонса. При моделировании вычислительная ячейка растягивалась вдоль умножением ее продольного размера на каждом шаге по времени на масштабный множитель (1+ε), где ε – малое число (0.00001), выбранное так, чтобы обеспечить требуемую скорость деформации. Координаты атомов при этом не менялись, т.е. при этом вводился зазор между атомами в смежных ячейках моделирования. При этом нагрузка прикладывалась к торцам кристалла, что лучше соответствует эксперименту. Поперечный размер системы контролировался с помощью динамического уравнения (12). Перед деформацией система приводилась в равновесное состояние с заданной температурой коротким прогоном с помощью МД. При вычислении кривых напряжение-деформация проводилось усреднение напряжения по атомам всей системы. С помощью закона сохранения энергии контролировалась правильность работы программы. Для наблюдения за временной эволюцией деформируемого кристалла использовалась анимация. Блок-схема программы приведена ниже.

 Скачать реферат