薄膜结构流固耦合效应的简化数值模拟方法.pdfVIP

薄膜结构流固耦合效应的简化数值模拟方法 孙晓颖武岳沈世钊 (哈尔滨工业大学土木工程学院 哈尔滨150090) 摘要：针对薄膜结构流固耦合问题，提出了一种更为实用的简化数值模拟方法。该方法的基本思想是，先将结构响应分 为平均响应、背景响应和共振响应三部分，再针对这三种响应分量的不同性质，采用不同的求解方法，以达到尽量减少CFD 运算时间的目的．应用上述简化数值模拟方法对单向柔性屋盖和鞍形膜结构屋盖进行了分析，并将计算结果与直接数值模拟 方法对比，验证了简化数值模拟方法的精确性和高效性。 。 关键词：薄膜结构；风振分析；流固耦合；CFD数值模拟；耦合影响因子 1引言 薄膜结构是典型的轻质、柔性结构体系，对脉动风荷载的作用十分敏感。膜结构在风荷载的作用下通 常会产生较大的变形和振动，这种大幅的变形和振动反过来也会影响到其表面风压分布，产生所谓的“流 固耦合”效应。因此，要对膜结构的风荷载以及风振响应做出准确估算，除了借助常规分析手段外，还需 要特别关注流固耦合的影响。但是由于膜结构的外形和自身力学性能的复杂性，特别是几何非线性的影响， 使得对膜结构流固耦合问题的研究异常复杂。 薄膜结构在风荷载作用下的耦合振动问题在理论上可归结为，不可压缩粘性流体与几何非线性弹性体 之间的非定常耦联振动问题。对这一问题的求解包括流体域、结构域和网格域三个计算模块：(1)流体域 主要是模拟近地面大气边界层风场，属于钝体空气动力学范畴；(2)结构域主要是模拟张拉膜结构的风致 动力响应，属于几何非线性弹性体的大位移、小应变受迫振动问题；(3)网格域主要是以ALE技术为基础 的动态网格计算问题以及流体一结构网格之间的数据传递问题。 基于上述求解策略，武岳【l’2】建立了适用于索膜结构流固耦合风振分析的数值风洞方澍31。应用该方法， 实现了对二维悬臂板和单向柔性屋盖的流固耦合效应模拟。在这一理论框架下，孙晓颖在文献[4]中综合 运用计算流体力学软件和自行编制的膜结构动力分析程序，搭建了薄膜结构流固耦合的数值模拟平台，并 基于该软件平台实现了对三维结构流固耦合风致振动的数值模拟。这种方法能给出耦合振动的具体过程， 从而得到此过程中的一些统计信息，如均值、方差等，但是这种方法的计算量十分巨大，在微型机上求解 一道题往往需要数十天时间，无法在工程实践中进一步推广。从工程角度看，人们所关心的不是耦合振动 的具体过程，而仅仅是此过程中的一些统计信息。因此不妨考虑换个思路，即放弃对流固耦合具体过程的 模拟，转而采用某种更为简化的方法，仅通过少量计算来确定此过程的统计信息；当然，这样做的前提是 募给出关于薄膜结构流固耦合效应机理的合理解释(或假定)。为此，本文首先基于对薄膜结构流固耦合效 应机理的探讨，提出了一种更具可操作性的简化数值模拟方法，并且对该方法中的若干具体问题进行了详 细探讨。 326 2数值计算方法 根据Davenport提出的结构风振理论协川，可将结构风振响应分为：平均响应，背景响应和共振响应三 部分，如图1所示。其中，平均响应F是结构在平均风荷载作用下的响应；背景响应白是由那些低于结构 基频的脉动风分量引起的，本质上是准静态的，不受结构动力特性的影响；共振响应咯是由接近结构基频 的脉动风分量引起的，具有明显的动力放大效应。这里假定薄膜结构的风振响应亦具有上述特征，则针对 不同响应分量的性质，可以考虑采用不同的求解方法。 在平均响应部分，主要考虑由于结构形状变化所导致的体型系数变化，这是一个静态过程；在背景响 应部分，主要考虑脉动风中的低频部分与结构之间的耦合作用，体现了脉动风的空间相关性对结构整体振 动的影响。这是一个慢变的平稳过程，可以看作是拟静态的；在共振响应部分，主要考虑脉动风中接近结 构基频的部分与结构之间的动力耦合作用。 划瞰f l lag．fi-euqueucy 图1结构风激动力响应的示意图 根据前面所述，可以将风与薄膜结构之间的流固耦合过程分为：静态耦合、拟静态耦合和瞬态耦合三 部分(示意图见图2)。其中，平均响应F属于静态耦合，背景响应瑶属于拟静态耦合，共振响应瑶属于瞬态 耦合。对于