Научная статья на тему 'Неоднородности полей деформаций в зернах поликристаллических материалов и задача Эшелби'

Неоднородности полей деформаций в зернах поликристаллических материалов и задача Эшелби Текст научной статьи по специальности «Физика»

CC BY
314
63
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕОДНОРОДНЫЕ ДЕФОРМАЦИИ / ПОЛИКРИСТАЛЛЫ / ИНТЕГРАЛЬНЫЕ УРАВНЕНИЯ / INHOMOGENEOUS STRAINS / POLYCRYSTALS / INTEGRAL EQUATIONS

Аннотация научной статьи по физике, автор научной работы — Ташкинов А.А., Шавшуков В.Е.

Предложен и реализован способ вычисления неоднородных полей деформаций в зернах поликристаллических материалов. Вычисления основаны на разработанном ранее методе решения краевой задачи механики неоднородных поликристаллических тел с помощью оригинального варианта теории возмущений, основанного на аналогиях c квантовой теорией поля. Краевая задача для неоднородных полей деформирования в дифференциальной форме преобразуется в интегральное уравнение для тензора деформаций. Решение интегрального уравнения строится в виде ряда по интенсивности взаимодействия деформаций. Это позволяет интерпретировать неоднородную деформацию в какой-либо точке зерна как суперпозицию макродеформации, обусловленной граничными условиями, и двух составляющих, обусловленных внутризеренным и межзеренным взаимодействием. Показано, что в нетекстурированных поликристаллах, несмотря на дальнодействующий характер упругого взаимодействия, при оценке влияния межзеренного взаимодействия на неоднородность деформаций в выделенном зерне можно ограничиться учетом взаимодействия только с ближайшими и вторыми по удалению зернами-соседями. Вклады от взаимодействия с более далекими зернами взаимно компенсируют друг друга. Неоднородное в пределах одного зерна поле деформаций аппроксимируется кусочно-постоянной функцией. Для этого каждое зерно разбивается на большое число малых субзерен, в пределах которых поля деформаций принимаются однородными. Такая аппроксимация сводит интегральные уравнения для локальных деформаций к линейным алгебраическим, которые решаются численно. Применение метода к классической задаче вычисления деформаций в сферическом включении, погруженном в неограниченную матрицу, дает решение Эшелби. На модельных поликристаллах цинка выполнена численная оценка неоднородных деформаций. Вблизи границ сферического зерна экстремальные значения деформаций, обусловленные межзеренным взаимодействием, на 30 % превосходят средние. В материалах с более низкой упругой симметрией зерен концентрация деформаций существенно выше.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

INHOMOGENEITIES IN GRAINS OF POLYCRYSTALLINE MATERIALS AND ESHELBY PROBLEM

The paper presents the method aimed at calculating inhomogeneous strain fields in grains of polycrystalline materials. The calculations are based on the earlier developed method of solving boundary values problem for inhomogeneous polycrystalline bodies by means of the original perturbation theory variant based on analogies with the quantum fields theory. The boundary value problem for inhomogeneous strain fields in a differential form transforms into the integral equation for strains tensor. The solution of the integral equation is formed as a series upon the intensity of strains interaction. This allows interpreting inhomogeneous strain at any point in a grain as a superposition of macrostrain, caused by boundary conditions and two components conditioned by intragrain and intergrain interaction. It is shown that in untextured polycrystals, despite the long range type of elastic interaction, one can take into account the interaction only with the nearest and second neighbor grains to evaluate how the intergrain interaction influences the inhomogeneity in the given grain. The contributions of interactions with farther grains mutually annihilate each other. The strain field inhomogeneous within one grain is approximated by the step-wise constant function. For that, each grain is divided into a great quantity of small subgrains, where subgrain strain fields are supposed to be homogeneous. This approximation reduces the integral equations for local strains into linear algebraic ones, which are solved numerically. The application of this method to a classical problem related to calculating strains in a spherical inclusion embedded into the infinite matrix gives Eshelby solution. The numerical evaluation of strain inhomogeneities is made using model zinc polycrystals. Close to boundaries in spherical grains the extreme strain values, caused by intergrain interaction, surpass mean strain values by 30 percent. The strains concentration is much higher in materials with a lower elastic symmetry of grains.

Текст научной работы на тему «Неоднородности полей деформаций в зернах поликристаллических материалов и задача Эшелби»

Ташкинов А.А., Шавшуков В.Е. Неоднородности полей деформаций в зернах поликристаллических материалов и задача Эшелби // Вестник Пермского национального исследовательского университета. Механика. - 2018. - № 1. - С. 58-72. DOI 10.15593/perm.mech/2018.1.05

Tashkinov A.A., Shavshukov V.E. Inhomogeneities in grains of poly crystalline materials and Eshelby problem. PNRPU Mechanics Bulletin, 2018, no. 1, pp. 58-72 DOI 10.15593/perm.mech/2018.1.05

ВЕСТНИК ПНИПУ. МЕХАНИКА № 1,2018 PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/ini7

DOI 10.15593/perm.mech/2018.1.05 УДК 539.3

НЕОДНОРОДНОСТИ ПОЛЕЙ ДЕФОРМАЦИЙ В ЗЕРНАХ ПОЛИКРИСТАЛЛИЧЕСКИХ МАТЕРИАЛОВ И ЗАДАЧА ЭШЕЛБИ

А.А. Ташкинов, В.Е. Шавшуков

Пермский национальный исследовательский политехнический университет, Пермь, Россия

О СТАТЬЕ

АННОТАЦИЯ

Получена: 11 декабря 2017 г. Принята: 23 января 2018 г. Опубликована: 30 марта 2018 г.

Ключевые слова:

неоднородные деформации, поликристаллы, интегральные уравнения.

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

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

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

© ПНИПУ

© Ташкинов Анатолий Александрович - доктор физико-математических наук, профессор, e-mail: tash@pstu.ru Шавшуков Вячеслав Евгеньевич - кандидат физико-математических наук, доцент, e-mail: shavshukov@pstu.ru

Anatolii A. Tashkinov - Doctor of Physical and Mathematical Sciences, Professor, e-mail: tash@pstu.ru Viacheslav E. Shavshukov - CSc in Physical and Mathematical Sciences, Assistant Professor, e-mail: shavshukov@pstu.ru

Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

INHOMOGENEITIES IN GRAINS OF POLYCRYSTALLINE MATERIALS AND ESHELBY PROBLEM

A.A. Tashkinov, V.E. Shavshukov

Perm National Research Polytechnic University, Perm, Russian Federation

ARTICLE INFO ABSTRACT

The paper presents the method aimed at calculating inhomogeneous strain fields in grains of pol-ycrystalline materials. The calculations are based on the earlier developed method of solving boundary values problem for inhomogeneous polycrystalline bodies by means of the original perturbation theory variant based on analogies with the quantum fields theory. The boundary value problem for inhomoge-neous strain fields in a differential form transforms into the integral equation for strains tensor. The solution of the integral equation is formed as a series upon the intensity of strains interaction. This allows interpreting inhomogeneous strain at any point in a grain as a superposition of macrostrain, caused by boundary conditions and two components conditioned by intragrain and intergrain interaction. It is shown that in untextured polycrystals, despite the long range type of elastic interaction, one can take into account the interaction only with the nearest and second neighbor grains to evaluate how the intergrain interaction influences the inhomogeneity in the given grain. The contributions of interactions with farther grains mutually annihilate each other.

The strain field inhomogeneous within one grain is approximated by the step-wise constant function. For that, each grain is divided into a great quantity of small subgrains, where subgrain strain fields are supposed to be homogeneous. This approximation reduces the integral equations for local strains into linear algebraic ones, which are solved numerically. The application of this method to a classical problem related to calculating strains in a spherical inclusion embedded into the infinite matrix gives Eshelby solution.

The numerical evaluation of strain inhomogeneities is made using model zinc polycrystals. Close to boundaries in spherical grains the extreme strain values, caused by intergrain interaction, surpass mean strain values by 30 percent. The strains concentration is much higher in materials with a lower elastic symmetry of grains.

© PNRPU

