《有限元数值求解微分方程原理及其约束条件的处理方法xx》由会员分享,可在线阅读,更多相关《有限元数值求解微分方程原理及其约束条件的处理方法xx(36页珍藏版)》请在金锄头文库上搜索。
1、有限元数值求解微分方程原理及其约束条件的处理方法xxStillwatersrundeep.流静水深流静水深,人静心深人静心深Wherethereislife,thereishope。有生命必有希望。有生命必有希望微课内容摘要 在讲述有限元数值求解微分方程过程原理基础上,针对各种初值与边界约束条件,重点提出了一种可行的处理方法。该方法简洁清晰,易学易用,对学习理解有限元法具有启发性。同时,这种方法具有可靠的依据,加以引申可拓展有限元对约束条件的适定范围,其理论推导值得深思。1.1 一般二阶线性微分方程形式1、预备知识、预备知识例1: 具体的二阶线性微分方程及约束条件改写1.2 最简积分型泛函1、
2、预备知识、预备知识例2:1、预备知识、预备知识1.2 最简积分型泛函有有限限元元法法数数值值求求解解微微分分方方程程过过程程原原理理伽辽金变分方法2、泛函方程近似求解、泛函方程近似求解近似解思路:就是将无限维函数解空间降维处理。 具体说,就是构造有限N维函数空间,在每一维函数坐标确定一个基函数 ,同时每个基函数满足初始或边界条件。N个基函数要求线性无关,由这N个基函数进行唯一的线性结合,作为连续泛函方程的近似解。 随着N维数增加,解的子空间扩大,近似解的精度提高。2.1 近似求解思想 如左图所示,用4维函数空间的4个基函数的线性组合去近似更高维(或无限维)的函数图1 4维函数空间选取的4 个基
3、函数 因为每个基函数满足约束条件,所以连续的Galerkin变分形式可不必考虑约束条件。例4:两点边值问题: 暂不必考虑边界两点约束条件不必考虑边界两点约束条件,上述定解问题对应的等价泛函形式: 2.2泛函方程近似求解实例微微分分方方程程有有限限元元法法数数值值求求解解过过程程原原理理3 有限元法求解微分方程及约束条件处理有限元法求解微分方程及约束条件处理xxxxx 3.1 限元法思想限元法思想 将解空间的区间分段,解函数在每段上用拉格朗格插值基函数进行拟合,最后进行约束处理。1.021.251.751.51.751.251.02.01.5y(1)y(1.25)y(1.5)y(1.75)y(2
4、)YX叠加形成叠加形成 “总体刚度矩阵总体刚度矩阵” 和和 “总体截荷向量总体截荷向量”:微微分分方方程程有有限限元元法法数数值值求求解解过过程程原原理理3.2、微分方程约束条件的处理、微分方程约束条件的处理例4(续)两点边值问题: 现考虑边界两点约束条件,即对总体刚度矩阵与截荷矩阵附加约束。但先要对所给的初始边界条件改写,使之与第三类约束条件在形式上对应相符:*初始及边界条件处理技巧 1)总结: a)本微课中的对约束条件的处理适用于如下4个式子的两两任意组合。 4 总结与思考总结与思考 b)本微课中的对约束条件的处理方法的理论依据是微分方程对应的变分形式。伽辽金变分方法 2)思考:a)微分方程的第三类约束条件是如何产生?它有何理论意义或含义? b)讨论与作业:下两点边值问题的约 束条件如何处理,并求其数值解。*附有相关的Matlab计算程序可用 谢谢大家欢迎交流E-mail: