Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Савёлова Методы решения некорректных задач 2012

.pdf
Скачиваний:
16
Добавлен:
12.11.2022
Размер:
6.26 Mб
Скачать

значим c

ni

 

– вес i -го зерна. В качестве ориентации i -го зерна

 

i

N

 

 

 

 

 

 

 

 

 

берем середину i -го интервала.

 

 

 

На рис. 8.13,ав, 8.14,ав изображены графики ЦНР

с

ε 1 2, ε 1 8

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

для

n 50 (случай

а),

n 100 (случай б), обычная выборка (случай в) для N 500 . Точ-

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

Из приведенных рисунков видим, что для ε 12 различия не видно, для ε 18 для случая 8.14,а график сильно изрезан в об-

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

чае острых текстур ( ε мало) требуется увеличивать объем выборки.

а

б

в

Рис. 8.13. Моделирование зеренной структуры ε 12

201

а

б

в

Рис. 8.14. Моделирование зеренной структуры ε 18

Контрольные вопросы

1. Ядерные методы вычисления ФРО и ПФ по набору отдель-

ных ориентаций на SO(3)

2. Проекционные методы вычисления ФРО и ПФ по набору от-

дельных ориентаций.

3.Вычисление математического ожидания и дисперсии для оценки ФРО и ПФ ядерными и проекционными методами.

4.Оценки точности вычисления ФРО и ПФ по набору отдель-

ных ориентаций.

5.Свойство состоятельности ядерных и проекционных оценок

ФРО.

6.Как влияет зависимость ориентаций на вычисление ФРО?

7.Выбор параметра сглаживания в ядерных методахвычисле-

ния ФРО.

8. Как влияет симметрия кристалла на точность вычисления

ФРО и ПФ?

202

9. Устойчивы ли ядерные и проекционные методы относитель-

но погрешностей измерения ориентаций?

203

Глава 9. Применение ФРО при исследовании физических

свойств

Многие физические свойства поликристаллических материалов анизотропны [2, 10, 13, 24, 25, 40, 47, 75, 108]. Свойство поликри-

сталла, измеренное в некотором направлении, – среднее значение свойств индивидуальных зерен, составляющих поликристалл. Та-

ким образом, данное свойство зависит от свойства монокристалла и распределения ориентаций кристаллитов в поликристалле. Если это распределение равномерное, то материал макроскопически изо-

тропен. Наличие преимущественных ориентаций (текстуры) при-

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

кристаллического материала: наличие пор или трещин, форма кри-

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

венной причиной является текстура. В таких случаях свойства по-

ликристалла могут быть найдены путем усреднения соответст-

вующих свойств образующих поликристалл зерен с функцией рас-

пределения зерен по ориентациям.

9.1. Вычисление средних значений свойств поликристалла

 

 

 

 

 

 

 

 

Среднее значение

E свойства поликристалла полностью опре-

деляется свойствами составляющих его зерен и ФРО:

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

E

 

E r dV

 

 

 

 

E g

dV E g f g dg, (9.1)

V

V

 

 

 

 

V

SO(3)

V g

SO(3)

 

 

 

 

 

 

 

 

 

204

 

E g
E g

где r – орт направления; – свойство кристаллита, имеюще-

го ориентацию g ; в системе координат образца (в зависимости от

того, каким образом математически описывается свойство кристал-

лита – с помощью тензора или указательной поверхности, величи-

ны E и в (9.1) являются тензорами или функциями направ-

лений соответственно).

Формулу (9.1) можно переписать в более наглядном виде. Пусть

EB – свойство отдельного монокристалла в системе координат кристаллита KB . Допустим, имеется образец с закрепленной с ним системой координат K A . Образец состоит из N кристаллитов с

объемами Vi , у каждого из которых – своя система координат KB ,

 

 

 

 

 

 

i

связанная

с

системой K A

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

gi

i 1, 2,..., N . При этом свойство

EA представляет собой супер-

позицию свойств отдельных монокристаллов:

 

 

 

 

N

 

 

 

 

 

 

EA ai

EB

,

(9.2)

 

 

 

i 1

 

gi

 

 

 

 

 

 

 

где ai Vi

V ,

EB

– свойство EB в системе K A .

 

 

 

 

gi

 

 

 

