第二章有限单元法的数学基础重点.pptVIP

第二章 有限单元法的数学基础 有限元最初是由结构工程师们基于物理概念提出来的，但这种被称为直接刚度法的方法仅适用于力学领域，很难将其推广到其它领域 。 自然界中的许多规律都是以微分方程的形式提出的，微分方程是关于场变量及其导数的一个方程，微分方程和相应的初边值条件一起构成了所谓的微分方程的定解问题。 数值方法，粗略的理解就是一种近似方法，但这种近似方法不是任意的，而是能够通过数值方法本身所提供的某种手段，通过加大运算量而使得近似解越来越逼近解析解。 有限元法的先行课程(按依赖程度)：线性代数、数理方程、数值分析 + 一门连续统力学，这门连续统力学对不同学科有不同的要求，对于土木工程学科，可以是弹性力学或结构力学。 引入目的： 书写和推导公式的简便 ； 阅读其它采用张量记法的有关文献 。 加权余量法示例 该问题的解析解为 什么是泛函 如果对某一类函数中的每一个函数v(x)，按照一定的变换规则 J，都有一个确定的实数，记为J(v)，与之对应，那么就称该 变换是一个泛函。 泛函就是由函数空间到一维欧氏空间R的一个变换，泛函的“自 变量”或者“原象”取自某一类函数，而变换的结果是一个实数。 例如 由于其满足了 所以可以称为线性泛函。 下面两个属于非线性泛函。 什么是泛函的变分 如果泛函J 的增量△J 可表示为 则称泛函J可微，并称 假设J 可微，则对于任一实数a ,有 其中 关于δv是线性的, 为一阶变分，简称变分，记为 等式两端除a 并令 泛函极值的充分条件和必要条件 如果函数u(x)∈M 使得对应的J(u)取极值，则 仿照函数的Taylor展开 的极值点为 函数取极值的必要条件： 泛函取极值的必要条件 一阶变分 二阶变分 若 且 ≥0(或≤0),则u(x)是J(v)的一个极小或极大函数 泛函取极值的充分条件 若 但 随δu有正有负，则u(x)是J(v)的一个驻点。 * 有限单元法方程的基本思想 按其推导方法分: (1)直接法 根据单元的物理意义，建立有关场变量表示的单元性质方程。 (2) 加权余量法 直接从控制方程中得到有限单元方程，是一种近似解法。 (3)变分法 直接从求解泛函的极值问题入手，把泛函的极植问题规划成线性代数方程组，然后求其近似解的一种计算方法。 有限单元法的数学基础 我们今天要面临的问题，不仅限于力学问题，要涉及到其它许多相关领域的问题，如，当我们运用数值方法模拟降雨型滑坡的成因和机理时，我们不只是考虑应力场，还需考虑地下渗流场、地表径流等诸多因素的共同作用。 变分法是建立有限元法的通用方法，同时变分法还是其它现代数值方法，如边界元法、数值流形法、有限体积法等的基础。 有限单元法的数学基础 三大类微分方程： 椭圆型方程用于描述与时间无关的稳态问题； 抛物型方程的典型应用是热传导问题； 双曲型方程的典型对象则是振动问题 由于问题的复杂性，能够获得解析解的问题是非常有限的。为了求解一般情况下的定解问题，就需要借助于某种数值方法。 有限单元法的数学基础 有限元法就是求解微分方程的一种数值方法，而且它特别适用求解椭圆型边值问题。 在有限元法问世之前所广泛采用的一种数值方法是有限差分法，有限差分法是一切微分方程数值方法的基础，至今仍在某些领域如流体力学问题中起着主导作用。有限差分法的基本思想就是用差商代替微商，而使得原来的定解问题的求解变为线性方程组的求解。所有的微分方程数值解法都是化微分方程定解问题为线性方程组的求解问题。 有限单元法的数学基础 Einstein求和约定 首先将三维空间的坐标x，y，z分别记为 Einstein求和约定：如果表达式中的乘积项中含有两个相同指标，则表示这个指标遍历1，2，…，n，n = 2或3，是空间的维数，而后将这n个同类型的表达式相加而得到的总和。如果相同的指标有好几对，则表示每一对指标做上述求和。 Einstein求和约定 相同的指标是求和指标，称为哑指标， 而其余的指标称为自由指标。 求和约定中求和指标遍历的是空间维数n，若在一个求和表达式中的求和指标不是遍历空间维数n，则我们不应省略求和号? 求和约定的好处示例：三维变形体的平衡方程 变分法的基本引理 子空间 如果存在 在区间[a,b]上有连续n 阶导数的全体函数空间为 假设 则在区间[a,b]上 变分法的基本引理 证明：采用反证法。 设存在一点 根据函数的连续性存在 构造连续函数 但 微分方程的等效积分形式 工程和物理学中的许多问题，都可以以未知