VOF法在斜拉索风雨激振数值模拟中应用.docVIP

VOF法在斜拉索风雨激振数值模拟中应用 　　摘要：将气液两相流理论与体积分数法（VOF法）相结合，提出一种模拟发生风雨激振现象时斜拉索表面水线形成和运动的新方法；应用这一方法采用计算流体软件CFX进行数值模拟得到拉索表面水线运动规律和拉索气动力；将气动升力施加到拉索上得到拉索的振动响应；通过与风洞试验结果及混合子结构方法的计算结果进行比对，验证了这一方法的准确性与合理性；通过分析水线运动、拉索气动力变化与拉索振动间的内在联系，揭示风雨激振现象的产生机理：上水线的环向周期性振荡导致拉索升力的周期性变化，引起拉索发生大幅度振动，上水线与拉索间的共振是产生风雨激振现象的主要原因。 　　关键词：风雨激振；斜拉索；水线；水膜；VOF法 　　中图分类号：TU973.2+13；U443.38 　　文献标识码：A 　　文章编号：1004-4523（2015）01-0130-09 　　引言 　　1984年Hikami和Shiraishi在日本的Meiko-Nishi桥上首次观测到斜拉索在一定风雨条件下会发生大幅度振动，即风雨激振现象，并发现水线在斜拉索表面上的形成与振荡是风雨激振现象的特殊标志和决定因素。科研人员首先通过现场观测和风洞实验来研究这一现象的产生原因，并在此基础上提出诸多理论模型。近年来，CFD软件被逐渐应用于风雨激振研究。Li等首先在风洞试验中应用超声波测厚系统对斜拉索风雨激振时表面的水线几何及动力特性进行了定量测量与分析，随后将风洞试验与CFD数值模拟相结合，假设水线在斜拉索表面做环向运动，应用混合子结构方法得到斜拉索风雨激振的气动力时程。毕继红等则将滑移理论与CFD数值模拟相结合，考虑斜拉索表面水膜形态对风压力系数和风摩擦力系数的影响，建立了模拟风雨条件下斜拉索表面水膜形态变化的气、液两相耦合二维断面模型，并推导出风雨作用下斜拉索气动力计算公式，通过数值模拟得到不同风速下的斜拉索表面水膜形态及气动力变化规律；随后又在此基础上，将滑移理论与单质点单自由度振动理论相结合，通过建立耦合的水膜运动方程和斜拉索振动方程，首次提出能够模拟水线形成及运动的气、液、固三相耦合二维断面模型。 　　由于应用已有的理论模型很难求得发生风雨激振现象时不同水膜形态下的拉索气动力，为分析水线对拉索气动力的影响，本文基于气液两相流理论和VOF法提出一种模拟水线产生过程的新方法；采用计算流体软件CFX模拟风场中斜拉索周围的气液自由界面，通过数值计算得到的拉索表面水膜（水线）形态及拉索气动力时程变化，将拉索气动力施加到斜拉索模型上得到斜拉索的振动响应；与已有的风洞试验及数值模拟结果进行比对验证，研究发生风雨激振现象时拉索表面水线运动、拉索气动力变化与拉索振动间的内在联系，初步揭示风雨激振现象的产生机理。 　　1VOF法简介 　　VOF法是一种在固定的欧拉网格下的表面跟踪方法，模型中的两种或多种流体没有互相穿插，可以形成互不相融的流体间界面。针对增加到模型中的每一个附加相，引进一个变量来计算单元里流体相的体积分数，在每个控制体内所有相的体积分数的和为1。如果某一单元内第q相流体的体积分数即为αq，则代表的含义为： 　　应用CFX软件进行数值模拟，为方便计算建模，选取的流场计算区域平行于UN，如图3所示。设拉索直径为D=2R，计算区域取边长为30D×20D的矩形，将拉索中心置于坐标原点，拉索中心距离人口10D，距离出口20D，距离上下边界10D。应用采用ICEM软件中的“O”Type划分流场网格，如图4所示。流场选择连续相气相空气（Air at 25℃）和连续相液相水（Water），水的表面张力系数为0.072N/m；采用均相流模型，气相和液相的切向速度与法向速度相同；流场参考压强为一个标准大气压，在gN方向设置重力加速度；连续相气相和连续相液相流场均采用剪切应力传输模型（SST模型）。 　　3模型验证 　　3.1无气流作用（U=0）时的水膜形状 　　为验证应用气液两相流理论和VOF法模拟水线产生过程这一方法的合理性，首先模拟水膜只受重力作用时的形态变化，如图5所示。为便于观察，将θ=360°置于坐标轴中心。与自然现象及已有的数值模拟结果相一致，拉索上表面处的水膜逐渐变薄，而拉索下表面处的水膜沿重力方向不断积聚的同时，最终形成水线。 　　3.2无重力作用（g=0）时的水膜形状 　　为验证这一方法的准确性，考察水膜只受气流作用时的形态变化。采用CFD软件计算时流场从不稳定状态发展到稳定状态需要一定的时间，本模型经过大约500个时间步（0.025s）的计算后流场已达到稳定状态，如图6所示，拉索附近以外的流场速度均为给定风速U=7.72m/s。 　　图7显示了无重力作用，风速U=7.72m/s时拉索表面水膜的形态变化。在拉索表面形成了