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

книги из ГПНТБ / Тепловые процессы при обработке металлов и сплавов давлением учеб. пособие для студентов металлург. спец. вузов

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

тельный коэффициент теплоотдачи;

Кг

• темпера­

турный критерии внешнего трения;

<*2

(То — Тм)

 

 

оо

 

 

/ " e r f c z ^ j І"-1

erfc yd у.

 

Заметим, что выражение (1.4.13) получено без учета конвектив­

ного массопереноса внутри очага деформации

и в силу этого может

служить лишь для ориентировочных расчетов температуры прокаты­ ваемого металла. На основании этого выражения авторы работы [64] предлагают график для определения температуры поверхности ме­

талла в момент выхода из очага деформации

(рис. 1.2).

Из формулы (1.4.13) несложно определить количество тепла, пе­

редаваемого через

1 м2поверхности

раската:

 

Q

 

к*

 

<*2 0

Тм)

hïat

Vnh У at

1

( 1 — exp h2at erfc h Y at)

 

hïat

 

 

 

0 , 7 52

 

(1.4.14)

 

h V at

h2at

 

 

На рис. 1.3 и 1.4 приведены графики зависимости количества теп­ ла, переданного от металла к валкам в очаге деформации, от вре­ мени при различных зна­ чениях температурного критерия внешнего тре­

ния Кт [64].

5. ВЛИЯНИЕ РАСПРЕДЕЛЕНИЯ ПЛАСТИЧЕСКИХ ДЕФОРМАЦИЙ НА ТЕМПЕРАТУРУ МЕТАЛЛА

0,02 0,04ОМ 0,0д 0,2 Ofi 0,8

hzat

Рис. 1.2. График для определения туры металла в момент выхода

деформации:

1-КТ=0,5; 2 — К т = 0 , 3 ; 3 К т = 0 , 2 ;

5 — К„-0

 

 

Во

время

нахождения

 

 

металла

в

очаге

 

дефор­

 

 

мации

в

 

деформируемом

 

 

объеме

возникает

 

поле

 

 

скоростей

 

сдвиговых

де­

 

 

формаций,

которое

вызы­

 

 

вает

 

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

мощ­

 

 

ности

тепловых

источни­

0,8 1 3 5

7 9ков

в

соответствии

с

вы­

 

 

ражением

 

(1.4.5).

В

об­

темпера­

щем

 

случае

распределе­

из

очага

ние

мощности

тепловых

4 К

т-0,1;

источников по высоте рас­

ката

 

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

Как

 

 

 

20

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

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

Рассматриваемую за­ дачу сформулируем сле­ дующим образом. Неогра­

ниченная

пластина

тол­

щиной 2R,

имеющая

по­

стоянную по объему

тем­

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

Перечисленные

 

условия

действуют

в течение

вре­

мени

0 ^ ^ 2 .

 

Прини­

мая, что теплообмен

меж­

ду

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

пласти­

ны

и

валком

происходит

по

закону

Ньютона

при

коэффициенте

теплоотда­

чи а, определим

темпера­

турное поле пластины в указанном интервале вре­ мени (т. е. в течение вре­ мени деформирования).

0,02 ОМ 0,0В 0,0В 0,1 ОЛ 0,80,81 3 5 7 3 hzat

Рис. 1.3. Количество тепла, переданное вал­ кам в очаге деформации п — 0):

/ - / С т = 0 ; 2 - К т = 0 , 1 ; 3 - К Т = 0 , 2 ; 4 - К т =0,3;

5 — Кт = 0,5

і,0

0,8

*Г 0,0

І-~

0,2

ь--

0,02 ОМ 0,0В 0,08 0,2 OA 0.S0,81 3 5 7 3 h'at

Рис. 1.4. Количество тепла, переданное вал­

кам

в очаге

деформации

(/Сд = 0,05) :•

1-КТ=0;

2 - І С І

= 0 Д ; 3-Kr=0,2;

4 - Кт=0,3;

21

Д л я этого дифференциальное уравнение теплопроводности

 

 

 

dFo

дХ2

 

 

 

 

(1.5.1)

 

 

 

 

 

 

 

 

 

решаем при следующих краевых условиях:

 

 

 

 

 

 

дѵ

 

=

0;

 

(1.5.2)

 

 

 

~дХ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дѵ

= В і [ х » в - г ; ( 1 , Fo)] f

KU

(1.5.3)

 

 

дХ

 

 

 

 

 

 

 

 

 

 

 

 

 

^ l F o = 0

= ^o =

 

const.

 

(1.-.4)

Все величины в системе уравнений

(1.5.1) — (1.5.4) выражены

в

