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

книги из ГПНТБ / Шичков А.Н. Температурный режим листопрокатных валков

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

Попутно следует отметить, что если вращение отсутствует, т. е. <»0 = 0, то уравнение (1.31) примет вид

&Ап

1 дАг

■Ап = 0.

(1.35)

öp2

др

 

 

Общим решением этого уравнения будет А п (р) =

Схр", где при

р = 1 Сі = Л^„; тогда

 

 

 

 

Л„(р) = ЛГпр".

(1.36)

Выражение (1.34) имеет комплексную форму записи. Выделим в нем мнимую и действительную части. Из теории Бесселевых функ­

ций известно [7, 8,

21, 36, 40, 54, 95], что

 

/„ (]/7 P d

p) = ber„(]/r Pdp) + ib e i„ (]/P d p ) ,

(1.37)

где Ьег„ и Ьеі„ — функции Кельвина (Томсона) n-го порядка. В то же время известно, что

еі PdFoJ = cos (Pd Fo) +

i sin (Pd Fo).

 

(1.38)

Подставляя (1.37) и (1.38) в (1.34),

получим

 

 

А ( о ) — N

Ьег” { У Pd р) + i bein ( У Pd р)

 

"

Ьегп ( К м ) + г-Ьеіл (У Й )

Х

 

X [cos (Pd Fo) i sin (Pd Fo)].

 

(1.39)

Выделим в выражении (1.39) действительную

и мнимую части

и в результате преобразований найдем

 

 

Л„(р) =

^ п[КеИп(р) + ; Im (р)].

 

(1.40)

Модуль этой функции будет равен

 

 

 

\Ап (Р) I = Nn Y ReM„(p) + ImMn(p) ,

(1.41)

где

ReH„(p) =

 

 

 

 

 

[ber„ { У Pd р) ber„ ( У Pd) + bei„ [ У Pd p) bei„ t y Pd)] X X cos (Pd Fo) — [ber„ ( ]/P d ) bei„ (K Pd p) —

____________— bein { У Pd) ber„ [ У Pd p)] sin (Pd Fo)________

Ьег^ [ У Pd) + beij] [ У Pd)

Im A n (p) =

[ber„ ( У Pd p) ber„ [ У Pd) + bein (]/p d p) bei„ ( У Pd)] X X sin (Pd Fo) — [bei„ ( V Pd) ber„ [ У Pd p) —

___________ -- bern [ У Pd) bein { У Pd p)] cos (Pd Fo)_________

ber^ І У Pd) + bei^ [ У Pd)

(1.42a)

(1.426)

20

Уравнение, которому удовлетворяет функция А*п (р, Fo), будет

д2К .

1

дК

п2 л *__ дА*а

(1.43)

dp'-

Р

’ др

р2 " д Fo

 

Общее решение этого уравнения имеет вид

