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

1185

.pdf
Скачиваний:
6
Добавлен:
15.11.2022
Размер:
11.27 Mб
Скачать

Таблица 7.10

Планирование опытов при изучении процесса хлорирования титановых шлаков

Основной уровень

16,6

750

3,5

2

5,65

1 , 0

8,3

1

6,65

5,7

4,7

7

7

8

1 : 1

Интервал варьиро-

5,2

50

1,5

1

1,35

0,5

1,7

1

1,35

1,3

1,3

3

3

2

0,5*

вання

 

800

 

3

7,0

 

 

 

 

 

7,0

 

 

 

 

1:0,5

Верхний уровень

2 1 , 8

5

1,5

1 0

, 0

2

8 , 0

6 , 0

1 0

1 0

1 0

Нижний уровень

11,4

700

2

1

4,3

0,5

6

, 6

0

5,3

4,4

3,4

4

4

6

1:1,5

Кодовые обозначения переменных

Помер опыта

3CQ

 

X;

 

х,

*5

 

 

 

Xs

 

1

+

+

 

 

 

+

+

+

+

2

+

+

+

3

+

 

4-

+

 

4

+

+

4-

+1

5

+

+

 

+

“Г

+

6

+

+

+

7

+

+

+

+

+

+

8

+

+

+

+

+

+1

9

+

 

+

 

+

+

+

1 0

+

+

+

 

И

+

4~

+

+

+

+

1 2

+

+

 

+

+

+

— .

13

+

+

+

+

+

14

+

+

+

+

+

15

+■

+

+

+

+

 

+

+

+

16

+

+

 

 

+

+

+

Ь

75,25

8,69

3,75

1 , 2 0

28,79

-5,09

8,38 7,40 1,29

 

 

 

 

 

 

 

 

 

X,

* 1 0

Хц

х12

х„

Xu

*15

Vi

V:

3

+

+

+

+

__

_

4

39,2

48,5

43,85

+

1 1

4

48,0

47,5

47,75

 

+

+

-1_

 

4

44,7

48,5

46,6

+

 

4

4

1

47,2

43,7

45,45

 

“Г

+

 

+

4

4

41,7

33,4

37,55

+

+

4

4

40,5

48,2

44,35

 

+

1

4

48,2

49,2

48,7

+

1

49,4

49,5

49,45

4

4

98,5

95,0

96,75

1

+

4-

4

4

115,0

92,5

103,75

+

4

1 1 0 , 0

102,5

106,25

 

+

+

4

4

50,4

47,0

48,7

+

 

4

151,0

151,1

151,05

+

+

41

93,5

103,7

98,6

+

4-

+

4

4

~Т~

4

127,0

123,5

125,25

 

 

4

4

85,0

103,0

94,0

1 , 6 6

1,65 -5,12

—3,71

7,78

0 , 6 8

-4,99

 

 

 

* Интервал относится ко второму члену пропорции.

факторов от 32 до 63 (здесь необходимы 64 опыта) и т. д. Однако для решения столь сложных задач рекомендуется применять методы отбора факторов, например метод случайного баланса [5].

Поэтому в «Ограничениях» мы и отметили, что будем рассмат­ ривать задачи с числом факторов от 2 до 15.

В качестве последнего примера рассмотрим применение плана 215-11.

Пример 5. Изучалось влияние следующих 15 факторов на скорость хлорирования титановых шлаков в расплаве [6 ]. Пусть — расход хлора; х2 — температура; х3 — концентрация углерода в расплаве; 2 4 — концен­ трация ТЮ2 в расплаве; хь—хц — компоненты, определяющие состав шлака; х12х1Ъ— компоненты расплава.

Столь большое число факторов было включено в программу исследова­ ния потому, что процесс хлорирования плохо изучен. Планы такого типа применяются в планировании эксперимента довольно редко, так как число степеней свободы равно нулю и невозможно провести последовательный статистический анализ. Здесь даже нет смысла выписывать совместные оценки, линейные эффекты смешаны со 105 парными взаимодействиями. Освободиться от них можно с помощью метода «перевала»: нужна дополнить первую реплику второй, изменив все знаки в матрице планирования на об­ ратные. В табл. 7.10 представлены уровни факторов и интервал варьирова­ ния, матрица планирования и результаты экспериментов.В последней строке помещены коэффициенты.

