流固耦合原理52457ppt课件

流固耦合定义：它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者相互作用的一门科学。

流固耦合力学的重要特征是两相介质之间的相互作用,变形固体在流体载荷作用下会产生变形或运动。

变形或运动又反过来影响流，从而改变流体载荷的分布和大小，正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。

(一)流固耦合动力学：求解方法与基本理论---张阿漫，戴绍仕●有限元法●边界元法●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) 出。

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一下。

流固耦合教材流固耦合是流体力学和固体力学相互作用的一种重要现象，广泛应用于工程领域中。

本教材旨在介绍流固耦合的基本概念、理论模型和数值计算方法，以及其在工程实践中的应用。

流体力学和固体力学是研究流体流动和固体结构力学特性的两个独立的学科。

然而，在实际工程中，往往需要考虑流体和结构之间的相互影响，这时就需要使用流固耦合分析方法。

二、流固耦合的基本理论1. 流动力学基本原理流体的运动可以通过质量守恒、动量守恒和能量守恒三个基本方程来描述。

在流固耦合分析中，需要将流体的力学行为与固体结构的行为相结合，建立适当的数学模型。

2. 固体力学基本原理固体力学是研究物质的应力与应变关系的学科，包括弹性力学、塑性力学、断裂力学等。

在流固耦合分析中，需要考虑固体结构在流体流动中的变形和应力响应。

3. 耦合理论的建立流固耦合分析中，需要将流体力学和固体力学的方程进行耦合，建立相应的耦合模型。

常见的耦合模型有一体化耦合模型和分离式耦合模型。

三、流固耦合的数值计算方法1. 有限元方法有限元方法是流固耦合分析中常用的数值计算方法之一。

通过将连续的物理域离散化为有限个单元，求解各单元上的位移、速度、压力等参数，并通过单元之间的边界条件进行耦合。

2. 边界元方法边界元方法是一种将物理问题转化为表面积分方程求解的数值方法。

在流固耦合分析中，边界元方法可以用于求解流体域和结构域的界面上的力学参数。

3. 网格less方法网格less方法是一种将物理问题空间离散化的方法。

在流固耦合分析中，可以利用网格less方法进行流体和结构的耦合求解。

四、流固耦合的应用1. 结构物的水动力问题流体对结构物的作用会导致结构物的振动、变形和破坏。

流固耦合分析可以用于预测结构物在水动力环境下的响应，并进行可靠性评估。

2. 水力机械的设计与优化水力机械的设计与优化需要考虑流体流动对结构的影响。

通过流固耦合分析，可以优化水力机械的结构形式和流道参数，提高水力机械的工作效率和安全性。

第五章 轴流泵的流固耦合5-1 流固耦合概论流固耦合问题一般分为两类，一类是流‐固单向耦合，一类是流‐固双向耦合。

单向耦合应用于流场对固体作用后，固体变形不大，即流场的边界形貌改变很小，不影响流场分布的，可以使用流固单向耦合。

先计算出流场分布，然后将其中的关键参数作为载荷加载到固体结构上。

典型应用比如小型飞机按刚性体设计的机翼，机翼有明显的应力受载，但是形变很小，对绕流不产生影响。

当固体结构变形比较大，导致流场的边界形貌发生改变后，流场分布会有明显变化时，单向耦合显然是不合适的，因此需要考虑固体变形对流场的影响，即双向耦合。

比如大型客机的机翼，上下跳动量可以达到5 米，以及一切机翼的气动弹性问题，都是因为两者相互影响产生的。

因此在解决这类问题时，需要进行流固双向耦合计算。

下面简单介绍其理论基础。

连续流体介质运动是由经典力学和动力学控制的，在固定产考坐标系下，它们可以被表达为质量、动量守恒形式:()0v tρρ∂+∇⋅=∂ (1) ()B v vv f tρρτ∂+∇⋅-=∂ (2) 式中，ρ为流体密度；v 为速度向量；Bf 流体介质的体力向量；τ为应力张量；在旋转的参考坐标系下，控制方程变为： ()0r v v tρρ∂+∇⋅=∂ （3） (-)+B r r c v v v f f tρρτ∂+∇⋅=∂ （4） 形式和固定坐标系下基本相同，只是速度变成了相对速度，另外就是增加了附加力项c f 。

固体有限元动力控制方程为：[]{}[]{}{}...[]{}M u C u K u F ++= （5）式中，[]M ，[]C ，[]K 分别是质量矩阵，阻尼矩阵以及刚度矩阵，{}F 为载荷矩阵。

流固耦合遵循最基本的守恒原则，所以在流固耦合交界面处，应满足流体与固体应力、位移、热流量、温度等变量的相等或守恒，即满足如下四方程：f f s s n n ττ⋅=⋅ (6)f s d d = (7)f s q q = (8)f s T T = (9)5-2 单向流固耦合思路分析：轴流泵的单向流固耦合仅仅考虑流场对结构的影响，并不考虑结构变形对流场的影响，所以其数据的传递是单向的，流场和结构的分开计算，完成流场计算之后将其作为结构的边界条件加载到结构域上。