流固耦合是流体力学与固体力学交叉而生成的一门独立的力学分支，它的研究对象是固体在流场作用下的各种行为以及固体变形或运动对流场影响。流固耦合力学的重要特征是两相介质之间的交互作用，固体在流体动载荷作用下会产生变形或运动，而固体的变形或运动又反过来影响流场，从而改变流体载荷的分布和大小，正是这种相互作用，将在不同条件下产生形形色色的流固耦合现象。

从20世纪80年代以来，流固耦合的研究便一直受到世界学术界的广泛关注，近年来流固耦合研究发展的3个标志为：

• 由线性流固耦合问题发展到非线性流固耦合问题；

• 由固体结构的变形和强度问题发展到固体的屈曲问题；

• 计算格式从单纯的固体有限元格式或流体的差分格式到混合或兼容的流固格式。

现已能在固体和结构中考虑材料非线性和几何非线性，在流体中也开始考虑有粘性和空化等效应的流体模型，从而得以模拟出晃动、空化、飞溅等流固耦合行为。在流体激发振动中也已经开始考虑复杂的结构阵列和流体流动，使其更加接近真实情况，从而可以更好的应用于实际情况中。

单步耦合应用频域法假设结构体以一个已知的频率和幅值进行运动，然后求解非定常气动力做功来判断稳定性。单步耦合往往需要先求解结构体的变形，然后通过将结构体的变形作用于流场，进而计算系统的阻尼和稳定性。单步耦合中对流场的求解经历了从线性到非线性的发展过程。

He先利用ANSYS计算出叶片模态振型，然后将模态振型以一定幅值耦合到流体边界，求解气动功和气动阻尼。北京航空航天大学张小伟等利用ANSYS计算了NASA67的弯曲振动阶模态，然后在流场中给定叶片振幅计算了气动力和气动阻尼。张正秋、邹正平等也利用单步耦合方法对叶轮机颤振预测和稳定性分析作了讨论。单步耦合研究叶片结构的稳定性，没有考虑到结构体和流体的相互作用，因此需要加以改进。

多步耦合方法与单步耦合方法相同之处在于，都需要对结构体和流体场进行分别求解；不同之处在于，单步耦合仅进行了一次数据交互，而多步耦合需要在多个时间点上进行交互计算，即每一次计算完成之后都需要在流体和结构体的交界面上进行载荷和位移等参数的传递。

多步耦合法的难点在于进行时间离散之后，结构体和流体场之间的数据交互总是存在滞后。Volker Carsterns介绍了多步耦合中使用的常规交错迭代法及其改进方法；S.Piperno对带预估的交错迭代方法进行了介绍；M.Sadeghi开发出叶栅颤振的多步耦合程序，研究了不同的数据传递方法在计算中的应用。西北工业大学徐敏等针对柔性大展弦比机翼发展了一种CFD/CSD的多步耦合方法。南京航空航天大学郭同庆、陆志良等用二级精度的龙格-库塔时间推进对结构运动方程进行求解，用非定常欧拉方程双时间有限体积推进对气动力进行求解，用多步耦合的方法计算了机翼的静气弹特性。

直接耦合法又称为整体积分法，该方法对结构体和流体场用统一的方程进行描述，按照统一的数值方法进行离散求解，从而在时间上实现了同步，不存在滞后现象。Bendiksen用一种混合欧拉-拉格朗日方程对流固耦合系统进行了求解，在耦合边界面实现了欧拉格式向拉格朗日格式的转换；Ge-Cheng Zha等利用直接耦合法对失谐叶盘进行了高周疲劳预测分析。由于直接耦合法涉及到不同模型和求解方法的转换，理论尚未完全成熟，开展的应用较少，国内尚处于起步阶段。

比较三种耦合方法可知，单步耦合法计算量较小，能较快得到结果，但因为没有考虑后续时间里流场对结构体的反作用，不能反映两种介质之间的能量传递。直接耦合方法准确直观，但是还需深入研究。多步耦合在目前的条件下比较容易开展研究。

在流固耦合问题的计算中，各国学者提出了不少的方法，经过归纳终结，基本可以概括为以下两个方面：

• 一类是结构部分和流体部分都按有限元法进行离散，建立流体与固体耦合的振动方程式；

• 另一类是结构部分仍按有限元法进行离散，而流体部分用边界元法离散，所谓边界元法与有限元法在连续体域内划分单元的基本思想不同，边界元法是指在定义域的边界上划分单元，用满足控制议程的函数去逼近边界条件。

边界元法与有限元相比具有单元的未知数少，数据准备简单等优点。然后建立流固耦合振动方程式。应用流体有限元和结构有限元结合的方法可以计算流体对复杂形状结构的影响，但这一方法一般要求电子计算机有较大的容量，并且计算机时较长，这给实际计算带来困难。边界元方法只对边界积分方程离散求解，计算量相对较小，在工程中得到广泛应用。

并且，当流体为无限域情况时，有限元法及差分法就显得力不从心，解决这种无限域困难的方法之一是Bettes等提出的无限元法。Bettes，Orsero等都用有限元法和无限元法结合起来处理流固耦合问题中的无限流场，但由于解的稳定性和衰减长度的不确定性，限制了无限元法的进一步应用。

相对而言，边界元法能十分有效地处理流体水动力计算，特别是在处理无限域流场时，更是得天独厚。大量学者在该领域进行了深入的研究，沈惠明、赵德有结合流体边界元和结构有限元求解流固耦合问题，采用迭代法求解流固耦合振动的特征，为了使迭代迅速、波动小，用结构在空气的振动模态（干模态）作为初始迭代向量，经过若干次迭代收敛于湿模态。安泽幸隆等人将结构部分用有限元离散，流体部分采用边界元，同时对结构和流体相互作用的界面模型做出假设，计算结果证明假设是合理的。