Ап(Р, Щ = ^ С п^ п ( ^ Р ) е

^ fcFo,

(1.44)

ft=i

 

 

где Jп (рлАр) — функция Бесселя 1-го рода,

п-го порядка.

 

Характеристическое уравнение в этом случае будет

 

К (н-nft) = 0.

 

(1.45)

а коэффициент Cnk, как и в предыдущем случае, находится по вы­

ражению

1г

 

Іп { Ѵ f Pdp) ßi PdFo

Jn (рлй p) P^P

о

In (KFPd \

(1.46)

Cnk —

 

 

 

\ p j l ( p nkp)dp

 

 

о

 

Окончательное выражение для функции Ап (р, Fo) с учетом (1.29), (1.34),- (1.44) и (1.46) будет иметь вид

2 ft=i

Л,«(р, Fo) = N neI PdFo

7„(K fPdp)_

 

 

 

In [ V i Pd)

 

p

/„ [ V i Pd p r ,

 

 

j P

-----J n (.PnkP)dp

 

n

Irn (F T p d )

~V-nk?°

• (1-47)

 

 

Jn (ilnfe P) e

J P ^ ( ^ p ) dP

Выражение для коэффициента Bn (р, Fo) окажется "абсолютно таким же, как (1.47), только вместо коэффициента Nn будет коэф­ фициент М п.

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

и (1.47):

Ѳ(р, Ф, Fo) = 0,5Л0 (р, Fo) + 2 А п (р, Fo) cos (пф) +

п—1

+ ß„(p, Fo) sin (wp) = Ѳ—

21'

-2j'[0 — Ѳ0] -MPoftP) Pdp

J 0 (pOjfcP) e - ^ F0+

J P^ot^oftP) dP

0

CO

 

p I nW £ Pd p) r ,

, j

In (іЛ-Pdp)

- - [ ■/■----- \ L Jn (P/»feP) pdp

i In [ V i Pd)

X

l n (lA'Pd)

1

 

 

J PJl (Pnk P) dP

0

X / „ ( |W ) e ^ Ро] e‘ PdF°[A£„cos(n(p) + M„sin(ncp)]. (1.48)

Из рассмотрения выражения (1.48) можно сделать следующие выводы:

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

пи k сходится быстрее, чем сумма по k при п — 0, так как jx0*<С

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

в ней зависит только от числа Fo и заканчивается при Fo äs 0,8, а в неосесимметричной составляющей— от двух комплексов: Pd и Fo. А именно, как будет показано ниже, с увеличением аргумента

V P d , определяющего при всех

равных условиях скорость враще-

ния валка, модуль

Іп (VTPd р)

*

 

уже при р, близких к единице,

Іп [У і Pd)

стремится к нулю, а следовательно, с увеличением скорости враще­ ния время неустановившегося процесса в неосесимметричной состав­ ляющей уменьшается. Даже в предельном случае, когда P d=0 время неустановившегося процесса в неосесимметричной составляющей будет меньше, чем в осесимметричной, в связи с тем, что Ро/г<М-л£> поэтому экспонента в степени быстрее стремится к нулю, чем в степени р^.

3. Наличие в выражении (1.48) разности между средней темпе­

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

кой к заданной средней температуре их поверхности, т. е. Ѳ0

Ѳ.

22

Интегралы знаменателей С0/г и Cnk равны

I рJо (иоар) Ф = - у Л (роа) ;

J* Р^п лаР) dp = — Jn+i (Р'па)•

о

2

Если при данном распределении, температур на поверхности валка вращения не будет, т. е. ш0 = 0, а следовательно, Pd = О,

то (1.48) с учетом (1,36) примет вид

оо1

{[Ѳ— Ѳ0] /0(ИоАР) РФ

Ѳ(р, ф, Fo) = Ѳ

 

о

 

X

 

■2к=1 1

 

 

 

 

 

І р^ооаР) Ф

 

 

' со

Г"

со

1

 

 

 

 

 

J Рn+lJn (Pnk Р) Ф

 

Х - / 0 (ц0*Р)в ^ Fo+

 

 

іо

ow)^

X

п= 1

L

А=1

 

X j n (рЛАр)<Г^Ро] [N„cos(ncp) + Mnsia(ncp)J.

(1.49)

На квазистационарном режиме (т. е. при Fo

оо), так как ис­

чезнут решения Ад (р, Fo) и А*п (р, Fo), выражения (1.48) и (1.49) примут более простой вид

со

Ѳ(Р, ф) = Ѳ+

\У і Pd р)

i PdFo г н

/ \

I

п ѵ

ѵ’- е

 

[ІѴпсоз(лф) +

 

і Л Ѵ і Pd

)

 

 

 

 

 

л=1

 

 

 

 

 

 

+ M„sin (лф)]

 

 

(1.50)

и при Pd = О

 

 

 

 

 

 

0 (р, Ф) =

Ѳ+ 2 Рп [Na cos (Лф) +

М пsin (лф)],

(1.51)

 

П= 1

 

 

 

 

 

т. е. осесимметричная составляющая на квазистационарном режиме

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

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

23

(1.50), но при этом следует оценивать поведение модуля функции. Согласно (1.37)

/„ \ У і Pd р) _

ber„ (V'Pd р) +

i bei„ ( Y Pd p)

 

/„ ( Y i Pd )

ber„ (y^Pd) +

i bein ( Y Pd)

 

Если учесть, что | In ( Y i Pd p) | =

] /' ber^ [ Y Pd p) + bei* [ Y Pd p) =

= bn ( Y Pd p), то

 

 

 

 

lg [ V * Pd p)

b a i V Pdp)

(1.52)

/„(K iP d )

6„(KPd) '

 

Исследования зависимости ( .52) упростятся, если принять во внимание^ что bn (]/P d )< b 0 (j/P d), а при P d > 4 bn (У"РЗ)

b0(]/P d ). Тогда при анализе достаточно рассматривать пове­ дение модуля модифицированной функции Бесселя комплексного аргумента нулевого порядка, т. е. при Pd>>4 имеет место равенство

Іп (K fP dp) _

bn {YPA_p) _

fro (V^Pd_p)

(1.53)

/„ (K /P d )

bn [YPd)

é*(K P d)

 

При этих же значениях аргумента функции Кельвина можно считать по асимптотическим формулам [21, 40, 95], т. е.

b e r ( l / p H

) « ^

i c o s

