Файл: Мастеров, В. А. Практика статистического планирования эксперимента в технологии биметаллов.pdf

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

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

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

Добавлен: 21.10.2024

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

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

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

Xi и параметрами модели Q,-. Одновременно возрастает вероятность одновременного существования нескольких моделей для одного процесса и одного г/*.

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

Итак, предположим, что такие модели заданы и вхо­ дящие в них параметры, факторы и функции описывают то, что инженер-технолог интуитивно считает «механиз­ мом» исследуемого процесса.

Уточним постановку задачи. Пусть сформулирована модель процесса в виде функции отклика:

У = У ( Х

1, * 2,

Ql t Q2,...,Q t,...)= y (X , Q),

где а-!, х2,

..., А',-,... — факторы;

Qi, Q2, ,..., Q,-,...— параметры модели.

Возможны три случая:

 

1.

Функция у(Х,

Q)

известна и притом единств

ная.

Требуется определить

пли уточнить неизвестные

параметры:

 

 

Qx

Q2

Q v

2 . Предложены две или более функций отклика:

У х ( Х , Q )

У(Х) =

. У а ( Х , 0 ) .

В общем случае набор параметров у них может быть различным. Требуется определить,' какая из функций является «истинной», и найти неизвестные параметры.

Для рассматриваемых случаев развиты эффективные методы статистического планирования эксперимента [3, 76]. Познакомимся с ним на примере из обработки металлов давлением.

96

ИССЛЕДОВАНИЕ СОПРОТИВЛЕНИЯ МЕТАЛЛА ГОРЯЧЕЙ ПЛАСТИЧЕСКОЙ ДЕФОРМАЦИИ [77]

В литературе предложено несколько уравнений связи сопротивления деформации os с температурой, степенью и скоростью деформации

 

 

ф; =

0's (Т, е, е),

(41)

например, уравнение (1)

из работы Ю. Н. Работнова [2]:

 

п

п

а

/

Tsn^\

 

os = a cre

s-exp (—КТ),

где

Т — температура, К;

 

е— интенсивность конечной (логарифми­ ческой) деформации;

е— интенсивность скорости деформации

(производная от е по времени), с-1; п, а, Ост (кгс/мм2), К (град) - 1 — парамет­ ры материала, которые находятся из опытов и характеризуют индивидуаль­

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

as— сопротивление деформации (интенсив­ ность напряжений) кгс/мм2.

Рис. 29. Сопротивление горячей деформации стали Ст.З [78]

7— 1193

97


Типичное представление опытных данных показано на рис. 29. Каждая точка на кривой — среднее арифме-. тическое нескольких параллельных измерений os, раз­ личающихся погрешностями эксперимента. Разброс опытных значений as затрудняет расчет параметров п, а, сгСт, К и проверку пригодности уравнений типа (41), например уравнения (1). Значения констант и пригод­ ность уравнения можно оценить лишь с некоторой веро­ ятностью (надежностью). Она приближается к 100% при возрастании числа параллельных опытов и точности аппаратуры.

На практике «на одну точку» расходуют 3—5 образ­ цов, относительная погрешность оценки отдельного из­ мерения составляет 68 % [78].

Функцию отклика (1) разложим в кратный ряд Тей­ лора и для технологических расчетов ограничимся первыми членами разложения в ряд, например, непол­ ным квадратичным полиномом

у ^ а0 + пу х,_ -|- а2 х2 -[- а3 х3 + а12лу х2 +

 

cii2x±x2 -\- a22x2xs а±22х^х2х2

(42)

или даже линейными членами разложения

 

у ^ а0 + оулу + а2х2 + а3 х3.

(43)

В случае уравнения (1) после логарифмирования и подобных преобразований получим

у = In су = In 0ГСТ -f — In е + —

In 8 ---- — Т.

(44)

 

11

п

 

п

 

Обозначим:

 

 

 

 

 

ду = In е, х3 — Т,

а2 = а/п

]

 

х2 = In е, а0 =

In аст,

а3 =

К/п

[

(45)

оу

л-1.

 

 

I

 

Тогда выражение (44) примет вид

у = In оу = а0 + а1х1 + а2х2 + а3х3,

совпадающий с линейной частью разложения уравнения

(42).

Факторный эксперимент типа 23 (см. гл. III) позво­ ляет статистически достоверно оценить коэффициенты уравнения (42) и (43) и адекватность уравнения (42)

98


или его линейной части (43),

или суммы любых членов

уравнения (42), опытным данным.

деформа­

Например, изучается [78]

сопротивление

ции стали СтЗ

при 9 0 0 ^ 0

^ 1100°С, 0,051 ^ е ^ О ,698

(т. е. 0,05^

относительное обжатие

^ 0,50),

0,5<е<!50 с-1. Для удобства

вычислений вводятся ко­

дированные значения переменных х*:

 

v

х( — 0,5 (max xi + minx,-)

 

 

0,5 (maxxi — minх{)

 

Эксперимент ставится в 8 точках (23) факторного прост­ ранства Х\, Х2, Х3, расположенных согласно планирова­ нию 23. В каждой точке проводится с > 1 параллельных опытов. Условия опытов указаны в табл. 26, матрица пла­ нирования п результаты опытов — в табл. 27. Планирова­ ние 23 позволяет строить модель вида

у = Ь0 Ь1Х1 + Ь2 Х2 + Ьа Х3 + Ь12 Х± Х2 +

 

+ Ь13 Хг Х3 +

Ь23 Х2 Х3 -|- Ь123 Хх Х2 Х3.

 

 

 

 

 

 

 

 

Таблица 26

 

 

Условия опытов на пластометре

 

 

 

 

Натуральные значения

 

 

Кодированные зна­

 

Интер­

 

чения переменных

Обозна­

переменных для уровней

 

уровней:

 

 

 

 

валы варь­

Код

 

 

 

чение

верх­

основного

 

ирования

верх­

ос­

ниж­

 

нижнего

Axt

 

 

 

нов­

 

него

*01

 

 

 

него

ного

него

лц= 1пе

+ 3 ,9 4

+ 1,62

— 0,70

2,32

 

+ i

0

— 1

х„ — 1пв

- 0 , 3 6

— 1,69

—3,02

1,33

х ,

+ i

0

— 1

* з = Т

1373

1273

1173

100

Х3

+ i

0

— 1

Коэффициенты модели рассчитывают, как описано выше:

bn =

Y u .. =

=

2,190,

Ь°

N ^ i y“

 

 

 

U= 1

 

 

 

8

3,07

 

 

 

 

&i = :t

I 1 ^ X i“ =

^8 = 0 ’384,

 

U=1

 

 

 

8

2,11

 

 

 

 

=

ы=1

T l =

0’264-

 

 

 

7*

99



8

63 =

T S

* 3" = ^

= ~

° ’334’

 

и=

1

 

 

 

8

 

 

 

ьа =

т

XiH А'2» = °~т =

° ’0137-

 

и=1

 

 

 

8

 

 

 

&13 = 7Г S

X l“ Хз“ = ^

=

° ’0762’

«= 1

8

Ью =

Y £

*■« *з« = ° ^

= 0,0162,

 

ы=1

 

 

 

8

 

 

= - у

2

*з„ =

= - 0,00375.

 

« = 1

 

 

Видно, что коэффициенты b0, b \, b2, bz значительно больше остальных, однако этого недостаточно для дока­ зательства адекватности ( 1 ) опытным данным.

Неясна точность использованных опытных данных. Из практики исследований на кулачковом пластометре

Таблица 27

Матрица планирования, результаты опытов на пластометре и расчета по модели

 

Кодиро--

 

 

 

 

 

 

 

 

ванные

Данные опыта

Данные

Отклонения расчетных

Номер

условия

 

 

данных от опытных

 

опыта

 

 

расчета

 

 

 

опыта

 

 

 

 

 

Л

 

 

 

и

 

 

 

 

Уи =

«и

Л

 

 

 

X,

х 2

X,

(5s)„

\»и - « и 1’ х

 

 

 

= In (°s)«

 

|"и-"н|

 

 

 

 

кгс/мм*

 

 

X 10<

1

 

 

 

7,3

2,00

1,90

0,10

 

100

2

+

11,4

2,46

2,40

0,06

 

36

3

+

+

13,0

2,58

2,64

0,06

 

36

4

+

____

22,0

3,11

3,17

0,06

 

36

 

 

 

 

 

 

 

 

 

5

 

 

+

3,1

1,14

1,21

0,07

 

47

6

+

+

5,3

1,68

1,73

0,05

 

25

7

+

____

+

7,6

2,04

1,98

0,06

 

36

8

+

+

13,4

2,62

2,54

0,08

 

64

N = 8

 

 

 

Сумма =- 17,63

 

Сумма =

380

100