ŠKODA AUTO a.s. Vysoká škola
DIPLOMOVÁ
2011
PRÁCE
Lukáš Kotfald
ŠKODA AUTO a.s. Vysoká škola
Studijní program: N6208 Ekonomika a management Studijní obor: 6208T088 Podniková ekonomika a management provozu
SYSTÉM ZAJIŠŤOVÁNÍ KVALITY NAKUPOVANÝCH DÍLŦ VE ŠKODA AUTO a.s.
Lukáš KOTFALD
Vedoucí práce: doc. Ing. Eva Jarošová, CS
ANOTAČNÍ ZÁZNAM AUTOR
Lukáš Kotfald
STUDIJNÍ OBOR
Podniková ekonomika a management provozu
NÁZEV PRÁCE
Systém zajišťování kvality nakupovaných dílŧ ve ŠKODA AUTO a.s.
VEDOUCÍ PRÁCE
doc. Ing. Eva Jarošová, CSc.
INSTITUT
IPT
POČET STRAN
62
POČET OBRÁZKŦ
29
POČET TABULEK
2
POČET PŘÍLOH
1
STRUČNÝ POPIS
KLÍČOVÁ SLOVA
ROK ODEVZDÁNÍ
2011
Diplomová práce se zaměřuje na aplikaci vybraných statistických nástrojŧ, v podobě regulačních diagramŧ a indexŧ zpŧsobilosti, při zajišťování kvality nakupovaných dílŧ ve společnosti ŠKODA AUTO a.s.. V první části je popsán soubor nástrojŧ na udržování stability procesu a zlepšování jeho zpŧsobilosti prostřednictvím redukce variability. Uvedené nástroje jsou využity v případových studiích, kde je ověřena stabilita a zpŧsobilost procesu u kapot a přitahovače. V případě kapot, se sledovaným znakem nanášená vrstva lepidla, je zjištěno, že v bodě č. 13 není proces považován na stabilní. Je doporučen okamžitý zásah do procesu. U přitahovače jsou sledovány hned tři regulované veličiny – síla, zdvih a čas. Regulovaná veličina zdvih vypovídá o nedostatečné přesnosti měření. V případě regulované veličiny čas je zjištěna nenormalita, nestabilní a nezpŧsobilý proces. Z dŧvodŧ značné variability je doporučen zásah do procesu.
Statistická regulace procesu; Zpŧsobilost; Systém měření.
ANNOTATION AUTHOR
Lukáš Kotfald
FIELD
Business economics and production management
THESIS TITLE
Quality assuarance system for purchased parts of SKODA AUTO a.s.
SUPERVISOR
doc. Ing. Eva Jarošová, CSc.
INSTITUTE
IPT
YEAR
NUMBER OF PAGES
62
NUMBER OF PICTURES
29
NUMBER OF TABLES
2
NUMBER OF APPENDICES
1
SUMMARY
KEY WORDS
2011
This thesis focuses on the application of selected statistical tools in the form of control charts and capability indices, while ensuring quality of purchased parts in SKODA AUTO a.s.. The first section describes a set of tools to maintain process stability and improve its capacity through the reduction of variability. These tools are used in case studies, where stability and capability of a bonnet and an adductor are tested. In the case of the bonnet, with the observed feature a silted layer of glue is found that the process at point No. 13 is not regarded as stable. An immediate action in the process is recommended. In the case of the adductor even three controlled variables are observed - the power stroke, and time. Controlled variable stroke indicates the lack of measurement system. Controlled variable time is detected abnormal, unstable and non-capability process. Because of its considerable variability an action to the process is recommended.
Statistical Process Control; Capability; Measurement Systems Analysis.
2
Prohlašuji, že jsem diplomovou práci vypracoval samostatně s použitím uvedené literatury pod odborným vedením vedoucího práce. Prohlašuji, že citace použitých pramenu je úplná a v práci jsem neporušil autorská práva (ve smyslu zákona č. 121/2000 Sb., o právu autorském a o právech souvisejících s právem autorským).
V Mladé Boleslavi, dne ...........................
...........................
3
Děkuji doc. Ing. Evě Jarošové, CSc. za odborné vedení diplomové práce, poskytování rad a informačních podkladů.
4
1. ÚVOD ………………………………………………………………………………....9 2. STATISTICKÁ REGULACE PROCESU………………………………………… 11 2.1. Etapy statistické regulace procesu……………………………………….. 12 2.1.1. Etapa analýz procesu…………………………………………………...12 2.1.2. Etapa udržování procesu……………………………………………….13 2.1.3. Etapa zlepšování procesu……………………………………………... 13 2.2. Nástroje statistické regulace procesu……………………………………. 13 2.2.1. Histogram………………………………………………………………... 14 2.2.2. Pravděpodobnostní graf……………………………………………...... 15 2.2.3. Ověření normality pomocí testů………………………………………. 16 2.2.4. Regulační diagramy…………………………………………………..... 16 2.3. Regulace měřením…………………………………………………………… 19 2.3.1. Diagramy pro aritmetické průměry a rozpětí………………………… 20 2.3.2. Diagramy pro aritmetické průměry a směrodatné odchylky……….. 20 2.3.3. Diagramy pro mediány a rozpětí……………………………………… 21 2.3.4. Diagramy pro individuální hodnoty a klouzavá rozpětí……………... 22 2.4. Zpŧsobilost procesu……………………………………………………….... 23 2.4.1. Výpočet indexů způsobilosti při normálním rozdělení……………… 25 2.4.1.1. Indexy způsobilosti………………………………………………25 2.4.1.2. Indexy výkonnosti………………………………………………. 27 2.4.2. Výpočet indexů způsobilosti při nenormálním rozdělení…………… 28 2.4.2.1. Box-Coxova transformace……………………………………... 28 2.4.2.2. Clementsova metoda……………………………………………29 3. SYSTÉM MĚŘENÍ…………………………………………………………………. 31 4. PŘEDSTAVENÍ SPOLEČNOSTI………………………………………………… 34 5. KAPOTY…………………………………………………………………………….. 36 5.1. Analýza procesu v bodě č. 13……………………………………………… 38 5.2. Analýza procesu v bodě č. 27……………………………………………… 42 6. PŘITAHOVAČ……………………………………………………………………… 47 6.1. Analýza procesu s regulovanou veličinou síla…………………………. 49 6.2. Analýza procesu s regulovanou veličinou zdvih……………………….. 52 6.3. Analýza procesu s regulovanou veličinou čas…………………………. 53 7. ZÁVĚR………………………………………………………………………………. 58
5
Seznam použitých zkratek apod. a podobně ISO
mezinárodní organizace pro standardizaci
např. například obr.
Obrázek
P/T
poměr přesnost/ tolerance
QMS systém řízení kvality resp. respektive R&R repeatability a reproducibility SPC statistická regulace procesu tab.
tabulka
tj.
to jest
tzv.
takzvaný
VW
Volkswagen AG
6
Seznam požitých symbolŧ Cp
index způsobilosti
C pk
index způsobilosti
C pm
index způsobilosti
CC pk index způsobilosti
CL
centrální přímka
CPL
dolní index způsobilosti
CPU horní index způsobilosti g1
šikmost
g2
špičatost
k
počet podskupin
LCL
dolní regulační mez
LL x
dolní přirozená mez
Lp
dolní kvantil
LSL
dolní specifikační mez
LWL dolní varovná mez
M
medián v podskupině
M
průměrná hodnota mediánů v podskupinách
MR
klouzavá rozpětí
n
rozsah podskupiny
Pp
index výkonnosti
Ppk
index výkonnosti je-li brána v úvahu poloha procesu
Ppm
index výkonnosti v případě asymetrických mezí
PPpk
index výkonnosti dle Kane
R
rozpětí v podskupině
R0
základní hodnota pro rozpětí
R
průměrná hodnota rozpětí R v podskupinách
s
výběrová směrodatná odchylka
s0
základní hodnota pro průměr směrodatných odchylek v podskupinách
T
cílová hodnota procesu
7
UCL horní regulační mez USL
horní specifikační mez
UL x
horní přirozená mez
Up
horní kvantil
UWL horní varovná mez. x
hodnota měřeného znaku jakosti
X0
základní hodnota pro průměr procesu
střední hodnota výrobního procesu
0
základní hodnota pro parametr polohy procesu
směrodatná odchylka procesu
0
základní hodnota pro průměr průměrů v podskupinách
ˆ
odhad skutečné hodnoty směrodatné odchylky výrobního procesu uvnitř podskupin
T2
celková variabilita
8
1. ÚVOD Podmínkou úspěchu výrobních podniků na trhu je nabídnutí bezchybného P Podmínkou úspěchu každého podniku je nabídnutí bezchybného produktu. I přes veškeré úsilí však může být zájem zákazníka ohrožen. Výrobě tzv. zmetků se bohužel nelze vyhnout. Je však důležité, aby se tyto výrobky nikdy nedostaly až k finálnímu spotřebiteli. K ověření správného chodu výrobního procesu podniky využívají různých kontrol. Tradiční způsob zajištění jakosti vyráběných produktů spočívá v kontrole dodatečné, která je provedena po ukončení výrobního procesu. Jedná se však o kontrolu neekonomickou, neboť v případě neshodného výrobku byly náklady vynaložené na výrobek zbytečné, a pasivní, neboť není schopna upozornit na produkci neshodných výrobků. Využitím dodatečných kontrol náklady podniku narůstají. Zbytečnému vynakládání finančních prostředků lze zabránit, věnuje-li společnost zvláštní pozornost analýze procesu, jeho regulaci a následnému zhodnocení. Výsledkem může být dosažení stabilního výrobního procesu, výhod v podobě
vyšší
produktivity,
nižších
výrobních
nákladů
či
eliminování
nehospodárnosti. Mezi nejčastější preventivní nástroje řízení jakosti patří statistická regulace procesu. Statistická regulace procesu slouží k včasnému odhalení jednotlivých odchylek od stanovené úrovně sledovaného znaku a udržení procesu v požadovaném stavu. Podstatou statistické regulace procesu je sledování, kontrola, ale i případná intervence, je-li odhalena odlišnost v chování procesu. Důvodem zásahu bývá působení řady vlivů majících negativní vliv na výrobní proces. Působením těchto vlivů dochází ke zvýšení variability výrobního procesu. Pokud jde o variabilitu v rámci určitých „malých“ mezí, neznamená to, že požadavky zákazníka jsou ohroženy. V opačném případě se produkt může stát nepřijatelný, což může vést až k samotné ztrátě zákazníka. Součástí statistické regulace procesu je hodnocení způsobilost procesu. Způsobilost procesu je schopnost procesu vytvářet produkty, které odpovídají požadavkům kvality. Znalost způsobilosti procesu je podmínkou zlepšování jakosti. Informace o způsobilosti jsou určeny převážně zákazníkovi, který je tímto způsobem informován o schopnosti procesu vyhovět předepsaným kritériím jakosti.
9
Statistická regulace procesu a hodnocení způsobilosti procesu je obsahem předkládané práce. Cílem práce je výběr vhodných postupů při regulaci procesu a hodnocení jeho způsobilosti a jejich aplikace na proces zajišťování kvality nakupovaných dílů ve společnosti ŠKODA AUTO a.s.. V první části je popsán soubor nástrojů na udržování stability procesu a zlepšování jeho způsobilosti prostřednictvím redukce variability. Uvedené nástroje jsou využity v druhé části práce, kde je ověřena stabilita a způsobilost procesu u vybraných nakupovaných dílů dodávaných do společnosti ŠKODA AUTO a.s..
10
2. Statistická regulace procesu Statistická regulace procesu je soubor nástrojů na udržování stability procesu a zlepšování jeho způsobilosti prostřednictvím redukce variability. 13 Cílem každé společnosti je úspěšný finální produkt. Toho lze dosáhnout tehdy, je-li společnost schopna uspokojit očekávání svých zákazníků. Tato očekávání však budou naplněna jen v případě stabilního výrobního procesu. Efektivní
způsoby
zabezpečení
jakosti
zabraňují
nadbytečnému
plýtvání
prostředků potřebných při výrobním procesu. Zabránit plýtvání lze cestou neustálého monitorování procesu a jeho pečlivou analýzou. Účelem regulace procesu je působit na výrobní proces takovým způsobem, aby výstupy dosahovaly požadovaných vlastností. Statistická regulace procesu umožňuje taková opatření, která zajistí udržení dlouhodobého stabilního výrobního procesu. Stabilním výrobním procesem lze nazvat takový proces, během kterého je vyráběn produkt s přijatelnou
variabilitou ukazatelů
kvality. Variabilita,
jejímž důsledkem je
neexistence dvou stejných výrobků, je považována za výsledek působení řady vlivů, ovlivňujících proces i za relativně standardních podmínek. Jelikož je možné variabilitu dlouhodobě sledovat, lze na základě získaných informací odhadnout chování výrobního procesu i v budoucnosti. Obecně existují dva druhy příčin, které ovlivňují variabilitu výstupu procesu: náhodné a vymezitelné příčiny Náhodné příčiny jsou stálou součástí procesu nebo systému a ovlivňují všechny složky procesu.14 Stabilní proces je takový systém, ve kterém se objevují jen náhodné příčiny způsobující pouze přirozenou variabilitu. Výskyt této variability v čase můžeme považovat za neměnný. V případě výstupu ze stabilního procesu mluvíme o minimální variabilitě. Tuto variabilitu je možné předpokládat v rámci určitých vymezených hranic. Vymezitelné příčiny nejsou stálou součástí procesu nebo systému, ale vznikají v důsledku specifických okolností.15 13
TEREK, M. HRNČIAROVÁ, L. Štatistické riadenie kvality. Bratislava: Iura Edition, 2004, ISBN 80-89047-97-1, s. 21 14 TEREK, M. HRNČIAROVÁ, L. Štatistické riadenie kvality. Bratislava: Iura Edition, 2004, ISBN 80-89047-97-1, s. 21 11
Opakem stabilního procesu je proces nestabilní. Na takovýto proces působí jak náhodné tak i vymezitelné příčiny. U nestabilního procesu se variabilita v čase mění. Aby se z nestabilního procesu stal proces stabilní, je nutné identifikovat a posléze odstranit vymezitelné příčiny. Výstup z každého výrobního procesu je charakterizován tzv. regulovanou veličinou. V praxi se můžeme setkat se dvěma typy regulací. Vždy záleží na tom, zdali je regulovanou veličinou kvantitativní či kvalitativní znak. V případě kvantitativního znaku hovoříme o regulaci měřením. Jedná-li se o znak kvalitativní, hovoříme o regulaci srovnáním.
2.1. Etapy statistické regulace procesu Statistická regulace procesu probíhá jako třístupňový cyklus skládající se z analýzy procesu – udržování procesu – zlepšování procesu. 2.1.1. Etapa analýz procesu Etapa analýz procesu usiluje o uvedení procesu do požadovaného stavu (proces je nastaven na požadovanou hodnotu a osciluje kolem ní s určitým vyhovujícím rozptylem) a stabilního stavu (kolísání vyvoláno pouze náhodnými příčinami). V první fázi etapy analýz jde zejména o diagnózu procesu, tedy stanovení, je-li proces schopen regulace či nikoliv. Fáze druhá je založena na napravování procesu, tedy odstraňování vymezitelných příčin ovlivňujících regulovanou veličinu, popřípadě centrování procesu. Výstupem etapy analýz procesu je stanovení odpovídajících mezí, přičemž vzdálenost mezí souvisí s velikostí kolísání daného procesu. Na počátku jsou stanoveny pokusné regulační meze, na jejichž základě se posoudí působení vymezitelných příčin. V případě, že je proces stabilní, leží body uvnitř regulačních mezí a pokusné meze lze považovat za platné i nadále. V opačném případě je nutné vymezitelnou příčinu odhalit, vyloučit podskupiny s vymezitelnými příčinami a přepočítat regulační meze.
15
TEREK, M. HRNČIAROVÁ, L. Štatistické riadenie kvality. Bratislava: Iura Edition, 2004, ISBN 80-89047-97-1, s. 21 12
2.1.2. Etapa udržování procesu Cílem etapy udržování procesu je udržení procesu ve stabilním stavu za pomoci monitorování regulačním diagramem, z důvodu včasného zásahu do procesu ještě před okamžikem produkce neshodných výrobků. Neustálé sledování je nezbytné, neboť i u momentálně stabilního procesu může časem dojít k překročení regulačních mezí, kdy některá z náhodných příčin se může stát příčinou vymezitelnou (například v případě jejího zesílení). Jelikož v dané fázi pracujeme s procesem regulovaným, je vhodné zavádět složitější typy regulačních diagramů, mezi které patří regulační diagramy pro aritmetický průměr s výstražnými mezemi, jež citlivěji zachycují posun úrovně procesu (ČSN ISO 7873) či regulační diagramy využívající kumulovaná data, například Cusum (ČSN 01 0266). 2.1.3. Etapa zlepšování procesu Cílem etapy zlepšování procesu je redukovat variabilitu a nadále zlepšovat jakost výstupu. Na základě hlubšího poznání jsou identifikována kritická místa procesu, na která je nezbytné zaměřit pozornost, a to z důvodu snížení kolísání vyvolaného náhodnými příčinami. Během etapy zlepšování procesu je proces hodnocen pomocí indexů způsobilosti, pomocí kterých se redukce variability projeví.
2.2. Nástroje statistické regulace procesu Nástroje statistické regulace procesu lze rozdělit na nástroje pro analýzu procesu a nástroje sloužící k samotné regulaci. Nástroje pro analýzu procesu mají za úkol odhalit, je-li proces schopen regulace či nikoliv. K samotné regulaci procesu se používají regulační diagramy. Jednou z podmínek regulace schopného procesu je ověření předpokladu normality. K ověření předpokladu normality slouží grafické metody a statistické testy.
Mezi
nejčastěji
pravděpodobnostní
graf.
používané V případě
grafické
metody
statistických
Shapirův-Wilkův test, vyznačující se dobrou účinností.
13
testů
patří se
histogram často
či
používá
2.2.1. Histogram Histogram (obr. 1) zobrazuje rozdělení četností sledovaného znaku, a informuje tak o charakteru proměnlivosti procesu, jeho přesnosti či centrování.
Obr. 1 Histogram regulované veličiny vrstva lepidla v bodě č. 27, Statgraphics Z důvodů reprezentativnosti výběru ze základního souboru se doporučuje při tvorbě histogramu vycházet nejméně ze 100 údajů. Pro volbu vhodného počtu intervalů se používá vzorec
1 3,3 * log n
(1)
kde n je rozsah podskupiny. Histogram poskytuje čtyři základní informace o výrobním procesu. 1. Na základě tvaru histogramu je často možné posoudit, zdali na proces působí pouze náhodné či i vymezitelné příčiny. 2. U regulovaných veličin s normálním rozdělením lze odhadnout, v jakém rozmezí se bude nacházet většina jejich hodnot. 3. Na základě polohy nejvyššího sloupce rozpoznáme, je-li proces centrován či nikoliv. 4. Zakreslíme-li do histogramu meze určené specifikací, můžeme posoudit způsobilost procesu.
14
2.2.2. Pravděpodobnostní graf Pravděpodobnostní graf (obr. 2) slouží k posouzení, zdali data pocházejí z normálního rozdělení či nikoliv. Na vodorovné ose jsou vyneseny naměřené hodnoty xi . Y –ová souřadnice bodů odpovídá distribuční funkci normovaného normálního rozdělení odvozené z pořadí hodnot x. V grafu je zachycena regresní přímka odpovídající rovnici
Exi u i
(2)
kde je skutečná střední hodnota výrobního procesu a skutečná hodnota směrodatné odchylky výrobního procesu uvnitř podskupin, u i značí kvantit normovaného normálního rozdělení. V případě, že se jedná o výběr z normálního rozdělení, měly by body ležet v přímce. Z důvodů ověření správnosti závěrů bývá graf doplněn některým ze statistických testů normality.
Obr. 2 Pravděpodobnostní graf regulované veličiny síla, Statgraphics
15
2.2.3. Ověření normality pomocí testŧ Nulovou hypotézou je tvrzení, že výběr pochází ze základního souboru s normálním rozdělením. Často používaným testem je Shapirův-Wilkův test s testovou statistikou W [20]
u x W u x x 2
i
i
(3)
2
2 i
i
která představuje korelační koeficient mezi naměřenými hodnotami xi a kvantily normovaného normálního rozdělení u , x a u i , kde i značí pořadí hodnoty x po seřazení podle velikosti. Platí, že čím blíže hodnotě jedna, tím spíše lze předpoklad normality považovat za splněný. Statistické programy uvádějí kromě hodnoty testové statistiky i tzv. phodnotu, která se porovnává s hladinou významnosti (volíme hodnoty 0,05, 0,01 a 0,001). Je-li p-hodnota nižší než zvolená hladina významnosti, hypotéza o normalitě se zamítá. 2.2.4. Regulační diagramy Regulační diagramy jsou považovány za nejdůležitější nástroj statistické regulace procesu, na jehož základě je možné operativně určovat, je-li proces stabilní
či
nikoliv.
Regulační
diagramy
vycházejí
z posloupností
výběrů
prováděných v pravidelných časových intervalech. Charakteristika zvoleného výběru, např. průměr či rozpětí, vypovídá o aktuálním stavu procesu vůči regulačním mezím. Regulační meze jsou určeny na základě variability procesu při působení náhodných příčin. Za tohoto předpokladu leží výběrové charakteristiky s největší pravděpodobností uvnitř zmíněných mezí. Výběrové charakteristiky tedy znázorňují průběh procesu v čase. Regulační diagram je zachycen na obr. 3. Vodorovná osa znázorňuje pořadí podskupin. Na svislé ose se vynášejí hodnoty sledované charakteristiky.
16
Obr. 3 Stabilní proces regulované veličiny vrstva lepidla v bodě č. 27, Statgraphics Centrální přímka CL zobrazuje požadovanou, tzv. referenční hodnotu regulované veličiny v případě, kdy je proces stabilní. Referenční hodnotou může být hodnota daná technickým předpisem či hodnota založená na minulé zkušenosti s daným procesem. Dolní regulační mez LCL a horní regulační mez UCL vymezují pásmo působení pouze náhodných příčin variability procesu. Zmíněné meze jsou považovány za základní rozhodovací kritérium, zdali zasáhnout do procesu či nikoliv. V některých případech se zakreslují i výstražné meze. Regulační diagram lze považovat za grafické zobrazení testu statistické hypotézy, týkající se hodnoty parametru rozdělení regulované veličiny. Za nulovou hypotézu lze považovat situaci, kdy je proces ve statisticky zvládnutelném stavu. Opakem nulové hypotézy je hypotéza alternativní, která odpovídá statisticky nezvládnutému procesu. Za obor přijetí nulové hypotézy považujeme pásmo v rozmezí akčních mezí LCL a UCL regulačního diagramu. Oblast mimo dolní a horní regulační mez odpovídá oboru kritickému. Pravděpodobnost chyby I. druhu, tzv. riziko zbytečného signálu, je dána vzdáleností akčních mezí. Jde o signál v případě, kdy vymezitelná příčina neexistuje, tj. parametry procesu se nezměnily. O zbytečném signálu hovoříme tehdy, zůstává-li proces ve statisticky zvládnutém stavu, avšak bod náhodně 17
padne mimo regulační meze. Důsledkem je zbytečný zásah, který může působit kontraproduktivně. V případě chyby II. druhu, tzv. chybějícího signálu, dochází k opačné situaci, kdy se změnily parametry procesu, avšak výběrová charakteristika se ještě nachází uvnitř regulačních mezí. V daném okamžiku regulační diagram včas neodhalí změnu v procesu vyvolanou vymezitelnými příčinami. Riziko chybějícího signálu můžeme zmenšit pomocí testů zvláštních příčin. K odhalení vymezitelných příčin se používá Ishikawův diagram (obr. 4). Ishikawův diagram, známý také jako diagram příčin a následků nebo fishbone diagram příčin a následků, zachycuje potenciální zdroje vymezitelných příčin, které lze následně analyzovat. Při tvorbě diagramu se řešený problém znázorní jako hlava rybí kostry. Následně se formou brainstormingu identifikují možné příčiny problému, které tvoří rybí kosti. Tyto kosti se dále dělí na jednotlivé kostičky.
Zdroj: Převzato z16 Obr. 4 Diagram příčin a následků
16
Interní materiály Škoda Auto a.s. 18
Při implementaci regulačních diagramů lze použít následující pravidla: Proces se považuje za statisticky zvládnutý za předpokladu, že všechny body leží uvnitř regulačních mezí LCL a UCL. Za dané situace není třeba zásahu, neboť nepůsobí žádná vymezitelná příčina. Pokud některý z bodů leží vně výstražných mezí, lze proces považovat za statisticky nezvládnutý. Za dané situace je třeba identifikovat vymezitelnou příčinu a učinit opatření vedoucí k .jejímu odstranění. V případě existence výstražných mezí dále platí: Pokud dva body za sebou leží mezi výstražnou mezí LWL a LCL, popřípadě UWL a UCL, je doporučeno zasáhnout do procesu následujícím způsobem: Je proveden další výběr (bez ohledu na kontrolní interval) a nastává jedna ze dvou možností., neboť Pokud nový bod leží mezi výstražnými mezemi, není třeba zasáhnout do procesu. V případě bodu vně výstražných mezí, je zásah nezbytný, neboť na proces pravděpodobně působí vymezitelná příčina.
2.3. Regulace měřením Regulace měřením plní funkci regulace jako metody prevence, neboť je schopna včas rozpoznat a upozornit na zhoršující se kvalitu, a to ještě před samotným okamžikem výroby neshodných výrobků. V případě regulace měřením se vždy konstruují dva diagramy - jeden pro polohu, druhý pro variabilitu. V praxi se využívají následující kombinace:
Diagramy pro aritmetické průměry a rozpětí.
Diagramy pro aritmetické průměry a směrodatné odchylky.
Diagramy pro mediány a rozpětí.
Diagramy pro individuální hodnoty a klouzavá rozpětí. U regulačních diagramů mohou být předem stanoveny základní hodnoty.
Základní hodnoty představují požadovanou střední hodnotu a kolísání procesu. Nejsou-li základní hodnoty předem známy, je nezbytné provést minimálně 20 výběrů. Stačí malý rozsah výběrů (n = 2 až 5). Interval mezi jednotlivými výběry závisí na charakteru procesu. Obecně platí, že rozsah intervalu v přípravné etapě bývá kratší. 19
Je třeba, aby výběry (podskupiny) byly tvořeny tak, aby v rámci nich nedocházelo k podstatným změnám procesu. To znamená, aby v rámci výběru nepůsobily vymezitelné příčiny. Z tohoto důvodu je doporučeno jednotlivé výběry skládat z prvků, vyrobených v chronologickém sledu. Jedná se tedy o vytváření již zmíněných logických podskupin. 2.3.1. Diagramy pro aritmetické prŧměry a rozpětí Centrální přímka a regulační meze diagramu pro rozpětí se určí podle vzorců [3] CLR R
(4)
UCLR R D4
(5)
LCL R R D3
(6)
V diagramu pro průměr platí [3] CL X
(7)
UCLX X A2 R
(8)
LCL X X A2 R
(9)
Součinitele A2, D3 a D4, závisející na rozsahu podskupin, jsou uvedeny v tabulce přílohy č. 1. 2.3.2. Diagramy pro aritmetické prŧměry a směrodatné odchylky Pro výběry o rozsahu 10 a větší je vhodnější místo rozpětí používat směrodatnou odchylku. Pro směrodatnou odchylku platí
s
1 k sj k j 1
kde k je počet podskupin.
20
(10)
Centrální přímka a regulační meze jsou dány vztahy [3] CLS s
(11)
UCLS B4 s
(12)
LCLS B3 s
(13)
Používáme-li k měření variability v podskupinách směrodatnou odchylku, jsou regulační meze v diagramu pro průměr dány vzorci [3]
UCL X X A3 s
(14)
LCL X X A3 s
(15)
2.3.3. Diagramy pro mediány a rozpětí Diagram je vhodné použít tehdy, provádíme-li výpočty provádíme ručně. Určení mediánu je jednodušší, ovšem za podmínky rozsahu výběru n = 3 či n = 5. Diagram pro rozpětí je stejný jako výše. Centrální přímka a regulační meze diagramu pro medián se určí podle vzorců [3] CLM M
(16)
UCLM M A4 R
(17)
LCL M M A4 R
(18)
Hodnoty konstanty A4 jsou zobrazeny v tabulce 1. Tab. 1 Hodnoty konstanty A4 n
2
3
4
5
6
A4
1,88
1,19
0,8
0,69
0,55
7
8
9
10
0,51
0,43
0,41
0,36
17
Zdroj: Upraveno z
17
ČSN ISO 8258 Shewhartovy regulační diagramy. Praha: Český normalizační institut, 1998, s. 13 21
2.3.4. Diagramy pro individuální hodnoty a klouzavá rozpětí Není-li v praxi možné vytvořit logické podskupiny, používají se diagramy pro individuální hodnoty a klouzavá rozpětí. Protože nemáme k dispozici logické podskupiny, nemůže se měřit kolísání uvnitř podskupin, a proto se tvoří umělé podskupiny ze dvou po sobě jdoucích hodnot a k měření variability se používají klouzavá rozpětí. Klouzavá rozpětí jsou absolutní hodnoty rozdílu mezi sousedními hodnotami, tj.
MR j x j x j 1
(19)
kde j 2,...,k. Horní regulační mez pro diagram rozpětí získáme podle vzorce UCLR MR D4
(20)
kde D4 = 3,267 odpovídá rozsahu podskupin n = 2. Dolní regulační mez je nulová viz tabulka přílohy č. 1. Centrální přímku a regulační meze pro diagram individuálních hodnot získáme podle vzorců [3] CL X X
(21)
UCL X X MR A3
(22)
LCL X X MR A3
(23)
Z důvodů prokázání stability by měl být diagram individuálních hodnot sestrojen na základě dostatečného počtu hodnot. V [6] se doporučuje 100 hodnot. Při konstrukci regulačních diagramů, u kterých předem známe základní hodnoty, jsou regulační meze odvozeny z těchto základních hodnot. Výpočet centrální přímky a regulačních mezí lze provést dle tabulky 2, kde základní hodnoty pro polohu procesu jsou označeny X 0 nebo 0 . Základní hodnoty pro kolísání procesu jsou označeny R0 , s 0 nebo 0 . Hodnoty součinitelů uvedených v tabulce lze nalézt v příloze č. 1.
22
Tab. 2 Vzorce pro výpočet regulačních mezí při předem stanovených základních hodnotách polohy střední hodnoty a kolísání procesu Výběrová charakteristika
Centrální přímka
UCL a LCL
X0 nebo µ0
X0 ± A*σ0
R
R0 nebo d2*σ0
D2*σ0, D1*σ0,
s
s0 nebo C4*σ0
B6*σ0, B5*σ0,
X
X0 nebo µ0
X0 ± 3*σ0
MR
R0 nebo d2*σ0
D2*σ0, D1*σ0,
Zdroj: Upraveno z18
2.4. Zpŧsobilost procesu Způsobilost procesu je schopnost procesu vyhovět daným specifikacím. Informace o způsobilosti jsou určeny převážně zákazníkovi, který je tímto způsobem informován o přesnosti procesu v budoucnosti. Přesnost procesu souvisí s variabilitou, která představuje kolísání způsobené pouze náhodnými příčinami. Při hodnocení způsobilosti pomocí nejznámějších indexů způsobilosti se předpokládá normální rozdělení sledované veličiny, avšak existují i postupy pro jiná než normální rozdělení. Technické specifikace korespondují s požadavky externích zákazníků, a proto jsou určujícím faktorem požadované přesnosti procesu. Technické specifikace bývají v praxi zachyceny různými způsoby. Cílová (nominální) hodnota může nebo nemusí být udána. Vždy jsou však udány toleranční meze (oboustranné nebo jednostranné). Nominální hodnota bývá často umístěna ve středu intervalu. Inherentní
variabilitu
měříme
pomocí
směrodatné
odchylky,
kterou
odhadneme pomocí některého z následujících vzorců podle toho, jaký regulační diagram je používán při předcházející regulaci procesu [3].
ˆ s
18
R d2
HŮLOVÁ, M. JAROŠOVÁ, E. Statistické metody v managementu kvality, enviromentu a bezpečnosti. Praha: Oeconomica, 2004, ISBN 80-245-0691-2, s. 26 23
(24)
ˆ
s C4
(25)
V prvním případě se vychází z rozpětí, v druhém případě ze směrodatných odchylek v podskupinách. Součinitele d2 a C4, které se používají při konstrukci mezí v regulačních diagramech, jsou uvedeny v tabulce přílohy č. 1. Směrodatnou odchylku lze také odhadnout na základě vzorce [3] k
ˆ
s i 1
2 i
(26)
k
Pro grafické vyjádření se určují tzv. přirozené meze dané vzorci [3]
ULx X 3 ˆ
(27)
LL x X 3 ˆ
(28)
Meze určené specifikací a přirozené meze jsou zakresleny v obr. 5, ve kterém platí, že skutečné rozmezí kolísání je menší než dovolené rozmezí kolísání. Může však nastat i situace, kdy skutečné rozmezí kolísání je rovno dovolenému rozmezí kolísání či větší než dovolené rozmezí kolísání.
LSL
LLx
ULx
USL
Obr. 5 Histogram znaku vrstva lepidla, Statgraphics
24
Poloha průměru procesu a přirozených mezí vůči mezím daným technickou specifikací graficky určuje způsobilost procesu. 2.4.1. Výpočet indexŧ zpŧsobilosti a výkonnosti při normálním rozdělení Rozdíl mezi indexy způsobilosti a výkonnosti spočívá ve způsobu odhadu směrodatné odchylky procesu. U indexů způsobilosti se směrodatná odchylka procesu odhaduje pomocí vzorce (28). Splňuje-li proces předpoklad normality, používají se k výpočtu způsobilosti procesu následující indexy. 2.4.1.1.
Indexy zpŧsobilosti
Indexy způsobilosti slouží k porovnání skutečných a požadovaných rozmezí kolísání v číselné formě. Index Cp
USL LSL 6
(29)
lze považovat za základní charakteristiku způsobilosti procesu [6]. V případě Cp < 1 je skutečné rozmezí kolísání je větší než dovolené rozmezí kolísání a proces považován za nezpůsobilý. Dospějeme-li k tomuto výsledku v praxi, je nezbytné provést 100% kontrolu vyrobených produktů, nalézt pro naše výrobky jiného zákazníka s „volnějšími“ technickými požadavky či docílit snížení variability pomocí úpravy procesu. Pokud Cp = 1, tedy skutečné rozmezí kolísání se rovná dovolenému rozmezí kolísání, je proces považován za nezpůsobilý. Je-li 1 < Cp < 1,33, hodnotíme proces jako podmíněně způsobilý. Je-li Cp ≥ 1,33, je proces považován za způsobilý. V případě indexu Cpk
je brána v
úvahu poloha procesu, a to
prostřednictvím vzdálenosti střední hodnoty od cílové hodnoty technických specifikací. Index Cpk získáme dle vzorce [6] C pk min CPU , CPL
kde
25
(30)
CPU
USL 3
(31)
CPL
LSL 3
(32)
Rovnají-li se obě vypočtené hodnoty, je zřejmé, že proces je centrován. V opačném případě není střední hodnota uprostřed tolerančního pole. Zatímco index Cp informuje o potenciální způsobilosti procesu, hovoříme v případě indexu Cpk o přesnosti aktuální. Jsou-li obě hodnoty totožné, je proces centrován. Není-li tomu tak, vždy platí, že Cpk < Cp, a tedy skutečná střední hodnota procesu je posunuta. Je-li v praxi stanovena pouze jedna z mezí, neboť stanovení meze druhé nemá reálný smysl, používá se k výpočtu indexu Cpk vždy jeden ze vzorců (33) a (34). Je-li v technické specifikaci cílová hodnota T odlišná od středu tolerančního pole, lze použít následující indexy [9]
T LSL USL T C pm min , 3 3
(33)
Pro výpočet indexu CC pk Kane [9] uvádí vzorec CC pk min CPU , CPL
(34)
kde
CPU
USL T 3
T 1 USL T
(35)
CPL
T T LSL 1 3 T LSL
(36)
Při výpočtu způsobilosti vycházíme z dat získaných z jednotlivých výběrů, nikoliv ze všech existujících hodnot, a proto pro vymezení neznámé hodnoty indexu způsobilosti konstruujeme konfidenční meze. Tyto meze závisejí na počtu hodnot vstupujících do výpočtu, přičemž platí, že se vzrůstajícím počtem použitých hodnot se konfidenční meze zužují. K výpočtu se používají vzorce [15]
26
USL LSL 6s
min CPU , CPL
2.4.1.2.
2 k n 1 2
k n 1 2 k n 1 2
k n 1
< Cp <
USL LSL 6s
2 k n 1 1
2
k n 1
< C pk < min CPU , CPL
(37)
2 k n 1 1
2
k n 1
(38)
Indexy výkonnosti
Indexy výkonnosti slouží ke zhodnocení procesu v určitém čase. Rozdíl ve výpočtu indexů výkonnosti oproti indexům způsobilosti spočívá ve způsobu odhadu směrodatné odchylky procesu [6].
( X i X )2 k * n 1 i 1 k *n
ˆ tot stot
(39)
Pro index Pp platí [6]
Pp
USL LSL 6 * ˆ tot
(40)
Pro výpočet indexu Ppk je třeba znát hodnoty CPU a CPL .
CPU
USL ˆ 3ˆ tot
(41)
CPL
ˆ LSL 3ˆ tot
(42)
Ppk min CPU , CPL
(43)
Výpočty indexů výkonnosti v případě jednostranné specifikace [6]
Ppk
USL ˆ 3ˆ tot
(44)
Ppk
ˆ LSL 3ˆ tot
(45)
27
Výpočty indexů výkonnosti v případě asymetrických mezí [6] T LSL USL T Ppm min , ˆ 3 3ˆ tot tot
(46)
Pro výpočet indexu PPpk je třeba znát hodnoty CPU a CPL . T X 1 USL T
(47)
T X T LSL CPL 1 3ˆ tot T LSL
(48)
PPpk min CPU , CPL
(49)
USL T CPU 3ˆ tot
Konfidenční meze pro indexy výkonnosti jsou dány vzorci [15]
USL LSL 6ˆ tot
min CPU , CPL
2 k n 1 2
k n 1 2 k n 1 2
k n 1
< Pp <
USL LSL 6ˆ tot
2 k n 1 1
2
k n 1
< Ppk < min CPU , CPL
(50)
2 k n 1 1
2
k n 1
(51)
2.4.2. Výpočet indexŧ zpŧsobilosti a výkonnosti při nenormálním rozdělení Neplňuje-li proces předpoklad normality, je vhodnější zvolit jiný přístup. První možností je transformovat data pomocí Boxovy-Coxovy transformace a a postupy aplikovat na transformovaná data. Druhou možností je využít speciální postup, tzv. Clementsovu metodu. 2.4.2.1.
Box-Coxova transformace
Transformaci dat obecně použijeme v okamžiku, zjistí-li se během analýzy dat, že rozdělení výběru neodpovídá rozdělení normálnímu. Boxova-Coxova transformace zlepší symetrii rozdělení a stabilizuje rozptyl. Boxova-Coxova transformace veličiny y se provede pomocí vzorce
28
Y ´ (Y 2 ) 1
(52)
kde parametry 1 a 2 se odhadují pomocí věrohodnostní funkce [13] numericky s využitím počítače. Parametr 2 je často roven nule. 2.4.2.2.
Clementsova metoda
Clements [2] navrhnul metodu výpočtu indexu C p a C pk v případě jiného než normálního rozdělení sledované veličiny. K transformaci se používá systém Pearsonových křivek. Vzhledem k tomu, že se využívají všechny hodnoty bez ohledu na podskupiny, počítá se naopak jen index výkonnosti. Indexy Pp a Ppk mají následující podobu
Pp
USL LSL U p Lp
(53)
Ppk min( CPU , CPL)
(54)
CPU
USL M Up M
(55)
CPL
M LSL M Lp
(56)
Up je 99,865% kvantil a Lp je 0,135% kvantil. M je medián. Abychom určili potřebné kvantily, musíme odhadnout šikmost a špičatost. Šikmost se určí pomocí vzorce [13] n
g1
n ( xi x ) 3 i 1
(n 1)( n 2) s 3
(57)
špičatost odhadneme pomocí vzorce [13] n
g2
n(n 1) ( xi x ) 4 i 1
(n 1)( n 2)( n 3) s 4
29
3(n 1) 2 (n 2)( n 3)
(58)
Na základě vypočtené šikmosti a špičatosti se v tabulkách převzaté z [2] naleznou hodnoty Lp´ , Up´ a M ´ následujícím postupem: Nabývá-li šikmost kladných hodnot, hledáme tabulkovou hodnotu Lp´ v tabulce 1a. Nabývá-li šikmost záporných hodnot, hledáme tabulkovou hodnotu
Lp´ v tabulce 1b. Nabývá-li šikmost kladných hodnot, hledáme tabulkovou hodnotu Up´ v tabulce 1b. Nabývá-li šikmost záporných hodnot, hledáme tabulkovou hodnotu
Up´ v tabulce 1a. Tabulkovou hodnotu M ´ hledáme v tabulce 2. Nabývá-li šikmost kladných hodnot, otočíme znaménko výsledné hodnoty. Nabývá-li šikmost záporných hodnot, znaménko neměníme. Potřebné kvantity Lp , Up a M určíme pomocí vzorců [2] Up X s Up´
(59)
Lp X s Lp´
(60)
M X s M´
(61)
a hodnoty pak dosadíme do vzorců (53) a (54).
30
3. Systém měření „Aby bylo možné efektivně řídit variabilitu libovolného procesu, je třeba proces měření poznat a ovládat.“19 Každému zhodnocení výrobního procesu musí předcházet počáteční sběr dat. Teprve na základě získaných hodnot lze provést analýzu, ověřit stabilitu či následně testovat způsobilost. I když při výpočtu postupujeme svědomitě, nemůžeme si být nikdy jistí, zdali výsledek koresponduje se skutečným stavem procesu. Z tohoto důvodu bychom na samém počátku hodnocení měli věnovat pozornost kvalitě systému měření. Cílem měření by mělo být snížit variabilitu naměřených hodnot. Celkovou variabilitu naměřených hodnot, skládající se z variability procesu a variability systému měření je možné vyjádřit pomocí vzorce [19], 2 2 T2 procesu mereni
(66)
kde pro variabilitu systému měření platí [12] 2 2 2 mereni opakovate ln ost reprodukovate ln ost
(67)
Během analýzy systému měření hodnotíme variabilitu systému měření a polohu vůči skutečné hodnotě. V případě variability systému měření hovoříme o tzv. R&R studii. Název studie je odvozen z počátečních písmen anglických slov
repeatability a
reproducibility. Ty se překládají jako opakovatelnost resp. reprodukovatelnost. Opakovatelnost je „variabilita měření získaných jedním měřicím přístrojem, pokud jej osoba provádějící hodnocení použije několikrát při měření identického znaku u téhož dílu.“20 Reprodukovatelnost je „variabilita průměrů měření provedených různými operátory při použití stejného měřidla při měření znaku u jednoho dílu.“ 21
Measurement Systems Analysis [online]. [cit. 20.9.2010]. Dostupné z: www.jakost.cz/metodiky/msa3/frame.htm 20 Measurement Systems Analysis [online]. [cit. 20.9.2010]. Dostupné z: www.jakost.cz/metodiky/msa3/frame.htm 21 Measurement Systems Analysis [online]. [cit. 20.9.2010]. Dostupné z: www.jakost.cz/metodiky/msa3/frame.htm 19
31
V případě polohy vůči skutečné hodnotě hovoříme o studii linearity a strannosti. Strannost, často označována jako přesnost, znázorňuje rozdíl mezi referenční hodnotou a vypočteným průměrem. Strannost ukazuje velikost systematické chyby systému měření. Linearita je „rozdíl strannosti v očekávaném měřícím rozsahu zařízení.“22
Cílem studie opakovatelnosti a reprodukovatelnosti je odhalit podíl velikosti variability systému měření na variabilitě celkové. Rozdělení celkové variability je patrné v obr. 6.
Zdroj: Upraveno dle23 Obr. 6 Rozdělení celkové variability
Measurement Systems Analysis [online]. [cit. 20.9.2010]. Dostupné z: www.jakost.cz/metodiky/msa3/frame.htm 23 MSA-Analýza systému měření [online]. [cit. 20.9.2010]. Dostupné z: www.statspol.cz/request/request2006/prezentace/bednar.pdf 22
32
Využívá se následujících charakteristik [19]. %R&R =
P /T
mereni 100% T
6 mereni 100% USL LSL
(68)
(69)
Výsledek můžeme interpretovat třemi způsoby:
Je-li R&R% či P/T < 10% lze systém měření považovat za způsobilý s ohledem na variabilitu procesu či uvedené toleranční meze.
Je-li 10% < R&R% či P/T < 30% lze systém měření považovat za podmíněně způsobilý s ohledem na proces či uvedené toleranční meze. V daném případě je nezbytné zvážit náklady spojené s uvedením měřidel do požadovaného stavu či zhodnotit význam monitorovaného znaku.
Je-li R&R% či P/T > 30%
považuje se systém měření za
nezpůsobilý s ohledem na proces či uvedené toleranční meze
33
4. Představení společnosti Společnost ŠKODA AUTO a.s. je považována za jednu z nejvýznamnějších společností v České Republice, která je spjata s více než stoletou tradicí ve výrobě automobilů, s podnikatelskou činností zaměřenou na vývoj, výrobu a prodej osobních vozů, komponentů, originálních dílů a příslušenství značky Škoda a poskytování servisních služeb. Značka Škoda je zároveň považována za jednu z nejstarších automobilových značek na světě. Od roku 2007 nastává zlom v historii firmy, kdy se jediným akcionářem stává společnost Volkswagen AG. Vedle mateřského závodu ŠKODA AUTO a.s. skupinu tvoří i její plně konsolidované dceřiné společnosti ŠKODA AUTO Deutschland GmbH, ŠKODA AUTO Slovensko, s.r.o., SKODA AUTO Polska S.A., SKODA AUTO India Private Ltd. a přidružený podnik OOO Volkswagen Rus. V současné době se skupina ŠKODA AUTO skládá z montážních závodů v České Republice, Slovensku, Ukrajině, Bosně a Hercegovině, Rusku, Kazachstánu, Indii a Číně.
Zdroj: Převzato z24 Obr. 6 Výrobní a montážní závody vozů značky Škoda 24
Škoda Auto - Výroční zpráva 2008 [online]. [cit. 22.8.2010]. Dostupné z: http://www.skodaauto.cz/company/CZE/Documents/Pro_investory/Vyrocni_zpravy/SkodaAuto_AnnualReport_200 8_CZ.pdf , s.56 34
Od okamžiku vstupu společnosti ŠKODA AUTO a.s. do koncernu Volkswagen sílí image značky, firma navyšuje výrobu, rozšiřuje produktové portfolio a buduje silné postavení v rámci celého světa. V současné době skupina Škoda Auto působí na více než 100 trzích, na kterých bylo v roce 2009 zákazníky poptáváno celkem 684 226 vozů značky Škoda.
Zdroj: Převzato z25 Obr. 7 Prodejní regiony skupiny ŠKODA AUTO V roce 2010 zahrnuje produktové portfolio značky pět modelových řad, které jsou vyvíjeny dle zásady „Simply Clever“. Tyto produkty nabízejí rozumný poměr hodnoty vozu a jeho ceny, a to vše při dodržení vysokého standardu kvality.
25
Škoda Auto - Výroční zpráva 2009 [online]. [cit. 22.8.2010]. Dostupné z: http://www.skodaauto.cz/company/CZE/Documents/Pro_investory/Vyrocni_zpravy/SkodaAuto_AnnualReport_200 9_CZ.pdf, s.41 35
5. Kapoty Popis problému Kapota je složena ze dvou dílů, které jsou k sobě fixovány pomocí vrstvy lepidla. Nanesená vrstva lepidla má zásadní vliv na schopnost výrobku vyhovět požadavkům na kvalitu. V případě nánosu velkého množství lepidla dochází k proznačení slepených ploch do svrchní části kapoty, což je neakceptovatelné ze strany konečného zákazníka. V případě nánosu příliš malého množství hrozí naopak odtržení svrchní části od spodní. Z důvodů vysokých nákladů v případě produkce zmetků bylo rozhodnuto o nezbytnosti regulace procesu nanášení vrstvy lepidla.
Zdroj: Firemní materiály Obr. 8 Kapota Vrstva lepidla je nanášena na více míst najednou. Na základě odlišných vrstev je třeba sledovat každý bod samostatně. V práci se budeme věnovat pouze vybraným bodům č. 13 a č. 27, znázorněných v obr. 8.
36
Zdroj: Firemní materiály Obr. 9 Nanesená vrstva lepidla Technické specifikace Technické
specifikace
jsou
definovány
pomocí
tolerančního
1,5 – 3,5 mm s cílovou hodnotou 2 mm.
Zdroj: Firemní materiály Obr. 10 Řez kapotou se zakreslenými technickými mezemi 37
pole
Analýza procesu lepení Regulovanou veličinou je vrstva naneseného lepidla. Logické podskupiny zahrnují výrobky vyrobené během jednoho pracovního dne. Důvodem je použití nového barelu nanášeného lepidla a seřízení výrobního zařízení vždy na počátku ranní směny. Rozsah logických podskupin n odpovídá třem po sobě jdoucím jednotkám, odebraným po pěti hodinách od zahájení výroby denní produkce. Zmíněný časový odstup ponechává prostor pro ustálení výrobního procesu a předpokládá, že jednotky vznikají za stejných či podobných podmínek.
5.1. Analýza procesu v bodě č. 13 Získaná data zachycují období březen až červenec roku 2010. Celkem bylo vytvořeno 94 podskupin o rozsahu 3. Histogram na obr. 11 znázorňuje rozdělení naměřených hodnot. V obrázku je zakreslena cílová hodnota 2 mm a dolní a horní toleranční mez (1,5 mm a 3,5 mm). Dále je zobrazeno rozpětí odpovídající přirozenému kolísání velikosti 6s, kde s je odhad sigma, znázorněné kratšími úsečkami.
Obr. 11 Histogram regulované veličiny vrstva lepidla v bodě č. 13, Statgraphics
38
Z obrázku je patrné, že tvar rozdělení neodpovídá zakreslenému modelu Gaussovy křivky. Proces je posunut od cílové hodnoty T doprava směrem k vyšším hodnotám. Pro objektivnější posouzení normality rozdělení hodnot byl použit normální pravděpodobnostní graf a Shapiro-Wilkův test, které nabízí Statgraphics.
Obr. 12 Pravděpodobnostní graf regulované veličiny vrstva lepidla v bodě č. 13, Statgraphics Protože body neleží v přímce, nelze předpokládat normální rozdělení. To potvrzuje
i Shapiro-Wilkův test. Na základě velmi malé p-hodnoty zamítáme
hypotézu o normalitě na hladině významnosti 0,001. Konstrukce regulačních diagramŧ Pro statistickou regulaci procesu lepení byl zvolen regulační diagram pro průměr a rozpětí. Na základě vzorce (24) byla odhadnuta směrodatná odchylka procesu a střední hodnota:
ˆ 2,89397
s 0,0967689
39
Pomocí
vzorců (4), (5), (6), (7), (8) a (9) byla vypočtena centrální přímka a
regulační meze diagramu pro rozpětí (obr. 13)
UCL 0,42173
CL 0,16383
LCL 0 ,
diagramu pro průměr (obr. 14)
UCL 3,06158
CL 2,89397
LCL 2,72636 .
Obr. 13 Regulační diagram pro rozpětí regulované veličiny vrstva lepidla v bodě č. 13, Statgraphics
40
Obr. 14 Regulační diagram pro průměr regulované veličiny vrstva lepidla v bodě č. 13, Statgraphics V diagramu pro rozpětí se vyskytují dva body odpovídající podskupinám č. 20 a 82 nad horní regulační mezí. Příčinu vychýlení se však ani po konzultaci s dodavatelem nepodařilo zjistit. V diagramu pro průměr je patrný především trend, tj. systematický růst střední hodnoty procesu, a to na počátku sledovaného období, tedy ve 3. – 4. měsíci roku 2010, a na konci sledovaného období, tedy v 7. měsíci roku 2010. Na základě výše popsaných skutečností proces nelze považovat za stabilní. Výpočet indexů způsobilosti nemá smysl. Na základě vyhodnocených dat lze současnou situaci považovat za kritickou. Je vhodný okamžitý zásah do procesu v podobě odhalení řady vymezitelných příčin působících na proces.
41
5.2. Analýza procesu v bodě č. 27 Použitá data zachycují období březen až červen roku 2010. Celkem bylo vytvořeno 67 podskupin o rozsahu 3. Histogram na obr. 15 znázorňuje rozdělení naměřených hodnot. V obrázku je zakreslena cílová hodnota 2 mm a dolní a horní toleranční mez (1,5 mm a 3,5 mm). Dále je zachyceno rozpětí odpovídající přirozenému kolísání velikosti 6s, znázorněné kratšími úsečkami.
Obr. 15 Histogram regulované veličiny vrstva lepidla v bodě č. 27, Statgraphics Z obrázku je patrné, že tvar rozdělení odpovídá zakreslenému modelu Gaussovy křivky. Proces je posunut od cílové hodnoty T doprava směrem k vyšším hodnotám. Následně je sestrojen pravděpodobnostní graf a vypočítán Shapiro-Wilkův test.
42
Obr. 16 Pravděpodobnostní graf regulované veličiny vrstva lepidla v bodě č. 27, Statgraphics Protože body leží v přímce, lze předpokládat normální rozdělení. To potvrzuje i Shapiro-Wilkův test. Na základě vysoké p-hodnoty rovné 0,963828 nezamítáme hypotézu o normalitě na hladině významnosti 0,05. Konstrukce regulačních diagramŧ Na základě vzorce (24) byla odhadnuta směrodatná odchylka procesu a střední hodnota:
ˆ 2,68245
s 0,148186
Výpočet centrální přímky a regulačních mezí je podobný jako v případě bodu č. 13. Postupuje se dle stejných vzorců. Diagram pro rozpětí (obr. 17)
UCL 0,644336
CL 0,249629
LCL 0
Diagram pro průměr (obr. 18)
UCL 2,93976
CL 2,68245
43
LCL 2,42514
Obr. 17 Regulační diagram pro rozpětí regulované veličiny vrstva lepidla v bodě č. 27, Statgraphics
Obr. 18 Regulační diagram pro průměr regulované veličiny vrstva lepidla v bodě č. 27, Statgraphics V obou diagramech se veškeré hodnoty nacházejí uvnitř regulační mezí. Proces lze považovat za stabilní. Následně lze přistoupit k výpočtu indexů způsobilosti.
44
Výpočet indexŧ zpŧsobilosti a výkonnosti Mezi nejpoužívanější indexy způsobilosti patří index C p , označován jako krátkodobý index způsobilosti. Vedle indexu C p počítáme i index výkonnosti PP . Pomocí vzorců (29), (30), (40) a (43) byly vypočteny indexy C p 2,24943
PP 2,23716
C pk 1,83902
PPk 1,82899
Na základě vypočítaných indexů C p a PP je zřejmé, že proces je způsobilý. Indexy C pk a Ppk udávají polohu vůči mezím daným specifikací. Z důvodů jejich rozdílných hodnot je patrné, že proces je posunut vůči středu tolerančního pole směrem k vyšším hodnotám. Aby bylo možné posoudit polohu procesu vůči cílové hodnotě T, použijeme vzorce (33), (34), (46) a (49). C pm 1,124713
Ppm 1,118583
CC pk -0,410408
PPpk -0,408171
Na základě vypočítaných indexů C p a PP ≥ 1 lze říci, že výrobní proces lze považovat za podmíněně způsobilý. Důvodem záporných hodnot CC pk a PPpk je skutečnost, že střední hodnota je posunuta značně vpravo od cílové hodnoty. Uvedené hodnoty indexů jsou pouze bodové odhady získané na základě naměřených hodnot. Abychom vymezili interval, v němž se může nacházet index popisující celý proces, konstruujeme konfidenční meze. K výpočtu byly použity vzorce (37), (38), (50) a (51). 2,029 < C p < 2,470 1,653 < C pk < 2,026 2,017 < Pp < 2,457 1,644 < Ppk < 2,015
45
V případě bodu č. 27 lze situaci považovat za uspokojivou. Proces je stabilní a považován za podmíněně způsobilý. Navrhuji však nespokojit se s danou situací a usilovat o přiblížení se k cílové hodnotě vrstvy lepidla 2 mm, což by vedlo ke snížení pořizovacích nákladů na lepidlo.
46
6. Přitahovač Popis problému Do modelu Škoda Super je montován zámek pátých dveří, tzv. přitahovač. Jedná se o dosti složitý díl, neboť pro jeho správnou funkci jsou důležité hned tři veličiny – síla, jenž musí být dostatečně velká, aby dokázala „docvaknout“ víko pátých dveří, zdvih, který musí urazit jistící trn, aby zamezil víku jeho opětovnému otevření a čas nezbytný pro uživatelský komfort. Teprve za předpokladu souladu všech tří veličin lze díl pokládat za funkční.
Zdroj: Firemní materiály Obr. 19 Přitahovač víka pátých dveří
47
Zdroj: Firemní materiály Obr. 20 Umístění přitahovače ve voze Technické specifikace Technické specifikace jsou definovány pomocí následujících tolerančních polí. Velikost síly, definované pomocí dolní přípustné meze, při které proces pracuje s požadovanou přesností, odpovídá velikosti 1000 N. Dráha zdvihu je definována pomocí tolerančního pole 8,5 mm ± 0,5 mm s požadovanou hodnotou ve středu intervalu. Časový interval, odpovídající požadavkům zákazníka, je stanoven pomocí tolerančního pole 2,8 s ± 0,3 s požadovanou hodnotou ve středu intervalu.
48
Analýza procesu Regulovanou veličinou jsou síla, zdvih a čas. Použitá data odpovídají 1357 jednotkám vyrobených během jednoho pracovního dne. Jelikož se měří všechny vyráběné jednotky, použije se diagram pro individuální hodnoty.
6.1. Analýza procesu s regulovanou veličinou síla Histogram na obr. 21 znázorňuje rozdělení naměřených hodnot. V obrázku je zakreslena dolní specifická mez 1000 N společně s dolní přirozenou mezí 1060 N, znázorněné pomocí kratší úsečky.
Obr. 21 Histogram regulované veličiny síla, Statgraphics Z obrázku je patrné, že tvar rozdělení odpovídá zakreslenému modelu Gaussovy křivky. Patrný je malý rozptyl naměřených hodnot. Stanovená toleranční mez není přesažena. Pro posouzení normality jsou použity stejné nástroje statistické regulace jako v případě kapot.
49
Obr. 22 Pravděpodobnostní graf regulované veličiny síla, Statgraphics Protože body leží v přímce, lze předpokládat normální rozdělení. To potvrzuje i Shapiro-Wilkův test. Na základě p-hodnoty rovné 0,186517 nezamítáme hypotézu o normalitě na hladině významnosti 0,05. Konstrukce regulačního diagramu Pro statistickou regulaci procesu s regulovanou veličinou síla byl zvolen regulační diagram pro individuální hodnoty. Na základě vzorce (24) byla odhadnuta směrodatná odchylka procesu a střední hodnota:
ˆ 1081,79
s 7,53938
Pomocí vzorců (21), (22) a (23) byla vypočtena centrální přímka a regulační meze diagramu pro individuální hodnoty (obr. 23).
UCL 1104,41
CL 1081,79
50
LCL 1059,17
Obr. 23 Regulační diagram pro individuální hodnoty regulované veličiny síla, Statgraphics V diagramu se vyskytují tři body vně regulačních mezí. Jelikož naměřená data hodnotíme zpětně, je nyní patrné, že se jedná o zbytečný signál. K žádnému zásahu do procesu nedošlo. Proces se sám vrátil do původního stavu. Proces je považován za stabilní. Následně lze přistoupit k výpočtu indexů způsobilosti. Výpočet ukazatelŧ zpŧsobilosti a výkonnosti Pomocí vzorců (30) a (43) byly vypočteny indexy C pk 3,61606
PPk 3,69284
Na základě odlišných hodnot obou indexů je zřejmé, že skutečný střed procesu je posunut směrem k vyšším hodnotám. Následnou konstrukcí konfidenčních mezí dle vzorců (38) a (51) určíme interval, ve kterém se nachází index popisující celý proces. 3,47881 < C pk < 3,7533 3,55273 < Ppk < 3,83295
51
6.2. Analýza procesu s regulovanou veličinou zdvih Histogram na obr. 24 znázorňuje rozdělení naměřených hodnot. V obrázku je zakreslena cílová hodnota 8,5 mm a dolní a horní toleranční mez (8 mm a 9 mm). Dále je zachyceno rozpětí odpovídající přirozenému kolísání velikosti 6s, znázorněné kratšími úsečkami.
Obr. 24 Histogram regulované veličiny zdvih, Statgraphics Proces je měřen nedostatečně přesně. Z obrázku je patrné, že tvar rozdělení neodpovídá zakreslenému modelu Gaussovy křivky. Výrobní proces je téměř centrovaný. Stanovené toleranční meze nepřesahuje. Nyní je sestrojen pravděpodobnostní graf a vypočítán Shapiro-Wilkův test.
52
Obr. 25 Pravděpodobnostní graf regulované veličiny zdvih, Statgraphics Protože body neleží v přímce, nelze předpokládat normální rozdělení. To potvrzuje i Shapiro-Wilkův test. Na základě velmi malé p-hodnoty zamítáme hypotézu o normalitě na hladině významnosti 0,001. Výsledkem je pouze osm různých obměn sledované veličiny. Měření bylo prováděno s přesností na jedno desetinné číslo. Přesnost měření je nedostatečná. Jelikož jak z histogramu tak pravděpodobnostního grafu je patrná nenormalita a nedostatečná přesnost měření, nemá smysl pokračovat regulací a indexy způsobilosti. Má-li mít regulace v tomto případě smysl, je nezbytné, aby přesnost měření byla charakterizována maximální hodnotou 0,02. O dané problematice pojednává kapitola 3. Výrobci doporučuji, aby se zaměřil na způsobilost používaných měřidel.
6.3. Analýza procesu s regulovanou veličinou čas Histogram na obr. 26 znázorňuje rozdělení naměřených hodnot. V obrázku je zakreslena cílová hodnota 2,8 s a dolní a horní toleranční mez (2,5 s a 3,1 s). Dále je zakresleno rozpětí odpovídající přirozenému kolísání velikosti 6s, znázorněné kratšími úsečkami.
53
Obr. 26 Histogram regulované veličiny čas, Statgraphics Z obrázku je patrné, že tvar rozdělení neodpovídá zakreslenému modelu Gaussovy křivky. Sloupce s vyššími hodnotami klesají pozvolněji. Proces je posunut od cílové hodnoty T doleva směrem k nižším hodnotám. Za pomoci pravděpodobnostního grafu a Shapiro-Wilkova testu je ověřena normalita.
Obr. 27 Pravděpodobnostní graf regulované veličiny čas, Statgraphics
54
Protože body neleží v přímce, nelze předpokládat normální rozdělení. To potvrzuje i Shapiro-Wilkův test. Na základě velmi malé p-hodnoty zamítáme hypotézu o normalitě na hladině významnosti 0,001. Konstrukce regulačních diagramŧ Pro statistickou regulaci procesu s regulovanou veličinou čas byl zvolen regulační diagram pro individuální hodnoty. Pomocí vzorce (24) byla odhadnuta směrodatná odchylka procesu a střední hodnota:
ˆ 2,07892
s 0,21104
Centrální přímka a regulační meze v obr. 28 jsou vypočteny dle stejných vzorců jako v případě regulované veličiny síla.
UCL 2,71204
CL 2,07892
LCL 1,4458
Obr. 28 Regulační diagram pro individuální hodnoty regulované veličiny čas, Statgraphics Na základě regulačního diagramu je patrné, že značný počet bodů spadá vně regulačních mezí. Proces není považován za stabilní.
55
Výpočet indexŧ zpŧsobilosti a výkonnosti při nenormálním rozdělení Nyní je bráno v úvahu nenormální rozdělení, kdy je vhodnější zvolit jiný přístup. První možností je transformovat data a použít regulaci procesu a indexy způsobilosti. Druhou možností je využít speciální postup při určování způsobilosti procesu, tzv. Clementsovu metodu. Box–Coxova transformace V případě Box-Coxovy transformace byla data transformována dle vzorce (52). Program sám nalezl optimální hodnotu parametru lambda.. Výsledkem je regulační diagram znázorněný na obr. 29.
Obr. 29 Regulační diagram pro individuální hodnoty regulované veličiny čas v případě Box-Coxovy transformace, Statgraphics Závěry týkající se překračování mezí jsou podobné. Proces není pod kontrolou. V případě ověření způsobilosti procesu byly nejprve aplikovány indexy na transformovaná data. Následně byl použit vzorec (24) k odhadu směrodatné odchylky procesu a střední hodnoty:
ˆ 0,504091
s 0,0498059
56
Pomocí vzorců (29), (30), (40) a (43)) byly vypočteny indexy C p 0,268935
PP 0,258996
C pk -0,593061
PPk -0,571143
Na základě uvedených indexů C p a PP je zřejmé, že proces je nezpůsobilý. Odlišné hodnoty indexů C pk a PPk vyjadřují, že skutečný střed procesu je posunut směrem k nižším hodnotám. Clementsova metoda K výpočtu Clementsovy metody bylo využito postupu blíže popsaného v teoretické části práce. Při výpočtu indexů způsobilosti a výkonnosti byly využity vzorce (53) a (54). PP 0,44099
PPk 1,05453
Na základě vypočítaného indexu PP je zřejmé, že proces je nezpůsobilý. Index PPk udává polohu vůči mezím daným specifikací. Skutečný střed procesu je posunut směrem k nižším hodnotám. Jak v případě Box-Coxovy transformace tak Clementsovy metody jsme dospěli k závěru, že proces je nezpůsobilý. Je patrné, že v případě regulované veličiny čas není proces zcela pod kontrolou. Doporučuji zásah do procesu za účelem odhalení vymezitelných příčin, které způsobují vzniklou variability.
57
7. Závěr Cílem diplomové
práce
bylo ukázat
možnosti
aplikace vybraných
statistických nástrojů, v podobě regulačních diagramů a indexů způsobilosti, při zajišťování kvality nakupovaných dílů ve společnosti ŠKODA AUTO a.s.. Pozornost jsem věnoval dodávaným dílům, u nichž se v minulosti objevily problémy v podobě neshodných výrobků, stejně tak jako dílům složitým, u nichž není jednoduché udržet proces v požadovaném stavu. Potřebná data, zachycující hodnoty sledované regulované veličiny, jsem získal od jednotlivých dodavatelů, stejně tak jako znalosti o průběhu výrobních procesů. Následně jsem aplikoval nástroje statistické regulace procesu, ověřující, zdali je proces stabilní či nikoliv. Nejen v případě prokázání stability jsem pokračoval zhodnocením způsobilosti procesu. V případě, že se prokázala nestabilita či nezpůsobilost procesu, navrhl jsem plán možného řešení. V průběhu tvorby diplomové práce byly zkoumány dva nakupované díly dodávané do společnosti ŠKODA AUTO a.s.. První případová studie se věnovala nanášení vrstvy lepidla mezi spodní a svrchní částí kapoty. Analyzovaná data zachycovala vrstvu nanášeného lepidla v bodě č. 13 a 27. Na základě vyhodnocených dat v bodě č. 13 byl zjištěn především patrný trend růstu střední hodnoty. Proces tedy nelze považovat za stabilní. K docílení stability procesu by bylo vhodné hledání příčiny. V případě bodu č. 27 lze situaci považovat za uspokojivou. Proces je stabilní a považován za podmíněně způsobilý. Navrhl jsem však nespokojit se s danou situací a usilovat o přiblížení se k cílové hodnotě vrstvy lepidla 2mm, což by vedlo ke snížení pořizovacích nákladů na lepidlo. V druhé případové studii jsem se zaměřil na díl, u něhož jsou sledovány hned tři regulované veličiny najednou (síla, zdvih a čas). V případě regulované veličiny síla je situace v pořádku. Proces je stabilní a způsobilý. Analyzovaná data zachycují regulovanou veličinu zdvih vypovídají o nedostatečné přesnosti měření. Výsledkem je pouze osm různých obměn sledované veličiny. Má-li mít regulace v tomto případě smysl, je nezbytné zabývat se způsobilostí měřidel.
58
V případě poslední sledované veličiny čas bylo zjištěno nenormální rozdělení sledované veličiny. Na základě značného počtu bodů vně regulačních mezí není proces považován za stabilní. Jelikož bylo bráno v úvahu nenormální rozdělení, byla ověřena způsobilost procesu využitím Box-Coxovy transformace a Clementsovy metody. Ani v případě Box-Coxovy transformace tak Clementsovy metody se situace nemění. Proces je považován za nezpůsobilý. Z důvodů značné variability doporučuji zásah do procesu. Na základě případových studií byly odhaleny jisté nedostatky výrobních procesů u zkoumaných dílů dodávaných do společnosti ŠKODA AUTO a.s.. Uvedené návrhy opatření by však měli pomoci k jejich odstranění. Využití statistické regulace procesu lze tedy považovat za klíčový nástroj k dosažení vysoké kvality výrobku.
59
Seznam literatury [1]
BISSELL, D. Statistical Methods for SPC and TQM. 1. vyd. London: Chapman & Hall, 1994, 373 s. ISBN 0-412-39440-5.
[2]
CLEMENTS, J. A. Process Capability Calculations for Non-Normal Distributions. Qualityware. Curitiba: 1989, september, s. 95-100.
[3]
ČSN
ISO
8258
Shewhartovy
regulační
diagramy.
Praha:
Český
normalizační institut, 1998, 36 s. [4]
ČSN ISO 7880 Regulační diagramy – všeobecné pokyny a úvod. Praha: Český normalizační institut, 1995, 12 s.
[5]
ČSN ISO 690 Bibliografické citace – obsah, forma a struktura. Praha: Český normalizační institut, 1996, 32 s.
[6]
HŮLOVÁ, M. JAROŠOVÁ, E. Statistické metody v managementu kvality, enviromentu a bezpečnosti. Praha: Oeconomica, 2004, 119 s. ISBN 80-245-0691-2.
[7]
Firemní materiály ŠKODA AUTO a.s.
[8]
JOHN, P. W. M. Statistical Methods in Engineering and Quality Assuarance. The United States of America: A Wiley-Interscience Publication, 1990, 373 s. ISBN 0-471-82986-2.
[9]
KANE, V. E. Process Capability Indices.Journay of Duality Technology. Milwaukee: 1986, roč. 18., č. 1., s. 41-52.
[10]
KŘEPELA, J. HORÁLEK, V. FABIAN, F. Statistické metody řízení jakosti. Praha: ČSJ, 2007, 390 s. ISBN 978-80-02-01897-1.
[11]
-, Malý metodický průvodce. 1. vyd. Mladá Boleslav: Škoda, automobilová, 2005
[12]
RYAN, T. P. Statistical Methods for Quality Improvement. 2. vyd. New York: A Wiley-Interscience Publication, 2000, 555 s. ISBN 0-471-19775-0.
[13]
Statgraphics Centurion XVI, User manual, StatPoint Technologies, Inc. 2009.
[14]
-, Školení ke kvalitě. 1. vyd. Mladá Boleslav: Škoda, automobilová, 2009
60
[15]
TEREK, M. HRNČIAROVÁ, L. Štatistické riadenie kvality. Bratislava: Iura Edition, 2004, 234 s. ISBN 80-89047-97-1.
Internetové zdroje [16]
Aplikace SPC metodiky pro oblast předvýroby a konečné montáže pro projekt X7. [cit. 17.11.2010]. Dostupné z: http://dspace.khihovna.utb.cz/handle/10563/6311
[17]
Measurement Systems Analysis [online]. [cit. 20.9.2010]. Dostupné z: www.jakost.cz/metodiky/msa3/frame.htm
[18]
Mocninná a Box-Coxova transformace v průzkumové exploratorní analýze dat
[online].
[cit.
11.10.2010].
Dostupné
z:
http://meloun.upce.cz/docs/publication/059.pdf [19]
MSA-Analýza systému měření [online]. [cit. 20.9.2010]. Dostupné z: www.statspol.cz/request/request2006/prezentace/bednar.pdf
[20]
Ověřování předpokladu normality [online]. [cit. 17.11.2010]. Dostupné z: http://isq.cz/npj/12 - Testy_normality.pdf
[21]
Statistická regulace procesu [online]. [cit. 20.9.2010]. Dostupné z: http://www.csq.cz/cs/nabidka-vydavatelstvi/prirucky-qs-9000/spc-statistickaregulace-vyrobniho-procesu.html
[22]
Škoda Auto - Výroční zpráva 2009 [online]. [cit. 22.8.2010]. Dostupné z: http://www.skodaauto.cz/company/CZE/Documents/Pro_investory/Vyrocni_zpravy/SkodaAut o_AnnualReport_2009_CZ.pdf
[23]
Škoda Auto - Výroční zpráva 2008 [online]. [cit. 22.8.2010]. Dostupné z: http://www.skodaauto.cz/company/CZE/Documents/Pro_investory/Vyrocni_zpravy/SkodaAut o_AnnualReport_2008_CZ.pdf
61
Seznam příloh Příloha č. 1 Součinitele pro výpočet přímek regulačních diagramů
62