流固耦合概述及应用研究进展摘要流固耦合力学是流体力学与固体力学交叉而生成的一门力学分支。
顾名思义，它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者交互作用的一门科学。
流固耦合力学的重要特征是两相介质之间的交互作用(fluid．solid interaction)：变形固体在流体载荷作用下会产生变形或运动，而变形或运动又反过来影响流场，从而改变流体载荷的分布和大小。
总体上 ,流固耦合问题按耦合机理可分为两大类:一类的特征是流固耦合作用仅仅发生在流、固两相交界面上 ,在方程上耦合是由两相耦合面的平衡及协调关系引入的;另一类的特征是流、固两相部分或全部重叠在一起 ,耦合效用通过描述问题的微分方程来实现。
1 流固耦合概述1．1引言历史上，人们对流固耦合现象的早期认识源于飞机工程中的气动弹性问题。
Wright兄弟和其它航空先驱者都曾遇到过气动弹性问题。
直到1939年二战前夕，由于飞机工业的迅猛发展，大量出现的飞机气动弹性问题的需要，有一大批科学家和工程师投入这一问题的研究。
从而，气动弹性力学开始发展成为一门独立的力学分支。
如果将与飞机颤振密切相关的气动弹性研究作为流固耦合的第一次高潮的话，则与风激振动及化工容器密切相关的研究可作为流固耦合研究的第二次高潮。
事实上，从美国ASME应用力学部召开的历次流固耦合研讨会上可以看出，流固耦合问题涉及到很多方面。
比如：空中爆炸及响应，噪声相互作用问题，气动弹性，水弹性问题，充液结构内的爆炸分析，管道中的水锤效应，充液容器的晃动及毛细流中血细胞的变形，沉浸结构的瞬态运动，流固相互冲击，板的颤振及流体引起的振动，圆柱由于热交换引起支持附件松动的非线性流固耦合系统，声音与结构的相互作用，涡流与结构的相互作用，机械工程中的机械气动弹性问题等等。
1.2流固耦合力学定义和特点流固耦合力学是流体力学与固体力学交叉而生成的--I'l力学分支。
顾名思义，它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者交互作用的一门科学。
流固耦合力学的重要特征是两相介质之间的交互作用（fluid-solid interaction)．变形固体在流体载荷作用下会产生变形或运动，而变形或运动又反过来影响流场，从而改变流体载荷的分布和大小。
正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。
流固耦合问题可由其耦合方程来定义，这组方程的定义域同时有流体域与固体域，而未知变量含有描述流体现象的变量及描述固体现象的变量，一般而言，具有以下两点特征：a)流体域或固体域均不可能单独地求解；b)无法显式地消去描述流体运动的独立变量或描述固体运动的独立变量。
1．3流固耦合力学涉及领域及分类流固耦合问题涉及到很多方面。
比如：工程实际中所涉及到的流固耦合问题，诸如地下水抽放和油气开采所引起的地表沉降的流固耦合，水库诱发地震的流固耦合，岩坡和坝基(含深基坑开挖)的稳定性，煤层瓦斯的耦合渗流和突出，地热开发利用中的流固耦合，核废料地下储存以及石油、天然气开发过程中的流固耦合问题，储层流固耦合渗流的研究方法和理论进展，流激振动及其控制，空中爆炸及响应，噪声相互作用问题，气动弹性，水弹性问题，充液结构内的爆炸分析，管道中的水锤效应，充液容器的晃动，沉浸结构的瞬态运动，流固相互冲击，板的颤振及流体引起的振动，圆柱由于热交换引起支持附件松动的非线性流固耦合系统，声音与结构的相互作用，涡流与结构的相互作用，机械工程中的机械气动弹性问题，生物力学中血流同血管，组织液同细胞组织等耦合问题，毛细流中血细胞的变形问题，飞行器等在水中坠落涉及的耦合问题。
此外，在近代科学技术的发展中，许多涉及物相变化、化学变化的两相耦合问题，例如，宇宙飞行器在星际航行中接近周围有大气的星球时，其环境发生剧烈变化，强光、强热的影响、粒子侵蚀、表面烧蚀、催化等物理化学作用使结构的气动外形变化，出现更为复杂的气固耦合问题，高温、高热下液态金属的固化过程中，或高速冲击、穿甲过程中同一物质液态、固态同时共存，且区域不停变化，其整个过程的研究遇到的液固耦合问题等，以及有关流固耦合振动的利用问题，如旋转叶片的气弹自适应弯矩、气弹变形制动等。
显然，流固耦合作用的研究在航空、航天、水利、建筑、石油、化工、海洋以及生物领域都有着十分重要的意义。
从总体上来看，流固耦合问题按其耦合机理可分为两大类。
第一大类问题的特征是两相域部分或全部重叠在一起，难以明显地分开，使描述物理现象的方程，特别是本构方程需要针对具体的物理现象来建立，其耦合效应通过描述问题的微分方程而体现。
第二大类问题的特征是耦合作用仅仅发生在两相交界面上，在方程上耦合是由两相耦合面的平衡及协调关系引入的。
问题a)是流体与固体结构之间有大的相对运动的问题。
其典型例子是机翼颤振或悬桥振荡中发生的气固相互作用，这被人们习惯称其为气动弹性力学问题。
问题b)是具有流体有限位移的短期问题。
这类问题由引起位形变化的流体中的爆炸或冲击引起。
其特点是：人们极其关心的相互作用是在瞬间完成的，总位移是有限的，但流体的压缩性是十分重要的。
问题c)是具有流体有限位移的长期问题，如近海结构对波或地震的响应、噪声振动的响应、充液容器的液固耦合振动、船水响应等都是这类问题的典型例子。
对这类问题，人们主要关心的是耦合系统对外加动力载荷的动态响应。
1．4流固耦合问题研究现状1903年Langley的单翼机首次有动力的飞行试验时因机翼断裂而坠入Potomac河中。
10年后，Brewer的研究指出这一事故是一典型的气动弹性静扭转发散问题。
第一次世界大战初期，Handley Page轰炸机的尾翼颤振坠毁导致TLanchester以及Bairstow与Fage进行了第一批有目的的气动弹性颤振研究。
Theodorsen系统地建立了非定常气动力理论，为气动弹性不稳定及颤振机理的研究奠定了基础。
直到1939年二战前夕，由于飞机工业的迅猛发展，大量出现的飞机气动弹性问题的需要，有一大批科学家和工程师投入这一问题的研究。
从而，气动弹性力学开始发展成为一门独立的力学分支，并出现了一些在历史上有重要影响的专著及综述，其中有重要贡献的科学家如Theordosen，Garrick，Bisplinghoff,Y．C．Fung以及将随机概念引入气动弹性的Liepman，Y．K．Lin，Davenport等。
在气动弹性的研究史上，发生过另一件划时代的事件，使气动弹性问题的重要性首次明显地表现于航空技术外的领域内。
这就是1940年Tacoma 悬桥在相当低的风速下(19米/秒)振动激励而倒塌。
此后桥梁颤振问题得到了建筑工程师的高度重视。
在土木工程中，一般认为Westergaard是研究坝水耦合的创始人。
他在1933年发表了他的著名论文，给出了刚性重力坝在水平地震载荷下的动水压力分布，其解至今仍为许多国家的坝工抗震设计规范所沿用。
然而，由于Westergaard解中，假定了刚性坝并给出了地面运动的规律，因而其解本质上并未计及两相的耦合作用，只是求解了一个给定边界条件的流体动力问题。
至于对这类问题，考虑耦合作用的研究可见美加利福尼亚大学Chopra的论文。
与坝水耦合问题类似的一类液固耦合问题是水中柱体的耦联振动问题。
对这类问题，国内以著名力学家郑哲敏院士为代表的许多专家作了不少的研究，如流体引起管道振动的研究源于横垮阿拉伯输油管道振动分析的需要。
Feodosev，Housner及Niordson是早期研究简支直输液管道稳定性的学者。
应用不同的方法，他们得到了相同的基本运动方程并得出有关稳定性的相同结论。
至于悬臂输流直管的开创性研究，可见Benjamin及Gregory和Paidoussis的论文。
又有记载，早在1939年Bourrieres 就对悬臂输流管作了最早的著名研究，但由于计算工具的落后，他未能确定出不稳定的临界条件。
有关流体引起管道振动的早期专著及综述文章可见，其中给出了各类问题的详细历史评述。
船水耦合动力学的早期研究，考虑船是一刚性浮体，探讨其由于船的刚体运动引起的扰动。
Haskind对这类问题作了深入的分析，他应用Green定理构造了由于船的存在及运动引起的速度位，并推导了点源Green函数的形式。
沿着Haskind的开创性工作，人们开始致力于通过速度位表示的线性边界值问题的构造和求解，考虑船舶的存在及运动对流场的影响，进行船舶刚性摇荡运动与水相互作用的研究。
在Mays的博士论文及Newman的影响甚大的综述文章中对有关历史作了详细叙述，并建立了船舶水动力学这类问题的统一理论。
至于考虑船体弹性变形的水动弹性理论的创始性工作应归功于Bishop和他的学生Price，他们在其专著中对其研究及理论作了详述。
这一专著的出版对船水动力学的研究有深远的影响，已被译成俄文。
储液容器的振荡问题的开创性研究归功于许多俄罗斯的专家。
在Moiseev的两篇著名综述报告中，详述了有关基本理论及文献。
由于当时计算工具的落后，Moiseev的研究主要局限于刚性不动容器或作简单运动刚体(如单摆)中液体的振荡，方程的导出应用了动力学的变分原理，求解采用了一些近似方法，如Ritz 法。
对于弹性结构与液体的耦合振荡，他指出这是一相当复杂的问题，仅仅讨论了梁中有液体腔的扭转及弯曲振荡问题。
沉浸于液体中的结构对爆炸波的瞬态响应短期流固耦合问题似乎最早为Taylor于1941年所研究，他在文中讨论了水中无穷大平板对一维爆炸波的响应问题。
Taylor所讨论的问题虽然简单，但他的开创性研究给出了一个重要的一般性研究方法，这就是将总噪声速度位分解为入射、反射和幅射速度位。
另一种极为重要的短期问题是变形固体的着水撞击问题。
这一问题早期研究的开创性工作归功于VonKarman，他于1929年提出了入水理论，在两相耦合面上引入基于试验结果的附加水质量，分析研究了水机降落过程中的冲击现象。
1932年Wagner 对VonKarman的方法作了修正，引入水波影响因子，使结果更符合实际。
VonKarman和Wagner的方法虽没有对流场进行仔细分析，但基于实验的假设使问题简单明了，在一定程度上也能解释撞水现象的机理，因而这一方法仍为工程设计的初期用于估算撞水力。
至于对这类问题的全场分析，则是有限元、边界元方法充分发展后的成果。
我国的流固耦合研究起自20世纪80年代，王世忠，王茹(1992)采用三维管道固液耦合分析，推出具有阻尼机制的反对称矩阵；王世忠，刘玉兰等(1998)讨论了流体的速度，压强变化以及固-液耦合阻尼和固-液耦合振动微分方程；张智勇等(2000，2001)对充液直管管系中固．液耦合对管系轴向振动相应的影响进行了研究，推导了低频情况下的充液直管轴向，横向振动传递矩阵与弯管单元的传递矩阵；薛景宏，欧进萍等(2001)研究了在沿轴向传递的剪切波和压缩波波动作用下管道的动力响应；王世忠，于石声等(2002)用有限元法推导流体输送管道固．液耦合振动方程，讨论了流体流速，压强变化以及固．液耦合阻尼和刚度对管道固有频率的影响；张敦福，王锡平等(2004)采用直接解法求出了输送管道自由振动的固有频率和极限流速；张进国，吕英民等(2004)针对储罐与管道连接这个抗震薄弱环节，研究了地震作用下储罐与管道连接波纹管的动力响应；孙建刚，滕振超等(2004)利用Ansys有限元分析软件建立了跨越管道有限元数值分析模型，对不同地震烈度下的地震响应进行了有限元数值计算，得到不同地震作用下管道的地震动响应。