哈工大传输原理半无限大物体温度场计算机模拟无限长铝棒.docxVIP

包套中铝棒端面受热过程中内部温度场数值模拟院系：哈尔滨工业大学材料学院姓名：日期：2014年11月30日指导教师：闫牧夫包套中铝棒端面受热过程中内部温度场数值模拟摘要：物体在加热过程中，其内部温度场不断变化，包含相变时，情况更复杂。本文通过解析方法对500度和700度铝的内部温度进行了计算机拟合，所得结果可以用于实际参考。关键词：相变；Neumann解；温度场问题提出实际生产过程中，物体温度分布是非常重要的状态特征。一方面它通过影响变形抗力而影响物体变形的力学性能参数；另一方面，它通过影响金属显微组织变化，而影响产品的组织和性能。在受热过程中，物体的温度对变形程度，内部显得组织变化及冷却后物体的力学性能都有重要影响。物体内部的应力，应变，金属的最终性能和外形质量都与温度分布有密切的关系。问题重述具有初始温度，半无限长包套中的纯铝棒，其中一端面受热瞬时达到和，并保持恒定。求温度场，熔化温度和总受热层层厚度。问题求解半无限大物体模型数学描述模型的解析解化简为温度分布情况将数据代入解析式利用origin可以画出如下图像受热层厚度工程上可辨别的温度为所以本文取此值作为温度上升的标志查误差函数表有时，所以受热层的厚度可以由以下公式计算问题求解数学模型：根据模型，并假定没有对流换热现象发生，可建立相应的数学模型。液相区固相区其中：--物质的导温系数--固相；--液相--相界面位置函数熔化过程中试样没有过热现象发生，则相界面上边界条件为相界面能量平衡方程式为其中：--物质的密度--物质的相变潜热分析解固态状态温度公式液态铝温度公式相界位置式中：为密度比，并且；为固体初始温度为升温后温度为物质的熔点为熔化常数，并可由下式确定式中：为定压热容为物质的相变潜热利用模型简化分析，只有固相区温度待定，两区域问题转化为一区域问题，确定凝固常数可简化为利用迭代法求解，在的误差范围内熔化层厚度相界位置即为熔化层厚度由前述相界位置公式及受热层厚度公式受热层厚度受热层厚度由两部分组成：①熔化层厚度；②未熔化但有温度上升的部分所以总受热层厚度温度场式(2-5)描述固态温度场分布式(2-6)描述液态温度分布场分布令为一定值，在时间内，该位置可能是固相区，也可能是液相区，分两种情况绘制温度场分布曲线相关数据初始温度液相导热系数熔化温度液相密度壁面温度相变潜热固相导热系数比热容固相密度模拟结果与分析模拟结果①铝棒无熔化层；其受热层厚度与时间存在式(1-7)所表示的)函数关系；其温度场分布为式(1-4)，温度分布与时间误差函数存在线性关系。②铝棒内熔化层随时间变化关系式(2-11)，熔化层厚度与时间的平方根成正比例关系；受热层厚度与时间变化关系式(2-12)，相比式(2-11)添加了有导热但未熔化部分的厚度；温度场分布式(2-5)、(2-6)，温度场是时间和距离的函数，同时距离小于相界，多个变量之间相互制约，不能明确地表达出每两个变量之间的关系。过程评价①半无限大物体模型的假设之一是与轴垂直的面上温度均相同，而本题中并不能证明这一条件成立，尤其是在阶段，因为固液态密度不同，相界不会是规则几何面；②熔化常数的方程从(2-8)简化到(2-9)用到了模型，而该模型的前提之一是液体与固体密度相等，在固液两相密度相差不大时，可以通过此模型简化；③熔化常数的方程(2-9)为超越方程，该方程没有确定解，利用迭代法，在误差小于时，其值为0.21；④本文中将铝的导热系数、密度、相变潜热、比热容等与状态有关的参数本文中按常数处理，简化了模型，但同时也影响了模型的精度；⑤熔化部分铝液内部存在导热和对流等多种情况，本文只考虑液体内部导热，不考虑液体的其他热传递规律。参考文献[1]张适阔.铝合金高温储热材料相变储热机理研究.武汉：武汉理工大学.2009[2]陆煜，程林.传热原理与分析.北京：科学出版社.1997.12[3]杨世铭，陶文铨.传热学（第四版）.北京：高等教育出版社.2006.8[4]吴孟怀.液态金属热扩散研究.西北工业大学学报.第10卷第1期.1992.1[5]李长庚，周孑民，刘健君.金属材料相变导热系数的研究――实验与数值模拟的结合.矿冶工程.2005年01期