基于第二代小波理论有限元自适应方法.pdfVIP

基于第二代小波理论的有限元自适应方法1 1 1 2 何育民 ，陈雪峰 ，何正嘉 1 西安交通大学机械工程学院 (710049) 2 西安交通大学机械制造系统工程国家重点实验室 (710049) E-mail：he_yumin@163.com 摘 要：基于第二代小波理论，建立了有限元多分辨分析方法，构造了相应的自适应算法。 采用第二代小波作为基函数构造了逐级嵌套的有限元逼近空间，并引入小波消失矩建立了有 限元求解方程刚度矩阵解耦的条件。通过消除刚度矩阵中不同尺度空间之间的耦合项，可以 实现有限元求解方程按空间尺度解耦，使得方程可以在各个空间内独立求解。低分辨空间可 以通过增加细节空间提升到高分辨空间，从而方程可以在低分辨空间求得初始解，然后通过 不断增加在细节空间的解，得到高分辨空间的解，从而满足精度要求。轴力杆的自适应分析 算例验证了方法的有效性。 关键词：第二代小波，多分辨分析，自适应算法，有限元 1．引言 Daubechies 小波、样条小波等常常作为基函数来构造有限元逼近空间[1,2]，低分辨空间 通过增加细节空间可以提升到高分辨空间，但是低分辨逼近空间和细节空间之间以及不同尺 度的细节空间之间均存在耦合项，为了消除耦合项，实现平滑逼近和细节信息的分离，本文 采用 Sweldens 于 1996 年在时域中采用提升方法(Lifting scheme)构造出一类小波，称为第二 代小波[3,4]，作为基函数来构造有限元逼近空间。通过设计预测器和更新器可以构造满足特 定要求的第二代小波小波，提供适应不同情况的多分辨小波基[3-5] 。第二代小波已经被应用 于很多领域，文献[6]研究了运用第二代小波配置法求解偏微分方程的问题，S. D Heedene 和 K Amaratunga 等人在文献[7]中提出了基于第二代小波的有限元分析思想。本文建立了基 于第二代小波的有限元多分辨分析方法，并给出了相应的自适应算法。 2 ．第二代小波变换的原理 第二代小波构造不依赖于傅立叶变换，其对信号的分解和重构在时域中完成[3-5] 。其分 解过程为： （1）剖分：将分辨率为j +1 的一组信号uj +1 分成分辨率为j 的偶样本uj 和奇样本 uj 。 e o （2 ）预测：用偶样本预测奇样本，预测估计的误差为 j () j 分辨率下的细节信号r ，P ⋅ 为预 测器。 r j uj =−P uj (1) ( ) o e j () （3 ）更新：用预测误差修正偶样本，得到j 分辨率下的逼近信号u ，U ⋅ 为更新器。 uj uj =+U rj (2) e ( ) 重构过程为分解过程的逆过程，包括恢复更新、恢复预测和合并，可直接由分解过程得到。 3 ．自适应有限元分析 { } 多分辨分析 V 提供一个互相嵌套，逐级扩大的子空