При переходе к непрерывному случаю сумма заменяется инте-

гралом на группе вращений трехмерного евклидова пространства

SO 3 :

EA EB g f g dg ,

(9.3)

SO(3)

205

f g

где в качестве весовой функции в интеграле берется нор-

мальное распределение канонического вида (2.3).

По формуле (9.3) проводится вычисление усредненных значе-

ний модулей упругости и модулей упругой податливости. Изло-

женный материал представляет также достаточно полное описание методов и результатов численного расчета упругих усредненных свойств поликристаллического магния (+4,5%Al+1%Zn) и титана

(экспериментально чистый образец).

9.2. Вычисление усредненных упругих свойств для

поликристаллов магния и титана

Для измеренных на спектрометре количественного анализа тек-

стуры (СКАТ) материалов магния (рис. 9.1,а) и титана (рис.

9.2,а) с экспериментальными полюсными фигурами, были числен-

но построены модельные ПФ для магния (рис. 9.1,б) и титана (рис.

9.2, б), которые с выбранными параметрами полуширин распреде-

ления могут аппроксимировать экспериментальные ПФ.

Рис. 9.1. Экспериментальная (а) и расчетная (б) ПФ для магния с параметрами полуширины α11 0,9; α22 0,8; α33 0, 2

206

Рис. 9.2. Экспериментальная (а) и расчетная (б) ПФ для титана с параметрами полуширины α11 0, 01, α22 0, 4, α33 0,8

Ниже проиллюстрирован метод расчета упругих характеристик,

в данном случае – модулей упругой податливости, так как для уп-

ругих свойств он наиболее интересен.

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

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

Метод расчета изложен на примере образцов, обработанных го-

рячей и холодной прокаткой. Для таких обработок (присутствие оси прокатки) находят как данный алгоритм, так и методы опреде-

ления упругих характеристик с помощью указательной поверхно-

сти.

Усреднение производится с помощью интегрирования произве-

дения функции распределения и тензора четвертого ранга:

207

 

 

 

 

 

ijkl A

 

1

2

 

ijkl

 

B

 

 

 

 

 

 

 

E

 

 

E

g

 

f

 

g dg.

(9.4)

 

 

 

 

