Научтруд
Войти
Сайт продаётся: mail@nauchtrud.com

РАЗРАБОТКА ТЕПЛОВОЙ МОДЕЛИ ШПИНДЕЛЬНОЙ ОПОРЫ МЕТАЛЛОРЕЖУЩЕГО СТАНКА

Автор: А.Ф. Денисенко

УДК 621.892

РАЗРАБОТКА ТЕПЛОВОЙ МОДЕЛИ ШПИНДЕЛЬНОЙ ОПОРЫ МЕТАЛЛОРЕЖУЩЕГО СТАНКА

© 2020 А.Ф. Денисенко, Л.Ю. Подкругляк

Самарский государственный технический университет

Статья поступила в редакцию 12.05.2020

На основе использования электротепловой аналогии при описании процесса теплообмена разработана инженерная методика определения теплового состояния опор шпиндельных узлов, не требующая «тяжелых» программных пакетов, основанных на применении метода конечных элементов. Особенности конструкции современных ШУ позволили обоснованно представить тепловую модель в виде плоской модели. Предложены принципы разбиения блоков плоской модели на укрупненные элементы, учитывающие особенности геометрии рассчитываемой конструкции и источников тепловыделения. Разработана методика определения термического сопротивления элемента расчетной тепловой модели с учетом граничных условий. Использование укрупненных элементов тепловой модели позволяет легко оценить влияние технологических особенностей обработки и сборки отдельных деталей конструкции на ее тепловое состояние. Представленная методика разработки тепловой модели может быть использована для оценки стационарных тепловых процессов любых осесимметричных конструкций при незначительных вариациях температур.

ВВЕДЕНИЕ

Работа металлорежущего станка, как правило, характеризуется достаточно интенсивным выделением тепла в электродвигателях, гидросистеме, зонах трения, зоне резания. Технические и эксплуатационные характеристики станка существенно зависят от того, в каком температурном режиме работают его подсистемы. Исходя из этого, существуют обоснованные рекомендации, ограничивающие значения температурного критерия. Так, например, для шпиндельных подшипников рекомендованы примерные величины температуры допустимого нагрева наружных колец подшипников, указанные в табл.1 [1].

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

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

Денисенко Александр Федорович, доктор технических наук, профессор кафедры «Технология машиностроения, станки и инструменты». E-mail: sammortor@yandex.ru Подкругляк Любовь Юрьевна, аспирант. E-mail: podkruglak@mail.ru

Во-первых, конструкции, используемые в станкостроении, имеют точечные или близкие к точечным источники тепловыделения (зубчатые передачи, зона резания, подшипники качения и т.д.)

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

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

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

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

Геометрические условия определяют форму и размеры тела, в котором протекает изучаемый процесс. Физические условия определяют теплофизические характеристики тела. Граничные условия определяют особенности протекания процесса на поверхности тела.

Таблица 1. Допустимый нагрев наружных колец шпиндельных подшипников

Класс точности станка Н П В А С

Допустимая избыточная температура наружного кольца, °С 50 30...35 20.25 15.20 8.10

ИСПОЛЬЗОВАНИЕ ЭЛЕКТРОТЕПЛОВОЙ АНАЛОГИИ ПРИ ОПИСАНИИ ПРОЦЕССА ТЕПЛООБМЕНА

Принимая, что в соответствии с законом Фурье плотность теплового потока в телах, когда температура отдельных его участков неодинакова, прямо пропорциональна градиенту температуры, можно записать значение для теплового потока [2]

