Статистические закономерности распределения напряжений в поликристаллах

Традиционно [1] моделирование поведения металла ведется для трех структурных уровней (рис.1).

а)                                         б)                                в)

Рис.1. Три структурных уровня металла:

а) макроуровень - инженерные модели (модель сплошной среды);

б) мезоуровень - инженерно-физические модели;

в) микроуровень - физические модели

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

Для каждого структурного уровня можно ввести свои параметры состояния: напряжения и деформации. При этом параметры макроуровня можно рассматривать как результат осреднения соответствующих параметров мезоуровня. Например, связь между макронапряжениями с и мезонапряжениями а может быть записана через операцию осреднения [2] по представительному объему V следующим образом:

С = ^ = ± j a ( X ) dV .                                  ⁽¹⁾

^V V V

Целью данной работы является изучение с помощью соответствующей математической модели статистических параметров распределения мезонапряжений а в различных поликристаллах и неоднородных материалах, а также эффективного напряжения ( a^V в зависимости от величины объема осреднения V.

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

Исследуемая среда занимает квадратную область Q (рис.2) с границей B = U B i , i = 1

( Q = Q U B) . Область Q разбивается на n квадратных ячеек Q i (i=1,2,...,n) (соответственно, сторона области разбивается на m равных отрезков: n=m² ). В пределах ячейки Q i свойства материала принимаются однородными. По сторонам области Q задаются однородные граничные условия. Напряженно- деформированное состояние в любой ячейке Q i удовлетворяет:

уравнениям равновесия

V- a = 0, x e Q i ,                                                       ⁽²⁾

определяющим уравнениям а = C “Р, x e Q.,

геометрическим соотношениям в = 1 (V u + V u T), x eQi,

Рис.2. Схема нагружения

где а - тензор мезонапряжений, в - тензор мезоде-

^*

формаций, C - тензор модулей упругости материала, и - вектор перемещений.

Граничные условия на В:

• ^а)    x ^e B i^: U x = ^u , U y = °;

• б)    x e B 2 U B 4- a xy = °, U y = °;

• ^в)    x ^e B 3- a xy = °, U x = °.

Численная реализация алгоритма проводилась с помощью метода конечных элементов. Все расчеты проводились для области единичного размера при m=30, т.е. объем получаемых выборок для оценки статистических параметров составлял 900 элементов. Для сопоставимости статистических параметров для различных материалов использовались относительные мезонапряжения:

∗ α

α

1 а (x) d 1.

Ω

Были исследованы однофазные металлы с кубическим типом решетки. Для

характеристики степени анизотропии кристалла использовался параметр анизотропии [3]

A =

2C 44

C 11 - C 12

который для изотропного материала равен единице. В табл.1 приведены данные о модулях упругости и параметре анизотропии исследованных материалов. Принимался равномерный закон распределения кристаллографических осей материала по объемам Ω i .

Обработка полученных в результате решения краевых задач выборок а * показала,

что гипотеза о нормальном законе распределения компонент α ^∗ _ij может быть принята с

уровнем значимости 90%. Математическое ожидание M =

1 jet *( x) d 1

Ω

принимало

единичное значение.

Таблица 1.

Название материала	Модули упругости,кбар			А	Среднеквадратическое отклонение
	С11	С12	С44		S( α∗ xx )	S( α ∗ yy )	S( α ∗ xy )
Натрий	74	62	42	7.00	16.46	26.39	3.07
Медь	1684	1214	754	3.21	3.66	6.09	26.04
Тантал	2650	1590	831	1.57	0.56	1.43	28.58
Алюминий	1073	609	283	1.22	0.12	0.50	9.37

В качестве характеристики разброса   α∗ij   определялась величина среднеквадратического отклонения

S

а *— м ) d ¹ ,

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

Исследование зависимости степени разброса эффективного напряжения ^ а ") от величины объема осреднения V выполнялось следующим образом. В качестве области осреднения выбирался квадрат размером m на m ячеек. Последовательным смещением выбранного квадрата по области Ω на одну ячейку по горизонтали или вертикали получали выборку из K объемов осреднения. В данном случае принимали K=100. Затем находили среднее для каждого V_i эффективное напряжение

