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

книги / Тепловые процессы в технологических системах

..pdf
Скачиваний:
3
Добавлен:
12.11.2023
Размер:
15.09 Mб
Скачать

деления f (х„, у„, г„, т)

и аппро­

 

 

 

 

ксимируют его той или иной

 

 

 

 

идеализированной

функцией.

 

 

 

 

По мере

 

совершенствования

 

 

 

 

техники

эксперимента

и разви­

 

 

 

 

тия

теории

контактных

про­

 

 

 

 

цессов

 

при

 

взаимодействии

 

 

 

 

твердых

тел в системах, формы

 

 

 

 

идеализированных законов f (хи,

 

 

 

 

г/и,

 

т)

все

с большей точно­

 

 

 

 

стью отражают

реальные физи­

 

 

 

 

ческие

явления.

 

 

 

 

 

 

 

 

Покажем

это

на

примере

 

 

 

 

эволюции

идеализированного

 

 

 

 

закона

распределения

плотно­

 

 

 

 

сти

тепловыделения,

вызван­

 

 

 

 

ного

 

трением

на

площадке

Рис. 1.19. Графическая интерпретация

контакта

 

между

стружкой

и

передней

 

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

резца.

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

те­

 

пловыделения (см. табл. 1 .1)

 

В ранних работах по тепло­

 

 

 

 

 

физике

резания

тепловыделе­

 

 

 

 

ние

на

этой

 

площадке

полагалось

равномерным

(вариант

7,

рис.

 

1.18).

В

дальнейшем

экспериментальным

путем было

установлено,

что силы трения на площадке контакта рас­

пределены

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

по

закону,

близкому к

нормальному.

ГВ связи с этим и тепловыделение стали описывать несимметрич­ ным законом нормального распределения (вариант 2). Впослед­

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

На рис. 1.19 приведены законы распределения плотности тепло­ выделения, часто используемые при теплофизическом анализе технологических подсистем. Для простоты изображения эти за­ коны даны для одномерных источников, поскольку из одномер­ ных всегда можно сконструировать двухили трехмерные законы. Математическое описание функций f {хи) приведено в табл. 1.1.

Приняв тот или иной вид функции / (л:„, */и, zM, т) и зная q0f

можем рассчитать количество теплоты по формуле (1.47). На прак­ тике, однако, при анализе тепловых процессов в технологических подсистемах приходится, как правило, решать обратную задачу. Обычно известны количество теплоты Q или средняя мощность тепловыделения W за время %и требуется при заданном законе распределения определить наибольшую плотность q0. Полагая,

что закон распределения не зависит от времени, и подставляя вы­ ражение (1.48) в уравнение (1.47), получаем

 

(1.49)

где

 

Рю ^и) dxB dya dzЦ.

(1.50)

Формула (1.49) справедлива для трех-, двух- и одномерных источников. Различие будет лишь в интеграле /, который для дву­

мерного источника

примет вид

 

 

 

/j = JJ / (хв,

уи) dxB dyB,

(1.51)

а для одномерного

 

 

 

 

h = \ f ( x B)dxB.

(1.52)

 

/

 

 

Для двумерного равномерно распределенного источника

(рис. 1.19) f (хв) = f (ув) = 1; h = \d yB\d xB = ab

и <7о = -^-*

Если двумерный

о

о

 

источник распределен по нормальному не­

симметричному закону в направлении оси хя и равномерно —

вдоль оси ув, то / (хи) =

exp [—kQxl ] и f (ув) = 1. Тогда

ь

I

l

