流固耦合理论研究述评3—流固耦合问题求解策略.pdfVIP

第七届全国现代结构工程学术研讨会 流固耦合理论研究述评(3) 一流固耦合问题求解策略 苏波1袁行飞2聂国隽3钱若军1 1同济太学土木工程学院．上海200092)(2浙江大学土术工程学院，杭州土木上程学院)(3同济大学航空航天与力学学院，上海2000鸵) 摘要：流固耦合(FsI)是目前很多领域研究的热点和也是难点之一。目前．国内研究对此相对滞后或比较片面。为在同一 软件平台上实现多种耦合算法，本文对流固耦合理论进行了研究。内容包括FsI研究进展；流体力学基奉理论和流体计算力 学力学(cFD)有限元法：耦合界面条件、边界追踪、网格更新、以及大型非线性方程的求解等诸多问题；特别是对目前广为 关注和讨论的强耦合和弱耦合算法作了理论上的探讨和区分。本文为流固耦合理论数值模拟提供理论和方法的指导。 关键词：流固耦台，强耦合．弱耦台，边界条件，^LE法 七、流固耦合求解策略 从本质上讲，流固耦合问题的求解，应该根据统一的耦合方程，在一个时间步内对流体域和固体域中 所有的未知量同时求解，即所隋强耦台d但由于目前的计算机计算能力和问题的求解精度不同，可以有相 弱耦合从理论上作了区分，并对比了不同方法的计算效率。在此基础上本文通从单一物理方程的出发．阐 明了强耦合、弱耦合及单向耦合的本质，这为我们在统一软件平台基础上，实现多种算法提供了可能。 (一)耦合方程的建立及其凝聚 l、统一耦台方程 采用有限元法来求解流固耦合问题，即对流体域、结构域均采用有限元方法来进行离散，边界条件 在单元的组合过程中强迫实现。假定△X’为耦合系统的变量增量，根据节点所在的物理域的不同，可 将△XA写成三部分： △x； △X扣=△x窖 (7．1) △X： 式中，△x；为流体域内的节点未知数向量。△x：为流固耦合边界上的节点来知数向量．△x；为固体域 内的节点未知数向量。下角标，为表示物理域内部，C为流固耦合边界。 线性方程化后的耦合方程可以采用如下简单形式表示： A矗△X^=R扫 (7．2) 式中，A’为耦合系统的等效质量矩阵，R。为耦合系统的外力矢量。上式可展开为如下形式： 工业建筑2007增刊 717 第七届全国现代结构工程学术研讨会 《 』4丢 o1f艋j1 《l《c+A‰ (7．3) 铭}{船于}亍 O 砭 筋j【继jJ ，●●●●●●●●●●●●●●●●l 吖硭彤 、●●●●●●●L，●●●●●J 式中，下标C和f分别表示为耦台界面处和内部节点上的变量。 2、统一耦合方程的凝聚 公式(7．1)是一个统一的耦合方程，我们可以两种方法对之进行静力凝聚： 1)消除流体内部变量，结合界面对结构进行求解．过程如下： a结合界面的结构求解： ． (7．4) I如+A舌[《r《+屹l缸兽+站△x；=硭一心[A玎1吖 健。丝：+雀l丝≯毯