где η = m , i = 1.. K .

Обрабатывая выборку из /^ \ , определяли математическое ожидание и дисперсию:

KK     ²   (9)

*                            *                    *

M ( \а /г^Ё \а )D^D ( \а /)^-^ Ё \а - MM ( \а /) • K j = 1       'J                 K j - 1         jj

С точки зрения математической статистики задачу определения дисперсии

D(\О6 /) можно сформулировать следующим образом. Имеется выборка n (не обязательно независимых) случайных величин аа 2, ••■, *„ с математическим ожиданием ц и дисперсией S2. Из данной выборки случайным образом выделяется K подвыборок объемом по п элементов и находится среднее [а 3 по каждой подвыборке. Требуется определить математическое ожидание и дисперсию найденного среднего а / • Из свойств математического ожидания и дисперсии можно записать:

M (ка > M (1 Еа *1 = 1Ё M (а *) - *• х ' k V J=1 j ) Л J=1      ' JJ

(^п А ( Л

^{D (} к«>^D ( п ^ к)- _? ^D | ^- j j .

Для дисперсии суммы случайных величин имеем[6]:

D

п *

^Е « j k j - 1 J

П П

- Е Е P ijSiSj , i - 1 j - 1

где p jj - коэффициент корреляции между элементами а * и а * , выражаемый через ковариацию X jj:

j                                                                        *               *          *                *

P j ^- P ji ^-7Т=Г’ P i ^- 1, где L j ^- L ji = M [ (б а , - M ⁽6 a J)(6 a. j - M (^ a j) ] ;

V ^L«^L jj

*        ~*

Si и Sj - среднеквадратические отклонения элементов a, и a j, равные в данном случае S. Учитывая симметрию матрицы моментов Pij, дисперсию [^^ / можно окончательно представить в следующем виде:

2 ( п -1 п А

D(о^■)-ST п+2ЕЕ Pij п k ,-1 j -,+1     )

Если величины a * ,( a 2 , ^ ,( Z*„ независимы, то P ij =0 при i ^ j В этом случае d ((< a’> SI-

Величина p ij по модулю не превышает 1, поэтому для дисперсии эффективных напряжений получаем следующее ограничение:

0 ^ D ( Oz > S 2 .                                              ⁽¹³⁾

Таким образом, из ограничения (13) следует, что в результате осреднения α∗ij по объему V величина дисперсии а *) не может превосходить дисперсию мезонапряжений. Если ввести понятие относительной дисперсии эффективных напряжений в виде

D ( а *))

S ²

D ⁼

то можно сравнить для различных материалов зависимость D* от величины объема осреднения мезонапряжений. При этом, по близости D* (m) к функции 1/η можно будет судить о независимости мезонапряжений.

Рис.3. Зависимость D*(m)

На рис.3 для рассмотренных выше металлов представлены зависимости D^*(m). Как можно видеть из приведенных результатов, величина относительной дисперсии при увеличении объема осреднения стремится к зависимости 1/ η . Данное обстоятельство позволяет сделать следующие выводы:

• 1.    С увеличением области осреднения дисперсия ^Об *) уменьшается, т.е. уменьшается случайность эффективных напряжений. Если задаться малой положительной величиной 5 и условиться, что при снижении дисперсии ^ *^ ниже выбранной 5 последние можно считать дисперсными (т.е. не случайными), то для любого материала можно подобрать такой объем осреднения V, при котором будет выполняться ограничение

• 2.    Так как D^*(m) при увеличении m стремится к 1/ η , то можно сделать вывод, что зависимость напряженного состояния в данной точке от напряжений в соседних убывает с увеличением расстояния между точками, т.е. чем точки дальше друг от друга, тем их состояния менее коррелированы друг с другом. Аналогично δ можно ввести малую положительную величину ε , являющуюся мерой независимости состояний в двух точках тела. Если принять, что ε равно δ , то условие (15) может служить критерием независимости напряженных состояний.

В работе [2] экспериментально исследована зависимость мезонапряжений в соседних точках. Показано, что в однофазных поликристаллах закон распределения мезонапряжений близок к нормальному. Поэтому для выяснения вопроса о независимости напряжений строились функции корреляции. Было установлено, что на расстоянии 4-5 средних диаметров зерна для алюминия и 6-7 средних диаметров для меди напряженно-деформированное состояние практически некоррелировано. Аналогичные выводы были сделаны в работах [4,5]. Расстояние, на котором напряженно-деформированное состояние в зернах можно считать независимым, предлагалось выбирать в качестве размера представительного объема на макроуровне. В этом случае, если взять представительный объем в виде куба, то для алюминия в нем должно содержаться более 50 зерен, а для меди - более 200. Для рассмотренного в данной работе плоского конечно-элементного автомата это соответствует разбиению 7 на 7 для алюминия и 14 на 14 для меди. Как можно видеть из рис.3, величины относительных дисперсий для данных материалов при соответствующих m близки друг к другу и равны примерно 0.1. Данное обстоятельство позволяет сделать вывод, что 5 в критерии (15) можно выбирать равным 0.1.

Вторая серия численных экспериментов с конечноэлементным автоматом была посвящена исследованию статистических свойств представительного объема для изотропных неоднородных материалов, состоящих из двух фаз. В табл.2 приведены cвойства фаз, отличающиеся значениями модуля С11 (для изотропного материала С44=(C11-C12)/2). В качестве характеристики материала использовался коэффициент неоднородности:

K =

Таблица 2

Номер фазы	Модули упругости,кбар			К	S( а ) при 50% 2-й фазы
	С11	С12	С44		S( a xx )	S( a yy)	S( a xy)
1	100	60	20	-	-	-	-
2	300	60	120	3	0.17	0.38	127.16
2	600	60	270	6	0.32	0.78	93.68
2	1000	60	470	10	0.39	1.23	23.14

Получаемый закон распределения мезонапряжений существенно зависит как от процентного содержания второй фазы, так и от граничных условий на стороне В 1 . При малом содержании одной из фаз и задании на В 1 граничного условия в перемещениях деформация осуществляется в основном за счет той фазы, объемная доля которой больше. В этом случае наибольший разброс напряжений наблюдается при примерно равном содержании фаз. Если на В 1 задать усилия, то наибольший разброс будет наблюдаться, если доля одной из фаз составляет примерно 10%. Кроме того, в этом случае получаемая выборка напряжений является не унимодальной, т.е. содержит 2 или даже 3 горба. Последнее обстоятельство можно объяснить резким разделением значений деформации жестких включений, прилегающих к ним областей и мягкой матрицы.

Влияние коэффициента неоднородности и процентного содержания второй фазы на величину среднеквадратического отклонения S(axx) показано на рис.4,а. Наибольший разброс наблюдается, если доля второй фазы составляет 50%. В табл.2 для 50%-ного содержания фаз приведены значения S( а), которые показывают, что с ростом коэффициента неоднородности материала величина среднеквадратического отклонения увеличивается.

а)

б)

Рис.4. Влияние коэффициента неоднородности и процентного содержания на S( α xx )

Зависимость относительной дисперсии D^*( α xx ) от величины объема осреднения m и от коэффициента неоднородности показана на рис.4,б. Приведенные результаты позволяют сделать вывод, что для подобной модели неоднородного материала зависимость относительной дисперсии от объема осреднения убывает практически обратно пропорционально числу элементов осреднения, т.е. напряженно-деформированное состояние в соседних ячейках слабо зависит друг от друга.

Итак, проведенные численные эксперименты и сравнение с экспериментальными результатами показали, что полученное соотношение (15) можно использовать в качестве критерия для определения размеров представительного объема. Данный критерий имеет следующие преимущества перед критерием [2], использующим функцию корреляции.

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

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

Авторы выражают свою благодарность к.т.н. М.Б.Гитману и проф. П.В.Трусову за участие в обсуждении данной работы и высказанные ими ценные замечания.