有限元第7章.pptVIP

第八章 流体流动问题的有限元方法 控制方程 伽辽金加权余量法: 有一个待定常数 代入控制方程 第八章 流体流动问题的有限元方法 两个未知数 代入控制方程 伽辽金加权余量法: 控制方程 根据边界条件,设 边界条件 取前2项 为形函数. 代入 代入可得到 * 泛函中各函数的确定 2.坐标变换（八节点等参数单元） 1)雅可比矩阵（八节点等参数单元） 泛函中各函数的确定 2)面积积分的变换（八节点等参数单元） 令 3)线积分的变换（八节点等参数单元） 泛函中各函数的确定 在1,8,7和3,4,5边 3)线积分的变换（八节点等参数单元） 在1,2,3和5,6,7边 温度刚度矩阵的形成 令 代入 温度刚度矩阵的形成 用高斯求积法 温度刚度矩阵的形成 求泛函极值,上式对温度Ti求偏导数 温度刚度矩阵的形成 求泛函极值,上式对温度Ti求偏导数 将上式展开,写成矩阵形式 温度刚度矩阵的形成 右端项的形成 泛函中各函数的确定 在1,8,7和3,4,5边 线积分的变换（八节点等参数单元） 在1,2,3和5,6,7边 右端项的形成 右端项的形成 右端项的形成 右端项的形成 三 总体合成 为单元温度刚阵的最后形式 求解区域内的总温度泛函等于各单元温度泛函之和,从而有 为按节点号叠加的总体温度刚度矩阵. 为按节点号叠加的总体节点温度向量. 为按节点号叠加的总右端向量 三 总体合成 根据泛函求极值的条件 泛函求极值问题变成多元函数的极值问题 总结: 1.将描述物理现象的微分方程及其边界条件通过变分原理转化为某一泛函求极值问题.再通过一定的推导得到有限元计算格式.这一方法可以应用到其它场问题的有限元求解中.重要我们找到与微分方程及其边界条件等价的泛函,就可以用变分原理. 2.泛函求极值问题与微分方程边值问题比较,有边界条件自动满足,导数阶次降低,推导有限元格式方便等优点. 3.推导方法也可用于不同单元节点数的二维或三维等参数元计算格式,也可用于其它场问题的有限元计算格式推导. 习题: 1.在稳定温度场有限元计算格式推导过程中,泛函求极值问题是如何转换成多元函数求极值问题的? 2.简述利用变分原理推导稳定温度场有限元计算的主要思路和基本过程. 第八章 流体流动问题的有限元方法 一 问题的提出 1. 采用变分原理可以推导某些偏微分方程所描述的非结构问题的有限元格式.但是泛函求极值问题和微分方程求解问题只是在一定的条件下等价,既要求求解函数连续且可导; 2. 给出一个泛函形式,通常可以找到相对应的微分方程;反之,若给出一个微分方程,并不一定能够找到与之对应泛函极值问题的具体形式. 要用有限元方法求解某一场问题的微分方程,如果找不到与之相对应的泛函,则无法用变分原理来推导相应的有限元计算格式. 第八章 流体流动问题的有限元方法 很多问题都不能用变分法 流体流动问题 高速车辆运行时的空气阻力问题 列车进入隧道时的压力波动问题 两高速运动的列车会车时的压力波动问题 桥梁的风致震动问题 汽车高速行驶时的气动稳定性问题 方程为非线性偏微分方程 第八章 流体流动问题的有限元方法 求解方法:加权余量法 思想: 某物理问题的控制微分方程及边界条件分别为 在域 内 在域 边界S上 为待求函数. 选择一试探函数 待求常数 试探函数项 第八章 流体流动问题的有限元方法 将试探函数代入控制方程和边界条件，在域内和边界上会产生误差 在域 内 在域 边界S上 为余量（残数、残差、残值）。 基本思想：在域内或边界上寻找n个线性无关和函数 使余量R和Rb在加权平均的意义上等于零 第八章 流体流动问题的有限元方法 尽管试探函数本身不能满足控制方程和边界条件。但是当其雨量与许多线性无关的权函数相乘并积分十，这个余量在总体上接近于零，也就是说试探函数在积分意义上满足微分方程式和边界条件。当n足够大时，试探函数就趋近于真实解。 基本方法：1.最小二乘加权余量法 首先找到一个已经满足边界条件式的试探函数 代入控制方程会产生余量 第八章 流体流动问题的有限元方法 基本方法：1.最小二乘加权余量法 代入控制方程会产生余量 简单地说，最小二乘的思想就是要使得观测点和估计点的距离的平方和达到最小.这里的“二乘”指的是用平方来度量观测点与估计点的远近（在古汉语中“平方”称为“二乘”），“最小”指的是参数的估计值要保证各个观测点与估计点的距离的平方和达到最小. 希望余量在最小二乘的意义下为最小. 第八章 流体流动问题的有限元方法 方差为: 由此可知道: 方差最小的条件: 通过求解线性方程组,可得到待定系数. 第八章 流体流动问题的有限元方法 控制方程 最小二乘加权余量法: 有一个待定常数 代入控制方程 对于复杂工程问题,不需要求精确解 第八章 流体流动