Определение комплексной оценки качества продукции

1. Собрать информацию о качестве данного вида продукции, выпускаемой на конкретном предприятии:

- объем выпускаемой продукции;

- виды и количество брака.

Выполнить пооперационный анализ с указанием видов дефектов и брака.

2. Собрать исходную информацию о фактических абсолютных значениях единичных показателей качества данного вида продукции (провести измерения, воспользоваться сертификатами качества и др.), оценить ее на соответствие требованиям НД и ТД.

3. Выбрать метод определения весомостей единичных показателей, определить коэффициенты весомости и их статус.

4. Определить базовые, номинальные и предельные значения по каждому показателю качества, выбрать формулу для определения оценки и произвести оценку качества продукции по единичным показателям.

Цель оценки состоит в определении степени соответствия фактических значений и показателей, характеризующих качество изготовленной продукции, установленным требованиям нормативной и конструкторской документации, стандартам, техническим условиям и другой документации.

5. Определить технический уровень качества заданного вида металлургической продукции по отношению к лучшему аналогу.

Аналогами являются изделия отечественного и зарубежного производства того же вида, что и сравниваемое изделие, обладающее общностью функционального назначения и условий применения. Если показатели аналогов ниже показателей качества, установленных соответствующими стандартами, то данные этих стандартов принимаются за показатели аналога.

6. Определить комплексную оценку качества данной продукции.

1.5.Определение причин возникновения дефектов
и рекомендации по повышению качества

При определении факторов, влияющих на какой-либо показатель качества, часто применяют причинно-следственную диаграмму (диаграмму Исикавы). Она позволяет выявить причины возникновения дефектов и сосредоточиться на устранении этих причин. При этом анализируют пять основных причинных фактора: человек, машина (оборудование), материал, метод работ (технология) и измерение. При анализе этих факторов выявляют вторичные, а может быть, и третичные причины, приводящие к дефектам и подлежащие устранению. Поэтому для анализа дефектов и построения диаграммы необходимо определить максимальное число причин, которые могут иметь отношение к допущенным дефектам.

Часто диаграмму Исикавы строят одновременно со столбиковыми диаграммами Парето. Последняя позволяет наглядно представить величину потерь в зависимости от различных дефектов. Благодаря этому можно сначала сосредоточить внимание на устранении тех дефектов, которые приводят к наибольшим потерям.

ГРАФИЧЕСКАЯ ЧАСТЬ

1. Построить схему технологического процесса производства конкретного вида продукции с перечнем видов брака по операциям и (или) построить контрольную карту «дата – доля дефектных изделий».

2. Построить гистограммы распределения фактических значений замеренных параметров продукции.

3. Построить диаграммы Парето и Исикавы с целью выявления факторов и условий, влияющих на качество продукции.

3. ПРИМЕР. ОПРЕДЕЛЕНИЕ КОМПЛЕКСНОЙ ОЦЕНКИ
КАЧЕСТВА НЕРАВНОПОЛОЧНОГО ГНУТОГО ШВЕЛЛЕРА

Определить комплексную оценку и уровень качества неравнополочного гнутого швеллера 67x65x35x3,0х9000 мм, выпускаемого по ГОСТ 8281-80 в ЛПЦ-7 ОАО «Магнитогорский металлургический комбинат» (рис. 2).

 

Рис. 2. Допуски на размеры неравнополочного швеллера по ГОСТ 8281-80 для профилей высокой точности

 

В данной задаче для определения номенклатуры, статуса показателей качества гнутых профилей и их весомостей использовали метод экспертного опроса специалистов (ГОСТ 23554-79), который отличается гибкостью и наглядностью.

Построим «дерево свойств» для неравнополочного гнутого швеллера (см. рис. 1).

Для примера возьмем только геометрические показатели.

Пределы изменения геометрических параметров неравнополочного гнутого швеллера 67x65x35x3,0х9000 мм по различным стандартам представлены в табл. 1 и рис. 2.

Таблица 1

Выписка из стандартов «Швеллеры гнутые неравнополочные. Сортамент»

