Файл: Баженов, Ю. М. Перспективы применения математических методов в технологии сборного железобетона.pdf

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

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

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

Добавлен: 29.10.2024

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

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

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

Матрица планирования (п кодовом обозначении) н результаты определения р (по Зг) приведены в табл. IV.7.

Т а б л и ц а

IV.7. Матрица планирования и результаты эксперимента

Номер опыта

Код

!/„=P

1 Номер опыта

Код

»«=p

1

abode

84,5

9

асе

85,6

2

abc

89

10

bde

123,9

3

abd

156,4

11

bee

88,2

4

acd

128

12

a

107,9

5

bed

125,8

13

b

122,5

6

abc

112,7

14

C

112,4

7

ade

109,7

15

d

105,8

8

ede

95,7

16

e

90,8

По формулам (1V.56—IV.58) при N = 1 6

рассчитываются

ко­

эффициенты регрессии, в результате чего модель имеет вид:

 

р = 108,68 — 0 ,54*і +

4,19*а — 7 ,53*3 +7,54.*,— 9 ,79*5 —

 

2,72х1хі -- 4,91*1*3 + 2 ,88*1*4 1 ,30*і*в — 8,46*г*3 +

 

+ 2 ,23*2*4 — 0 ,75*2*6— 0,19*3*4 — 2 ,85*3*5 — 2 ,98*4*s .

 

При а = 5°/о и / —40 определяем по (IV.60)

значение 6КР;

 

з( р)

8,1

 

(б)

Ь*р = і

= 2, 02 — ^ 1 = 4,09,

V N

V 16

 

 

если Ьі или bij меньше йкр, то такие эффекты можно удалить и по­

лучить модель:

 

р = 108,7 +

4,2*2 — 7,5*з + 7,5*4 — 9,8*5 —

 

 

— 4 ,9*і*з — 8,5*2*з •

(в)

Адекватность модели проверяется по F-критерню. Для этого пред­

варительно

рассчитаем

по (IV.61) сумму квадратов

неадекватности

S S HA иа

(суммирование только по значимым эффектам!):

ssHа =

S у і - ь0 (ОУ)- S bt m -

S b (ÜY) =

 

 

U=1

1=1

t-fj

 

=

194 796 — 188 984 — 2 485 — 2 680 = 647;

(r)

 

■S-Sh a ^ ha _647:(16 — 7)

1,096.

(Д)

 

s" (p)

65,61

 

 

 

Поскольку F ф <C Ftа б л === 2,1

(a = 5%, fi= 9 ,

/2= 4 0 ), модель

(в)

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

149



Переход к четырем моделям (е—и) для

каждой пары

добавок

осуществляется при подстановке в модель (в) значений *3

и х4

па

уровнях

±1 (например, для

добавок NH4C l+N a2S 0 4 подставляем

хз— +1

и х4 = + 1):

 

 

 

 

 

р {NH4C1 +

Na2S04} =

108,7 — 4,9*,;— 12,6*3 — 9,8х6;

(е)

 

p{NH4Cl +

K2C03) = 9 3 ,6 — 4,9*1— 12,6*3 — 9 ,8хв;

 

(ж)

 

р {КОИ +

Na2S04] =

123,8 +

4,9*! +

4,3*2 — 9,8*6;

 

(з)

 

р {КОН +

К2С 03) =

108,7 +

4,9*! +

4,3*2 — 9,8*6.

 

(и)

По моделям (е—и) для каждой пары добавок можно уточнить соотношение между добавками (такое *і, чтобы р увеличивалось, например, для p{NH4Cl+Na2S 0 4} * і = — 1) и уровень концентрации

—х2. Результаты оптимизации представлены в табл. ІѴ.8.

Т а б л и ц а IV.S. Оптимальные уровни добавок в примере ІѴ.10

Добавка

Уровень добавки

Прирост р,

(%) при Ц / В

А,: А,

С. %

1,6

2,2

 

KOH+NaoS04 ..................

0,7

2,5

143

128

NH4C l+N a2S 0 4 . . . .

0,3

0,5

136

116

КОН+КгСОз . . .

0,7

2,5

128

108

МН4СІ+Н2СОз . . . .

0,8

0,5

122

102

Из табл. ІѴ.8 можно сделать, в частности, такие выводы:

при

а)

эффективной является добавка KOH+Na2S 0 4, которая

любых

Ц/В (в пределах эксперимента) увеличивает прочность

не

менее чем на 2s{р} = 16,2%;

 

 

б)

максимальное повышение прочности р=

143% после двух су­

ток твердения достигается при Ц/В = 1,6 для

комбинированной

до­

бавки КОН : Na2SO4= 0 ,7 : 0,3 при общей концентрации 2,5%;

 

в)

все четыре парные добавки повышают эффективность при пе­

реходе к меньшим Ц/В.

 

 

ІѴ.8. Планы для построения квадратичных моделей.

Для определения коэффициентов регрессии Ьц при квад­

ратичных членах в модели (ІѴ.44) нужен такой план эксперимента, в котором каждая переменная будет варьироваться хотя бы на трех уровнях. Такое плани­ рование может быть получено путем добавления некото­ рого количества точек к «ядру», в котором уже были проведены опыты при линейном приближении. Эти но­ вые точки располагаются некоторым оптимальным об­ разом (с точки зрения целей данного эксперимента [56]). Последовательно строящиеся планы называются композиционными.

