数值流形方法在流固耦合谐振分析中的应用.pdfVIP

VoL 第24卷第6期 计算力学学报 24，No．6 Chinese of Mechanics December2007 2007年12月 Journal Computational 数值流形方法在流固耦合谐振分析中的应用 苏海东¨“， 黄玉盈1 (1．华中科技大学土木工程与力学学院武汉430074；2．长江科学院材料与结构研究所．武汉430010) 摘 要：数值流形方法(流形法)是石根华博士利用现代数学中流形分析的有限覆盖技术建立起来的新的数值分 析方法。统一解决了连续和非连续变形的力学问题，具有广阔的应用前景。本文将流形法应用于交界面耦合的 流固振动分析．采用平面矩形数学网格。针对无粘、无旋、不可压缩流体和无阻尼的固体结构。提出分析流固桶合 系统筒谐振动的高阶流形法公式，其中．采用拉格朗日乘子法引入流场的已知边界条件。本文还初步研究了在 特殊的无限远沈场中采用解析解覆盖函数的实现技术。文中算例体现了流形法网格划分的方便性和计算的高 精度，显示出流形法在数值解和解析解联合运用上的优势。 关键词：数值流形方法；漉固耦合振动I高阶覆盖函数；解析解覆盖I数值解和解析解的结合 中图分类号：035；0241文献标识码：A 升高时的P型特性有利于在不改变网格的情况下 1 引 言 提高计算精度。 数值流形方法[1](简称流形法)是我国留美学 (3)流场中的某些特殊区域，如无限远流场， 者石根华博士利用现代数学中流形分析的有限覆 采用解析解级数的特殊覆盖函数可以减少数值解 盖技术建立起来的一种新的数值分析方法，该方法 的离散区域和待求解的未知量，提高计算效率。 统一解决了连续和非连续变形的力学问题，具有广 关于上述第(3)点，中科院数学所梁国平教授 阔的应用前景。 在文献[1]的序中指出，在处理无穷区域问题和奇 流形法通过两套相互独立的覆盖体系——反 点问题时，经典的解析法(即级数展开法)往往有很 映数值解精度的数学网格和表示几何边界和材料 高的收敛速度，但一般难以满足复杂求解域上的全 分区的物理网格，将整个研究区域划分成有限个相 部边界条件，因而通常对其收敛域有一定的限制。 互重叠的集合(称为物理覆盖)，在各个覆盖上独立 由此，人们考虑在某些局部区域采用解析法，其他 定义局部覆盖场函数，通过权函数加权平均得到整 区域采用数值方法。可是这样做又会引起新的困 个求解域上的总体场函数。局部覆盖函数可以是 难，由于不同区域待求未知量的性质不同而产生非 多项式或其他类型的级数甚至是解析解级数。物 协调问题，处理起来比较麻烦。例如，文献[4]采用 理覆盖的公共区域称为流形元，相当于有限元法的 Ele