Параметр (показатель качества) неравнополочного гнутого швеллера 67х65х35х3,0 мм (длина 9 м)	Предельные отклонения
Россия – ГОСТ 8281-80	Болгария – БДС 13403-76	Чехия – CSN 426973 jk 154	Япония JIS G 3350 (1977)
Точность профилирования	Поштучное профилирование	Непрерывное профилирование
высокая	повышенная	обычная
1. Высота стенки (67 мм)	+0,5 -1,0	±1,5	±2,0	±1,5	±1,3	±1,0	±1,5
2. Ширина меньшей полки (35 мм)	±0,5	±1,0	±1,5	±1,5	±1,0	±0,5	±1,5
3. Ширина большей полки (65 мм)	+0,5 -1,0	±1,5	±2,5	±2,0	±1,5	±1,0	±1,5
4. Длина профиля (9000 мм)	+50	+50	+50	+50	+50	+50	+50
5. Толщина профиля (3,0 мм)	±0,21	Нет сведений	±0,28
6. Предельные отклонения от угла подгибки (90о), не более, град.	±1,0	±1,5	±1,5	±1,5	k2	±1,5
7. Волнистость, мм/1 м	<2,0	<2,0	2,0	-
8. Кривизна, %, от длины	<0,1	<2,5 мм/м	-	<0,2 %
9. Скручивание, град/м	1,0 (<10 град.)	1,0(<10 град.)	2,0	-
10. Радиус подгибки, Rв, мм	£5,0	-	-	-

Из таблицы видно, что самые жесткие требования предъявляются по ГОСТ 8281-80

к профилям высокой точности.

 

Определяем доверительный интервал фактических значений по стенке профиля.

Номинальное значение h = 67,0 мм.

Производим n замеров фактических значений (х_i) по высоте стенки профиля. Всего сделано 48 замеров. Результаты замеров сводим в табл. 2.

Таблица 2

К расчету `х и s по высоте стенки профиля

хi	n	хi n	n`(х- хi)2
65,5 66,0 66,5 67,0 67,5 68,0		65,5 132,0 1396,5 1340,0 202,5 68,0	1,5876 1,1552 1,4196 1,1520 1,6428 1,5376
N = 48	S = 8,4948

 

Строим гистограмму по замеренным значениям, на которую наносим пределы изменения h по ГОСТ 8281-80 для швеллеров нормальной, повышенной и высокой точности (рис. 3).

Определяем среднее арифметическое значение замеренных параметров (математическое ожидание) ` х и среднеквадратическое отклонение s.

; ;

где х_i - замеренные значения параметра (случайная величина);

n - объем выборки.

`х = 66,76 мм; s = 0,4251 мм.

 

 

Рис. 3. Гистограмма фактических значений замеров по высоте стенки швеллера

*) - допуск по ГОСТ 8281-80 по стенке швеллера h

для профилей нормальной точности;

**) - то же для профилей повышенной точности;

***) - то же для профилей высокой точности;

67 – номинальное значение, которое принимаем за базовое.

 

Доверительный интервал фактических значений с надежностью Р =0,99 определяем с помощью правила трех сигм: отклонение истинного значения измеряемой величины от среднего арифметического значения результатов измерений не превосходит утроенной средней квадратической ошибки этого среднего значения

