薄膜结构风致耦合作用数值初探.pdfVIP

薄膜结构风致耦合作用数值初探 朱伟亮杨庆山 (jE京交通大学土木建筑工程学院，jE京100044) 提要本文采用数值方法，对典型马鞍面封闭式薄膜结构的风致耦合作用进行了模拟研究，并对结构风致振动特性以及流场 对结构振动的影响进行探讨。数值结果表明，结构的平均位移与结构形式、边界条件及风荷载分布有关；而脉动位移则与风 场的涡形、演化规律以及结构形式的关系明显。由于空气与结构的相互作用，薄膜上由于气流作用出现“行波”。在低预应 力条件下，薄膜上各部分的振动特性有明显区别，气流分离、漩涡脱落强烈以及“行波”尾端鞭梢效应明显的区域，脉动位 移幅值及其振动频率较大。 关键词薄膜结构；流固耦合；CFD；有限元 一、引言 膜结构重量轻，刚度小，在风荷载作用下会产生较大位移、速度及加速度，并有可能显著改变结构周 围空气的绕流情况。因此，分析其风致动力效应时，应当考虑流固耦合的作用，这也一直是众多学者努力 的方向。 目前，研究膜结构与风环境流固耦合作用(FluidStructure 武岳等(2003)采用频域方法，生成符合给定功率谱的随机风速时程，根据膜结构体形系数将其转化为外 力施加于膜结构，并在时域内求解结构动力方程，求解过程中通过风与结构的相对速度不断修正风荷载， 以此来考虑膜结构与风环境的“耦合”作用。杨庆山等(2001)提出了一种简化的气弹模型，在结构运动 方程中，通过引入耦合参数(附加质量、气动阻尼、气承刚度等)修正结构运动方程，来考虑膜结构与风 Fluid Structure Dynamics)和CSD(ComputationalDynamics)，建立一种数值风洞模型，通过中间数据交换 平台传递CFD与CSD耦合数据，间接地实现了膜结构与风环境的数值耦合作用。随着CFD理论及计算机 硬件的发展，该方法无疑具有良好的研究前景，正在为众多研究人员所认识，成为一种前沿的研究领域。 本文采用数值方法，对典型马鞍面封闭式薄膜结构的风致耦合作用进行了数值模拟，并对结构风致振 动特性以及流场对结构振动的影响进行探讨。 二、数值模拟方法 一般而言，薄膜结构风致耦合作用的数值模拟方法分为直接耦合与迭代耦合。本文采用迭代耦合方法 作为数值研究的手段。一个迭代耦合下FSI求解系统有三个主要部分组成：流体求解系统，结构求解系 统和一个用于以上两个求解系统的信息交换平台。采用迭代算法直至整个耦合方程的达到收敛解为止。 (一)流体域 在流固耦合作用问题中，耦合界面是可变形的，一般用ALE坐标系描述流体流动，流体流动的未知 量既包括通常的流体变量(压力、速度)，也包括结构运动引起的网格运动速度。不可压缩流体ALE坐标 321 下的控制方程为： 连续方程 u．f=0 动量方程： p(鲁+(_一叶)廿锄w (2) 本构方程： 互口=一p6口+2At％ (3) 式中：p为流体密度；V，w分别为薯方向流体流动速度与网格运动速度；q，为应力张量；．f丑为体力； 1， 、 P为压力；岛为Krenecker符号；At为流体粘度；勺=寺【M√+匕，f)为速度应变张量。 二 利用上述流体控制方程求解流体流动的算法称为直接模拟(DirectNumerical 解高雷诺数湍流流动需要十分强大的计算资源，目前阶段较难实行。因此有必要引入相关的湍流模型。湍 simulation) 流模型主要有时间平均的RANS(ReynoldsaveragedNS)模型和空间平均的LES(Largeeddy 模型。而薄膜结构体型的多样性与