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

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

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

Выполняя в выражении (2.2.35) интегрирование во внутренних интегралах, а затем интегрируя полученные функции по частям или подстановкой, находим

^ erfc

X

м

 

 

 

 

2V

 

 

 

 

 

 

t

 

 

 

 

 

 

2_У± exp

X

M

y+

 

 

 

YnM

2 Ѵ Ч

 

 

 

 

e x p ( - A ! X ) e r f c / _ X

 

М у -

(2.2.36)

 

 

 

 

 

 

Подставляя^ (2.2.36) в (2.2.33), получаем

X

 

 

V*

 

exp(-MX)cdc

I X

M

(2.2.37)

 

 

 

Ш

 

2

 

Дл я неподвижной системы отсчета эта функция запишется сле­ дующим образом:

Ф (Xv

t ) =

?^4:exp

 

X,

2Уг

 

 

У л

— (— 4-Х, —Мл)

erfc/

 

921

м

s

\ 2

у

ехр [ — М(Х1

М[х)] erfc

X:А_

-і-(м[+м)Ух . (2.2.38)

 

2 ] / ,

2

Если предположить, что деформации по сечению раската рас­ пределяются равномерно, то можно принять М = 0. Раскрывая не­ определенности, появляющиеся в этом случае в выражениях (2.2.37) и (2.2.38), находим:

Н т Л . 0

[ Ф ( * . 1 > ) ] = ^ е

х Р (

- - ^

X erfc

X

; (2.2.39)

 

 

2 ^

 

 

 

 

 

 

Ііш М-+0

[ф №. *)]=1/я

ехр

 

Ух

 

 

 

(Л^—Л*іт)егіс

X,

м\

1/t

 

(2.2.40).

 

 

 

2]А

 

 

 

40

sfc= 5,0 6,0 г
.0,40
0,20

Функция (2.2.39) представляет собой решение дифференциаль­ ного уравнения теплопроводности для полупространства с непод­

вижной

границей

и граничными

условиями

I I рода. Табулирован­

ные значения

этой

функции

2,00

 

приведены

в

приложении 1.

 

 

На

рис. 2.3 и 2.4

представле­

1,80

 

ны

графики

функций

Ф(Хи

1,60

 

т)

для

подвижной

границы

 

 

0,1

(т. е. Х\ = М/х)

при

различ­

IAO

 

ных значениях M и времени

1,20

 

(принято Мі = 0).

 

 

 

1,0

 

Д л я

окончательного про­

V w

 

 

ведения

обратного

преобра­

 

 

зования

выражения

(2.2.29)

0,80

1,5

используем

теорему

запаз­

 

0,50

 

дывания, которая

состоит в

2,0

следующем:

 

 

 

 

 

 

 

-Л.П-

XX*-У,

(2.2.41)

где г)(т — Ь)

единичная

функция, определяемая со­ отношением (1.6.5).

і.во

CD 1I

 

1

1,40

i l

 

 

сэіІІ

',20

1,5 2,0

 

0,80

0,8 1,0

 

 

