复合材料特征值的高阶多尺度Rayleigh商校正.pdfVIP

2013年 11月 计 算 数 学 第 35卷第 4期 NOV．，2013 MATHEMATICA NUM ERICA SINICA Vo1．35，No．4 复合材料特征值的高阶多尺度 Rayleigh商校正术1) 张 磊 (解放军后勤学院后勤管理系，北京 100858； 中国科学院数学与系统科学研究院，计算数学与科学工程计算研究所，北京100190) 曹礼群 (LSEC，NCMIS，中国科学院数学与系统科学研究院，计算数学与科学工程计算研究所，北京100190) 摘 要 本文讨论了周期结构复合材料特征值的多尺度计算，提出了高阶多尺度 Rayleigh商校正算 法，并给出了收敛性分析．最后，通过大量数值实验结果表明，新算法是有效且必要的． 关键词：复合材料；特征值；高阶多尺度校正；Rayleigh商 MR (2000)主题分类：65D17 1．引 言 复合材料周期结构声波和振动问题可转化为以下具有周期振荡间断系数椭圆特征值问 题 三 一 鑫((詈) ) ( (脚( ()inQ， {乱x)=0onFo， (1．1) L ￡(。)=0 OilFi， 其中Q cR ，佗 1是有界区域，边界为 aQ=FouF1，FonFi= ，aQ是分片光滑的 Lipschitz连续边界．这里E0是周期性微结构的相对尺寸，即0 =争1，其中f，L分 别是复合材料结构周期单胞和整体区域的尺寸，如图1(1)和 (2)所示．这里P(吾)是质量密 度函数， ()三 毗j(罟) ．本文中采用Einstein求和约定． 我们对a (吾)，ao(罟)和P(吾)作如下假设： (A1)设 ∈=C--1X，0 (∈)，00(∈)和P(∈)关于 ∈是 1一周期函数． (A2)系数矩阵 (∈)=(aij(∈))对称正定． (As)aid，ao，P∈ (Q)，P(罟) Po0，ao(x) 0，其中Po是常数． 由于复合材料结构包含大量的夹杂或者孔洞，直接进行数值计算需要非常细的网格，从 而产生难以承受的计算规模和时间．为此，物理学家和力学家提出了各类平均化方法，数学家 发展了均匀化方法，其核心思想是寻找等效参量．如Kesavan[121]很早就研究了椭圆特征值 2013年 6月 4日收到． )基金项 目：国家 自然科学基金 90916027)和国家重点基础研究发展计划 (2010CB832702)项 目 资助． 432 计 算 数 学 问题的均匀化方法并给出了收敛性结果．大量数值模拟结果告诉我们当 ￡不是特别小的时候， 均匀化方法的数值精度远远无法满足实际工程计算的需要，见文献 3[-5，8，14]，于是很多多尺 度算法被提出．Oleinik等 [10．]讨论了具有快速振荡系数的一维 Sturm—Liouville方程，并给 出了特征值和特征函数的完全渐近展开．曹礼群和崔俊芝等 3[4]对一般区域及孔洞区域的椭 圆特征值问题给出了高阶渐近分析结果． 糕 麟 鞣 缀 瓣 麟 嬲 瓣 麟 麟 嬲 黼 黼 麟 黼 — — 麟 貉 麟