T_uo_ot_tau_cylinder_1
Определяет температуру конца ускоренного охлаждения для круглого проката

function T_uo = T_uo_ot_tau_cylinder_1(T_0, d, tau)
% Определяет температуру конца ускоренного охлаждения T_uo [град. C]
% круглого проката диаметром d [мм] с начальной температуры T_0 [град. C]
% при длительности охлаждения tau [сек].

% При вычислениях будем использовать сплайн-интерполяцию трех зависимостей
load Cylinder_1_M_Spline.mat i_ot_T Tetta_ot_Fo A_ot_i

% Поскольку величина приведенного коэффициента температуропроводности a_pr
% зависит от конечной температуры, будем использовать метод последователь-
% ного приближения. Начинаем расчеты с величины приведенного коэффициента
% температуропроводности равного 5 мм.кв./с
a_pr_1=5;
a_pr=0.1;

while any(abs(a_pr_1-a_pr)./a_pr.*100 >= 0.1)
% Расчеты производятся до тех пор, пока ошибка в определении a_pr
% не будет ниже 0.1 %
a_pr = a_pr_1;
% Число Фурье при принятом нами a_pr
Fo=a_pr.*tau./(d/2).^2;

% Относительное теплосодержание в конце ускоренного охлаждения
% как функция, зависящая от времени охлаждения (числа Фурье)
Tetta_i = spline(Tetta_ot_Fo(1,:), Tetta_ot_Fo(2,:), Fo);

% Теплосодержание проката в начале ускоренного охлаждения
i_0 = spline(i_ot_T(1,:), i_ot_T(2,:), T_0);

% Теплосодержание стали, соответствующее температуре поверхности
% раската в процессе ускоренного охлаждения
T_p = 40; % Принимаем температуру охлаждающей воды равной 40 град. C
i_p = spline(i_ot_T(1,:), i_ot_T(2,:), T_p);

% Теплосодержание раската в конце ускоренного охлаждения
i_1 = Tetta_i.*(i_0 - i_p) + i_p;

% Интегральная величина коэффициента температуропроводности для
% теплосодержаний i_0 и i_1 составит
A_0 = spline(A_ot_i(1,:), A_ot_i(2,:), i_0);
A_1 = spline(A_ot_i(1,:), A_ot_i(2,:), i_1);

% Приведенный коэффициент температуропроводности
a_pr_1=(A_1-A_0)./(i_1-i_0);
end

% Температура конца ускоренного охлаждения как функция от теплосодержания
% стали
T_uo = spline(i_ot_T(2,:), i_ot_T(1,:), i_1);

К началу страницы

Web-сайт “Термист” (termist.com)
Термомеханическое упрочнение арматурного проката

Отсутствие ссылки на использованный материал является нарушением заповеди "Не укради"

Редактор сайта: Гунькин И.А. (termist.com@gmail.com)