7.6. Резюме

Дробные реплики находят широкое применение при получе­ нии линейных моделей. Целесообразность их применения воз­ растает с ростом количества факторов. В табл. 7.10 показано, что при исследовании влияния 15 факторов можно в 2048 раз сокра­ тить число опытов, применяя реплику большой дробности (16 опытов вместо 32768). Эффективность применения дробных реплик зависит от удачного выбора системы смешивания линейных эф­ фектов с эффектами взаимодействия, а также от умелой-стратегии экспериментирования в случае значимости некоторых взаимодей­ ствий. Априорные сведения о взаимодействиях могут оказать большую услугу экспериментатору.

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

Реплики, которые используются для сокращения опытов в раз, где т = 1, 2, 3, 4, . . ., называются регулярными. Они пользуются большой популярностью, так как позволяют про­ изводить расчет коэффициентов уравнения так же просто, как и в случае полного факторного эксперимента.

111

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

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

^ = £ 22:3, х2= х 1х3

и х3= х гх2,

то определяющий контраст будет

1 = х 1х2х3.

реплики

зависит от системы смешивания.

Эффективность

Реплики, у которых линейные эффекты смешаны с взаимодей­ ствиями наивысшего порядка, являются наиболее эффективными,

так

как обладают наибольшей разрешающей способностью.

Г

Для освобождения линейных эффектов от взаимодействий

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

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

трастов.

Мы научились строить полные и дробные факторные экспери­

менты. Давайте теперь посмотрим,

как их реализовать.

Л и т е р а т у р а

1.W. S. Connor, S. Young. Fractional Factorial * Designs’’ for^ Experiments with Factors of two and three Levels. National Bureau of Standards, Applied

Mathematics Series, 58, 1961.

2.0. L. Davies, W. A. Hay. The construction and uses of fractional factorial designs in industrial research. Biometrics, 6 , 1950.

3.B. 3. Бродский. Многофакторные регулярные планы. М., изд-во МГУ, 1972.

4.Е. В. Маркова, П. Фарка, И. П. Борисова и др. Метод планирования эк­

сперимента при поиске оптимальных условий получения полимерного серусодершащего антиоксиданта. — Пластические массы, 1969, № 1 .

5. F. Е. Satterthwaite. Random Balance Experimentation. — Technometrics, 1959, 1, N 2.

6 . П. П. Хомяков, Ю. П. Адлер, В. В. Налимов. Выявление фактиров, влияю­ щих на скорость хлорирования титановых шлаков в расплаве. — Завод­ ская лаборатория, 1963, 29, № 1.

Глава восьмая

ПРОВЕДЕНИЕ ЭКСПЕРИМЕНТА

Перед ошибками захлопываем дверь.

Всмятенье Истина: «Как я войду теперь?»

Р. Тагор. Искри

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

тесь

на землю,

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

все

детали постановки опыта.

В этой главе

рассказывается, как нужно готовиться к опыту,

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

8.1. Анкета для сбора априорной информации

Постановка задачи, выбор параметров оптимизации

1.Краткое описание процесса, объекта.

2.Формулировка цели исследования (если задач несколько — проранжировать их по степени важности).

*3. Выбор параметров оптимизации (откликов).

Заполните следующую таблицу, включив в нее все возможные отклики

Номер

Название

Размерность

Область

Точность

Примечание

отклика

определения

8 Заказ. М 588

из

4.Желаемый результат. Число и точность.

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

Выбор факторов

1.Список всех «подозреваемых»: факторов, которые могут влиять на процесс.

2.Список факторов, включаемых в реальный эксперимент.

Номер

Название

Размер­

Область

Область

Точность Примечание

фактора

ность

определения

интереса

3.Существуют ли возможности установления значения фак­ тора на любом заданном уровне?

4.Сохраняются ли заданные значения уровней в течение опыта?

5.Могут ли некоторые комбинации уровней факторов при­ вести к остановке процесса (например, взрыв, нетехнологичность

ит. д.)?

