CC BY
13
УДК 517 : (53+57/59)
А.И. Мельников*, Е.А. Митюшов*
* - Уральский государственный университет путей сообщения
* - Уральский государственнй технический университет (УПИ)
МОДЕЛЬНАЯ ЗАДАЧА ПОТЕРИ УСТОЙЧИВОСТИ КРИСТАЛЛА С ГПУ- СТРУКТУРОЙ ПРИ РАСТЯЖЕНИИ
Abstract
A model of a single crystal with hexagonal close packed structure has been performed which described its behaviour. According the theories of bifurcation and nonlinear analysis the equation of stability (instability) has been obtained. The grafical solution of the equation and the value of the critical force have been found.
Согласно экспериментальным данным [1], [2] разрушение идеального кристалла при нагружении происходит за счёт сдвиговой деформации, которая характеризуется образованием полос скольжения. Как известно, скольжение происходит в случаях, когда фактор Шмида отличен от нуля. Однако следует ожидать, что при определённых условиях сдвиговая деформация будет развиваться и в случае, когда фактор Шмида равен нулю, то есть когда направление растяжения перпендикулярно плоскости скольжения.
Теоретические исследования такой ситуации были проведены в работах [3], [4], [5]. Используя методы нелинейного анализа, изучалась потеря устойчивости двумерного кристалла при растяжении, атомная решётка которого соответствовала расположению атомов в плоскости (0001) металла, имеющего ГПУ-структуру. Применяя потенциал межатомного взаимодействия Леннарда-Джонса в работе [4] и самого общего вида в работе [5], получены критерии потери устойчивости атомного слоя и определена критическая нагрузка.
Предлагаемая работа является продолжением исследований, проведённых в [3], [5], в которой строится модель, характеризующая потерю устойчивости трёхмерного кристалла при растяжении.
Рассмотрим кристаллическую решётку с ГПУ-структурой и трансляцией решётки а, которая растягивается нагрузкой интенсивности q = const, перпендикулярной базисной плоскости (рис. 1) При малых значениях q кристалл ведёт себя как упругое тело и деформация будет происходить за счёт смещения базисных плоскостей относительно друг друга в направлении приложения нагрузки. В этом случае состояние равновесия кристалла устойчивое. При достижении нагрузки некоторого значения q = qc кристалл может перейти в неустойчивое положение равновесия, при котором базисные плоскости повернутся и составят некоторый угол с направлением нагрузки, отличный от 90°. При дальнейшем нагружении деформация кристалла будет развиваться за счёт сдвиговых деформаций путём скольжения базисных плоскостей.
Построим модель, которая характеризует потерю устойчивости монокристалла и на основании которой определяется предельное значение нагрузки qc.
Ч
Рис.1. Плотноупакованные слои атомов в напряжённом состоянии
Рис 2. Потенциал межатомного взаимодействия г) и сила при растяжении соседних атомов
Предположим, что при нагружении атомы, лежащие в базисных плоскостях, относительно друг друга не смещаются, то есть базисные плоскости представляют собой абсолютно твёрдые тела. В процессе деформации взаимодействуют только соседние атомы, силы межатомного взаимодействия являются потенциальными (сила взаимодействия а двух атомов и межатомный потенциал представлены на рис.2 ). При приложении нагрузки ц деформацию кристалла можно представить в виде вертикального смещения базисных плоскостей и их поворота.
Нетрудно видеть, что в качестве элементарной ячейки, описывающей поведение всего кристалла, можно выбрать четырёхатомную ячейку, представляющую тетраэдр, один атом которой закреплён, а три других ведут себя как абсолютно твёрдое тело
(рис.З) Атом В4 в начале координат, а ось 02 проходит через центр тяжести треугольника#,#,.^. Распределённую нагрузку заменим приведённой силой Р. которая приложена в центре тяжести треугольника ВХВ..В.К, в качестве обощённых координат выберем г, ф, 0. Координата г характеризует смещение центра тяжести треугольника ВхВгВ-. вдоль оси 02, угол 0 - поворот треугольника вокруг оси ОТ, угол 9- поворот треугольника вокруг биссектрисы угла В,В1В2.
X
Рис. 3. Четырёхатомная ячейка в недеформированном и деформированном состоянии
Используя методик}', описанную в работах [2], [6], [7], выведем условия потери устойчивости кристалла.
Запишем потенциальную энергию ячейки:
з
п = X “»(р.) -рг , (1>
Л = 1
а2 л 2аг
\а~ т аг рг = В2В4 = Г + ^/=5ІП 6 + аг СОБ #8Н1 ф ,
где и (ра)- потенциальная энергия взаимодействия атома Вл с атомами Д (л-= 1,2,3).
Условия равновесия имеют вид:
дг ор1 дг
(2)
п = ап - £ 50(9дРх - о
5ф др, аф
Значение потенциальной энергии П в исследуемом положении равновесия определяет его устойчивость (неустойчивость), которое в свою очередь определяется через вторые производные П Составим матрицу
где индексы, стоящие у функций П, показывают, что берутся частные производные по соответствующим переменным.
Если матрица А - знакопределённая положительная, то соответствующее положение равновесия устойчивое; точка, в которой нарушается знакоопределённость, является точкой бифуркации [6], [7].
При вычислении коэффициентов устойчивости П^, Пге, Пг#, П8в>, ПгЭ, ГТ,„ , П берём их значени при 0 = ^ = 0, так как исследуется устойчивость соответствующего положения равновесия. Опуская промежуточные выкладки, запишем окончательный результат
П,ТП,9П^ А = П П т П 9ф
|_п „п „п т
Зг2 а+а2-Р)
\
Р К
(3)
п = п „ = п. = о
Гф ГІ/ ву
где
Из соотношений (3) следует, что при выполнении неравенств
Ъг2а+—а > О Р
(4)
Р
матрица А будет знакоопределённой положительной. Знакоопределённость нарушится, когда одно из неравенств (4) переходит в равенство.
При обращении в ноль величины Пп происходит потеря устойчивости по г Теоретически такой случай рассмотреть возможно, но практически это означает разрыв атомных связей, то есть разрушение кристалла предшествует выполнению условий
Величины П м и Г1обращаются в ноль при достижении расстояния между атомами рс, где рс решение уравнения
Вид потенциала межатомного взаимодействия металлов мало изучен [1], [4] и в явном виде найти решение уравнения невозможно. Но качественное определение решения (5) можно получить графически. На рис.4 приведены графики, поясняющие смену знака величин П9в и П и определяющие значение рс. Так как рс < рт, то потеря устойчивости положения равновесия 9 = <р = 0 произойдёт при силе взаимодействия между атомами ос, которая меньше временного сопротивления.
пп=о
(5)
Р
СТс
а
2а
Рис. 4. Графическое решение уравнения (5)
Критическое значение силы Р определяется из уравнения (2).
Зная геометрическую форму кристалла, при помощи несложных преобразований всегда можно перейти от силы Р к нагрузке q и наоборот
Вышеприведённые теоретические выкладки показывают, что сдвиговые деформации могут развиваться и в случае, когда направление растяжения перпендикулярно базисной плоскости. При медленном нагружении кристалла нагрузкой с] деформация развивается за счёт увеличения расстояния между базисными
плоскостями, при достижении нагрузки q значения qc происходит потеря устойчивости положения равновесия Q = (p = 0, при которой плоскости повернутся и
составят некоторый угол с направлением растяжения, отличный от -5-, и при дальнейшем нагружении деформация будет развиваться за счёт сдвиговых смещений
Библиографический список
1. Carrard М., Martin J.L. A study of (001) glide Al single crystals // Philosophical magazine. Ser. A. - 1987. - Vol. 5,6. - № 3.
2. Хоникомб P. Пластическая деформация металлов. - М.: Мир, 1972. 408 с.
3. Мельников А.И., Митюшов Е А. Об одной модели образования полос сдвига в идеальном монокристалле. // Ш Всесоюз. конфер. по нелинейной теории упругости Тез. докл.; - Сыктыквар, 1969. - С. 202-209.
4. Томпсон Дж. М.Т. Неустойчивости и катастрофы в науке и технике - М.: Мир, 1985. -256 с.
5. Мельников А.Н. Об одной модели потери устойчивости идеального плоского кристалла при растяжении // Прикладная механика и техническая физика - АН СССР, Сибир. отдел. Новосибирск.-1992.-№ 1. - С.141-144.
6. Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний. - М.. Физматгиз. 1959. -916 с,
7. Малкин И.Г. Теория устойчивости движения. М.: Наука, 1966. - 530 с.
Получено 10 03 2000