Received: 11 December 2017 Accepted: 23 January 2018 Published: 30 March 2018

Keywords:

inhomogeneous strains, polycrystals, integral equations.

Введение

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

Проблема неоднородности полей деформаций и напряжений в поликристаллических материалах имеет давнюю историю. Первые модели поликристаллов, с успехом применяемые и в настоящее время, - модели Фойгта [5] и Ройсса [6]. В модели Фойгта деформации принимаются однородными в пределах всего поликристаллического тела, а напряжения однородны в пределах зерна, но испытывают скачок при переходе от зерна к зерну. В модели Ройсса напряжения однородны по всему поликристаллическому телу, а деформации однородны в пределах зерна, но скачкообразно изменяются при переходе границ зерен. Более поздние модели, такие как корреляционное, сингу-

лярное и самосогласованное приближения в теории неоднородных материалов [7], уточняют модели Фойгта и Ройсса в количественном отношении, но оставляют неизменным положение об однородности полей внутри индивидуальных зерен. Эти модели дают достаточно точные оценки средних деформаций и их дисперсий в отдельных зернах. Поэтому они приводят к хорошим результатам в вычислении эффективных упругих свойств материалов, связывающих средние напряжения и деформации, и успешно применяются во многих прикладных вопросах механики. Однако даже из общих соображений уравновешенности напряжений на границах зерен и совместности деформаций следует [8], что поля напряжений и деформаций должны быть переменны внутри зерен и необходима количественная оценка степени их неоднородности, а также влияния на эту неоднородность различных параметров микроструктуры поликристаллов.

Неоднородность деформаций в поликристаллах в рамках континуальной механики деформируемого твердого тела можно рассматривать на нескольких масштабных уровнях - от макроскопического (масштаб размеров тела) до микроскопического (масштаб размеров дислокаций). Объединение моделей поведения материалов на разных масштабах составляет суть многоуровнего подхода в механике материалов [9-14]. Многие явления хорошо описываются двухуровневыми моделями, оперирующими только с макро- и мезодеформациями. Настоящая работа посвящена исследованию неоднородностей мезодеформа-ций в пределах индивидуальных зерен, влияния окру-

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

Гипотетически исчерпывающее описание неоднородных деформаций в поликристаллах дает метод прямого компьютерного моделирования. Метод широко применяется в современном вычислительном материаловедении [15, 16]. Основным дискуссионным моментом в этом подходе является размер компьютерной модели материала (измеряемый количеством структурных элементов-зерен в модели), позволяющий адекватно описывать механическое поведение макроскопических образцов поликристаллических материалов, и статистическая достаточность метода. Доступные для реализации на современных компьютерах прямые модели содержат несколько сотен зерен. Реальные макрообразцы содержат миллионы зерен. Размер компьютерной модели должен быть во всяком случае намного больше размеров представительного элемента объема поликристалла. Дискуссия о размерах представительного элемента объема продолжается с момента введения концепции Р. Хиллом [17] до настоящего времени [18, 19]. Для вычисления эффективных упругих модулей считается достаточным представительный объем в несколько десятков зерен. Однако для анализа локальных явлений, для которых существен учет неоднородностей деформаций в отдельных зернах, требуемый размер представительного элемента существенно возрастает. В работе [19] отмечается, что в этом случае адекватный размер представительного элемента должен быть в 8-10 раз больше. Соответственно на порядок и более возрастает требуемый размер копьютерной модели, что не доступно для современных компьютеров. Поэтому актуальны альтернативные методы вычисления неоднородных полей деформаций в зернах. В настоящей работе оценки неоднородности упругих деформаций в зернах сделаны с помощью разработанного ранее метода решения краевых задач механики поликристаллов, основанного на аналогиях с квантовой теорией поля [20, 21].

Обычно принимается, что неоднородность полей внутри зерна невелика и ей пренебрегают. Впервые эту гипотезу применил И.М. Лифшиц [22] при вычислении эффективных модулей упругости поликристаллов со слабо анизотропными зернами. В дальнейшем гипотеза однородности получила некоторую теоретическую поддержку в работах Эшелби [23], показавшего, что в сферическом (и в более общем случае - эллипсоидном) зерне, погруженном в однородную среду, поле деформаций однородно. Аналогичные результаты были независимо получены М.А. Кривоглазом [24] и применены для вычисления упругих модулей твердой смеси, где форма отдельных зерен принималась сферической. Несмотря на геометрическую невозможность упаковать сферические зерна в сплошную среду, модель сферических зерен широко применяется в вычислительном материаловедении [25]. Гипотеза однородности в пределах зерен хорошо работает в задачах, где требуется найти только макро-

скопические поля деформирования, то есть осредненные по некоторому представительному объему материала.

Однако даже в сферическом зерне, погруженном в неоднородную среду окружающих зерен, деформации будут неоднородны. Это обусловлено упругим взаимодействием зерен.

1. Теоретико-полевой метод нахождения полей деформаций в зернах поликристалла

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

и1 (?) = (и,> + и;(г), е,(г) = <е,> + е,(г),

С« (?) = <с„ >+ст (г),

после чего уравнения равновесия в отсутствие объемных сил

дх,

Ст (Г)

дик(г)

дх,

= 0

принимают вид

(С,, >

д \ (г)

дх, дх1

дх,

Ск (г)

дик(г)

дх,

(1)

(2)

где (С > - осредненный по объему тела тензор модулей УпPУгости, а Ст (г) = С,,к1(г) -<Сцк1 > - его флуктуации. Решение (2) представляют в виде

и (г ) = и* (г ) + и; (г ), (3)

где и*(г) - решение уравнения (2) с нулевой правой

частью, то есть решение соответствующей задачи для однородного тела с осредненными модулями. Граничные условия выбираем в перемещениях

и (г)|г = и,(0)(?)|г

(4)

где и,(0)(г)|г - заданная функция на поверхности тела Г .

В силу сказанного выше все функции считаются непрерывными и дифференцируемыми. После этого осуществляется переход к интегральному уравнению для деформаций, впервые предложенный И.А.Куниным [26] и впоследствии ставший общепринятым [7, 27]

е,(?) = £*.(?) +1^(г -Г)С'т(гХп(г-), (5)

V

где 8,,Ы = ^л,Л + ) - тензоР ГРина сРеДы с осРеД-

ненным по объему тела тензором модулей упругости (С ,ы); Ол ^ - вторая производная тензора Кельвина-

Сомильяны; е*. (г) - решение краевой задачи для однородного тела с осредненными модулями (в общем случае неоднородное). Тензор Грина имеет известный вид [7]

g цы(r) = ■

1

3<Ц>

+ (1 - 5 X) DiJU

(1 -X) Vki +

fw(r)

8(r) +-

<Ц>

(6)

3( K )+(ц) й , K)

где x = 3^к)—) - материальный параметр; (K)

и (ц) - осредненные модули объемного сжатия и сдвига поликристалла; 8;j - символ Кронекера; VijU и DijU -шаровая и девиаторная части единичного тензора

1 .и = ^С8и8ji +8Й 8jk X а функция j (г) зависит от ? следующим образом:

j(r) = Г3 {""i8 .и + "j"i 8ik)"

8л |r|

" 2 j + X {{ + hw +1. ) " 3X -{ ""j 8ki + ""и 8 ji + "t",8 jk + "j" 8n +i

+ njni8л + ntni8,j - 5«,^«)],

(7)

где = ^ - направляющие косинусы радиуса-вектора ? . г

Если в краевой задаче выбрать граничные условия в перемещениях, соответствующие макроскопически однородной глобальной деформаций ц (г)|г = е*jxj, где

е*. - макродеформация, то первый член в правой части (6) - решение соответствующей однородной задачи -будет равен просто константе е*.. В дальнейшем будет

рассматриваться именно этот случай.

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

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

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

и нулю вне зерна (Ы - полное число зерен в поликристаллическом теле):

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

е. (г) = £Ч (г^ (г), С„т„ (г) = £ (г). (г) . (8)

Тогда уравнение для глобальных полей (5) сводится к системам уравнений для локальных взаимодействующих полей е^г) [21],

