1 of 27

Кацусика Хокусай (1760 –1849)

2 of 27

3 of 27

ВВЕДЕНИЕ в

компьютерное

материаловедение

5. СТРУКТУРИРОВАНИЕ СИСТЕМЫ ФАКТОРОВ.

СТРУКТУРИРОВАННЫЕ ЭС- МОДЕЛИ

Противоречие между «простотой»

интегральной полиномиальной ЭС-модели,

описывающей полное поле Y(x), в полной факторной области Ωx,

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

Декомпозиция системы факторов,

начиная с содержательной (технологической) постановки задачи

4 of 27

Априорное структурирование – до эксперимента,

при его «предпланировании» и планировании

РТ-факторы разделяют на подсистемы,

прежде всего, в соответствии с характером и степенью

их влияния на структуру и свойства материала.

Это могут быть:

группа факторов, определяющих свойства матрицы

(водо-вяжущее отношение, концентрации добавок и т.п.)

факторы, задающие количество наполнителя,

его дисперсионный и (или) минеральный состав

состав комплексной добавки

еще

Априорное структурирование – до эксперимента,

при его «предпланировании» и планировании

5 of 27

etc.

эксплуатационные факторы, задающие агрессивность среды

совокупности факторов из разных подобных групп,

к варьированию которых

наиболее или наименее чувствительны отклики Y

(по априорным данным)

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

Апостериорное структурирование – когда

ЭС-модели уже построены

может потребоваться при анализе

влияния РТ-факторов на свойства – тогда

могут быть выделены другие РТ-подсистемы,

любые полезные комбинации градиентных и изменяющих факторов

в конце темы

Апостериорное структурирование – когда

ЭС-модели уже построены

6 of 27

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

факторного пространства:

p-мерные кубы (отрезки, квадраты, кубы…)

взаимонезависимых нормализованных факторов xi

−1 ≤ xi ≤ +1, x = (x1, x2, …, xi, …, xp)

-1

0

+1

∙ ∙ ∙

7 of 27

0

+1

-1

x1

7

10

4

8

6

14

15

11

12

13

x3

3

9

1

5

2

-1

-1

+1

+1

x2

«технологическое» структурирование

8 of 27

(q-1)-мерные симплексы (отрезки, треугольники, тетраэдры…)

q линейно связанных смесевых факторов vi,

0 ≤ vi ≤ 1, Σvi = 1, v = (v1, v2, …, vi, …, vq);

в системе может быть выделено две смеси

(w – второе обозначение)

1

0

1

0.5

0

вырезанные дополнительными ограничениями на x, v, w

подобласти кубов и симплексов,

а также призмы и другие «произведения» полных

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

9 of 27

Точки плана эксперимента при технологических ограничениях на содержание суперпластификатора (x1) и микрокремнезема (x2)

и гипотезе

об эффектах 4-го порядка в ЭС-модели

x1

x2

-1

1

1

-1

∙ 1

∙ 2

∙ 4

∙ 7

∙ 5

∙2

8

∙ 9

∙ 6

∙ 3

10 of 27

Отображения

факторных областей

в 3-мерном пространстве «структура вектора факторов»

См. стр. 69-71 в книге

«ЭВМ и оптимизация композиционных

материалов»

/ В.А. Вознесенский,

Т.В. Ляшенко,

Я.П. Иванов,

И.И. Николов.

К.: Будивэльнык,

1989. – 240 с.

https://drive.google.com/file/d/0BzKYSjvwhyieS3g0NFd5RjUyOEE/view

11 of 27

12 of 27

13 of 27

14 of 27

Подсистемы факторов, определяющих набор свойств материала – его качество Q (от quality), обозначены:

T (от technology) – «технология»,

Tp – подсистема p взаимонезависимых x,

из Tp могут быть выделены несколько подсистем «технология» (в частности, Tp1 и Tp2, p1 + p2 = p);

M (от mixture) – «смесь»,

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

Mq – смесь q компонентов

Это и факторы рецептуры

(например, дозировки добавок

относительно базового компонента),

и параметры процессов производства и эксплуатации

15 of 27

Особенности пространства факторов

в задачах разработки материалов

могут приводить к рассмотрению систем

«смеси, технологии – свойства» (MMTTQ).

Вектор факторов для них в общем случае

