по |
При |
цене |
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 . |
|
|
|
|
|
|
П Р И Л О Ж Е Н И Е 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 - |
Д л я безразмерных температур и потенциалов получаем следующие отношения:
|
|
|
|
t—t, |
а— |
и2 |
|
|
|
|
|
|
Wi — и. |
|
в |
которых |
за |
нулевой уровень отсчета |
приняты температура h и потенциал |
и%, |
а |
в |
качестве |
собственных масштабов |
для температур — максимальная |
раз |
ность |
t\ — |
h |
и для потенциалов — соответственно |
разность «i — иг. |
|
Рис. I . Принципиальная схема интегратора ЭГДА9/60 с моделью валка из электро проводной бумаги или Л-сетки:
ПУ — блок питающего устройства; ИУ — блок |
измерительного устройства; |
Р — реохорд; Г — гальванометр-индикатор; Д |
— декада сопротивлений |
Итак, математическое описание явления теплопроводности в безразмерном виде при граничных условиях первого рода для тела А может быть представлено следую щим уравнением:
Граничные условия этого уравнения можно пояснить на примере рис. I . На по верхности Ft 6 = 100%, на поверхности F2 6 = 0%.
Аналогично для явления электропроводности в модели В имеем:
|
dp2 |
_ L |
Ё1 |
дг2 |
= 0. |
|
(252) |
|
р |
др |
|
Граничные |
условия этого |
уравнения |
на |
соответствующей |
поверхности Fi |
и = 100%, на |
поверхности F 2 |
и = 0%. |
|
|
|
полную |
Сравнивая обе системы уравнений для детали и модели мы констатируем |
аналогичность основных уравнений и одинаковые числовые |
значения |
краевых |
условий. |
|
|
|
|
|
|
|
|
|
6Ср'0[Чср"0 |
|
|
|
|
z
Г "
acp-i "ср-1
Рис. П. Двумерная электрическая сетка
Отсюда, как следствие, появляется возможность выразить интересующую нас функциональную зависимость в форме (6 = и = / (р, г), что утверждает формально одинаковое решение процессов различного физического содержания.
Для проведения исследований двумерного стационарного осевого температур ного поля валка был использован интегратор ЭГДА9/60 (электрогидродинамических аналогий) и интегратор из сетки сопротивлений.
Интегратор ЭГДА 9/60 (см. рис. I) использовали в том случае, если модель валка была выполнена из электропроводной бумаги, а сетчатый — для модели, представ ленной в виде дискретной сетки сопротивления.
Принципиально оба интегратора ничем не различаются.
Модель валка, выполненная из электропроводной бумаги или сетки сопротивле
ний, подключается |
к интегратору по схеме, показанной |
на рис. |
I . В частности, |
на нем схематично |
изображена модель четверти рабочего |
валка |
листопрокатного |
стана с внутренним охлаждением. Для наглядности модели, выполненные из элек тропроводной бумаги и сетки сопротивлений, совмещены.
Моделирование стационарного температурного поля в осевом сечении валка на электропроводной бумаге заключается в следующем. Модель валка в требуемом масштабе вырезают из электропроводной бумаги с припуском по 3—5 мм на стороны
для подключения шин с соответствующими граничными условиями первого рода. На поверхностях, где заданы граничные условия третьего рода (коэффициент тепло обмена а ) , при вырезании модели добавляют дополнительный слой бумаги шириной 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рг + Др |
|
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]. |
|
Этот метод заключается в использовании аналогии уравнения теплопровод ности, записанного в конечных разностях для произвольной узловой точки тела (рис. 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
Это значит, что параметры /?-сетки можно рассчитывать из соотношений:
Формулы (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), |
|