е^г;) = е* + |g¡jkl(г, -фск®е!Р(г:)+

+ 11 (г - г )сжи)(гл'), (9)

где ст=с<т -(сышп) - постоянный тензор в пределах ^-го зерна.

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

е«Чгр) = е<°>« >(г,) + е Г «(г,) + + е(2)(; )(г;) + е(3'|(; ;(г,) +..., (10)

получим следующую бесконечную систему уравнений для поправок:

е^)- | .(г, -ет0)С;)(г') =е*, е(1)С,)(г|)^ gjU (г, - г^ст =

=11 (г, - г; )с^ етгчг;),

)- ет^е?)=

= У Глр е (г - г )с;(;) е(1)(;) (г;)

} .кАЧ к1тп тп V;''

(^)ёг,; - гТ)^ етп(?)(г:) = =Е1 <кш (г, - осгтп е т2п)(л)(г;)

(11)

Л*? о

и т.д.

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

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

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

2. Нулевое приближение по межзеренному взаимодействию. Задача Эшелби

Исходное интегральное уравнение (9) для деформаций в ^-м зерне в нулевом приближении по межзерен-ному взаимодействию принимает вид

,(0)«)

(п ) = S* +

+ J drl g1]kl(Г, - S^)- (12)

Полевую переменную ) в этом уравнении

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

вГЧ?; )=S^ a (Г ^mr ).

(13)

При этом объем ^-го зерна равен сумме объемов суб-

n

зерен ю, = Sroa, где n - число субзерен в зерне. После

a=1

подстановки индикаторных разложений в интегральное уравнение (12) получим систему интегральных уравнений для деформаций s(0)(ç)(a)(r, ) в субзернах ^-го зерна, которая по структуре аналогична вышеприведенной системе (9) для деформаций в зернах без разбиения на субзерна,

вГ)(a)(ra ) = £* + J dr: gji(ra -T^ZS™® +

+ SJ Î,gm (?a - КСЩ(14)

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

,(0)(,)(a)(0)

(Г ) -

- [dr ' g (Г - Г ')C'(,) S(0)(,)(a)(0)(Î') = S*

J" 'a&ijklVa 'a^klmn^mn У,> bij'

,(0)(,)(a)(1)

(Г ) -

-J dr' gl]U (ra - r^cz smo)(,)(a)(1)(ra)= = S J dr g (r - r')C s(0)(,)(a)(0)(î')

J b ijkl a b klmn mn b

,(0)(,)(a)(2)

(Г ) -

- J dr' gjl (ra - OC'm sl0n)(Ç)(a)(2)(r:) = = SJ d?g„ (ra - VClmn S m0n(Ç)(i)(1)(r'),

(15)

s(0)(,)(a)(3) ij

(ra ) -

- J dr' gk(î: -r'a)C'mn s00n)(ç)(a)(3)(ra') =

= SJ ^ (Îa - r')C'Umn S 00n)(Ç)(i)(2)(r'},

и т.д.,

где с, Ь = 1, 2, ... п.

Полное решение исходного уравнения (15) равно сумме всех поправок,

s

(0)(, )(a )

(Îa ) =S

_„(0)(, )(a )(0)

(Îa ) +s;

(0)(,)( a )(1)

(Îa ) +

+ S

(0)(,)(a )(2)

(Îa ) + S

(0) (,)(a)(3)

(Îa ) + ■

(16)

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

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

a=1

,(0)(У(а)(0) = е* "Ы у'

|Т' ук, В(к, е,

Г Т - В(ааХ1 Р(0Ж)(а)(-) = V В("Ь)К')р С0)С,)СЛ)С0)

_ ук1 к1 ^ к1 ~ ¿^^уЫ к1

Ь *а

Г Т - Я(аа)(,) 1 _ С0)С;)Са)С2) = V В(аЬ)(?)р(0)«ХЬ)(1)

_ ук1 ук1 ^ к1 £^иук1 к,

Ь *а

Г Т - В(аа)(О 1 р С0)С;)Сa)С3) = V Я(аЬ)«) 10Ш,Ь)(2)

\Jiikl ыук1 ^ к1 ук1 к,

Ь* а

и т.д., а = 1,2, ... п,

(17)

где

В,

(аа)(,) =

| ^ЗуИ (га - Г'а) • с;

к,тп '

В.

(аЬ)Й) _

I ^и (К - г,') • ск^тп (а * Ь ) (18)

есть постоянные тензорные коэффициенты, описывающие взаимодействие деформаций внутри одного субзерна и между субзернами соответственно. Всего (17) содержит п • N систем линейных уравнений размерностью 6^6 для шести компонент деформаций во всех субзернах всего поликристалла.

Первые уравнения в (17) для поправок нулевого порядка для всех субзерен независимы друг от друга, так что поправки нулевого порядка определяются параметрами только рассматриваемого субзерна.

Коэффициенты В^ безразмерны и уменьшаются

с увеличением расстояния между субзернами. Они вычисляются численным интегрированием для каждой конкретной микроструктуры поликристаллического материала. Вычисления показывают, что даже для соседних субзерен их величина по модулю порядка 0,02 и меньше. Поэтому итерационная процедура (17) по взаимодействию быстро сходится. Ограничимся поправками первого порядка по межсубзеренному взаимодействию, то есть примем

^0)©^) ,^(0X0^X0) + р(0)0(а)(1) . . .

(19)

Тогда в бесконечной цепочке (18) останется два первых уравнения

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Г . - В.аа) >2*0(аХ0) =е*, (20)

Г Т - В(аа) 1 р(0)й)Сa)С1) = V С0)С;)СЬ)С0)

^ук, иук1 ^ к1 ¿^ ук1 к1

Ь *а

Выберем форму и размеры всех субзерен одинаковыми. В этом случае решения нулевого порядка для всех субзерен будут совпадать е(;0)(О(а)(0) = е(0Х®(ЬХ0),

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

Г т - в(аа) 1е

| ук1 ук1 JЧ

С0)c;)Сa)С0) = к1 ьу

Т - в(аа) 1 ^«Ха)(1) = В(a)еС0)С;)Са)(0)

_ укI ук, J кI укI кI >

содержащих неизвестные только для одного субзерна и новый коэффициент взаимодействия выделенного субзерна со всем остальным объемом зерна

Вт, - I (га - г;) • с'кк

(22)

где область интегрирования, обозначенная как (ю - а), равна (ю, -юа).

Коэффициенты в уравнениях (21) определяются параметрами а-го субзерна, его положением в зерне и параметрами зерна в целом. Особенностью (21) является то, что дискретным образом выступает только а-е субзерно. Остальные субзерна входят интегрально. Разбиение остальной части зерна на субзерна становится несущественным, так как все остальные субзерна включены в интеграл (22) в качестве единой области интегрирования. Это позволяет выбрать объем а-го субзерна в виде шара. Ниже будет показано, что для шарового субзерна коэффициенты В^ не зависят от его

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

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

Коэффициент взаимодействия В^ для сферического субзерна вычисляется аналитически в общем ви-

де. В интеграле (18) для В1

( аа)

В':""'1 - I (¥' е.к1 (г - г') • О,,

утп } а&ук ^ а а' кт

начало координат поместим в центр субзерна (га = 0). Подставляя выражение для тензора Грина и производя интегрирование по угловым координатам, убеждаемся, что интегралы от функции . (г) в (7)

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

В (аа) =

утп

3(ц)

(1 -Х). + (1 -тХ)

укI

с '

к,тп'

(23)

Коэффициенты взаимодействия Вутп какого-либо

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

ю-а

Грина через производные тензора Кельвина-Сомилианы, перепишем выражение (22) в следующем виде:

B (a) = 1 — iJmn 2 cX(a)

j dr'Gkki(ra -r') • Cm +

J l,

1 д

+ 2 аХ^" j drGJkI (r - r) • Ck'mn -

-j dA-

3<Ц>

(1 -X)V,« +

+(1 - 5 X) DjU

5ft -Г')^ Q

(24)

где 0,к1 и 0,1 - первые производные тензора Кельви-

