流固耦合的研究与发展综述.doc

流固耦合的研究与发展综述

目 录

1.引言 - 1 -
2.流固耦合的分类与发展 - 1 -
3.流固耦合的研究方法 - 2 -
4.流固耦合计算法 - 4 -
5.软件应用方法 - 6 -
6.总结与展望 - 14 -
参考文献 - 15 -

流固耦合的研究与发展
1.引言
近来，航空航天工业在世界上发展迅速，而作为“飞机心脏”的航空发动机是限制其发展的主要因素。目前，航空发动机日益向高负荷、高效率和高可靠性的趋势发展，高负荷导致的高你压力梯度容易引起流动分离，同时随着科技的发展，航空发动机的设计使得材料越来越轻，越来越薄，这就使得发动机内部的不稳定流动对叶片的影响大大增加，成为发动机气动及结构设计要考虑的关键问题之一。而以往单单考虑气动或结构因素不能满足实际的需求，必须将气动设计和结构设计相结合，考虑其相互作用的影响，因此流固耦合的研究应运而生。
流固耦合是流体力学与固体力学交叉而生成的一门独立的力学分支,它的研究对象是固体在流场作用下的各种行为以及固体变形或运动对流场影响。流固耦合力学的重要特征是两相介质之间的交互作用,固体在流体动载荷作用下会产生变形或运动，而固体的变形或运动又反过来影响流场，从而改变流体载荷的分布和大小，正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。
2.流固耦合的分类与发展
总体上，从流固耦合的机理上可以分为两大类：第一类，耦合作用仅仅发生在两相交界面上，在方程上的耦合是由两相耦合面上的平衡及协调来引入的如气动弹性、水动弹性等；第二类，两相部分或全部重叠在一起，难以明显地分开，使描述物理现象的方程，特别是本构方程需要针对具体的物理现象来建立,其耦合效应通过描述问题的微分方程来体现。
从20世纪80年代以来, 流固耦合的研究便一直受到世界学术界的广泛关注，近年来流固耦合研究发展的3个标志为：①由线性流固耦合问题发展到非线性流固耦合问题；②由固体结构的变形和强度问题发展到固体的屈曲问题；③计算格式从单纯的固体有限元格式或流体的差分格式到混合或兼容的流固格式。现已能在固体和结构中考虑材料非线性和几何非线性；在流体中也开始考虑有粘性和空化等效应的流体模型,从而得以模拟出晃动、空化、飞溅等流固耦合行为。在流体激发振动中也已经开始考虑复杂的结构阵列和流体流动。使其更加接近真实情况，从而可以更好的应用于实际情况中。
3.流固耦合的研究方法
流固耦合的研究经历了持续的发展，按照发展的先后顺序，可以分为单步耦合、多步耦合、直接耦合三个阶段。
1) 单步耦合
单步耦合应用频域法假设结构体以一个已知的频率和幅值进行运动，然后求解非定常气动力做功来判断稳定性。单步耦合往往需要先求解结构体的变形，然后通过将结构体的变形作用于流场，进而计算系统的阻尼和稳定性。单步耦合中对流场的求解经历了从线性到非线性的发展过程。Stuart Moffatt 和 Li He先利用ANSYS计算出叶片模态振型，然后将模态振型以一定幅值耦合到流体边界，求解气动功和气动阻尼。北京航空航天大学张小伟等利用ANSYS计算了NASA67的弯曲振动阶模态，然后在流场中给定叶片振幅计算了气动力和气动阻尼。张正秋、邹正平等也利用单步耦合方法对叶轮机颤振预测和稳定性分析作了讨论。单步耦合研究叶片结构的稳定性，没有考虑到结构体和流体的相互作用，因此需要加以改进。
2) 多步耦合
多步耦合方法与单步耦合方法相同之处在于都需要对结构体和流体场进行分别求解，不同之处在于单步耦合仅进行了一次数据交互，而多步耦合需要在多个时间点上进行交互计算，即每一次计算完成之后都需要在流体和结构体的交界面上进行载荷和位移等参数的传递。多步耦合法的难点在于进行时间离散之后，结构体和流体场之间的数据交互总是存在滞后。Volker Carsterns介绍了多步耦合中使用的常规交错迭代法及其改进方法；S.Piperno对带预估的交错迭代方法进行了介绍；M.Sadeghi开发出叶栅颤振的多步耦合程序，研究了不同的数据传递方法在计算中的应用。西北工业大学徐敏等针对柔性大展弦比机翼发展了一种CFD/CSD的多步耦合方法。南京航空航天大学郭同庆、陆志良等用二级精度的龙格-库塔时间推进对结构运动方程进行求解，用非定常欧拉方程双时间有限体积推进对气动力进行求解，用多步耦合的方法计算了机翼的静气弹特性。
3) 直接耦合
直接耦合法又称为整体积分法，该方法对结构体和流体场用统一的方程进行描述，按照统一的数值方法进行离散求解，从而在时间上实现了同步，不存在滞后现象。Bendiksen用一种混合欧拉-拉格朗日方程对流固耦合系统进行了求解，在耦合边界面实现了欧拉格式向朗格朗日格式的转换；Ge-Cheng Zha等利用直接耦合法对失谐叶盘进行了高周疲劳预测分析。由于直接耦合法涉及到不同模型和求解方法的转换，理论尚未完全成熟，开展的应用较少，国内尚处于起步阶段。
比较三种耦合方法可知，单步耦合法计算量较小，能较快得到结果，但因为没有考虑后续时间里流场对结构体的反作用，不能反映两种介质之间的能量传递；直接耦合方法准确直观，但是还需深入研究；多步耦合在目前的条件下比较容易开展研究。
4.流固耦合计算法
流固耦合的数值计算问题，早期是从航空领域的气动弹性问题开始的，这也就是通过界面耦合的情况，只要满足耦合界面力平衡，界面相容就可以。求解气动弹性问题的耦合方法通常可以分为两种:强耦合和弱耦合，强耦合方法需要对CFD和CSD方程同时进行求解，弱耦合方法是模块化的形式。其耦合通过CFD网格点上的载荷转换到CSD节点上和CSD节点上的位移插值到CFD网格点上数据交换实现。在这种弱耦合方法中,CSD和CFD网格位移可保持高精度。Guruswamy通过在动网格上建立带有欧拉/纳维-斯托克斯方程模型的方法证明了弱耦合技术。Guruswamy和Byun提出了求解二维翼型的气动弹性的一种弱耦合方法。并证明了这种松耦合方法是有效和精确的。
在流固耦合问题的计算中，各国学者提出了不少的方法，经过归纳终结，基本可以概括为以下两个方面：一类是结构部分和流体部分都按有限元法进行离散，建立流体与固体耦合的振动方程式；另一类是结构部分仍按有限元法进行离散，而流体部分用边界元法离散，所谓边界元法与有限元法在连续体域内划分单元的基本思想不同，边界元法是指在定义域的边界上划分单元，用满足控制议程的函数去逼近边界条件。所以边界元法与有限元相比具有单元的未知数少，数据准备简单等优点。然后建立流固耦合振动方程式。应用流体有限元和结构有限元结合的方法可以计算流体对复杂形状结构的影响,但这一方法一般要求电子计算机有较大的容量，并且计算机时较长，这给实际计算带来困难。边界元方法只对边界积分方程离散求解，计算量相对较小，在工程中得到广泛应用。并且，当流体为无限域情况时，有限元法及差分法就显得力不从心，解决这种无限域困难的方法之一是Bettes等提出的无限元法，Bettes，Orsero等都用有限元法和无限元法结合起来处理流固耦合问题中的无限流场，但由于解的稳定性和衰减长度的不确定性，限制了无限元法的进一步应用。相对而言,边界元法能十分有效地处理流体水动力计算，特别是在处理无限域流场时，更是得天独厚。大量学者在该领域进行了深入的研究，沈惠明、赵德有结合流体边界元和结构有限元求解流固耦合问题，采用迭代法求解流固耦合振动的特征，为了使迭代迅速、波动小，用结构在空气的振动模态(干模态)作为初始迭代向量，经过若干次迭代收敛于湿模态。安泽幸隆等人将结构部分用有限元离散，流体部分采用边界元，同时对结构和流体相互作用的界面模型做出假设，计算结果证明假设是合理的。
5.软件应用方法
ANSYS是目前十分常用的典型的流固耦合分析软件。分析机理为流体与固体部分分开进行，第一个分析作为第二个分析的荷载，如果分析是完全耦合的，那么第二个分析的结果又会影响或成为第一个分析的荷载，如此将流体与固体场耦合起来。复杂的几何图形建模可以通过UG、CATIA、PROE等专业软件完成，他们与有限元分析软件都有很好的接口，可以方便的传送文件。流固耦合的软件分析大致分为以下几个步骤：
1.首先要做好固体.CDB文件和流场.CAS文件,这个在HyperMesh里面可以分别导出。流体部分采用HyperMesh9.0分网，按照流体分网步骤即可，没有特殊要求。HyperMesh9.0划分固体网格。设置边界条件，载荷选项，求解控制，导出.cdb文件
2.导入流体网格
3.设置分析类型（ANALYSIS TYPE）-ANSYS MULTIFIELD，输入固体网格文件，设置瞬态分析，时间设置
4.建流体材料，设置属性
5.设置默认域（default domain）流场的一些特性
6.添加边界条件，与网格中的边界相对应
7.初始化
8.求解控制设置
9.输出控制设置
10.监视变量设置
11.求解
12.后处理
一般来说，CAE分析工程师80%的时间都花在了有限元模型的建立、修改和网格划分上，而真正的分析求解时间是消耗在计算机工作站上，所以采用一个功能强大，使用方便灵活，并能够与众多CAD系统和有限元求解器进行方便的数据交换的有限元前后处理工具，对于提高有限元分析工作的质量和效率具有十分重要的意义。下面就提供了一些常见的前处理器软件，下表中显示的是一些常用的前处理器软件及它们各自的工作环境、特点、优缺点等。
表：常用前处理器软件
ICEM-CFD
1.主要有四个模块：Tetra（最高水平）、Hexa（方便使用）、Global（笛卡尔网格划分软件）、AutoHexa（应用不多）。
2. 拥有强大的CAD模型修复能力、自动中面抽取、独特的网格“雕塑”技术、网格编辑技术以及广泛的求解器支持能力。
3.接口多，几乎支持所有流行的CFD软件（包括CATIA, CADDS5, ICEM Surf/DDN, I-DEAS, SolidWorks, Solid Edge, Pro/ENGINEER and Unigraphics等）。
4.使用方便，一个月内可以学会，两个月就可以针对课题努力了。同时这个软件还有后处理模块Visual3，但应用相对较少。
Gridgen
1. Gridgen是专业的网格生成器，结构网格划分很好，可以生成多块结构网格、非结构网格和混合网格，可以引进CAD的输出文件作为网格生成的基础。
2.生成的网格可以输出十几种常用商业流体软件的数据格式，直接让商业流体软件使用。对用户自编的CFD软件，可选用公开格式（Generic）。
3. Gridgen网格生成主要分为传统法和各种新网格生成方法，形成了各种现代网格生成技术。传统方法的思路是由线到面、由面到体的装配式生成方法。各种新网格生成法，如推进方式可以高速的由线推出面，由面推出体。另外还采用了转动、平移、缩放、复制、投影等多种技术。
Pointwise
1.新的用户界面：新的界面给用户崭新的视觉享受，同时也充分考虑了用户的使用习惯，增添了许多的操作功能。
　　2.可靠性：无论是结构、非结构、混合网格，Pointwise使用其高质量的网格技术配合其强大的网格生成控制功能，使用户在最小的计算机资源的情况下，可以得到最精确，最可靠的网格。
3.灵活性：首先是自动化技术，可以生成和人工干预及控制下生成的同样高质量的网格。其次是可以从不太完美的CAD数据生成符合你的求解器要求的高质量网格的技术。
Gambit
1.目前最常用的CFD前处理器，ACIS内核基础上的全面三维几何建模能力，通过多种方式直接建立点、线、面、体，而且具有强大的布尔运算能力，ACIS内核已提高为ACIS R12。该功能大大领先于其它CAE软件的前处理器。
2. 可以导入PRO/E、UG、CATIA、SOLIDWORKS、ANSYS、PATRAN等大多数CAD/CAE软件所建立的几何和网格。导入过程新增自动公差修补几何功能，以保证GAMBIT与CAD软件接口的稳定性和保真性，使得几何质量高，并大大减轻工程师的工作量。
3. 强大的几何修正功能，在导入几何时会自动合并重合的点、线、面；新增几何修正工具条，在消除短边、缝合缺口、修补尖角、去除小面、去除单独辅助线和修补倒角时更加快速、自动、灵活，而且准确保证几何体的精度。
4. 居于行业领先地位的尺寸函数（Size function）功能可使用户能自主控制网格的生成过程以及在空间上的分布规律，使得网格的过渡与分布更加合理，最大限度地满足CFD分析的需要。
5.可以说GAMBIT是目前最有优势的CFD网格软件，功能十分强大，但其使用要在Exceed环境下，占用内存比较多，常常会跑死机（不是个别的问题）。
CFX-build
1. CFX是世界上唯一采用全隐式耦合算法的大型商业软件。算法上的先进性，丰富的物理模型和前后处理的完善性使ANSYS CFX在结果精确性，计算稳定性，计算速度和灵活性上都有优异的表现。
2. CFX-BUILD是一种以结构分析软件MSC/PATRAN为基础的图形处理系统，会用Patran就会它。可以直接访问各种CAD软件, 可以从任一CAD系统以IGES格式直接读入CAD图形。
3. 具有很强的操作功能，具有出色的几何造形能力，具有高度自动的曲面和体网格划分能力, 以保证生成高质量的网格。
CFD-Geom
1.传统的生成有限元网格过程乏味而且复杂，GEOM很好的提供了解决这些问题的办法，并有效的促进模型的生成和网格的建立。
2.GEOM具有大量的几何结构，丰富的网格生成方法，支持多种几何体、网格以及边界条件的输出等。
3.目前CFD-GEOM V2009增加了许多新功能，包括表面网格离散，表面三角网格推进式生成等。这些新功能大大增加了建立网格的可靠性，减小了与实际物体的差异。
HyperMesh
1. Hypermesh最著名的特点是它具有强大的、高性能的有限元网格划分前处理功能和后处理功能。
2.在处理几何模型和有限元网格的效率和质量方面，HyperMesh具有很好的速度，适应性和可定制性，并且模型规模没有软件限制，其强大的几何处理能力可以很快的读取那些结构非常复杂，规模非常大的模型数据，从而大大提高工作效率，也使得很多应用其他前后处理软件很难或者不能解决的问题变得迎刃而解。
3. Hypermesh具有很高的有限元网格划分和处理效率，可以大大提高CAE分析工程师的效率。
4. HyperMesh具有工业界主要的CAD数据格式接口。可以直接把已经生成的三维实体模型导入到Hypermesh中，而且一般导入的模型的质量都很高，基本上不需要对模型进行修复。
5.在建立和编辑模型方面，HyperMesh提供用户一整套高度先进、完善的、易于使用的工具包。对于2D和3D建模，用户可以使用各种网格生成模板以及强大的自动网格划分模块。
TrueGrid
1.用户可以完全控制网格设计，所有的网格由块结构化六面体或四边形网格构成，
2.与当前流行的模拟软件完全兼容，除了简单的生成网格，还可以进行预处理操作，生成控制参数、选项、载荷、接触面，以及条件等，还可以指定单元类型、剖面以及材料属性等。
3. TrueGrid适用于流固耦合的分析，其块结构化设计和投影方法可用于创建用于流体或结构力学分析的网格，此外，这种技术可以很轻松完美的构建结构与流体界面，这种方法可以将一个网格嵌入到其他网格中。
4. TrueGrid可以交互式或通过脚本来操作，脚本语言可用于持续地重定义模型，可使用参数、代数公式、条件语句、循环语句、数组以及用户自定义方程等。
在流固耦合分析中，流体部分网格划分的好坏对分析结果的准确性至关重要，同时也是相对复杂的部分，因此选择适当的网格划分软件十分必要，Gambit、HyperMesh都是目前应用最为广泛的软件，它们在复杂结构上具有强大的网格效率与准确性。
此外，ADINA 也是当今最为可靠的结构非线性、流固耦合计算系统。ADINA-2F中使用的程序是基于有限元和有限体积离散图，带有非常全面和高效的解决方法，可解决任意几何学中的全部流动问题。一旦计算区域的任何一部分发生变形，对流体的Eulerian描述就不再可用了。因此，ADINA求解流体的控制方程使用Arbitrary-Lagrangian-Eulerian(ALE)表示。
1.质量守恒方程

