УДК 65.018
ОЦЕНКА ЭФФЕКТА ОТ ВНЕДРЕНИЯ ДВУХКОНТУРНОЙ СИСТЕМЫ АКТИВНОГО КОНТРОЛЯ И РАННЕГО ОБНАРУЖЕНИЯ «ОСОБОЙ» ПРИЧИНЫ ИЗМЕНЧИВОСТИ НА ОПЕРАЦИЯХ МЕХАНООБРАБОТКИ
© 2018 С.И. Клейменов1, А.В. Заятров1, В.Н. Козловский1, А.Н. Чекмарев2, А.В. Седельников2
1 Самарский технический университет 2 Самарский национальный исследовательский университет имени академика С.П. Королёва
Статья поступила в редакцию 06.12.2018
В работе представлены результаты анализа эффективности внедрения двухконтурной системы активного контроля и раннего обнаружения «особой» причины изменчивости, на операциях ме-ханообоработки.
Ключевые слова: качество; эффективность; механообработка; автомобилестроение.
Для оценки эффективности от внедрения новых средств контроля или модернизации уже имеющегося необходимо выполнение расчётов, позволяющих сравнить преимущества, достигаемые за счёт использования средств контроля [7-8]. Применительно к продукции механоо-брабатывающего производства можно провести оценку по характеру распределения контролируемых параметров для партии обработанных на данном оборудовании средств контроля [910]. Для сравнения можно воспользоваться моментами распределения, описывающих свойства распределений [1-3].
Первый момент - математическое ожидание или координата центра распределения. В общем случае центр распределения определяется как точка на оси абсцисс, слева и справа от которой вероятность появления различных значений случайной величины равны между собой.
X = J xp(x )dx,
(1)
Второй момент - дисперсия случайной величины Б, характеризующая рассеяние отдельных значений от центра распределения. Действующее значение рассеяния или среднее квадрати-ческое отклонение является более удобной для оценки характеристикой и имеет размерность случайной величины:
I & - X )2
a = ^J D =11 i=1
(3)
Для получения несмещённой оценки дисперсии стандартного отклонения пользуются формулой:
I & - X )2
S = ll ^-
n -1
(4)
где p(x) - плотность распределения вероятностей.
При дискретных отсчётах, что имеет место в данном случае вместо вычисления интеграла допустимо перейти к вычислению среднего арифметического:
_ n X
X = IX. (2)
i=1 n
Клейменов Сергей Иванович, соискатель ученой степени кандидата технических наук, заместитель директора ООО «Бизнес-Консалт», г. Тольятти Заятров Алексей Викторович, кандидат технических наук.
Козловский Владимир Николаевич, доктор технических наук, профессор. E-mail: kozlovskiy-76@mail.ru Чекмарев Анатолий Николаевич, доктор технических наук, профессор.
Седельников Андрей Валерьевич, доктор технических наук, доцент.
В этой формуле вместо количества образцов п используется количество степеней свободы п-1. Таким образом, первая формула используется при расчёте по всей генеральной совокупности, во втором случае - по выборке. При этом разница в расчёте не заметна при большом объёме рассматриваемой выборки (больше 100).
Третий момент характеризует асимметрию или скошенность распределения относительно его среднего значения:
Л = 11Г Xl-X
n i=11 S
(5)
На практике для суждения о степени симметричности большее применение имеет безразмерный коэффициент асимметрии, дающий несмещённую оценку:
A = I
(n - i).(n - 2 ))=1
-л3
X - X S
(6)
Четвёртый момент характеризует протяжённость распределения и позволяет оценить величину остроконечности распределения по отношению к нормальному:
n
3
На рисунках 1 и 2 представлены гистограммы размеров обработанных на шлифовальном станке деталей с наложенным графиком функции плотности распределения вероятностей нормального закона со старым средством контроля и с использованием двухконтурной структуры активного контроля соответственно. Помимо замены самого средства активного контроля дополнительно проводились работы по модернизации шлифовального станка. Дополнительно в таблицу 2 сведены статистические показатели полученных результатов. Все значения размера обработанных деталей представлены в относительном виде относительно центра поля допуска на обрабатываемую деталь - 26,975 мм.
Анализ рисунков 1 и 2, а также данных в таблице 2 показывает положительный эффект от результатов модернизации контрольного оборудования и модернизации оборудования. В
Таблица 1. Результаты испытаний средств контроля
№ детали Время' ¿2 , МКМ Т' ОБ, сек Время Т ОБ, сек ^ык , мкм/сек Ь2 , МКМ
1 12:29 8 25,5 10:18 27,4 13,2 0
2 12:31 11 25,9 10:19 29 9,4 4
3 12:34 7 25,6 10:21 28,3 10,3 4
4 12:35 9 25,5 10:22 28,3 11,6 4
5 12:37 3 23,4 10:23 27,9 11,9 2
6 12:38 5 25,6 10:24 28 12,2 2
7 12:39 5 25,7 10:32 28,3 11,9 3
8 12:40 6 26,2 10:33 28,1 11,0 2
9 12:41 8 26 10:34 28,2 11,0 4
10 12:42 7 25,8 10:35 28 12,2 4
11 12:43 6 25,7 10:36 27,8 12,7 1
12 12:44 7 26 10:37 28,7 8,2 4
13 12:46 4 25,5 10:38 28,4 10,8 3
14 12:47 5 25,3 10:40 25,8 12,4 4
15 12:49 1 23,2 10:41 28,3 10,4 2
16 12:50 4 26,1 10:42 28,6 9,7 5
17 12:51 5 25,8 10:43 28,2 11,7 3
18 12:52 3 26 10:44 28,5 8,7 4
19 12:53 7 26,2 10:45 28 11,9 2
20 12:54 1 26,2 10:46 28,2 11,3 3
21 12:55 4 26,3 10:47 28 11,9 2
22 12:56 7 26,2 10:52 28,1 11,6 3
23 12:57 6 26 10:53 28,3 10,0 4
24 12:58 5 26 10:54 28,3 10,5 2
25 13:00 2 23,8 10:55 27,6 13,5 1
26 13:01 0 25,5 10:57 28,3 10,8 4
27 13:02 5 26 10:58 28 11,8 2
28 13:03 3 25,4 10:59 28 13,7 0
29 13:04 6 26,3 11:00 27,9 11,9 1
30 13:05 5 26,3 11:10 28,3 11,4 4
Ч1Ш- (7)
Обычно оценку проводят по коэффициенту эксцесса, который позволяет судить о протяжённости рассматриваемого распределения по отношению к нормальному (для которого в =3):
и•(„+1) »(х-хУ _ 3-(п-1)2
в (п _ \)-(п _ 2)(п _ 3)Ь I 5 ) (п _ 2)(п _ 3). (8)
Испытания выполнялись на оборудовании одного из ведущих российских предприятий по выпуску автомобилей, в частности круглошли-фовальном станке производства завода Пупка Юктёа, для наружной шлифовки. На указанном оборудовании выполняется обработка детали с допуском: 027-0,015 - 27-0,036 мм. Результаты замеров до модернизации и после сведены в таблицу 1. Показатели отмеченные знаком «'» относятся к данным, полученным на оборудовании до модернизации.
-4 -2 0 2 4 6 8 10 12 14
Рис. 1. Гистограмма распределения замеров на этапе послеоперационного контроля до модернизации оборудования
Таблица 2. Статистические показатели
Наименование показателя Значение до модернизации Значение после модернизации
Максимальное значение в выборке 11,0 5,0
Минимальное значение в выборке 0,0 0,0
Среднее арифметическое выборки 5,17 2,77
Стандартное отклонение 2,45 1,33
Коэффициент асимметрии -0,04 -0,48
Коэффициент эксцесса 0,25 -0,64
Рис. 2. Гистограмма распределения замеров на этапе послеоперационного контроля после модернизации оборудования
частности, более чем в 2 раза сократилось величина рассеяния итогового размера детали от центра распределения, а также сам центр распределения сместился ближе центру поля допу-
ска контролируемого параметра. Коэффициент асимметрии говорит о большем, чем до модернизации смещении результатов относительно среднего значения влево. Коэффициент эксцес-
са также увеличился и сменил знак по отношению к первоначальному, что свидетельствует о снижении величины остроконечности по отношению к нормальному закону распределения. Таким образом, для дальнейшего смещения центра распределения к центру поля допуска контролируемого параметра требуется провести ряд дополнительных настроек.
Выявив положительный эффект от модернизации используемого оборудования определим возможность раннего выявления момента разладки технологического процесса. Для этого сопоставим данные, получаемые на этапе активного контроля основного контура в темпе технологического процесса - конечную скорость снятия припуска в момент окончания обработки, с реальным размером детали, получаемым на этапе послеоперационного контроля дополнительного контура [4 - 6].
Для этого воспользуемся данными таблицы 1 и рассчитаем отклонение конечной скорости снятия припуска ЛУМК , в качестве номинального значения выбрав среднее:
^ (9)
А также отклонение размера детали ЛТ2 для каждой детали:
(10)
детали в выборке, результаты проиллюстрированы на рисунке 4.
Анализ контрольных карт, представленных на рисунке 4, уже достаточно информативен и позволяет отследить общие тенденции на обоих графиках, в особенности для карты средних. Коэффициент корреляции, рассчитанный для значений отклонения конечной скорости снятия припуска и отклонения размера детали, можно определить по следующему выражению:
Я =
±((г -АЬг)(( -АУмк)
I=1
£ (г -А Ь 2) £ ( -АУмк )
. (11)
I=1
I=1
АУмк, = Умк, - Умк.
АТ2 г - Т2 г _ Т .
Результаты представлены на рисунке 3.
В подобном виде диаграммы малоинформативные и по ним сложно судить о степени соответствия между данными, получаемыми в темпе технологического процесса и в статических условиях дополнительного контура. Проведём простейшую статистическую обработку методов контрольных карт средних и размахов по 3
Коэффициент корреляции получился .=-0,8987, что подтверждает высокую степень связи между анализируемыми параметрами. В этой связи можно рассчитать коэффициенты аппроксимирующего выражения для отклонения размера через показатели отклонения конечной скорости съёма припуска, которое будет иметь вид:
АЬ = Л -АУМК + Л2, (12)
где А1 и А2 - постоянные коэффициенты, для конкретного оборудования и параметров обработки.
В рассматриваемом случае выражение (12) примет вид:
АЬ = -0,769- АУМК - 0,01,
(13)
Представим на одном графике (рисунок 5) контрольную карту средних значений отклонения размеров деталей, полученных непосредственно в дополнительном контуре ЛТ2 с отклонением размеров деталей, рассчитанных по
Рис. 3. Точечные диаграммы конечной скорости снятия припуска и размеров деталей
Карта средних -т-
5 6 7
Карта размахов
9 10
№ выборки
-^Осн. конт. ^^Доп. КОНТ.^- А
Г/
1 1 Т 1
2 3 * 5 < 5 7 ; 8 9 1
№ выборки
Рис. 4. Контрольные карты средних и размахов для отклонения скорости и размера 1.5
а 0.5
-0,5
-1
-1.5
л ж Ж 41 ж ^ /с
п \ я * Ж С1 ^ J Л/1 1 ^^ ■*
■ N г г ■ ^
1 \ V.!
1 Доп. когг. ■Осн. конт, I I
10
№ выборки
Рис. 5. Карта средних значений отклонения размера деталей для основного ДЬ1 и дополнительного ЛГ2 контуров
выражению (13) на основе отклонения конечной скорости снятия припуска в момент окончания обработки в основном контуре ДЬХ.
Как видно из сравнения представленных графиков рассчитанные значения отклонения размера детали в основном контуре хорошо повторяют результаты статического контроля дополнительного контура. Максимальные расхождения между отклонением размера, полученном в динамическом режиме и статическом режимах прослеживаются только на первых двух выборках, достигая максимальной разницы в 0,4 мкм.
Таким образом, данные получаемые в темпе технологического процесса могут быть использованы для оценки качества готовых деталей,
являясь довольно точным показателем отклонения геометрических параметров детали от номинального значения. Также это позволит использовать данные параметры для раннего обнаружения причин изменчивости, выводящих технологический процесс из устойчивого состояния.
Стоит отметить, что величины погрешности обработки ЛЬ, даже при постоянстве времени отвода шлифовального круга в момент окончания обработки, близки к результатам, полученным во втором контуре системы контроля только в случае больших значений конечной скорости снятия припуска УМК. Если значения УМК находятся в диапазоне от 5 мкм/с и выше,
то совпадение динамической погрешности размера, полученное в основном контуре с результатами статических измерений в дополнительном контуре велико. В случае, когда значения УМК находятся в диапазоне до 5 мкм/с, что часто имеет место в производственной практике, увеличивается доля влияния прочих возмущающих факторов - температурный режим обработки, состояние шлифовального круга и периодичность его правки, качество СОЖ и её количество, поступающее в зону обработки. В подобных ситуациях требуется контроль дополнительных факторов для более точных вычислений.
Помимо проблемы повышения качества выпускаемой продукции в направлении увеличения точности обработки, особенно при модернизации уже имеющегося оборудования возникает вопрос об экономической целесообразности подобных улучшений. Современные приборы как активного, так и послеоперационного контроля в совокупности с измерительной оснасткой (измерительное приспособление, датчики контроля, эталоны и т.п.) стоят довольно дорого. Помимо стоимости самих приборов необходимо учитывать затраты, связанные с их обслуживанием в период эксплуатации (ремонт, настройка и калибровка). Если вопрос о целесообразности установки на новом оборудовании, впервые вводимом в эксплуатацию, средств активного и послеоперационного контроля не встаёт, то при модернизации он чрезвычайно актуален. Для этого проводится экономический расчёт, причём оцениваются показатели как для уже функционирующего оборудования, так и для нового с тем, чтобы сопоставить разницу между текущими показателями и величиной положительного эффекта, который может быть достигнут после внедрения новых средств контроля.
СПИСОК ЛИТЕРАТУРЫ
1. Годлевский В.Е. Статистические методы анализа качества. Самара, НВФ «Сенсоры. Модули. Системы», 1998. 103 с.
2. Заятров А.В., Козловский В.Н. Комплексная оценка качества и надёжности электрооборудования транспортных средств: монография. Самара: СамНЦ РАН. 2014. 176 с.
3. Заятров А.В. Использование двухконтурных систем активного контроля для повышения качества продукции / / Молодежная наука в развитии регионов: материалы Всерос. (с междунар. участием) науч.-практ. конф. студентов и молодых ученых (Березники, 25 апреля 2018). Пермь: Изд-во Перм. нац. исслед. политех. ун-та, 2018. С. 84 - 87.
4. Кобзарь А.И. Прикладная математическая статистика. Для инженеров и научных работников. М.: ФИЗМАТЛИТ, 2006. 816 с.
5. Новицкий П.В., Зограф И.А. Оценка погрешностей результатов измерений. Л.: Энергоатомиздат. Ле-нингр. отд-ние, 1991. 304 с.
6. Решетов А.Г. Автоматизация шлифования и размерного контроля деталей. СПб.: Политехника, 2003. 193 с.
7. Козловский В.Н., Горбачевский Н.И., Сорокин А.Г. Аналитический комплекс прогнозирования надежности электромобилей и автомобилей с комбинированной силовой установкой // Вестник Казанского технологического университета. 2014. Т. 17. № 3. С. 227-229.
8. Козловский В.Н. Обеспечение качества и надежности системы электрооборудования автомобилей: Дис. ... докт. техн. наук. Моск. гос. автомобил.-до-рож. ин-т (техн. ун-т). Тольятти, 2010
9. Гречников Ф.В., Яковишин А.С., Захаров О.В. Минимизация объема измерений плоских поверхностей деталей при сборке // Сборка в машиностроении, приборостроении. 2018. № 2. С. 56-58.
10. Гречников Ф.В., Яковишин А.С., Захаров О.В. Минимизация объема измерений при контроле цилиндрических поверхностей на основе статистического моделирования // Вестник Пермского национального исследовательского политехнического университета. Машиностроение, материаловедение. 2017. Т. 19. № 4. С. 101-110.
ASSESSMENT OF THE EFFECT OF THE INTRODUCTION OF A TWO-CIRCUIT SYSTEM OF ACTIVE CONTROL AND EARLY DETECTION OF A "SPECIAL" CAUSE OF VARIABILITY IN MACHINING OPERATIONS
© 2018 S.I. Kleimenov1, A.V. Zayatrov1, V.N. Kozlovsky1, A.N. Chekmarev2, A.V. Sedelnikov2
1 Samara Technical University 2 Samara National Research University named after Academician S.P. Korolyov
The paper presents the results of an analysis of the effectiveness of the implementation of a two-circuit active control system and the early detection of a "special" cause of variability in mechanical processing operations.
Keywords: quality; efficiency; mechanical processing; automotive industry.
Sergey Kleimenov, Candidate for a Scientific Degree of Candidate of Technical Sciences, Deputy Director of Business Consult LLC, Togliatt
Alexey Zayatrov, Candidate of Technics, Professor. Vladimir Kozlovsky, Doctor of Technics, Professor. E-mail: kozlovskiy-76@mail.ru Anatoly Chekmarev, Doctor of Technics, Professor. Andrei Sedelnikov, Doctor of Technics, Associate Professor.