Файл: Совершенствование теплового процесса листовой прокатки..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 78

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

 

2Ra

Ш

R

(253)

откуда

 

 

 

 

Д л я того чтобы

направление тока в слое h было по

нормали к поверхности,

необходимо разрезать

этот слой на полоски

и переложить

их электроизолирующей

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

Так как электрическое сопротивление на модели изменяется по линейному закону, а в цилиндрическом теле равенству тепловых сопротивлений соответствует равенство:

R

1 п - ^ = 1 п - ^ - =

=

 

'3

то полученное на модели температурное поле необходимо перестроить, т. е. разбить модель по радиусу на участки r(- соответствующие одинаковым теп­ ловым сопротивлениям.

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

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

Рассмотрим случай, когда узлы сопротивлений сосредоточены в центре объемов, как указано на рис. I I . Уравнение (251) для узла с центром 6 i будет иметь в этом случае вид:

1

 

 

 

 

 

Ар L Ар

Ар

Pi

L

2 Ар

2 Др

 

 

9 1

6 2

0.

AZ

 

AZ

AZ

 

 

 

 

Если привести подобные члены, умножить все члены на Др и преобразовать это выражение, то получим:

e 6 - e t

+

е , - в 4

 

з 2 - е х

2 Др

2 Др

AZ2

AZ2

2рг + Др

 

г — Др

Pi Др

Pi Др

19 А. В. Третьяков

 

 

 

289


 

2Др

 

 

2Др

 

 

 

 

 

 

 

 

 

 

 

выражения ^

и

^

^

можно

приближенно

представить

логарифми­

ческими

функциями:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l n

Pi

 

2 Д Р

 

Q

, pt- + Др

 

2 Др .

 

 

р,- — Др

2pi

Др'

 

р{

 

 

2pi

+

Др '

 

тогда окончательно

данное

выражение

будет иметь

вид:

 

 

 

 

 

е 5 — e

i , e i — е 4 , 9 з — 9 1 , е 2 — 9

i _ п

 

 

 

In

Pi

 

in Pi-i

 

Р/Др

 

Р/Др

 

 

 

 

 

 

tl±L.

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

теперь >уравнение

(252)

подобным

образом

привести

к

виду

(254), то по-

лучим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

» 5 ~ » 1

. »! — и 4

 

» 3

— % ,

ы2

« 1

 

 

 

 

* рi+l. 1

 

1\Г\1 +

J

^

i + ! f

L

_

J i

L

=

o.

(255)

 

 

*pl

 

*z p l

 

 

/?zp <

 

 

 

 

Выражение (255) есть не что иное, как закон Кирхгоффа для узла эквивалент­ ной электрической цепи, показанной на рис. I I , где 9i соответствует и\ и т. д.

Чтобы относительные потенциалы эквивалентной электрической цепи соот­ ветствовали относительным температурам в сходственных точках, необходимо выпол­ нение следующих равенств:

 

 

 

 

щ 9J±L R n

=

R p { + 1 ,

m A

-

RN =

R P I

^

R

N

= R z p { ,

 

(256)

где

/?л? — номинальное

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

интегратора.

 

 

 

 

 

Если

на

каких-либо

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

заданы

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

первого

рода

( 9 п о в ~

относительная

температура

поверхности),

то

к

соответствующим

узлам

модели

подводится

относительный

потенциал

и П 0 В )

 

а если заданы граничные усло­

вия

третьего

рода,

например, в

радиальном

направлении

 

 

 

 

 

 

 

 

 

 

 

 

 

дд

_

а

(9 среды 9 п о в )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dp

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

или

по

оси Z

 

 

 

 

 

ае =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Череды — 9пов).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dZ

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

%

 

 

 

 

 

 

 

 

 

 

 

то,

преобразовав

эти

выражения

к

виду

 

 

 

 

 

 

 

 

 

 

 

 

 

9 п о в

 

9 л

 

 

Q

 

 

 

о

 

бпов

 

Од

_

а

 

а

 

 

 

 

а

 

 

р 0

 

" ~ "

с р е

д ы

— Опов,

а

 

 

Г"^52

— "средыд ы

«п

 

 

 

 

 

 

 

 

 

~ ° '

 

AZ

~

 

 

 

 

 

 

X

"

 

 

Др

 

 

 

 

 

 

к

 

4 р , Д р

 

 

 

 

 

 

 

 

 

Р о +

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

формулы

для

расчета

электрических

сопротивлений,

соответствующих

тепловым сопротивлениям:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

«

i „

 

Р

о

 

п...

о

_ _

?

 

 

A

z \

 

 

 

 

 

 

 

ч

=

х

1

п

P -^V

 

 

=

 

-

- TPT V RN-

(257)

 

Подключение модели к питающему и измерительному устройствам

интегратора

и методика

определения

потенциальных — температурных

полей

аналогичны

опи­

санному

выше (см. рис. I) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исследование нестационарного радиального температурного поля в активной

зоне

валка

проводили

на интеграторе

с

J?-сеткой

 

методом

Либмана

[79, 80].

 

290


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

В частности, уравнение нестационарной теплопроводности в безразмерном виде для полуограниченного тела имеет вид:

а ' дх ~~ дХ2 '

( 2 5 8 )

где X— безразмерная координата.

Рис. I I I . Электрическая модель для исследования неста­ ционарного температурного поля полуограниченного стержня

Для t-той узловой точки в и-ный период времени это уравнение в конечных разностях будет иметь вид:

 

Ji

( n - i ) — %

 

ДХ

 

 

L\X

 

 

 

 

а Дт

 

 

 

 

ДХ

 

 

 

 

 

 

 

 

 

„ — в ,

 

e ( t - i ) п — е(-

>

 

(259)

 

 

 

 

 

ДХ2

 

 

 

д у 2

"

 

 

 

 

 

 

 

г

 

ДА2

 

 

 

 

или

 

0 г - ( п - 1 ) - 9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

«ДТ

~ =

( 9 ( ' - 1 >

» ~

+

( в «+1>

» -

9 < J -

 

 

 

ДХ2

 

 

 

 

 

 

 

 

 

 

Уравнение Кирхгоффа для этой же точки эквивалентной электрической цепи

имеет

вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

us

« о

» 2

— " о

, Щ — " о

 

 

(260)

 

 

 

 

R%

 

Ri

 

 

 

 

 

 

где

Rx — электрическое

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

моделирующее

время

процесса;

Ri

и # 2 — электрические сопротивления, моделирующие тепловые сопротивле­

 

 

ния

(см. рис. I I I ) .

 

 

 

 

 

 

 

 

Уравнения

(259)

и

(260) будут

идентичными,

если

выполняется

условие

 

 

 

 

 

1 -

1

^ т

— р . ,

 

 

 

 

 

 

 

 

 

 

-abx-«N-

 

 

 

 

ДХ2

Это значит, что параметры /?-сетки можно рассчитывать из соотношений:

19*

291


Формулы (261) были использованы для расчета /?-сетки при исследованиях нестационарного температурного поля активной зоны валка. Расчет сетки сопротив­ лений для моделирования температурных полей в более сложных системах приведен в работе [59].

Нестационарное радиальное температурное поле валка в его основной зоне может быть получено на электрической модели из R— С-сетки [72]. Этот метод состоит в использовании аналогии дифференциального уравнения нестационарной теплопроводности и электрического поля, а именно:

^ а д д ^

j _ я .

( 2 6 2 )

где

г — текущий

радиус

валка;

 

R

— максимальный радиус валка;

гэ

— текущее

электрическое

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

R3

— максимальное электрическое сопротивление системы;

RQ3

— удельное

электрическое

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

с о э

— удельная

электрическая

емкость;

с о т = с рР — удельная

объемная тепловая емкость;

т э

— время

электрического процесса;

т х

— время

теплового

процесса;

К — теплопроводность

или

— R0T — тепловое сопротивление;

 

 

 

 

 

А

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

и— безразмерный потенциал; 6 — безразмерная температура.

При R — С-сетке применяется аналогия, согласно которой температуре, как основной величине системы, поставлено в соответствие электрическое напряжение, тепловому сопротивлению — электрическое сопротивление и тепловой емкости — электрическая емкость. Если обозначить

 

&р =

О

Т

 

сек/ч;

kR

=

Ух*

ом-м/[ккал/(м-ч-град)];

 

— ;

kx = —,

 

А;о х

 

 

Рт

T x

 

 

К

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

& с =

ф'град/ккал;

k, = —т

,

 

 

 

 

 

'-от

 

 

 

 

 

 

 

 

и

подставить

их в

уравнения

(262) и (263),

то

указанные уравнения

идентичны,

а

описываемые ими

явления

подобны

при

 

выполнении

условия

 

 

 

 

 

 

ф 2 * ^ в =

1,

 

 

(264)

 

 

 

 

 

 

"•х

 

 

 

 

 

 

Расчет параметров R — С-сетки

производят

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

 

 

Аналогично решению задачи стационарного

температурного поля

на 7?-сетке

правая часть уравнения (262) представляется системой конечно-разностных уравне­ ний, а производная температуры по времени (или ди/дх) моделируется на R — С-сетке непрерывно и не представляется в конечных разностях.

292

I


Сопротивления и емкости для электрической сетки можно найти по формулам:

 

 

 

 

Ро

Д Р Л

ом;'

R2 z

=

Р 0

A?2kR

 

 

 

 

 

7?(2p0

+ A p i H

(2р0 -

Др2 ) X

 

и т. д.

 

 

 

 

 

 

 

 

 

 

где

V — объем элемента тела.

 

 

 

 

 

 

 

Радиусы

центров

одинаковых объемов р п

можно определить по

формуле

 

 

 

 

Р Ч

= ) / Р О

+

 

О-Ро) .

 

 

где

л = 1,

2,

3, 4, . . ., k;

 

 

 

 

 

 

 

 

k — число элементов, на которые разбито моделируемое тело.

 

 

Применение R — С-сетки для расчета более

сложных систем

описано в ра­

боте

[72].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П Р И Л О Ж Е Н И Е I I I

РЕШЕНИЕ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ

(6)

 

 

 

Если на внешней поверхности валка задано граничное условие первого рода (7),

то с учетом

вращения

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

представить в

виде;

 

 

 

 

 

 

F(ty)

exp(iPdFo),

 

(265)

где

i -

- мнимая

единица;

 

 

 

 

 

 

Pd

=

-критерий

Предводителева;

 

 

 

 

 

 

 

со — угловая

скорость вращения

валка.

 

 

 

 

Уравнение (6) при краевых условиях (265) и (8) интегрируется приемом, пред­

ложенным Г., А. Гринбергом

[58]: искомое решение представляется

в виде триго­

нометрического ряда по полярному углу tf> с коэффициентами, являющимися функ­ циями полярного радиуса р и времени Fo.

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

аппарат бесселевых функций

нулевого

0 ),

я-ного

(со„) порядков

и их

модифика­

цию (/„).

 

 

 

 

 

 

 

 

 

 

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

(265),

представленное рядом

Фурье,

будет

иметь

вид:

F(\]3)exp (iPdFo)

=

-^а0

+ X

й п c o s

("^

+ 6 " s i n ("^

е х р ( i P d F o )

=

 

 

L

п=\

 

 

 

 

л

 

 

 

 

со

 

 

 

 

 

 

 

(266)

9 +

S Nnco%(n^)+Mn

sin

(т|э)

ехр

(iPdFo),

 

293