CN111898263B - 工程含热源结构热传导情况求解方法以及装置 - Google Patents

工程含热源结构热传导情况求解方法以及装置 Download PDF

Info

Publication number
CN111898263B
CN111898263B CN202010720362.3A CN202010720362A CN111898263B CN 111898263 B CN111898263 B CN 111898263B CN 202010720362 A CN202010720362 A CN 202010720362A CN 111898263 B CN111898263 B CN 111898263B
Authority
CN
China
Prior art keywords
boundary
heat conduction
temperature
heat
equation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010720362.3A
Other languages
English (en)
Other versions
CN111898263A (zh
Inventor
王桥
刘彪
周伟
岳强
黄诚斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202010720362.3A priority Critical patent/CN111898263B/zh
Publication of CN111898263A publication Critical patent/CN111898263A/zh
Application granted granted Critical
Publication of CN111898263B publication Critical patent/CN111898263B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Geometry (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computer Graphics (AREA)
  • Control Of Temperature (AREA)

Abstract

本发明提供工程含热源结构热传导情况求解方法以及装置,能够高效、精确的获得工程结构的温度变化和热传导情况。该方法包括:步骤1.利用温度探测器探测待求解含热源工程结构的外部边界温度值;步骤2.在结构温度热传导求解器中,利用建模软件建立待求解模型,输入待求解含热源工程结构的材料参数、网格化分数、网格类型以及温度探测器探测得到的边界温度值等信息,然后输出模型数据信息文件;步骤3.建立热传导控制方程;步骤4.建立温度积分方程;步骤5.将域积分转化为边界积分;步骤6.建立遍及积分方程的求解矩阵;步骤7.基于步骤6的求解矩阵,得到模型边界和内部各点的温度值,进而得到工程结构内部的温度变化和热传导情况。

Description

工程含热源结构热传导情况求解方法以及装置
技术领域
本发明属于固体热传导技术领域,具体涉及一种工程含热源结构热传导情况求解方法以及装置。
背景技术
含热源的工程结构可具体理解为内部含有热源加热或者由自热材料构成的结构。此类工程结构在日常生活中普遍存在,例如:在火灾时,房屋结构会被高温失火源加温,致使内部结构因高温发生变形。此外,自热材料构成的结构在各大工程结构中也广泛存在。由于自热材料在使用过程中,由于物理、化学作用,会自身产热,而温度过高则会影响工程结构的安全性,因此分析材料结构由于温度变化带来的影响非常关键,该问题也受到了国内外学者的关注。以工程中被广泛应用并且极具代表性混凝土为例,大体积混凝土工程结构的温度控制问题一直是受到生产单位极度重视的工程问题,如若处理不当,温度过高,则会产生裂缝对工程结构带来巨大的安全隐患。
边界元法是工程中最常用的方法之一。与有限元法相比,边界元法的一个优点是可以将问题降低一维。边界元法可以节省大量的计算资源和时间,因为只需要对边界进行网格划分,并且模型重建过程也较为简单。然而,由于边界积分方程(BIE)无法知道域中的信息,并且对于非齐次问题存在域积分,这一优点也是其缺点。同样的缺点也存在于许多其他的边界类型方法中,如边界节点法(BNM)、伽辽金边界节点法(GBNM)和混合边界节点法(混合BNM)。
由于需要考虑结构内部的热源,在传统的热传导边界积分方程中会出现域积分,由于域积分的存在,需要对问题与内进行离散计算,故而大大降低了计算效率,边界元法将失去只用单元进行边界离散化的优势。因此寻找一种精确并且高效的域积分处理方法是十分必要的。
许多研究者都将注意力集中在避免体单元的方法上,目前已经提出了几种方法,如双互易法(DRM)。在DRM中,非齐次项可以用径向基函数(RBF)等一系列函数来逼近,并应用第二个互易性将域积分转化为边界积分。只有域或边界上的点需要提供由非齐次项表述的信息。然而,DRM的精度在很大程度上取决于域点的分布和位置,以及用于近似非齐次项的函数类型。此外,复杂域中的点的排列可能不容易,特别是对于三维问题。类似于DRM的方法是多重互易法(MRM),其中互易定理通过一系列高阶基本解反复应用,将域积分转化为边界。另一种方法是径向积分法(PSM),该方法可以将域积分转化为边界积分和径向界积分,但是当采用高斯积分进行计算径向积分时,径向积分法则比较耗时。
发明内容
本发明是为了解决上述问题而进行的,目的在于提供工程含热源结构热传导情况求解方法以及装置,能够高效、精确的获得结构的工程结构内部任意一点任意时刻的温度,进而可以准确地反映工程结构的温度变化和热传导情况。
本发明为了实现上述目的,采用了以下方案:
<方法>
本发明提供一种工程含热源结构热传导情况求解方法,其特征在于,包括:
步骤1:利用温度探测器探测待求解含热源工程结构的外部边界温度值;
步骤2:在结构温度热传导求解器中,利用建模软件建立待求解模型,输入待求解含热源工程结构的材料参数、网格化分数、网格类型以及温度探测器探测得到的边界温度值等信息,然后输出模型数据信息文件;
步骤3:建立热传导控制方程
含热源结构的热传导控制方程式为:
Figure BDA0002599741860000021
式中,k为热传导系数;θ为温度;B(x)为含热源结构热源项;x为研究域内的点;i,j,k为取值1-3的整数;
步骤4:建立温度积分方程
位移边界方程为:
Figure BDA0002599741860000022
式中,Γ为问题研究域Ω的边界;
Figure BDA0002599741860000023
为正则化的温度值,具体表达为
Figure BDA0002599741860000024
q为热通量;x和y代表问题域Ω内的源点和场点;u*(x,y)和q*(x,y)为基本解,具体为:
Figure BDA0002599741860000031
Figure BDA0002599741860000032
式中,r代表源点和场点之间的距离;n为边界Γ的单位外法线向量;
步骤5:将域积分转化为边界积分
对于公式(2)中的域积分
D1=∫Ωu*(x,y)B(y)dΩ(y) (5)
基于直线积分法,利用如下公式转化为等效边界积分:
Figure BDA0002599741860000033
式中,nl1是边界Γ上外法向量n在yl1-axis方向上的元素,所有起始点可以定义在同一个直线{yl1=m}上,a=a(yl2,yl3)=m为任意常数值,l1和l2分别取2和3中不同值,把边界离散为N个单元;此外
Figure BDA0002599741860000034
式中,
Figure BDA0002599741860000035
代表第i个边界积分点,t{t,y2,y3}∈BM代表积分点的个数,wi
Figure BDA0002599741860000036
分别代表第i个边界积分点在yl1方向的权重和外法向量,k(x,y)是定义在域内,对任意x,y∈Ω0,x≠y都连续的微奇异核函数,并且对于x,y∈Ω0,x≠y,都存在一个常数C1使得
|k(x,y)|≤C1r,0≤β<3 (8)
步骤6:建立遍及积分方程的求解矩阵
经过离散,除去域积分的位移边界积分方程可写为如下形式
Figure BDA0002599741860000037
式中,Γe为e个单元的边界;对于光滑边界,Ci=0.5;对于内部点,Ci=0;Ne为离散单元数;并且
Figure BDA0002599741860000041
通过整合,可得以下矩阵方程:
Hu=Gt (11)
将已知的边界条件代入上式,即可得
Ax=y (12)
步骤7:基于步骤6的求解矩阵,利用高斯消去求解法即可得到模型边界和内部各点的温度值,进而得到工程结构内部的温度变化和热传导情况。
<装置>
另外,本发明还提供了一种工程含热源结构热传导情况求解装置,其特征在于,包括:热传导求解器,根据上文<方法>中所描述的工程含热源结构热传导情况求解方法得到待求解含热源工程结构的温度变化和热传导数据;输入显示器,与热传导求解器通信相连,让操作者输入操作指令,并显示相应的信息;以及控制器,与热传导求解器、输入显示器通信相连,控制热传导求解器和输入显示器的运行。
进一步,本发明所提供的工程含热源结构热传导情况求解装置,还可以具有这样的特征,热传导求解器包括:
探测单元,通过温度探测器探测待求解含热源工程结构的外部边界温度值;
建模单元,与探测单元通信相连,在结构温度热传导求解器中,建立待求解模型,并基于待求解含热源工程结构的材料参数、网格化分数、网格类型以及温度探测器探测得到的边界温度值信息,输出模型数据信息文件;
热传导方程建立单元,与建模单元通信相连,建立热传导控制方程;含热源结构的热传导控制方程式为:
Figure BDA0002599741860000042
式中,k为热传导系数;θ为温度;B(x)为含热源结构热源项;x为研究域内的点;i,j,k为取值1-3的整数;
温度积分方程建立单元,与热传导方程建立单元通信相连,建立温度积分方程;
Figure BDA0002599741860000051
式中,Γ为问题研究域Ω的边界;
Figure BDA0002599741860000052
为正则化的温度值,具体表达为
Figure BDA0002599741860000053
q为热通量;x和y代表问题域Ω内的源点和场点,u*(x,y)和q*(x,y)为基本解,具体为:
Figure BDA0002599741860000054
Figure BDA0002599741860000055
式中,r代表源点和场点之间的距离;n为边界Γ的单位外法线向量;
转化单元,与温度积分方程建立单元通信相连,将域积分转化为边界积分;对于公式(2)中的域积分:
D1=∫Ωu*(x,y)B(y)dΩ(y) (5)
基于直线积分法,利用如下公式转化为等效边界积分:
Figure BDA0002599741860000056
式中,nl1是边界Γ上外法向量n在yl1-axis方向上的元素,所有起始点定义在同一个直线{yl1=m}上,a=a(yl2,yl3)=m为任意常数值,l1和l2分别取2和3中不同值,把边界离散为N个单元;此外
Figure BDA0002599741860000057
式中,
Figure BDA0002599741860000058
代表第i个边界积分点,t{t,y2,y3}∈BM代表积分点的个数,wi
Figure BDA0002599741860000059
分别代表第i个边界积分点在yl1方向的权重和外法向量,k(x,y)是定义在域内,对任意x,y∈Ω0,x≠y都连续的微奇异核函数,并且对于x,y∈Ω0,x≠y,都存在一个常数C1使得
|k(x,y)|≤C1r,0≤β<3 (8)
求解矩阵建立单元,与转化单元通信相连,建立遍及积分方程的求解矩阵;
Figure BDA0002599741860000061
式中,Γe为e个单元的边界,;对于光滑边界,Ci=0.5;对于内部点,Ci=0;Ne为离散单元数;并且
Figure BDA0002599741860000062
进一步可得以下矩阵方程:
Hu=Gt (11)
将已知的边界条件代入上式,即可得
Ax=y (12)
求解单元,与求解矩阵建立单元通信相连,基于求解矩阵,求解得到模型边界和内部各点的温度值,进而得到工程结构内部的温度变化和热传导情况。
进一步,本发明所提供的工程含热源结构热传导情况求解装置,还可以包括:图像生成器,与热传导求解器、输入显示器、控制器均通信相连,根据热传导求解器建立的待求解模型生成含热源工程结构的三维模型结构图,并根据热传导求解器求解得到的数据生成反映温度高低和变化情况的热图像;其中,输入显示器对图像生成器生成的图像进行显示,并根据操作指令对热传导求解器求解得到的具体数值进行显示。
进一步,本发明所提供的工程含热源结构热传导情况求解装置,还可以包括:判断预警器,与热传导求解器、输入显示器、控制器均通信相连,存储有温度阈值,并基于工程结构内部的温度变化和热传导情况对该工程结构各部位温度是否超出温度阈值进行判断,并在判断为是的情况下,发出预警信息;其中,输入显示器还将预警信息显示在三维模型结构图的相应位置处。
发明的作用与效果
本发明所提供的工程含热源结构热传导情况求解方法以及装置,可以在探测到工程结构温度边界条件以及结构内部热源的情况下,高效、准确地求解得到相应工程结构内部任意一点的温度,进而能够随时掌握该工程结构的温度变化状况和热传导情况,可以有效地用于实际工程中分析材料结构由于温度变化带来的影响,时刻把控工程结构温度分布情况,以确保工程结构的安全性和可靠性。
附图说明
图1为本发明实施例中涉及的工程含热源结构热传导情况求解方法的流程图;
图2为本发明实施例中涉及的模型的结构示意图;
图3为本发明实施例中涉及的本方案计算所得温度与精确解对比图;
图4为本发明实施例中涉及的工程含热源结构热传导情况求解装置的结构示意图。
具体实施方式
以下结合附图对本发明涉及的工程含热源结构热传导情况求解方法以及装置的具体实施方案进行详细地说明。
<实施例>
如图1所示,本实施例所提供的工程含热源结构热传导情况求解方法,具体包括以下步骤:
步骤1:利用温度探测器探测待评估含热源结构的外部边界温度值。本实施例中为验证准确性及有效性,选取一个具有解析解的模型,其示意图见图2,该模型为一个混凝土正方体结构,设其热传导系数为k=0.01W/(m·K),正方体的边长为1m,其温度边界条件为上下表面温度为0,其他表面进行了绝热处理。该结构自身产热热源具体表达式为Q(x)=e-x
温度解析解为:θ=[(1-e-x)-x(1-e-L)]/k。
步骤2:在结构温度求解器中,利用建模软件建立待求解模型,输入模型的材料参数、网格化分数、网格类型以及边界条件等信息,然后输出模型数据信息文件。
步骤3:建立热传导控制方程
含热源结构的热传导控制方程式为:
Figure BDA0002599741860000081
式中,k为热传导系数;θ为温度;B(x)为含热源结构热源项;x为研究域内的点;i,j,k为取值1-3的整数。
步骤4:建立温度积分方程
位移边界方程为:
Figure BDA0002599741860000082
式中,Γ为问题研究域Ω的边界;
Figure BDA0002599741860000083
为正则化的温度值,具体表达为
Figure BDA0002599741860000084
q为热通量;x和y代表问题域Ω内的源点和场点u*(x,y)和q*(x,y)为基本解,可写为:
Figure BDA0002599741860000085
Figure BDA0002599741860000086
式中,r代表源点和场点之间的距离;n为边界Γ的单位外法线向量。
步骤5:将域积分转化为边界积分
对于公式(2)中的域积分
D1=∫Ωu*(x,y)B(y)dΩ(y) (5)
基于直线积分法,利用如下公式转化为等效边界积分:
Figure BDA0002599741860000087
式中,nl1是边界Γ上外法向量n在yl1-axis方向上的元素,为了简化问题,所有起始点可以定义在同一个直线{yl1=m}上,a=a(yl2,yl3)=m为任意常数值,l1和l2分别取2和3中不同值,把边界离散为N个单元。此外
Figure BDA0002599741860000088
式中,
Figure BDA0002599741860000091
代表第i个边界积分点,t{t,y2,y3}∈BM代表积分点的个数,wi
Figure BDA0002599741860000094
分别代表第i个边界积分点在yl1方向的权重和外法向量,k(x,y)是定义在域内,对任意x,y∈Ω0,x≠y都连续的微奇异核函数,并且对于x,y∈Ω0,x≠y,都存在一个常数C1使得
|k(x,y)|≤C1r,0≤β<3 (8)
步骤6:建立遍及积分方程的求解矩阵
经过离散,除去域积分的位移边界积分方程可写为如下形式
Figure BDA0002599741860000092
其中,Γe为e个单元的边界;对于光滑边界,Ci=0.5;对于内部点,Ci=0;Ne为离散单元数;并且
Figure BDA0002599741860000093
将cij(P)整合到公式右边的第二项中去即可的以下矩阵方程:
Hu=Gt (11)
将已知的边界条件代入上式,即可得
Ax=y (12)
然后,利用高斯消去求解法即可得到模型的边界和内部点的温度值,进而可以掌握工程结构内部的温度情况,并结合设计标准进行对比以确保该工程的温度处在可控范围之内。
请见图4,图4为计算模型的温度在x方向的变化,可见计算结果与解析解高精度吻合。
基于以上计算得到数据,可以得到工程结构内部的稳温度分布情况,然后根据结构设计标准的温度耐受范围,判断是否需要对该结构采取相应的温度控制措施以保证工程结构的稳定性和安全性。
另外,以上过程还可以采用工程含热源结构热传导情况求解装置来实现,该装置包括:热传导求解器、图像生成器、判断预警器、输入显示器以及控制器。
热传导求解器能够执行上述工程含热源结构热传导情况求解方法,得到待求解含热源工程结构的温度变化和热传导数据。热传导求解器具体包括探测单元、建模单元、热传导方程建立单元、温度积分方程建立单元、转化单元、求解矩阵建立单元、求解单元。
探测单元通过温度探测器探测待求解含热源工程结构的外部边界温度值。
建模单元与探测单元通信相连,在结构温度热传导求解器中,建立待求解模型,并基于待求解含热源工程结构的材料参数、网格化分数、网格类型以及温度探测器探测得到的边界温度值信息,输出模型数据信息文件。
热传导方程建立单元与建模单元通信相连,建立热传导控制方程;含热源结构的热传导控制方程式为:
Figure BDA0002599741860000101
式中,k为热传导系数;θ为温度;B(x)为含热源结构热源项;x为研究域内的点;i,j,k为取值1-3的整数;
温度积分方程建立单元与热传导方程建立单元通信相连,建立温度积分方程;
Figure BDA0002599741860000102
式中,Γ为问题研究域Ω的边界;
Figure BDA0002599741860000103
为正则化的温度值,具体表达为
Figure BDA0002599741860000104
q为热通量;x和y代表问题域Ω内的源点和场点,u*(x,y)和q*(x,y)为基本解,具体为:
Figure BDA0002599741860000105
Figure BDA0002599741860000106
式中,r代表源点和场点之间的距离;n为边界Γ的单位外法线向量;
转化单元,与温度积分方程建立单元通信相连,将域积分转化为边界积分;对于公式(2)中的域积分:
D1=∫Ωu*(x,y)B(y)dΩ(y) (5)
基于直线积分法,利用如下公式转化为等效边界积分:
Figure BDA0002599741860000111
式中,nl1是边界Γ上外法向量n在yl1-axis方向上的元素,所有起始点定义在同一个直线{yl1=m}上,a=a(yl2,yl3)=m为任意常数值,l1和l2分别取2和3中不同值,把边界离散为N个单元;此外
Figure BDA0002599741860000112
式中,
Figure BDA0002599741860000113
代表第i个边界积分点,t{t,y2,y3}∈BM代表积分点的个数,wi
Figure BDA0002599741860000114
分别代表第i个边界积分点在yl1方向的权重和外法向量,k(x,y)是定义在域内,对任意x,y∈Ω0,x≠y都连续的微奇异核函数,并且对于x,y∈Ω0,x≠y,都存在一个常数C1使得
|k(x,y)|≤C1r,0≤β<3 (8)
求解矩阵建立单元与转化单元通信相连,建立遍及积分方程的求解矩阵;
Figure BDA0002599741860000115
式中,Γe为e个单元的边界,;对于光滑边界,Ci=0.5;对于内部点,Ci=0;Ne为离散单元数;并且
Figure BDA0002599741860000116
进一步可得以下矩阵方程:
Hu=Gt (11)
将已知的边界条件代入上式,即可得
Ax=y (12)
求解单元与求解矩阵建立单元通信相连,基于求解矩阵,求解得到模型边界和内部各点的温度值,进而得到工程结构内部的温度变化和热传导情况。
图像生成器与热传导求解器通信相连,根据热传导求解器建立的待求解模型生成含热源工程结构的三维模型结构图,并根据热传导求解器求解得到的数据生成反映温度高低和变化情况的热图像。
判断预警器与热传导求解器通信相连,存储有温度阈值,并基于工程结构内部的温度变化和热传导情况对该工程结构各部位温度是否超出温度阈值进行判断,并在判断为是的情况下,发出预警信息。
输入显示器与热传导求解器、图像生成器、判断预警器均通信相连,让操作者输入操作指令,并显示相应的信息。例如,输入显示器让操作者输入待求解工程结构的材料参数、网格化分数、网格类型以及边界条件等信息,并对求解结果进行显示。具体地,输入显示器对能够图像生成器生成的图像进行显示,并根据操作指令对热传导求解器求解得到的具体数值进行显示;进一步,输入显示器还能将预警信息显示在三维模型结构图的相应位置处。
控制器与热传导求解器、图像生成器、判断预警器、输入显示器均通信相连,控制它们的运行。
以上实施例仅仅是对本发明技术方案所做的举例说明。本发明所涉及的工程含热源结构热传导情况求解方法以及装置并不仅仅限定于在以上实施例中所描述的内容,而是以权利要求所限定的范围为准。本发明所属领域技术人员在该实施例的基础上所做的任何修改或补充或等效替换,都在本发明的权利要求所要求保护的范围内。

Claims (5)

1.一种工程含热源结构热传导情况求解方法,其特征在于,包括如下步骤:
步骤1:利用温度探测器探测待求解含热源工程结构的外部边界温度值;
步骤2:在结构温度热传导求解器中,利用建模软件建立待求解模型,输入待求解含热源工程结构的材料参数、网格化分数、网格类型以及温度探测器探测得到的边界温度值信息,然后输出模型数据信息文件;
步骤3:建立热传导控制方程
含热源结构的热传导控制方程式为:
Figure FDA0003796802520000011
式中,k为热传导系数;θ为温度;B(x)为含热源结构热源项;x为研究域内的源点;s为取值1-3的整数;
步骤4:建立温度积分方程
位移边界方程为:
Figure FDA0003796802520000012
式中,Γ为研究域Ω的边界;
Figure FDA0003796802520000013
为正则化的温度值,具体表达为
Figure FDA0003796802520000014
q为热通量;x和y代表研究域Ω内的源点和场点;u*(x,y)和q*(x,y)为基本解,具体为:
Figure FDA0003796802520000015
Figure FDA0003796802520000016
式中,r代表源点和场点之间的距离;n为边界Γ的单位外法线向量;
步骤5:将域积分转化为边界积分
对于公式(2)中的域积分
D1=∫Ωu*(x,y)B(y)dΩ(y) (5)
基于直线积分法,利用如下公式转化为等效边界积分:
Figure FDA0003796802520000017
式中,nl1是边界Γ上外法向量n在yl1-axis方向上的元素,所有起始点定义在同一个直线{yl1=m}上,a=a(yl2,yl3)=m为任意常数值,l1和l2分别取2和3中不同值,把边界离散为N个单元;此外
Figure FDA0003796802520000021
式中,
Figure FDA0003796802520000022
代表第i个边界积分点,t{t,y2,y3}∈B,M代表积分点的个数,wi
Figure FDA0003796802520000023
分别代表第i个边界积分点在yl1方向的权重和外法向量,k(x,y)是定义在域内,对任意x,y∈Ω0,x≠y都连续的微奇异核函数,并且对于x,y∈Ω0,x≠y,都存在一个常数C1使得
|k(x,y)|≤C1r,0≤β<3 (8)
步骤6:建立遍及积分方程的求解矩阵
经过离散,除去域积分的位移边界积分方程可写为如下形式
Figure FDA0003796802520000024
式中,Γe为e个单元的边界;对于光滑边界,Ci=0.5;对于内部点,Ci=0;Ne为离散单元数;并且
Figure FDA0003796802520000025
进一步可得以下矩阵方程:
Hu=GQ (11)
将已知的边界条件代入上式,即可得
AX=Y (12)
步骤7:基于步骤6的求解矩阵,得到模型边界和内部各点的温度值,进而得到工程结构内部的温度变化和热传导情况。
2.一种工程含热源结构热传导情况求解装置,其特征在于,包括:
热传导求解器,根据权利要求1所述的工程含热源结构热传导情况求解方法,得到待求解含热源工程结构的温度变化和热传导数据;
输入显示器,与所述热传导求解器通信相连,让操作者输入操作指令,并显示相应的信息;以及
控制器,与所述热传导求解器、所述输入显示器通信相连,控制所述热传导求解器和所述输入显示器的运行。
3.根据权利要求2所述的工程含热源结构热传导情况求解装置,其特征在于:
其中,所述热传导求解器包括:
探测单元,通过温度探测器探测待求解含热源工程结构的外部边界温度值;
建模单元,与所述探测单元通信相连,在结构温度热传导求解器中,建立待求解模型,并基于待求解含热源工程结构的材料参数、网格化分数、网格类型以及温度探测器探测得到的边界温度值信息,输出模型数据信息文件;
热传导方程建立单元,与所述建模单元通信相连,建立热传导控制方程;含热源结构的热传导控制方程式为:
Figure FDA0003796802520000031
式中,k为热传导系数;θ为温度;B(x)为含热源结构热源项;x为研究域内的源点;s为取值1-3的整数;
温度积分方程建立单元,与所述热传导方程建立单元通信相连,建立温度积分方程;
Figure FDA0003796802520000032
式中,Γ为研究域Ω的边界;
Figure FDA0003796802520000033
为正则化的温度值,具体表达为
Figure FDA0003796802520000034
q为热通量;x和y代表研究域Ω内的源点和场点,u*(x,y)和q*(x,y)为基本解,具体为:
Figure FDA0003796802520000035
Figure FDA0003796802520000036
式中,r代表源点和场点之间的距离;n为边界Γ的单位外法线向量;
转化单元,与所述温度积分方程建立单元通信相连,将域积分转化为边界积分;对于公式(2)中的域积分:
D1=∫Ωu*(x,y)B(y)dΩ(y) (5)
基于直线积分法,利用如下公式转化为等效边界积分:
Figure FDA0003796802520000041
式中,nl1是边界Γ上外法向量n在yl1-axis方向上的元素,所有起始点定义在同一个直线{yl1=m}上,a=a(yl2,yl3)=m为任意常数值,l1和l2分别取2和3中不同值,把边界离散为N个单元;此外
Figure FDA0003796802520000042
式中,
Figure FDA0003796802520000043
代表第i个边界积分点,t{t,y2,y3}∈B,M代表积分点的个数,wi
Figure FDA0003796802520000044
分别代表第i个边界积分点在yl1方向的权重和外法向量,k(x,y)是定义在域内,对任意x,y∈Ω0,x≠y都连续的微奇异核函数,并且对于x,y∈Ω0,x≠y,都存在一个常数C1使得
|k(x,y)|≤C1r,0≤β<3 (8)
求解矩阵建立单元,与所述转化单元通信相连,建立遍及积分方程的求解矩阵;
Figure FDA0003796802520000045
式中,Γe为e个单元的边界,;对于光滑边界,Ci=0.5;对于内部点,Ci=0;Ne为离散单元数;并且
Figure FDA0003796802520000046
进一步可得以下矩阵方程:
Hu=GQ (11)
将已知的边界条件代入上式,即可得
AX=Y (12)
求解单元,与所述求解矩阵建立单元通信相连,基于求解矩阵,求解得到模型边界和内部各点的温度值,进而得到工程结构内部的温度变化和热传导情况。
4.根据权利要求2所述的工程含热源结构热传导情况求解装置,其特征在于,还包括:
图像生成器,与所述热传导求解器、所述输入显示器、所述控制器均通信相连,根据所述热传导求解器建立的待求解模型生成所述含热源工程结构的三维模型结构图,并根据所述热传导求解器求解得到的数据生成反映温度高低和变化情况的热图像;
其中,所述输入显示器对所述图像生成器生成的图像进行显示,并根据操作指令对所述热传导求解器求解得到的具体数值进行显示。
5.根据权利要求4所述的工程含热源结构热传导情况求解装置,其特征在于,还包括:
判断预警器,与所述热传导求解器、所述输入显示器、所述控制器均通信相连,存储有温度阈值,并基于工程结构内部的温度变化和热传导情况对该工程结构各部位温度是否超出所述温度阈值进行判断,并在判断为是的情况下,发出预警信息;
其中,所述输入显示器还将所述预警信息显示在所述三维模型结构图的相应位置处。
CN202010720362.3A 2020-07-24 2020-07-24 工程含热源结构热传导情况求解方法以及装置 Active CN111898263B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010720362.3A CN111898263B (zh) 2020-07-24 2020-07-24 工程含热源结构热传导情况求解方法以及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010720362.3A CN111898263B (zh) 2020-07-24 2020-07-24 工程含热源结构热传导情况求解方法以及装置

Publications (2)

Publication Number Publication Date
CN111898263A CN111898263A (zh) 2020-11-06
CN111898263B true CN111898263B (zh) 2022-10-04

Family

ID=73190574

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010720362.3A Active CN111898263B (zh) 2020-07-24 2020-07-24 工程含热源结构热传导情况求解方法以及装置

Country Status (1)

Country Link
CN (1) CN111898263B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116796569B (zh) * 2023-07-10 2024-03-29 锦州佳润农业科技有限公司 一种温室环境调控系统及方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6513000B1 (en) * 1998-03-11 2003-01-28 Nec Corporation Simulation method of wiring temperature rise
CN104050307A (zh) * 2013-12-23 2014-09-17 南京理工大学 基于时域谱元法的微波无源电路电磁热一体化分析方法
CN110705057A (zh) * 2019-09-19 2020-01-17 武汉大学 各向同性固体材料的静态热弹性问题求解方法以及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6513000B1 (en) * 1998-03-11 2003-01-28 Nec Corporation Simulation method of wiring temperature rise
CN104050307A (zh) * 2013-12-23 2014-09-17 南京理工大学 基于时域谱元法的微波无源电路电磁热一体化分析方法
CN110705057A (zh) * 2019-09-19 2020-01-17 武汉大学 各向同性固体材料的静态热弹性问题求解方法以及装置

Also Published As

Publication number Publication date
CN111898263A (zh) 2020-11-06

Similar Documents

Publication Publication Date Title
Hartmann et al. An optimal order interior penalty discontinuous Galerkin discretization of the compressible Navier–Stokes equations
Gomez-Iradi et al. Development and validation of a CFD technique for the aerodynamic analysis of HAWT
Constantinescu et al. Numerical investigation of flow past a prolate spheroid
Feng et al. Fast and efficient analysis of transient nonlinear heat conduction problems using combined approximations (CA) method
Medic et al. Prediction of transition and losses in compressor cascades using large-eddy simulation
Alibrandi et al. First-order reliability method for structural reliability analysis in the presence of random and interval variables
Zhang et al. Numerical evaluation of single fiber motion for short-fiber-reinforced composite materials processing
Weng et al. The complex variable reproducing kernel particle method for two-dimensional inverse heat conduction problems
He et al. An Element-free Galerkin (EFG) scaled boundary method
Zhu et al. An improved direct-forcing immersed boundary method for fluid-structure interaction simulations
CN111898263B (zh) 工程含热源结构热传导情况求解方法以及装置
Hassanzadeh-Aghdam et al. Thermal expanding behavior of carbon nanotube-shape memory polymer nanocomposites
Galán-Marín et al. A new multivalued neural network for isomorphism identification of kinematic chains
Wu et al. A novel classification method to random samples for efficient reliability sensitivity analysis
Patil et al. Robust deep learning for emulating turbulent viscosities
Zhou et al. Fast transonic flow prediction enables efficient aerodynamic design
Acton et al. Effect of volume element geometry on convergence to a representative volume
Almuzaiqer et al. Effect of the aspect ratio and tilt angle on the free convection heat transfer coefficient inside Al2O3–water-filled square cuboid enclosures
Wang et al. Finite element analysis of the thermal characteristics and parametric study of steady rolling tires
Bidar et al. Sensor placement for data assimilation of turbulence models using eigenspace perturbations
CN112213638A (zh) 一种电机线圈的发热状态监测方法、相关设备和介质
Zhu et al. A multiple regression convolutional neural network for estimating multi-parameters based on overall data in the inverse heat transfer problem
Zeng et al. An explicit phase field material point method for modeling dynamic fracture problems
Kumar et al. Finite element analysis using nonconforming mesh
Li et al. A local meshless collocation method for solving certain inverse problems

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant