传热与接触两类问题耦合作用的有限元分析Ξ.PDFVIP

Vol. 21 　No. 3 　　 　　　　　　 　 　　　　　　　　　 　 第21卷 第3期 固 体 力 学 学 报 2000年 　9 月 ACTA MECHANICA SOLIDA SINICA September 2000 传热与接触两类问题耦合作用的有限元分析 张洪武 　　顾元宪 　　钟万勰 (大连理工大学工程力学系、工业装备与结构分析国家重点实验室 ,大连 ,116024) 摘 　要 　考虑传热接触耦合作用的热力学分析问题大量存在于工程中 ,分析的难点是必须 考虑热与可移动接触边界间的耦合作用. 针对这类问题的求解 ,该文给出了接触边界上热交换 与温度边界条件 ,并在此基础上建立了两类变分方程 ,一类是热力学变分泛函 ,其考虑了接触区 域对结构热传导的影响 ;另一类是二维热弹性接触分析的参数变分原理 ,可以方便地对接触问 题进行求解. 文中给出有限元分析的离散公式 ,并进一步给出两类问题耦合分析的迭代算法 ,其 中接触分析的惩罚因子是可以消除的 ,数值结果验证了该文的理论与算法. 关键词 　传热 ,接触分析 ,迭代法 ,规划法 ,有限元法 1 　引言 将热与接触两类问题相结合进行耦合分析近来有代表性的工作主要由文[1～4 ] 给出. 文[1 ]主要侧重研究界面的热力学摩擦特性及本构关系 ,文[2 ]给出热接触问题模型的数学 描述 ,且自认为是进行该领域有限元分析研究的较早的系统性工作 ,文[3 ]给出热摩擦接触 问题的有限元求解模型,所采用的计算方法是对两类问题均采用迭代分析技术进行求解 ,文 [4 ]则以互补关系进行热接触问题的求解 ,特点是对热接触面给出以热流与温度差为变量的 互补关系 ,但这在另一方面使接触面的传热条件过于严格与保守 ,无法考虑接触间隙介质可 能进行的热交换性态. [5 ] 本文基于时域的逐步温度增量分析与参变量变分原理的二次规划算法 进行热与接触 问题两类问题的求解,为达到耦合分析的目的 , 引入迭代技术进行两类问题的交替求解过 程. 对热传导分析采用常规的分析技术 ,但为模拟接触面与接触区域的传热特性 ,文中建立 了接触面热交换与温度关系模型 ,这一模型包含了文[4 ] 的接触面模型 ,且较文[4 ] 的模型更 为灵活 ,可以模拟间隙内介质的导热特性. 所建立的接触面传热关系可由文中给出的能量泛 函的变分过程导出. 另外 ,对于采用由参变量变分原理导出的参数二次规划算法 ,所导出的 接触刚度阵中的惩罚因子在热接触分析中同样可以消除 ,从而保证了接触问题的求解精度. 文中采用迭代过程进行耦合问题求解 ,实际计算表明这种迭代过程在通常情况下有很好的 收敛特性. 文中给出了两个有代表性的数值算例 ,计算结果在验证本文计算模型与算法的同 时 ,进一步表明耦合计算对传热接触分析的重要性. 作为研究工作的第一步 ,本文只讨论二 维问题. ( ) ( ) 国家自然科学基金、国家重点基础研究专题经费 、国家教育部优秀青年教师基金 资助课题. 收到第 1 稿 收到修改稿. © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 　　　　　　　　　　　　　　　　　固体力学学报 　　　　 　　　　　　　　　　　2000 年 第 21 卷 2·18 · 2 　热接触耦合问题的