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

книги из ГПНТБ / Совершенствование теплового процесса листовой прокатки

..pdf
Скачиваний:
16
Добавлен:
22.10.2023
Размер:
14.01 Mб
Скачать

по

При

цене

1

опорного валка Цайв

 

=

15 ООО руб. величина

Э3

формуле (245) составит:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Э3

=

15 000-1,76 =

24 600

руб.

 

 

 

 

 

Экономическую

эффективность

от

уменьшения

потерь

 

металла

при

 

сокращении

количества

 

обрывов полосы

можно

рассчитать

по

формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

4

=

 

 

 

 

 

 

 

'

 

 

 

 

( 2 4 ? )

где

II1

и Цт

— цены

годного

и

отправляемого

в шихту

металла,

 

 

d±

н d2

 

руб/т;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— количество порванных

полос до и после

внедрения

 

 

 

 

 

данного

мероприятия,

%;

 

 

 

 

 

 

 

 

 

 

L n o p

средняя

величина потерь

металла

при

одном

по­

 

 

 

 

 

рыве

полосы,

т;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qx — масса одного рулона, т.

 

 

стойкости

валков

 

Эффект от уменьшения

простоев

и увеличения

при

сокращении

количества

обрывов можно

подсчитать

по формулам

(241),

(243—245).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Экономическую

эффективность

от увеличения

 

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

стана

при освоении

прокатки

более

широких

полос

можно

опреде­

лить

по

формуле (241)

или

(243), в

которой

процент

увеличения

годовой производительности следует вместо формулы (242) вычислить'

по

формуле

 

 

 

 

 

 

а= А*-А>

-100%,

(248)

где

А г и А 2 — годовая

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

до и после

вне­

 

дрения

данного

мероприятия.

 

 

 

Приведенные формулы и примеры расчетов

показывают,

что

совершенствование теплового процесса листовой прокатки способ­ ствует улучшению качественных показателей работы листопрокатных цехов и снижению себестоимости продукции.

П Р И Л О Ж Е Н И Е I

ИСХОДНЫЕ Д А Н Н Ы Е ДЛ Я РАСЧЕТА ЭНЕРГЕТИЧЕСКОГО Э К В И В А Л Е Н Т А СОПРОТИВЛЕНИЯ МЕТАЛЛА ДЕФОРМАЦИИ [13]

 

X

*

 

 

 

X

 

 

 

 

 

 

 

 

 

Марка сплава

о

 

 

Марка

сплава

о

 

 

в

 

 

в

 

 

 

, 6

 

с

 

 

оо

 

с

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

12,0

2,92

0,63

Х18Н9Т (ЭН1Т)

41,0

2,40

0,91

железо

 

 

 

Х22Н5АГ9 (ЭП20)

55,0

5,10

0,72

Хромансиль

14,5

2,49

0,58

Динамная

сталь

37,5

3,40

0,63

Низкоуглероди­

25,0

0,91

0,96

Трансформаторная

50,0

4,80

0,62

стая сталь (особо

 

 

 

сталь

 

 

 

 

мягкая)

 

 

 

ЭЗА

 

40,0

11,2

0,38

Малоуглеродистая

29,0

1,76

0,74

ЭИ496

 

32,5

7,20

0,45

сталь

 

 

 

Х23Н12 (ЭИ417)

36,0

17,0

0,46

Ст.О

25,0

5,62

0,46

ЭИ409

 

48,0

16,0

0,47

Ст.1

26,0

1,33

0,73

ЭИ602

 

54,0

7,30

0,65

Ст.2

30,0

3,06

0,62

ЭИ659

 

70,0

0,38

1,10

08пс

30,0

7,70

0,48

12Х5МА

 

48,0

2,30

0,72

08кп

23,0

3,46

0,60

ЭИ712

 

37,5

3,60

0,62

10сп

30,0

2,95

0,64

ЭИ763

 

47,0

1,40

0,80

20сп

37,5

3,16

0,64

ЭИ811

 

67,0

2,70

0,68

20А

35,0

6,45

0,50

Э91

 

38,0

2,40

1,0

40

35,0

8,36

0,48

Р9

 

27,0

7,90

0,61

45

35,0

8,66

0,48

ВЖ98

 

58,0

8,90

0,62

50

40,0

10,0

0,47

Ж В П

 

38,0

1,10

0,95

85

50,0

14,7

0,43

50ХФА

 

45,0

3,2

0,66

09Г2

32,0

5,90

0,46

ЭИ846

 

33,0

5,40

0,74

10Г2

35,0

4,40

0,59

ЭИ852

 

33,0

19,4

0,26

12Г2А

50,0

10,0

0,34

ЭИ962

 

60,0

4,0

0,64

65Г

40,0

17,6

0,35

ЭИ435

 

50,0

7,10

0,58

У8А

39,0

1,80

0,84

 

 

 

 

 

У9А

30,0

12,7

0,41

Титановые

сплавы

 

 

 

У10

45,0

2,5

0,79

 

 

 

 

 

У10А

40,0

6,40

0,60

Тех ническийтитан

50,0

6,70

0,36

У12

55,0

2,70

0,76

ВТ-1

 

42,0

5,20

0,48

У12А

62,0

0,76

1,025

ВТ-2

 

44,0

4,20

0,54

У8ГА

40,0

11,4

0,44

ВТ-6

 

00,0

1,10

0,76

85ХФ

35,0

9,0

0,49

ОТ4

 

64,0

2,40

0,68

25ХГСА

38,0

5,70

0,57

ИМП-7

 

83,5

1,20

0,86

30ХГСА

47,5

8,60

0,45

ИМП-1А

 

62,5

4,0

