有限元流体流动问题的有限元法.pptx

《有限元流体流动问题的有限元法.pptx》由会员分享，可在线阅读，更多相关《有限元流体流动问题的有限元法.pptx（17页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、描述稳态不可压缩流体流动的微分方程就没有与之对应的泛函表达式。而流体流动问题却是工程中经常遇到的问题（1）车辆高速运行时的气动稳定性；（2）两列高速运行的列车会车时的压力波动；（3）列车进入隧道时的压力波动；（4）建筑物的风载荷；（5）室内的通风与空调；（6）桥梁的风致振动；（7）船舶的运行阻力；（8）飞机的升力、阻力。第1页/共17页不论是那种原因，如果找不到与微分方程对应的泛函表达式，那么就无法利用变分原理建立有限元的计算格式。这时我们只有寻求另外的途径。这个途径就是：加权余量法。二 加权余量法加权余量法的基本思想：通过使试探函数与真值的加权误差在求解域内的总和为零，以求得满足微分方程的近

2、似解。设某物理问题的控制微分方程及其边界条件分别为第2页/共17页为待求函数。如果无法或不易直接求解，可选一个试探函数式中 ci待定常数；i试探函数项。将试探函数带入控制微分方程及其边界条件，一般来讲不可能正好满足方程，在域内和边界S上会产生误差，即式中R和Rb称为余量（或残数，残差，残值）。第3页/共17页加权余量法的基本思想：在域内 和/或 边界S上寻找n个线性无关的函数Wi(i=1,2,n)，使余量R和Rb在加权求和的意义上等于零，即这里Wi称为权函数。加权余量法所假设的试探函数并不能满足微分方程及其边界条件，但是当加权的试探函数与真值的误差（余量）在求解域上积分为零时，那么试探函数就在

3、总体上满足微分方程及其边界条件。当n足够大时，试探函数就趋近于真解。第4页/共17页介绍两种常用的权函数。1 最小二乘加权余量法设有满足边界条件的试探函数带入控制微分方程将产生余量如果希望余量R在最小二乘的意义下为最小（即令R的平方和为最小），则构造第5页/共17页使比较可知，权函数通过求解可求出ci，进而得到 。第6页/共17页2 伽辽金加权余量法如果选用试探函数中的试探函数项i作为权函数Wi，就成为伽辽金加权余量法。即在许多物理问题控制微分方程的有限元法求解过程中，都采用伽辽金加权余量法推导有限元计算格式。第7页/共17页7-2 二维流体流动的有限元计算格式二维稳态可不压缩流体流动方程由连

4、续方程和动量方程描述方程中的待求变量为流体速度u,v和压力p。第8页/共17页根据有限元法的计算思路，首先选取插值函数来近似描述速度u,v和压力p在单元内的变化情况。式中N单元形状函数；ui,vi,pi单元节点处的速度和压力值。第9页/共17页这里的插值函数，就作为加权余量法中的试探函数；其中的形状函数，就作为加权余量法中的权函数。第10页/共17页经过推导和简化，可得单元方程为或其中分别为单元刚度矩阵、单元节点列向量和单元节点受到的来自“铰链”的“节点力”。第11页/共17页第12页/共17页将 在求解区域内分别按节点号叠加，就可以构成整个流场的有限元计算的总体方程第13页/共17页7-3

5、流场有限元分析的几个特殊问题1 速度和压力插值函数的阶次。速度插值函数要高于压力一阶，否则方程会出现“病态”。2 主对角线元素为零。采用罚函数法，将压力用速度表示。求出速度后，再计算压力。3 刚度矩阵不对称。原来介绍的压缩存储方法全部没有用。4 非线性方程组需要迭代计算。第14页/共17页小结：（1）本章讨论了利用有限元法求解非结构问题的又一个例子流体流动问题的计算。所用方法为加权余量法，通过将试探函数带入控制微分方程，基于使所产生的误差（余量）在加权平均的意义上等于零的思想，来推导该控制方程的有限元计算格式。（2）本章简要介绍了加权余量法的基本概念，最小二乘加权余量法中权函数的计算，以及伽辽金加权余量法中权函数的确定。对于无法利用变分原理，即找不到等价的泛函极值问题的控制微分方程有限元求解问题，一般来讲，都可以利用加权余量法推导其有限元计算格式。第15页/共17页（3）本章简述了二维稳定流场问题的有限元计算格式的推导思路，同时涉及了流场有限元计算的一些特殊问题，计有流场速度插值函数与压力插值函数不同阶问题；流场有限元方程系数矩阵的“病态”问题；流场有限元方程系数矩阵的非对称性问题；流场有限元方程的非线性问题。第16页/共17页感谢您的观看！第17页/共17页