总结归纳有限元步骤.docxVIP

总结归纳有限元步骤 　　有限元基础理论学习总结报告　　中国矿业大学14级硕士王涛　　通过课上和课下的学习，对有限元基础理论有了一定的了解和认识。经过学习，更加深刻的理解了有限元的离散、单元类型、插值函数构造和等参变换等知识，现对有限元的基本理论和用法做了如下学习和报告。　　已经发展的偏微分方程数值分析方法可以分为两大类。一类是有限差分法，其特点是直接求解基本方程和相应定解条件的近似解，求解步骤归纳为：首先将求解域划分为网格，然后在网格的节点上用差分方程来近似微分方程。借助于有限差分法能够求解相当复杂的问题，特别是求解方程建立于固结在空间的坐标系的流体力学问题，有限差分法有自身的优势，因此在流体力学领域内，至今仍占支配地位。但是对于固体结构问题，由于方程通常建立于固结的物体上的坐标系和形状复杂，另一类数值分析方法——有限元法则更为合适。　　有限差分法：　　特点：以差分方程近似微分方程，直接数值求解原问题的微分方程，　　在流体力学，岩土力学领域占重要地位。　　有限元法：　　特点：区别于有限差分法，即不是直接从问题的微分方程和相应的定　　解条件出发，而是从等效的积分形式出发，数值求解原问题的　　等效积分方程。　　基本思想：1将求解域离散为有限个子域的集合　　2分片逼近待求函数　　分析过程：1单元特性分析，单元节点位移与节点力之间的关系2系统特性分析，将单元刚度矩阵集成整体刚度方程　　1.有限元法的理论基础——加权余量法和变分原理　　微分方程的等效积分形式和加权余量法　　微分方程的等效积分形式　　工程或物理学中的许多问题，通常是以未知场函数应满足的微分方程和边界　　条件形式提出来的，可以一般地表示为未知函数应满足微分方程组　　A(u)?0　　域Ω可以是体积域、面积域等。同时未知函数还应满足边界条件　　B(u)?0　　Г是域Ω的边界。　　由于微分方程组在域Ω中每一点都必须为零，因此就有　　?　　其中??TA(?)d????(?1A1(?)??2A2(?)?...)d??0　　是函数向量，它是一组和微分方程个数相等的任意函数。式与微分方程组式是完全等效的积分形式。同理，加入边界条件也同时在边界上每一点都得到满足，则其等效积分形式为　　T???A(?)d????B(?)d??0()T　　对分部积分得到等到另一种形式??CT(?)D(?)d????E()F(?)d??0()　　其中C、D、E、F是微分算子，它们中包含的阶数较式的A低，这样对函数只需要求较低阶的连续性就可以了。在式中降低的连续性要求是以提高?和的连续性要求为代价的。这种通过适当提高对任意函数?和的连续性要求，以降低对微分方程场函数的连续性要求所建立的等效积分形式称为微分方程的等效积分“弱”形式。　　基于等效积分形式的近似方法——加权余量法　　对微分方程式和边界条件式所表达的物理问题，假设未知场函数可以采用近似函数来表示。近似函数是一族带有待定参数的已知函数，一般形式是　　?????Niai?Na　　i?1n　　其中ai是待定参数；Ni是称之为试探函数的已知函　　数，它取自完全的函数序列，是线性独立的。　　显然，近似解不能精确满足微分方程式和全部边界条件式，它们将产生残差R和R，即A(Na)?R;B(Na)?R。残差R和R亦称为余量。在式中用n个规定的函数来代替任意函数?和，即　　??Wj；?Wj(j?1~n)　　Wj和Wj称为权函数。　　对应等效积分“弱”形式式，同样可以得到它的近似形式为?C?T(Wj)D(Na)d???ET(Wj)F(Na)d??0(j?1,...,n)?　　采用使余量的加权积分为零来求得微分方程近似解得方法称为加权余量法。对于权函数不同的选择可分为配点法，子域法，最小二乘法，力矩法和伽辽金法。　　变分原理　　如果微分方程具有线性和自伴随的性质，则不仅可以建立它的等效积分形式，并利用加权余量法求其近似解，还可以建立与之相等效的变分原理，并进而得到基于它的另一种近似求解方法，即里兹方法。　　线性、自伴随微分方程变分原理的建立　　1.线性、自伴随微分算子　　若有微分方程　　L(u)?b?0　　　　其中微分算子L具有如下性质　　L(au1??u2)?aL(u1)??L(u2)　　则称L为线性算子，方程为线性微分方程。其中a和?是两个常数。　　现定义L(u)和任意函数的内积为　　?　　随有边界项。结果可表示如下：?L(u)vd?对上式进行分部积分直至u的倒数消失，这样就可以得到转化后的内积并伴　　??L(u)vd???L*(v)d??(u,v)?　　L*称为L的伴(u,v)表示在?的边界Г上由u和v及其导数组成的积分项。　　随算子。若L*=L，则称算子是