Расчет температур методом конечных разностей при ленточном шлифовании с постоянным усилием прижима при косоугольной схеме ленточного шлифования Текст научной статьи по специальности «Механика и машиностроение»

МАШИНОСТРОЕНИЕ

УДК 621.923

Д. И. Волков, А. А. Коряжкин

РАСЧЕТ ТЕМПЕРАТУР МЕТОДОМ КОНЕЧНЫХ РАЗНОСТЕЙ ПРИ ЛЕНТОЧНОМ ШЛИФОВАНИИ С ПОСТОЯННЫМ УСИЛИЕМ ПРИЖИМА ПРИ КОСОУГОЛЬНОЙ СХЕМЕ ЛЕНТОЧНОГО ШЛИФОВАНИЯ

Исследования направлены на изучение распределения температур при ленточном шлифовании с использованием схемы с разворотом ролика к направлению подачи. Ленточное шлифование; метод конечных разностей; максимальная контактная температура; остаточные напряжения

В настоящее время в качестве окончательной обработки проточной части крупногабаритных лопаток вентилятора распространение получил метод строчечного шлифования на станках с ЧПУ. Одним из основных параметров, влияющих на величину съема и производительность при ленточном шлифовании, является величина подачи на строку Нс [1-5]. Но необходимо отметить, что с увеличением величины подачи на строку Нс увеличивается такой параметр качества поверхности, как высота межстрочеч-ных гребешков Нг, внешний вид такой поверхности приведен на рис. 1.

Рис. 1. Внешний вид обработанной поверхности

Допускаемое значение высоты межстрочеч-ных гребешков Нг невелико, в пределах 0,005-

0,2 мм [1, 2], и в ряде случаев наблюдается выход этого параметра за допустимый предел.

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

Контактная информация: [email protected]

на угол а. Данная схема, по аналогии с известной схемой, применяемой при плоском шлифовании, может быть названа «косоугольной». При традиционной схеме обработки остаточный гребешок формируется всеми режущими зернами, участвующими в снятии припуска (рис. 2, а). В случае применения косоугольной схемы зона контакта между инструментом и заготовкой также разворачивается (рис. 2, б), и остаточный гребешок формируется меньшим числом зерен, что способствует уменьшению его высоты и снижению параметра Яа. При одинаковом усилии прижима съем металла в первом и втором случае одинаков. Поэтому в первом и втором случаях площадь оставляемой канавки будет также одинакова. Во втором случае срезаемый слой имеет большую ширину, но меньшую глубину, чем в первом случае. Это означает, что минимальная величина остаточного гребешка будет достигаться тогда, когда выполняется зависимость:

а = аг^—, (1)

Ь

где /к - длина площадки контакта абразивной ленты с деталью, м; Ь - ширина ролика и абразивной ленты, м.

С целью определения влияния разворота ролика на шероховатость обрабатываемой поверхности и величину остаточного гребешка был поставлен эксперимент - обрабатывались лопатки из материала ВТ3-1 на режимах Уд = = 4,8 м/мин; Нс = 1 мм; Р = 250 Н, лента ТЙЕаС 953БА А16; ¥р = 10 м/с, ролик (рифленый) диаметром 150 мм, 32 выступа, соотношение размеров выступов к впадинам 3:1, твердость эластичного покрытия ролика 70 ед. по Шору, Ь = = 25 мм (ширина ролика).

Результаты экспериментов по влиянию угла разворота оси ролика а на параметр обработанной поверхности Яа (рис. 3) показали снижение

шероховатости поверхности по параметру Яа в 1,2-1,4 раза и снижение величины «остаточного гребешка» на 30-55% (рис. 4), что позволяет рекомендовать данную схему при чистовой обработке.

а б

Рис. 2. Схема формирования обработанной поверхности: а - традиционная схема; б - косоугольная схема

2,0

мкм

1,5

1,0

Яа

0,5

0,0

сс= 0° -0-а= 10° —Д— а = 30°

50

100

Рис. 3. Влияние размера зерна и угла разворота оси ролика на параметр шероховатости поверхности Яа