0,61

Х05

50,0

0,75

0,99

Алюминиевые

 

 

 

Х18Н9 (ЭЯ1)

25,0

1,90

1,0

 

 

 

2X18Н9 (ЭЯ2)

60,0

3,80

0,70

сплавы

 

 

 

 

Х18Н25С2 (ЭЯЗС)

60,0

3,60

0,76

 

 

 

 

 

Х13Н4Г9 (ЭИ100)

34,0

1,80

0,84

Алюминий (99,5%)

6,00

0,64

0,62

ЭИ401

35,0

7,20

0,69

Алюминий

 

1,80

0,28

0,74

0Х20Н4АГ10

50,0

2,5

0,92

(99,99%)

 

 

 

0,71

(ННЗ)

 

 

 

АМц

 

5,0

0,60

Х14П4НЗ Т

34,0

3,30

0,78

АМг

 

10,0

1,30

0,59

(ЭИ711)

 

 

 

АСм

 

4,0

1,60

0,45

Х14Г14Н (ЭП212)

30,0

9,30

0,62

АМЗ

 

7,50

6,40

0,30

Х18АГ14 (ЭП213)

45,0

3,60

0,86

АМг58

 

16,0

4,30

0,42

0Х17Н5Г9АБ

60,0

5,90

0,69

Д1

 

8,80

3,50

0,41

(ЭП55)

 

 

 

Дюралюминий М-2

10,0

3,20

0,39

* 1 кГ/мм 2 равен

9,8 Мн/м2 .

 

 

 

 

 

 

285

П Р И Л О Ж Е Н И Е II

МЕТОД ЭЛЕКТРОТЕПЛОВОЙ АНАЛОГИИ

В гл. I I и V I I были представлены результаты исследования температурных полей валков, полученные методом электротепловой аналогии.

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

Теоретические основы метода электротепловой аналогии достаточно полно разработаны и широко освещены в специальной литературе (59; 72—80].

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

ных в

гл. I I и V I I .

Общие

положения

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

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

Моделирование стационарных тепловых процессов

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

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

Применительно к стационарному процессу в электрическом поле в двумерном пространстве это уравнение будет иметь вид:

^

+ _ L | L +

! ! £ = <),

(250)

дг2

г

дг

дг2

 

где t — температура; и — потенциал.

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

Переход к безразмерной форме уравнений Лапласа (249) и (250) и краевых усло­ вий производится посредством метода собственных масштабов.

С этой целью геометрические размеры изучаемого тела А выражают в виде отно­ шения линейных размеров тела к условно принятому за единицу измерения раз­ меру R, а для моделирующего тела В аналогичное отношение сходственных разме­ ров выражают через размер Ri. Таким образом, получаем следующие безразмерные линейные величины для тел А и В:

= p = i d e m ;

x = ^ = z = i d e m -

286

Д л я безразмерных температур и потенциалов получаем следующие отношения:

 

 

 

 

t—t,

а—

и2

 

 

 

 

 

 

Wi — и.

 

в

которых

за

нулевой уровень отсчета

приняты температура h и потенциал

и%,

а

в

качестве

собственных масштабов

для температур — максимальная

раз­

ность

t\ —

h

и для потенциалов — соответственно

разность «i — иг.

 

Рис. I . Принципиальная схема интегратора ЭГДА9/60 с моделью валка из электро­ проводной бумаги или Л-сетки:

ПУ — блок питающего устройства; ИУ — блок

измерительного устройства;

Р — реохорд; Г — гальванометр-индикатор; Д

— декада сопротивлений

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

Граничные условия этого уравнения можно пояснить на примере рис. I . На по­ верхности Ft 6 = 100%, на поверхности F2 6 = 0%.

287

Аналогично для явления электропроводности в модели В имеем:

 

dp2

_ L

Ё1

дг2

= 0.

 

(252)

 

р

др

 

Граничные

условия этого

уравнения

на

соответствующей

поверхности Fi

и = 100%, на

поверхности F 2

и = 0%.

 

 

 

полную

Сравнивая обе системы уравнений для детали и модели мы констатируем

аналогичность основных уравнений и одинаковые числовые

значения

краевых

условий.

 

 

 

 

 

 

 

 

 

6Ср'0[Чср"0

 

 

 

 

 

 

V 7

-4-

4—

-J

 

 

1пов

z

Г "

acp-i "ср-1

Рис. П. Двумерная электрическая сетка

Отсюда, как следствие, появляется возможность выразить интересующую нас функциональную зависимость в форме (6 = и = / (р, г), что утверждает формально одинаковое решение процессов различного физического содержания.

Для проведения исследований двумерного стационарного осевого температур­ ного поля валка был использован интегратор ЭГДА9/60 (электрогидродинамических аналогий) и интегратор из сетки сопротивлений.

Интегратор ЭГДА 9/60 (см. рис. I) использовали в том случае, если модель валка была выполнена из электропроводной бумаги, а сетчатый — для модели, представ­ ленной в виде дискретной сетки сопротивления.

Принципиально оба интегратора ничем не различаются.

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

ний, подключается

к интегратору по схеме, показанной

на рис.

I . В частности,

на нем схематично

изображена модель четверти рабочего

валка

листопрокатного

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

Моделирование стационарного температурного поля в осевом сечении валка на электропроводной бумаге заключается в следующем. Модель валка в требуемом масштабе вырезают из электропроводной бумаги с припуском по 3—5 мм на стороны

288

для подключения шин с соответствующими граничными условиями первого рода. На поверхностях, где заданы граничные условия третьего рода (коэффициент тепло­ обмена а ) , при вырезании модели добавляют дополнительный слой бумаги шириной 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

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