Рефетека.ру / Промышленность и пр-во

Лабораторная работа: Охлаждение изолированного провода

Кафедра КТЭИ


Переработка полимеров


Лабораторная работа


"Охлаждение изолированного провода"


Специальность – электроизоляционная, конденсаторная и кабельная техника


2009

Цель лабораторной работы – изучение процесса охлаждения изолированной жилы. Задачей лабораторной работы является исследование влияния параметров технологического режима охлаждения изолированной жилы на процесс охлаждения с применением метода математического моделирования и численных методов.

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


Vo, м/мин Ттпж, ˚С Тиз, ˚С Rтпж, мм Rиз, мм Т1, ˚С Т2, ˚С Т3, ˚С
45 20 250 1,825 3,3 80 60 20

Λи, Λж, Ρи, кг/м3 Ρж, кг/м3 Си, Дж/кгЧ°С Сж, Дж/кгЧ°С
370 0,18 770 9210 1940 305

 1, Охлаждение изолированного провода

 2, Охлаждение изолированного провода

3, Охлаждение изолированного провода

1000 1500 2000

Теоретические сведения


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

Модель №1

Геометрия изолированной жилы представлена на рис. 1.


Охлаждение изолированного провода

Рис. 1. Геометрия изолированной жилы


Допущения данной модели:

процесс стационарный;

теплофизические параметры постоянны;

диффузией тепла в направлении оси z можно пренебречь;

температура жилы не изменяется.

Уравнение энергии будет иметь вид:


Охлаждение изолированного провода (1)


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


Охлаждение изолированного провода

где Тпр – температура жилы провода; Тср – температура воды в охлаждающей ванне.

Начальные условия:


Охлаждение изолированного проводаОхлаждение изолированного провода


где Твых – температура изоляции на выходе из кабельной головки.

Для решения данной задачи воспользуемся методом конечных разностей. Выберем равномерную сетку по z и r:


Охлаждение изолированного провода


Аппроксимируем производные и получим:


Охлаждение изолированного провода (2)


С учётом выражений (2) выражение (1), граничные и начальные условия будут иметь вид:


Охлаждение изолированного провода (3)

Преобразуем выражения (3):


Охлаждение изолированного провода, (4)

Охлаждение изолированного провода. (5)


Алгоритм решения данной задачи:

с учётом начальных условий рассчитываем температуру Тi,1;

по выражению (4) рассчитываем поле температур для шага по длине j+1;

по выражению (5) рассчитываем температуры в точка на границе изоляции;

переходим к следующему шагу по длине z и повторяем расчёт сначала.

Модель №2

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

Уравнения энергии для изоляции и для жилы соответственно будут иметь вид:


Охлаждение изолированного провода (6)

Охлаждение изолированного провода

Рис. 2. Схема разбиения


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

Охлаждение изолированного провода

Начальные условия:


Охлаждение изолированного провода


Также как и в предыдущем случае для решения системы дифференциальных уравнений (4) используем метод конечных разностей. Разностные уравнения (использована явная разностная схема) имеют вид:


Охлаждение изолированного провода (7)

Охлаждение изолированного провода (8)

Охлаждение изолированного провода (9)

Алгоритм решения данной задачи:

с учётом начальных условий рассчитываем температурное поле по выражениям (7) и (8) во внутренних точках областей I и II;

по выражениям (9), рассчитывается температура в граничных точках;

переходим к следующему шагу по длине z и повторяем расчёт сначала.

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

Модель №3

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


Охлаждение изолированного провода

Рис. 3. Зависимость теплофизических характеристик от температуры


Теплофизические характеристики для текущей точки можно определить (ниже приведён пример для теплоемкости с):


Охлаждение изолированного провода

При решении задачи методом конечных разностей можно воспользоваться постановкой задачи для модели 2.

Алгоритм решения данной задачи:

с учётом начальных условий рассчитываем температуру Тi,1;

по выражениям (5) и (6) рассчитываем температурное поле во внутренних точках областей I и II;

по выражениям (7), (8) и (9) рассчитывается температура в точках на границе изоляции и между изоляцией и жилой а также в точке, лежащей в центре жилы;

пересчитываем значения теплофизических характеристик для каждой точки сечения;

пересчитываем поле температур в соответствии с выражениями (5), (6), (7), (8) и (9) для этого же сечения и пересчитываем теплофизических характеристики, повторяем расчёт до тех пор пока разность между предыдущими значениями и значениями полученными в текущем расчёте не будет меньше некоторой величины, взятой за погрешность;

переходим к следующему шагу по длине z и повторяем расчёт сначала.

Расчётная программа

uses crt;

const rg=0.00185; ri=0.0033;

l1=10; l2=10; l3=7;

tv1=80; tv2=60; tv3=20;

alpha1=1000; alpha2=1500; alpha3=2000; alphav=40;

v0=45; lamdai=0.18; lamdap=370;

roi=770; rop=9210;

ci=1940; cp=305;

tp=20; ti=250;

var t:array [1..50] of double;

t1:array [1..50] of double;

f, f1, f2:text;

n, m, n1, i:integer;

j, l:longint;

q, r, hr1, hr2, hz1, hz2, hz3:double;

procedure vanna (var f:text; l:longint; alpha, tv:real);

begin

for j:=1 to l do

begin

r:=hr1;

for i:=2 to n1–1 do

begin

t1 [i]:=t[i]+lamdap/(cp*rop*v0)*hz1*(1/r*((t [i+1] – t [i 1])/(2*hr1))+(t [i+1] – 2*t[i]+t [i 1])/(hr1*hr1));

r:=r+hr1;

