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

книги / Управление большими системами. УБС-2017

.pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
17.48 Mб
Скачать

Информационные технологии в управлении техническими системами и технологическими процессами

Вкачестве основных параметров зёренно-фазовой структуры исследуемого материала выбраны: объёмная доля фаз, соответственно ФКС, М и Б, а также параметры зёренной структуры – средний размер зерна, коэффициент вариации размера зерна, степень анизотропии зерен и ОДЗ.

Вкачестве эксплуатационных характеристик исследуемого материала выбраны, например, такие прочностные характеристи-

ки, как твердость по Виккерсу (HV30) и ударная вязкость (КС). Отметим, что, как указано в работе [2], в отличие от обыч-

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

значать как μi для микросвойств и ν i – для макросвойств ис-

следуемого материала) соответствует функции принадлежности элемента в обычном нечетком множестве.

Отметим, что общее количество пар элементов в нечетком множестве Ai равно значению d, а общее количество пар эле-

ментов в нечетком множестве Вi равно значению w, которые определяются следующим образом:

 

k

r

(2)

d = pl ;

w = tl .

 

i=1

i=1

Таким образом, для каждого i-го эксперимента (i = 1, …, l) общий вид специальных нечетких множеств Ai и Вi, которых можно представить следующим образом:

Ai = (

ФКС

,

М

, ... ,

ОДЗ

) ,

p1 пар элементов

 

p2 пар элементов

 

pk пар элементов

Ci = (

c1

,

c2

, ... ,

cr

) .

 

t1 пар элементов

 

t2 пар элементов

 

tr пар элементов

 

Задача оценки близости исследуемого материала к определенному классу, характеризующемуся заданным комплексом физико-механических характеристик C*, теперь определится мерой близости параметров зёренно-фазовой структуры исследуемого материала к параметрам зёренно-фазовой структуры

463

501

Управление большими системами. Выпуск XX

всех тестовых образцов (число которых равно l), для каждого из которых известно, принадлежит ли этот образец к классу с комплексом физико-механических характеристик C* или нет.

Для определенности будем считать, что образцы с номерами от 1 до s, (s l) принадлежат этому классу, а образцы с номерами от s + 1 до l – нет. Очевидно, что количество образцов, которые не принадлежат к классу с комплексом физико-механических характеристик C*, равно m, где m = l – s. Отметим, что если все образцы принадлежат (или не принадлежат) классу с характеристиками C*, то задача принадлежности исследуемого образца этому классу сразу решается однозначно.

Построим процедуру оценки близости исследуемого материала к классу с комплексом физико-механических характеристик C*.

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

Теперь мера близости между специальными нечеткими множествами A и В (расстояние между этими множествами),

которую обозначим ρi (так как A соответствует i-му эксперименту), может быть определена следующим образом:

 

d

 

(аi

bi )

 

 

 

 

 

 

(3)

ρi = μi

 

 

 

 

,

max

(аi , bi )

 

i=1

 

 

 

где d, как указано выше, определяется количеством пар элементов μi / ai и μi / bi в нечетких множествах A и В соответственно.

Пример. В качестве демонстрационного примера формирования функционального материала рассмотрим получение арматурной стали с заданными прочностными свойствами путем термообработки с различными скоростями охлаждения (исследования проведены в НИИ наносталей при Магнитогорском государственном техническом университете [3]. Для выявления качественных и количественных характеристик формирующейся структуры использовался исследовательский комплекс GLEEBLE

464

502

Информационные технологии в управлении техническими системами и технологическими процессами

3500 с оптическим микроскопом Meiji Techno с применением системы компьютерного анализа изображений Thixomet PRO, а также сканирующий электронный микроскоп JSM 6490 LV).

Результаты обработки экспериментальных данных анализа зё- ренно-фазовой структуры и эксплуатационных характеристик исследуемой стали, в том числе с использованием интеллектуальных технологий [7, 11], для 7 образцов, полученных при различных режимахобработки(скоростяхохлаждения), приведенывтабл. 1.

Таблица1. Результатыобработкиданныханализазёренно-фазовой структурыиэксплуатационныххарактеристикстали

Скорость

Параметры

Параметры

Твердость

п/п

охлажде-

зёренной структуры

фазового

материала

 

ния, °C/с

(размерызернавмкм)

состава, %

HV30

1

10

Среднийразмерзерна: 0,10

ФКС: 100

400

 

 

Коэффициентвариациираз-

 

 

 

 

меразерна: 0,11

 

 

 

 

Степеньанизотропии: 0,98

 

 

 

 

Объемное содержание зе-

 

 

 

 

рен: 0,87

 

 

2

20

Среднийразмерзерна: 0,10

ФКС: 60

412

 

 

Коэффициентвариациираз-

Бейнит: 10

 

 

 

меразерна: 0,07

Мартенсит: 30

 

 

 

Степеньанизотропии: 1,14

 

 

 

 

Объемное содержание зе-

 

 

 

 

рен: 0,26

 

 

3

25

Среднийразмерзерна: 0,16

ФКС: 70

436

 

 

Коэффициентвариациираз-

Бейнит: 20

 

 

 

меразерна: 0,06

Мартенсит: 10

 

 

 

Степеньанизотропии: 1,09

 

 

 

 

Объемное содержание зе-

 

 

 

 

рен: 0,44

 

 

4

30

Среднийразмерзерна: 0,26

ФКС: 65

469

 

 

Коэффициентвариациираз-

Бейнит: 20

 

 

 

меразерна: 0,17

Мартенсит: 15

 

 

 

Степеньанизотропии: 1,18

 

 

 

 

Объемное содержание зе-

 

 

 

 

рен: 0,98

 

 

 

 

 

 

465

 

 

 

 

503

Управление большими системами. Выпуск XX

Окончание табл. 1

Скорость

Параметры

Параметры

Твердость

п/п

охлажде-

зёренной структуры

фазового

материала

 

ния, °C/с

(размерызернавмкм)

состава, %

HV30

5

40

Среднийразмерзерна: 0,25

ФКС: 50

554

 

 

Коэффициентвариациираз-

Бейнит: 30

 

 

 

меразерна: 0,03

Мартенсит: 20

 

 

 

Степеньанизотропии: 1,25

 

 

 

 

Объемное содержание зе-

 

 

 

 

рен: 0,37

 

 

6

50

Среднийразмерзерна: 0,21

ФКС: 40

617

 

 

Коэффициентвариациираз-

Бейнит: 30

 

 

 

меразерна: 0,06

Мартенсит: 30

 

 

 

Степеньанизотропии: 1,1

 

 

 

 

Объемное содержание зе-

 

 

 

 

рен: 0,79

 

 

7

60

Среднийразмерзерна: 0,13

ФКС: 10

800

 

 

Коэффициентвариациираз-

Мартенсит: 90

 

 

 

меразерна: 0,07

 

 

 

 

Степеньанизотропии: 1,08

 

 

 

 

Объемное содержание зе-

 

 

 

 

рен: 0,91

 

 

Значения важности, использованные при построении функции расстояния, указаны в табл. 2.

Таблица 2. Значения важности при построении функции расстояния

 

 

 

 

 

Коэффи-

 

 

Название

 

 

Мар-

Средний

циент

Степень

Объемное

параметра

ФКС

Бейнит

тенсит

размер

вариации

анизо-

содержание

 

 

 

 

зерна

размера

тропии

зерен

 

 

 

 

 

зерна

 

 

Важность

0,3

0,4

0,5

0,7

0,6

0,6

0,8

Для указанных значений важности по формуле (3) получим матрицу расстояний между элементами:

466

504

Информационные технологии в управлении техническими системами и технологическими процессами

 

 

0

1,883

1,981

1,839

2,495

1,858

1,24

 

 

1,883

0

1,278

1,865

1,536

1,376

1,747

 

 

 

 

(4)

1,981

1,278

0

1,332

1,225

1,121

1,737

 

1,839

1,865

1,332

0

1,38

1,217

1,881

 

 

 

2,495

1,536

1,225

1,38

0

1,136

2,264

 

 

 

1,858

1,376

1,121

1,217

1,136

0

1,427

 

 

 

 

 

 

1,24

1,747

1,737

1,881

2,264

1,427

0

 

 

 

 

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

(5)

 

 

ρ 0 (x, y) =

 

 

x y

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

225

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисляя расстояния между элементами, получим

 

 

 

 

0

0,0533

0,16

 

0,306

0,684

0,965

1,77

 

 

 

0,0533

0

0,106

 

0,253

0,631

0,911

1,72

 

 

 

 

 

(6)

 

0,16

0,106

0

 

0,146

0,524

0,805

1,61

 

0,306

0,253

0,146

 

0

 

 

 

0,378

0,658

1,47

 

 

 

0,684

0,631

0,524

 

0,378

0

0,28

1,09

 

 

 