Число опытов

1.Желаемое число опытов, ограничения на число опытов.

2.Желаемый срок проведения исследования.

3.Примерная длительность одного опыта.

4.Стоимость и затраты труда при проведении одного опыта

серии.

5.Желаемое число уровней для одного фактора.

6. Возможность выполнения параллельных опытов и их же­ лаемое число.

7. Возможность проведения параллельных, измерений.

8. Желаемая стратегия проведения опытов (например, по од­ ному в день и т. д.).

Учет априорной информации

1.Условия и результаты, достигнутые при изучении аналогич­ ных процессов.

2.Результаты предварительного эксперимента и данные (ли­ тературные и собственные) о величине ошибки эксперимента.

3.Взаимодействия факторов.

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

114

8.2. Реализация плана эксперимента

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

Чтобы облегчить работу лаборанта и исключить ошибки при

выборе условий опыта, в рабочей матрице

планирования целе­

сообразно проставлять не только кодовые

значения факторов,

но и натуральные.

 

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

результатов

эксперимента.

Все

расчеты должны сохраняться

до окончания

работы.

 

 

Пример 1. В качестве примера приведем оформление журнала при опти­

мизации процесса

получения сульфадимизина.

 

 

 

 

Страница 1

Планирование

эксперимента

при

оптимизации процесса получения

сульфадимизина

 

 

Цель исследования: определение оптимальных условии процесса

конден­

сации сульгина о ацетилацетоном в присутствии уксусной кислоты.

 

Параметры,

характеризующие процесс: ух — выход •уЛьфадммизжна

по сульгину, %,

уа — качество сульфадимизина.

 

 

8*

115

Формулировка аадачи оптимизации. Достижение максимального выхода сульфадимизина, качество которого удовлетворяет требованиям фармакопеи: Pi — параметр оптимизации (ух -> 1 0 0 %), — ограничение.

Качество продукта определяется по процентному содержанию сульфа­ димизина в получаемом продукте и по его температуре плавления. Согласно требованиям фармакопеи содержание основного вещества в получаемом про­ дукте должно быть не менее 99%, а температура плавления должна нахо­

диться

в пределах 196—200° С,

т. е. у£ > 99%, 196* С <

у\ <

200е С.

 

 

 

 

 

 

 

Страница 2

Факторы, определяющие процесс: 5Х— время реакции,

час; 52

— содрр

жание

ацетилацетона в

реакционной массе, %; Ж9 — содержание

уксусной

кислоты в реакционной

массе,

%; 54

— температура реакционной

массы,

°С;

— качество ацетилацетона,

%;

5 , — качество

сульгина,

%.

 

Выбор варьируемых

факторов.

Принято решение

изменять

 

в

опытах

первые три фактора. Качество ацетилацетона и сульгина решено поддержи­ вать постоянным (таким, как на действующем производстве) х5= 9 0 % ; 5а=98% .

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

Страница 3

Выбор технологии. Сульгин загружается одновременно с ацетилацетоном и уксусной кислотой. Реакция проводится при перемешивании реакцион­ ной смеси и непрерывном отгоне воды.

Необходимые анализы. Анализ исходного сырья: ацетилацетона, суль­ гина, уксусной кислоты (следует описание методик). Анализы получаемых продуктов: сульфадимизина в осадке, сульфадимизина в фильтрате, сульгина в фильтрате (следует описание методик).

Описание экспериментальной установки. Опыты проводятся на лабора­ торной установке, состоящей из стеклянной конической колбы емкостью 250 мл, снабженной металлической якорной мешалкой и обратным холодиль­ ником. Температура реакционной массы измеряется термопарой, подключен­ ной к электрическому потенциометру, и непрерывно записывается на карто­ грамму. Колба обогревается электрической баней, наполненной вазелиновым маслом. Температура в бане автоматически регулируется с помощью реле и контактного термометра и поддерживается около 160° С.

Страница 4

Выбор основного уровня и интервалов варьирования. Для того чтобы выбрать уровни факторов, следуем собрать и проанализировать литератур­ ные и заводские данные. По заводскому регламенту процесс проводится при следующих условиях: 2Х=27 час, £а= 4 % , 2Я=16% . При этом уг= 84%.

