第10章 有限元法基础.ppt

10.1 有限元法概要 参考书： 《传热的有限元计算方法》孔祥谦 例1 最速下降线问题 A、B两点之间可以画许多条曲线，质点沿着哪条曲线运动所需的时间最短？ 显然不是直线段！！！ 最速降线问题的数学描述 曲线： 速度： 弧长： 时间： 总时间： 端点条件： 函数与泛函 函数：对于自变量 x，有一个值与其对应 泛函：对于自变量 y(x)，有一个值与其对应 微分 与 变分 函数微分：让自变量 x 略微变动一点，看函数 y(x) 相应地变化多少。 泛函变分：让自变量 y(x) 略微变动一点，看泛函 J(y) 相应地变化多少 泛函 与 变分 自变量的变分： 泛函的变分： 泛函变分 多元泛函 变分： 2 变分的运算法则 1）和差的变分 2）积的变分 3）商的变分 4）微分的变分 5）积分的变分 变分的计算 对于泛函 考虑自变量的小增量： 泛函取极值的条件 10.2.2 微分问题与变分问题的等价关系 1、函数内积与算子 1）内积 2）算子 ? 一种运算 内积空间中某一集合 D 和另一集合 R 建立了某种一一对应的关系。 3） 线性算子 （4）伴随算子 对内积进行分部积分 5） 自伴算子 （6） 本质边界条件和自然边界条件 本质边界条件：规定 F(u) 自然边界条件：规定G(u) 二阶偏微分方程： 本质边界条件：是规定函数值(即第1类BC)， 自然边界条件：是规定导数值(即第2、3类BC) （7） 对称算子 （8） 正算子 （9） 算子方程 例2 负Laplace算子 线性算子？ 内积分部积分 正算子？ 结论：线性对称正算子！ 变分原理 2 对称正算子方程的变分原理 3. 自伴算子方程非齐次边界条件的变分原理 4. 变分问题转化为微分问题——Euler方程 同一物理问题 微分方程描述 变分描述 此泛函取极值 此泛函取极值 微分方程 泛函 第10章 有限元法基础 胡 茂 彬 /~humaobin/ humaobin@ 10.1 有限元法概要 离散化：“分块逼近” ? “总体合成” 代数化：基于变分原理、加权余量法，将微分方程转化为代数方程 固体力学计算中用有限元法较多，而流体计算一般用有限容积法。流固耦合计算中，流体部分用有限容积法，而固体部分用有限元 10.2 变分原理和 Ritz 法 10.2.1 变分原理 如果存在于某个函数空间 D 中的物理状态函数 u，决定了一个依赖它的泛函 J(u)，则在一切可能的物理状态中，真实的物理状态是 D 中使泛函 J(u) 取极值的那个函数。 泛函：函数的函数，以函数为自变量 变分问题：求泛函极值的问题 变分问题 光线沿着用时最短的路径传播 - 费马原理 Lagrange 最小作用量原理 Hamilton原理、最小势能原理 最速下降线问题 同一物理问题 微分方程描述 变分描述 定义域 D 值域 R 算子L 边界 对称算子！