风洞试验测量覆冰导线气动力特性的数值模拟

风洞试验测量覆冰导线气动力特性的数值模拟

磁线的旋转是信号安全运行中必须解决的难题。舞动问题的研究主要包括舞动机理、数值模拟和覆冰导线气动特性及舞动试验研究等。舞动机理的研究,最早出现的有DenHartog垂直舞动机理和Nigol扭转舞动机理,此外还有多种耦合舞动理论和稳定性理论。近年来,随着计算机硬件和软件的发展,采用数值方法模拟研究舞动的方法倍受关注,并取得了不少的研究成果。覆冰导线的气动特性是诱发舞动的主要因素,因此,对覆冰导线气动特性的风洞试验研究十分重要,国内外已有不少覆冰导线气动特性风洞试验的研究成果。此外,对于各种舞动理论和数值分析方法均需要试验验证。国际上日本和加拿大等国较早建立了试验线路用于观测和研究输电线路的舞动。2008年,我国国家电网公司在河南省郑州尖山建立了我国第一条真型输电线试验线路,并于2011年6月在导线上安装D型人造覆冰模型后诱发了舞动。利用试验线路研究舞动问题存在重复性差、效率低、在现场不易准确获取各种数据的缺点,不便于理论和数值方法的验证以及参数研究等。利用风洞试验模拟舞动是一种途径。文献提出了建立输电导线风洞试验气动弹性模型的方法,并对双分裂导线的气动弹性问题进行了风洞试验研究。然而,对于覆冰导线,要设计满足相似性原理的气动弹性模型却极为困难,特别是要同时满足阻力、升力和扭矩的相似性几乎不可能。因此,基于输电导线大档距小弧垂的特点,文献采用导线节段模型,在风洞中模拟了一种实际冰形覆冰导线的驰振过程。本文针对新月形和扇形两种冰形,首先利用风洞试验测量覆冰导线的气动特性,进而实施覆冰导线节段模型驰振风洞模拟试验,并利用作者提出的舞动数值模拟方法模拟风洞驰振试验,对数值方法进行验证。试验结果对输电线路舞动问题的研究具有重要的价值。1空气动力系数测量1.1覆冰导线测力试验新月形和扇形冰形是输电线路中两种常见的覆冰形式,该两种冰形覆冰导线的截面示意图如图1所示。本文试验中研究的新月形覆冰的厚度分别为12mm和20mm;扇形覆冰厚度取18mm,外圆弧张角为120°,内圆弧张角为140°。试验在中国空气动力研究与发展中心低速空气动力研究所1.4m×1.4m风洞中完成。该风洞为直流式低速风洞,截面形状为切角矩形,试验段长2.8m,风速范围为0～65m/s。使用TG0151A天平测量导线模型的阻力、升力、扭矩三分力。覆冰导线测力试验数据采集系统选用的是PXI系统。角度控制、速压控制由相应的工控机系统实现。设备之间由网络通讯传递指令。1.2攻角范围试验时通过同步转动上下端板改变风向角,每隔5°测试一组升力、阻力和扭矩数据,测试攻角范围为0～180°。12mm厚新月型覆冰导线在不同风速下的升力系数、阻力系数和扭矩系数随风攻角的变化曲线如图3所示,18mm厚扇形覆冰导线在12m/s下的气动系数随风攻角变化曲线如图4所示。2覆冰导线模型对于覆冰导线,要设计同时满足阻力、升力和扭矩相似性的气动弹性模型目前还难以实现。尽管如此,根据输电导线大档距小弧垂的特点,可以采用覆冰导线节段模型,在风洞中模拟其驰振过程,其试验结果可以用于验证覆冰导线舞动方法的正确性。利用节段模型,覆冰导线模型可以采用和原型相同的截面尺寸,即长度缩尺比为1;使单位长度覆冰导线模型的质量与原型一致,则质量缩尺比为1;而导线发生舞动时的风速范围一般为4m/s～20m/s,驰振风洞试验的风速可选择实际风速,即速度缩尺比为1,由此可知其气动力和运动速度的缩尺比也为1。2.1导线节段模型的数字模拟覆冰导线节段模型驰振风洞模拟试验原理参见图5。覆冰导线模型的两端固定在两个竖直圆盘上,导线模型由四根竖直弹簧悬挂起来,另用四根弹簧约束其在水平方向的运动。可通过改变弹簧的刚度和悬挂点调整模型的垂直、水平和扭转刚度,还可通过调整配重块在吊架上的位置改变模型的转动惯量,从而改变系统的固有频率。系统的固有频率可用模态分析系统测得。在导线节段模型的一端部安装三个均匀排列的LED点光源,分别用三个CCD摄像机记录该三个光点的运动轨迹,通过数据采集系统,由计算机记录导线节段模型的垂直、水平位移随时间的变化以及导线模型的驰振运动轨迹。驰振试验仍在1.4m×1.4m风洞中完成。图6所示为试验装置,用于驰振试验的覆冰导线模型长度为1300mm。各竖直弹簧的刚度为40N/m,各水平弹簧的刚度为36N/m。2.2导线风攻角分别对新月形和扇形两种冰形覆冰导线节段模型进行了驰振模拟试验。现分别给出两种冰形情况下的试验结果。对于新月形覆冰导线,由图3可见,当风攻角在40°～100°范围内其升力系数曲线斜率为负,根据DenHartog舞动机理,当升力系数CL和阻力系数CD满足∂CL/∂α+CD<0时可能诱发垂直舞动,故选取初始攻角60°进行试验。测得系统的1阶竖直固有频率和扭转固有频率均为0.81Hz。图7所示为两种厚度新月形覆冰导线节段模型驰振风洞模拟试验结果,图中给出了12mm冰厚18m/s风速和20mm冰厚三种风速下导线的垂直位移时程曲线和驰振运动轨迹。可见,对于相同冰厚的覆冰导线模型,驰振的振幅随风速的增大而增大。此外,各种情况下导线模型的驰振均主要为垂直运动,运动轨迹近似为椭圆,和实际观测的DenHartog垂直舞动现象相符。对于扇形覆冰导线,由图4可见,风攻角在110°～160°范围内其升力系数曲线斜率为负,可能诱发DenHartog垂直舞动,故选取初始攻角120°进行试验。测得系统竖直方向的1阶固有频率和扭转方向的固有频率均为0.82Hz。图8所示为风速为13m/s时扇形覆冰导线节段模型驰振风洞模拟试验结果,图中给出了导线的垂直位移时程曲线和驰振运动轨迹。与新月形覆冰情况类似,运动轨迹仍然近似为椭圆,类似于实际输电线路的垂直舞动现象。3假定误差模型的直线段模型3.1覆冰导线风洞试验模型本文作者提出了一种求解覆冰分裂导线舞动的数值模拟方法,在此利用该方法模拟覆冰导线节段模型风洞试验的驰振过程,进而对该数值方法进行验证。将覆冰导线风洞试验模型简化为如图9所示的数值模型。在数值模型中,用具有扭转自由度的索单元模拟覆冰导线,用空间梁单元模拟水平吊架,弹簧用弹簧单元模拟。利用前述风洞试验测得的气动系数曲线可计算节段模型驰振过程中的气动载荷,气动荷载的计算和驰振的数值模拟方法详见文献。值得一提的是,风洞试验仅给出了典型风速下覆冰导线的气动系数,其它风速下的气动系数可利用试验风速下的测试值插值得到。3.2试验结果对比数值模拟得到的新月形覆冰导线节段模型的驰振结果如图10所示。驰振过程数值模拟中,覆冰导线在风速为18m/s和14m/s的气动系数直接利用风洞试验结果,风速13m/s的气动系数由12m/s和14m/s线性插值得到,15m/s的值由14m/s和18m/s的值插值得到。比较图10和图7的结果可见,数值模拟结果与试验结果一致性较好。数值模拟得到的扇形覆冰导线节段模型的驰振结果如图11所示。比较图11和图8的结果可见,数值模拟得到的垂直位移和运动轨迹均与试验结果吻合。由数值模和试验结果的一致性,验证了本文提出的覆冰导线舞动数值方法的正确性。4模拟分析试验结果本文利用风洞试验测试了两种典型冰形覆冰导线的气动特性,并实现了覆冰导线节段模型驰振风洞模拟试验,得到如下结论:①得到了典型厚度新月形和扇形覆冰导线的升力系数、阻力系数和扭矩系数随风攻角的变化曲线,结果表明,气动系数受风速的影响很小;②在风洞中实现了新月形和扇形覆冰导线节段模型的驰振,得到的驰振运动轨迹符合DenHartog垂直舞动模式;③利用作者研究得到的舞动数值