Рис. 4. Влияние размера зерна и угла разворота оси ролика на величину остаточного гребешка Нг

Для определения зависимости величины остаточных напряжений от скорости шлифования Ур и угла разворота ролика а был поставлен эксперимент - обрабатывались лопатки из материала ВТ3-1 на режимах Уд = 6 м/мин; Нс = 1 мм; Р = 100 Н, лента ТЙЕаС; 953БЛ Л16; Ур = 618 м/с, ролик (рифленый) диаметром 150 мм,

32 выступа, соотношение размеров выступов к впадинам 3:1, твердость эластичного покрытия ролика 70 ед. по Шору, Ь = 25 мм (ширина ролика).

-50

МПа

-100

- 1Ьи

- 2ии -2Ы1 -300

аост-з50

-400

-450

-Л- -1

-о- -2

3

І XI

/

Л

в к V -

12 14 м/с 18 20

Рис. 5. Исследование остаточных напряжений в поверхностном слое (рентгеноструктурный метод):

1 - а = 30°; 2 - а = 15°; 3 - а = 0°

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

1. СОСТОЯНИЕ ВОПРОСА

Изучению тепловых процессов при шлифовании уделялось достаточное внимание [5-7]. В то же время исследований тепловых процессов при ленточном шлифовании выполнено недостаточно. В работе [7] двумерное установившееся распределение температуры определяется уравнением

Т =

2да' кХу

| в-иКо (х2 + и2)0

(2)

где X - коэффициент теплопроводности, Вт/(мх х оС); а - коэффициент температуропроводности, м2/с; д - интенсивность теплового источника, Вт/м2; х - безразмерная глубина; V - скорость движения теплового источника, м/с; и - переменная интегрирования; Н - безразмерная половина ширины источника тепла.

Вместе с тем процесс ленточного шлифования остается одним из самых сложных с точки зрения моделирования тепловых явлений в процессе обработки и прогнозирования максимальных контактных температур в зоне обработки. Результаты работ [5-7] не позволяют учитывать особенности шлифования криволинейных поверхностей эластичным абразивным инструментом.

2. ПОСТАНОВКА ЗАДАЧИ

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

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

В качестве примера для апробации предложенной модели был выбран процесс ленточного шлифования на шестикоординатном лентошлифовальном станке МеаЬо 6КС - 1000. Данный станок позволяет осуществлять обработку криволинейных поверхностей, в том числе и пера лопаток ГТД.

Исходным уравнением для моделирования распределения температур при ленточном шлифовании являлось дифференциальное уравнение теплопроводности (3).

дв 1

дв ,. дв

д 1^7 + ••

дт ер ^ до до

дв .,дв дв .,дв +—(X—)+—(X— ду ду дг дг

(3)

Граничные и начальные условия, необходимые для однозначного решения уравнения теплопроводности (рис. 6), в этом случае представляют систему следующих уравнений:

дв

дУ

дв

дУ

д( х, г) а к (х, г)

У=0

х^(0+умст'т;Н+уис

ге(-1 /2,1 /2)

У=0

хє(0+уист ' т;Н+уист' т

г<-1 /2, г>1 /2)

X

=00 0, X

X

в;

(4)

где д(х, г) - интенсивность теплового источника в точке контакта с координатами х, у, Вт/м2; а0 - коэффициент теплообмена на свободной поверхности, Вт/(м2-К); ае - коэффициент теплообмена в зоне контакта, Вт/(м2- К); Н -ширина ролика, мм; I - длина контактной площадки, мм.

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

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

3. МЕТОДИКА ИССЛЕДОВАНИЯ

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

При расчете на исследуемое сечение тела наносится сетка вертикальных и горизонтальных линий (рис. 3), соответственно і = 0,1,2,3.... Ік; ] = 0,1,2,3....Л; р = 0,1,2,3..РК.

Р=1 Р~2 р=3 £5=4 р=Рк Рис. 7. Схема сеточной области

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

Частные производные при численном методе расчета заменялись приближенными конечно-разностными эквивалентами, а уравнение (3) записывалось в виде выражения (5). Значения I, 7, р в уравнении (5) обозначают номера интервалов пространства вдоль осей Х, У, 2, а индекс к - номер интервала времени т. Уравнение в конечных разностях позволяет определить температуры во всех внутренних узлах сеточной области.

т)

