CC BY
28
УДК 621.9.048.6
Д.Н. МАЛЕЕВ, В.Х. АЛЬ-ТИББИ, Д.А. ЧИЛИКИН ОПТИМИЗАЦИЯ ПРОЦЕССА ЭЛЕКТРОАКУСТИЧЕСКОГО НАПЫЛЕНИЯ ПО КРИТЕРИЮ МИКРОТВЕРДОСТИ
Исследуется влияние акустических параметров (амплитуда ультразвуковых колебаний А, частота ультразвуковых колебаний / и угол между осью волновода и направлением вектора продольно-крутильных колебаний а ) процесса электроакустического напыления защитных покрытий на качество покрытия (микротвер-досгь). На базе проведенных исследований предложена математическая модель зависимости микротвердости покрытия от акустических параметров процесса и проверена адекватность полученной модели. Ключевые слова: акустические параметры, амплитуда ультразвуковых колебаний, частота ультразвуковых колебаний, продольно-крутильные колебания, математическая модель, адекватность модели.
Введение. Одним из перспективных методов нанесения токопроводящих износожаростойких защитных покрытий является метод электроакустического напыления - ЭЛАН [1,2]. Задача формирования качественного покрытия с высокими требованиями к шероховатости и микротвердости поверхностного слоя в процессе ЭЛАН является одной из важнейших при дальнейшем развитии этого метода для поверхностно-упрочняющей обработки деталей машин и формообразующего инструмента. Многовариативность выбора параметров процесса и их несомненное влияние на физико-механические свойства покрытия, необходимость исследования физической картины процесса требует проведения статистических исследований с целью обоснования режимов обработки и оптимизации показателей качества покрытия. Комплекс таких параметров существенно влияет на износостойкость и долговечность деталей машин и на другие их эксплуатационные характеристики. В задачу исследований входило изучение влияния акустических параметров процесса ЭЛАН на качество наносимого покрытия.
Методика исследований и некоторые априорные сведения. Исследование процесса ЭЛАН, безусловно, требует изучения влияния не только акустических, но и электрических параметров процесса на качество покрытия. Тем не менее, такая задача является слишком объемной в рамках одного исследования, поэтому в данной работе мы решили сосредоточиться для начала только на акустических параметрах процесса. Напряжение разрядного импульса как электрический параметр не претерпело изменений. По результатам предыдущих исследований [3,4] было выявлено, что удовлетворительную шероховатость покрытия (сопоставимую по классу с шероховатостью обрабатываемой поверхности) можно получить при работе на следующих усредненных по величине режимах: напряжение разрядного импульса и=32 В (при напылении твердого сплава ВК6М на сталь 45), емкость разрядного конденсатора С=12000 мкФ, частота подачи импульсов /г = 20 кГц. При этом скорость перемещения электрода составляла V = 0,04 м/мин.
Для разработки математической модели процесса был применен полнофакторный эксперимент (ПФЭ), позволяющий получить математическую модель, выбрать контролируемые параметры и определить движение по градиенту, т.е. направление, при движении по которому параметр оптимизации уменьшается или увеличивается быстрее, чем в любом другом направлении. Были запланированы три вида статистических исследований, в каждом из которых фиксировалось три важнейшие характеристики качества покрытия: отношение микротвердости покрытия к микротвердости подложки Нмп /Нмм, шероховатость покрытия Ra, мкм, толщина покрытия Н, мкм.
В качестве параметра оптимизации у1, было взято отношение НцП /НцМ, которое вычислялось после замеров на микротвердомере модели ПМТ-3, в соответствии с требованиями ГОСТ 9450-73.
В качестве переменных (факторов) были выбраны следующие акустические параметры ЭЛАН (рисунок): х1 - амплитуда ультразвуковых колебаний А (мкм), х2 - частота ультразвуковых колебаний / (кГц), х3 - угол между осью волновода и направлением вектора продольнокрутильных колебаний а (градусов).
Схематическое изображение угла а между осью волновода и направлением вектора продольно-крутильных колебаний
Для описания процесса была выбрана линейная модель уравнения регрессии
у = Ь0 + Ь1 х1 + Ь2 х2 + Ь3 х3 + Ь12 х1 х2 + Ь13 х1 х3 + Ь23 х2 х3 + Ь123 х1 х2 х3, (1)
где у - параметр оптимизации; bi - параметры модели; х; - факторы процесса ЭЛАН; / - номер параметра оптимизации, / =1,2,3..
Границы изменения факторов х1, х2, х3 определены из установочных исследований [3,4] и соответствуют данным, приведенным в табл.1.
Таблица 1
Значения исследуемых факторов
Уровни Факторы
А, мкм I, кГц а , град
Верхний 17 22,8 82
Нижний 4 16 10
Основной 10,5 19,4 46
Интервал варьирования 6,5 3,4 36
Кодовое обозначение *1 *2 *3
Планирование и проведение эксперимента. Был построен план матрицы планирования полнофакторного эксперимента типа 23, где 3 - число факторов. Кроме того, для обеспечения расчетов ввели фиктивную переменную х0, которая принимает во всех опытах значения +1, а также графы взаимодействия факторов [5]. Данные сведены в табл.2.
Таблица 2
Матрица планирования ПФЭ 23
Номер точки плана Значения факторов Комбинации произведений с эакторов
*0 *1 *2 *3 *1 *2 *1 *3 *2 *3 *1 *2 *3
1 + - - - + + + -
2 + + - - - - + +
3 + - + - - + - +
4 + + + - + - - -
5 + - - + + - - +
6 + + - + - + - -
7 + - + + - - + -
8 + + + + + + + +
Перед реализацией плана эксперимента опыты, предусмотренные матрицей, проводились в случайной последовательности. Порядок реализации опытов, выбранный по таблице равномерно-распределенных случайных чисел, представлен в табл.3.
Таблица 3
Рандомизация опытов
Порядок реализации Номер опыта
1 2 З 4 5 б 7 8
Серия 1 1 б 5 8 2 7 4 З
Серия 2 4 2 7 5 З 1 б 8
Серия 3 2 7 4 З 1 8 5 б
Серия 4 8 1 б 7 5 2 4 З
Серия 5 З 8 1 б 4 5 2 7
Нужное значение факторов устанавливали поочередно в каждой точке плана матрицы в порядке реализации текущей серии эксперимента в соответствии с матрицей планирования и
поддерживали эти значения постоянными в течение опыта. Параметр оптимизации у1(——)
Н фЛ
фиксировали по окончании опыта и заносили в табл.4.
Таблица 4
—
Результаты измерений параметра у1(——)
—фЛ
Номер опыта Результаты эксперимента в сериях
Y11 Y12 Y^ Y14
1 1.79 1.91 1.7З 1.8б
2 1.б1 1.54 1.58 1.4б
З 1.77 1.81 1.79 1.84
4 1.42 1.б9 1.З9 1.б4
5 1.58 1.41 1.б0 1.42
б 1.З2 1.17 1.24 1.21
7 1.49 1.27 1.41 1.З5
8 1.З4 1.17 1.28 1.21
Последующая обработка экспериментальных данных проводилась на ЭВМ, для чего была реализована программа в среде математического программирования Math Cad.
H П
Исследование зависимости —— = f (A,v,a). Воспроизводимость измерений проверяли пу-
H м
тем подсчета среднего показателя параметра оптимизации процесса при реализации параллельных наблюдений, также была рассчитана выборочная дисперсия воспроизводимости параллельных наблюдений для каждой точки плана матрицы. Однородность дисперсий проверили по критерию Кохрена.
Определили значимые коэффициенты регрессии и их значения с использованием критерия Стьюдента. Расчетные значения коэффициентов представлены в табл.5. В математическую модель процесса включили только значимые коэффициенты. Получили уравнение регрессии в виде
Y = 1.50625 - 0.11437 • X1 - 0.16437 • X3 . (2)
З41
Таблица 5
Коэффициенты уравнения регрессии
Коэффициент Значение t - критерий Значимость
Ъо 1,50625 96,55 Да
Ъ -0,11437 7,275 Да
Ъ2 -0,02062 1,28 Нет
Ъз -0,16437 10,51 Да
Ъ12 0,02125 1,35 Нет
Ъ13 0,015 1 Нет
Ъ23 -0,00625 0,4 Нет
Ъ123 0,01312 0,83 Нет
Проверка адекватности модели. Адекватность модели проверили по формуле Фишера. Табличное значение F =2,62 (при принятом количестве степеней свободы и уровне значимости). Поскольку Fpac4(Fma6n, то гипотеза адекватности принимается. Были найдены коэффициенты уравнения регрессии в действительных значениях. После простейших преобразований получили:
= 1,901 - 0,017A - 0,0046а
Н,м
Так как в уравнении регрессии все коэффициенты, стоящие при х2, незначимые, то данный фактор можно не учитывать. Фактор хх (амплитуда ультразвуковых колебаний А (мкм)) оказывает большее влияние на показатель параметра оптимизации по сравнению с фактором х3 (угол между осью волновода и направлением вектора продольно-крутильных колебаний а (градусов)). Это следует учитывать при контроле процесса. Следовательно, на микротвердость покрытия наибольшее влияние оказывает амплитуда ультразвуковых колебаний. Существенен также параметр угла наклона вектора продольно-крутильных колебаний. Частота колебаний влияет незначительно. Знаки коэффициентов уравнения регрессии отрицательны, что говорит о том, что увеличение микротвердости можно достичь уменьшением параметров A и а .
Методика оптимизации. Оптимизация проводилась методом крутого восхождения и начиналась от основных уравнений значимых факторов х, и х3 .Факторы изменялись пропорционально
1осн зосн
величинам коэффициентов регрессии с учетом их знаков. План крутого восхождения учтен при расчете в среде математического программирования Math Cad.
С целью проверки полученной модели оценили показатель параметра оптимизации по полученному уравнению регрессии, используя рассчитанные значения факторов в плане крутого восхождения. Это необходимо не для предсказания выхода процесса, а для определения направления движения. Данные "мысленных опытов " представлены в табл.6.
Таблица 6
"Мысленный эксперимент" для уравнения = 1,901 - 0,017A - 0,0046а
Нм
Номер шага крутого восхождения А, мкм а0 «Мысленное значение» НПn /НfM
1 9.5 38.2 1.558
2 8.6 31.4 1.605
3 7.7 24.6 1.652
4 6.8 17.2 1.702
5 5.9 10 1.751
6 5.07 2.8 1.8
7 4.24 4.4 1.846
Полученный "выход" показывает, что направление движения выбрано верно - микротвердость покрытия возрастает. Вместе с тем показано, что нет необходимости в проведении опытов 1-3, так как эффект варьирования параметров в указанных областях незначителен.
Реализация оптимальных режимов для микротвердости покрытия. В соответствии с планом крутого восхождения были проведены дополнительные эксперименты по вышеуказанной методике. Результаты экспериментов представлены в табл.7. Приведенные данные показывают, что экстремальное значение параметр Н/Нпринимает в точке 5 при следующих оптимальных
режимах: А = 5,9 мкм, а = 10° .
Таблица 7
Результаты дополнительных экспериментов
Номер шага А, мкм а0 H /лП / H фМ
3 7,7 24,6 1,91
4 6,8 17,2 2,08
5 5,9 10 2,15
6 5,07 2,8 2,03
Выводы. При постоянных параметрах напряжения разрядного импульса, емкости разрядного конденсатора и скорости перемещения электрода:
1) на микротвердость покрытия наибольшее влияние оказывает амплитуда ультразвуковых колебаний;
2) существенен также параметр угла наклона вектора продольно-крутильных колебаний;
3) частота колебаний влияет незначительно;
4) знаки коэффициентов уравнения регрессии отрицательны, т.е. увеличение микротвердости можно достичь уменьшением параметров A и а .
5) оптимальные значения акустических параметров ЭЛАН по критерию твердости: амплитуда ультразвуковых колебаний A =12 мкм, частота ультразвуковых колебаний f =20 кГц, угол между осью волновода и направлением вектора продольно-крутильных колебаний а =45°.
Библиографический список
1. Минаков В.С. Электроакустическое напыление / В.С. Минаков, А.Н. Кочетов // СТИН.
- 2003. - № 4. - С.32-35.
2. Аль-Тибби В.Х. Получение наноструктурных материалов методом электроакустического напыления: журнал / В.Х. Аль-Тибби, В.С. Минаков, Д.Д. Дымочкин // СТИН. - 2007. - № 4.
- С.28-32.
3. Минаков В.С. Разработка комплексных механических и электро-физических процессов обработки на основе использования энергии трансформируемых ультразвуковых колебаний: дис...д-ра техн. наук: 05.03.01 / - Ростов н/Д, 1989. - 516 с.
4. Кочетов А.Н. Барьерно-дислокационный механизм упрочнения деталей машин методом электроакустического напыления: дис. ... к-та техн. наук: 05.03.01 / - Ростов н/Д, 1998. - 241 с.
5. Кохановский В.А. Организация планирования эксперимента / В.А. Кохановский, М.Х. Сергеева. - Ростов н/Д: Издательский центр ДГТУ, 2008.
Материал поступил в редакцию 15.03.10.
D.N. MALEYEV, W.H. AL-TIBBI, D.A. CHILIKIN OPTIMIZATION OF THE ELECTRO-ACOUSTIC SPUTTERING PROCESS BY THE MICROHARDNESS CRITERION
Influence of acoustic parameters (amplitude of ultrasonic oscillations A (mcm), frequency of ultrasonic oscillations (Khz), and a corner between the axis of the wave guide and the vector direction of longitudinal-twisting oscillations (degrees)) process of an electro-acoustic dusting of the protective coatings on quality of coating (micro hardness) is investigated. On the basis of the executed research the mathematical model of dependence of micro hardness of a coating from acoustic parameters of process is developed and adequacy of the received model is checked.
Keywords: acoustic parameters, amplitude of ultrasonic oscillations, frequency of ultrasonic oscillations, longitudinal-twisting oscillations, mathematical model, adequacy of the model.
МАЛЕЕВ Данил Николаевич (р. 1985), аспирант первого года обучения кафедры «Автоматизация производственных процессов» Донского государственного технического университета. Окончил Донской государственный технический университет (2009).
Область научных интересов: ультразвуковые процессы и технологии.
Автор 7 публикаций.
АЛЬ-ТИББИ Висам Хусамович (р. 1979), кандидат технических наук (2006), преподаватель кафедры «Автоматизация производственных процессов» ДГТУ. Окончил Донской государственный технический университет (2002).
Область научных интересов: ультразвуковые процессы и технологии.
Автор 20 публикаций.
ЧИЛИКИН Денис Алексеевич (р. 1983), аспирант третьего года обучения кафедры «Автоматизация производственных процессов» Донского государственного технического университета. Окончил Донской государственный технический университет (2007).
Область научных интересов: ультразвуковые процессы и технологии.
Автор 8 публикаций.
