流固耦合问题的数值分析研究.pdfVIP

第八届全国现代结构工程学术研讨会 流固耦合问题的数值分析 陈进于钱若军 (同济大学建筑工程系。上海200092) 提要：在本文中采用有限单元法来模拟三维流固耦合问题。采用二阶显式Taylor-Galerkin方法和四面体单元分析高压缩流体流 动．引入ALE表达，对不可压缩流体采用速度修正的分形方法．固体采用三节点三角单元，其中每个节点包含6个自由度 (三个线位移向量和三个角位移向量)。采用Nemark方法对时间进行积分并且对动态平衡方程采用修正的Lagrangean表 达方法．采用共轭梯度法求解代数方程，采用限制线位移和角位移的递增迭代方法来解决非线性问题。对节点采用向量处 理进行最优化。 关键词#流固耦合，ALE．网格运动，Taylcr--Galerkin方法 一、引言 在土木工程之中。会遇到很多由于流体和固体相互作用引起的物理现象，例如桥梁和高耸建筑的风振响应等。 流崮耦合问题是目前广为关注和具有挑战性的课题之一，涉及多个学科和多方面的技术，解决这一问题需要对流 固耦合理论体系做广泛和深入的研究。 流固耦合全称为流体．结构(同体)的耦合作用(Huid-Slructures／SolidIntel粥ti∞)，英文简称为FSl，它是指 在流体(空气，水流)作用下，结构特别是柔性结构(如索结构，膜结构)的响应(位移，速度和加速度)对周 围流场会产生较大影响，而流场的改变会进一步改变作用在结构表面上的流体力大小，从而形成的流体与结构之 间的相互耦合作用。 流固耦合问题的研究主要有三种方法和手段。一是试验方法，如风洞试验；二是解析方法；三是数值方法． 由于在试验中刚性模型居多，所以不能完全模拟流固耦合现象，而解析方法则是对流体模型进行简化，在应用上 具有较大的局限性。而随着计算机技术的发展，流固耦合的数值模拟技术正越来越广泛地得到应用。从20世纪 70年代以来，国外很多的研究团体的流固耦合理论及流固耦合数值求解方法取得了很多成果。例如美国麻省理 Bathe教授领导开发的ADINA软件系统，其FSl求解模块可以较好地模拟流体和涉 工学院(MIT)Klaus-Jurgcn 及大变形、接触等具有非常明显非线性响应的结构之间完全耦合的物理现象。国内一些专家如沈士钊、武岳就采 用了分块迭代法来计算流固耦合作用。但与国外相比，目前国内对流固耦合理论体系及其数值模拟的研究还处于 初步阶段。由同济大学钱若军教授领导的流固耦合模拟研究小组，目前正在对流固耦合理论体系及其数值模拟做 系统的研究，探寻一种能够在单一软件平台上实现流固耦合的强、弱耦合算法。 流固耦合问题的模拟方法主要被分为两类：分离方法和集成方法。分离方法是指流体和固体的状态控制方程 分别对时间进行积分；而在集成方法中流体域和固体域考虑为单一整体，允许同时对其进行时问积分。在由 力边界条件在子域间保持协调变化。在这些方法中每一个子域都可以用离散方法和有效运算法则分别进行单独求 解。在新方法和模型可以引入模块方法以得到更大的柔性。另外每一个子域网格划分的建立可以不考虑边界面流 体和固体节点的协调问题。采用这些模数特性和柔性可以使计算变得方便。 本文介绍了一个采用分离方法来模拟流固耦合问题的运算法则。采用线性四面体有限单元的二阶显式 等人广泛使用。对固体的分析采用包含旋转自由度的广义协调三角板壳单元，考虑几何非线性，采用递增迭代方 法和共轭梯度法分别解决非线性问题和代数问题，对节点则采用向量处理。 工业建筑2008增刊 617 第八届全国现代结构工程学术研讨会 二、流体动力学求解 (一)高压缩流体的二阶显式Taylor-Galerkin方法 在ALE表达中表示质量、动量和能量守恒定律的方程可以分别表达为： 补嚣刊笔．(f都’3) 其中jI=兄+％且。 p p’)I O 肌 p_Iivt+p4l —气