х _факт= ` х ± 3 s = 66,76 ± 0,4251 × 3 = (65,48...68,04) мм

или х _факт = мм.

То есть доверительный интервал по стенке h швеллера

D_h = 2,56 мм, а Р_i^пр =Р_i^номин ± DР_max = 67 – 1,52 = 65,48 мм » 65,5 мм.

Аналогично находим доверительные интервалы фактических значений по остальным геометрическим параметрам и записываем их в таблицу 6.

 

Далее определяем весомости показателей качества.

Весомость показателей качества определяем экспертным методом. Пусть пять независимых экспертов оценивают степень важности показателей качества по пятибалльной шкале. Исходные данные заносим в таблицу 3.

Таблица 3

Исходные данные по степени важности единичных показателей

Показатель (отклонение от размеров) j =1,2,.., n	Весомость (степень важности)	D/К
Номер эксперта (е =1,2,…, m)
1. Высота стенки, h	4 D	4 D	5 D	5 D	5 D	D	Rj1= 23
2. Ширина малой полки, b 1	5 D	5 D	5 D	5 D	5 D	D	Rj2= 25
3. Ширина большей полки, b 2	4 D	5 D	5 D	5 D	4 D	D	Rj3= 23
4. Толщина профиля, S	2 D	3 D	2 D	3 D	2 D	D	Rj4= 12
5. Угол подгибки, a	5 D	4 D	4 D	5 D	5 D	D	Rj5= 23
6. Длина профиля, l	5 D	5 K	5 D	5 D	4 D	D	Rj6= 24
7. Скручивание	5 D	5 D	5 D	4 D	5 D	D	Rj7= 24
8. Кривизна	5 D	4 D	3 K	4 K	5 D	D	Rj8= 21
9. Волнистость	5 D	4 D	3 K	4 K	4 D	D	Rj9= 20
10. Радиус подгибки, R в	3 K	2 K	4 D	3 K	3 K	K	Rj10= 15

Примечание: D – доминирующие показатели, К – компенсируемые показатели. Нулевое значение любого из доминирующих показателей вызывает нулевое значение комплексного показателя, нулевое значение какого-либо единичного компенсируемого показателя не влечет снижения до нуля комплексного показателя. Эти свойства позволяют учитывать не только количественное различие в значимости единичных показателей, но и качественное или иначе – статус доминирующего или компенсируемого показателя. Такое разделение параметров качества наиболее полно отражает специфику продукции различных отраслей.

Дальнейшую обработку данных ведем по формулам:

- параметр весомости j -го показателя у е -го эксперта:

; (1)

- параметр весомости единичного j- го показателя:

, (2)

где R_je – оценка весомости j –го свойства у е – го эксперта, баллы;

m – количество экспертов;

n – количество свойств.

Весомость j -го показателя у е -го эксперта для первого значения: 4 / 43 = 0,093; весомость единичного j- го показателя для первого значения: 0,548 / 5 = 0,1096.

Проводим пересчет весомостей с учетом, что

S М_1-9^D = 1 - 0,0716 = 0,9284 – принимаем за единицу, тогда для первого доминирующего показателя справедлива пропорция:

0,1096 – х

0,9284 – 1 Отсюда х = 0,1181.

Результаты расчетов заносим в таблицу 4.

Таблица 4

Определение весомости единичных показателей

Показатель качества, j= 1 ,…,n	Mje	S Mje	Mj	Mj п
e = 1,…, m
1. h	0,093	0,090	0,182	0,116	0,119	0,548	0,1096	0,1181
2. b1	0,116	0,122	0,122	0,116	0,119	0,595	0,1190	0,1282
3. b2	0,093	0,122	0,122	0,116	0,095	0,548	0,1096	0,1181
4. S	0,047	0,073	0,049	0,070	0,048	0,287	0,0574	0,0618
5. a	0,116	0,098	0,098	0,116	0,119	0,547	0,1094	0,1178
6. l	0,116	0,122	0,122	0,116	0,095	0,571	0,1142	0,1230
7. Скручив.	0,116	0,122	0,122	0,093	0,119	0,572	0,1144	0,1232
8. Кривизна	0,116	0,098	0,073	0,093	0,119	0,499	0,0998	0,1075
9. Волнист.	0,116	0,098	0,073	0,093	0,095	0,475	0,0950	0,1023
10. R в	0,070	0,049	0,098	0,070	0,071	0,358	0,0716	(K) 1,0
	=1	=1	=1	=1	=1	SS=5	S Mj =1	S D =1 S К =1

 

Определяем степень согласованности назначенных экспертами индивидуальных коэффициентов весомости. Результаты расчетов заносим в табл. 5.

Степень согласованности назначенных экспертами индивидуальных коэффициентов весомости определяют по коэффициенту вариации (S_j – среднее квадратическое отклонение):

Таблица 5

Определение степени согласованности мнений экспертов k_v

Показатель качества	Мjср	(Мjср-Мjе)2	S(Мjср- -Мjе)2	Sj	kv
1. h	0,1096	0,00027556	0,00013456	0,00015376	0,00004096	0,00008836	0,00069296	0,013162	0,12
2. b1	0,1190	0,00000900	0,00000900	0,00000900	0,00000900		0,00003600	0,003000	0,03
3. b2	0,1096	0,00027556	0,00015376	0,00015376	0,00004096	0,00021316	0,00083696	0,014500	0,13
4. S	0,0574	0,00010816	0,00024336	0,00007056	0,00015876	0,00008836	0,00066896	0,0129321	0,23
5. a	0,1094	0,00004356	0,00012996	0,00012996	0,00004356	0,00009216	0,00043890	0,0104749	0,10
6. l	0,1142	0,00000324	0,00006084	0,00006084	0,00000324	0,00036864	0,00055504	0,0117796	0,10
7. Скручив.	0,1144	0,00000256	0,00005776	0,00005776	0,00045796	0,00002116	0,00059736	0,012200	0,11
8. Кривизна	0,0998	0,00026244	0,00000324	0,00071824	0,00004624	0,00036864	0,00139864	0,0186991	0,19
9. Волнист.	0,0950	0,00044100	0,00000900	0,00048400	0,00000400		0,00094400	0,0153622	0,16
10. R в	0,0716	0,00000256	0,00051076	0,00069696	0,00000026	0,00000036	0,00121326	0,0174159	0,24

 

;

;

.

Согласованность индивидуальных коэффициентов весомости считают достаточной, если k_vj < [ k_v ] = 0,25.

В случае недостаточной согласованности коэффициенты весомости определяют групповым методом с взаимодействием (когда с индивидуальными мнениями знакомят всех экспертов и находят общее решение) [ РД 50-553-85 ].

Для первого значения k_v1 = 0,013162 / 0,1096 = 0,12.

Все коэффициенты вариации имеют достаточную согласованность.

Определяем оценку единичных показателей качества неравнополочного гнутого швеллера.

Относительные показатели по доверительному интервалу определяем по формулам (3) и табл. 6:

, (3)

где i = 1,..., n;

P_i - значение i- го показателя качества оцениваемой продукции;

P_i^б - значение i -го базового показателя;

n - количество показателей качества продукции, используемых для оценки.

Формулу выбирают в зависимости от того, увеличению или уменьшению относительно базового показателя отвечает улучшение качества продукции.

За базовые значения показателей Р_i^баз при определении относительных показателей качества принимаем допустимый интервал значений по нормативной документации Δ Р_i^ГОСТ (графы 2-5, табл. 6). Так, для стенки швеллера допустимый интервал по ГОСТ 8281-80 для профилей высокой точности - Р_i^баз =1,5 мм (см. табл. 1); повышенной точности - Р_i^баз = 3,0 мм; обычной точности Р_i^баз =4,0 мм.

Чем меньше доверительный интервал Р_i^факт, тем выше оценка, поэтому при определении относительных единичных показателей качества пользуемся формулой q_i = Р_i^баз / Р_i^факт. Тогда для стенки швеллера (профили обычной точности) q_i =4,0/2,56=1,56 (табл. 6, графа 9). Аналогично определяем значения q остальных единичных показателей профилей относительно базовых значений обычной, повышенной и высокой точности по ГОСТ 8281 (Россия), и относительно JIS G 3350 (Япония). Результаты заносим соответственно в графы 7-10

табл. 6.

Из таблицы видно, что оцениваемый профиль превышает требования российского стандарта для профилей обычной точности (все единичные показатели больше единицы), не соответствует требованиям японского стандарта по одному единичному показателю – ширине большей полки (b₂), по двум показателям – требованиям российского стандарта для профилей повышенной точности и по трем – для профилей высокой точности (значения их относительных показателей меньше единицы).

 

Таблица 6

Оценка относительных показателей качества неравнополочного гнутого швеллера
67х65х35х3,0 мм по доверительному интервалу [формула (3)]

 

Показатель качества	Допустимый интервал значений Δ РiГОСТ	Доверительный интервал Δ Piфакт	Оценка единичных показателей, qi
ГОСТ 8281-80	JIS G 3350 Япония	высокая точность	повышенная точность	обычная точность	Япония
высокая точность	повышенная точность	обычная точность
h, мм	1,5	3,0	4,0	3,0	2,56	0,59	1,17	1,56	1,17
b1, мм	1,0	2,0	3,0	3,0	2,90	0,34	0,69	1,03	1,03
b2, мм	1,5	3,0	5,0	3,0	3,20	0,47	0,94	1,56	0,94
S, мм	0,42	0,42	0,42	0,56	0,40	1,05	1,05	1,05	1,40
a, град	2,0	3,0	3,0	3,0	2,5	0,80	1,20	1,20	1,20
l, мм						1,67	1,67	1,67	1,67
Скруч, град/м	1,0	1,0	1,0	-	0,98	1,02	1,02	1,02	-
Крив., % / м	0,1	0,1	0,1	0,2	0,09	1,11	1,11	1,11	2,22
Волн., мм / м	2,0	2,0	2,0	-	1,8	1,11	1,11	1,11	-
R в, мм	5,0	5,0	5,0	-	3,5	1,43	1,43	1,43	-

 

Однако даже если продукция полностью соответствует требованиям нормативной документации по доверительному интервалу, возможны случаи производства несоответствующей продукции. Это может быть связано с неправильной настройкой оборудования, на котором производят данную продукцию. При этом доверительный интервал может быть значительно ΄уже требуемого по НД, но смещен, когда одна из его границ выходит за пределы допустимых значений. В этом случае следует проводить оценку качества продукции с учетом предельных значений показателей качества (Р_i^пр), например, по следующей формуле:

 

.

Оценку единичных показателей качества неравнополочного гнутого швеллера по предельным значениям произведем по формуле (4), предложенной Г.С. Гуном и Е.А. Пудовым (табл.7):

. (4)

Формула (4) позволяет оценивать качество продукции по свойствам, увеличение значений которых повышает или снижает качество продукции (технологии), а также для свойств, выход абсолютных значений которых из определенного интервала понижает уровень качества. Следует обратить внимание на тот факт, что при использовании данной формулы величина оценки не будет превышать единицы, но чем выше оценка, тем лучше качество оцениваемого объекта.

При определении оценки единичных показателей по формуле (4) в качестве базового значения единичных показателей примем номинальные значения этих показателей (графа 2 табл. 7), что более корректно, чем доверительный интервал. Тогда предельные значения показателей определим по формуле Р_i^пр =Р_i^номин ± DР_max.

(DР_max – это наибольшее из фактических значений, составляющих доверительный интервал). Д ля стенки швеллера Р_i^пр =67 – 1,52 = 65,48 мм » 65,5 мм. Аналогично следует определять предельные значения остальных показателей. Результаты заносим в графу 3 табл. 7. А фактическое значение показателя равно математическому ожиданию: P_i^факт=`х (см. ранее). Для стенки швеллера P_i^факт = 66,76 мм » 66,8 мм. Аналогично определяем фактические значения остальных показателей и заносим их в табл. 7.

