利用瑞利里兹方法计算固有频率的程序_虎门大桥异常抖动背后：流固耦合仿真方法解析...

原创 wei_lewis 仿真秀

导读：5月5日，广东虎门大桥发生异常抖动，全桥路段已实施双向全封闭，禁止通行。网友调侃，马上就要收费了，大桥激动地不能自已。

一、虎门大桥突发异常晃动

作为虎门大桥的设计者之一，西南交通大学土木工程学院教授郑凯锋，通过网络表达了自己的看法。他认为此次异常振动，原因在于近期桥梁维修施工，桥面两侧靠近护栏处摆放了隔离带(俗称“水马”)，梭形风嘴的流线型梁体多了突兀的流场边界(一排水马)，特定风场作用下，桥梁产生涡激振动。

前不久，我在仿真秀官网发布《流固耦合专题之一：什么是流固耦合》原创文章(点击文尾阅读原文可以查看我的专栏)谈到过：流固耦合影响工程安全的标志性事件，即Tacoma大桥在低风速下(19米/秒)下发生剧烈振动并倒塌的事故，这一事件使得颤振成为桥梁风工程的重要课题。而此次虎门大桥的涡振，使得流固耦合问题再次得到大家广泛的关注，甚至成为当天CCTV的热点新闻事件。

热闹之余，作为一名CAEer心中肯定会有疑问，这类抖动需如何仿真分析？

本文将从自然风特性、结构自身的动力特性和风与结构的相互作用三个方面展开介绍风对结构的作用这一复杂的流固耦合现象。

二、自然风特性

风是由于大气中热力和动力现象的时空不均匀性，使相同高度上的两点之间产生压差所造成的。下图为实时更新的地球风场在线可视化网站(https://earth.nullschool.net/)提供的东亚半球上的地球瞬时风场。

在靠近地表的某一范围内，由于摩阻力的存在，使风速的平均值会随着高度的降低而减小，当接近地表时风速接近于零。当高度在1000~2000m之间时，地表摩阻力的影响基本消失，这时风速趋近于常数，在该高度以下即为大气边界层，大气边界层的厚度由风力、地形粗糙度以及纬度而定。

边界层内的风速随高度增大而增大，边界层顶的风速值称为梯度风速。

在边界层内，风受到结构物的阻碍时，其部分动能将转化为作用在结构物上的力，这个力就是我们通常所说的风荷载。

结构物上的风荷载，随空间位置变化，同时受时间历程的影响，因此它是非定常的随机荷载。

大量的实测数据表明，风速时程曲线中包含两种成分，即周期在10min以上的长周期成分和周期仅有几秒钟的短周期成分。因此&#x