一种有效的广义特征值分析方法.pdfVIP

Vol. 24 No. 4 　　 　　　　　　 　 　　　　　　　　　 　 第24卷 第4期 固 体 力 学 学 报 2003年 　12月 ACTA MECHANICA SOLIDA SINICA December 　2003 一种有效的广义特征值分析方法 刘寒冰 　　龚国庆 　　刘建设 ( 吉林大学南岭校区交通学院 ,长春 ,130025) 摘 　要 　提出了一种适合于自适应有限元分析中求解广义特征值问题的多重网格方法. 这 种方法充分利用了初始网格下的结果 ,通过插值或最小二乘拟合技术来得到网格变化后的新的 近似特征向量 ,然后由多重网格迭代过程实现对结构广义特征值问题的求解. 在多重网格迭代 的光滑步中 ,选择了收敛梯度法以提高其收敛率 ;在粗网格校正步中 ,则导出了一种近似求解特 征向量误差的方程. 这种方法将网格离散过程和数值求解过程很好地相结合 ,建立了一个网格 细分后广义特征值问题的快速重分析方法 ,与传统有限元方法相比较 ,具有计算简便、计算量少 等特点 ,可以作为结构动力问题自适应有限元分析的一种十分有效的工具. 关键词 　多重网格迭代 ,广义特征值问题 , 自适应有限元 0 　引言 用自适应有限元方法分析结构动力问题时 ,有限元网格经常会发生变化 ,若用传统的有 ( ) 限元方法 如子空间迭代法 去分析自适应过程则存在一定的缺陷 , 因为传统的有限元方法 ( ) 包含的两个过程 即离散过程和数值求解过程 是相互独立的 ,两者无相互作用 ,这就造成了 ( ) ( ) 很大的浪费 : 1 没有充分利用初始网格下的结果 ; 2 当网格逐步加密时 ,传统有限元方法 的计算量会很快增大 ,当网格很多 , 自由度数很大时,计算量会更大 ,但收敛效果不一定好 , 这种缺陷对于动力问题会更明显. 尽管目前对于特征值问题发展了一些有效的算法 , 如 Lanczos 方法[1 ,2 ] 、Arnoldi [3 ] 方法等 , 以上的缺陷仍然存在. 多重网格方法则弥补了传统有限元方法之不足 ,它将离散过程和数值求解过程相结合 , 利用初始网格下的结果 ,通过多重网格迭代 ,可以达到简化计算过程 ,节省计算量的目的. 目 前 ,多重网格技术已经发展到了可以应用于实际工程问题 ,商用软件也已出现 ,线性和非线 性多重网格求解器建立了高效的求解偏微分方程的方法. 用多重网格方法求解特征值问题 也引起了广泛的关注. 下面回顾一下多重网格技术在特征值问题方面的研究状况. 传统的多重网格方法求解特征值问题是将此问题转换成非线性问题 ,然后利用非线性 ( ) [4 ] 多重网格求解器 ,如“全近似方法”full approximation scheme 进行求解 . 这种求解器应用于 非线性问题通常是有效的[5 ,6 ] ,但对于线性特征值问题 ,这种技术未必总是有效的 ,因为它没 有充分利用特征值问题的一些特性 ,而且 ,非线性多重网格特征求解器通常要对特征簇进行 特殊的处理 ,这又进一步使程序复杂化了. 此外 ,此技术的另一个缺陷就是对粗网格解的近 似程度要求较高. ( ) 国家自然科学基金 资助. 收到第 1 稿 收到修改稿.