Об одном критерии потери устойчивости в атомарных моделях

УДК 539.3

ОБ ОДНОМ КРИТЕРИИ ПОТЕРИ УСТОЙЧИВОСТИ В АТОМАРНЫХ МОДЕЛЯХ

В.В. Стружанов

Abstract

The concept of stability (instability) ofphysical state is introduced for the elastoplastic lattices. The criterion is developed for detection of such states and determination of the moment when the instability of equilibrium states of the whole system occurs. The criterion has been applied studying the stability of special model deforming.

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

В рамках структурного подхода эффективным методическим приемом для анализа качественных закономерностей процесса разрушения служит исследование различных решеточных атомарных моделей материала [1 - 7]. В ряде работ данного направления [2,3,7] отмечается, что скачкообразное (катастрофическое) возникновение разрывов в решетке связано с потерей устойчивости положения равновесия отдельных атомов или атомных рядов, происходящее в процессе нагружения. При этом изучение данных явлений, как правило, опирается на общую теорию устойчивости дискретных консервативных систем, которая, к сожалению, не может быть распространена на континуальные модели.

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

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

Максимальное значение функция q = q(x) принимает при х = хс . Тогда

при х < хс, к р < 0, х > хс . Кроме того X р (хс) = О, X р (0) = X,

Xр <Х. Здесь величины Xр имеют смысл касательных модулей, X - модуля продольной упругости.

В целях моделирования пластических (неупругих) свойств материала полагаем, что при снятии нагрузки сила взаимодействия изменяется по линейному закону

Х(х- хр) при фиксированном значении хр, которое представляет собой неупругую составляющую смещения х Считая далее, что деформация системы равняется х , получаем выражение для потенциальной энергии упругих деформаций Пе =0,5А,(х-хр)2.

Определим теперь закон приращения неупругих составляющих смещений Имеем в общем случае q = X (х - хр) . Отсюда ё q = X (с! х - ё хр). С другой стороны с! q - X х Приравнивая оба выражения, находим

dxp

ч

р

d х =--------d q .

Х^р

И наконец отметим, что мощность диссипации

d D = q dxp = q(X - X p)d q/X X p , при растяжении двух атомов - величина положительная

2. Опираясь на теорему Лагранжа дадим следующее определение Физическое состояние системы (ФСС) в некоторой конфшурации, характеризуемой совокупностью параметров состояния и управления (нагрузок), устойчивое, если при произвольном малом изменении этих параметров на любом активном пути нагружения потенциальная энергия упругих деформаций увеличивается Если существуют пути нагружения, на которых потенциальная энергия упругих деформаций уменьшается, то исходное физическое состояние системы является неустойчивым

При разгрузке, очевидно, состояние системы всегда устойчивое.

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

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

3. Для произвольного положения системы двух атомов приращение потенциальной энергии упругих деформаций dne = Х(х - xp)(dx - dxp) Отсюда dne>0, когда dx>dxp, и dne<0, когда dx<dxp Вспоминая формулу (1) и определение устойчивог о физического состояния, получаем, что при Xр (х) > О

состояние системы устойчивое, а при X р (х) < 0 - неустойчивое.

Введем теперь следующий функционал

R = X dxdx - X dxpdxp ,

где dxp и dx связаны выражением (1). Будем называть его R - суммой данной системы. Очевидно, что R > 0 отвечает устойчивым ФСС, R < 0 - неустойчивым.

Далее из постулата устойчивости Друккера о положительности работы внешних сил в случае возникновения необратимой деформации на замкнутом цикле их приложения и снятия следует, что при устойчивом физическом состоянии системы атомов выполняются неравенства сЦс1х>0, сЦ<1хР >0 Когда система находится в неустойчивом физическом состоянии, то сЦс!х <0, < 0.

Преобразуем выражение <1цс1х с использованием равенства ёх = с1хе+с1хр, с!ц = X (]х с . Имеем

с1цс1х = Цс1х - с!хр )(с!хс + ёх1’) - X ёхёх - ёяёхр - X сЗхрскр .

Отсюда, если постулат Друккера выполняется, то Я>0. В противном случае

Я<0.

С другой стороны, подставляя выражение (1) в формулу (2) и используя соотношение - X рёх, получаем

И. = А, _1Л Р(2Х - X р)сЫх = \ ~1(2Х - X р)с1Чс1х .

Следовательно, сЦёх > 0, когда Я > 0, и с^<1х < 0 при К < 0. Кроме того

сЦёх = X с1хис!хе + dqdxp .

Так как X ёх^хе > 0, то из условия Я < 0 следует dqdxp < 0 и наоборот.

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

4.Рассмотрим теперь решетку, состоящую из трех параллельных рядов атомов. Один из крайних рядов закреплен, к атомам второго приложены параллельные растягивающие усилия, равнодействующая которых равна р. Ряды могут перемещаться только параллельно друг другу. То есть атомы крайнего ряда имеют одно и то же перемещение и, а атомы среднего ряда - перемещение х .

Взаимодействие атомов среднего ряда с соответствующими им закрепленными атомами определяется функциями ц, (х) () - 1,..., п, п - число

атомов в рядах). Таким образом, имеем известный разброс свойств. Характеристики

1 Р

этих связей описаны выше, а именно: их касательные модули равны X! , модули

продольной упругости - А,; , формула (1) связывает приращение смещений с их

неупругими составляющими

Взаимодействие же атомов подвижного крайнего ряда с атомами среднего ряда полагаем линейным с коэффициентом пропорциональности с.

Сначала исследуем устойчивость процесса растяжения, используя аппарат, изложенный в [3]. Запишем выражение полной энергии системы:

г пс ^ г

П = |час1х + —

О О

П

где ца = ^ Ч ] Критические точки функции П определяются уравнениями

