A Boundary Element Method for Steady-State Heat Transfer Problems Based on a Novel Type of Interpolation Elements
-
摘要: 为了提高边界元法在求解稳态热问题时的计算精度,通过使用一种新型单元插值方法(称为扩展单元插值法)，实现对稳态传热问题的求解。扩展单元是在传统不连续单元的边界配置虚拟节点，把原非连续单元变成高阶的连续单元，并将其作为新型的插值单元。利用虚拟节点和内部源节点构造出的插值函数，可以精确插值边界上的连续和不连续物理场，插值精度要比原始不连续单元高两阶。另外，边界积分方程只在传统的不连续单元的内部节点处建立，只包含内部源节点的自由度，而虚拟节点的自由度可通过与内部源节点之间的关系消除掉，因此最终系统方程的求解规模不会增加。这种新型的插值单元继承了传统连续和不连续单元的优点，克服了它们的缺点。数值结果表明，此种单元插值方法用于求解稳态传热问题时可获得较高的计算精度和收敛性。Abstract: To improve the computational accuracy of the boundary element method (BEM) for solving the steady-state heat problems, a new method was studied with a new element interpolation method (called expansion element interpolation method). The novel expansion interpolation element was obtained through configuration of virtual nodes at the boundary of the traditional discontinuous element, to transform the original discontinuous element into a higher-order continuous element. The interpolation function constructed with the virtual node and the internal source node can accurately work in the continuous and discontinuous physical fields on the boundary, and the interpolation accuracy increases by 2 orders compared with that of the original discontinuous element. In addition, the boundary integral equation is only valid at the internal source nodes of the traditional discontinuous element and contains only the degrees of freedom of the source nodes, while the degrees of freedom of the virtual node can be eliminated through the relationship between the virtual node and the internal source node, so the size of the final linear system equations will not increase. This new interpolation element inherits the advantages of the traditional continuous and discontinuous elements and overcomes their disadvantages. The numerical results show that, the proposed method helps solve the steady-state heat transfer problem with high computational accuracy and convergence.
-
图 2 新型单元用于插值已知边界变量
Figure 2. The novel elements used for interpolation of the known boundary variables
图 3 新型单元用于插值未知边界变量
Figure 3. The novel elements used for interpolation of the unknown boundary variables
图 4 新型单元用于插值非连续边界变量
Figure 4. The novel elements used for interpolation of the discontinuous boundary variables
图 5 带孔平板的几何模型（单位：mm）
Figure 5. The geometric model for the plate with a hole (unit: mm)
图 7 内孔边界上热流量q的数值结果
Figure 7. The numerical results of thermal flux on the boundary of the internal hole
图 8 底边和侧边界上热流量q的相对误差和收敛性图
Figure 8. The relative errors and convergence rates of q on the bottom and side edges
图 10 复合支架所施加的混合边界条件
Figure 10. The mixed boundary conditions applied for the composite bracket
图 11 左侧边界和圆弧边界（R＝0.5 m）上温度的计算结果
Figure 11. The computational results of temperature on left and arc (R＝0.5 m) boundaries
图 12 本文方法计算得到的温度云图
注 为了解释图中的颜色，读者可以参考本文的电子网页版本。
Figure 12. The temperature nephogram obtained with the proposed method
-
[1] 姚振汉, 王海涛. 边界元法[M]. 北京: 高等教育出版社, 2010.YAO Zhenhan, WANG Haitao. Boundary Element Method[M]. Beijing: Higher Education Press, 2010. （in Chinese） [2] BREBBIA C A. Boundary Element Methods[M]. Berlin: Springer-Verlag, 1981. [3] ZHANG Y M, GONG Y P, GAO X W. Calculation of 2D nearly singular integrals over high-order geometry elements using the sinh transformation[J]. Engineering Analysis With Boundary Elements, 2015, 60: 144-153. doi: 10.1016/j.enganabound.2014.12.006 [4] 郭钊, 郭子涛, 易玲艳. 多裂纹问题计算分析的本征COD边界积分方程方法[J]. 应用数学和力学, 2019, 40(2): 200-209. (GOU Zhao, GOU Zitao, YI Lingyan. Analysis of multicrack problems with eigen cod boundary integral equations[J]. Applied Mathematics and Mechanics, 2019, 40(2): 200-209.(in Chinese) [5] ZHANG J M, ZHONG Y D, DONG Y Q, et al. Expanding element interpolation method for analysis of thin-walled structures[J]. Engineering Analysis With Boundary Elements, 2018, 86: 82-88. doi: 10.1016/j.enganabound.2017.10.014 [6] 董荣荣, 张超, 张耀明. 三维位势问题的梯度边界积分方程的新解法[J]. 力学学报, 2020, 52(2): 472-479. (DONG Rongrong, ZHANG Chao, ZHANG Yaoming. A new method for solving the gradient boundary integral equation for three dimensional potential problems[J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(2): 472-479.(in Chinese) doi: 10.6052/0459-1879-19-308 [7] 窦智峰, 晋玉祥, 郭新飞. IGBT瞬态短路失效分析及其有限元热电耦合模型研究[J]. 轻工学报, 2018, 33(6): 101-108. (DOU Zhifeng, JIN Yuxiang, GUO Xinfei. Research on transient short-circuit failure analysis and finite element thermoelectric coupling mode of IGBT[J]. Journal of Light Industry, 2018, 33(6): 101-108.(in Chinese) doi: 10.3969/j.issn.2096-1553.2018.06.012 [8] 王勖成. 有限单元法[M]. 北京: 清华大学出版社, 2003.WANG Xucheng. Finite Element Method[M]. Beijing: Tsinghua University Press, 2003. （in Chinese） [9] 李聪, 牛忠荣, 胡宗军, 等. 三维切口/裂纹结构的扩展边界元法分析[J]. 力学学报, 2020, 52(5): 1394-1408. (LI Chong, NIU Zhongrong, HU Zongjun, et al. Analysis of 3-D notched/cracked structures by using extended boundary element method[J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(5): 1394-1408.(in Chinese) doi: 10.6052/0459-1879-20-129 [10] YAO Z H. A new type of high-accuracy BEM and local stress analysis of real beam, plate and shell structures[J]. Engineering Analysis With Boundary Elements, 2016, 65: 1-17. doi: 10.1016/j.enganabound.2015.12.011 [11] 周枫林, 谢贵重, 张见明, 等. 角度-距离复合变换法消除边界积分方程近奇异性[J]. 应用数学和力学, 2020, 41(5): 530-540. (ZHOU Fenglin, XIE Guizhong, ZHANG Jianming, et al. Near-singularity cancellation with the angle-distance transformation method for boundary integral equations[J]. Applied Mathematics and Mechanics, 2020, 41(5): 530-540.(in Chinese) [12] GU Y, SUN L L. Electroelastic analysis of 2D ultra-thin layered piezoelectric films by an advanced boundary element method[J]. International Journal for Numerical Methods in Engineering, 2021, 122(6): 1-19. [13] 周焕林, 严俊, 余波, 等. 识别含热源瞬态热传导问题的热扩散系数[J]. 应用数学和力学, 2018, 39(2): 160-169. (ZHOU Huanlin, YAN Jun, YU Bo, et al. Identification of thermal diffusion coefficients for transient heat conduction problems with heat sources[J]. Applied Mathematics and Mechanics, 2018, 39(2): 160-169.(in Chinese) [14] 高效伟, 曾文浩, 崔苗. 等参管单元及其在热传导问题边界元法中的应用[J]. 计算力学学报, 2016, 33(3): 328-334. (GAO Xiaowei, ZENG Wenhao, CUI Miao. Isoparametric tube elements and their application in heat conduction BEM analysis[J]. Chinese Journal of Computational Mechanics, 2016, 33(3): 328-334.(in Chinese) [15] 张驰, 校金友. 用边界元分析功能梯度材料的稳态热传导[J]. 科学技术与工程, 2013, 13(13): 3563-3565, 3571. (ZHANG Chi, XIAO Jinyou. Boundary element method for FGM in steady-state heat conduction[J]. Science Technology and Engineering, 2013, 13(13): 3563-3565, 3571.(in Chinese) doi: 10.3969/j.issn.1671-1815.2013.13.008 [16] PARREIRA P. On the accuracy of continuous and discontinuous boundary elements[J]. Engineering Analysis, 1988, 5(4): 205-211. doi: 10.1016/0264-682X(88)90018-4 [17] MUKHOPADHYAY S, MAJUMDAR N. A study of three-dimensional edge and corner problems using the neBEM solver[J]. Engineering Analysis With Boundary Elements, 2009, 33(2): 105-119. doi: 10.1016/j.enganabound.2008.06.003 [18] HAMZAH K B, LONG N M A N, SENU N, et al. Stress intensity factor for multiple cracks in bonded dissimilar materials using hypersingular integral equations[J]. Applied Mathematical Modelling, 2019, 73: 95-108. doi: 10.1016/j.apm.2019.04.002 -
下载:
点击查看大图
计量
- 文章访问数: 576
- HTML全文浏览量: 224
- PDF下载量: 71
- 被引次数: 0