Опубликованные данные и сведения из отчетов (априорная информация).

Влияние времени реакции (хх). Данные об оптимальном времени реакции в лабораторных условиях противоречивы. Так, в отчете № 1 указано, что опыты проводились при 2Х=21 час, затем время уменьшили до 12 час. Умень­ шение времени не снизило существенно выход реакции. В отчете № 2 описы­

116

ваются опыты с различным временем: 18, 24 и 30 час. Наилучший выход получен при * х= 2 4 час.

Влияние избытка ацетилацетона (х2). Данные о влиянии избытка аце­ тилацетона от стехиометрического соотношения также противоречивы. В одном отчете указано, что содержание ацетилацетона сверх 10% является нецелесообразным, в другом оптимальным считается 40% избытка ацетил­ ацетона.

Влияние уксусной кислоты (х3). Вопрос о влиянии процентного содержа­ ния уксусной кислоты специально не исследован. Считается, что х.л целесо­ образно поддерживать около 16—17%. Предполагается, что это раствори­ тель и его концентрация может изменяться в широких пределах.

На основании анализа имеющихся сведений решено выбрать следующие уровни и интервалы варьирования факторов (табл. 8 .1 )

Таблица 8.1 Уровни факторов и интервалы варьирования

 

 

 

Уровни

 

I

 

Факторы

 

—1

0

+1

Размерность

 

 

 

 

хх

 

16

18

20

2

час

 

 

20

24

28

4

%

*3

 

12

15

18

3

%

Расчет

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

 

 

Страница 5

 

 

 

Ацетилацетон

 

Сульгин

 

Уксусная кислота

сна- с о

 

NHa

 

 

О

I

 

+

С—NHSOa— С„Н4—NHa

Н3С-С<^он

сн

----►

I

 

I

 

 

 

СНз-СОН

 

HN

 

 

 

100

 

 

2 14

 

 

60

Сульфадимизин

 

 

 

 

 

СНа—C =N

 

 

 

 

 

СН C—NHSOa—CeH4—NHa+2H aO.

 

 

СН8 —С—N

 

 

 

 

 

2 78

 

 

36

 

 

 

Расчет по стехиометрическим соотношениям на 12,5 г ацетилацетона:

100—214

 

 

214 • 12,5 = 26,7 г;

 

 

 

12,5 — z

 

z —

 

 

 

 

100

 

 

 

ацетилацетона — 12,5

а, сульгина — 26,7

а.

 

 

Расчет

ацетилацетона с избытком 20%:

12,5+12,5-0,2=15 г и т. д.

117

Страница 6

Таблица 8.2 Матрица планирования и результаты эксперимента

Номер

Случай­

Дата

ный по­

опыта

рядок

проведе­

в мат­

реализа­

ния опы­

рице

ции опы­

тов *

 

тов

 

1

2

10/Ш

 

9

5/IV

2

6

2 1

 

13

14/IV

 

 

 

*3

 

 

St

 

St

 

 

К

О

*

О

 

 

о

к

к

 

 

X

 

16

 

2

0

 

+

2 0

___

2 0

__ _

3

1

8/1II

16

+

28

 

15

21/IV

 

 

 

 

 

4

7

27/Ш

+

2 0

+

28

 

1 0

6 /IV

 

 

 

 

 

5

3

1 2

16

_

2 0

+

 

16

27/IV

 

 

 

 

 

6

8

1/IV

+

2 0

___

2 0

+

 

14

16/IV

 

 

 

 

 

 

14'

17/IV

 

 

 

 

 

7

4

15/III

___

16

+

28

+

 

1 2

12/IV

 

 

 

 

 

8

5

20/1II

+

2 0

+

28

+

 

1 1

9/IV

 

 

 

 

 

 

 

У\,

0/

 

 

/0

 

 

повтор­ ные опыты

средний результат

1

2

80,23

81,08

 

 

81,93

 

1

2

86,50

85,65

 

 

84,80

 

1

2

82,45

82,27

 

 

82,10

 

1

2

89,50

90,40

 

 

91,30

 