0у,р,к+1 ЛЯ0х_4/',р,к(0,-1,/,р,к + ©г+1 ,/,р,к) +

+ Лр0у_у,р,к(0у-1,р,к+1 + 0у+1,р,к)+

+ ЛЯ

0г_у,р,к(0 у,р+ 1,к + 0,+ 1 ,/,р+1

,к)+

(5)

+ ЛР,

0у_1,},р,к

+

+ (1 (ЛМ0х_17,р,к

+ ЛЯ0г i,j,p,kj)0i,j,p,k,

где ЛГ0х_/7,р,к, Лр0у_7р,к, Лр0г_17,р,к - безразмерные критерии Фурье, определяемые в соответствии с параметрами разбиения рассматриваемой области и текущим значением температуропроводности материала.

Определялись критерии Фурье по системе формул (6)

(в,, 7, р ,к )Лт.

ЛЯ

а ,у

_ ■,7, р,к'

0 х _ ,, 7, р ,к

ЛЯ

0 У _ /, 7, р,к

ЛЯ

Лх2

^, ], р ,к(в,, ], р ,к)Л^. Лу2 !

*,,7,р,к (в,,],р,к ^

(6)

0 г _ ,, 7, р,к

Л

где Г, к (в, 7 рк_х) - коэффициенты температу-

проводности материала определялись по выражению (7), м2/с

(в,, ,■, р,к) =

/,7,р,к ' /,7',р,к

1■■ к(в-■ к)

/,7,р,кУ /,7,р ,к / (7)

,7,р,к (в-,7,р,к )Г

где Хц,к,р(0у,к,р) - коэффициент теплопроводности материала, соответствующий интервалу времени к, определялся для титановых сплавов при помощи выражения (8), м2/с; с(0у,кр) - коэффициент теплоемкости материала, соответствующий интервалу времени к, определялся для титановых сплавов при помощи выражения (7), Дж/(кг- К); р - плотность материала, кг/м3.

1 ■ 7 , (в- - к) = 2,239 10_5 в2 _

1,7,р,кУ 1,7,р,к' ’ , ,7,р,к

_ 2,646 10_3 -в,^ + 9;

(8)

с^хр^Хр) = 0,130у,к,р + 590. (9)

Выражения (7) и (8) были получены при помощи математической обработки данных из работы [8]. Адекватность зависимостей (7) и (8) оценивалась коэффициентом корреляции, равным 0,988 и 0,98 соответственно. Из анализа литературных данных видно, что зависимость Х(0) имеет сложный вид. В одних случаях происходит рост величины Х(0) с увеличением 0, в других уменьшение, причем темп изменения Х(0) разный для разных интервалов температур.

Вычислительный эксперимент, проведенный для следующих исходных данных: Рг = =150 Н, £ = 400 мм2, Н = 20 мм, гладкий ролик диаметром 150 мм, Ур =10 м/с, Уд = 5 м/мин, а0 = 5000 Вт/(м2-К), ак = 2000 Вт/(м -К), показал, что разность максимальных температур, рассчитанных для случая с постоянными теплофизиче-

скими величинами а(320) = 4-10-6 м2/с, Х(320) = = 11 Вт/(м • К), и изменяемыми, рассчитанными с использованием формул (6), (7) и (8), составляет 10%. На практике при обработке лопаток с постоянным усилием прижима на многокоординатных станках с ЧПУ часто возникает проблема с возникновением прижогов.

С использованием метода конечных разностей удалось получить изменение температур в зоне для следующих исходных данных: Рг = = 150 Н, I = 20 мм, Н = 20 мм, Ур = 10 м/с, гладкий ролик диаметром 150 мм, Уд = 1 м/мин, а0 = = 5000 Вт/(м2 • К), ак = 2000 Вт/(м2 • оС), рис. 8, 9 - с разворотом ролика к направлению подачи (а = 45о), рис. 10, 11 - без разворота к направлению подачи.

