求解瞬时温度场的有限元显式算法
A Finite Element Explicit Algorithm for Solving the Temporal Temperature Fields
-
摘要: 在一公共节点为中心的各单元中,对于线性形函数,实际计算和数字实验表明,温度在单元各节点上的时间导数用它在中心节点上的时间导数表示是可取和合理的。由此可在对微分方程用有限元法进行空间离散的基础上得到单个节点温度的时间导数与其周围节点温度的关系,建立温度场的显式计算格式。它具有计算简捷的特点。用最大值原理对稳定性的分析导出了与稳式算法类似的稳定性条件。Abstract: Practical calculations and numerical experiments in this paper have shown that in elements relating to a common node it is acceptable and reasonable for derivaties of temperature with respect to time on nodes of those elements to be presented with one on common node,if linear interpolation shape function is taken.The relation between the derivative of temperature to time on a certain node and the temperature on other nodes around that node may therefore be established after discretization of the differential equation is made in space by the finite element method.Then an explicit scheme for calculating the temperature fields may be constructed.The obtained algebraic equations,being simple and the procedure being straight will be its two tangible advantages and its calculating will,therefore,be fast.The stability analysis by the maximum principle,as in the example quoted,proves that the stability condition is similar to that in implicit algorithms.
-
[1] Wilson,E.L.and R.E.Nickell,Application of the finite element method to heat conduction analysis,Nuclear Engineering and Design,4(1966),276-286. [2] Zienkiewicz,O.C.,The Finite,Element Method in Engineering Science(1971). [3] 王润富、陈和群、李克敌,在有限单元法中灾用平衡原理求解不稳定温度场,全国计算力学会议论文(1980). [4] 515科研组,瞬变温度场问题的有限元解法和最大模原理,计算数学,4.4(1982). [5] 周顺兴,解抛物型方程的Dufort-Frankal格式和有限单元法,计算数学,4(1978),53-64. [6] 朱伯芳,‘有限单元法原理与应用》,水利电力出版社(1979). [7] 刘润广,有限元法求解不稳定温度场及温度应力,全国计算力学会议论文(1980). [8] 孔徉谦,有限单元法求解不稳定导热问题中的时间差分格式及其变步长计算,工程热物理学报,3,3(1982). [9] 南京大学计算数学专业,嵘偏微分方程数值解法》,科学出版(1919). [10] Nogotov,E.F.,Application of Numerical Heat Transfer(1979). [11] [ll]黄玉霞,带第三类边界条件的热传导方程的守恒型格式和最大模原理,计算数学,5.2(1983). [12] [l2] Ozisik,M.N.,Heat Conduction(1980).
计量
- 文章访问数: 1842
- HTML全文浏览量: 73
- PDF下载量: 644
- 被引次数: 0