Файл: Горюшко, В. Е. Планирование эксперимента в бытовой химии [обзор].pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 01.11.2024

Просмотров: 40

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Г Л А В А IV

ПЛАНИРОВАНИЕ ЭКСПЕРИМЕНТА ПРИ ИССЛЕДОВАНИИ ДИАГРАММ СОСТАВ — СВОЙСТВО

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

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

Планирование эксперимента применительно к рецептурам раз­ вивается в двух направлениях. Первое из них предполагает опти­ мизацию состава смеси, второе — изучение диаграмм состав — свойство.

Для первого направления характерно исследование свойств сме­ си с помощью ортогональных факторных планов [7, 9]. Поскольку сумма компонентов равна 10 0 %, это условие приводит к коррелированности факторов. При математическом описании системы обычно опускают один компонент, а для остальных применяют полные и дробные факторные планы (см., например, [10]).

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

Факторное пространство представляет собой правильный симп­ лекс с

2 ^ = 1 , x i> °> * '= 1 , 2 , . . . , ? i=i

и q вершинами .в (q 1 ) -мерном пространстве, где х* — относитель­ ное содержание t-ro компонента в смеси [12—15].

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

63

// = *0+ 2

2 b n x i x i :

2

b i j b x i x i x b> (Iv- 1)

1 < i < q

K J < j < q

1< / < / < S < 7

 

а приведенными полиномами, получаемыми из (IV. I) с учетом ус­ ловия

ч

 

 

2

Х ;= 1 ’

 

/=1

 

со значительно меньшим числом коэффициентов.

факторов

Так, например, полином

второй степени для трех

имеет вид

 

 

УК \- blx l-f- b2x 2ф b3x3-j- Ьпх гх 2-(- b13x vx3-j- bn x 2x 34-

+-bnx\j -)- Ь22х\ Ь33х\.

(IV. 2)

С учетом уравнения

 

(IV. 3)

x ]-L-x.z-i-x3= l

этот полином можно преобразовать к приведенной форме

 

У — $1Х14~РаХ24" Рз-*3 4~ФиХ1Х24" ?13Х 1Х 3 4- $23Х 2Х 3-

(IV. 4)

При переходе к приведенной форме постоянный член Ь0 был ис­ ключен из (IV.2) путем умножения обеих частей (IV.3) на Ь0:

Кхг+ Ь0Х24- Ь0х 3= К

(IV. 5)

и подстановкой полученных результатов в (IV.2):

 

у = (^о4" )x i 4" (^о4" х 2 4" 4~^з) х з 4" ьпх 1 х з4" ^з-^’Л

b23x 2x 3-j- bnх\ Ь22х\ b33x\.

(IV.6 )

Квадратичные члены исключали, подставляя в уравнение зна­ чения X]2, х22 и Хз2:

xi --х1 ЧЛ2

Х ‘ = Х 2 Х хХ 2 х 2х 3,

*а = *з- -х,х.

1Х 3 Х 2Х 31

образованные умножением (IV.3) на Xj, х2 и х3 соответственно. Таким образом было получено выражение

