Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Корчак, С. Н. Производительность процесса шлифования стальных деталей

.pdf
Скачиваний:
12
Добавлен:
23.10.2023
Размер:
20.79 Mб
Скачать

Если положить условие ф (лг, у) = U 0 = const, что вполне реально, то первый интеграл в уравнении (46) существенно упрощается. Непосредственные вычисления показывают что

°° .

г

( * - * ') 2+(у-г/')г

+

J

J -ш г

^

О— со

+ е

( Х+х') 2+ ( у —у') 1

4><*

J d / = [/0.

Таким образом, при вполне реальном начальном усло­ вии U0 = const интегральное решение задачи (39)—(41) имеет вид

 

U (х , у, 0 =

 

 

_2_

 

Я {У', ?)

(У—у')2+х2

 

^ U - n

dy',

Ф

4ях — Р) е

где Dl — заштрихованные

области,

лежащие

ниже пря­

мой

t{ = t (см.

рис.

49).

 

расчетной

Полученное

решение можно

обобщить для

объемной схемы на случай пространственной задачи тепло­ проводности, описываемой уравнением (31). Плоская схема действия источников (во времени и одном направле­ нии пространства Y) становится трехмерной, так как добавляется координата Z (по ширине источника — зерна и всей ширине шлифования). Параллелограммы на рис. 49 превращаются в цилиндры (рис. 50), направляющими ко­ торых являются контуры источников — зерен Фг (z, х) = 0, а образующими — прямые, параллельные прямой:

( У=

К + «к)

1 z =

0.

Эти цилиндры, являющиеся следами во времени кон­ туров плоских YZ источников — зерен, заполняют парал­ лелепипед 0 10 2Н 1Н 2М 1М 2Р 1Р 2 (для длины L зоны кон­ такта) в соответствии с частотой расположения зерен на круге (в плоскости YZ) и временем вступления зерен в ра­ боту друг за другом (по координате Т). На рис. 50 условно показан один источник и его след. В общем случае эти цилиндры-следы заключены между плоскостями

у = 0; у = L; у = vRt\ у = vAt — L0.

148

Интенсивность источника для объемной схемы распре­ делится так:

q (у, Z, t) (внутри цилиндров);

Р (У, 2, t) = О

(вне их).

Рис. 50. Объемная схема прохождения участка детали L через. зону кон­ такта с действием тепловых импуль­ сов (на схеме изоб­ ражен один им­

пульс)

Для полупространства функция Грина имеет вид

G(x, у, z,

t, х 'у у',

z', t') =

1

( х - х ' У + ( У - У ' ) * + ( г - г ’ )*

 

4 и ( t - V )

[4ях (t — /']3/2

 

 

 

( * - |- * ') г-Н г/—j/ 'P + f z - z ') 1

+ е

4к i t - t ' )

 

Проведя совершенно

ту же

цепочку преобразова­

ний, что и в плоском случае, получим интегральное реше-

149

ние для пространственной задачи:

U (х,

у, z,

t) =

J

СС

Я (у',

г', ? )

v

JJ Y n n K ( t - t ' ) Y

 

 

 

о

 

 

 

 

 

 

 

 

 

_

*г+ (у —г/')2+ (г—г ')2

 

 

t/0,

 

 

е

 

« -п

ей/' dz' +

 

где D(- — область

действия

г-го

источника — зерна

в

системе

YZT,

ограниченная

поверхностью

Ф * (У. z)

=

о.

с учетом

ряда

источников

Интегрирование .ведется

тепла, попавших в рассматриваемую зону. Если интере­ сующая нас точка с координатами (х, у, z, t) попадает внутрь цилиндра (см. рис. 50), т. е. в зону действия тепло­ вого источника р 0 {у, z, t), то происходит нагревание и температура повышается; если точка оказывается вне

действия источника и р (у, z, t)

=

0, то происходит про­

цесс остывания.

 

 

 

 

Анализ расчетной

схемы

 

и интегрального решения.

 

Вывод расчетных формул

Полученное

интегральное решение применимо лишь

в зоне шлифования (полоса шириной OL, см. рис. 49),

так как только

в этом случае

правомочно допущение

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

дачу.

Применим

интегральное решение только при

0 <

г/ < L, так

как для проверки основной гипотезы

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

Прежде чем переходить к рассмотрению теплоотвода по направлениям Z и Y, отметим еще раз масштабные соот­ ношения размеров источников тепла (зерен) и глубину поверхностных слоев детали в направлении + Х , прогрев которых необходимо определить. Размер источника не-

150

сколько больше круга диаметром /3 = 0,1 мм (так как надо еще учесть зону пластического сдвига), а средне­ статистическая толщина среза одним зерном составляет, как было показано, порядка 0,005 мм, т. е. в 20 раз (ли­ нейно) меньше.

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

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

Таким образом, в зоне контакта (зона ОМ PH, см. рис. 48) можно полагать преимущественно одномерное распределение тепла в глубь металла (по крайней мере, до глубины, равной толщине среза), исключая поверхност­ ные ( YZ) границы зоны контакта, где отвод тепла будет двумерным. Однако ширина «запретной» зоны по пери­

метру зоны

контакта будет

небольшой, соизмеримой

с толщиной среза

(а — 0,005

мм), т. е. с глубиной рас­

сматриваемого

нагреваемого

слоя.

Таким образом, для решения локальной задачи, ка­

сающейся тонких

поверхностных слоев детали в зоне

151

контакта с кругом, целесообразно упростить принятую расчетную схему. Вырежем в полупространстве х ^ О полубесконечный стержень (рис. 51) с осью, перпендику­ лярной границе полупространства, и весьма малым попе­ речным сечением. Размеры сечения стержня будем счи­ тать существенно меньшими по сравнению с размерами источника — зерна. Исследуем распределение темпера­ туры по выделенному стержню.

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

1---о

О г t

Рис. 52. Схема интенсивности источников и их действия во времени

На рис. 51 показана графическая расчетная схема с уче­ том многократного действия источников. Схема действия источников во времени может быть получена (рис. 52) из схемы для плоской задачи (см. рис. 49) как сечение ее любой прямой (в любой точке длины линии контакта), например АА, параллельной оси Т. Таким образом, след источника в виде параллелограмма выродился в линейные

отрезки оси

Т .

 

Интенсивность источника может быть задана функцией

 

q (t)

t£ (Tp, T t -1-т)

p(t) =

.0

t (Гt + t; T /+i) i = 0, 1, 2 ,..., n. (46)

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

G (х,

t, х', f) =

1

X

 

/ 4 яу. (/—/')

 

г

(х-х'У

(*+*')*

1

(47)

X

4* <<-П

4 х (<-(') J .

152

Учитывая

выражение (47)

по аналогии с решением

на стр. 150,

сразу получим

 

 

t

J£2

U(x,t) = U0 + - ^ \ -

e ~ ~ ^ = n d t ' . (48)

ФJ у 4ях ( t t )

0

Вэтом случае интегрирование практически ведется по всем интервалам (Г,-; Тг + т) с отличной от нуля интен­

сивностью источников,

лежащим

левее точки

t' = t

(см. рис. 52).

 

 

 

Учитывая проведенное выше осреднение интенсив­

ности источника, можем положить q (t) = q =

const.

В этом случае интеграл в решении уравнения (48)

может быть выражен через специальные функции.

Рассмотрим интервал

времени

0 •< t <f т. Для

этого

интервала интегральное решение уравнения (48) примет вид

U (.х,

 

 

 

t

 

 

 

X2

 

t) = Ua + % - [

r

J .

. - Г

4* (*-*'> dt'.

(49)

у

'

ср

J

Y4nv.(t — t')

 

v

Используя

введенные

Хартри

функции

 

 

 

 

 

 

 

СО

 

 

 

 

 

i n+ i e r fc х _ J (п e r fc х d x ,

 

(5 0 )

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

со

 

 

 

где t° erfc л: =

erfc х =

 

%

| е~хг dx

функция

оши-

бок, выражение (49)

 

У 71

X

 

 

 

 

преобразуется:

 

 

 

U(x, t) = U0

- ~ = V ti erfc —р = .

(51)

 

 

 

 

 

уДср

Y 4v.t

 

Для интервала т <

 

t

«< Т г решение

(48) будет

 

 

 

 

 

t

 

l

■e

dt,

 

U(x,

f) = U0 + %-

J

 

 

 

 

 

 

 

Ф

 

}Дзтх (t t')

 

 

 

 

 

 

о

 

 

 

 

 

153

которое с использованием первообразных от функции ошибок (50) приводится к виду

U(x, t) = u , + ^ { y t i ^ k ^ = - ~

- l / f — ieric V i x

(52)

Таким образом, выражение (51) описывает нагрев металла от воздействия первого импульса, а (52) — осты­ вание после ухода первого зерна-источника.

По мере прохождения через выбранное сечение каждого последующего источника в выражении для U (х, t) до­ бавляется только новое слагаемое, структурно совпадаю­ щее с имеющимися. Используя единичную функцию Хеви­ сайда, можно построить аналитическое выражение U (х , t) для п источников, справедливое в любой момент времени:

 

 

 

 

П

 

 

 

 

 

и (х,

t) = и о +

-Ш=г V

 

erfc ..

*

X

V ’

 

Vkcp

 

к

/ 4 x ( t - T k)

 

 

 

 

k—0

 

 

 

 

 

 

 

x H ( t - T k) - V t

— x - T kiX

 

 

 

 

X erfc _Л,

 

H(t — x — Tk),

(53)

 

 

 

/ 4 х ( / t Tk)

 

 

 

 

 

 

 

0 a: <

0 — единичная

функция,

введен-

где Н (х)

[ 1 x ^ 0

ная Хевисайдом.

 

 

Покажем, как выражается решение (53) на примере

интервала

Т г <

t <

7 \ + т.

В

указанном

интервале

разности

t Tk

отрицательны

 

для

всех

k,

начиная

с k =

2,

поэтому

Н ( t Tk) =

0 для k ^ 2 ,

т. е. эти

слагаемые не будут присутствовать в выражении (53).

Аналогично

исчезнут слагаемые, содержащие

разности

(t — т — Тк)

для

k p s l . Таким образом, для

интервала

Т г < t < . Т 1 + т

решение (53) примет вид

 

U(x, t) = U0 + 2q

YUerfc У ш

 

 

 

УХс(>

 

 

Y t — xi erfc

У 4k (t t) +

 

 

+ l/V — T-J, erfc

X _______

(54)

 

 

У4я ( t - T J

154

Функции

in erfc x

изучены

и затабулированы.

Таб­

лица первых

шести

функций

приведена в

работе

[23].

В сочетании

с таблицами выражение (53)

вполне

при­

годно для инженерных расчетов и позволяет строить температурное поле детали в зоне шлифования как ре­ зультат совокупного действия зерен-источников с учетом их движения по дуге контакта.

Для анализа обрабатываемости сталей шлифованием наиболее существенна температура остывания; это тем­ пература, до которой успеет остыть металл в рассматри­ ваемом сечении (нагретый действием предыдущих зерен) к моменту подхода очередного зерна.

Принимая, что зерна на круге расположены равно­ мерно, получим Т { = iT (см. рис. 52); Tt — время входа в рассматриваемое сечение i + 1 зерна. Температура ме­ талла в этот момент и определяет его теплофизическое и механическое состояние. Из выражения (53) можно полу­ чить формулу, определяющую температуру остывания в зависимости от номера (т) последнего прошедшего через сечение источника. Для этого вычислим значение U (х , t)

для t =

т Т :

m

 

 

 

 

 

 

 

 

 

 

U (х , тТ) = U0-

-%=■ V V{k — i) Ti erfc -7__ *

 

0

Y t e p - j i i J

 

'

Y i x T { k - i )

 

 

 

г=о

 

 

 

 

 

У (k i) Т — т erfc ______ X______

(55)

 

 

 

/

4к (kT — Д=7Т)

 

Меняя индекс суммирования,

из

выражения (55)

по

лучим

 

U (х , тТ) U„

 

 

 

 

 

 

 

m

 

 

 

 

 

2дУт

erfc

У k ai erfc

(56)

УХер

H V i i

Yk

 

 

Yk-

 

 

k=i

 

 

 

 

 

где 6

X T

(57), a -f- (58) — безразмерные крите­

рии процесса.

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

155

 

 

 

 

 

Т а б л и ц а 9

 

Критерии а

и 6 для различных количеств зерен

 

 

 

 

 

б ДЛ Я

 

п

Т в мкс

а

х — 0 мкм

х = 2 мкм

х — 5 мкм

* = 1 0 мкм

 

 

 

8

500

0,006

0

0,0155

0,0386

0,077

1 2

333

0,009

0

0,0190

0,0475

0,095

2 0

2 0 0

0,015

0

0,024

0,0612

0 , 1 2

Для различных условий шлифования при оценке обра­ батываемости а ^ 0,015 и б 0,12, т. е. эти величины весьма малы. Воспользовавшись этими оценками, можно существенно упростить формулу (56). Разлагая функцию

Рис. 53. Соотношение интен­ сивностей тепловыделения от пластического деформи­ рования металла (а) и от трения зерна о металл (б)

ттттгтпт

х t

Б)

под оператором суммирования в ряд Тейлора в окре­ стности (а = 0, б = 0) и ограничиваясь в силу быстрой сходимости получающегося ряда двумя первыми членами разложения, получим с учетом выражения (58):

_ т 1 _ _ 62_

у(*'тГ)=и”+7Щ

«

где т — 1, 2, . . п,

k—1

 

Погрешность формулы (59) для сформулированных выше условий не превышает 1,5%. Формула (59) позволяет без построения всего температурного поля сразу находить интересующую нас температуру остывания.

При интегрировании для упрощения было сделано допущение, что q (t) = q — const, т. e. была усреднена

156

Т а б л и ц а 10

Изменение температуры от интенсивности и длительности действия источника

Глубина

 

Щтих (нагрев)

 

от поверх­

 

 

 

 

ности X

Для источника 1 (10д; 0,1т)

Для источника 2 (^т)

в мкм

0

^шах (0.01 г ) -3,16

 

 

2

 

 

Umax<2-

t ) - 0-8

 

 

 

5

" - " №

 

 

 

 

 

10

 

 

^ т а х '10'

* ) - ° . 136

 

 

 

 

 

£/0 (остывание й(=100т)

 

0

и'

(0; 100-с) = 0.05652-4^1-

П20СТ(»1 100Т) - 0,05650

 

0l-T

УХер

 

 

2

,

аУ Т

Ц2ост (2; Ю0Т) -0,05645 ^

Щ СТ (2; 100т) =0,05645

 

ост

УУГср

 

 

5

Щ

(5; 100т) =0,05640 - ? - =

U0CT (5; W0X)-0,06640

 

0СТ

с р

 

 

10

и '

(10; 100т) = 0,05610—™~Lr

Н2СТ(10; 100т)-0,05610

 

ост

УХ ер

 

 

интенсивность источника. Но из рис. 44 видно, что зави­ симость q (t) имеет большие колебания около ее среднего значения. Простота полученных решений (53) и (59) по­ зволяет оценить эффект усреднения интенсивности источ­ ника. Так как зависимость q (t) кусочно-постоянна, цнтеграл для нее вычисляется по формулам (53) и (59), которые содержат только больше слагаемых. Для упро­ щения оценок отбросим в q (t) постоянную составляю­ щую и рассмотрим различие температурных полей от двух источников (рис. 53).

Суммарные количества тепла, выделяемого этими источниками, примерно одинаково;

Qi = 10<7«0, 1т = qx\ Q2 = qx.

157

Соседние файлы в папке книги из ГПНТБ