(C: 4 \

0,60

 

0Л0

/>>

 

CSJ

 

11

0,20

 

'

О 2

S

8

10

12 14- Т

Рис. 2.3. Зависимость

функции

Ф ( М і ' т, X) от времени

при

различных

 

значениях

M

 

 

а

Рис. 2.4. Зависимость функции

ФЩ/

т, т)

от величины M в различные

моменты

вре­

мени

 

 

41

1-я пауза

( О ^ г ^ ^ т і ) :

 

 

 

 

 

решение в области изображений

 

 

 

 

Ѵ(Х,

/>) =

КііФ(Л", Р)\

 

(2-2.42)

решение в области действительной переменной

 

 

 

 

Ѵ(Х,

ф) =

Кі 1 Ф(Л', Ф);

 

(2.2.43)

/-е обжатие ( т і ^ а і ^ т г о ) :

 

 

 

 

решение в области изображений

 

 

 

 

X®(X,

р)

^ і - е х р ( - / л г 1 ) ;

 

(2.2.44)

 

 

 

 

 

Р

 

 

решение в области действительной переменной

 

 

Ѵ(Х,

^) =

Ki^(X,

 

ф) +

( К і 2 - К і 1 ) Ф ( А ' , Ф - t O X

 

X

(Ф -

xj -

Р о 2 (ф - 1 , ) Tj (ф - 1 , ) .

 

(2.2.45)

2-я пауза

( т о ^ г | ) ^ т 0 + Т і ) :

 

 

 

решение в области изображений

 

 

Ѵ(Х, р)=Кі1Ф(Х,

/?) +

( K i 2 - K i 1 ) [ e x p ( - / ? t 1

) X

 

X®(X,

р)-ехр(-рх0)Ф(Х,

 

pïï-^expi-pxj;

 

(2.2.46)

 

 

 

 

 

р

 

 

решение в области действительной переменной

 

 

Ѵ(Х, ф) = К І ! Ф ( ^ , Ф Ж К І 2 - К і , ) [ Ф ( Л \ ф _ ^ ) т і ( Ф — r t ) -

_ Ф ( Л \

ф - г 0 ) 7 ] ( ф - г 0

) ] - Р о 2 ( ф - г 1 ) 7 1 ( ф - г 1

)

(2.2.47)

ит. д.

Вобщем случае будем иметь:

для (п + 1)-й

паузы

 

 

 

 

Ѵ(Х, ф) = К 1 х Ф ( ^ ,

Ф) + (КІ2 -Кі1 )2 < ф ^ '

ф - ( « - і ) х

. X t o - t J ТІ[Ф—(m— l)tr 0 — ti] — Ф (^Г,

* - m t 0

) i j ( * - m t 0 ) } -

- Po2 2

( m - 1) * o - * i l 1

( m - - ) r o - r i l -

m =

l

 

 

 

 

 

-

(ф -

mt0 ) -q (Ф— mt0 )} ;

(2.2.48)

для (п + ^)-го

обжатия

 

 

 

^(Х,Ф) = К і 1 Ф ( ^ , Ф) +

( К і а - К І ! ) 2

Ф ( ^ ,

Ф - о т Ѵ - t j x

 

 

 

m-0

 

 

Х т і ( Ф - т * 0 - ^ ) - 2 ф ( ^ > Ф - ^ 0 ) і ) ( Ф - т г 0 ) ] - т - 1

42

— Po,

л

(2.2.49)

m-l

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

0,32

0,28

0,24

0,20

z.0,10

^ 0,12

t

0,08

ОМ

0,00

 

-0,0t

 

 

 

 

 

 

л

 

 

 

 

 

 

 

ІБ IB 20 T

 

0

2

U

6

W

12

m

 

Рис. 2.5. Динамика

изменения вс

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

 

 

 

верхности / и центра раската

/ /

для

(п+1)-й

паузы

 

 

 

 

 

ѵ(Хѵ

т) = Кі1 Ф(А'1 ,

т)4-(Кі2 -КІ!)2

{ ® [ ^ 1 - c - ( m - l ) T 0 - t 1 ] X

 

 

 

 

 

 

m =

l

 

(Зля v(Xv

X Ч [t (m — 1) t 0 — t j — Ф (Xv

x - mt0 ) Kj (t mt0 )} —

- P o 2 2 { [ t - ( m — l J t o - ^ T j l t - î m - l ) t 0 — t j —

nt-1

 

 

 

 

— (t — wt 0 )7](t — m t 0 ) } ;

(2.2.50)

(п + 1)-го обжатия

2 ф ( * і «

f - m t o - t J X

x) = Kh®(Xv

т) + ( К 1 2 - К ч )

 

m - 0

X ч (t m t 0 — t j ) - 2 Ф

t - mt 0 ) 7) (t — mx0)

m = l

 

43

"

ft-

 

 

— Po2 2 (t~mx0

— х1)уі{х— тхй— xj — ^

(t — mt 0 )т](t — mt0 )

_m=0

m =

l

_

 

 

 

(2.2.51)

На рис. 2.5 приведена динамика изменения во времени темпера­ туры на поверхности / и в центре / / раската в течение 7 пропусков. Д л я расчета приняты следующие значения параметров теплооб­ мена:

В і 0 = 1 0 ;

Л =

330 \\м;

К ч =

0,02;

Кд2 =

0,70;

t 1 = 2 ,

5 0 ;

* а =

0,05;

t 0 = 2 , 5 5 ;

Л1І =

Л11 =

3,64;

Р о 2 =

0,04.

 

Пример

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

Безразмерное

время,

соответствующее концу третьей

паузы,

равно 2то+Ті =

=2(2,50+ 0,05)

+2,50=

7,60.

 

 

 

 

 

 

 

Д л я расчета

используем решение

(2.2.50),

которое

в данном случае

запи­

шется так:

 

 

 

 

 

 

 

 

 

 

 

V

[M[z,

Т ) =

0.02Ф (7,60) + 0,68 [Ф (5,10) — Ф (5,05)

+

 

 

 

 

 

 

Ф (2,55) — Ф (2,50)] — 2 Р о 2 т 2 .

 

 

 

Функция

Ф(т)

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

выражением

(2.2.40) после подстановки

значе­

ния координаты

расчетной

точки

Хі =

М\'х:г

 

 

 

 

2 " ^ т

V я

Вычисляя необходимые значения функции Ф ( т ) , получаем

v(M'z, т = 7,60) = 0,02-3,110 + 0,68(2,258 — 2,247 +

+

1,597— 1,581)—2 0,040-0,05 = 0,082 —0,004 =

0,078.

 

Таким образом, безразмерная температура поверхности

раската

в конце

третьей паузы

составляет 0,078.

 

 

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

44

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

дифференциального

уравнения

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

/

dv(Xh

Yh

х)

д2ѵь

Yu

x) ,

дѢ(Хх, Yh

x)

 

 

 

 

дХ\

f

 

 

 

 

 

 

ÔY\

 

• M ,

àv(Xb

Ki,

T)

•Mo

dv(Xh

Y1

•Po(t)

 

(М[х <

 

< оо,

A f 2 t < r ! < o o ) ,

(2.2.52)

при следующих краевых

условиях:

 

 

 

 

 

дѵ

 

 

- K i ' ( t ) ;

(2.2.53)

 

 

дХі

 

 

 

 

X^Mf

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дѵ

 

 

- K i " ( t ) ;

(2.2.54)

 

 

ÔY

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.2.55)

 

 

дѵ(Хъ

Ylt

x)

 

0;

(2.2.56)

 

 

 

 

дХг

 

 

 

 

 

 

 

X,=oo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.2.57)

 

 

 

 

 

 

 

= 0:

(2.2.58)

 

 

 

 

 

 

К, = oo

 

 

 

 

 

 

 

 

 

(2.2.59)

где Ро(т) — критерий Померанцева, определенный как суммарный эффект диссипации механической энергии формоизменения по на­ правлениям Хі и Уі; Кі'(т) —критерий Кирпичева, определенный по

тепловому потоку, который проходит через

поверхность Х\ — М\'х;

Кі"(т) — критерий Кирпичева, определенный

по тепловому

потоку,

который

проходит через

поверхность Уі = М2 Ѵ, M] — компонент

безразмерной скорости течения частиц металла

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

Х\\

М2 — то же в направлении

Уь М/—безразмерная

скорость

движе­

ния поверхности раската в направлении Х\\

М2'

— то же в направ­

лении У].

 

 

 

 

 

Под

функцией ѵ(Хи У ь

т) в выражениях

(2.2.52) — (2.2.59)

сле­

дует понимать температурное поле области,

образованной

пересе­

чением двух полуограниченных тел. Начало координат этой области

расположено

в вершине двугранного

угла

XiQY]

