流体中两弹性板自振荡运动的数值模拟.docVIP

流体中两弹性板自振荡运动的数值模拟 余钊圣，邵雪明 浙江大学航空航天学院力学系，杭州，310027, yuzhaosheng@zju.edu.cn 摘要：本文采用分布式拉格朗日乘子/虚拟区域(DLM/FD)方法，对两平行放置的弹性板在绕流时的自振荡运动进行了数值模拟。方法考虑了弹性板的体积效应和非线性效应。计算结果证实了两弹性板在间距较小时作同相位振荡，在间距较大时作反相位振荡。我们没有发现明显的过渡区，一个很小的间距增加量可导致振荡从同相位到反相位，不过，在刚进入反相位振荡的间距区时，两弹性板的振荡关于中线（面）是不对称的，相应地尾涡结构也相对不规则，因而可以称为‘拟过渡区’。同相位振荡时的振幅比反相位振荡时大，除‘拟过渡区’，反相位振荡时的振幅对于间距的改变不敏感。同相振荡频率小于单板（独立）振荡的频率，而反相振荡频率大于单板振荡的频率。从同相振荡过渡到反相振荡时，频率有较大地跳跃增长，然后随着间距增加，频率逐渐减小，直至单板振荡频率。 关键词：弹性板，流固耦合，自振荡，虚拟区域方法 1. 引言 流体与弹性体的相互作用现象广泛存在于人体生理活动、鱼类游动、鸟类飞行以及工业生产中。目前适用于流体弹性体相互作用的数值方法可分为两大类：一类是基于贴体网格的任意拉格朗日－欧拉有限元方法（ALE）[1,2]，该类方法的特点是流场网格随着弹性体的运动而运动；另一类则是基于非贴体网格的方法[3-9]，其特点是固体内部也作为流场的一部分，因而流场网格是固定的，不受弹性体运动影响。贴体网格方法的优点是计算精度比较高，但算法复杂，计算效率不高。非贴体网格方法由于流场网格是固定的，算法简单，计算效率高，特别是很容易实现高效的并行计算，因而尤其适合复杂的大问题。 非贴体网格方法大致可分为三种：基于非体力的方法[8,9]，基于拟体力为非拉格朗日乘子的浸没边界方法(Immersed Boundary Method)[3-5] 和基于拟体力为拉格朗日乘子的虚拟区域方法[6,7](Distributed Lagrange Multiplier /Fictitious Domain Method)。目前基于拉格朗日乘子的针对流体弹性体相互作用的虚拟区域方法有两个版本：一个是Baaijens 在2001 年根据经验提出的[7]，另一个是余钊圣于2005 年通过推广Glowinski 等人提出的针对流体刚性体相互作用的虚拟区域方法[10]得到的[6]，两个版本的主要区别在于前者的拉格朗日乘子只定义在弹性体边界上，而后者的定义在弹性体内部及边界上。Baaijens 的方法已被成功的应用于心脏瓣膜的运动，但还没有被应用到其它问题。余钊圣提出的虚拟区域方法和其它方法相比具有以下优点：１）该方法具有严格的数学基础，而且通用性好，比如可应用于用任何本构方程描述的流体和固体。２） 在Farnell 等人[8,9]的非体力方法和Peskin 的浸没边界法中，只考虑没有体积的线弹性材料，但实际弹性体材料在大变形时一般都是非线性的，余钊圣的虚拟区域方法考虑了具有体积的不可压缩非线性弹性材料；３）在Baaijens 的虚拟区域方法中，为了避免数值失稳，所有方程耦合在一起求解，因此计算量非常大，余钊圣的虚拟区域方法成功地把原问题分解成几个子问题，大大简化了计算，提高了计算效率。余钊圣对其提出的虚拟区域方法已在两维情况下做了数值实现[6]，并成功模拟了一弹性板在振荡流及绕流中与流体的相互作用。 本文采用余钊圣的虚拟区域方法来研究两块平行放置的弹性板在流体中的自振荡。Zhang等人[11]用实验方法分别研究了一根和两根弹性丝在肥皂膜中的自振荡。在一根弹性丝情形中，在流速超过一定值时，弹性丝可以自发作振荡运动。在两根弹性丝的实验中，发现在两弹性丝间距较小时，它们作同相位振荡，当间距较大时，作反相位振荡。Zhu 和 peskin 用浸没边界法模拟了一弹性丝[4]和两弹性丝的自振荡[5], Farnell等人用非体力法也模拟了一弹性丝[8]和两弹性丝的自振荡[9]。但如前所述，这两个方法都无法考虑材料的体积效应和非线性效应，而且给出的结果也不多。 2. 数值方法 我们假定流体为不可压缩牛顿流体，固体为Neo-Hookean非线性材料。用表示固体材料占据的区域，用表示包含流体和固体的整个区域。弱解形式的无量纲控制方程如下[6]： （1） 国家自然科学基金资助项目(No.. （2） （3） （4） 以上方程中：分别为流体速度，流体压力，拉