на-Сомилианы [7],

0к^=- ^х

х [ ппп + 28,кп1 - Х(п,8и + П 8а + П 8 л)]. (25)

Интегрирование в первых двух интегралах (24) производится по всей области ^-го зерна, а в третьем интеграле только по области а-го субзерна. Третий интеграл устраняет слагаемое с 5-функцией в (6), которое возникает при выполнении дифференцирования в (24), но не

входит в исходное выражение (22) для Б^,, так как в нем всюду га Ф г' и 5-функция обращается в ноль.

Рассмотрим первый интеграл в (24) и сделаем в нем замену переменных интегрирования га - г' = г ,

j dr'Gk j (Fa - r') = -j drGik j (r ).

(26)

Он вычисляется методом, использованным Эшелби для подобного интеграла в работе 1957 г. [23]. Запишем его в сферических координатах г, 6, ф. Начало координат

поместим в центр с-го субзерна, то есть га = 0. Область интегрирования по г есть вся шаровая область на рис. 1.

Рис. 1. Сферическое зерно и область интегрирования в (26) Fig. 1. Spherical grain and integration area

На рисунке L = (L1,L2,L3) - радиус-вектор, проведенный из начала координат в центр зерна. Верхний

предел интегрирования по радиусу г (п), зависящий от

направления П (от углов 6, ф), есть точка пересечения

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

- Г \2 ( - J V / _ т \2 XX Li I I Хъ L I I Xo L-,

R

R

R

= 1.

Легко видеть, что верхний предел интегрирования находится из квадратного уравнения

г2 - 2(пЬ)г - Я2 + Ь2 = 0.

Из условия положительности г > 0 оставляем один корень уравнения

r (П) = (nL) + yl(nL)2 + R2 - L2.

Производя в (26) сначала интегрирование по радиусу, получим

j dFGikl(r) =—— jdо[(nL) + V(nLf + R2 -L2 l, [

x [n^n, + 25лп, -x(n5и + nk5a + n,5л^

где о - телесный угол, Со = 8ш6с!6с!ф, а интегрирование производится по всему телесному углу вокруг начала координат. Слагаемое с квадратным корнем исчезает в силу нечетности функции в правых квадратных скобках, и остаются только табличные интегралы от произведений направляющих косинусов

| !г01к11(г ) Iс(пЬ)[3Хп,пкп1 +

»-а ' (27)

+ 28,кп1 - Х(п, 8к1 + пк 8,1 + п, 8,к)].

Теперь в (27) возвращаемся к переменной гс и, замечая, что

д д

dxja) dLJ

получаем

jd/Gik I (Fa ) =-¿^ d ^

dx(a >

j i,

<[3ХпЛП +25ikni -X(n5ki + nk5ii + n5ik)] = =__

=-<ц>

1 - 3 x }vj -15 xD,kji

(28)

аналогично для второго слагаемого в (24)

д

r - ' J J 1 r

I dr Gjkl(ra -r ) =--I dоч

J Jk,i a 8rc<^>J

dx(a 5 i i,

([3Xnjnkni + 25Jkni -X(nj5ki + nk5Ji + n5Jk)] = 1

<ц>

1 -3AVm -15XDiiJk

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(29)

1

д

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

В(а> = 0.

утп

(30)

Из (21) следует, что поправки 1-го порядка по межсуб-зеренному взаимодействию обращаются в ноль, и далее из (17) следует, что все поправки высших порядков также обращаются в ноль. Деформации во всех субзернах зерна оказываются одинаковыми е. = е(j0)(a) = и определяются из одной системы линейных уравнений

Тутп + 3(ц)

+ (1 - 5 х) оук,

(1 -Х). +

с' [е(а) = е*..

катп | тп у

(31)

Таким образом в сферическом зерне в нулевом приближении по межзеренному взаимодействию деформации однородны. При переходе к соседнему зерну деформации испытывают скачок вследствие скачкообразной смены в (31) сктп - сШп -(ск1тп) , ВШИИШ^

поворотом кристаллографических осей зерен.

Нулевое приближение по межзеренному взаимодействию имеет особое значение в связи с задачей Эшелби. Применим изложенный формализм к одиночному включению сферической формы, погруженному в неограниченную матрицу. В качестве включения возьмем зерно с тензором модулей упругости скС). Все остальные (т|-е) зерна выберем одинаковыми с тензором модулей упругости мат-

С( таЗТ )

.к , так что они составят неограниченную однородную матрицу. Тогда нужно положить (сук1) = Ст^ ,

так как объемная доля одиночного включения ничтожно мала. Легко видеть, что в правых частях точных уравнений в (11), начиная со второго, остаются нули, так как в них

. - с; - с,) = стт) - суг " = 0. После

"уй

-'ук

ук

"ук,

этого, сло-

жив все уравнения (11), получим, что поле деформаций

внутри включения е('пс1) (г) определяется из одного точно-

го уравнения

^(г) - I ёг; (г - г') х

х(с<т) - скта,г))е|>с1)(г/) = е*..

V катп катп ' тп V , у

(32)

Таким образом, уравнение (32) есть уравнение задачи Эшелби для поля деформаций внутри сферического включения в матрице.

Разобьем включение на большое число малых суб-

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

все деформации в субвключениях будут одинаковыми,

а значит, е(.'пе) (г) однородна внутри включения. Поэтому уравнение (32) сводится к линейному алгебраическому уравнению, совпадающему по форме с (31),

1уп + 3(ц)

(1 -X). +

(с(пс,) - с(та1г)) {Лтс,) = * /"З-ЗЧ

у к,тп к,тп ) Гьтп Ьу • (33)

+(1 - 5 х) в,

Уравнение (33) является точным для одиночного включения. Его легко решить численно для любых типов анизотропии включения и матрицы и произвольной взаимной ориентации их кристаллографических осей.

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

с) = 3 К

ук,

тс^уЫ + 2ц пс^уЫ,

(34)

с (та!г) = V + 2и В (35)

ук -"^тайг* ук, ^ГтаГг^уМ' V3-V

где К и ц - модули объемного сжатия и сдвига включения (шс1) и матрицы (шаП).

Подставляя (34), (35) в (33), получим

^^ + 3Ц

3(1 -Х)(Кш - Кта,г Уш +

+ 2(1 ^X)(Цincl таг )ВуЫ

пС та&' ук

сincl) _ 0*

и его решение

'пс,) _

1

'к,

(36)

1 +

(1 -Х)(КЩ -Кта,г )

Цта!г

V.к, +

2 Вуы

2(1 -ТХ)(Цпс, таг )

1+

3Цт

(37)

Выражение (37) совпадает с решением Эшелби для деформаций внутри изотропного сферического включения [23, 27].

Уравнение (33), решаемое численно, является обобщением классической задачи Эшелби на случай анизотропного сферического включения в анизотропной матрице.

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

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

+

е

к,

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

Рис. 2. Модельный поликристалл со сферическими зернами Fig. 2. Model polycrystal with spherical grains

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

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

Необходимо учесть взаимодействие всех субзерен в выделенном 4-м зерне и взаимодействие этих субзерен со всеми субзернами остальных "Л^-х зерен. В точном уравнении (9) разложим полевые переменные во всех зернах на переменные в субзернах,

) = (Г )е(Гх:)(Г)

(38)

^(Г, ) = e (Г, К

(е)(л) ij

(r,),

и получим точную систему интегральных уравнений для деформаций во всех субзернах всех зерен

еГЧга) -1 ск .(г -ОС^ЧО =

»с

= 8* + £ | !гь Ер, (г - гь )с<ц8тчп) +

ЬтЬ

+ X X I ^м(га -г.). (39)

Л*? «=1 т.

В (38) и (39) индексы с, Ь нумеруют субзерна в выделенном ^-м зерне, индекс . нумерует субзерна во всех остальных зернах поликристалла.

Теперь применяем к уравнению (39) теорию возмущений. Все межсубзеренные взаимодействия, то есть оба члена с суммами в правой части, принимаем за возмущение и представляем неизвестные переменные в виде рядов, суммирующих поправки всех порядков,

е

( аАе)«,л)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(r ) = е(0)(«■,)(?) +

+ е

(1)(аАе)«,л)