end;

r:=r+hr2;

t1 [1]:=t1 [2];

for i:=n1+1 to n 1 do

begin

t1 [i]:=t[i]+lamdai/(ci*roi*v0)*hz1*(1/r*((t [i+1] – t [i 1])/(2*hr2))+(t [i+1] – 2*t[i]+t [i 1])/(hr2*hr2));

r:=r+hr2;

end;

t1 [n]:=(t1 [n 1]/hr2+alpha/lamdai*tv)/(1/hr2+alpha1/lamdai);

t1 [n1]:=(t1 [n1–1]*lamdap+t1 [n1+1]*lamdai)/(lamdap+lamdai);

for i:=1 to n do

t[i]:=t1 [i];

if j mod 10000=0 then

begin

for i:=1 to n do

begin

write (f, t[i]:2:2);

write (f, ' ');

end;

writeln(f);

end;

end;

end;

begin

clrscr;

n:=30;

n1:=18;

for i:=1 to n1 do

t[i]:=20;

for i:=n1+1 to n do

t[i]:=250;

hr1:=rg/(n1–1);

hr2:=(ri-rg)/(n-n1–1);

hz1:=0.00001;

assign (f, 'tabl.txt');

rewrite(f);

for i:=1 to n do

begin

write (f, t[i]:2:2);

write (f, ' ');

end;

writeln(f);

l:=round (l1/hz1);

vanna (f, l, alpha1, tv1);

close(f);

assign (f1,'tabl1.txt');

rewrite(f1);

for i:=1 to n do

begin

write (f1, t[i]:2:2);

write (f1,' ');

end;

writeln(f1);

l:=round (l2/hz1);

vanna (f1, l, alpha2, tv2);

close(f1);

assign (f2,'tabl2.txt');

rewrite(f2);

for i:=1 to n do

begin

write (f2, t[i]:2:2);

write (f2,' ');

end;

writeln(f2);

l:=round (l3/hz1);

vanna (f2, l, alpha3, tv3);

close(f2);

end.

Результаты расчётов


Охлаждение изолированного провода

Рис. 4. Распределение температуры по длине охлаждающих ванн


Охлаждение изолированного провода

Рис. 5. Распределение температуры по радиусу кабеля


Охлаждение изолированного провода

Рис. 6. Зависимость температуры от скорости изолирования жилы

Охлаждение изолированного провода

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


Охлаждение изолированного провода

Рис. 8. Зависимость средней арифметической температуры от скорости изолирования жилы


Охлаждение изолированного провода

Рис. 9. Зависимость температуры от коэффициента теплопроводности

Охлаждение изолированного провода

Рис. 10. Распределение температуры по радиусу для различных коэффициентов теплопроводности


Охлаждение изолированного провода

Рис. 11. Зависимость средней арифметической температуры от коэффициента теплопроводности


Охлаждение изолированного провода

Рис. 12. Зависимость температуры от удельной теплоёмкости

Охлаждение изолированного провода

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


Охлаждение изолированного провода

Рис. 14. Зависимость средней арифметической температуры от удельной теплоемкости


Охлаждение изолированного провода

Рис. 15. Зависимость температуры от температуры первой охлаждающей ванны

Охлаждение изолированного провода

Рис. 16. Распределение температуры по радиусу для различных температур первой ванны


Охлаждение изолированного провода

Рис. 17. Зависимость средней арифметической температуры от температуры первой ванны


Охлаждение изолированного провода

Рис. 18. Зависимость температуры от температуры второй охлаждающей ванны

Охлаждение изолированного провода

Рис. 19. Распределение температуры по радиусу для различных температур второй ванны


Охлаждение изолированного провода

Рис. 20. Зависимость средней арифметической температуры от температуры второй ванны


Охлаждение изолированного провода

Рис. 21. Зависимость температуры от температуры третьей охлаждающей ванны

Охлаждение изолированного провода

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


Охлаждение изолированного провода

Рис. 23. Зависимость средней арифметической температуры от температуры третьей ванны


Охлаждение изолированного провода

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

Вывод: в данной лабораторной работе была использована вторая модель для охлаждения изолированного провода. Были построены распределения температур по длине и радиусу провода. Из Рис. 4 видно, что внешние слои охлаждаются быстрей, чем внутренние, т. к. они имеют лучшие условия для охлаждения. Из Рис. 5. видно, что внутренние слои изоляции охлаждаются плохо, это происходит из-за того, что полимер имеет низкую теплопроводность, а следовательно отводить тепло из внутренних слоёв довольно сложно.

Из Рис. 6 и 7 видно, что при увеличении скорости изолирования жилы температура жилы уменьшается, а температура изоляции увеличивается. Это происходит потому, что изоляция не успевает отвести тепло и жила не успевает нагреться от воздействия на неё тепла, отводимого от изоляции.

По графикам 8 и 9 видно, что при увеличении коэффициента теплопроводности температура изоляции падает, т. к. при увеличении теплопроводности полимер быстрее отдаёт тепло.

Из Рис. 11 и 12 видно, что при увеличении теплоемкости температура увеличивается. Это объясняется тем, что с увеличением теплоемкости полимер может дольше сохранять в себе тепло.

По графикам зависимостей температуры от температур в охлаждающих ваннах видно, что при увеличении температур в ваннах, температура на выходе увеличивается, т. к. ухудшаются условия охлаждения.

По рис. 24. видно, что при увеличении толщины изоляции температура на выходе из охлаждающих ванн увеличивается. Это происходит из-за того, что требуется охлаждать больший объём полимера.

Рефетека ру refoteka@gmail.com