УДК 621.923
РЕШЕНИЕ ТЕПЛОВОЙ ЗАДАЧИ ПРИ ШЛИФОВАНИИ С УЧЕТОМ ИМПУЛЬСНОГО ХАРАКТЕРА ТЕПЛОВОГО ПОТОКА
В.А. Смирнов
Ижевский государственный технический университет, Воткинский филиал E-mail: smivial@mail.ru
Приведено решение нестационарного нелинейного двумерного уравнения теплопроводности при плоском шлифовании периферией круга с учетом импульсного характера теплового потока. Предложено представлять тепловой поток в зоне контакта шлифовального круга и заготовки в виде суммы распределенной и импульсной частей. Полученные результаты позволяют прогнозировать изменение температуры заготовки при шлифовании.
Ключевые слова:
Шлифование, температурное поле, тепловой поток Key words:
Grinding, temperature field, heat-flow rate.
Исследование температурного поля заготовки является одной из центральных проблем теории шлифования. Такие ученые как А.И. Исаев, С.С. Силин, С.Г Редько, А.Н. Резников, В.И. Островский, В.А. Сипайлов, А.В. Якимов, Е.П. Калинин, А.Ш. Хусаинов, Е.С. Киселев и др. посвятили данной проблеме множество своих работ [1-5 и др.]. В.А. Сипайловым и А.В. Якимовым получено аналитическое решение двумерного нестационарного линейного уравнения теплопроводности при действии равномерно распределенного полосового источника теплоты постоянной интенсивности, рассмотрены случаи шлифования сплошным и прерывистым кругом [3, 5]. Позже в работах [1, 4] получено численное решение нелинейной теплофизической задачи при шлифовании для тел различной конфигурации. В работах А.Н. Резникова [2] из теоретико-вероятностных соображений получены основные соотношения для расчета локальной и средней температуры при шлифовании. Результаты, полученные перечисленными учеными, внесли большой вклад в развитие теории теплофизики шлифования, однако остался ряд нерешенных проблем.
Существенный недостаток разработанных математических моделей - недостаточное внимание уделяется граничным условиям в зоне контакта шлифовального круга и заготовки. Как правило, в известных математических моделях проводится усреднение теплового потока по площадке контакта круга и заготовки. В результате не учитывается импульсная природа теплового потока из-за действия локальных источников теплоты (отдельных режущих зерен). С другой стороны в вероятностных моделях для расчета локальной температуры принимаются допущения, справедливость которых проверить не представляется возможным, вводятся переменные, о которых практически ничего не известно. В результате существующие расчетные схемы лишь частично соответствуют реальному процессу шлифования. В настоящей статье мы попытались устранить существующие недостатки.
Обзор проведенных исследований показал, что недостаточно внимания уделяется колебаниям температуры заготовки. Экспериментальные данные, полученные нами с помощью метода полуискус-ственной термопары, показывают, что колебания контактной температуры могут достигать 25...30 % от усредненного значения. Амплитуда колебаний температуры нарастает с увеличением времени шлифования. Как отмечает А.Н. Резников, локальная температура часто ограничивает режимы шлифования, поскольку именно под воздействием этих температур возникают структурные изменения и прижоги в поверхностных слоях заготовки. Локальная температура влияет на изнашивание и разрушение режущих зерен, а значит на работоспособность и стойкость режущего инструмента [2].
Колебания температуры обусловлены действием единичных зерен и изменением мгновенной производительности вследствие вибраций или неравномерного износа круга. Влияние колебаний силы резания на изменение контактной температуры были рассмотрены нами в [7]. В рамках данной статьи рассмотрим влияние единичных зерен на колебания температуры заготовки.
Рассмотрим плоское шлифование периферией круга. Обрабатываемую заготовку представим в виде пластины размерами ЬУ^Ь1 (рис. 1). Обозначим 0(1,у,т) - температурное поле заготовки. Рассматривается нелинейное нестационарное двумерное уравнение теплопроводности:
сЛ<Э)Рз = Шу(Аз(©)У©( у, у),
от
где Яз(0), сз(0) - соответственно, коэффициент теплопроводности и удельная теплоемкость материала заготовки, зависящие от температуры 0; т -время.
На рис. 1 р - тангенциальная составляющая силы резания; q - плотность теплового потока, поступающего в заготовку; У3 - скорость подачи заготовки; Гкр - окружная скорость вращения шлифовального круга.
N = Nо +
(1)
Рис. 1. Расчетная схема для теплофизической задачи
Нелинейное уравнение теплопроводности выбрано потому, что при шлифовании возникают значительные контактные температуры, при этом теплофизические свойства материала заготовки изменяются в значительной мере (до двух раз).
На верхней границе расчетной области зададим граничные условия. В зоне контакта круга с заготовкой задается граничное условие 2-го рода - известен тепловой поток через границу:
Я3(©) д0 = Ч( г,т).
ду
Предполагаем, что смазывающе-охлаждающая жидкость (СОЖ) практически не проникает в зону контакта шлифовального круга и заготовки. Подобная картина наблюдается при подаче СОЖ поливом и при сухом шлифовании.
Вне зоны контакта задаем граничное условие 3-го рода - охлаждение потоком СОЖ с заданным коэффициентом теплоотдачи аж и температурой 0ж:
¿3(0) ^ = «.(0-0 ж).
ду
Закон распределения плотности теплового потока на площадке контакта шлифовального круга с заготовкой представим в следующем виде (рис. 2). Тепловой поток состоит из распределенной и импульсной частей (д0 и д^. Распределенная часть д0 обусловлена работой резания, выполненной активными зернами, уже прошедшими через зону резания, а также работой трения в межзерновом пространстве между связкой круга и материалом заготовки. д0 обеспечивает общее повышение температуры на площадке контакта шлифовального круга и заготовки. Распределенная часть д0 перемещается в направлении оси I со скоростью подачи заготовки Уз. Импульсная часть теплового потока д1 возникает вследствие резания единичными зернами и перемещается относительно д0 с окружной скоростью вращения круга Гкр. Величина импульсной части д1 является случайной величиной, так как положение режущих зерен и величина их выступания из связки являются случайными величинами.
На рис. 2 йъ- средний диаметр площадки износа зерен круга; Ьк - длина контакта шлифовального круга с заготовкой; В - ширина круга.
Для расчета величин q0 и д1 используем следующие условия. Во-первых, суммарная мощность, выделяемая источниками q0и д1, должна быть равна общей мощности теплового потока в зоне резания:
Рис. 2. Закон распределения плотности теплового потока по площадке контакта
Величину N рассчитаем как:
N = 40бщ БЬк>
где 4общ
РУ
2 кр
БЬ
плотность теплового потока,
поступающего в заготовку; Р1 - тангенциальная составляющая силы шлифования (рис. 1); ав - коэффициент ввода, показывающий долю тепла, уходящей в заготовку (может быть рассчитан по формуле В.А. Сипайлова [3]).
Величину Иф находим аналогично N
о = Чо ВЬк.
Для расчета величины используем формулу:
^і = 4і
ті
где д 1 - среднее значение дх\ по.р. - количество одновременно режущих зерен:
Ь
где преж - количество режущих зерен на периферии круга (может быть найдено по формуле А.Н. Резникова [2]), Я - радиус шлифовального круга.
Средний диаметр площадок износа зерен круга йи найдем по формуле:
- I 8Аа2
= Ха
'3К (1 -£)’
где -=1,0176-10^ - средневероятный размер зерна [2]; N - зернистость шлифовального круга по ГОСТ 3647-80; А - относительная опорная площадь поверхности круга, %; а - коэффициент формы зерна [2]; К=4(60-2С) - концентрация режущего материала [2]; С - номер структуры шлифовального круга по ГОСТ 3647-80; є - относительная критическая глубина заделки зерен.
в
После преобразований условия (1) получим:
Я ~1яЯг + д°Б = д°бщБ •
8К
(2)
Во-вторых, введем безразмерный параметр кн -наибольший коэффициент неравномерности распределения плотности теплового потока: шах(^1)
к =-
Яо
По определению кн>0. Параметр кн предлагается определять на основании экспериментальных исследований.
Учитывая случайный характер теплового потока дь запишем:
Яі = ЯоК х„. (3)
где хч - случайная величина, имеющая равномерное распределение в интервале [0;1]. Среднее значение 2г=0,5.
Для расчета д0 используем среднее значение ^=0,5 в уравнениях (2) и (3):
= 16 ЯБ
Яо = Яобщ 16 ЯБ + Кн і>реж '
Далее из уравнения (3) находим д1.
Найдем частоту появления режущих зерен на площадке контакта. Рассматривать будем площадку шириной Необходимо определить среднее время между контактами режущих зерен с поверхностью заготовки. Среднее количество режущих зерен на ширине іи:
4,
И . = П —-•
4И реж Б
Время одного оборота круга составляет '2тЯ • Следовательно, среднее время между
У„
контактами зерен на ширине определяется следующим образом:
2пЯВ
—
п У і
реж кр и
тод имеет высокую точность при наличии больших градиентов температуры, характерных для шлифования. На базе языков программирования С# и С++ нами разработана прикладная вычислительная система, реализующая используемый численный метод.
Результаты расчетов температуры на поверхности заготовки при различных значениях коэффициента кн представлены на рис. 3. Исходные данные для расчета: обрабатываемый материал - 9ХС; шлифовальный круг - электрокорунд белый; Гкр=35 м/с; К=0,05 м/с; глубина резания /=0,01-10-3 м; N=40; С=6; ^ =0,2-10-3 м.
Фактическое время между контактами зерен с заготовкой - величина случайная, поэтому для его расчета используем формулу:
т1 = ХЛ’
где Хт - случайная величина, имеющая равномерное распределение плотности вероятности в интервале [0;2].
Для решения поставленной задачи используем численный метод конечных разностей, который подробно описан в [6]. Для аппроксимации уравнения теплопроводности используем явную двухслойную схему, которая является линейной комбинацией двух известных схем: явной шеститочечной и Алена-Чена. Производная д0/ду в граничных условиях заменяется конечной разностью третьего порядка точности. Используемый численный ме-
Рис. 3. Изменение температуры на поверхности заготовки при к: а) 0; б) 10; в) 20
В результате расчета получена одна из реализаций случайного процесса. По рис. 3 видно, что размах колебаний температуры нарастает с увеличением коэффициента кн. Следует отметить, что график изменения средней контактной температуры
Т. С
Т. С
*1 =
X. С
в целом соответствует результатам, полученным В.А. Сипайловым в работе [3].
В качестве критерия оценки колебаний температуры вследствие работы активных зерен используем величину Rв - размах колебаний температуры. На рис. 4 показана типичная зависимость Rв от коэффициента ^
к н
Рис. 4. Типичная зависимость И0(к„)
Рис. 4 позволяет на основании экспериментальных данных определить значение ^ характерное для рассматриваемой шлифовальной операции. Можно также проследить изменение ^ на протяжении периода стойкости круга.
Пример экспериментального графика температуры заготовки из стали 9ХС представлен на рис. 5.
Рис. 5. Экспериментальные исследования температуры на поверхности заготовки
На экспериментальном графике, рис. 5, видны относительно низкочастотные колебания температуры вследствие изменения силы резания [7], а так-
же относительно высокочастотные колебания температуры вследствие работы активных зерен. На рис. 6 представлен теоретический график изменения температуры, полученный для тех же исходных данных (с учетом результатов [7]).
Рис. 6. Теоретический график изменения температуры на поверхности заготовки
Таким образом, математические модели, разработанные нами в [7] и данной статье, подтверждаются экспериментальными данными. Расхождение результатов экспериментальных и теоретических исследований не превышает 3 %.
Выводы
На основе нестационарного нелинейного двумерного уравнения теплопроводности разработана математическая модель для расчета температурного поля заготовки при шлифовании с учетом импульсного характера теплового потока. Показана возможность представлять тепловой поток в зоне контакта шлифовального круга и заготовки в виде суммы распределенной и импульсной частей. Математическая модель позволяет прогнозировать среднюю контактную и мгновенную температуру заготовки при шлифовании с учетом действия импульсных источников теплоты от активных зерен, благодаря чему обеспечивается требуемое качество шлифованных поверхностей.
Статьярекомендована Оргкомитетом V Международной научно-технической конференции «Современные проблемы машиностроения», г. Томск, 2010.
СПИСОК ЛИТЕРАТУРЫ
Киселев Е.С., Ковальногов В.Н. Теплофизический анализ концентрированных операций шлифования. - Ульяновск: УлГТУ, 2002. - 140 с.
Резников А.Н. Теплофизика процессов механической обработки материалов. - М.: Машиностроение, 1981. - 279 с., ил. Сипайлов В.А. Тепловые процессы при шлифовании и управление качеством поверхности. - М.: Машиностроение, 1978. - 167 с., ил.
Хусаинов А.Ш. Повышение производительности бездефектного шлифования заготовок клиновидных изделий: дис. ... докт. техн. наук: 05.03.01. - Ульяновск, 2006. - 425 с.
Якимов А.В. Прерывистое шлифование. - Киев: Вища школа, 1986. - 175 с.
Смирнов В.А. Численный метод решения нелинейной двумерной задачи теплопроводности при шлифовании // Современные технологии. Системный анализ. Моделирование. -2010. - № 3. - С. 23-26.
Смирнов В.А. Влияние колебаний силы резания на контактную температуру при шлифовании // Современные проблемы машиностроения: Труды V Междунар. научно-техн. конф.; ТПУ. - Томск: Изд-во ТПУ, 2010. - С. 625-630.
Поступила 21.02.2011 г.
т. С