0,965

0,911

0,805

 

0,658

0,28

0

0,813

 

 

 

 

 

 

 

1,77

1,72

1,61

 

1,47

1,09

0,813

0

 

 

 

 

 

Сравнивая полученные матрицы, можно заметить, что полученные значения плохо согласуются.

Для сравнения вычислим расстояние Минковского по формуле [8]:

 

 

(ρ 0

N 1

N

ρij0 ρij1

p

1 p

(7)

lp

, ρ1 ) =

.

 

 

 

i=1

j=i+1

 

 

 

Получаем значение l1 = 19,54. Несоответствие мер близости

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

467

505

Управление большими системами. Выпуск XX

и возможной недостаточной точностью экспертного выбора значений важности.

Для того чтобы улучшить согласованность функций расстояния в пространствах значений твердости образцов и параметров микроструктуры, найдем значения важностей, как решение задачи минимизации расстояния Минковского l1. Данная задача является задачей оптимизации функции нескольких переменных с ограничениями типа неравенств. Наличие ограничений обусловлено требованиями, предъявляемыми к значениям функций принадлежности для специальных нечетких множеств. Решение задачи с ограничениями сводится к решению задачи без ограничений методом штрафных функций [3], которая,

всвою очередь, решается методом Нелдера–Мида [1].

Втабл. 3 представлены значения важности, полученные при таком подходе.

Таблица 3. Значения важности

 

 

 

 

 

 

Средний

Коэффи-

Степень

 

Объемное

 

 

 

 

 

 

циент

 

Название

 

 

 

Мар-

 

содер-

параметра

ФКС

Бейнит

тенсит

 

размер

вариации

анизо-

 

 

жание

 

 

 

 

 

 

 

зерна

размера

тропии

 

 

зерен

 

 

 

 

 

 

 

 

зерна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Важность

1,0

0,15

0

 

0,16

0,11

0

 

 

0

 

Для полученных значений важности расстояния между

элементами будут равны

 

 

 

 

 

 

 

 

 

 

0

0,594

0,564

0,641

0,83

0,888

0,977

 

 

0,594

0

0,295

0,316

0,427

0,535

 

1,02

 

 

 

 

 

(8)

0,564

0,295

0

 

0,203

0,449

0,518

 

1,05

 

0,641

0,316

0,203

0

0,378

0,537

 

1,14

 

 

 

 

0,83

0,427

0,449

0,378

0

0,28

 

1,09

 

 

 

0,888

0,535

0,518

0,537

0,28

0

0,981

 

 

 

 

 

 

0,977

1,02

1,05

 

1,14

1,09

0,981

 

0

 

 

 

 

 

 

Сравнивая с (6), можно заметить, что значения функции расстояния с оптимизированными значениями важности согла-

468

506

Информационные технологии в управлении техническими системами и технологическими процессами

суются лучше, чем в (7). Расстояние Минковского для метрик в данном случае l1 = 5, 43. Однако нет возможности принять,

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

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

Чтобы решить такое несоответствие, необходимо изменить вид функций принадлежности специальных нечетких множеств

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

феррито-карбидной

смеси a1 .

Вид зависимости

представлен

в табл. 4.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 4. Вид зависимости

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Мар-

 

 

 

 

 

Коэффи-

 

 

 

 

 

 

 

 

 

 

 

 

Бей

 

Средний

 

циентва-

Степень

Объемное

ФКС

 

нит

тен-

 

размерзерна

риации

 

анизотро-

содержа-

 

 

 

 

 

сит

 

 

 

 

 

размера

 

 

пии

 

ниезерен

 

 

 

 

 

 

 

 

 

 

 

 

 

зерна

 

 

 

 

 

 

 

 

 

0,9×

 

 

 

 

0, 2×

 

0, 2(1+ q5 ×

0,1(1+ q6 ×

0,1(1+ q7 ×

 

 

 

 

 

 

 

a1

 

 

 

 

a1

 

 

a1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

×(1q1

 

a1

)

0,2

0,1

 

×(1q4

 

a1

)

×

 

 

1 )

×

 

 

1 )

×

 

1 )

100

 

 

 

100

 

100

 

 

 

100

 

 

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для определения расстояния при таком виде значений функций принадлежности запишем функцию расстояния в виде

(9) ρ (A, B) = μi (а1 ) + μi (b1 )

Exp(2

μi (а1 ) μi (b1 ) )

аi bi

.

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

2

 

 

 

 