(рис. 2.6).

При

этом считаем,

что призма /—2—3—4 —4'V2'—3'

с точки

зре­

ния теории теплопроводности ведет себя в

направлениях Х\

и Y\

как полуограниченное

тело. Это положение вытекает из обсуждав­

шегося выше

вопроса

об инерционном

времени, которое при про-

45

катке на обжимных станах больше времени прокатки. Зная темпе­ ратурное поле призмы /—234 —4'—/'—2 3', мы имеем всю информацию о температурном поле прокатываемой заготовки (вследствие симметрии).

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

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

(штриховыми линиями показаны контуры реальной заготовки, имеющей высоту 2 R и ширину 2 s)

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

 

Х = Х1 М1х,

 

 

У=У1—М2х,

 

ф = т .

 

(2.2.60)

Принимая во внимание, что

 

 

 

 

 

 

 

 

 

v(Xv

Ѵъ

 

х) = Ѵ(Х,

V,

ф),

 

(2.2.61)

а также следующие соотношения [63]:

 

 

 

 

 

 

дѵ

дѴ

дХ

, дѴ

дѴ

1

 

 

д\>

дѴ

(2.2.62)

дХг

дХ

дХх

 

dY

дХх

д Ѵ

 

дХ

дХ

1

1

аф •

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д^ѵ