(r) + е

(2)(а,4,е)(?,л) j

(Г ) +.

(40)

Поправки находятся из бесконечной цепочки систем уравнений

8] ^ (гс ) - I ^с 8]к1 (гс - гс)Ск1^п8™ с {гс) =8*.,

»с

¿«(г) -1 сг; Е,]к1 (гв - ^тп ¿тсю=

»с

= X | (Т. -Ч)Ст8«\7Ь) +

_ n

ЬЕ S j dregw (Га _ re )Ct

Л*? e=1 ,

(Л)e(0)(e)(,) (Г )

klmn mn ^ e '

и т.д.

(41)

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

е(0)(а)(0 _ в(аа)(Эе(0)(а)(:) = * ij ijk kl ij'

е(1)(а)(:) _ в(аа))(?)е(1)(а)(:) =

ij ijkl mn

= S B

( аЬЖ)^(0)(Ьт ijkl

kl

-XX в:

л*? e=1

( аe)(л^(0)(e)(л) jkl

kl

. и т.д.,

(42)

где коэффициенты

в.

(аЬ)С) =

j drbg jkl (Г _ Г') • C

klmn

B (ae)(л) =

j drbgjk (Га _ Г,') • C

'(л)

klmn

(43)

описывают взаимодействие субзерен в одном зерне и в разных зернах соответственно.

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

а=1

и

e=1

уравнениями в (42) и (30). Кроме того, для сферических зерен решения нулевого порядка е(г0)(а)® однородны, поэтому в них можно опустить индекс субзерна е(о)(а)(э - е(Ш), а член с одинарной суммой во втором

уравнении (42) обращается в нуль. Остается система из двух уравнений

р(0)(э _ в(»»)(9р(0)(9 = р* а . к у >

р(1Хя)«) - о(ои)«)p(l)(a)(0 = Y B(ол)р(0)(л)

ij ijki m" / . ijki ki '

Л*?

а новый коэффициент

B (ал) =

ijm"

j <8ijk,

(Г - г' с,(л)

\'a 'л^ kim"

(44)

(45)

определяет интенсивность взаимодействия субзерна в ^-м зерне с п-м зерном в целом. Этот коэффициент зависит от положения субзерна в £-м зерне, от взаимного расположения ^-го и п-го зерен и от ориентации кристаллографических осей этих двух зерен. Таким образом вновь взаимодействие с субзернами других зерен входит в теорию интегрально, то есть фактически разбивать на субзерна достаточно только рассматриваемое конкретное зерно. Неоднородность деформаций в зерне обусловливается только межзеренным взаимодействием.

Заметим, что е(,0Хт|) в правой части второго уравне-

ния (43) определяется из уравнения р(°)(л) - bi

э(аа)(л)р(0)(л) _

_р*. , в котором коэффициенты Bjaif(Л) ся выражением B^ - j dr'agijkl х (ra - DC

определяют-

(л)

Umn > то естЬ

по форме совпадают с ранее найденными коэффициентами (18).

Из линейности системы (44) следует, что межзерен-ное взаимодействие, описываемое членом с суммой, дает аддитивную добавку к решению предыдущего параграфа. Эта добавка является стохастической. Ясно, что она изменяется в некоторых конечных пределах, и для нее можно построить функцию распределения методом, изложенным в [28]. Для оценки величины неоднородности деформаций, вызванной межзеренным взаимодействием, приведем здесь только маргинальные значения этих добавок на ансамбле всевозможных конфигураций окружения выделенного зерна. Число возможных конфигураций для макроскопического поликристаллического образца чрезвычайно велико. Маргинальные значения для решений (44) реализуются только для некоторого небольшого набора конфигураций из всего ансамбля. Поэтому для поиска маргинальных значений пришлось бы решать систему (44) для очень большой статистической выборки конфигураций ансамбля. Однако аддитивность вкладов взаимодействия с отдельными зернами кардинально сокращает объем вычислительной работы, делая ее доступной для однопроцессорных компьютеров.

Следуя работе [29], разбиваем сумму во втором уравнении (44) по всем п-м зернам на сумму по ближайшим зернам-соседям (1-я координационная сфера содержит 26 зерен), следующим по удалению зернам (2-я координационная сфера содержит 98 зерен), 3-я координационная сфера (216 зерен) и т.д.

p(1)(a)(?) - в(aa)р(1)(a)«) _

ij ijki m"

ZB(aл)р(0)(л) + V B(ол)р(0)(л) + ijki ki ijki ki

Л*? Л*?

(46)

Вклады зерен всех координационных сфер аддитивны и могут быть рассмотрены отдельно. Поэтому для нахождения маргинального значения деформации в какой-либо заданной точке (субзерне) рассматриваемого зерна достаточно найти ориентацию кристаллографических осей, дающую максимальный вклад для каждого зерна-соседа по отдельности, а затем просуммировать все вклады. Для каждого зерна-соседа рассматривается статистическая выборка из 105-106 его случайных ориента-ций. Уравнения (46) для поправок первого порядка решаются численно для каждой реализации случайного окружения выделенного зерна из этой выборки. Стохас-тичность для модельного поликристалла сводится к случайной ориентации кристаллографических осей зерен. Случайная ориентация осей каждого зерна задается генератором случайных чисел и специальным распределением углов Эйлера, обеспечивающим равновероятность всех ориентаций [см., например, 30]. По массиву 105-106 полученных решений в каждой точке ^-го зерна выбираются маргинальные значения, обусловленные каждым зерном-соседом по отдельности. Затем маргинальные вклады от всех зерен складываются, что дает экстремальное значение деформации в рассматриваемой точке.

В качестве примера рассмотрим нетекстурированный поликристалл цинка, подвергнутый макроскопической деформации чистого сдвига. Зерна цинка имеют гексагональную упругую симметрию. Упругие свойства монокристалла зерна взяты из [31]. Кристаллографическую ось симметрии 6-го порядка выбранного (^-го) зерна направим вдоль оси Х3 глобальной системы координат поликристалла. Диаметр зерна принят равным единице.

xA k

/ 4 /______j— 3 ■ _ A

Рл 2

j . ¡i X2

— j- _

Рис. 3. Характерные точки в зерне Fig. 3. Characteristic points in grain

Рассмотрим макросдвиг в двух плоскостях XX и ХХ3, соответствующих в выбранной системе координат максимальной и минимальной жесткости на сдвиг выделенного зерна; е12 = 3 • 10 4 и е23 = 3 • 10 4.

Сначала найдем точки возможных максимумов деформаций. Рассмотрим характерные точки, показанные на рис. 3. Максимально удаленные от центра зерна точки находятся на расстоянии 0,45 от центра зерна (для наглядности рисунок сделан не в масштабе). Точки на границе зерна не рассматриваются, так как в упругой постановке в них возможны сингулярности. Предлагаемый метод работает в области, ограниченной штриховыми линиями.

В силу упомянутой аддитивности для нахождения координат точек максимумов достаточно рассмотреть влияние только 26 ближайших соседей. Учет взаимодействия с зернами более удаленных координационных сфер будет увеличивать значения максимальных деформаций, не меняя координаты точек максимумов. Максимальному (или минимальному) значению деформации в каждой точке на рис. 3 будет соответствовать вполне определенная ориентация кристаллографических осей соседних зерен, причем для разных точек такие «экстремальные» ориентации осей соседей будут разные.

Результаты вычислений максимальных положительных деформаций сдвига в выбранных точках зерна в12)(° (г) = е12)(° (г) + е12)(а)(° (г) 126 при макросдвиге

s*2 = 3 -10 4 в плоскости ХХ2 и sg*5)(r ) _s(u3,(y(r ) +

_о(0)(5),

+ s:

(1)(a)( 5)

(r) 126 при макросдвиге e23 = 3 -10 4 в плоскости

Х2Х3 и учете взаимодействия только с ближайшими соседями приведены соответственно в табл. 1 и табл. 2. Для информативности отдельно приведены поправки нулевого порядка и 1-го порядка по взаимодействию.

Таблица 1

Максимальные деформации б(2)(?) в характерных точках зерна (*104)

Table 1

Maximum strains s(2)(?) at characteristic points of grain (*104)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Поправка Точка

1 2 3 4 5 6 7

s(0)(5) 12 2,469 2,469 2,469 2,469 2,469 2,469 2,469

s(1)(a)©| 12 126 0,309 0,342 0,335 0,343 0,327 0,362 0,363

s(a)(5) _ s(0)(5) + 12 _ 12 "r + s(1)(a)(5)| 12 126 2,778 2,811 2,804 2,812 2,796 2,831 2,832

Для обоих случаев наибольшие деформации наблюдаются в точках вблизи границы зерна с абсолютными максимумами в точках 7 и 6 (зенит и надир). Однако различие максимальных значений в разных граничных точках порядка 1 %. Учитывая модельность задачи, можно сказать, что в любой точке зерна вблизи границы может наблюдаться максимально возможное значение деформации. Каждому такому максимальному значе-

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

Таблица 2

Максимальные деформации е2з)№) в характерных

точках зерна (*10 )

Table 2

Maximum strains s2a)(5) at characteristic points of grain (*104)

Поправка Точка

1 2 3 4 5 6 7

s(0)(5) 23 3,241 3,241 3,241 3,241 3,241 3,241 3,241

s(1)(a)(5)| 23 126 0,364 0,381 0,415 0,383 0,413 0,427 0,428

s(a)(5) _ s(0)(5) + 23 _ 23 + s(1)(a)(5)| 23 126 3,605 3,622 3,666 3,624 3,654 3,668 3,669

Для нахождения абсолютного максимума деформации в какой-либо точке нужно учесть взаимодействие с зернами 2-й и последующих координационных сфер. Рассмотрим деформацию в точке 7 на рис. 3. Максимальные вклады от взаимодействия с 98 зернами 2-й координационной сферы в этой точке для двух рассмотренных случаев деформирования, найденные путем раздельного поиска наибольших аддитивных добавок на выборке 105-106 случайных ориентаций всех зерен по

.,(!)(точка7)(0 I =

отдельности,

равны соответственно

(1)(точка7)(5) I _п <-><">« 1П-4

= 0,191 • 104 и е^™198 = 0,225 • 10 4. В обоих случаях вклад взаимодействия с 98 зернами 2-й координационной сферы составляет около 52 % от взаимодействия с 26 зернами 1-й координационной сферы. Это соответствует сделанному в работе [29] выводу о медленном убывании с расстоянием упругого взаимодействия в поликристаллах. Максимальные вклады от взаимодействия с 218 зернами 3-й координационной сферы и всех последующих также легко вычисляются, но для неупорядоченных нетекстурированных поликристаллов в этом нет практического смысла по следующим соображениям. Вычисленным экстремальным добав-

26

кам s

(1)(точка7)(5) I

и s

(1)(точка7)(5) I

у | 26 — -у 198 соответствуют

и 98 троек вполне определенных углов Эйлера для зерен 1-й и 2-й координационных сфер. Если обозначить вероятность экстремальной ориентации осей отдельного зерна-соседа через Р, то вероятность одновременной экстремальной ориентации осей зерен 1-й координационной сферы равна Р26, 1-й и 2-й координационных сфер - Р26Р98, 1-й, 2-й и 3-й координационных сфер -

26 98 218

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

12

98

для компоненты е23 в точке 7 и равен е™точка7)№) = = (3,241 + 0,428 + 0,225)•10-4 = 3,894•10-4. Это значение превосходит макродеформацию е23 = 3 • 10-4 на 29,8 %. Учитывая вышеизложенное, такой вывод относится к любой точке вблизи границы зерна. Таким образом, флуктуации деформаций в неупорядоченном поликристалле значительны и не могут полагаться малыми.

4. Оценка вероятностей реализации экстремальных конфигураций локальной микроструктуры

Произведем оценку вероятности реализации экстремальных значений деформаций в зерне. Эта вероятность определяется вероятностью реализации экстремальной конфигурации окружения выделенного зерна. Ограничимся анализом экстремальной конфигурации 26 ближайших зерен-соседей. Структура поликристалла принимается полностью неупорядоченной. В рассматриваемой модели это означает статистическую независимость ориентаций осей всех зерен. Однако вклады взаимодействия в неоднородность деформаций от разных соседей одной координационной сферы сильно различаются. В табл. 3 приведены вычисленные аддитивные вклады ближайших зерен-соседей в деформацию в точке 7 на рис. 3. Номера зерен-соседей нумеруются по рис. 2. Зерно № 1 - в зените, зерна № 2-8 - первый сверху пояс зерен, № 9-18 - зерна экваториального пояса, № 19-25 -нижний пояс и зерно № 26 - в надире.

Вклад зерна-соседа в зените (33,4 10-6) более чем в 30 раз превышает вклад зерна-соседа в надире и дает около четверти общего вклада в неоднородность. Суммарный вклад зерна в зените и шести (из семи) зерен верхнего пояса составляет более 60 % суммарного эффекта взаимодействия со всеми 26 соседями 1-й координационной сферы. Таким образом, семь ближайших

зерен дают основной вклад в появление экстремальных деформаций. Поэтому вероятность экстремальной конфигурации микроструктуры практически определяется вероятностью совместной экстремальной ориентации только этих зерен. Вычислительные эксперименты показали, что зависимость вклада в неоднородность от ориентации отдельного зерна-соседа вблизи максимума достаточно пологая. Вероятность попадания случайно ориентированного зерна в положение, дающее вклад 85 % и более от возможного абсолютного максимума для всех ближайших зерен-соседей колеблется около Р = 0,05. Вероятность реализации экстремальной локальной микроструктуры из п зерен вокруг какого-либо определенного зерна имеет порядок Рп, то есть весьма мала. Эта локальная микроструктура представляет собой кластер небольшого количества согласованно ориентированных зерен. Макроскопический образец поликристаллического материала содержит очень много зерен, поэтому вероятность появления в нем таких кластеров может составлять конечную величину. Оценим порядок такой вероятности. В поликристаллическом образце каждое зерно имеет случайное окружение, поэтому макроскопический образец, содержащий N зерен, представляет собой статистическую выборку из N случайных локальных микроструктур. Легко видеть, что вероятность появления в такой выборке кластера экстремальной микроструктуры равна

П = 1 - (1 - Рп).

В макроскопическом образце размерами 10*10x10 мм и средним диаметром зерна 20 мкм содержится N = 1,25-108 зерен. Для п = 7 искомая вероятность П = 0,093. Это заметная величина. Следовательно, при измерении некоторых механических свойств на партии из порядка 10 образцов весьма вероятно появление на единичных образцах больших отклонений от средних значений по партии, обусловленных не наличием каких-либо

Таблица 3

3

Вклады от взаимодействия с ближайшими соседями (*106)

Contributions from interactions with nearest neighbors (*10 )

Table

Зенит

№ = 1

33,4329

Верхний пояс, зерна № = 2-8

2 3 4 5 6 7 8

3,6359 9,1634 8,3775 8,7724 8,7264 8,3823 9,0074

Экваториальный пояс, зерна № = 9-18

9 10 11 12 13 14 15 16 17 18

2,9284 4,3017 4,1246 4,1107 4,2894 2,9289 4,3077 4,1071 4,1031 4,3017

Нижний пояс, зерна № = 19-25

19 20 21 22 23 24 25

1,5001 2,2656 1,8092 2,0778 2,0695 1,7133 2,2509

Надир N = 26 1,0598

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

Заключение

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

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

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

1. Antolovich S.D., Armstrong R.W. Plastic strain localization in metals: origins and consequences // Progress in Materials Science. -2014. - Vol. 59. - P. 1-160. DOI: org/l0.l0l6/j.pmatsci.20l3.06.00l

2. Pineau A., Benzerga A.A., Pardoen T. Failure of metals I: Brittle and ductile fracture // Acta Materialia. - 2016. - Vol. 107. -P. 424-483. DOI: org/l0.l0l6/j.actamat.20l5.l2.034

3. Failure of metals II: Fatigue / A. Pineau, D.L. McDowell, E.P. Busso, S.D. Antolovich // Acta Materialia. - 2016. -Vol. 107. - P. 484-507. DOI: org/l0.l0l6/j.actamat.20l5.05.050

4. Pineau A., Benzerga A.A., Pardoen T. Failure of metals III. Fracture and fatigue of nanostructured metallic materials //

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

Вклад взаимодействия зерен определяется введенными коэффициентами взаимодействия (18). В работе авторов [21] было показано, что в зависимости от типа анизотропии зерен величина коэффициентов взаимодействия изменяется более чем на порядок. Поэтому в поликристаллах с более низкой упругой симметрией зерен (каковыми являются, например, сложные многокомпонентные керамики) вклад взаимодействия в неоднородность деформаций может достигать 100 % и более от средних значений.

Прогнозирование возможных конфигураций экстремальных кластеров зерен может иметь практическое значение для относительно новой области материаловедения (годом рождения считается 1984-й [32]) - инженерии зерен (см. недавние обзоры [33, 34]), занимающейся поиском и формированием оптимальных микроструктур конструкционных материалов. Например, упрочнение металлов твердыми частицами-зернами (карбидами и т.п.) будет более эффективно, если окружающие зерна матричной фазы будут ориентированы так, что дают максимум напряжений в зернах упрочняющей фазы. Тогда по теореме о среднем напряжении основная нагрузка ляжет на упрочняющие зерна, что приведет к дополнительному упрочняющему эффекту. Предлагаемый нами метод позволяет прогнозировать такие оптимальные микроструктуры.

Благодарности

Работа выполнена при финансовой поддержке РФФИ (гранты №15-01-06912, 16-01-00682, 17-41-590433 р_а).

Acknowledgments

The work has been carried out with the financial support of the Russian Foundation for Basic Research (Grants No. 15-01-06912, 16-01-00682, 17-41-590433 р_а).

Acta Materialia. - 2016. - Vol. 107. - P. 508-544. DOI: org/10.1016/j.actamat.2015.07.049

5. Voight W. Lerbuch der Kristallphysuk. - Leipzig und Berlin:Teubner, 1928. - 978 p.

6. Reuss A., Berechnunug der Fliebgrenze von Mischkristallen auf Grund der plastizitatsbending fur einkristalle // Z. Angew. Math. und Mech. - 1929. - Vol. 9. - No. 1. - P. 49-56.

7. Шермергор Т.Д. Теория упругости микронеоднородных сред. - М.: Наука, 1977. - 400 с.

8. Штремель М.А. Прочность сплавов. Ч. II. Деформация. -М: Изд-во МИСИС, 1997.-527 с.

9. Multiscale Modeling and Simulation of Composite Materials and Structures. Eds Y.W. Kwon, D.H. Allen, R. Talreija - Springer Science+Business Media, LLC, 2008. - 630 p.

10. Overview of constitutive laws, kinematics, homogenization and multiscale methods in crystal plasticity finite-element modeling: Theory, experiments, applications / F. Roters, P. Eisenlohr, L. Hant-cherli, D.D. Tjahjanto, T.R. Bieler, D. Raabe // Acta Materialia. -2010. - Vol. 58. - P. 1152-1211. DOI: 10.1016/j.actamat.2009.10.058

11. Geers M.G.D., Kouznetsova V.G., Brekelmans W.A.M. Multi-scale computational homogenization: Trends and challenges // Journal of Computational and Applied Mathematics. - 2010. -Vol. 234. - P. 2175-2182. DOI: 10.1016/j.cam.2009.08.077

12. Multiscale modelling of plasticity and fracture by means of dislocation mechanics. Eds Pippan R., Gumbsch P. - Springer Wien New York, 2010. - 394 p.

13. Трусов П.В., Швейкин А.И.. Многоуровневые физические модели моно- и поликристаллов. Статистические модели // Физическая мезомеханика. - 2011. - Т. 14, № 4. - С. 17-28.

14. Трусов П.В., Швейкин А.И. Многоуровневые физические модели моно- и поликристаллов. Прямые модели // Физическая мезомеханика. - 2011. - Т. 14, № 5. - С. 5-30.

15. Raabe D. Computational materials science. - WILEY-VCH Verlag Gmbh, 1998. - 382 p.

16. Computational Methods for Microstructure-Property Relationships. Eds. Ghosh S., Dimiduk D. - Springer Science+ Business Media, LLC, 2011. - 658 p. DOI: 10.1007/978-1-4419-0643-4

17. Hill R. Elastic properties of reinforced solids - some theoretical principles // Journal of Mechanics and Physics of Solids. - 1963. - Vol. 11. - No. 5. - P. 357-372.

18. Moussaddy H., Therriault D., Levesque M. Assessment of existing and introduction of a new and robust efficient definition of the representative volume element // International Journal of Solids and Structures. - 2013. - Vol. 16. - P. 3817-3828. DOI: org/10.1016/j.ijsolstr.2013.07.016

19. Determination of the size of the representative volume element for random quasi-brittle composites / Pelissou, Baccou, Monerie, F. Perales // International Journal of Solids and Structures. -2009. - Vol. 46. - P. 2842-2855. DOI: 10.1016/j.ijsolstr.2009.03.015

20. Shavshukov V., Tashkinov A. Quantum Field Theory Approach to Mechanics of Polycrystals // Solid State Phenomena. - 2016. -Vol. 243. - P. 131-138. DOI: 10.4028/www.scientific.net/SSP.243.131

References

1. Antolovich S.D., Armstrong R.W.. Plastic strain localization in metals: origins and consequences, Progress in Materials Science, 2014, vol. 59, pp. 1-160. DOI: 10.1016/j.pmatsci.2013.06.001

2. Pineau A., Benzerga A.A., Pardoen T.. Failure of metals I: Brittle and ductile fracture, Acta Materialia, 2016, vol. 107, pp. 424-483. DOI: 10.1016/j.actamat.2015.12.034.

3. Pineau A., McDowell D.L., Busso E.P., Antolovich S.D. Failure of metals II: Fatigue, Acta Materialia, 2016, vol. 107, pp. 484-507. DOI: 10.1016/j.actamat.2015.05.050.

4. Pineau A., Benzerga A.A., Pardoen T.. Failure of metals Ш. Fracture and fatigue of nanostructured metallic materials, Acta Materialia, 2016, vol. 107, pp. 508-544. DOI: 10.1016/j.actamat.2015.07.049.

5. Voight W. Lerbuch der Kristallphysuk. Leipzig und Berlin, Teubner, 1928, 978 p.

6. Reuss A., Berechnunug der Fliebgrenze von Mischkristallen auf Grund der plastizitatsbending fur einkristalle. Z. Angew. Math. und Mech., 1929, vol. 9, no. 1, pp. 49-56.

7. Shermergor T.D. Teorija uprugosti mikroneodnorodnykh sred [Theory of elasticity of micro inhomogeneous media]. Moscow, Nauka, 1977, 400 p.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

21. Ташкинов А.А., Шавшуков В.Е. Решение задач механики деформирования поликристаллических материалов на основе теории возмущений // Вычислительная механика сплошных сред. -2016. - Т. 9, № 4. - С. 486-497. DOI: 10.7242/1999-6691/2016.9.4.41

22. Лифшиц И.М., Розенцвейг Л.Н. К теории упругих свойств поликристаллов // Журн. эксперимент. теор. физ. -1946. - Т. 16, № 11. - С. 967-980.

23. Eshelby J.D. The determination of the elastic field of an ellipsoidal inclusion and related problems // Proc. R. Soc. A 241. -1957. - P. 376-396. DOI: 10.1098/rspa.1957.0133

24. Кривоглаз М.А., Черевко А.С. Об упругих модулях твердой смеси // Физика металлов и металловедение. - 1959. -Т. 8, № 2. - С. 164-164.

25. Christensen R. Mechanics of Composite Materials. - New York: John Wiley & Sons, 1979. - 336 p.

26. Кунин И.А., Соснина Э.Г. Эллипсоидальная неоднородность в упругой среде // Докл. Акад. наук СССР. - 1971. -Т. 9, № 3. - С. 571-574.

27. Mura T. Micromechanics of Defects in Solids. -Dordrecht. - Martinus Nijhoff Publishers, 1987. - 587 p.

28. Шавшуков В.Е. Распределение полей напряжений в поликристаллических материалах // Физическая мезомехани-ка. - 2012. - Т. 15, № 6. - С. 85-91.

29. Шавшуков В.Е. Упругое взаимодействие зерен в поликристаллических материалах // Вестник Пермского национального исследовательского политехнического университета. Механика. -2014. - № 4. - С. 197-220. DOI: 10.15593/perm.mech/2014.4.08

30. Богачев И.Н., Вайнштейн А.А., Волков С.Д. Статистическое металловедение. - М.: Металлургия, 1984. - 176 с.

31. Францевич И.Н., Воронов Ф.Ф., Бакута С.А. Упругие постоянные и модули упругости металлов и неметаллов. -Киев: Наукова думка, 1982. - 286 с

32. Tadao Watanabe. An approach to grain boundary design for strong and ductile polycrystals // Res Mechanica. - 1984. -Vol. 11. - No. 1. - P. 47-84.

33. Randle V. Grain boundary engineering: an overview after 25 years // Materials Science and Technology. - 2010. - Vol. 26. -No. 3. - P. 253-261.

34. Watanabe T. Grain boundary engineering: historical perspective and future prospects // Journal of Materials Science. -2011. - Vol. 46. - P. 4095-4115. DOI: 10.1007/s10853-011-5393-z

8. Shtremel M.A. Prochnost splavov. Chast II. Defoimatsia [Strength of alloys. Part II. Deformation]. Moscow, MISIS, 1997, 527 p.

9. Kwon Y.W., Allen D.H., Talreija R. Multiscale Modeling and Simulation of Composite Materials and Structures. Springer Science+Business Media, LLC, 2008, 630 p.

10. Roters F., Eisenlohr P., Hantcherli L., Tjahjanto D.D., Bieler T.R., Raabe D. Overview of constitutive laws, kinematics, homogenization and multiscale methods in crystal plasticity finite-element modeling: Theory, experiments, applications, Acta Materialia, 2010, vol. 58, pp. 1152-1211. DOI: 10.1016/j.actamat.2009.10.058.

11. Geers M.G.D., Kouznetsova V.G., Brekelmans W.A.M. Multi-scale computational homogenization: Trends and challenges, Journal of Computational and Applied Mathematics, 2010, vol. 234, pp. 2175-2182. DOI: 10.1016/j.cam.2009.08.077.

12. Pippan R., Gumbsch P Multiscale modelling of plasticity and fracture by means of dislocation mechanics. SpringerWienNewYork, 2010, 394 p.

13. Trusov P.V., Shveikin A.I. Mnogoyrovnevye phisicheskie modeli mono- i polikristallov. Statisticheskie modeli. [Multilevel physical modeling of single- and polycrystals. Statistical models].

Physicheskaya mesomechanika - Physical mesomechanics, 2011, vol. 14, no. 4, pp. 17-28.

14. Trusov P.V., Shveikin A.I. Mnogoyrovnevye phisicheskie modeli mono- i polikristallov. Pryamyie modeli. [Multilevel physical modeling of single- and poly crystals. Direct models]. Physicheskaya mesomechanika - Physical mesomechanics, 2011, vol.14, no. 5, pp. 5-30.

15. Raabe D.. Computational materials science. WILEY-VCH Verlag Gmbh. 1998. 382 p.

16. . Ghosh S., Dimiduk D. Computational methods for microstructure-property relationshipseds. Springer Science+Business Media, LLC, 2011, 658 p. DOI: 10.1007/978-1-4419-0643-4.

17. Hill R. Elastic properties of reinforced solids - some theoretical principles. Journal of Mechanics and Physics of Solids, 1963, vol. 11, no. 5, pp. 357-372.

18. Moussaddy H., Therriault D., Levesque M.. Assessment of existing and introduction of a new and robust efficient definition of the representative volume element. International Journal of Solids and Structures, 2013, vol. 16. pp.3817-3828. DOI: 10.1016/j.ijsolstr.2013.07.016.

19. Pelissou C., Baccou J., Monerie Y., Perales F.. Determination of the size of the representative volume element for random quasi-brittle composites. International Journal of Solids and Structures, 2009, vol. 46, pp 2842-2855. DOI: 10.1016/j.ijsolstr.2009.03.015

20. Shavshukov V. and Tashkinov A.. Quantum Field Theory Approach to Mechanics of Polycrystals. Solid State Phenomena, 2016, vol. 243, pp 131-138. DOI: 10.4028/www.scientific.net/SSP.243.131

21. Tashkinov A.A., Shavshukov V.E. Reshenie zadach mekhaniki deformirovaniya polikristallicheskikh materialov na osnove teorii vozmushchenii [Solving the problems of deformation mechanics of polycrystalline materials on the basis of perturbation theory]. Vychislitelnaia mekhanika sploshnykh sred -Computational continuum mechanics, 2016, vol.16, no. 4, pp. 486497. DOI: 10.7242/1999-6691/2016.9.4.41

22. Lifshitz I.M., Rosenzweig L.N. K teorii uprugikh svoistv polikristallov [On the theory of elastic properties of polycrystals]. Jurnal eksper. i teoreticheskoi fiziki - Journal of Experimental and Theoretical Physics, 1946, vol.16, no. 11, pp.967-980.

23. Eshelby, J. D. The determination of the elastic field of an ellipsoidal inclusion and related problems. Proc. R. Soc. A 241, 1957, pp. 376-396. DOI: 10.1098/rspa.1957.0133.

24. Krivoglaz M.A., Cherevko A.S. Ob uprugikh modulyakh tverdoi smesi [On elastic moduli of solid mixtures]. Fizika metallov I metallovedenie - Physics of metals and metallography,1959, vol. 8, no.2, pp. 164-164.

25. Christensen R. Mechanics of Composite Materials. New York: John Wiley & Sons, 1979, 336 p.

26. Kunin I.A., Sosnina E.G. Ellipsoidalnaia neodnorodnost v uprugoi srede [ Ellipsoidal inhomogeneity in elastic medium]. Doklady Akademii nauk SSSR -Proceedings of USSR academy of sciences, 1971, vol. 9, no. 3, pp. 571-574.

27. T.Mura. Micromechanics of Defects in Solids. Dordrecht: Martinus Nijhoff Publishers, 1987, 587 p.

28. Shavshukov V.E. Raspredeleniye polei napryagzenii v polikristallicheskikh materialakh [Distribution of stress fields in poly crystalline materials]. Fizicheskaya mezomekhanika -Physical Mesomechanics, 2012, vol. 15, No.6, pp.85-91.

29. Shavshukov V.E. Elastic interaction of grains in polycrystalline materials. PNRPU Mechanics Bulletin. 2014, no. 4, pp.197-220. D0I:10.15593/perm.mech/2014.4.08.

30. Bogachev I.N., Vainshtein A.A., Volkov S.D. Statisticheskoje metallovedenije [Statistical metallography]. Moscow, Metallurgiya, 1984, 176 p.

31. Franchevich I.N., Voronov F.F., Bakuta S.A. Uprugie postoyannye I moduli uprugosti metallov I nemetallov [Elastic constants and moduli of merals and nonmetals]. Kiev, Naukova dumka, 1982, 286 p.

32. Tadao Watanabe. An approach to grain boundary design for strong and ductile polycrystals. Res Mechanica, 1984, vol.11, no.1, pp. 47-84.

33. V.Randle. Grain boundary engineering: an overview after 25 years. Materials Science and Technology, 2010, vol.26, no.3, pp 253-261.

34. Tadao Watanabe. Grain boundary engineering: historical perspective and future prospects. Journal of Materials Science, 2011, vol. 46, pp. 4095-4115. DOI: 10.1007/s10853-011-5393-z

i Надоели баннеры? Вы всегда можете отключить рекламу.