科学及工程计算有限元.ppt

矩阵的形式为 总刚度阵 位移向量 总荷载向量 1：插值函数 当精确函数 y = f(x) 非常复杂或未知时，但能测量 或者统计出在一列节点 处测得函数值 ， ， 根据这些已知条件，构造一个简单易算的近似函数g(x)， 满足条件 求一个至多n次多项式函数 利用待定系数法 因为方程组的系数矩阵为范德蒙特行列式， 存在而且唯一 插值多项式是一次多项式 求解方程组得 记 Lagrange基函数 个互异节点时, 构造出每个节点的Lagrange基函数 它们满足 由上面的条件 可以写为 求得常数 每个节点的基函数 插值多项式为 两点边值问题的有限元方法 考虑两点边值问题 第一步 将数学物理问题化为等价的变分问题 变分问题 求 第二步 定解区域剖分 对[a,b]进行剖分，使 每一个小区间称为单元。 第三步 有限元空间的构造,线性有限元空间 构造 的有限维子空间 应该满足如下要求 任何元素（函数）限制在任意一个单元上 是一次函数 任何元素（函数）满足一定的光滑性，在 是连续函数 在端点满足 因为在单元上是一次函数 根据插值多项式的结论. 几何图形 选取的函数满足 称为整体基函数。它们的函数图像分别为 函数列是线性无关的 因为，令 利用基函数的性质得 函数列可以成为子空间 的基底 称此空间为线性有限元空间。 存在实数列 做变换 则此变换将单元 变换到区间 称为参考单元 形函数 第四步 生成有限元方程 离散变分问题 写成矩阵形式 ， 因此得到有限元的系数矩阵每行至多只有三个非零元素 如果左边值条件为非齐次的，即 则得到有限元方程 计算单元刚度阵。它是 的二次型，写成 是单元刚度阵。 于是得到在整个定解区域上的有限元方程 总刚度矩阵 总荷载向量。 一维高次有限元 区间上构造二次Lagrange插值基函数 和 考虑 满足插值条件 表达式为 变换 单元的Lagrange基函数 下面构造在整个定解区域上的整体基函数 考虑节点 是单元 的共同的端点 将相邻两个单元上的Lagrange基函数拼接起来得到在节点 的整体基函数 单元中点的整体基函数就是在此单元上是Lagrange基函数 其形状为 二次有限元的近似解可以表示为 在每个单元上，近似解是二次多项式，在节点处连续。 因此，在整个定解区域近似解是连续函数。 二维椭圆边值问题的有限元方法 考虑二维椭圆边值问题 第一步 将数学物理问题化为等价的变分问题 变分问题 求 ， 第二步 定解区域剖分 三角形剖分 定解区域的边界分片光滑 ,如果也不是由直折线段组成， 那么截弯取直的办法。用适当的折线段逼近. 三角形剖分就是将定解区域分割为有限个小三角形单元。 三角形的顶点称为节点，记为 每一个三角形称为单元，记为 三角形剖分必须注意几点，即三角形剖分必须是正则的。 1:每个单元的顶点必须是另外一个单元的顶点，不能是 相邻单元的边上的点。原因是方便构造插值函数。 2:三角形之间内部无重叠，即 3:在剖分过程中尽量避免大钝角三角形，要求满足不等式 4:如果方程的系数在区域内部是间断的，则用折线近似间 断线，而且这些折线必须是三角形单元的边界。 5:如果定解区域的边界是由有限条光滑曲线组成的，则 这些光滑曲线的交点也必须是节点，如果在边界条件中 的系数也是间断的。则间断点也必须是节点。 6:对于任何一个单元，单元最多有两个顶点落在边界上。 注：（1）三角剖分具有很好的灵活性和适应性。 （2）区域剖分好后，必须对节点和单元编号 (3)子空间不是真解空间的子空间 第三步 有限元空间的构造 在上述剖分的基础上，构造有限维子空间 . 任何元素（函数）限制在任意一个单元上 是一次函数 任何元素（函数）满足一定的光滑性 是连续函数 在边界上满足 构造形函数 满足以下条件 在单元是一次函数 单元之外，每个形函数都等于0。 整体基函数，它是以节点 为顶点的所有单元上的形函数 拼接而成的，它是分片线性的，而且局部非零的。 它满足条件 有限元空间 在网格步长 有限元空间收敛到广义解空间。 第四步 形成有限元方程 ，