1885,10 84,95

84,80

18

90,30

89,95

 

 

89,60

1885,60 85,25

84,90

18

8 8 , 0 2 88,25

 

188,48

* Начало проведения всех опытов в 8 ч. В таблице необходимо указать фамилии ла­ борантов по сменам. Для каждого опыта дать картограмму температуры.

Имея такую таблицу, можно приступить к обработке экспе­ риментальных данных^

Каждая горизонтальная строка матрицы — это условия опыта. Условия опытов чрезвычайно разнообразны. Ведь мы занимаемся планированием многофакторного эксперимента, когда все фак­ торы изменяются одновременно. Приступая к планированию экс­ перимента, мы должны отказаться от привычного однофакторного эксперимента, который проводится по принципу «изменяй один фактор, а прочие держи постоянными».

Мы хотим заниматься исследованием сложных многофактор­ ных систем и понимаем, что однофакторный эксперимент нам не поможет. И это не только наше мнение. У. Р. Эшби во «Вве­ дении в кибернетику» писал: «Тот факт, что в течение столетий могли принимать такую догму, нак «изменяйте факторы по од­ ному», показывает, что ученые занимались в основном исследова-

118

Страница 7

 

9

°/.

-------------тг

° с

 

 

 

Примечания

 

 

1/2.

 

1/2.

 

 

 

 

(

 

5

(

 

Н

 

 

 

 

 

 

 

— *

 

 

 

 

 

Р*

ч

« н

 

 

« н

 

 

 

 

 

S л

О

g

S А

 

начало

 

начало

выпадение

О

g

и (=:

Нч

 

 

о 3 С

м >■

5

© з

н > »

 

кипения

 

помутнения

осадка

QiO

О 3 С

SS

 

 

 

 

 

С И о

о п

с т о

О Р .

 

 

 

 

 

94,2

97,05

196.

196,5

8

ч 45 мин

2 0

ч 2 0 мин

20 ч 40 мин

99,9

 

197

 

8 ч 38 мин

2 0

ч 18 мин

20 ч 35 мин

99,6

99,80

197

196,5

8

ч 38 мин

2 0

ч 2 0 мин

20 ч 40 мин

1 0 0 , 0

 

196

 

8

ч 40 мин

20

ч 23 мин

20 ч 43 мин

99,2

99,60

195

195,5

8

ч 32 мин

20

ч 40 мин

20 ч 55 мин

1 0 0 , 0

 

196

 

8

ч 30 мин

20

ч 45 мин

2 1 ч

99,8

99,85

196

196,5

8

ч 2 0 мин

2 0

ч

20 ч 32 мин

99,9

 

197

 

8

ч 25 мин

20

ч 15 мин

2 1 ч

99,5

99,65

197

197,5

8

ч 25 мин

19

ч

19 ч 15 мин

99,8

 

198

 

8

ч 25 мин

19

ч 05 мин

19 ч 18 мин

99,6

99,65

196

196,0

8

ч 30 мин

2 0

ч 1 0 мин

20 ч 35 мин

— При кипении из колбы сильно выбрасывало

99,7

 

196

 

8

ч 2 0

мин

99,5

99,55

197

196,5

8

ч 25 мин

99,6

 

196

 

8

ч 28

мин

99,5

99,60

197

197,0

8

ч 35

мин

99,7

 

197

 

8

ч 32

мин

жидкость * *

to о Л

18 ч 15 мин

18 ч 45 мин

20 ч 45 мин

2 1 ч

2 0

ч 2 0 мин

18

ч 35 мин

19

ч

2 1

ч

21 ч 45 мин

* * Опыт повторить.

нием систем, допускающих этот метод, ибо в сложных системах он часто неприменим по существу».

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

рименте,

варьируя одним фактором и стабилизируя все прочие

на произвольно выбранных уровнях,

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

зависимость

параметра

оптимизации

только

от

одного фак­

тора и

определяет локальный оптимум. Далее

он повторяет

аналогичную

процедуру

для второго, третьего

и к-то фактора.

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

За время эксперимента могут происходить изменения в аппа­ ратуре, сырье и т. д. Все это вносит изменения в результаты

119

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