具有小周期参数抛物型方程一种新的双尺度有限元分析.pdfVIP

第27卷 第3期 河 南 科 学 Vo1．27 No．3 2009牟3月 HENAN SCIENCE Mar．2009 文章编号：1004—3918(2009)03—0253—05 具有小周期参数抛物型方程一种 新的双尺度有限元分析 孙艳萍 ， 王 自强 ， 曹俊英 (1．河南工程学院数理科学系，郑州 451191； 2．贵州民族学院数学与计算机科学学院，贵阳 550025) 摘 要：利用均匀化，渐进展开式和双尺度有限元方法，对具有高阶震荡系数的抛物型方程给出了一种新的半离散 有限元格式，即在时间上保持连续，空间上利用双尺度有限元进行离散．并分析了该格式的收敛性． 关键词：均匀化；双尺度有限元；复合材料；抛物型方程 中图分类号：0242．2 文献标识码：A 自从7O年代起，随着计算机软硬件的发展，出现了许多涉及具体物理问题的算法与理论．但对于小 周期复合材料的计算，若el非常小时，口()变化非常频繁，差分或有限元等数值求解时，网格剖分要求非 常细，这样最终求解线性方程组的计算量非常大．针对类似这样的问题，渐渐发展起来了均匀化和多尺度 渐进展开方法，在宏观尺度上求解均匀化解，在微观尺度上增加调整项，均匀化解和调整项构成原始解的渐 进展开形式，从而降低了计算量．文献[4]讨论二阶椭圆问题的均匀化和渐进展开，对抛物型方程仅给出了 均匀化结果．在文献 [2．3，7]之中也从不同方面对 [4]中的结果进行了很大的改进，但也只对椭圆型方程进 行了讨论．对抛物型方程多尺度渐进展开的讨论还至今没见到报道．我们在文献 [9]的基础上再给出一种 新的半离散格式的双尺度有限元格式． 平面复合材料的热传导问题可描述为 去(∽ )， ×0( u ，)：0， ，￡)∈aD×(0， ， (1) ， 0)：M。 )， ∈Q， 其中：“表示 ；表示温度分布；／是热源；。)为初始温度；。；()(，=1，2)为热传导系数，是以2为周期 的周期函数．对任意2ex28的单胞K，口()=。，(， Q=[一I，1]×[_i，1]，Q为参考单胞． = =×+：， ∈Q)，XK为K中心，考虑到不同材料问的粘结 比较紧密，可 以认为n()∈C(fD，在数学上理解为函数的光 滑化．另外，假设 (n())2x2∈ ，卢)，即存在常数Oa／3，使得 2 2 ∑77≤aij) ≤卢∑叼，V(，)≠(0，0)． (2) i=1 =l 文 中采用 Einstein求和记号，即带相同下标量的乘积表示从 1到 2求和．为了以后使用方便，记 (Q) 为 ( )的子空间，其中的函数是以Q为周期，即每个分量是以2e为周期的， )为 (在空间H ) 中的完备化．文 中不同地方出现的常数 C大小可能不相等，但都是与 无关，我们假定feL(0，；H一(Q)) 和 ／t~ )eL (YD ． 问题 (1)的弱形式为 ∈ =(；veL(0， (Q))，121eL(0， ( (Q))))， 收稿 日期：2008一l1—01 基金项 目：河南工程学院校级青年基金项 目(Y2007002)：贵州民族学院校级科研项 目(GzNC2008JS01) 作者简介：孙艳萍(1977一)，女，河南郑州人，讲师，硕士，研究方向为有限元、多尺度方法； 王 自强 (1981一)，男，贵州贵阳人，讲师，硕士，研究方向为有限元、多尺度方法． 河 南 科 学 第27卷 第3期