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

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

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

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

Добавлен: 29.10.2024

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

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

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

Номера

1

2

3

4

 

5

 

G

7

точек

 

 

 

—1

+ і

1/2

1/2

- 1

/

2

- 1 /2

0

дга

0

0

і Т / 2

- Т Т /2

V T

/2

~ Ѵ г 12

0

Т а б л и ц а

IV.IO.

Матрицы

планов

Бокса — Бенкина

 

 

Номер опыта

 

 

 

Фактор

 

 

 

 

*1

хг

х3

*1

 

 

 

 

 

1-4

- ы

+

1

 

0

0

 

0

5-8

+ 1

 

0

-1 1

0

 

0

9-12

0

+

1

+

1

0

 

0

13-16

0

 

0

+

1

±1

 

0

17-20

+ 1

 

0

 

0

± 1

 

0

21-24

0

± 1

 

0

± 1

 

0

25-28

0

± 1

 

0

0

± 1

29-32

0

 

0

 

0

+ 1

±1

33-36

0

 

0

± 1

0

і

1

39-40

± 1

 

0

 

0

0

±

1

 

0

 

0

 

0

0

 

0

П р и м е ч а н и е .

Нулевую

строку

повторяют три раза

для К = 3 н- 4 и '"есть

раз для К=5.

 

 

 

 

 

 

 

 

154


bo='T1( O Y ) - T i %(UY) = y0,

(IV.78)

i=i

 

в которой i/й — среднее по результатам опытов в центре плана.

Хорошими статистическими характеристиками обла­ дают [57] композиционные, близкие к Д-оп'тималы-іым

планы на кубе, минимизирующие обобщенную диспер-

л

сию s2 {Y}. Наиболее приемлемы для решения задач бетоноведения и технологии близкие к Д-оптимальиым планы на кубе при Д=4-^5 типа В4, В5 и На5, обладаю­ щие симметричным расположением точек в факторном пространстве и дающие независимые оценки Ьі и bij коэффициентов модели. Эти планы в кодовой записи приведены в табл. IV.11. Расчетные формулы подобны

(IV.70) — (IV.77).