безразмерном виде, а именно:

 

 

 

 

 

 

 

-v(X, Fo) — T(X'

— безразмерная

температурная функция; Р о Л

=

 

^о — Тв

 

 

 

 

 

 

 

 

 

— критерий Померанцева;

X коэффициент теплопро-

водности;

В і = — R

—критерий

Био;

К і ^

 

критерии

Кирпичева;

Fo —

 

критерий Фѵрье; а—

коэффициент тем-

пературопрозодноещ;

с — удельная

теплоемкость;

р —плотность;

Х=-^—относительная

 

координата (— 1 -< X -< 1); А; координата.

К дифференциальному уравнению (1.5.1) применим интеграль­ ное преобразование Лапласа — Карсона [32—40], после чего полу­ чим неоднородное дифференциальное уравнение, лишенное диффе­ ренциальных операций по времени:

ю"{Х,

р)-рѵ(Х,

p)=-pv0-PoRX2,

;i.5.5)

где

оо

 

 

 

 

 

v{p)

= P J г> (Fo) exp ( — p Fo) d Fo,

; 1.5.6)

 

p — комплексный параметр.

Граничные условия (1.5.2) и (1.5.3) запишутся после преобра­

зования:

 

 

 

 

 

dv

=

0;

(1.5.7)

 

dX |x=o

 

 

 

 

dv

= Bi[vB-v(\,

 

/7)] + Кі.

(1.5.8)

dX

 

x-i

 

 

 

22

Общее решение дифференциального уравнения (1.5.5), удовлет­ воряющее условию (1.5.7), запишется следующим образом:

v(X, p ) - V o = Ä(p)chVpX

+ ^ - X 2

+ ^ - . (1.5.9)

 

Р

Р2

Величину Ä(p), постоянную относительно переменной X, опре­ делим из граничного условия (1.5.8):

Кі -

Ві - fÜ* (Bi + 2) -

~

Рор

 

А(Р):

Р

Р2

к

(1.5.9а)

Vpsh Ѵр + Ві сьУ р

Отсюда решение поставленной задачи в области изображений можно записать как

ѵ(Х,

р)-ѵ0--

chVpX

( K l — B l )

Po (2 + Bi)

 

5

 

 

1/^/7 sh V p + Bi c\iV

P

 

 

 

2 P 0 / ? B i

2Por

(1.5.10)

 

 

 

 

Принимая во внимание то обстоятельство, что прокатываемый металл находится в очаге деформации доли секунды, имеет смысл искать решение поставленной задачи, удобное для вычислений при малых значениях критерия Fo. Следуя [32], получим в этом случае решение в области изображений:

ѵ(х,

р ^ ^ ^ - У р ^

- ^ +

^ Л - Ѵ р і ^

X)]

х

 

X

(Ki — Bi)

Po D

 

2Po D Bi

P o „

 

2Po D

(1.5.11)

5. (2 + BI)

S—

+ _ * A"2_|

 

 

 

P

 

P2

P

 

P2

 

Переходя обычным методом от решения в области

изображений

к решению в области действительной переменной, получим

выраже­

ние для функции

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

температуры по толщине

пласти­

ны [41]:

 

 

 

 

 

 

 

 

v0-v(X,

F o ) « ( l - | j - ) < I > ( * ,

Р ° ) + Р

0 * в , +

