凸域上温度分布的有理插值近似.pdf

计算物理 第z擘鬻期年 11月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS V磊01．22N0v．j2005 [文It编号】1001．246x{2005)06．0507—08 凸域上温度分布的有理插值近似 王兆清1’2， 冯 伟1 (1．上海市应用数学和力学研究所，上海大学，上海200072；2．山东建筑工程学院工程力学研究所，山东济南250101) [摘要】 构造出凸多边形上的一个新型有理插值，并将其推广到含有边节点的情况．给出了多边形上有理函 数插值形函数的计算表达式．利用构造出的多边形有理插值，采用凸多边形逼近任意凸域，通过区域边界温度 离散值，插值近似区域内的温度场分布．给出了圆形区域和正方形区域温度分布插值近似的算例．与现有的插 值方法相比，具有程序实现简单、精度高等优点． [关键词】 温度分布；有理插值；凸区域；多边形单元；形函数；边节点 [中图分类号】0241．3 [文献标识码] A O 引言 许多工程问题的研究需要知道求解域内的温度场分布，如力学中的热应力计算等…．稳态温度场分布的 界条件，难以得到温度分布的解析表达式．并且在实际工作中，我们能知道的一般是区域边界温度分布的离 散数据，如何通过这些边界上的离散数据近似得到区域内的温度分布，是工程技术人员所关心的问题． 近似区域内温度分布的传统方法是有限元法．将求解区域划分成三角形或四边形单元，在每个单元上用 多项式插值单元节点温度值，得到单元上的近似温度分布，进而得到整个区域的近似温度分布．有限元法需 要在区域内部布置节点，所得到的温度场的梯度是不连续的． 采用多边形划分区域，可以减少单元数量，进而改善区域内温度场梯度的不连续性．对于凸区域可以采 用一个多边形单元插值，如对圆形区域采用一个圆内接正20边形单元插值，由此得到的近似温度场在区域 内梯度是连续的．对于边数大于4的多边形单元，难以构造满足协调性要求的多项式形式插值函数b1．既然 多项式形式的插值函数难以构造，设想能否采用有理函数构造多边形单元的插值函数．因为有理函数也是一 种比较简单的函数类，并且在有些情况下，有理函数的近似精度比多项式函数更好H1． (Projective 构造出含有边节点的四边形单元上的无理函数插值，并应用于近似正方形区域的温度分布口1．范亚玲等采用 正多边形等参变换法构造出任意多边形单元上的有理函数插值u…． 本文采用几何度量的方法直接构造任意凸多边形上的有理函数插值，在插值过程中既考虑插值点与插 值节点的距离，又照顾到插值节点的分布情况，并应用于凸区域的温度分布插值近似．所构造的有理函数插 [收稿日期】2004—06—02；[修回日期】2004—10—19 【作者简介】王兆清(1965一)，男，山东荣成，副教授，博士，从事计算力学及其应用软件、工程中的数值方法方面的研究，山东省济南市临 港开发区泉港路山东建筑工程学院科技处工程力学研究所250101． 万方数据万方数据 计 算 物 理 第22卷 式形式是相同的．并且可以推广到含有边节点的任意多边形． 1 凸多边形上有理函数插值的构造 (i=1，2，…，11,)．函数“(髫，y)称作点(算，，Y。)，(菇：，Y：)，…，(石。，Y。)的插值函数． 在不同的函数类中可以构造出不同的插值函数．例如在多项式函数类中构造的为多项式插值函数；在有 理函数类中构造的即为有理函数插值．插值函数一般采用形函数(数学上称作基函数)表达， M(石，Y)=：≯；(戈，Y)“；， (1) 膏 其中声；(戈，Y)为对应节点(戈i，Y；)的形函数．有了节点处的形函数即可得到插值函数． 对顶点为(Ⅳ。，Y。)，(z：，Y：)，…，(戈。，Y。)的多边形，为推导 方便起见，记z；=(并i，Y；)，x=(菇，Y)．给定多边形内的一个点 P，记P与多边形第i个顶点Pi的距离为h；(z)，依次作线段 PP；的垂直平分线，得到一个凸多边形，如图1所