1. 首先要做好固体.CDB文件和流场.CAS文件，这个在HyperMesh里面可以分别导出。流体部分采用HyperMesh9.0分网，按照流体分网步骤即可，没有特殊要求。HyperMesh9.0划分固体网格。设置边界条件，载荷选项，求解控制，导出.cdb文件；

2. 建流体材料，设置属性；

3. 添加边界条件，与网格中的边界相对应；

一般来说，CAE分析工程师80%的时间都花在了有限元模型的建立、修改和网格划分上，而真正的分析求解时间是消耗在计算机工作站上，所以采用一个功能强大，使用方便灵活，并能够与众多CAD系统和有限元求解器进行方便的数据交换的有限元前后处理工具，对于提高有限元分析工作的质量和效率具有十分重要的意义。

下面就提供了一些常见的前处理器软件，下表中显示的是一些常用的前处理器软件及它们各自的工作环境、特点、优缺点等。

表1 常用前处理器软件

在流固耦合分析中，流体部分网格划分的好坏对分析结果的准确性至关重要，同时也是相对复杂的部分，因此选择适当的网格划分软件十分必要，Gambit、HyperMesh都是目前应用最为广泛的软件，它们在复杂结构上具有强大的网格效率与准确性。

也是当今最为可靠的结构非线性、流固耦合计算系统。ADINA-2F中使用的程序是基于有限元和有限体积离散图，带有非常全面和高效的解决方法，可解决任意几何学中的全部流动问题。一旦计算区域的任何一部分发生变形，对流体的Eulerian描述就不再可用了。因此，ADINA求解流体的控制方程使用Arbitrary-Lagrangian-Eulerian(ALE)

此外，可以针对流固耦合数学模型以及其有限元数值模型等，在这些模型的基础上，应用目前的编程软件，如C++、Matlab等，开发基于流固耦合下的航空发动机叶片颤振的数值模拟软件等，加入非线性的影响，使之尽量贴近实际情况，并最终应用于工程的具体问题中。总之，目前国内流固耦合的技术还处于初级阶段，还需要不断地完善与发展。

《流固耦合教学》由会员分享，可在线阅读，更多相关《流固耦合教学（20页珍藏版）》请在人人文库网上搜索。

4，调为正视，Generate一下新建一个Sketch：单击，显示 ，在此Sketch中建立液舱模型 草图单击约束（Constrains），将草图中的“水平线”调整为水平，“垂直线”调整为垂直：事实上仅用Horizontal（水平）和Vert

3、ical（垂直）就OK了。以水平约束为例，先单击Horizontal，再依次单击草图中的水平线段。调整后如下图所示：定义尺寸：左下角空缺的部分是预留贴“应变片”的部分，需要单独建模单击Extrude（拉伸），设置Operation（下拉列表中改选为Add Frozen）与拉伸尺寸（0.1m）： 然后Generate一下将第三个Solid重命名为Fluid，拉伸后的效果如下： 再新建一个Sketch，显示，在空缺处画一个长方形，然后拉伸0.1m，（其中Operation属性同样选为Add Frozen），Generate一下，同样把第四个Solid重命名为Fluid 建立舱壁：在Concept中

Shell选中Fluid（内流场），将其属性改为Fluid，（Fluent中默认均为Solid）选中“内流场”，右击，选择Form New part，并重命名为Fluid再选中舱壁（Fluid Shell） 也组成一个part，并重命名为Fluid Shell到此，液舱（内流场与舱壁）就建完了，

5、然后将二者都执行Hide body（3） 在ZX Plane内建外流场选择，调整为正视，旋转坐标系先确定外部尺寸，再确定内部尺寸：外部流场关于坐标轴（横轴）对称，两边各距离横轴0.3m，前后距离纵轴距离分别为：0.3m、3.14m. 内部为船体位置，横向（船宽）为0.14m、纵向（型深）为0.11m拉伸（Extrude）一下，拉伸长度为船长，即0.4m ，其中Operation选择Add Frozen，Generate一下图中显示外流场把船体的位置给空了出来，将重命名为Out Fluid，同时将属性改为Fluid 接下来进行流场切分（Slice）：在Tools中选择Freeze，产生透明效果单

Bodies，点击一下Bodies，选中外流场，Apply一下，再Generate，此时外流场就被分为两部分。再分别以船体所在位置内侧（以逆时针）第二、第三、第四个面为Target Face，同样进行切分操作。最终结果如下所示：外流场被切分成8个体 将这8个体组成一个part，并重命名为Out Fluid

7、：将Ship Solid、Ship Shell、Fluid及Fluid Shell部分同时show body，可以看到整个模型（包括外流场、船体、液舱及内流场）至此，几何模型建立完毕，退出，在project中保存，命名为try-、 材料设置点击打开Engineering Data，单击，选择，在表中找到结构钢（Structural），右击，选择Add to B2，单击返回设置Structural Steel 1的材料参数（这里只对Density（密度）和Youngs Modulus（弹性模量）进行设置）：表中找到以下两个参数，并照此设置这代表刚性同样对Structural

Shell，把子项中全部选上，进行如下设置Hide一下Ship Shell，先划液舱的网格右击mesh，在Insert中选sizing，先选中除与“应变片”共边的其他10条边这样，网格单位长度为0.01m再选中剩下的两条边（除去了与“应变片”共部分），每条边撒9个点（由于