CN104964808A - 一种风雨致输电导线振动的非定常气动力系数计算方法 - Google Patents
一种风雨致输电导线振动的非定常气动力系数计算方法 Download PDFInfo
- Publication number
- CN104964808A CN104964808A CN201510379303.3A CN201510379303A CN104964808A CN 104964808 A CN104964808 A CN 104964808A CN 201510379303 A CN201510379303 A CN 201510379303A CN 104964808 A CN104964808 A CN 104964808A
- Authority
- CN
- China
- Prior art keywords
- vibration
- wind
- rain
- model
- power transmission
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供一种风雨致输电导线振动的非定常气动力系数计算方法,即实验与CFD数值模拟混合子结构方法。包括以下步骤:(1)研制模拟风洞实验台;(2)构建实验用输电导线模型;(3)测试装置获取模型面内外振幅和频率;(4)将实验得到的模型振动响应作为动边界施加到CFD流场子结构;(5)仅对动边界的绕流场进行CFD数值模拟计算得出导线振动的非定常气动力系数。由于该方法不必进行结构的有限元计算和数据的交换,因此可以极大地提高结构风致振动的计算效率;同时,由于结构的振动是物理测试得到的,因此绕流场特性的数值计算结果具有较高的精度。
Description
技术领域
本发明涉及一种计算结构振动的气动力系数的方法,特别是一种计算风雨致输电导线振动非定常气动力系数的方法。
背景技术
输电导线的风雨激振是固体,液体,气体之间的耦合振动。风雨激振具有振幅大,频率低、偶发性等特点,易造成金具剧烈磨损及导线伤股、断股现象,严重影响输电线路的安全稳定运行。因此,对于输电导线发生风雨激振机理的研究是亟待解决的关键问题。为揭示风雨激振诱发机理,必须获得输电导线发生振动的非定常气动力系数。
目前,针对输电风雨激振的研究方法主要包括:现场观测、风洞试验、理论分析以及CFD数值模拟计算。由于在风雨条件下输电导线会发生面内和面外振动,引起风偏角,雨线位置不断变化,因此难以通过风洞实验方法得到风雨致输电导线振动非定常气动力系数的变化规律。采用CFD数值模拟方法研究输电导线风雨激振时可以获得导线上的气动力,但通常情况下,随着雷诺数的增大,网格数和计算量以指数级增长,使计算效率大大降低,甚至对某些问题难以计算得到。
发明内容
针对上述现有技术的不足,本发明提供了一种风雨致输电导线振动的非定常气动力系数计算方法。该方法不必进行结构的有限元计算和数据交换,可以精准,高效的计算出输电导线发生风雨激振时非定常气动力系数。
为实现上述发明目的,本发明采用的技术方案如下:
一种风雨致输电导线振动的非定常气动力系数计算方法,即实验与CFD数值模拟混合子结构方法。包括以下步骤:(1)研制模拟风洞实验台;(2)构建实验用输电导线模型;(3)测试装置获取模型面内外振幅和频率;(4)将实验得到的模型振动响应作为动边界施加到CFD流场子结构;(5)仅对动边界的绕流场进行CFD数值模拟计算得出导线振动的非定常气动力系数。
前述风雨致输电导线振动的非定常气动力系数计算方法中,所述子结构包括输电导线模型子结构、上雨线子结构和流场子结构;其中上雨线子结构截面形状为类抛物线;宽度为8mm,高度为5mm。
前述风雨致输电导线振动的非定常气动力系数计算方法中,步骤(1)中所述的风洞实验台是由斜流变频风机和矩形风洞构成。矩形风洞整体由80×80mm2的方钢管焊接而成;实验台总长4.5m,高2m,占地面积约9m2。通过变频器改变风机的转速,从而模拟不同的风速。
前述风雨致输电导线振动的非定常气动力系数计算方法中,步骤(2)中所述的实验用模型是由输电导线模型和人工雨线模型组成。输电导线模型是用空心铝管模拟输电导线,铝管两端分别焊接两根实心轴,并在实心轴紧靠铝管处安装滚动轴承;在实心轴的两端分别悬挂两组弹簧阻尼系统,弹簧阻尼系统固定于两个矩形框;可实现水平、竖直方向的移动,以实现输电导线倾角和风向角的变化。人工雨线模型是依附于模拟导线外表面上,其两端固联于滚动轴承的外径上。
前述风雨致输电导线振动的非定常气动力系数计算方法中,步骤(3)中所述的测试装置是由加速度传感器、位移传感器以及无线动态信号采集器组成。传感器分别固定于实心轴的两端,所测得的振动信号被其旁的无线动态信号采集器接收。
前述风雨致输电导线振动的非定常气动力系数计算方法中,步骤(4)中所述的动边界设置如下:输电导线边界按照风洞实验实测的振动做横向风振动;上雨线在跟随输电导线模型边界做横向风振动的同时,沿着导线模型边界做环向振荡。
前述风雨致输电导线振动的非定常气动力系数计算方法中,步骤(5)中所述的CFD数值模拟计算是仅对动边界绕流场的计算。
与现有技术相比,本发明的有益效果在于本发明不需要进行结构的有限元计算和数据交换,可以极大地提高结构风致振的计算效率;同时由于结构振动是在风洞实验台中测到的,测到的水线和导线振动响应真实地反映了它们与流场之间的相互作用,因此,绕流场特性的数值计算结果具有较高的精度。
附图说明
图1是本发明的数据处理流程图;
图2是风洞实验台结构示意图;
图3是带有人工雨线的输电导线结构示意图;
图4是CFD模拟结构示意图;
图5是气动力升力系数时程变化曲线;
图6是气动力阻力系数时程变化曲线;
其中,1-变频器,2-斜流式风机,3-矩形洞体,4-空心铝管,5-人工雨线,6-无线动态信号采集器,7-弹簧阻尼系统,8-矩形框,9-轴承,10-实心轴,11-位移传感器,12-加速度传感器,13-来流风,14-上水线子结构,15-输电导线子结构,16-流场子结构,17-面内振动阻尼系数,18-面内振动刚度,19-面外振动刚度,20-面外振动阻尼系数。
具体实施方式
实施例:
一种风雨致输电导线振动的非定常气动力系数计算方法,即实验与CFD数值模拟混合子结构方法。图(1)所示包括以下步骤:(1)研制模拟风洞实验台;(2)构建实验用输电导线模型;(3)测试装置获取模型面内外振幅和频率;(4)将实验得到的模型振动响应作为动边界施加到CFD流场子结构;(5)仅对动边界的绕流场进行CFD数值模拟计算得出导线振动的非定常气动力系数。
具体计算过程如下:
(1)如图(2)所示采用80×80mm2的方钢管焊接成长度为4.5m,高度为2m的矩形风洞体,将斜流变频式风机(1,2)安装在矩形风洞体3的一端,中心高度与输电导线模型4中心高度持平。
(2)如图(3)所示采用长为1.2m,直径为30mm的空心铝管4两端分别焊接两根直径为10mm的实心轴10,并实心轴紧靠铝管处安装两个滚动轴承9,轴承9内圈与实心轴10过盈配合。人工雨线5依附于模拟导线4外表面上,其两端固联于滚动轴承9的外径上。在实心轴10的两端分别悬挂两组弹簧阻尼系统7,弹簧阻尼系统7固定于两个矩形框8。通过矩形框8的水平,竖直方向相对移动,实现不同倾角和风向角下的模拟输电导线气动力分析。
(3)将加速度传感器12和位移传感器11分别固定于实心轴10两端。通过控制变频器1和矩形框8的移动分别实现不同风速,不同风向角下的模拟输电导线和上雨线模型的振动响应;传感器(11,12)将测得的振动响应传给无线信号采集器6。
(4)如图(3)所示将CFD模拟系统分为流场子结构16、输电导线模型子结构15和上雨线子结构14;具体参数:流场子结构16的流场域尺寸为1200×600(mm),网格数250000;导线的CFD模型15以LGJ-500/35导线为原型,导线模型D=30mm,密度为7000kg/m3;网格数为88;上水线14采用宽8mm,高5mm的类抛物线形的截面形状,水线的相位角为60°;导线和雨线模型距离流场入口边界和上边界各为300mm。
(5)将步骤(3)测得的振动响应施加在流场子结构16:输电导线15边界按照实测的振动做横向风振动;上雨线14在跟随输电导线模型边界做横向风振动的同时,沿着导线模型边界做环向振荡。
(6)进行CFD数值模拟计算,计算步数设置为3000;计算得到风雨致输电导线振动的非定常气动力升力系数图(5)和阻力系数图(6)。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受所述实施例的限制,其他的任何在未脱离本发明宗旨的前提下做出的各种变化均属于本发明的保护范围之内。
Claims (6)
1.一种风雨致输电导线振动的非定常气动力系数计算方法,即实验与CFD数值模拟混合子结构方法。包括以下步骤:(1)研制模拟风洞实验台;(2)构建实验用输电导线模型;(3)测试装置获取模型面内外振幅和频率;(4)将实验得到的模型振动响应作为动边界施加到CFD流场子结构;(5)仅对动边界的绕流场进行CFD数值模拟计算得出导线振动的非定常气动力系数。
2.根据权利要求1所述实验与CFD数值模拟混合子结构方法,其特征在于所述子结构包括输电导线模型子结构、上雨线子结构和流场子结构;其中上雨线子结构截面形状为类抛物线。
3.根据权利要求1所述一种风雨致输电导线振动的非定常气动力系数计算方法,其特征在于步骤(1)中风洞实验台是由斜流变频风机和矩形风洞构成。
4.根据权利要求1所述一种风雨致输电导线振动的非定常气动力系数计算方法,其特征在于步骤(2)中所述的实验用模型是由输电导线模型和人工雨线模型组成。
5.根据权利要求1所述一种风雨致输电导线振动的非定常气动力系数计算方法,其特征在于步骤(3)中所述的测试装置是由加速度传感器、位移传感器以及无线动态信号采集器组成。
6.根据权利要求1所述一种风雨致输电导线振动的非定常气动力系数计算方法,其特征在于步骤(4)中所述的动边界设置如下:输电导线边界按照风洞实验实测的振动做横向风振动;上雨线在跟随输电导线模型边界做横向风振动的同时,沿着导线模型边界做环向振荡。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510379303.3A CN104964808A (zh) | 2015-07-02 | 2015-07-02 | 一种风雨致输电导线振动的非定常气动力系数计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510379303.3A CN104964808A (zh) | 2015-07-02 | 2015-07-02 | 一种风雨致输电导线振动的非定常气动力系数计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104964808A true CN104964808A (zh) | 2015-10-07 |
Family
ID=54218858
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510379303.3A Pending CN104964808A (zh) | 2015-07-02 | 2015-07-02 | 一种风雨致输电导线振动的非定常气动力系数计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104964808A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110909512A (zh) * | 2019-11-26 | 2020-03-24 | 国网新疆电力有限公司电力科学研究院 | 一种基于非定常方法的双分裂背风子导线气动荷载仿真方法 |
CN111027261A (zh) * | 2019-11-15 | 2020-04-17 | 四川大学 | 一种用于研究结构风激励响应的混合模拟试验方法 |
CN112464589A (zh) * | 2020-12-11 | 2021-03-09 | 国网新疆电力有限公司电力科学研究院 | 输电导线气动阻力简化数值模拟方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101221553A (zh) * | 2008-01-04 | 2008-07-16 | 东南大学 | 斜风作用下大跨度桥梁抖振响应分析时域方法 |
JP4284688B2 (ja) * | 2005-08-25 | 2009-06-24 | 株式会社フジクラ | 電線防護管の品質評価方法 |
CN102254068A (zh) * | 2010-12-01 | 2011-11-23 | 东南大学 | 大跨度桥梁抖振响应的多尺度分析方法 |
CN104599568A (zh) * | 2014-12-08 | 2015-05-06 | 华北电力大学 | 一种输电导线用风雨致振模拟实验装置 |
CN104596727A (zh) * | 2014-12-08 | 2015-05-06 | 华北电力大学 | 一种带人工雨线的高压输电导线气动力特性实验装置 |
CN104596783A (zh) * | 2014-12-08 | 2015-05-06 | 华北电力大学 | 一种用于模拟下击暴流作用下高压输电导线动力特性的实验装置 |
-
2015
- 2015-07-02 CN CN201510379303.3A patent/CN104964808A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4284688B2 (ja) * | 2005-08-25 | 2009-06-24 | 株式会社フジクラ | 電線防護管の品質評価方法 |
CN101221553A (zh) * | 2008-01-04 | 2008-07-16 | 东南大学 | 斜风作用下大跨度桥梁抖振响应分析时域方法 |
CN102254068A (zh) * | 2010-12-01 | 2011-11-23 | 东南大学 | 大跨度桥梁抖振响应的多尺度分析方法 |
CN104599568A (zh) * | 2014-12-08 | 2015-05-06 | 华北电力大学 | 一种输电导线用风雨致振模拟实验装置 |
CN104596727A (zh) * | 2014-12-08 | 2015-05-06 | 华北电力大学 | 一种带人工雨线的高压输电导线气动力特性实验装置 |
CN104596783A (zh) * | 2014-12-08 | 2015-05-06 | 华北电力大学 | 一种用于模拟下击暴流作用下高压输电导线动力特性的实验装置 |
Non-Patent Citations (1)
Title |
---|
陈文礼 等: "斜拉索风雨振-实验与数值模拟混合子结构方法", 《哈尔滨工业大学学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111027261A (zh) * | 2019-11-15 | 2020-04-17 | 四川大学 | 一种用于研究结构风激励响应的混合模拟试验方法 |
CN111027261B (zh) * | 2019-11-15 | 2023-09-26 | 四川大学 | 一种用于研究结构风激励响应的混合模拟试验方法 |
CN110909512A (zh) * | 2019-11-26 | 2020-03-24 | 国网新疆电力有限公司电力科学研究院 | 一种基于非定常方法的双分裂背风子导线气动荷载仿真方法 |
CN112464589A (zh) * | 2020-12-11 | 2021-03-09 | 国网新疆电力有限公司电力科学研究院 | 输电导线气动阻力简化数值模拟方法 |
CN112464589B (zh) * | 2020-12-11 | 2023-02-03 | 国网新疆电力有限公司电力科学研究院 | 输电导线气动阻力简化数值模拟方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103969022B (zh) | 一种高超声速风洞湍流度间接测量方法 | |
CN202075304U (zh) | 一种超声波风速仪 | |
CN104596727A (zh) | 一种带人工雨线的高压输电导线气动力特性实验装置 | |
CN104964808A (zh) | 一种风雨致输电导线振动的非定常气动力系数计算方法 | |
Lueck | Calculating the rate of dissipation of turbulent kinetic energy | |
CN105004506A (zh) | 自升式海洋平台桩腿波流载荷系数测试实验系统 | |
CN109799053B (zh) | 一种输电设备动力特性分析方法 | |
CN104792493A (zh) | 一种喷流反射速度场测试方法 | |
RU179254U1 (ru) | Электромеханический стенд | |
CN103048153A (zh) | 一种用于圆筒形结构声学性能测试的试验系统 | |
CN209841331U (zh) | 移动模架支撑系统风洞试验装置 | |
CN104732071A (zh) | 一种动量轮与航天器结构的耦合动响应获取方法 | |
RU147970U1 (ru) | Конструкция несущей арматуры электроакустических преобразователей ультразвукового 3d-анемометра | |
CN202928632U (zh) | V型内锥流量计 | |
CN106560685A (zh) | 一种针对于垂直轴风力发电机的翼型动态气动特性试验台 | |
CN205749173U (zh) | 一种减振材料的振动测试装置 | |
CN202075303U (zh) | 一种微差压式高精度测风装置 | |
CN111577516B (zh) | 一种海洋温差能发电平台海水管振动监测系统 | |
CN104807413A (zh) | 自定位道路监测设备 | |
Li et al. | Investigation on passive simulation method and factors influencing the type-C-terrain wind profile of a structural wind-resistant moving-vehicle tester | |
CN104180926B (zh) | 一种喷流反射温度场测试方法 | |
CN203798531U (zh) | 沙粒启动风速观测装置 | |
CN109001038B (zh) | 管片加载试验装置及方法 | |
CN211085467U (zh) | 一种用于带冠叶片测振的标定系统 | |
CN203948229U (zh) | 一种风电机组塔筒倾斜测量装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151007 |
|
RJ01 | Rejection of invention patent application after publication |