у—{Ьо4~ )х1~\~{Ьз~\-^2 +^22)-^24“(А4 ^“г^ззАзЧ-

4 ( ^ 1 2 ^11

^ 22) Х \ Х 2 4 ~ А з

Ь \1 ^33> Х 1Х а + ( Ь 23 ^22 ^33 ) х 2-^з'

 

 

(IV. 7;

Используя

обозначения

= b0 + bi + 6 И, p2 = b0 + b2 + b22 и р3

= Ь0 + Ьг+ Ь-с3, придем к приведенной форме (IV.4).

64


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

Число экспериментов для построения приведенного полинома степени я от q факторов (компонентов) можно рассчитать по фор­ муле

К = Сп =

+ п ~ Ч!

q + n - 1

n \ ( q — 1)!

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

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

К _ тя{д + Ч

,

д(д — !)(<? —2)

2

1

6

Применив спмплекс-решетчатые планы с приведенными полино­ мами вместо обычных, получим экономию на числе экспериментов, определяемую по формуле [14]

 

С"

Сп

= Cq, ..

 

q + n

q + n — 1

q + n — 1

Например,

для трех

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

(у = 3, я = 2} экономия составит

 

 

С?:

4!

опыта.

 

3! (4 — 3)!

 

 

 

Шеффе [16]

предложил планы,

обеспечивающие равномерное

распределение экспериментальных точек по —1)-мерному симп­ лексу. Идея метода состоит в том, что в {д, я}-симплексной решет­

ке используются я+1 равнорасположенных значений для

каждого

фактора ( Хг = 0, 1/я, 2/я, ..., 1) и их комбинации.

 

Так, для приближения поверхности отклика полиномами второй

степени (я = 2) необходимо использовать = 0, 7г, 1, для

прибли­

жения поверхности отклика полиномами третьей степени

(я = 3) —

соответственно Xj== 0, 7з, 7г, 1.

Примеры симплексных {3, я}-решеток представлены на рис. 4. Анализ расположения экспериментальных точек на планах по­ казывает, что они являются последовательно строящимися (компо­ зиционными). Например, добавляя к квадратичной решетке одну точку в центре с координатами (7з, 7з, 7з), можно получить непол­ нокубическую решетку. Большая часть экспериментальных точек этой решетки может быть использована для решетки четвертой сте­

пени.

Удобно ввести специальные обозначения для откликов в экспе­

риментальных точках. Отклик для чистого

компонента в вершине

симплекса обозначим через г/г, отклик для

бинарной

(1:1) смеси

компонентов — через уц, отклик для тройной (1:1:1)

смеси — че­

рез г/гjh (i< j< k ), отклик для бинарных (2:1

и 1: 2) смесей — соот­

ветственно через г/ггj и г/зд (i<j) и т. п.

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


маковых индексов — на количество компонентов смеси, число оди­ наковых индексов — на повторение дозировки.

Так, отклик #12зз соответствует сочетанию рецептурных факто­ ров, имеющих координаты решетки четвертой или неполной четвер­ той степени.

Три различных индекса указывают на трехкомпонентную систе­ му с пропорциями */4 + lU= '/г Для компонента 3 и — Для каждо­ го из компонентов 1 и 2.

Рис. 4. Симплексные {3, п)-решетки.

Рассмотрим один из планов, в частности план, принятый нами для исследования моющей способности СМС:

опыта

М

х 2

-«3

У

1

1

0

0

>’1

2

0

1

0

У2

3

0

0

1

Уз

4

0,5

0,5

0

У12

5

0,5

0

0,5

У13

6

0

0,5

0,5

У23

7

0,5

0,25

0,25

У 1123

8

0,25

0,5

0,25

У 1223

9

0,25

0,25

0,5

У1233

10

0,333

0,333

0,333

У123

66


План составлен таким образом, чтобы использовать свойство композиционное™. По этому плану могут быть построены в трех­ компонентной системе модели состав — свойство второй, неполной третьей и неполной четвертой степени.

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

Подставляя в уравнение (IV.4) координаты точек х и х2 и Хз. лежащих в вершинах концентрационного треугольника, получим

соответственно у i = pb у2=р2 и Уз= Рз-

Подстановкой в это же уравнение координат точек, расположен­

ных на серединах сторон треугольника, можно получить

 

У12— У1 -7,—b&2 — + ?12— .

(IV. 8)

откуда с учетом yi~=$i

 

гЧ2 —4^12 ^У\

 

и т. д.

Аналогично определяют остальные коэффициенты модели. Модель первого порядка для трехкомпонентной смеси имеет вид

 

 

У=

+ Р-2-^2 +

 

 

 

(IV. 9}

Модель второго порядка:

 

 

 

 

 

 

 

 

y=i>1X1J\-$2X2JC$3X3~\~$nXlXZJC\J\3XlX3~ir$ZZX‘lXZ'

(IV- Ш-

где

 

1^12= fyiz

 

 

2?/2,

 

 

 

 

 

Из~ ^У\з ^Уг ^Уз'

 

 

 

 

 

А23= 4^/23

 

2 г/2

2у3.

 

 

 

Неполная кубичная модель:

 

 

 

 

 

 

 

У — $1Х1^$2Х2~\~$3Х3~\~$12Х1Х2~\~$13Х1Х3~\~ i'23X2X3JT$lZ3XlX‘2X3'

(IV- П)

где

?12з—27«/ш— 12 {Уп^У1 з ^ У 2з)^гЗ{у1-\-у2^ у 3).

 

Модель неполной четвертой степени:

 

 

 

 

У= ?1-^1 4~~ ^2-^2 4“ РзХз4"~ Pl2-^l-^2 4“ ИзМ^З iJ23"^2-^3 ~Ь

 

 

+

14123-^^3 4- ?122зМЛ:2Д:з + ^ гззМ ^ з’

 

(IV12)

где

2з~ 32 (Зг/1123

У1223

 

У1233)

3 (Зг/12-)- 3у53

 

 

 

 

У23 4~У\

У2

Уз)’

 

 

 

^1223“

3 2 (З г /122з

У1123

У1233)

 

' 8

( 3 //124 - 3 г/23

У1з~\~ У-г

У1

Уз)’

iJi233=

32 (Зг/1233

Уц2з

У1223)

 

®(3^1з4“ ЗУ2 3 '

У124”Уз

У1

Уг)-

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

67


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

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

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

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

При числе параллельных опытов в кажой точке, равном п, дис­ персия среднего значения может быть рассчитана по формуле [12]

(IV. 13)

Множитель | появляется в результате неортогональности плана и зависит от местоположения точки на плане.

Например, для модели второго порядка

 

 

5=-- 2

*?+ 2

 

4 ,

где

ai= x i {2xl— 1);

 

 

 

 

aij — AxiXj.

 

 

 

 

Изолинии |

в виде «контурных карт»

для тройных систем приво­

дятся в литературе [13—15].

 

 

 

Адекватность, как правило, проверяется по критерию Стьюдента

 

 

< V ;/ =

АУу П

>

(IV. 14)

 

 

--------------

 

 

 

г,

 

 

 

 

 

si/ У 1+ $

 

 

где

а — уровень значимости;

 

 

 

 

/ — число контрольных точек;

 

 

 

А= | jfaccn

У расч ] ■

модели принимается в случае tBксп<

 

Гипотеза об адекватности

< ^табл-

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

делью, можно рассчитать по формуле

 

Д = + 4 /* ;/4 М 1/2,

(IV. 15)

У П

 

где k — число коэффициентов модели.

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

Обозначим факторы (содержание ПАВ) через Х\, Х2 и х3. Матри­ ца и результаты эксперимента представлены в таблице:

68