Файл: Горюшко, В. Е. Планирование эксперимента в бытовой химии [обзор].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