Файл: Баженов, Ю. М. Перспективы применения математических методов в технологии сборного железобетона.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 29.10.2024
Просмотров: 78
Скачиваний: 0
Матрица планирования (п кодовом обозначении) н результаты определения р (по R»Зг) приведены в табл. 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