一种基于热辐射边界条件的三维有限元模拟方法技术

本发明专利技术属于三维热分析数值求解技术领域，涉及一种基于热辐射边界条件的三维有限元模拟方法。本发明专利技术首先对要进行热分析的器件进行建模，然后将辐射边界条件引入热传导问题，采用伽辽金残数加权的方法，得到热辐射边界条件的有限元弱形式。接着采用四面体网格剖分模型，选择二阶叠层基函数，离散有限元弱形式方程，配合Newton‑Raphson迭代方法得到有限元单元矩阵和右端向量，集成最终的方程组，最后运用科学的非线性收敛判据，经过不断迭代，快速准确地得到最终的数值计算结果。

【技术实现步骤摘要】
一种基于热辐射边界条件的三维有限元模拟方法
本专利技术属于三维热分析数值求解
，涉及一种基于热辐射边界条件的三维有限元模拟方法。
技术介绍
作为卫星大功率核心器件的空间行波管，在工作过程中不仅要受到真空高低温的影响，还要受到空间太阳辐射、红外辐射和真空冷黑背景温度的长期影响，这些热环境将影响空间行波管的电参数和热稳定性，对空间行波管的正常工作构成严重威胁。为了保证空间行波管各个部件保持在正常的工作温度，散热设计变得尤为重要。由于空间行波管暴露在外界真空环境，辐射式散热几乎是其唯一散热方式，因此有必要对空间行波管热辐射问题进行研究。目前，在数值计算领域，有关的热传导差分和有限元法著作[如数值传热学，作者：陶文铨]以及其他相关文献，大都只涉及线性边界条件(狄利克雷边界、诺曼边界、罗宾边界)的处理方法，而对于非线性边界条件如热辐射只是大略提及，没有阐明具体的解决方法。西安交通大学的张博提出了一种间接迭代的方式来处理，但是没有从有限元的本质上处理辐射边界问题，只是将非线性辐射问题转化成了诺曼边界条件来求解，只能应用在特定领域，而且针对复杂模型也缺乏必要的验证。
技术实现思路
针对上述存在问题或不足，为解决热辐射边界条件在有限元方法中的应用难题，本专利技术提供了一种基于热辐射边界条件的三维有限元模拟方法，该方法采用二阶叠层基函数以及Newton-Raphson迭代方法来进行辐射边界条件的有限元求解，可以很快地求得高精度的数值模拟结果。其具体技术方案，包括以下步骤：A.对目标器件进行建模，建立对应的几何结构模型；B.采用伽辽金加权残数法得到热辐射边界条件的有限元弱形式；C.采用四面体网格剖分求解域；D.选择基函数，离散步骤B中得到的有限元弱形式，得到热辐射边界条件的有限元方程；E.对步骤D中的有限元非线性方程组不断地进行迭代，直到其温度值满足收敛规则,收敛规则是||F-S(q+1)A(q+1)||＜ε或者ε为设定的收敛精度值。进一步优选，所述的步骤D中，在基本的三维四面体有限元基函数基础上，构造了二阶叠层基函数，相比于插值高阶基函数，叠层基函数构造方法更加简便，而且对于后面的有限元处理过程也有极大的好处，提高了有限元求解的精度。此外，步骤D中的迭代方法选用Newton-Raphson非线性迭代方法，通过这种方法的使用，使得辐射边界的高度非线性难题迎刃而解，并且可以很快地达到收敛，极大地提高了求解速度。Newton-Raphson方法最终的非线性方程组迭代形式如下：J(q)ΔA(q)＝F-S(q)A(q)(1)其中J(q)是雅可比矩阵，是本专利技术需要求解的最重要矩阵，ΔA(q)是前后两次迭代的温度差值，F是右端向量，q是迭代次数，S(q)是采用Newton-Raphson方法之前的有限元初始矩阵项，A(q)是前一次迭代的温度值。因为辐射的高度非线性，矩阵的积分求解过程很难完成，本专利技术采用了高斯数值积分的方式，采用12阶高斯积分，很方便地求解了非线性积分矩阵项，并且保证了矩阵计算的精度。本专利技术针对热辐射边界在有限元方法中的应用问题，提出了一种通用的处理方法，将热辐射边界条件与热传导控制方程联立，得到有限元弱形式，并采用高阶叠层基函数离散成非线性方程组，与Newton-Raphson迭代方法结合，通过收敛判据(||F-S(q+1)A(q+1)||＜ε或不断地迭代，直到求得最终的温度值。与现有技术相比，本专利技术可以准确、快速地求解高度非线性热辐射边界问题，并且解决了现有技术关于这方面描述不足以及一些特定解决方法的适用性不足的难题。附图说明图1是本专利技术的流程图；图2是四节点四面体单元示意图；图3是Newton-Raphson非线性迭代方法的流程图；图4是下标对应关系图。具体实施方式下面结合附图和具体实施例来详细描述本专利技术的技术方案。参照图1，一种基于热辐射边界条件的三维有限元模拟方法，包括以下步骤：A.对目标器件进行建模，建立对应的几何结构模型。B.采用伽辽金加权残数法得到热辐射边界条件的有限元弱形式。要对热辐射边界条件进行有限元求解，就必须先构建边值问题，包括热传导的控制微分方程和热辐射边界条件两部分，具体如下方程所示：其中为拉普拉斯算子，u是求解域内的温度值，k是热传导系数，Q是内部产热量，ρ是密度，c是比热容，t是时间，ε是热辐射发射率，n是法向方向，σ是Stefan-Boltzmann常数，uf是外界环境温度值。采用伽辽金残数加权法，可以得到(2)式、(3)式的加权残数表达式如下其中表示残差项，Ω表示求解域，Γ表示边界条件。从而进一步可以得到其中v1,v2是权函数，c1是任意实数。这里定义面积分和体积分如下(u,v)Ω＝∫Ω(u,v)dV(7)<u,v>Γ＝∫Γ(u,v)dS(8)其中u、v表示任意两个函数，V表示体积，S表示面积。对于(6)式中的通过格林定理，可以展开为由于v1,v2的任意性，令c1＝-1,v1＝v2，可得为了利用弱形式得到问题的近似解，首先应该选择试函数Ni(x,y,z)来代替真实解，如下式，且必须满足必要的边界条件。其中C0、Ci是任意实数，Ni为简单函数，例如低阶多项式。在伽辽金方法中，直接采用试函数自身作为权函数，即v1(x,y,z)＝Ni(x,y,z)(12)后面的有限元过程考虑稳态热传导，不考虑控制方程中的时间项，所以三维热传导方程弱形式写成C.采用四面体网格剖分求解域；采用四面体网格剖分求解域，剖分后的求解域被分割为三维四面体网格，从而将连续的几何结构空间转化为离散的网格空间。D.选择叠层基函数，离散B中得到的有限元弱形式，得到热辐射边界条件的有限元方程；如图2所示四面体单元中i,j,k,l代表四个顶点的编号，我们首先得到四个最基本的基函数：式中将(18)式、(19)式、(20)式和(21)式中的i,j,k,l轮换，得到aj,ak,al，bj,bk,bl，cj,ck,cl，dj,dk,dl。V为四面体的体积。标量二阶叠层基函数，选择N1,N2,N3,N1N2,N1N3,N2N3这六个基函数。对于有限元过程来说，把域Ω离散为M个单元之后，如同(13)式所示的弱形式定积分，可以通过将每个单元的积分贡献简单相加，即对于每一个单元来说，跟有限元方程组右端项有关系的和的求解在很多有限元基础材料中都有介绍，在此不再赘述。本专利技术的核心是用Newton-Raphson迭代方法来进行非线性辐射边界的处理，因此重点考虑的是有限元方程组左端矩阵的求解。假设非线性方程组的形式为：SA＝F(24)其中S为左端矩阵，A为待求解向量，F为右端向量。对于Newton-Raphson法是一种梯度算法，将(24)式的有限元方程组写成如下形式：f(A)＝SA-F＝0(25)其中f(A)为非线性函数。应用Newton-Raphson方法则有A(q+1)＝A(q)-[f′(A(q))]-1f(A(q))(26)其中上标q表示迭代次数，令J(q)＝f′(A(q)),并化简有J(q)ΔA(q)＝F-S(q)A(q)(27)其中J(q)为非线性函数f(A)对A的导数矩阵[称为雅可比(Jacobi)矩阵]的第q次迭代值，ΔA(q)前后两次迭代的温度差值，满足ΔA(q)＝A(q+1)-A(q)(28)其本文档来自技高网...

【技术保护点】
一种基于热辐射边界条件的三维有限元模拟方法，包括以下步骤：A.对目标器件进行建模，建立对应的几何结构模型；B.采用伽辽金加权残数法得到热辐射边界条件的有限元弱形式；C.采用四面体网格剖分求解域；D.选择基函数，离散步骤B中得到的有限元弱形式，得到热辐射边界条件的有限元方程；E.对步骤D中的有限元非线性方程组不断地进行迭代，直到其温度值满足收敛规则,收敛规则是||F‑S

【技术特征摘要】
1.一种基于热辐射边界条件的三维有限元模拟方法，包括以下步骤：A.对目标器件进行建模，建立对应的几何结构模型；B.采用伽辽金加权残数法得到热辐射边界条件的有限元弱形式；C.采用四面体网格剖分求解域；D.选择基函数，离散步骤B中得到的有限元弱形式，得到热辐射边界条件的有限元方程；E.对步骤D中的有限元非线性方程组不断地进行迭代，直到其温度值满足收敛规则,收敛规则是||F-S(q+1)A(q+1)||＜ε或者其中ε为设定的收...

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