流固耦合问题的数值模拟方法.docxVIP

流固耦合问题的数值模拟方法 目前，流弹性力学或形状变形力学是连续介质力学的一个非常流行的研究内容。它是由水流力学和固体力学的交叉作用产生的力学分支。同液体、气体相互作用的薄壁结构,多具有壳、壁板和薄板的形式。它们既具有大的抗拉刚度,在载荷作用下也很容易发生弯曲变形,这将导致物体附近的流体速度场及压力场的改变,因而,作用它的载荷本身也要发生变化。由于在大多数情况下,弹性薄壁物体的变形为几何非线性,因此,问题多归结为对结构变形引起的流体非线性运动问题的研究。 流固耦合问题由于其交叉性质,从学科上涉及流体力学、固体力学、动力学、计算力学等学科的知识;从技术上与不同工程领域,如土木、航空航天、船舶、动力、海洋、石化、机械、核动力、地质、生物工程等均有关系。研究和应用的领域都是非常宽广的。 近几十年来,国内外学者对流固动力耦合的理论和计算方法开展了广泛研究,取得了一些成果。但是,由于流固耦合问题的复杂性,无论在理论分析还是数值计算方面,还存在许多假设,还远远没有达到理论与实践的统一,有必要深入研究。对于流固动力耦合系统的求解,比较简单的问题,可以采用解析法和半解析法,而具有复杂边界条件的实际工程问题,很难给出其解析解答。在应用有限元法解流固耦合问题中,明显地存在欧拉坐标和拉格朗日坐标在耦合界面上的变化问题。对于解决小扰动问题和非线性问题,与保持耦合界面上的协调条件和平衡条件是大不一样的。若将流固耦合问题明确分类,以便针对不同类型的问题采用不同的解决办法,无论是求解解析解还是数值解都是大有裨益的。分类的目的是合理简化运动方程和接触条件,正确选择解决问题的简便方法。 从总体上来看,流固耦合问题按其耦合机理可分为两大类。第一大类问题的特征是两相域部分或全部重叠在一起,难以明显地分开,使描述物理现象的方程,特别是本构方程需要针对具体的物理现象来建立,其耦合效应通过描述问题的微分方程而体现。第二大类问题的特征是耦合作用发生在两相交界面上,在方程上耦合是由两相耦合面的平衡及协调关系引入的。通过耦合界面,流体动力影响固体运动,而固体的运动又影响流场。在耦合界面上,流体动力及固体的运动事先都不知道,只有在系统地求解了整个耦合系统后,才可给出它们的解答,这正是相互作用的特征所在。若没有这一特征,其问题将失去耦合作用的性质。大多数情况下,弹性薄壁构件的变形为几何非线性,再加上流体方程非线性,导致在界面上强非线性。因此有必要专门讨论流固耦合界面的描述方法和其分类简化准则。本文所谈论的是产生相互作用的接触面条件问题,属于上述第二大类流固耦合问题。 流体弹性力学的分类准则是基于壳体位移的法向分量及线性元素的转角、壳体位移场的可变形性及流体速度、压力的数值估值的基础上的。在法线位移呈线性变化的剧烈弯曲情况下、在中等弯曲(包括浅拱形壳)及小弯曲变形(线性理论)的情况下,可以采用参数构成的方法进行讨论。在具体问题中,还要考虑到接触面附近的流体状态等等。所引用流体弹性问题的分类,可进一步分类细化和将边界条件表达为其它形式。以便在应用中进一步完善并建立新的数学模型,将流固耦合理论研究提高到一个更高的水平。 当前,无论是大变形流固耦合问题,还是中变形或者小变形的流固耦合问题的解析解还很少见,随着计算技术的发展,数值分析方法中的有限元法有了长足的发展。为了促进理论分析及数值计算方法更快地完善起来,对流固耦合问题的界面相互作用的描述方法进行分类、给出简化准则,可为合理简化不同状态下的流固耦合问题,采用适当的理论分析和计算方法提供依据。 本文介绍了文献给出的界面相互作用的描述方法,在引入描述弹性体变形特征的数值m、n、k和描述流体运动特征的数值λ、ν的基础上,根据诺沃日洛夫B B在非线性弹性力学中,从几何非线性方面的分类方法出发,对流体弹性力学的流固耦合问题进行分类,明确了分类简化准则。由此,可为按类别对运动学条件、动力学条件及界面上的接触条件进行相应的简化提供可靠的依据。 1 拉格朗日变换 在流体弹性力学范围内研究流体和固体的共同运动及其相互作用时,对它们的描述是建立在连续介质力学基本原理上的。连续介质力学的基本内容,是详细地研究用欧拉观点或用拉格朗日观点描述运动之方法的转换。这是因为在传统的研究方法中,解决固态变形体的力学问题时,只使用拉格朗日变数法;而在流体力学中,主要使用的是欧拉法;但在流体弹性力学中,却又出现这两种方法都应用的情况。由于在流体和弹性物体接触面这一条件的表达方法中,包含有这两种变数系统,因此,要求有比经典力学文献中更详细描述运动特点的方法及它们之间相互转换的知识。 1.1 流固耦合问题的解决方法 对于相互接触的两种介质,根据守恒原理和受力平衡,在其接触面便可以结合拉格朗日法和欧拉法建立相互作用方程。其主要有以下四种方法: 1)单一拉格朗日法。 两种相互作