2.动量守恒方程
控制固定域上的牛顿流体流动问题的增强形式由控制方程和对应的初始边界条件组成，控制流体问题的方程是Navier-Stokes方程的ALE描述，其中是物质速度与网格速度u之差，称为相对速度。

3.能量守恒方程

直接耦合求解的办法中，流体方程和结构方程是组合起来在一个方程组(一个刚度阵)中处理的，线性化和求解使用Newton-Raphson迭代算法。迭代耦合方法比直接计算占用的内存要小，因此可以用来求解大规模问题。
6.总结与展望
本文总结了当今国内外流固耦合相关方面的研究现状与成果、发展方向以及各种研究方法等，着重介绍了目前最常用的基于计算机软件技术的流固耦合问题的求解方法。其中，对流体以及构件的网格划分部分进行了详细的说明，在流固耦合的计算中，网格的好坏对计算结果的精度极为重要，因此本文详细介绍了现阶段常用划分网格软件的使用情况、优缺点等。但由于实际状态下流固耦合的情况是十分复杂的，目前还没有很好的划分出能够十分准确的表示实际状态的网格，尤其是在流体部分，值得进一步研究。
此外，可以针对流固耦合数学模型以及其有限元数值模型等，在这些模型的基础上，应用目前的编程软件，如C++、Matlab等，开发基于流固耦合下的航空发动机叶片颤振的数值模拟软件等，加入非线性的影响，使之尽量贴近实际情况，并最终应用于工程的具体问题中。总之，目前国内流固耦合的技术还处于初级阶段，还需要不断地完善与发展。
参考文献
[1]王勖成.有限单元法.北京:清华大学出版社.2007.
[2]杨新华,陈　雷.用ANSYS软件分析流固耦合问题.科技信息.2006
[3]余　浩,曾永清.流固耦合作用综述.西部探矿.2008
[4] A Dermudez,R Rodriguez.Finite element computation of the
vibration modes of a fluid2solid system.Computer Methods in
Applied Mechanics and Engineering.1994
[5]王昊森.CFX流固双向耦合的实现.2010
[6]增强.压气机转子叶片流固耦合计算及软件集成研究.2006
[7]王福军.《计算流体动力学分析—CFD软件原理与应用》.清华大学出版社
[8]郭术义,陈举华.流固耦合应用研究进展.济南大学学报.2004
[9]邢景棠,周盛,崔尔杰.流固耦合力学概述.力学进展.1997
[10]娄 涛.基于ANSYS的流固耦合问题数值模拟.2008
[11]董平川,徐小荷,何顺利.流固耦合问题及研究进展.1999
[12]李雪梅,陈文元.流固耦合问题计算方法综述.2009
[13]戴大农,王勖成.流固耦合系统动力响应的模态分析理论.固体力学学报.1990
[14]王　超,李红云.流固耦合系统动力响应分析的精细时程积分法.上海交通大学学报.2002
[15]崔　鹏,韩景龙.一种局部形式的流固耦合界面插值方法.2009
[16]强大的网格生成工具——TrueGrid.中仿科技公司
[17]CFD-GEOM V2009.0 User Manual
[18]Xu M,Chen S L.CFD/CSD coupling calculation research. Chinese Journal of Applied Mechanics.2004
[19]徐　敏,安效民,陈士橹.一种CFD/CSD耦合计算方法.航空学报.2006
[20]Xu M,Shi Z J,Chen S L.A data exchange method for fluid-structure interaction problems.Journal of North western Polytec hnical University.2003
[21] 周盛.叶轮机械流固耦合有关模型综述.燃气涡轮试验与研究.1995
[22] 张小伟1涡轮机械叶片的流固耦合数值计算方法。航空动力学报。2009
[23]Miller G D. Active flexible wing(AFW)technology. AFWAL-TR-87-3096.1987
[24]Guruswamy G P.Coupled finite difference/finite element approach for wing body aeroelasticity. AIAA Paper 92-4680. 1992
[25]Guruswamy G P. Byun C. Fluid-structural interactions using Navier-Stokes flow equations coupled with shell finite element structure. AIAA Paper 93-3087.1993
[26] Sadeghi M,Liu F. Coupled fluid-structure simulation for
turbo machinery blade rows.AIAA Paper.2005
[27] W. Dettmer, D. Peric. A computational framework for fluid–structure interaction: Finite element formulation and applications. computer methods in applied mechanics and engineering. 2006
[28] Matthew F. Barone,Irina Kalashnikova,Matthew R.Brake, and Daniel J.Segalman.Reduced Order Modeling of Fluid-Structure Interaction.SANDIA REPORT.2009
[29] Joris Degroote,Robby Haelterman,Sebastiaan Annerel, Peter Bruggeman,Jan Vierendeels. Performance of partitioned procedures in fluid-structure interaction.Computers and Structures.2010
[30] Folco Casadei，Serguei Potapov. Permanent fluid-structure interaction with non-conforming interfaces in fast transient dynamics. computer methods in applied mechanics and engineering. 2004
[31] Irena Lasiecka , Yongjin Lua. Interface feedback control stabilization of a nonlinear fluid-structure interaction. Nonlinear Analysis.2012
[32] Patrick Le Tallec, Jean-Frédéric Gerbeau , Patrice Hauret , Marina Vidrascu . Fluid-structure interaction problems in large deformation. comptes rendus mecanique.2005
[33] Suncica Canic, Andro Mikelic, Josip Tambaca. A two-dimensional effective model describing fluid–structure interaction in blood flow:analysis, simulation and experimental validation. comptes rendus mecanique.2005
[34] R.Jaiman, P.Geubelle, E.Loth, X.Jiao. Combined interface boundary condition method for unsteady fluid–structure Interaction. computer methods in applied mechanics and engineering. 2011
[35] R.Jaiman, P.Geubelle, E.Loth, X.Jiao. Transient fluid–structure interaction with non-matching spatial
and temporal discretizations. Computers&Fluids.2011
[36]纪兵兵,陈金瓶.ANSYS ICEM CFD网格划分技术实例详解.中国水利水电出版社.2012
[37]刘志远,郑源,张文佳,司佳钧.ANSYS-CFX单向流固耦合分析的方法.水利水电工程设计.2009
[38]贺李平.ANSYS13.0与HyperMesh11.0联合仿真有限元分析.机械工业出版社.2012