книги из ГПНТБ / Мастеров, В. А. Практика статистического планирования эксперимента в технологии биметаллов
.pdf8
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
известно, что суммарная относительная погрешность б оценки 0Эравна примерно 7%, а результаты параллель ных опытов подчиняются [78] нормальному распределе нию. Это позволяет оцепить точность взятых с рис. 29 данных в виде дисперсии воспроизводимости:
|
|
|
|
r |
Vс |
W=8_ |
|
|
|
|
||||
|
|
|
|
б |
|
S |
|
уи |
|
|
|
|
||
|
|
S2M |
= |
|
|
и=1 |
|
|
|
|
|
|
||
|
|
|
^v; а |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
где 6 = 0,07; с = 3; |
N = 8 ; |
8 _ |
|
|
|
|
/v,a—коэффицн- |
|||||||
Ег/„= 17,53; |
||||||||||||||
ент Стьюдента '. |
|
|
1 |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|||
В рассматриваемом примере число степеней свободы |
||||||||||||||
v= c— 1 = 2 |
и |
для |
уровня |
значимости |
0,05 |
величина |
||||||||
£=4,30. Тогда s{i/} = 0,0617, |
a s2{ y } = 38 •10-4. |
|
|
|||||||||||
Теперь можно проверить значимость |
коэффициентов |
|||||||||||||
модели: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
s |Ш |
|
. |
|
|
|
|
|||
|
|
|
|
[ - |
|
Г ) |
|
v; |
а |
|
|
|
|
|
Здесь дисперсия s2{6} ошибки определения коэффи |
||||||||||||||
циента одинакова для всех коэффициентов модели: |
|
|||||||||||||
|
|
S2 {</} |
38,2-10 ~ 1 |
|
1,59-10_ \ |
|
|
|||||||
|
|
N-c |
|
8-3 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
s{b} = 1,26-10-2. |
|
|
|
|
|||||||
Значение коэффициента |
Стьюдента |
£v,a = |
2,12 |
при |
||||||||||
N (с—1) = 8 |
■2 = |
16 степенях свободы и уровне значимо |
||||||||||||
сти 0,05; соответственно |
^знач^О)0267. |
Поскольку- |
b 12, |
|||||||||||
£?2з и Ь123 меньше 0,0267, они приравниваются нулю. |
|
|||||||||||||
Далее рассматривается уравнение-модель |
|
|
||||||||||||
у = 2,190 + |
0,364 |
+ |
0,264 Х2— 0,334 Х 3 + |
0,0762 Хг Х3. |
||||||||||
Если уравнение |
(41) |
верно, |
|
то вклад |
слагаемого |
|||||||||
0,0762 Х\ХЪстатистически должен быть незначим. |
|
|||||||||||||
Проверка адекватности модели |
|
|
|
|
|
|
|
|
||||||
у = |
2,190 + |
0,384 |
+ |
0,264 Х2 - |
0,334 X , |
(46) |
1 Если действительно ставится эксперимент 23 с с параллельны ми измерениями в 8 точках факторного пространства, то проверка воспроизводимости проводится, как описано в гл. III.
101
осуществляется |
с помощью дисперсии неадекватности |
и /•'-отношения, |
как описано в гл. III: |
s9: = —i— 382КГ'1= 95,5. КГ4
ад 8 — 4
счислом степеней свободы
v= N — d = 8 — 4 = 4.
О
F = |
95,5-10~4 = 2,50. |
s2{y) |
38,2.10—' |
Табличное значение FViy„]a— FA- 16; 0,os=3,01. Посколь ку оно больше экспериментального, гипотеза об адек ватности уравнения (1 ) опытным данным не отвергает ся с надежностью 95%.
В заключение определим параметры п, а, К и аст ин дивидуального материала. Из (1), (45), (46) следует:
|
|
Ах, |
2,32 |
|
„ п . |
|
|
|
|
п = —- = —— = 6,04, |
|
|
|||||
|
|
Ьх |
0,384 |
|
|
|
||
|
_ |
nb2 _ |
6,04-0,264 = 1,38, |
|
||||
|
|
Ах, |
1,16 |
|
|
|
||
— пЬз |
6,04-0,334 |
= |
n |
_ i |
, |
|||
к = |
|
|
100 |
|
0,0202 град |
|||
Ах. |
|
|
|
|
|
|
||
In0 СТ— Ь0 |
Д |
А- 0 1 |
Д |
А'02- |
*03 — |
|||
|
|
|
|
|
|
|
Да*я |
|
= 2t190 — 0^384(1 62) — 0i£64 (— 1 69)------М З?(1273) = |
||||||||
2,32 |
v |
|
1,33 |
v |
|
|
100 |
v |
= |
6,507; |
crCT= 630 кгс/мм2. |
|
|||||
Для иллюстрации применения модели (46) построе |
||||||||
ны графики as= a s(e) для е = 5 |
с- 1 и температур 900, |
|||||||
1000, 1100° С, рассчитанные |
по |
(46) |
линии |
показаны |
пунктиром. Следует отметить удовлетворительное сов падение расчетных и экспериментальных кривых. Одна ко для построения графиков с помощью планирования эксперимента потребовалось всего 8-3, т. е. 24 опыта вместо 108 при традиционной методике экспериментиро вания.
Описанный метод можно применить и для обработ ки литературных данных.
102
СПОСОБ ПОСТРОЕНИЯ ПОЛИНОМИАЛЬНОЙ МОДЕЛИ ПРОЦЕССА ПО ЛИТЕРАТУРНЫМ ДАННЫМ [7]
В ряде случаев несложно построить модель техно логического процесса по литературным данным в виде таблиц или графиков. Так, если выполнены следующие условия:
1 ) факторы Х\, х2...... хп варьируются на числе уров ней gi, g 2, ..., g n соответственно через равные интерва лы Ахй
2 ) в наличии полный перебор условий опытов, т. е. каждый уровень фактора ху сочетается со всеми уровня
ми остальных факторов (всего выполнено |
g iX g ^ X -'X |
||||||||
Xgn = N опытов); |
|
|
|
|
|
|
|||
3) для всех N опытов измерено значение параметра |
|||||||||
оптимизации уи, |
|
|
|
|
Т а б л и ц а 28 |
||||
|
|
|
|
|
|
|
|
||
|
|
|
Матрица планирования |
|
|
|
|||
Номер |
Х„ |
х, |
X, |
X? |
х,х. |
Х1 |
От |
||
опыта |
клик |
||||||||
и |
|
|
|
|
|
|
|
|
V % |
i |
+ 1 |
+ 1 |
+ 1 |
+ 1 |
+ 1 |
+ 1 |
41 |
||
2 |
+ 1 |
+ 1 |
+ 0 ,6 |
+ 1 |
+ 0 ,6 |
+ 0 ,3 6 |
42 |
||
3 |
+ 1 |
+ 1 |
+ 0 ,2 |
+1 |
+ 0 ,2 |
+ 0 ,0 4 |
44 |
||
4 |
+ 1 |
+ 1 |
—0,2 |
+ 1 |
—0,2 |
+ 0 ,0 4 |
56 |
||
5 |
+ |
1 |
- И |
—0,6 |
+ 1 |
- 0 , 6 |
+ 0 ,3 6 |
52 |
|
6 |
+ 1 |
- и |
— 1 |
+ 1 |
— 1 |
+ Г |
54 |
||
7 |
+ 1 |
0 |
+ 1 |
0 |
|
0 |
+ 1 |
61 |
|
8 |
+ 1 |
0 |
+ 0 ,6 |
0 |
|
0 |
+ 0 ,3 6 |
76 |
|
9 |
+ 1 |
0 |
+ 0 ,2 |
0 |
|
0 |
+ 0 ,0 4 |
67 |
|
10 |
+1 |
0 |
—0,2 |
0 |
|
0 |
+ 0 ,0 4 |
93 |
|
11 |
+ |
1 |
0 |
—0,6 |
0 |
|
0 |
+ 0 ,3 6 |
109 |
12 |
+ 1 |
0 |
— 1 |
0 |
|
0 |
+ 1 |
77 |
|
13 |
+ 1 |
— 1 |
+ 1 |
+ 1 |
—1 |
+ 1 |
63 |
||
14 |
+ 1 |
— 1 |
+ 0 ,6 |
+ 1 |
- 0 , 6 |
+ 0 ,3 6 |
69 |
||
15 |
+ 1 |
— 1 |
+ 0 ,2 |
+ 1 |
—0,2 |
+ 0 ,0 4 |
99 |
||
16 |
+ |
1 |
— 1 |
—0,2 |
+ 1 |
+ 0 ,2 |
+ 0 ,0 4 |
138 |
|
17 |
+ |
1 |
— 1 |
—0,6 |
+ 1 |
+ |
1 |
+ 0 ,3 6 |
201 |
18 |
-Ы |
— 1 |
+ 1 |
+ 1 |
+ |
1 |
. + 1 |
221 |
|
N= 18 {07} = |
{17} = |
{27} = |
{117} = |
{127} = |
{227} = |
|
|||
|
= 1563 |
= —502 = |
—307,4 |
= 1080 |
= |
223,6 |
= 734,52 |
|
ЮЗ
то без затрат на новый эксперимент можно построить модель процесса в виде полинома второй степени:
У = Ь0 + T ,bi x i + |
II bu |
X , + |
t bn X ? |
(47) |
|
i = i |
t = i |
|
; = i |
|
|
|
»'</ |
|
|
|
|
или неполного квадратичного полинома |
(при g — 2 ). |
|
|||
Решение задачи рассмотрим |
на примере из |
работы |
|||
[8]. Исследовали на gi = |
6 уровням влияние |
фактора |
х { и на g 2— 2 уровням фактора х2 на параметр у (отно сительная прочность, %) строительного материала. По сле перехода к безразмерным значениям факторов по формуле
v |
ду — 0,5 (max .\y-|-min *,-) |
|
|
Л; — ------------------- :---------- |
' |
||
|
0 ,5 (шах лу — min лу) |
||
была получена матрица эксперимента (табл. |
28). |
||
Коэффициенты модели (47) определяются по фор |
|||
мулам: |
|
|
|
b0 = N - ^ 0{OY} + N -1 £ г [у (аГ }; |
|
||
b{ = N~' срт-1 {/К}; |
|
||
b ^ N |
- ' |
{ijY } |
|
Ьи = Х - 1Ци {ИУ} + М -] ^ {О У }, |
|
||
|
{OK} = |
£ Уи, |
(48) |
|
|
«=i |
т = s х -jj,-
и=i
{//K}= £ xinx!uyn-
{ay}
«=1 |
; |
где /— число определяемых данным планом квадратич ных эффектов.
104
Остальные необходимые параметры в случае двух факторных экспериментов считываются из табл. 28, а при большем числе факторов определяются с помощью уравнений (48), причем:
Число |
|
|
|
уровней |
|
■47 |
47т |
фактора |
|
|
|
2 |
1 |
— 3 |
9/2 |
3 |
2/3 |
||
4 |
5/9 |
—45/16 |
81/16 |
5 |
1/2 |
—20/7 |
40/7 |
6 |
7/15 |
—375/128 |
5625/896 |
7 |
4/9 |
—3 |
27/4 |
Если в литературном источнике сообщается ошибка опыта s{t/} (или ее можно вычислить по результатам параллельных опытов), представляется возможным определить дисперсии коэффициентов модели и оценить адекватность модели:
s2 (&0) = |
ЛГ- 1 -ф0.52 [у\, |
s2 {b.\ = |
N - ' v f s2 {у}, |
s2 \ЬЦ) = |
N~'cp- 1 s2 {г/}, |
s2 (bu\= |
N ~ % t s2 {у}, |
cov {bQ, bu} = N ~ 1$ { s‘l {y}.
Оценка значимости коэффициентов модели произво дится по /-критерию Стьюдента:
, _ |
\bj\ |
1 |
S{*,} |
с числом степеней свободы |
v — N (c—1 ), где с — число |
повторений каждого опыта. Табличные значения t счи тываются из приложения II.
Проверка адекватности модели состоит в сравнении с табличным значением Е-отношения дисперсии неадек
ватности s2fl к дисперсии |
воспроизводимости s2{«/}: |
р ^ |
2 |
San |
^ S ^ y} ’
где дисперсия неадекватности
105
—
о
П л а н g i X f i 2
2 X 2
2 X 3
2 X 4
2 X 5
2 X 6
2 X 7
3 X 3
3 X 4
3 X 5
3 X 6
3 X 7
4 X 4
4 X 5
4 X 6
4 X 7
5 X 5
5 X 6
5 X 7
6 X 6
6 X 7
7 X 7
N Ч 0
V *
V 2
4l / l 2 8
17 / 70
3 0 3 / l 5 3 0
ч»
б / о
7 3 / l 9 2
3 1 / 106
-5-6 9 / 23 0 4
1 3 / б З
3 3 / 128
4 4 7 /2240 603/3072 187/l3«
2 7 /l76
3101/г0880 79/ 13Б
239 / 2304
7421/ 16128 U / 147
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а 29 |
|
Значения параметров для двухфакторной модели |
|
|
|
||||||||||
w |
1 Ф , 1 |
" |
‘ ф г 1 |
i V |
|
1 1| )j |
w |
4 2 |
w - 4 n |
W _ 1 Ф о 2 |
W 1 Ф ; 1 ф 2 1 |
|||
1 / 4 |
V 4 |
|
_ |
|
|
__ |
|
_ |
__ |
v « |
||||
V o |
V |
* |
' |
----- |
|
- |
V |
, |
— |
3 / 4 |
V 4 |
|||
V s |
9 / 40 |
|
— |
|
|
- |
4 5 / l 2 8 |
— |
8 1 / i 2 8 |
9 / 4 0 |
||||
|
|
|
|
|||||||||||
1 / 10 |
V o |
|
— |
|
|
- 2 / 7 |
— |
4 7 |
|
V o |
||||
* / 1 2 |
V |
2 8 |
, |
----- |
|
- |
1 2 ? / з 1 2 . |
■— |
1 8 7 5 / з 0 8 9 |
V 28 |
||||
V l 4 |
9 / в о |
|
— |
|
|
----- 3 / l 4 |
— |
27 / 50 |
9 / о о |
|||||
|
|
|
|
|||||||||||
V o |
V o |
- v |
|
3 |
- v |
3 |
V , |
Чг |
v « |
|||||
V s |
3 / 20 |
- V |
|
4 |
|
V / 04 |
3 / e |
2 7 / 0 4 |
V 40 |
|||||
1 / l 0 |
2 / ю |
- |
V |
o |
|
-----4 / 2 1 |
3 / l 0 |
8 / |
21 |
V o |
||||
V |
12 |
? |
/ « |
- |
V |
o |
|
- |
1 2 5 / |
708 |
V4 |
1 8 7 6 / |
0380 |
6 / 28 |
V14 |
3 / 28 |
- Ч 7 |
—V» |
3/l4 |
9/28 |
9 / о о |
||||||||
9 / s o |
9 / з о |
|
4 V |
25 0 |
-«/200 |
81/2oo |
8 1 / г оо |
81 / 400 |
||||||
9 / |
100 |
1 / i o |
— |
V04 |
-----V7 |
8 1 / з 20 |
2/7 |
9 / о о |
||||||
3 / |
40 • |
3 / о 0 |
— |
1Э/ |
128 |
----- / 1024 |
27 / 128 |
1876/7108 |
9 // 280 |
|||||
9/l40 |
9/112 |
|
4 - / 448 |
|
V |
28 |
8V448 |
27 / 112 |
8 1 / о о о |
|||||
°/25 |
V |
26 |
— |
4 /35 |
- |
4 / з о |
8 / 35 |
8 / з о |
4/2S |
|||||
V |
l 5 |
х / и |
- |
2 / |
2i |
- |
25/25 o |
4/г1 |
3 7 8 / 1792 |
Ч7 |
||||
2 / 35 |
9 / 140 |
— |
V |
|
49 |
— 3/35 |
8/49 |
27 / 140 |
9 /70 |
|||||
V |
84 |
V84 |
1 9 V |
|
1036 |
— |
125 / ю з е |
02 6 W |
8 2 8 / 3084 |
26/l90 |
||||
S/08 |
3 / о о |
— 12 ° / 1792 |
|
V |
14 |
18,5/l2544 |
9 / о о |
4 S / 392 |
||||||
9 / |
100 |
° |
/ 190 |
|
3 / |
49 |
|
3 / |
49 |
27 / 190 |
27/ 190 |
8V784 |
с числом степеней свободы v\ = N—d [d — число членов полинома (47)], а дисперсия воспроизводимости
s,fs}=s1{M + SMs! M +
|
|
1=1 |
|
|
+ |
+ г Е * ? с о » М » | |
|
|
<</ |
|
|
с числом |
степеней свободы V2 = N ( c — 1 ). |
||
Если |
полученное |
значение A ^ l, |
модель заведомо |
адекватна, в остальных случаях сравниваем F с таблич |
|||
ным значением (см. |
приложение I). |
При 77< / гтабл ги |
потеза об адекватности не отвергается с достоверностью (1—а ). Здесь а — уровень значимости критерия F, обычно а = 0,05 .
Для приведенного примера модель имеет вид:
у'— 79,67 — 41,83Xj — 36,60 Х 2 + 9,50 + 39,93 X tX2 + 1,79 Х|.
Глава V
ТЕОРИЯ И ПРАКТИКА ИССЛЕДОВАНИЯ ПРОЧНОСТИ БИМЕТАЛЛИЧЕСКИХ СОЕДИНЕНИЙ
В данной главе на примере биметалла сталь+алюминий рассмотрены: влияние отжига на сопротивление отрыву слоев а и ударную вязкость; статистическое пла нирование эксперимента для построения функции рас пределения а; вероятностная оценка минимальной проч ности и коэффициента запаса прочности биметалличес ких соединений; физическая модель падения прочности биметалла после нагрева и статистическое планирование эксперимента для ее проверки.
Листовой биметалл Х18Н10Т+АД1+АМг6 получен прокаткой нагретого до 450° С пакета с общим обжати ем 35—40% [46, 47]. Режим термической обработки го
107
тового листа на заводе-нзготовнтеле биметалла — отжиг при 480° С, 6 ч [48]. В готовом виде биметалл представ ляет собой трехслойный материал толщиной 10 мм с со отношением толщин слоев стали и сплава АМгб при мерно 1 : 1 и с толщиной алюминиевого подслоя марки АД1 0,35—0,40 мм.
При освоении промышленного производства и внед рении биметаллических изделий на основе композиции сталь+алюминий встретились трудности: нестабиль ность структуры и свойств при нагреве биметалла, уве личение разброса механических свойств после ' повтор ных нагревов.
Сходные проблемы возникают при освоении других композиций разнородных металлов и сплавов. Поэтому биметаллические и сварные соединения разнородных ме таллов стали объектом исследований в ряде отечествен ных и зарубежных организаций: ИЭС им. Патона. Гипроцветметобработка, ВИЛС, УИЧМ, NASA (США) идр.
ЛИТЕРАТУРНЫЕ ДАННЫЕ
Прочность сцепления оценивается значением а, вели чину которого большинство авторов связывают прежде
|
|
всего |
с |
наличием |
интерме- |
|||||
|
|
таллидов, |
образование |
и |
||||||
|
|
рост которых стимулируется |
||||||||
|
|
нагревом биметалла при его |
||||||||
|
|
изготовлении, сварке с ос |
||||||||
|
|
тальной |
конструкцией |
или |
||||||
|
|
при |
эксплуатации |
изделий |
||||||
|
|
из него. |
|
|
|
|
|
|||
|
|
На рис. 30 представлены |
||||||||
|
|
результаты Р. М. Сизовой, |
||||||||
|
|
В. К. Короля и Н. Д. Лу- |
||||||||
|
|
кашкина |
[45], |
из |
которых |
|||||
|
|
следует, |
что |
после отжига |
||||||
|
|
вначале о увеличивается, а |
||||||||
|
|
затем |
|
происходит |
падение |
|||||
|
|
прочности, которое исследо |
||||||||
|
|
ватели |
связывают |
с низкой |
||||||
Рис. 30. Влияние продолжительно |
прочностью |
и |
хрупкостью |
|||||||
интерметаллических |
соеди |
|||||||||
сти отжига (2, 6, |
12 ч) на сопротив |
|||||||||
ление отрыву |
слоев биметалла |
нений. |
Наиболее обстоятель |
|||||||
Х18Н10Т+АД1+АМг6 [45] |
|
|
|
|
|
|
|
|
но
ное исследование выполнено в ИЭС нм. Патона [49, 50], полученные результаты иллюстрируются рис. 31, 32.
В. Р. Рябов, Л. Г. Гвинчевская и А. В. Лозовская [50] отмечают, что прочность сцепления слоев биметалла оп-
Рис. 31. Влияние температуры отжига на сопротивление отрыву слоев биметалла Х18Н10Т+АД1+АМг6 после 10 (а) и 20 (б) мни отжига [49]
Рис. 32. Влияние длительности отжига при 520' (а) и 550 (б) °С на со противление отрыву слоев и толщину интерметаллндного слоя биметал
ла Х18НЮТ+АД1+АМг6 [50]
ределяется физико-химическим состоянием переходного слоя, прежде всего наличием интерметаллидов. Установ ленные авторами работы [50] (путем просмотра шли фов) температурные границы появления интерметалли дов показаны на рис. 31 вертикальной линией. В то же
109