Таблица 7

Оценка единичных показателей качества неравнополочного гнутого швеллера 67х65х35х3,0 мм по предельным значениям [формула (4)]

 

Пока-затель	Базовое значение показателя Рiбаз	Предельное значение показателя Рiпр	Фактическое значение показателя Piфакт	Оценка единичных показателей
h, мм	67,0	65,5	66,76	1-[(67-66,76)/(67-65,5)]2=0,9822
b1, мм	35,0	36,4	36,1	0,38
b2, мм	65,0	63,9	64,7	0,93
S, мм	3,0	2,85	2,95	0,89
a, град	90,0	90,5	89,8	0,84
l, м	9,0	9,03	9,02	0,56
Скруч, град/м		0,98	0,95	0,06
Крив., %/м		0,09	0,055	0,63
Волн., мм/м		1,8	0,9	0,75
R в, мм	3,0	3,5	3,2	0,84

 

Далее определяем уровень качества профиля по отношению к требованиям стандартов по формуле:

Тогда по отношению к требованиям российского стандарта уровень качества неравнополочного гнутого швеллера производства ЛПЦ-7 ОАО «ММК» для профилей обычной, повышенной и высокой точности соответственно равен:

К_у^высок = 0,93; К_у^обычн = 1,28; К_у^повыш = 1,13.

