Моделирование теплофизической задачи, возникающей при обработке глубокого отверстия

УДК 621.787.4

Л.Р. Милованова, Я.И. Барац, М.В. Стекольников, В.С. Разуменко МОДЕЛИРОВАНИЕ ТЕПЛОФИЗИЧЕСКОЙ ЗАДАЧИ, ВОЗНИКАЮЩЕЙ ПРИ ОБРАБОТКЕ ГЛУБОКОГО ОТВЕРСТИЯ

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

Тепловой источник, глубокое сверление, теплопроводность, метод источников тепла, принцип пространственно-временного соответствия.

L.R. Milovanova, Ya.I. Baratz, M.V. Stekolnikov, V.S. Razumenko MODELING OF THERMOPHYSICAL PROBLEMS ARISING IN THE DEEP-HOLES PROCESSING

The axisymmetric thermalphysic problem of deep-hole drilling has been solved in this artical, with the help of mathematical model calculation of temperature based on heat source method with spatio-temporal correspondence principal and local influence principal. The temperature in contact zone ‘tool-component’, and also the temperature field of work stock is calculated. Adequacy of the offered mathematical model is proved.

Heat source, deep-hole, thermo-conduction, heat sources method, principle space-temporary correspondence.

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

Отверстие считается глубоким, если отношение его длины к диаметру более 10. Однако на практике уже при обработке отверстия на глубину превышающую три диаметра наблюдаются трудности с выводом стружки, снижение стойкости сверла и производительности обработки. Технологии сверления глубоких отверстий отличаются от обработки обычных отверстий и требуют применения специфического оборудования и технологического оснащения, а именно, специальных режущего, вспомогательного, измерительного инструментов и приспособлений. При глубоком сверлении к инструменту подается смазочноохлаждающая жидкость (СОЖ), это также определяет специфику не только процесса, но и применяемого оборудования и оснастки [1].

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

Изображение теплового процесса позволяет использовать фундаментальное решение уравнения теплопроводности для осесимметричных теплофизических задач, например, соотношение для расчета температурного поля в неограниченном теле с мгновенным источником мощностью д Дж/м в виде окружности радиусом г’, действующего в плоскости г’=0 в момент времени t=0 [2]:

а = Х/(ср) - коэффициент температуропроводности; ср - объемная теплоемкость; X - коэффициент теплопроводности.

В соответствии с правилом изображения тепловой режим в теле с отверстием радиуса г0, на адиабатической поверхности которого действует кольцевой источник мощностью Q = 2пг0 д , будет адекватен тепловому режиму в неограниченном теле в области г > г0, если

в последнем на окружности г = г0 разместить кольцевой источник удвоенной мощности

Q0 = 2Q, а в точке г = 0 поместить сток теплоты мощностью Q с = —Q.

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

В соответствии с принципом местного влияния [3] допустимо рассматривать температурное поле в детали как сумму двух полей: общего поля вдали от источника и местного поля непосредственно в пространстве источника. При этом, определяя общее поле, можно всемерно схематизировать источник теплоты, учитывая в то же время более точно условия теплоотвода. При расчете же местного поля условия теплоотвода могут быть схематизированы, в то время как распределения интенсивности местного источника должны быть наиболее полно учтены.

В связи с этим, рассчитывая местное температурное поле в детали, можно допустить следующую схематизацию. Обрабатываемую деталь условно представить неограниченным телом с цилиндрическим отверстием радиуса г0. Далее, имея в виду, что окружная скорость источника на несколько порядков больше скорости движения подачи, то по значению критерия Пекле Ре = 2УЯ / а > 10 , такой быстродвижущийся местный источник в пределах одного оборота детали можно считать пространственным мгновенным кольцевым [4] с распределением интенсивности в направлении радиуса г и координаты г по нормальному закону Гаусса:

В соответствии с принципом пространственно-временного соответствия [5], чтобы перейти от решения, описывающего температурное поле в неограниченном теле с мгновенным кольцевым источником (1) к соответствующему решению с нормально-тороидальным источником (2) достаточно в соотношении (1) к временной координате і добавить постоянную времени