ф ( ^ ' W

+

 

 

 

 

 

 

о

 

 

 

Fo

 

 

 

 

 

 

 

+

2 Р о Л J (Fo-t)O(X,

t) dt - Pop* 2 Fo -

Po^Fo2 ,

(1.5.12)

23-

где

Ф(Х,

F o ) - e r f c - b ^ L - e x p [ ( l - A ' ) B i + ß i 2 F o ]

X

 

 

2 ]/FO

 

- e x p [ ( l f

 

X e r î c (

^ ^

+ ВіІ^РсЛ -I-erfc

1

Х)Ш±

\2"|/FO

/

2|/FO

 

 

 

+

Bi2 Fo] erfc / 1 + _ f - f Bi Wo) .

(1.5.13)

 

 

\ 2 K F O

 

/

 

Область практических значений критерия Кд при прокатке на обжимных станах не превышает 1-=-2. Учитывая высокие в этом случае значения критерия Био, можно считать, что

^ - « 1 .

(1-5.14)

Ві

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

зультате работы против внешних сил

трения, можно

пренебречь.

В то же время при прокатке металла

небольшой толщины (напри­

мер, листа) тепловой поток дтр следует

учитывать.

 

Используя решение (1.5.12), рассчитали безразмерную темпера­

туру поверхности раската.

Некоторые

из полученных

результатов

представлены на графиках

(рис. 1.5—1.7). На этих графиках по вер­

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

 

ѵ0-ѵ(Х,

Fo) = 0(A", Fo, Ві, Р о д , КО,

(1.5.15)

а по горизонтальной — значение критерия Fo. Сплошной линией изображается температура раската при учете неравномерности рас-

24

пределения тепловых источников по высоте. Штриховой линией — температура поверхности, рассчитанная при условии, что тепловые

0,700

0,800

 

 

/

 

 

-- —

0,500

 

.

^ 0,400

 

 

г _

S.

— • —

 

 

 

— "

 

£ 0,300

 

 

 

в 0,200

te

1

 

11'/

s**

 

 

5

0,100W

 

 

 

 

Fa-/0 s

/0

 

 

 

Рис. 1.6. Изменение во времени безразмерной темпе­

ратуры поверхности

раската ( Р о я = 5 0 0 ) :

/ — Ві = 100;

2 — Ві=80; 3 — Bi=60; 4 — Ві=40;

S — Ві=20

источники распределены равномерно по высоте. Анализируя приве­ денные графики, можно заключить, что разность между значениями

этих функций может

достигать

50%. Далее, из анализа решения

0.700• 1

1 , 1

1 1 1 1 1 ,

!

Z 3

4

5 В

7

в

9 W

 

 

 

 

Ft

ІО5

 

Рис. 1.7. Изменение во времени безразмерной темпе­

ратуры

поверхности

раската

(Род = 1000):

/ _ В І = 100;

2 — Ві - 80;

3 — Ві=60; 4 — Ві=40;

5 — В і = 2 0

(1.5.13) следует, что учитывать

в расчетах

температуры неравно­

мерность распределения мощности тепловых источников по высоте раската следует в тех случаях, когда Р о д > 2 0 0 .

25

6.ПОСТАНОВКА ЗАДАЧИ ТЕПЛОПРОВОДНОСТИ ДЛЯ ДЕФОРМИРУЕМОГО МЕТАЛЛА

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

описываться

дифференциальным

уравнением

Фурье — Кирхгофа

[6,

42]

 

 

 

 

 

 

 

 

 

i r =

a ( ^ + ^ + ^ i - ) - ( s

^ + S 2

^ - + 5

3 ^ - ) + 7 p - '

( 1 - 6 Л )

где

 

Т(х, у, z,

t)

— функция

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

температуры

по

объему

деформируемого

металла; su

s2, s3

— компоненты

скорости

течения

частиц металла

по координатным

осям х,

у и

г;

W(x,

у,

z, t) —

функция распределения мощности тепловых источников по объему металла.

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

Прежде чем написать условия однозначности для дифференци­

ального уравнения

(1.6.1), остановимся на способах

аналитической

записи функции коэффициента

теплоотдачи

a(t)

(на примере про­

дольной прокатки

металла).

 

 

 

 

 

 

Как упоминалось выше, функция a(t)

имеет различные значе­

ния для пауз (ai)

и обжатий

(аг). Помимо

этого, коэффициент

аі

 

 

 

 

 

 

Т а б л и ц а

1.1

Значения коэффициента гхі в чистовой группе стана

 

1700 при прокатке полосы 2,5x720

мм

 

 

 

GCJ

для начала

az

для

конца

Ошибка в резуль­

Номер клети

тате усреднения

 

паузы,

 

паузы,

коэффициента

а,

 

ккал/м2 'Ч-град

ккал/м'-ч- град по пропускам,

%

VI

1,04

116,8

 

112,6

9,55

 

VII

 

 

2,00

111,4

 

109,4

5,45

 

V I I I

 

 

3,40

108,7

 

105,7

2,38

 

IX

 

 

4,89

103,0

 

101,0

2,58

 

X

 

 

6,72

99,5

 

97,5

5,92

 

XI

 

 

8,22

95,6

 

93,6

9,65

 

 

 

 

26

изменяется в течение прокатки (поскольку изменяются темпера­ тура поверхности и скорость движения раската). С другой сторо­ ны, коэффициент а2 изменяется в течение времени обжатия и зави­ сит от номера пропуска. Для иллюстрации сказанного приведем численные значения коэффициента си при прокатке полосы сечени­ ем 2,5X720 мм в чистовой группе стана 1700. Температурный ре­ жим раската принят по результатам работы [11]. Значения ai полу­ чены расчетным путем и приведены в табл. 1.1.

Как следует из табл. 1.1, коэффициент а\ зависит от номера про­ пуска и изменяется в течение каждой паузы. В то же время усред­ нение коэффициента ai по пропускам обусловливает погрешность, не превышающую 10%- Поскольку точность определения парамет­

ров

внешнего и внутреннего

теплообмена

обычно

не

превышает

10—15%, целесообразно

 

использовать

в расчетах

усредненное

по

пропускам значение коэффициента оц.

 

 

 

 

 

 

Аналогичный анализ, проведенный для коэффициента а2,

пока­

зывает, что и в этом случае возможно применение в расчетах

усред­

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

а2.

 

 

 

 

 

 

 

 

 

 

Таким образом, для описания функции a(t) можно использовать

единичную функцию [1, 32, 43—45]:

 

 

 

 

 

 

 

 

а (/) = с^-Наа о ^ М О ,

 

(1.6.2)

 

где q>(t) имеет следующие значения:

 

 

 

 

 

 

для (п+ 1)-й паузы

 

 

 

 

 

 

 

 

 

 

 

 

т = 1

 

 

 

 

 

 

 

 

 

 

для

(п-\-\)-то обжатия

 

 

 

 

 

 

 

 

 

 

 

9(0

= 2

^

-

^

0 - ^ -

2

4 ( ' - ' " ' o ) ,

 

(1-6.4)

 

 

т=0

 

 

 

 

т — 1

 

 

 

 

а единичная функция ц (z)

определяется

как

 

 

 

 

 

 

 

 

 

(0

при

z<T0

 

 

 

 

 

 

 

ѵ

 

(1

при z

>

0,

 

 

ѵ

 

причем t\ — средняя за

время

прокатки длительность

паузы;

t2

средняя длительность обжатия *;

 

 

 

 

 

 

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

27

Выражение (1.6.2) будет правильно описывать функцию

a (t) в

тех случаях, когда, с одной стороны, коэффициенты

сц и а2

в тече­

ние прокатки изменяются

незначительно, а с другой стороны, дли­

тельность пауз и обжатий

(т. е. величины ti и t2)

одинакова для

каждого пропуска. Если перечисленные условия не выполняются, то

математическое аппроксимирование

функции

a(t) принципиально

не усложнится, в чем и заключается

большое удобство применения

единичной функции r\(z).

Например, изображенная на рис. 1.8 кри­

вая представится следующим образом:

 

o.(t) = ali-{a2~al)ri(t~tl)

— (a2 — a[)rl(t

— t1 — t2) - f

+ (а'2 аі)т (

(t tx —12—t\) — (a2

al) X

X y l ( t - t l - t i - t ' l - t l ) + . . . .

(1-6.6)

oc(t)\

 

 

 

CM

 

T3

T3

 

* CM

TS \?

TS

T5

 

 

 

 

"raj

"T?

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Г/

 

t;

t;

tf

t"

t;

.///,

 

 

1 R 2

 

Рис. 1.8. Зависимость коэффициента теплоотдачи от времени

 

Аналогично запишем функцию

теплового

потока через

поверх­

ность прокатываемого

металла:

 

 

 

 

 

 

 

?(0 =

?і +

( ? 2 - ? і М ' ) ,

 

(1-6.7)

где qi — среднее за прокатку

значение удельного теплового

потока

через поверхность раската в течение пауз; q2 — то ж е в течение вре­ мени обжатия.

В течение времени обжатия в деформируемом объеме действуют тепловые источники мощностью W2. В течение пауз мощность ис­ точников равна нулю. Следовательно, функцию W(t), входящую в дифференциальное уравнение (1.6.1), можно записать следующим образом:

W{t) = W&{t).

(1.6.8)

На протяжении процесса прокатки поверхность раската отдает тепло в окружающую атмосферу (в течение паузы) и в валки (в те-

28

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

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

Ось раската

Обжатие

Пауза

Рис. 1.9.

Изменение во времени положения поверхности раската:

Re — половина

начальной толщины заготовки;

/?к —половина конечной толщины

 

раската; і, — деформация

за один пропуск

его поверхности. Если ограничиться рассмотрением изменения по­

ложения поверхности раската только в одном направлении

(напри­

мер, в вертикальном), то схематически процесс может быть

изобра­

жен, как показано на рис. 1.9. Функция f(t) представляет

собой

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

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

поверх­

ности раската. При таком законе перемещения поверхности

решение

дифференциального уравнения (1.6.1) связано с большими

трудно­

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

Используя изложенные допущения, запишем краевые условия для дифференциального уравнения (1.6.1) в случае граничных усло­ вий второго рода [1, 6, 35, 42]:

. дТ дх

дТ

ду

дТ dz

= Ял (*).

(1.6.10)

-Яь С),

(1.6.11)

-Я ad),

(1.6.12)

29