|
11.1. Úvod
Měření rychlosti šíření ultrazvuku v matrici kovových
materiálů, jako charakteristika jejich vnitřního uspořádání,
je známá a používaná metoda již řadu let, a to i v
oblasti kontroly litinových odlitků. Zprvu jen pro odlitky z
litiny s lupínkovým grafitem (LLG) a dnes, s postupným nárůstem
používáním litiny s kuličkovým grafitem (LKG), i pro
kontrolu vlastností této matrice.
Měření rychlosti podélných ultrazvukových vln ve vztahu
k mechanickým vlastnostem (jakými jsou například hodnoty
meze pevnosti, tvrdosti, tažnosti, E-modulu a další) se osvědčilo
zejména u litiny s lupínkovým grafitem a byly nalezeny
odpovídající korelační vztahy.
U skupiny LKG se však ukazuje, že nelze vždy věrohodně s
vysokým koeficientem korelace aplikování zjištěné funkční
závislosti na druhý soubor dat. To je zapříčiněno tím,
že litiny vykazují větší nehomogenitu a celkovou proměnlivost
své struktury než například oceli
Teorie měření rychlosti ultrazvuku předpokládá, že
plocha, která je ve styku s měřící sondou by měla být
geometricky rovná, tj. broušená, protilehlé plochy rovnoběžné,
jejich plocha trojnásobkem průměru sondy a tloušťka
prozvučovaného materiálu kolem 20 mm. Podíváme-li se však
na odlitky, zjistíme, že prozvučované stěny neodpovídají
často jak geometrickým charakteristikám, tak především
požadavkům na kvalitu povrchu (úpravou by se ztratila jedna
z výhod této metody měření - možnost 100% kontroly v sériové
výrobě) [98].
Přesto, na základě experimentálních prácí, například
v [99] byly nalezeny odpovídající vztahy s relativně vysokými
korelačními koeficienty i pro LKG. Zde měření probíhala
na válečcích vypracovaných z přilitých tyčí. Pokusili
jsme se proto využít této metody měření i u dalších
odlitků z litiny s kuličkovým grafitem a to i přesto, že
vzhledem k nízkým korelačním koeficientům u tohoto materiálu
se nedoporučuje používat tuto metodu pro zjišťování tažnosti
a tvrdosti, naopak je vhodná pro kontrolu pevnosti v tahu.
Jedná se rovněž o metodu, jak hodnotit strukturní
charakteristiky matrice - např. podíl zrnitého grafitu,
neboť vlivem nevyhovujícího tvaru grafitu může dojít ke
značnému zhoršení mechanických vlastností.
Měření ultrazvukovou metodou tak představuje rychlou nepřímou
metodu zkoušení, která dovoluje stanovit na základě
rychlosti průchodu ultrazvuku stěnou odlitku její
charakteristiky.
11.2. Experimentální měření
Za účelem ověření, zda-li je možné zodpovědně převzít
korelační vztahy dříve publikované v literatuře a zda-li
je vůbec možné takové vztahy, které by obecně popisovaly
chování rychlosti ultrazvuku v závislosti na materiálových
vlastnostech matrice odlitků z LKG, nalézt, jsme na našem
pracovišti VUT FSI - Odbor slévárenství provedli, ve
spolupráci se slévárnou Motor JIKOV, České Budějovice,
slévárna litiny a.s., následující experiment, při kterém
bylo prováděno měření na několika typech odlitků.
Přestože se ve všech případech jednalo o stejnou litinu s
kuličkovým grafitem, odlitky byly odlévané v jedné slévárně
(jednou technologii) a do jednotné bentonitové formovací směsi
(vyjma lití do kovové kokily), pro každý typ vznikly
samostatné statistické soubory. Díky tomu můžeme sledovat
změny parametrů, které jsou ovlivněny z velké části jen
konstrukci odlitku, resp. návrhem slévárenské formy. Cílem
je rovněž i ověřit homogennost statistických souborů a
testovat hypotézu o normálním rozdělení.
Pro měření byl použit přístroj ECHOMETR 1060, sonda HB2
- 10 MHz o průměru 10 mm, akustická vazba olej. Úprava
povrchu pro měření u tyčí - zarovnány čela jemným
soustružením, u válce, třmenů, držáků a pák byly
protilehlé plochy broušeny. Odlitky přilitých tyčí, třmenů,
držáků a pák byly odlévány do bentonitové směsi,
odlitek válce do litinové kokily.
Ve slévárně Motor JIKOV, slévárna litiny a.s. bylo odlito
těchto pět typů vzorků z taveb o předepsaném chemickém
složení pro litinu ČSN 42 23 05.:
1. 154 přilitých tyčí určených pro tahovou zkoušku
[99]
2. 33 komerčně litých odlitků třmenů, kde jsme
prováděli měření na místech označených jako čelo
3. Třmen - ramena (dva vzorky na každém odlitku, po
třech měřeních na jednom rameni)
4. Třmeny v provozu
5. Experimentální odlitek válce litý do kokily - měření
bylo prováděno ve čtyřech řezech a na každém z nich na
36 místech rovnoměrné rozmístěných po ploše
6. Držák
7. Páka
Jako veličinu, kterou jsme hodnotili u všech odlitků, byla
absolutní rychlost šíření ultrazvuku. Dále byla měřena
tvrdost HB a to vždy na stejných místech jak měření cL.
Na odlitku válce (jediný experimentálně navržený
odlitek, neboť ostatní odlitky jsou komerčně lité) a třmenů
bylo změřeno teplotní pole a provedena analýza
mikrostruktury a simulace teplotního pole. Na odlitcích přilitých
tyčí byla vyhodnocena jak rychlost ultrazvuku, tak i pevnost
Rm, tažnost A5, tvrdost HB, chemické složení a struktura.
Na zbylých odlitcích byla vyhodnocena jen rychlost
ultrazvuku, popřípadě tvrdost HB.
Vztah mezi rychlosti šíření ultrazvuku v závislosti na
tvrdosti materiálu HB jsme zvolili nejen pro jednoduchost a
rychlost jejího určení, ale především proto, že nám
podává obrázek o chování materiálu, kde byla měřena i
rychlost cL sondou a nebylo pro ní třeba vzorek dále
mechanicky opracovávat, plochy byl jen obrobeny na strojních
bruskách.
Cílem takto rozsáhlé analýzy bylo přiřadit mechanické
vlastnosti jednotlivým místům a tak zjistit jejích
citlivost na tvar odlitku, teplotní pole a strukturu. A také
možnost věrohodně zjišťovat výsledné mechanické
vlastnosti pomocí měření průchodu ultrazvuku. Za tím účelem
byla u každého typu odlitku provedena základní popisná
statistika.
11.3. Přilité tyče
V analýze bylo navázáno na práci [99], kde bylo prováděno
měření na přilitých tyčích určených pro zkoušku
tahem - litina s kuličkovým grafitem ČSN 42 2305. Před měřením
rychlosti ultrazvuku, tvrdosti HB, určení meze pevnosti v
tahu Rm, tažnosti a také metalografických parametrů byly
tyče obrobeny.
Nejprve jsme ověřili, zda nedocházelo při měření
vlastností k systematické chybě s postupujícím měřením.
Za tím účelem jsme zhotovili: graf 1 - Rychlost šíření
ultrazvuku v závislosti na pořadí měření, graf 2 - Mez
pevnosti v tahu vs. pořadí měření, graf 3 - Tvrdost HB vs.
pořadí měření a graf 4 - Tažnost A5 vs. pořadí měření.
Hodnotili jsme i chemické složení vs. pořadí měření.
Grafy 1-4 jsou v příloze na CD-ROM k této kapitole.
Výsledek analýzy
Ze všech sledovaných průběhů je patrná mírná závislost
na postupujícím měření. Ve všech případech dochází,
při aproximaci hodnot polynomem druhého stupně, k poklesu
(nárůstu) s postupujícím měřením, s extrémem kolem 80
měření a opětovnému nárůstu (poklesu). Například u
rychlosti ultrazvuku činí tento rozdíl 50 ms-1, u pevnosti
cca 50 MPa.
Následně byla proto provedena i analýza chemického složení.
V grafu 5 v příloze na CD-ROM je to obsah uhlíku, v grafu 6
- obsah křemíku, graf 7 - obsah manganu, graf 8 - obsah
fosforu až graf 9 - obsah síry.
Z grafů 5 - 9 je vidět, že obsahy prvků uhlíku a křemíku
"kopírují" průběh mechanických vlastností
(extrém je kolem 80 měření), obsahy dalších prvků již
mají trend buď trvale (při aproximaci polynomem druhého
stupně) klesající - obsah síry, nebo trvale rostoucí -
obsah fosforu a manganu.
V grafu 10, resp. 11 je zobrazen procentuální podíl složky
perlitu, resp. feritu. I tyto grafy ukazují na klesající
(stoupající) trend s extrémem kolem 60 měření a s následným
poklesem (nárůstem).
Polynom druhého stupně byl zvolen s ohledem na zachycení
trendu celého souboru dat a zároveň omezení vlivu místních
výchylek. Lze soudit, že tento trend není zapříčiněn
subjektivními vlivy, neboť měření provádělo několik
osob na několika pracovištích, ale vlivem technologie výroby,
neboť vzorky byly lity průběžně s výrobou.
První série vzorků "přilitých tyčí" pochází
z práce [99], kde bylo provedeno základní měření
statistického souboru. V tabulce 36 je výsledek základní
popisné statistiky proměnných rychlost ultrazvuku cL [ms-1],
pevnosti v tahu Rm [MPa], tvrdosti HB a tažnost A5 [%].
Tab.
36 Popisná statistika pro odlitek přilitých tyčí
|
Proměnná
|
Rychlost
cL
|
Pevnost
Rm
|
Tvrdost
|
A5
|
|
|
[ms-1]
|
[MPa]
|
[HB]
|
[%]
|
|
Velikost
souboru
|
154
|
154
|
154
|
154
|
|
Průměrná
hodnota
|
5569,6
|
562.9
|
188.5
|
9.2
|
|
Medián
|
5570
|
560.5
|
187
|
9
|
|
Modus
|
5610
|
538
|
187
|
7
|
|
Geometrický
průměr
|
5569,3
|
560.3
|
188.2
|
8.8
|
|
Rozptyl
|
3220,1
|
3091.9
|
135.0
|
8.0
|
|
Směrodatná
odchylka
|
56,8
|
55.6
|
11,6
|
2.8
|
|
Standardní
chyba
|
4,573
|
4.481
|
0.936
|
0.228
|
|
Minimum
|
5410
|
417
|
170
|
2
|
|
Maximum
|
5700
|
758
|
224
|
20
|
|
Rozpětí
|
290
|
341
|
54
|
18
|
|
Dolní
kvartil
|
5540
|
525
|
180
|
7,3
|
|
Horní
kvartil
|
5610
|
591
|
195
|
11
|
|
Šikmost
|
-0,366
|
0.523
|
0.820
|
0.527
|
|
Standardizovaná
šikmost
|
-1,855
|
2.649
|
4.154
|
2.672
|
|
Špičatost
|
0,384
|
1.471
|
0.220
|
1.159
|
|
Standardizovaná
špičatost
|
0,973
|
3.727
|
0.559
|
2.936
|
|
Koef.
rozptylu
|
1,019
|
9.877
|
6.165
|
30.651
|
|
Suma
|
857710
|
86691
|
29029
|
1423
|
V tabulce
37 je základní statistika pro odlitek přilitých tyčí –
obsah perlitu a feritu.
Tab. 37 Popisná statistika pro odlitek přilitých tyčí
|
Proměnná
|
Perlit
|
Ferit
|
|
|
[%]
|
[%]
|
|
Velikost
souboru
|
154
|
154
|
|
Průměrná
hodnota
|
28,8
|
70,9
|
|
Medián
|
25
|
75
|
|
Modus
|
20
|
80
|
|
Geometrický
průměr
|
-
|
69,26
|
|
Rozptyl
|
218,6
|
218,9
|
|
Směrodatná
odchylka
|
14,78
|
14,80
|
|
Standardní
chyba
|
1,19
|
1,19
|
|
Minimum
|
0
|
30
|
|
Maximum
|
70
|
100
|
|
Rozpětí
|
70
|
70
|
|
Dolní
kvartil
|
20
|
60
|
|
Horní
kvartil
|
40
|
80
|
|
Šikmost
|
0,34
|
-0,32
|
|
Standardizovaná
šikmost
|
1,28
|
-1,60
|
|
Špičatost
|
-0,78
|
-0,76
|
|
Standardizovaná
špičatost
|
-1,94
|
-1,93
|
|
Koef.
rozptylu
|
51,28
|
20,86
|
|
Suma
|
4440
|
10924
|
V tab.
38 je statistika pro odlitky přilitých tyčí – chemické
složení
Tab.
38 Popisná statistika pro odlitek přilitých tyčí
|
Proměnná
|
C
|
Si
|
Mn
|
P
|
S
|
Cu
|
|
|
[%]
|
[%]
|
[%]
|
[%]
|
[%]
|
[%]
|
|
Velikost
souboru
|
150
|
150
|
150
|
150
|
150
|
99
|
|
Průměrná
hodnota
|
3,57
|
2,70
|
0,23
|
0,20
|
0.02
|
0,23
|
|
Medián
|
3,57
|
2,70
|
0,22
|
0,20
|
0.02
|
0,26
|
|
Modus
|
3,60
|
2,60
|
0,20
|
0,20
|
0.01
|
0,27
|
|
Geometrický
průměr
|
3,57
|
2,70
|
0,23
|
0,20
|
0.02
|
0,21
|
|
Rozptyl
|
0,00658
|
0,0145
|
0,001919
|
2,41
E-6
|
4.08
E-6
|
0,00707
|
|
Směrodatná
odchylka
|
0,08117
|
0,1204
|
0.043807
|
1,55
E-3
|
2.02
E-3
|
0,08407
|
|
Standardní
chyba
|
0,00662
|
0,0098
|
0,003577
|
0,00013
|
0,000170
|
0,00845
|
|
Minimum
|
3,35
|
2,14
|
0,16
|
0,013
|
0,011
|
0,08
|
|
Maximum
|
3,88
|
2,8
|
0,35
|
0,022
|
0,021
|
0,45
|
|
Rozpětí
|
0,53
|
0,76
|
0,19
|
0,009
|
0,01
|
0,37
|
|
Dolní
kvartil
|
3,51
|
2,6
|
0,19
|
0,015
|
0,014
|
0,15
|
|
Horní
kvartil
|
3,63
|
2,81
|
0,26
|
0,017
|
0,017
|
0,27
|
|
Šikmost
|
0,6181
|
-0,5508
|
0,6041
|
0,3997
|
0,11026
|
0,2365
|
|
Standardizovaná
šikmost
|
3,091
|
-2,7542
|
3,0206
|
1,998
|
0,55130
|
0,9607
|
|
Špičatost
|
0,6984
|
1,544
|
-0,3939
|
0,6770
|
-0,425
|
-0,1169
|
|
Standardizovaná
špičatost
|
1,74601
|
3,861
|
-0,9846
|
1,25
|
-1,064
|
-0,2373
|
|
Koef.
rozptylu
|
2,27
|
4.56
|
18,99
|
9,67
|
13,51
|
36,54
|
|
Suma
|
536,020
|
405.380
|
34,600
|
2.410
|
2.245
|
22.780
|
Analýza
grafů
V příloze k této kapitole jsou uvedeny grafy jak
popisující naměřené hodnoty chování matrice a jejího
chemického složení v závislosti na pořadí měření
(graf 1 - 11), tak i histogramy (graf 12-18), a korelační
grafy (graf 19-21). Pro prvotní náhled na naměřená data
slouží histogramy pro rychlé vizuální posouzení naměřených
dat, jejich homogenity a normality.
Hodnocení
histogramů četnosti
Graf 12 Histogram četnosti pro pevnost Rm [MPa]
Z grafu se dá soudit na možnou normalitu naměřených dat.
Graf 13 Histogram četnosti pro tvrdost HB
U tohoto grafu jsou patrné dva vrcholy četnosti, což by
ukazovalo na nehomogenní soubor dat
Graf 14 Histogram četnosti pro tažnost A5 [%]
I u tohoto grafu jsou patrné rovněž dva vrcholy, jeden
kolem 7% a druhý ve 12%.
Graf 15 Histogram četnosti pro absolutní rychlost ultrazvuku
UV [ms-1]
Lze soudit na možnost homogenních dat, počáteční naměřené
hodnoty mezi 5400 až 5500 ms-1 vypovídají o horších
vlastnostech odlitků
Graf 16 Histogram četnosti pro obsah uhlíku C [%]
Zde stojí za povšimnutí odlehlé body (3,3 % až 3,4 %) a
(3,75 % až 3,9%), kterých je však velmi malý počet a
ukazují na nevyrovnanost ve složení taveniny
Graf 17 Histogram četnosti pro obsah křemíku Si [%]
Pomineme-li výrazně odlehlý bod (2,1%), který je zapříčiněn
nejspíše chybou měření, je patrné, že se jedná o
nehomogenní soubor dat se dvěma extrémy.
Graf 18 Histogram četnosti pro obsah manganu Mn [%]
Co platí o histogramu v grafu 17, platí i zde. Dva extrémy
a kdybychom neznaly fyzikální podstatu, soudilo by se na dva
sloučené statistické soubory v jeden
Takřka u všech histogramů se ukazují dva vrcholy v grafu.
Nejpatrnější je to u grafu rozdělení obsahu křemíku a
manganu.
Hodnocení regresní závislosti
Graf 19 Regresní závislost mezi absolutní rychlosti
ultrazvuku cL [ms-1] a pevnosti v tahu Rm [MPa] spolu s pásy
spolehlivosti na hladině 95 a 99 %
Z tohoto grafu je patrné, že mezi těmito dvěma veličinami
lze získat velmi dobrý koeficient korelace. To dokazuje i
graf 48 v příloze, s koeficientem korelace je 0,7473 (Rm [MPa]
= 0,8471 cL [ms-1] - 4154,9).
Graf 20 Regresní závislost mezi absolutní rychlosti cL [ms-1]
a tvrdosti HB spolu s pásy spolehlivosti na hladině 95% a 99
%
Co platilo v předcházejícím bodě, nebude platit pro závislosti
mezi tvrdostí HB a rychlosti cL. Pásy spolehlivosti jsou značně
od sebe, a koeficient korelace je 0,22 (HB = 0,0965 UV [ms-1]
- 349,1).
Graf 21 Regresní závislost mezi absolutní rychlosti
ultrazvuku UV [ms-1] a tažnosti Á5 [%]
Pásy spolehlivosti jsou výrazně od sebe Ţ velmi nízký
koeficient korelace
Analýza hlavních komponent
Za účelem komplexní analýzy pro soubor dat z odlitků
přilitých tyčí jsme použili metodu hlavních komponent,
tab. 39. Analýza hlavních komponent patří mezi vícerozměrné
statistické metody užívané k redukci výchozího počtu
proměnných. Jestliže součet prvních (nejvyšších) k těchto
podílu je dostatečně blízký jedné (ze zkušenosti se požaduje
kolem 0,85 - 0,9), postačí nám brát v úvahu právě těchto
prvních k hlavních komponent. [100]
Tab. 39. Výsledky analýzy hlavních komponent
|
Komponenta
|
Procentuální
podíl [%]
|
Kumulativní
podíl [%]
|
|
Rm
[MPa]
|
48,023
|
48,023
|
|
Tvrdost
[HB]
|
16,775
|
64,798
|
|
A5
[%]
|
14,185
|
78,984
|
|
Obsah
C [%]
|
9,055
|
88,039
|
|
Obsah
Si
[%]
|
4,708
|
92,748
|
|
Obsah
Mn [%]
|
3,389
|
96,137
|
|
Obsah
P
[%]
|
2,300
|
98,438
|
|
Obsah
S
[%]
|
1,319
|
99,757
|
|
Obsah
Cu [%]
|
0,242
|
100,000
|
Z
tabulky 39 je zřejmé, že pro vysvětlení chování
rychlosti ultrazvuku u tohoto typu odlitku je nutné brát do
úvahy první čtyři komponenty tj. vedle meze pevnosti Rm a
tvrdosti HB také hodnotu tažnosti A5 a obsah uhlíku. [101,
102]
Vyhodnocení statistické analýzy přilitých tyčí
Vzhledem k dostatečnému rozsahu měření u souboru dat
"Přilité tyče" jsme, jako u jediného souboru v této
práci, přistoupili ke komplexnější statistice za použití
profesionálního programového vybavení - QC Expert [101].
Ostatní měření byla vyhodnocována pomocí software
Stagraphics.
Veškeré tabulky zjištěných statistických charakteristik
jsou v příloze k této kapitole. Kvůli odlišnostem v
statistických charakteristikách obou programů, jsou v příloze
stručně popsány významy jednotlivých charakteristik
program QC Expert. Popis je převzat z literatury [101].
Cílem analýzy bylo zjistit, zda-li naměřené soubory dat
pocházejí z normálního rozdělení a zda-li jsou data
homogenní.
Hodnocení odchylky od nuly je testem, kdy zkoumáme to,
zda-li je hodnota šikmosti statisticky významně odlišná
od 0. V takovém případě nelze data považovat za symetrická,
neboť normální a každé symetrické rozdělení má šikmost
nulovou. Spolehlivější je však test normality [102].
Podobným testem je hodnocení rozdílu od tří, kdy normální
rozdělení má špičatost 3. Je-li hodnota špičatosti
statisticky významně odlišná od 3, lze předpokládat, že
data neodpovídají normálnímu rozdělení. Spolehlivější
je však opět test normality [102]. U statistiky dělané
pomocí programu Statgraphics je tato hodnota porovnávána k
nule.
Výstupem programu QC Expert je slovní výrok o tom, zda-li
jsou tyto odchylky významné či nevýznamné.
Z provedené analýzy v tab. 1 v příloze je vidět, že
soubor dat "Absolutní rychlosti šíření
ultrazvuku" je významně odlišný od nuly i od tří. A
proto, podle těchto testů, zkoumaný soubor dat nejspíš
nepochází z normálního rozdělení a není ani symetrický.
Ke stejnému závěru jsme dospěli i u souboru dat
"Pevnosti v tahu", a souboru měření
"tvrdosti".
Šikmost u dat rychlosti cL je -12,04, což je významně odlišné
od nuly. Hodnota špičatosti je dokonce 148,23, na rozdíl od
hodnoty tří.
Šikmost u souboru dat Pevnost Rm je -2,78, což je odlišné
od nuly. Hodnota špičatosti je dokonce 26,08, na rozdíl od
hodnoty tří.
Šikmost u souboru Tvrdost HB je -5,92, což je odlišné od
nuly. Hodnota špičatosti je dokonce 112, na rozdíl od
hodnoty tří.
Všechny tří soubory dat nepocházejí, podle tohoto testu,
z normálního rozdělení.
U souboru dat Tažnost je šikmost 0,33, což je hodnota blízká
nule a software QC Expert ne zamítl hypotézu o normálním
rozdělení. Avšak podle testu špičatosti je její hodnota
4,28, což je hodnota vzdálená již od tří a normalita
souboru je zamítnuta.
Po provedeném testu normality - tab. 3 v příloze, program
QC Expert u souborů dat Abs. rychlost ultrazvuku, Pevnost Rm
a Tvrdost HB zamítl hypotézu o normálním rozdělení.
Hypotéza byla opět přijata jen u souboru Tažnost.
Dalším provedeným testem byl test na vybočující body -
tab. 4 v příloze, resp. na homogenitu dat. Ve všech čtyřech
souborech byla homogenita zamítnuta. U souboru Tažnost byla
vypočtena hodnota "Spodní mez" pod nulou, což je
z fyzikálního modelu nesprávné.
Autokorelační funkce slouží jako nástroj k posouzení závislosti
mezi jednotlivými hodnotami řady. Tato procedura odhaduje
korelační koeficienty pro jednotlivé dvojice stejně vzdálených
hodnot v řadě, přičemž její řád (tj. vzdálenost
dvojic) je specifikován. V našem případě zařazujeme tuto
statistickou proceduru pro možnost odhalení vzájemně závislých
měření - například chyb při měření.
Mimo souboru dat Abs. rychlost ultrazvuku, se u dalších
souborů vyskytuje významná autokorelace, která se však s
postupujícím řádem autokorelace stává nevýznamnou. U
použití čtvrtého řádu je autokorelace u všech čtyř
souborů dat nevýznamná - tab. 5 v příloze.
Nyní budeme podrobovat analýze obsahy chemických prvků:
uhlíku C, křemíku Si, manganu Mn, fosforu P, síry S a mědi
Cu na vzorcích přilitých tyčí.
Z provedené analýzy šikmosti a špičatosti - tab. 6 v příloze,
je vidět, že soubory dat popisující chemické složení C,
Si, Mn, P a S jsou významně odlišné od nuly i od tří. A
proto, podle těchto testů, zkoumané soubory dat nepocházejí
z normálního rozdělení a nejsou ani symetrické.
Jen u souboru popisující obsah mědi, jsou oba sledované
parametry vyhodnoceny jako statisticky nevýznamně odlišné
od porovnávacích hodnot.
Šikmost u obsahu uhlíku je -10,85, což je statisticky významně
odlišné od nuly. Hodnota špičatosti je dokonce 128,39, na
rozdíl od hodnoty tří.
Šikmost u obsahu křemíku je -8,24, což je statisticky významně
odlišné od nuly. Hodnota špičatosti je dokonce 88,69, na
rozdíl od hodnoty tří.
Šikmost u obsahu manganu je 8,93, což je statisticky významně
odlišné od nuly. Hodnota špičatosti je dokonce 83,05, na
rozdíl od hodnoty tří.
Šikmost u obsahu fosforu je 6,78, což je statisticky významně
odlišné od nuly. Hodnota špičatosti je dokonce 47,57, na
rozdíl od hodnoty tří.
Šikmost u obsahu síry je -1,52, což je statisticky odlišné
od nuly. Hodnota špičatosti je 12,07, na rozdíl od hodnoty
tří.
Šikmost u obsahu mědi má hodnotu 0,10, což je statisticky
velmi blízko nule a hodnota špičatosti je 2,98. Na základě
těchto hodnot byl soubor popisující obsah mědi, vyhodnocen
jako soubor dat pocházející z normálního rozdělení.
Pro potvrzení závěru z předcházejícího testu slouží
test normality. Zde program QC Expert u souborů dat popisující
obsah uhlíku, křemíku, manganu, fosforu a síry zamítl
hypotézu o normálním rozdělení. Hypotéza byla opět přijata
jen pro soubor popisující obsah mědi - tab. 8 v příloze.
Dalším provedeným testem byl test na vybočující body,
resp. na homogenitu dat tab. 9 v příloze. A zde, u prvních
pěti souborů dat, byla homogenita zamítnuta. Nezamítnuta
byla jen u souboru dat Cu.
U souboru popisující obsah mědi byla vypočtena hodnota
"Spodní mez" pod nulou, což je z fyzikálního
modelu nesprávné.
Posouzení závislosti mezi jednotlivými hodnotami časové
řady je v tab. 10 v příloze. U souborů dat obsah C, Si a
Mn nebyla, ani při jednom řádu autokorelace, detekována
pravděpodobnost autokorelace. U souboru Cu byla zjištěna
autokorelace do třetího řádu. U řádu autokorelace 4, již
ani pro tento soubor, nebyla odhadnuta autokorelace. Avšak
pro soubor dat popisující obsah síry se vyskytuje
autokorelace i při řádu 4.
Nyní budeme analyzovat procentuální podíl jednotlivých
struktur matrice: obsah perlitu, feritu a cementitu na vzorcích
tyčí.
Z provedené analýzy šikmosti a špičatosti je vidět, že
soubory dat popisující procentuální podíl feritu a
cementitu jsou významně odlišné od nuly a od tří - tab.
11 v příloze. A proto, podle těchto testů, zkoumané
soubory dat nepocházejí z normálního rozdělení a nejsou
ani symetrické. Jen u souboru popisující obsah perlitu,
jsou oba sledované parametry vyhodnoceny jako statisticky nevýznamně
odlišné od porovnávacích hodnot Ţ je naděje, že
testovaný soubor pochází z normálního rozdělení.
Šikmost u obsahu feritu je -0,76, což je statisticky odlišné
od nuly. Hodnota špičatosti je 4,25, na rozdíl od hodnoty tří.
Šikmost u obsahu cementitu je 8,22, což je statisticky významně
odlišné od nuly. Hodnota špičatosti je dokonce 74,78, na
rozdíl od hodnoty tří.
Šikmost u obsahu perlitu je hodnota 0,31, což je statisticky
velmi blízko nule a hodnota špičatosti je 2,23. Na základě
těchto hodnot byl soubor popisující obsah perlitu v
matrici, vyhodnocen jako soubor dat pocházející z normálního
rozdělení.
Pro potvrzení závěru z předcházejícího testů slouží
test normality. Zde program QC Expert u souborů dat popisující
obsah feritu a cementitu opět zamítl hypotézu o normálním
rozdělení. Hypotéza byla přijata jen pro soubor popisující
obsah perlitu - tab. 13 v příloze.
Dalším provedeným testem byl test na vybočující body,
resp. na homogenitu dat. A zde, u souborů dat feritu a
cementitu, byla homogenita zamítnuta. Nezamítnuta byla u
souboru dat perlitu.
U souboru perlitu byla vypočtena hodnota "Spodní
mez" pod nulou, což je z fyzikálního modelu nesprávné
- tab.14 v příloze.
U souboru dat popisující obsah cementit nebyla ani při
jednom řádu autokorelace detekována pravděpodobnost
autokorelace. U souborů periltu a feritu byla zjištěna
autokorelace do třetího řádu. U řádu autokorelace 4, již
ani pro tyto soubory, nebyla odhadnuta autokorelace - tab. 15
v příloze.
Závěr z analýzy přilitých tyčí
Pozornost, kterou jsme věnovali kontrole normálního
rozdělení a homogenitě naměřených souborů dat, byla
založena na snaze analyzovat jednotlivé procesy a jejich
ovlivnění dalšími faktory, které nejsou podchyceny. Výsledky
několikerých testů jsou hodnoceny jak numericky, tak slovně
(slovní výsledek je výsledek matematického algoritmu
programu QC Expert 2.). Pokud hodnoty signalizují Nenormální
rozdělení jsou označeny červeně,
je-li výsledek pro normální rozdělení pozitivní - označeno
zeleně - tab.40.
Testy robustních charakteristiky, přestože je nekomentují,
mohou rovněž mnohé napovědět o homogenitě dat. Jsou-li
hodnoty průměrů či mediánu blízko sebe i při odstranění
10 %, 20 % a 40 % hodnot, svědčí to o stabilitě naměřených
hodnot a jejich nevychýlenosti.
Pro kontrolu jakosti v průmyslu jsou podstatné i hodnoty
dolní a horní hranice, pod níž je možno data považovat
za vybočující.
Tab. 40 Testování hypotéz o normálním rozdělení a
homogenitě souborů
|
Proměnná
|
Normální
rozdělení
|
Homogenita
|
|
Abs.
rychlost ultrazvuku
|
Zamítnuto
|
Zamítnuta
|
|
Pevnost
Rm
|
Zamítnuto
|
Zamítnuta
|
|
Tvrdost
HB
|
Zamítnuto
|
Zamítnuta
|
|
Tažnost
A5
|
Přijato
|
Zamítnuta
|
|
Obsah
C
|
Zamítnuto
|
Zamítnuta
|
|
Obsah
Si
|
Zamítnuto
|
Zamítnuta
|
|
Obsah
Mn
|
Zamítnuto
|
Zamítnuta
|
|
Obsah
P
|
Zamítnuto
|
Zamítnuta
|
|
Obsah
S
|
Zamítnuto
|
Zamítnuta
|
|
Obsah
Cu
|
Přijato
|
Přijata
|
|
Zastoupení
perlitu
|
Zamítnuto
|
Zamítnuto
|
|
Zastoupení
feritu
|
Zamítnuto
|
Zamítnuta
|
|
Zastoupení
cementitu
|
Zamítnuto
|
Zamítnuta
|
Přestože
soubor dat z měření na přilitých tyčí pochází z
jednoho typu odlitku, litého ze stejné slitiny, ve stejné
slévárně Motor JIKOV ČB a.s., do stejné bentonitové
formy, bylo použito stejné technologie obrobení vzorků,
stejné metodiky měření veličin, vždy ve stejných
laboratořích, z testů je vidět, že většina proměnných
pochází z nehomogenních souborů a nejsou normálního rozdělení.
To svědčí o vlivu dalších nepodchycených parametrů,
které ovlivňují samotné měření a to i v laboratorních
podmínkách.
S tímto závěrem je nutné uvažovat při dalších analýzách
naměřených hodnot, tvorbě korelačních vztahů a analýze
jejich, většinou, nízkých korelací. Pro vizuální ověření
těchto závěrů, jsou v grafech 12 - 18 v příloze
histogramy četnosti naměřených hodnot.
11.4.
Odlitek Třmene
Druhým typem odlitku byly třmeny (výrobce - Motor JIKOV, slévárna
litiny a.s., tavby o předepsaném chemickém složení pro
litinu ČSN 42 23 04, bentonitová směs, podrobněji v
kapitole 10). Na těchto odlitcích jsme provedli měření cL
a HB na čelní ploše a to na osmi a na dvanácti místech -
obr. 48 a 49, dále na místech ramen - obr. 50. Tyto
experimenty byly prováděny v laboratořích VUT FSI, odbor
slévárenství. Měření rychlosti ultrazvuku cL bylo rovněž
prováděno přímo ve výrobním závodě - třetí soubor
dat. Vzhledem k tomu, že se jedná o rozdílné statistické
soubory, budou analyzovány samostatně.
Popisná statistika byla dělána softwarem Statgraphics. Výsledné
hodnoty analýzy jsou v tabulkách 41 - měření na osmi místech,
42 - měření na dvanácti místech a 43 - statistika při
zahrnutí obou typů měření. Histogramy a grafy hodnot jsou
v příloze k této kapitole.
Analýza dat
Nejprve jsme přistoupili k analýze dat prostřednictvím
histogramů hodnot.
Graf 22 v příloze. Histogram rychlosti ultrazvuku pro
odlitek třmene - měření na 12 místech
Soubor 184 měření se jeví jako homogenní data pocházející
z mírně vychýleného normálního rozdělení.
Graf 23. Histogram tvrdosti HB pro odlitek třmene - měření
na 12 místech
Co platilo pro měření rychlosti ultrazvuku nelze říci pro
tyto data, která jsou zjevně nehomogenní a nepocházejí z
normálního rozdělení.
Graf 24. Histogram tvrdosti HB pro odlitek třmene - měření
na 8 místech
Je vidět, že tyto data pocházejí nejspíše ze dvou rozdělení
sloučených v jeden soubor dat.
Graf 25. Histogram rychlosti ultrazvuku pro odlitek třmene -
měření na 8 místech
Tento graf hodnot rychlosti ultrazvuku kopíruje předcházející
graf dat získaný měřením tvrdosti. Z toho se dá soudit,
že se nejedná o náhodné vlivy měření, ale popis skutečnosti.
Budeme-li testovat normální rozdělení u souboru dat -
rychlost ultrazvuku - tab. 41, pak podle hodnoty šikmosti
tento soubor dat není symetrického rozdělení (šikmost =
0,309). Podle hodnoty špičatosti (0,001 " 0) soubor
pochází z normálního rozdělení.
U souboru tvrdostí HB měřených na osmi místech - tab. 42,
je hodnota šikmosti je 0,505 (nesymetrické) a hodnota špičatosti
je 0,25 Ţ což je hodnota signalizující .na normální rozdělení
Budeme-li testovat normální rozdělení u souborů hodnot
pro měření na dvanácti místech pak zjistíme, že šikmosti
u souboru rychlostí ultrazvuku na dvanácti místech je
-0,78. Podle špičatosti (1,710) se jeví jako soubor, který
nepochází z normálního rozdělení. Podobně se to ukazuje
i u souborů měření tvrdostí na dvanácti místech - graf
23. Hodnota šikmosti je -0,21 (nesymetrické) a hodnota špičatosti
je -0,91, což je významně odlišné od nuly, vůči které
testujeme na normální rozdělení.
Když sloučíme předcházející analyzovaná data v jeden
soubor - tab. 43, zjistíme, že hodnota šikmost je 0,25 pro
měření rychlost ultrazvuku a 0,08 pro měření tvrdost.
Podle tohoto typu testu je měření tvrdostí symetrické. U
hodnot špičatosti jsou odchylky od hodnoty platné pro normální
rozdělení větší, pro soubor rychlosti cL je to hodnota
-0,24 a pro hodnoty tvrdosti HB je to -0,44 Ţ data nepocházejí
z normálního rozdělení.
11.5. Třmen - ramena
Druhým místem měření na odlitcích třmenů byly
plochy označované jako ramena. Zde jsme prováděli měření
na třech místech označených 1, 2 a 3 dle obr. 45 (viz. také
obr. 43 a 44). Na stejných místech bylo prováděno jak měření
rychlosti ultrazvuku, tak i tvrdosti HB, abychom mohli tyto
dva statistické soubory vzájemně porovnávat. Celkem tak
bylo provedeno na 42x3 měření rychlosti cL a 29x3 měření
tvrdostí HB.
Ramena byly před měřením broušena na strojní brusce v dílnách
odboru slévárenství. Cílem měření je posoudit, zda-li
se i na tak malé ploše jakou jsou ramena (rozměr 20 x 30 x
8 mm) projeví vliv polohy na výsledný statistický soubor a
dále možnost přiřazení závislosti mezi měřením
rychlosti ultrazvuku a tvrdosti HB.
Sumární statistika pro odlitek třmene - ramen je v tab. 45
a 46 a v tab. 44, kde jsou seskupeny všechny tři soubory dat
pro jednotlivá místa měření v jeden soubor.
Analýza dat
Budeme-li testovat hypotézu o normálním rozdělení pro
data pocházející z odlitků třmene - ramena (tab. 44) pro
měření cL, nepochází podle hodnoty šikmosti = 1,07 a špičatosti
0,90 tento soubor dat z normálního rozdělení. U souborů
tvrdosti HB je hodnota šikmosti -0,81, což ukazuje na
nesymetrické rozdělení. Hodnota špičatosti je dokonce
6,17, což je výrazně mimo oblast pro kterou by byla přijatelná
hypotéza o normálním rozdělení.
Když budeme analyzovat samostatně jednotlivé místa na
ramenech (tab. 45), je u rychlosti šíření ultrazvuku
hodnota šikmosti, postupně od místa 1 do místa 3, rovna
-1,59 , 0,75 a 0,62. U hodnoty špičatosti je to 7,69 , 0,44
a 0,097 a hypotéza o normálním rozdělení by byla přijatelná
jen pro soubor z místa 3.
U analýzy souborů dat popisující tvrdost HB na jednotlivých
místech zjistíme, že hodnota šikmosti je pro místo 1 =
1,31, pro místo 2 = 0,61 a pro místo 3 = 0,66. Z toho lze
soudit, že data nejsou symetrické. Hodnota špičatosti je u
naměřeného souboru pro místo 1 = 1,94, pro místo
2 = -0,50 a pro místo 3 = -0,12. Jen místo 3 by mohlo, podle
tohoto testu, pocházet z normálního rozdělení, avšak
test na symetrii tento závěr nepodporuje.
Při porovnání výsledků popisné
statistiky pro sobory dat z jednotlivých míst na ramenech,
zjistíme, že pocházejí ze tří statistických souborů s
"posunutou" jak průměrnou hodnotou (místo 1 -
5556 ms-1, místo 2 - 5550 ms-1 a místo 3 - 5549 ms-1), tak i
mediánem (místo 1 - 5556 ms-1, místo 2 - 5548 ms-1 a místo
3-5545 ms-1). Průměrná hodnota ze všech třech míst je
5551,2 ms-1 a rozdíl mezi místy 1 a 3 je 7 ms-1, rozdíl v
mediánu je 11 ms-1 [110].
Přestože výsledky maximálních rozdílů mezi jednotlivými
místy nejsou velké, je patrný vliv měřeného místa na
vzorku. To bylo utvářeno teplotním polem při tuhnutí a
chladnutí odlitku, a statisticky významným způsobem ovlivňuje
strukturu matrice i pro tak omezenou oblast jako plocha vzorku
ramen (20 x 30 mm) při lití do bentonitu. Výsledky jsou
zachyceny v tab. 47.
Tab.
47 Výběr ze statistiky pro rychlosti ultrazvuku a tvrdosti
HB
|
Místo
|
1
|
2
|
3
|
Rozdíl
|
|
Průměr
HB
|
5555.33
|
5549.45
|
5548.64
|
6,69
|
|
Medián
HB
|
5554
|
5547.5
|
5544.5
|
9,50
|
|
Geometrický
průměr HB
|
5555.2
|
5549.4
|
5548.58
|
6,62
|
|
Průměr
cL
|
203.4
|
200.5
|
195.6
|
7,8
|
|
Medián
cL
|
195
|
198
|
191
|
7
|
|
Geometrický
průměr cL
|
202
|
200
|
195
|
7
|
11.6.
Třmeny v provozu
Třetí hodnocení odlitku třmene pomocí měření cL,
bylo prováděno přímo v závodě Motor JIKOV České Budějovice
a.s., slévárna litiny. Místo a podmínky měření odpovídaly
používané metodě a postupu, včetně výběru místa, jaký
se používá jako provozní výstupní kontrola. Podrobnosti
v kapitole 10.7. - obr. 50.
Analýza dat
Sumární statistka je v tab. 48. Hodnota šikmosti je
0,23 a hodnota špičatosti 1,057. Výsledky, i podle přílohy
ke kapitole 11 - graf 31 a 32, se jeví pro zamítnutí hypotézy
o normální rozdělení. Na rozdíl od ostatních měření,
kde probíhalo měření na více místech po ploše odlitku,
zde se měřilo vždy na jednom místě. To se projevilo ve
shodě hodnoty aritmetického a geometrického průměru, mediánu
a modusu = 5729 ms-1. Rovněž rozpětí minimální a maximální
hodnoty je jen 115 ms-1, zatímco u měření na ramenech je
rozpětí cL 247ms-1. Přesto podle provedeného testu tento
soubor dat pochází spíše z nenormálního rozdělení, což
je zapříčiněno neobrobennosti plochy před měřením. Zajímavé
je srovnání s obrobenou plochou na ramenech třmenů - zde v
místech 2 a 3 je rozptyl hodnot ještě nižší, a to jen
106, resp. 113 ms-1.
11.7. Obsah perlitu na třmenech
Posledním měřením na odlitku třmene, které bylo
statisticky vyhodnocováno, je obsah perlitu na čelech třmene.
V tab. 49 je souhrnná statistika
Z výsledků hodnoty šikmosti (0,12) se dá soudit, že data
pocházejí ze symetricky rozděleného souboru, hodnota špičatosti
(-0,80) však nepotvrzuje hypotézu o normálním rozdělení
dat. U měření obsahu perlitu se tak projevil vliv místa,
resp. teplotního pole, jak bylo prokázáno v kap. 10.
Hodnota kvartilového rozpětí je 15, což je takřka 50 % z
průměrné hodnoty, zatímco například u měření
rychlosti ultrazvuku dle tab. 48 je kvartilové rozpětí jen
0,5 % z průměrné hodnoty.
Graf průběhu obsahu perlitu po ploše vzorku je na obr. 35 v
příloze ke kap. 10.
11.8. Experimentální válec
Odlitek válce byl podroben nejkomplexnějšímu
vyhodnocování tvrdosti a rychlosti ultrazvuku ze všech
zkoumaných odlitků. Důvodem bylo nejen korektní změření
teplotního pole a to na pěti místech na jednom odlitku, ale
především "dostatek" vhodných ploch po obrobení,
které bylo možné použít pro měření jak tvrdosti, tak i
měření rychlosti ultrazvuku. Na obr. 73 je rozdělení válce
na jednotlivé vzorky. Červená šipka označuje směr válce
v kokile. Obrobení bylo provedeno v mechanické dílny VUT
FSI odbor slévárenství na čtyři válce. Plochy byly broušeny
na mechanické strojní brusce.
Každá plocha byla rozdělena na 37 míst a na nich bylo změřeno
jak tvrdost HB, tak i rychlost ultrazvuku.
V tab. 50 je statistika souboru dat tvrdosti pro jednotlivé
řezy na válci a v posledním sloupci výsledky po sloučení
všech hodnot tvrdosti ze všech řezů
V tab. 51 je statistika souboru dat hodnot rychlosti
ultrazvuku pro jednotlivé řezy na válci a v posledním
sloupci po sloučení všech hodnot ze všech řezů.
Analýza dat
Při porovnání výsledků z jednotlivých řezů dospějeme
k jednoznačnému závěru, že co řez to samostatný
statistický soubor. Hodnoty aritmetického, geometrického průměru
či mediánu nám charakterizují polohu křivky rozložení
hodnot vůči sobě a jsou závislé na poloze teplotního
pole odlitku. Tento závěr platí jak pro hodnoty tvrdosti,
tak i pro rychlosti ultrazvuku.
Testování normálního rozdělení ukazuje, že ani pro
jeden řez nelze potvrdit hypotézu o normálním rozdělení.
To je dáno tím, že teplotní pole formující výsledné
vlastnosti odlitku bylo proměnlivé jak po výšce kokily,
tak také od středu odlitku k jeho povrchu Ţ vlastnosti
odlitku jsou v jednotlivých řezech ovlivněny gradientem
teplotního pole a nelze, jak plyne z technické představy,
mluvit o jednom homogenním statistickém souboru dat. Naše závěry
potvrzuje pohled na příslušné histogramy 33 - 42 v příloze.
V tabulkách 52 a 53 jsou uvedeny hodnoty šikmosti a špičatosti
ze statistických testů.
Tab.
52 Testování hypotézy o normálním rozdělení pro soubor
hodnot tvrdosti
|
Typ
souboru
|
HB_A
|
HB_B
|
HB_C
|
HB_D
|
|
Šikmost
|
1,361
|
0,552
|
0,463
|
0,4667
|
|
Špičatost
|
2,301
|
-0,876
|
-1,216
|
-1,408
|
Tab.
53 Testování hypotézy o normálním rozdělení pro soubor
hodnot rychlosti ultrazvuku
|
Typ
souboru
|
UV_A
|
UV_B
|
UV_C
|
UV_D
|
|
Šikmost
|
0,654
|
0,374
|
0,609
|
0,097
|
|
Špičatost
|
-0,258
|
-1,254
|
-0,762
|
0,912
|
Rozložení
tvrdosti po ploše jednotlivých řezů 1 až 4 je na obrázcích
88 až 91 v příloze ke kap. 10. Při pohledu na obrázek 82
v příloze ke kap. 10 vidíme teplotní pole na válci. Je zřejmé,
že větší rychlost chladnutí je směrem od osy válce ke
stěně - styku odlitku s chladnou kokilou, což způsobuje nárůst
tvrdosti a také rychlosti cL.
Z obrázků 105 až 108 v příloze ke kap. 10 vidíme, že
nejrychleji chladnoucí plochou byl řez 4, následován řezem
3, 1 a 2. Když toto pořadí porovnáme s hodnotami tvrdosti
z tab. 50, pak nejvyšší průměrné tvrdosti - 329 HB je
dosaženo na vzorku 4 (plocha D). Následováno tvrdosti 322
HB na vzorku 3 (plocha C), tvrdosti 289 HB na vzorku 1 (plocha
A). Nejnižší průměrná tvrdost je na vzorku 2 (plocha B)
- 286 HB. Stejného pořadí získáme i pro geometrické průměry.
Když porovnání průměrné rychlosti ultrazvuku na vzorcích
A až D - řezy 1 až 4, zjistíme, že nejvyšší průměrné
hodnoty 5794ms-1 je dosaženo opět v řezu 4 (plocha D). Dále
však již pořadí neodpovídá, neboť na ploše A je
hodnota 5711, na ploše B 5710 a na ploše C 5669 ms-1. Příčinu
lze hledat v komplexním působení chemických prvků a
jejich nehomogenného zastoupení v jednotlivých místech.
K podobnému paradoxu dospějeme, když hodnotíme velikost
grafitu na ploše. V řezu 4, na ploše D je největší
grafit, zatímco na ploše B je nejmenší. Přitom, podle
teorie, by měl být v místě rychleji chladnoucím větší
počet grafitu o menší velikosti. Zde je však tomu naopak
(obr. 94 v příloze ke kap. 10). Možné vysvětlení je v
obsahu uhlíku v jednotlivých řezech - obr. 85 v příloze
ke kap. 10.
11.9. Odlitek držáku
Jako další odlitek pro experimentální měření
rychlosti ultrazvuku a tvrdosti dle HB byl vzat odlitek držáku,
jehož podobu s vyznačeným směrem výřezem vzorků vidíte
na obr. 100 a obr. 101, v kapitole 10.
Měření rychlosti ultrazvuku bylo prováděno na vyříznutých
a obrobených vzorcích buď na 8-mi nebo na 12-ti místech. Měření
tvrdosti dle HB bylo prováděno na stejných vzorcích a to
na 8-mi místech. Celkem bylo změřeno 176 měření
rychlosti šíření ultrazvuku a 160 měření tvrdosti.
Statistické vyhodnocení následuje v tab. 54.
Vyhodnocení
popisné statistiky
Testovat hypotézy o normálním rozdělení pro odlitek
držáku, podobně jako u odlitku experimentálního válce,
není dost dobře možné, neboť data nepocházejí z jednoho
místa, ale jedná se o sběr dat po délce vzorků, kde bylo
proměnlivé teplotní pole.
Naše závěry dokumentují i hodnoty šikmosti, resp. špičatosti.
U souborů pro tvrdost je hodnota šikmosti 1,2 a hodnota špičatosti
0,73, což jest odlišné od hodnot pro které by byla přijatelná
hypotéza o symetrickém normálním rozdělení. U hodnot
rychlosti ultrazvuku, by bylo možné dělat opatrnou domněnku
na symetrické rozdělení, neboť šikmost je blízká nule
(-0,1), ale špičatost svědčí o zamítnutí hypotézy o
normalitě (0,40).
Na grafu 43 44 v příloze k této kapitole jsou histogramy
rozdělení naměřených hodnot rychlosti ultrazvukových vln
a tvrdosti HB. I zde je patrné, že soubor dat pocházející
z měření cL je homogennější a bližší normálnímu
rozdělení, než-li soubor dat tvrdosti.
11.10. Páka
Na odlitcích páky - obr. 59 byla analyzována rychlost
ultrazvuku po dvou vzorcích na každém odlitku. Čela
experimentálních válečků byla broušeny na mechanické
strojní brusce pro zajištění hladkého povrchu a dosažení
lepší akustické vazby. Celkem bylo takto zpracováno na 17
odlitků. Výsledek popisné statistiky je v tab. 55.
Vyhodnocení popisné statistiky:
Jak plyne z výsledků popisné statistiky, je hodnota šikmosti
1,83, špičatost 3,23 a nelze potvrdit hypotézu o symetrickém
normálním rozdělení. Ke stejnému výsledků dospějeme i
při hodnocením histogramu - graf 45 v příloze.
11.11.Závěr
Oblast tyčí pro tahovou zkoušku se vyznačuje nejnižší
tvrdosti (od 170 HB) a nejnižší rychlosti ultrazvuku - od
5410 ms-1. Avšak rychlosti šíření ultrazvuku odpovídající
nejnižším hodnotám tvrdosti "zasahují" až do
oblasti hodnot rychlosti ultrazvuku, která by spíše odpovídala
kokilovému lití válce (max. 5710 ms-1) a tudíž i vysokým
hodnotám tvrdosti. Rozpětí hodnot cL je 290 ms-1. Pokud
budeme hodnotit kvartilové rozpětí, je 70 ms-1. To
dokumentuje graf 46 v příloze.
Na odlitku třmene jsme prováděli v laboratořích měření
na dvou oblastech - čelo a ramena a výsledkem jsou dva
statistické soubory. Oblasti Odlitek třmene - čelo je nejnižší
tvrdost od 174 HB a rychlost ultrazvuku je od 5529 do 5686 ms-1.
Zatímco oblasti "Odlitek třmene - ramena" odpovídá
vyššímu rozptylu rychlosti cL (od 5387 do 5634 ms-1) a vyšším
tvrdostem (až do hodnoty 285 HB). Je zřejmé, že se jedná
o samostatné statistické soubory, které byly převážně
utvářeny rozdílným teplotním polem.
Na obr. 18 a 19 v příloze ke kap. 10 je teplotní pole na čele
třmene. Ukazuje se, že utvářející podmínky se měnily
bod od bodu. To dokazuje i statistická analýza, která
nepotvrdila, že data pocházejí z normálního rozdělení.
Na obrázku 24 v příloze ke kap. 10 je zobrazen průběh průměrné
hodnoty tvrdosti po obvodu čela třmene a na obrázku 27 tamtéž
jsou zobrazeny průměrné hodnoty cL po obvodě čela třmene.
Potvrzení, že hodnoty tvrdosti a cL jsou závislé na poloze
na čele, resp. na rozdílném teplotním poli, dokazuje i
obr. 36 v příloze ke kap. 10, kde je zobrazen průběh
obsahu perlitu po obvodě po čele třmene a obr. 38, kde je
zobrazen průměrný počet grafitu po obvodě čela třmene.
Tento obrázek koresponduje s obr. 40, kde je zobrazen průměrný
počet pravidelného grafitu.
Určit souvislost mezi teplotním polem a výslednou
strukturou byl experiment na části ramen třmene. Bylo prokázáno,
že pro jednotlivá místa měření tvrdosti a rychlosti
ultrazvuku 1 až 3 vznikají samostatné statistické soubory
- tab. 47, ke kterým je možné přiřadit teplotní pole -
obr. 13 a 14 v příloze ke kap. 10. Teplotní pole ovlivňuje
samozřejmě i strukturu (obr. 74 v příloze ke kap 10, kde
je zobrazena simulace obsahu feritu na ramenech třmene), která
následně působí na tvrdost a rychlost průchodu
ultrazvukových vln.
Při porovnání výsledků popisné statistiky pro sobory dat
ze tří míst na ramenech zjistíme, že pocházejí ze tří
statistických souborů se vzájemně posunutou jak průměrnou
hodnotou (místo 1 - 5556 ms-1, místo 2 - 5550 ms-1 a místo
3 - 5549 ms-1), tak i mediánem (místo 1 - 5556 ms-1, místo
2 - 5548 ms-1 a místo 3 - 5545 ms-1). Přestože rozdíly
mezi jednotlivými místy nejsou velké, je patrný vliv měřeného
místa. Teplotní pole při tuhnutí a chladnutí odlitku tak
statisticky významným způsobem ovlivňuje strukturu matrice
i pro tak omezenou oblast jakou je plocha vzorku ramen při
lití do bentonitu. Z měření je rovněž patrné, že
rychlost cL popisující chování uvnitř stěny odlitku
vykazuje obecně menší rozptyl hodnot, než měření
tvrdosti HB. Například hodnoty tvrdosti pro místo 1 na
ramenech třmene jsou mezi 170 až 285 HB, což je rozdíl 115
HB Ţ procentuální podíl z průměrné hodnoty je přes 50
%! Rovněž v místě měření 2 je tento poměr přes 30 %,
podobně jako v místě 3.
U měření rychlosti ultrazvuku je poměr rozdílu maximální
a minimální hodnoty ku průměrné hodnotě mnohem menší -
pro místo 1 je rozdíl 247 ms-1, což je jen 4,5 %. V místě
2 je to 2 % a v místě 3 opět přes 2 %. K podobným závěrům
dospějeme i u dalších měření. Z toho vyplývá, že měření
rychlosti ultrazvuku vykazuje řádově menší rozptyl
hodnot, a soubory dat jsou blíže k normálnímu rozdělení.
Mnohem korektnější by samozřejmě bylo vztáhnout tyto úvahy
na soubor dat s vyloučenými krajními hodnotami. Tento krok
by však snížil rozptyl jen nepatrně.
Rychlosti cL u měření na třmenech v provozu je od 5677 ms-1
do 5790 ms-1, s průměrnou hodnotou 5729 ms-1. Průběh je na
obrázku 33 v příloze ke kap. 10. Zde byla všechna měření
provedena jen na jednom místě odlitku Ţ omezený rozptyl
hodnot. Ten je jen 247 ms-1, což je D4,3 % z průměrné
hodnoty. Oproti obrobené ploše na třmenech je tento rozptyl
větší o cca 2 %.
Oblast odlitku držáku je v rozmezí hodnot tvrdosti 179 až
219 HB (D 39 HB, tj. 20% z průměrné hodnoty), rychlost
ultrazvuku je od 5379 po 5721 ms-1 (D 342 ms-1, tj. 6,2% z průměrné
hodnoty) - graf 47 v příloze k této kapitole. Větší
rozptyl hodnot u tohoto měřením, než na třmenech v
provozu, je vlivem měření na několika místech na odlitku.
Ukazuje se tedy, že vliv místa měření má podstatnější
vliv na rozptyl hodnot, než vliv obrobené/neobrobené
plochy.
Rychlosti cL u měření na odlitku páky je od 5530 ms-1 po
5890 ms-1, s průměrnou hodnotou 5646 ms-1. Předcházející
závěr dokazují i tyto výsledky, kde je procentuální
rozptyl hodnot přes 6 % z průměrné hodnoty. Měření zde
probíhala na dvou místech na odlitku.
Další oblasti je Odlitek válce. Vzhledem k tomu, že se
jednalo o lití do kokily, dosahuje nejvyšších tvrdosti (až
415 HB) a tomu i odpovídající rychlostí cL (max. 5853 ms-1).
Ale jak svoji dolní hranici tvrdosti (min. 225 HB), tak
rychlosti ultrazvuku (min. 5625 ms-1) "zasahuje" do
oblastí pro odlitky lité do bentonitu. I zde platí, že co
oblast na odlitku to samostatný statistický soubor, lišící
se jak svým teplotním polem, chemickým složením, výslednou
strukturou, tvrdosti a rychlosti cL.
U válce je rozptyl hodnot v řezu A přes 45 % z průměrné
hodnoty. V řezu B je to 48 %, v řezu C je to 46 % a v dolním
vzorků řezu D je to přes 62 %. Tento značný rozptyl je
jednoznačně dán značně nehomogenním teplotním polem při
chladnutí odlitku v kokile a rovněž chemickou
nestejnorodosti, zvláště ve vzorku řezu D.
U měření rychlosti ultrazvuku je maximální rozdíl hodnot
1,8 % v řezu A, 1,5 % v řezu B, 1,6 % v řezu C a 2,5 % v řezu
D. Z toho vyplývá, že hodnoty změřené v povrchové vrstvě
vzorků (HB, chemické složení, . . .), mají větší
rozptyl hodnot, neboť vypovídají o slabé vrstvě v místě
měření. Zatímco měření vlastností přes stěnu odlitku
(rychlost ultrazvuku) zahrnují velké množství těchto
"povrchových vrstev" a jedná se o součet hodnot z
mnoha těchto vrstev, Z toho vyplývá jejich mnohem menší
rozptyl a menší citlivost na lokální odlišnosti.
Na grafu 46 a 47 v příloze jsou porovnány oblasti naměřených
hodnot tvrdosti a rychlosti ultrazvukových vln. Stručné
porovnání vybraných statistických charakteristik na všech
typech vzorků je v tabulce 16 a-d) v příloze.
Z grafu 46 je patrné, že přestože zde publikována data
pocházejí z odlitků ČSN 42 2305 o přibližně stejném
chemickém složení, není možné nalézt funkci, která by,
byť jen pro tento jeden typ materiálu a formovací směs (a
tudíž i rychlost odvodu tepla), popisovala dostatečně
spolehlivě (s vysokým korelačním koeficientem) chování
rychlosti podélných ultrazvukových vln v závislosti na
tvrdosti.
Větší rozptyly tvrdostí u litiny s kuličkovým grafitem
(než například u litiny s lupínkovým grafitem) spočívají
také v tom, že množství perlitu, který je přítomno ve
struktuře těchto odlitků, závisí na rychlosti procesu
tuhnutí, která je různá v různých tloušťkách stěn
(tyče o průměru 6 mm, třmen - čelo o tloušťce 8 mm a průměru
55 mm, třmen - ramena o tloušťce 8 mm, válec o průměru
65 mm).
Je známo, že například v tloušťkách 6 - 50 mm kolísá
obsah perlitu od 30 až do 70 %, což samozřejmě ovlivňuje
například i náklady na obrábění. Množství perlitu ale
také závisí na vzdálenosti mezi modulemi grafitu, větší
množství modulí znamená menší množství perlitu a toto
všechno podstatně ovlivňuje rozdíly ve tvrdosti.
U odlitku přilitých tyčí jsme měli k dispozici dostatek
dat, abychom mohli porovnat závislost mezi jednotlivými veličinami
a rychlosti průchodu vln. Na grafu 48 v příloze je závislost
mezi tvrdostí a rychlostí cL. Korelační vztah je HB =
0,0965 cL [ms-1] - 349,1 [HB] pro interval rychlostí od 5400
po 5700 ms-1, se slabým koeficientem 0,22. Závislost pro
pevnost Rm a rychlost vln je silnější: Rm = 0,8471 cL [ms-1]
- 4154,9 [MPa] s vysokým koeficientem 0,75 - grafu 53.
Pro ostatní veličiny jsou koeficienty korelace velmi slabé.
Schopnost popsat rychlost šíření ultrazvukových vln na základě
více vlastností byl proveden v tab. 39. Za použití analýzy
hlavních komponent jsme přišli k závěru, že pro dostatečný
popis cL je potřeba použít alespoň čtyř veličin:
Pevnosti v tahu, tvrdosti, tažnosti a obsahu uhlíku.
Před započetím měření nelze věrohodně určit velikosti
rychlosti ultrazvuku a jí odpovídající hodnoty tvrdosti,
pevnosti, tažnosti či dalších vlastností matrice. Vždy před
zavedením této metody do výstupní kontroly odlitků je třeba
provést přímé měření těchto vlastností na statisticky
významném počtu vzorků a teprve následně považovat tuto
metodu za věrohodnou a průkaznou pro zjišťování stavu
kvality odlitků, s periodicky se opakujícími přímými
kontrolami pro ověření trendu vlastností.
Zároveň nelze přenášet korelační vztahy zjištěné na
jednom typu odlitku na druhý, ba ani v rámci jednoho odlitku
z jedné oblasti, pro kterou byly splněny výše zmíněné
podmínky měření, na oblast druhou, jak dokázal případ třmenů
- ramen, kdy i v rámci velmi blízkého okolí získáváme
odlišné statistické soubory. Neboť by bylo nutné dodržet
stejnou tloušťku materiálu, stejné rychlosti chladnutí ve
formě a přibližně stejné chemické složení, což v
provozních podmínkách není vždy možné.
Zároveň se prokázalo, že statistické soubory většinou
nepocházejí z normálního rozdělení, což signalizuje na
další nepodchycené parametry ovlivňující měření a to
i v laboratorních podmínkách.
Stáhnout
PDF soubor s kap 11 (1,2 MB)
Předcházející kapitola: 10.
Experimentální část
Následující kapitola: 12.
Simulace TP a mikrostruktury
Upozornění: Pokud
použijete část z moji disertační práce, dodržujte
Autorský zákon a dbejte na správnost citací
|