Рис. 8. Распределение температуры на поверхности лопатки при косоугольной схеме

Тк

500

400 300

100 0

Рис. 9. Распределение температуры по ширине зоны контакта, расстояние между сечениями 5 = 1,6 мм

Рис. 10. Распределение температуры на поверхности лопатки при традиционной схеме

2

500

300

200

Тн

100

г

\

Л

1

40

60

80

100

мм

Рис. 11. Распределение температуры по ширине зоны контакта, расстояние между сечениями 5 = 1,6 мм

4. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ

С целью апробации предложенной модели были произведены эксперименты по замеру максимальных контактных температур по ширине зоны контакта рис. 12.

*

600

1°

500

400

300

Тк

200

100

0/7

[/ ] А

/ 1 /

? /.

Л У

'у<Л > —■— 1

10

15

20

25

Рис. 12. Сравнение теоретических и экспериментальных данных: без разворота ролика: 1 - экспериментальные данные;

2 - теоретические данные; с разворотом ролика (а = 45о); 3 - экспериментальные данные; 4 - теоретические данные

Сравнение экспериментальных данных с данными вычислительного эксперимента показало расхождение не более 15%. Причем расхождение имеет тенденцию уменьшения по краям зоны контакта и увеличения в середине. Это обстоятельство можно объяснить неравномерным распределением теплового потока интенсивности теплового источника вдоль ширины контакта.

ВЫВОДЫ

1. Применение косоугольной схемы ленточного шлифования позволило снизить шероховатости поверхности по параметру Яа в 1,21,4 раза и снизить величину «остаточного гре-

бешка» на 30-55%, повысить долю сжимающих остаточных напряжений

2. Учет теплофизических величин позволяет уточнить величину максимальной контактной температуры на 10%.

3. Применение метода конечных разностей позволило рассчитать изменение максимальной контактной температуры при развороте ролика к направлению подачи.

4. За счет разворота ролика к направлению подачи на угол а возможно снизить температуру на окончательно формообразующих участках зоны контакта на 10-15%.

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

1. Фельдман Е. Ю., Юнусов Ф. С. Бескопир-ное шлифование крупногабаритных фасонных деталей. Казань: изд-во Казанск. гос. техн. ун-та, 1999. 155 с.

2. Юнусов Ф. С., Фельдман Е. Ю. Бескопи-ровальное шлифование крупногабаритных фасонных деталей. Казань, 1999. 155 с.

3. Коряжкин А. А., Михрютин В. В. Повышение точности ленточного шлифования широкохорд-ных лопаток компрессора ГТД на станках с ЧПУ // Справочник. Инженерный журнал. 2008. № 7. С. 7-11.

4. Полетаев В. А. Михрютин В. В., Коряжкин А. А. Ленточное шлифование крупногабаритных лопаток гаотурбинных двигателей на станках с ЧПУ // Справочник. Инженерный журнал. 2005. № 12. С. 7-11.

5. Бабошкин А. Ф. Теория и методы повышения шлифования абразивными лентами: автореферат дисс. на соиск. уч. ст. д-ра наук. С.-Петербург, 2005.

6. Сипайлов В. А. Тепловые процессы при шлифовании и управление качеством поверхности. М.: Машиностроение, 1978. 167 с.

7. Иегер Дж. К. Движущиеся источники тепла и температура трения // Прикладная механика и машиностроение. 1952. № 6. С. 22-39.

8. Резников А. Н. Теплофизика процессов механической обработки. М.: Машиностроение, 1981. 279 с.

ОБ АВТОРАХ

Волков Дмитрий Иванович, проф., зав. каф. резания материалов, станков и инструментов им. С. С. Силина РГАТА имени П. А. Соловьева. Д-р техн. наук. Иссл. в обл. глубинного шлифования.

Коряжкин Андрей Александрович, доц. той же каф. Канд. техн. наук. Иссл. в обл. ленточного шлифования.