ЕСТЕСТВЕННЫЕ И ТОЧНЫЕ НАУКИ
УДК 538.3
МЕТОДИКА РАСЧЕТА НА ОСНОВАНИИ МАКРОСКОПИЧЕСКОЙ ДИСЛОКАЦИОННОЙ МОДЕЛИ ПОЛЕЙ НАПРЯЖЕНИЙ У ПОЛИСИНТЕТИЧЕСКОГО ДВОЙНИКА, НАХОДЯЩЕГОСЯ ВДАЛИ ОТ ПОВЕРХНОСТИ ДВОЙНИКУЮЩЕГОСЯ МАТЕРИАЛА
Канд. физ.-мат. наук ОСТРИКОВ О. М.
УО «Гомельский государственный технический университет имени П. О. Сухого»
Разработанная в [1] дислокационная макроскопическая модель клиновидного двойника позволила решить проблему, имеющуюся в существующей в настоящее время дислокационной модели тонкого двойника [2], которая не давала возможности рассчитывать поля напряжений внутри двойника. Это существенно ограничивало область применения дислокационных моделей двойников, так как в большинстве случаев на практике приходится иметь дело с двойниками, ширина которых соизмерима с их длиной [3-5].
При деформировании двойникующихся материалов двойники в них, как правило, формируются группами [6, 7]. Поэтому практический интерес представляет решение задачи по математическому моделированию полисинтетического двойникования, являющегося результатом развития системы параллельных двойников. Полисинтетическому двойникованию подвержены материалы при их деформировании сжатием, растяжением или кручением [6], т. е. таким видам деформаций, которые часто встречаются при эксплуатации двойникующих-ся материалов.
Целью данной работы стала разработка методики расчета полей напряжений у полисинтетического двойника, находящегося вдали от поверхности кристалла.
Решение такой задачи интересно для случая крупнозернистых поликристаллических материалов, когда при внешнем деформировании система параллельных двойников зарождается
на границах зерен. При этом отвлечемся от напряжений границ зерен и рассмотрим только те напряжения и эффекты, которые связаны только с двойникованием.
Для расчета полей напряжений с. (х, у) у
единичного клиновидного двойника с произвольной формой границ в [1] были предложены следующие соотношения:
с. (х, у) = с(1)(х, у) + с<2)( х, у); (1)
1^
(х, у) = {V1 + САК))2 Р1(х0) X
0
хс.-0)(х, у, Хо)^Хо;
ь _
(х, у) = {V1 + (/г(Хо))2р2(Хо) X
хс(2,0)(х, у, хо)ёхо;
сао) = -
мАс
2 п(1 -V)
(у - /1(хо))[3(х - хо)2 + (у - /1(хо))2].
[(х - хо)2 + (у - /1(хо))2]2 ;
с°'0) = МАкр х
2 п(1 -V)
(у - /1(хо))[(х - хо)2 - (у - /(хо))2].
2П2
[(х - хо)2 + (у - /1(хо))2]
(2)
(3)
X
сао = х
^ 2 п(1 -v)
х (x - x0)[(x - Xq)2 - (y - /1(xp))2]; [(x - Xq)2 + (y - f1(Xo))2]2 '
аао) = y - f1(Xo)
" n(1 -v)(x - Xq)2 + (y - f1(Xo))2
C(1,0) = - Mb. У - f1(X0) ;
ZX 2л (X - Xq)2 + (y - _/1(Xq))2'
_(1,0) = Mb.
zy /„ „ \2 , \\2
2л (X - Xq)2 + (y - /1(Xq))2 '
с**) = х
" 2 n(1 -v)
. (y - /2(Xq))[3(x - Xq)2 + (y - /2(Xq))2] ;
■ [(X - Xq)2 + (y - /2(Xq))2]2 '
c(,0) = х
yy 2 n(1 -v)
х (y - /2 (Xq))[(X - Xq)2 - (y - /2(Xq))2] ; [(X - Xq)2 + (y - /2(Xq))2]2 '
C™ =х ^ 2 n(1 -v)
х (x - Xq)[(x - Xq)2 - (y - /2(Xq))2] ; [(X - Xq)2 + (y - /2(xq))T '
c(2,q) = му у - /2(xq)
n(1 -v) (X - Xq)2 + (y - /2(Xq))2'
(2,0) = - мЬ. У - /2(xq) ;
ZX (x - Xq)2 + (y - /2(Xq))2 '
c(2,0) =МЬв
(x - Xq) + (y - /2(Xq))
В данных соотношениях х, у) и
с(2) (х, у) - напряжения на каждой из двойниковых границ; Ь - длина двойника; /1(х0), /2( х0) - функции, определяющие профиль каждой из двойниковых границ; р1( х0), р2( х0) -плотности дислокаций на каждой двойниковой границе; с(1,0)(х, у, х0), с(2'0)(х,у, х0) - напря-
(4)
(5)
жения, создаваемые отдельными дислокациями, находящимися на разных двойниковых границах; ц - модуль сдвига; V - коэффициент Пуассона; Ькр и Ьв - краевая и винтовая составляющие вектора Бюргерса частичных двойникующих дислокаций.
Полученные соотношения могут быть положены в основу макроскопической дислокационной модели полисинтетического двойника в следующей модификации:
N L
c®(x, y) = 1 + WW)2 P1(Xq) х
n=0 о
хс, '0)(x, У - nD, Xq
N L
(6)
c(2)(x,y) = 1 + (/2'(Xq))2P2(Xq) х
n=0 0
(7)
хС2,01^,y - nD, yo)dxQ.
Здесь N связано с числом двойников в полисинтетическом двойнике. Так как суммирование ведется от нуля, рассчитываемое количество двойников будет на единицу больше, чем N.
Для использования (6) и (7) в (4) и (5) необходимо произвести замену у ^ у - пО.
На рис. 1 схематически представлен полисинтетический двойник, для которого справедливы соотношения (6) и (7). Пусть направления краевой и винтовой составляющих вектора Бюргерса двойникующих дислокаций по отношению к полисинтетическому двойнику будут такими, как это показано на рис. 1. В случае прямолинейных границ двойников в полисинтетическом двойнике имеют место следующие соотношения:
/1( XQ) = H\1 -^1;
2
L
f2( X0) = - ^ l1 - ^
(8)
(9)
При p(x0) = p2(x0) = const = p для скалывающих напряжений с^ получим результат,
представленный на рис. 2. Принималось: L = = 100 мкм; H = 21 мкм; D = 31 мкм; N = 5. Для
удобства рассчитывалась приведенная величина с*ху (л; у) = с^(x, у)/РА, где А = 2ПА- v).
н
о
о-
ЛЫ /
/2(уо)
В
X
Рис. 1. Схематическое изображение полисинтетического двойника, находящегося вдали от поверхности кристалла
Из рис. 2 видно, что напряжения локализуются на границах двойников, входящих в состав полисинтетического двойника, и у их вершин. У вершин двойников напряжения меняют знак на противоположный. С удалением от вершин двойников в сторону их роста также имеется область локализации положительных сдвиговых напряжений.
Рис. 2. Конфигурация полей сдвиговых напряжений Схх у полисинтетического двойника, находящегося вдали от поверхности кристалла
Аналогичным образом можно рассчитать поля напряжений у полисинтетического двойника с различной формой двойниковых границ входящих в его состав двойников. Это осуществимо путем задания функций /1(хо ) и /2 (хо). Плотность же двойникующих дислокаций на двойниковых границах задается функЦиями рДхо) и Р2(хо).
В Ы В О Д
На основании дислокационной макроскопической модели предложена методика расчета полей напряжений у полисинтетического двойника, находящегося вдали от поверхности кристалла. Методика позволяет рассматривать напряженное состояние внутри двойников, не прибегая к приближению тонких двойников. В качестве примера приведен расчет скалывающих напряжений у полисинтетического двойника с двойниками, имеющими прямолинейные границы при равномерном распределении двойникующих дислокаций на границах.
Л И Т Е Р А Т У Р А
1. Остриков, О. М. Дислокационная макроскопическая модель клиновидного двойника / О. М. Остриков // Вестник ГГТУ им. П. О. Сухого. - 2ооб. - № 2. - С. 10-18.
2. Косевич, А. М. Дислокационная теория упругого двойникования кристаллов / А. М. Косевич, В. С. Бойко // Успехи физических наук. - 1971. - Т. 104, № 2. -С. 101 -255.
3. Остриков, О. М. Форма клиновидных двойников в локально деформируемых ионноимплантированных монокристаллах висмута / О. М. Остриков // Изв. высш. учеб. заведений. Черная металлургия. - 2ооб. - № 9. -С. 5-7.
4. Остриков, О. М. Использование полипараксилиле-новых тонких пленок при исследовании пластической деформации монокристаллов висмута/ О. М. Остриков // Прикладная механика и техническая физика. - 2ооб. -Т. 47, № 4. - С. 162-166.
5. Остриков, О. М. Некоторые особенности формы клиновидных двойников в монокристаллах висмута, деформированных сосредоточенной нагрузкой / О. М. Остриков // Физика металлов и металловедение. - 2ооо. -Т. 90, № 1. - С. 91-95.
6. Классен-Неклюдова, М. В. Механическое двойни-кование кристаллов / М. В. Классен-Неклюдова. - М.: АН СССР, 1960. - 266 с.
7. Остриков, О. М. Исследование механического двойникования монокристаллов сурьмы методом наноин-дентирования / О. М. Остриков, С. Н. Дуб // Инженерно-физический журнал. - 2оо3. - Т. 76, № 1. - С. 170-172.
Поступила 30.10.2007
У
Ь
кр
Ь
ь