一种新型的接触热阻有限元求解方法技术

本发明专利技术属于三维热传导有限元数值求解技术领域，涉及一种新型的接触热阻有限元求解方法。本发明专利技术首先对存在接触问题的器件进行建模，然后开创性地将接触热阻问题作为一种边界条件引入热传导问题，并采用伽辽金残数加权的方法，得到接触热阻问题的有限元弱形式。接着采用四面体网格剖分模型，选择叠层基函数，离散有限元弱形式方程，得到有限元单元矩阵和右端向量，然后集成总的有限元线性方程组，经过直接法或者迭代法进行求解，从而快速准确地得到最终的数值计算结果。

【技术实现步骤摘要】
一种新型的接触热阻有限元求解方法
本专利技术属于三维热传导有限元数值求解
，涉及一种新型的接触热阻有限元求解方法。
技术介绍
接触热阻(ThermalContactResistance)问题已经涉及到航天、机械制造、微电子、化工、低温超导、生物医学、核反应堆以及仪器仪表等众多科学与工程领域，其产生机理，广大学者普遍认为是由于粗糙表面间的不完全接触所造成的热流收缩而导致的。从理论上讲，完全接触的接触面应该保持同一温度，而在实际工程中，任何表面在微观上都是粗糙的，故任何固体表面之间都不可能完全接触，接触的地方直接导热，在不接触处存在空隙，产生热流收缩，存在传热阻力，即接触热阻。实际中，由于接触热阻的影响，对物体进行热分析时会造成一定的误差，甚至得到完全错误的结果。近些年来随着计算机的不断发展，人们开始利用计算机进行数值分析研究，并且将会有越来越多的研究人员使用有限元方法进行接触热阻的数值模拟分析，数值分析可以为接触热阻实验提供很好的借鉴意义。到目前为止，接触热阻在使用有限元方法进行数值模拟分析方面发表的文献不多，主要是把界面微凸体等效成长方体、圆柱体、圆锥体等，或者是在两个物体的接触面上施加一个薄层来进行过渡。以上的这些方法因为要改变原模型的结构，效率会很低下，误差也难以控制，甚至可能会得到完全错误的结果，在实际应用中会很难实施，尤其是对于形状不规则以及规模大的模型。
技术实现思路
针对上述存在问题或不足，为解决现有方法在实际应用中效率低下，误差难以控制，甚至得到完全错误的结果和难以实施的问题，本专利技术提供了一种新型的接触热阻有限元求解方法，该方法将接触热阻问题转化为有限元边界条件进行处理，高效简便地解决了接触热阻难题。其具体技术方案，包括以下步骤：A.针对热分析的对象进行建模，建立对应的几何结构模型；B.采用四面体网格对模型进行剖分；C.将接触热阻转化为边界条件，用伽辽金残数加权法得到有限元的弱形式；对于单位面积的交界面，接触热阻定义如下：其中R表示接触热阻，uA、uB表示接触面两侧温度，q”表示平均热流密度，文字表述为：接触热阻等于两个接触面温度之差除以平均热流密度。通过伽辽金方法，将接触热阻问题转化为边界条件之后，可以得到最终的有限元弱形式为其中为拉普拉斯算子，k为热传导系数，v、v1和v2为权函数，u、u1和u2为温度，δc为接触热导(接触热阻的倒数)，Ω为整个求解域，Γ12和Γ21为接触边界。D.用叠层基函数进行目标离散，得到最终的有限元方程组。E.求解步骤D中矩阵和右端项形成的线性方程组，得到最终的温度解。步骤C中，将接触热阻问题转化为热传导问题的边界条件来进行处理是本专利技术的最优之处，通过考虑接触热阻的定义以及在接触面上的热通量对于接触面两侧来说是相同的(接触面两侧的温度值是不同的)这两点，开创性地将其转化为有限元边界条件，进而用有限元方法进行处理。综上所述，本专利技术从接触热阻本质出发，虽然接触导致接触面两侧的温度不一致，但是接触面上的热通量是连续的，这里将接触热阻问题等效为热传导边值问题，用有限元方法进行求解。相比于以前的接触热阻处理方法，本专利技术最大的优势就是在不改变原模型的基础上，非常规范、高效地对接触热阻问题进行了求解，在实际工程领域非常易于实施。附图说明图1是本专利技术的流程图；图2是接触热阻作为有限元边界条件进行处理的示意图；图3是四节点四面体单元示意图；图4是接触单元计算示意图。具体实施方式下面结合附图和具体实施例来详细描述本专利技术的技术方案。一种新型的接触热阻有限元计算方法，包括以下步骤：A.针对热分析的对象进行建模，建立对应的几何结构模型。B.采用四面体网格对模型进行剖分；剖分后的求解域被分割成三维四面体网格，从而将连续的几何结构空间转化为离散的网格空间。C.将接触热阻转化为边界条件，用伽辽金残数加权法得到有限元的弱形式。考虑一个域(R3表示三维空间)，热传导的控制微分方程由下式给出:其中，u代表在有限区域上随时间变化的温度分布，为拉普拉斯算子，k是热传导系数，Q是内部产热量，ρ是密度，c是比热容，t是时间，Ω表示求解域。为了说明问题，这里将求解域分为两个区域，进而考虑两个区域接触面上的接触热阻。就如图2中所示的那样：Ω:＝Ω1∪Ω2，将表示子域的接触面表示为外部边界定义为另外，表示单位法向量，方向是在边界Γ12上指向区域Ωi外部。为了更好地描述，这里定义体积分和面积分:(u,v)Ω＝∫Ω(u,v)dV(4)<u,v>Γ＝∫Γ(u,v)dS(5)其中u、v表示任意两个函数，V表示体积，S表示面积。由接触热阻的定义，对两个区域分别有：其中δc是两个区域交界面上的接触热导值(接触热阻的倒数)，k1,k2为两个区域的热传导系数，u1,u2为接触面上的温度，n1,n2为法线。考虑稳态热分析过程，即去掉(3)中的时间项，也不讨论内部产热项Q。由伽辽金方法可以得到稳态热分析控制微分方程的残差表达式：其中RΩ表示稳态热分析控制微分方程的残差。同样地，两个区域的残差表达式可以写为：其中Rc1,Rc2表示两个区域上接触热导表达式的残差。由伽辽金残数加权法，根据(8)(9)(10)式，可以写为如下的表达式其中c1,c2是任意实数，v,v1,v2为权函数。由格林公式，可以展开为同样地，在接触面上，(11)式后两项可以展开为由于c1,c2,v1,v2的任意性，令c1＝c2＝-1，v＝v1＝v2，可以将(11)式写为如下的表达式，也就是最终得到的有限元伽辽金弱形式。D.用叠层基函数进行离散，得到最终的有限元方程组。如图3所示四面体单元中i,j,k,l代表四个顶点的编号，我们得到四个最基本的基函数，具体推导过程为一种公知过程，这里不再阐述：式中将(20)式、(21)式、(22)式、(23)式中的i,j,k,l轮换，得到aj,ak,al，bj,bk,bl，cj,ck,cl，dj,dk,dl。V为四面体的体积。为了使四面体的体积不为负值，单元节点编号i,j,k,l必须依照一定的顺序。在右手坐标系中，当i,j,k的方向转动时，右手螺旋应指向l的方向。确定了基函数为N1,N2,N3,N4，对于有限元过程来说，把域Ω离散为M个单元之后，如同(15)所示的弱形式定积分，可以通过将每个单元的积分贡献简单相加，即由(25)式，定义单元矩阵为Sij＝Kij+Cij，其中Kij为矩阵K的计算在很多基础的热传导有限元书中都有介绍，本专利技术不再赘述，主要介绍接触热阻矩阵项C，以图4所示的两个四面体为例进行说明。两个四面体为abcd和efgh，三角面bcd和efg是同一个面，即其接触面。C矩阵可以写为如下的表达式其中Γbcd、Γefg分别表示三角面bcd和三角面efg，代表三角面bcd的基函数，代表三角面efg的基函数。按照上面的公式就可以建立接触面两侧之间的关系，形成最终的矩阵。E.求解步骤D中矩阵和右端项形成的线性方程组，得到最终的温度解。在步骤D中已经求得有限元的矩阵和右端向量，可以形成最终的线性方程组，有很多方法可以进行方程组的求解，直接法和迭代法都在很多文献资料里有详细的论述，本专利技术不再进行具体描述。求解完有限元方程组，就得到了最终的温度解，可以发现接触面两侧温度值出现了跳跃，接触热阻问题得以本文档来自技高网...

【技术保护点】
一种新型的接触热阻有限元求解方法，包括以下步骤：A.针对热分析的对象进行建模，建立对应的几何结构模型；B.采用四面体网格对模型进行剖分；C.将接触热阻转化为边界条件，用伽辽金残数加权法得到有限元的弱形式；对于单位面积的交界面，接触热阻定义如下：

【技术特征摘要】
1.一种新型的接触热阻有限元求解方法，包括以下步骤：A.针对热分析的对象进行建模，建立对应的几何结构模型；B.采用四面体网格对模型进行剖分；C.将接触热阻转化为边界条件，用伽辽金残数加权法得到有限元的弱形式；对于单位面积的交界面，接触热阻定义如下：其中R表示接触热阻，uA、uB表示接触面两侧温度，q”表示平均热流密度，文字表述为：接触热阻等于两个接触面温度之差除以平均热流密度；通过伽辽金方法，将接触热阻问题转化为边界条件之后，得到最终的有限元弱形式为

【专利技术属性】
技术研发人员：徐立，谢鹏，李斌，杨中海，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51