Y Pd__ n

(1.54)

Y l

8

 

V2nVPâ

 

 

"Pd

eVvAh'2

i ' Y Pd

тс

(1.55)

r -

г

- S in

— -------------

 

Y 2 n V P d

 

[ Y 2

8

 

§ 3. Температурное поле валка с наружным и внутренним (комбинированным) охлаждением

В этом и в последующих параграфах будет рассмотрена эффек­ тивность применения наружного, внутреннего и комбинированного охлаждения валков и роликов УНРС. Решение уравнений темпе­ ратурного поля для валка с комбинированным охлаждением ана­ логично решению, приведенному выше. Отличие заключается в сле­ дующем. При решении уравнения (1.20) с граничными условиями (1.136) получим выражение

Ло(Рі. р) = 2

Ѳ—Ѳх ln p + Ѳ

(1.56)

 

ln 1/рх

 

Общим решением (1.22) для полого цилиндра будет

^ о ( Р і> Р> F ° ) — ' S C okV ok{ \ i о /,р )е —^0feFo

(1.57)

fc=i

24 t

Здесь Vok (р.олр) — комбинированная

функция

Бесселя

нулевого

порядка, определяемая уравнением

 

 

 

 

 

 

 

FoА(PoaP) — Jo(PoaP) Y о (Poa) — Y oa(PoaP) J о(Poa)>

 

 

где

Y 0 (p0A)

и

Y 0 (j.i0/;P) — Бесселевы

функции

нулевого

порядка

2-го

рода; |xoft — корень характеристического уравнения,

получен­

ный из граничных условий. Так как

по

условию

(рі, р =

1,

Fo)

= 0 и

Л;

(р = Pl, Fo) = 0, то

в

(1.57)

Vok ( p j

= 0

и

Foa(PoaPi) =

0.

Откуда

 

 

 

 

 

 

 

 

Jо (PoaPi) Y о (Роа) — J о(Poa) Y о (РоаРі) — 0-

 

(1.58)

Выражение (1.58) является характеристическим уравнением. Аналогично предыдущему получим выражение для Cok — коэффи­ циентов разложения функции А0 (рі, р, 0) в ряд по комбинирован­ ным функциям Бесселя, а именно с учетом (1.56) и (1.19)

 

Ѳ—0!

ln p + Ѳ

— Ѳ0 Voa(PoAP) РФ

 

 

-2.f

 

 

 

 

 

С,OA

О М- ln Upi

 

 

 

(1.59)

 

 

 

 

 

 

 

 

IPFoa (PoaP) dP

 

 

 

Результирующая зависимость для

коэффициента А 0 (р,

Fo) бу­

 

дет иметь вид

 

 

 

 

 

 

А о(рі, р, Fo) = A0(p) + i4o(p, Fo) = 2(jb=-^lnp +

0 —

 

 

 

 

ln l/Pi

 

 

 

-S

lnp + 0

— Ѳо VoA (PoaP) pdp

-^ oa1*

 

0

 

 

 

 

1

 

VOA(pOA p) e

 

 

 

 

 

 

 

 

 

JpFoA(^0AP)dP

 

 

 

 

 

 

 

 

(1.60)

 

Нахождение коэффициента A n (р,- Fo) аналогично случаю сплош­

 

ного валка, т. е.

Ап (pi, р,

Fo) = А п (рі, р) + А п (pi,

р,

Fo). Об­

 

щим решением уравнения (1.31) является выражение

 

 

 

а л р і.

p ) = c 1i n ( V m P) +

c 2K n ( V i P d p ) ,

 

(1-61)

,

где Кп [ V t'Pdp) — модифицированная функция Бесселя (функция Макдональда) комплексного аргумента n-го порядка.

Произвольные постоянные найдутся из граничных условий

(1.156), а именно

А п (Рі, р = 1) = i V / pdFo= C1/ n { Ѵ Щ + С2К П(K iP d );

An(P= Pi) = 0 = C J n (jA P d pi) + C2K n ( ] / i Pd Pi);

25

откуда

 

 

 

 

 

 

 

 

 

 

 

Ci = _____________ NneiVÜT°_Kn { V »Pd Pl)

 

 

.

 

- I n i V T P d Pl) Kn [ViPd) + In {ViPd) KniViPdpx)

^

 

- Nneipdp0In { V T PdPl)

 

 

 

 

2 '

/ л [ V i Pd) Kn [ V i Pd Pl) -

In [ V t Pd p j Kn [ V i Pd)

 

Подставляя Cx и C2 в (1.61), получим

 

 

 

 

 

 

 

 

 

А Л р і. p) =

 

 

 

 

 

 

 

 

Kn { У ’■Pd Pi)

/ J K t' P d p ) _

/,t (l/ 'tPdp1)

Kn(lA-Pdp)

pj g£ PdFo

Kn ('KtPd)

/ „ ( і Л -Pd)

i n {V7pâ)

Kn ( K t P d ) .

 

 

jCB ( V 7 P d Pl)

 

/ „ ( K t P d p J

 

 

 

 

 

 

Kn{ VTpd)

 

i n {ViPâ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.62)

Обозначим

в (1.62) выражения ^ ( V ^ P d f j j __т ^ и !п (V^t'Pd pj

_

 

 

 

 

ff„ (K /P d )

 

/„ (K iP d )

 

= пи.

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_

/„(KiPdp)

 

... ^„(KtPdp) .

 

 

 

 

 

iil\

у—г-

—~Г~

 

ГЛо — ' v

——

 

 

А п (рх, p) = J V /PdFo - 2*SYj

F*\______ U K y P d )

 

(1.63)

 

 

 

 

 

m-i то

 

 

 

 

Если вращение отсутствует, т. е. ю0

— 0 и Pel — 0,

то

общее

решение уравнения (1.35)

для полого цилиндра будет

иметь

вид

 

 

А п (Рі>

Р) — £іРп + СгР

п>

 

 

(1-64)

где Сх и С2

найдутся из граничных условий, а именно А п (Рі,

р =

= 1) = Nп =

-j- С2; А а (р == рх) = 0 =

0 1р”-(-02р1

”.

Следо­

вательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

Ап (Рі> рЧ_

(Р/Рі)" — (Рі/Р)

” ,,

 

 

(1.65)

 

 

р)

 

 

 

 

 

 

 

 

 

Так как зависимость (1.63) получена в комплексной форме, то выделим в ней действительную и мнимую части.

Из предыдущего параграфа нам известно представление моди­ фицированных функций Бесселя 1-го рода n-го порядка, комплекс­ ного аргумента в виде функций Кельвина. Аналогичные выраже­ ния [21, 36, 95] для функций Макдональда имеют вид

К п [ V i Pd) = ker„ [ у Pd) + i kei„ ( V pd) ■

(1.66)

26

Чтобы сократить форму записи, введем обозначения

 

 

а — кег„ (]/"Pd рх);

b =

kei„ (]/P d рх);

с = кег„ ( ] /

Pd);

 

d = k ei„(]/P d );

m = ber„ (]/P d рг);

 

n bei„ (]/P d pj);

 

u — ber„ (j/^Pd);

v ~ bein (j/^Pd];

 

x ber„ (j/^Pd p);

 

y = be\n (]/P d p );

