Размещено на http://www.allbest.ru/

Институт гидродинамики им. М.А. Лаврентьева СО РАН

Модель динамического деформирования кварца

В.В. Литвенко, Л.А. Мержиевский

Новосибирск, Россия

Основное содержание исследования

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

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

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

Рис.1 Приближенная фазовая диаграмма диокисида кремния

За основу при моделировании ударно-волновых процессов в диоксиде кремния взята модель вязкоупругого тела максвелловского типа [1]. Особенностью этой модели является тот факт, что она учитывает процесс релаксации касательных напряжений в материале при его деформировании. Такой подход не требует формулирования дополнительных феноменологических условий пластичности и позволяет единообразно описывать все состояния среды от упругого до гидродинамического. Соответствующая система уравнений состоит из законов сохранения массы, импульса, энергии и дифференциальных уравнений, описывающих эволюцию компонент тензора упругих деформаций. Для плоского одномерного случая система уравнений имеет вид:

где , . Здесь , - время и пространственная координата, , - начальная и текущая плотности вещества, - массовая скорость вещества, - удельная внутренняя энергия, - главные значения тензора деформация Генки, - главные значения тензора напряжений.

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

Зависимость для уравнения состояния строилась на основе широко известных традиционных уравнений Ми-Грюнайзена [2]. При этом предполагается, что вклад девиаторной составляющей в изменение энергии можно учесть с помощью дополнительного слагаемого . В этом случае уравнение состояния приобретает следующий вид:

,

где , - холодная и тепловая составляющие соответственно.

Входящие в уравнение слагаемые записываются в виде:

,

,

.

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

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

В работе построены все необходимые для замыкания модели зависимости. На рис.2 в координатах D,u (скорость ударной волны - массовая скорость) показана рассчитанная ударная адиабата кварца (сплошная линия) в сравнении с экспериментальными данными [3-9] (точки). Сравнение показывает хорошее соответствие расчетных и экспериментальных данных.

Рис.2 Сравнение расчетной ударной адиабаты с экспериментальными данными

Эти же данные в координатах плотность - давление показаны на рис.3. Другой характеристикой свойств веществ в условиях ударно-волнового нагружения являются адиабаты разгрузки из сжатого состояния. Немногочисленные экспериментальные данные об адиабатах разгрузки сравниваются с расчетом на рис.4. Определенное отличие в ходе адиабат разгрузки может объясняться как недостаточно аккуратным подбором использованных в расчете параметров замыкающих соотношений, так и неточностями в определении самих экспериментальных данных, связанными с особенностями использованных методик. В целом можно считать описание адиабат разгрузки приемлемым.

деформирование кварц диоксид кремний

Рис.3 Сравнение расчетной ударной адиабаты с экспериментальными данными

Рис.4 Сравнение расчетных адиабат разгрузки с экспериментальными

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

Работа выполнена при поддержке гранта РФФИ № 06-02-17335 и Интеграционного проекта СО РАН № 115.

Литература