№3
2006
621.923
ТЕПЛОФИЗИЧЕСКИЙ АНАЛИЗ КАК ОСНОВА ПРОЕКТИРОВАНИЯ
КОМПОЗИЦИОННЫХ ШЛИФОВАЛЬНЫХ КРУГОВ
Докторант, канд. техн. наук В. И. КОВАЛЬНОГОВ, докторант, канд. техн. наукС М. МИХАИЛИН
Предложены математические модели, расчетные методики и программное обеспечение для численного моделирования теплового состояния заготовки и композиционного шлифовального круга в процессе обработки а также для расчета расхода смазочных элементов такого круга за счет их выплавления. Приведены результаты численного моделирования. Установлены закономерности расхода смазочных элементов композиционного круга, свидетельствующие о возлюэюности рационального выбора значений геометрических параметров такого круга, обеспечивающих согласование расхода режущих и смазочных его элементов в конкретных технологических условиях.
The mathematical models, settlement procedures and software for a numerical modeling of a thermal state of preform and composition wheel at grinding are offered and also for calculation of the rate of flux of lubricating elements of such wheel at the expense of them fusion. The results of a numerical modeling are given. The legitimacies of the rate of flux of lubricating elements of a composition wheel showing an opportunity of rational choose of parameters of such wheel ensuring the coordination of the rate af flux of cutting and lubricating elements in particular technological requirements are set.
Композиционные шлифовальные круги (КШК) [1], сочетая преимущества прерывистых и импрегнированных шлифовальных кругов, позволяют снижать тепловую напряженность в зоне обработки и, как следствие, получать детали с качественными поверхностными слоями при более высокой производительности. Особенно эффективны такие круги в условиях затрудненного транспортирования смазочно-охлаждающей жидкости (СОЖ) в зону обработки (плоское торцовое, торцекруглое, фасонное, внутреннее шли-
I
фование).
Конструктивно КШК — это круг, на рабочей поверхности которого чередуются режущие выступы и заполненные твердым смазочным материалом (ТСМ) вырезы, выполненные в виде пазов, прорезей или каналов (рис. 1). Среди комплекса вопросов, требующих решения при проектировании КШК, в затем и в процессе применения, особое место занимает синхронизация линейного износа режущих (абразивных) и смазочных элементов при одновременном надежном обеспечении заданного температурного режима в зоне обработки в конкретных технологических условиях.
Если расход ТСМ вследствие интенсивного его плавления превышает расход абразивных элементов, то происходит переход к схеме шлифования обычным прерывистым шлифовальным кругом, увеличивается ударная нагрузка на обрабатываемую поверхность, что сопровождается ухудшением ее макро- и микрогеометрии. Напротив, если скорость износа абразивных элементов превышает скорость износа смазочных элементов, то избыток твердой смазки переносится на поверхности заготовки и круга, а также попадает в СОЖ и на поверхности технологического оборудования. В результате ТСМ (как правило, химически активный и дорогостоящий) используется неэффективно, а санитарно-гигиенические условия на рабочем месте ухудшаются.
Расходы режущих и смазочных элементов КШК можно синхронизировать за счет изменения их числа, соотношения между их длинами и подбора ТСМ.
№3
2006
Рис. I. Схема плоского торцового шлифования композиционным кругом: 1
твердый смазочный материал; 3 — заготовка
шлифовальный круг; 2
Массовый расход в плавящегося ТСМ без учета парообразования его жидкой ф
азы
определяется зависимостью
а-
гсм
СМ + СТСМ I- пл
(1)
где <2ТСМ —тепловой поток, поглощенный смазочными элементами КШК; ггсм, сТСМ , Т — соответственно удельная теплота плавления, удельная теплоемкость и температура плавления ТСМ; Г0— начальная температура смазочных элементов (температура окру-
жающей среды).
Мощность теплового потока <2ТСМ, в свою очередь, может быть определена по изменению энтальпии в объеме V смазочных элементов круга
2тсм ~~ /Ртсмстсм
(г,-г0)
с/У,
где ртсм — плотность ТСМ; Т3 — местная мгновенная температура в анализируемой точке смазочного элемента, полученная в предположении, что весь тепловой поток затрачивается на его нагрев; т — продолжительность теплового воздействия на смазочный элемент.
Из анализа выражений (1) и (2) следует, что для расчета массового расхода плавящегося ТСМ необходимо, помимо его теплофизических характеристик, знать распределение в нем нестационарного температурного поля, соответствующего конкретным технологическим условиям обработки. Для получения этой информации разработана математическая модель, методика расчета и программное обеспечение, основанные на совместном решении дифференциальных уравнений теплопроводности композиционного шлифовального круга и заготовки с общим граничным условием в зоне контакта. Это позволяет исключить использование в расчетах информации о распределении тепловых потоков между контактирующими объектами, что вносит существенную погрешность в результаты анализа [2].
Генерируемая в зоне плоского торцового шлифования теплота распределяется между заготовкой, шлифовальным кругом (в том числе его смазочными элементами), стружкой и СОЖ (рис, 1). Кроме того, имеет место теплообмен между нагретыми поверхностями взаимодействующих объектов и окружающей средой. В результате этих процессов изменяется температурное поле всех взаимодействующих при шлифовании объектов.
Температурное поле заготовки в любой момент времени т определяется дифференциальным уравнением теплопроводности, которое в декартовой системе координат
^зЬ связанной с заготовкой (неподвижной относительно нее), имеет вид
9 А
Эт Эх
/
X.
эт
Эх
\
+
V
Эу,
X.
\
дт Эк
\
+
_э
дг.
г
X.
дт
\
дг
(3)
/
\
1 У
№3
2006
где — соответственно плотность, удельная теплоемкость и теплопроводность ма-
териала заготовки; Т — температура; х,у,г — пространственные декартовы координаты.
Для КШК дифференциальное уравнение теплопроводности удобно записать в цилиндрической системе координат {гк, <рк, хк], неподвижной относительно него (вращающейся вместе с ним),
где р , с %Х
Г КШК' кшк кшк
ность КШК;
о с
г кшк кшк
Э Т д
{
дт д/
к
X
кшк
V
дт
д г
\
"Xl-I.ii- 9/ кшк__
г дг
к к
+
1 Э
'"к2 Эф»
I
КШК
\
дт_
Эф,
Л
+
д
д>
К
к.
Э7
(4)
КШК
V
дх
к
/
соответственно плотность, удельная теплоемкость и теплопровод-
Р*
кшк
рк? в точках абразивного элемента Ртсм» в точках смазочного элемента
(5)
С
кшк
<
ск, в точках абразивного элемента стсм, в точках смазочного элемента
(6)
^кшк
Хк, в точках абразивного элемента Хтш, в точках смазочного элемента
(7)
где рк, сК, — соответственно плотность, удельная теплоемкость и теплопроводность материала абразивных элементов круга; гК, фк, хк — соответственно радиальная, окружная и осевая координаты цилиндрической системы координат.
Кроме системы дифференциальных уравнений теплопроводности, математическая модель должна содержать геометрические, физические, начальные и граничные условия однозначности. Геометрические условия задают форму и размеры контактирующих объектов и зоны контакта, а также траектории их перемещения в процессе обработки. Физические условия однозначности определяют теплофизические свойства материалов взаимодействующих объектов и их зависимость от температуры. Начальные условия задают распределение температур во взаимодействующих объектах в начальный момент времени (т = 0). Полагаем это распределение однородным и равным температуре окружающей среды Т^
т -Т 1 ./ •
(8)
Граничные условия однозначности определяют условия теплообмена на свободных поверхностях взаимодействующих объектов и в зоне контакта. В соответствии с избранным подходом граничное условие в зоне контакта имеет вид
X
' дТ 4
кшк
дп
/
+ х
\
к А =0
V
дт
дп
\
я
выд
<7с ~ Ч
ст
ч
пл
(9)
Т = Т
КШК * '
где и , и
нормаль к поверхности соответственно шлифовального круга и заготовки на
поверхностная плотность теплового потока, поверхностная плотность теплового потока, отводимого из зоны контакта соответственно вместе с СОЖ, стружкой и плавящимся сма-
анализируемом участке зоны контакта; ¿/8Ь выделившегося в зоне контакта; а,, д ,
• С ** I
выд 9пл
№3
2006
зочным элементом (при расчете температуры Тэ плотность теплового потока дпл полагаем равной нулю); 71 , Т3 — местная контактная температура круга и заготовки соответственно.
Для плоского торцового шлифования композиционным кругом (далее результаты численного моделирования), с учетом данных работы [3] и поправки на прерывистость круга, принято, что поверхностная плотность теплового потока дШ1{ в зоне шлифования
определяется зависимостью
,г/выд= 383,Зу(°-7(^Г4Л^-°-6, (Ю)
где у — скорость стола; £ = ¿¿р/(яО) — отношение совокупной длины режущих элементов круга к длине его рабочей поверхности; г — число смазочных элементов круга; 1р — длина режущего элемента; О — диаметр круга; 5 — фактическая площадь контакта круга с заготовкой; Г — снимаемый за один ход стола припуск; Е — длина обрабатываемой поверхности (ход стола).
Поверхностная плотность дс теплового потока, отводимого в СОЖ,
N11 , Хг / ч
Яс—^Рп-Т,), (11)
А
оЛ,.
где N11, =—- — число Нуссельта (для расчета можно использовать эмпирические урав-
с1
нения подобия [4]); а — коэффициент теплоотдачи; Хс— теплопроводность СОЖ; с1 — эффективный зазор между поверхностями заготовки и круга (примерно равный половине номера зернистости круга).
Поверхностная плотность д теплового потока, отводимого в стружку, будет
с пг Т
г, — 3 ст з
5
где т — масса материала, удаляемого с заготовки за единицу времени.
(12)
Поверхностная плотность д теплового потока, затрачиваемого на плавление ТСМ,
Ч
0. при < Тпя
пл
*ТСМ ^
5
, при Т > Т 03)
На свободных (внеконтактных) поверхностях заготовки и круга задаем граничное условие 3-го рода теплообмена с окружающей средой (воздухом или СОЖ), Для расчета
соответствующих коэффициентов теплоотдачи целесообразно воспользоваться эмпирическими уравнениями подобия [4].
Систему дифференциальных уравнений (3) и (4) с учетом условий однозначности решали методом баланса на основе явной шеститочечной разностной схемы. Для этого анализируемые объекты разбивали на расчетные элементы в форме сегментов цилиндра для шлифовального круга (рис. 2, а) и в форме прямоугольных параллелепипедов для заготовки (рис. 2, б). Далее для каждого расчетного элемента составляли тепловой баланс, представляющий собой дискретный (алгебраический) аналог дифференциального уравнения или граничного условия.
В результате получили следующую систему алгебраических уравнений.
1. Для расчетных элементов заготовки вне зоны контакта
№3
2006
а)
б)
Рис. 2. Пример разбиения на расчетные элементы, системы координат и направления нумерации элементов
КШК (я) и заготовки (б)
ГП - Т * 1 _ л
1 1 >\].к 1 Лл
/=1
л
\
+ А. х
)
+сЛц.к + ътъ_и+Е,хик+
М+1
(14)
где — температура (/,7, -го расчетного элемента заготовки соответственно в
рассматриваемый и предыдущий момент времени; А^В^С^й^Е^Е^ — коэффициенты;
№3
а
. I А
'"''г
ЬУгН
• Е =
. I А
Ау.
' /Г =
I >} Л-—
7
и+1
Д
Дг.
Дт — время между рассматриваемым и предыдущим моментами времени (шаг по времени); — объем (/, у, £)-го расчетного элемента заготовки; 50, 5,, 5Э, 53, 54, 55 — площадь общей грани (?, у, к)-то расчетного элемента соответственно с (/ - 1,у, к)-мл (/' + 1, у, £)-м, (/,у - 1, к)-м9 (/, у + 1, £)-м, (/, у\ к - 1)-м и (/,у, Ж) -м элементом;
Дг3, Ду3, Дг3 — расстояние между соседними расчетными элементами вдоль соответствующих координатных осей; индексы у теплофизических свойств (ри Хг означают, что эти свойства следует брать по температуре соответствующего расчетного элемента (или элементов), например, \
• 1
означает, что теплопроводность заготовки
взята при средней температуре элементов (/ - 1, у, к) и (/, у, к)\ I
код свободной
поверхности заготовки; а3/ —коэффициент теплоотдачи к /-й свободной поверхности
— площадь общей поверхности (/,у, А)-го расчетного элемента с /-й сво-
заготовки; — бодной поверхностью заготовки; Т(] ную поверхность заготовки.
температура среды, окружающей /-ю свобод-
2. Для расчетных элементов круга вне зоны контакта:
/
тК = тК
1-4
V
10
//
/=5
\
V
+ А х
(15)
К)
1=5
где —температура (/, у, &)-го расчетного элемента круга соответственно в рас-
сматриваемый и предыдущий момент времени; А ,Д.,С.С
коэффициенты:
А. | 50 А. 1
Дт кшк/—-¿А шк/+-,/.Л
Л ---; Д =-^-; С
к (о с ) УК ' * Дг ' К Дг
^РкшкЧшк "'к/Ч
^ I ^ , X X
ЛФЛ ' К= АФЛ, ^ Дх/ Дд,'
Км — объем (/, у, к)-го расчетного элемента круга; Дхк, Дгк, Дфк — расстояние между соседними расчетными элементами круга вдоль соответствующих координатных осей; ак/ — коэффициент теплоотдачи к I -й свободной поверхности круга. 3. Для расчетных элементов заготовки в зоне контакта:
№3
2006
rpKi _ гр КЗ J 1 _ Л
1 i,J.k 1 i,j.k ^1 л
4 Л/
/=1 I
10 Л
/=5
/ -J
(16)
/=1
M
х W 17? Т*'" 4- f Тк/" 4- Г) Ткш 4- Ткш -Г ^JS«, L MtyiA ^ /+1,7,А т ^к/И ^ /,/+U "t"
/»=1
10
+ ]+Я/к [?Выд -ch-q,T~q
1=5
пл
где Глм, —температура (/, у, А:) -го расчетного элемента заготовки, находящегося в контакте, соответственно в рассматриваемый и предыдущий момент времени; А — коэффициент
2Дт
д = --—---
(РзСз )/j(jt ^ ХД» (РкшАшк )„,/,/,£
Ш=1
— площадь общей поверхности (г, у, к)-го расчетного элемента заготовки с рабочей поверхностью круга; = — весовые коэффициенты; Бт — площадь контакта /, &)-го
5/к
расчетного элемента заготовки с т-м элементом круга; М — число расчетных элементов круга, находящихся в контакте с (/,7, к)-и расчетным элементом заготовки в анализируемый момент времени.
4. Для расчетных элементов круга в зоне контакта:
N 9
т^кк — X1 —ÎLT*3
1 i,№ ~ ZJ n 1 » '
/»=1 J ß
(17)
где T™tk —температура (¿,у, к)-го расчетного элемента круга, находящегося в контакте в рассматриваемый момент времени; Sn — площадь контакта (г,у', к)-го расчетного элемента круга с п-м элементом заготовки; N — число расчетных элементов заготовки, находящихся в контакте с (/,у, к)-м расчетным элементом круга в анализируемый момент времени; Sfi —■ площадь общей поверхности (/, у, к) -го расчетного элемента круга с обрабатываемой поверхностью заготовки.
Зная температуру этого и соседних с ним элементов в предыдущий момент времени, можно по (14)—(17) рассчитать температуру любого расчетного элемента в любой момент времени, Последовательность расчета температурных полей в заготовке и шлифовальном круге, которую иллюстрирует блок-схема (рис. 3), реализована в пакете прикладных программ для персонального компьютера. Программный пакет разработан в среде Microsoft Visual С++ и состоит из двух модулей: интегрированного в систему трехмерного моделирования SolidWorks модуля ripe- и постпроцессора, а также рас-
№3
1. Задание геометрических, физических, начальных
и граничных условий однозначности
2. Генерация расчетных элементов заготовки и инструмента
3. Задание начальной температуры всех расчетных элементов т
4. Задание условия окончания расчета: продолжительности теплового процесса ^ или параметра установления процесса £
5. Переход к новому моменту времени т
6. Расчет по уравнениям (14) температур внеконтактных элементов заготовки
нет
11. Визуализация теплового состояния
12. Переприсвоение 7* = Т для всех
расчетных элементов
Рис. 3. Укрупненная блок-схема расчета температурных полей в заготовке и шлифовальном круге при
шлифовании композиционными кругами
четного модуля (решателя). Первый модуль обеспечивает возможность создания геометрических моделей инструмента и заготовки, генерации расчетных элементов и визуализации результатов расчета. Второй модуль реализует изложенную выше методику численного решения системы дифференциальных уравнений, а также позволяет изменять исходные данные (физические, начальные и граничные условия однозначности).
N3 3 2006
Таблица I
и
Максимальная контактная температура Тк в зоне плоского торцового шлифования композиционным кругом
Число г Отношение
смазочных
элементов 1,0 0,78 0,55 0,33 0,11
круга
0 1019,48
2 | 873,12 718,24 578,81 381,13
4 I 873,04 750,35 540,98 382,10
б | 873,97 750,74 542,70 552,33
8 874,69 751,67 616,11 709,74
12 852,32 791,96 621,17 842,12
16 736,16 1014,44 810,41 645,97 1222,0
24 III | 1018,97 912,30 1254,81
32 1018,51 1019,41 997,18 1256,45
Использование пакета 8оМ\¥огкз обеспечило возможность, с одной стороны, вычислять геометрические параметры (объемы, площади поверхностей) трехмерных расчетных элементов и контактирующих объектов, а с другой стороны, способствовало увеличению быстродействия решателя, поскольку все вычисления, связанные с определением объемов и площадей поверхностей расчетных элементов (зависимости (14)—(15)), вынесены из него в модуль препроцессора. Кроме того, пакет ЗоИсГ^огкв является мощным средством визуализации трехмерных объектов и расчетных элементов, а значит, и температурных полей.
Адекватность методики расчета теплового взаимодействия круга и заготовки на основе совместного решения уравнений их теплопроводности с общим граничным условием (9) в зоне контакта, а также соответствующего программного обеспечения, проверена и подтверждена ранее применительно к операциям шлифования (в том числе круглого наружного, торцекруглого, фасонного и многокругового) обычными шлифовальными кругами [2]. При этом отмечено хорошее согласование результатов расчета с экспериментальными данными, расхождение между которыми не превышало погрешности эксперимента. Полномасштабные исследования адекватности разработанной математической модели, методики и программного обеспечения реальным условиям плоского торцового шлифования заготовок композиционными кругами проводятся в настоящее время. Тем не менее, уже предварительные результаты численного моделирования теплового процесса, сопровождающего шлифование композиционными кругами, полученные с их использованием, позволили выявить некоторые общие его закономерности. В качестве иллюстрации в табл. 1 приведены некоторые результаты численного моделирования тепловой напряженности плоского торцового шлифования модельной заготовки в виде прямоугольного параллелепипеда размерами 45x15x15 мм из стали 45 (ИКС 40...42) композиционным кругом 1 — 250x76x25 характеристики 24А16НС18К26 со смазочными элементами из парафина с подачей в зону шлифования 3 %-ной эмульсии Укринол-1.
При численном моделировании варьировали количеством г, смазочных элементов круга и отношением Режим шлифования: рабочая скорость круга ук = 30 м/с, скорость подачи стола \> = 10 м/мин, врезная подача / = 0,03 мм/ход. Приведенные в табл. 1 резуль-
№3
таты соответствуют десяти оборотам круга, когда происходит «квази-установление» теплового режима обработки.
Как видно из табл. 1, зависимость контактной температуры немонотонна как относительно числа г, так и относительно параметра При любом фиксированном значении г с уменьшением параметра £ увеличивается длительность фазы охлаждения (контакта заготовки со смазочными элементами круга) и сокращается длительность фазы нагрева (непосредственно шлифования). При этом увеличивается плотность теплового потока в зоне контакта абразивного элемента с обрабатываемой поверхностью за счет уменьшения его площади, (9) и (10). В результате с уменьшением параметра £ контактная температура сначала уменьшается, а затем увеличивается. Как показали расчеты (табл. 1), в анализируемых технологических условиях обработки при значениях г, больших или равных 6, минимальные значения контактной температуры соответствуют равному 0,33, при меньших значениях г минимум контактной температуры смещается в область меньших £(£ = 0,11).
С увеличением г при фиксированном значении отношения £ контактная температура сначала уменьшается, а затем увеличивается, причем значение г, при котором обеспечивается минимум контактной температуры, для разных £ различно: в анализируемых условиях при £ = 0,78 оно равно шестнадцати, при £ = 0,55 — двум, а при £ = 0,33 — четырем. Повышение контактной температуры при увеличении г сверх этих значений также обусловлено увеличением плотности теплового потока в зоне контакта за счет уменьшения площади абразивного элемента. При этом при значениях £ > 0,11 с увеличением числа смазочных элементов контактная температура выходит на уровень, примерно соответствующий шлифованию обычным кругом. При меньших значениях £ контактная температура с ростом г увеличивается даже сверх этого уровня.
Таким образом, область рационального применения КШК ограниченна: эти круги эффективны только при определенных для конкретных технологических условий обработки сочетаниях их конструктивных параметров. При неправильно выбранных сочетаниях этих параметров переход на шлифование композиционными кругами взамен обычных может оказаться не только нецелесообразным, но и контр-продуктив-ным. Вместе с тем, большинство проанализированных вариантов исполнения композиционного шлифовального круга (табл. 1) позволяет заметно снизить контактную температуру в зоне шлифования по сравнению с обработкой обычным кругом. Так, температура в зоне шлифования обычным кругом составила 1019 К, что превышает температуру эвтектоидного превращения в стали 45 (1003 К). Последнее свидетельствует о высокой вероятности образования шлифовочных дефектов в поверхностном слое заготовки при обработке в анализируемых технологических условиях. Применение в этих условиях КШК с параметрами, ограниченными областью, обведенной в табл. 1, исключает такую вероятность. Окончательный выбор того или иного варианта исполнения КШК должен быть основан на анализе его расхода как в сравнительном аспекте по отношению к обычному кругу, так и в аспекте соотношения между расходом режущих и смазочных элементов. Поэтому на втором этапе исследований для тех же условий подсчитывали расход смазочных элементов с использованием зависимостей (1) и (2) по найденному с помощью изложенной методики распределению температуры в КШК и в зоне обработки (табл. 2).
При анализе данных табл. 2 обращает на себя внимание следующая особенность,! Несмотря на то, что массовый расход смазочных элементов, зависимости (1) и (2), прямо пропорционален локальной температуре, повышение контактной температуры не всегда
№ 3 2006
Таблица 2
Расход смазочной части композиционного круга
Число г Отношение
1,0 0,78 0,55 0,33 0,11
0 0
2 0,292/0,018 0,473/0,014 0,442 / 0,009 0,137 /0,002
4 0,292/0,018 0,454/0,013 0,405 / 0,008 0,139/0,002
6 0,293/0,018 0,455/0,013 0,408/0,008 0,486/0,007
8 0,293/0,018 0,456 /0,013 0,467 / 0,009 0,822/0,012
12 0,281 /0,017 0,498/0,015 0,475 / 0,009 1,104/0,016
16 0,219/0,013 0,518/0,015 0,515/0,010 1,437/0,021
24 0,368/0,022 0,740 / 0,022 0,941/0,019 1,983/0,029
32 0,370 / 0,022 0,741 / 0,022 1,077/0,021 1,986/0,029
Примечание. В числителе массовый расход (мг/с), в знаменателе — размерный (мкм/с).
приводит к увеличению этого расхода: например, при любом фиксированном числе г из диапазона от 2 до 12 максимальный массовый расход смазочных элементов соответствует £ = 0,55, при котором контактная температура не является максимальной (табл. 2). Это свидетельствует о значимом влиянии на массовый расход смазочных элементов параметра определяющего объем V смазочного элемента, зависимость (2). Что касается размерного расхода, то его распределение в целом соответствует распределению контактной температуры.
Широкий диапазон изменения как массового, так и размерного расхода смазочных элементов (табл. 2) при одновременном обеспечении заданного уровня контактных температур (табл. 1) позволяет согласовывать расход режущих и смазочных элементов КШК,
Результаты исследований будут использованы при проектировании композиционных кругов для торцового шлифования, производимых ОАО «Димитровградхим-маш».
Выводы
1. Для обеспечения технологической эффективности КШК на стадии их проектирования следует проводить теплофизический анализ соответствующей операции шлифования с учетом конкретных технологических условий обработки.
2. Предложены математические модели, расчетные методики и программное обеспечение для численного моделирования теплового состояния заготовки и КШК в процессе обработки, а также для расчета расхода смазочных элементов такого круга за счет их выплавления.
3. Выявлены закономерности расхода смазочных элементов КШК, свидетельствующие о возможности рационального назначения параметров КШК, а следовательно, согласованного расхода режущих и смазочных его элементов в конкретных технологических условиях.
№3
2006
СПИСОК ЛИТЕРАТУРЫ
1. Худобин Л, В., В е т к а с о в Н.И. Шлифование композиционными кругами. — Ульяновск: УлГТУ, 2004.
2. Киселев Е, С,, К о в а л ь н о г о в В. Н. Теплофизический анализ концентрированных операций шлифования. — Ульяновск: УлГТУ, 2002. — 140 с.
3. Терган В. С., Л и б е р м ан Б. С., Генис Б. М. Плоское шлифование. — М.: Высшая школа, 1964. — 320 с.
4. Болгарский А. В., М у х а ч е в Г. А., Щ у к и н В. К. Термодинамика и теплопередача. — М.: Высшая школа, 1975. — 495 с.
256 с.