一种通用含热阻温度场计算方法.pdfVIP

第 6卷 第 1期 现 代 应 用 物 理 Vo1．6，NO．i 2015年 3月 M0DERN APPLIED PHYSICS M ar．2O15 一 种通用的含热阻温度场计算方法 ， 张相华 ，彭国良 ，杜太焦 ，闫 辉 高银军 (1l西jE核技术研究所 ，西安 710024； 2．激光与物质相互作用国家重点实验室，西安 710024) 摘 要：针对含 热阻温度场 的计算问题，提 出了一种通用的有 限元求解方法。由于含 热阻界 面热流密度表达式具有对流换热的形式，该方法将热阻条件视为第三类边界条件。推导了热 阻条件下的有限元方程，通过计算含热阻界面两侧单元形函数在单元相交 区域的局部积分，实 现 了区域界面热阻条件的约束。算例表明，该方法不受区域界面数量的限制，能同时适用于对 齐 网格和非对齐 网格 ，具有 良好 的通用性。 关键词 ：有限元方法；热阻；温度场；对齐网格 ；非对齐网格；局部积分 中图分类号 ：TN249；0343．6 文献标志码 ：A 文章编号 ：2095—6223(2015)01—060—06 A GeneralM ethodforTemperatureFieldCalculationwithThermalResistance ZHANGXiang—hua，PENGGuo—liang，DU Tai—jiao一，YAN Hui，GAO Yin—jun (1．NorthwestInstituteofNuclearTechnology，Xi’an 710024，China； 2．StateKeyLaboratoryofLaserInteractionwithMatter， Xi’an 710024，China) Abstract：A generalfiniteelementmethod (FEM )wasdevelopedforthecalculationoftem- peraturefieldwiththermalresistance，andtheequationswerededucedinthispaper．Asthe expressionofheatflow ratedensitythrough theinterfacewith thermalresistancetakesthe similarform totheexpressionofconvectionboundarycondition，thermalresistancecondition isconsideredasthe3rdtypeboundarycondition．Bycalculatinglocalintegralofinterfaceel— ementshapefunctions，thethermalresistanceconditionconstraintcanbeimplemented．The examplesshow thatthemethodworkswellforanynumberofinterfaceswithalignedmeshor unalignedmesh． Keywords： finiteelementmethod9‘thermalresistance；temperaturefield；alignedmesh； unalignedmesh；localintegral 含热阻温度场计算问题在工程应用 中十分常 时，需要的界面单元信息过大，导致实现较为困难等 见 ，国内外学者已提出了一系列算法 ，典型的如热阻 不足 。 单元法u ]、节点分裂法_5]、间断有 限元法 。和分段 本文将