1、第 44 卷第 2 期2023 年 2 月哈 尔 滨 工 程 大 学 学 报Journal of Harbin Engineering UniversityVol.44.2Feb.2023二维功能梯度材料热传导格点型有限体积法研究龚京风,徐宗著,宣领宽,江致远,郑文利(武汉科技大学 汽车与交通工程学院,湖北 武汉 430065)摘 要:为了分析二维功能梯度材料热传导问题,发展格点型有限体积法,本文利用交错网格技术处理任意空间分布的材料物性,并使用 C+语言开发数值求解程序。利用该程序分析不同网格类型下功能梯度材料圆环和不同边界条件下功能梯度材料板的热传导问题。结果表明:发展的格点型有限体积法对
2、于线性三角形单元及双线性四边形单元的适用性良好,不同边界条件下格点型有限体积法解与解析解的误差不超过0.12%,验证了求解程序的正确性。与有限元及解析解结果对比发现,当结构角点处存在剧烈温度变化且网格较为粗糙时,有限元结果出现数值振荡的现象,而格点型有限体积法结果仍然合理。分析表明:对于复杂的功能梯度材料热传导问题,格点型有限体积法比有限元方法具有更好的数值稳定性。关键词:交错网格;格点型有限体积法;二维功能梯度材料;热传导;数值振荡;C+程序语言;误差分析;CAE 软件DOI:10.11990/jheu.202105015网络出版地址:https:/ 文献标志码:A 文章编号:1006-70
3、43(2023)02-0321-08Investigation on cell-vertex finite volume method for the heat conduction of two-dimensional functionally graded materialsGONG Jingfeng,XU Zongzhu,XUAN Lingkuan,JIANG Zhiyuan,ZHENG Wenli(School of Automobile and Traffic Engineering,Wuhan University of Science and Technology,Wuhan 4
4、30065,China)Abstract:To analyze the heat conduction problem of two-dimensional functionally graded materials(FGMs),the staggered cell-vertex finite volume method(CV-FVM)is developed in this research.Staggered grid technology is ap-plied to process the material properties of arbitrary spatial distrib
5、ution,and the C+language is used to develop the numerical solution program.This program is used to analyze the heat conduction of FGM rings under different grid types and FGM plates under different boundary conditions.The results show that the developed CV-FVM method has good applicability for linea
6、r triangular elements and bilinear quadrilateral elements.The error between the CV-FVM and analytical solutions under different boundary conditions does not exceed 0.12%,which verifies the correctness of the solution procedure.Comparing the finite element and analytical results shows that when sever
7、e temperature changes occur at the corners of a structure with a coarse grid,the finite element results exhibit numerical oscillations,but the CV-FVM results are still reasonable.The analyses indicate that compared with the complicated finite element method,CV-FVM has better numerical stability for
8、complicated heat conduction analysis of FGMs.Keywords:staggered grid;cell-vertex finite volume method;functionally graded materials;heat conduction;nu-merical oscillation;C+programming language;deviation analysis;CAE software收稿日期:2021-05-01.网络出版日期:2022-10-28.基金项目:国家自然科学基金青年基金项目(51909197);湖北省自然科学基金项目
9、(2019CFB136);武汉科技大学预研基金项目(GF201911).作者简介:龚京风,女,副教授,硕士生导师;宣领宽,男,高级工程师.通信作者:宣领宽,E-mail:nsmpf_wust .功能梯度材料(functionally graded materials,FGM)被广泛应用于航空航天、核反应堆、内燃机等领域。由于 FGM 常被应用于高温高压的恶劣工作环境,FGM 结构内的温度场、应力场变化剧烈,其热力性能得到了广泛关注。刘文光等1研究了热流密度、陶瓷体积分数指数和热环境等参数对功能梯度圆柱壳热传导的影响;Fu 等2分析了在弹性基础上,考虑非线性热环境的多孔功能梯度材料圆柱壳的热声响
10、应;许新等3基于 Eluer-Bernoulli 梁理论和单向耦合的热传导理论,研究了 FGM 微梁的热弹性阻尼;梅靖等4研究了考虑厚度方向稳态温度分布的石墨烯增强功能梯度复合材料板条热弹性问题。Steinberg 等5提出为承受复杂的温度场,要求材料参数性能在 2 个甚至 3 个方向变化,因此,研究哈 尔 滨 工 程 大 学 学 报第 44 卷二维 FGM 热传导问题是有必要的。Rahul 等6在一阶剪切变形理论基础上,对二维温度变化作用下的双向功能梯度圆板进行了自由振动分析。Tang等7研究了具有横向和纵向位移耦合的双向功能梯度梁的非线性湿热动力学。这些研究都考虑多方向功能梯度材料的温度场
11、影响,研究热环境下的二维 FGM 结构性能,准确预测其温度场是前提。许杨健等8推导了复杂边界条件下的二维FGM 板的稳态温度场解析解;蒋科9在离散域采用数值法,非离散域采用解析法的混合方法,研究了FGM 结构的热传导特性。解析方法虽然能获得精确解,但难以应用于复杂工程问题,而数值模拟方法具有更广泛的适用性。Sladek 等10采用无网格局部 Petrov-Galerkin(meshless local Petrov-Galerkin,MLPG)方法分析了非均匀各向异性 FGM 三维轴对称结构的瞬态热传导问题;Ching 等11采用 MLPG方法求解了功能梯度梁的瞬态热弹性变形。MLPG方法计算
12、量较大,处理结构复杂网格量大的热机械问题存在一定困难。仝国军等12采用有限元方法(finite element methed,FEM)研究了非均匀温度场下变物性二维 FGM 板的瞬态热传导分布问题。针对 FGM 热传导问题,Gong 等13-14提出非结构时域有限体积法(finite volume methed,FVM),该方法将四边形单元处理为双线性单元,有效避免由于将四边形单元作为线性单元处理的数值振荡问题;随后进一步发展该方法,提出交错格点型有限体积法(cell-vetex finite volume methed,CV-FVM),用于研究活塞、功能梯度多孔材料等热传导问题15。当前,基
13、于 FEM 的数值计算方法在固体热传导分析领域占据主导地位。考虑到 FVM 在流体力学里的广泛应用及其在固体力学中的成功应用,为了发展统一的多物理场数值求解方法,拟基于 FVM 进一步探究其在 FGM 热传导问题中的应用。为了研究二维 FGM 材料热传导问题,本文进一步发展交错 CV-FVM,将材料物性作为空间坐标的函数,考虑全域的物性参数变化。基于线性三角形单元和双线性四边形单元,建立二维 FGM 热传导数值求解方法,并验证其正确性,讨论其适用性;同时,在处理复杂 FGM 热传导问题,考虑大温度梯度区域难以预知情况时,探究 CV-FVM 在不加密该区域网格情况下的数值稳定性。本文通过研究求解
14、二维FGM 问题的 CV-FVM 解法,为开发拥有自主知识产权的 CAE 分析软件奠定基础。1 热传导数值方法1.1 控制方程 在控制体 M 上建立导热方程,控制体体积为 V,边界面积为 S。根据能量守恒定律建立热传导方程:V(cT)tdV=S-qndS+VQdV(1)式中:、c、T 分别是密度、比热容和温度;q 为热流密度;n 为控制体边界面的外法线矢量;Q 为热源在单位体积内产生的热量。由傅里叶定律可知-q=kT,则可将式(1)写成:V(cT)tdV=SkTndS+VQdV(2)式中 k 为导热系数,考虑各项同性问题,k 为标量。1.2 数值离散1.2.1 控制体的建立 对于二维问题,CV
15、-FVM 依次将围绕节点的相邻单元中心、边中点连接,构造控制体。采用三角形单元或四边形单元离散计算域。图 1 为局部网格示意图,空心点表示单元中心或边中点,实心点为单元节点。图 1 二维 CV-FVM 离散单元示意Fig.1Schematic diagram of two-dimensional CV-FVM discrete unit以节点 1 为例,构造的控制体为图 1 中虚线围成的区域 O1-L2-O2-L3-O3-L4-O4-L1-O1,该区域即围绕节点 1 的控制体 M,其中,L1-O1-L2表示节点 1 在相邻单元 1 内的控制体边界面,节点 27 皆为节点1 的相邻节点。采用交错
16、网格的思想,将待解温度定义在网格节点,并假设温度在控制体内均匀分布;将材料参数、体积热源定义在单元中心,假设其在网格单元内均匀分布。由于材料物性的空间分布,在控制体内,物性参数非均匀,从而将物性参数的空间分布引入离散方程。1.2.2 导热方程的离散 为了方便方程的离散,将式(2)改写为:VcTtdV=SkTxnx+Tyny()dS+VQdV(3)式中 nx、ny分别表示控制体边界面的法向矢量在 x、223第 2 期龚京风,等:二维功能梯度材料热传导格点型有限体积法研究y 方向的分量。方程(3)右端第 1 项借用 FEM 的思想,采用形函数对该项进行离散:SkTxnx+Tyny()dS=ni=1kimj=1TijSiNijxnxdS+mj=1TijSiNijynydS()(4)式中:n 表示当前节点的相邻单元总数;m 表示第 i个相邻单元的节点总数;Si表示当前节点在第 i 个相邻单元内的控制体边界面积;Si(Nij/x)nxdS 以及Si(Nij/y)nydS 分别表示第 i 个单元内第 j 个节点上的形函数对 x、y 的偏导与积分路径的乘积,该项只与网格有关,具体计算方法参考文献16