daulechies条件小波有限元法.docxVIP

daulechies条件小波有限元法 1 小波有限元法 传统的基于有限值法在计算工程值的领域发挥着重要作用。随着这一理论的不断发展和更新，该理论在工程中的应用日益深入。但是,由于传统有限元法一般采用幂函数作为位移场函数的试探函数,在某些工程问题中当位移场函数不是幂函数时,采用传统有限元法进行求解将产生较大的误差,而利用小波函数插值可以获得新的有限元逼近空间 小波有限元法吸收了传统有限元法离散逼近的优点,可以方便地处理复杂的边界条件;同时又拥有小波函数特有的多分辨特性,可以提供另一种提高精度的细化算法,即在不改变网格剖分的前提下提高其分辨率。小波函数最大的特点就在于它具有多尺度、多分辨、紧支性的特点。因此,小波有限元法可以根据实际需要任意改变分析尺度:在变化梯度小的求解域用大的分析尺度;在变化梯度大的求解域则采用小的分析尺度。在这方面,小波有限元法可以弥补传统有限元法的不足;并且小波有限元法算法的数值稳定性好、运算速度快、求解精度高,因此有必要对其展开研究。 国内外针对小波有限元法进行了大量研究工作 本文主要研究Daubechies条件小波有限元法。该法区别于其他参考文献中的常规小波有限元法,是以已有研究成果Daubechies条件小波Ritz法 2 有限元法的构建 本文以结构中常见的受拉(压)弹性杆、平面受弯梁为基本研究对象,阐述Daubechies条件小波有限元法的构建过程。 2.1 拉紧弹性套管 图1为两端固结且受均布荷载的弹性拉压杆,f(x)为所受均布荷载,σ 2.1.1 基于lagrange函数的修正泛函 首先,建立弹性拉压杆的局部坐标系。设:杆的长度为l;左右端点的总体坐标分别为 显然有,局部坐标ξ的取值范围为0≤ξ≤1。 其次,在局部坐标中预设Daubechies小波基位移场函数。 其中:u 根据p阶Daubechies小波基函数的支撑域特性,令N=2p-1(注:若无特殊说明,文后所用字母N均有N=2p-1成立,其中p为小波函数消失矩的阶数),则对于 设λ为Lagrange乘子,为约束边界Γ域中的一组独立坐标的函数向量。由于图1中的约束边界Γ是两个端点,则可简化为 其中b为边界条件矩阵。 采用Lagrange乘子法将需满足的边界条件E(u)引入自然变分条件下的泛函,可以形成新的修正泛函,即 其中:λ 其中:E为弹性模量;A为构件的横截面积。为体现荷载的普遍性,将f(x)设为幂函数,即令 其中:q为荷载集度系数;r为幂函数指数。若构件受均布荷载,则r=0。E(u)为必须满足的边界条件,且有 其中:u为位移函数; 由修正泛函变分为零可得到下列一组方程 得 其中:Φ′为小波基函数一阶导数;(Φ′) 这里 2)左端固结、右端自由 3)左端自由、右端固结 为方便组装整体求解矩阵,可将式(10)中的b 式(13)即为所构建的单元刚度矩阵。 另外,在上述求解过程中涉及了Daubechies小波尺度函数导数的求解问题,由于Daubechies小波为离散小波,其尺度函数和小波函数不存在显式表达式,因此其导数计算也必须采用特殊方式进行。基本思路借助两尺度方程进行隐式求导,具体求解方法可参考文献 2.1.2 单元求解矩阵 进行整体求解矩阵的组装时,首先构造中间单元求解矩阵。以两端固结边界条件为例,为方便构造整体求解矩阵,将第一行及第一列取负值,即 其次,构造端部单元求解矩阵。左端单元求解矩阵为 右端单元求解矩阵为 再次,根据弹性拉压杆的特点,按公共节点位移相等的原则组装整体求解矩阵,基本过程如式(17)所示 式中(0-0)表示矩阵合成时要组合的两个矩阵在相应位置处元素互相叠加。最后,由整体求解矩阵解得各单元的小波基待定系数a 2.2 竖向位移函数 图2为两端固结且受线性荷载作用的平面弯曲梁。图中:I为截面惯性矩;v(x)为竖向位移函数。同样,实际工程中的边界条件及荷载可在求解过程中灵活变换。 2.2.1 dawellchies小波基函数的列阵 建立局部坐标系的过程与弹性拉压杆完全一致。根据能量法建立平面弯曲梁(注:此处指长度远远大于其截面尺寸的细长梁,忽略剪切应力对变形的影响)的泛函表达式见式(18)。同理,为体现荷载的普遍性,将f(x)设为幂函数,即f(x)=ql 其中:Φ′为Daubechies小波基函数Φ(ξ)对局部坐标ξ求二阶导数所组成的列阵;(Φ′) 两端竖向支撑的位移边界条件可描述为 若梁的两端固结,则两端转角θ为0,即 以边界条件两端固结为例,将式(18)～式(21)代入式(5)中,得 式中,C=[(v-v)(v 由修正泛函变分为零可得到下列一组方程 式中,D=[φφ 进而得到 式中 则有 其中: 1)两端固结 2)两端铰结 3)左端固结、右端铰结 4)左端固结、右端自由 同样,考虑到条件小波有限元组装整体求解矩阵的方便