基于CFD仿真的平板非线性气动力系统特征研究摘要：为研究平板气动力系统的非线性特征，基于非定常雷诺时均NavierStokes（RANS）方程和SSTkω湍流模型，数值模拟了在单位位移激励下平板非定常运动的绕流场，获得了作用在平板上的气动力时程，并基于Volterra理论开展了平板非线性气动力系统识别.研究表明，本文建立的平板非线性气动力模型能对一定?l率带宽和一定幅值范围的激励产生合理的响应；在本文研究的强迫运动位移幅值和频率范围内，平板非线性气动力模型响应没有表现出对振动幅值和频率的明显相关性，且其气动力的非线性效应并不明显，因而可以认为小攻角下的平板绕流属于气动力弱非线性系统.本文研究证明了CFD模拟在桥梁主梁气动力系统识别上的明显优势.关键词：气动力非线性；平板；Volterra理论；CFD：U448.213文献标志码：A不断突破桥梁跨径记录的跨江跨海工程需求和新材料、新结构形式的应用，可能会使得大跨度桥梁的频率和阻尼不断降低，因而导致结构或其构件对风的敏感程度大大增加.合理评价大跨度桥梁气动性能，准确预测桥梁风致振动稳定性，对保障大跨度桥梁设计合理和运营安全具有深远意义.然而，现有桥梁抗风研究涉及的风荷载，均是建立在气动力的线性模型之上.以气弹分析广泛采用Scanlan自激力模型为例[1]，对二自由度桥梁断面而言，其气动自激力表达为位移、速度与颤振导数的线形组合形式，也即气弹自激力是通过线性叠加多项而得到，其中每一项与位移或速度有关，实际反映的是不同广义位移对应的相对攻角所产生的效应，这是一种典型的线性颤振理论.因此，Scanlan这一桥梁颤振理论也是基于线性和小扰动理论，没有考虑桥梁气动外形对气流场可能产生的大扰动，以及大跨度桥梁大幅运动对气流场的强烈干扰，如流动分离、再附和涡系的相互作用.这些流动现象是桥梁气流场的非线性行为，将导致气动力响应呈现非线性特性，这在Scanlan的早期研究中就有阐述，后续研究者也有类似报导[2].流动非线性对大跨度桥梁气弹响应和颤振稳定性的影响目前并不清楚.因此，对气动外形明显钝化的桥梁主梁，或大幅度振动的桥梁结构，研究气动力的非线性，有助于给出大跨度桥梁风致振动稳定性的合理评价.3.2计算域网格及流动条件包含平板主梁模型的计算域及分区如图4所示.采用圆形计算域，计算域外圆到平面中心的距离均为16B，对应的模型堵塞度不到0.2%，因此无需考虑模型堵塞度带来的影响.采用计算域分区划分以控制网格的正交性和网格缩放比.在平板模型外第一个分区为椭圆形“刚性网格区”，CFD模拟模型运动时，该区域网格与模型刚性固定，在每一时间步上与模型同步运动.该刚性网格区外边界为椭圆，对该区域作四边形结构网格划分.计算域的绝大部分区域为“静止网格区”，该区域外边界是计算域外边界，内边界为离开刚性网格区外椭圆一定距离，且包围刚性网格区的圆形.静止网格区采用四边形单元剖分，从内到外采用合适的网格放大比例.静止网格区和刚性网格区在CFD模拟过程中一直使用计算开始时的网格系统.在“静止网格区”和“刚性网格区”之间为“动网格区”，“动网格区”采用三角形单元剖分.在每一时间步上，“动网格区”根据桥梁断面的运动位置并由设定的网格系统质量要求可能需重新进行网格划分.紧靠桥梁断面的区域流场变量变化剧烈，特别是断面迎风侧和断面法向，因此网格划分适应流场变量的变化程度，并沿各个方向采用适度的网格放大率，实现与动网格区域网格的平顺过渡，如图5所示.计算域的网格总数为96226.边界条件设置为：计算域入口水平均匀速度边界，湍流度为零；下游出口施加流动出口边界条件；模型表面使用无滑移壁面条件；初始场采用入口速度初始化.数值计算采用非定常二阶隐式格式，采用速度压力解耦的SIMPLE算法，二阶格式离散压力方程，动量、湍动能和湍流耗散率方程均采用二阶迎风格式，数值模拟的Re=3.08×105，本文所有数值模拟均基于CFD专用程序Fluent6.3.26开展.3.3单位采样激励的实现在CFD模拟中，连续型Volterra级数无法使用，需采用离散型的Volterra级数形式.对小攻角下平板模型，因其有较好的气动外形，气动力中的非线性可能并不显著.对此类弱非线性气动力系统，...