ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ И
ПРОЦЕССЫ УПРАВЛЕНИЯ N.1, 2021 Электронный журнал, рег. Эл № ФС77-39410 от 15.04.2010 ISSN 1817-2172
http://diffjournal.spbu.ru/ e-mail: _ jodiff@mail.ru
Численное моделирование
Численно - аналитическое исследование турбулентного течения смазочно-охлаждающей жидкости при взаимодействии с деталью и инструментом в процессе обработки глубоких отверстий
1 * 2 ** 2 *** 2 **** Игумнов Л.А , , Грезина А.В. , , Метрикин В.С. , , Панасенко А.Г. ,
1Институт Механики Нижегородского государственного университета
им. Н.И. Лобачевского Нижегородский государственный университет им. Н.И. Лобачевского
e-mail:
*
igumnov@mech.unn.ru
**
al eksandra-grezina@yandex .ru
***
v.s.metrikin@mail.ru
****
a.g.panasenko@ yandex.ru
Аннотация. В работе приведена методика моделирования взаимодействия смазочно-охлаждающей жидкости (СОЖ) с обрабатываемой деталью и инструментом при растачивании глубоких отверстий. Моделирование основывается на численном решении системы нелинейных дифференциальных уравнений в частных производных, описывающих турбулентное течение (СОЖ) и перенос тепловой энергии в каналах детали и инструмента. Жидкость считается несжимаемой, а для замыкания системы дифференциальных уравнений используется стандартная к — £ модель турбулентности. Геометрическое моделирование объема проточной части проводилось в системе автоматизированного проектирования SolidWorks. Численное моделирование распределения давления, скоростей, температуры в потоке СОЖ, температуры обрабатываемой детали и режущего инструмента - с помощью программного комплекса FlowVision HPC. Получены распределения давления, скоростей, температуры в потоке СОЖ, а также распределение температур в зоне взаимодействия обрабатываемой детали и режущего инструмента в процессе теплообмена, и интегральные гидродинамические характеристики. Разработаны расчетные модели, позволяющие проводить вычислительные эксперименты при различных входных данных и способах подвода и отвода СОЖ. Сравнение результатов численного моделирования и
экспериментальных данных показало хорошую сходимость. Предложенная методика численного моделирования взаимодействия СОЖ с деталью и инструментом может быть использована при модернизации технологических процессов, используемых при обработке глубоких отверстий на токарных станках.
Ключевые слова: турбулентный поток, диссипативные системы, математическая модель, численное моделирование, смазочно-охлаждающая жидкость, несжимаемая жидкость, борштанга, глубокое отверстие.
1. Введение
При обработке металлов резанием на токарных станках приходится решать задачи повышения точности обработки глубоких отверстий путем поддержания температуры в оптимальном режиме в зоне резания и уменьшения силы трения между режущим инструментом и деталью [1]. В большинстве случаев это достигается с помощью СОЖ, которая не только охлаждает инструмент и деталь, но и способствует отводу стружки из зоны резания. Эта задача многократно усложняется при растачивании глубоких отверстий, так как характеристики потока СОЖ, взаимодействующего с инструментом (борштангой) и деталью существенно меняются вдоль потока. Технически регулировать процесс теплообмена в зоне резания можно лишь за счет изменения скорости и температуры СОЖ на входе.
Качественная токарная обработка металла в большинстве случаях невозможна без применения СОЖ, поэтому исследованию ее свойств уделяется большое внимание [2-4]. Так в [2] дана оценка охлаждающим воздействиям ряда смазочно-охлаждающих жидкостей при обработке металлов резанием, в [3] приведена оценка влияния СОЖ на точность и чистоту обработки, а в [4] ее соответствие обрабатываемому материалу. В работах [5-7] представлены методы повышения точности обработки глубоких отверстий и основные положения, связанные с разработкой технологических процессов глубокого сверления отверстий малых диаметров с подводом СОЖ. В [8,9] приведены основные причины появления дефектов и методы повышения точности обработки глубоких отверстий.
Известно [1-7], что технически регулировать процесс теплообмена и удаление стружки в зоне резания можно лишь за счет оптимального способа подвода на входе и отвода СОЖ в системе. Из-за сложности объекта исследования установление связи между характеристиками потока в разных областях возможно только с применением численного моделирования, которое начинается с создания физической модели. Для построения адекватной физической модели рассматриваемого технологического процесса нужно определить физико-технологические характеристики всех взаимодействующих объектов: СОЖ, деталь и режущий инструмент.
В работе приведены математическая модель, описывающая распределения температуры, скорости и давления, обрабатываемой детали и режущего инструмента при теплообмене, представляющая собой систему дифференциальных уравнений в частных производных (уравнения Навье - Стокса, уравнение энергии и уравнения турбулентности) с граничными условиями и методика численного моделирования взаимодействия смазочно-охлаждающей жидкости с обрабатываемой деталью и инструментом при растачивании глубоких отверстий. Следует отметить, что численное моделирование распределения давления, скорости, температуры в потоке СОЖ, температуры детали и режущего инструмента проведено с использованием программного комплекса FlowVision HPC и аналитических результатов, полученных в настоящей работе. В результате, определены: давление, скорость, температура в потоке СОЖ, температуры обрабатываемой детали и режущего инструмента при теплообмене, а также интегральные гидродинамические характеристики. Кроме того, разработаны расчетные модели, позволяющие проводить вычислительные эксперименты при
различных входных данных и способах подвода и отвода СОЖ, не прибегая к дорогостоящим натурным экспериментам.
2. Математическая модель
Рассмотрим эквивалентную механическую модель токарного станка, предназначенного для растачивания глубоких отверстий (рис.1) [10]. На рисунке изображены элементы рабочего инструмента - стебель борштанги с расточной головкой, а также элементы установки - суппорт, шпиндельная бабка, люнеты и станина. Обрабатываемая деталь крепится к патрону, вращающемуся вместе со шпинделем с п (об/мин), а борштанга - к суппорту, перемещающемуся поступательно с подачей Я (мм/мин).
Рис. 1. Эквивалентная механическая модель токарного станка. Обозначения: 1 - стебель борштанги, 2 - расточная головка, 3 - обрабатываемая деталь, 4 - суппорт, 5 - шпиндельная
бабка, 6 - люнеты, 7 - станина
При растачивании глубоких отверстий инструмент имеет большую длину, поэтому его выполняют составным по длине, используя разъемное соединение двух основных частей -режущего и вспомогательного инструментов [10, 11]. Режущий инструмент выполняется в виде так называемой расточной головки, устанавливаемой на конце вспомогательного инструмента - стебля борштанги.
Для численного моделирования турбулентного взаимодействия СОЖ с деталью и инструментом при составлении математической модели, описывающей распределения температуры, скорости и давления, приняты следующие предположения:
1.СОЖ интерпретируется, как вязкая несжимаемая жидкость.
2. Течение вязких жидкостей может быть ламинарным или турбулентным [12-14]. Переход с ламинарного режима на турбулентный происходит при повышении скорости потока при достижении определенных значений числа Рейнольдса. Проведенные оценки для нашего случая взаимодействия СОЖ с деталью и инструментом при обработке глубоких отверстий показывают на турбулентный характер течения.
3. Характеристики турбулентного течения СОЖ можно разделить на две составляющие: осредненную и пульсирующую. При осреднении по времени (большому по сравнению с периодом пульсации) остается только «медленно» меняющаяся составляющая, которая и представляет главный интерес для практических приложений.
Так как уравнение Навье-Стокса, описывающее ламинарное движение, нелинейное, то при переходе к осредненным переменным, получаем уравнения Рейнольдса в виде
дтц
Щ дУ; д1 дхк рдХ] рдх]■
Ч х дУ1
-риж
(1)
где VI - декартовая компонента скорости, д^ - ускорение свободного падения, р - плотность жидкости.
Уравнение неразрывности в силу его линейности остается неизменным и для осредненных величин, а уравнение Навье-Стокса изменилось, в нем кроме тензора напряжений
о ^ (дУгдУЛ (2)
где р - давление, д - динамическая вязкость жидкости, появляется тензор турбулентных напряжений (тензор напряжений Рейнольдса) т^. Тензор турбулентных напряжений т^ выражается через корреляцию второго порядка пульсации скорости и\, которые являются мгновенными добавками к осредненной скорости
*и = -р*Щ, (3)
ту - симметричный тензор второго ранга содержит 6 компонент, которые надо находить вместе с ^ и р, в силу чего система уравнений Рейнольдса становится незамкнутой. Простейший способ замыкания уравнений Рейнольдса заключается во введении турбулентной вязкости
_ дуЛ
+ (4)
Турбулентная вязкость не является константой. Она вычисляется в рамках стандартной к — £ модели по формуле [15]
к2 к2 (5)
^т = СцР—>$т = Сц —
где т9т - турбулентная кинематическая вязкость.
Удельная энергия турбулентных пульсаций к и скорость ее диссипации £ подчиняются уравнениям переноса[15]
дк дк _ д д1 1 дХ[ дХ[
( дк
^+—)-;—
V ок) дх1.
+ — (6)
де де д
д1 1 дХ[ дХ[
/ дт\ де' ( аЕ) дх(.
£ £2 + Ст р ~ С], Спр—;—, 1Е к к 2Е к
(7)
где
дУ1 дЦ
Си = 0.09, С1£ = 1.44, С2е = 1.92, ак = 1,а£ = 1.3.
(8) (9)
Распределение температуры в турбулентном потоке определяется уравнением
дТ дТ _ д
д1 1 дХ[ дХ[
(а + ат) —
дТ
дхь
(10)
где а = Х/Ср - коэффициент температуропроводности среды, X - коэффициент теплопроводности, с размерностью , С- удельная теплоемкость среды с размерностью
ат- турбулентная температуропроводность, которая для инженерных расчетов
определяется соотношением ат = дт/Ргт. Здесь Ргт - число Прандтля (близкое к единице).
Решения системы уравнений (1), (6), (7), (10) ищутся при численном моделировании, предполагая, что они удовлетворяют следующим граничным условиям:
1. На входе во внешний канал борштанги задается нормальная составляющая массовой скорости течения СОЖ (рУсож), температура СОЖ (ТП и начальная степень турбулизации потока рабочего вещества (к — г)
_ Я _ (0.03Усож)2
рУп1Ь = рУсож =^~>Т1Ь = Ты; к1Ь = 5 , ^Ъ =
С^(к1ъ)3/2
Ъ
0.001
(11)
2. Температура СОЖ на входе считается равной 30°С .
3. На выходе из внутреннего канала стебля борштанги задается условие нулевых возмущений по статическому давлению
Р1ъ = 0;
йУп
йп
й Т й к й
= 0; — = 0; — = 0; —
Ъ й п Ъ й п Ъ й п
= 0.
(12)
Ъ
4. На всех стенках конструкции задаются: условие непротекания Уп1ъ = 0 с турбулентным пограничным слоем, характеризующимся логарифмическим законом изменения касательной компоненты скорости У^ъ, а также к1ъ и е1ъ - через пристеночные функции с граничными условиями
йк йп
йе
Ъ й п
= 0.
(13)
Ъ
3. Аналитическая оценка теплообмена между подаваемой СОЖ и обрабатываемой деталью
Известно [16], что теплопередача от жидкости к твердому телу определяется явлениями конвекции и теплопроводности. Каждое из этих явлений зависит от формы стенки, отделяющей горячую среду (жидкость) от холодной (твердого тела) (рис. 2).
Для вывода аналитического выражения оценки теплообмена между подаваемой СОЖ и обрабатываемой деталью рассмотрим задачу, в которой стенка обрабатываемой детали имеет цилиндрическую форму с внутренним радиусом г3 и внешним радиусом г4. Подача СОЖ осуществляется по внешнему каналу, находящемуся между борштангой и обрабатываемой деталью. С целью снижения энергетических затрат СОЖ подается в нагретом (до температуры 25-30°) состоянии [10]. Поэтому внутри обрабатываемой детали находится горячая среда с температурой Тп. Снаружи поддерживается холодная среда с температурой Тм,. Температура внутренней поверхности обрабатываемой детали Т изменяется под действием горячей и холодной среды. Значение данной температуры и является искомым.
Рис. 2. Теплопередача через цилиндрическую стенку детали
Изменение температуры между горячей и холодной средой происходит только в радиальном направлении г. Выделив внутри стенки кольцевой слой радиусом г и толщиной йг, ограниченный изотермическими поверхностями, получим, что, согласно закону Фурье, количество теплоты, проходящее в единицу времени через этот слой, равно
Q = -ЛF?í = -2Anrl^ , (14)
йг йг
где F- поверхность обрабатываемой детали, Л - коэффициент теплопроводности, I - длина цилиндра, t - температура обрабатываемой детали. Из уравнения (14) получим t = — + С. Подставляя значения переменных на границах стенки (при г = г3,1 = и при
т = г4,Ь = 12 и исключая постоянную С , находим следующую расчетную формулу
2Лп1
(15)
¿3
Количество теплоты, переданное в единицу времени через стенку трубы, прямо пропорционально коэффициенту теплопроводности Я, длине I и температурному напору № = — и обратно пропорционально натуральному логарифму отношения внешнего диаметра детали й4 к внутреннему й3. Количество теплоты, проходящее через стенку трубы, может быть отнесено к единице длины I . При этом расчетная формула принимает вид
_ 2ХпМ (16)
Ч1 = - = ЦГ-
¿3
При установившемся тепловом состоянии системы количество теплоты, отданное горячей и воспринятое холодной средой, одно и тоже. Следовательно, на основе формулы (16) для определения температуры внутренней поверхности Л обрабатываемой детали Т можно написать соотношение
2пЛж 2пЛгт (17)
¿2 й3
где й2 - внешний диаметр борштанги, й3- внутренний диаметр обрабатываемой детали, й4 -внешний диаметр обрабатываемой детали, Лж- коэффициент теплопроводности внутренней (горячей) среды и Лгт - коэффициент теплопроводности стенки обрабатываемой детали. Из полученного соотношения температура внутренней поверхности обрабатываемой детали равна
ЛжТ1п1 ~Т + ЛстТ\н1п~Г „ ч
Т =-Ц-(18)
Л^Ы — + ЛГт1п —
3
2
Для получения аналитического выражения для оценки температуры на связанной границе используется формула (18). Данная формула предназначается для расчета температуры на внутренней поверхности обрабатываемой детали в случае, когда деталь находится между горячей и холодной средой. В рассматриваемой задаче горячая среда находится как между цилиндрическими поверхностями стебля борштанги и обрабатываемой детали, так и между цилиндрическими поверхностями расточной головки и обрабатываемой детали. В первом случае аналитическое решение примет вид
ЛжТ1п1 ~Г + ЛпТы\п ~г
т - аз_(19)
1 апаШ1с_1 ~ а а > 4 '
Лж1п—+ Лгт1п —
а во втором случае
ЛжТШ ~Г + ЛпТы\п ~Т
т __¿3_(20)
Т апа1Шс_2 = а а . у !
Лж1п—+ Лг-Лп —
При соотношении длин стебля борштанги и расточной головки (1:5) выражение аналитической оценки температуры на сопряженной стенке имеет вид
5 15 ЛжТ1ш ~т + Лп^1п~ Т — _гг _|__т__¿3_
1 апаШгс ~ (-1апаШ1с_1^(-1апа1Шс_2=^ ¿л й
б^апаиис 1 ■ /■ ~ апаимс 2= г п* й,
_ 6 _2 6 а ^nй± + лcт1nйз
Чз ст ¿2
+ ЛгтТу^1п^Г з_
Лж1nй4 + ЛгT1nйз (21)
1 ЛжТШ ~Г + ЛгтТлл,\П~Т
+1 ¿3
6 Л^п^ + Л^Ы^
Ж ¿3 с
з
2
з
6
4. Численное моделирование
Численное моделирование проводилось с помощью программного комплекса FlowVision НРС, в котором применяется метод конечных объёмов.
Метод конечных объемов предполагает интегрирование уравнений движения жидкости и переноса скалярных величин по объемам ячеек расчетной сетки. По теореме Гаусса для произвольной векторной или тензорной величины Б
/
(V^F)dü=I(Fi-ni)-ASi, (22)
где П- объем ячейки, АЯ^- площадь /-ой грани ячейки. При интегрировании решаемых уравнений в ячейке производится суммирование потоков массы, импульса, энергии и турбулентных величин, вычисленных на гранях ячеек. Поскольку каждая грань разделяет две соседние ячейки, соответствующий поток входит в дискретные уравнения для обеих ячеек. Этим обеспечивается консервативность массы, импульса, энергии и других искомых величин в расчетной области [13].
Интегрирование системы уравнений (1),(6),(7),(10) в FlowVision НРС проводится методом расщепления по физическим процессам. Для интегрирования уравнений по времени используется конечно-разностная схема 2-го порядка точности, тип схемы - неявная.
Для иллюстрации метода расщепления по физическим процессам опишем общую расчетную схему интегрирования уравнений на одном шаге по времени.
Рассмотрим уравнения Эйлера (вязкими членами пренебрегаем). Проинтегрируем уравнение неразрывности и уравнение импульсов по расчетной ячейке - многограннику
рП+1 — рп I (23)
--^- + -1(рп+1Уп+1-п-А5) = 0.
рп+1уп+1 — „пуп 1 (24)
---—--+ -!(рпУп ■п-АБ)Уп+1 = —Чрп+1 + рп¥п
Перепишем уравнение (24) в виде
„П + 1уП+1 — ХП + 1 уП + 1 ХП+1 уП+1 — 0пуп 1
-----+ ---—--+ ■ п ■△ Б)^п+1 = (25)
= -ЧрП+1 + рп¥п - Чрп + Чрп
Разделим (расщепим) его на два уравнения:
-1 уп+1 — рпуп 1 ^
-----+ -!(рпУп ■ п ■ АБ)^1 = —Чрп + рпТп, (26)
0п+1уп+1 — дп+1 УП+1 (27)
г г „ п+1 , уу п V '
Vpn+L + Vpn.
т
Знаком «тильда» обозначим промежуточные решения. Далее выразим массовую скорость на п+1 временном слое из уравнения (27)
рп+^уп+1 = РП+1 УП+1 + Т(-Чрп+1 + Чрп). (28)
Подставляя выражение (28) в дискретное уравнение неразрывности (20), получим уравнение Пуассона для давления
рП+1 — рП 1
--— + -1(рп+1Уп+1 • п • Д5) = т(Арп+1 - Арп).
тП
(29)
Далее выполняются следующие действия:
1.Находится промежуточное (неконсервативное) поле скорости \п+1 из уравнения (26).
2.Решается уравнение для давления (29) и по давлению определяется плотность рп+1.
3.Рассчитывается окончательное (консервативное) поле скорости \п+1, используя уравнение (28).
4.Решается уравнение энергии, в котором используется поле \п+г.
5.На п+1 временном слое рассчитывается температура
Ьп+1 = к(Тп+1). (30)
б.Решаются уравнения стандартной к- е модели турбулентности, в которых используется поле Уп+1.
Описанный выше численный алгоритм, реализуемый в НРС, был апробирован
на большом числе прикладных задач и показал высокую устойчивость, стабильную сеточную сходимость и согласованность численных результатов с экспериментальными данными [13].
Исследование сходимости по сетке для задачи, представленной в статье, описывается ниже в п. 3.2.
4.1. Создание расчетных моделей
Для трехмерного моделирования в программном комплексе FlowVision НРС течения СОЖ и ее теплообмена с деталью и борштангой в системе SolidWorks созданы две геометрии объема проточной части, определяемые тремя элементами: стеблем борштанги, расточной головкой и обрабатываемой деталью при внешнем и внутреннем подводе СОЖ [17]. Для численного моделирования выбрана расточная головка фирмы Бо1ек (Германия) одностороннего резания с определенностью базирования [11], предназначенная для растачивания отверстий диаметром от 45 мм до 250 мм.
Геометрическая модель объема проточной части с тремя подобластями ("деталь", "инструмент" и "течение") для численного моделирования при внешнем подводе СОЖ представлена на Рис. 3.
к
!
1
Рис. 3. Геометрическая модель объема проточной части при внешнем подводе СОЖ
Рис. 4. Расточная головка фирмы Botek (Германия) одностороннего резания с определенностью базирования для растачивания глубоких отверстий
Ниже представлены физические и геометрические параметры модели рис.4 , используемые для численного моделирования. Физические и геометрические параметры рабочего инструмента и обрабатываемой детали при обработке глубоких отверстий взяты из работ [10, 11]. В качестве СОЖ выбран продукт ЕСОСЦТ HFN 5 LE [17]. Рабочий инструмент и обрабатываемая деталь изготовлены из стали 45. Параметры борштанги:
- длина борштанги Ь = 1.3 м,
- внутренний радиус борштанги г1 = 0.03 м,
- внешний радиус борштанги г2 = 0.06 м,
- теплопроводность борштанги А^^ = 46Вт ■м-1 К-1 м Параметры обрабатываемой детали:
- внутренний радиус обрабатываемой детали г3 = 0.0855 м,
- внешний радиус обрабатываемой детали г4 = 0.1155 м,
- теплопроводность обрабатываемой детали ■ ЛаеС = 46 Вт • м-1 • К-1 Параметры расточной головки (рис.3, рис. 4):
- длина расточной головки I = 0.24 м,
- внутренний радиус расточной головки г5 = 0.03 м,
- внешний радиус расточной головки г6 = 0.08 м,
- длина выступа направляющих элементов п = 0.055 м,
- длина направляющих элементов к = 0.16 м,
- длина выступа режущей пластины т = 0.055 м,
- теплопроводность расточной головки Л1Пз1г = 46 Вт • м-1 • К-1 Параметры смазочно-охлаждающей жидкости:
- плотность СОЖ рцц = 830 кг • м-3,
- вязкость СОЖ = 0.00332 кг • м-1 • с-1,
- теплопроводность СОЖ Лцч = 8.04104 Вт • м-1 • К-1,
- теплоемкость СОЖ СрП({ = 2422Дж • кг-1 • К-1. Входные параметры:
- массовая скорость СОЖ на входе У^ = 415кг • м-2 • с-1,
- скорость вращения обрабатываемой детали шае1 = 62.8 рад/с,
- температура подаваемой на входе СОЖ Т^ = 30°С,
- температура обрабатываемой детали Тйе1 = 22°°,
- температура рабочего инструмента (борштанги и расточной головки) Т1Пз1г = 22°, температура в зоне резания Тг = 800°.
Для численного моделирования подключались следующие физические процессы программного комплекса: Движение = Ньютоновская жидкость, Теплоперенос = Конвекция и теплопроводность, Турбулентность = Стандартная к — £ модель турбулентности.
4.2. Исследование сходимости по сетке
В ходе проведения численного моделирования одной из сложнейших задач оказывается поиск оптимальной расчетной сетки, позволяющей при минимальной размерности задачи отследить тонкие физические явления. В исследуемой задаче к таким явлениям относятся вихреобразование, распространение тепловых потоков из области резания и образование пограничных слоев вблизи стенок конструкции [19].
Поиск оптимальных методов построения сетки осуществлен путем проведения исследования сходимости. Исследование сходимости -- один из принципов оценки точности получаемого решения. Он заключается в проведении серии расчетов одной и той же задачи на сетке, которая последовательно сгущается во всей области расчета. При уменьшении расчетных ячеек точность решения исходных уравнений увеличивается пропорционально Л-п, где Л - размер расчетной ячейки, п - порядок аппроксимации расчетной схемы. Моделируемые параметры задачи (например, тепловой поток, проходящий внутри детали, температура СОЖ на выходе из внутреннего канала борштанги) при этом сходятся к некоторому значению, соответствующему бесконечно мелкой сетке.
На рис. 5 приведены графики изменения температуры СОЖ на выходе из внутреннего канала борштанги с течением времени при разных уровнях адаптации. В основе построения расчетных моделей использована равномерная начальная расчетная сетка с размером ячейки 0.8^0.6x0.6 см. Адаптация расчетной сетки проведена последовательно с 1-го до 3-го уровня на внутренней поверхности обрабатываемой детали и поверхностях инструмента. При этом стационарное решение, полученное для предыдущего уровня адаптации, было использовано в качестве начального приближения.
Относительное изменение значения температуры СОЖ на выходе при переходе с 1-го на 2-й уровень адаптации составило 10.2%, со 2-го на 3-й - 2%. При этом число ячеек начальной расчетной сетки с 1-ым уровнем адаптации составило 1 млн, со 2-ым - 3.3 млн, с 3-им - 12 млн.
т.ч:
1
}
—-1 уровжь адоптации —2 урогемь адаптации Я уровень ПТЛЦИИ Т-1-1—J
С ДО 40 Ы5 30 ДО0 дго Д40 1« 190
1,1
Рис. 5. Изменение значения температуры СОЖ на выходе из внутреннего канала борштанги с течением времени при разных уровнях адаптации
Из приведенного рисунка и оценки относительного изменения значения выходного параметра следует, что исследование сходимости может быть завершено. Таким образом, наиболее точный результат получен при адаптации 3-ого уровня на поверхностях рабочего инструмента и обрабатываемой детали.
Для меньшего ущерба масштабируемости задачи (с точки зрения распараллеливания вычислений) принято решение проводить итоговые расчеты на начальной сетке с размером ячейки в 2 раза меньше и со 2-м максимальным уровнем адаптации на поверхностях конструкции. Таким
образом, минимальный размер расчетной ячейки, как и в третьем варианте задачи, составит 0.1x0.075x0.075 см.
4.3. Результаты численного моделирования
В результате проведения вычислительных экспериментов на выбранной сетке получены картины распределения поля температур вблизи зоны резания, распределение поля скоростей и давления, гидродинамические характеристики, необходимые для проверки корректности численного решения. Критерием установления решения выбрано изменение значения температуры СОЖ на выходе из отводного канала борштанги.
На рисунке 6 представлены способы подвода СОЖ и отвода стружки при растачивании глубоких отверстий, используемые при проведении вычислительных экспериментов. На рисунке 6а приведен способ наружного подвода СОЖ и внутреннего отвода стружки. СОЖ подводится по наружному каналу Н между поверхностью инструмента 2 и стенками отверстия в заготовке, а отводится вместе со стружкой через окно С и далее по внутреннему каналу В в инструменте. Сквозное отверстие в заготовке закрывают заглушкой 4 с уплотнением 3, поджимаемым плотно к торцу заготовки. Перед включением рабочей подачи инструмента полость Д в заготовке полностью заполняется СОЖ. На рисунке 6б приведен внутренний способ подвода СОЖ. СОЖ подводится через внутренний канал В в инструменте, через окно С попадает в зону резания и выводится вместе со стружкой через необработанное отверстие в заготовке. Более детальное описание режущей головки, приведенной на рисунке, дано в [11].
Наружный подвод СОЖ (а) Внутренний подвод СОЖ (б)
Рис. 6. Схемы подвода и отвода СОЖ
На рисунке 7 представлены в виде цветовой заливки картины рассчитанного распределения давления СОЖ в подобласти течения в плоскости ХУ и распределение давления в подобласти течения в плоскости У2 (сечение) соответственно.
Наружный подвод СОЖ (а)
■
ММ
■I
12150
10500
8850
щя
ия
■
Внутренний подвод СОЖ (б)
Рис. 7. (а) Распределение давления (Па) в подобласти течения в плоскости ХУ, (б) распределение давления в подобласти течения в плоскости У2
Анализ результатов расчетов (рис. 7 (а)) показал, что наибольше давление наблюдается между обрабатываемой деталью и внешней поверхностью стебля борштанги. Давление СОЖ падает в зоне растачивания и в проточной части борштанги при приближении к свободному выходу. На рис. 7 (б) наибольшее давление наблюдается в проточной части борштанги. В зоне растачивания и при приближении к свободному выходу через отверстие детали давление уменьшается. Наименьшее давление наблюдается между обработанной поверхностью детали и наружной поверхностью стебля борштанги. Наибольшее давление в целом наблюдается при наружном подводе СОЖ. На рисунке 8 представлено рассчитанное распределение векторного поля скоростей потока СОЖ в плоскости XY при наружном подводе СОЖ.
Рис. 8. Распределение векторного поля скоростей (м/с) потока СОЖ в плоскости XY
Из рисунка видно, что поток ускоряется в зоне расточной головки и замедляется, попадая в проточную часть обрабатываемой детали. На рисунке 9 представлено распределение температуры в расточной головке в плоскости ХУ. В зоне резания наблюдается максимальная температура. При внутреннем подводе СОЖ расточная головка нагревается сильнее, чем при внешнем подводе.
Наружный подвод СОЖ (а)
Внутренний подвод СОЖ (б)
Рис. 9. Распределение температуры (°С) в расточной головке в плоскости XY
При растачивании глубоких отверстий могут быть использованы оба способа обработки, но большее предпочтение отдается внутреннему подводу СОЖ, так как тепло из зоны резания отводится лучше, в результате расточная головка нагревается меньше.
Из представленных рисунков видно, что при течении по внешнему каналу за счет теплообмена с обрабатываемой деталью и инструментом происходит незначительное охлаждение СОЖ. Область максимальных температур находится в зоне резания. Наибольшее изменение температур наблюдается в расточной головке на расстоянии 0.04 м от зоны резания. Далее температура падает и становится близкой к начальной температуре инструмента, равной 22°. Это обусловлено тем, что в процессе растачивания глубоких отверстий СОЖ отводит значительное количество теплоты из зоны резания, а значит, её использование является неотъемлемой частью технологического процесса.
Векторное поле скоростей СОЖ в области расточной головки, приведенное на рисунке 8, свидетельствует о нескольких характерных зонах течения. Во внешнем канале, образованном наружной поверхностью инструмента и поверхностью обработанного отверстия в заготовке, образуется зона ламинарного течения. Значительное изменение скорости наблюдается при попадании СОЖ в область между направляющими элементами и обрабатываемой деталью: поток ускоряется. Далее при попадании в проточную часть обрабатываемой детали происходит замедление движения жидкости. Течение СОЖ в проточной части обрабатываемой детали характеризуется наличием развитых вихревых зон.
Проверка корректности численного результата была проведена с использованием формулы (21), представляющей собой теоретическую оценку значения температуры на внутренней поверхности обрабатываемой детали. Значение температуры поверхности детали с обрабатываемой стороны, подсчитанное по приведенной формуле, составляет 25.6° . В результате численного моделирования получено значение температуры Тру = 28°. Таким образом, относительная погрешность численного решения, подсчитанная по формуле
составляет 9%.
Решение, полученное численным путем, на 2.4 градуса выше теоретического. Результат объясняется тем, что в теоретической оценке не учитывается передача тепла от трех направляющих элементов, расположенных на расточной головке. Направляющие элементы являются базирующими и находятся на близком расстоянии (1 мм) от внутренней поверхности обрабатываемой детали. В процессе растачивания данные элементы расточной головки также нагреваются и передают тепло обрабатываемой детали.
8 =
100%
(30)
Температура циркулирующей СОЖ, замеряемая при выходе из отводного канала перед входом в бак, в соответствии с экспериментальными данными должна составлять от 30° до 50° [10]. Численное значение температуры СОЖ на выходе, полученное в программном комплексе, равняется 48.9°С, что удовлетворяет рекомендуемому специалистами диапазону.
5. Заключение
• В работе приведена методика численного моделирования взаимодействия смазочно-охлаждающей жидкости (СОЖ) с обрабатываемой деталью и инструментом при обработке глубоких отверстий.
• Построены расчетные модели, включающие в себя: систему из расчетной области, дифференциальных уравнений, описывающих физические процессы, протекающие в исследуемой системе, граничных условий и сетки.
• Описан процесс теплообмена между подаваемой СОЖ и обрабатываемой деталью.
• Проведены теоретические оценки для значения температуры внутренней поверхности обрабатываемой детали в зависимости от гидродинамических характеристик.
• Результаты численного моделирования иллюстрируются картинами течения СОЖ и теплообмена).
• Корректность численного решения установлена на основе сравнения теоретических результатов и экспериментальных данных.
• Анализ результатов численного моделирования показал, что расчетные модели адекватно и достаточно точно отображают процессы взаимодействия смазочно-охлаждающей жидкости (СОЖ) с обрабатываемой деталью и инструментом при обработке глубоких отверстий.
Разработанные расчетные модели позволяют многократно проводить вычислительные эксперименты при различных входных данных, не прибегая к дорогостоящим натурным экспериментам.
6. Благодарности
Работа выполнена при поддержке Министерства науки и высшего образования Российской Федерации, соглашение Ко 075-15-2020-808.
Список литературы
[1] Бобров, В.Ф. Основы теории резания металлов. М.: Машиностроение, 1975. 344 с.
[2] Кисель, А.Г., Пуртов, Е.Д., Дейлова, А.В., Кочура, Н.Н. Оценка охлаждающих свойств смазочно-охлаждающих жидкостей / Омский научный вестник. 2017. № 1(151). С. 2729.
[3] Кисель, А.Г., Реченко, Д.С., Титов, Ю.В., Пуртов, Е.Д., Петров, И.В. Подбор смазочно-охлаждающей жидкости для чистовой обработки/ Системы. Методы. Технологии. 2015. № 3(27). С. 39-43.
[4] Кисель, А.Г., Ражковский, А.А., Реченко, Д.С., Попов, А.Ю. Повышение точности токарной обработки за счет применения смазочно-охлаждающих жидкостей / Технология машиностроения. 2014. № 2. С. 18-20.
[5] Немцев, Б.А., Яковлев, П.Д., Яковлев С.П. Технология глубокого сверления отверстий малых диаметров с наружным подводом СОЖ // Металлообработка. 2015. № 4(88). С. 19-24.
[6] Biermann, D., Sachrow, A., Wohlgemuth, K. Simulation of the BTA deep-hole drilling process // Prodaction Engineering Res. Devel. 2009. P. 339-346.
[7] Weinert, K., Weihs, C., Webber, Oliver, Raabe, Nils. Varying bending eigenfrequencies in BTA deep hole drilling: mechanical modeling using statistical parameter estimation // Production Engineering Res. Devel. 2007. P. 127-134.
[8] Горелова, А.Ю., Плешаков, А.А., Кристаль, М.Г. Методы повышения точности обработки глубоких отверстий // Известия Тульского государственного университета. 2013. Вып. 7. Ч. 2. С. 363-370.
[9] Novakova, T., Jackson, M. J. Chatter problems in microand macrocutting operations, existing models, and influential parameters - a review // Manuf. Technol. 2010. P. 597-620.
[10] Уткин, Н.Ф. Обработка глубоких отверстий, Л.: Машиностроение, 1988. 269 с.
[11] Кожевников, Д.В., Гречишников, В.А., Кирсанов, С.В., Кокарев, В.И., Схиртладзе, А.Г. — Под редакцией С.В. Кирсанова. Режущий инструмент. 3-е. изд. — М.: Машиностроение, 2007. — 528 с.
[12] Комаров, В.Н., Грезина, А.В., Артемьева, С.А. Моделирование процесса теплообмена при растачивании глубоких отверстий //Вестник Нижегородского университета им. Н И. Лобачевского. 2011. № 4-2. С. 87-91.
[13] FlowVision HPC. Руководство пользователя. Версия 3.08.04. ООО ТЕСИС. М. 19992013.
[14] Бэтчелор, Дж. Введение в динамику жидкости. М.: Мир, 1973. — 760 с.
[15] Мазо, А.Б. Модель турбулентных течений несжимаемой жидкости - КГУ, 2007,106с.
[16] Михеев, М.А., Михеева, И.М. Основы теплопередачи. Изд. 2-е., стереотип. М., «Энергия», 1977, 344 с.
[17] Мюррей, Д. SolidWorks. Издательство "ЛОРИ", 2009 г., 604 с.
[18] http://uniks42.ru/upload/products/6670/ECOCUT HFN 5 LE ru.pdf
[19] Роуч, П. Вычислительная гидродинамика. М.: Мир, 1980 г., 616 с.
Numerical - analytical study of the turbulent flow of lubricoolant when interacting with a part and a tool in the process of machining
deep holes
LA. Igumnov1' , A.V. Grezina2 , V.S. Metrikin2' , AG. Panasenko2
e-mail: *igumnov@mech.unn.ru aleksandra- grezina@yandex .ru v.s.metrikin@mail.ru a. g.panasenko@yandex.ru
Abstract. The paper presents an approach to modeling of the interaction of lubricoolant with the workpiece and the tool in the process of deep hole boring. The modeling is based on the numerical solving of the system of nonlinear partial differential equations describing the turbulent flow of the lubricoolant and the transfer of thermal energy in the channels of the part and tool. The fluid is assumed to be incompressible, and the standard turbulence k -s model is used to close the system of differential equations. SolidWorks is used to model the geometry of the flow volume. The distribution of pressure' velocities' the temperature of the lubricoolant flow as well as the temperature of the workpiece and the cutting tool is numerically simulated by using FlowVision HPC. The dynamics of the temperature balance between the lubricoolant, tool and workpiece is analyzed. The integral hydrodynamic characteristics of the flow are given in the analytical form. The developed numerical models allow us to carry out numerical experiments for various combinations of in- and out- flow of the lubricoolant. The comparison between numerical results and experimental data testifies to their good convergence. The proposed approach can be applied to the optimization of technological processes of machining deep holes on lathes.
Keywords: turbulent flow, dissipative system, incompressible fluid, boring bar, deep hole, mathematical model, numerical modeling, lubricoolant.
Acknowledgements The work was supported by the Ministry of Science and Higher Education of the Russian Federation, agreement No 075-15-2020-808.