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

При решении задачи методом конечных разностей можно воспользоваться постановкой задачи для модели 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;