流体流动问题的有限元法优秀课件.ppt

流体流动问题的有限元法第1 页，本讲稿共17 页描述稳态不可压缩流体流动的微分方程就没有与之对应的泛函表达式。而流体流动问题却是工程中经常遇到的问题（1）车辆高速运行时的气动稳定性；（2）两列高速运行的列车会车时的压力波动；（3）列车进入隧道时的压力波动；（4）建筑物的风载荷；（5）室内的通风与空调；（6）桥梁的风致振动；（7）船舶的运行阻力；（8）飞机的升力、阻力。第2 页，本讲稿共17 页不论是那种原因，如果找不到与微分方程对应的泛函表达式，那么就无法利用变分原理建立有限元的计算格式。这时我们只有寻求另外的途径。这个途径就是：加权余量法。二 加权余量法加权余量法的基本思想：通过使试探函数与真值的加权误差在求解域内的总和为零，以求得满足微分方程的近似解。设某物理问题的控制微分方程及其边界条件分别为第3 页，本讲稿共17 页 为待求函数。如果 无法或不易直接求解，可选一个试探函数式中 ci待定常数；i试探函数项。将试探函数带入控制微分方程及其边界条件，一般来讲不可能正好满足方程，在域 内和边界S 上会产生误差，即式中R 和Rb称为余量（或残数，残差，残值）。第4 页，本讲稿共17 页加权余量法的基本思想：在域 内 和/或 边界S 上寻找n 个线性无关的函数Wi(i=1,2,n)，使余量R 和Rb在加权求和的意义上等于零，即这里Wi称为权函数。加权余量法所假设的试探函数并不能满足微分方程及其边界条件，但是当加权的试探函数与真值的误差（余量）在求解域上积分为零时，那么试探函数就在总体上满足微分方程及其边界条件。当n 足够大时，试探函数就趋近于真解。第5 页，本讲稿共17 页介绍两种常用的权函数。1 最小二乘加权余量法设有满足边界条件的试探函数带入控制微分方程将产生余量如果希望余量R 在最小二乘的意义下为最小（即令R的平方和为最小），则构造第6 页，本讲稿共17 页使比较可知，权函数通过求解可求出ci，进而得到。第7 页，本讲稿共17 页2 伽辽金加权余量法如果选用试探函数中的试探函数项i作为权函数Wi，就成为伽辽金加权余量法。即在许多物理问题控制微分方程的有限元法求解过程中，都采用伽辽金加权余量法推导有限元计算格式。第8 页，本讲稿共17 页7-2 二维流体流动的有限元计算格式二维稳态可不压缩流体流动方程由连续方程和动量方程描述方程中的待求变量为流体速度u,v和压力p。第9 页，本讲稿共17 页根据有限元法的计算思路，首先选取插值函数来近似描述速度u,v和压力p 在单元内的变化情况。式中N 单元形状函数；ui,vi,pi单元节点处的速度和压力值。第10 页，本讲稿共17 页这里的插值函数，就作为加权余量法中的试探函数；其中的形状函数，就作为加权余量法中的权函数。第11 页，本讲稿共17 页经过推导和简化，可得单元方程为或其中分别为单元刚度矩阵、单元节点列向量和单元节点受到的来自“铰链”的“节点力”。第12 页，本讲稿共17 页第13 页，本讲稿共17 页将 在求解区域内分别按节点号叠加，就可以构成整个流场的有限元计算的总体方程第14 页，本讲稿共17 页7-3 流场有限元分析的几个特殊问题1 速度和压力插值函数的阶次。速度插值函数要高于压力一阶，否则方程会出现“病态”。2 主对角线元素为零。采用罚函数法，将压力用速度表示。求出速度后，再计算压力。3 刚度矩阵不对称。原来介绍的压缩存储方法全部没有用。4 非线性方程组需要迭代计算。第15 页，本讲稿共17 页小结：（1）本章讨论了利用有限元法求解非结构问题的又一个例子流体流动问题的计算。所用方法为加权余量法，通过将试探函数带入控制微分方程，基于使所产生的误差（余量）在加权平均的意义上等于零的思想，来推导该控制方程的有限元计算格式。（2）本章简要介绍了加权余量法的基本概念，最小二乘加权余量法中权函数的计算，以及伽辽金加权余量法中权函数的确定。对于无法利用变分原理，即找不到等价的泛函极值问题的控制微分方程有限元求解问题，一般来讲，都可以利用加权余量法推导其有限元计算格式。第16 页，本讲稿共17 页（3）本章简述了二维稳定流场问题的有限元计算格式的推导思路，同时涉及了流场有限元计算的一些特殊问题，计有 流场速度插值函数与压力插值函数不同阶问题；流场有限元方程系数矩阵的“病态”问题；流场有限元方程系数矩阵的非对称性问题；流场有限元方程的非线性问题。第17 页，本讲稿共17 页