I I

сШ

—- = Яге-пс

О х

г5 П

=пс

о и

Очевидно, что при заданных значениях р все положения равновесия

получаются из решения уравнений (4), (5).

Для определения вырожденных критических точек, где происходит смена

устойчивых положений равновесия на неустойчивые, к уравнениям (4), (5) необходимо

добавить уравнение, получающееся приравниванием нулю детерминанта матрицы Гессе Имеем

дгП/дк[ дх^ = пс^~~—+ пс]-п2с2 - 0 (/,/=1,2).

п

Здесь X \ ’ X} = х, Х2 ~ и. Отсюда при X (Х>0 процесс квазистатического

1=]

деформирования устойчив, при А, ^<0 положения равновесия системы неустойчивы Потеря устойчивости реализуется для тех значений (х, и, р) - решений уравнений (4),

(5), для которых X £= 0. В этом случае образуется разрыв между закрепленным и

средним рядами атомов.

Если реализуется кинематическое нагружение, при котором задается перемещение и, а состояние системы определяет один параметр х, то в выражении (3) исчезает последний член и, следовательно, остается только одно уравнение равновесия

(4) Проделывая аналогичные операции, находим, что при А,^>-пс процесс

растяжения устойчив, при X £<-пс - неустойчив. Потеря устойчивости происходит,

1 Р

когда X а~ -пс.

5. Запишем N -сумму данной системы, суммируя выражения (2) для каждой связи Сначала рассмотрим случай, когда ёи > 0, с!х > 0, с!и>с1х, то есть все связи работают в условиях активного нагружения Имеем

Я = X ас1хс1х - X ас1хрс1хр + пс

п

Здесь с1и,с1х - малое возмущение параметров состояния системы,

¡-1

находящейся в некотором положении равновесия, определяемом параметрами (х,и,р). В силу уравнений равновесия (4),(5) с учетом равенства с1яа=\(*с1х получаем

¿рффа)

пс с!и ёр А 5 + пс

ёх ---------, ёи--=——--------.

Х^+ПС ПС ^

Очевидно, что данный случай реализуется при ёр > О, А £> 0. Подставляя выражения

(7) в формулу (6) и используя соотношение (1), находим квадратичную форму

\ 2

я = яр = в

где В^Аа-Ха1(Ха~хЦ)2 + (пс)~](А(*)2.

р

Отсюда Яр > 0 при Аа>0 и физическое состояние атомной решегки устойчиво Кроме того, устойчиво и исходное (невозмущенное) положение

р

равновесия. Потеря устойчивости происходит, когда Аа=0 и Яр = сс

Пусть теперь ёи < О, ёх > 0, или ёи > 0, ёх > 0, ёх > ёи, то есть упругая связь

Р Р

разгружается. В первом случае ёр < О, X а< -пс, во втором ёр < 0, 0 > X а> - пс Так

как упругая связь теряет энергию, то Я -сумма в данном случае принимает вид Я] = А (Хёхёх - А аёхрёхр - пс(ёи - ёх)(ёи- ёх).

Соответственно

Я1р=В1ёрёр/(А(Р'2

__ 1 п 9 __ | р о

гдеВ| = Аа-Ха (Аа~ А(Х) - (пс) (Аа) .Отсюда Я ]р < 0 и физическое состояние

системы неустойчиво. Очевидно, что и исходное положение равновесия системы неустойчивое.

Для кинематического нагружения в первом варианте получаем Я = Яи = В(пс)2ёи ёи / |а £+ пс|

Когда А £>0, то Яи>0 и физическое состояние системы устойчиво.

Устойчиво также исходное положение равновесия. Кроме того, при А(Х=0 потери устойчивости не происходит.

Если 0>А^>-пс (ёх>ёи>0), то Я[и = В1(пс)2ёиёи /|а^+ пс| <0.

Физическое состояние системы атомов неустойчиво, деформирование может идти только с потерей энергии упругих деформаций. Однако исходное (невозмущенное) положение равновесия системы сохраняет устойчивость. Потеря устойчивости

равновесия происходит при А -пс и Я1и = -ж .

Когда А^<-пс (ёи<0,ёх>0), то Я1и<0 и физическое состояние системы неустойчиво Кроме того, неустойчиво и исходное (невозмущенное) положение равновесия.

Очевидно, что приведенный анализ с помощью Я -сумм полностью совпадает с анализом устойчивости процесса растяжения атомной решетки, изложенным в пункте 4.

6. Заключение

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

Библиографический сиисок

1 Новожилов В.В. О необходимом и достаточном условии хрупкой прочности // ПММ. - 1969 - Т.ЗЗ, Вып.2. - С.212-222.

2 Новожилов В.В. К основам теории равновесных трещин в упругих телах /У ПММ - 1969. - Т.ЗЗ, Вып. 5. - С.797-812.

3 3 Стружанов В.В Об одном подходе к изучению механизма зарождения трещин // ПМТФ - 1986 - № 6 - С 118-123.

4 Томпсон Д М.Т. Неустойчивость и катастрофы в науке и технике,- М. Мир, 1985 - 254с,

5 Дайнис Г , Пэскин А Моделирование трещин с помощью вычислительных машин// Атомистика разрушения. - М.: Мир, 1987. -С.177-212.

6 Корнев В.М., Тихомиров Ю В. О критерии хрупкого разрушения тел с трещиной при наличии дефекта атомной решетки // Изв. РАН МТТ. - 1994. -№ 2,- С. 185 -193.

7 Корнев В.М., Тихомиров Ю.В. Потеря устойчивости цепочки атомов при наличии примеси Снижение прочности хрупких трещиноватых тел // ПМТФ -1996. - Т.37, №3. - С. 160-173.