Значения Т{ приведены в прил. IX.

Т а б л и ц а IV.11.

Близкие к Д-оптимальным планы на кубе

План

 

 

 

Матрица

 

 

 

 

Общее число

 

 

 

 

 

 

 

точек

 

 

 

 

 

 

 

 

 

 

В4

ПФЭ-24

плюс восемь точек: (+1,

0,

0, 0),

24

 

(0,±1,

0,

0),

(0,

0,±1,

0),

(0,

0,

0,±1)

 

в5

ПФЭ 25

плюс десять точек:

(± 1, 0, 0,

0, 0),

42

 

(0, ± 1, 0, 0,

0),

(0, 0,

± 1,

0,

0),

(0, 0,

 

 

0, ± 1 ,

0),

(0, 0, 0, 0±1)

 

 

 

 

 

Нй5

Полуреплика 25

(см. прил. VIII)

плюс де­

27

 

сять периферийных точек по плану Bs,

 

 

плюс одна центральная точка

 

 

 

 

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

Действительно, при изучении влияния тепловлажно­ стной обработки при атмосферном давлении на механи­ ческие свойства бетонов температура среды как фактор не может быть более 100° С (это дает ограничение в ги­ перпространстве в виде плоскости), а в технологических

15.5


задачах часто приходится исследовать совместное влия­ ние факторов Хі именно в окрестностях этой плоскости.

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

Для решения основного вопроса об адекватности опи­ сания полиномом экспериментальных данных необходи­ мо рассчитать /"'-отношение (IV.18) и сравнить его с F.табл (обычно а = 5%). Оценка дисперсии s23v, характе­

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

плана (you) от их средних значений (уаи) при числе сте­ пеней свободы центральных опытов '»о—1 :

4 = t (£/„„— ,У: (« о - !)■

<!Ѵ-7Э>

U„=l

 

Оценка дисперсии неадекватности spIA для полинома

с числом значимых эффектов L рассчитывается по соот­ ношению

s h a ~ S V , '•[н л ~ [ 5 5 ОСт

S [УОн

Уои) ] •

 

и„=1

 

: ( N - L - n 0)t

(IV.80)

где остаточная сумма квадратов S S 0Cт определяется или

л

прямым расчетом по Уи и уи (IV.16), или по коэффи­ циентам регрессии:

550СТ= £

у і - 0 (0 Y) + V b i (LY)+

 

«=1

і=і

 

+ £ b l[(üY) + Z b l1(ijY)\.

(IV.8 Г)

і= 1

ІФІ

 

Следует отметить, что такой подход к проверке адек­ ватности возможен только в том случае, если измерения в - каждом из N опытов не дублируются. В противном случае следует проверять адекватность с учетомпарал­ лельных измерений (см. п. ІІІ.5 и [33]).

156


Для определения значимости коэффициентов регрес­ сии рассчитывается bKp — tTiS3 (IV.25). Если часть. Ьц оказывается статистически незначимой, то их не следует удалять из модели, так как все Ьц связаны между собой

и с Ь0 (или нужно полностью пересчитать модель по ме­

тоду

последовательного регрессионного анализа —

п. IV.4). От технологической интерпретации ßü = 0 сле­

дует

воздержаться.

Пример IV.11. При изучении на физических моделях [69] свойств ячеистых бетонов оценивалось влияние на R„av четырех показателей,

характеризующих мембраны ячеистой структуры (табл. IV. 12). В ка­ честве параметра выхода Y рассматривалось (а) относительное из­ менение прочности на изгиб любого u-образца по сравнению с нан-

лучшнм результатом (Дизг) макс, достигнутым во всей совокупности опытов.

Т а б л и ц а IV. 12. Варьирование характеристик мембран пористой структуры

Хі

A',

X,

X*

Характеристики

 

—1

0

- и

Капиллярная

пористость

мембран

10

32,5

55

в %

 

 

 

75

65

55

Плотность мембран в % плотности

смеси

н

дыр в

шт. на

10

30

50

Число трещин

100 пор. образующихся в мембра­

 

 

 

нах в процессе формирования по­

 

 

 

ристой структуры

 

 

250

450

650

Число мнкротрещнн в шт. на I л

материала, образующихся в мем­

 

 

 

бранах во время

тепловлажиост-

 

 

 

ной обработки

 

 

 

 

 

Выбор_ такого критерия оптимизации для данной

задачи был

весьма удобен, так как

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

собой степень влияния Х і

на разные свойства материала,

но и срав­

нить степень влияния различных элементов структуры Х і

и Z; (кро­

ме системы «Мембраны»

с факторами Х і моделировалась система

типа «Поры» с факторами Z ,):

 

100%;

(а)

У и ~

І^нзг.макс

 

 

—о I {R u}

S2 {Ru

(б)

[Уи) = Уи

Г>±

Rl

Амакс.

 

Для планирования эксперимента использован план типа В4. В табл. IV. 13 показана' матрица плана в ко­ довой записи и результаты эксперимента. По формулам

(IV.70) — (IV,75) рассчитана модель:

157


 

Т а б л и ц а

1Ѵ.13. План

В4 и

результаты

опытов

 

 

 

 

 

Результат

 

 

 

 

 

 

Результат

 

Код

 

 

 

2

 

 

Код

 

 

2

 

опыта

 

»н

 

опыта

 

 

»„

 

 

 

 

SIIW

 

 

 

 

 

 

suw

 

1

abed

 

26

2,39

 

13

 

cd

 

76

20,45

2

abc

 

48

8,16

 

14

 

с

 

 

90

28,67

3

abd

 

30

3,19

 

15

 

d

 

 

85

25,58

4

ab

 

52

9,57

 

16

 

а )

 

100

35,4

 

5

aed

 

35

4,34

 

17

a (0)

 

50

8,85

6

ас

 

57

11,5

 

18

-

a

(0)

 

86

26,18

7

ad

 

46

7,49

 

19

 

6(0)

 

53

9,94

8

а

 

62

13,61

 

20

- b

(0)

 

68

16,37

9

bed

 

69

16,85

 

21

 

e (0)

 

56

11,1

 

10

be

 

82

23,8

 

22

—c (0)

 

66

15,42

11

bd

 

73

18,86

 

23

d ( 0)

 

49

8,5

 

12

b

 

88

27,41

 

24

- d

(0)

 

76

20,45

 

Л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У — 63,28 — 19,22*! — 5,39*2 — З,5*3 — 9,22*., +

 

+ 6,24*^ — 2,79*о — 2,29*з — 0,79*^ — 0,31*, *2 +

0,31*, *3 —

 

 

1,56*,*! +

1,06*2*3 0,31*2*4 0, 19*3*4,

 

(в)

По 5зш= 3,95

рассчитаны Ь,ір

при

/=1,677

для а=10%

(г)

и получена

(д)

окончательная

модель

(незначимые

Ьц взяты

в скобки), адекватная по /^-критерию

при

а = 5%:

 

 

 

(^о)кр == 3,17;

(&;)кр =

1,56;

(і>т/)кр =

1 >66;

(ö,-,-)Kp =

4,12.

(г)

Для сравнительного анализа н технологической интерпретации

результатов использована

модель

 

 

 

 

 

 

 

 

У =

63 —

19*, —

5*2— 4*з —

9*4+ 6*^ —

( 3*2+

2*з +

*4) .

(д)

В частности, отмечалось:

а) среди четырех анализируемых факторов структуры пористого материала капиллярная пористость влияет на предел прочности при изгибе весьма нелинейно; при малой пористости ее изменение ска­ зывается на R„3t сильнее, чем при большой; влияние всех остальных

факторов на 7?НЗг почти прямо пропорционально (в исследуемом диа­

пазоне

значимы лишь линейные эффекты модели)

их величине;

б)

увеличение величины всех кодированных факторов * і , что со­

ответствует увеличению капиллярной пористости— Хі

и росту числа

дефектов в мембранах А'3 и Х4, но одновременно снижению плотно­ сти перегородок Х 2 (см. табл. Ѵ.12), приводит к уменьшению проч­

ности на

изгиб

(все

й ,-< 0);

в)

по

абсолютной

величине коэффициентов ) 6» | можно оценить

степень

влияния

каждого из структурных факторов; поскольку

158