где Я — характерный размер источника тепла, м; к — коэффициент сосредоточенности интенсивности источника.

где Q = 2пг'д ; 10(гг'/(2at)) - функция Бесселя нулевого порядка от мнимого аргумента;

і0 = Я2 /4-ак ,

Известно, что для малых промежутков времени, то есть для больших значений аргумента функции Бесселя, в соотношениях (1) допустима замена

гг' гг' /I

~ ехр Л

_ 2аг _ _ 2аг _ / V

2п

гг

2аг

После преобразований для малых промежутков времени решение, описывающее тепловой режим в неограниченном теле с отверстием радиуса г0 с нормально-тороидальным ис-

точником, при

г' = г, г' = 0

получим в виде соотношения

в

О

8 Ап2

1ггп

г +

Я

2 Л

ехр

4ак

(г - го )2 +

4а

г + -

Я

2

4ак

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

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

Предположим, что теплоисточник /о (рис. 1), определяющий местное поле в области около диаметрального сечения г=0, возникает в момент времени 1=0. Следовательно, предыдущий виток источника при угловой скорости С и величиной подачи $ за один виток источника (на схеме - мгновенный кольцевой источник /1) завершился раньше в момент времени г1 = 2п/ С на расстоянии г1 = $ от источника /0.

Еще раньше в момент времени г2 = 2г1 и

на расстоянии г2 = 2г1 действовал кольцевой источник /2, и так далее вплоть до первого мгновенного кольцевого импульса /т, возникшего в момент времени гт = 2пт /с в положении гт = т$. В соответствии с принятой схе-

Рис. 1. Схема к расчету аккумуляции теплоты в теле с отверстием

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

т

во = Ев (/).

Поскольку повышение температуры местного поля от каждого источника определяется за достаточно короткий промежуток времени, то для расчета температуры общего поля можно

использовать соотношение (3), положив г = г0,г{ = 2пт/с, = т$ , получим выражение

С2

в0 =------Е—ехр

16пАг ы т

8па

т

2

Влияние стока теплоты, который учитывается при изображении процесса в неограниченном теле, можно учесть с использованием соотношения (1). В этом случае следует в последнем положить г' = 0, г = г0, г = г,, г = г,,10 (0) = 1. В результате получим, что снижение температуры местного поля за т витков источника теплоты (при г, = 2лт/С, ) определяется соотношением

в =-

- Qc

(

8ср(ла)/

-Z

о

2пп

exp

-о-

2 + (Sm. )2 8пат

.

Таким образом, в соответствии с приемами изображения для тела, ограниченного изнутри отверстием, местное поле в области действия нормально-

тороидального источника, рассчитывается по формуле (3); повышение температуры местного поля за счет аккумуляции теплоты определяется соотношением для общего температурного поля (4) и снижение температуры местного поля, вызванное стоком теплоты для каждого из т, кольцевых источников рассчитывается по формуле (5).

300

253

203

1)

>-103

53

Без охлаждения СОЖ

С охлажден іем СОЖ

0,01

0,02 0,03

BpeiYffl, СЕК

С,04

0 05

Рис. 2. Зависимость температуры на поверхности отверстия от времени при сверлении стали 45

(А = 60-^; а = 0,08  10-4 —;

м-°С с

ср = 5,02 Л06-ДЖ; я = 2,5  10-3 м; м -°С

Ро = 2174Н; к = 2п; г = г0; г = 0)

V2

r

Рис. 3. Влияние аккумуляции теплоты на температуру на поверхности отверстия

200 180 160 и 140 | 120 S. 100

ф

| 80 £ 60 40 20 О

0123456789 10

Количество учитываемых источи и ков тепла

Рис. 4. Влияние стока теплоты на температуру в центре нормально-тороидального источника в зависимости от радиуса кольцевого источника

-о.о:

о.ос )6

На рис. 3 и рис. 4 показаны уровень аккумуляции теплоты в детали и действие стока теплоты для следующих параметров процесса: г0 = 0,01 м; V = 25 м/с; $ = 0,25 об/мм при обработке стали 45.

