流固耦合模拟分析之探讨.docx

MSC.Dytran 流固耦合模拟分析之探讨 摘要：流场会驱动固体；而固体的运动也可能驱动流体，甚至引起流场振动。这称为流固互制或流固耦合。MSC.Dytran软件的流固耦合功能包含一般耦合、任意偶合，且采用拉格朗日法与欧拉法分别描述固体与流体的运动。Dytran中流固耦合计算种类拉格朗日的元素节点依附在材料上，节点随着材料质点作运动，故各物理量也作用在节点上随材料流动而变化。相反，除任意耦合外，欧拉元素网格与节点不随时间而变，其物理量虽也作用在欧拉元素节点上，但对于通过欧拉元素面的各时间的质量、动量与能量的进与出，加之模拟，即模拟元素面的材料流，而不模拟各材料质点的时间历程。因为对一般固体材料，要模拟各材料质点的时间历程，因此大多用拉格朗日法。而对于流体不需要模拟材料质点的时间历程，故采用欧拉法，MSC.Dytran的欧拉法需用三维的计算域、三维的体元素与DMAT通用材料。此外，欧拉法容许一个元素内含有两种以上的材料，这就是模拟计算材料流的扩散与混合行为欧拉法模拟计算材料流的扩散与混合行为。欲推广应用型的计算软件，需有充分的应用范例。关于流固耦合的模拟，除MSC.Dytran中文范例手册（1999b）与Example Problem Manual(2001a)论述到一般耦合与任意耦合的应用范例之外，本文进一步探讨其他的应用。而搭配的MSC.Patran软件除有Results的后处理工具外；该软件的Insight工具，能透明地看到体元素所构成的欧拉域内部，因此，更需用Insight，去展示欧拉域内部的流固耦合计算结果。 Dytran中流固耦合计算种类 欧拉法模拟计算材料流的扩散与混合行为 前言 MSC.Dytran之流固耦合计算功能大致上包括一般耦合(general coupling)和任意耦合(ALE, Arbitrary Lagrangian-Eulerian)两大类。欲推动一项泛用型计算软件被广泛地应用，须有可供参考的文件及丰富的应用案例。而MSC.Dytran之中范例手册包括一个一般耦合和两个任意耦合的应用范例。MSC.Dytran 之Example Problem Manual (2001)也包括两个一般耦合、两个任意耦合的范例。 本文之要旨是经由案例，探讨MSC.Dytran软件对流固耦合（欧拉域）与固体（拉格朗日域）的一般耦合与任意耦合的模拟计算，以及MSC.Dytran软件对流固耦合计算之前后处理与模型结果展示。本文所叙述的流固耦合计算应用场合如下： ? 水下之固体物的告诉移动? 造波板与之波浪水槽 ? 海面上之高速物体撞击混凝土墙? 上游侧与下游侧水深相等情况之潜堰? 上游侧与下游侧水深不等情况之固定堰? 上游侧与下游侧水深不等情况之阶梯式渠道? 固定开度之水利闸门? 隧洞内之气爆压作用在钢板上? 隧洞内之水流推动固体物 对运动的数学描述法 拉格朗日法与欧拉法是对运动现象的两类不同的数学描述，可说是分别对材料质点流与空间流之描述。拉格朗日法与欧拉法之元素网格可在同一计算模型内，但拉格朗日法的元素与欧拉法的元素分别拥有节点，只采用介面(interface)，称为耦合面，才能将两者连结在一起；否则，纵使两者在空间内相互重叠，也彼此不相干，即忽视对方之存在。 2.1 拉格朗日法 对固定的坐标系而言，拉格朗日元素的节点可相对地运动。因节点系附在材料上，故材料连续体之节点系一起随着材料质点流而运动。各拉格朗日元素的质量是不变量(invariant)，但其元素体积可随时间而改变。此外，速度、压力强度或质量密度等物理量系作用在拉格朗日元素的节点上，因此，各物理量系随着材料流(material flow)而改变。因对固体材料之行为， 较须追踪各材料质点之时间历程，故适宜采用拉格朗日法。拉格朗日法也适宜用以分析材料破坏(failure)或应变硬化(strain hardening)问题。 2.2欧拉法 除任意耦合(ALE)外，所有欧拉元素的网格与节点均保持固定，不随着时间或其他运动或变形。换言之，欧拉元素不随时间而变。各时间之速度、 压力强度或密度等物理量也是作用在欧拉节点上，系计算经欧拉元素面的各时间的质量、动量与能量等的进与出之量，而不追踪各材料质点之时间历程。 因对流体较不须追踪各材料质点之时间历程，故一般系采用欧拉法。MSC.Dytran之欧拉法必定使用三维的计算域及三维的体元素(solid elements)，且限于通用材料 DMAT。欧拉法之特点是：须采用较大之计算域，计算结果才不会受到计算域的边界之影响.欧拉法的特点和要求 欧拉法的特点和要求 由于欧拉法系仿真经过欧拉元素面的材料流；且在一个元素之内，容许两种以上之材料，因此，应用欧拉法可模拟计算空气或水等材料流之扩散与混合现象。欧拉法计算混合现象