【传输原理】半无限大铝棒加热讲义.docx

Harbin Institute of Technology传输原理论文·一课程名称：传输原理设计题目：包套中纯铝棒端面受热过程中内部温度场数值模拟院系：材料科学与工程班级：1219001班设计者：缪克松学号：1121900133哈尔滨工业大学摘要：包套中金属棒定壁温端面加热是典型的半无限大问题，“定壁温”符合激光加热的特点，包套隔绝外部环境，其内部温度场的数值分布直接关系该类材料诸如比热容、导热系数、密度等一系列本征物性参数，其温度场的情况对于材料导热或隔热的情况下的设计具有指导意义，本文以铝棒为对象，进行了有无相变两类情况的模拟分析。关键词：半无限大空间温度场相变模拟设计题目具有初始温度℃，半无限长包套中的纯铝棒，其中一个端面受热瞬时打到500℃和700℃，并保持恒定。求温度场，熔化层厚度和总受热层层厚度。问题分析及解决步骤设计此题目相当于半无限大问题的非稳态导热问题，其中铝的熔点为660℃，即500℃和700℃分别涉及铝棒有无相变“s → l” 过程，可以建立数学模型求解，然后利用matlab编程实现数值分析并输出图像，最后分析模拟数据，得出结论。铝的物性参数为方便计算，取液态物性参数和固态物性参数为定值。参数比热容J/( kg·℃ )密度kg/m3导热系数W/(m·℃)固态0.990*1032.638*103236.6液态1.076*1032.35*10362.8熔点℃660融化热J/kg3.93*105500℃端面温度求解数学模型：引入过余温度，一维温度场方程： 0 x∞,τ 0初始条件：,x ≥0,τ= 0边界条件： , x= 0, τ 0运用相似性解法，引入相似性变量：变换方程为常微分得到：满足以上条件的方程解为：其中：编程生成图像实验材料：铝初始设置：初始温度20℃，边界温度瞬时达到500℃实验条件：加热10s、100s、1000s、10000s。图像如下：此时受热层厚度和时间关系表格存在：时间/s10100100010000厚度/m0.361.133.5711.28程序如下：clcclearclft0=20; %初始温度t1=500; %边界温度p=2.638E3; %密度c=0.99E3; %比热容l=236.6; %导温系数a=l/(p*c);T=10;X=0:0.01:12;x=[10 100 1000 10000;0 0 0 0;];%存储受热层厚度t=t0+(t1-t0)*(1-erf(0.5*X/sqrt(a*T)));plot(X,t,b-);x(2,1)=min(X(find(t==20)));hold onxlabel(位置/m)ylabel(温度/℃)title(不同时间铝加热过程中温度场分部（无相变）)grid on;text(x(2,1),30,{\downarrow},Color,b)T=100;t=t0+(t1-t0)*(1-erf(0.5*X/sqrt(a*T)));plot(X,t,g-);x(2,2)=min(X(find(t==20)));text(x(2,2),30,{\downarrow},Color,g)T=1000;t=t0+(t1-t0)*(1-erf(0.5*X/sqrt(a*T)));plot(X,t,r-);x(2,3)=min(X(find(t==20)));text(x(2,3),30,{\downarrow},Color,r)T=10000;t=t0+(t1-t0)*(1-erf(0.5*X/sqrt(a*T)));plot(X,t,k-);x(2,4)=min(X(find(t==20)));text(x(2,4),30,{\downarrow},Color,k);legend(10s,100s,1000s,10000s)小结：随着时间增加，受热层厚度增加，在单次加热过程中（时间一定），温度场分布越远离加热端面越低，直到等于初始温度，降温幅度逐渐减小。700℃端面温度求解数学模型：在液相区内：0xX(τ),τ0tl = t0 x=0, τ0在固相区内：0xX(τ),τ0ts → t∞ x→∞, τ0ts = t∞ x0, τ=0分界面x=X(τ)处的耦合条件为：ts = tl = tm x = X(τ), τ0x=X(τ), τ0引入待定常数η：得到液相区和固相区温度分布：以及计算参数η的超越方程：编程生成图像实验材料：铝初始设置：初始温度20℃，边界温度瞬时达到700℃实验条件：加热10s、100s、1000s、10000s。图像如下：此时融化层厚度关系表格：时间/s10100100010000融化层厚度/mm2476238753受热层厚度过大不予以计算，理论上计算过程与500℃中相同。程序如下：clea