Как показывают расчеты, на некоторое повышение температуры местного поля оказывают влияние не более 6 витков источника (рис. 3). Как видно из рис. 4 сосредоточенный сток теплоты не оказывает влияния при обработке отверстия диаметром 10 мм, влияние происходит когда г0< 5 мм.

В рамках выполнения научной работы проводилось сверление глубоких отверстий диаметрами 10 и 20 мм специальными свёрлами для глубокого сверления. Обработка велась на универсальном токарном станке модели 16К20, модернизованном для глубокого сверления. Во время обработки сверлением проводилось экспериментальное исследование температуры в зоне контакта инструмент-деталь, а также распределения тепла в обрабатываемой детали, так как, по данным литературных источников, при сверлении около 75% теплоты переходит в обрабатываемую деталь.

Были проанализированы особенности различных методов экспериментального исследования тепловых явлений (метод искусственных и естественных термопар, пирометрический метод, метод термочувствительных красок, метод оплавления чистых металлов и другие) и предложено измерять температуру методом искусственной термопары.

Так как теоретическое исследование показало, что максимальная температура не должна превышать 300 оС (без охлаждения СОЖ) и 180 оС (при охлаждении СОЖ), то для измерения температуры использовались стандартные термопары типа хромель-алюмель. Термопары были подведены через отверстия диаметрами 2 мм в зону обработки в соответствии с разработанной схемой (рис. 5), которая позволила отследить изменение температуры в материале вдоль оси отверстия, а также в плоскости перпендикулярной оси отверстия. Определить значение температуры непосредственно в зоне контакта инструмент-деталь стало возможно теоретически на основании экспериментально полученного градиента температур (рис. 6 и рис. 7).

ОтЬ под установку

ПтЬ ппг! пгтпнпЬм]

Рис. 5. Схема расположения каналов под установку термопар в обрабатываемой детали

т, с

400т

т, с

400т

100

200

300 С

Рис. 6. Графики изменения температур для точек, расположенных на разном расстоянии от торца заготовки

Рис. 7. Графики изменения температур для точек, расположенных на разной глубине от поверхности отверстия

Выводы

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

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

Экспериментально доказана адекватность предложенной математической модели. Работа выполнена при финансовой поддержке гранта РФФИ 10-08-00669-а.

ЛИТЕРАТУРА

1. Подураев В.Н. Резание труднообрабатываемых материалов / В.Н. Подураев. М.: Высшая Школа, 1974. 590 с.

2. Резников А.Н. Теплофизика резания / А.Н. Резников. М.: Машиностроение, 1969. 288 с.

3. Карслоу Г. Теплопроводность твердых тел / Г. Карслоу, Д. Егер. М.: Наука, 1964. 488 с.

4. Рыкалин Н.Н. Теория нагрева металла местными источниками теплоты / Н.Н. Рыкалин // Тепловые явления при обработке металлов резанием: сб. науч. тр. М.: НТО Машпром, 1959. С.14-45.

5. Барац Я.И., Милованова Л.Р. Метод расщепления уравнения теплопроводности при решении теплофизических задач, возникающих при механической обработке / Я.И. Барац, Л.Р. Милованова / Милованова Л.Р., Барац Я.И. // Вестник Саратовского государственного технического университета. 2010. №1 (44). С. 28-34.

Милованова Людмила Руслановна -

кандидат технических наук, доцент кафедры «Материаловедение», декан механикомашиностроительного факультета Энгельсского технологического института филиала Саратовского государственного технического университета им. Гагарина Ю.А.

Барац Яков Ильич -

доктор технических наук, профессор кафедры «Технологии и оборудование электрофизических, электрохимических методов обработки» Энгельсского технологического института филиала Саратовского государственного технического университета им. Г агарина Ю.А.

Стекольников Максим Владимирович -

кандидат технических наук, доцент кафедры «Технологии и оборудование электрофизических, электрохимических методов обработки» Энгельсского технологического института филиала Саратовского государственного технического университета им. Г агарина Ю.А.

Разуменко Виктория Сергеевна -

студентка 5 курса специальности «Технология машиностроения» Энгельсского технологического института филиала Саратовского государственного технического университета им. Гагарина Ю.А.

Статья поступила в редакцию 27.07.11, принята к опубликованию 9.11.11