SO(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Под EB

g

Eijkl g

подразумеваются

компоненты

тензора

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

рицы поворота по закону преобразования тензора четвертого ранга:

Eijkl g Aim Ajn Ako Alp Eijkl .

(9.5)

Если в качестве входных параметров используются модули же-

сткости или модули упругой податливости, то можно записать

Eijkl cijkl

или

Eijkl sijkl .

В работе усреднение производится с помощью интегрирования методом Гаусса.

Интеграл (9.3) приводится к виду

 

 

ijkl

1

 

 

 

ijkl g B

N aii , 0 (g)dg,

 

E

 

E

(9.6)

2

 

 

 

 

SO(3)

 

 

где Eijkl g

 

– свойство отдельного

кристаллита. В

качестве

 

B

 

 

 

 

 

 

 

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

Упругие свойства гексагональных кристаллов магния и титана описываются пятью независимыми параметрами. Такими парамет-

рами являются, например, модули упругости: C11, C12 , C13 , C33 ,

208

C44 или модули упругой податливости S jk . Модули C jk , ГПа, и

S

jk

10 2 , ГПа-1, монокристаллов магния и титана представлены в

 

 

табл. 9.1.

 

На рис. 9.1,а и 9.2,а изображены экспериментальные ПФ магния

и титана соответственно, на рис. 9.1,б и 9.2,б – расчетные ПФ, по-

лученные аппроксимацией ПФ от КНР. Для магния с найденными

параметрами распределения α11 0,9,

α22 0,8, α33 0, 2 в цен-

тре g0 α0 0 , γ0 0, 0, 0 и титана

с параметрами α11 0, 01,

α22 0, 4, α33 0,8 . В табл. 9.2

приведены результаты расчетов

эффективного тензора S

jk

10 2 ,

ГПа-1, коэффициентов упругой

 

 

 

 

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

приближении Фойгта, Ройсса и Хилла: S

 

, S

P

,

S X

[24, 68, 94,

 

 

 

 

 

 

 

 

 

 

ij

 

 

ij

 

 

ij

 

 

 

109].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 9.1. Модули ( C

jk

, ГПа) и (

S

jk

10 2 , ГПа-1) для

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

монокристаллов магния и титана

Магний, C jk , ГПа

 

 

 

Магний, S jk 10 2 , ГПа-1

 

 

59,3

25,7

21,4

 

0

0

0

2,20

0,78

 

0,50

 

0

0

0

 

25,7

59,3

21,4

 

0

0

0

0,78

 

2,20

 

0,50

 

0

0

0

 

21,4

21,4

61,5

 

0

0

0

0,50

 

0,50

 

1,97

 

0

0

0

 

0

0

0

 

16,4

0

0

0

 

 

0

 

0

 

 

 

6,09

0

0

 

0

0

0

 

0

16,4

0

0

 

 

0

 

0

 

 

 

0

6,09

0

 

0

0

0

 

0

0

16,8

0

 

 

0

 

0

 

 

 

0

0

5,96

 

Титан, C jk , ГПа

 

 

 

Титан, S jk

10 2 ,

 

ГПа-1

 

 

 

160,0

90,0

66,0

 

0

0

0

0,97

0,47

 

0,18

 

0

0

0

 

90,0

160,0

66,0

 

0

0

0

0,47

 

0,97

 

0,18

 

0

0

0

 

66,0

66,0

181,0

 

0

0

0

0,18

 

0,18

 

0,67

 

0

0

0

 

 

 

 

 

 

 

209

 

 

 

 

 

 

 

 

 

 

 

 

Продолжение табл. 9.1

Титан, C jk , ГПа

 

 

 

Титан,

S jk 10 2 ,

ГПа-1

 

 

0

0

0

46,5

0

0

0

0

 

0

2,15

 

0

0

0

0

0

0

46,5

0

0

0

 

0

0

 

2,15

0

0

0

0

0

0

35,0

0

0

 

0

0

 

0

2,88

Таблица 9.2. Значения вычисленных коэффициентов упругой по-

 

 

 

 

 

 

 

 

 

датливости

S

для поликристаллических магния

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Mg+4,5%Al+1%Zn) и титана (Ti)

Приближение

 

 

 

S

 

 

S

S

 

S

S

 

S

 

S

 

S

 

 

 

 

 

11

 

 

12

13

22

33

 

 

44

 

 

55

 

66

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ройсс (Mg)

 

 

 

2,25

 

 

–0,66

–0,65

2,30

2,21

 

5,79

 

 

5,79

 

5,82

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Фойгт (Mg)

 

 

 

2,25

 

 

–0,66

–0,64

2,25

2,21

 

5,78

 

 

5,78

 

5,82

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Хилл (Mg)

 

 

 

2,25

 

 

–0,66

–0,65

2,25

2,24

 

5,79

 

 

5,78

 

5,82

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sij (Mg)

 

 

 

0,01

 

 

0,01

0,01

0,01

0,01

 

0,01

 

 

0,01

 

0,01

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δSij , % (Mg)

 

 

0,07

 

 

0,89

1,18

0,10

0,15

 

0,19

 

 

0,19

 

0,11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ройсс (Ti)

 

 

 

0,88

 

 

–0,27

–0,26

0,90

0,88

 

2,36

 

 

2,29

 

2,36

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Фойгт (Ti)

 

 

 

0,88

 

 

–0,30

–0,27

0,91

0,88

 

2,38

 

 

2,31

 

2,38

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Хилл (Ti)

 

 

 

0,88

 

 

–0,29

–0,26

0,91

0,88

 

2,37

 

 

2,30

 

2,37

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Sij ( Ti)

 

 

 

0,01

 

 

0,02

0,02

0,01

0,01

 

0,02

 

 

0,02

 

0,03

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δSij , % (Ti)

 

 

 

0,87

 

 

7,76

6,15

0,70

0,95

 

1,04

 

 

0,84

 

1,05

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

также

 

приведены

значения

S

 

S P S

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ij

 

ij

ij

 

δSij

 

 

 

Sij

 

 

 

 

 

100 %

соответственно,

абсолютные и

 

 

 

 

 

 

 

 

 

 

 

 

min

 

Sij

 

,

 

Sij P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

относительные погрешности вычисления эффективных коэффици-

ентов упругой податливости.

210

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]