Р = (( - Хг )/Е , (1)

где Е = ( — х2 )/(Х£) - тепловой коэффициент;

X - коэффициент теплопроводности;

£ - площадь изотермической поверхности;

х1; х2 - координаты расположения изотермической поверхности;

1Х ;(2 - температура в точках, соответствующих координатам хх; х2.

Сравнивая зависимость электрического тока на некотором участке цепи от напряжения и сопротивления этого участка, описываемую законом Ома,

I = (( — и2 )), (2)

с выражением (1), легко установить аналогии между электрическими и тепловыми параметрами: электрический ток - тепловой поток (мощность теплового потока); электрическое напряжение - температура; электрическое сопротивление - тепловой коэффициент.

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

По аналогии с электрическими цепями величина стт = — х2 ) называется термической (тепловой) проводимостью с размерностью (Вт/ (м2К)), а обратная ей величина #т = ( — х2 )/Х - термическим (тепловым) сопротивлением с размерностью (м2К/Вт).

Наряду с указанными значениями стт и Ят используются значения полной (абсолютной) термической проводимости стТП с размерностью (Вт/К ) и полного (абсолютного) термического сопротивления #ТП с размерностью (К/Вт) [4]:

°ТП = ®ТБ = (( - Х2 ); (3)

Ятп = Ят/Б = ( - х2)/().

Тогда уравнение (1) можно записать в интегральном виде [4]:

Р = сш ( - гг ) = ((1 - гг )/Ктп. (4)

Таким образом, на основании электротепловой аналогии процесс теплообмена может быть

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

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

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

МЕТОДИКА РАЗРАБОТКИ ТЕПЛОВОЙ МОДЕЛИ ОПОРЫ ШУ

Рассмотрим методику построения тепловой модели на примере шпиндельной опоры металлорежущего станка (рис. 1).

Рис. 1. Конструкция передней опоры шпиндельного узла

Для шпиндельного узла (ШУ) разработка тепловой модели, используемой при проектном анализе, может быть упрощена в связи с тем, что конструкция является, как правило, осесимме-тричной, что особенно характерно для гильзованных ШУ. Это позволяет с некоторым допущением привести 3Э-модель к плоской, рассматривая тепловую модель в виде тонкой платины постоянной толщины Н, расположенной по диаметральному сечению шпинделя по одну сторону от его оси. Учитывая одинаковые условия распространения тепла в пластинах можно считать, что теплообмена между пластинами не происходит.

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

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

Я-„ = оХ-„- * + + 0.5ЛГ;

Я -. = 0,5Я[_1 - _ + Я- р + 0,5<;

Я - * = ^Я-п - * + Я- * + 0,5ЯТя ;

-п-* - р

Я- = 0,5^ - р + Я- + 0,5ЯТ.

Тепловое сопротивление элемента Я по аналогии с электрическим сопротивлением определяется зависимостью [5]

Рис. 2. Блоки плоской тепловой модели наружного кольца подшипника передней опоры шпиндельного узла

4-" 14\\ .15 + " + 1 16 + м 4&

+4 1 4 1 . 2 + 1 +

______ +* +7

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

Узлы разрабатываемой тепловой модели удобно разместить в центре прямоугольного элемента (рис. 3). Тогда для любого г -го элемента (рис. 4), имеющего тепловое сопротивление , тепловые сопротивления его сторон Я&С_1 ,ЯС_2 ,Я<г_3 4 (соединения с сопряженными деталями и окружающей средой) и тепло-сопряженных элементов

вые сопротивления Я^; Ятр; Я^; Я^ ячейка тепловой модели будет представлена в виде, показанном на рис. 5,а, где

г_1_р - термические сопротивления тела элемента по соответствующей траектории.

Рис. 4. К определению теплового сопротивления ячейки тепловой модели

Рис. 5. Ячейка тепловой модели

где - элемент координатной линии теплового потока;

X - коэффициент теплопроводности материала;

£(() - площадь сечения, перпендикулярного тепловому потоку в точке I координатной линии;

1 и 2 - изотермические зоны, между которыми расположен участок тела, характеризующийся тепловым сопротивлением.

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

л Ь Ь

(() ХаИ&

где Н - толщина пластины;

Ь - размеры элемента вдоль распространения теплового потока. Аналогично

„т а ¿1 а Ка = Ь

о () ХЪИ& Тогда для центрального элемента на рис. 4

Я^ = Ь/ХаН; Я--р = а/ХЬН. (12)

При составлении тепловой модели наиболее сложной задачей является определение контактного термического сопротивления (КТС) в местах, соответствующих соединениям деталей сборочной единицы.

Тепло с поверхности одной детали к поверхности сопряженной может передаваться следующими путями:

- теплопроводностью через места непосредственного контакта;

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

- конвективным переносом тепла средой, заполняющей это пространство;

- лучистым теплообменом между поверхностями.

Однако, при предварительных проектных оценках можно считать, что условия контактирования в соединениях шпиндельной опоры одинаковы и при оценке контактного термического сопротивления толщина зазора в контакте в среднем вдвое меньше максимального расстояния между впадинами шероховатостей 5тах = 0,03 мм. Поскольку твердые тела касаются друг друга только вершинами профилей шероховатостей, площадь контакта которых составляет 0,01...0,1% от номинальной [6], практически весь тепловой поток передается через воздушный зазор. Тогда считая, что для воздуха X = 0,025 Вт/(м-К), контактное термическое сопротивление соединения деталей подшипникового узла равно 0,0006 м2К/Вт. Тогда полное (абсолютное) термическое сопротивление стыка в соответствии с выражением (3) Яс = 0,0006/ £ К/ Вт, где £ - площадь стыка, м2, и для рис. 4 имеем

Я<-„=Я--8 = о,0006/аН К/Вт и Я?-р = Я£1 =

=0,0006/ ЬН К/Вт. Для более полного учета факторов, влияющих на формирование КТС, можно

воспользоваться исследованиями, приведенными в работах [7, 8].

Если элемент контактирует не с твердым телом, а с газообразной средой (воздухом), то воспользовавшись законом Ньютона для описания конвективного теплообмена между изотермической поверхностью Б с температурой ТБ и омывающей эту поверхность газообразной средой с температурой ТСР , говорящим о том, что количество теплоты, отданное нагретой поверхностью окружающей среде в единицу времени (то есть мощность теплоотдачи РТ ), пропорционально разности температур поверхности и среды и величине теплоотдающей поверхности [5]

Рт ~а{Т8 Тср )£,

где а - коэффициент конвективной теплоотдачи с поверхности в среду, Вт/(м2*К); (для воздуха при естественной конвекция а=5...25 Вт/(м2*К);

ТБ - ТСР - температурный напор, обусловливающий процесс теплообмена;

Б - площадь теплоотдающей поверхности,

можно записать выражение для теплового сопротивления процессу конвективной теплоотдачи с поверхности Б тела в окружающую среду [5], К/Вт:

ЯТа= 1/ И ). (14)

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

Мощность тепловыделения, Вт, в подшипниках качения [9]

б = 1,047 -Ю-4 Мтп, (15)

где п - частота вращения, об/мин;

МТ - момент трения, Н*мм:

Мт = 10-7/ (при > 2000);

V - кинематическая вязкость смазочного материала, мм2/с;

/0 - коэффициент, зависящий от типа подшипника и условий смазывания;

йср - средний диаметр подшипника, мм.

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

РЕЗУЛЬТАТЫ РАСЧЕТА

Воспользовавшись электротепловой аналогией и законом Кирхгофа для определения температур в элементах тепловой модели, составим уравнения равновесия тепловых потоков в ее узлах (табл. 2).

Значения термических проводимостей <_п; <_р ; < ; <_; ячеек тепловой модели, размеры блоков которой указаны на рис. 2, приведены в табл. 3.

Для среднего диаметра подшипника йср= 100 мм; кинематической вязкости масла V = 20 мм2/с и коэффициента, зависящего от типа подшипника и условий смазывания, /0=1 значения удельной мощности тепловыделения приведены в табл. 4.

Результаты расчета температур для ряда узлов тепловой модели при различных частотах вращения подшипника приведены на рис. 6.

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

Таблица 2. Уравнения равновесия тепловых потоков в узлах тепловой модели

№ узла Уравнение равновесия тепловых потоков

1 °1-2 (Т - ¿2 )= Р
2 а2-1(Т2 - ¿х )+о2-4 (т - и )+°2-з (Т - ¿з )+°2-15 (Т - ¿15) = 0
3 °3-2 (Т - ¿2 ) + °3-5 (Т - ¿В ) + °3-17 ( - ¿17 ) + ^3-16 ( - ¿16 ) = 0
4 °4-2( -¿2) + °4-6(Т -¿6 ) + °4-5 (Т -¿5 )+°4-14(Т -¿14) = 0
19 °19-18 ((19 - ¿18 ) + °19-51 (Т19 - ¿В ) + °19-52 ((19 - ¿В ) + °19-53 ((19 - ¿В ) = ^
20 °20-18 ((20 - ¿18 ) + °20-21 ((20 - ¿21) + °20-51 ((20 - ¿В ) + °20-52 ((20 - ¿В ) = ^
21 ^21-17Т21 -¿17) + ^21-51Т21 -¿В )+^21-52Т21 -¿В )+^21-20((21 -¿20) = ^

Таблица 3. Термические проводимости

№ узла Термическая проводимость а • 105 , Вт/К

1 а!-2 = 11,99
2 а 2-3 = 5000 а з-5 = 8,99 а3-17 = 729,2 а 2-15 = 919,5
3 а3-2 = 5000 а 3-5 = 8,99 а 3_17 = 729,2 а3-16 = 689,7
20 а 20-18 = 1300 а 20-21 = 138 а 20-51 = 179,5 а 20-Я 2 = 0
21 а 21-17 = 466>7 а 21-51 = 179,8 а21-52 = 179,8 а21-20 = 13800

Таблица 4.Удельная мощность тепловыделения в наружное кольцо подшипника

П, об/мин 500 1000 1500 2000 3000 4000 6000 8000 10000 15000

Р, Вт/мм 0,0039 0,012 0,024 0,039 0,077 0,124 0,243 0,393 0,57 1,121

1 *С 120

> / х

) 1 а 15 • ■ ■ У&

> И5 (17 \\ А/& Га&

) 1 Л к*

2x10 4x10 6хЮ3

ги об/мин

8хЮ
1хЮ

Рис. 6. Зависимость температуры в узлах № 2, № 5, № 15 и № 17 от частоты вращения шпинделя

ШУ (аз 17=1458,4.10-5 Вт/К; аз 16=1379,4.10-5 Вт/К; а2_15=1839,0.10-5 Вт/К; а4-14=1379,4.10-5 Вт/К; а4-5=1971,8.10-5 Вт/К) обеспечивает снижение температуры £2 (наружное кольцо подшипника) на 2,997°.

ЗАКЛЮЧЕНИЕ

Разработанная инженерная методика определения теплового состояния опор шпиндельных узлов не требует применения «тяжелых» программных пакетов, основанных на использовании метода конечных элементов. Укруп-ненные элементы тепловой модели позволяют

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

СПИСОК ЛИТЕРАТУРЫ

1. Фигатнер А.М. Расчет и конструирование шпиндельных узлов с подшипниками качения металлорежущих станков. Обзор. - М.: НИИМаш, 1971. 193 с.
2. Удалов А.И. Тепловое проектирование радиоэлектронных средств: Учебное пособие / Московский государственный институт радиотехники, электроники и автоматики (технический университет) - М., 2007. 184 с.
3. Лыков А.В. Теория теплопроводности. - М. : Высшая школа, 1967. 600 с.
4. Измайлов В.В., Чаплыгин С.А. Электротепловая аналогия и расчет проводимости дискретного контакта деталей машин // Интернет-журнал «НАУКОВЕДЕНИЕ» Том 8, №2 (2016). URL: http ://naukovedenie.ru/PDF/26TVN216.pdf. DOI: 10.15862/26TVN216 (дата обращения 01.10.19)
5. Овчинников С.В. Введение в теорию теплообмена: теплопроводность в твердых телах: учеб.-метод пособие [Электронное издание]. - Саратов: СГУ имени Н.Г.Чернышевского, 2015. URL: http:// elibrary.sgu.ru/uch_lit/1486.pdf (дата обращения
01.10.19)
6. Демкин Н.Б. Контактирование шероховатых поверхностей.- М.: Наука, 1970. 227с.
7. Денисенко А.Ф., Назаров Н.С. Формирование регрессионной модели контактного термического сопротивления плоских соединений шпиндельных узлов/ Сборка в машиностроении, приборостроении, 2017. №7. С.325-329.
8. Меснянкин С.Ю., Викулов А. Г., Викулов Д.Г. Современный взгляд на проблемы теплового контактирования твердых тел // Успехи физических наук. 2009. Т. 179. № 9. С. 945-970. Б01: ИИрз://^. ог^10.3367/иБ№.0179.200909с.0945
9. Проектирование металлорежущих станков и станочных систем: Справочник-учебник. В 3-х т. Т. I: Проектирование станков / А.С. Проников, О.И. Аверьянов, Ю.С. Аполлонов и др.; Под общ. ред. А.С. Проникова. - М.: Изд-во МТТУ им. Н.Э. Баумана: Машиностроение, 1994. 444 с.

DEVELOPMENT OF THE HEAT MODEL OF THE SPINDLE SUPPORT METAL CUTTING MACHINE

© 2020 А.F. Denisenko, LY. Podkruglyak

Samara State Technical University

Based on the use of the electrothermal analogy in describing the heat transfer process, an engineering technique has been developed for determining the thermal state of the supports of spindle assemblies, which does not require «heavy» software packages based on the application of the finite element method. The design features of modern control rooms made it possible to reasonably present the thermal model in the form of a flat model. The principles of dividing the blocks of a flat model into enlarged elements are proposed, taking into account the geometry features of the calculated structure and heat sources. A technique has been developed for determining the thermal resistance of an element of a calculated thermal model taking into account boundary conditions. The use of enlarged elements of the thermal model makes it easy to assess the influence of technological features of processing and assembly of individual structural parts on its thermal state. The presented methodology for developing a thermal model can be used to estimate stationary thermal processes of any axisymmetric structures with insignificant temperature variations.

DOI: 10.37313/1990-5378-2020-22-3-49-55

Alexander Denisenko, Doctor of Technics, Professor of the Department «Engineering Technology, Machines and Tools». Е-mail: sammortor@yandex.ru Lyubov Podkruglyak, Graduate Student of the Department «Engineering Technology, Machines and Tools». Е-mail: podkruglak@mail.ru

опора шпиндельного узла температура электротепловая аналогия тепловая модель термическое контактное сопротивление spindle unit support temperature electrothermal analogy thermal model thermal contact resistance.
Другие работы в данной теме:
Научтруд |