УДК 539.21
РАСЧЕТ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ В НАНОЗЕРНЕ, ОБУСЛОВЛЕННОГО НАНОДВОЙНИКОВАНИЕМ
Т.В. ДРОБЫШЕВСКАЯ, канд. физ.-мат. наук О.М. ОСТРИКОВ (Гомельский государственный технический университет им. П.О. Сухого)
Разработана методика расчета полей смещений и напряжений в нанозерне при наличии в нем нанодвойника. Изучено напряженное состояние в нанозерне, имеющем форму правильного многоугольника и содержащем нанодвойник. Определены точки концентрации напряжений в зерне наноматериала. Проанализирована зависимость напряжений от параметров двойника.
Ключевые слова: наноматериалы, нанозерно, зернограничные дислокации, двойникующие дислокации, нанодвойникование.
Введение. Наноматериалы - ультрамелкозернистые поликристаллические материалы [1]. Широкое использование наноструктурных металлов и сплавов в качестве конструкционных и функциональных материалов в настоящее время сдерживается из-за их особенных свойств. Так, наряду с высокой прочностью и твердостью наноструктурные металлы облают очень низкой пластичностью, что создает существенные препятствия для их применения в конструкциях деталей машин [2]. Повышение пластичности наноструктурных материалов может быть достигнуто за счет искусственного нанодвойникования. Благодаря чему наноматериалы (металлы) приобретают высокие показатели механических свойств, позволяющие существенно расширить их область применения [1]. При этом проанализировать приобретаемые в результате нанодвойникования структуры и свойства нанометаллов можно путем изучения напряженно-деформированного состояния зерна наноматериала (металла), содержащего нанодвойники.
Цель данной работы - изучение напряженно-деформированного состояния в зерне нанострук-турного материала, обусловленного нанодвойникованием.
Рассмотрим представленное на рисунке 1 единичное зерно наноматериала. Данное наноразмер-ное зерно имеет вид я-угольника и содержит в своем объеме единичный нанодвойник. В рассматриваемом случае границы зерна смоделированы в виде стенки полных дислокаций, а границы двойника в виде цепочки частичных дислокаций. Зернограничные и двойникующие дислокации в рамках изучаемой задачи в пределах каждой отдельной границы параллельны оси 02 (рисунок 2).
В рамках проводимого исследования примем, что рассматриваемое нанозерно расположено на некотором расстоянии от поверхности материала. Кроме того, при расчете напряженно-деформированного состояния зерна не будем учитывать никакие напряжения кроме напряжений, обусловленных зерногра-ничными и двойникующими дислокациями [3].
В соответствии с поставленной задачей и принятыми условиями, с учетом методик, описанных в [4, 5], смещения и напряжения, создаваемые двойниковыми и зеренными границами определяем путем суммирования
Граница
ГРЕ зер
Рисунок 1. - Схематическое изображение единичного зерна наноматериала, содержащего единичный нанодвойник
смещений и напряжений, создаваемых отдельными дислокациями (
Рисунок 2. - Взаимное расположение зернограничных и двойникующих дислокаций и их компонент вектора Бюргерса
Щ = Ы;
(1) (2)
\щ ) (х, у) + (щ ) (х, у),
х 'Тм ^ 'Ь
(Щ ) (х у)= Е (щР) (х- гйгм>у + ^гм^м -гкм)+ Е (и(2)) (х- гйгм>У - + гкм)
у 'м 'гм г=0
(и), (х у)= Е Ет Е (Щк)) (хк, Ук - гкъ);
^ /Ь к=1 I г=0К 1
= (°У I(х,у ИО* )Ь(х,у),
(Оу ) (х, у )= Е (°//}) (х - гйгм>у + К^Ьгм - гкм)+ ЕЕ- (о(2) (х - Г^>У - + гкм), (5)
\ п-ю -1 'м г=0 м
(Оу)и(х,У) = Е Етит]м Е №) (хк,ук -гкъ);
Ь к=11,м г=0 'Ъ
И*
г=0
(3)
(4)
(6)
где I иу - индексы, принимающие значения х, у или г; I им - индексы, принимающие значения хк , ук или Ъу-; к - номер зеренной границы; т - номер двойниковой границы; (щ ) (х,у), (Оу ) (х,у) -
смещения и напряжения, создаваемые двойникующими дислокациями; (щ ) (х, у), (оу ) (х, у) - смещения и напряжения, создаваемые зернограничными дислокациями; , Иь - количество двойникующих и зернограничных дислокаций, расположенных на отдельной границе (двойника либо зерна); йгм, Игм, кь - параметры, характеризующие взаимное расположение дислокаций (см. рисунок 2); системы координат ХУ2 и Хк Ук X (см. рисунок 2) - системы координат, используемые в расчетах; Тц , Тум - элементы матрицы преобразования [6].
В данном случае рассматриваем поворот правосторонней системы координат ХУХ вокруг оси X по часовой стрелке на угол ак (см. рисунок 2). При этом матрица преобразования принимает вид [6]
Т =
С08 (ак ) !5Ш (ак ) 0 -^П(ак ) СШ(ак ) 0 0 0 1
(7)
В соответствии с представленными на рисунке 2 системами координат и ориентировками векторов Бюргерса, а также с [3, 6], дальнейшие расчеты проводим следующим образом:
{ \
(их) (х у)= X
1 г=0
- ьп> 2я
агс1ё У+^Ьп - гЛ>п + х -
(у+Ntw^w - гк,)(х -г(>п) 2 (1 -У)(( У + Ntwhiw - гЛп )2 +( X - )2)
^-1 ье + X ь/п
г=0 2р
агс1ё У - + гЛ>п +
х - г(ил>
(иу ) ( x, У ) = - X
1 Нп г=0
^ ье 2я
1 - 2У 2р
(У - Ntwhw + гЛп)(х- г(/п)
2 (1 -у)(( У - + гК )2 +(х - Чп )2 )
1п ((У + Ntwhtw - гЛп)2 +(х- г(/п)2) +
(х- г(п)2-(У + - гК)2 4(1 -у)((У + ^пЛп - гЛ^)2 + (х- (п)2)
^ -1 ье _ X ь>п
г=0 2р
1 - 2У 2Р
1п ((У - ^^ + гЛ>п) +( х г(/п) )
2\ (х- г(/п) -(У- + гЛп)
4 (1 -у)(( У - + гЛп )2 +( х - (п )2 )
N>4, ий
/п 1л-
и..) (х, У )= агсге
пп г=0 2я
г ьп агс,„ гЛп , ^1 _____У - ^„Лт + гЛ/п
2 я х - (
-+ X ^агсге
г=0
2я
х - г(,
(8)
/ \ я ^-1 Ьи
(их I (x, У)= X с°8(ак ) X 7Т-
1 'Ь к=1 г=0 2р
хк
(У^ - гЛь)
Ь \хк
п ^ -1 Ьь
+X К) X 2Р
к=1 г=0 2Р
1 - 2У , 2р
1п 1( Ук - гЛЬ ) + хк
2(1 -у) ((Ук - гЛь) + хк:
-(Ук -гЛь) 4(1 -у)((Ук - гЛь )2 + хк:
хк
, > я NЬ-1 ие
(иУ ) (x, У)= X Нп (ак )) X
к=1
V %
£1 ы
2р
г=0
Ук - гЛь (Ук - гЛь) хк
хк
2(1 -у)((Ук - гЛь) + хк:
+ X с°8(ак) X
к=1 г=0 2Р
1 - 2у , и ; , \2 ; 21 хк
1п I ( Ук - гЛь ) + хк I +
2
( Ук; - гЛь )2
4(1 -У) ((У^ - гЛь) + хк:
я ьь Ук - гЛь
(и. ) (х, У)= X X ь-атс1§^^ ь
К /ь 1—1 —п 2Р хк
( \ Ут (°хх )/п =" X
2 2 ^ иье (У + Ntwhtw - гЛ>п) 3(х- г(/п) +(У + ^пк/п - гЛп) X Г-ьРп___
=02р(1 У) ^(х-г(/п)2 + ( У + ^пЬ/п - гЛп?
+
+
+
я
2 2 Им-1 ИЪе (у-Игмкгм + гкм) 3(х-) +(у-+ гК) - е г-Ъгм__1__.
г=0 2р(1 V) |-( х - )2 + ( у - + гкм )2
(о уу I =Е
гм г=02я(1 -V)
2
Им цЪм (у + - ^) 1_(х- ) -(у + - ^) цъм (у-Игм'кгм + ^) (х-^)2-(у-Имкм + гкгм)
22 (х- чм) +(у + ЫгмК - гкм)
+ Е0 М1 -V)
22 (х - Чм) +(у - Игмкгм + гкгм)
Им, ,.ие
(О ) =- Е ^мУ__у + Игмкгм - гкгм__тЪмп__у - "гмкм + гкгм .
1 ЪЪ'гм г=0 р(1 х - гйм )2 + ( у + Имкм - гкм )2 г=0 р(1 -п)( х - гйм )2 +(у - Имкм, + гкм )2;
Им -1 .. ив
(о ху )г = Е
1 г=0
Чы
м 1 МЬ}^__у-Имкм + гк^
-п)(х - Чм)2 + 2 / л, , , \2"
Ым цЪ^" (х- ^) (х- ) -(у + - гкм)
гм г=0 2я(1 -V)
22 (х-^гм) +(у + Игмкгм -гкгм)
Игм-1 мЬв (х-) (х-)2-(у-Игмкгм + гкм)2
+ Е м
г=0 2Р(1 -V)
22 (х - ) +( у - ИгмК + гкм)
(ох) =-Е 0
у + Игмкгм - гкм
-1 пЬ5 - е М-Ьг№
у - Игмкгм + гкм
г=0 2р (х- )2 +(у + - гкм)2 г=0 2р (х- )2 + (у- + гкм)
2
(Оу) = Е о
Им . . >,4
МЬл
х - гйи
Им -1.. ии
+Е
МЪп
х - гйи
г=0 2р (х- )2 + (у + Игмк1м - гкгм)2 г=0 2р (х- Мм)2 +(у- + ^)
(10)
/ \ п г. "Ъ-1 мЪв
(Охх )Ь =-Е cos2 (ак ) Е ^ТТГ^
v 'Ь к=1 г=0 2р(1 ^
"ь -1
г=0
(ук - гкь) 3хк2 + (ук - гкь)
хк +(ук - гкь)
^N^-1 МЬв
+Е 81П2 (ак) Е , ч к=1 г=0 2р(1 -V)
(ук-гкь) х'к2-(ук-гкь)
хк +(ук - гкь)
+ Е ^ (2ак) Е
"Ь-1 мЪЬв хк
хк -(у к - гкь)
к =1 г-
0 2^(!-V)
хк +(ук -гкь)
2 '
п МЬв
(о уу )Ь=- Е §1п2 (ак) Е
(ук - гкь) 3хк2 + (ук - гкь)
к=1 г
=0 2Р(1 -V)
хк +(ук -гкь)
2 ( ) N,-1 мЬьв (- гкЬ )х^2-(- гкЪ )
+Е С082 (ак) Е —, ч---
1 к> г=0 2р(1 -V)
хГ + (Ук - гкь)
- Е (2ак) Е
"ъ--1 мъъ хк
хк -(ук -гкь)
к=1
0 2*(1 -V)
+ (ук - гкь)
2
2
+
2
+
2
2
+
2
п
2
+
2
п
г
n.%1 mbebn yk - rhb
(szz )b =- Z Z P
V 'b k=l r=0 p
k=1 r=0 '"(1 v) xk2 +( yk - rhb
-(yk - rhb)
n ,e ( yk - rhb ) 3xk2 +( yk - rhb )
(Sxy 1 = Z C0S (ak ) Sin (ak ) b v)-"-
b k=1 2p(1 -V) Г '2 I ' , \2
xk +(yk -rhb)
+ Zcos(ak)sin(ak) Z v)
k=1 r=0 2p(1-v)
Nb-1 b (yk -rhb) xk2-(yk - rhb)
xk +(yk -rhb)
n Nb-1 Ubf xk
+ Z cos (2ak) Z m1b k=1
xk -(yk - rhb)
r=0 2p(1-V)
xk2 + (yk - rhb)
(°zx 1 =- Zcos (ak) Z
v 'b k=1 r=0 2p
N¿_-1 me (yk- rhb)
i) , ,N¿-1 me xk
+ Z Sin (ak) Z k
xr + (yk - rhb)
k=1
r=0
2p
x*2 + (yk - rhb)
(° zy 1= Zsin (ak) Z
v y'b k=1 r=0 2p
Nb_-1 ^ (- rhb)
i) n / míe xk + Z cos (ak) Z k
xr + (yk - rhb)
k=1
r=0
2p
xk +(yk -rhb)
(11)
Здесь V - коэффициент Пуассона; m - модуль сдвига; bfw - модуль краевой составляющей вектора Бюргерса двойникующей дислокации; b?w - модуль винтовой составляющей вектора Бюргерса двой-
никующей дислокации; b| - модуль вектора Бюргерса полной краевой дислокации; bí - модуль вектора Бюргерса полной винтовой дислокации.
Границы двойника и зерна принимаем прямолинейными, а дислокации на рассматриваемых границах равноудаленными друг от друга. Тогда координаты в системах Xk Yk Z и координаты границ зерен в данных системах координат определяем в соответствии с выражением (4):
xk = x • cos (ak)-y • sin (ak), (12)
y'k = x • sin (ak) + y • cos (ak).
(13)
Результаты расчетов и их обсуждение. В качестве базового варианта для исследования было принято шестиугольное наноразмерное зерно, содержащее в своем теле единичный нанодвойник, материал - железо (^ё). Для расчетов принимались следующие параметры: Ьё = ЬЬ = 0,248 нм; Ьё„ = Ь= 0,124 нм [7]; т = 81 ГПа [6]; у = 0,29 [7]. Базовые параметры зерна и двойника: = 35,
Н^ = 2Ъё^ = 0,248 нм; = 1,736 нм; Иь = 32 , Нь = 1,736 нм.
Анализ напряженного состояния был проведен на базе анализа гидростатических Он (рисунок 3, а),
нормальных охх, О уу, О, комбинированных касательных напряжений т (рисунок 3, б), а также касательных напряжений о^ (рисунок 3, в), оХ2 и оуХ . Дополнительные расчеты проводили в соответствии с [7]:
0Н = Е °«73, (14)
t = -v/s yz 2 +sxz 2
(15)
+
n
2
2
В рассматриваемом случае гидростатические (о к) и сдвиговые (о^ , т) локализованы на границах элементов системы двойник - зерно. Максимальные гидростатические Ок и сдвиговые о^ напряжения наблюдаются у устья двойника, а максимальные комбинированные касательные напряжения т -в вершинах зерна.
Концентрация гидростатических напряжений Ок внутри зерна несколько выше нежели за его пределами. В анализируемом поле напряжений имеют место как сжимающие, так и растягивающие напряжения. Однако выявить четкую характеристику знакопеременности напряжений не представляется возможным (см. рисунок 3, а).
Локализация нормальных напряжений о хх внутри зерна значительно выше, нежели за его пределами. При этом выше оси ОХ преобладают сжимающие напряжения, а ниже - растягивающие. На участках максимальной концентрации (зеренные и двойниковые границы) напряжения о хх носят сжимающий характер. При этом максимальные значения о хх отмечены у устья двойника.
Нормальные напряжения оуу локализованы на зеренных и двойниковых границах и имеют относительно равномерное распределение внутри зерна и за его пределами. В пределах зерна имеют место как растягивающие (в I и III четвертях относительно осей, проходящих через вершину двойника), так и сжимающие напряжения (соответственно во II и IV четвертях). За пределами зерна преобладают отрицательные значения о уу. Напряжения на границах двойника и зерна носят сжимающий характер
и имеют максимальные значения на границе зерна, к которой прилегает устье двойника.
Концентрация нормальных напряжений оЪЪ внутри зерна и за его пределами отличается несущественно. При этом выше оси ОХ напряжения отрицательны, а ниже данной оси - положительны. На границах же зерна и двойника напряжения оЪЪ, аналогично напряжениям охх, о уу , являются сжимающими, а максимальные значения оЪЪ соответствуют устью двойника. При этом максимальные напряжения оЪЪ на порядок ниже максимальных значений охх , о уу.
Комбинированные касательные напряжения т можно охарактеризовать высокой концентрацией напряжений на зеренных и двойниковых границах, а также в узловых точках. Концентрация напряжений т за пределами зерна выше концентрации внутри него. При этом максимальные численные значения т имеют место в вершинах зерна (см. рисунок 3, б).
Распределение напряжений охъ характеризуется максимальными значениями, соответствующими устью двойника. Преобладающими являются сжимающие напряжения охъ, они расположены на границах двойника и зерна (места наибольшей локализации охъ), а также выше оси ОХ. Ниже оси ОХ имеют место растягивающие напряжения (кроме границ).
Как было отмечено выше сдвиговые напряжения о уЪ имеют максимальную концентрацию на
границах двойника и зерна, при этом напряжения носят сжимающий характер. Концентрация данной компоненты тензора напряжений внутри зерна превышает концентрацию за его пределами. При удалении от границ зерна концентрация напряжений оуЪ существенно снижается, также снижаются и численные значения напряжений. Максимальные значения оуЪ отмечены на границе зерна, к которой прилегает устье двойника.
Сдвиговые напряжения оху знакопеременны относительно оси, проходящей параллельно оси ОУ
через середину двойника. Слева от данной оси напряжения оху отрицательны, а справа - положительны.
В местах локализации оху (зеренные и двойниковые границы) преобладают сжимающие напряжения,
максимальные их значения можно отметить на границе, к которой прилегает устье двойника (см. рисунок 3, в).
В рамках данного исследования также были проведены расчеты полей напряжений для пятиугольного и семиугольного зерен. Результаты свидетельствуют о том, что форма зерна не существенно влияет на характер распределения полей напряжений. В случае пятиугольного и семиугольного зерна аналогично варианту с шестиугольным зерном основные элементы рассматриваемой системы (зеренные и двойниковые границы), являются не только концентраторами напряжений, но соответственно и местом локализации напряжений.
х, нм
а - вл(х,у); б - а^у); в - Сху(х,у) Рисунок 3. - Распределение напряжений в шестиугольном зерне нанокристалла, обусловленных наличием единичного клиновидного нанодвойника
а
б
в
Кроме того, было рассчитано влияние параметров двойника на напряженное состояние в зерне наноматериала. В качестве варьируемых параметров приняты количество дислокаций на границах двойника (Ntw) и расстояние между двойникующими дислокациями, измеренное вдоль оси OX (dtw).
Так, при увеличении количества двойникующих дислокаций (Ntw = 2, 7,15, 25, 35,45) увеличивается концентрация гидростатических напряжений о^ в пределах зерна, а за его пределами концентрация не изменяется. В отношении гидростатического напряжения о^ на основных элементах системы можно отметить, что количество двойникующих дислокаций Ntw оказывает влияние на о^ только на границе зерна, к которой прилегает устье двойника. При этом характер распределения напряжений на данной границе не изменяется. Имеет место лишь незначительное изменение численных значений о^ при изменении Ntw. На остальных границах зерна количество двойникующих дислокаций Ntw не оказывает влияние на гидростатическое давление. Следует также отметить отсутствие зависимости напряжений о^ от количества двойникующих дислокаций Ntw , измеренное вдоль оси OX.
Анализируя распределение гидростатических напряжений о^ вдоль границ двойника, можно отметить следующее. Характер распределения анализируемых напряжений не зависит от Ntw, лишь размерные характеристики распределения обусловлены длиной двойника и соответственно Ntw .
В отношении комбинированных касательных напряжений t и сдвиговых напряжений о^ можно
отметить, что количество двойникующих дислокаций Ntw не оказывает существенного влияния на характер распределения полей напряжений. При этом численные значения анализируемых компонент тензора напряжений (оxy и t) имеют тенденцию к увеличению при увеличении Ntw.
По результатам исследования влияния расстояния между двойникующими дислокациями, измеренного вдоль оси OX (dtw = 0,248; 0,744; 1,240; 1,736; 1,984; 2,232 нм), можно отметить, что на анализируемые характеристики напряженного состояния (о^, t, оxy) параметр dtw не оказывает существенного влияния.
Таким образом, изучено напряженно-деформированное состояние в зерне наноструктурного материала, обусловленного нанодвойникованием. Разработана методика расчета полей смещений и напряжений в нанозерне, имеющем форму многоугольника и содержащем в своем теле единичный нанодвойник. Выявлены места концентрации напряжений в нанозерне. Проведен анализ зависимости распределения полей напряжений от параметров нанодвойника.
ЛИТЕРАТУРА
1. Овидько, И.А. Влияние миграции двойниковых границ на трещиностойкость нанодвойникованных металлов / И.А. Овидько, А.Г. Шейнерман // Materials Physics and Mechanics. - 2014. - № 21. - С. 248-258.
2. Объемные наноструктурные металлы и сплавы с уникальными механическими свойствами для перспективных применений / Р.З. Валиев [и др.] // Вестник УГАТУ. - 2006. - № 3. - С. 23-34.
3. Дробышевская, Т.В. Задача о расчете напряженно-деформированного состояния, обусловленного единичным некогерентным двойником в пяти-, шести- и семигранном зерне поликристалла с гексагональной структурой / Т.В. Дробышевская, О.М. Остриков // Физ. мезомех. - 2015. - Т. 18. - № 2. - С. 79-88.
4. Остриков, О.М. Механика двойникования твердых тел : моногр. / О.М. Остриков. - Гомель : ГГТУ им. П.О. Сухого, 2008. - 301 с.
5. Остриков, О.М. Поле упругих деформаций и смещений в системе «Клиновидный нанодвойник - цепочка полных дислокаций» / О.М. Остриков // Вестник ГГТУ им. П.О. Сухого. - 2011. - №4 (47).
6. Хирт, Дж. Теория дислокаций / Дж. Хирт, И. Лоте. - М. : Атомиздат, 1972. - 600 с.
7. Миркин, Л.И. Физические основы прочности и пластичности (Введение в теорию дислокаций) / Л.И. Миркин. -М. : Изд-во МГУ, 1968. - 538 с.
Поступила 13.10.2017
THE CALCULATION OF THE STRESS-STRAIN STATE IN NANOZINE DUE TO NEODOLIMAYA
Т. DROBYSHEVSKAYA, O. OSTRIKOV
The developed method of calculation of fields of displacements and stresses in nanozine in the presence of nanovoice. The stress state in a nanosecond having the shape of a regular polygon and containing a nanodvoir is studied. The stress concentration points in the grain of nanomaterial are determined. The dependence of stresses on the parameters of the double is analyzed.
Keywords: nanomaterials, nanozero, grain boundary dislocations, double dislocations, nanodvoynikovanie.