coeffp.m

用于计算脉冲激光照射下物体的温度分布

function f=CoeffP(i,j,l)
%    计算控制容积结点上的扩散系数，是温度（和浓度）的函数
%    输入结点位置（i，j）,时层序号l
%    返回该结点上的扩散系数值
%    此公式结果单位为um^2/s
global r TT tn;
%load D_Dcoeff.mat;
%load N_Dcoeff.mat;
%load T_Dcoeff.mat;
load pp.mat
rdd=50;                    %热斑半径，认为近似与激光焦斑半径相等%单位：um
TTT=TT+273.15;


%k=1.33e-23;
%f=0.602153e-5*exp(-1.56e-19/(k*((TT(l)*exp(-(r(j)/rdd).^2)))))*1.0e8;%um^2/s
%f=0.602153e-5*exp(-1.56e-19/(k*((TT(l)*exp(-(r(j)/rdd).^2)))));%cm^2/s

 m=1;    

 f=10.^(polyval(pp(m,:),((1./(TTT(l)*exp(-(r(j)/rdd).^2))))))*1.0e8;%um^2/s
 %f=10.^(polyval(pp(m,:),((TT(l)*exp(-(r(j)/rdd).^2)))))*1.0e8;%um^2/s
 if f==0,f=eps;end;