粘性不可压缩流的摄动--拉氏乘子杂交%2f混合有限元方法.pdfVIP

14—2l 4 力季与工杜一杜庆华院士八十寿辰庆贺丈集北京：清华大学出版社．1999 粘性不可压缩流的摄动一拉氏乘子 杂交／混合有限元方法 陈大鹏1柯尊平1卞学镘2 (1西南交通大学计算工程科学研究所，610031，成都 2麻省理工学院宇航系，美国) 摘要本文由杂交／混台变分原理出发．引用小参数撮动．把不可压缩约束化为轻微可压 缩条件，构造了一类9节点四边彤杂交／混合有限单元。数值实验表明，该单元性能优于 常见的各种住移型单元，包括选择性减缩积分元、加罚一选择性战缩积分单元等等，也优 干“杂交应力’单元，不出现自锁，且由于求解未知量只有一类，故计算效率高．同时， 由于引用了多变量列式。敌应力解和速度解皆具有高精度． 关蕾词杂交／混合变分原理i杂交／混舍有限单元 一、引言 的有限元分析工作，一直是具有多种背景的研究工作者集中探讨的问题之一。由于引 用了不同的方法来处理不可压缩约束，出现了多种列式方法。 流函数和流函数一涡度列式【l删通过引入流函数而恒等满足不可压缩条件，故不存 在与连续方程相关的问题。然而，此类方法所固有的缺陷，却对边界条件的实施工作 造成了障碍，并且使向三维问题的延拓工作遇到极大困难。 速度一压力列式吲【4】采用拉氏乘子法，以静水压力作为拉氏乘子，将不可压缩约束 引入列式工作。人们较早地认识到，速度和压力近似试解必须满足一定的一致性条件， 即BB条件口Ⅱ6】，才能使解答收敛、稳定。不幸的是，许多简单自然的单元，如采用等阶 插值单元，并不能满足BB条件，从而导致伪压力模式m．尽管也有一些收敛的单元，但 它们在实施时并不方便。 减缩积分一罚函数法口嗍，采用外罚函数法来处理不可压缩约束，并引用减缩积分来 计算罚项，以求避免“自锁”。该方法的突出优点是能在单元一级消去压力，从而使动 量平衡与连续方程脱离耦合，使有限元方程系统得以缩减。为了保证解答收敛和稳定， 速度与相关的压力空间仍需满足BB条件。根据文[10～12】，许多单元是不稳定的。为 了保证稳定性。人们不得不进一步降低减缩积分的阶次。由此招致过分的数值积分误 差，又破坏了解答的精度和收敛速率。 杂交有限元法【13】【I”，同样以压力作为拉氏乘子，借以引入不可压缩约束。与此同时， 陈大鹏等：粘性不可压缩流的援动一拉氏来子挛变，混合有限元方法 15 又令偏应力张量和静水压力“事先”满足线动量平衡约束。研究工作表明，该方法能 够提供良好的速度解和压力解。然而，由于不可压缩约束的非耦合特性，仍然无法避 免节点速度和(在单元内任意变化的)压力场的常数项同时作为方程系统的求解未知 量的局面， 本文采用摄动连续方程代替非耦合的不可压缩约束，建议了一种用于分析slokes 流动的罚一杂交变分原理。其中，偏应力张量和压力满足线动量平衡约束。在由此建 立的有限元模型中，完全可以在单元一级消去压力，仅以节点速度作为方程系统的求 解未知量。推导了几种9节点四边形单元，并通过数值算例，检验了单元特性。 二、Stokes流动与分析 令Q表示∥的一个有界区域，d是空间维数。删的边界由a．0表示，并假设是分片 光滑的。进一步假设将a口分成两个不相交的子区域墨和s。。stokes流动方程与边界 条件可写为： E力 一·V一即+PF=D (1) V·Ⅳ=O E妇 (2) D=(HV+VⅣ)／2E力 (3) 盯=一∥+盯’， 盯’=2，Ⅲ