z =

ker„ (|/P d p );

 

щ =

kei„ (]/P d p);

 

и далее

 

 

 

 

 

 

 

 

 

 

 

 

 

(ax+ nwbytnz) —Ф;

(bx + aynztnw) = F\

 

 

(<0.11 + ndbvmc) = D\

(bu-\-va mdnc) = K-

 

 

Тогда (1.63) примет вид

 

 

 

 

 

 

 

 

А п(pi,

р) =

Nn (0D +

FK)A~1

~ ° K) X [cos (PdFo) + i sin (PdFo)],

 

 

 

 

К 2 + D2

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

An (pi,

P) = N n |R e^„(pi, p) +

i l m A n (p1, p)),

(1.67)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

И е Л Л Р ь

P) = N j

-<*£ +

FK) cos (P d F o )

-

(F D -

■* «> s i!L g d F a . I

;

n v

n \

 

 

K2Jr D 2

 

J

 

 

 

 

 

 

 

 

 

 

 

 

 

(I.68a)

Іш Л

(p

p) =

N

f (Ф° + FI<) sin (PdFo) +

(FD — ФАЗ cos (PdFo))

 

 

n ^ 1’

n \

 

 

K2+ Ö2

 

* }T

'

 

 

 

 

 

 

 

 

 

 

 

 

( 1. 686)

Модуль (1.67) равен

 

 

 

 

 

 

 

 

 

 

 

 

I A n (Pi. P) I =

V Re2 A

(Pi. P) +

 

Im2 А

(Рь P) •

(L69)

Далее найдем решение A (P n p, Fo). Общим решением уравнения (1.43) для полого цилиндра будет выражение вида

Ап (рі>

р> Fo)

Сп^Уп (pnft р) е-^feFo

(1.70)

 

k = \

 

 

 

Аналогично предыдущему

Vn (p,nÈp) есть комбинированная

функция Бесселя n-го порядка, т. е.

 

 