max(

 

аi

 

,

 

bi

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

469

 

 

 

507

Управление большими системами. Выпуск XX

Отметим, что в случае равенства коэффициентов q1 , q4 , q5 , q6 ,q7 нулю формула (9) аналогична (3). Значения коэффициен-

тов вычислим также из решения задачи оптимизации. Вычисленные коэффициенты представлены в табл. 5.

Таблица 5. Вычисленные коэффициенты

q1

q4

q5

q6

q7

0,8

1,0

1,0

1,0

0,37

Полученная в результате матрица расстояний примет вид:

 

 

0

0,679

0,59

0,595

0,913

1,0

1,85

 

 

0,679

0

0,323

0,375

0,415

0,503

1,47

 

 

 

 

(10)

0,59

0,323

0

0,233

0,413

0,529

1,61

 

0,595

0,375

0,233

0

0,391

0,52

1,6

 

 

 

0,913

0,415

0,413

0,391

0

0,28

1,42

 

 

 

1,0

0,503

0,529

0,52

0,28

0

1,18

 

 

 

 

 

 

1,85

1,47

1,61

1,6

1,42

1,18

0

 

 

 

 

Можно утверждать, что согласованность полученной функции расстояния с метрикой в пространстве значений твердости образцов стала выше, чем в случае (9). В частности, упомянутое нарушение порядка следования элементов для образца № 7 удалось исправить. Значение расстояния Минковского в данном случае равно l1 = 4,35. Таким образом, использование более

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

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

470

508

Информационные технологии в управлении техническими системами и технологическими процессами

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

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

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

Исследование проводились при финансовой поддержке Министерства образования и науки РФ в рамках проекта RFMEFI58617X0055. Название совместного проекта с ЕС в рамках программы «Горизонт 2020»: «Fracture across Scales and Materials, Processes and Disciplines» Акроним проекта: «FRAMED».

Литература

1.БАНДИ Б. Методы оптимизации. Вводный курс: пер. с англ. –

М.: Радиоисвязь, 1988. – 128 с.

2.ВОЖАКОВ А.В., ГИТМАН М.Б., ФЕДОСЕЕВ С.А. Ком-

плексное оценивание при выборе оптимального плана производства на тактическом уровне с учетом нечетких крите-

471

509

Управление большими системами. Выпуск XX

риев и ограничений // Управление большими системами.

Вып. 30. – М.: Изд-во ИПУ РАН, 2010. – С. 164–179.

3.ЗАНГВИЛЛ У. Нелинейное программирование. Единый под-

ход; пер. с англ., под ред. Е.Г. Гольдштейна. – М.: Сов. радио, 1973. – 312 с.

4.Исследование процесса механотермической обработки арматуры из высокоуглеродистых сталей / М.П. БАРЫШ-

НИКОВ, Д.К. ДОЛГИЙ, К.Ю. КУРАНОВ, М.В. ЗАЙЦЕВА //

Сталь. – 2012. – № 2. – С. 89–97.

5.КУРДЮМОВ Г.В., УТЕВСКИЙ Л.М., ЭНТИН Р.И. Превра-

щения в железе и стали. – М.: Наука, 1977. – 238 с.

6.ПОТАПОВ А.С. Распознавание образов и машинное восприятие: общий подход на основе принципа минимальной длины описания. – СПб.: Политехника. 2007. – 548 с.

7.Разработка интеллектуальной системы распознавания сложных микроструктур на шлифах металлов и сплавов/ С.И. ША-

РЫБИН, В.Ю. СТОЛБОВ, М.Б. ГИТМАН, М.П. БАРЫШНИКОВ // Нейрокомпьютеры: разработка, применение. – 2014. – №12. – С. 50–56.

8.РАУШЕНБАХ Г.В. Меры близости и сходства // Анализ не-

числовой информации в социологических исследованиях. –

М.: Наука, 1985. – С. 169–203.

9.STOLBOV V.Y., GITMAN M.B., SHARYBIN S.I. Application of Intelligent Technology in Functional Materials Quality Control // Materials Science Forum. Materials Engineering and Technologies for Production and Processing II. – Vol. 870. – P. 717–724.

10.Shape Memory Effects in alloys / Ed. J. Perkins. – NY, London: Plenum Press, 1975. – 480 p.

11.SHARYBIN S.I., KLYUEV A.V., STOLBOV V.YU. Complex grain structure of metals and alloys // Steel in translation. – 2013. – Vol. 43. – No 5. – P. 245–248.

472

510