h = | dyB ( exp [— Ло*и] dxB — b | exp [— koxX\ dxB.

оо

Известно, что

 

f exp [— p V ] du= ^Q - (erf [pu2] — erf [,puJ),

(1.53)

«»

 

где erf [%] — функция, называемая интегралом вероятностей или интегралом ошибок (от французского erreur fonciion — функция

ошибок). Функция erf (х 1» значения которой в зависимости от аргумента %табулированы [4 ], имеет следующие свойства: erf [0] =

= 0; erf [оо ] == 1; erf [—

—erf [%]. С учетом свойств функ­

ции получаем

 

' * =

w r e r,[' 1/1:1

и далее

2WVk^ _

 

b 1 / я erf [/ V *ol

Обратим внимание на то, что коэффициент kQимеет размерность м~2.

С точностью, достаточной для ин­ женерных расчетов, функция erf [%1

может быть аппроксимирована выра­ жением

erf Ш « У 1 — ехр[— 1,26х2]. (1.55)

Для осесимметричных источников за­ кон распределения плотности тепловыде­ ления удобно задавать в полярной, ци­ линдрической или сферической системах координат. Рассмотрим, например, дву­ мерный круговой источник с нормальным законом распределения (рис. 1.20). Для него f (г) = ехр [—k0r2]. Выделим на пятне контакта диаметром 2R элементар­ ный участок площадью г dtp dr. Тогда

2Я R

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

W = <7„ | dqp | гехр [— ktf2]dr.

ао

Положив k0r2 = и, после интегрирования получаем

а-

”° я ( 1 — ехр[— Л*/?а])

Функция ехр [—k0/?я], описывающая закон распределения тепловых потоков, обращается в нуль при R -*• оо. Следова­

тельно, теоретически пятно контакта нормально-кругового источ­ ника имеет бесконечные размеры. Однако уже при k0R2 3 ехр [—k0R2] 0,05, т. е. ордината кривой, описывающей закон

распределения, на окружности радиусом R = Y 3/й0 составляет только 5 % ординаты в центре источника. Полагая k0R2& 3,

получаем

<7о » W/R2.

(1.56)

Аналогично рассчитывают величину q0 для источников с другими

законами распределения.

Скорость перемещения и длительность функционирования

источников. По скорости перемещения источники разделяются на неподвижные (о = 0), движущиеся (о Ф 0) и быстродвижущиеся (о > У). Быстродвижущиеся — это источники, скорость v перемещения которых превышает скорость V распространения

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

где / — длина источника в направлении перемещения, м» v

скорость перемещения источника, м/с; со — коэффициент темпе­ ратуропроводности материала, в котором (или по которому) пе­ ремещается источник, м2/с.

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

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

Длительность функционирования источника

характеризуют

безразмерным критерием

 

Fo = шт/Р,

(1.58)

где и — время, в течение которого действует источник. Безразмерный комплекс Fo носит название критерия Фурье

или безразмерного времени. По длительности функционирования источники можно разделить на мгновенные (Fo -► 0), действующие в течение конечного промежутка времени (Fo > 0) и действующие столь длительное время, что процесс теплообмена под влиянием источника можно полагать устанбвившимся (Fo -*■оо). Предель­

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

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

Начальные и граничные условия. Начальные условия отвечают на вопрос о том, каково было температурное поле в момент вре­ мени, принятый за начало отсчета. Они описываются выраже­ нием 0|т=о = f0 (х, у, г). Очень часто температура компонентов

технологических подсистем в начальный момент времени может быть принята равной температуре 0О окружающей среды, т. е. fo (х>Уу *) = В этом случае удобно, как отмечалось выше,

вести расчет в так называемых избыточных температурах, условно считая, что /о (х, у, z) = 0, а затем по окончании расчета к ре­ зультату прибавляя 0О.

Граничными называются условия взаимодействия поверх­ ностей тел с окружающей средой или другими телами. Различают несколько разновидностей граничных условий. При граничных условиях первого рода (ГУ1) предполагают, что известен закон распределения температур на граничных поверхностях тела 0а = = fa (х, у, г, т). Пусть, например, требуется определить темпе­

ратурное поле внутри какой-либо детали или инструмента. Сде­ лать это экспериментальным путем, не разрушая объект измерения, довольно трудно, измерить же температуру на поверхности де­ тали, инструмента или другого твердого тела экспериментальным путем значительно проще, это может быть выполнено без повреж­ дения объекта. Если мы знаем ГУ1 в виде закона распределения температур на поверхностях тела, то, решая дифференциальное

уравнение теплопроводности,

можем рассчитать

поле температур

внутри

детали,

инструмента

и

т. д. Частным

случаем ГУ1

яв­

ляется

условие

изотермичности

поверхностей

тела, т. е. 0а

=

= 0! =

const.

 

 

 

 

 

Граничные условия второго рода (ГУ2) предусматривают, что известен закон распределения плотности тепловых потоков qa =

ф (х, у, г, т), следующих через граничные поверхности. В част­ ном случае qa = 0. Это означает, что рассматриваемая поверх­

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

Граничные условия третьего рода (ГУЗ) используют в том случае, когда теплообменом поверхности с окружающей средой пренебречь нельзя. В этом случае должны быть заданы темпе­ ратура 0Осреды, с Которой соприкасается данное тело, и так на­ зываемый коэффициент теплоотдачи а, Вт/(м2-°С), характери­

зующий теплообмен между средой и поверхностью (см. гл. 3). Согласно закону Ньютона—Рихмана плотность теплового по­ тока пропорциональна разности температур поверхности 0S и

окружающей ее среды 0О, т. е.

 

Яа = а ({я — 0О).

(1.59)

Формула (1.59) Дает возможность определить количество теп­ лоты <7s, Вт/м2, Которое в единицу времени с единицы поверхности

отводится в окружающую среду. Как следует из закона Фурье [см. (1.9)], к Поверхности тела подводится поток

Яв = — bgrad0s ------

Следовательно,

Выражение (1.60)' представляет собой математическое описание граничных условий третьего рода.

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

0S1 = 0S2. (1-61)

С целью упрощения расчетов часто вместо равенства темпе­ ратур в каждой точке контакта в качестве ГУ4 принимают ра­ венство средних температур на поверхности контакта, т. е. вместо формулы (1.61) полагают

©si = 0S2-

(1.62)

Граничные условия четвертого рода,

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

в гл. 2, используют при решении балансовых задач, т. е. при ана­

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

Заканчивая рассмотрение вопроса о граничных условиях, отметим, что на разных участках реальных тел могут иметь место различные граничные условия. Рассмотрим, например, процёсс плоского шлифования заготовки торцом чашечного круга (см. рис. 1.16). Если решена задача о распределении теплоты шлифо­ вания между кругом и заготовкой, то по отношению к заготовке имеем следующие граничные условия: ГУЗ — на поверхности соприкосновения с жидкостью; ГУ2 — на контактной поверхности с кругом, где известна плотность теплового потока, и на торце заготовки, который можно считать адиабатическим, если пре­ небречь его теплоотдачей в воздух; ГУ4 — на поверхности, где заготовка соприкасается с магнитным столом станка.

1 .в . К О Д И Р О В А Н И Е Т Е П Л О В Ы Х З А Д А Ч

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

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

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

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

мко TV

р. сд 1у

Первая группа символов содержит информацию о мерности (М), конфигурации зоны тепловыделения (К) и ограниченности

(О) источника. Вторая группа информирует о законах распреде­ ления плотности тепловыделения по осям координат (Р), ско­ рости движения источника (С) и длительности его функциониро­ вания (Д). Третья группа, состоящая из двух символов, фикси­ рует форму тела (Т), на котором действует источник, и род гра­

ничных условий (У).

Каждому из символов в конкретных условиях соответствует одна цифра (табл. 1.2). Исключение составляет символ Р, который при постановке в код тепловой задачи должен содержать три цифры, поскольку должно быть указано распределение плот­ ности теплообразующих потоков по всем трем осям координат. Если мерность источника такова, что плотность потоков распре­ деляется по двум или только по одному направлению, то на соот­ ветствующих местах В символе Р должны быть поставлены нули.

Проиллюстрируем кодирование на примере условий задачи, с которой мы начали этот параграф. Примем М =2 (источник дву­ мерный); К = 1 (площадка плоская прямоугольная); 0 = 2 (источник ограничен по двум направлениям), Р =710 (по оси х„ распределе­ ние симметричное нормальное, по оси уя — равномерное, по оси 2И— нет, так как источник двумерный). Далее С =1 (источник движется); Д = 2 (действует длительно); Т = 1 (на полупространстве) и У = 3 (при ГУЗ). Итак, полный код рассматриваемой задачи

имеет вид

212

7 1 0 .1 2

1 .2 . К оди рован и е

теп л о в ы х з а д а ч . С труктура кода.

j y

 

 

Р. ВД

Символ

Признак источника и тела

Коды

м

к

О

р

(см. рис. 1.19, табл. 1 .1)

С

д

т

(рис. 1.15)

У

Точечный; одно-, двух- и трехмерный

0;

1;

2;

3

Прямой, плоский

(полосовой,

прямоугольный),

 

 

1

 

призматический

 

 

 

 

2

 

Кольцевой, круговой, цилиндрический

 

 

 

Шаровой

 

 

 

 

3

 

Неограниченный;

ограниченный

по направлению

0; 1 2; 3

одной, двух или трех осей координат

 

 

 

 

 

Распределенный равномерно

 

 

 

1

 

Распределенный линейно

 

 

2;

3

 

Распределенный по экспоненте

 

 

 

4

в

 

Нормально распределенный несимметричный

 

5;

 

Нормально распределенный симметричный

 

 

7

 

 

Комбинированный

 

 

 

 

8

 

 

Неподвижный; движущийся; быстродвижущийся

0;

1;

2

Мгновенный, действующий некоторое время; дей­

0;

1;

2

ствующий длительно (процесс установился)

 

 

 

 

 

Неограниченное тело; полупространство; пластина;

0;

1;

2;

3

параллелепипед

 

 

 

4

 

 

Стержень неограниченный

 

 

 

 

Стержень, ограниченный с одной стороны

 

5

 

 

Стержень конечной длины

 

 

6

 

 

Цилиндр

 

 

 

7

 

 

Клин

 

 

 

8

 

 

Шар

 

 

 

9

 

 

Граничные условия

1, 2, 3, 4 рода

1;

2;

3;

4

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

скорости движения.

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

+

212

83

212

82.

8 1 0 .0 2

110.02

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

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

Вопросы для самопроверки к пп. 1.5 и 1 . 6

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

2.Перечислите основные условия однозначности, дополняющие дифферен­ циальное уравнение теплопроводности при решении конкретных задач.

3.Чем отличаются пассивные граничные поверхности от активных? Какие манипуляции можно выполнять с пассивными граничными поверхностями при схематизации формы тел?

4.Как определить количество теплоты, выделенной источником за некоторое время, если известны максимальная плотность тепловыделения и закон распреде­ ления потоков на площадке контакта между источником и нагреваемым телом?

5.Перечислите основные свойства функции, называемой интегралом веро­

ятности.

6 . Каковы особенности быстродвижущихся источников теплоты?

7. Назовите различные виды граничных условий и дайте каждому из них краткую характеристику.

8 . Из каких групп символов состоит код, описывающий особенности тепло­ вых задач?

Задачи к пп. 1.5 и 1.6

13.При изготовлении отверстия D =

а= 2 0 мм в стальной заготовке (рис. 1 .2 1 )

к

сверлу, вращающемуся с частотой п =

=

300 мин-*1, приложен момент М =

=6 6 Н -м. Работа сверления полностью

переходит в теплоту. В первом прибли­

 

 

жении

можно

принять, что теплота вы­

 

 

деляется

непосредственно

на

режущих

 

 

кромках

инструмента. Рассчитать плот­

 

 

ность потока а точке О на наружном диа­

 

 

метре сверла при следующих вариантах

 

 

законов тепловыделения: 1) равномерно

 

 

вдоль режущих кромок (длиной попереч­

 

 

ной кромки Д пренебречь);

2 )

пропорцио­

 

 

нально

скорости резания

в каждой точ­

 

 

ке кромки; 3) по экспоненциальному

 

 

закону.

 

 

решения и

комментарии

 

 

Алгоритм

Рис. 1.21. Распределение плотности

к нему:

 

 

 

 

 

 

 

мощность

тепловыделения

на кромке сверла:

а)

рассчитать тепловую

1 — равномерное;

2 — по линейному

процесса,

приходящуюся

на

одну кром-

закону; 3 — по

экспоненциальному

ку (W =

1036

Вт);

 

 

закону

 

Рис. 1.22. Схема процесса бесцентро­ вого шлифования:

J — ведущий круг; 2 — заготовка; 8 ■- режущий круг; 4 —>опорный нож

Рис. 1.23. Закон распределения плот­ ности источника тепловыделения при нагревании поверхности заготовки ши­ рокой газовой горелкой

б)

положить начало координат в точке О, ось Х ш направить вдоль режущей

кромки

(см. рис.

1 .2 1 );

 

 

в)

рассчитать плотность тепловыделения по варианту 1 (foi « 8,98 *104 Вт/м);

г) полагая, что тепловыделение в центре сверла близко к нулю, написать

линейный

закон

распределения в пределах

0 ^

хв ^ R/cos 30° [/ (*и) =

= 1 — хя cos 30°//?];

J (хи)

рассчитать значение

д)

по

формуле (1.52) при известном законе

для второго варианта и далее по формуле (1.49) — плотность потока q02 (q0a «

«

18X104 Вт/м);

0 , написать, пользуясь табл. 1 .1 ,

 

е) имея в виду, что значения ехр [—3] «

экспоненциальный закон распределения

плотности потоков

[/ (хи)] =

=

ехр (—Зхи cos 30°//?)];

 

 

«

ж) по формуле (1.52) рассчитать значения /* для третьего варианта (/* «

0,385) и по формуле (1.49) — плотность потока q03 (q03 ж 26,9*104 Вт/м).

 

14. При врезном бесцентровом шлифовании (рис. 1 .2 2 ) 1 1 ,8 %

мощности

W = 5 кВт расходуется на трение в механизмах станка, а 1 ,5 % — на преодоле­ ние трения между заготовкой и опорным ножом. Окружная скорость вращения заготовки 0 ! = 24 м/мин, сила трения между заготовкой и ведущим кругом Рт= 215 Н. Полагая все источники тепловыделения одномерными, равномерно распределенными по ширине b = 80 мм обрабатываемой поверхности (размер, перпендикулярный к плоскости чертежа), выполнить кодирование тепловой за­ дачи для заготовки и определить плотность теплообразующих потоков в местах

ее соприкосновения с кругами и ножом.

Вт/м2;

Ответ: qoi «

5,31 - 10е Вт/м2; ?0а «0.11-10®

<7оз «

0,094*10® Вт/м2.

широкой газовой го­

15. При нагревании

массивной заготовки пламенем

релки плотность тепловыделения по поверхности пятна нагрева распределена по закону, показанному на рис. 1.23. Кривые нормального распределения имеют коэффициент сосредоточенности k0 = З/ / 2 (см. табл. 1.1). Пятно нагрева медленно движется в направлении v; за пределами этого пятна теплоотдача заготовки в окру­ жающую среду пренебрежимо мала. Представляя источник теплоты в виде трех самостоятельных частей, написать код тепловой задачи для заготовки. Опреде­ лить наибольшую плотность q0 потока, если известно, что за т секунд источник выделил Q джоулей теплоты.

Л 1,95*104Q о / •

0твш : ь ~ Щ ь + оЖ У

Соседние файлы в папке книги