Чтобы определить уровень качества швеллера по отношению к требованиям стандарта Японии, следует пересчитать коэффициенты весомости на 7 показателей, которые нормируются стандартом JIS G 3350 (см. табл. 1).

Отсутствуют сведения по показателям № 7, 9, 10, непересчитанная суммарная весомость которых равна 0,1144+0,095+0,0716=0,281. Тогда весомость остальных показателей составит:

1–0,281=0,719 ® принимаем за единицу.

Тогда коэффициенты весомости для показателей по Японскому стандарту будут (для первого показателя «высота стенки h »: 0,1096/0,719=0,1526):

 

1 – h	2 – b1	3 – b2	4 – S	5 – a	6 – l	8 – Крив
0,1526	0,1655	0,1526	0,080	0,1524	0,1590	0,1390

К_у^Яп =1,17х0,1526 + 1,03х0,1655 + 0,94х0,1526 + 1,4х0,08+
+ 1,2х0,1524 + 1,67х0,159 + 2,22х0,139 = 1,36.

То есть, если мы правильно выбрали номенклатуру показателей качества, можно сказать, что производимые профили в 1,36 раза превышают требования Японского стандарта JIS G 3350 (1977), в 1,28 раза – Российского ГОСТ 8281-80 для профилей обычной точности, в 1,13 раза – повышенной точности и уступают требованиям ГОСТ 8281-80 для профилей высокой точности в 1,075 раза (1,0 / 0,93 = 1,075).