-д

 

х

(—]

 

dW

 

 

 

(2.2.62а)

 

 

дХ\

дХ

 

[дх)

 

 

 

 

 

 

дѵ

дѴ

дХ

, дѴ dY

 

дѴ

 

 

дѴ

(2.2.626)

 

дХ

дГ1

1 дг

 

дУі

 

дь

дУі

дГ

 

 

1

 

 

 

д2ѵ

 

д

(дѴ\

_

д?Ѵ

 

 

(2.2.62B)

 

 

дѴ\

 

 

 

[дУ

) ~

дУ2

 

 

 

 

 

 

 

 

 

 

46

дѵ

дѴ

дХ

,

дѴ

дГ .

дѴ

сіф

дг

дХ

dz

'

дГ

дг

1

d>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дХ

 

•^2 —

+ —

,

(2.2.63)

 

 

 

 

 

 

 

 

 

запишем систему

уравнений

(2.2.52) —(2.2.59)

следующим обра­

зом:

 

 

 

 

 

 

 

 

 

дѴ(Х,

Y, ф) _

dW(X,

Y, ф) t

dW(X,

 

Y, ф) |

<Эф _ dV(X, Y, ф)

dX

0 < X^oo,

dX

дѴ dY

dX2

 

dY2

 

•My

dV(X,

Y, ф)

•Ро(ф),

 

 

'dY

 

0 <

Г < о о .

(2.2.64)

=

-

K i ' (Ф),

(2.2.65)

=

- К і " ( Ф ) .

(2.2.66)

 

 

 

(2.2.67)

дѴ

= 0,

(2.2.68)

дХ

Х= (

 

 

 

Ѵ{Х,

Y, ф ) | г = с о # с х э ,

(2.2.69)

дѴ

- О ,

(2.2.70)

dY

 

 

^

Ф)Іф_о = ° .

(2.2.71)

 

где МХ = М[ — М{, МУ — М'іМ2.

Подвергая уравнение (2.2.64) интегральному преобразованию Лапласа — Карсона, получаем дифференциальное уравнение в част­ ных производных по X и Y, лишенное дифференциальных операций по времени:

д2Ѵ(Х, Г , р) •

dW(X,

 

Y, p)

•M,

dV(X,

Y, p)

 

дХ2

 

 

dY2

 

dX

 

 

 

 

 

 

-My

дѴ(Х,

Y,

p)

•pV(X,

V, p) = Po(p).

(2.2.72)

 

dY

 

 

 

 

 

 

 

 

 

Граничные условия в этом случае получают следующий вид:

дѴ

= - К і ' (р),

(2.2.73)

дХ

х=о

 

 

 

47

 

дѴ

 

 

 

 

(2.2.74)

 

дГ

іу=о

 

 

 

 

 

 

 

 

 

 

 

 

Ѵ(Х,

У,

 

p)\x,ao^œ,

(2.2.75)

 

 

дѴ

 

=

0,

(2.2.76)

 

 

дХ

 

 

 

ІХ= оо

 

 

 

 

 

 

 

 

Ѵ(Х,

У,

 

рПу^фоо,

(2.2.77)

 

 

дѴ

 

•0.

(2.2.78)

 

 

дѴ

 

 

 

 

 

 

 

Общим решением дифференциального уравнения (2.2.72) яв­

ляется функция

 

 

 

 

 

Ѵ(Х, У,

p)=Ä(p)exp

 

 

 

+

 

 

 

 

+

ß ( p ) e x p

 

 

 

+

 

 

 

 

 

 

+

С(/?)ехр

У

M,

 

 

 

 

 

 

 

+ 7 3 ( / ? ) е х р [ г |

 

 

Mi

PO (p) .

(2.2.79)

 

 

 

Удовлетворяя общее решение (2.2.79) граничным условиям

(2.2.73) — (2.2.78), находим:

 

 

 

 

 

 

В(р)

= 0,

D(p)

= 0;

 

 

А(р) =

 

 

K î ' (/>)

 

 

 

M X

+

 

 

 

 

 

2

У

4

 

 

 

'

 

С{р):

КГ" (/>)

 

 

Y

48

Таким образом, решение задачи в области изображений запи­ шется следующим образом:

ехр

X

Af,

V

Mi

+

р

 

 

 

 

 

 

 

 

 

 

 

Ѵ(Х, V, р) = КѴ{р)

 

+ / 4

 

 

 

 

 

M,

+

р

 

 

ехр

 

мі

+ р

 

 

 

 

 

 

 

Po

(p)

(2.2.80)

- fKi"(/>)

 

 

 

 

 

 

 

 

 

 

Af,.

Ml

 

 

 

 

 

 

+

 

+ P

 

 

 

 

 

 

 

 

 

 

 

 

Обозначив

 

 

 

 

 

 

 

ехр — X

AT*

 

 

 

 

 

 

 

 

 

+

P

 

 

Ali

 

 

Mi

 

 

 

 

 

 

 

+

P

 

 

 

M,,

/ 4

 

 

 

 

 

 

 

 

 

 

exp

 

 

 

 

 

 

 

< W P)-

 

 

 

 

 

 

 

A l ,

 

Mi

 

 

 

 

 

+

 

+ p

 

 

 

 

 

 

 

 

 

 

 

получим после обратного преобразования решение в области дейст­

вительной

переменной:

 

для (п +1 )-й

паузы

 

Ѵ(Х,

У,

ф) = КіІФ,(А-> ф) + (КІ2—Кіі) 2

< ф і [ * ' Ф -

 

 

m =

l

- ( т - 1 ) т 0 - т 1 ] т ) [ 1 ) ; - ( т - 1 ) т ( ) - т 1 ] - Ф 1 ( Х , ф - mt0 ) т) ( ф - т т 0 ) ) +

+ К І ; Ф 2 ( К , ф) + (КІ2 — К І І " ) 2 l ) t 0 - t j X

m-l

Х ч [ ф — С " — 1 ) ^ о - ^ і ] - ф 2 ( ^ . ф—/пг0 |(ф mt0 )> —

- ( ф - т т 0 ) т ) ( ф - т т 0 ) } ;

(2.2.81)

49