можно записать в виде z = (v, w, (x1, x2)

Система факторов «смеси, технологии» (MMTT) обобщает все рассматриваемые

в строительном материаловедении

комбинации линейно связанных

и взаимонезависимых переменных,

включаемых в ЭС-модели

16 of 27

В таблице ниже обозначен

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

(целесообразно до 7 факторов)

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

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

Некоторые обозначения: M3T2, M3M3T1, T1T3

17 of 27

Обозначения факторных систем «смеси, технологии»

с вариантами отображения факторных областей

См. книгу [1], с. 49

18 of 27

Структурированным системам факторов соответствуют и структурированные ЭС-модели

Для систем TT это обычные полиномы, записанные блоками, соответствующими отдельным подсистемам и их синергизму – как модель для вязкости η1 (Па⋅с) фиброполимерцементной композиции

(при скорости сдвига γ′ = 1 с-1)

в зависимости от содержания диспергируемого полимера (Vx1), молекулярной массы (MMCx2) и дозировки (MCx3) метилцеллюлозы, количества полипропиленовой фибры (Fx4)

Vinnapas effects at median MMC, MC, F

Effects of MMC & MC

at median V & F

Synergism

Fibre effects at median V, MMC, MC

19 of 27

Линейная связь факторов приводит к специальному виду полиномиальных моделей.

Это – приведенные полиномы Шеффе (H. Sheffe, 1958, 1963)

Число эффектов и соответствующих коэффициентов

в приведенном полиноме степени m от q факторов v

равно Cq+m-1m,

меньше, чем Ck+mm в обычном полиноме,

т к = q взаимонезависимых x)

Так приведенный полином 2-й степени (с коэффициентами A)

для тройной системы имеет следующий вид

(6 коэффициентов вместо 10-ти)

20 of 27

В общем случае –

приведенный полином 2-й степени от q перемнных

Коэффициенты имеют физический смысл:

Ai равны Y (уровню свойства) для «чистого» i-го компонента,

нелинейная часть – синергизм

(положительный, если смешивание компонентов вызывает

увеличение отклика, Aij > 0,

отрицательный – антагонизм, при уменьшении отклика, Aij < 0)

21 of 27

Приведенные полиномы Шеффе неполной 3-й и 3-ей степени

Коэффициенты приведенных полиномов также определяются методом наименьших квадратов (МНК),

по значениям факторов и откликов

в спланированном эксперименте

22 of 27

Для систем со смесями предложены «приведенные» полиномы

Например, для системы M3T2

приведенный полином 2-го порядка

описывает влияние фракционного состава серпентинитового заполнителя (w1, w2, w3 – доли «мелкой», «средней» и «крупной» фракций), концентрации СП и В/Ц (взаимонезависимые x1 и x2)

на прочность (МПа) серпентинитового бетона

(а) полином Sheffe, влияние многофракционного заполнителя при x4 = x5 = 0 (в центре диапазонов факторов матрицы), коэффициенты при wi уровни R для отдельных фракций, другие оценивают эффект от их смешивания, их синергизм

(b) ) учитывает нелинейность влияния факторов матрицы при любом фракционном составе серпентинита

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

23 of 27

«Произведения полиномов»

если анализ априорных данных не позволяет выдвинуть гипотезу о соответствии зависимости Y (z) приведенному полиному,

например, об отсутствии эффектов vivjxk

В частности,

произведение полиномов 2-го порядка (Шеффе и обычного) для системы M3T1

Или

Но слишком дорогое удовольствие

(18 коэффициентов, не менее 18 опытов)

24 of 27

10 коэффициентов, не менее 10 опытов

Но не позволяет, в частности, выделить эффекты vivjx – охарактеризовать связь синергизма внутри смеси с уровнем «технологического» фактора

Через модель-произведение можно соединить с ЭС-моделью «квазифундаментальную» модель

Параметры K и m модели Оствальда-Виля, зависимости вязкости от скорости сдвига,

(где K – вязкость η1, при единичной скорости, m – темп разрушения)

могут быть выражены функциями состава

Слайд 18

25 of 27

При апостериорном структурировании,

при анализе результатов моделирования

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

Могут быть рассмотрены любые полезные комбинации выделения

градиентных и изменяющих факторов из полного вектора факторов.

Перемещение локального поля Y(xgr) по области изменяющих факторов

Впереди анализ и оптимизация – примеры решения задач

Может оказаться полезным:

сравнение локальных полей Y(xgr)

для «контрастных» уровней xch,

с наибольшим перепадом свойств

в диапазоне изменяющих факторов;

графики и «карты» (изолинии, поверхности) локальных полей при разных уровнях изменяющих факторов, etc.

26 of 27

27 of 27

The End

The End