边界元法在膜结构和风耦合的研究中的应用.pdfVIP

边界元法在膜结构与风耦合研究中的应用 徐闻 叶继红 单建 （东南大学土木工程学院，南京 210096 ） 摘 要：膜结构是一种风敏感性结构，结构与风的耦合作用不容忽视。本文将膜结构所处的风场简化为不可压缩势流，运用边界元 法求解绕流后作用在膜结构上的风压，然后计算在风压作用下膜结构产生的变形，修正膜结构形状，再次求解变形后膜结构周围流 场及作用在其上的风压，这样依次迭代，直至得到收敛的解。这种方法的有效性通过典型算例和其他方法计算结果的对比得到验证。 变换膜结构的矢跨比、初始预张力及风场平均风速，研究参数变化对流固耦合的影响，为膜结构抗风研究与工程实践提供参考。 关键词：张拉膜结构，流固耦合，边界元法 中图分类号：TU 393.3,O241,O35 文献标识码：A 严格满足相似比准则且风洞试验技术的难度也较 0 引 言 大，因此膜结构的气弹模型试验较少，有重要参考 近年来，膜结构在大跨度空间结构应用日益广 价值的文献则更少。对于张拉膜结构风洞试验相关 泛，这与它质地轻柔、经济、建造安装简易、抗震 文献[3][4][5] 中对膜结构的振动形态和可能引起气弹 防火性能优越等特点密不可分。然而，薄膜结构对 失稳的原因还存在不同的认识，尚未得到可以明确 风的作用十分敏感，风荷载是该类结构设计中的主 地用于实际结构设计的规律。 要控制荷载。已有大型膜结构不乏风毁实例：1995 简化气弹性模型的方法是通过引入一些必要 年，佐治亚穹顶在强风大雨作用下四片膜材被撕 假设对气弹力学模型进行简化，引入附加气动项， 裂；韩国为 2002 年世界杯修建的济洲岛体育馆， 叠加到静止结构表面的风荷载上以考虑结构与气 其顶部膜结构屋盖先后两次在台风下撕裂；2005 流的耦合作用，然后在时域内利用非线性有限元方 年，台风维达在我国登录，三亚的标志性建筑“美 法求解运动方程，最终得到结构的风振响应。关于 丽之冠”顶部严重撕裂，而博鳌亚洲论坛主会场的 附加气动项的描述，一般常用附加质量、气动阻尼 膜材在台风过后所剩无几。与大跨桥梁、高层建筑 和气动刚度。其中气动刚度一般在充气膜结构中才 及高耸结构相比，膜结构的风振响应研究相对滞 予以考虑，在张拉膜结构中可以忽略。附加质量可 [6][7][8] 。膜结构的气动阻尼很 后，尚未形成一套成熟的方法来满足工程设计和使 以通过薄翼理论来确定 用需要。 难确定，文献[3]、[9]以声致阻尼的形式来处理气动 膜结构通常为典型的大变形柔性体系，不具有 阻尼。而利用而这些处理方法还缺乏理论和试验的 弯曲刚度，其刚度主要由预张力和负高斯曲面所构 比较验证，并且不能在根本上解释膜结构风致振动 成的几何刚度来提供。在风载作用下会产生较大的 气动弹性效应的机制[2] 。文献[10]提出了气弹模型 变形和振动，而膜结构几何形状的改变和振动会影 风洞试验获得附加质量与气动阻尼的方法。文献[11] 响周围流场的变化，进而改变作用在薄膜结构表面 借鉴这一方法由风洞试验获得附加质量与气动阻 的风压，形成风与结构的流固耦合效应（气动弹 尼。 性）。风－结构的流固耦合效应估算一直是结构风 数值风洞也是一个极有前景的方法。气动弹性 工程研究的热点问题之一。由于其具有较大的理论 的数值模