Уп (М’лйР)

Jп (pnfcP) Уп (рлй)

Уп лАР) /„ (pnft) >

 

где |хпА — корень характеристического

уравнения

 

 

Уп лаРі) = 0.

 

(1.71)

27

Соответственно с учетом (1.63) получим зависимость для Cnk\

 

 

 

ті

и (У~iі Pdp)

■m.

/Cn(]/ 7Pdp)

 

 

 

Мпёі PdFo

 

i)

K n {V <Pd)

 

 

 

 

 

 

PKn(PnftP) dP

 

 

 

 

 

тг m2

 

 

 

 

-'nk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 'P ^ ( ^ P ) dP

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.72)

Окончательное

выражение для

коэффициента Ап (рх,

р, Fo)

с учетом (1.29), (1.63),

(1.70)

и (1.72) будет

 

 

 

 

 

 

 

L /п (K lPdp) _

Kn (l^Tpd р)

Л„ (рі, р, Fo) =

N

j

MF"

‘ ' ^

P d ) -------- ' к Л У Ы

_

 

 

 

 

 

 

 

rnt т 2

 

 

 

 

со

 

/„(K t Pdp)

Kn (VTpdp)

 

 

 

 

 

1

/ n(K7pd)

' КА У Т Ъ )

pVn(PnkP)dp

 

 

 

 

 

/ті! — m2

 

 

 

 

 

 

 

 

 

X

Е

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

*=і

 

 

 

JpK (P„fcP) dP

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.73)

Зависимость для коэффициента Вп (рх, р, Fo) будет аналогична (1.73) , только вместо Nn будет М п. Если в (1.6) подставить (1.60) (1.73) и выражение для Bn (р, Fo), то получим решение уравнения теплопроводности для цилиндра с комбинированным охлаждением

 

Ѳ(р,

ф,

Fo) =

^

i i n p

+ 0-

 

 

1

Ѳ— Ѳг lnp +

 

 

 

СО

I

0 - 0 0

Ko (pofeP) pdp

 

ln 1/pj

 

 

 

 

■S

 

 

 

 

 

X

 

 

 

 

 

 

 

fc=i

 

 

J PKo(PoftP) dP

 

 

 

 

 

 

 

 

 

 

 

CO

 

___

___

 

 

 

L

 

^ (K ip d p )

K n i V i Pdp)

X VOA(po,P) e ^ ‘Fo +

7

j [

1 , л ѵ ш

)

* К П( Ѵ m )

 

 

 

 

 

 

in* m9

28

I n { V i P d p )

Kn ( K t P d p )

 

1/„ (ѴШ)

2 (vra)

PVn (PnftP? dp

m1 — m2

 

 

 

X

 

J P ^ O W ) dP

 

 

0

 

X Vn (М-плР) e ^ 0] e PdFo [Nacos (ncp) +

sin (ncp)]. (1.74)

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

1. Осесимметричная составляющая представляет собой неуста­ новившееся температурное поле с температурами: на внешней по­ верхности Ѳ и на внутренней поверхности

2.Время неустановившегося процесса определяется также осе­ симметричным полем, но в этом случае оно зависит от рх = R J R q. Например, при рх = 0,2 Fo « 0,4, т. е. практически в два раза меньше времени неустановившегося процесса в сплошном цилиндре.

3.Неустановившийся процесс в неосесимметричной составляю­ щей поля также зависит от чисел Pd и Fo, но критерий Pd является

аргументом

комбинации функций

Іп { у i Pd) и К п І У i Pd). Интег­

ралы знаменателей Cok и Cnk в этом случае равны:

 

 

1

 

2 [dp (pQfr)

( W

l) ] .

 

 

J РѴІ (цоаР) Ф

 

 

•n'2tx0&^0 (tl0fe Pi)

 

 

 

о

 

 

 

 

IpV* (Fnfep) Ф

2 [Jl (Pnk)~j2n OW l)] _

 

 

 

nYnkJn (^nfePl)

 

 

Если Pd =

0, T. e. вращение

отсутствует,

то (1.74)

с учетом

(1.64) примет вид

 

 

 

 

 

 

 

Ѳ(р, ф, Fo) = y~ - p ln p + Ѳ—

 

 

 

 

 

ln I/Pi

 

 

 

03

1

9 — 0!

 

Fo(PofeP) Pdp

 

 

 

 

ln p -j- Ѳ — 90

 

 

 

 

ln 1/Pi

 

 

 

 

^F0+

 

 

 

 

 

Vo(Foftp)e

&=і

 

Jp^o (Po*P)dP

 

 

 

 

0

 

 

 

 

 

29

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