有限元方法读书报告.docVIP

有限元读书笔记 在2009年暑假期间，在汤老师的指导下，我学习了有限元这门课程，了解了有限元的基本概念和基本方法。利用课外书籍、资料也了解了有限元方法的相关知识和应用情况，成为一种丰富多彩、应用广泛且实用高效的数值分析方法，对我们方程组迭代法研究过程中有着极为重要的影响。 一、基础知识 在课程学习中，我们学习了变分原理与变分法，通过最速降线问题、短程线问题、等周问题引入变分法的概念，即研究泛函极值问题就是变分法。研究变分问题的求解，变分问题旨在讨论相应的约束条件和边界条件下的泛函极值问题和多个独立变量的变分问题，而对于这类变分问题已有许多求泛函极值的必要条件和充分条件，如Euler方程组等，就此分类讨论一元（）以及二元（）状态下的求解方法，在学习中，我掌握了变分法基本引理，求泛函变分的参数导数定理以及泛函取极值的必要条件，也掌握了泛函求极值的直接算法与过程，也了解了具有一个自变量，一个自变函数和一阶导数的泛函的Euler方程；具有一个自变量，一个自变函数和高阶导数的泛函的Euler方程；具有多个自变量，一个自变函数泛函的Euler方程，以及参数变分问题。在变分问题求解学习中，对于Dirichlet原理的学习，求解，在已知，时，得到Euler方程，同时找到序列，取，再引入坐标函数（满足线性无关和完备性），使得，代入，求得参数为取得极值的关键所在解是唯一的，且等价于取极小值解。 定理1的应用非常重要，在研究微分方程边值问题等价的变分过程中十分关键，例：，存在齐次边界条件，首先对算子进行研究，验证其是否是线性正算子，应用定理1可解得，若边界条件为非齐次，则可化解成齐次边界条件，方法十分巧妙。 对于偏微分方程边值问题等价的变分，则对其边值问题进行系统分析，要考虑三种基本边值问题的边界条件： （1）Dirichlet问题： （2）Neumann问题： （3）Robin问题： 在三类边值问题的条件下，应用定理1，使得研究变分问题，而Robin问题又可划分为三类变分问题，对于的取值不同，则分解为几类问题，其中一类涉及到非齐次边界条件时，应采取引入坐标函数的方法求解，但区别的是此时引入的坐标函数需满足边界条件，区别于Ritz法引入坐标函数的概念。由此得出Galerkin法和Ritz法的优劣性；Galerkin法与Ritz法得到的结果是一致的，但是Galerkin法适用于更广泛的一类微分算子；Ritz法所选坐标系，在所满足的边界条件中不必考虑自然边界条件，而Galerkin法需满足自然边界条件；两个方法的背景不相同。 在前期大量的基础学习中，第二章正式进入有限元理论的理论基础学习，首先对泛函分析里的线性空间、逆算子、Sobolev空间以及弱收敛性、紧致性概念进行了回顾，对等价模定理进行了更深一层次的分析。最后研究了有限元方法，对出Galerkin法和Ritz法进行了改进，采取分片多项式的方法，进行单元插值运算，因此总结出有限元求解问题的基本步骤： 第一步：问题及求解域定义：根据实际问题近似确定求解域的物理性质和几何区域。 第二步：求解域离散化：将求解域近似为具有不同有限大小和形状且彼此相连的有限个单元组成的离散域，习惯上称为有限元网络划分。显然单元越小（网络越细）则离散域的近似程度越好，计算结果也越精确，但计算量及误差都将增大，因此求解域的离散化是有限元法的核心技术之一。 第三步：确定状态变量及控制方法：一个具体的物理问题通常可以用一组包含问题状态变量边界条件的微分方程式表示，为适合有限元求解，通常将微分方程化为等价的泛函形式。 第四步：单元推导：对单元构造一个适合的近似解，即推导有限单元的列式，其中包括选择合理的单元坐标系，建立单元试函数，以某种方法给出单元各状态变量的离散关系，从而形成单元矩阵（结构力学中称刚度阵或柔度阵）。 为保证问题求解的收敛性，单元推导有许多原则要遵循。对工程应用而言，重要的是应注意每一种单元的解题性能与约束。例如，单元形状应以规则为好，畸形时不仅精度低，而且有缺秩的危险，将导致无法求解。 第五步：总装求解：将单元总装形成离散域的总矩阵方程（联合方程组），反映对近似求解域的离散域的要求，即单元函数的连续性要满足一定的连续条件。总装是在相邻单元结点进行，状态变量及其导数（可能的话）连续性建立在结点处。 第六步：联立方程组求解和结果解释有限元法最终导致联立方程组。联立方程组的求解可用直接法、选代法和随机法。求解结果是单元结点处状态变量的近似值。对于计算结果的质量，将通过与设计准则提供的允许值比较来评价并确定是否需要重复计算Galerkin法和Ritz法出发，研究有限元方法，并初步掌握了有限元方法的核心。 二、知识的强化 在学习变分法的过程当中，我们一般情况下研究的是传统的变分问题，传统的变分法所讨论的泛函指被积函数只含