金属凝固过程界面换热系数的Tikhnonv正则化辨识.pdfVIP

cHINEsE oFcoMPuTATIoNALPHYsIcsPHⅧcs 2008年7月瑟姜焉堋 cH。NEsE，总NAL磊c。M韶队叽罴LJouRNAL V鼍苗j糍Iul，j2(x培 [文章编号】 1001．246x(2008)04．0463一07 ’ 隋大山， 崔振山 (上海交通大学国家模具cAD工程研究中心，上海200030) [摘 要】 金属凝固过程中，铸件与铸型间的界面换热系数是数值模拟所必需的边界条件之一．充分考虑非线 和Moro肋v偏差原理确定正则参数，采用灵敏度系数和Newton．Raph∞n迭代法求解该泛函．具体的案例分析表明， 该方法克服了热传导反问题的不适定性，能够保证辨识结果的稳定性和精度．与其他方法相比，随着测温误差 的增大，该方法具有较高的辨识精度． [关键词】 凝固；正则化；反问题；参数辨识；灵敏度系数；界面换热系数 [中图分类号】 0482．2 [文献标识码】 A O 引言 数值模拟技术正逐渐成为铸造、热处理等热加工工艺设计和优化的主要工具，准确完整的材料参数、初 始条件和边界条件是数值模拟所必需的已知条件，直接影响模拟精度．而缺少完备的模拟参数是数值模拟面 临的主要问题之一．与边界条件相关的重要参数——铸件与铸型间的界面换热系数，是界面接触状况、物理 条件、化学条件以及界面温度、压力等多种因素综合作用的体现，由于凝固过程中铸件和铸型之间存在气隙， 界面处同时存在热传导、热对流和热辐射等多种传热方式，目前还无法对其进行直接测量，也没有成熟完善 的数学模型进行定量描述．通常根据具体铸造工艺，假设该系数为某些经验数值或者一定的分布函数形式进 行处理…． 实验与理论相结合的热传导反算法是辨识界面换热系数的方法之一【2。]．Stolz是最早研究热传导反问题 热传导反问题的方法中旧’30．随后，一些学者在Beck提出的参数估计方法的基础上做了许多研究工作H“J． 近年来，随着计算技术的发展，不少学者采用正则化方法、动态规划技术、平滑技术、卡尔曼滤波技术、共轭梯 度法等求解热传导反问题【1¨墉1．由于热传导反问题是一个不适定问题，辨识结果对输入数据(测量温度)的 扰动非常敏感，并且金属凝固过程释放结晶潜热是一个非线性瞬态热传导问题，热电偶在测温过程的延迟效 应和衰减效应，加大了测量误差，导致不适定性更为突出口。J． 本文针对金属凝固过程数值模拟技术的发展需求，充分考虑非线性瞬态热传导过程测温数据的误差特 点，根据Tikhonov正则化理论构造了正则化泛函o 的精度和稳定性． 1 金属凝固过程热传导正问题的数学模型 由已知的材料物性参数、初始条件和边界条件求解温度场的问题称为热传导正问题．金属凝固过程温度 场的数值分析就是求解Fourier导热微分方程，并且金属凝固要释放结晶潜热，因此，可将该过程视为一个有 内热源的瞬态热传导问题，其数学模型为 [收稿日期】200r7一Ol一22；[修回日期】20昕一09一ll 学国家模具cAD工程研究中心20∞30． 万方数据 计 算 物 理 第25卷 IDc，瓦2赢l^瓦J+丽I^丽J+瓦＼^瓦J+V’．oc，警=未(A警)+未(A警)+未(A警)+o， …(1) 式中lD，．=I，郇分别为材料的密度、导热系数和比热，O为内热源，r为温度，t为时间． 金属凝固过程中，释放结晶潜热是其显著特点．常用的处理结晶潜热方法有等价比热容法、温度回升法 和热焓法等…．本文采用等价比热容法处理结晶潜热，即 o：以誓：以普警， (2) 式中，L为合金的结晶潜热，^为温度r时的固相率，是温度的函数． 求解Fourier导热微分方程的初始条件为 丁(菇，)，，z，￡)I。：o=死， (3) 式中死为初始温度． 第三类边界条件为 。蕃l。“∽-一