150


На рис. IV.6 показана схема планирования для трех факторов. После постановки ПФЭ в вершинах куба к плану добавляется шесть точек С, вынесенных за грани куба на некоторое расстояние а ^ І от центра эксперимента, где также поставлены дополнительные опыты. Точки С называются звездными точками.

Исследования в области математической теории экс­ перимента показали, что критерий ортогональности пла­ на для получения модели второго порядка не является

п

л

Рис. ІѴ.6. Трехфакториын план второго порядка, полученный добавлением шести звездных точек С к ПФЭ 2Э (точки А)

достаточно эффективным, поскольку коэффициенты Ь0, Ьі, Ьц и bij определяются с разными дисперсиями, вели­ чина которых меняется в зависимости от поворота ко­ ординатных осей факторного пространства. Было пред­ ложено (Бокс и Хантер, 1957) считать оптимальным планирование, при котором количество информации, со­ держащееся в уравнении регрессии, будет одинаково для всех точек, равностоящих (радиус г) от центра эксперимента. Такие планы были названы ротатабель­ ными [55]; для них в табл. IV.9 показаны значения звездного плеча а и другие необходимые характеристики.

Ковариационная матрица для планов второго порядка в общем случае недиагональна, но ее можно предста­ вить в виде четырех подматриц:

[Дх]

(IV.67)

|0] [Да!

 

из которых две нулевые [0 ], одна — [Дг] для оценки ли­ нейных эффектов и взаимодействий диагональна с эле­ ментами Гз и Т6 и одна [Ді] для оценки bQи квадра­ тичных эффектов содержит внедиагональные элементы

151

Г2 и Т5, определяющие корреляцию между Ь0 и Ьц, а также между è,-,- и &<>

 

7\ Го • • • т 3

 

[Ді] =

т,тп .. • т б

(IV. 68)

 

 

 

т ат 8 •• ' Т п

 

 

длл

-----------■--------

 

 

оценки

для оценки

 

Т3

0

[Д.

0

0

0

• • • 00

т3

о о

(IV.69)

0

 

•••То 0

0 . . . ° т е

для для оценки

оценки ьа

ь,

Отсюда следуют расчетные формулы для определе­ ния оценок коэффициентов регрессии:

Ь0 = Тх(ОУ) - Та£ («Г), s

 

 

 

(ІѴ.70)

 

 

 

і= і

 

 

 

 

 

 

ь, =

т3 (iY),

s {b,} = ) / T 3.s3 =

T8-s3;

(IV.71)

Т а б л и ц а

IV.9.

Ротатабельные

композиционные планы

( п а — число

звездных

точек,

/і0 — число опытов

в

центре

плана)

Число факторов

 

2

3

'1

 

5

б

Первый этап

 

ПФЭ-22

ПФЭ-23

ПФЭ-24

Полу-

Полу-

 

 

 

 

 

 

реплика

реплика

 

 

 

 

 

 

25—1*

20—г*

Второй этап

( а

1,414

1,682

2

 

2

2,378

< па

4

 

6

8

 

10

12

 

 

5

 

6

7

 

6

9

Общее число

точек

13

 

20

31

 

32

53

* Из прнл. ѴШ.

152


к

b„ = Т, (HY) Г5 У) (HY) - T2 (07); i=l

 

Ti = Tn

T5;

5 {^(i}

~ У'т'и S3 — T'os3;

Ьц = 70 (t/7);

S {ö;/} -

K 70s3 = 7\0s3;

P{Mn} = / 7 V / 7 y ,

p { ^ £-

= V

n -•

(IV.72)

(IV.73)

(IV.74)

(IV.75)

(IV.76)

(IV.77)

Для планов в табл. IV.9 значения расчетных коэф­ фициентов Ті приведены в прил. IX.

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

Для К = 2 рекомендуются [63] наименьшие рота­ табельные планы типа правильного пятиугольника и шестиугольника с центральными точками (/г0 ^ 1 ) . План ' на шестиугольнике показан на рис. IV.7 (расчетные ко­ эффициенты Ті в формулах приведены в прил. IX). Как показано в [28, 49], такой план весьма эффективен, на­ пример, при изучении действия индивидуальных химиче­ ских добавок на бетой при разных Ц/В (или других технологических факторов). Особенность плана заклю­ чается в том, что один фактор варьируется на трех уров­ нях, а другой на пяти — это очень удобно в технологи­ ческих экспериментах.

Ротатабельными (или почти ротатабелы-іыми [63]) являются трехуровневые планы на кубе, предложенные Боксом и Бенкиным. Особенностью этих планов являет­ ся то, что во всех строках отличны от нуля только два (при /( = 3-f-5), три (при К= 6 , 7, 9) или четыре (при /(=10, 11, 12, 16) фактора. Эта особенность планов весь­ ма важна для решения рецептурно-технологических за­ дач, так как позволяет стабилизировать ряд Хі в много­ факторной ситуации в течение некоторой группы опытов. Для /С=Зч-5 матрицы планирования приведены в табл. IV. 10 (с разрезкой на блоки они даны в [63]). Расчет­ ные формулы подобны (ІѴ.7І) —(ІѴ.77), кроме расчета

Ь0 по (ІѴ.78).

153