Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Мастеров, В. А. Практика статистического планирования эксперимента в технологии биметаллов

.pdf
Скачиваний:
4
Добавлен:
20.10.2023
Размер:
7.86 Mб
Скачать

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

известно, что суммарная относительная погрешность б оценки 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 (c1 ), где с — число

повторений каждого опыта. Табличные значения 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

Соседние файлы в папке книги из ГПНТБ