При необходимости сравнения качества швеллеров аналогичного сортамента, выпускаемых на различных агрегатах или по различным технологиям, можно определить комплексную оценку качества по формуле (5):

, (5)

,

где a_i и b_y – коэффициенты весомости доминирующих и компенсируемых показателей;

D_i – оценка качества по i доминирующему показателю;

K_y - оценка качества по y компенсируемому показателю.

Тогда (см. табл. 7):

К_к = (0,98^0,1181´0,38^0,1282´0,93^0,1181´0,89^0,0618´0,84^0,1178´0,56^0,123´
´0,06^0,1232´0,63^0,1075´0,75^0,1023+0,84´1,0)^0,5= 0,66.

С целью эффективного воздействия на качество продукции в первую очередь следует обратить внимание на единичные показатели, которые имеют наибольшие коэффициенты весомости и наименьшие значения единичных оценок. Из табл. 7 видно, что такими единичными показателями являются ширина меньшей полки b₁ (наибольший коэффициент весомости 0,1282 и низкая оценка - 0,38)и показатель скручивания (соответственно, 0,1232 и 0,06). Воздействие на технологический процесс с целью повышения оценки этих показателей позволит максимально улучшить качество неравнополочных гнутых швеллеров в целом.

Для определения основных причин низких оценок по выше указанным единичным показателям можно воспользоваться одним из методов выявления факторов и условий, влияющих на качество продукции.

Например, построим диаграмму Исикавы (рис. 4). Для построения причинно-следственной диаграммы обычно используют метод «мозгового штурма». При этом важно на стадии подготовки диаграммы привлечь к обсуждению как можно больше специалистов, чтобы диаграмма была полной, и в ней ничего не было упущено, поскольку фактор, не обозначенный на первоначальной стадии еще до построения диаграммы, на более поздних стадиях уже не появится.

 

 

 

Анализ причинно-следственной диаграммы (диаграммы Исикавы) показывает, что основной причиной невысоких оценок при существенной значимости таких геометрических показателей неравнополочного холодногнутого швеллера, как скручивание профиля и ширина меньшей полки, является неправильно рассчитанная калибровка валков, которая создает разные по величине крутящие моменты при подгибке разноширинных полок. Устранение данной причины позволит существенно повысить оценки указанных единичных показателей и качество холодногнутого неравнополочного швеллера в целом.

РЕКОМЕНДУЕМАЯ ЛИТЕРАТУРА

1. Рашников В.Ф., Салганик В.М., Шемшурова Н.Г. Квалиметрия и управление качеством продукции: Учебн. пособ. – Магнитогорск: МГТУ, 2000. – 184 с.
2. Фомин В.Н. Квалиметрия. Управление качеством. Сертификация: Курс лекций. – М.: Экмос, 2000. – 320 с.
3. Всеобщее управление качеством: Учебник / О.П. Глудкин, Н.М. Горбунов, А.И. Гуров и др. – М.: Радио и связь, 1999. – 600 с.
4. Управление качеством. Часть 1: Семь простых методов: Уч. пособие / Ю.П. Адлер, Т.М. Полховская, В.Л. Шпер, П.А. Нестеренко. – М.: МИСИС, 2001. – 138 с.
5. Федюкин В.К., Дурнев В.Д., Лебедев В.Г. Методы оценки и управления качеством промышленной продукции: Учебник. – М.: Филинъ - Рилант, 2000. – 328 с.
6. Гун Г.С. Управление качеством высокоточных профилей. – М.: Металлургия, 1984. – 151 с.
7. Анализ нормативных документов: Методич. указ. / Х.Н. Белалов, В.А. Кулеша, Н.А. Галкина. – Магнитогорск: МГТУ, 1999. – 10 с.
8. Шемшурова Н.Г., Моллер А.Б. Методы сбора и предварительного анализа данных в управлении качеством: Методич. указ. – Магнитогорск: МГТУ, 2002, – 38 с.