Файл: Мастеров, В. А. Практика статистического планирования эксперимента в технологии биметаллов.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 образ цов, относительная погрешность оценки отдельного из мерения составляет 6—8 % [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