下载文档原格式
流固耦合定义：它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者相互作用的一门科学。
流固耦合力学的重要特征是两相介质之间的相互作用,变形固体在流体载荷作用下会产生变形或运动。
变形或运动又反过来影响流，从而改变流体载荷的分布和大小，正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。
(一)流固耦合动力学：求解方法与基本理论---张阿漫，戴绍仕●有限元法●边界元法●SPH法与谱单元法●瞬态载荷作用下流固耦合分析方法●小尺度物体的流固耦合振动●水下气泡与边界的耦合效应按耦合机理分两大类：1 耦合作用只发生在两相交界面---界面耦合（场间不相互重叠与渗透），耦合作用通过界面力（包括多相流的相间作用力等）起作用。
它的计算只要满足耦合界面力平衡，界面相容就可以了（其耦合效应是通过在方程中引入两相耦合面边界条件的平衡及协调关系来实现的）。
如气动弹性，水动弹性等。
按照两相间相对运动的大小及相互作用分为三类：(1)流体和固体结构之间有大的相对运动问题"最典型的例子是飞机机翼颤振和悬索桥振荡中存在的气固相互作用问题,一般习惯称为气动弹性力学问题"(2)具有流体有限位移的短期问题"这类问题由引起位形变化的流体中的爆炸或冲击引起"其特点是:我们极其关心的相互作用是在瞬间完成的,总位移是有限的,但流体的压缩性是十分重要的"(3)具有流体有限位移的长期问题"如近海结构对波或地震的响应!噪声振动的响应!充液容器的液固耦合振动!船水响应等都是这类问题的典型例子"对这类问题,主要关心的是耦合系统对外加动力荷载的动态响应"2 两域部分或全部重叠在一起，难以明显的分开，使描述物理现象的方程，特别是本构方程需要针对具体的物理现象来建立，其耦合效应应通过建立与不同单相介质的本构方程等微分方程来体现。
按耦合求解方法分两大类：1 直接耦合求解：直接耦合是在一个求解器中同时求解不同物理场的所有变量，需要针对具体的物理现象来建立本构方程，其耦合效应通过描述问题的微分方程来体现。
关键词流固耦合; 气动弹性; 水动弹性; 非线性动力学; 计算力学1 定义和特点流固耦合力学是流体力学与固体力学交叉而生成的一门力学分支. 顾名思义, 它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者交互作用的一门科学. 流固耦合力学的重要特征是两相介质之间的交互作用(f lu id2so lidin teract ion) : 变形固体在流体载荷作用下会产生变形或运动, 而变形或运动又反过来影响流场, 从而改变流体载荷的分布和大小. 正是这种相互作用将在不同条件下产生形形色色的流固耦合现象.流固耦合问题可由其耦合方程来定义[ 1 ]. 这组方程的定义域同时有流体域与固体域, 而未知变量含有描述流体现象的变量及描述固体现象的变量, 一般而言, 具有以下两点特征:a) 流体域或固体域均不可能单独地求解;b) 无法显式地消去描述流体运动的独立变量或描述固体运动的独立变量.从总体上来看, 流固耦合问题按其耦合机理可分为两大类. 第一大类问题的特征是两相域部分或全部重叠在一起, 难以明显地分开, 使描述物理现象的方程, 特别是本构方程需要针对具体的物理现象来建立, 其耦合效应通过描述问题的微分方程而体现. 图1给出的渗流问题是这类问题的典型例子, 描述其现象的微分方程如下[ 2 ]:其中, ý表示梯度算子, u 表示土壤骨架的位移矢量, p 为渗流压力, R 是应力张量(用矢量式) , L 是相应于应变的微分算子, D 是弹性矩阵, b 是体力矢量, k 是渗透率, K f 是流体的体积模量, n 是空隙率, m = [1, 1, 1, 0, 0, 0 ]T . 这里, 由于耦合效应, 固体的本构关系中出现了压力项.土壤渗流相互作用第二大类问题的特征是耦合作用仅仅发生在两相交界面上, 在方程上耦合是由两相耦合面的平衡及协调关系引入的. 本文中, 我们主要讨论这一类问题.对于第二大类问题, Zienk iew icz与其合作者Bet tess 在文[3 ] 中按两相间相对运动的大小及相互作用性质将其分为三小类. 图2中示出了这三种问题.问题a) 是流体与固体结构之间有大的相对运动的问题. 其典型例子是机翼颤振或悬桥振荡中发生的气固相互作用, 这被人们习惯称其为气动弹性力学问题. 在这类问题中的基本物理关系和物理过程可用易于理解的所谓方块图加以描述, 这种方法由著名力学家冯元桢(Y. C. Fung) 教授[ 425 ]引用到气动弹性力学中来, 特别是对于气动弹性稳定问题中的反馈过程, 用这种方法说明是很有启发性的. 图3示出了机翼颤振的这一方框图. 图中三个方框表示了机翼(结构) 在这类问题中执行的三种不同功能: 首先它产生空气动力, 其次是产生惯性力, 再就是它产生弹性变形. 机翼按空气动力学规律产生升力A , 而机翼振动时则引起惯性力I . 这两种力A + I 使弹性机翼产生变形H, 从而又产生新的作用力A 和I , 这样, 以反馈过程的形式构成一条闭合回路, 如果出现变形的振幅随时间不断增大的现象, 则称为颤振.图4 流固耦合问题中各种力之间的相互关系图问题b) 是具有流体有限位移的短期问题. 这类问题由引起位形变化的流体中的爆炸或冲击引起.其特点是: 人们极其关心的相互作用是在瞬间完成的, 总位移是有限的, 但流体的压缩性是十分重要的.问题c) 是具有流体有限位移的长期问题, 如近海结构对波或地震的响应、噪声振动的响应、充液容器的液固耦合振动、船水响应等都是这类问题的典型例子. 对这类问题, 人们主要关心的是耦合系统对外加动力载荷的动态响应.图4 中示出了流固耦合中各种力之间的相互影响关系. 其中, 两个虚线描绘的大圆周分别划出了流体与固体. 在这两个圆周相切的地方, 用一个小圆表示了两相耦合界面. 通过耦合界面, 流体动力影响固体运动, 而固体的运动又影响流场. 在耦合界面上, 流体动力及固体的运动事先都不知道, 只有在系统地求解了整个耦合系统后, 才可给出它们的解答, 这正是相互作用的特征所在. 若没有这一特征, 其问题将失去耦合作用的性质. 例如, 若给定流固交界面上的流体动力或交界面上固体结构的运动规律, 耦合机理将会消失, 原来的耦合系统将被解耦而成为单一固体在给定表面力下的动力问题及单一流体在给定边界条件下的流体力学边值或初边值问题.在最一般情况下, 流体与固体通过两相交界面的相互作用同时受流体及固体各自的弹性力和惯性力影响, 这_______就是两个大圆周中间方框中表示的一般流固耦合问题. 随着研究问题的目的不同, 可将着眼点放在流场或固体结构上进行研究. 流体力学工作者多着眼于流场, 而固体力学工作者则注重结构. 在工程实际问题中, 可针对不同性质的问题, 作相应的简化, 从而便有简化后的耦合问题. 例如, 研究水同结构相互作用的非短期问题时, 水的可压性可以不计, 这就构成不可压流体同固体的耦合问题. 类似地, 若忽略结构的弹性变形, 就有刚体同流体的相互作用问题. 在航空中, 独成一个学科的刚体飞机飞行力学问题就是重要的例子.也可以在某些问题中忽略流体或固体的惯性效应, 从而有忽略流体惯性的耦合问题及忽略固体惯性的耦合问题. 在空气弹性力学中的静力发散, 舵面效率等问题即是重要的忽略结构惯性的流固耦合问题. 至于忽略流体惯性时的耦合问题, 其本质就是将流体(通常为气体) 视为一弹簧, 如空气弹簧, 这在工程中也常常见到. 所有这些简化后的耦合问题, 包括非耦合性质的可压流体动力学及变形固体动力学问题, 在图4中用虚线圆周上的方框表示出来. 于是,每种流固耦合问题可以按该问题中诸力所处的相互关系而进行直观的区分.2 发展简史流固耦合问题由于其交叉性质, 从学科上涉及流体力学、固体力学、动力学、计算力学等学科的知识; 从技术上与不同工程领域, 如土木、航空航天、船舶、动力、海洋、石化、机械、核动力、地震地质、生物工程等均有关系. 其研究问题甚广, 难以确定合适的研究分类,而且随着科学技术的发展, 其分类也在不停的变化, 这里以美国机械工程师学会(A SM E) 出。
流固耦合分析（FSI）流固耦合分析（FSI）是涉及流体和固体之间相互作用的问题研究，其理论包括了几个主要方面：流体力学、固体力学、耦合边界条件、求解器等。
以下是流固耦合分析的详细理论讲解，带有相关公式和尽量详细的说明。
一、流体力学1. 守恒定律质量守恒定律：$$ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{u}) = 0 $$动量守恒定律：$$ \rho \frac{\partial \mathbf{u}}{\partial t} + \rho (\mathbf{u} \cdot \nabla) \mathbf{u} = \nabla \cdot \tau + \mathbf{f} $$其中，$\rho$是流体密度，$\mathbf{u}$是流体速度，$\tau$是应力张量，$\mathbf{f}$是体力。
2. 纳维-斯托克斯方程$$ \rho \frac{\partial \mathbf{u}}{\partial t} + \rho (\mathbf{u} \cdot \nabla) \mathbf{u} = \nabla \cdot (-p\mathbf{I} + \tau) + \mathbf{f} $$其中，$p$是静压力，$\mathbf{I}$是单位张量。
3. 边界条件（1）速度边界条件：$\mathbf{u} = \mathbf{u}_b$，其中$\mathbf{u}_b$是边界上的速度。
（2）压力边界条件：$p = p_b$，其中$p_b$是边界上的压力。
4. 流体力学求解器常用的流体力学求解器有OpenFOAM、ANSYS Fluent等。
二、固体力学1. 力学基本方程$$ \tau = \sigma\cdot \mathbf{n} $$其中，$\tau$是表面上的接触力，$\sigma$是固体的应力张量，$\mathbf{n}$是表面的单位法向量。
1、打开ANSYS Workbench, 拖动各模块到空白区，并照此连接各模块。
2 2、打开第一个模块当中的Geometry，建立几何模型：（1）在XY Plane内建立Ship Shell船长：0.4、船宽：0.14、型深0.11将第一个Solid重命名为Ship Solid在Concept中选择Surfaces From Faces，选中模型的六个面，然后Apply、Generate。
重命名第二个Ship Solid为Ship Shell右击Ship Solid, 选择Hide Body，显示Ship Shell, 然后对Ship Shell执行同样操作（即隐去）（2）在YZ Plane内建立液舱单击（New Plane），选择YZ plane，，Apply一下将YZ Plane 向X正方（图中为法向，即Z）向偏移0.02mGenerate一下，然后Show body 一下Ship Solid 与Ship Shell可以看到YZ Plane已平移到Body内了再将Ship Solid 与Ship Shell 都Hide，选择Plane 4，调为正视，Generate一下新建一个Sketch：单击，显示，在此Sketch中建立液舱模型草图单击约束（Constrains），将草图中的“水平线”调整为水平，“垂直线”调整为垂直：事实上仅用Horizontal（水平）和Vertical（垂直）就OK了。
以水平约束为例，先单击Horizontal，再依次单击草图中的水平线段。
调整后如下图所示：定义尺寸：左下角空缺的部分是预留贴“应变片”的部分，需要单独建模单击Extrude（拉伸），设置Operation（下拉列表中改选为Add Frozen）与拉伸尺寸（0.1m）：然后Generate一下将第三个Solid重命名为Fluid，拉伸后的效果如下：再新建一个Sketch，显示，在空缺处画一个长方形，然后拉伸0.1m，（其中Operation属性同样选为Add Frozen），Generate一下，同样把第四个Solid重命名为Fluid建立舱壁：在Concept中选择Surfaces From Faces，选中除“应变片”外的其余9个面图中“应变片”显示为未着色，即不选中，然后Generate一下将第五个Fluid重命名为Fluid Shell再Surfaces From Faces一次，选中“应变片”，Apply，Generate，同样将其重命名为Fluid Shell选中Fluid（内流场），将其属性改为Fluid，（Fluent中默认均为Solid）选中“内流场”，右击，选择Form New part，并重命名为Fluid再选中舱壁（Fluid Shell）也组成一个part，并重命名为Fluid Shell到此，液舱（内流场与舱壁）就建完了，然后将二者都执行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，产生透明效果单击Slice（或者在Create中单击Slice），在Slice Type中选择Slice by surface，点击Target Face，选中船体所在位置（即图中外流场所空出来的位置）内侧某一个面（以左侧面为第一个面为例），Appy一下。
