chuanre_12.m

这是一套完整的求解连铸二冷区喷嘴模拟试验传热系数的程序

function [c,q,c0,q0,f0,t0,a] =chuanre_12
%chuanre Summary of this function goes here
%  t1 是1点的温度
%  t2 是2点的温度
%  t3 是3点的温度
%  c  是t1点的比热
%  q  是t1点的导热系数
%  c0  是试件比热
%  q0  是导热系数
%  f0 是t0点的傅立叶准则数
%  t0 是表面温度
%  a  是换热系数
[t1,t2,t3]=textread('t.txt','%f%f%f','delimiter','');
T=[t1,t2];
m=size(T);
t0=zeros(m(1)-1,1);
f0=zeros(m(1)-1,1);
a=zeros(m(1)-2,1);
q=11.77034+0.00833*t1+5.13709*10^(-6)*t1.^2-1.03044*10^(-8)*t1.^3+9.51184*10^(-12)*t1.^4;
c=474.73669+0.07218*t1+8.40877*10^(-5)*t1.^2-2.17759*10^(-7)*t1.^3+2.12546*10^(-10)*t1.^4;
for i=1:(m(1)-1)
    t0(i)=0.3144*c(i)*(t1(i+1)-t1(i))/q(i)+2*t1(i)-t2(i);
end 
q0=11.77034+0.00833*t0+5.13709*10^(-6)*t0.^2-1.03044*10^(-8)*t0.^3+9.51184*10^(-12)*t0.^4;
c0=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)-1)
    f0(i)=2.5*10^4*q0(i)/(7.86*10^3*c0(i));
end
for i=1:(m(1)-2)
    a(i)=157.2*(0.5*(t0(i+1)-t0(i))+f0(i)*(t0(i)-t1(i)))*c0(i)/(10-t0(i));
end