李志青, 李远飞, 张文彬

(广州华商学院,广东 广州 511300)

近年来,在科学和工程的计算领域中,复合材料的出现对科学技术的发展产生了极大的推动作用.在新型复合材料的研发过程中,经常需要对一些复合材料进行等效性能评价分析,双尺度方法和均匀化方法因此得到了发展[1-11].以上研究主要对复合材料进行了多尺度分析、均匀化分析、渐近分析以及双尺度解和有限元解的误差估计.本文在以上文献的基础上进一步研究周期结构带阻尼项椭圆边值问题,就其有效算法和近似算法展开讨论，利用均匀化方法及高阶双尺度方法对周期结构带阻尼项椭圆边值问题进行了分析,给出了这类方程的双尺度渐近展开式并分析双尺度有限元解的误差估计,最后设计对应的有限元算法.

周期结构带阻尼项椭圆边值问题可以表示为

其中：

1)i、j、h、k=1,2,…,n.

2)∂Ωε=Γε满足Lipschitz边界条件,并且Γε=∂Ω∩εω,Ωε为与ε有关的小周期状区域,其中ω为以1为周期的无界区域,Ω为Rn中的有界闭区域.

(2)

(3)

1 θε(x)的双尺度渐近展开式及均匀化方程

假设问题(2)有如下的双尺度形式渐近展开式：

(4)

其中Hα1(ξ)、Hα1α2(ξ)、Hα1α2α3(ξ)、…在整个空间Rn上定义,它们在各个方向是关于ξ为周期的标量函数,且可以在单位胞体Q上定解.

联合(4)和(2),再由ε的任意性,通过计算并比较ε-1、ε0和ε1的两边系数及边界条件可知

(5)

而θ0(x)是问题(2)的均匀化解,满足

(6)

其中：

(7)

(8)

2 uε(x)的双尺度渐近展开式及均匀化方程

问题(3)的解uε(x)是与x、ξ和ε有关的函数,假设uε(x)有如下的形式渐近展开式：

(9)

其中，u0(x)是待定的充分……