chuanre.m

这个程序是关于连铸二冷区喷嘴模拟实验中传热系数计算的程序。

function [c,q,f0,f0m,t0,a,am] =chuanre
%chuanre Summary of this function goes here
%  t1 是1点的温度
%  t2 是2点的温度
%  t3 是3点的温度
%  c  是t0点的比热
%  q  是t0点的导热系数
%  f0 是t0点的傅立叶准则数
%  f0m是另一个方法算得的傅立叶准则数
%  t0 是表面温度
%  a  是换热系数
%  am  是另一个方法算得的换热系数
[t1,t2,t3]=textread('t.txt','%f%f%f','delimiter','');
T=[t1,t2,t3];
m=size(T);
t0=zeros(m(1)-1,1);
f0=zeros(m(1)-1,1);
f0m=zeros(m(1)-1,1);
a=zeros(m(1)-2,1);
am=zeros(m(1)-2,1);
for i=1:(m(1)-1)
    if t1(i)-2*t2(i)+t3(i)==0
        f0(i)=f0(i-1);
        t0(i)=t0(i-1);
    else f0(i)=(t2(i+1)-t2(i))/(t1(i)-2*t2(i)+t3(i));
         t0(i)=(t1(i+1)-t1(i))/f0(i)+2*t1(i)-t2(i);
    end
end 
q=11.77034+0.00833*t0+5.13709*10^(-6)*t0.^2-1.03044*10^(-8)*t0.^3+9.51184*10^(-12)*t0.^4;
c=474.73669+0.07218*t0+8.40877*10^(-5)*t0.^2-2.17759*10^(-7)*t0.^3+2.12546*10^(-10)*t0.^4;
for i=1:(m(1)-2)
    a(i)=157.2*(0.5*(t0(i+1)-t0(i))+f0(i)*(t0(i)-t1(i)))*c(i)/(10-t0(i));
    am(i)=(78.6*c(i)*(t0(i+1)-t0(i))+500*q(i))/(10-t0(i));
end
for i=1:(m(1)-1)
    f0m(i)=2.5*10^4*q(i)/(7.86*10^3*c(i));
end