第11章 流固耦合分析.pdfVIP

第十一章第十一章 流固耦合分析流固耦合分析 第十一章第十一章 流固耦合分析流固耦合分析 LS-DYNA 除了具有强大的结构动力分析模块外，还具有强大的流体与固体相互耦合的 功能，广泛应用于各种爆炸 （水下、空中、建筑物和土壤中）、气囊展开、体积成型，罐内 液体晃动等分析中。 下面具体讨论怎样应 LS-DYNA 进行流固耦合分析及应该注意的一些问题，也希望通 过本章的学习对LS-DYNA 流固耦合分析方法有深入的了解和掌握。 11.1 Lagrangian、、Eulerisn 和和 ALE 算法算法 、、 和和 算法算法 在LS-DYNA 中三维单元有三种基本的算法：Lagrangian 、Eulerisn 和ALE （任意拉格 朗日欧拉算法）算法，是由关键字*SECTION_SOLID 中的ELFORM 控制的: 从上面可知ELFORM 的取值不一样导致单元算法的不一样，可用于流固耦合分析的单 元公式如下 （不包括1）： 1 常应力实体单元 （纯粹的Lagrangian 算法） 5 中心单点积分的ALE 算法 （单元内为单一材料） 6 中心单点积分的Eulerisn 单元 （单元内为单一材料） 7 中心单点积分的环境Eulerisn 单元 （用于Eulerisn 计算的进出口边界） 11 中心单点积分的ALE 多物质单元 （一个单元内可以包含多种物质） 12 中心单点积分的带空白材料的单物质ALE 单元 149 Lagrangian、、Eulerisn 和和ALE 算法的区别算法的区别：： 、、 和和 算法的区别算法的区别：： LS-DYNA 程序具有Lagrange 、Euler 和ALE 算法，Lagrange 算法的单元网格附着在材 料上，随着材料的流动而产生单元网格的变形。但是在结构变形过于巨大时，有可能使有限 元网格造成严重畸变，引起数值计算的困难，甚至程序终止运算。 ALE 算法和Euler 算法可以克服单元严重畸变引起的数值计算困难，并实现流体- 固体 耦合的动态分析。ALE 算法先执行一个或几个Lagrange 时步计算，此时单元网格随材料流 动而产生变形，然后执行ALE 时步计算：（1）保持变形后的物体边界条件，对内部单元进 行重分网格，网格的拓扑关系保持不变，称为Smooth Step ；（2 ）将变形网格中的单元变量 （密度、能量、应力张量等）和节点速度矢量输运到重分后的新网格中，称为Advection Step 。 用户可以选择ALE 时步的开始和终止时间，以及其频率。Euler 算法则是材料在一个固定的 网格中流动，在LS-DYNA 中只要将有关实体单元标志Euler 算法，并选择输运 advection 算法。LS-DYNA 还可将Euler 网格与全Lagrange 有限元网格方便地耦合，以处理流体与结 构在各种复杂载荷条件下的相互作用问题。 为了说明问题，以一个2D 的长方形变形为例来说明。对于同一个物理过程，可以用不 同的方式来描述：A. Lagrangian，，B. Eulerisn，，C. ALE，分别以上面三种方式来分别描述 ，， ，， 该物体的变形，分析他们的差别，如图所示： Lagrangian 网格平动网格平动 网格平动网格平动 材料变形材料变形 材料变形材料变形 A Eulerian 网 格 在 空 间固定 B C ALE 网格平动网格平动 网格平动网格平动 实体材料实体材料 空的或空气空的或空气 实体材料实体材料 空的或空气空的或空气 ALE 网网 格格 + 网网 格格 t- dt t 可以移动可以移动 可以移动可以移动 经过一个dt 的时间变化后，比较三种描述的构形变化： A. Lagrangian ：： ：： 对于拉格朗日描述，空间网格的节点与假想的材料点是一致的，也就是说，网格变形， 材料也跟着网格变形，如图中A 所示，所以对于大变形情况，网格可能发生严重畸变。 B. Eulerisn ：： ：： 对于欧拉描述，可以这样理解，即有两层网格重迭在一起，一个是空间网格固定在空间 中不动，另一层附着在材料上随材料在固定的空间网格中流动，并通过下面两步来实现： 首先，材料网格以一个拉格朗日步变形 （象A 所描述的那样），然后拉格朗日单元的状 态变量被映射或输送回到固定的空间网格中去。这样网格总是不动和不变形的，相当于 材料在网格中流动，如图中B 所示，从而可以处理流体流动等大变形问题。 150 C. ALE 对于ALE （任意拉格朗日欧拉）描述，与欧拉描述一样，可以理解有两层网格重迭在一 起，但空间网格可以在空间任意运动，其 与欧拉描述一样，有物质的输送在两层网格 中发生。如图C