识别材料导热系数与导温系数的温度场逆分析.pdfVIP

中国工奢撼■理学套 传热传魔学擎拳会嵌 王登刚刘迎曦李守巨李风志 (大连理工大学工程力学系工业袭备结构分析田素重点实验室116眈4) 擅蔓为了识别材料热传导系数和导温系数，采用观测温度同时识别材料导热系数和 比热与密度的乘积。在建立同时确定热传导系数和比热与密度乘积的瞬态温度场逆分析 j I 数值计算模型的基础上，引入混沌优化方法求解该优化模型．以克服阻尼牛顿等方法求 ， 、 解该模型所遇到的困难．在瞬态温度场正问题求解中采用了精细积分方法以提高计算精 ● 度．算倒验证了本文方法的有效性． 关t司热传导逆问题，导热系数，导温系数，混沌优化 l引言 导热系数A和导温系数口是材料的重要热物理性质．文【l】曾采用不稳定导热理论 拟定了同时测定热绝缘材料导温系数和寻热系数的浏试方法与技术．本文从逆问题的角 度，对由温度观测信息同时确定工作结构材料的导热系鼓和导温系数的计算方法进行研 究．基于导熟系数、导温系数之间的公式口=州(卢c)，把本文问题转化为同时确定导 热系数且和比热与密度乘积卢c．在建立同时确定丑和卢c的反演计算模型基础上．把 混沌优化方法引入到逆问题求解领域，以克撮采用阻尼牛顿等方法求解同类阔晨所遇到 的困难。为提高计算精度，采用精细积分法求解瞬态温度场有限元平衡方程． 2瞬态温度场逆问题计算模型 考虑具有实际物理意义的特求变量五和．印的上、下界限． ∥≤丑蔓Z．∥‘≤肛≤肛。 (1) 采用有限元等数值方法求解热传导定解方程．并考虑到实际中所能观测到的只能是一些 离数点处的温度值时，导热逆问题的数值计算模型可建立为 ‘E p，，‘】 nlinJ(五，肛，)=∑∑F(而，f，)一L(一，f／)，a!＆‘r』E sJ． 【c【蒯{n+k以粕)={P} (2) ∥s五≤∥，肛。≤卢S伊。 式中： r“，0)和乙(t，0)分别是在0时刻观测点而处的温度计算制和班测值； 砌。是观测点所在区域：‘、‘分别是起始观测时刻和终止观测时刻；等式约束为瞬 态温度场有限元平衡方程；【c(肛)】、晖∞)】、{P}分别是比热容矩阵、热传导矩阵和 右端向量． 3．逆问题计算模型求解的混沌优化方法 7 ———————————1F■_i百蕊■—————————～一…—J 由于模型(2)中等式约束即为瞬态温度场有限元平衡方程，对于某一组工和‘℃，可 以由该式求出观测点处的温度值r(‘，r，)，式(2)实质上是一个只有变量上下界约束的优 化问题．笔者曾经试图采用阻尼牛顿方法求解该问题．但是往往会遇到求解计算失败的 情况。这不仅与逆问题目标函数的非凸性有关，温度对五和弗这两个变量的灵敏度相 差较大，二者在计算中耦合作用也是引起该困难的原因之～．为此本文把混沌优化方法 引入反问题求解领域来求解该问题。 利用混沌搜索求解问题(2)时，首先采用式(3)建立特求变量工与混沌变量f的一一对 应关系．式(3)中：岛和4；分别是变量x的第f个分量的时、下界限，变量个数一=2， 后采用Logistic映射式(4)在[0，1]中产生一个轨迹不同的混沌变量f：进行优化搜索．式 “)中“=4．已经证明越=4时是“单片”混沌，‘在[0，1]之间历追． ‘=(t一口；)，(6I一口J) f=1，2，…，一